Reconstruction of Photon Conversions in the MPD Experiment

: Spectra of thermal photons carry important information on the temperature of the hot and dense medium produced in heavy ion collisions. Photons can be measured via their conversion into electron-positron pairs in the detector material. In this contribution, challenges in the photon reconstruction are discussed and feasibility studies on photon conversion measurements in the future multipurpose detector (MPD) experiment at NICA are presented. The obtained results indicate good prospects for thermal photon measurements.


Introduction
The future multipurpose detector (MPD) at the Nuclotron-based Ion Collider fAcility (NICA) is aimed towards the studies of hot and dense medium produced in heavy ion collisions in the range of center-of-mass energies √ s NN = 4 to 11 GeV [1]. In contrast to heavy ion collisions at the Relativistic Heavy Ion Collider (RHIC) and the Large Hadron Collider (LHC) corresponding to vanishing net baryonic densities, the matter produced in the NICA energy range is expected to reach large baryonic densities far exceeding the normal nuclear density ρ 0 . According to several model calculations described in [2], the nuclear fireball created in central heavy ion collisions will be compressed to a density above 8ρ 0 already at the lowest NICA energy of √ s NN = 4 GeV. At higher NICA energies, the maximum net baryon density is expected to increase further. Thus, heavy ion collisions at NICA complement experimental efforts at RHIC and LHC colliders in the studies of the phase diagram of nuclear matter and searches for the deconfinement phase transition. This energy range will be also studied by the STAR experiment in the RHIC beam energy scan II program (down to √ s NN = 7.7 GeV in the collider mode and down to √ s NN = 3 GeV in the fixed target mode) [3]. Even lower center-of-mass energies ( √ s NN from 2 to 5 GeV) will be accessed by the High Acceptance DiElectron Spectrometer (HADES) and the Compressed Baryonic Matter (CBM) experiments at the Facility for Antiproton and Ion Research (FAIR) [4,5]. This broad and complementary research program provides an opportunity to prob the nuclear equation-of-state at high baryonic densities which is necessary for an understanding of neutron star interiors [6]. The construction of the MPD experiment will be carried out in two stages. In the first stage, the time projection chamber (TPC), the time-of-flight (TOF) detector, the electromagnetic calorimeter (ECal), the forward hadron calorimeter (FHCal), and the fast forward detector (FFD) will be installed. The TPC will be placed in a superconducting solenoid providing efficient tracking, particle identification (PID) and momentum measurements in a wide pseudorapidity range |η| < 1.2. The TOF detector will be based on multigap resistive plate counters and will be used to perform particle identification up to a momenta of about 2 GeV/c. The ECal will be used to identify electrons, photons, and neutral hadrons in the central rapidity region while FHCal is aimed for centrality determination purposes. The FFD will be used for trigger purposes. In the second stage, the inner tracking system (ITS) located between the beam pipe and the TPC will be installed for a precise reconstruction of primary and secondary vertices.
Among various observables studied in heavy ion collisions, direct photons play a special role. They are emitted at all stages of the evolution of the produced hadronic matter and carry information on the temperature of the medium at a time of their emission [7]. One usually distinguishes high-energy prompt photons produced in hard scatterings at early stages and low-energy thermal photons radiated by the expanding fireball. The inverse slope of the thermal photon momentum distribution reflects the effective temperature averaged over space-time evolution of the medium. The effective photon temperatures measured at RHIC [8] and LHC energies [9] appear to be much larger than the temperature of about 170 MeV predicted for the transition from the hadron gas to the deconfined quarkgluon plasma phase [10]. According to the RHIC beam energy scan results [11], the effective temperature tends to decrease towards lower collision energies. However, according to theoretical predictions [12,13], the initial temperatures in central Au-Au collisions at NICA energies still raise above 200 MeV, while time-averaged effective temperatures are expected to be in the range from 150 to 170 MeV. Thus the MPD experiment has a unique opportunity to perform precision direct photon measurements and trace the transition from hadronic to partonic degrees of freedom in the NICA energy range.
One of the challenges in the direct photon measurements is related to the fact that the measured photon spectra are dominated by decay photons originating mainly from neutral π and η meson decays. In order to get the direct photon spectrum, the decay contribution has to be subtracted from the inclusive photon spectrum. Usually one constructs a ratio R γ of inclusive and decay photon yields where decay photon yields are determined from the measured spectra of π 0 and η mesons [7].
The MPD experiment will explore two methods to measure photons. The first method will utilize the ECal capabilities. The second method will be based on the reconstruction of photons converted in the detector material into electron-positron pairs measured in the TPC. The latter method has an advantage of much better energy resolution at low energies and smaller combinatorial backgrounds in neutral meson measurements. However the conversion method has a disadvantage of a much smaller efficiency that is proportional to the effective radiation length of the detector material in front of the tracking volume. In the first stage of the MPD experiment, photons will mainly be converted in the beam pipe (0.3% X 0 ) and in the inner TPC vessel structures (2.4% X 0 ) therefore the conversion reconstruction efficiency will hardly exceed a few percent.
In this contribution, feasibility studies on photon conversion measurements in the future MPD experiment at NICA are presented.

Photon Reconstruction
The studies of the photon conversion probability were performed with the GEANT4 package coupled to the official MpdRoot framework used for Monte Carlo simulations of the MPD detector response, track reconstruction, and pattern recognition. The photon conversion probability is shown in Figure 1 as a function of the photon transverse momentum p T for photons in the central pseudorapidity range |η| < 1. The conversion probability rapidly grows towards a higher momenta and reaches a plateau of about 3.5% at p T > 0.5 GeV/c. As shown in Figure 1, an application of fiducial selection criteria on the transverse momentum (p e T > 50 MeV/c) and pseudorapidity (|η e | < 1) of conversion electrons and positrons results in a reduction of the fraction of accepted photons, mainly at a low transverse momenta. Reconstruction of converted photons requires the efficient identification of electrons. The MPD experiment provides efficient electron separation using energy loss measurements in the TPC in the momentum range 0.2 < p e < 0.5 GeV/c while TOF measurements allow one to access the momentum region p e < 0.3 GeV/c. A dedicated procedure to select electron candidates has been developed. Tracks with energy loss and effective squared mass measurements are selected in ±4 standard deviations around their mean values at a given momentum. Distributions of energy losses in the TPC and reconstructed squared mass in the TOF for electrons are shown in Figure 2 together with corresponding selection limits. Secondary vertices of electron-positron pairs from conversions were reconstructed using a dedicated MpdParticle tool that allows one to build a vertex and determine its parameters for any combination of positive and negative tracks. The tool is based on the Kalman filter algorithm [14]. Figure 3 shows typical distributions of the χ 2 parameter representing the quality of secondary vertices for signal electron-positron pairs and random background pairs as a function of pair p T . A requirement χ 2 < 10, used to reject random pairs of tracks, is also shown in the figure. The χ 2 selection does not remove background pairs of tracks coming from the primary vertex, such as electron-positron pairs from Dalitz decays, since they are characterized by small distances of the closest approach. In order to remove these pairs from the sample of conversion candidates, one can reject secondary vertices at small distances from the primary vertex. Distributions of signal and background pairs are shown in Figure 4 as a function of the radial position of the vertex (R vtx ) and transverse momentum of the reconstructed pair. Signal conversion pairs originate from radii close to the beam pipe and inner TPC structures while random background pairs mainly come from radii below 20 cm. The radial resolution of secondary vertices is rather poor therefore a tight requirement on R vtx > 10 cm is used to suppress the background from random pairs of primary tracks. Background pairs of tracks can be further suppressed using the selection on the invariant mass of the pair assuming electron mass hypothesis for daughter tracks. Distributions of signal and background pairs as a function of invariant mass and transverse momentum of the reconstructed pair are shown in Figure 5. The invariant masses for conversion pairs are close to zero. As for the background pairs, the distribution is much broader though the ridge at small masses is also visible. The ridge corresponds to electron-positron pairs coming from a primary vertex (e.g., from Dalitz decays) and can be suppressed with R vtx selection. Since the invariant mass resolution depends on daughter track momenta, we use a p T -dependent selection: m pair < 0.0226 + 0.0174 p pair T where m pair and p pair T are measured in GeV units.  Finally the important difference between track pairs from photon conversions and weak decays is that the opening angle between electron and positron tracks from conversions is small resulting in a preferred direction of the dielectron plane with respect to the magnetic field vector. Distributions of the angle Ψ pair between the plane perpendicular to the magnetic field and the plane which is spanned by the e + and e − momentum vectors is shown in Figure 7 for signal and background pairs. A selection |Ψ pair | < 0.1 helps to improve the purity of the converted photon sample. All the selection criteria were carefully optimized using a multivariate technique with the aim to maximize significance of the signal. The reconstruction efficiency after optimized single electron track cuts and after pair cuts is shown in Figure 8. The conversion reconstruction efficiency reaches a maximum of about 1.5%. Contribution of the nonphoton background does not exceed 10-20% and can be further improved with tighter selection criteria.

Direct and Decay Photon Production Estimates
Widely used general purpose event generators do not simulate the production of direct photons and hence the direct photon feasibility studies require external inputs. A data-driven method was developed to estimate the spectra of direct photons in Au-Au collisions at NICA energies. The method relies on the universal scaling of direct photon yields with the (pseudo)rapidity density of final state charged particles dN ch /dη observed at RHIC and LHC energies [11]. The PHENIX experiment found that the integrated yields of direct photons above p T > 1 GeV/c scale as a power law (dN ch /dη) 1.25 in a wide range of centralities and collision energies. This scaling works well for collisions with dN ch /dη > 20 and is broken only at lower multiplicities. Assuming the scaling still holds at NICA energies for relatively high multiplicities, one can expect its validity for Au-Au collisions in the range of (0-60)% centrality at √ s NN = 11 GeV and (0-40)% centrality at √ s NN = 4 GeV.
Transverse momentum spectra of direct photons can be also predicted assuming a universal p T scaling observed by PHENIX in the region p T > 0.4 GeV/c [11,15]. Differential direct photon yields in Au-Au collisions at √ s NN = 4 and 11 GeV were calculated in the region of p T > 0.6 GeV/c assuming a d 3 N d 2 p T dy = 2.755 × 10 −4 · p −4.5 T × (dN ch /dη) 1.25 dependence as shown in Figure 9, left. In the region of low p T < 0.6 GeV/c, direct photon spectra were assumed to follow the thermal spectrum dN/dp T ∼ p T exp(−p T /T eff ) with a conservative effective temperature T eff = 150 MeV [12,13].
The decay photon spectrum coming from π 0 → γγ and η → γγ decays was estimated in accordance with projections on p T -differential π 0 and η meson measurements described in [16][17][18]. The inclusive photon yields were obtained as a sum of direct and decay photon yields. The excess photon ratio R γ was calculated as a ratio of an inclusive photon spectrum to decay photon spectrum and is shown in Figure 9, right, in three centrality classes. The performed study shows that R γ significantly exceeds unity starting from 1.03 for semi-central collisions and increases with p T exceeding 1.2 at 1.4 GeV similar to previous measurements at higher energies. Assuming a sample of about 10 9 minimum bias events, expected statistical uncertainties on the inclusive photon spectrum will not exceed 1% in the considered p T range. The precision of the R γ ratio and the direct photon spectrum measurement will be limited mainly by the systematic uncertainty related to the knowledge of the decay photon contribution in the inclusive photon yield. According to the experience of other experiments, typical systematic uncertainties on R γ can be reduced to about 2% [11]. Thus the obtained R γ estimates justify the possibility to measure direct photon production with the MPD detector down to low transverse momenta.  Figure 9. The expected direct photon spectra (left) and excess photon ratio R γ (right) as a function of the photon transverse momentrum in Au-Au collisions at NICA energies.

Conclusions
In this contribution, we summarized that photons serve as valuable probes of the hot and dense medium produced in heavy ion collisions. The procedure of photon reconstruction with the photon conversion method was described. The feasibility of measuring thermal photon production spectra in Au-Au collisions at √ s NN = 11 GeV with the MPD detector at NICA was shown.