We study the influence of geometric and dynamical anisotropies on the development of flow harmonics and, simultaneously, on the second- and third-order oscillations of femtoscopy radii. The analysis is done within the Monte Carlo event generator HYDJET++, which was extended to dynamical triangular deformations. It is shown that the merely geometric anisotropy provides the results which anticorrelate with the experimental observations of either $v_2$ (or $v_3$) or second-order (or third-order) oscillations of the femtoscopy radii. Decays of resonances significantly increase the emitting areas but do not change the phases of the radii oscillations. In contrast to the spatial deformations, the dynamical anisotropy alone provides the correct qualitative description of the flow and the femtoscopy observables simultaneously. However, one needs both types of the anisotropy to match quantitatively the experimental data.

Triangular flow $v_3$ of identified and inclusive particles in Pb+Pb collisions at $\sqrt{s_{NN}} = 2.76$ TeV is studied as a function of centrality and transverse momentum within the \textsc{hydjet++} model. The model enables one to investigate the influence of both hard processes and final-state interactions on the harmonics of particle anisotropic flow. Decays of resonances are found to increase the magnitude of the $v_3(p_{\rm T})$ distributions at $p_{\rm T} \geq 2$ GeV/$c$ and shift their maxima to higher transverse momenta. The $p_{\rm T}$-integrated triangular flow, however, becomes slightly weakened for all centralities studied. The resonance decays also modify the spectra towards the number-of-constituent-quark scaling fulfillment for the triangular flow, whereas jets are the main source of the scaling violation at the energies available at the CERN Large Hadron Collider (LHC). Comparison with the corresponding spectra of elliptic flow reveals that resonance decays and jets act in a similar manner on both $v_3(p_{\rm T})$ and $v_2(p_{\rm T})$ behavior. Obtained results are also confronted with the experimental data on differential triangular flow of identified hadrons, ratio $v_3^{1/3}(p_{\rm T}) / v_2^{1/2}(p_{\rm T})$ and $p_{\rm T}$-integrated triangular flow of charged hadrons.

The LHC data on event-by-event harmonic flow coefficients measured in PbPb collisions at center-of-mass energy 2.76 TeV per nucleon pair are analyzed and interpreted within the HYDJET++ model. To compare the model results with the experimental data the unfolding procedure is employed. The essentially dynamical origin of the flow fluctuations in hydro-inspired freeze-out approach has been established. It is shown that the simple modification of the model via introducing the distribution over spatial anisotropy parameters permits HYDJET++ to reproduce both elliptic and triangular flow fluctuations and related to it eccentricity fluctuations of the initial state at the LHC energy.

The LHC data on azimuthal anisotropy harmonics from PbPb collisions at center-of-mass energy 2.76 TeV per nucleon pair are analyzed and interpreted in the framework of the HYDJET++ model. The cross-talk of elliptic $v_2$ and triangular $v_3$ flow in the model generates both even and odd harmonics of higher order. Comparison with the experimental data shows that this mechanism is able to reproduce the $p_{\rm T}$ and centrality dependencies of quadrangular flow $v_4$, and also the basic trends for pentagonal $v_5$ and hexagonal $v_6$ flows.

We estimate the amount of collective "elliptic flow" expected at mid-rapidity in proton-proton (p-p) collisions at the CERN Large Hadron Collider (LHC), assuming that any possible azimuthal anisotropy of the produced hadrons with respect to the plane of the reaction follows the same overlap-eccentricity and particle-density scalings as found in high-energy heavy ion collisions. Using a Glauber eikonal model, we compute the p-p eccentricities, transverse areas and particle-multiplicities for various phenomenological parametrisations of the proton spatial density. For realistic proton transverse profiles, we find integrated elliptic flow v2 parameters below 3% in p-p collisions at sqrt(s) = 14 TeV.

The azimuthal anisotropy of charged particles in heavy ion collisions is an important probe of quark-gluon plasma evolution at early stages. In the present paper the elliptic flow pattern in Pb+Pb collisions at sqrt{s}=5.5 TeV is analyzed for different hadron species in the frameworks of HYDJET++ Monte-Carlo model. The influence of resonance decays on particle flow is investigated. The different methods of elliptic flow reconstruction are compared under LHC conditions.

HYDJET++ is a Monte-Carlo event generator for simulation of relativistic heavy ion AA collisions considered as a superposition of the soft, hydro-type state and the hard state resulting from multi-parton fragmentation. This model is the development and continuation of HYDJET event generator (Lokhtin & Snigirev, 2006, EPJC, 45, 211). The main program is written in the object-oriented C++ language under the ROOT environment. The hard part of HYDJET++ is identical to the hard part of Fortran-written HYDJET and it is included in the generator structure as a separate directory. The soft part of HYDJET++ event is the "thermal" hadronic state generated on the chemical and thermal freeze-out hypersurfaces obtained from the parameterization of relativistic hydrodynamics with preset freeze-out conditions. It includes the longitudinal, radial and elliptic flow effects and the decays of hadronic resonances. The corresponding fast Monte-Carlo simulation procedure, C++ code FAST MC (Amelin et al., 2006, PRC, 74, 064901; 2008, PRC, 77, 014903) is adapted to HYDJET++. It is designed for studying the multi-particle production in a wide energy range of heavy ion experimental facilities: from FAIR and NICA to RHIC and LHC.

The model to simulate rescattering and partonic energy loss in ultrarelativistic heavy ion collisions is presented. The full heavy ion event is obtained as a superposition of a soft hydro-type state and hard multi-jets. This model is capable of reproducing main features of the jet quenching pattern at RHIC, and is applied to probe jet quenching in various novel channels at LHC.

We analyze the azimuthal anisotropy of jet spectra due to energy loss of hard partons in quark-gluon plasma, created initially in nuclear overlap zone in collisions with non-zero impact parameter. The calculations are performed for semi-central Pb-Pb collisions at LHC energy.