Next Article in Journal
Extracting Hypernuclear Properties from the (e, eK+) Cross Section
Next Article in Special Issue
Probing Dense QCD Matter: Muon Measurements with the CBM Experiment at FAIR
Previous Article in Journal / Special Issue
Techniques for Reconstruction of Strange Objects at MPD
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A Vector Finder Toolkit for Track Reconstruction in MPD ITS

by
Dmitry Zinchenko
*,†,
Eduard Nikonov
,
Veronika Vasendina
and
Alexander Zinchenko
Joint Institute for Nuclear Research, Joliot-Curie 6, 141980 Dubna, Moscow Region, Russia
*
Author to whom correspondence should be addressed.
These authors contributed equally to this work.
Submission received: 30 March 2021 / Revised: 24 April 2021 / Accepted: 26 April 2021 / Published: 29 April 2021
(This article belongs to the Special Issue Analysis Techniques and Physics Performance Studies for FAIR and NICA)

Abstract

:
As a part of the future upgrade program of the Multi-Purpose Detector (MPD) experiment at the Nuclotron-Based Ion Collider Facility (NICA) complex, an Inner Tracking System (ITS) made of Monolitic Active Pixel Sensors (MAPSs) is proposed between the beam pipe and the Time Projection Chamber (TPC). It is expected that the new detector will enhance the experimental potential for the reconstruction of short-lived particles—in particular, those containing the open charm particle. To study the detector performance and select its best configuration, a track reconstruction approach based on a constrained combinatorial search was developed and implemented as a software toolkit called Vector Finder. This paper describes the proposed approach and demonstrates its characteristics for primary and secondary track finding in ITS, ITS-to-TPC track matching and hyperon reconstruction within the MPD software framework. The results were obtained on a set of simulated central gold–gold collision events at s N N = 9 GeV with an average multiplicity of ∼1000 charged particles in the detector acceptance produced with the Ultra-Relativistic Quantum Molecular Dynamics (UrQMD) generator.

1. Introduction

The Nuclotron-Based Ion Collider Facility (NICA) and the Multi-Purpose Detector (MPD) at the NICA collider, which are currently under construction at the Joint Institute for Nuclear Research (JINR) in Russia [1,2], are expected to strengthen the experimental investigation of the high-density baryonic matter produced in heavy-ion collisions. New experiments will be performed to collect measurements of the bulk properties of the nuclear matter and study different probes with relatively high production rates. These are also planned during the upgrade of the detector to enhance its capabilities to study rear probes such as multi-strange (anti)hyperons and charmed particles.
The investigation of the formation and propagation of particles containing charm and/or anticharm quarks in heavy-ion collisions opens the possibility of exploring the early phase of the fireball due to the large mass of the charm quarks, which can be produced in hard processes only. The elliptic flow of charm quarks is an important diagnostic aspect of the early fireball and its prevailing degrees of freedom. The propagation of open charm particles in the fireball provides information about the transport properties of hot and dense matter [3,4].
The NICA energy region of s N N = 4 11 GeV for heavy-ion collisions is at the kinematic threshold for the production of charmed mesons D ± , D 0 , D s and baryons Λ c , Ξ c [5]. Therefore, the study of charmed particles in the MPD experiment could provide valuable information on their production mechanisms (such as subthreshold production, for instance).
The MPD upgrade plan foresees the installation of the Inner Tracking System (ITS) based on the next-generation silicon pixel detectors known as Monolitic Active Pixel Sensors (MAPSs) [6] between the beam pipe and the Time Projection Chamber (TPC). Such a detector will increase the research potential of the experiment in terms of both the high-luminosity proton–proton and high-multiplicity nucleus–nucleus interactions. The main purpose of the ITS is to improve the quality and precision of track, primary and secondary vertex reconstruction in MPD in the region close to the interaction point.
In order to optimize the ITS configuration and study the performance of the upgraded detector, a dedicated ITS track reconstruction method based on a constrained combinatorial search was developed and implemented as a Vector Finder toolkit [7,8]. The approach is described below, and some results of its application to track and hyperon decay reconstruction in the MPD ITS are presented.

2. MPD Tracking System Configuration

The MPD detector at the NICA collider has been designed to detect hadrons, electrons and photons over a large phase-space at the high event rate achieved at NICA. At the designed luminosity of L = 1027 cm−2s−1 for Au + Au, the collision rate will be ∼7 kHz.
Figure 1a shows a three-dimensional view of the MPD setup. All the subdetectors are located inside a superconducting solenoid, which produces a magnetic field along the beam axis with a nominal strength of 0.5 T. The detailed description of the MPD detector can be found in [9]. Here, only the tracking system is presented.
The MPD TPC is the main tracking device for the reconstruction of charged particle trajectories over the pseudorapidity range of | η | < 1.5 . It provides the momentum measurement for charged particles with a precision better than 3.5% at a transverse momentum p T below 2 GeV/c and particle identification via the specific energy loss measurement ( d E / d x ) in the TPC gas with a resolution better than 8%. The track reconstruction quality outside the midrapidity region is enhanced by a straw-tube End Cap Tracker (ECT). The Inner Tracker (IT) system, made of five layers of silicon pixel detectors (Figure 1b,c), is intended for the very precise determination of the position of primary and secondary vertices and improvement of the track reconstruction performance for low- p T tracks. These goals are achieved due to the higher acceptance of the ITS at pseudorapidity ( | η | < 2.0 ) and transverse momentum of tracks because of its close proximity to the beam line (∼3 cm instead of ∼40 cm in the TPC) and excellent spatial resolution (<10 μm in each direction as compared with 0.5 mm and 1 mm in transverse and longitudinal directions, respectively, for the TPC).

3. Vector Finder Toolkit Algorithms

The track finding method is based on the combinatorial search for detector hit combinations that can potentially belong to the same particle. Here, a hit is a reconstructed point of a particle trajectory crossing of a detector layer. Possible hit combinations are selected on detector layers using prior constraints defined by the detector geometry and event characteristics, such as the momentum and angular spectra of particles to be reconstructed.
The collected hit combinations, forming sequences of joint vectors connecting pairs of hits on different detector layers, are further fitted by the Kalman filter to verify hit-to-track associations via the track quality values defined by the χ 2 -metrics and to obtain track parameters.

3.1. Primary Track Reconstruction

Primary tracks are produced in the collision point of two particle beams, called the primary vertex. Therefore, the acceptance window determination for hits on each detector layer can utilize hit azimuthal and polar angles in the coordinate system with the Z-axis directed along the beam line, the Y-axis vertically upward and X-axis horizontally.
In the longitudinal projection, the particle track is not affected by the magnetic field, so it can be considered as a straight line. Therefore, the track candidate can be extended on each detector layer by adding hits with sufficiently close polar angle theta (Figure 2).
In the transverse projection, the charged particle track is close to a circle arc because of the magnetic field. Therefore, hits added to a track candidate on each layer should correspond to the estimated track curvature. Thus, the area of interest can be defined by two angles epsphi lower and epsphi higher shifted from the track candidate’s last accepted hit azimuthal angle phi in the direction defined by the particle charge.
The size of the hit acceptance window on each detector layer can be evaluated from the Monte Carlo-simulated event samples. For example, Figure 3 shows the distributions of the angular position differences of hits on adjacent layers from the same particle as functions of p T . From the plots, one can select optimal cuts which are usually taken to be dependent on the processing iteration. Running the reconstruction procedure in several passes allows one to relax the acceptance criteria, with the number of remaining hits decreasing after previous iterations, resulting in the acceleration of the processing. From Figure 3b, one can also see that to determine the proper shift of the transverse angle window, it is necessary to evaluate p T (by taking two first hits and the primary vertex nominal position, for example). More details can be found in [11].

3.2. Secondary Track Reconstruction

Secondary tracks can be produced by decay products of short-lived particles; therefore, they do not pass through the primary vertex. Thus, the primary vertex cannot be used to define good constraints, and the utilization of longitudinal and transverse angles theta and phi is no longer optimal. Thus, an alternative solution should be used; e.g., based on estimated track parameters at the earlier steps of its propagation.
In the longitudinal projection, the particle track is still close to a straight line, because it is not affected by the magnetic field. Thus, linear extrapolation can be used to find possible hits to extend track candidates to detector layers (Figure 4). In the transverse projection, the particle track is close to a circle arc, so the circle arc propagation is used. Initial acceptance windows and p T estimates still can be taken as for the case of primary tracks with some adjustment to account for the approximation.

3.3. Track Matching

Separate track reconstruction procedures in ITS and TPC require a method to combine tracks from both detectors to achieve optimal track reconstruction quality. The implemented track-matching procedure propagates tracks from both detectors to a cylinder surface lying between them to find a set of corresponding TPC tracks for each ITS track, as shown in Figure 5. The best possible match is chosen using a special quality function that takes into account both the number of hits and χ 2 of the resulting track, defined as
q u a l i t y = N h i t s + ( 100.0 m i n ( χ 2 , 100.0 ) ) / 101.0 ,
where N h i t s is the number of hits in the combined track and χ 2 is the total χ 2 value.

4. Results

The algorithms described were implemented as a Vector Finder toolkit within the MpdRoot software framework [12] and tested on a simulated 200-event set of central Au + Au collisions produced with the Ultra-Relativistic Quantum Molecular Dynamics (UrQMD) generator [13]. Produced particles were transported through the detector setup with the Geant3 code [14]. To have better defined conditions for track reconstruction (namely, a cleaner track matching environment), TPC response simulation was performed using a simplified approach based on a Gaussian-smeared-hit approximation; i.e., effects of track multiplicity in the TPC on the results were mostly excluded.
Results for track reconstruction efficiencies as functions of pseudorapidity η and transverse momentum p T can be seen in Figure 6. The efficiency is defined with respect to particles having produced hits in at least three layers of the ITS. One can see that the efficiency remained close to 100% up to the geometric acceptance limit of | η | = 2 and down to ∼0.05 GeV/c in p T . Some decrease for secondary tracks below 0.25 GeV/c can be explained by the particle specie composition; i.e., there was a relatively large fraction of protons in this momentum region. They experienced non-negligible energy loss in the outer two silicon layers with relatively large thickness and changed their curvature in the magnetic field, which was not accounted for in the algorithm, since there was no particle identification at this reconstruction stage. This explanation is supported by earlier results obtained for a detector configuration with all 50 μm layers [11], where no such efficiency decrease was observed.
The track matching efficiency remained above ∼90% up to | η | = 2 for p T -values above ∼0.1 GeV/c (Figure 7).
To demonstrate the effect of the ITS on the short-lived particle decay reconstruction quality, a Monte Carlo event set of 250,000 central Au + Au collisions at s N N = 9 GeV was produced and processed to reconstruct and select Λ and Ω hyperons via their weak decays into p + π and Λ + K , respectively. The simulation and reconstruction were done for two detector configurations—with and without the ITS. For this study, the TPC simulation package included the part realistically describing the detector response to properly take into account high-multiplicity effects [15], although decay products were identified using Monte Carlo truth information. To ensure a fair comparison of the results, the hyperon selection criteria were chosen to obtain the maximum significance of the decay product invariant mass peak [16]. The significance was defined as S / S + B , where the signal S and background B were calculated as integrals within the ± 2 σ interval of the peak position of the signal and background distributions, obtained from the invariant mass fit to a sum of gaussian (signal) and polynomial (background) functions. From Figure 8, one can see that the ITS gives a significant improvement in hyperon reconstruction quality, especially for objects with a shorter decay length (decay length c τ = 7.9 cm for Λ and 2.5 cm for Ω ), and the presented track reconstruction approach improves the results as compared with the method based on a Kalman filter extension of tracks from the TPC to the ITS (Figure 8a).

5. Conclusions

A track reconstruction approach based on a constrained combinatorial search was developed and implemented as a software toolkit called Vector Finder for the optimization of studies of the future inner tracking system of the MPD experiment.
It demonstrated good performance results for simulated high-multiplicity events of Au + Au collisions, showing improvements in the reconstruction quality of short-lived particles over the existing method.
The role of the ITS for future studies of short-lived particles has been also demonstrated, and this research direction will be further followed by evaluating the reconstruction feasibility of open charm particles.

Author Contributions

Conceptualization, A.Z.; Investigation, V.V.; Software, D.Z.; Supervision, E.N. 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-40060.

Acknowledgments

The authors would like to thank V. Kondratev (SPbSU) and Yu. Murin (JINR) for providing the description of the MPD ITS geometry.

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. Kekelidze, V.D.; Lednicky, R.; Matveev, V.A.; Meshkov, I.N.; Sorin, A.S.; Trubnikov, V.G. Three stages of the NICA accelerator complex. Eur. Phys. J. A 2016, 52, 211. [Google Scholar] [CrossRef]
  2. 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]
  3. Buss, O.; Gaitanos, T.; Gallmeister, K.; van Hees, H.; Kaskulov, M.; Lalakulich, O.; Larionov, A.B.; Leitner, T.; Weil, J.; Mosel, U. Transport-theoretical Description of Nuclear Reactions. Phys. Rep. 2012, 512, 1. [Google Scholar] [CrossRef] [Green Version]
  4. Friman, B.; Hohne, C.; Knoll, J.; Leupold, S.; Randrup, J.; Rapp, R.; Senger, P. The CBM Physics Book: Compressed Baryonic Matter in Laboratory Experiments; Springer: Berlin/Heidelberg, Germany, 2011; Volume 814, p. 1. [Google Scholar]
  5. Chen, H.X.; Chen, W.; Liu, X.; Liu, Y.R.; Zhu, S.L. A review of the open charm and open bottom systems. Rep. Prog. Phys. 2017, 80, 076201. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  6. Abelev, B.; Adam, J.; Adamová, D.; Aggarwal, M.M.; Rinella, G.A.; Agnello, M.; Agostinelli, A.; Agrawal, N.; Ahammed, Z.; Ahmad, N.; et al. Technical Design Report for the Upgrade of the ALICE Inner Tracking System. J. Phys. G 2014, 41, 087002. [Google Scholar] [CrossRef] [Green Version]
  7. Zinchenko, D.; Nikonov, E.; Zinchenko, A. A track finding algorithm for the inner tracking system of MPD/NICA. EPJ Web Conf. 2019, 204, 07006. [Google Scholar] [CrossRef] [Green Version]
  8. Zinchenko, D.; Zinchenko, A.; Nikonov, E. A “vector finder” approach to track reconstruction in the inner tracking system of MPD/NICA. AIP Conf. Proc. 2019, 2163, 060006. [Google Scholar]
  9. Abraamyan, K.U.; Afanasiev, S.V.; Alfeev, V.S.; Anfimov, N.; Arkhipkin, D.; Aslanyan, P.Z. The MPD detector at the NICA heavy-ion collider at JINR. Nucl. Instrum. Methods Phys. Res. Sect. A 2011, 628, 99. [Google Scholar] [CrossRef]
  10. Zinchenko, A.I.; Igolkin, S.N.; Kondratiev, V.P.; Murin, Y.A. NICA-MPD Vertex Tracking Detector Identification Capability for Reconstructing Strange and Charmed Particle Decays. Phys. Part. Nucl. Lett. 2020, 17, 856. [Google Scholar] [CrossRef]
  11. Zinchenko, D.; Zinchenko, A.; Nikonov, E. Vector Finder—A Toolkit for Track Finding in the MPD Experiment. Phys. Part. Nucl. Lett. 2021, 18, 107. [Google Scholar] [CrossRef]
  12. Gertsenberger, K.; Merts, S.; Rogachevsky, O.; Zinchenko, A. Simulation and analysis software for the NICA experiments. Eur. Phys. J. A 2016, 52, 214. [Google Scholar] [CrossRef]
  13. Bleicher, M.; Zabrodin, E.; Spieles, C.; Bass, S.A.; Ernst, C.; Soff, S.; Bravina, L.; Belkacem, M.; Weber, H.; Stoecker, H.; et al. Relativistic hadron hadron collisions in the ultrarelativistic quantum molecular dynamics model. J. Phys. G 1999, 25, 1859. [Google Scholar] [CrossRef]
  14. Brun, R.; Bruyant, F.; Maire, M.; McPherson, A.C.; Zanarini, P. GEANT3. CERN 1987. Available online: https://cds.cern.ch/record/1119728 (accessed on 30 March 2021).
  15. Kolesnikov, V.; Mudrokh, A.; Vasendina, V.; Zinchenko, A. Towards a Realistic Monte Carlo Simulation of the MPD Detector at NICA. Phys. Part. Nucl. Lett. 2019, 16, 6. [Google Scholar] [CrossRef]
  16. Ilieva, M.; Kolesnikov, V.; Murin, Y.; Suvarieva, D.; Vasendina, V.; Zinchenko, A.; Litvinenko, E.; Gudima, K. Evaluation of the MPD Detector Capabilities for the Study of the Strangeness Production at the NICA Collider. Phys. Part. Nucl. Lett. 2015, 12, 100. [Google Scholar] [CrossRef] [Green Version]
Figure 1. (a) Three-dimensional view of the MPD detector. (b) Beam pipe, TPC and ITS geometry representation in the simulation package. (c) Technical drawing of the ITS cut view [10]: shown are the beam pipe, three inner layers with an effective thickness of 50 μm of Si and two outer layers with an effective thickness of 700 μm of Si. Average layer radii are 24, 43, 63, 146 and 196 mm.
Figure 1. (a) Three-dimensional view of the MPD detector. (b) Beam pipe, TPC and ITS geometry representation in the simulation package. (c) Technical drawing of the ITS cut view [10]: shown are the beam pipe, three inner layers with an effective thickness of 50 μm of Si and two outer layers with an effective thickness of 700 μm of Si. Average layer radii are 24, 43, 63, 146 and 196 mm.
Particles 04 00017 g001
Figure 2. (a) Longitudinal projection scheme for a primary track. (b) Transverse projection scheme for a primary track.
Figure 2. (a) Longitudinal projection scheme for a primary track. (b) Transverse projection scheme for a primary track.
Particles 04 00017 g002
Figure 3. (a) Longitudinal angle difference of hits on sequential layers dtheta as a function of p T for primary tracks. (b) Transverse angle difference of hits on sequential layers dphi as a function of signed p T for primary tracks. Different colors correspond to different distances between detector layers.
Figure 3. (a) Longitudinal angle difference of hits on sequential layers dtheta as a function of p T for primary tracks. (b) Transverse angle difference of hits on sequential layers dphi as a function of signed p T for primary tracks. Different colors correspond to different distances between detector layers.
Particles 04 00017 g003
Figure 4. (a) Longitudinal projection scheme for a secondary track. (b) Transverse projection scheme for a secondary track.
Figure 4. (a) Longitudinal projection scheme for a secondary track. (b) Transverse projection scheme for a secondary track.
Particles 04 00017 g004
Figure 5. TPC and ITS track matching scheme.
Figure 5. TPC and ITS track matching scheme.
Particles 04 00017 g005
Figure 6. (a) Track reconstruction efficiency as a function of | η | for primary and secondary tracks. (b) Efficiency vs p T for primary and secondary tracks. Symbols and solid lines correspond to primary and secondary tracks, respectively.
Figure 6. (a) Track reconstruction efficiency as a function of | η | for primary and secondary tracks. (b) Efficiency vs p T for primary and secondary tracks. Symbols and solid lines correspond to primary and secondary tracks, respectively.
Particles 04 00017 g006
Figure 7. (a) Track matching efficiency vs | η | for primary and secondary tracks. (b) Matching efficiency vs p T for primary and secondary tracks. Notations as in Figure 6.
Figure 7. (a) Track matching efficiency vs | η | for primary and secondary tracks. (b) Matching efficiency vs p T for primary and secondary tracks. Notations as in Figure 6.
Particles 04 00017 g007
Figure 8. (a) Reconstructed invariant mass spectra of Λ hyperon decay products: Λ p + π . Red, blue and hatched histograms show results from the Vector Finder, Kalman filter track follower and TPC-only reconstruction, respectively. (b) Reconstructed invariant mass spectra of Ω hyperon decay products: Ω Λ + K . Filled and empty symbols present results for TPC and TPC–ITS configurations. Histograms are fitted to a sum of a gaussian and a polynomial functions. The numbers in the legends present the mean and sigma of the gaussian as well as the signal/background ratio, significance and signal reconstruction efficiency, and efficiency is taken with respect to all hyperons produced within 50 cm of the interaction point. Left and right legends are, respectively, for TPC and TPC–ITS.
Figure 8. (a) Reconstructed invariant mass spectra of Λ hyperon decay products: Λ p + π . Red, blue and hatched histograms show results from the Vector Finder, Kalman filter track follower and TPC-only reconstruction, respectively. (b) Reconstructed invariant mass spectra of Ω hyperon decay products: Ω Λ + K . Filled and empty symbols present results for TPC and TPC–ITS configurations. Histograms are fitted to a sum of a gaussian and a polynomial functions. The numbers in the legends present the mean and sigma of the gaussian as well as the signal/background ratio, significance and signal reconstruction efficiency, and efficiency is taken with respect to all hyperons produced within 50 cm of the interaction point. Left and right legends are, respectively, for TPC and TPC–ITS.
Particles 04 00017 g008
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Zinchenko, D.; Nikonov, E.; Vasendina, V.; Zinchenko, A. A Vector Finder Toolkit for Track Reconstruction in MPD ITS. Particles 2021, 4, 186-193. https://0-doi-org.brum.beds.ac.uk/10.3390/particles4020017

AMA Style

Zinchenko D, Nikonov E, Vasendina V, Zinchenko A. A Vector Finder Toolkit for Track Reconstruction in MPD ITS. Particles. 2021; 4(2):186-193. https://0-doi-org.brum.beds.ac.uk/10.3390/particles4020017

Chicago/Turabian Style

Zinchenko, Dmitry, Eduard Nikonov, Veronika Vasendina, and Alexander Zinchenko. 2021. "A Vector Finder Toolkit for Track Reconstruction in MPD ITS" Particles 4, no. 2: 186-193. https://0-doi-org.brum.beds.ac.uk/10.3390/particles4020017

Article Metrics

Back to TopTop