# Baryon Stopping and Charged Particle Distributions in Central Pb+Pb Collisions at 158 GeV per Nucleon

## Abstract

Net proton and negative hadron spectra for central Pb+Pb collisions at 158 GeV per nucleon at the CERN SPS were measured and compared to spectra from lighter systems. Net baryon distributions were derived from those of net protons, utilizing model calculations of isospin contributions as well as data and model calculations of strange baryon distributions. Stopping (rapidity shift with respect to the beam) and mean transverse momentum $`p_T`$ of net baryons increase with system size. The rapidity density of negative hadrons scales with the number of participant nucleons for nuclear collisions, whereas their $`p_T`$ is independent of system size. The $`p_T`$ dependence upon particle mass and system size is consistent with larger transverse flow velocity at midrapidity for Pb+Pb compared to SS central collisions.

PACS numbers: 25.75-q,25.75.Dw

Lattice QCD predicts that strongly interacting matter at an energy density greater than 1-2 GeV/fm<sup>3</sup> attains a deconfined and approximately chirally restored state known as the quark-gluon plasma (for an overview, see ). This state of matter existed in the early universe, and it may influence the dynamics of rotating neutron stars . The collision of nuclei at ultrarelativistic energies offers the possibility in the laboratory of creating strongly interacting matter at sufficiently high energy density to form a quark gluon plasma . Hadronic spectra from these reactions reflect the dynamics of the hot and dense zone formed in the collision. The baryon density, established early in the reaction, is an important factor governing the evolution of the system . Comparison of model predictions with measured rapidity and transverse momentum distributions and correlation functions constrains the possible dynamical scenarios of the reaction , such as those for longitudinal and transverse flow . In addition, the mechanism by which the incoming nucleons lose momentum during the collision (baryon stopping ) is an important theoretical problem , and the measurement of baryon stopping in heavy ion collisions provides essential data on this question.

In this letter, we present measurements by the NA49 collaboration of rapidity and transverse momentum distributions of participating baryons and negative hadrons over a large fraction of phase space for central Pb+Pb collisions at 158 GeV per nucleon ($`\sqrt{s_{NN}}=17.2`$ GeV). Hadronic spectra from SS collisions at 200 GeV per nucleon  and N+N (nucleon-nucleon) collisions at 200 and 400 GeV  serve as important references, helping to identify effects that depart from those expected from the linear superposition of many N+N collisions.

The NA49 apparatus is described in . A beam of <sup>208</sup>Pb ions at 158 GeV per nucleon ($`y_{lab}=5.8`$) from the CERN SPS struck a 224 mg/cm<sup>2</sup> thick natural Pb target. The 5% most central collisions were selected by measurement of the forward-going energy in the phase space occupied by the projectile spectator nucleons. The reaction products passed through a dipole magnetic field, and tracks from charged particles used for this analysis were measured in two large Main Time Projection Chambers (MTPCs) placed downstream of the magnets on either side of the beam axis. Over 1000 tracks were measured in each event. Identification of protons utilized the specific ionization ($`dE/dx`$) of the gas of the MTPC . A relative $`dE/dx`$ resolution of 6% was achieved.

The evolution of the incoming baryons (baryon stopping) was studied through measurement of the difference of the proton and antiproton distributions (net protons, denoted $`p\overline{p}`$) to eliminate the effect of baryon-antibaryon pair production. Meson production was studied through the yield of negative hadrons ($`h^{}`$), comprising primarily $`\pi ^{}`$, with an admixture of $`K^{}`$ and $`\overline{p}`$. $`5\times 10^4`$ central events were used for the $`p\overline{p}`$ analysis and $`7\times 10^3`$ central events were used for the $`h^{}`$ analysis.

The net proton yield was determined by an identification technique designed to minimize the systematic errors over wide acceptance . For each phase space bin ($`y`$,$`p_T`$), the distribution of $`dE/dx`$ from negatively charged particles was subtracted from that of positively charged particles, resulting in a distribution that contained a peak with positive amplitude (more positive than negative tracks, principally due to the net proton yield), and a peak at higher $`dE/dx`$ with negative amplitude (net pions). In order to extract the particle yields, Gaussian functions were fitted to the two extrema in each difference distribution. A correction to the net proton yield for the difference between $`K^+`$ and $`K^{}`$ yields, based upon NA49 data , was 15% at $`y_{cm}<0`$ and was negligible for $`y_{cm}>1`$. The negative hadron yield was determined from the yield of all tracks from negatively charged particles excluding electrons, identified via $`dE/dx`$.

Corrections for detector acceptance and track reconstruction efficiency were calculated by embedding and reconstructing simulated tracks in real events . The tracking efficiency was greater than 95% over most of the acceptance, falling below 80% in regions of high track density and near the edges of the acceptance. Yield as a function of rapidity was determined only for bins having acceptance to $`p_T`$=0, with the exception of the two lowest rapidity bins for $`h^{}`$where extrapolation was based upon the adjacent bin having full coverage. At high $`p_T`$, the acceptance was limited to the region in which detector effects were well understood in this analysis, extending at least to $`p_T`$=2.0 GeV/c for all rapidity bins. Instrumental background due to secondary interactions with detector material was estimated to be 5% of the total measured yield, using a GEANT-based Monte Carlo simulation and the VENUS event generator , which reproduces transverse and forward energy distributions in central Pb+Pb collisions . Excluding the decay corrections discussed below, the systematic errors of all rapidity distributions were less than 10%.

The measured yields contain contributions from the products of weak decays that were incorrectly reconstructed as primary vertex tracks. The background correction to the $`h^{}`$ yield due to the decay of $`K_S^0`$ was estimated using NA49 data  to be 5% at midrapidity, decreasing strongly at higher rapidity. Background to the $`h^{}`$ yield due to $`\mathrm{\Lambda }`$ decay is less than 1%. The background correction to the $`p\overline{p}`$ yield is due to the decays of $`\mathrm{\Lambda }`$, $`\overline{\mathrm{\Lambda }}`$, $`\mathrm{\Sigma }^+`$, and $`\overline{\mathrm{\Sigma }}^{}`$, and was assessed using a GEANT-based simulation of the decay of $`\mathrm{\Lambda }`$ and $`\overline{\mathrm{\Lambda }}`$ and reconstruction of their charged decay products. ($`\mathrm{\Lambda }`$ and $`\overline{\mathrm{\Lambda }}`$ should be understood in this paper to include the contributions of $`\mathrm{\Sigma }^0`$ or $`\overline{\mathrm{\Sigma }^0}`$ as well as feeddown from weak decays, which are not distinguishable from primary $`\mathrm{\Lambda }`$ or $`\overline{\mathrm{\Lambda }}`$ in this analysis.) To investigate the dependence of this correction on the phase space distribution of the decaying particles, three different shapes of rapidity distribution for $`\mathrm{\Lambda }\overline{\mathrm{\Lambda }}`$ were used, based upon (I) predictions of the RQMD model  and (II) the VENUS model , and (III) a preliminary NA49 measurement . In each case the $`\mathrm{\Sigma }^+\overline{\mathrm{\Sigma }}^{}`$ yield was assumed to be 30% of the $`\mathrm{\Lambda }\overline{\mathrm{\Lambda }}`$ yield  and to have the same rapidity distribution. Strangeness is not conserved if both the hyperon yields in I or II and the measured kaon distribution are used. The hyperon yields in these cases were therefore scaled by the ratio of total charged kaon yields in the data and the model, thereby imposing strangeness conservation. In order to compare corrections due to $`\mathrm{\Lambda }\overline{\mathrm{\Lambda }}`$ distributions differing only in shape but not total yield, the distribution in III was scaled so that the ratio of $`\mathrm{\Lambda }`$+$`\mathrm{\Sigma }^0`$ to charged kaon total yield agreed with that from the models. The $`p\overline{p}`$ distribution corrected by other $`\mathrm{\Lambda }\overline{\mathrm{\Lambda }}`$ distributions can be derived from the figure.

Fig. 1, upper panel, shows the event normalized net proton yield as a function of rapidity for central Pb+Pb collisions, incorporating the three $`\mathrm{\Lambda }\overline{\mathrm{\Lambda }}`$ corrections and corresponding $`\mathrm{\Lambda }\overline{\mathrm{\Lambda }}`$ rapidity distributions. Also shown is the proton rapidity distribution for p+p collisions at 400 GeV , which is qualitatively different. The comparison highlights the importance of multiple collisions to baryon stopping in nuclear collisions.

To determine the net baryon ($`B\overline{B}`$) rapidity distribution, the contribution of the remaining net baryons, in addition to $`p\overline{p}`$ and $`\mathrm{\Lambda }\overline{\mathrm{\Lambda }}`$, must be estimated. Model calculations  indicate that the rapidity distribution of net neutrons follows that of net protons over most of phase space, with a $`7\%`$ larger yield. The contributions of $`\mathrm{\Sigma }^+\overline{\mathrm{\Sigma }}^{}`$ and $`\mathrm{\Sigma }^{}\overline{\mathrm{\Sigma }}^+`$ were accounted for by assuming the same rapidity distributions as $`\mathrm{\Lambda }\overline{\mathrm{\Lambda }}`$, and scaling the $`\mathrm{\Lambda }\overline{\mathrm{\Lambda }}`$ distribution by an empirical factor $`1.6\pm 0.1`$ derived from hadronic reactions . Multistrange baryons contribute less than 2% to $`B\overline{B}`$ , and were not accounted for. The net baryon yield $`B\overline{B}`$ was then calculated as

$$B\overline{B}=(2.07\pm 0.05)(p\overline{p})+(1.6\pm 0.1)(\mathrm{\Lambda }\overline{\mathrm{\Lambda }}).$$
(1)   
Fig. 1, lower panel, shows the $`B\overline{B}`$ rapidity distribution for Pb+Pb. The data points and errors are for $`\mathrm{\Lambda }\overline{\mathrm{\Lambda }}`$ correction I. The variation in $`B\overline{B}`$ due to the use of different $`\mathrm{\Lambda }\overline{\mathrm{\Lambda }}`$ corrections is smaller than that due to other sources of systematic error. The $`B\overline{B}`$ yield is less sensitive than the $`p\overline{p}`$ yield to the assumed $`\mathrm{\Lambda }\overline{\mathrm{\Lambda }}`$ distribution, because the $`\mathrm{\Lambda }\overline{\mathrm{\Lambda }}`$ distribution is removed from the measured $`p\overline{p}`$ distribution by the decay correction but added again in Eq. 1 to calculate $`B\overline{B}`$, with the two contributions approximately cancelling. Conclusions on baryon stopping are therefore based upon the net baryon rather than net proton rapidity distribution. Also shown in Fig. 1, lower panel, is the $`B\overline{B}`$ rapidity distribution for the 3% most central SS collisions at 200 GeV per nucleon , using a coefficient of $`2.0`$ for $`p\overline{p}`$ in Eq. 1 and with integral normalized to the number of nucleon participants in Pb+Pb. Within $`|y_{cm}|<2.5`$, there are $`352\pm 12`$ participant baryons for Pb+Pb and $`52\pm 3`$ for SS central collisions .

The $`B\overline{B}`$ rapidity distribution is narrower for Pb+Pb than for SS collisions, indicating increased baryon stopping for Pb+Pb collisions. (In the CM frame, $`y_{beam}=2.9`$ for Pb+Pb and $`3.0`$ for SS.) The median rapidity shift of the leading nucleon with respect to the beam in high energy p+Pb collisions is 2-2.5 units . Due to the symmetry of the nuclear reactions studied in this work, distributions of target and projectile nucleon rapidity shifts in full phase space cannot be determined separately, limiting the magnitude of a calculated rapidity shift even in the case of large baryon stopping. The mean rapidity shift relative to $`y_{beam}`$ for participant baryons within $`0<y_{cm}<2.5`$ is $`\mathrm{\Delta }y=1.76\pm 0.05`$ for Pb+Pb and $`\mathrm{\Delta }y=1.63\pm 0.16`$ for SS collisions.

Fig. 2 shows the event normalized rapidity distribution (assuming pion mass) of $`h^{}`$ in central Pb+Pb and SS  collisions, and isoscalar inelastic N+N  collisions. An independent analysis of Pb+Pb collisions using the NA49 Vertex TPCs found a distribution consistent with these data . The SS and N+N distributions were scaled relative to that for Pb+Pb by the ratio of the number of participant nucleons and a factor 0.96 to account for the energy dependence of average multiplicities measured in hadronic collisions . No correction for the net isospin difference between Pb+Pb and the other systems was applied. Under the assumption that the net isospin in the final state of a Pb+Pb collision is carried entirely by mesons, the additional scaling factor needed for comparison of the SS and N+N distributions to that of Pb+Pb is 1.12.

The enhancement of meson yield at midrapidity for nuclear collisions relative to N+N collisions has been noted previously . Taking into account the net isospin difference, the agreement between the Pb+Pb and scaled SS distributions is striking. This scaling of the yield of produced particles with number of participants, first observed in p+A collisions , is also observed in the collision of very heavy ions. Extrapolation of the $`h^{}`$ yield to full phase space was performed by taking into account the asymmetry in the $`h^{}`$ rapidity distribution (assuming pion mass) due to the contribution of $`K^{}`$ , resulting in an estimated total $`h^{}`$ yield of $`695\pm 30`$ particles. The number of $`h^{}`$ per participant nucleon pair, not adjusted for isospin, is $`4.0\pm 0.2`$ for Pb+Pb collisions, compared to $`3.6\pm 0.2`$ for SS collisions and $`3.22\pm 0.06`$ for isoscalar N+N collisions .

Fig. 3, upper panel, shows transverse mass spectra ($`m_T=\sqrt{m^2+p_{T}^{}{}_{}{}^{2}}`$) near midrapidity for $`h^{}`$ and $`p\overline{p}`$ for central Pb+Pb collisions. NA49 has previously reported the fit of an expanding hadronic source model  to midrapidity $`h^{}`$ and deuteron $`m_T`$ spectra and $`h^{}`$ correlation functions , giving a freezeout temperature of $`T=(120\pm 12)`$ MeV and transverse expansion velocity $`\beta _{}=(0.55\pm 0.12)`$. Fig. 3, upper panel, shows the fit of this model with fixed $`\beta _{}=0.55`$ to the $`h^{}`$ and $`p\overline{p}`$ spectra reported here. The resulting freezeout temperatures are $`T=(126\pm 2)`$ MeV for $`h^{}`$ and $`T=(118\pm 5)`$ MeV for $`p\overline{p}`$.

Fig. 3, lower panel, shows the rapidity dependence of $`p_T`$ calculated within $`0<p_T<2.5`$ GeV. For $`p\overline{p}`$ for central Pb+Pb collisions, $`p_T=(825\pm 37)`$ MeV at $`y_{cm}=0`$ and $`(600\pm 17)`$ MeV at $`y_{cm}=2.4`$, compared to the rapidity-averaged value $`p_T=(622\pm 26)`$ MeV for central SS collisions . For $`h^{}`$ for central Pb+Pb collisions, $`p_T=(385\pm 18)`$ MeV at $`y_{cm}=0`$, compatible with $`(377\pm 4)`$ MeV for central SS collisions . (In considering $`p_T`$ for $`h^{}`$ at high rapidity, note that the kinematic limit for pions for N+N collisions at 158 GeV falls below $`p_T=1`$ GeV near beam rapidity).

A similar characterization of $`p_T`$ distributions results from fitting $`(1/m_T)dn/dm_T`$ with a function $`A\mathrm{exp}(m_T/T)`$. A fit within $`0<m_Tm_0<0.8`$ GeV to the $`p\overline{p}`$ data near midrapidity for central Pb+Pb collisions yielded $`T=(308\pm 15)`$ (NA44 reported $`(289\pm 7)`$ MeV for protons ). NA35 found rapidity–averaged $`T=(235\pm 9)`$ MeV for central SS collisions . The increase in $`T`$ and $`p_T`$ with particle mass is consistent with a larger transverse radial flow velocity at midrapidity in Pb+Pb than SS collisions .

In summary, baryon stopping for central collisions at ultrarelativistic energies increases with system size. The integrated yield and rapidity density of negative hadrons exhibit scaling with the number of participant nucleons for nuclear collisions, and a small enhancement with respect to N+N collisions. At midrapidity, a large increase is observed in $`p_T`$ for the stopped baryons in Pb+Pb relative to SS collisions, with no significant increase in $`p_T`$ for negative hadrons. The increase in $`p_T`$ with particle mass is consistent with an increase in transverse radial flow velocity for heavier colliding systems.

This work was supported by the Director, Office of Energy Research, Division of Nuclear Physics of the Office of High Energy and Nuclear Physics of the U.S. Department of Energy under Contracts DE-AC03-76SF00098 and DE-FG02-91ER40609, the U.S. National Science Foundation, the Bundesministerium für Bildung und Forschung, Germany, the Alexander von Humboldt Foundation, the U.K. Engineering and Physical Sciences Research Council, the Polish State Committee for Scientific Research (2 P03B 02615 and 2 P03B 9913), the EC Marie Curie Foundation, and the Polish-German Foundation.

deceased <sup>$</sup>present address: Kent State Univ., Kent, OH, USA <sup>#</sup>present address: Physikalisches Institut, Universitaet Heidelberg, Germany <sup>+</sup>present address: Max-Planck-Institut für Physik, Munich, Germany