A high-resolution pixel silicon Vertex Detector for open charm measurements with the NA61/SHINE spectrometer at the CERN SPS

The study of open charm meson production provides an efficient tool for the investigation of the properties of hot and dense matter formed in nucleus–nucleus collisions. The interpretation of the existing di-muon data from the CERN SPS suffers from a lack of knowledge on the mechanism and properties of the open charm particle production. Due to this, the heavy-ion programme of the NA61/SHINE experiment at the CERN SPS has been extended by precise measurements of charm hadrons with short lifetimes. A new Vertex Detector for measurements of the rare processes of open charm production in nucleus–nucleus collisions was designed to meet the challenges of track registration and high resolution in primary and secondary vertex reconstruction. A small-acceptance version of the vertex detector was installed in 2016 and tested with Pb + Pb collisions at 150AGeV/c\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$A\,\text{ GeV }\!/\!c$$\end{document}. It was also operating during the physics data taking on Xe + La and Pb + Pb collisions at 150AGeV/c\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$A\,\text{ GeV }\!/\!c$$\end{document} conducted in 2017 and 2018. This paper presents the detector design and construction, data calibration, event reconstruction, and analysis procedure.


Introduction
The charm production mechanism is one of the important questions in relativistic heavy-ion physics.Several models were introduced to describe charm production.Some are based on dynamical and others -on statistical approaches.Predictions of these models on the mean number of produced cc pairs ( cc ) for central Pb+Pb collisions at 158A GeV/c differ by up to a factor of 50 [1,2].Moreover, the system size dependence is different in these approaches and the the predictions suffer from large systematic uncertainties [3,4].Precise data on cc will allow to disena Corresponding author: pawel piotr.staszel@uj.edu.pl.
tangle between theoretical predictions and learn about the charm quark and hadron production mechanism.Obtaining good estimate of cc requires measurements of D 0 , D + and their antiparticles.This is because these mesons carry about 85% of the total produced charm in Pb+Pb collisions at the top SPS energy [5,6].
Besides this, a study of open charm meson production was proposed as a sensitive tool for detailed investigations of the properties of hot and dense matter formed in nucleus-nucleus collisions at ultra-relativistic energies [7,8,9].In particular, charm mesons are of vivid interest when studying the phase transition between confined hadronic matter and the quark-gluon plasma (QGP).The cc pairs produced in the collisions are converted into open charm mesons and charmonia (J/ψ mesons and their excited states).The charm production is expected to be different in the confined and deconfined matter because of the different properties of charm carriers in these phases.In confined matter, the lightest charm carriers are D mesons, whereas, in deconfined matter, the carriers are charm quarks.The production of a DD pair (2m D = 3.7 GeV) requires more energy than the production of a cc pair (2m c = 2.6 GeV).Since the effective degrees of freedom of charm, hadrons and charm quarks are similar [10], more abundant charm production is expected in deconfined than confined matter.Consequently, in analogy to strangeness [3,11], a change in collision energy dependence of cc production may indicate an onset of deconfinement.
Finally, systematic measurements of open charm production are urgently needed to interpret existing results on J/ψ.Such measurements would allow disentangling between initial and final state effects, revealing hidden and open charm transport properties through the dense medium created in nucleus-nucleus collisions and testing the validity of theoretical models [8].
Measurements of open charm mesons are challenging since the yields of D mesons are low, and their lifetimes are relatively short (cτ = 122 µm).The measurements require precise tracking and high primary and secondary vertex resolutions.To meet these challenges, a novel highresolution Small Acceptance Vertex Detector (SAVD) was designed and built under the leadership of the Jagiellonian University group participating in the NA61/SHINE experiment.SAVD was installed as a part of the NA61/SHINE facility in December 2016.Test data on Pb+Pb collisions at 150A GeV/c beam momenta were collected and analyzed.The main goal of the test was to prove the feasibility of precise tracking in the large track multiplicity environment and demonstrate the ability of precise primary and secondary vertex reconstruction.In 2017 and 2018, data on Xe+La and Pb+Pb collisions at the beam momenta of 150A GeV/c were recorded with SAVD included in the detector setup.The data quality and statistics were sufficient for the first direct observation of a D 0 + D 0 signal in the π + K decay channel in nucleus-nucleus collisions at the SPS energy.This paper presents SAVD design and construction, data calibration, event reconstruction, and analysis procedure.
It is foreseen that the NA61/SHINE Collaboration will perform large statistics measurements after 2022.These data will allow for the first insight into the centrality dependence of open charm [4].
The following variables and definitions are used in this paper.The particle rapidity y is calculated in the nucleonnucleon collision center of the mass system (c.m.s.) with where E and p L are the particle energy and longitudinal momentum, respectively, the transverse momentum is denoted as p T , and m is the particle mass.The quantities are given either in GeV or in MeV.The results shown in this paper were obtained for Xe+La collisions at the beam momenta of 150A GeV/c.

NA61/SHINE experimental facility
The SPS Heavy Ion and Neutrino Experiment (NA61/ SHINE) [26] at CERN was designed to study the properties of the onset of deconfinement and search for the critical point of the strongly interacting matter.These goals are being pursued by investigating p+p, p+A and A+A collisions at different beam momenta from 13A to 158A GeV/c for ions and up to 400 GeV/c for protons.
The layout of the experimental setup is shown in Fig. 1.The setup includes the beam position detectors (BPD), Cherenkov counters and the scintillator detectors located upstream of the target.They provide information on the timing, charge and position of beam particles.Further, the experiment includes two Vertex Time Projection Chambers (VTPC-1 and VTPC-2) located inside the vertex magnets, two main TPCs (MTPC-L and MTPC-R) for dE / dx measurements and Gap TPC and Forward TPCs that complete the coverage between MTPCs.These TPCs provide acceptance in the full forward hemisphere, down to p T = 0.The TPCs allow tracking, momentum and charge determination, and measuring the mean energy loss per unit path length.The time-of-flight (ToF) walls used for additional particle identification are located behind the main TPCs.The projectile spectator detector (PSD) measures the energy of the projectile spectator and delivers information on the collision centrality.

Vertex Detector rationale
For open charm measurements in nucleus-nucleus collisions, NA61/SHINE was upgraded with SAVD.As was already mentioned, open charm mesons are difficult to measure because of their low yields and short lifetime.They can be measured in their decay channels into pions and kaons.However, in heavy-ion collisions, pions and kaons are produced in large numbers in other processes giving huge combinatorial background.To distinguish the daughter particles of D 0 mesons from hadrons produced directly in the nucleus-nucleus interaction, one selects hadron pairs created in secondary vertices.The vertex reconstruction is done by extrapolating the track trajectories back to the target and identifying intersection points.The primary vertex will appear as the intersection point of multiple tracks while the tracks originating from selected decays will intersect at the displaced point (secondary vertex), see Fig. 2.
Until the development of silicon sensors for particle tracking, it was not possible to perform secondary vertex reconstruction with resolution sufficient to measure open charm.Consequently, an open charm meson production

SAVD hardware
SAVD is positioned between the target and VTPC-1 (see Fig. 1) in the in-homogeneous and weak (0.13 -0.25T) fringe field of the VTPC-1 magnet.A photograph of the device is shown in Fig. 3.It consists of two arms called Jura and Saleve arm.This naming follows the NA61/ SHINE convention for the left and right partition of the experiment in the direction of the beam, respectively, and corresponds to the location of the nearby mountains.SAVD is composed of four detection planes (stations) equipped with the position-sensitive MIMOSA-26AHR CMOS Monolithic Active Pixel Sensors (MAPS) [28,29,30] provided by the PICSEL group of the IPHC Strasbourg.The arms are horizontally movable, allowing the sensors to be placed safely during beam tuning.The stations, called Vds1, Vds2, Vds3 and Vds4, are located 5, 10, 15 and 20 cm downstream the target, respectively.The sensors are held and water-cooled by vertically oriented ALICE ITS carbon fibre support "ladders" [31] developed by St. Petersburg State University and CERN.The ladders are mounted in C-frames made from aluminum.The four C-frames of each arm share a movable support plate.The first (Vds1) and second station (Vds2) consist of two ladders, each holding one sensor only, the third station consists of two ladders, each holding two sensors, and the last station is composed of four ladders, each hosting two sensors (see Fig. 8).A holder for targets is placed on an additional, movable support.
The whole structure is installed on a thick aluminum base plate, which provides mechanical stability.Four brass screws serve as legs for the plate and enable fine adjustment of the vertical position when installed on the beamline.The pink color box structure in the photograph is made of plexiglass covered with conducting paint.The base plate, together with the plexiglass structure and front and back mylar windows (dismounted on the photograph) served as a gas-tight detector box.During data taking, the detector box is filled with helium gas at atmospheric pressure, which reduces beam-gas interactions and unwanted multiple Coulomb scatterings between the target and sensors.
The readout of the sensors was done via 20 cm long, copper-based single-layer Flex Print Cables (FPC).The non-shielded cables were chosen to minimize the material in the acceptance of the TPC, knowing that they may inject pick-up noise into the sensors.

Sensor technology and integration
The MIMOSA-26AHR sensors have a 1.06 × 2.13 cm 2 sensitive area, which is covered by 1156 columns made of 576 pixels giving 663.5k pixels per chip.The pixel pitch is 18.4 µm in each direction, which leads to an excellent spatial resolution of 4.5 µm.The sensor readout is done with a column-parallel rolling shutter.The readout time is equivalent to the time resolution of the device and amounts to 115.2 µs.The slow control of the sensors is done via a JTAG interface, and the most relevant voltages are generated with internal DACs.A prominent exception to this rule is the so-called clamping voltage, which has to be provided from an external source and sets the dark output signal of the pixels.The sensor performs internal sig-  The 50 µm thin sensors are flexible and initially slightly bent.Their integration was carried out at the Institute of Nuclear Physics (IKF) of the Goethe-University Frankfurt am Main.The sensors were first glued together with the flex print cable to a base plate made from carbon fibre.This base plate is used as a mechanical adapter.It is needed as the sensor, and cable size exceeds the ITS ladder's width.After gluing, the bending of the sensor was eliminated, and it was wire bonded to the FPC.Finally, the base plate was glued on the ladder structure.A photograph of the module obtained is shown in Fig. 4. The estimated average material budget of the module in its active area amounts ∼ 0.3% X 0 .

The DAQ system of SAVD
A schematic diagram of the local SAVD DAQ is depicted in Fig. 5.It relies on hardware and software modules, which were initially developed for the prototype of the CBM Micro Vertex Detector [32] and adapted to the needs of SAVD.
The sensors are connected with the FPCs to a Front End Boards (FEB) are located outside of the acceptance on the C-frames.The FEB boards perform noise filtering.A conventional flat cable connects the FEBs with the socalled converter boards located at the outer side of the box.The converter boards host remote-controlled voltage regulators.Moreover, the boards host a latch-up protection system.This system monitors the bias currents of the sensors and can detect possible over-currents as caused by a latch-up.If a latch-up is detected, a rapid power cycle on a given sensor is enforced to extinguish the related meta stable short circuit.
The sensors are steered and read out by two TRBv3 FPGA boards [34].The standard TDC firmware of these boards was replaced by a dedicated code for steering MIMOSA-26AHR sensors.Hereafter, each board serves a readout of eight sensors (data produced in each arm).During the 2016 test run, the two boards were operated with independent clocks.Consequently, the data was synchronized based on the global trigger of NA61/SHINE only.Starting from 2017, the boards operated on a common clock, and the sensors remained also synchronized in hardware.
The sensors and the TRBv3 boards operate continuously and stream out their data with the UDP protocol through the gigabit-Ethernet interface to a DAQ-PC.To synchronize the data with the trigger of NA61/SHINE, the TRBv3 boards receive the trigger signal via the converter board.Information on the arrival time of the trigger is added in real-time to the data stream, but for the sake of simplicity, the data selection is performed in software on the DAQ-PC.Five sensor frames per trigger were forwarded to the central DAQ after the selection was performed, all other data was rejected.The DAQ-PC also performs basic checks on data integrity.In the case of inconsistencies suggesting sensor malfunctioning, a sensor reset is scheduled and the necessary reprogramming of the sensors via the JTAG interface is performed during the next spill break.
The central NA61/SHINE DAQ runs in a data push mode.To prevent mixing events with different trigger numbers, each subsystem must deliver a busy logic signal.If any of the detector's busy logic lines are asserted, the whole system is halted.If this waiting time surpasses the delay limit, data acquisition is stopped, and all subsystems run through a restart procedure.The SAVD busy signal is generated by its local DAQ program using an external Arduino board.

Detector performance and event reconstruction
The Vertex Detector was designed for high-efficiency tracking and finding of primary and secondary vertices with high resolution.The detector concept was developed based on simulations [35,36,37].The goal was to keep the number of sensors low while requiring the system covers most of the produced open charm mesons.
For studying the detector efficiency and acceptance, the simulations were performed using the Geant4 package (for more details, see Ref. [38]).The background was described using the AMPT model [39] and for the parametrization of the open charm meson spectra, the AMPT and the PHSD models were used.Figure 6 presents the distribution of transverse momentum -rapidity of all generated D 0 and D 0 and those D 0 and D 0 , that pass the detector acceptance, i.e. when both of the daughter tracks have sufficient for reconstruction number of SAVD (3 or more) and TPC (10 or more) hits.The simulation for Xe+La collisions at 150A GeV/c shows that about 7.8 % and 5.9 % acceptance of D 0 +D 0 in π and K decay channel for AMPT and PHSD phase space distributions, respectively.

Sensor operation and efficiencies
In SAVD, the sensors are located as close as 3 mm from the beam center.Thus they are exposed to primary beam ions from the beam halo and nuclear beam fragments.It was considered that the related impacts would create latchup and do severe damage to the sensors.Fortunately, although the beam halo ranged to 1 cm from the beam axis, this was not the case.The ion impacts were observed to create clusters. of the size up to 200 pixels, but no sensor was destroyed by the radiation during the detector operation.This was certainly a success of the related protection system and reflected the unexpectedly good robustness of the sensors.
A dedicated radiation test has shown that 30A GeV/c Pb ions created an integrated, non-ionizing radiation damage of 300 n eq /cm 2 (upper limit).As expected by our radiation dose estimates, the radiation damage in the sensors remained below the radiation tolerance of the sensor, which amounts to ∼ 150 krad and 10 13 n eq /cm 2 at modest cooling (typically the coolant temperature was chosen with 10 • C).
Due to a lack of resources, no near-time monitoring providing a sensor detection efficiency was available during the data taking.The thresholds of the sensors were thus lowered until the highest reasonable dark occupancy of ∼ 10 −4 was reached.Based on the sensor's known efficiency/dark occupancy curve, we expected to reach a good efficiency.However, disappointing efficiencies of 10 − 94 % were observed in the 2016 Pb test run, and two sensors did not work.This lack of efficiency was dominantly caused by a bad synchronization of the data selected by the trigger, which rejected valid data in some cases.This was corrected for the 2017 Xe-La run.Moreover, the biasing voltages were adapted for the nominal settings to account for the ohmic losses in the FPCs.Still, the impact on the clamping voltage had not been considered properly.This issue generated a saturation of the pre-amplifiers of multiple pixels.Once identified, it was corrected by adapting a reference voltage of the pre-amplifiers by slow control.
Thanks to the modifications and sensor repairs, all sensors were operational in the 2017 Xe+La run.Unfortunately, the above-mentioned coarse approach for threshold tuning had to be used again.Still, an efficiency between 84 % and the nominal > 99 % was observed and most sensors showed an efficiency significantly above 90 %.

SAVD internal geometry calibration
The alignment of SAVD was done using track candidates found by the combinatorial method with data taken with zero magnetic fields.The purpose of geometry tuning is to  find the corrections for the sensor positions (each sensor has 6 degrees of freedom: offsets from the nominal geometry in x, y and z position and rotation along x, y and z axes).For correct geometry alignment, hits produced by the same particle should lie in a straight line.To define the collinearity of three hits, the variable "dev", which represents the deviation of the position of the middle cluster from the straight line connecting the other two clusters, was introduced: where the variables are explained in Fig. 7.For properly calibrated internal geometry, the distribution of the "dev" variables should show a narrow correlation peak centered at zero.The positions resolutions in x and y directions can be then determined from the obtained distributions, which are approximately equal to 2  3 σ x,y , were σ x (σ y ) represent the width of the dev x (dev y ) distribution.The factor 2  3 refers to the equal Vds1 to Vds2 and Vds2 to Vds3 distances in z coordinate (see Fig. 7).
The calibration algorithm uses the MIGRAD function of the MINUIT [40] package.The Variable Metric method was used to minimize the "dev" function to find the optimal alignment parameters.
A detailed description of the applied geometry reconstruction procedure is provided [45].It is seen from the plot presented in Fig. 9, that the obtained position resolution provided by sensors is on the level of the nominal 4.5 µm in both x and y coordinates.

Cluster, track and vertex reconstruction
The first step of data reconstruction is cluster recognition.A particle passing through a sensor may fire more  than one pixel in a given sensor.These pixels should thus not be considered to indicate independent particle hits but rather together constituting a particle hit.Such a composite particle hit is called a "cluster".A computer algorithm, the so-called "clusteriser", identifies such It takes each pixel as a starting point and searches neighboring pixels containing signals in both dimensions.The search is repeated recursively for neighboring fired pixels until no more neighboring fired pixels can be found.The set of fired pixels is used to calculate the center of gravity, taken as the center of the resulting cluster.
The tracks registered in SAVD are slightly curved because of the magnetic field.This curvature is small enough to use a straight line to identify clusters in different stations on the same track.Consequently, a straight line was chosen to describe the tracks: x(z) = Az + x 0 , y(z) = Bz + y 0 . (2) Using this parametrization, a combinatorial track identification procedure based on checking the combinations of all hits from different stations was introduced.If the hits detected on different SAVD stations lie on a straight line according to a χ 2 criterion, the combination is accepted as a reconstructed track.
The track reconstruction procedure was first implemented for the field-off data set.From the distributions of the residuals of hits from the reconstructed and fitted with a straight line tracks, the spatial sensor resolution was determined to be on the level of 5 µm, as was expected.
It turned out that the same straight-line combinatorial method could also be applied to reconstruct the tracks for physics data sets with the field on.However, if the straight-line track model is applied, the hits on the third and fourth stations of SAVD visibly deviate from the fitted straight line.The result of this is a double-peak structure in the distribution of cluster deviations for the x-direction rather than a Gaussian distribution.This effect is caused by the vertical B y component of the magnetic field in the SAVD volume.Therefore in the next steps of the reconstruction, the positions of hits are fitted using a secondorder polynomial function for x and linear for y coordinate: The distribution of the ∆x/∆z ratios for the reconstructed tracks is shown in Fig. 10.The ratios are calculated for track lines reconstructed in the target region, referring to tracks emission angles in the xz-plane.The distribution reflects a clear three-peak structure for each arm.Firstly, the narrow inner-most peak (green peak at small angles) is associated with particles produced far upstream and traveling parallel to the beam for a long distance.Next, the middle structure (gray histogram) corresponds to particles produced upstream of the target.Finally, the outer peak (brown color histogram) is generated by particles produced in the target -these tracks are selected for further analysis.
The primary vertex is 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 squares 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 x prim and y prim coordinates of the primary vertex are afterwards calculated as the average of x and y positions of tracks at z = z prim .
To support the interpretation of components from Fig. 10, the primary vertex reconstruction was performed on the event by event basis separately for tracks within the |∆x/∆z| the interval from 0.01 to 0.025 (gray histogram) and with |∆x/∆z| > 0.025 (brown histogram).By looking at the longitudinal distribution of the primary vertex for these samples of tracks (see Fig. 11) it can be seen that, indeed, the tracks associated with the the most outer peak in Fig. 10 (brown) originates from the target, which is located 47 mm upstream from the first VD station.The primary vertices associated with tracks from the middle peak (gray in Fig. 10) are relatively smoothly distributed upstream of the target in the range from -1200 mm (exit from the beam-line) to -50 mm (near the target).At -190 mm, the distribution has a sharp peak related to interactions in the aluminized Mylar front window of the SAVD box.One can also see that between the window and the target, the frequency of interaction drops by a factor of 5 due to the presence of helium gas in the SAVD vessel.
A target segmented to three 1 mm thick La layers was used for Xe+La data taking.The target structure can be well seen in the z prim distribution shown in Fig. 12.It is seen that the precision of the primary vertex reconstruction allows for determining on which particular target segment the Xe+La collision occurred.
To determine the spatial resolution of the primary vertex reconstruction, the SAVD tracks from an event were split into two non-overlapping sub-events, namely every second track from Jura and Saleve arms, were assigned to sub-event 1, whereas the remaining tracks were assigned to sub-event 2. In this way, one obtains two equivalent track samples.The primary vertex spatial resolutions obtained with sub-event 1 and sub-event 2 are expected to be identical since the opening angle range for both samples is the same.The distributions of differences between x, y and z coordinates of the primary vertices reconstructed using sub-event 1 and sub-event 2 tracks are shown in Fig.    for the Xe+La data.The red lines correspond to Gaussian fits of the distributions.The observed widths of the peaks can be converted to the spatial resolution of the primary vertex, namely σ x = 1.3 µm, σ y = 1 µm and σ z = 15 µm, for x, y and z coordinate, respectively.
After the primary vertex is found, the next step of track reconstruction searches for tracks using the Hough transform (HT) method (for details, see Ref. [41]).It is a global method of track reconstruction where each cluster is processed only once.Thus, the computation time of this method is proportional to the number of all detected hits and is much faster than the combinatorial method, which accesses clusters in the nested loops over clusters grouped according to the station of their detection.However, the HT method requires information about the origin point thus, it is implemented as a second step of the SAVD track reconstruction chain.The HT procedure is based on representing the track as a set of two slope parameters (a x , a y ), which can be used to describe straight track lines according to the following parametrization: where x, y, z are cluster coordinates with respect to the primary vertex position.Then, for each hit its position in coordinate space (x, y, z) are transformed to so-called Hough space of parameters (a x , a y ).Further, hits left by the same particle would have the same track parameters and appear as peaks in the Hough space presented as a 2-dimensional histogram.The algorithm searches for such local peaks which correspond to tracks.However, due to multiple scattering and track curvature, hits that belong to the same track might appear in different bins of the Hough space histogram.Thus, the algorithm performs the clusterisation procedure: combining neighboring bins into one cluster.

Track reconstruction efficiency
To test the reconstruction efficiency of SAVD, the Geant4based simulation study was performed (the effect of the sensor inefficiency was excluded).The efficiency was determined as the ratio between the number of the reconstructed SAVD tracks and the number of the simulated SAVD tracks with three and four hits.Fig. 14 shows the dependence of the efficiency versus track momenta.It is seen that the efficiency is close to 100 % for higher track multiplicity.However, it starts to drop for tracks with momentum < 1 GeV/c.Low momentum tracks have large curvature in the SAVD region (the magnetic field is low but not zero).Thus such tracks can neither be reconstructed within the straight-line model of the combinatorial reconstruction, nor during the Hough Transform stage as the hits belonging to these tracks are transformed into the different Hough space regions.

VD-TPC global geometry calibration
The track multiplicity correlation between tracks reconstructed (all collected events, no trigger selection) in SAVD and TPCs is shown in Fig. 15.As one can see, the multiplicities of SAVD and TPC tracks are well correlated, proving that the tracking procedures described above are correct.It may be observed that for some events, tracks were reconstructed in either VD or TPCs, but not both.These cases are related to low-track multiplicity events selected by the minimum bias trigger.Merging the track fragments measured by SAVD and TPCs requires the SAVD alignment relative to the TPCs.By observing the difference between the positions of reconstructed primary vertices in the SAVD and the TPCs in a given event, the SAVD position was calibrated with an accuracy of 16 µm, 6 µm and 100 µm in the x, y and z coordinate, respectively.

Global tracking
The merging of SAVD and TPC track fragments is done in three steps: (i) Since tracks are not affected by the magnetic field in the y direction, all SAVD tracks are combined with VTPC tracks, and for each SAVD-VTPC track pair, the difference between the tracks slopes in the y coordinate, ∆a y , is calculated.The distribution of ∆a y shows a sharp peak on a large combinatorial background.A ± 5σ cut around this peak is applied to pre-select SAVD and TPC track pairs that potentially match.(ii) For a given track pair, the TPC momentum is assigned to the SAVD track.This allows extrapolating the SAVD track to the VTPC front surface where both are matched in x, y (z is matched by construction as it defines the merging plane) coordinates and the difference of the track positions ∆x and ∆y are calculated.Fig. 16 (left) shows the distribution of ∆y versus y of Saleve side SAVD tracks matched to Jura side tracks of VTPC1.Because the average value of ∆y depends on y, narrow ranges of y of the distribution are projected onto ∆y.The projected distributions (slices) are then fitted with a sum of a second-order polynomial which describes the background related to false-merging cases and a Gaussian peak that accounts for the true ones.An example of a single slice is shown in Fig. 16 (right).
The dependence of the fitted mean ( ∆y ) and variance (σ ∆y ) on y are then fitted with a third-order polynomial function.The results of these fits are shown as red ( ∆y (y)) and blue lines (± σ ∆y (y)) in Fig. 16 (left).A similar procedure was used for ∆x versus z merging.Both ∆y versus y and ∆x versus z distributions were constructed for Jura -Jura, Jura -Saleve, Saleve -Saleve and Saleve -Jura track combinations, separately for VTPC1 and VPTC2.(iii) The values of ∆y , σ ∆y and ∆x , σ ∆x obtained from the fits are used to apply elliptic cuts to select the best merge candidate.
Fig. 17 shows the distribution of the difference between SAVD and TPC momentum components ∆p x and ∆p z calculated at the merging plane for SAVD and TPC track combinations that passed the cut on ∆a y (blue) and with the elliptical 4σ cuts on ∆x and ∆y (red).It can be seen that after the ∆x and ∆y cuts, the distributions are practically free of background.About 75 % of the SAVD tracks are merged with the VTPC tracks.This result corresponds to the performed Geant4-based simulations.The remaining tracks either miss the VTPC acceptance, decay before reaching the VTPC or are not merged due to the SAVD-TPC merging inefficiency, which is about 5 %.
Finally, the global track, which has hits in both SAVD and TPCs, is refitted using a method based on Kalman Filter [43] and used for further analysis.

Secondary vertex resolution
The position resolution of the reconstructed secondary vertices related to open charm mesons decays was determined in the Geant4-based simulations by comparing the simulated and reconstructed positions of the vertices.The differences ∆x, ∆y, ∆z of the coordinates of the reconstructed secondary vertex position and the one defined in the Geant4-based simulations are shown in Fig. 18.The sigma of these distributions determines the primaryvertex resolution to be 20 µm, 11 µm, 170 µm for x, y and z coordinates, respectively.

Invariant mass spectra in Xe+La data
The performance results are based on the 2017 Xe+La data since it is currently the most thoroughly investigated data set.The SAVD tracks matched to TPC tracks are used to search for the D 0 + D 0 signal.The particle identification (PID) information was not used in the analysis.Each SAVD track is paired with another SAVD track and is assumed to be either a kaon or a pion.Thus each pair contributes twice in the combinatorial invariant mass distribution.The combinatorial background is several orders of magnitude higher than the D 0 + D 0 signal due to the low yield of charm particles.Five cuts were applied to reduce the large background.The cut parameters were chosen to maximize the signal-to-noise ratio (SNR) of the reconstructed D 0 + D 0 peak and were determined from the Geant4-based simulations.These cuts are: The d and D parameters are defined as the shortest distance between the primary vertex and the track line of a single track and D 0 candidate, respectively.Note that the last four cuts are based on information delivered by the SAVD.Fig. 19 shows the invariant mass distribution of unlike charge daughter candidates with the applied cuts for 1.86M 0-20% central Xe+La events.One observes a peak emerging at 1.86 GeV/c 2 , consistent with a D 0 + D 0 production.The invariant mass distribution was fitted using an exponential function to describe the background and a Gaussian to describe the D 0 + D 0 signal contribution.Both lines representing signal plus background and background alone are drawn on the plot in red.The indicated errors are statistical only.From the fit, one finds the width of the peak to be 12 ± 3.5 MeV/c 2 , consistent with the value obtained in simulations taking into account instrumental effects.The total yield amounts to 80 ± 28 with a ±3σ integrated SNR of 3.4.The feasibility of D + and D − measurements has been demonstrated so far only by simulations.However, these measurements are foreseen not to be more difficult than these of D 0 and D 0 .4.9 K 0 S and Λ in the Xe+La data The same strategy of background suppression as that described in the previous section can be applied for the reconstruction of K 0 S and Λ particles.Fig. 20 shows the invariant mass distribution in the regions of the K 0 S mass of the unlike sign pairs assigning π mass to both tracks in the pair.The results are drawn for 1.1 × 10 6 collisions of Xe+La at the beam momentum of 150A GeV/c.No event selection was applied.A clear K 0 S peak is seen at 0.498 GeV/c 2 .For the same data Fig.21 presents invariant mass distribution in the regions of the Λ mass for the unlike sign pairs assigning the proton mass to positively charged track and the π − mass to negatively charged track in the pair.As in the case of K 0 S , a clear Λ peak appears at the mass of 1.1156 GeV/c 2 .In both figures, the red line represents a fit with the Gaussian function to account for the signal plus the second-order polynomial to account for the remaining background.The cut parameters were not optimized to maximize the signal significance.In this analysis, we used rather arbitrary cuts to demonstrate the ability of K 0 S and Λ reconstruction.As expected, the Λ peak width is significantly smaller than the width of the K 0 S peak.Our reconstruction over-predicts masses of K 0 S and Λ by 2 MeV/c 2 and 0.7 MeV/c 2 , respectively.Although the shifts are small, they are much larger than the statistical uncertainty and are related to the limited control of the absolute value of the magnetic field.The observed discrepancy can be used to calibrate the absolute strength of the magnetic field.

Summary and outlook
This paper presents the design and construction of a Small Acceptance Vertex Detector developed within NA61/SHINE at the CERN SPS for pioneering measurements of open charm production.Moreover, the SAVD data calibration, event reconstruction and analysis procedure are also presented.

Fig. 1 .
Fig. 1.The layout of the NA61/SHINE experimental setup (top view, not to scale).

Fig. 2 .
Fig. 2. Schematics of reconstruction strategy of D 0 → π + + K − decay channel with the help of the Vertex Detector.

Fig. 3 .
Fig. 3. Photograph of SAVD before closing the detector with the front and exit windows.The detector elements are indicated.For detail, see the text.

Fig. 4 .
Fig. 4. Single SAVD unit composed of two MIMOSA-26AHR sensors, carbon fibre extension plate, flex-print cable and a supporting ladder.The right-bottom part of the figure shows an enlarged view of the sensors, which are also visible in the central part of the ladder.

Fig. 6 .
Fig. 6.Rapidity -transverse momentum spectra of D 0 + D 0 mesons for 20% of the most central Xe+La collisions at 150A GeV/c according to AMPT (top plots) and PHSD (bottom plots) predictions.The left plots show the generated phase space, and the right plots show geometrically accepted phase space.The plots are obtained for 1M generated D 0 + D 0 decaying in the π and K channel.

Fig. 7 .
Fig. 7.The graphical representation of the "dev" variable used for geometry tuning.

Fig. 8 .
Fig. 8.The naming convention of the SAVD sensors.The first number following "Vds" denotes the station, while the second gives the sensor number in a given station.

Fig. 10 .
Fig. 10.∆x/∆z (x-slope) distribution tracks reconstructed in Jura (positive values) and Saleve (negative values) arms.The plot is done for Xe+La at 150A GeV/c data set taken in 2017.Different colors refer to different production components explained in the text.

Fig. 11 .
Fig. 11.Distribution zprim of primary vertices for tracks produced on target (brown color histogram) and production out of the target (dark color histogram).Seed text for more explanation.

Fig. 12 .
Fig. 12. Distribution of of the longitudinal coordinate zprim for the Xe+La data at 150 A GeV/c recorded in 2017.

Fig. 13 .
Fig.13.Distributions of differences between x, y and z coordinates of the primary vertexes reconstructed using sub-event 1 and sub-event 2 tracks (see text for details) for the Xe+La data at 150A GeV/c recorded in 2017.

Fig. 14 .
Fig. 14.The SAVD reconstruction efficiency versus track momenta assuming fully efficient sensors for Xe+La at 150A GeV/c.

Fig. 16 .
Fig. 16.Left: Difference in y coordinate of SAVD and TPC tracks (∆y) versus y at the merging plane.Right: example of the projection of distribution of ∆y versus y onto the ∆y coordinate for −7 mm < y < −2.5 mm (single slice).

Fig. 17 .
Fig.17.Difference of momentum components ∆px and ∆pz calculated at the merging plane for SAVD -TPC track combinations that passed the cut on ∆ay (blue) and after additional elliptical 4σ cuts on ∆x and ∆y (red).

Fig. 18 .
Fig. 18.Distributions of differences between coordinates of the reconstructed secondary vertices and simulated vertices for Xe+La at 150A GeV/c data.Lines correspond to Gaussian fits of the distributions, and the σ parameters resulting from the fits are indicated.
(i) cut on the track transverse momentum, p T > 0.34 GeV/c; (ii) cut on the track impact parameter, d > 37 µm; (iii) cut on the longitudinal distance between the D 0 decay vertex candidate and the primary vertex, V z > 1050 µm; (iv) cut on the impact parameter D of the back extrapolated D 0 candidate momentum vector, D < 18 µm; (v) cut on daughter tracks distance at the closest proximity, DCA < 36 µm.

Fig. 19 .
Fig. 19.Invariant mass distribution of unlike charge sign π, K decay track candidates for Xe+La collisions at 150A GeV/c taken in 2017.The presented results refer to 1.86M 0-20 % central events.

Fig. 20 .
Fig. 20.Invariant mass distribution of unlike charge sign π, π decay track candidates for Xe+La collisions at 150A GeV/c.The plot was done for 1.86M 0-20 % central events.

Fig. 21 .
Fig. 21.Invariant mass distribution of unlike charge sign π, p decay track candidates for Xe+La collisions at 150A GeV/c.The plot was done for 1.86M events of 0-20 % central collisions.