Anisotropic Flow Measurements of Identiﬁed Hadrons with MPD Detector at NICA

: The primary scientiﬁc mission of the Multi-Purpose Detector (MPD) at the accelerator Nuclotron-based Ion Collider facility (NICA) (Dubna) is to investigate the properties of strongly interacting matter at high net-baryon densities. The goal of this work is to study the performance of the MPD detector for directed and elliptic ﬂow measurements of identiﬁed hadrons by using the realistic Monte Carlo simulations of heavy-ion collisions at energies √ s NN = 4.5–11 GeV.


Introduction
The Nuclotron-based Ion Collider facility (NICA) is under construction at the Joint Institute for Nuclear Research (JINR), Dubna, Russia. NICA will collide heavy-ions ( 198 Au + 198 Au, 209 Bi + 209 Bi) at energies in the range of √ s NN = 4 to 11 GeV per nucleon pair in the center-of-mass system to provide an opportunity of studying the matter in the region of high net-baryon density [1]. The Multi-Purpose Detector (MPD) experiment at NICA will measure various prominent diagnostic probes sensitive to the Equation-of-State (EOS) and transport properties of the strongly interacting matter [2,3]. Among them, the most prominent one is the azimuthal collective flow of the produced hadrons relative to the collision symmetry plane [4]. It can be quantified by Fourier coefficients v n in the expansion of the particle azimuthal distribution as: where n is the order of the harmonic, ϕ is the azimuthal angle of particle of the given type, and Ψ n is the azimuthal angle of the n th order symmetry plane. The n th order flow coefficients v n can be calculated as v n = cos[n(ϕ − Ψ n )] , where the brackets denote the average over the particles and events. Directed (v 1 ) and elliptic (v 2 ) flows are the dominant flow signals and have been studied by experiments carried out at several facilities such as the Schwer-Ionen-Synchrotron (SIS), the Alternating Gradient Synchrotron (AGS), the Super Proton Synchrotron (SPS), the Relativistic Heavy-Ion Collider (RHIC) and the Large Hadron Collider (LHC) [4,5]. The collision energy dependence of flow coefficients is of particular importance. The v 1 and v 2 results from the STAR Beam Energy Scan program at RHIC [6][7][8][9], NA61/SHINE program at SPS [10], E895 program at AGS [11,12] and HADES program at SIS [13] have indicated strong energy dependence at NICA energies [14]. In this work, we discuss the anticipated performance of the MPD detector for differential anisotropic flow measurements of identified hadrons at NICA energies.

MPD Detector System and Simulation Setup
The MPD has been designed as a 4π spectrometer for detecting charged hadrons, electrons and photons in heavy-ion collisions at high luminosity. In the first stage of operation in 2023, the MPD will consist of the Time Projection Chamber (TPC), the Time-of-Flight (TOF) detector, the electromagnetic calorimeter (ECal), the forward hadron calorimeter (FH-Cal) and the fast forward detector (FFD); see left panel of Figure 1. The TPC will provide 3D tracking of charged particles, as well as measuring the specific ionization energy loss dE/dx to identify the particles with |η| < 1.2. The TPC will be surrounded by a cylindrical barrel of the Time-of-Flight (TOF) detector. The TOF system of the MPD developed to identify the charged hadrons is based on the technology of Multigap Resistive Plate Chambers (MRPC). The detector is designed to provide both the time and coordinate measurements with the accuracy of the order of 80 ps and 0.5 cm, respectively. Two arms of hadron calorimeters (FHCal), made of 44 modules of lead-scintillator sampling calorimeters and covering the pseudorapidity range of 2.0 < |η| < 5, will measure the forward going energy distribution. The information from FHCal about the energy deposit will be used to reconstruct the event plane, see right panel of Figure 1. In this work, we used the cascade version of the Ultra-relativistic Quantum Molecular Dynamics (UrQMD) model (version 3.4) [15,16] to simulate the heavy-ion collisions at NICA energies: √ s NN = 4.5, 7.7 and 11 GeV. In total, the sample of 100 M of minimum bias Au + Au and Bi + Bi events was used to analyze the directed and elliptic flow signals of the identified hadrons. We applied the term "true" for these v 1 and v 2 results, obtained from the direct analysis of the generated events. At the next step, the same sample of UrQMD minimum bias events was made as an input for the full chain of realistic simulations of the MPD detector subsystems based on the GEANT4 [17] platform and reconstruction algorithms built in the MPDROOT. The v n results obtained from the flow analysis of these fully reconstructed events are termed as the "reco" data.

Analysis Procedure
The track reconstruction in MPD is based on the Kalman filter technique [18] and the minimal requirement of 16 TPC hits ensures a low momentum error. We introduced a 3D distance of the closest approach (DCA) between the track and the reconstructed primary vertex. The primary tracks were selected with the 2σ cut on the DCA. The relative transverse momentum resolution (∆p T /p T ) for primary tracks as a function of p T is shown in the left panel of Figure 2. ∆p T /p T is less than 3% for tracks with 0.1 < p T < 1.8 GeV/c and for the pseudorapidity range |η| < 1.5, see right panel of Figure 2. The resolution degrades rapidly above η = 1.5 due to the decrease of the number of TPC space points [18].
The analysis was performed for tracks from the kinematic regions of TPC with the higher tracking efficiency: 0.2 < p T < 2.5 GeV/c and |η| < 1.5. The centrality classes are defined based on the uncorrected charged particle multiplicity (N ch ) distribution in the TPC for pseudorapidity |η| < 0.5 and full azimuth. As an example, Figure 3 shows the track multiplicity distribution for the Au + Au (left panel) and Bi + Bi (right panel) collisions at √ s NN = 7.7 GeV compared to those from Monte Carlo (MC)-Glauber simulations (red line) [19]; 10% centrality classes are indicated with black vertical lines in the Figure 3. Identification of the charged hadrons in the MPD experiment is based on a combination of momentum information, the specific energy loss (dE/dx) in the Time-Projection Chamber (TPC) and time-of-flight measurements from the TOF detector [18]. Time-of-flight measurements from the TOF detector were used in conjunction with the measured particle momentum and flight-path length to generate a mass-squared (m 2 ) distribution, see left panel of Figure 4. The identified candidates (hadrons and light nuclei) can be selected within the predefined elliptical ranges around a nominal position in dE/dx and m 2 axes fixed for each type of particles, see right panel of Figure 4. The widths of these distributions provide an additional criterion of identification. Charged pions and kaons can be easily distinguished up to 1.5 GeV/c in transverse momentum, whereas at higher momenta, the particle species start to significantly overlap. Protons and mesons can be separated up to 2.5 GeV/c.
Short-lived weakly-decaying particles, such as K 0 s and Λ, have been reconstructed using the invariant mass technique. The combinatorial background from uncorrelated particles has been reduced by the selection criteria based on the topology of the specific decay. The topological information about the primary and secondary decay vertex positions, the distance of the closest approach (dca) of the daughter particles to the primary vertex, the dca of the mother particle to the primary vertex, and the dca between the daughter tracks have been obtained by the Kalman filtering algorithm based on the MpdParticle tool [20]. The cuts have been applied to optimize the signal of K 0 s and Λ particles. As an example, Figure 5 shows the invariant mass distributions for K 0 S (left panel) and Λ (right panel) particles with 0.2 < p T < 2.4 GeV/c for 10-40% central Au + Au collisions at √ s NN = 11 GeV. The event plane method is used to obtain the present results for the directed (v 1 ) and elliptic (v 2 ) flow coefficients of the produced particles. The event plane method correlates the azimuthal angle ϕ of each particle with the azimuthal angle Ψ n of the event plane determined from the anisotropic flow itself [4,21]. Directed flow v 1 is large at NICA energies compared to other flow harmonics [14]. It is the strongest in the forward rapidity region (i.e., in FHCal acceptance area: 2< |η| < 5). For these reasons the first harmonic event plane Ψ 1,FHCal is used to study the present flow performance. The event plane angle Ψ 1,FHCal has been calculated from the energy deposition in a given module of the FHCal by constructing the so-called flow Q-vector Q 1,FHCal (two-dimensional vector in the transverse to the beam plane): where ϕ i is the azimuthal angle of the center of the i th FHCal module in the transverse plane, and E i is the energy deposition in the i th module of FHCal (weight to improve the event plane resolution). The weights E i have opposite signs for backward and forward rapidities due to the anti-symmetry of the v 1 as a function of rapidity y. The reconstructed Ψ 1,FHCal can be used to measure the differential directed flow v .
where R n (Ψ 1,FHCal ) represents the event plane resolution factor. The two-subevent method with extrapolation algorithm was applied to estimate the event plane resolution factors [21]. Figure 6 shows the centrality dependence of the event plane resolution factors R n (Ψ 1,FHCal ) for directed v 1 (left panel) and elliptic v 2 (right panel) flow measurements for Au + Au collisions at √ s NN = 4.5, 7.7 and 11 GeV. For the mid-central Au + Au events at √ s NN =11 GeV the resolution factor is as high as 0.9-for v 1 and 0.65-for v 2 measurements. The event plane resolution degrades slowly while decreasing the collision energy.    interval from fits to the K 0 s and Λ invariant mass distributions, see Figure 5. Values for v S n signal for K 0 S and Λ particles were extracted via a direct fit to the v SB n (M inv ) for each p pair T selection by Equation (4). That is, the background v B n (M inv ) was parametrized as a linear or quadratic function of M inv (depending on the p pair T selection) and v S n is taken as a fit parameter. The accuracy of the extraction procedure was verified by checking that the invariant mass (M inv ) dependencies of the sine coefficients v SB sin,n = sin[n(ϕ pair − Ψ 1,FHCal )] , were all equal to zero within statistical errors. As an example, Figure 8 shows the demonstration of the invariant mass fit method to extract the v S 1 and v S 2 signals for K 0 S (upper part) and Λ (lower part) particles for 10-40% central Au + Au collisions at √ s NN = 11 GeV.
The v n results obtained by the event plane analysis can be affected by non-flow and flow fluctuations [4]. The non-flow effects are mainly due to few particle correlations not associated with the reaction plane: Bose-Einstein correlations, resonance decays, momentum conservation and di-jets. In the present study, the large rapidity gap ∆η > 0.  The open symbols correspond to v n results from the analysis of the fully reconstructed events "reco" and closed symbols and to the results from the generated "true" UrQMD events. The available statistics of 20 M minimum-bias events have allowed us to perform the detailed p T -differential measurements of charged pions, kaons and protons up to 1.5 GeV/c. The more detailed p T -differential studies as a function of centrality and rapidity will require a larger data sample of up to 300 M of minimum-bias events. The open symbols correspond to v n results from the analysis of the fully reconstructed events ("reco") and closed symbols to the results from generated ("true") UrQMD events. Figures 11 and 12 illustrate the MPD detector system performance for the differential directed v 1 (y) and elliptic flow measurements v 2 (p T ) of K 0 S and Λ particles for 10-40% central Au + Au collisions at √ s NN = 11 GeV. The results were obtained from the analysis of 25 M minimum-bias fully reconstructed UrQMD events. A good agreement is observed between the v n results of K 0 S and Λ particles extracted by invariant-mass fit method illustrated in the Figure 8 from the fully reconstructed data "reco" and generated "true" UrQMD events. The final part of this study is related to the comparison of the differential v 1 and v 2 results of identified hadrons from two colliding systems Au + Au and Bi + Bi at √ s NN = 7.7 GeV. The NICA collider is planned to start with first beams of 209 Bi ions in 2023. The delivery of 198 Au ions will be accomplished after this phase of NICA operation. Figure 13 presents the detailed comparison of the differential v 1 (y) results for charged pions and protons for 0-10% central (left panels) and 10-40% central (right panels) Au + Au and Bi + Bi collisions at √ s NN = 7.7 GeV. The v 1 (y) results have been obtained from the event plane analysis of 30 M fully reconstructed minimum-bias UrQMD events. The results received for Au + Au (Bi + Bi) collisions are marked as open (filled) symbols. Figure 14 shows the performance of measuring the p T dependence of directed v 1 (left) and

Results
elliptic v 2 (right) flow coefficients of charged pions and protons for 10-40% central Au + Au (open symbols) and Bi + Bi (filled symbols) collisions. The expected small difference is observed in the v n results for two colliding systems.

Conclusions
In this work, the performance of the MPD experiment for directed (v 1 ) and elliptic (v 2 ) flow measurements was studied with Monte Carlo simulations using 198 Au + 198 Au and 209 Bi + 209 Bi collisions at NICA energies. A large sample of generated UrQMD minimum bias events has been used as an input for the full chain of realistic simulations of the MPD detector subsystems based on the GEANT4 platform and reconstruction algorithms built in the MPDROOT. Realistic procedures for centrality determination, particle identification and event plane reconstruction have been used in the analysis. The resulting performance of the MPD has been verified for v 1 and v 2 measurements of identified charged pions, kaons, protons, K 0 s and Λ particles as a function of rapidity and transverse momentum in different centrality classes. The detailed comparison of the results obtained from the analysis of the fully reconstructed data and generator-level data has allowed us to conclude that MPD system will provide the detailed differential measurements of directed and elliptic flows with high efficiency. In future we plan to include the data from other transport models [23] and extend the study to other colliding systems.