$P_{\rm H}$ measurements are shown as a function of hyperon $p_{\rm T}$ at $\sqrt{s_{\rm NN}}$=19.6 and 27 GeV. Statistical uncertainties are represented as lines while systematic uncertainties are represented as boxes. There is no observed dependence of $P_{\rm H}$ on $p_{\rm T}$ at $\sqrt{s_{\rm NN}}$=19.6 or 27 GeV, consistent with previous observations.

$P_{\rm H}$ measurements are shown as a function of hyperon $y$ at $\sqrt{s_{\rm NN}}$=19.6 and 27 GeV. Statistical uncertainties are represented as lines while systematic uncertainties are represented as boxes. The data set at $\sqrt{s_{\rm NN}}$=19.6 GeV takes advantage of STAR upgrades to reach larger $|y|$.

Example of combinatorial background subtracted invariant mass distributions and the extracted yields as a function of $\cos \theta^*$ for $\phi$ and $K^{*0}$ mesons. \textbf{a)} example of $\phi \rightarrow K^+ + K^-$ invariant mass distributions, with combinatorial background subtracted, integrated over $\cos \theta^*$; \textbf{b)} example of $K^{*0} (\overline{K^{*0}}) \rightarrow K^{-} \pi^{+} (K^{+} \pi^{-})$ invariant mass distributions, with combinatorial background subtracted, integrated over $\cos \theta^*$; \textbf{c)} extracted yields of $\phi$ as a function of $\cos \theta^*$; \textbf{d)} extracted yields of $K^{*0}$ as a function of $\cos \theta^*$.

Example of combinatorial background subtracted invariant mass distributions and the extracted yields as a function of $\cos \theta^*$ for $\phi$ and $K^{*0}$ mesons. \textbf{a)} example of $\phi \rightarrow K^+ + K^-$ invariant mass distributions, with combinatorial background subtracted, integrated over $\cos \theta^*$; \textbf{b)} example of $K^{*0} (\overline{K^{*0}}) \rightarrow K^{-} \pi^{+} (K^{+} \pi^{-})$ invariant mass distributions, with combinatorial background subtracted, integrated over $\cos \theta^*$; \textbf{c)} extracted yields of $\phi$ as a function of $\cos \theta^*$; \textbf{d)} extracted yields of $K^{*0}$ as a function of $\cos \theta^*$.

Example of combinatorial background subtracted invariant mass distributions and the extracted yields as a function of $\cos \theta^*$ for $\phi$ and $K^{*0}$ mesons. \textbf{a)} example of $\phi \rightarrow K^+ + K^-$ invariant mass distributions, with combinatorial background subtracted, integrated over $\cos \theta^*$; \textbf{b)} example of $K^{*0} (\overline{K^{*0}}) \rightarrow K^{-} \pi^{+} (K^{+} \pi^{-})$ invariant mass distributions, with combinatorial background subtracted, integrated over $\cos \theta^*$; \textbf{c)} extracted yields of $\phi$ as a function of $\cos \theta^*$; \textbf{d)} extracted yields of $K^{*0}$ as a function of $\cos \theta^*$.

Efficiency corrected $\phi$-meson yields as a function of cos$\theta$* and corresponding fits with Eq.1 in the method section.

Efficiency and acceptance corrected $K^{*0}$-meson yields as a function of cos$\theta$* and corresponding fits with Eq.4 in the method section.

$\phi$-meson $\rho_{00}$ obtained from 1st- and 2nd-order event planes. The red stars (gray squares) show the $\phi$-meson $\rho_{00}$ as a function of beam energy, obtained with the 2nd-order (1st-order) EP.

$\phi$-meson $\rho_{00}$ with respect to different quantization axes. $\phi$-meson $\rho_{00}$ as a function of beam energy, for the out-of-plane direction (stars) and the in-plane direction (diamonds). Curves are fits based on theoretical calculations with a $\phi$-meson field. The corresponding $G_s^{(y)}$ values obtained from the fits are shown in the legend.

$\rho_{00}$ as a function of transverse momentum for $\phi$ for different collision energies. The gray squares and red stars are results obtained with the 1st- and 2nd-order EP, respectively.

$\rho_{00}$ as a function of transverse momentum for $\phi$ for different collision energies. The gray squares and red stars are results obtained with the 1st- and 2nd-order EP, respectively.

$\rho_{00}$ as a function of transverse momentum for $\phi$ for different collision energies. The gray squares and red stars are results obtained with the 1st- and 2nd-order EP, respectively.

$\rho_{00}$ as a function of transverse momentum for $\phi$ for different collision energies. The gray squares and red stars are results obtained with the 1st- and 2nd-order EP, respectively.

$\rho_{00}$ as a function of transverse momentum for $\phi$ for different collision energies. The gray squares and red stars are results obtained with the 1st- and 2nd-order EP, respectively.

$\rho_{00}$ as a function of transverse momentum for $\phi$ for different collision energies. The gray squares and red stars are results obtained with the 1st- and 2nd-order EP, respectively.

$\rho_{00}$ as a function of transverse momentum for $\phi$ for different collision energies. The gray squares and red stars are results obtained with the 1st- and 2nd-order EP, respectively.

$\rho_{00}$ as a function of transverse momentum for $\phi$ for different collision energies. The gray squares and red stars are results obtained with the 1st- and 2nd-order EP, respectively.

$\rho_{00}$ as a function of transverse momentum for $\phi$ for different collision energies. The gray squares and red stars are results obtained with the 1st- and 2nd-order EP, respectively.

$\rho_{00}$ as a function of transverse momentum for $\phi$ for different collision energies. The gray squares and red stars are results obtained with the 1st- and 2nd-order EP, respectively.

$\rho_{00}$ as a function of transverse momentum for $\phi$ for different collision energies. The gray squares and red stars are results obtained with the 1st- and 2nd-order EP, respectively.

$\rho_{00}$ as a function of transverse momentum for $\phi$ for different collision energies. The gray squares and red stars are results obtained with the 1st- and 2nd-order EP, respectively.

$\rho_{00}$ as a function of transverse momentum for $K^{*0}$ for different collision energies.

$\rho_{00}$ as a function of transverse momentum for $K^{*0}$ for different collision energies.

$\rho_{00}$ as a function of transverse momentum for $K^{*0}$ for different collision energies.

$\rho_{00}$ as a function of transverse momentum for $K^{*0}$ for different collision energies.

$\rho_{00}$ as a function of transverse momentum for $K^{*0}$ for different collision energies.

$\rho_{00}$ as a function of transverse momentum for $K^{*0}$ for different collision energies.

$\rho_{00}$ as a function of transverse momentum for $K^{*0}$ for different collision energies.

$\rho_{00}$ as a function of centrality for $\phi$ (upper panels) and $K^{*0}$ (lower panels). The solid squares and stars are results for the $\phi$ meson, obtained with the 1st- and 2nd-order EP, respectively. The solid circles are results for the $K^{*0}$ meson, obtained with the 2nd-order EP.}

$\rho_{00}$ as a function of centrality for $\phi$ (upper panels) and $K^{*0}$ (lower panels). The solid squares and stars are results for the $\phi$ meson, obtained with the 1st- and 2nd-order EP, respectively. The solid circles are results for the $K^{*0}$ meson, obtained with the 2nd-order EP.}

$\rho_{00}$ as a function of centrality for $\phi$ (upper panels) and $K^{*0}$ (lower panels). The solid squares and stars are results for the $\phi$ meson, obtained with the 1st- and 2nd-order EP, respectively. The solid circles are results for the $K^{*0}$ meson, obtained with the 2nd-order EP.}

$\rho_{00}$ as a function of centrality for $\phi$ (upper panels) and $K^{*0}$ (lower panels). The solid squares and stars are results for the $\phi$ meson, obtained with the 1st- and 2nd-order EP, respectively. The solid circles are results for the $K^{*0}$ meson, obtained with the 2nd-order EP.}

$\rho_{00}$ as a function of centrality for $\phi$ (upper panels) and $K^{*0}$ (lower panels). The solid squares and stars are results for the $\phi$ meson, obtained with the 1st- and 2nd-order EP, respectively. The solid circles are results for the $K^{*0}$ meson, obtained with the 2nd-order EP.}

$\rho_{00}$ as a function of centrality for $\phi$ (upper panels) and $K^{*0}$ (lower panels). The solid squares and stars are results for the $\phi$ meson, obtained with the 1st- and 2nd-order EP, respectively. The solid circles are results for the $K^{*0}$ meson, obtained with the 2nd-order EP.}

$\rho_{00}$ as a function of centrality for $\phi$ (upper panels) and $K^{*0}$ (lower panels). The solid squares and stars are results for the $\phi$ meson, obtained with the 1st- and 2nd-order EP, respectively. The solid circles are results for the $K^{*0}$ meson, obtained with the 2nd-order EP.}

Global spin alignment measurement of $\phi$ and $K^{*0}$ vector mesons in Au+Au collisions at 0-20\% centrality. The solid squares and stars are results for the $\phi$ meson, obtained with the 1st- and 2nd-order EP, respectively. The solid circles are results for $K^{*0}$-meson, obtained with the 2nd-order EP.

Global spin alignment measurement of $\phi$ and $K^{*0}$ vector mesons in Au+Au collisions at 0-20\% centrality. The solid squares and stars are results for the $\phi$ meson, obtained with the 1st- and 2nd-order EP, respectively. The solid circles are results for $K^{*0}$-meson, obtained with the 2nd-order EP.

Numerical values of double-differential yields $d^{2}n/dydp_{T}$ presented in Fig. 10, given in units of $10^{−3} (GeV/c)^{−1}$. The first uncertainty is statistical, while the second one is systematic

Numerical values of double-differential yields $1/m_T d^{2}n/dm_{T}dy$ given in units of $10^{−3}(GeV)^{−2}$ and presented in Fig. 11; the values of $m_T−m_0$ specify the bin centers. The first uncertainty is statistical, while the second one is systematic

Numerical values of double-differential yields $1/m_T d^{2}n/dm_{T}dy$ given in units of $10^{−3}(GeV)^{−2}$ and presented in Fig. 11; the values of $m_T−m_0$ specify the bin centers. The first uncertainty is statistical, while the second one is systematic

Numerical values of double-differential yields $1/m_T d^{2}n/dm_{T}dy$ given in units of $10^{−3}(GeV)^{−2}$ and presented in Fig. 11; the values of $m_T−m_0$ specify the bin centers. The first uncertainty is statistical, while the second one is systematic

Numerical values of the $p_T$-integrated and extrapolated $dn/dy$ distribution presented in Fig. 12. The first uncertainty is statistical, while the second one is systematic. Additionally, the width of the Gaussian fit to the $dn/dy$ distribution, as well as the mean multiplicity of $K^{*}(892)^0$ mesons are shown (see the text for details)

The mean multiplicities $\langle K^{*}(892)^0\rangle$ and the widths of the rapidity distributions $\sigma_y$ obtained from $dn/dy$ distributions (see the text for details). The first uncertainty is statistical and the second systematic

The mean multiplicity of $K^{*}(892)^0$ mesons for 158 GeV/$c$ inelastic p+p interactions compared to theoretical multiplicities obtained within Hadron Gas Models

The mean multiplicities of different particle species measured in nucleus-nucleus collisions at 158$A$ GeV/$c$ by NA49 and NA61/SHINE. The total uncertainties of $\langle K^{*}(892)^0\rangle$, $\langle K^{+}\rangle$ and $\langle K^{-}\rangle$ were taken as the square roots of the sums of squares of statistical and systematic uncertainties. For NA49 p+p data, the $\langle K^{+}\rangle$ and $\langle K^{-}\rangle$ results include statistical uncertainties only($\langle K^{+}\rangle$ =0.2267$\pm$0.0006 and $\langle K^{-}\rangle$=0.1303$\pm$0.0004), whereas systematic uncertainties for total yields were not reported. Therefore, NA61/SHINE $\langle K^{+}\rangle$ and $\langle K^{-}\rangle$ values were used in the $\langle K^{*}(892)^0 \rangle / \langle K^{+} \rangle$ and $\langle K^{*}(892)^{0}\rangle / \langle K^{-} \rangle$ ratios. The numbers of $\langle K^{+} \rangle$ and $\langle K^{−}\rangle$ and their uncertainties for the 5% most central Pb+Pb collisions were multiplied by a factor 262/362 in orderto estimate charged kaon multiplicities in the 23.5% most central Pb+Pb reactions

<K+>/<PI>+.

Inverse slope parameter T.

Inverse slope parameter T.

Inverse slope parameter T.

Inverse slope parameter T.

K-/PI- at y=0.

K-/PI- at y=0.

<K->/<PI->.

<K->/<PI->.

Transverse mass $m_T$ spectrum at midrapidity of $\phi$ mesons produced in minimum bias p+p collisions at beam momentum of 158 GeV/c. $m_0$ is the PDG mass of the $\phi$ meson.

Transverse mass $m_T$ spectrum at midrapidity of $\phi$ mesons produced in minimum bias p+p collisions at beam momentum of 80 GeV/c. $m_0$ is the PDG mass of the $\phi$ meson.

Dependence of the slope parameter $T$ of the transverse momentum spectrum on rapidity $y$, of $\phi$ mesons produced in minimum bias p+p collisions at beam momentum of 158 GeV/c.

Dependence of the slope parameter $T$ of the transverse momentum spectrum on rapidity $y$, of $\phi$ mesons produced in minimum bias p+p collisions at beam momentum of 80 GeV/c.

Dependence of the width $\sigma_y$ of the rapidity distributions on $p_T$, of $\phi$ mesons produced in minimum bias p+p collisions at beam momentum of 158 GeV/c.

Dependence of the width $\sigma_y$ of the rapidity distributions on $p_T$, of $\phi$ mesons produced in minimum bias p+p collisions at beam momentum of 80 GeV/c.

Rapidity spectrum of $\phi$ mesons produced in minimum bias p+p collisions at beam momentum of 158 GeV/c.

Rapidity spectrum of $\phi$ mesons produced in minimum bias p+p collisions at beam momentum of 80 GeV/c.

Rapidity spectrum of $\phi$ mesons produced in minimum bias p+p collisions at beam momentum of 40 GeV/c.

Energy dependence of ratios of total yields of $\phi$ mesons to mean total yields of pions produced in in minimum bias p+p collisions at CERN SPS energies.

Energy dependence of double ratios of total yields of $\phi$ mesons to mean total yields of pions in central Pb+Pb collisions over minimum bias p+p collisions at CERN SPS energies.

Energy dependence of total yields of $\phi$ mesons produced in minimum bias p+p collisions at CERN SPS energies.

Energy dependence of midrapidity yields of $\phi$ mesons produced in minimum bias p+p collisions at CERN SPS energies.

Widths of rapidity distributions of $\phi$ mesons produced in minimum bias p+p collisions at CERN SPS energies, as a function of beam rapidity.

Raw $\Delta N_k$ distributions in Au+Au collisions at 19.6 GeV for 0–5%, 30–40%, and 70–80% collision centralities at midrapidity. The distributions are not corrected for the finite centrality bin width effect nor the reconstruction efficiency.

Raw $\Delta N_k$ distributions in Au+Au collisions at 27 GeV for 0–5%, 30–40%, and 70–80% collision centralities at midrapidity. The distributions are not corrected for the finite centrality bin width effect nor the reconstruction efficiency.

Raw $\Delta N_k$ distributions in Au+Au collisions at 39 GeV for 0–5%, 30–40%, and 70–80% collision centralities at midrapidity. The distributions are not corrected for the finite centrality bin width effect nor the reconstruction efficiency.

Raw $\Delta N_k$ distributions in Au+Au collisions at 62.4 GeV for 0–5%, 30–40%, and 70–80% collision centralities at midrapidity. The distributions are not corrected for the finite centrality bin width effect nor the reconstruction efficiency.

Raw $\Delta N_k$ distributions in Au+Au collisions at 200 GeV for 0–5%, 30–40%, and 70–80% collision centralities at midrapidity. The distributions are not corrected for the finite centrality bin width effect nor the reconstruction efficiency.

Collision centrality dependence of cumulants (C1, C2, C3, and C4) of $\Delta N_k$ distributions in Au+Au collisions at 7.7 GeV. The error bars are statistical uncertainties and the caps represent systematic uncertainties.

Collision centrality dependence of cumulants (C1, C2, C3, and C4) of $\Delta N_k$ distributions in Au+Au collisions at 11.5 GeV. The error bars are statistical uncertainties and the caps represent systematic uncertainties.

Collision centrality dependence of cumulants (C1, C2, C3, and C4) of $\Delta N_k$ distributions in Au+Au collisions at 14.5 GeV. The error bars are statistical uncertainties and the caps represent systematic uncertainties.

Collision centrality dependence of cumulants (C1, C2, C3, and C4) of $\Delta N_k$ distributions in Au+Au collisions at 19.6 GeV. The error bars are statistical uncertainties and the caps represent systematic uncertainties.

Collision centrality dependence of cumulants (C1, C2, C3, and C4) of $\Delta N_k$ distributions in Au+Au collisions at 27 GeV. The error bars are statistical uncertainties and the caps represent systematic uncertainties.

Collision centrality dependence of cumulants (C1, C2, C3, and C4) of $\Delta N_k$ distributions in Au+Au collisions at 39 GeV. The error bars are statistical uncertainties and the caps represent systematic uncertainties.

Collision centrality dependence of cumulants (C1, C2, C3, and C4) of $\Delta N_k$ distributions in Au+Au collisions at 62.4 GeV. The error bars are statistical uncertainties and the caps represent systematic uncertainties.

Collision centrality dependence of cumulants (C1, C2, C3, and C4) of $\Delta N_k$ distributions in Au+Au collisions at 200 GeV. The error bars are statistical uncertainties and the caps represent systematic uncertainties.

Collision centrality dependence of the $M/\sigma^2$ for $\Delta N_k$ distributions in Au+Au collisions at 7.7 GeV. The error bars are statistical uncertainties and the caps represent systematic uncertainties.

Collision centrality dependence of the $M/\sigma^2$ for $\Delta N_k$ distributions in Au+Au collisions at 11.5 GeV. The error bars are statistical uncertainties and the caps represent systematic uncertainties.

Collision centrality dependence of the $M/\sigma^2$ for $\Delta N_k$ distributions in Au+Au collisions at 14.5 GeV. The error bars are statistical uncertainties and the caps represent systematic uncertainties.

Collision centrality dependence of the $M/\sigma^2$ for $\Delta N_k$ distributions in Au+Au collisions at 19.6 GeV. The error bars are statistical uncertainties and the caps represent systematic uncertainties.

Collision centrality dependence of the $M/\sigma^2$ for $\Delta N_k$ distributions in Au+Au collisions at 27 GeV. The error bars are statistical uncertainties and the caps represent systematic uncertainties.

Collision centrality dependence of the $M/\sigma^2$ for $\Delta N_k$ distributions in Au+Au collisions at 39 GeV. The error bars are statistical uncertainties and the caps represent systematic uncertainties.

Collision centrality dependence of the $M/\sigma^2$ for $\Delta N_k$ distributions in Au+Au collisions at 62.4 GeV. The error bars are statistical uncertainties and the caps represent systematic uncertainties.

Collision centrality dependence of the $M/\sigma^2$ for $\Delta N_k$ distributions in Au+Au collisions at 200 GeV. The error bars are statistical uncertainties and the caps represent systematic uncertainties.

Collision centrality dependence of the $S\sigma$ for $\Delta N_k$ distributions in Au+Au collisions at 7.7 GeV. The error bars are statistical uncertainties and the caps represent systematic uncertainties.

Collision centrality dependence of the $S\sigma$ for $\Delta N_k$ distributions in Au+Au collisions at 11.5 GeV. The error bars are statistical uncertainties and the caps represent systematic uncertainties.

Collision centrality dependence of the $S\sigma$ for $\Delta N_k$ distributions in Au+Au collisions at 14.5 GeV. The error bars are statistical uncertainties and the caps represent systematic uncertainties.

Collision centrality dependence of the $S\sigma$ for $\Delta N_k$ distributions in Au+Au collisions at 19.6 GeV. The error bars are statistical uncertainties and the caps represent systematic uncertainties.

Collision centrality dependence of the $S\sigma$ for $\Delta N_k$ distributions in Au+Au collisions at 27 GeV. The error bars are statistical uncertainties and the caps represent systematic uncertainties.

Collision centrality dependence of the $S\sigma$ for $\Delta N_k$ distributions in Au+Au collisions at 39 GeV. The error bars are statistical uncertainties and the caps represent systematic uncertainties.

Collision centrality dependence of the $S\sigma$ for $\Delta N_k$ distributions in Au+Au collisions at 62.4 GeV. The error bars are statistical uncertainties and the caps represent systematic uncertainties.

Collision centrality dependence of the $S\sigma$ for $\Delta N_k$ distributions in Au+Au collisions at 200 GeV. The error bars are statistical uncertainties and the caps represent systematic uncertainties.

Collision centrality dependence of the $\kappa\sigma^2$ for $\Delta N_k$ distributions in Au+Au collisions at 7.7 GeV. The error bars are statistical uncertainties and the caps represent systematic uncertainties.

Collision centrality dependence of the $\kappa\sigma^2$ for $\Delta N_k$ distributions in Au+Au collisions at 11.5 GeV. The error bars are statistical uncertainties and the caps represent systematic uncertainties.

Collision centrality dependence of the $\kappa\sigma^2$ for $\Delta N_k$ distributions in Au+Au collisions at 14.5 GeV. The error bars are statistical uncertainties and the caps represent systematic uncertainties.

Collision centrality dependence of the $\kappa\sigma^2$ for $\Delta N_k$ distributions in Au+Au collisions at 19.6 GeV. The error bars are statistical uncertainties and the caps represent systematic uncertainties.

Collision centrality dependence of the $\kappa\sigma^2$ for $\Delta N_k$ distributions in Au+Au collisions at 27 GeV. The error bars are statistical uncertainties and the caps represent systematic uncertainties.

Collision centrality dependence of the $\kappa\sigma^2$ for $\Delta N_k$ distributions in Au+Au collisions at 39 GeV. The error bars are statistical uncertainties and the caps represent systematic uncertainties.

Collision centrality dependence of the $\kappa\sigma^2$ for $\Delta N_k$ distributions in Au+Au collisions at 62.4 GeV. The error bars are statistical uncertainties and the caps represent systematic uncertainties.

Collision centrality dependence of the $\kappa\sigma^2$ for $\Delta N_k$ distributions in Au+Au collisions at 200 GeV. The error bars are statistical uncertainties and the caps represent systematic uncertainties.

Collisions energy dependence of $M/\sigma^2$ for $\Delta N_k$ multiplicity distributions from 0–5% most central and 70–80% peripheral collisions in Au+Au collisions at \sqrt{s_{NN}} = 7.7, 11.5, 14.5, 19.6, 27, 39, 62.4 and 200 GeV. The error bars are statistical uncertainties and the caps represent systematic uncertainties.

Collisions energy dependence of $S\sigma$ for $\Delta N_k$ multiplicity distributions from 0–5% most central and 70–80% peripheral collisions in Au+Au collisions at \sqrt{s_{NN}} = 7.7, 11.5, 14.5, 19.6, 27, 39, 62.4 and 200 GeV. The error bars are statistical uncertainties and the caps represent systematic uncertainties.

Collisions energy dependence of $\kappa\sigma^2$ for $\Delta N_k$ multiplicity distributions from 0–5% most central and 70–80% peripheral collisions in Au+Au collisions at \sqrt{s_{NN}} = 7.7, 11.5, 14.5, 19.6, 27, 39, 62.4 and 200 GeV. The error bars are statistical uncertainties and the caps represent systematic uncertainties.

Identified particle (Kaon Plus) RCP for RHIC BES energies. The colored shaded boxes describe the point-to-point systematic uncertainties. The uncertainty bands at unity on the right side of the plot correspond to the pT -independent uncertainty in Ncoll scaling with the color in the band corresponding to the color of the data points for that energy.

Identified particle (Kaon Minus) RCP for RHIC BES energies. The colored shaded boxes describe the point-to-point systematic uncertainties. The uncertainty bands at unity on the right side of the plot correspond to the pT -independent uncertainty in Ncoll scaling with the color in the band corresponding to the color of the data points for that energy.

Identified particle (Proton) RCP for RHIC BES energies. The colored shaded boxes describe the point-to-point systematic uncertainties. The uncertainty bands at unity on the right side of the plot correspond to the pT -independent uncertainty in Ncoll scaling with the color in the band corresponding to the color of the data points for that energy.

Identified particle (Antiproton) RCP for RHIC BES energies. The colored shaded boxes describe the point-to-point systematic uncertainties. The uncertainty bands at unity on the right side of the plot correspond to the pT -independent uncertainty in Ncoll scaling with the color in the band corresponding to the color of the data points for that energy.

Charged hadron Y(<Npart>) for two ranges of pT (pT 3.0 - 3.5 GeV/c). Statistical uncertainty bars are included, mostly smaller than point size, as well as shaded bands to indicate systematic uncertainties.

Charged hadron Y(<Npart>) for two ranges of pT (pT 4.0 - 4.5 GeV/c). Statistical uncertainty bars are included, mostly smaller than point size, as well as shaded bands to indicate systematic uncertainties.

Glauber Fit Parameters

Nch at each Collision Energy (GeV)

Ncoll at each Collision Energy (GeV)

Npart at each Collision Energy (GeV)

The value of $\sigma^{NN}_{inel}$ used in the Monte Carlo Glauber simulation at each collision energy

Charged hadron $\frac{1}{2\pi p_{T}}$ * $\frac{d^{2}N}{d\eta dp_{T}}$ $\pm$ stat. $\pm$ sys. $(GeV/c)^{-2}$ for $\sqrt{s_{NN}}$ = 7.7 GeV/c

Charged hadron $\frac{1}{2\pi p_{T}}$ * $\frac{d^{2}N}{d\eta dp_{T}}$ $\pm$ stat. $\pm$ sys. $(GeV/c)^{-2}$ for $\sqrt{s_{NN}}$ = 7.7 GeV/c

Charged hadron $\frac{1}{2\pi p_{T}}$ * $\frac{d^{2}N}{d\eta dp_{T}}$ $\pm$ stat. $\pm$ sys. $(GeV/c)^{-2}$ for $\sqrt{s_{NN}}$ = 11.5 GeV/c

Charged hadron $\frac{1}{2\pi p_{T}}$ * $\frac{d^{2}N}{d\eta dp_{T}}$ $\pm$ stat. $\pm$ sys. $(GeV/c)^{-2}$ for $\sqrt{s_{NN}}$ = 11.5 GeV/c

Charged hadron $\frac{1}{2\pi p_{T}}$ * $\frac{d^{2}N}{d\eta dp_{T}}$ $\pm$ stat. $\pm$ sys. $(GeV/c)^{-2}$ for $\sqrt{s_{NN}}$ = 14.5 GeV/c

Charged hadron $\frac{1}{2\pi p_{T}}$ * $\frac{d^{2}N}{d\eta dp_{T}}$ $\pm$ stat. $\pm$ sys. $(GeV/c)^{-2}$ for $\sqrt{s_{NN}}$ = 14.5 GeV/c

Charged hadron $\frac{1}{2\pi p_{T}}$ * $\frac{d^{2}N}{d\eta dp_{T}}$ $\pm$ stat. $\pm$ sys. $(GeV/c)^{-2}$ for $\sqrt{s_{NN}}$ = 19.6 GeV/c

Charged hadron $\frac{1}{2\pi p_{T}}$ * $\frac{d^{2}N}{d\eta dp_{T}}$ $\pm$ stat. $\pm$ sys. $(GeV/c)^{-2}$ for $\sqrt{s_{NN}}$ = 19.6 GeV/c

Charged hadron $\frac{1}{2\pi p_{T}}$ * $\frac{d^{2}N}{d\eta dp_{T}}$ $\pm$ stat. $\pm$ sys. $(GeV/c)^{-2}$ for $\sqrt{s_{NN}}$ = 27 GeV/c

Charged hadron $\frac{1}{2\pi p_{T}}$ * $\frac{d^{2}N}{d\eta dp_{T}}$ $\pm$ stat. $\pm$ sys. $(GeV/c)^{-2}$ for $\sqrt{s_{NN}}$ = 27 GeV/c

Charged hadron $\frac{1}{2\pi p_{T}}$ * $\frac{d^{2}N}{d\eta dp_{T}}$ $\pm$ stat. $\pm$ sys. $(GeV/c)^{-2}$ for $\sqrt{s_{NN}}$ = 39 GeV/c

Charged hadron $\frac{1}{2\pi p_{T}}$ * $\frac{d^{2}N}{d\eta dp_{T}}$ $\pm$ stat. $\pm$ sys. $(GeV/c)^{-2}$ for $\sqrt{s_{NN}}$ = 39 GeV/c

Charged hadron $\frac{1}{2\pi p_{T}}$ * $\frac{d^{2}N}{d\eta dp_{T}}$ $\pm$ stat. $\pm$ sys. $(GeV/c)^{-2}$ for $\sqrt{s_{NN}}$ = 62.4 GeV/c

Charged hadron $\frac{1}{2\pi p_{T}}$ * $\frac{d^{2}N}{d\eta dp_{T}}$ $\pm$ stat. $\pm$ sys. $(GeV/c)^{-2}$ for $\sqrt{s_{NN}}$ = 62.4 GeV/c

$\\p$ $\frac{1}{2\pi p_{T}}$ * $\frac{d^{2}N}{d\eta dp_{T}}$ $\pm$ stat. $\pm$ sys. $(GeV/c)^{-2}$ for $\sqrt{s_{NN}}$ = 7.7 GeV/c

$\\p$ $\frac{1}{2\pi p_{T}}$ * $\frac{d^{2}N}{d\eta dp_{T}}$ $\pm$ stat. $\pm$ sys. $(GeV/c)^{-2}$ for $\sqrt{s_{NN}}$ = 7.7 GeV/c

$\\p$ $\frac{1}{2\pi p_{T}}$ * $\frac{d^{2}N}{d\eta dp_{T}}$ $\pm$ stat. $\pm$ sys. $(GeV/c)^{-2}$ for $\sqrt{s_{NN}}$ = 11.5 GeV/c

$\\p$ $\frac{1}{2\pi p_{T}}$ * $\frac{d^{2}N}{d\eta dp_{T}}$ $\pm$ stat. $\pm$ sys. $(GeV/c)^{-2}$ for $\sqrt{s_{NN}}$ = 11.5 GeV/c

$\\p$ $\frac{1}{2\pi p_{T}}$ * $\frac{d^{2}N}{d\eta dp_{T}}$ $\pm$ stat. $\pm$ sys. $(GeV/c)^{-2}$ for $\sqrt{s_{NN}}$ = 14.5 GeV/c

$\\p$ $\frac{1}{2\pi p_{T}}$ * $\frac{d^{2}N}{d\eta dp_{T}}$ $\pm$ stat. $\pm$ sys. $(GeV/c)^{-2}$ for $\sqrt{s_{NN}}$ = 14.5 GeV/c

$\\p$ $\frac{1}{2\pi p_{T}}$ * $\frac{d^{2}N}{d\eta dp_{T}}$ $\pm$ stat. $\pm$ sys. $(GeV/c)^{-2}$ for $\sqrt{s_{NN}}$ = 14.5 GeV/c

$\\p$ $\frac{1}{2\pi p_{T}}$ * $\frac{d^{2}N}{d\eta dp_{T}}$ $\pm$ stat. $\pm$ sys. $(GeV/c)^{-2}$ for $\sqrt{s_{NN}}$ = 19.6 GeV/c

$\\p$ $\frac{1}{2\pi p_{T}}$ * $\frac{d^{2}N}{d\eta dp_{T}}$ $\pm$ stat. $\pm$ sys. $(GeV/c)^{-2}$ for $\sqrt{s_{NN}}$ = 19.6 GeV/c

$\\p$ $\frac{1}{2\pi p_{T}}$ * $\frac{d^{2}N}{d\eta dp_{T}}$ $\pm$ stat. $\pm$ sys. $(GeV/c)^{-2}$ for $\sqrt{s_{NN}}$ = 27 GeV/c

$\\p$ $\frac{1}{2\pi p_{T}}$ * $\frac{d^{2}N}{d\eta dp_{T}}$ $\pm$ stat. $\pm$ sys. $(GeV/c)^{-2}$ for $\sqrt{s_{NN}}$ = 27 GeV/c

$\\p$ $\frac{1}{2\pi p_{T}}$ * $\frac{d^{2}N}{d\eta dp_{T}}$ $\pm$ stat. $\pm$ sys. $(GeV/c)^{-2}$ for $\sqrt{s_{NN}}$ = 39 GeV/c

$\\p$ $\frac{1}{2\pi p_{T}}$ * $\frac{d^{2}N}{d\eta dp_{T}}$ $\pm$ stat. $\pm$ sys. $(GeV/c)^{-2}$ for $\sqrt{s_{NN}}$ = 39 GeV/c

$\\p$ $\frac{1}{2\pi p_{T}}$ * $\frac{d^{2}N}{d\eta dp_{T}}$ $\pm$ stat. $\pm$ sys. $(GeV/c)^{-2}$ for $\sqrt{s_{NN}}$ = 62.4 GeV/c

$\\p$ $\frac{1}{2\pi p_{T}}$ * $\frac{d^{2}N}{d\eta dp_{T}}$ $\pm$ stat. $\pm$ sys. $(GeV/c)^{-2}$ for $\sqrt{s_{NN}}$ = 62.4 GeV/c

$\overline{p}$ $\frac{1}{2\pi p_{T}}$ * $\frac{d^{2}N}{d\eta dp_{T}}$ $\pm$ stat. $\pm$ sys. $(GeV/c)^{-2}$ for $\sqrt{s_{NN}}$ = 7.7 GeV/c

$\overline{p}$ $\frac{1}{2\pi p_{T}}$ * $\frac{d^{2}N}{d\eta dp_{T}}$ $\pm$ stat. $\pm$ sys. $(GeV/c)^{-2}$ for $\sqrt{s_{NN}}$ = 7.7 GeV/c

$\overline{p}$ $\frac{1}{2\pi p_{T}}$ * $\frac{d^{2}N}{d\eta dp_{T}}$ $\pm$ stat. $\pm$ sys. $(GeV/c)^{-2}$ for $\sqrt{s_{NN}}$ = 11.5 GeV/c

$\overline{p}$ $\frac{1}{2\pi p_{T}}$ * $\frac{d^{2}N}{d\eta dp_{T}}$ $\pm$ stat. $\pm$ sys. $(GeV/c)^{-2}$ for $\sqrt{s_{NN}}$ = 11.5 GeV/c

$\overline{p}$ $\frac{1}{2\pi p_{T}}$ * $\frac{d^{2}N}{d\eta dp_{T}}$ $\pm$ stat. $\pm$ sys. $(GeV/c)^{-2}$ for $\sqrt{s_{NN}}$ = 14.5 GeV/c

$\overline{p}$ $\frac{1}{2\pi p_{T}}$ * $\frac{d^{2}N}{d\eta dp_{T}}$ $\pm$ stat. $\pm$ sys. $(GeV/c)^{-2}$ for $\sqrt{s_{NN}}$ = 14.5 GeV/c

$\overline{p}$ $\frac{1}{2\pi p_{T}}$ * $\frac{d^{2}N}{d\eta dp_{T}}$ $\pm$ stat. $\pm$ sys. $(GeV/c)^{-2}$ for $\sqrt{s_{NN}}$ = 14.5 GeV/c

$\overline{p}$ $\frac{1}{2\pi p_{T}}$ * $\frac{d^{2}N}{d\eta dp_{T}}$ $\pm$ stat. $\pm$ sys. $(GeV/c)^{-2}$ for $\sqrt{s_{NN}}$ = 19.6 GeV/c

$\overline{p}$ $\frac{1}{2\pi p_{T}}$ * $\frac{d^{2}N}{d\eta dp_{T}}$ $\pm$ stat. $\pm$ sys. $(GeV/c)^{-2}$ for $\sqrt{s_{NN}}$ = 19.6 GeV/c

$\overline{p}$ $\frac{1}{2\pi p_{T}}$ * $\frac{d^{2}N}{d\eta dp_{T}}$ $\pm$ stat. $\pm$ sys. $(GeV/c)^{-2}$ for $\sqrt{s_{NN}}$ = 27 GeV/c

$\overline{p}$ $\frac{1}{2\pi p_{T}}$ * $\frac{d^{2}N}{d\eta dp_{T}}$ $\pm$ stat. $\pm$ sys. $(GeV/c)^{-2}$ for $\sqrt{s_{NN}}$ = 27 GeV/c

$\overline{p}$ $\frac{1}{2\pi p_{T}}$ * $\frac{d^{2}N}{d\eta dp_{T}}$ $\pm$ stat. $\pm$ sys. $(GeV/c)^{-2}$ for $\sqrt{s_{NN}}$ = 39 GeV/c

$\overline{p}$ $\frac{1}{2\pi p_{T}}$ * $\frac{d^{2}N}{d\eta dp_{T}}$ $\pm$ stat. $\pm$ sys. $(GeV/c)^{-2}$ for $\sqrt{s_{NN}}$ = 39 GeV/c

$\overline{p}$ $\frac{1}{2\pi p_{T}}$ * $\frac{d^{2}N}{d\eta dp_{T}}$ $\pm$ stat. $\pm$ sys. $(GeV/c)^{-2}$ for $\sqrt{s_{NN}}$ = 62.4 GeV/c

$\overline{p}$ $\frac{1}{2\pi p_{T}}$ * $\frac{d^{2}N}{d\eta dp_{T}}$ $\pm$ stat. $\pm$ sys. $(GeV/c)^{-2}$ for $\sqrt{s_{NN}}$ = 62.4 GeV/c

$K^{+}$ $\frac{1}{2\pi p_{T}}$ * $\frac{d^{2}N}{d\eta dp_{T}}$ $\pm$ stat. $\pm$ sys. $(GeV/c)^{-2}$ for $\sqrt{s_{NN}}$ = 7.7 GeV/c

$K^{+}$ $\frac{1}{2\pi p_{T}}$ * $\frac{d^{2}N}{d\eta dp_{T}}$ $\pm$ stat. $\pm$ sys. $(GeV/c)^{-2}$ for $\sqrt{s_{NN}}$ = 7.7 GeV/c

$K^{+}$ $\frac{1}{2\pi p_{T}}$ * $\frac{d^{2}N}{d\eta dp_{T}}$ $\pm$ stat. $\pm$ sys. $(GeV/c)^{-2}$ for $\sqrt{s_{NN}}$ = 11.5 GeV/c

$K^{+}$ $\frac{1}{2\pi p_{T}}$ * $\frac{d^{2}N}{d\eta dp_{T}}$ $\pm$ stat. $\pm$ sys. $(GeV/c)^{-2}$ for $\sqrt{s_{NN}}$ = 11.5 GeV/c

$K^{+}$ $\frac{1}{2\pi p_{T}}$ * $\frac{d^{2}N}{d\eta dp_{T}}$ $\pm$ stat. $\pm$ sys. $(GeV/c)^{-2}$ for $\sqrt{s_{NN}}$ = 14.5 GeV/c

$K^{+}$ $\frac{1}{2\pi p_{T}}$ * $\frac{d^{2}N}{d\eta dp_{T}}$ $\pm$ stat. $\pm$ sys. $(GeV/c)^{-2}$ for $\sqrt{s_{NN}}$ = 14.5 GeV/c

$K^{+}$ $\frac{1}{2\pi p_{T}}$ * $\frac{d^{2}N}{d\eta dp_{T}}$ $\pm$ stat. $\pm$ sys. $(GeV/c)^{-2}$ for $\sqrt{s_{NN}}$ = 14.5 GeV/c

$K^{+}$ $\frac{1}{2\pi p_{T}}$ * $\frac{d^{2}N}{d\eta dp_{T}}$ $\pm$ stat. $\pm$ sys. $(GeV/c)^{-2}$ for $\sqrt{s_{NN}}$ = 19.6 GeV/c

$K^{+}$ $\frac{1}{2\pi p_{T}}$ * $\frac{d^{2}N}{d\eta dp_{T}}$ $\pm$ stat. $\pm$ sys. $(GeV/c)^{-2}$ for $\sqrt{s_{NN}}$ = 19.6 GeV/c

$K^{+}$ $\frac{1}{2\pi p_{T}}$ * $\frac{d^{2}N}{d\eta dp_{T}}$ $\pm$ stat. $\pm$ sys. $(GeV/c)^{-2}$ for $\sqrt{s_{NN}}$ = 27 GeV/c

$K^{+}$ $\frac{1}{2\pi p_{T}}$ * $\frac{d^{2}N}{d\eta dp_{T}}$ $\pm$ stat. $\pm$ sys. $(GeV/c)^{-2}$ for $\sqrt{s_{NN}}$ = 27 GeV/c

$K^{+}$ $\frac{1}{2\pi p_{T}}$ * $\frac{d^{2}N}{d\eta dp_{T}}$ $\pm$ stat. $\pm$ sys. $(GeV/c)^{-2}$ for $\sqrt{s_{NN}}$ = 39 GeV/c

$K^{+}$ $\frac{1}{2\pi p_{T}}$ * $\frac{d^{2}N}{d\eta dp_{T}}$ $\pm$ stat. $\pm$ sys. $(GeV/c)^{-2}$ for $\sqrt{s_{NN}}$ = 39 GeV/c

$K^{+}$ $\frac{1}{2\pi p_{T}}$ * $\frac{d^{2}N}{d\eta dp_{T}}$ $\pm$ stat. $\pm$ sys. $(GeV/c)^{-2}$ for $\sqrt{s_{NN}}$ = 62.4 GeV/c

$K^{+}$ $\frac{1}{2\pi p_{T}}$ * $\frac{d^{2}N}{d\eta dp_{T}}$ $\pm$ stat. $\pm$ sys. $(GeV/c)^{-2}$ for $\sqrt{s_{NN}}$ = 62.4 GeV/c

$K^{-}$ $\frac{1}{2\pi p_{T}}$ * $\frac{d^{2}N}{d\eta dp_{T}}$ $\pm$ stat. $\pm$ sys. $(GeV/c)^{-2}$ for $\sqrt{s_{NN}}$ = 7.7 GeV/c

$K^{-}$ $\frac{1}{2\pi p_{T}}$ * $\frac{d^{2}N}{d\eta dp_{T}}$ $\pm$ stat. $\pm$ sys. $(GeV/c)^{-2}$ for $\sqrt{s_{NN}}$ = 7.7 GeV/c

$K^{-}$ $\frac{1}{2\pi p_{T}}$ * $\frac{d^{2}N}{d\eta dp_{T}}$ $\pm$ stat. $\pm$ sys. $(GeV/c)^{-2}$ for $\sqrt{s_{NN}}$ = 11.5 GeV/c

$K^{-}$ $\frac{1}{2\pi p_{T}}$ * $\frac{d^{2}N}{d\eta dp_{T}}$ $\pm$ stat. $\pm$ sys. $(GeV/c)^{-2}$ for $\sqrt{s_{NN}}$ = 11.5 GeV/c

$K^{-}$ $\frac{1}{2\pi p_{T}}$ * $\frac{d^{2}N}{d\eta dp_{T}}$ $\pm$ stat. $\pm$ sys. $(GeV/c)^{-2}$ for $\sqrt{s_{NN}}$ = 14.5 GeV/c

$K^{-}$ $\frac{1}{2\pi p_{T}}$ * $\frac{d^{2}N}{d\eta dp_{T}}$ $\pm$ stat. $\pm$ sys. $(GeV/c)^{-2}$ for $\sqrt{s_{NN}}$ = 14.5 GeV/c

$K^{-}$ $\frac{1}{2\pi p_{T}}$ * $\frac{d^{2}N}{d\eta dp_{T}}$ $\pm$ stat. $\pm$ sys. $(GeV/c)^{-2}$ for $\sqrt{s_{NN}}$ = 14.5 GeV/c

$K^{-}$ $\frac{1}{2\pi p_{T}}$ * $\frac{d^{2}N}{d\eta dp_{T}}$ $\pm$ stat. $\pm$ sys. $(GeV/c)^{-2}$ for $\sqrt{s_{NN}}$ = 19.6 GeV/c

$K^{-}$ $\frac{1}{2\pi p_{T}}$ * $\frac{d^{2}N}{d\eta dp_{T}}$ $\pm$ stat. $\pm$ sys. $(GeV/c)^{-2}$ for $\sqrt{s_{NN}}$ = 19.6 GeV/c

$K^{-}$ $\frac{1}{2\pi p_{T}}$ * $\frac{d^{2}N}{d\eta dp_{T}}$ $\pm$ stat. $\pm$ sys. $(GeV/c)^{-2}$ for $\sqrt{s_{NN}}$ = 27 GeV/c

$K^{-}$ $\frac{1}{2\pi p_{T}}$ * $\frac{d^{2}N}{d\eta dp_{T}}$ $\pm$ stat. $\pm$ sys. $(GeV/c)^{-2}$ for $\sqrt{s_{NN}}$ = 27 GeV/c

$K^{-}$ $\frac{1}{2\pi p_{T}}$ * $\frac{d^{2}N}{d\eta dp_{T}}$ $\pm$ stat. $\pm$ sys. $(GeV/c)^{-2}$ for $\sqrt{s_{NN}}$ = 39 GeV/c

$K^{-}$ $\frac{1}{2\pi p_{T}}$ * $\frac{d^{2}N}{d\eta dp_{T}}$ $\pm$ stat. $\pm$ sys. $(GeV/c)^{-2}$ for $\sqrt{s_{NN}}$ = 39 GeV/c

$K^{-}$ $\frac{1}{2\pi p_{T}}$ * $\frac{d^{2}N}{d\eta dp_{T}}$ $\pm$ stat. $\pm$ sys. $(GeV/c)^{-2}$ for $\sqrt{s_{NN}}$ = 62.4 GeV/c

$K^{-}$ $\frac{1}{2\pi p_{T}}$ * $\frac{d^{2}N}{d\eta dp_{T}}$ $\pm$ stat. $\pm$ sys. $(GeV/c)^{-2}$ for $\sqrt{s_{NN}}$ = 62.4 GeV/c

$\pi^{+}$ $\frac{1}{2\pi p_{T}}$ * $\frac{d^{2}N}{d\eta dp_{T}}$ $\pm$ stat. $\pm$ sys. $(GeV/c)^{-2}$ for $\sqrt{s_{NN}}$ = 7.7 GeV/c

$\pi^{+}$ $\frac{1}{2\pi p_{T}}$ * $\frac{d^{2}N}{d\eta dp_{T}}$ $\pm$ stat. $\pm$ sys. $(GeV/c)^{-2}$ for $\sqrt{s_{NN}}$ = 7.7 GeV/c

$\pi^{+}$ $\frac{1}{2\pi p_{T}}$ * $\frac{d^{2}N}{d\eta dp_{T}}$ $\pm$ stat. $\pm$ sys. $(GeV/c)^{-2}$ for $\sqrt{s_{NN}}$ = 11.5 GeV/c

$\pi^{+}$ $\frac{1}{2\pi p_{T}}$ * $\frac{d^{2}N}{d\eta dp_{T}}$ $\pm$ stat. $\pm$ sys. $(GeV/c)^{-2}$ for $\sqrt{s_{NN}}$ = 11.5 GeV/c

$\pi^{+}$ $\frac{1}{2\pi p_{T}}$ * $\frac{d^{2}N}{d\eta dp_{T}}$ $\pm$ stat. $\pm$ sys. $(GeV/c)^{-2}$ for $\sqrt{s_{NN}}$ = 14.5 GeV/c

$\pi^{+}$ $\frac{1}{2\pi p_{T}}$ * $\frac{d^{2}N}{d\eta dp_{T}}$ $\pm$ stat. $\pm$ sys. $(GeV/c)^{-2}$ for $\sqrt{s_{NN}}$ = 14.5 GeV/c

$\pi^{+}$ $\frac{1}{2\pi p_{T}}$ * $\frac{d^{2}N}{d\eta dp_{T}}$ $\pm$ stat. $\pm$ sys. $(GeV/c)^{-2}$ for $\sqrt{s_{NN}}$ = 14.5 GeV/c

$\pi^{+}$ $\frac{1}{2\pi p_{T}}$ * $\frac{d^{2}N}{d\eta dp_{T}}$ $\pm$ stat. $\pm$ sys. $(GeV/c)^{-2}$ for $\sqrt{s_{NN}}$ = 19.6 GeV/c

$\pi^{+}$ $\frac{1}{2\pi p_{T}}$ * $\frac{d^{2}N}{d\eta dp_{T}}$ $\pm$ stat. $\pm$ sys. $(GeV/c)^{-2}$ for $\sqrt{s_{NN}}$ = 19.6 GeV/c

$\pi^{+}$ $\frac{1}{2\pi p_{T}}$ * $\frac{d^{2}N}{d\eta dp_{T}}$ $\pm$ stat. $\pm$ sys. $(GeV/c)^{-2}$ for $\sqrt{s_{NN}}$ = 27 GeV/c

$\pi^{+}$ $\frac{1}{2\pi p_{T}}$ * $\frac{d^{2}N}{d\eta dp_{T}}$ $\pm$ stat. $\pm$ sys. $(GeV/c)^{-2}$ for $\sqrt{s_{NN}}$ = 27 GeV/c

$\pi^{+}$ $\frac{1}{2\pi p_{T}}$ * $\frac{d^{2}N}{d\eta dp_{T}}$ $\pm$ stat. $\pm$ sys. $(GeV/c)^{-2}$ for $\sqrt{s_{NN}}$ = 39 GeV/c

$\pi^{+}$ $\frac{1}{2\pi p_{T}}$ * $\frac{d^{2}N}{d\eta dp_{T}}$ $\pm$ stat. $\pm$ sys. $(GeV/c)^{-2}$ for $\sqrt{s_{NN}}$ = 39 GeV/c

$\pi^{+}$ $\frac{1}{2\pi p_{T}}$ * $\frac{d^{2}N}{d\eta dp_{T}}$ $\pm$ stat. $\pm$ sys. $(GeV/c)^{-2}$ for $\sqrt{s_{NN}}$ = 62.4 GeV/c

$\pi^{+}$ $\frac{1}{2\pi p_{T}}$ * $\frac{d^{2}N}{d\eta dp_{T}}$ $\pm$ stat. $\pm$ sys. $(GeV/c)^{-2}$ for $\sqrt{s_{NN}}$ = 62.4 GeV/c

$\pi^{-}$ $\frac{1}{2\pi p_{T}}$ * $\frac{d^{2}N}{d\eta dp_{T}}$ $\pm$ stat. $\pm$ sys. $(GeV/c)^{-2}$ for $\sqrt{s_{NN}}$ = 7.7 GeV/c

$\pi^{-}$ $\frac{1}{2\pi p_{T}}$ * $\frac{d^{2}N}{d\eta dp_{T}}$ $\pm$ stat. $\pm$ sys. $(GeV/c)^{-2}$ for $\sqrt{s_{NN}}$ = 7.7 GeV/c

$\pi^{-}$ $\frac{1}{2\pi p_{T}}$ * $\frac{d^{2}N}{d\eta dp_{T}}$ $\pm$ stat. $\pm$ sys. $(GeV/c)^{-2}$ for $\sqrt{s_{NN}}$ = 11.5 GeV/c

$\pi^{-}$ $\frac{1}{2\pi p_{T}}$ * $\frac{d^{2}N}{d\eta dp_{T}}$ $\pm$ stat. $\pm$ sys. $(GeV/c)^{-2}$ for $\sqrt{s_{NN}}$ = 11.5 GeV/c

$\pi^{-}$ $\frac{1}{2\pi p_{T}}$ * $\frac{d^{2}N}{d\eta dp_{T}}$ $\pm$ stat. $\pm$ sys. $(GeV/c)^{-2}$ for $\sqrt{s_{NN}}$ = 14.5 GeV/c

$\pi^{-}$ $\frac{1}{2\pi p_{T}}$ * $\frac{d^{2}N}{d\eta dp_{T}}$ $\pm$ stat. $\pm$ sys. $(GeV/c)^{-2}$ for $\sqrt{s_{NN}}$ = 14.5 GeV/c

$\pi^{-}$ $\frac{1}{2\pi p_{T}}$ * $\frac{d^{2}N}{d\eta dp_{T}}$ $\pm$ stat. $\pm$ sys. $(GeV/c)^{-2}$ for $\sqrt{s_{NN}}$ = 14.5 GeV/c

$\pi^{-}$ $\frac{1}{2\pi p_{T}}$ * $\frac{d^{2}N}{d\eta dp_{T}}$ $\pm$ stat. $\pm$ sys. $(GeV/c)^{-2}$ for $\sqrt{s_{NN}}$ = 19.6 GeV/c

$\pi^{-}$ $\frac{1}{2\pi p_{T}}$ * $\frac{d^{2}N}{d\eta dp_{T}}$ $\pm$ stat. $\pm$ sys. $(GeV/c)^{-2}$ for $\sqrt{s_{NN}}$ = 19.6 GeV/c

$\pi^{-}$ $\frac{1}{2\pi p_{T}}$ * $\frac{d^{2}N}{d\eta dp_{T}}$ $\pm$ stat. $\pm$ sys. $(GeV/c)^{-2}$ for $\sqrt{s_{NN}}$ = 27 GeV/c

$\pi^{-}$ $\frac{1}{2\pi p_{T}}$ * $\frac{d^{2}N}{d\eta dp_{T}}$ $\pm$ stat. $\pm$ sys. $(GeV/c)^{-2}$ for $\sqrt{s_{NN}}$ = 27 GeV/c

$\pi^{-}$ $\frac{1}{2\pi p_{T}}$ * $\frac{d^{2}N}{d\eta dp_{T}}$ $\pm$ stat. $\pm$ sys. $(GeV/c)^{-2}$ for $\sqrt{s_{NN}}$ = 39 GeV/c

$\pi^{-}$ $\frac{1}{2\pi p_{T}}$ * $\frac{d^{2}N}{d\eta dp_{T}}$ $\pm$ stat. $\pm$ sys. $(GeV/c)^{-2}$ for $\sqrt{s_{NN}}$ = 39 GeV/c

$\pi^{-}$ $\frac{1}{2\pi p_{T}}$ * $\frac{d^{2}N}{d\eta dp_{T}}$ $\pm$ stat. $\pm$ sys. $(GeV/c)^{-2}$ for $\sqrt{s_{NN}}$ = 62.4 GeV/c

$\pi^{-}$ $\frac{1}{2\pi p_{T}}$ * $\frac{d^{2}N}{d\eta dp_{T}}$ $\pm$ stat. $\pm$ sys. $(GeV/c)^{-2}$ for $\sqrt{s_{NN}}$ = 62.4 GeV/c

Transverse momentum-rapidity spectrum of K+ produced in inelastic p+p interactions at 20 GeV/c with systematic uncertainties.

Transverse momentum-rapidity spectrum of p produced in inelastic p+p interactions at 20 GeV/c with statistical uncertainties.

Transverse momentum-rapidity spectrum of p produced in inelastic p+p interactions at 20 GeV/c with systematic uncertainties.

Transverse momentum-rapidity spectrum of pi+ produced in inelastic p+p interactions at 20 GeV/c with statistical uncertainties.

Transverse momentum-rapidity spectrum of pi+ produced in inelastic p+p interactions at 20 GeV/c with ssystematic uncertainties.

Transverse momentum-rapidity spectrum of pi− produced in inelastic p+p interactions at 20 GeV/c with statistical uncertainties.

Transverse momentum-rapidity spectrum of pi− produced in inelastic p+p interactions at 20 GeV/c with systematic uncertainties.

Transverse momentum-rapidity spectrum of K− produced in inelastic p+p interactions at 31 GeV/c with statistical uncertainties.

Transverse momentum-rapidity spectrum of K− produced in inelastic p+p interactions at 31 GeV/c with systematic uncertainties.

Transverse momentum-rapidity spectrum of K− produced in inelastic p+p interactions at 31 GeV/c with statistical uncertainties.

Transverse momentum-rapidity spectrum of K− produced in inelastic p+p interactions at 31 GeV/c with systematic uncertainties.

Transverse momentum-rapidity spectrum of pi− produced in inelastic p+p interactions at 31 GeV/c with statistical uncertainties.

Transverse momentum-rapidity spectrum of pi− produced in inelastic p+p interactions at 31 GeV/c with systematic uncertainties.

Transverse momentum-rapidity spectrum of pi− produced in inelastic p+p interactions at 31 GeV/c with statistical uncertainties.

Transverse momentum-rapidity spectrum of pi− produced in inelastic p+p interactions at 31 GeV/c with systematic uncertainties.

Transverse momentum-rapidity spectrum of pi− produced in inelastic p+p interactions at 31 GeV/c with statistical uncertainties.

Transverse momentum-rapidity spectrum of pi− produced in inelastic p+p interactions at 31 GeV/c with systematic uncertainties.

Transverse momentum-rapidity spectrum of pi− produced in inelastic p+p interactions at 31 GeV/c with statistical uncertainties.

Transverse momentum-rapidity spectrum of pi− produced in inelastic p+p interactions at 31 GeV/c with systematic uncertainties.

Transverse momentum-rapidity spectrum of K− produced in inelastic p+p interactions at 40 GeV/c with statistical uncertainties.

Transverse momentum-rapidity spectrum of K− produced in inelastic p+p interactions at 40 GeV/c with systematic uncertainties.

Transverse momentum-rapidity spectrum of K− produced in inelastic p+p interactions at 40 GeV/c with statistical uncertainties.

Transverse momentum-rapidity spectrum of K− produced in inelastic p+p interactions at 40 GeV/c with systematic uncertainties.

Transverse momentum-rapidity spectrum of pi− produced in inelastic p+p interactions at 40 GeV/c with statistical uncertainties.

Transverse momentum-rapidity spectrum of pi− produced in inelastic p+p interactions at 40 GeV/c with systematic uncertainties.

Transverse momentum-rapidity spectrum of pi− produced in inelastic p+p interactions at 40 GeV/c with statistical uncertainties.

Transverse momentum-rapidity spectrum of pi− produced in inelastic p+p interactions at 40 GeV/c with systematic uncertainties.

Transverse momentum-rapidity spectrum of pi− produced in inelastic p+p interactions at 40 GeV/c with statistical uncertainties.

Transverse momentum-rapidity spectrum of pi− produced in inelastic p+p interactions at 40 GeV/c with systematic uncertainties.

Transverse momentum-rapidity spectrum of pi− produced in inelastic p+p interactions at 40 GeV/c with statistical uncertainties.

Transverse momentum-rapidity spectrum of pi− produced in inelastic p+p interactions at 40 GeV/c with systematic uncertainties.

Transverse momentum-rapidity spectrum of K− produced in inelastic p+p interactions at 80 GeV/c with statistical uncertainties.

Transverse momentum-rapidity spectrum of K− produced in inelastic p+p interactions at 80 GeV/c with systematic uncertainties.

Transverse momentum-rapidity spectrum of K− produced in inelastic p+p interactions at 80 GeV/c with statistical uncertainties.

Transverse momentum-rapidity spectrum of K− produced in inelastic p+p interactions at 80 GeV/c with systematic uncertainties.

Transverse momentum-rapidity spectrum of pi− produced in inelastic p+p interactions at 80 GeV/c with statistical uncertainties.

Transverse momentum-rapidity spectrum of pi− produced in inelastic p+p interactions at 80 GeV/c with systematic uncertainties.

Transverse momentum-rapidity spectrum of pi− produced in inelastic p+p interactions at 80 GeV/c with statistical uncertainties.

Transverse momentum-rapidity spectrum of pi− produced in inelastic p+p interactions at 80 GeV/c with systematic uncertainties.

Transverse momentum-rapidity spectrum of pi− produced in inelastic p+p interactions at 80 GeV/c with statistical uncertainties.

Transverse momentum-rapidity spectrum of pi− produced in inelastic p+p interactions at 80 GeV/c with systematic uncertainties.

Transverse momentum-rapidity spectrum of pi− produced in inelastic p+p interactions at 80 GeV/c with statistical uncertainties.

Transverse momentum-rapidity spectrum of pi− produced in inelastic p+p interactions at 80 GeV/c with systematic uncertainties.

Transverse momentum-rapidity spectrum of K− produced in inelastic p+p interactions at 158 GeV/c with statistical uncertainties.

Transverse momentum-rapidity spectrum of K− produced in inelastic p+p interactions at 158 GeV/c with systematic uncertainties.

Transverse momentum-rapidity spectrum of K− produced in inelastic p+p interactions at 158 GeV/c with statistical uncertainties.

Transverse momentum-rapidity spectrum of K− produced in inelastic p+p interactions at 158 GeV/c with systematic uncertainties.

Transverse momentum-rapidity spectrum of pi− produced in inelastic p+p interactions at 158 GeV/c with statistical uncertainties.

Transverse momentum-rapidity spectrum of pi− produced in inelastic p+p interactions at 158 GeV/c with systematic uncertainties.

Transverse momentum-rapidity spectrum of pi− produced in inelastic p+p interactions at 158 GeV/c with statistical uncertainties.

Transverse momentum-rapidity spectrum of pi− produced in inelastic p+p interactions at 158 GeV/c with systematic uncertainties.

Transverse momentum-rapidity spectrum of pi− produced in inelastic p+p interactions at 158 GeV/c with statistical uncertainties.

Transverse momentum-rapidity spectrum of pi− produced in inelastic p+p interactions at 158 GeV/c with systematic uncertainties.

Transverse momentum-rapidity spectrum of pi− produced in inelastic p+p interactions at 158 GeV/c with statistical uncertainties.

Transverse momentum-rapidity spectrum of pi− produced in inelastic p+p interactions at 158 GeV/c with systematic uncertainties.

Rapidity spectrum of K- from p+p->K-+X at 158 GeV/c with statistical uncertainties.

Rapidity spectrum of K- from p+p->K-+X at 158 GeV/c with systematic uncertainties.

Rapidity spectrum of K+ from p+p->K++X at 158 GeV/c with statistical uncertainties.

Rapidity spectrum of K+ from p+p->K++X at 158 GeV/c with systematic uncertainties.

Rapidity spectrum of pi- from p+p->pi-+X at 158 GeV/c with statistical uncertainties.

Rapidity spectrum of pi- from p+p->pi-+X at 158 GeV/c with systematic uncertainties.

Rapidity spectrum of pi+ from p+p->pi++X at 158 GeV/c with statistical uncertainties.

Rapidity spectrum of pi+ from p+p->pi++X at 158 GeV/c with systematic uncertainties.

Rapidity spectrum of p from p+p->p+X at 158 GeV/c with statistical uncertainties.

Rapidity spectrum of p from p+p->p+X at 158 GeV/c with systematic uncertainties.

Rapidity spectrum of barp from p+p->barp+X at 158 GeV/c with statistical uncertainties.

Rapidity spectrum of barp from p+p->barp+X at 158 GeV/c with systematic uncertainties.

Rapidity spectrum of K- from p+p->K-+X at 80 GeV/c with statistical uncertainties.

Rapidity spectrum of K- from p+p->K-+X at 80 GeV/c with systematic uncertainties.

Rapidity spectrum of K+ from p+p->K++X at 80 GeV/c with statistical uncertainties.

Rapidity spectrum of K+ from p+p->K++X at 80 GeV/c with systematic uncertainties.

Rapidity spectrum of pi- from p+p->pi-+X at 80 GeV/c with statistical uncertainties.

Rapidity spectrum of pi- from p+p->pi-+X at 80 GeV/c with systematic uncertainties.

Rapidity spectrum of pi+ from p+p->pi++X at 80 GeV/c with statistical uncertainties.

Rapidity spectrum of pi+ from p+p->pi++X at 80 GeV/c with systematic uncertainties.

Rapidity spectrum of p from p+p->p+X at 80 GeV/c with statistical uncertainties.

Rapidity spectrum of p from p+p->p+X at 80 GeV/c with systematic uncertainties.

Rapidity spectrum of barp from p+p->barp+X at 80 GeV/c with statistical uncertainties.

Rapidity spectrum of barp from p+p->barp+X at 80 GeV/c with systematic uncertainties.

Rapidity spectrum of K- from p+p->K-+X at 40 GeV/c with statistical uncertainties.

Rapidity spectrum of K- from p+p->K-+X at 40 GeV/c with systematic uncertainties.

Rapidity spectrum of K+ from p+p->K++X at 40 GeV/c with statistical uncertainties.

Rapidity spectrum of K+ from p+p->K++X at 40 GeV/c with systematic uncertainties.

Rapidity spectrum of pi- from p+p->pi-+X at 40 GeV/c with statistical uncertainties.

Rapidity spectrum of pi- from p+p->pi-+X at 40 GeV/c with systematic uncertainties.

Rapidity spectrum of pi+ from p+p->pi++X at 40 GeV/c with statistical uncertainties.

Rapidity spectrum of pi+ from p+p->pi++X at 40 GeV/c with systematic uncertainties.

Rapidity spectrum of p from p+p->p+X at 40 GeV/c with statistical uncertainties.

Rapidity spectrum of p from p+p->p+X at 40 GeV/c with systematic uncertainties.

Rapidity spectrum of barp from p+p->barp+X at 40 GeV/c with statistical uncertainties.

Rapidity spectrum of barp from p+p->barp+X at 40 GeV/c with systematic uncertainties.

Rapidity spectrum of K- from p+p->K-+X at 31 GeV/c with statistical uncertainties.

Rapidity spectrum of K- from p+p->K-+X at 31 GeV/c with systematic uncertainties.

Rapidity spectrum of K+ from p+p->K++X at 31 GeV/c with statistical uncertainties.

Rapidity spectrum of K+ from p+p->K++X at 31 GeV/c with systematic uncertainties.

Rapidity spectrum of pi- from p+p->pi-+X at 31 GeV/c with statistical uncertainties.

Rapidity spectrum of pi- from p+p->pi-+X at 31 GeV/c with systematic uncertainties.

Rapidity spectrum of pi+ from p+p->pi++X at 31 GeV/c with statistical uncertainties.

Rapidity spectrum of pi+ from p+p->pi++X at 31 GeV/c with systematic uncertainties.

Rapidity spectrum of p from p+p->p+X at 31 GeV/c with statistical uncertainties.

Rapidity spectrum of p from p+p->p+X at 31 GeV/c with systematic uncertainties.

Rapidity spectrum of barp from p+p->barp+X at 31 GeV/c with statistical uncertainties.

Rapidity spectrum of barp from p+p->barp+X at 31 GeV/c with systematic uncertainties.

Rapidity spectrum of K- from p+p->K-+X at 20 GeV/c with statistical uncertainties.

Rapidity spectrum of K- from p+p->K-+X at 20 GeV/c with systematic uncertainties.

Rapidity spectrum of K+ from p+p->K++X at 20 GeV/c with statistical uncertainties.

Rapidity spectrum of K+ from p+p->K++X at 20 GeV/c with systematic uncertainties.

Rapidity spectrum of pi- from p+p->pi-+X at 20 GeV/c with statistical uncertainties.

Rapidity spectrum of pi- from p+p->pi-+X at 20 GeV/c with systematic uncertainties.

Rapidity spectrum of pi+ from p+p->pi++X at 20 GeV/c with statistical uncertainties.

Rapidity spectrum of pi+ from p+p->pi++X at 20 GeV/c with systematic uncertainties.

Rapidity spectrum of p from p+p->p+X at 20 GeV/c with statistical uncertainties.

Rapidity spectrum of p from p+p->p+X at 20 GeV/c with systematic uncertainties.

Midrapidity (|y| < 0.1) transverse momentum spectra for (d) K− in Au+Au collisions at √sNN = 7.7 GeV for different centralities. The spectra for centralities other than 0–5% are scaled for clarity as shown in the figure. The curves represent the Bose-Einstein, mT -exponential, and double-exponential function fits to 0–5% central data for pions, kaons, and (anti)protons, respectively. The uncertainties are statistical and systematic added in quadrature.

Midrapidity (|y| < 0.1) transverse momentum spectra for (c) K+ in Au+Au collisions at √sNN = 7.7 GeV for different centralities. The spectra for centralities other than 0–5% are scaled for clarity as shown in the figure. The curves represent the Bose-Einstein, mT -exponential, and double-exponential function fits to 0–5% central data for pions, kaons, and (anti)protons, respectively. The uncertainties are statistical and systematic added in quadrature.

Midrapidity (|y| < 0.1) transverse momentum spectra for (f) p¯ in Au+Au collisions at √sNN = 7.7 GeV for different centralities. The spectra for centralities other than 0–5% are scaled for clarity as shown in the figure. The curves represent the Bose-Einstein, mT -exponential, and double-exponential function fits to 0–5% central data for pions, kaons, and (anti)protons, respectively. The uncertainties are statistical and systematic added in quadrature.

Midrapidity (|y| < 0.1) transverse momentum spectra for (e) p in Au+Au collisions at √sNN = 7.7 GeV for different centralities. The spectra for centralities other than 0–5% are scaled for clarity as shown in the figure. The curves represent the Bose-Einstein, mT -exponential, and double-exponential function fits to 0–5% central data for pions, kaons, and (anti)protons, respectively. The uncertainties are statistical and systematic added in quadrature.

Midrapidity (|y| < 0.1) transverse momentum spectra for (b) π− in Au+Au collisions at √sNN = 11.5 GeV for different centralities. The spectra for centralities other than 0–5% are scaled for clarity as shown in the figure. The curves represent the Bose-Einstein, mT -exponential, and double-exponential function fits to 0–5% central data for pions, kaons, and (anti)protons, respectively. The uncertainties are statistical and systematic added in quadrature.

Midrapidity (|y| < 0.1) transverse momentum spectra for (a) π+ in Au+Au collisions at √sNN = 11.5 GeV for different centralities. The spectra for centralities other than 0–5% are scaled for clarity as shown in the figure. The curves represent the Bose-Einstein, mT -exponential, and double-exponential function fits to 0–5% central data for pions, kaons, and (anti)protons, respectively. The uncertainties are statistical and systematic added in quadrature.

Midrapidity (|y| < 0.1) transverse momentum spectra for (d) K− in Au+Au collisions at √sNN = 11.5 GeV for different centralities. The spectra for centralities other than 0–5% are scaled for clarity as shown in the figure. The curves represent the Bose-Einstein, mT -exponential, and double-exponential function fits to 0–5% central data for pions, kaons, and (anti)protons, respectively. The uncertainties are statistical and systematic added in quadrature.

Midrapidity (|y| < 0.1) transverse momentum spectra for (c) K+ in Au+Au collisions at √sNN = 11.5 GeV for different centralities. The spectra for centralities other than 0–5% are scaled for clarity as shown in the figure. The curves represent the Bose-Einstein, mT -exponential, and double-exponential function fits to 0–5% central data for pions, kaons, and (anti)protons, respectively. The uncertainties are statistical and systematic added in quadrature.

Midrapidity (|y| < 0.1) transverse momentum spectra for (f) p¯ in Au+Au collisions at √sNN = 11.5 GeV for different centralities. The spectra for centralities other than 0–5% are scaled for clarity as shown in the figure. The curves represent the Bose-Einstein, mT -exponential, and double-exponential function fits to 0–5% central data for pions, kaons, and (anti)protons, respectively. The uncertainties are statistical and systematic added in quadrature.

Midrapidity (|y| < 0.1) transverse momentum spectra for (e) p in Au+Au collisions at √sNN = 11.5 GeV for different centralities. The spectra for centralities other than 0–5% are scaled for clarity as shown in the figure. The curves represent the Bose-Einstein, mT -exponential, and double-exponential function fits to 0–5% central data for pions, kaons, and (anti)protons, respectively. The uncertainties are statistical and systematic added in quadrature.

Midrapidity (|y| < 0.1) transverse momentum spectra for (b) π− in Au+Au collisions at √sNN = 19.6 GeV for different centralities. The spectra for centralities other than 0–5% are scaled for clarity as shown in the figure. The curves represent the Bose-Einstein, mT -exponential, and double-exponential function fits to 0–5% central data for pions, kaons, and (anti)protons, respectively. The uncertainties are statistical and systematic added in quadrature.

Midrapidity (|y| < 0.1) transverse momentum spectra for (a) π+ in Au+Au collisions at √sNN = 19.6 GeV for different centralities. The spectra for centralities other than 0–5% are scaled for clarity as shown in the figure. The curves represent the Bose-Einstein, mT -exponential, and double-exponential function fits to 0–5% central data for pions, kaons, and (anti)protons, respectively. The uncertainties are statistical and systematic added in quadrature.

Midrapidity (|y| < 0.1) transverse momentum spectra for (d) K− in Au+Au collisions at √sNN = 19.6 GeV for different centralities. The spectra for centralities other than 0–5% are scaled for clarity as shown in the figure. The curves represent the Bose-Einstein, mT -exponential, and double-exponential function fits to 0–5% central data for pions, kaons, and (anti)protons, respectively. The uncertainties are statistical and systematic added in quadrature.

Midrapidity (|y| < 0.1) transverse momentum spectra for (c) K+ in Au+Au collisions at √sNN = 19.6 GeV for different centralities. The spectra for centralities other than 0–5% are scaled for clarity as shown in the figure. The curves represent the Bose-Einstein, mT -exponential, and double-exponential function fits to 0–5% central data for pions, kaons, and (anti)protons, respectively. The uncertainties are statistical and systematic added in quadrature.

Midrapidity (|y| < 0.1) transverse momentum spectra for (f) p¯ in Au+Au collisions at √sNN = 19.6 GeV for different centralities. The spectra for centralities other than 0–5% are scaled for clarity as shown in the figure. The curves represent the Bose-Einstein, mT -exponential, and double-exponential function fits to 0–5% central data for pions, kaons, and (anti)protons, respectively. The uncertainties are statistical and systematic added in quadrature.

Midrapidity (|y| < 0.1) transverse momentum spectra for (e) p in Au+Au collisions at √sNN = 19.6 GeV for different centralities. The spectra for centralities other than 0–5% are scaled for clarity as shown in the figure. The curves represent the Bose-Einstein, mT -exponential, and double-exponential function fits to 0–5% central data for pions, kaons, and (anti)protons, respectively. The uncertainties are statistical and systematic added in quadrature.

Midrapidity (|y| < 0.1) transverse momentum spectra for (b) π− in Au+Au collisions at √sNN = 27 GeV for different centralities. The spectra for centralities other than 0–5% are scaled for clarity as shown in the figure. The curves represent the Bose-Einstein, mT -exponential, and double-exponential function fits to 0–5% central data for pions, kaons, and (anti)protons, respectively. The uncertainties are statistical and systematic added in quadrature.

Midrapidity (|y| < 0.1) transverse momentum spectra for (a) π+ in Au+Au collisions at √sNN = 27 GeV for different centralities. The spectra for centralities other than 0–5% are scaled for clarity as shown in the figure. The curves represent the Bose-Einstein, mT -exponential, and double-exponential function fits to 0–5% central data for pions, kaons, and (anti)protons, respectively. The uncertainties are statistical and systematic added in quadrature.

Midrapidity (|y| < 0.1) transverse momentum spectra for (d) K− in Au+Au collisions at √sNN = 27 GeV for different centralities. The spectra for centralities other than 0–5% are scaled for clarity as shown in the figure. The curves represent the Bose-Einstein, mT -exponential, and double-exponential function fits to 0–5% central data for pions, kaons, and (anti)protons, respectively. The uncertainties are statistical and systematic added in quadrature.

Midrapidity (|y| < 0.1) transverse momentum spectra for (c) K+ in Au+Au collisions at √sNN = 27 GeV for different centralities. The spectra for centralities other than 0–5% are scaled for clarity as shown in the figure. The curves represent the Bose-Einstein, mT -exponential, and double-exponential function fits to 0–5% central data for pions, kaons, and (anti)protons, respectively. The uncertainties are statistical and systematic added in quadrature.

Midrapidity (|y| < 0.1) transverse momentum spectra for (f) p¯ in Au+Au collisions at √sNN = 27 GeV for different centralities. The spectra for centralities other than 0–5% are scaled for clarity as shown in the figure. The curves represent the Bose-Einstein, mT -exponential, and double-exponential function fits to 0–5% central data for pions, kaons, and (anti)protons, respectively. The uncertainties are statistical and systematic added in quadrature.

Midrapidity (|y| < 0.1) transverse momentum spectra for (e) p in Au+Au collisions at √sNN = 27 GeV for different centralities. The spectra for centralities other than 0–5% are scaled for clarity as shown in the figure. The curves represent the Bose-Einstein, mT -exponential, and double-exponential function fits to 0–5% central data for pions, kaons, and (anti)protons, respectively. The uncertainties are statistical and systematic added in quadrature.

Midrapidity (|y| < 0.1) transverse momentum spectra for (b) π− in Au+Au collisions at √sNN = 39 GeV for different centralities. The spectra for centralities other than 0–5% are scaled for clarity as shown in the figure. The curves represent the Bose-Einstein, mT -exponential, and double-exponential function fits to 0–5% central data for pions, kaons, and (anti)protons, respectively. The uncertainties are statistical and systematic added in quadrature.

Midrapidity (|y| < 0.1) transverse momentum spectra for (a) π+ in Au+Au collisions at √sNN = 39 GeV for different centralities. The spectra for centralities other than 0–5% are scaled for clarity as shown in the figure. The curves represent the Bose-Einstein, mT -exponential, and double-exponential function fits to 0–5% central data for pions, kaons, and (anti)protons, respectively. The uncertainties are statistical and systematic added in quadrature.

Midrapidity (|y| < 0.1) transverse momentum spectra for (d) k- in Au+Au collisions at √sNN = 39 GeV for different centralities. The spectra for centralities other than 0–5% are scaled for clarity as shown in the figure. The curves represent the Bose-Einstein, mT -exponential, and double-exponential function fits to 0–5% central data for pions, kaons, and (anti)protons, respectively. The uncertainties are statistical and systematic added in quadrature.

Midrapidity (|y| < 0.1) transverse momentum spectra for (c) k+ in Au+Au collisions at √sNN = 39 GeV for different centralities. The spectra for centralities other than 0–5% are scaled for clarity as shown in the figure. The curves represent the Bose-Einstein, mT -exponential, and double-exponential function fits to 0–5% central data for pions, kaons, and (anti)protons, respectively. The uncertainties are statistical and systematic added in quadrature.

Midrapidity (|y| < 0.1) transverse momentum spectra for (f) pbar in Au+Au collisions at √sNN = 39 GeV for different centralities. The spectra for centralities other than 0–5% are scaled for clarity as shown in the figure. The curves represent the Bose-Einstein, mT -exponential, and double-exponential function fits to 0–5% central data for pions, kaons, and (anti)protons, respectively. The uncertainties are statistical and systematic added in quadrature.

Midrapidity (|y| < 0.1) transverse momentum spectra for (e) p in Au+Au collisions at √sNN = 39 GeV for different centralities. The spectra for centralities other than 0–5% are scaled for clarity as shown in the figure. The curves represent the Bose-Einstein, mT -exponential, and double-exponential function fits to 0–5% central data for pions, kaons, and (anti)protons, respectively. The uncertainties are statistical and systematic added in quadrature.

Centrality dependence of dN/dy normalized by ⟨Npart⟩/2 for π+, π−, K+, K−, p, and p ̄ at midrapidity (|y|<0.1) in Au+Au collisions at √sNN = 7.7 GeV. Errors shown are the quadrature sum of statistical and systematic uncertainties. For clarity, ⟨Npart⟩ uncertainties are not added in quadrature.

Centrality dependence of dN/dy normalized by ⟨Npart⟩/2 for π+, π−, K+, K−, p, and p ̄ at midrapidity (|y|<0.1) in Au+Au collisions at √sNN = 11.5 GeV. Errors shown are the quadrature sum of statistical and systematic uncertainties. For clarity, ⟨Npart⟩ uncertainties are not added in quadrature.

Centrality dependence of dN/dy normalized by ⟨Npart⟩/2 for π+, π−, K+, K−, p, and p ̄ at midrapidity (|y|<0.1) in Au+Au collisions at √sNN = 19.6 GeV. Errors shown are the quadrature sum of statistical and systematic uncertainties. For clarity, ⟨Npart⟩ uncertainties are not added in quadrature.

Centrality dependence of dN/dy normalized by ⟨Npart⟩/2 for π+, π−, K+, K−, p, and p ̄ at midrapidity (|y|<0.1) in Au+Au collisions at √sNN = 27 GeV. Errors shown are the quadrature sum of statistical and systematic uncertainties. For clarity, ⟨Npart⟩ uncertainties are not added in quadrature.

Centrality dependence of dN/dy normalized by ⟨Npart⟩/2 for π+, π−, K+, K−, p, and p ̄ at midrapidity (|y|<0.1) in Au+Au collisions at √sNN = 39 GeV. Errors shown are the quadrature sum of statistical and systematic uncertainties. For clarity, ⟨Npart⟩ uncertainties are not added in quadrature.

Centrality dependences of <pT> for π+, π−, K+, K−, p, and p ̄ at midrapidity (|y|<0.1) in Au+Au collisions at √sNN = 7.7 GeV. Errors shown are quadrature sum of statistical and systematic uncertainties where the latter dominates.

Centrality dependences of <pT> for π+, π−, K+, K−, p, and p ̄ at midrapidity (|y|<0.1) in Au+Au collisions at √sNN = 11.5 GeV. Errors shown are quadrature sum of statistical and systematic uncertainties where the latter dominates.

Centrality dependences of <pT> for π+, π−, K+, K−, p, and p ̄ at midrapidity (|y|<0.1) in Au+Au collisions at √sNN = 19.6 GeV. Errors shown are quadrature sum of statistical and systematic uncertainties where the latter dominates.

Centrality dependences of <pT> for π+, π−, K+, K−, p, and p ̄ at midrapidity (|y|<0.1) in Au+Au collisions at √sNN = 27 GeV. Errors shown are quadrature sum of statistical and systematic uncertainties where the latter dominates.

Centrality dependences of <pT> for π+, π−, K+, K−, p, and p ̄ at midrapidity (|y|<0.1) in Au+Au collisions at √sNN = 39 GeV. Errors shown are quadrature sum of statistical and systematic uncertainties where the latter dominates.

Variation of π−/π+, K−/K+, and p ̄/p ratios as a function of ⟨Npart⟩ at midrapidity (|y| < 0.1) in Au+Au collisions at 7.7 GeV. Errors shown are the quadrature sum of statistical and systematic uncertainties where the latter dominates.

Variation of π−/π+, K−/K+, and p ̄/p ratios as a function of ⟨Npart⟩ at midrapidity (|y| < 0.1) in Au+Au collisions at 11.5 GeV. Errors shown are the quadrature sum of statistical and systematic uncertainties where the latter dominates.

Variation of π−/π+, K−/K+, and p ̄/p ratios as a function of ⟨Npart⟩ at midrapidity (|y| < 0.1) in Au+Au collisions at 19.6 GeV. Errors shown are the quadrature sum of statistical and systematic uncertainties where the latter dominates.

Variation of π−/π+, K−/K+, and p ̄/p ratios as a function of ⟨Npart⟩ at midrapidity (|y| < 0.1) in Au+Au collisions at 27 GeV. Errors shown are the quadrature sum of statistical and systematic uncertainties where the latter dominates.

Variation of π−/π+, K−/K+, and p ̄/p ratios as a function of ⟨Npart⟩ at midrapidity (|y| < 0.1) in Au+Au collisions at 39 GeV. Errors shown are the quadrature sum of statistical and systematic uncertainties where the latter dominates.

Variation of K−/π−, p ̄/π−, K+/π+, and p/π+ ratios as a function of ⟨Npart⟩ at midrapidity (|y| < 0.1) in Au+Au collisions at 7.7 GeV. Errors shown are the quadrature sum of statistical and systematic uncertainties where the latter dominates.

Variation of K−/π−, p ̄/π−, K+/π+, and p/π+ ratios as a function of ⟨Npart⟩ at midrapidity (|y| < 0.1) in Au+Au collisions at 11.5 GeV. Errors shown are the quadrature sum of statistical and systematic uncertainties where the latter dominates.

Variation of K−/π−, p ̄/π−, K+/π+, and p/π+ ratios as a function of ⟨Npart⟩ at midrapidity (|y| < 0.1) in Au+Au collisions at 19.6 GeV. Errors shown are the quadrature sum of statistical and systematic uncertainties where the latter dominates.

Variation of K−/π−, p ̄/π−, K+/π+, and p/π+ ratios as a function of ⟨Npart⟩ at midrapidity (|y| < 0.1) in Au+Au collisions at 27 GeV. Errors shown are the quadrature sum of statistical and systematic uncertainties where the latter dominates.

Variation of K−/π−, p ̄/π−, K+/π+, and p/π+ ratios as a function of ⟨Npart⟩ at midrapidity (|y| < 0.1) in Au+Au collisions at 39 GeV. Errors shown are the quadrature sum of statistical and systematic uncertainties where the latter dominates.

The midrapidity (|y| < 0.1) dN/dy normalized by ⟨Npart⟩/2 as a function of √sNN for π±, K±, and p and p ̄ in 0–5% Au+Au collisions at BES energies. Errors shown are the quadrature sum of statistical and systematic uncertainties where the latter dominates.

⟨mT⟩ − m of π±, K±, and p and p ̄ as a function of √sNN . Midrapidity (|y| < 0.1) results are shown for 0–5% central Au+Au collisions at BES energies. The errors shown are the quadrature sum of statistical and systematic uncertainties where the latter dominates.

π−/π+, K−/K+, and p ̄/p ratios at midrapidity (|y| < 0.1) in central 0–5% Au+Au collisions at √sNN = 7.7, 11.5, 19.6, 27, and 39 GeV. Errors shown are the quadrature sum of statistical and systematic uncertainties where the latter dominates.

K/π ratio at midrapidity (|y| < 0.1) for central 0–5% Au+Au collisions at √sNN = 7.7, 11.5, 19.6, 27, and 39 GeV. Errors shown are the quadrature sum of statistical and systematic uncertainties where the latter dominates.

The GCE model particle yields fits shown along with standard deviations for Au+Au 7.7 and Au+Au 39 GeV in 0–5% central collisions. Uncertainties on experimental data represent statistical and systematic uncertainties added in quadrature.

The GCE model particle ratios fits shown along with standard deviations for Au+Au 7.7 and Au+Au 39 GeV in 0–5% central collisions. Uncertainties on experimental data represent statistical and systematic uncertainties added in quadrature.

The SCE model particle yields fits shown along with standard deviations for Au+Au 7.7 and Au+Au 39 GeV in 0–5% central collisions. Uncertainties on experimental data represent statistical and systematic uncertainties added in quadrature.

The SCE model particle ratios fits shown along with standard deviations for Au+Au 7.7 and Au+Au 39 GeV in 0–5% central collisions. Uncertainties on experimental data represent statistical and systematic uncertainties added in quadrature.

Chemical freeze-out parameter γS plotted vs ⟨Npart⟩ in GCE for particle yields fit. Uncertainties represent systematic errors.

Chemical freeze-out parameter μB plotted vs ⟨Npart⟩ in GCE for particle yields fit. Uncertainties represent systematic errors.

Chemical freeze-out parameter μS plotted vs ⟨Npart⟩ in GCE for particle yields fit. Uncertainties represent systematic errors.

Chemical freeze-out parameter Tch plotted vs ⟨Npart⟩ in GCE for particle yields fit. Uncertainties represent systematic errors.

Chemical freeze-out parameter R plotted vs ⟨Npart⟩ in GCE for particle yields fit. Uncertainties represent systematic errors.

Ratio of chemical freeze-out parameter γS between results from particle yield fits to particle ratio fits in GCE plotted vs ⟨Npart⟩. Uncertainties represent systematic errors.

Ratio of chemical freeze-out parameter μB between results from particle yield fits to particle ratio fits in GCE plotted vs ⟨Npart⟩. Uncertainties represent systematic errors.

Ratio of chemical freeze-out parameter μS between results from particle yield fits to particle ratio fits in GCE plotted vs ⟨Npart⟩. Uncertainties represent systematic errors.

Ratio of chemical freeze-out parameter Tch between results from particle yield fits to particle ratio fits in GCE plotted vs ⟨Npart⟩. Uncertainties represent systematic errors.

Chemical freeze-out parameter γS plotted vs ⟨Npart⟩ in SCE for particle yields fit. Uncertainties represent systematic errors.

Chemical freeze-out parameter μB plotted vs ⟨Npart⟩ in SCE for particle yields fit. Uncertainties represent systematic errors.

Chemical freeze-out parameter Tch plotted vs ⟨Npart⟩ in SCE for particle yields fit. Uncertainties represent systematic errors.

Chemical freeze-out parameter R plotted vs ⟨Npart⟩ in SCE for particle yields fit. Uncertainties represent systematic errors.

Ratio of chemical freeze-out parameter γS between yield and ratio fits in SCE plotted vs ⟨Npart⟩. Uncertainties represent systematic errors.

Ratio of chemical freeze-out parameter μB between yield and ratio fits in SCE plotted vs ⟨Npart⟩. Uncertainties represent systematic errors.

Ratio of chemical freeze-out parameter Tch between yield and ratio fits in SCE plotted vs ⟨Npart⟩. Uncertainties represent systematic errors.

Ratio of chemical freeze-out parameter γS between GCE and SCE results using particle ratios in fits plotted vs ⟨Npart⟩. Uncertainties represent systematic errors.

Ratio of chemical freeze-out parameter μB between GCE and SCE results using particle ratios in fits plotted vs ⟨Npart⟩. Uncertainties represent systematic errors.

Ratio of chemical freeze-out parameter Tch between GCE and SCE results using particle ratios in fits plotted vs ⟨Npart⟩. Uncertainties represent systematic errors.

Ratio of chemical freeze-out parameter γS between GCE and SCE results using particle yields in fits plotted vs ⟨Npart⟩. Uncertainties represent systematic errors.

Ratio of chemical freeze-out parameter μB between GCE and SCE results using particle yields in fits plotted vs ⟨Npart⟩. Uncertainties represent systematic errors.

Ratio of chemical freeze-out parameter Tch between GCE and SCE results using particle yields in fits plotted vs ⟨Npart⟩. Uncertainties represent systematic errors.

Ratio of chemical freeze-out parameter R between GCE and SCE results using particle yields in fits plotted vs ⟨Npart⟩. Uncertainties represent systematic errors.

Extracted chemical freeze-out temperature vs baryon chemical potential for (a) GCE and (b) SCE cases using particle yields as input for fitting. Curves represent two model predictions [81,82]. The gray bands represent the theoretical prediction ranges of the Cleymans et al. model [81]. Uncertainties represent systematic errors.

Extracted chemical freeze-out temperature vs baryon chemical potential for (a) GCE and (b) SCE cases using particle yields as input for fitting. Curves represent two model predictions [81,82]. The gray bands represent the theoretical prediction ranges of the Cleymans et al. model [81]. Uncertainties represent systematic errors.

Extracted chemical freeze-out temperature vs baryon chemical potential for (a) GCE and (b) SCE cases using particle yields as input for fitting. Curves represent two model predictions [81,82]. The gray bands represent the theoretical prediction ranges of the Cleymans et al. model [81]. Uncertainties represent systematic errors.

"Choice on constraints: Extracted chemical freeze-out temperatures shown in panels (a), (c), and (e) and baryon chemical potentials shown in panels (b), (d), and (f) for GCE using particle yields as input for fitting, respectively, for Au+Au collisions at √sNN = 7.7, 19.6, and 39 GeV. Results are compared for three initial conditions: μQ = 0, μQ constrained to B/2Q value, and μQ constrained to B/2Q along with μS constrained to 0. Uncertainties represent systematic errors."

"Choice on constraints: Extracted chemical freeze-out temperatures shown in panels (a), (c), and (e) and baryon chemical potentials shown in panels (b), (d), and (f) for GCE using particle yields as input for fitting, respectively, for Au+Au collisions at √sNN = 7.7, 19.6, and 39 GeV. Results are compared for three initial conditions: μQ = 0, μQ constrained to B/2Q value, and μQ constrained to B/2Q along with μS constrained to 0. Uncertainties represent systematic errors."

"Choice on constraints: Extracted chemical freeze-out temperatures shown in panels (a), (c), and (e) and baryon chemical potentials shown in panels (b), (d), and (f) for GCE using particle yields as input for fitting, respectively, for Au+Au collisions at √sNN = 7.7, 19.6, and 39 GeV. Results are compared for three initial conditions: μQ = 0, μQ constrained to B/2Q value, and μQ constrained to B/2Q along with μS constrained to 0. Uncertainties represent systematic errors."

"Choice on constraints: Extracted chemical freeze-out temperatures shown in panels (a), (c), and (e) and baryon chemical potentials shown in panels (b), (d), and (f) for GCE using particle yields as input for fitting, respectively, for Au+Au collisions at √sNN = 7.7, 19.6, and 39 GeV. Results are compared for three initial conditions: μQ = 0, μQ constrained to B/2Q value, and μQ constrained to B/2Q along with μS constrained to 0. Uncertainties represent systematic errors."

"Choice on constraints: Extracted chemical freeze-out temperatures shown in panels (a), (c), and (e) and baryon chemical potentials shown in panels (b), (d), and (f) for GCE using particle yields as input for fitting, respectively, for Au+Au collisions at √sNN = 7.7, 19.6, and 39 GeV. Results are compared for three initial conditions: μQ = 0, μQ constrained to B/2Q value, and μQ constrained to B/2Q along with μS constrained to 0. Uncertainties represent systematic errors."

"Choice on constraints: Extracted chemical freeze-out temperatures shown in panels (a), (c), and (e) and baryon chemical potentials shown in panels (b), (d), and (f) for GCE using particle yields as input for fitting, respectively, for Au+Au collisions at √sNN = 7.7, 19.6, and 39 GeV. Results are compared for three initial conditions: μQ = 0, μQ constrained to B/2Q value, and μQ constrained to B/2Q along with μS constrained to 0. Uncertainties represent systematic errors."

"Choice on including more particles: Extracted chemical freeze-out parameters (a) Tch, (b) μB, and (c) γS along with (d) χ2/ndf for GCE using particle yields as input for fitting. Results are compared for Au+Au collisions at √sNN = 39 GeV for four different sets of particle yields used in fitting. Uncertainties represent systematic errors."

"Choice on including more particles: Extracted chemical freeze-out parameters (a) Tch, (b) μB, and (c) γS along with (d) χ2/ndf for GCE using particle yields as input for fitting. Results are compared for Au+Au collisions at √sNN = 39 GeV for four different sets of particle yields used in fitting. Uncertainties represent systematic errors."

"Choice on including more particles: Extracted chemical freeze-out parameters (a) Tch, (b) μB, and (c) γS along with (d) χ2/ndf for GCE using particle yields as input for fitting. Results are compared for Au+Au collisions at √sNN = 39 GeV for four different sets of particle yields used in fitting. Uncertainties represent systematic errors."

"Choice on including more particles: Extracted chemical freeze-out parameters (a) Tch, (b) μB, and (c) γS along with (d) χ2/ndf for GCE using particle yields as input for fitting. Results are compared for Au+Au collisions at √sNN = 39 GeV for four different sets of particle yields used in fitting. Uncertainties represent systematic errors."

"Blast wave model fits of π±, K±, p and p p¯ T spectra in 0–5% central Au+Au collisions at √sNN = (a) 7.7, (b) 11.5, (c) 19.6, (d) 27, and (e) 39 GeV. Uncertainties on experimental data represent statistical and systematic uncertainties added in quadrature. Here, the uncertainties are smaller than the symbol size."

"Blast wave model fits of π±, K±, p and p p¯ T spectra in 0–5% central Au+Au collisions at √sNN = (a) 7.7, (b) 11.5, (c) 19.6, (d) 27, and (e) 39 GeV. Uncertainties on experimental data represent statistical and systematic uncertainties added in quadrature. Here, the uncertainties are smaller than the symbol size."

"Blast wave model fits of π±, K±, p and p p¯ T spectra in 0–5% central Au+Au collisions at √sNN = (a) 7.7, (b) 11.5, (c) 19.6, (d) 27, and (e) 39 GeV. Uncertainties on experimental data represent statistical and systematic uncertainties added in quadrature. Here, the uncertainties are smaller than the symbol size."

"Blast wave model fits of π±, K±, p and p p¯ T spectra in 0–5% central Au+Au collisions at √sNN = (a) 7.7, (b) 11.5, (c) 19.6, (d) 27, and (e) 39 GeV. Uncertainties on experimental data represent statistical and systematic uncertainties added in quadrature. Here, the uncertainties are smaller than the symbol size."

"Blast wave model fits of π±, K±, p and p p¯ T spectra in 0–5% central Au+Au collisions at √sNN = (a) 7.7, (b) 11.5, (c) 19.6, (d) 27, and (e) 39 GeV. Uncertainties on experimental data represent statistical and systematic uncertainties added in quadrature. Here, the uncertainties are smaller than the symbol size."

"Variation of Tkin with <β> for different energies and centralities. The centrality increases from left to right for a given energy. The data points other than BES energies are taken from Refs. [43,66]. Uncertainties represent systematic uncertainties."

"Variation of Tkin with <β> for different energies and centralities. The centrality increases from left to right for a given energy. The data points other than BES energies are taken from Refs. [43,66]. Uncertainties represent systematic uncertainties."

"Variation of Tkin with <β> for different energies and centralities. The centrality increases from left to right for a given energy. The data points other than BES energies are taken from Refs. [43,66]. Uncertainties represent systematic uncertainties."

"Variation of Tkin with <β> for different energies and centralities. The centrality increases from left to right for a given energy. The data points other than BES energies are taken from Refs. [43,66]. Uncertainties represent systematic uncertainties."

"Variation of Tkin with <β> for different energies and centralities. The centrality increases from left to right for a given energy. The data points other than BES energies are taken from Refs. [43,66]. Uncertainties represent systematic uncertainties."

" (a) Energy dependence of kinetic and chemical freezeout temperatures for central heavy-ion collisions. The curves represent various theoretical predictions [81,82]. (b) Energy dependence of average transverse radial flow velocity for central heavy-ion collisions. The data points other than BES energies are taken from Refs. [43,53–64,66] and references therein. The BES data points are for 0–5% central collisions, AGS energies are mostly for 0–5%, SPS energies are for mostly 0–7%, and top RHIC and LHC energies are for 0–5% central collisions. Uncertainties represent systematic uncertainties."