HYDJET++ model combines the description of soft processes with the treatment of hard partons propagating hot and dense nuclear medium. The model is employed for the analysis of PbPb collisions at LHC energies. The cross-talk of elliptic v2 and triangular v3 flow in the model generates both even and odd harmonics of higher order, as well as long-range azimuthal dihadron correlations (“ridge” structure). The simple modification of the model via introducing the distribution over spatial anisotropy parameters permits HYDJET++ to reproduce both v2 and v3 flow fluctuations and related to it eccentricity fluctuations of the initial state. The data on momentum spectra and elliptic flow of charmed mesons are reproduced by HYDJET++ including thermal and non-thermal charm production mechanisms. A tolerable agreement of HYDJET++ predictions with the data for hard probes (including such observables as high-pT hadron and jet nuclear modification factors, and specific modification of longitudinal and radial jet profiles) testifies in favour of successful treatment of the hard component based on PYQUEN event generator.