Next Article in Journal
Transport Properties in Magnetized Compact Stars
Next Article in Special Issue
Anisotropic Flow Measurements of Identified Hadrons with MPD Detector at NICA
Previous Article in Journal
Particles —First 100 Papers Milestone Achieved
Previous Article in Special Issue
Deep Learning Based Impact Parameter Determination for the CBM Experiment
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Conference Report

Reconstruction of Photon Conversions in the MPD Experiment

1
Petersburg Nuclear Physics Institute Named by B. P. Konstantinov of NRC “Kurchatov Institute”, 1 mkr. Orlova roshcha, 188300 Gatchina, Russia
2
Peter the Great St. Petersburg Polytechnic University, Polytechnicheskaya 29, 195251 St. Petersburg, Russia
3
MEPhI (Moscow Engineering Physics Institute), National Research Nuclear University, 31 Kashirskoe Shosse, 115409 Moscow, Russia
*
Author to whom correspondence should be addressed.
For the MPD Collaboration.
Submission received: 13 January 2021 / Revised: 5 February 2021 / Accepted: 7 February 2021 / Published: 17 February 2021
(This article belongs to the Special Issue Analysis Techniques and Physics Performance Studies for FAIR and NICA)

Abstract

:
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.

1. 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 quark-gluon 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.

2. 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 T e > 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 T pair where m pair and p T pair are measured in GeV units.
The reconstructed total momentum of electron-positron pairs from converted direct photons should point to the primary vertex. The selection on the pointing angle, i.e., the angle between pair momentum vector and direction to the primary vertex, can be used to suppress the background further. Distributions of signal and background pairs as a function of pointing angle and transverse momentum of the reconstructed pair are shown in Figure 6. The pointing angle resolution improves towards higher transverse momentum of the pair. Therefore the following p T pair -dependent selection on the pointing angle was found to be optimal: θ < exp ( 2.777 2.798 p T pair ) + 0.0175 where p T pair is 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 non-photon background does not exceed 10–20% and can be further improved with tighter selection criteria.

3. 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 d N 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 ( d N ch / d η ) 1.25 in a wide range of centralities and collision energies. This scaling works well for collisions with d N 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 d y = 2.755 × 10 4 · p T 4.5 × ( d N 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 d N / d p 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.

4. 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.

Author Contributions

Conceptualization, E.K., V.R., Y.R.; methodology, E.K. and V.R.; software, M.M., Y.R. and D.I.; validation, D.K. and D.I.; formal analysis, E.K.; investigation, E.K. and V.R.; writing—original draft preparation, E.K.; writing—review and editing, V.R.; visualization, E.K.; supervision, D.I. and E.K.; project administration, D.I.; funding acquisition, E.K. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by RFBR according to the research project # 18-02-40045 and partially supported by the National Research Nuclear University MEPhI in the framework of the Russian Academic Excellence Project (contract No. 02.a03.21.0005, 27.08.2013).

Conflicts of Interest

The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript, or in the decision to publish the results.

References

  1. Golovatyuk, V.; Kekelidze, V.; Kolesnikov, V.; Rogachevsky, O.; Sorin, A. The Multi-Purpose Detector (MPD) of the collider experiment. Eur. Phys. J. A 2016, 52, 212. [Google Scholar] [CrossRef]
  2. Arsene, I.C.; Bravina, L.V.; Cassing, W.; Ivanov, Y.B.; Larionov, A.; Randrup, J.; Russkikh, V.N.; Toneev, V.D.; Zeeb, G.; Zschiesche, D. Dynamical phase trajectories for relativistic nuclear collisions. Phys. Rev. C 2007, 75, 034902. [Google Scholar] [CrossRef] [Green Version]
  3. Odyniec, G. Beam Energy Scan Program at RHIC (BES I and BES II)–Probing QCD Phase Diagram with Heavy-Ion Collisions. PoS 2019, CORFU2018, 151. [Google Scholar] [CrossRef] [Green Version]
  4. Senger, P. Probing dense QCD matter in the laboratory—The CBM experiment at FAIR. Phys. Scr. 2020, 95, 074003. [Google Scholar] [CrossRef]
  5. Salabura, P. Studying Baryonic Matter with HADES at GSI/FAIR. Acta Phys. Polon. B 2019, 50, 1205–1216. [Google Scholar] [CrossRef]
  6. Klähn, T.; Blaschke, D.; Typel, S.; Van Dalen, E.N.E.; Faessler, A.; Fuchs, C.; Wolter, H.H. Constraints on the high-density nuclear equation of state from the phenomenology of compact stars and heavy-ion collisions. Phys. Rev. C 2006, 74, 035802. [Google Scholar] [CrossRef] [Green Version]
  7. David, G. Direct real photons in relativistic heavy ion collisions. Rept. Prog. Phys. 2020, 83, 046301. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  8. Adare, A.; Afanasiev, S.; Aidala, C.; Ajitanand, N.N.; Akiba, Y.; Al-Bataineh, H.; Alexander, J.; Al-Jamel, A.; Aoki, K.; Aphecetche, L.; et al. Enhanced production of direct photons in Au+Au collisions at s NN = 200 GeV and implications for the initial temperature. Phys. Rev. Lett. 2010, 104, 132301. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  9. Adam, J.; Adamova, D.; Aggarwal, M.M.; Rinella, G.A.; Agnello, M.; Agrawal, N.; Ahammed, Z.; Ahn, S.U.; Aiola, S.; Akindinov, A.; et al. Direct photon production in Pb-Pb collisions at s NN = 2.76 TeV. Phys. Lett. B 2016, 754, 235–248. [Google Scholar] [CrossRef]
  10. Laermann, E.; Philipsen, O. The Status of lattice QCD at finite temperature. Ann. Rev. Nucl. Part. Sci. 2003, 53, 163–198. [Google Scholar] [CrossRef] [Green Version]
  11. Adare, A.; Afanasiev, S.; Aidala, C.; Ajitan, N.N.; Akiba, Y.; Akimoto, R.; Fleuret, F. Beam Energy and Centrality Dependence of Direct-Photon Emission from Ultrarelativistic Heavy-Ion Collisions. Phys. Rev. Lett. 2019, 123, 022301. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  12. Endres, S.; van Hees, H.; Bleicher, M. Photon and dilepton production at the Facility for Antiproton and Ion Research and beam-energy scan at the Relativistic Heavy-Ion Collider using coarse-grained microscopic transport simulations. Phys. Rev. C 2016, 93, 054901. [Google Scholar] [CrossRef] [Green Version]
  13. Rapp, R.; van Hees, H. Thermal Dileptons as Fireball Thermometer and Chronometer. Phys. Lett. B 2016, 753, 586–590. [Google Scholar] [CrossRef] [Green Version]
  14. Luchsinger, R.; Grab, C. Vertex reconstruction by means of the method of Kalman filter. Comput. Phys. Commun. 1993, 76, 263–280. [Google Scholar] [CrossRef]
  15. Khachatryan, V.; Praszałowicz, M. Scaling properties of direct photon yields in heavy ion collisions. Eur. Phys. J. C 2020, 80, 670. [Google Scholar] [CrossRef]
  16. Kryshen, E.; Ivanishchev, D.; Kotov, D.; Malaev, M.; Riabov, V.; Ryabov, Y. Perspectives of thermal photon measurements in heavy ion collisions at NICA. EPJ Web Conf. 2019, 222, 02006. [Google Scholar] [CrossRef]
  17. Kryshen, E.; Ivanishchev, D.; Kotov, D.; Malaev, M.; Riabov, V.; Ryabov, Y. Study of neutral meson production with photon conversions in the MPD experiment at NICA. J. Phys. Conf. Ser. 2019, 1400, 055055. [Google Scholar] [CrossRef]
  18. Ivanishchev, D.; Kryshen, E.; Kotov, D.; Malaev, M.; Riabov, V.; Riabov, Y. Feasibility of thermal photon measurements in heavy ion collisions at NICA energies. J. Phys. Conf. Ser. 2020, 1690, 012109. [Google Scholar] [CrossRef]
Figure 1. Photon conversion probability (black) and acceptance (red) as a function of the photon transverse momentum.
Figure 1. Photon conversion probability (black) and acceptance (red) as a function of the photon transverse momentum.
Particles 04 00008 g001
Figure 2. (Left) Distribution of energy losses in the time projection chamber (TPC) for electrons as a function of momentum. (Right) Distribution of reconstructed squared mass in the time-of-flight (TOF) for electrons as a function of momentum. Selection limits are shown with magenta lines.
Figure 2. (Left) Distribution of energy losses in the time projection chamber (TPC) for electrons as a function of momentum. (Right) Distribution of reconstructed squared mass in the time-of-flight (TOF) for electrons as a function of momentum. Selection limits are shown with magenta lines.
Particles 04 00008 g002
Figure 3. Distributions of signal (left) and background (right) pairs as a function of the vertex χ 2 and transverse momentum of the reconstructed pair. Selection limits are shown with magenta lines.
Figure 3. Distributions of signal (left) and background (right) pairs as a function of the vertex χ 2 and transverse momentum of the reconstructed pair. Selection limits are shown with magenta lines.
Particles 04 00008 g003
Figure 4. Distributions of signal (left) and background (right) pairs as a function of radial position of the vertex and transverse momentum of the reconstucted pair. Selection limits are shown with magenta lines.
Figure 4. Distributions of signal (left) and background (right) pairs as a function of radial position of the vertex and transverse momentum of the reconstucted pair. Selection limits are shown with magenta lines.
Particles 04 00008 g004
Figure 5. Distributions of signal (left) and background (right) pairs as a function of invariant mass and transverse momentum of the reconstructed pair. Selection limits are shown with magenta lines.
Figure 5. Distributions of signal (left) and background (right) pairs as a function of invariant mass and transverse momentum of the reconstructed pair. Selection limits are shown with magenta lines.
Particles 04 00008 g005
Figure 6. Distributions of signal (left) and background (right) events as a function of the pointing angle and transverse momentum of the photon pair. Selection limits are shown with magenta lines.
Figure 6. Distributions of signal (left) and background (right) events as a function of the pointing angle and transverse momentum of the photon pair. Selection limits are shown with magenta lines.
Particles 04 00008 g006
Figure 7. Distributions of signal (left) and background (right) pairs as a function of the Ψ pair angle and transverse momentum of the reconstructed pair. Selection limits are shown with magenta lines.
Figure 7. Distributions of signal (left) and background (right) pairs as a function of the Ψ pair angle and transverse momentum of the reconstructed pair. Selection limits are shown with magenta lines.
Particles 04 00008 g007
Figure 8. Photon reconstruction efficiency after single electron (magenta) and electron pair (blue) selection criteria.
Figure 8. Photon reconstruction efficiency after single electron (magenta) and electron pair (blue) selection criteria.
Particles 04 00008 g008
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.
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.
Particles 04 00008 g009
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Kryshen, E.; Ivanishchev, D.; Kotov, D.; Malaev, M.; Riabov, V.; Ryabov, Y. Reconstruction of Photon Conversions in the MPD Experiment. Particles 2021, 4, 55-62. https://0-doi-org.brum.beds.ac.uk/10.3390/particles4010008

AMA Style

Kryshen E, Ivanishchev D, Kotov D, Malaev M, Riabov V, Ryabov Y. Reconstruction of Photon Conversions in the MPD Experiment. Particles. 2021; 4(1):55-62. https://0-doi-org.brum.beds.ac.uk/10.3390/particles4010008

Chicago/Turabian Style

Kryshen, Evgeny, Dmitry Ivanishchev, Dmitry Kotov, Mikhail Malaev, Victor Riabov, and Yuriy Ryabov. 2021. "Reconstruction of Photon Conversions in the MPD Experiment" Particles 4, no. 1: 55-62. https://0-doi-org.brum.beds.ac.uk/10.3390/particles4010008

Article Metrics

Back to TopTop