Abstract
The NA61/SHINE experiment at the CERN Super Proton Synchrotron (SPS) studies the onset of deconfinement in hadron matter by a scan of particle production in collisions of nuclei with various sizes at a set of energies covering the SPS energy range. This paper presents results on inclusive double-differential spectra, transverse momentum and rapidity distributions and mean multiplicities of \(\pi ^\pm \), \(K^\pm \), p and \(\bar{p}\) produced in the 20% most central \(^7\)Be+\(^9\)Be collisions at beam momenta of 19A, 30A, 40A, 75A and 150A \({\mathrm{Ge} \mathrm{V}}\!/\!c\). The energy dependence of the \(K^\pm \)/\(\pi ^\pm \) ratios as well as of inverse slope parameters of the \(K^\pm \) transverse mass distributions are close to those found in inelastic p+p reactions. The new results are compared to the world data on p+p and Pb+Pb collisions as well as to predictions of the Epos, Urqmd, Ampt, Phsd and Smash models.
Similar content being viewed by others
Explore related subjects
Discover the latest articles, news and stories from top researchers in related subjects.Avoid common mistakes on your manuscript.
1 Introduction
This paper presents experimental results on inclusive spectra and mean multiplicities of \(\pi ^\pm , K^\pm , p\) and \(\bar{p}\) produced in the 20% most central \(^7\)Be+\(^9\)Be collisions at beam momenta of 19A, 30A, 40A, 75A and 150A \({\mathrm{Ge} \mathrm{V}}\!/\!c\) (\(\sqrt{s_{NN}}\)= 6.1, 7.6, 8.8, 11.9 and 16.8 GeV). These studies form part of the strong interactions programme of NA61/SHINE [1] investigating the properties of the onset of deconfinement and searching for the possible existence of a critical point. This requires a two dimensional scan in collision energy and nuclear mass number of the colliding nuclei. Such a scan allows to explore systematically the phase diagram of strongly interacting matter [1]. An increase of collision energy causes an increase of temperature and a decrease of baryon chemical potential of strongly interacting matter at freeze-out , whereas increasing the nuclear mass number of the colliding nuclei decreases the temperature [2].
The schematic layout of the NA61/SHINE experiment at the CERN SPS [13] showing the components used for the Be+Be energy scan (horizontal cut, not to scale). The beam instrumentation is sketched in the inset (see also Fig. 2). Alignment of the chosen coordinate system as shown in the figure: its origin lies in the middle of VTPC-2, on the beam axis. The nominal beam direction is along the z-axis. The magnetic field bends charged particle trajectories in the x–z (horizontal) plane. The drift direction in the TPCs is along the y (vertical) axis
Pursuing this programme NA61/SHINE recorded data on p+p, Be+Be, Ar+Sc, Xe+La and Pb+Pb collisions. Moreover, further measurements of Pb+Pb interactions are planned with an upgraded detector [3] starting in 2021.
The \(^7\)Be+\(^9\)Be collisions (see Ref. [4] for results on \(\pi ^-\) production) play a special role in the NA61/SHINE scan programme. First, it was predicted within the statistical models [5, 6] that the yield ratio of strange hadrons to pions in these collisions should be close to those in central Pb+Pb collisions and significantly higher than in p+p interactions. Second, the collision system composed of a \(^7\)Be and a \(^9\)Be nucleus has eight protons and eight neutrons, and thus is isospin symmetric. Within the NA61/SHINE scan programme the \(^7\)Be+\(^9\)Be collisions serve as the lowest mass isospin symmetric reference needed to study collisions of medium and large mass nuclei. This is of particular importance when data on proton-proton, neutron-proton and neutron-neutron are not available to construct the nucleon-nucleon reference [7].
The paper is organized as follows: after this introduction the experiment is briefly presented in Sect. 2. The analysis procedure, as well as statistical and systematic uncertainties are discussed in Sect. 3. Section 4 presents experimental results and compares them with measurements of NA61/SHINE in inelastic p+p interactions [8,9,10] and NA49 in Pb+Pb collisions [11, 12]. Section 5 discusses model predictions. A summary in Sect. 6 closes the paper.
The following variables and definitions are used in this paper. The particle rapidity y is calculated in the collision center of mass system (cms), \(y=0.5 \cdot ln{[(E+p_{L})/(E-p_{L})]}\), where E and \(p_{L}\) are the particle energy and longitudinal momentum, respectively. The transverse component of the momentum is denoted as \(p_{T}\) and the transverse mass \(m_{T}\) is defined as \(m_{T} = \sqrt{m^2 + (cp_{T})^2}\) where m is the particle mass in GeV. The momentum in the laboratory frame is denoted \(p_\text {lab}\) and the collision energy per nucleon pair in the center of mass by \(\sqrt{s_{NN}}\).
Results of the measurements correspond to collisions with low energy emitted into the forward beam spectator region. For \(^7\)Be+\(^9\)Be collisions this energy is not tightly correlated with geometric parameters of the interaction such as the collision impact parameter of the collision (see Sect. 3.1). This is caused by the small number of nucleons and the cluster structure of the Be nucleus. Nevertheless, following the convention widely used in the analysis of nucleus-nucleus collisions, the term central is used for events selected by imposing an upper limit on this energy.
2 Experimental setup of NA61/SHINE
2.1 Detector
The NA61/SHINE experiment is a multi-purpose facility designed to measure particle production in nucleus+nucleus, hadron+nucleus and p+p interactions [13]. The detector is situated at the CERN Super Proton Synchrotron (SPS) in the H2 beamline of the North experimental area. A schematic diagram of the setup is shown in Fig. 1.
The main components of the produced particle detection system are four large volume Time Projection Chambers (TPC). Two of them, called Vertex TPCs (VTPC), are located downstream of the target inside superconducting magnets with maximum combined bending power of 9 Tm. The magnetic field was scaled down in proportion to the beam momentum in order to obtain similar phase space acceptance at all energies. The main TPCs (MTPC) and two walls of pixel Time-of-Flight (ToF-L/R) detectors are placed symmetrically to the beamline downstream of the magnets. The fifth small TPC (GAP-TPC) is placed between VTPC1 and VTPC2 directly on the beam line. The TPCs are filled with Ar:CO\(_{2}\) gas mixtures in proportions 90:10 for the VTPCs and the GAP-TPC, and 95:5 for the MTPCs.
The Projectile Spectator Detector (PSD), which measures mainly the energy in the forward region of projectile spectators, is positioned 20.5 m (16.7 m) downstream of the target during measurements at 75A and 150\(A\,{\mathrm{Ge} \mathrm{V}}\!/\!c\) (19A, 30A, 40\(A\,{\mathrm{Ge} \mathrm{V}}\!/\!c\)) centered in the transverse plane on the position of the deflected beam. The PSD is used as a part of the trigger system (see Sect. 2.2) to accept collisions by imposing an upper limit on the energy measured in the 16 central modules and also to select central events in analysis procedure (see Sect. 3.1).
The beamline instrumentation is schematically depicted in Fig. 2. It is designed for obtaining high beam purity with secondary ion beams [14] produced by fragmenting the primary Pb\(^{82+}\) ions extracted from the SPS. A detailed discussion of the properties of selected \(^7\)Be ions can be found in Ref. [4].
A set of scintillation counters as well as beam position detectors (BPDs) upstream of the spectrometer provide timing reference, selection, identification and precise measurement of the position and direction of individual beam particles.
The schematic of the placement of the beam and trigger detectors in high-momentum (top) and low-momentum (bottom) data taking configurations showing beam counters S, veto counters V and beam position and charge detectors BPD, as well as a Cerenkov detector Z. Note, that the PSD calorimeter was almost 4 m closer for low momentum data taking
The target was a plate of \(^9\)Be of 12 mm thickness placed \(\approx \) 80 cm upstream of VTPC1. Mass concentrations of impurities in the target were measured at 0.3%, resulting in an estimated increase of the produced pion multiplicity by less than 0.5% due to the small admixture of heavier elements [15]. No correction was applied for this negligible contamination. Data were taken with target inserted (denoted I, 90%) and target removed (denoted R, 10%).
2.2 Trigger
The schematic of the placement of the beam and trigger detectors can be seen in Fig. 2. The trigger detectors consist of a set of scintillation counters recording the presence of the beam particle (S1, S2), a set of veto scintillation counters with a hole used to reject beam particles passing far from the centre of the beamline (V0, V1), and a Cherenkov charge detector (Z). Beam particles were defined by the coincidence T1 = \(\text {S1} \cdot \text {S2} \cdot \overline{\text {V1}} \cdot \text {Z(Be)}\) and T1 = \(\text {S1} \cdot \overline{\text {V0}} \cdot \overline{\text {V1}} \cdot \overline{\text {V1'}} \cdot \text {Z(Be)}\) for low and high momentum data taking respectively. An interaction trigger detector (S4) was used to check whether the beam particle changed charge after passing through the target. In addition, collisions were selected by requiring an energy signal below a set threshold from the 16 central modules of the PSD. The event trigger condition thus was T2 = T1\(\cdot \overline{\text {S4}} \cdot \overline{\text {PSD}}\) or T2 = T1\(\cdot \overline{\text {PSD}}\) for low and high beam momenta, respectively. The PSD threshold was set to retain from \(\approx \) 70% to \(\approx \) 40% of inelastic collisions at low and high beam momenta, respectively. The statistics of recorded events is summarised in Table 1.
3 Analysis procedure
This section starts with a brief overview of the data analysis procedure and the applied corrections. It also defines to which class of particles the final results correspond. A description of the calibration and the track and vertex reconstruction procedures can be found in Ref. [8].
The analysis procedure consists of the following steps:
-
(i)
application of event and track selection criteria,
-
(ii)
determination of raw spectra of identified charged hadrons using the selected events and tracks,
-
(iii)
evaluation of corrections to the raw spectra based on experimental data and simulations,
-
(iv)
calculation of the corrected spectra and mean multiplicities,
-
(v)
calculation of statistical and systematic uncertainties.
Corrections for the following biases were evaluated:
-
(a)
contribution from off-target interactions,
-
(b)
losses of in-target interactions due to the event selection criteria,
-
(c)
geometrical acceptance,
-
(d)
reconstruction and detector inefficiency,
-
(e)
losses of tracks due to track selection criteria,
-
(f)
contribution of particles other than primary (see below) charged particles produced in Be+Be interactions,
-
(g)
losses of primary charged particles due to their decays and secondary interactions.
Correction (a) was not applied due to insufficient statistics of the target removed data. The contamination of the target inserted data was estimated from the z distribution of fitted vertices to amount to \(\approx \) 0.35%.
Corrections (b)–(g) were estimated by data and simulations. MC events were generated with the Epos1.99 model (version CRMC 1.5.3) [16], passed through detector simulation employing the Geant 3.21 package [17] and then reconstructed by the standard program chain.
The final results refer to particles produced in central Be+Be collisions by strong interaction processes and in electromagnetic decays of produced hadrons. Such hadrons are referred to as primary hadrons. Central collisions refer to events selected by a cut on the total energy emitted into the forward direction as defined by the acceptance maps for the PSD given in Ref. [18].
The analysis was performed in (\(y\), \(p_T\)) bins. The bin size was chosen taking into account the statistical uncertainties and the resolution of the momentum reconstruction [8]. Corrections as well as statistical and systematic uncertainties were calculated for each bin.
3.1 Central collisions
A short description of the procedure for defining central collisions is given below. For more details see Refs. [4, 19].
Final results presented in this paper refer to Be+Be collisions with the 20% lowest values of the forward energy \(E_F\) (central collisions). The quantity \(E_F\) is defined as the total energy in the laboratory system of all particles produced in a Be+Be collision via strong and electromagnetic processes in the forward momentum region defined by the acceptance map in Ref. [18]. Final results on central collisions, derived using this procedure, allow a precise comparison with predictions of models without any additional information about the NA61/SHINE setup and used magnetic field.
For analysis of the data the event selection was based on the \(\approx \) 20% of collisions with the lowest value of the energy \(E_{PSD}\) measured by a subset of PSD modules (see Fig. 3) in order to optimize the sensitivity to projectile spectators. The forward momentum acceptance in the definition of \(E_F\) corresponds closely to the acceptance of this subset of PSD modules.
Online event selection by the hardware trigger (T2) used a threshold on the electronic sum of energies over the 16 central modules of the PSD set to accept \(\approx \) 40% of the inelastic interactions. The minimum-bias distribution was obtained using the data from the beam trigger T1 with offline selection of events by requiring an event vertex in the target region and a cut on the ionisation energy detected in the GTPC to exclude Be beams. The spectrum of \(E_{PSD}\) was calculated for the subset of modules and a properly normalized spectrum for target removed events was subtracted. Events for further analysis were selected by applying a cut in this distribution for the 20% of events with the smallest values of \(E_{PSD}\). More details on the PSD modules selection optimisation and \(E_{PSD}\) determination can be found in [19].
The forward energy \(E_F\) cannot be measured directly. However, both \(E_F\) and \(E_{PSD}\) can be obtained from simulations using the Epos1.99 (version CRMC 1.5.3) [16] model. A global factor \(c_{cent}\) (listed in Table 2) was then calculated as the ratio of mean negatively charged pion multiplicities obtained with the two selection procedures for 20% of all inelastic collisions. A possible dependence of the scaling factor on rapidity and transverse momentum was neglected. The resulting factors \(c_{cent}\) range from 1.00 and 1.04 which is only a small correction compared to the systematic uncertainties of the measured particle multiplicities. The correction was therefore not applied, but instead included in the systematic uncertainty.
Finally, the average number of wounded nucleons \(\langle W \rangle \) and the average collision impact parameter \(\langle b \rangle \) were calculated within the Wounded Nucleon Model [20] implemented in Epos for events with the 20% smallest values of \(E_F\). Results are listed in Table 2. Example distributions for the top beam momentum are shown in Fig. 4. As the Be nucleus consists of few nucleons these distributions are quite broad. For comparison \(\langle W \rangle \) and \(\langle b \rangle \) were also calculated from the GLISSANDO model [21] which uses a different Glauber model calculation. The results, also listed in Table 2, differ by about 10% for \(\langle W \rangle \). This discrepancy was included in the systematic error estimate (see Sect. 3.5.2) .
Examples of the distribution of the number of wounded nucleons W (left) and collision impact parameter b (right) for events with the 20% smallest forward energies \(E_F\) (central) at beam momentum of 150\(A\,{\mathrm{Ge} \mathrm{V}}\!/\!c\) simulated with the Epos model using the acceptance map provided in Ref. [18]
3.2 Event and track selection
3.2.1 Event selection
For further analysis Be+Be events were selected using the following criteria:
-
(i)
four units of charge measured in S1, S2, and Z counters as well as BPD3 (this requirement also rejects most interactions upstream of the Be target),
-
(ii)
no off-time beam particle detected within a time window of ±4.5\(~\upmu \)s around the trigger particle,
-
(iii)
no other event trigger detected within a time window of ±25\(~\upmu \)s around the trigger particle,
-
(iv)
beam particle detected in at least two planes out of four of BPD-1 and BPD-2 and in both planes of BPD-3,
-
(v)
a well reconstructed interaction vertex with z position (fitted using the beam trajectory and TPC tracks) not farther away than 15 cm from the center of the Be target (the cut removes less than 0.4% of T2 trigger (\(E_{PSD}\)) selected interactions),
-
(vi)
an upper cut on the measured energy \(E_{PSD}\) which selects 20% of all inelastic collisions.
The event statistics after applying the selection criteria is summarized in Table 1.
3.2.2 Track selection
In order to select tracks of primary charged hadrons and to reduce the contamination by particles from secondary interactions, weak decays and off-time interactions, the following track selection criteria were applied:
-
(i)
track momentum fit including the interaction vertex should have converged,
-
(ii)
fitted x component of particle rigidity \(q \cdot p_\text {lab} \) is positive. This selection minimizes the angle between the track trajectory and the TPC pad direction for the chosen magnetic field direction, reducing uncertainties of the reconstructed cluster position, energy deposition and track parameters,
-
(iii)
total number of reconstructed points on the track should be greater than 30,
-
(iv)
sum of the number of reconstructed points in VTPC-1 and VTPC-2 should be greater than 15 or greater than 4 in the GTPC,
-
(v)
the distance between the track extrapolated to the interaction plane and the (track impact parameter) should be smaller than 4 cm in the horizontal (bending) plane and 2 cm in the vertical (drift) plane.
3.3 Identification techniques
Charged particle identification in the NA61/SHINE experiment is based on the ionization energy loss, \(\text {d}E/\text {d}x\) , in the gas of the TPCs and the time of flight, tof, obtained from the ToF-L and ToF-R walls. In the region of the relativistic rise of the ionization at large momenta the measurement of \(\text {d}E/\text {d}x\) alone allows identification. At lower momenta the \(\text {d}E/\text {d}x\) bands for different particle species overlap and additional measurement of tof is required to remove the ambiguity. These two methods allow to cover most of the phase space in rapidity and transverse momentum which is of interest for the strong interaction programme of NA61/SHINE. The acceptance of the two methods is shown in Figs. 5 and 6 for the 20% most central Be+Be interactions at 30A and 150\(A\,{\mathrm{Ge} \mathrm{V}}\!/\!c\), respectively. At low beam energies the tof-\(\text {d}E/\text {d}x\) method extends the identification acceptance, while at top SPS energy it overlaps with the \(\text {d}E/\text {d}x\) method (for more details see Ref. [22]).
3.3.1 Identification based on energy loss measurement (\(\text {d}E/\text {d}x\) )
Time projection chambers provide measurements of energy loss \(\text {d}E/\text {d}x\) of charged particles in the chamber gas along their trajectories. Simultaneous measurements of \(\text {d}E/\text {d}x\) and \(p_\text {lab}\) allow to extract information on particle mass. The mass assignment follows the procedure which was developed for the analysis of p+p reactions as described in Ref. [9]. Values of \(\text {d}E/\text {d}x\) are calculated as the truncated mean (smallest 50%) of ionisation energy loss measurements along the track trajectory. As an example, \(\text {d}E/\text {d}x\) measured in the 20% most central Be+Be interactions at 75\(A\,{\mathrm{Ge} \mathrm{V}}\!/\!c\) is presented in Fig. 7, for positively and negatively charged particles, as a function of \(q \cdot p_\text {lab} \).
Distribution of charged particles in the \(\text {d}E/\text {d}x\) – \(q\cdot \) \(p_\text {lab}\) plane. The energy loss in the TPCs for different charged particles for events and tracks selected for the analysis of Be+Be collisions at 75\(A\,{\mathrm{Ge} \mathrm{V}}\!/\!c\) (the target inserted configuration). Expectations for the dependence of the mean \(\text {d}E/\text {d}x\) on \(p_\text {lab}\) for the considered particle types are shown by the curves calculated based on the Bethe–Bloch function
The contributions of \(e^{+}\), \(e^{-}\), \(\pi ^{+}\), \(\pi ^{-}\), \(K^{+}\), \(K^{-}\), p and \(\bar{p}\) are obtained by fitting the \(\text {d}E/\text {d}x\) distributions separately for positively and negatively charged particles in bins of \(p_\text {lab}\) and \(p_{T}\) with a sum of four functions [23, 24] each corresponding to the expected \(\text {d}E/\text {d}x\) distribution for the corresponding particle type. The small contribution of light (anti-)nuclei was neglected.
In order to ensure similar particle multiplicities in each bin, 20 logarithmic bins are chosen in \(p_\text {lab}\) in the range \(1-100\) \({\mathrm{Ge} \mathrm{V}}\!/\!c\) to cover the full detector acceptance. Furthermore, the data are binned in 20 equal \(p_{T}\) intervals in the range 0-2 \({\mathrm{Ge} \mathrm{V}}\!/\!c\).
The \(\text {d}E/\text {d}x\) distributions for negatively (top, left) and positively (top, right) charged particles in the bin 12.6 \(\le \) \(p_\text {lab}\) \(\le \) 15.8 \({\mathrm{Ge} \mathrm{V}}\!/\!c\) and 0.2 \(\le \) \(p_{T}\) \(\le \) 0.3 \({\mathrm{Ge} \mathrm{V}}\!/\!c\) produced in PSD selected Be+Be collisions at 75\(A\,{\mathrm{Ge} \mathrm{V}}\!/\!c\). The fit by a sum of contributions from different particle types is shown by solid lines. The corresponding residuals (the difference between the data and fit divided by the statistical uncertainty of the data) is shown in the bottom plots
Mass squared derived using time-of-flight measured by ToF-R (right) and ToF-L (left) versus laboratory momentum for particles produced in PSD selected Be+Be collisions at 75\(A\,{\mathrm{Ge} \mathrm{V}}\!/\!c\) (target inserted configuration). The lines show the expected mass squared values for different hadrons species
The distribution of \(\text {d}E/\text {d}x\) for tracks of a given particle type i is parameterised as the sum of Gaussians with widths \(\sigma _{i,l}\) depending on the particle type i and the number of points l measured in the TPCs. Simplifying the notation in the fit formulae, the peak position of the \(\text {d}E/\text {d}x\) distribution for particle type i is denoted as \(x_{i}\). The contribution of a reconstructed particle track to the fit function reads:
where x is the \(\text {d}E/\text {d}x\) of the particle, \(n_{l}\) is the number of tracks with number of points l in the sample and \(A_{i}\) is the amplitude of the contribution of particles of type i. The second sum is the weighted average of the line-shapes from the different numbers of measured points (proportional to track-length) in the sample. The quantity \(\sigma _{l}\) is written as:
where the width parameter \(\sigma _{0}\) is assumed to be common for all particle types and bins. A \(1/\sqrt{l}\) dependence on number of points is assumed. The Gaussian peaks could in principle be asymmetric if the tail of the Landau distribution persists to some extent even after truncation (for detail see [25]). However, no significant effect was found.
The fit function has nine parameters (4 amplitudes, 4 peak positions and width) which are difficult to fit in each bin independently. Therefore the following constraints on the fitting parameters were adopted:
-
(i)
positions of electrons, kaons and protons relative to pions were assumed to be \(p_{T}\)-independent,
-
(ii)
the fitted amplitudes were required to be greater than or equal to 0,
-
(iii)
the electron amplitude was set to zero for total momentum \(p_\text {lab}\) above 23.4 \({\mathrm{Ge} \mathrm{V}}\!/\!c\) (i.e. starting from the 13\(^{th}\) bin), as the electron contribution vanishes at high \(p_\text {lab}\),
-
(iv)
if possible, the relative position of the positively charged kaon peak was taken to be the same as that of negatively charged kaons determined from the negatively charged particles in the bin of the same \(p_\text {lab}\) and \(p_{T}\). This procedure helps to overcome the problem of the large overlap between \(K^+\) and protons in the \(\text {d}E/\text {d}x\) distributions.
The constraints reduce the number of independently fitted parameters in each bin from 9 to 6, i.e. the amplitudes of the four particle types, the pion peak position and the width parameter \(\sigma _{0}\).
Examples of fits are shown in Fig. 8 and the values of the fitted peak positions \(x_i\) are plotted in Fig. 9 versus momentum for different particle types i in selected Be+Be interactions at 150 \(A\,{\mathrm{Ge} \mathrm{V}}\!/\!c\). As expected, the values of \(x_i\) increase with \(p_\text {lab}\) but do not depend on \(p_{T}\).
In order to ensure good fit quality, only bins with total number of tracks grater than 1000 (500 in the 19\(A\,{\mathrm{Ge} \mathrm{V}}\!/\!c\) sample) are used for further analysis. The Bethe–Bloch curves for different particle types cross each other at low values of the total momentum. Thus, the technique is not sufficient for particle identification at low \(p_\text {lab}\) and bins with \(p_\text {lab}\) \(< 3.98\) \({\mathrm{Ge} \mathrm{V}}\!/\!c\) (bins 1–5) are excluded from the analysis based solely on \(\text {d}E/\text {d}x\) measurements.
3.3.2 Identification based on time of flight and energy loss measurements (tof-\(\text {d}E/\text {d}x\) )
Identification of \(\pi ^{+}\), \(\pi ^{-}\), \(K^{+}\), \(K^{-}\), p and \(\bar{p}\) at low momenta (from 2-8 \({\mathrm{Ge} \mathrm{V}}\!/\!c\)) is possible when measurement of \(\text {d}E/\text {d}x\) is combined with time of flight information tof. Timing signals from the constant-fraction discriminators and signal amplitude information are recorded for each tile of the ToF-L/R walls. Only hits which satisfy quality criteria (see Ref. [26] for detail) are selected for the analysis. The coordinates of the track intersection with the front face are used to match the track to tiles with valid tof hits. The position of the extrapolation point on the scintillator tile is used to correct the measured value of tof for the propagation time of the light signal. The distribution of the difference between the corrected tof measurement and the value calculated from the extrapolated track trajectory length with the assumed mass hypothesis can be well described by a Gaussian with standard deviation of 80 ps for ToF-R and 100 ps for ToF-L. These values represent the tof resolution including all detector effects.
Momentum phase space is subdivided into bins of 1 \({\mathrm{Ge} \mathrm{V}}\!/\!c\) in \(p_\text {lab}\) and 0.1 \({\mathrm{Ge} \mathrm{V}}\!/\!c\) in \(p_{T}\). Only bins with more than 200 (500 in the 150\(A\,{\mathrm{Ge} \mathrm{V}}\!/\!c\) sample) entries were used for extracting yields with the tof-\(\text {d}E/\text {d}x\) method.
The square of the particle mass \(m^2\) is obtained from tof, the momentum \(p\) and the fitted trajectory length l:
For illustration distributions of \(m^2\) versus \(p_\text {lab}\) are plotted in Fig. 10 for positively (left) and negatively (right) charged hadrons produced in Be+Be interactions at 75\(A\,{\mathrm{Ge} \mathrm{V}}\!/\!c\). Bands which correspond to different particle types are visible. Separation between pions and kaons is possible up to momenta of about 5 \({\mathrm{Ge} \mathrm{V}}\!/\!c\), between pions and protons up to about 8 \({\mathrm{Ge} \mathrm{V}}\!/\!c\).
Example distributions of particles in the \(m^2\)–\(\text {d}E/\text {d}x\) plane for the selected Be+Be interactions at 40\(A\,{\mathrm{Ge} \mathrm{V}}\!/\!c\) are presented in Fig. 11. Simultaneous \(\text {d}E/\text {d}x\) and tof measurements lead to improved separation between different hadron types. In this case a simple Gaussian parametrization of the \(\text {d}E/\text {d}x\) distribution for a given hadron type can be used.
Particle number distribution in the \(m^2\)-\(\text {d}E/\text {d}x\) plane for negatively (left) and positively (right) charged particles with momenta 3 < \(p_\text {lab}\) < 4 \({\mathrm{Ge} \mathrm{V}}\!/\!c\) and 0.3 < \(p_{T}\) < 0.4 \({\mathrm{Ge} \mathrm{V}}\!/\!c\)) for PSD selected Be+Be collisions at 40\(A\,{\mathrm{Ge} \mathrm{V}}\!/\!c\)
The tof-\(\text {d}E/\text {d}x\) identification method proceeds by fitting the 2-dimensional distribution of particles in the \(\text {d}E/\text {d}x\) -\(m^2\) plane. Fits were performed in the momentum range from 1–8 \({\mathrm{Ge} \mathrm{V}}\!/\!c\) and transverse momentum range 0–1 \({\mathrm{Ge} \mathrm{V}}\!/\!c\). For positively charged particles the fit function included contributions of p, \(K^+\), \(\pi ^+\) and \(e^+\), and for negatively charged particles the corresponding anti-particles were considered. The fit function for a given particle type was assumed to be a product of a Gauss function in \(\text {d}E/\text {d}x\) and a sum of two Gauss functions in \(m^2\) (in order to describe the broadening of the \(m^2\) distributions with momentum). In order to simplify the notation in the fit formulae, the peak positions of the \(\text {d}E/\text {d}x\) and \(m^2\) Gaussians for particle type j are denoted as \(x_{j}\) and \(y_{j}\), respectively. The contribution of a reconstructed particle track to the fit function reads:
where \(N_{j}\) and f are amplitude parameters, \(x_j\), \(\sigma _{x}\) are means and width of the \(\text {d}E/\text {d}x\) Gaussians and \(y_j\), \(\sigma _{y1}\), \(\sigma _{y2}\) are mean and width of the \(m^2\) Gaussians, respectively. The total number of parameters in Eq. (4) is 16. Imposing the constraint of normalisation to the total number of tracks N in the kinematic bin
the number of parameters is reduced to 15. Two additional assumptions were adopted:
-
(i)
the fitted amplitudes were required to be greater than or equal to 0,
-
(ii)
\(\sigma _{y1} < \sigma _{y2} \) and \(f > 0.7 \), the ”core” distribution dominates the \(m^2\) fit.
An example of the tof-\(\text {d}E/\text {d}x\) fit obtained in a single phase-space bin for positively charged particles in PSD selected Be+Be collisions at 40\(A\,{\mathrm{Ge} \mathrm{V}}\!/\!c\) is shown in Fig. 12.
Example of the tof-\(\text {d}E/\text {d}x\) fit (Eq. 4) obtained in a single phase-space bin (3 < \(p_\text {lab}\) < 4 \({\mathrm{Ge} \mathrm{V}}\!/\!c\) and 0.3 < \(p_{T}\) < 0.4 \({\mathrm{Ge} \mathrm{V}}\!/\!c\)) for positively charged particles in PSD selected Be+Be collisions at 40\(A\,{\mathrm{Ge} \mathrm{V}}\!/\!c\). Lines show projections of the fits for pions (red), kaons (green), protons (blue) and electrons (magenta). Bottom right panel shows the fit residuals
The tof-\(\text {d}E/\text {d}x\) method allows to fit the kaon yield close to mid-rapidity. This is not possible using the \(\text {d}E/\text {d}x\) method. Moreover, the kinematic domain in which pion and proton yields can be fitted is enlarged. The results from both methods partly overlap at the highest beam momenta. In these regions the results from the \(\text {d}E/\text {d}x\) method were selected since they have smaller uncertainties.
3.3.3 Probability method
The fit results allow to calculate the probability \(P_i\) that a measured particle is of a given type i = \(\pi , K, p, e\). For the \(\text {d}E/\text {d}x\) fits (see Eq. 1) one gets:
where \(\rho _{i}\) is the value of the fitted function in a given (\(p_\text {lab}\), \(p_{T}\)) bin calculated for \(\text {d}E/\text {d}x\) of the particle.
Similarly the tof-\(\text {d}E/\text {d}x\) fits (see Eq. 4) give the particle type probability as
For illustration, particle type probability distributions for positively and negatively charged particles produced in PSD selected Be+Be collisions at 19A and 150\(A\,{\mathrm{Ge} \mathrm{V}}\!/\!c\) are presented in Fig. 13 for the \(\text {d}E/\text {d}x\) fits and in Fig. 14 for the tof-\(\text {d}E/\text {d}x\) fits. In the case of perfect particle type identification the probability distributions in Figs. 13 and 14 will show entries at 0 and 1 only. In the case of incomplete particle identification (overlapping \(\text {d}E/\text {d}x\) or tof-\(\text {d}E/\text {d}x\) distributions) values between these extremes will also be populated.
The probability method allows to transform fit results performed in (\(p_\text {lab}\), \(p_{T}\)) bins to results in (\(y\), \(p_{T}\)) bins. Hence, for the probability method the mean number of identified particles in a given kinematical bin (e.g. (\(p_\text {lab}\), \(p_{T}\))) is given by [27]:
for the \(\text {d}E/\text {d}x\) identification method and:
for the tof-\(\text {d}E/\text {d}x\) procedure, where \(P_i\) is the probability of particle type i given by Eqs. (6) and (7), j the summation index running over all entries \(N_{trk}\) in the bin, \(N_{ev}\) is the number of selected events.
Statistical uncertainties of multiplicities calculated with probability method were derived from the variance of the distribution of \(P_i\) in the (\(p_\text {lab}\),\(p_{T}\)) bin:
3.4 Corrections and uncertainties
In order to estimate the true number of each type of identified particle produced in Be+Be interactions a set of corrections was applied to the extracted raw results. These were obtained from a simulation of the NA61/SHINE detector followed by event reconstruction using the standard reconstruction chain. Only inelastic Be+Be interactions were simulated in the target material. The Epos1.99 model (version CRMC 1.5.3) [16] was selected to generate primary inelastic interactions as it best describes the NA61/SHINE measurements. A Geant3 based program chain was used to track particles through the spectrometer, generate decays and secondary interactions and simulate the detector response (for details see Ref. [8]). Simulated events were then processed using the standard NA61/SHINE reconstruction chain and reconstructed tracks were matched to the simulated particles based on the cluster positions. The event selection was based on a dedicated simulation of the energy recorded by the PSD (see Sect. 3.1). Corrections depend on the particle identification technique (i.e. \(\text {d}E/\text {d}x\) or tof-\(\text {d}E/\text {d}x\) ). Hadrons which were not produced in the primary interaction can amount to a significant fraction of the selected tracks. Thus a special effort was undertaken to evaluate and subtract this contribution. The correction factors were calculated in the same bins of \(y\) and \(p_{T}\) as the particle spectra. Bins with a correction factor lower than 1.5 and higher than 4 are caused by low acceptance or high contamination of non primary particles and were rejected.
3.4.1 Corrections for the \(\text {d}E/\text {d}x\) method
The correction factor \(c_i^{dEdx}\)(\(y\),\(p_{T}\)) for biasing effects listed in Sect. 3.2.1 items (b)–(g) was calculated as:
where \(n[i]^{MC}_{gen}\) is the number primary particles in the (\(y\),\(p_{T}\)) bin for simulated events and \(n[i]^{MC}_{sel}\) the number of reconstructed tracks passing all event and track selection cuts. The uncertainty of \(c_i^{dEdx}\)(\(y\), \(p_{T}\)) was calculated assuming that the denominator \(n[i]^{MC}_{sel}(y, p_{T})\) is a subset of the numerator \(n[i]^{MC}_{gen}(y, p_{T})\) and thus has a binomial distribution. The uncertainty of \(c_i^{dEdx}\)(\(y\), \(p_{T}\)) is thus given by:
The mean multiplicity for particle type i in the PSD selected events in a (\(y\),\(p_{T}\)) is calculated as:
3.4.2 Corrections for the tof-\(\text {d}E/\text {d}x\) method
The corrections for the tof-\(\text {d}E/\text {d}x\) method were calculated based on simulation and data. The ToF tile efficiency \(\epsilon _{pixel}(p,p_{T})\) was calculated from the data in the following way. Each reconstructed track was extrapolated to the ToF walls and if it crossed one of the ToF tiles it was classified as a hit and summed in \(n[i]_{tof}\) . Moreover, it was accepted as a valid ToF hit, summed in \(n[i]_{hit}\), if the signal satisfied quality criteria given in Ref. [26]. Finally, the ToF tile efficiency \(\epsilon _{pixel}(p,p_{T})\) was calculated as the ratio of the number of tracks \(n[i]_{tof}\) crossing a particular tile to the number of tracks \(n[i]_{hit}\) with valid ToF hits. The corresponding efficiency factor \(\epsilon _{pixel}(p,p_{T})\) is given by:
This ToF pixel efficiency factor was used in the MC simulation by weighting each reconstructed MC track passing all event and track selection cuts by the efficiency factor of the corresponding tile. Then, the number of selected MC tracks \(n[i]^{MC}_{sel}\) becomes a sum of weights:
Only hits in working tiles with efficiency higher than 50% were taken into account in the identification and correction procedures. This results in the following correction factor for biasing effects listed in Sect. 3.2.1 items (b)–(g) as well as ToF efficiency, tof-\(\text {d}E/\text {d}x\) method acceptance, secondary interactions and contribution of particles other than primary:
where \(n[i]^{MC}_{gen}\) is the number primary particles in the (\(y\),\(p_{T}\)) bin for simulated events and \(n[i]^{MC}_{sel}\) the number of reconstructed tracks passing all event and track selection cuts weighted by the tile efficiency factor \(\epsilon _{pixel}(p,p_{T})\) given by Eg. 15.
The uncertainty of \(c_i^{dEdx,m^2}(y,p_{T})\) was calculated assuming that the denominator \(n[i]^{MC}_{sel}(y,pT)\) is a subset of the nominator \(n[i]^{MC}_{gen}(y,pT)\) and thus has a binomial distribution. The uncertainty of \(c_i^{dEdx, m^2}\)(\(y\), \(p_{T}\)) was calculated as follows:
The mean multiplicity for particle type i in the PSD selected events in a (\(y\),\(p_{T}\)) bin for the tof-\(\text {d}E/\text {d}x\) method is defined by:
3.5 Corrected spectra
Final spectra of different types of hadrons produced in Be+Be interactions are calculated as:
where \(\Delta y \) and \(\Delta p_{T} \) are the bin sizes, \(N_{ev}\) is the total number of accepted events, \(n[i]^{raw}\) represents the mean multiplicity for particle type i in the \(y\),\(p_{T}\) bin obtained as \(n[i]^{dEdx}\) and \(n[i]^{dEdx,m^2}\) for the \(\text {d}E/\text {d}x\) and tof-\(\text {d}E/\text {d}x\) identification method, respectively and finally the correction factor \(c_i\) stands for \(c_i^{dEdx}(y,p_{T})\) (Eq. 11) or \(c_i^{dEdx,m^2}(y,p_{T})\) (Eq. 16) for the \(\text {d}E/\text {d}x\) and tof-\(\text {d}E/\text {d}x\) identification method, respectively.
Resulting two-dimensional distributions \(\frac{d^{2}n}{dydp_T}\) of \(\pi ^{-}, \pi ^{+}, K^{-}, K^{+}, p\) and \(\bar{p}\) produced in the 20% most central Be+Be collisions at different SPS energies are presented in Fig. 15.
3.5.1 Statistical uncertainties
Statistical uncertainties of yields (Eq. 20) were calculated as:
where \(\sigma _{n[i]^{raw}}\) (Eq. 10) is the uncertainty of the uncorrected particle multiplicity \(n[i]^{raw}\) (Eqs. 8 and 9), \(\sigma _{c_i(y, p_T)}\) (Eqs. 12 and 17) denotes the uncertainty of the correction factors \(c^{dEdx}_i(y,pT)\) (Eq. 11) or \(c^{dEdx, m^2}_i(y,pT)\) (Eq. 16) for the \(\text {d}E/\text {d}x\) or tof-\(\text {d}E/\text {d}x\) identification method, i is the particle type and \(N_{ev}\) the total number of accepted events.
3.5.2 Systematic uncertainties
The contributions to the systematic uncertainty for the \(\text {d}E/\text {d}x\) and tof-\(\text {d}E/\text {d}x\) methods in selected \(p_{T}\) intervals at beam momenta of 19A (30A) and 150 \(A\,{\mathrm{Ge} \mathrm{V}}\!/\!c\) are presented in Figs. 16, 17, 18, 19. Assuming they are uncorrelated, the total systematic uncertainty, also shown in the plots, was calculated as the square root of the sum of squares of the described components.
Contributions to the systematic uncertainty of particle spectra obtained from the \(\text {d}E/\text {d}x\) method at 19\(A\,{\mathrm{Ge} \mathrm{V}}\!/\!c\) as a function of rapidity for the transverse momentum interval between 0.2–0.3 \({\mathrm{Ge} \mathrm{V}}\!/\!c\). \(\sigma _{\text {(i)}}\) (red lines) refers to event selection \(\sigma _{\text {(ii)}}\) (magenta lines) to the track selection procedure, \(\sigma _{\text {(iii)}}\) (blue lines) to the identification technique and \(\sigma _{\text {(iv)}}\) (orange lines) to the contamination by feeddown from weak decays of strange particles. Black lines, \(\sigma _{\text {(sys)}}\), show the total systematic uncertainty calculated as the square root of the sum of squares of the components
Contributions to the systematic uncertainty of particle spectra obtained from the \(\text {d}E/\text {d}x\) method at 150\(A\,{\mathrm{Ge} \mathrm{V}}\!/\!c\) as a function of rapidity for the transverse momentum interval between 0.2–0.3 \({\mathrm{Ge} \mathrm{V}}\!/\!c\). \(\sigma _{\text {(i)}}\) (red lines) refers to event selection \(\sigma _{\text {(ii)}}\) (magenta lines) to the track selection procedure, \(\sigma _{\text {(iii)}}\) (blue lines) to the identification technique and \(\sigma _{\text {(iv)}}\) (orange lines) to the contamination by feeddown from weak decays of strange particles. Black lines, \(\sigma _{\text {(sys)}}\), show the total systematic uncertainty calculated as the square root of the sum of squares of the components
Contributions to the systematic uncertainty of particle spectra obtained from the tof-\(\text {d}E/\text {d}x\) method at 30\(A\,{\mathrm{Ge} \mathrm{V}}\!/\!c\) as a function of \(p_{T}\) for the rapidity interval − 0.2 to 0.0 \({\mathrm{Ge} \mathrm{V}}\!/\!c\) (\(K^+\) and \(K^-\)) and − 0.4 to − 0.2 \({\mathrm{Ge} \mathrm{V}}\!/\!c\) (protons). \(\sigma _{\text {(i)}}\) (red lines) refers to event selection \(\sigma _{\text {(ii)}}\) (green lines) to the track selection procedure, \(\sigma _{\text {(iii)}}\) (blue lines) to the identification technique and \(\sigma _{\text {(iv)}}\) (orange lines) to the contamination by feeddown from weak decays of strange particles. Black lines, \(\sigma _{\text {(sys)}}\), show the total systematic uncertainty calculated as the square root of the sum of squares of the components
Contributions to the systematic uncertainty of particle spectra obtained from the tof-\(\text {d}E/\text {d}x\) method at 150\(A\,{\mathrm{Ge} \mathrm{V}}\!/\!c\) as a function of \(p_{T}\) for the rapidity interval 0.0 to 0.2 \({\mathrm{Ge} \mathrm{V}}\!/\!c\) (\(K^+\) and \(K^-\)) and − 0.4 to − 0.2 \({\mathrm{Ge} \mathrm{V}}\!/\!c\) (protons). \(\sigma _{\text {(i)}}\) (red lines) refers to event selection \(\sigma _{\text {(ii)}}\) (green lines) to the track selection procedure, \(\sigma _{\text {(iii)}}\) (blue lines) to the identification technique and \(\sigma _{\text {(iv)}}\) (orange lines) to the contamination by feeddown from weak decays of strange particles. Black lines, \(\sigma _{\text {(sys)}}\), show the total systematic uncertainty calculated as the square root of the sum of squares of the components
The considered contributions to systematic uncertainties are listed below:
-
(i)
Event selection:
Systematic uncertainty of final particle yields due to the slightly different procedures of event selection used for data (\(E_{PSD}\)) and simulated events (\(E_F\)) (see Sect. 3.1).
Systematic uncertainty related to the rejection of events with additional tracks from off-time particles was estimated by changing the width of the time window in which no second beam particle is allowed by \(\pm ~1~\upmu s\) with respect to the nominal value of ±4.5\(~\upmu \)s. The maximal difference of the results was assigned as the systematic uncertainty of the selection. This contribution does not affect the results of the \(tof-\) \(\text {d}E/\text {d}x\) identification method.
Systematic uncertainty due to the choice of selection window for the z-position of the fitted vertex was estimated by varying the selection criteria for the data and the Epos1.99 model in the range of ±25 \(\text{ cm }\) around the nominal position of the target.
-
(ii)
Track selection:
Systematic uncertainty from the value of the cut on the track impact parameter at the primary vertex was estimated by varying the cut for the data and the Epos1.99 model by ±50% around the nominal value.
Systematic uncertainty originating from the requirement on the number of measured points in the magnetic field (minimum number of points in VTPCs and GTPC) was estimated by changing the nominal requirement on the number of measured points by \(\pm ~5\) (33% of the standard selection) and \(\pm ~10\) (66% of the standard selection) for the \(\text {d}E/\text {d}x\) and \(tof-\) \(\text {d}E/\text {d}x\) identification methods, respectively.
-
(iii)
Particle identification:
Uncertainties of the \(\text {d}E/\text {d}x\) identification method were studied and estimated by a 10% variation of the parameter constraints for Eq. (1) fitted to the \(\text {d}E/\text {d}x\) spectra.
In case of tof-\(\text {d}E/\text {d}x\) identification, systematic uncertainties were estimated by shifting the mean (\(x_j\) and \(y_j\)) of the two-dimensional Gaussians (Eq. 4) fitted to the \(m^2-dE/dx\) distributions by \(\pm 1\%\).
An additional systematic uncertainty arises for the \(tof-\) \(\text {d}E/\text {d}x\) method from the quality requirements on the signals registered in the ToF pixels. This systematic uncertainty was estimated by changing the nominal thresholds by ±10%.
-
(iv)
Feeddown correction:
The determination of the feeddown correction is based on the Epos1.99 model which reasonably describes the available cross section data for strange particles in p+p collisions (see e.g. for \(K^{+}\), \(K^{-}\) Ref. [9] and for \(\Lambda \) at 158 \({\mathrm{Ge} \mathrm{V}}\!/\!c\) Ref. [28]). Systematic uncertainty comes from the lack of precise knowledge of the production cross section in Be+Be collisions of \(K^{+}\), \(K^{-}\), \(\Lambda \), \(\Sigma ^{+}\), \(\Sigma ^{-}\), \(K^{0}_{s}\) and \(\bar{\Lambda }\) in case of pions, and in addition of \(\Sigma ^{+}\) in case of protons, and \(\bar{\Lambda }\) in case of \(\bar{p}\).
4 Results
Two dimensional distributions \(\frac{d^{2}n}{dydp_T}\) of \(\pi ^{-}\), \(\pi ^{+}\), \(K^{-}\), \(K^{+}\), p and \(\bar{p}\) produced in the 20% most central Be+Be interactions at different SPS energies are presented in Fig. 15. Where available, results from the \(\text {d}E/\text {d}x\) method were used because of their smaller statistical uncertainties. Results from the tof-\(\text {d}E/\text {d}x\) method were taken to extend the momentum space coverage. Reflection symmetry around \(y\) = 0 was used for \(\frac{d^{2}n}{dydp_T}\) near mid-rapidity (see discussion in Sect. 4.2). Empty bins in momentum space (mostly for lower energies) are caused by insufficient acceptance for the identification methods used in the analysis.
4.1 Transverse momentum and transverse mass spectra
Transverse momentum spectra in rapidity slices of \(K^{+}\) produced in the 20% most central Be+Be collisions. Rapidity values given in the legends correspond to the middle of the corresponding interval. Black dots (blue squares) show results of the \(\text {d}E/\text {d}x\) (tof-\(\text {d}E/\text {d}x\) ) analysis, respectively. Shaded bands show systematic uncertainties
Transverse momentum spectra in rapidity slices of \(K^{-}\) produced in the 20% most central Be+Be collisions. Rapidity values given in the legends correspond to the middle of the corresponding interval. Black dots (blue squares) show results of the \(\text {d}E/\text {d}x\) (tof-\(\text {d}E/\text {d}x\) ) analysis, respectively. Shaded bands show systematic uncertainties
Transverse momentum spectra in rapidity slices of \(\pi ^{+}\) produced in the 20% most central Be+Be collisions. Rapidity values given in the legends correspond to the middle of the corresponding interval. Presented results were obtained with the \(\text {d}E/\text {d}x\) analysis method. Shaded bands show systematic uncertainties
Transverse momentum spectra in rapidity slices of \(\pi ^{-}\) produced in the 20% most central Be+Be collisions. Rapidity values given in the legends correspond to the middle of the corresponding interval. Presented results were obtained with the \(\text {d}E/\text {d}x\) analysis method. Shaded bands show systematic uncertainties
Transverse momentum spectra in rapidity slices of protons produced in the 20% most central Be+Be collisions. Rapidity values given in the legends correspond to the middle of the corresponding interval. Presented results were obtained with the \(\text {d}E/\text {d}x\) analysis method. Shaded bands show systematic uncertainties
Transverse momentum spectra in rapidity slices of antiprotons produced in the 20% most central Be+Be collisions. Rapidity values given in the legends correspond to the middle of the corresponding interval. Presented results were obtained with the \(\text {d}E/\text {d}x\) analysis method. Shaded bands show systematic uncertainties
Resulting double differential spectra as a function of transverse momentum \(p_{T}\) in intervals of rapidity \(y\) of \(K^{+}\), \(K^{-}\), \(\pi ^{+}\), \(\pi ^{-}\), p and \(\bar{p}\) produced at 19A, 30A, 40A, 75A, 150\(A\,{\mathrm{Ge} \mathrm{V}}\!/\!c\) beam momentum are shown in Figs. 20, 21, 22, 23, 24, 25. Entries in Fig. 15 below mid-rapidity were reflected in order to fill the gaps in acceptance as far as possible. Spectra in successive rapidity intervals were scaled by appropriate factors for better visibility. Vertical bars on data points correspond to statistical, shaded bands to systematic uncertainties. Systematic uncertainties plotted in the logarithmic scale look small.
Transverse momentum spectra of shown in Figs. 20, 21, 22, 23, 24, 25 were parametrized by an exponential function [29, 30]:
where m is the particle mass and S and T are the yield integral and the inverse slope parameter, respectively. The functions (Eq. 21) fitted at mid-rapidity are shown together with the data points for \(K^{+}\) and \(K^{-}\) in Fig. 26 and for protons and antiprotons in Fig. 27.
Transverse momentum spectra of \(K^{+}\) (left) and \(K^{-}\) (right) mesons produced in \(0<y< 0.2 (-0.2<y < 0.0\) for 30\(A\,{\mathrm{Ge} \mathrm{V}}\!/\!c\)) in the 20% most central Be+Be collisions. Colored lines represent the fitted function (Eq. 21)
Transverse momentum spectra of protons (left) and antiprotons (right) produced in \(0<y < 0.2\) in the 20% most central Be+Be collisions. Colored lines represent the fitted function (Eq. 21)
The fitted inverse slope parameter is plotted in Fig. 28 as a function of the rapidity. The value of T results from an interplay of the kinetic freeze-out temperature when final state rescattering stops as well as of the radial expansion flow. Thus the \(p_T\) spectra need not show a strictly exponential decrease. Note that results are only plotted for those rapidity intervals for which there were more than 6 data points in the \(p_{T}\)-distribution.
Inverse slope parameter T of the \(p_T\) distributions as function of rapidity y fitted with Eq. (21) in the full available \(p_T\) range. Results for the 20% most central Be+Be collisions at beam momenta of 19A, 30A 40A, 75A and 150\(A\,{\mathrm{Ge} \mathrm{V}}\!/\!c\). Only statistical uncertainties are shown, systematic uncertainties are estimated at about 10%. Dashed curves show predictions of the Epos [16, 31] model for comparison
Rapidity spectra dn/dy of \(K^{+}\) and \(K^{-}\) were obtained by integration of the transverse momentum distributions shown in Fig. 26. Both dn/dy and the corresponding inverse slope parameter T at mid-rapidity are tabulated in Table 3.
A step-like structure in the energy dependence of the inverse slope parameter T of kaons at mid-rapidity was predicted [32] at the onset of deconfinement (step). In this scenario it is caused by the softness of the equation of state of the mixed phase of hadrons and partons stalling the expansion of the initial state with rising energy density. Kaon \(p_{T}\) spectra are well described by a simple exponential because in contrast to pion \(p_{T}\) spectra they are not affected strongly by resonance decay products. As seen in Fig. 29 from a compilation of published results such a step is observed in central collisions of heavy nuclei (Au and Pb) in the SPS energy range which is consistent with the onset of the phase transition [11, 12]. The new NA61/SHINE measurements from central Be+Be collisions are similar to published results from inelastic p+p interactions. This indicates that not much expansion flow is created in the small Be+Be collision system at SPS energies. Although there appears to be a similar step feature, the values of T are much smaller than the ones from collisions of heavy nuclei. The intriguing similarity between the energy dependence of T between p+p and Pb+Pb collisions is discussed in detail in Ref. [10].
4.2 Rapidity spectra and mean multiplicities
Figure 30 presents rapidity spectra for all studied particle species at all available collision energies. Measurements of the \(\pi ^-\) rapidity distributions employing the \(h^-\) method indicated slight forward-backward asymmetry [4]. The acceptance for particle identification does not extend into the backward hemisphere making such a test impossible. The forward-backward asymmetry of \(\pi ^-\) rapidity spectra reported in Ref. [4] was at the level of 5%. Since this is smaller than the systematic uncertainty no correction was applied. In fact, Gaussian functions provide a good description for all particle types except for protons. The latter exhibit a strong leading particle effect. Mean forward multiplicities (protons excepted), were obtained by summing the measured points and adding an extrapolation estimated from the fitted Gaussian functions which are shown by the curves in Fig. 30. Assuming forward-backward symmetry total mean multiplicities were calculated and are summarized in Table 5.
The simplest model of nucleus-nucleus collisions, the Wounded Nucleon Model [20], suggests that total produced particle multiplicities scale approximately with the respective ratios of wounded nucleons. This ratio was derived using the Epos model and the NA61/SHINE centrality selection procedure for Be+Be collisions with the 20% smallest number of wounded nucleons (see Sect. 3.1) and inelastic p+p collisions Ref. [10]. The resulting ratio of about 4 is close to the experimental ratios calculated from the total multiplicities listed in Table 5 for Be+Be and Ref. [10] for p+p interactions.
The determination of the mean multiplicity is more complicated for protons due to the rapid rise of the yield towards beam rapidity and the lack of measurements in this region. A comparison of the rapidity distributions obtained in this analysis with predictions of the Urqmd [33, 34], Epos [16, 31], Ampt [35,36,37], Phsd 4.0 [38, 39] and Smash 1.6 [40, 41] models at 40A and 150\(A\,{\mathrm{Ge} \mathrm{V}}\!/\!c\) beam momentum is shown in Fig. 31. Note that protons with \(p_{T} \approx 0\) and beam rapidity are assumed to be spectator protons and rejected in the model calculations. To determine the mean proton multiplicity one may calculate extrapolation factors from the models as the ratio of mean multiplicity to that in the region covered by the measurements. The multiplicity of protons in the region of measurements and the extrapolation factor obtained from the Epos model are given in Table 4. Mean proton multiplicities are not provided since the values will be strongly model dependent.
The total multiplicities for all studied particle species are plotted as a function of collision energy in Fig. 32 and compared to the corresponding results from inelastic p+p interactions. Particle yields in Be+Be collisions are higher by approximately a factor of four consistent with expectations from the wounded nucleon model.
4.3 K/\(\pi \) ratio
The K/\(\pi \) ratio at SPS energies was shown to be a good measure of the strangeness to entropy ratio [42] which is different in the confined phase (hadrons) and the QGP (quarks, anti-quarks and gluons). A maximum (horn) and a subsequent plateau in the energy dependence of the \(K^+\)/\(\pi ^+\) ratio was observed by the NA49 experiment in the SPS energy range (see Fig. 33) in central Pb+Pb collisions [11, 12] and interpreted as one of the indications of the onset of deconfinement [32].
The NA61/SHINE experiment studies the evolution of this signal with respect to the size of the collision system in order to find when deconfinement starts occurring. In p+p interactions and Be+Be collisions the \(K^{\pm }\) and \(\pi ^-\) yield can be measured by NA61/SHINE over all of phase space. This is not possible for the \(\pi ^+\) yield which has to be derived from the \(\pi ^-\) yield using an isospin correction. Be+Be collisions, however, are isospin symmetric and therefore the mean multiplicities of \(\pi ^+\) and \(\pi ^-\) have to be the same. The \(\pi ^-\) yields for the 20% most central Be+Be collisions were calculated by scaling the \(\pi ^-\) multiplicity published in Ref. [4] by the ratio of \(\langle W \rangle \) for the 20% (see Table 2) and the 5% most central Be+Be collisions Ref. [4]. The results for the energy dependence of the \(K^+\)/\(\pi ^+\) ratio from the 20% most central Be+Be collisions are shown in Fig. 33 together with measurements in inelastic p+p, central Pb+Pb collisions and other reactions. The results from Be+Be collisions do not yet deviate from those in p+p reactions and are about a factor two lower for the total yield ratio than the values found in central Pb+Pb and Au+Au collisions. In particular, no horn structure is observed for the small collision systems.
5 Comparison with models
This subsection compares observables, sensitive to the onset of deconfinement, between measurements and model predictions. The Epos 1.99 [43], Urqmd 3.4 [33, 34], Ampt 1.26 [35,36,37], Phsd 4.0 [38, 39] and Smash 1.6 [40, 41] models were chosen for this study. In Epos the reaction proceeds from the excitation of strings according to Gribov–Regge theory to string fragmentation into hadrons. Urqmd starts with a hadron cascade based on elementary cross sections for resonance production which either decay (mostly at low energies) or are converted into strings which fragment into hadrons (mostly at high energies). Ampt uses the heavy ion jet interaction generator (Hijing) for generating the initial conditions, Zhang’s parton cascade for modeling partonic scatterings and the Lund string fragmentation model or a quark coalescence model for hadronization. Phsd is a microscopic offshell transport approach that describes the evolution of a relativistic heavy-ion collision from the initial hard scatterings and string formation through the dynamical deconfinement phase transition to the quark-gluon plasma as well as hadronization and the subsequent interactions in the hadronic phase. Smash uses the hadronic transport approach where the free parameters of the string excitation and decay are tuned to match the experimental measurements in inelastic p+p collisions. Selection of events in all model calculations follows the procedure for central collisions to which experimental results correspond to, see Sect. 3.1. This is particularly important when comparisons of yields with measurements are to be performed.
Comparisons of the \(p_T\) spectra at midrapidity and rapidity spectra of \(K^+\), \(K^-\) and \(\pi ^-\) in the 20% most central Be+Be collisions at 150\(A\,{\mathrm{Ge} \mathrm{V}}\!/\!c\) are shown in Figs. 34 and 35 respectively. Positively charged kaon \(p_{T}\) distributions are overpredicted by Ampt, Phsd and Epos, whereas Urqmd and Smash predictions are closer to the measurements (Fig. 34 (top-left)). In the case of negatively charged K mesons (Fig. 34 (top-right)), Epos next to Smash describes \(p_{T}\) the distribution well, while Urqmd, Phsd and Ampt overestimate the measured spectrum. For \(\pi ^-\) one finds a similar behaviour (Fig. 34 (bottom-left)). In the case of protons (Fig. 34 (bottom-right)), again all models overestimate the \(p_{T}\) spectrum. Moreover, Epos, Smash and Phsd predict narrower \(p_{T}\) spectra than measured, which result in smaller value of inverse slope parameter than the one obtained by the data fit (see Figs. 28 and 34).
Predictions for the rapidity distributions are too high for all models (Figs. 35 and 31).
The energy dependence of the inverse slope parameter of \(p_T\) spectra at mid-rapidity of positively (left) and negatively (right) charged K mesons for central Be+Be, Pb+Pb and Au+Au collisions as well as inelastic p+p interactions. Both statistical (vertical bars) and systematic uncertainties (shaded bands) are shown
Rapidity spectra of \(K^{+}\), \(K^{-}\), \(\pi ^{+}\), \(\pi ^{-}\), protons and antiprotons produced in the 20% most central Be+Be collisions. Curves depict Gaussian fits used to determine mean multiplicities. Spectra for different beam momenta were scaled by the following factors for better visibility: 150\(A\,{\mathrm{Ge} \mathrm{V}}\!/\!c\) by factor 8, 75\(A\,{\mathrm{Ge} \mathrm{V}}\!/\!c\) by factor 6, 40\(A\,{\mathrm{Ge} \mathrm{V}}\!/\!c\) by factor 4, 30\(A\,{\mathrm{Ge} \mathrm{V}}\!/\!c\) by factor 2 and 19\(A\,{\mathrm{Ge} \mathrm{V}}\!/\!c\) by factor 4
Proton rapidity distribution in the 20% most central Be+Be collisions at 40A and 150\(A\,{\mathrm{Ge} \mathrm{V}}\!/\!c\) compared with predictions of the Epos 1.99 [16, 31] (blue dashed line), Urqmd 3.4 [33, 34] (black solid line), Ampt 1.26 [35,36,37] (violet dotted line), Phsd 4.0 [38, 39] (brown dashed-dotted line) and Smash 1.6 (red dashed-double dotted line) [40, 41] models
Collision energy dependence of mean multiplicities of \(K^{+}\), \(K^{-}\), \(\pi ^{+}\), \(\pi ^{-}\), protons and antiprotons produced in the 20% most central Be+Be collisions. Results from inelastic p+p interactions [8, 9] are plotted for comparison. Both statistical (vertical bars) and systematic uncertainties (shaded bands) are shown. Gray points and band correspond to \(\pi ^{-}\) results obtained via \(h^{-}\) method. Results for \(\langle \pi ^{-} \rangle \) obtained using the \(h^{-}\) method [4] taking into account backward-forward asymmetry of the rapidity distribution. They were scaled from 5% centrality to 20% based on the ratio of wounded nucleons estimated using the Epos model
The energy dependence of the \(K^+\)/\(\pi ^+\) particle yields ratio at mid-rapidity (left) and full acceptance (right) for the 20% most central Be+Be, central Pb+Pb and Au+Au collisions, as well as inelastic p+p interactions. Both statistical (vertical bars) and systematic uncertainties (shaded bands) are shown
Comparison of the \(p_T\) spectra of \(K^+\) (top-left), \(K^-\) (top-right), \(\pi ^-\) (bottom-left) and p (bottom-right) at mid-rapidity for the 20% most central Be+Be collisions at 150\(A\,{\mathrm{Ge} \mathrm{V}}\!/\!c\) with models: Epos 1.99 (blue dashed line), Urqmd 3.4 (black solid line), Ampt 1.26 (violet dotted line), Phsd 4.0 (brown dashed-dotted line) and Smash 1.6 (red dashed-double dotted line)
Comparison of the \(K^+\) (left), \(K^-\) (right) and \(\pi ^-\) (bottom) rapidity spectra for the 20% most central Be+Be collisions at 150\(A\,{\mathrm{Ge} \mathrm{V}}\!/\!c\) with models: Epos 1.99 (blue dashed line), Urqmd 3.4 (black solid line), Ampt 1.26 (violet dotted line), Phsd 4.0 (brown dashed-dotted line) and Smash 1.6 (red dashed-double dotted line)
The measurements of the inverse slope parameter T in midrapidity in the 20% most central Be+Be collisions are shown versus collision energy for \(K^+\) and \(K^-\) in Fig. 36 (top) together with model predictions. Except for Ampt the predictions cluster around the measurements. Urqmd and Phsd feature a hadron rescattering phase as the last step of the system’s evolution. The values of T for \(\pi ^-\) shown in Fig. 36 (bottom-left) are well reproduced by Epos, Urqmd and Smash whereas the predictions of Ampt and Phsd are too high by 30–50 \({\hbox {Me} \hbox {V}}\). The limited range of studied beam momentum does not allow a definite conclusion on a possible step structure. Urqmd predictions for inverse slope parameter of protons reproduce the experimental data quite well. Ampt model significantly overestimates T value for p, while Phsd, Epos and Smash present opposite tendency. Comparison of T as a function of y with Epos predictions is presented in Fig. 28 for \(K^+\), \(K^-\), \(\pi ^+\), \(\pi ^-\), p and \(\overline{p}\). In general, Epos underestimates T for all types of identified hadrons except pions, for which simulation follows the experimental results.
Comparison of the energy dependence of the inverse slope parameter T of \(K^+\) (top-left), \(K^-\) (top-right), \(\pi ^-\) (bottom-left) and p (bottom-right) spectra at mid-rapidity for the 20% most central Be+Be collisions with models: Epos 1.99 (blue dashed line), Urqmd 3.4 (black solid line), Ampt 1.26 (violet dotted line), Phsd 4.0 (brown dashed-dotted line) and Smash 1.6 (red dashed-double dotted line)
Comparison of the energy dependence of \(K^+\)/\(\pi ^+\) (left) and \(K^-\)/\(\pi ^-\) (right) yields ratio at mid-rapidity for the 20% most central Be+Be collisions with models: Epos 1.99 (blue dashed line), Urqmd 3.4 (black solid line), Ampt 1.26 (violet dotted line), Phsd 4.0 (brown dashed-dotted line) and Smash 1.6 (red dashed-double dotted line)
Comparison of the energy dependence of \(K^+\)/\(\pi ^+\) (left) and \(K^-\)/\(\pi ^-\) (right) mean multiplicities ratio for the 20% most central Be+Be collisions with models: Epos 1.99 (blue dashed line), Urqmd 3.4 (black solid line), Ampt 1.26 (violet dotted line), Phsd 4.0 (brown dashed-dotted line) and Smash 1.6 (red dashed-double dotted line)
Finally, the energy dependence of the ratios of kaon to pion yields are compared to model predictions. Figure 37 shows the mid-rapidity results for \(K^+\)/\(\pi ^+\) (left) and \(K^-\)/\(\pi ^-\) (right), respectively. Figure 38 displays the ratios of total mean multiplicities. Unlike particle yields, particle ratios are not sensitive to the details of the event selection assuming that the shapes of the spectra do not change significantly in the studied centrality range. For \(K^+\)/\(\pi ^+\) Epos, Urqmd and Smash provide a good description whereas Ampt and Phsd overpredict the data. For \(K^-\)/\(\pi ^-\) all models are close to the experimental results.
A detailed investigation of the effects included in the models and their impact on the predictions of the experimental results is beyond the scope of this paper. However, the new NA61/SHINE measurements provide useful input for future refinements of the models.
6 Summary and conclusions
This paper reports measurements by the NA61/SHINE experiment at the CERN SPS of spectra and mean multiplicities of \(\pi ^\pm , K^\pm , p\) and \(\bar{p}\) produced in the 20 % most central \(^7\)Be+\(^9\)Be collisions at beam momenta of 19A, 30A, 40A, 75A and 150\(A\,{\mathrm{Ge} \mathrm{V}}\!/\!c\). This is the lightest nucleus-nucleus system investigated in the system size scan of NA61/SHINE. In this program data were also recorded from Ar+Sc, Xe+La and Pb+Pb collisions for which the analysis is ongoing. Results on central \(^7\)Be+\(^9\)Be collisions were found to be similar to those from inelastic p+p interactions for shapes of transverse momentum and rapidity spectra as well as the \(K^+\)/\(\pi ^+\) ratio. However, particle yields are higher by approximately a factor of four consistent with expectations from the wounded nucleon model. Summarising, neither measurements nor models show indications of a horn structure at low SPS energy for small collision systems in contrast to the results from central Pb+Pb interactions, for detailed discussion on possible indications of onset of deconfinement in small systems see Ref. [10].
The results were compared with predictions of the models: Epos 1.99, Urqmd 3.4, Ampt 1.26, Phsd 4.0 and Smash 1.6. None of the models reproduces all features of the presented results.
Data Availability Statement
This manuscript has no associated data or the data will not be deposited. [Authors’ comment: There are no external data associated with the manuscript.]
Change history
30 January 2023
An Erratum to this paper has been published: https://doi.org/10.1140/epjc/s10052-023-11234-2
References
N. Antoniou et al., [NA61/SHINE Collab.], Study of hadron production in hadron nucleus and nucleus nucleus collisions at the CERN SPS, Tech. Tep., CERN, CERN-SPSC-2006-034 (2006)
M. Gazdzicki, M. Gorenstein, P. Seyboth, Recent developments in the study of deconfinement in nucleus–nucleus collisions. Int. J. Mod. Phys. E 23, 1430008 (2014). https://doi.org/10.1142/S0218301314300082. arXiv:1404.3567 [nucl-ex]
A. Aduszkiewicz, [NA61/SHINE Collab.], Beam momentum scan with Pb+Pb collisions, Tech. Rep. CERN-SPSC-2015-038. SPSC-P-330-ADD-8, CERN, Geneva (2015). https://cds.cern.ch/record/2059811
A. Acharya et al., [NA61/SHINE Collab.], Measurements of \(\pi ^-\) production in \(^7\)Be+\(^9\)Be collisions at beam momenta from 19\(A\) to 150\(A\)GeV/\(c\) in the NA61/SHINE experiment at the CERN SPS. arXiv:2008.06277 [nucl-ex]
R. Poberezhnyuk, M. Gazdzicki, M. Gorenstein, Statistical model of the early stage of nucleus–nucleus collisions with exact strangeness conservation. Acta Phys. Polon. B 46(10), 1991 (2015). https://doi.org/10.5506/APhysPolB.46.1991. arXiv:1502.05650 [nucl-th]
A. Motornenko, V. Begun, V. Vovchenko, M. Gorenstein, H. Stoecker, Hadron yields and fluctuations at energies available at the CERN Super Proton Synchrotron: system-size dependence from Pb + Pb to \(p+p\) collisions. Phys. Rev. C 99(3), 034909 (2019). https://doi.org/10.1103/PhysRevC.99.034909. arXiv:1811.10645 [nucl-th]
M. Gazdzicki, O. Hansen, Hadron production in nucleon–nucleon collisions at 200-GeV/c: a compilation. Nucl. Phys. A 528, 754–770 (1991). https://doi.org/10.1016/0375-9474(91)90261-4
N. Abgrall et al., [NA61/SHINE Collab.], Measurement of negatively charged pion spectra in inelastic p+p interactions at \(p_{lab}\) = 20, 31, 40, 80 and 158 GeV/c. Eur. Phys. J. C 74, 2794 (2014). https://doi.org/10.1140/epjc/s10052-014-2794-6. arXiv:1310.2417 [hep-ex]
A. Aduszkiewicz et al., [NA61/SHINE Collab.], Measurements of \(\pi ^\pm \) , K\(^\pm \) , p and \({\bar{\text{p}}}\) spectra in proton–proton interactions at 20, 31, 40, 80 and 158 \( \text{ GeV }/c\) with the NA61/SHINE spectrometer at the CERN SPS. Eur. Phys. J. C 77(10), 671 (2017). https://doi.org/10.1140/epjc/s10052-017-5260-4. arXiv:1705.02467 [nucl-ex]
A. Aduszkiewicz et al., [NA61/SHINE Collab.], Proton–proton interactions and onset of deconfinement. Phys. Rev. C 102, 011901 (2020). https://doi.org/10.1103/PhysRevC.102.011901. arXiv:1912.10871 [hep-ex]
S. Afanasiev et al., [NA49 Collab.], Energy dependence of pion and kaon production in central Pb + Pb collisions. Phys. Rev. C 66, 054902 (2002). https://doi.org/10.1103/PhysRevC.66.054902
C. Alt et al., [NA49 Collab.], Pion and kaon production in central Pb + Pb collisions at 20-A and 30-A-GeV: evidence for the onset of deconfinement. Phys. Rev. C 77, 024903 (2008). https://doi.org/10.1103/PhysRevC.77.024903
N. Abgrall et al., [NA61 Collab.], NA61/SHINE facility at the CERN SPS: beams and detector system. JINST 9, P06005 (2014). https://doi.org/10.1088/1748-0221/9/06/P06005. arXiv:1401.4699 [physics.ins-det]
O. Berrig et al., [NA61/SHINE Collab.], The 2010 test of secondary light ion beams, Tech. Rep., CERN, CERN-SPSC-2011-005 (2011)
D. Banas, A. Kubala-Kukus, M. Rybczynski, I. Stabrawa, G. Stefanek, Influence of target material impurities on physical results in relativistic heavy-ion collisions. Eur. Phys. J. Plus 134(1), 44 (2019). https://doi.org/10.1140/epjp/i2019-12465-9. arXiv:1808.10377 [nucl-ex]
K. Werner, F.-M. Liu, T. Pierog, Parton ladder splitting and the rapidity dependence of transverse momentum spectra in deuteron-gold collisions at RHIC. Phys. Rev. C 74, 044902 (2006). https://doi.org/10.1103/PhysRevC.74.044902
R. Brun, F. Carminati, Geant detector description and simulation tool, CERN program library long writeup w5013 (1993). http://wwwasdoc.web.cern.ch/wwwasdoc/geant/geantall.html
A. Seryakov, [NA61/SHINE Collab.], PSD acceptance maps for event selection. CERN EDMS (2017). https://edms.cern.ch/document/1867336/1
E. Kaptur, [NA61/SHINE Collab.], Analysis of collision centrality and negative pion spectra in \(^7\)Be+\(^9\)Be interactions at CERN SPS energy range, Ph.D thesis (2017). https://edms.cern.ch/document/2004086/1
A. Bialas, M. Bleszynski, W. Czyz, Multiplicity distributions in nucleus–nucleus collisions at high-energies. Nucl. Phys. B 111, 461 (1976). https://doi.org/10.1016/0550-3213(76)90329-1
W. Broniowski, M. Rybczynski, P. Bozek, GLISSANDO: Glauber initial-state simulation and more. Comput. Phys. Commun. 180, 69 (2009). https://doi.org/10.1016/j.cpc.2008.07.016
M. Kuich, [NA61/SHINE Collab.], Kaon production in mid-rapidity in Be+Be collisions at the CERN SPS, Ph.D. thesis (2019). https://edms.cern.ch/document/2150851/1
M. van Leeuwen, [NA49 Collab.], Energy dependence of particle production in nucleus nucleus collisions at the CERN SPS, in Proceedings, 38th Rencontres de Moriond on QCD and High-Energy Hadronic Interactions (2003). arXiv:nucl-ex/0306004
M. van Leeuwen, A practical guide to de/dx analysis in na49 (2008)
S. Afanasev et al., [NA49 Collab.], The NA49 large acceptance hadron detector. Nucl. Instrum. Methods A 430, 210–244 (1999). https://doi.org/10.1016/S0168-9002(99)00239-9
T. Anticic et al., [NA49 Collab.], Antideuteron and deuteron production in mid-central Pb+Pb collisions at 158\(A\) GeV. Phys. Rev. C 85, 044913 (2012). arXiv:1111.2588 [nucl-ex]. https://doi.org/10.1103/PhysRevC.85.044913
A. Rustamov, M. Gorenstein, Identity method for moments of multiplicity distribution. Phys. Rev. C 86, 044906 (2012). https://doi.org/10.1103/PhysRevC.86.044906. arXiv:1204.6632 [nucl-th]
A. Aduszkiewicz et al., [NA61/SHINE Collab.], Production of \(\Lambda \)-hyperons in inelastic p+p interactions at 158 \({\rm GeV}\!/\!c\). Eur. Phys. J. C 76(4) 198, (2016). https://doi.org/10.1140/epjc/s10052-016-4003-2. arXiv:1510.03720 [hep-ex]
R. Hagedorn, Statistical thermodynamics of strong interactions at high energies. 3. Heavy-pair (quark) production rates. Nuovo Cim. Suppl. 6, 311–354 (1968)
W. Broniowski, W. Florkowski, L.Y. Glozman, Update of the Hagedorn mass spectrum. Phys. Rev. D 70, 117503 (2004). https://doi.org/10.1103/PhysRevD.70.117503. arXiv:hep-ph/0407290
T. Pierog, R. Ulrich, and private communication, EPOS 1.99 version CRMC (2018). https://web.ikp.kit.edu/rulrich/crmc/html
M. Gazdzicki, M.I. Gorenstein, On the early stage of nucleus–nucleus collisions. Acta Phys. Polon. B 30, 2705 (1999). arXiv:hep-ph/9803462
S. Bass et al., Microscopic models for ultrarelativistic heavy ion collisions. Prog. Part. Nucl. Phys. 41, 255–369 (1998). https://doi.org/10.1016/S0146-6410(98)00058-1, arXiv:nucl-th/9803035
M. Bleicher et al., Relativistic hadron hadron collisions in the ultrarelativistic quantum molecular dynamics model. J. Phys. G 25, 1859–1896 (1999). https://doi.org/10.1088/0954-3899/25/9/308. arXiv:hep-ph/9909407
Z.-W. Lin, C.M. Ko, B.-A. Li, B. Zhang, S. Pal, Multiphase transport model for relativistic heavy ion collisions. Phys. Rev. C 72, 064901 (2005). https://doi.org/10.1103/PhysRevC.72.064901
Z.-W. Lin, Evolution of transverse flow and effective temperatures in the parton phase from a multiphase transport model. Phys. Rev. C 90, 014904 (2014). https://doi.org/10.1103/PhysRevC.90.014904
B. Zhang, C.M. Ko, B.-A. Li, Z. Lin, Multiphase transport model for relativistic nuclear collisions. Phys. Rev. C 61, 067901 (2000). https://doi.org/10.1103/PhysRevC.61.067901
W. Cassing, E.L. Bratkovskaya, Parton transport and hadronization from the dynamical quasiparticle point of view. Phys. Rev. C 78, 034919 (2008). https://doi.org/10.1103/PhysRevC.78.034919
W. Cassing, E. Bratkovskaya, Parton-hadron-string dynamics: an off-shell transport approach for relativistic energies. Nucl. Phys. A 831(3), 215–242 (2009). https://doi.org/10.1016/j.nuclphysa.2009.09.007. http://www.sciencedirect.com/science/article/pii/S0375947409005995
J. Mohs, S. Ryu, H. Elfner, Particle production via strings and baryon stopping within a hadronic transport approach. J. Phys. G 47(6), 065101 (2020). https://doi.org/10.1088/1361-6471/ab7bd1. arXiv:1909.05586 [nucl-th]
J. Weil, V. Steinberg, J. Staudenmaier, L.G. Pang, D. Oliinychenko, J. Mohs, M. Kretz, T. Kehrenberg, A. Goldschmidt, B. Bäuchle, J. Auvinen, M. Attems, H. Petersen, Particle production and equilibrium properties within a new hadron transport approach for heavy-ion collisions. Phys. Rev. C 94, 054905 (2016). https://doi.org/10.1103/PhysRevC.94.054905
M. Gazdzicki, M. Gorenstein, P. Seyboth, Onset of deconfinement in nucleus-nucleus collisions: review for pedestrians and experts. Acta Phys. Polon. B 42, 307–351 (2011). https://doi.org/10.5506/APhysPolB.42.307. arXiv:1006.1765 [hep-ph]
K. Werner, The hadronic interaction model EPOS. Nucl. Phys. Proc. Suppl. 175–176, 81–87 (2008). https://doi.org/10.1016/j.nuclphysbps.2007.10.012
Acknowledgements
We would like to thank the CERN EP, BE, HSE and EN Departments for the strong support of NA61/SHINE. This work was supported by the Hungarian Scientific Research Fund (Grant NKFIH 123842/123959), the Polish Ministry of Science and Higher Education (Grants 667/N-CERN/2010/0, NN 202 48 4339, NN 202 23 1837 and DIR/WK/2016/2017/10-1), the National Science Centre Poland (Grants 2014/14/E/ST2/00018, 2014/15/B/ST2 /02537 and 2015/18/M/ST2/00125, 2015/19/N/ST2 /01689, 2016/23/B/ST2/00692, 2017/25/N/ST2/02575, 2018/30/A/ST2/00226, 2018/31/G/ST2/03910), the Russian Science Foundation, Grants 16-12-10176 and 17-72-20045, the Russian Academy of Science and the Russian Foundation for Basic Research (Grants 08-02-00018, 09-02-00664 and 12-02-91503-CERN), the Russian Foundation for Basic Research (RFBR) funding within the research project no. 18-02-40086, the National Research Nuclear University MEPhI in the framework of the Russian Academic Excellence Project (contract No. 02.a03.21.0005, 27.08.2013), the Ministry of Science and Higher Education of the Russian Federation, Project "Fundamental properties of elementary particles and cosmology" No 0723-2020-0041, the European Union’s Horizon 2020 research and innovation programme under Grant Agreement No. 871072, the Ministry of Education, Culture, Sports, Science and Technology, Japan, Grant-in-Aid for Scientific Research (Grants 18071005, 19034011, 19740162, 20740160 and 20039012), the German Research Foundation (Grant GA 1480/8-1), the Bulgarian Nuclear Regulatory Agency and the Joint Institute for Nuclear Research, Dubna (bilateral contract No. 4799-1-18/20), Bulgarian National Science Fund (Grant DN08/11), Ministry of Education and Science of the Republic of Serbia (Grant OI171002), Swiss Nationalfonds Foundation (Grant 200020117913/1), ETH Research Grant TH-01 07-3 and the Fermi National Accelerator Laboratory (Fermilab), a U.S. Department of Energy, Office of Science, HEP User Facility managed by Fermi Research Alliance, LLC (FRA), acting under Contract No. DE-AC02-07CH11359 and the IN2P3-CNRS (France).
Author information
Authors and Affiliations
Consortia
Corresponding author
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.
Funded by SCOAP3
About this article
Cite this article
Acharya, A., Adhikary, H., Aduszkiewicz, A. et al. Measurements of \(\pi ^\pm \), \(K^\pm \), p and \(\bar{p}\) spectra in \(^7\)Be+\(^9\)Be collisions at beam momenta from 19A to 150A \({\mathrm{Ge} \mathrm{V}}\!/\!c\) with the NA61/SHINE spectrometer at the CERN SPS. Eur. Phys. J. C 81, 73 (2021). https://doi.org/10.1140/epjc/s10052-020-08733-x
Received:
Accepted:
Published:
DOI: https://doi.org/10.1140/epjc/s10052-020-08733-x