KnE Energy | The 3rd International Conference on Particle Physics and Astrophysics (ICPPA) | pages: 154–161

1. Open charm measurements in NA61/SHINE

The SPS Heavy Ion and Neutrino Experiment (NA61/SHINE) [1] at CERN was designed for studies of the properties of the onset of deconfinement and search for the critical point of strongly interacting matter. These goals are being pursued by investigating p+p, p+A and A+A collisions at different beam momenta from 13 A to 158 A GeV/c for ions and up to 400 GeV/c for protons.

Heavy flavor quarks are produced at the early stage of the hadron collision and may serve as sensitive probes for properties of the QGP. However, only indirect measurements of open charm production in nucleus-nucleus collisions at the SPS energies exist, but precise data on both open and hidden charm production are required to provide information on the in-medium behavior of quarkonia in a model independent way [2].

For open charm measurements in nucleus-nucleus collisions at SPS energies, the experiment NA61/SHINE was upgraded with the new Small Acceptance Vertex Detector (SAVD) [3,4,5]. The low yields of D mesons and their short lifetime require detectors located close to the primary vertex that are capable of providing high tracking efficiency and high primary and secondary vertex resolution.

Since D0 mesons typically decay tens to hundreds of microns downstream of their production point, one needs to reconstruct the decay vertex with a precision of about 50 μ m in order to recognize daughter particles that originate from D0 decays.

Figure 1

The layout of the experimental setup of NA61/SHINE.

Figure 2

Photograph of the SAVD before closing the detector with the front and exit windows.


2. Small Acceptance Vertex Detector SAVD for NA61/SHINE

The requirements of high precision and resolution can only be reached by adding Vertex Detector to the NA61/SHINE set-up. The constructed SAVD is positioned between the target and the VTPC-1 (see Fig. figFigure1). It consists of two arms (called Jura and Saleve) with four planes (stations) of coordinate-sensitive detectors located at 5, 10, 15 and 20 cm distance from the target. High coordinate resolution MIMOSA-26 sensors [6] based on the CMOS Pixel Sensor technology were chosen as the basic detection element of the SAVD stations. The sensors provide a spatial resolution of 3.5 μ m, have very low material budget (50 μ m thickness), and their readout time is 115.2 μ s.

In order to reduce beam-gas interactions, the SAVD and the target were placed in a helium enclosure. The magnetic field in the SAVD volume is small and inhomogeneous (0.13 - 0.25T) because the detector is located close to the edge of the VTPC-1 magnet.

Figure figFigure2 shows the fully integrated detector in the helium vessel installed in December 2016 in the NA61/SHINE experiment for data taking with Pb+Pb collisions at 150 A GeV/c [7]. The main goal of the test was to prove precise tracking in the large track multiplicity environment, to demonstrate the ability of precise primary vertex reconstruction and eventually to reconstruct a D0 signal from the collected data.

3. Data reconstruction

Track and vertex reconstruction

A straight line was chosen as the track model: x(z)=Az+x0,y(z)=Bz+y0 .

The first step of track reconstruction uses a combinatorial method of track recognition based on checking of combinations of all hits from different stations. If the hits lie on a straight line according to a χ2 criterion, the combination is accepted as a reconstructed track.

From the distributions of the residuals of hits from the reconstructed tracks the spatial sensor resolution was determined to have values of σx = 4.7 μ m and σy = 5.0 μ m.

The angular distribution of the reconstructed tracks in the x,z plane is shown in Fig. figFigure3. The distribution for each arm a has clear three-peak structure for each arm. The narrow inner most peak (small angles) was associated with particles produced far upstream and travelling parallel to the beam for a long distance. The middle peak corresponds to particles produced upstream of the target. The outer peak is created by particles produced on the Pb target; these tracks are used for further analysis.

Figure 3

Angular distribution of reconstructed tracks in the x,z plane (x-slope) for Jura (positive values) and Saleve (negative values) arms.


The primary vertex is defined as the point of the closest convergence of all reconstructed tracks. Thus, the longitudinal coordinate of the primary vertex is found by minimizing the expression:


which describes the sum of the relative distances of all track pairs reconstructed in a single event at the given transverse plane defined by the longitudinal coordinate z . The xprim and yprim coordinates are afterwards calculated as the average of x and y positions of tracks at z=zprim .

Primary vertex reconstruction was done separately for two samples of tracks with x-slope in the interval 0.02 to 0.04 and with x-slope > 0.04. By looking at the longitudinal distribution of the primary vertex for these samples of tracks (the left picture in Fig. figFigure4) one may see that indeed the tracks associated with the most outer peak in Fig. figFigure3 (blue) originate from the target which is located about 49 mm upstream the fist VD station; the primary vertices associated with tracks from the middle peak (red in Fig. figFigure3) are relatively smoothly distributed upstream of the target in the range from -1200 mm to -55 mm. At -190 mm the distribution has a sharp peak which is related to interactions in the aluminized Mylar front window of the SAVD. One can also see that between the window and the target the frequency of interaction drops due to the presence of helium gas in the SAVD vessel.

The distribution of zprim from the interactions in the target is shown by the right panel of Fig. figFigure4. The near rectangular shape of the distribution allows to reconstruct the thickness of the target (= 1 mm) and its location with respect to the SAVD. The primary vertex resolution is σx = 5 μ m, σy = 1.5 μ m and σz = 30 μ m (for more details see [7]). The difference between σx and σy is caused by the presence of a vertical component of the magnetic field in the SAVD volume.

Figure 4

Left: Distribution of the longitudinal coordinate zprim of primary vertices (see text). Right: Distribution of zprim for tracks produced in the target.


After the primary vertex is found, the next step of track reconstruction searches for 3-hit tracks using the Hough transform method (for more details see [8]). It uses a parametric description of the track by a set of parameters. Once a parametric description of the track in the coordinate space is given, the hit position measurements are transformed into the track parameter space. Then, in this so-called accumulation space track candidates are obtained as local maxima using the voting procedure: the most popular candidates are accepted as real tracks.

The magnetic field in the Vertex Detector volume allows momentum reconstruction of the SAVD-tracks [8]. The following formula was used to calculate momentum:


where α1 , α2 are angles of the track in the first and last stations respectively. Momentum reconstruction is imprecise due to the fact that the field is weak and inhomogeneous ( Bydl = 0.04 T · m). Nevertheless, the obtained values are useful for verifying track matching between SAVD-tracks and TPC-tracks.

Track matching

The track matching algorithm tries all possible combinations of pairs of SAVD with TPC tracks, assigning to each SAVD track the total momentum of the paired TPC track (which is precise). Next, the SAVD tracks and TPC tracks are extrapolated to a common surface. The algorithm then searches for the best match among all combinations based on the difference in the position and direction of the tracks and uses charge and total momentum to verify the match (since the SAVD has small reconstruction power for charge and momentum). About 80% of SAVD tracks are in the acceptance of the TPCs according to GEANT4 simulations. The track matching efficiency is determined by number of matched SAVD tracks divided by number of SAVD tracks in the TPC acceptance, and is on the level of 90%.

4. Reconstruction of D0 signal

SAVD tracks matched to TPC tracks are used to search for the D0 signal. In the current analysis PID information was not used. Each SAVD track is paired with another SAVD track and is assumed to be either a kaon or pion. Thus each pair contributes twice in the combinatorial invariant mass distribution. The combinatorial background is several orders of magnitude higher than the D0+D0¯ signal due to the low yield of charm particles. In order to reduce the large background and maximize the signal to noise ratio (SNR) of the reconstructed D0 peak, four cuts are applied, which were determined from the simulations [4,5]:

  • A cut on the track transverse momentum pT>0.34 GeV/c;

  • A cut on the track impact parameter d>34μ m;

  • A cut on the longitudinal distance between the D0 decay candidate and the interaction point Vz>475μ m;

  • A cut on the impact parameter D of the back-extrapolated D0 candidate momentum vector D<21μ m.

The left panel of Fig. figFigure5 shows the invariant mass distribution of unlike charge daughter candidates with applied cuts. One observes a small peak emerging at 1.8 GeV/ c2 , which could be an indication of D0 production. The invariant mass distribution was fitted (red line) using of a third order polynomial function to describe the background and a Gaussian peak for the D0 signal contribution. The indicated errors are statistical only. From the fit one finds the width of the peak to be 40 ± 15 MeV and the total yield to amount to 55 ± 20 with a SNR of 3.

To test the detector capabilities and validate the D0 peak extraction procedure, it was also attempted to reconstruct the Ks0 , which is much more abundantly produced. However, one can reconstruct only a small fraction of all Ks0 with the SAVD since the detector design and the tracking algorithm were optimized for decays which have much shorter decay length. A clear peak corresponding to Ks0 decays is seen in the right panel of Fig. figFigure5 with a width of 18 ± 5 MeV and a total yield of 130 ± 30 with a SNR of about 4.5.

Figure 5

Left: Invariant mass distribution of unlike charge π,K decay track candidates. Right: Invariant mass distribution of π-,π+ decay track candidates. See text for details.


5. Summary

A small-acceptance version of the NA61/SHINE Vertex Detector (SAVD) was installed in 2016. A modest set of data of Pb+Pb collisions at 150 A GeV/c was collected and analyzed. A weak indication of a D0 signal was observed in the π,K decay channel. However, the data set is not finally calibrated, and further optimization of tracking, track matching and analysis algorithms, as well as of cuts is still ongoing. It may be premature to draw a final conclusion from the analysis result, though the result appears promising.

The SAVD will be used in p+Pb and Xe+La data taking in the autumn of 2017. Also, high statistics measurements in Pb-Pb collisions at 150 A GeV/c are planned in 2018. Looking forward, an upgraded version the so-called Large Acceptance Vertex Detector (LAVD) with more sensors is being planned; the exact design of this detector is currently under investigation.


The work related to instrumentation was supported by the Polish National Center for Science grant 2014/15/B/ST2/02537 and the work related to reconstruction and data analysis was partially supported by the Russian Science Foundation research grant 16-12-10176.



Abgrall N et al. (NA61 Collaboration), 2014, JINST 9, P06005.


H. Satz Calibrating the in-medium behavior of quarkoniaAdvances in High Energy Physics201320132-s2.0-8488148577010.1155/2013/242918242918


NA61/SHINE Collaboration, 2015, Addendum (Proposal), CERN-SPSC-2015-038, SPSC-P-330-ADD-8.


Y. Ali P. Staszel A. Marcinek J. Brzychczyk R. Płaneta Feasibility studies of open charm measurements with the NA61/shine experiment at cern-SPSActa Physica Polonica B20134410201920342-s2.0-8489001369410.5506/APhysPolB.44.2019


Y. Ali P. Staszel NA61/SHINE experiment upgrade with vertex detector for open charm measurementsJournal of Physics: Conference Series201450912-s2.0-8490162617210.1088/1742-6596/509/1/012083012083


Deveaux M et al., 2011, Journal of Instrumentation 6(02), C02004.


NA61/SHINE Collaboration, 2017, CERN-SPSC-2017-038 (SPSC-SR-221).


A. Merzlaya Track reconstruction in the inhomogeneous magnetic field for Vertex Detector of NA61/SHINE experiment at CERN SPSJournal of Physics: Conference Series201779812-s2.0-8501628773910.1088/1742-6596/798/1/012072012072



  • Downloads 17
  • Views 63



ISSN: 2413-5453