Alexander Snigirev

Heavy ion event generator HYDJET++ (HYDrodynamics plus JETs)

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 and Snigirev, EPJC 45 (2006) 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., PRC 74 (2006) 064901; PRC 77 (2008) 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. Program summary Program title: HYDJET++, version 2 Catalogue identifier: AECR_v1_0 Program summary URL:http://cpc.cs.qub.ac.uk/summaries/AECR_v1_0.html Program obtainable from: CPC Program Library, Queens University, Belfast, N. Ireland Licensing provisions: Standard CPC licence, http://cpc.cs.qub.ac.uk/licence/licence.html No. of lines in distributed program, including test data, etc.: 100 387 No. of bytes in distributed program, including test data, etc.: 797 019 Distribution format: tar.gz Programming language: C++ (however there is a Fortran-written part which is included in the generator structure as a separate directory) Computer: Hardware independent (both C++ and Fortran compilers and ROOT environment [1] (http://root.cern.ch/) should be installed) Operating system: Linux (Scientific Linux, Red Hat Enterprise, FEDORA, etc.) RAM: 50 MBytes (determined by ROOT requirements) Classification: 11.2 External routines: ROOT [1] (http://root.cern.ch/) Nature of problem: The experimental and phenomenological study of multi-particle production in relativistic heavy ion collisions is expected to provide valuable information on the dynamical behavior of strongly-interacting matter in the form of quark–gluon plasma (QGP) [2–4], as predicted by lattice Quantum Chromodynamics (QCD) calculations. Ongoing and future experimental studies in a wide range of heavy ion beam energies require the development of new Monte Carlo (MC) event generators and improvement of existing ones. Especially for experiments at the CERN Large Hadron Collider (LHC), implying very high parton and hadron multiplicities, one needs fast (but realistic) MC tools for heavy ion event simulations [5–7]. The main advantage of MC technique for the simulation of high-multiplicity hadroproduction is that it allows a visual comparison of theory and data, including if necessary the detailed detector acceptances, responses and resolutions. The realistic MC event generator has to include maximum possible number of observable physical effects, which are important to determine the event topology: from the bulk properties of soft hadroproduction (domain of low transverse momenta pT≲1 GeV/cpT≲1 GeV/c) such as collective flows, to hard multi-parton production in hot and dense QCD-matter, which reveals itself in the spectra of high-pTpT particles and hadronic jets. Moreover, the role of hard and semi-hard particle production at LHC can be significant even for the bulk properties of created matter, and hard probes of QGP became clearly observable in various new channels [8–11]. In the majority of the available MC heavy ion event generators, the simultaneous treatment of collective flow effects for soft hadroproduction and hard multi-parton in-medium production (medium-induced partonic rescattering and energy loss, so-called “jet quenching”) is lacking. Thus, in order to analyze existing data on low and high-pTpT hadron production, test the sensitivity of physical observables at the upcoming LHC experiments (and other future heavy ion facilities) to the QGP formation, and study the experimental capabilities of constructed detectors, the development of adequate and fast MC models for simultaneous collective flow and jet quenching simulations is necessary. HYDJET++ event generator includes detailed treatment of soft hadroproduction as well as hard multi-parton production, and takes into account known medium effects. Solution method: A heavy ion event in HYDJET++ is a superposition of the soft, hydro-type state and the hard state resulting from multi-parton fragmentation. Both states are treated independently. HYDJET++ is the development and continuation of HYDJET MC model [12]. The main program is written in the object-oriented C++ language under the ROOT environment [1]. The hard part of HYDJET++ is identical to the hard part of Fortran-written HYDJET [13] (version 1.5) and is included in the generator structure as a separate directory. The routine for generation of single hard NN collision, generator PYQUEN [12,14], modifies the “standard” jet event obtained with the generator PYTHIA 6.4 [15]. The event-by-event simulation procedure in PYQUEN includes 1. generation of initial parton spectra with PYTHIA and production vertexes at given impact parameter; 2. rescattering-by-rescattering simulation of the parton path in a dense zone and its radiative and collisional energy loss; 3. final hadronization according to the Lund string model for hard partons and in-medium emitted gluons. Then the PYQUEN multi-jets generated according to the binomial distribution are included in the hard part of the event. The mean number of jets produced in an AA event is the product of the number of binary NN subcollisions at a given impact parameter and the integral cross section of the hard process in NN collisions with the minimum transverse momentum transfer pTmin. In order to take into account the effect of nuclear shadowing on parton distribution functions, the impact parameter dependent parameterization obtained in the framework of Glauber–Gribov theory [16] is used. 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 (the adapted C++ code FAST MC [17,18]). Hadron multiplicities are calculated using the effective thermal volume approximation and Poisson multiplicity distribution around its mean value, which is supposed to be proportional to the number of participating nucleons at a given impact parameter of AA collision. The fast soft hadron simulation procedure includes 1. generation of the 4-momentum of a hadron in the rest frame of a liquid element in accordance with the equilibrium distribution function; 2. generation of the spatial position of a liquid element and its local 4-velocity in accordance with phase space and the character of motion of the fluid; 3. the standard von Neumann rejection/acceptance procedure to account for the difference between the true and generated probabilities; 4. boost of the hadron 4-momentum in the center mass frame of the event; 5. the two- and three-body decays of resonances with branching ratios taken from the SHARE particle decay table [19]. The high generation speed in HYDJET++ is achieved due to almost 100% generation efficiency of the “soft” part because of the nearly uniform residual invariant weights which appear in the freeze-out momentum and coordinate simulation. Although HYDJET++ is optimized for very high energies of RHIC and LHC colliders (c.m.s. energies of heavy ion beams s=200 and 5500 GeV per nucleon pair, respectively), in practice it can also be used for studying the particle production in a wider energy range down to s∼10 GeV per nucleon pair at other heavy ion experimental facilities. As one moves from very high to moderately high energies, the contribution of the hard part of the event becomes smaller, while the soft part turns into just a multi-parameter fit to the data. Restrictions: HYDJET++ is only applicable for symmetric AA collisions of heavy (A≳40A≳40) ions at high energies (c.m.s. energy s≳10 GeV per nucleon pair). The results obtained for very peripheral collisions (with the impact parameter of the order of two nucleus radii, b∼2RAb∼2RA) and very forward rapidities may be not adequate. Additional comments: Accessibility http://cern.ch/lokhtin/hydjet++ Running time: The generation of 100 central (0–5%) Au+Au events at s=200A GeV (Pb+Pb events at s=5500A GeV) with default input parameters takes about 7 (85) minutes on a PC 64 bit Intel Core Duo CPU @ 3 GHz with 8 GB of RAM memory under Red Hat Enterprise. References: [1] I.P. Lokhtin, A.M. Snigirev, Eur. Phys. J. C 46 (2006) 211. [2] N.S. Amelin, R. Lednicky, T.A. Pocheptsov, I.P. Lokhtin, L.V. Malinina, A.M. Snigirev, Iu.A. Karpenko, Yu.M. Sinyukov, Phys. Rev. C 74 (2006) 064901. [3] N.S. Amelin, I. Arsene, L. Bravina, Iu.A. Karpenko, R. Lednicky, I.P. Lokhtin, L.V. Malinina, A.M. Snigirev, Yu.M. Sinyukov, Phys. Rev. C 77 (2008) 014903.

QCD status of factorization ansatz for double parton distributions

Recent CDF measurements of the inclusive cross section for a double parton scattering attach a great importance to any theoretical calculations of two-particle distribution functions. Using a parton interpretation of the leading logarithm diagrams of perturbative QCD theory, generalized Lipatov-Altarelli-Parisi-Dokshitzer equations for the two-parton distributions are re-obtained. The solutions of these equations are not at all the product of two single-parton distributions what is usually applied to the current analysis as ansatz.

Double heavy meson production through double parton scattering in hadronic collisions

Abstract It is shown that the contribution from double parton scattering to the inclusive double heavy meson yield is quite comparable with the usually considered mechanism of their production at the LHC energy. For some pairs of heavy flavored quarks in the final state the double parton scattering will be a dominant mode of their production.

Nuclear geometry of jet quenching

The most suitable way to study the jet quenching as a function of distance traversed is varying the impact parameter b of ultrarelativistic nucleus-nucleus collision (initial energy density in nuclear overlapping zone is almost independent of b up to b R_A). It is shown that b-dependences of medium-induced radiative and collisional energy losses of a hard parton jet propagating through dense QCD-matter are very different. The experimental verification of this phenomenon could be performed for a jet with non-zero cone size basing on essential difference between angular distributions of collisional and radiative energy losses.

Angular structure of energy losses of hard jet in dense QCD-matter

Abstract Angular structure of radiative and collisional energy losses of a hard parton jet propagating through dense QCD-matter is investigated. For small angular jet cone sizes, θ 0 ≲5°, the radiative energy loss is shown to dominate over the collisional energy loss due to final state elastic rescattering of the hard projectile on thermal particles in the medium. Due to coherent effects, the radiative energy loss decreases with increasing the angular size the jet. It becomes comparable with the collisional energy loss for θ 0 ≳5°−10°.Angular structure of radiative and collisional energy losses of a hard parton jet propagating through dense QCD-matter is investigated. For small angular jet cone sizes,

Fast hadron freeze-out generator. II. Noncentral collisions

\theta_0<5^0

Enhanced J/ψ-pair production from double-parton scatterings in nucleus–nucleus collisions at the Large Hadron Collider

, the radiative energy loss is shown to dominate over the collisional energy loss due to final state elastic rescattering of the hard projectile on thermal particles in the medium. Due to coherent effects, the radiative energy loss decreases with increasing the angular size the jet. It becomes comparable with the collisional energy loss for

Probing the medium-induced energy loss of bottom quarks by dimuon production in heavy ion collisions at LHC

\theta_0>5^0-10^0

Same-sign WW production in proton–nucleus collisions at the LHC as a signal for double parton scattering

.

Towards Azimuthal Anisotropy of Direct Photons

The fast Monte Carlo procedure of hadron generation developed in our previous work is extended to describe noncentral collisions of nuclei. We consider different possibilities to introduce appropriate asymmetry of the freeze-out hypersurface and flow velocity profile. For comparison with other models and experimental data, we demonstrate the results based on the standard parametrizations of the hadron freeze-out hypersurface and flow velocity profile assuming either a common chemical and thermal freeze-out or the chemically frozen evolution from chemical to thermal freeze-out. The C++ generator code is written under the ROOT framework and is available for public use at http://uhkm.jinr.ru/.