
We study the influence of geometric and dynamical anisotropies on the
development of flow harmonics and, simultaneously, on the second and
thirdorder 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 secondorder (or thirdorder) 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 finalstate
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 numberofconstituentquark 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 eventbyevent harmonic flow coefficients measured in PbPb
collisions at centerofmass 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 hydroinspired freezeout 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
centerofmass energy 2.76 TeV per nucleon pair are analyzed and interpreted in
the framework of the HYDJET++ model. The crosstalk 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 midrapidity
in protonproton (pp) 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 overlapeccentricity and
particledensity scalings as found in highenergy heavy ion collisions. Using a
Glauber eikonal model, we compute the pp eccentricities, transverse areas and
particlemultiplicities for various phenomenological parametrisations of the
proton spatial density. For realistic proton transverse profiles, we find
integrated elliptic flow v2 parameters below 3% in pp collisions at sqrt(s) =
14 TeV.

The azimuthal anisotropy of charged particles in heavy ion collisions is an
important probe of quarkgluon 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++ MonteCarlo
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 MonteCarlo event generator for simulation of relativistic
heavy ion AA collisions considered as a superposition of the soft, hydrotype
state and the hard state resulting from multiparton 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
objectoriented C++ language under the ROOT environment. The hard part of
HYDJET++ is identical to the hard part of Fortranwritten 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 freezeout hypersurfaces obtained from the parameterization of
relativistic hydrodynamics with preset freezeout conditions. It includes the
longitudinal, radial and elliptic flow effects and the decays of hadronic
resonances. The corresponding fast MonteCarlo 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 multiparticle 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 hydrotype state and hard multijets.
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 quarkgluon plasma, created initially in nuclear overlap zone in
collisions with nonzero impact parameter. The calculations are performed for
semicentral PbPb collisions at LHC energy.