(Color online.) $J/\psi$ RAA versus pT for several centrality bins for Au+Au collisions at $\sqrt{s_{NN}}$ = 200 GeV. The statistical (systematic) uncertainties are shown with vertical bars (open boxes). The filled boxes about unity on the right show the size of the normalization uncertainty. PHENIX low-$p_{T}$ $J/\psi$ results [6] and STAR high-$p_{T}$ results in Cu+Cu collisions [15] are shown for comparison. The curves are from the predictions by Model I (Liu et al.) [11] and Model II (Zhao et al.) [12].

(Color online) $R_{AA}$ versus $N_{part}$ for high-$p_{T}$ $J/\psi$ in comparison with a) low-$p_{T}$ $J/\psi$ data from PHENIX [6], and high-$p_{T}$ $\pi^{\pm}$ from STAR [41, 42], and b) high-$p_{T}$ $J/\psi$ from CMS [8]. The statistical (systematic) uncertainties are shown in vertical bars (boxes). The shaded green band about unity shows the systematic uncertainties from $N_{bin}$ and the box about unity on the right shows the $R_{AA}$ normalization uncertainty from the statistical and global systematic uncertainties of the $p+p$ reference data. The percentage in the top horizontal axis of (a) and (b) refer to the collision centrality for high-$p_{T}$ $J/\psi$ data at the RHIC and LHC, respectively.

$A_N$ vs. $x_F$ at average pseudorapidity of 3.68 for $\pi^0$ and $\eta$. Also shown are the previously published results for $\pi^0$ at lower $x_F$ , derived from the same data set but without the center cut [12]. The error bars are statistical uncertainties only. The error boxes indicate the systematic uncertainties.

(Color online) Distribution of the product of the TPC reconstructed charge sign and 1/$p_{T}$ for candidates satisfying all the W signal selection criteria and $E^{e}_{T}$ > 25 GeV.

(Color online) Distributions of the invariant mass of pairs of oppositely charged $e^{\pm}$ candidates after sequentially applying the selection criteria described in Sec. IV A to both $e^{\pm}$ candidates.

(Color online) Distributions of the invariant mass of $Z/\gamma^{*} \to e^{+}e^{−}$ candidate events satisfying all selection criteria described in Sec. IVC. The $Z/\gamma^{*} \to e^{+}e^{−}$ MC distribution (dashed line) is shown for comparison. Note the larger bin widths relative to Fig. 7

(Color online) $E^{e}_{T}$ distribution of $W^{+}$ (top) and $W^{−}$ (bottom) candidate events, background components, and $W \to e\nu$ MC signal for comparison. Note the factor of two difference in the vertical scales.

(Color online) $E^{e}_{T}$ distribution of $W^{+}$ (top) and $W^{−}$ (bottom) candidate events, background components, and $W \to e\nu$ MC signal for comparison. Note the factor of two difference in the vertical scales.

(Color online) Lepton pseudorapidity distribution of $W^{+}$ (left) and $W^{−}$ (right) candidate events, background components, and $W \to e\nu$ MC signal for comparison.

(Color online) Lepton pseudorapidity distribution of $W^{+}$ (left) and $W^{−}$ (right) candidate events, background components, and $W \to e\nu$ MC signal for comparison.

(Color online) W cross section ratio, $R_{W}$, for the two $e^{\pm}$ pseudorapidity bins. Theory calculations at NLO from the fewz program using the MSTW08 and CTEQ 6.6 PDF sets (with 90% confidence level error eigenvector uncertainties) are shown for comparison.

The same as the left but for Cu+Cu collisions. The systematic errors are shown as thin lines with wide caps at the ends and statistical errors are shown as thick lines with small caps at the end. Statistical and systematic errors are very small.

The difference of charge-independent (CI) v2{2} and like-sign (LS) $v_2\{2\}$ for Au+Au and Cu+Cu collisions at 200 (top panel) and 62.4 (bottom panel) GeV vs. the log of $\langle dN_{ch}/d\eta\rangle$.The statistical errors are smaller than the marker size and not visible for most of the data.

The difference of charge-independent (CI) v2{2} and like-sign (LS) $v_2\{2\}$ for Au+Au and Cu+Cu collisions at 200 (top panel) and 62.4 (bottom panel) GeV vs. the log of $\langle dN_{ch}/d\eta\rangle$.The statistical errors are smaller than the marker size and not visible for most of the data.

The difference of charge-independent (CI) v2{2} and like-sign (LS) $v_2\{2\}$ for Au+Au and Cu+Cu collisions at 200 (top panel) and 62.4 (bottom panel) GeV vs. the log of $\langle dN_{ch}/d\eta\rangle$.The statistical errors are smaller than the marker size and not visible for most of the data.

The difference of charge-independent (CI) v2{2} and like-sign (LS) $v_2\{2\}$ for Au+Au and Cu+Cu collisions at 200 (top panel) and 62.4 (bottom panel) GeV vs. the log of $\langle dN_{ch}/d\eta\rangle$.The statistical errors are smaller than the marker size and not visible for most of the data.

The LS and CI four-particle cumulant $v_2\{4\}^4$ for Au+Au collisions at 200 and 62.4 GeV for $0.15 < pT < 2.0$ GeV/$c$. The systematic errors are shown as narrow lines with wide caps at the end and statistical errors are shown as thick lines with narrow caps at the end. Statistical errors are not visible for most of the points.

The LS and CI four-particle cumulant $v_2\{4\}^4$ for Au+Au collisions at 200 and 62.4 GeV for $0.15 < pT < 2.0$ GeV/$c$. The systematic errors are shown as narrow lines with wide caps at the end and statistical errors are shown as thick lines with narrow caps at the end. Statistical errors are not visible for most of the points.

The LS and CI four-particle cumulant $v_2\{4\}^4$ for Cu+Cu collisions at 200 and 62.4 GeV for $0.15 < p_T < 2.0$ GeV/c. The most central points (two points for Cu+Cu 62.4 GeV) gives $v_2\{4\}^4 < 0$ for all the data sets. The negative values are probably due to large fluctuations in agreement with Eq. (1). These may include contributions from impact parameter spread and finite multiplicity bin width.

The LS and CI four-particle cumulant $v_2\{4\}^4$ for Cu+Cu collisions at 200 and 62.4 GeV for $0.15 < p_T < 2.0$ GeV/c. The most central points (two points for Cu+Cu 62.4 GeV) gives $v_2\{4\}^4 < 0$ for all the data sets. The negative values are probably due to large fluctuations in agreement with Eq. (1). These may include contributions from impact parameter spread and finite multiplicity bin width.

The difference of charge-independent (CI) $v_2\{4\}$ and like-sign (LS) $v_2\{4\}$ for Au+Au collisions at 200 and 62.4 GeV vs. the log of $\langle dN_{ch}/d\eta\rangle$.

The difference of charge-independent (CI) $v_2\{4\}$ and like-sign (LS) $v_2\{4\}$ for Au+Au collisions at 200 and 62.4 GeV vs. the log of $\langle dN_{ch}/d\eta\rangle$.

(Left) The difference between $v_2\{2\}^2$ and $v_2\{4\}^2$ for 200 GeV Au+Au and Cu+Cu collisions for both LS and CI combinations.

(Left) The difference between $v_2\{2\}^2$ and $v_2\{4\}^2$ for 200 GeV Au+Au and Cu+Cu collisions for both LS and CI combinations.

(Right) The difference between $v_2\{2\}^2$ and $v_2\{4\}^2$ for 62.4 GeV Au+Au and Cu+Cu collisions for both LS and CI combinations. The statistical and systematic errors are shown as in previous figures.

(Right) The difference between $v_2\{2\}^2$ and $v_2\{4\}^2$ for 62.4 GeV Au+Au and Cu+Cu collisions for both LS and CI combinations. The statistical and systematic errors are shown as in previous figures.

The upper limit on $\sigma_{v_2}/\langle v_2 \rangle$ for 200 GeV (left) and 62.4 GeV (right) Au+Au collisions from Eq. (9) compared to $\sigma_\varepsilon/\varepsilon$ from Eq. (10) for three different models. The upper limit is found using the LS results for $v_2\{2\}$. Data are from the range $0.15 < p_T < 2.0$ GeV/$c$. The shaded bands reflect the uncertainties on the models which are dominated by uncertainty on the distribution of nucleons inside the nucleus. The uncertainty is only shown for the MCG-N and fKLN-CGC models. The uncertainty on the MCG-Q model is the same as for the MCG-N model but is not shown for the visual clarity.

The upper limit on $\sigma_{v_2}/\langle v_2 \rangle$ for 200 GeV (left) and 62.4 GeV (right) Au+Au collisions from Eq. (9) compared to $\sigma_\varepsilon/\varepsilon$ from Eq. (10) for three different models. The upper limit is found using the LS results for $v_2\{2\}$. Data are from the range $0.15 < p_T < 2.0$ GeV/$c$. The shaded bands reflect the uncertainties on the models which are dominated by uncertainty on the distribution of nucleons inside the nucleus. The uncertainty is only shown for the MCG-N and fKLN-CGC models. The uncertainty on the MCG-Q model is the same as for the MCG-N model but is not shown for the visual clarity.

The upper limit on $\sigma_{v_2}/\langle v_2 \rangle$ for 200 GeV (left) and 62.4 GeV (right) Au+Au collisions from Eq. (9) compared to $\sigma_\varepsilon/\varepsilon$ from Eq. (10) for three different models. The upper limit is found using the LS results for $v_2\{2\}$. Data are from the range $0.15 < p_T < 2.0$ GeV/$c$. The shaded bands reflect the uncertainties on the models which are dominated by uncertainty on the distribution of nucleons inside the nucleus. The uncertainty is only shown for the MCG-N and fKLN-CGC models. The uncertainty on the MCG-Q model is the same as for the MCG-N model but is not shown for the visual clarity.

The upper limit on $\sigma_{v_2}/\langle v_2 \rangle$ for 200 GeV (left) and 62.4 GeV (right) Au+Au collisions from Eq. (9) compared to $\sigma_\varepsilon/\varepsilon$ from Eq. (10) for three different models. The upper limit is found using the LS results for $v_2\{2\}$. Data are from the range $0.15 < p_T < 2.0$ GeV/$c$. The shaded bands reflect the uncertainties on the models which are dominated by uncertainty on the distribution of nucleons inside the nucleus. The uncertainty is only shown for the MCG-N and fKLN-CGC models. The uncertainty on the MCG-Q model is the same as for the MCG-N model but is not shown for the visual clarity.

The STAR data compared to PHOBOS data [34] on $\sigma_{v_2}/\langle v_2 \rangle$ with $\delta_2$ for $\Delta\eta > 2$ taken to be zero (see Fig. 6 from Ref. [34]). The shaded band shows the errors quoted from Ref. [34].

The STAR data compared to PHOBOS data [34] on $\sigma_{v_2}/\langle v_2 \rangle$ with $\delta_2$ for $\Delta\eta > 2$ taken to be zero (see Fig. 6 from Ref. [34]). The shaded band shows the errors quoted from Ref. [34].

The upper limit on $\sigma_{v_2}/\langle v_2 \rangle$ for 200 GeV (left) and 62.4 GeV (right) Cu+Cu collisions from Eq. (9) compared to $\sigma_\varepsilon/\varepsilon$ from Eq. (10) for three different models.

The upper limit on $\sigma_{v_2}/\langle v_2 \rangle$ for 200 GeV (left) and 62.4 GeV (right) Cu+Cu collisions from Eq. (9) compared to $\sigma_\varepsilon/\varepsilon$ from Eq. (10) for three different models.

The upper limit on $\sigma_{v_2}/\langle v_2 \rangle$ for 200 GeV (left) and 62.4 GeV (right) Cu+Cu collisions from Eq. (9) compared to $\sigma_\varepsilon/\varepsilon$ from Eq. (10) for three different models.

The upper limit on $\sigma_{v_2}/\langle v_2 \rangle$ for 200 GeV (left) and 62.4 GeV (right) Cu+Cu collisions from Eq. (9) compared to $\sigma_\varepsilon/\varepsilon$ from Eq. (10) for three different models.

$\rho^0$ production cross section determined with the data set collected with trigger B, as a function of momentum transfer $t$, fitted with a double exponential fit function in the range $t < 0.1$. The statistical errors are shown. The ratio of incoherent to coherent cross sections is measured to be $0.20 \pm 0.8$. The fit parameters are shown in Table I.

Comparison of theoretical predictions to the measured total cross section for coherent $\rho^0$ production as a functionof $\sqrt{s_{NN}}$ . The measured cross section is based on the dataset collected with trigger B. The error bars show the sum of the statistical and systematic uncertainties. See text for details.

Efficiency corrected pT spectra for the different centrality bins and for the various particles. Note that 7 centrality bins have been used for the KS0 and the Λ while only 6 and 3 have been used for the Ξ and Ω, respectively. Errors are statistical only. The Λ spectra are corrected for the feed-down of the Ξ decay.

Efficiency corrected pT spectra for the different centrality bins and for the various particles. Note that 7 centrality bins have been used for the KS0 and the Λ while only 6 and 3 have been used for the Ξ and Ω, respectively. Errors are statistical only. The Λ spectra are corrected for the feed-down of the Ξ decay.

Efficiency corrected pT spectra for the different centrality bins and for the various particles. Note that 7 centrality bins have been used for the KS0 and the Λ while only 6 and 3 have been used for the Ξ and Ω, respectively. Errors are statistical only. The Λ spectra are corrected for the feed-down of the Ξ decay.

Efficiency corrected pT spectra for the different centrality bins and for the various particles. Note that 7 centrality bins have been used for the KS0 and the Λ while only 6 and 3 have been used for the Ξ and Ω, respectively. Errors are statistical only. The Λ spectra are corrected for the feed-down of the Ξ decay.

Efficiency corrected pT spectra for the different centrality bins and for the various particles. Note that 7 centrality bins have been used for the KS0 and the Λ while only 6 and 3 have been used for the Ξ and Ω, respectively. Errors are statistical only. The Λ spectra are corrected for the feed-down of the Ξ decay.

Efficiency corrected pT spectra for the different centrality bins and for the various particles. Note that 7 centrality bins have been used for the KS0 and the Λ while only 6 and 3 have been used for the Ξ and Ω, respectively. Errors are statistical only. The Λ spectra are corrected for the feed-down of the Ξ decay.

Efficiency corrected pT spectra for the different centrality bins and for the various particles. Note that 7 centrality bins have been used for the KS0 and the Λ while only 6 and 3 have been used for the Ξ and Ω, respectively. Errors are statistical only. The Λ spectra are corrected for the feed-down of the Ξ decay.

Efficiency corrected pT spectra for the different centrality bins and for the various particles. Note that 7 centrality bins have been used for the KS0 and the Λ while only 6 and 3 have been used for the Ξ and Ω, respectively. Errors are statistical only. The Λ spectra are corrected for the feed-down of the Ξ decay.

Efficiency corrected pT spectra for the different centrality bins and for the various particles. Note that 7 centrality bins have been used for the KS0 and the Λ while only 6 and 3 have been used for the Ξ and Ω, respectively. Errors are statistical only. The Λ spectra are corrected for the feed-down of the Ξ decay.

Efficiency corrected pT spectra for the different centrality bins and for the various particles. Note that 7 centrality bins have been used for the KS0 and the Λ while only 6 and 3 have been used for the Ξ and Ω, respectively. Errors are statistical only. The Λ spectra are corrected for the feed-down of the Ξ decay.

Efficiency corrected pT spectra for the different centrality bins and for the various particles. Note that 7 centrality bins have been used for the KS0 and the Λ while only 6 and 3 have been used for the Ξ and Ω, respectively. Errors are statistical only. The Λ spectra are corrected for the feed-down of the Ξ decay.

Efficiency corrected pT spectra for the different centrality bins and for the various particles. Note that 7 centrality bins have been used for the KS0 and the Λ while only 6 and 3 have been used for the Ξ and Ω, respectively. Errors are statistical only. The Λ spectra are corrected for the feed-down of the Ξ decay.

Efficiency corrected pT spectra for the different centrality bins and for the various particles. Note that 7 centrality bins have been used for the KS0 and the Λ while only 6 and 3 have been used for the Ξ and Ω, respectively. Errors are statistical only. The Λ spectra are corrected for the feed-down of the Ξ decay.

Extrapolated average transverse momenta ⟨pT ⟩ as a function of dNch/dy for different particle species in Au+Au collisions at 62.4 GeV. Statistical uncertainties are represented by the error bars at the points while the systematic uncertainties are represented by the gray bars. The π, charged K and p data were extracted from Ref. [14].

Extrapolated average transverse momenta ⟨pT ⟩ as a function of dNch/dy for different particle species in Au+Au collisions at 62.4 GeV. Statistical uncertainties are represented by the error bars at the points while the systematic uncertainties are represented by the gray bars. The π, charged K and p data were extracted from Ref. [14].

KS0 dN/dpT spectra compared to the charged Kaon spectra for the event centrality of 0-5% and 30-40%. The charged Kaons data points are for rapidity range of |y| < 0.1 and were extracted from Ref. [14].

KS0 dN/dpT spectra compared to the charged Kaon spectra for the event centrality of 0-5% and 30-40%. The charged Kaons data points are for rapidity range of |y| < 0.1 and were extracted from Ref. [14].

KS0 dN/dpT spectra compared to the charged Kaon spectra for the event centrality of 0-5% and 30-40%. The charged Kaons data points are for rapidity range of |y| < 0.1 and were extracted from Ref. [14].

KS0 dN/dpT spectra compared to the charged Kaon spectra for the event centrality of 0-5% and 30-40%. The charged Kaons data points are for rapidity range of |y| < 0.1 and were extracted from Ref. [14].

Strange particle production yields at mid-rapidity in central Au+Au and Pb+Pb collisions versus the center of mass energy √sNN. The top panel shows results for K0S and Λ. The AGS values are from E896 [1] (centrality 0 − 5 %). The SPS values are from NA49 [20] (centrality 0 − 7 %) and the RHIC values are from STAR [4, 15] (centrality 0 − 5 %). For the multi-strange baryons Ξ and Ω (bottom panel), the SPS results are from NA57 [2] (centrality 0 − 11 %) and the RHIC values are from STAR [15, 21] (centrality 0 − 20 %).

Strange particle production yields at mid-rapidity in central Au+Au and Pb+Pb collisions versus the center of mass energy √sNN. The top panel shows results for K0S and Λ. The AGS values are from E896 [1] (centrality 0 − 5 %). The SPS values are from NA49 [20] (centrality 0 − 7 %) and the RHIC values are from STAR [4, 15] (centrality 0 − 5 %). For the multi-strange baryons Ξ and Ω (bottom panel), the SPS results are from NA57 [2] (centrality 0 − 11 %) and the RHIC values are from STAR [15, 21] (centrality 0 − 20 %).

Strange particle production yields at mid-rapidity in central Au+Au and Pb+Pb collisions versus the center of mass energy √sNN. The top panel shows results for K0S and Λ. The AGS values are from E896 [1] (centrality 0 − 5 %). The SPS values are from NA49 [20] (centrality 0 − 7 %) and the RHIC values are from STAR [4, 15] (centrality 0 − 5 %). For the multi-strange baryons Ξ and Ω (bottom panel), the SPS results are from NA57 [2] (centrality 0 − 11 %) and the RHIC values are from STAR [15, 21] (centrality 0 − 20 %).

Strange particle production yields at mid-rapidity in central Au+Au and Pb+Pb collisions versus the center of mass energy √sNN. The top panel shows results for K0S and Λ. The AGS values are from E896 [1] (centrality 0 − 5 %). The SPS values are from NA49 [20] (centrality 0 − 7 %) and the RHIC values are from STAR [4, 15] (centrality 0 − 5 %). For the multi-strange baryons Ξ and Ω (bottom panel), the SPS results are from NA57 [2] (centrality 0 − 11 %) and the RHIC values are from STAR [15, 21] (centrality 0 − 20 %).

Anti-baryon to baryon yield ratios for strange baryons versus the center of mass energy √sNN. Λ/Λ is shown in the top panel while the multi-strange baryons are on the bottom panel. The data from AGS are not corrected for the weak decay feed-down from the multistrange baryons while the data from SPS and RHIC are corrected. The lines are the results of a thermal model calculation (see text section IV A). The AGS values are from E896 [1] (centrality 0 − 5 %). The SPS values are from NA49 [20] (centrality 0 − 7 %) and the RHIC values are from STAR [4, 15] (centrality 0 − 5 %). For the multi- strange baryons Ξ and Ω (bottom panel), the SPS results are from NA57 [2] (centrality 0 − 11 %) and the RHIC values are from STAR [15, 21] (centrality 0 − 20 %).

Anti-baryon to baryon yield ratios for strange baryons versus the center of mass energy √sNN. Λ/Λ is shown in the top panel while the multi-strange baryons are on the bottom panel. The data from AGS are not corrected for the weak decay feed-down from the multistrange baryons while the data from SPS and RHIC are corrected. The lines are the results of a thermal model calculation (see text section IV A). The AGS values are from E896 [1] (centrality 0 − 5 %). The SPS values are from NA49 [20] (centrality 0 − 7 %) and the RHIC values are from STAR [4, 15] (centrality 0 − 5 %). For the multi- strange baryons Ξ and Ω (bottom panel), the SPS results are from NA57 [2] (centrality 0 − 11 %) and the RHIC values are from STAR [15, 21] (centrality 0 − 20 %).

Antibaryon-to-baryon yield ratios for strange particles and protons as a function of dNch/dy at √sNN=62.4 and 200 GeV. The p data were extracted from Ref. [14]. The √sNN=200 GeV strange hadron data were extracted from Ref. [15].

Antibaryon-to-baryon yield ratios for strange particles and protons as a function of dNch/dy at √sNN=62.4 and 200 GeV. The p data were extracted from Ref. [14]. The √sNN=200 GeV strange hadron data were extracted from Ref. [15].

Particle-yield ratios as obtained by measurements (black dots) for the most central (0–5%) Au+Au collisions at 62.4 GeV and statistical model predictions (lines). The ratios indicated by the dashed lines (blue) were obtained by using only π, K, and protons, whereas the ratios indicated by the full lines (green) were obtained by also using the hyperons in the fit.

Particle-yield ratios as obtained by measurements (black dots) for the most central (0–5%) Au+Au collisions at 62.4 GeV and statistical model predictions (lines). The ratios indicated by the dashed lines (blue) were obtained by using only π, K, and protons, whereas the ratios indicated by the full lines (green) were obtained by also using the hyperons in the fit.

Chemical freeze-out temperature Tch (a) and strangeness saturation factor γs (b) as a function of the mean number of participants.

Chemical freeze-out temperature Tch (a) and strangeness saturation factor γs (b) as a function of the mean number of participants.

Chemical freeze-out temperature Tch (a) and strangeness saturation factor γs (b) as a function of the mean number of participants.

Chemical freeze-out temperature Tch (a) and strangeness saturation factor γs (b) as a function of the mean number of participants.

Temperature and baryon chemical potential obtained from thermal model fits as a function of √sNN (see Ref. [22]). The dashed lines correspond to the parametrizations given in Ref. [22]. The solid stars show the result for √sNN=62.4 and 200 GeV.

Temperature and baryon chemical potential obtained from thermal model fits as a function of √sNN (see Ref. [22]). The dashed lines correspond to the parametrizations given in Ref. [22]. The solid stars show the result for √sNN=62.4 and 200 GeV.

Temperature and baryon chemical potential obtained from thermal model fits as a function of √sNN (see Ref. [22]). The dashed lines correspond to the parametrizations given in Ref. [22]. The solid stars show the result for √sNN=62.4 and 200 GeV.

Temperature and baryon chemical potential obtained from thermal model fits as a function of √sNN (see Ref. [22]). The dashed lines correspond to the parametrizations given in Ref. [22]. The solid stars show the result for √sNN=62.4 and 200 GeV.

Ratio of baryon (solid symbols) and antibaryon (open symbols) to π+ as a function of dNch/dy for √sNN=62.4 GeV (left) and √sNN=200 GeV (right). The π and p data were extracted from Ref. [14].

Ratio of baryon (solid symbols) and antibaryon (open symbols) to π- as a function of dNch/dy for √sNN=62.4 GeV (left) and √sNN=200 GeV (right). The π and p data were extracted from Ref. [14].

Ratio of baryon (solid symbols) and antibaryon (open symbols) to π+ as a function of dNch/dy for √sNN=62.4 GeV (left) and √sNN=200 GeV (right). The π and p data were extracted from Ref. [14].

Ratio of baryon (solid symbols) and antibaryon (open symbols) to π- as a function of dNch/dy for √sNN=62.4 GeV (left) and √sNN=200 GeV (right). The π and p data were extracted from Ref. [14].

Ratio of baryon (solid symbols) and antibaryon (open symbols) to π+ as a function of dNch/dy for √sNN=62.4 GeV (left) and √sNN=200 GeV (right). The π and p data were extracted from Ref. [14].

Ratio of baryon (solid symbols) and antibaryon (open symbols) to π- as a function of dNch/dy for √sNN=62.4 GeV (left) and √sNN=200 GeV (right). The π and p data were extracted from Ref. [14].

Ratio of baryon (solid symbols) and antibaryon (open symbols) to π+ as a function of dNch/dy for √sNN=62.4 GeV (left) and √sNN=200 GeV (right). The π and p data were extracted from Ref. [14].

Ratio of baryon (solid symbols) and antibaryon (open symbols) to π- as a function of dNch/dy for √sNN=62.4 GeV (left) and √sNN=200 GeV (right). The π and p data were extracted from Ref. [14].

Ratio of baryon (solid symbols) and antibaryon (open symbols) to π− as a function of √sNN. The lines are the results of the thermal model calculation (see text Sec. 4a). The SPS values are from NA49 [20] (centrality 0–7%) and the RHIC values are from STAR [4, 15] (centrality 0–5%). For the multistrange baryons Ξ and Ω (bottom), the SPS results are from NA57 [2] (centrality 0–11%) and the RHIC values are from STAR [15, 21] (centrality 0–20%).

Ratio of baryon (solid symbols) and antibaryon (open symbols) to π− as a function of √sNN. The lines are the results of the thermal model calculation (see text Sec. 4a). The SPS values are from NA49 [20] (centrality 0–7%) and the RHIC values are from STAR [4, 15] (centrality 0–5%). For the multistrange baryons Ξ and Ω (bottom), the SPS results are from NA57 [2] (centrality 0–11%) and the RHIC values are from STAR [15, 21] (centrality 0–20%).

Nuclear modification factor RCP, calculated as the ratio between 0–10% central spectra and 40–80% peripheral spectra, for π, K0S, Λ, and Ξ particles in Au+Au collisions at 62.4 GeV. The π RCP values were extracted from Ref. [10]. The gray band on the right side of the plot shows the uncertainties on the estimation of the number of binary collisions and the gray band on the lower left side indicates the uncertainties on the number of participants.

Nuclear modification factor RCP, calculated as the ratio between 0–10% central spectra and 40–80% peripheral spectra, for π, K0S, Λ, and Ξ particles in Au+Au collisions at 62.4 GeV. The π RCP values were extracted from Ref. [10]. The gray band on the right side of the plot shows the uncertainties on the estimation of the number of binary collisions and the gray band on the lower left side indicates the uncertainties on the number of participants.

Nuclear modification factor RCP, calculated as the ratio between 0–5% central spectra and 40–60% peripheral spectra, for Λ and Ξ particles measured in Au + Au collisions at 62.4 GeV. The gray band corresponds to the equivalent RCP curve for the Λ particles measured in Au+Au collisions at 200 GeV [15].

Nuclear modification factor RCP, calculated as the ratio between 0–5% central spectra and 40–60% peripheral spectra, for Λ and Ξ particles measured in Au + Au collisions at 62.4 GeV. The gray band corresponds to the equivalent RCP curve for the Λ particles measured in Au+Au collisions at 200 GeV [15].

Λ/K0S ratio as a function of transverse momentum for different centrality classes. 0–5% (solid circles), 40–60% (open squares), and 60–80% (solid triangles) in Au+Au collisions at 62.4 GeV.

Λ/K0S ratio as a function of transverse momentum for different centrality classes. 0–5% (solid circles), 40–60% (open squares), and 60–80% (solid triangles) in Au+Au collisions at 62.4 GeV.

Maximum value of the Λ/K0S ratio from Au+Au collisions at 62.4 GeV (solid circles) and 200 GeV (open circles) [11] as a function of ⟨Npart⟩ for different centrality classes. The lowest ⟨Npart⟩ point corresponds to p+p collisions at 200 GeV [44]. The maximum of the Λ––/K0S from Au+Au collisions at 62.4 GeV is shown as solid triangles.

Maximum value of the Λ/K0S ratio from Au+Au collisions at 62.4 GeV (solid circles) and 200 GeV (open circles) [11] as a function of ⟨Npart⟩ for different centrality classes. The lowest ⟨Npart⟩ point corresponds to p+p collisions at 200 GeV [44]. The maximum of the Λ––/K0S from Au+Au collisions at 62.4 GeV is shown as solid triangles.

Longitudinal single-spin asymmetry, AL, for W± events as a function of the leptonic pseudorapidity, $\eta_e$, for 25 < $E^e_T$ < 50 GeV in comparison to theory predictions

Longitudinal single-spin asymmetry, AL, for W± events as a function of the leptonic pseudorapidity, $\eta_e$, for 25 < $E^e_T$ < 50 GeV in comparison to theory predictions

The K$\pi$ pair invariant mass distribution integrated over the $K^{*0}$ $p_T$ for minimum bias Cu+Cu collisions at $\sqrt{s_{NN}}$ =62.4 GeV after mixed-event background subtraction.

The Kπ pair invariant mass distribution for various pT bins (top left) pT = 0.4–0.6 GeV/c in Au+Au collisions at √sNN = 200 GeV after the mixed-event background subtraction.

The Kπ pair invariant mass distribution for various pT bins (top right) pT = 0.6–0.8 GeV/c in Au+Au collisions at √sNN = 62.4 GeV after the mixed-event background subtraction.

The Kπ pair invariant mass distribution for various pT bins (bottom left) pT = 0.8–1.0 GeV/c in Au+Au collisions at √sNN = 200 GeV after the mixed-event background subtraction.

The Kπ pair invariant mass distribution for various pT bins (bottom right) pT = 1.0–1.2 GeV/c in Au+Au collisions at √sNN = 62.4 GeV after the mixed-event background subtraction.

The signal-to-background ratio for $K^{*0}$ measurements as a function of $p_T$ for different collision centrality bins (0-10%, 10-40%, 40-60%, 60-80%) in Au+Au collisions at 200 GeV.

$K^{*0}$ mass as a function of $p_T$ for minimum bias Au+Au collisions at $\sqrt{s_{NN}}$ = 62.4 GeV.

$K^{*0}$ mass as a function of $p_T$ for minimum bias Au+Au collisions at $\sqrt{s_{NN}}$ = 200 GeV.

$K^{*0}$ mass as a function of $p_T$ for minimum bias Cu+Cu collisions at $\sqrt{s_{NN}}$ = 62.4 GeV

$K^{*0}$ mass as a function of $p_T$ for minimum bias Cu+Cu collisions at $\sqrt{s_{NN}}$ = 200 GeV

$K^{*0}$ width as a function of $p_T$ for minimum bias Au+Au collisions at $\sqrt{s_{NN}}$ = 62.4 GeV

$K^{*0}$ width as a function of $p_T$ for minimum bias Au+Au collisions at $\sqrt{s_{NN}}$ = 200 GeV

$K^{*0}$ width as a function of $p_T$ for minimum bias Cu+Cu collisions at $\sqrt{s_{NN}}$ = 62.4 GeV

$K^{*0}$ width as a function of $p_T$ for minimum bias Cu+Cu collisions at $\sqrt{s_{NN}}$ = 200 GeV

The $K^{*0}$ reconstruction efficiency multiplied by the detector acceptance as a function of $p_T$ in Au+Au (|$\eta$| < 0.8) collisions at 200 GeV for different collision centrality bins (0-20% ,20-40% , 40-60%)

The $K^{*0}$ reconstruction efficiency multiplied by the detector acceptance as a function of $p_T$ in Cu+Cu (|$\eta$| < 1.0) collisions at 200 GeV for different collision centrality bins (0-20% ,20-40% , 40-60%)

Mid-rapidity $K^{*0}$ $p_T$ spectra for various collision centrality bins (0-20%, 20-40%, 40-60%, 60-80%) in Au+Au collisions at $\sqrt{s_{NN}}$ = 62.4 GeV

Mid-rapidity $K^{*0}$ $p_T$ spectra for various collision centrality bins (0-20%, 20-40%, 40-60%) in Cu+Cu collisions at $\sqrt{s_{NN}}$ = 62.4 GeV

Mid-rapidity $K^{*0}$ $p_T$ spectra for various collision centrality bins (0-20%, 20-40%, 40-60%, 60-80%) in Au+Au collisions at $\sqrt{s_{NN}}$ = 200 GeV

Mid-rapidity $K^{*0}$ $p_T$ spectra for various collision centrality bins (0-20%, 20-40%, 40-60%) in Cu+Cu collisions at $\sqrt{s_{NN}}$ = 200 GeV

The mid-rapidity yields dN/dy of $K^{*0}$ as a function of the average number of participating nucleons, $⟨N_{part}⟩$, for Au+Au collisions at $\sqrt{s_{NN}}$ = 62.4 GeV

The mid-rapidity yields dN/dy of $K^{*0}$ as a function of the average number of participating nucleons, $⟨N_{part}⟩$, for Cu+Cu collisions at $\sqrt{s_{NN}}$ = 62.4 GeV

The mid-rapidity yields dN/dy of $K^{*0}$ as a function of the average number of participating nucleons, $⟨N_{part}⟩$, for Au+Au collisions at $\sqrt{s_{NN}}$ = 200 GeV

The mid-rapidity yields dN/dy of $K^{*0}$ as a function of the average number of participating nucleons, $⟨N_{part}⟩$, for Cu+Cu collisions at $\sqrt{s_{NN}}$ = 200 GeV

The mid-rapidity $K^{*0}$ $⟨p_T⟩$ as a function $⟨N_{part}⟩$ for Au+Au collisions at $\sqrt{s_{NN}}$ = 62.4 GeV

The mid-rapidity $K^{*0}$ $⟨p_T⟩$ as a function $⟨N_{part}⟩$ for Cu+Cu collisions at $\sqrt{s_{NN}}$ = 62.4 GeV

The mid-rapidity $K^{*0}$ $⟨p_T⟩$ as a function $⟨N_{part}⟩$ for Au+Au collisions at $\sqrt{s_{NN}}$ = 200 GeV

The mid-rapidity $K^{*0}$ $⟨p_T⟩$ as a function $⟨N_{part}⟩$ for Cu+Cu collisions at $\sqrt{s_{NN}}$ = 200 GeV

The mid-rapidity $⟨p_T⟩$ of $\pi$, K, p and $K^{*0}$ as a function of $⟨N_{part}⟩$ for Au+Au collisions at $\sqrt{s_{NN}}$ = 62.4 GeV.

Mid-rapidity $N(K^{*0})/N(K^-)$ ratio for Au+Au at $\sqrt{s_{NN}}$ = 62.4 GeV as a function of $⟨N_{part}⟩$

Mid-rapidity $N(K^{*0})/N(K^-)$ ratio for Cu+Cu at $\sqrt{s_{NN}}$ = 62.4 GeV as a function of $⟨N_{part}⟩$

Mid-rapidity $N(K^{*0})/N(K^-)$ ratio for Au+Au at $\sqrt{s_{NN}}$ = 200 GeV as a function of $⟨N_{part}⟩$

Mid-rapidity $N(K^{*0})/N(K^-)$ ratio for Cu+Cu at $\sqrt{s_{NN}}$ = 200 GeV as a function of $⟨N_{part}⟩$

Mid-rapidity $N(K^{*0})N(K^-)$ in Au+Au collisions divided by $N(K^{*0})N(K^-)$ ratio in p+p collisions at $\sqrt{s_{NN}}$=200 GeV as a function of $⟨N_{part}⟩$.

Mid-rapidity $N(K^{*0})N(K^-)$ in Cu+Cu collisions divided by $N(K^{*0})N(K^-)$ ratio in p+p collisions at $\sqrt{s_{NN}}$=200 GeV as a function of $⟨N_{part}⟩$

Mid-rapidity $N(K^{*0})N(K^-)$ in d+Au collisions divided by $N(K^{*0})N(K^-)$ ratio in d+Au collisions at $\sqrt{s_{NN}}$=200 GeV as a function of $⟨N_{part}⟩$

Mid-rapidity $N(K^{*0})/N(K^-)$ ratio in minimum bias Au+Au collisions as a function of $\sqrt{s_{NN}}.

Mid-rapidity $N(K^{*0})/N(K^-)$ ratio in minimum bias Cu+Cu collisions as a function of $\sqrt{s_{NN}}.

Mid-rapidity $N(K^{*0})/N(K^-)$ ratio in minimum bias p+p collisions as a function of $\sqrt{s_{NN}}.

Mid-rapidity $N(K^{*0})/N(K^-)$ ratio in minimum bias Au+Au collisions as a function of $\sqrt{s_{NN}}.

Mid-rapidity $N(K^{*0})/N(K^-)$ ratio in minimum bias Cu+Cu collisions as a function of $\sqrt{s_{NN}}.

Mid-rapidity $N(K^{*0})/N(K^-)$ ratio in minimum bias p+p collisions as a function of $\sqrt{s_{NN}}.

Mid-rapidity $N(\phi)/N(K^{*0})$ ratio for Au+Au at $\sqrt{s_{NN}}$ = 62.4 GeV as a function of $⟨N_{part}⟩$

Mid-rapidity $N(\phi)/N(K^{*0})$ ratio for Cu+Cu at $\sqrt{s_{NN}}$ = 62.4 GeV as a function of $⟨N_{part}⟩$

Mid-rapidity $N(\phi)/N(K^{*0})$ ratio for Au+Au at $\sqrt{s_{NN}}$ = 200 GeV as a function of $⟨N_{part}⟩$

Mid-rapidity $N(\phi)/N(K^{*0})$ ratio for Cu+Cu at $\sqrt{s_{NN}}$ = 200 GeV as a function of $⟨N_{part}⟩$

Mid-rapidity $[N(\phi)/N(K^{*0})]$ in Au+Au collisions divided by $[N(\phi)/N(K^{*0})]$ ratio in p+p collisions at $\sqrt{s_{NN}}$=200 GeV as a function of $⟨N_{part}⟩$

Mid-rapidity $[N(\phi)/N(K^{*0})]$ in Cu+Cu collisions divided by $[N(\phi)/N(K^{*0})]$ ratio in p+p collisions at $\sqrt{s_{NN}}$=200 GeV as a function of $⟨N_{part}⟩$

Mid-rapidity $[N(\phi)/N(K^{*0})]$ in d+Au collisions divided by $[N(\phi)/N(K^{*0})]$ ratio in p+p collisions at $\sqrt{s_{NN}}$=200 GeV as a function of $⟨N_{part}⟩$

Mid-rapidity $N(\phi)/N(K^{*0})$ ratio in minimum bias Au+Au collisions as a function of $\sqrt{s_{NN}}$.

Mid-rapidity $N(\phi)/N(K^{*0})$ ratio in minimum bias Cu+Cu collisions as a function of $\sqrt{s_{NN}}$.

Mid-rapidity $N(\phi)/N(K^{*0})$ ratio in minimum bias p+p collisions as a function of $\sqrt{s_{NN}}$.

Mid-rapidity $N(\phi)/N(K^{*0})$ ratio in minimum bias Au+Au collisions as a function of $\sqrt{s_{NN}}$.

Mid-rapidity $N(\phi)/N(K^{*0})$ ratio in minimum bias Cu+Cu collisions as a function of $\sqrt{s_{NN}}$.

Mid-rapidity $N(\phi)/N(K^{*0})$ ratio in minimum bias p+p collisions as a function of $\sqrt{s_{NN}}$.

The $K^{*0}$ $v_2$ (Run IV) as a function of $p_T$ in minimum bias Au+Au collisions at $\sqrt{s_{NN}}$ = 200 GeV.

The $K^{*0}$ $v_2$ (Run II) as a function of $p_T$ in minimum bias Au+Au collisions at $\sqrt{s_{NN}}$ = 200 GeV.

The $K^{*0}$ $R_{CP}$ as a function of $p_T$ in Au+Au collisions at 62.4 and 200 GeV compared to the $R_{CP}$ of $K^0_S$ and $\Lambda$ at 200 GeV.

The $K^{*0}$ $R_{CP}$ as a function of $p_T$ in Au+Au collisions at 62.4 and 200 GeV compared to the $R_{CP}$ of $K^0_S$ and $\Lambda$ at 200 GeV.

The $K^{*0}$ $R_{CP}$ as a function of $p_T$ in Au+Au collisions at 62.4 and 200 GeV compared to the $R_{CP}$ of $K^0_S$ and $\Lambda$ at 200 GeV.

The $K^{*0}$ ~$R_{CP}$~ as a function of $p_T$ in Au+Au collisions at 62.4 and 200 GeV compared to the $R_{CP}$ of $K^0_S$ and $\Lambda$ at 200 GeV.

Centrality dependence of moments of $\Delta N_p$ distributions for Au+Au collisions at $\sqrt{s_{NN}}$ = 19.6 GeV.

Centrality dependence of moments of $\Delta N_p$ distributions for Au+Au collisions at $\sqrt{s_{NN}}$ = 19.6 GeV.

Centrality dependence of moments of $\Delta N_p$ distributions for Au+Au collisions at $\sqrt{s_{NN}}$ = 19.6 GeV.

Centrality dependence of moments of $\Delta N_p$ distributions for Au+Au collisions at $\sqrt{s_{NN}}$ = 19.6 GeV.

Centrality dependence of moments of $\Delta N_p$ distributions for Au+Au collisions at $\sqrt{s_{NN}}$ = 62.4 GeV.

Centrality dependence of moments of $\Delta N_p$ distributions for Au+Au collisions at $\sqrt{s_{NN}}$ = 62.4 GeV.

Centrality dependence of moments of $\Delta N_p$ distributions for Au+Au collisions at $\sqrt{s_{NN}}$ = 62.4 GeV.

Centrality dependence of moments of $\Delta N_p$ distributions for Au+Au collisions at $\sqrt{s_{NN}}$ = 62.4 GeV.

Centrality dependence of moments of $\Delta N_p$ distributions for Au+Au collisions at $\sqrt{s_{NN}}$ = 200 GeV.

Centrality dependence of moments of $\Delta N_p$ distributions for Au+Au collisions at $\sqrt{s_{NN}}$ = 200 GeV.

Centrality dependence of moments of $\Delta N_p$ distributions for Au+Au collisions at $\sqrt{s_{NN}}$ = 200 GeV.

Centrality dependence of moments of $\Delta N_p$ distributions for Au+Au collisions at $\sqrt{s_{NN}}$ = 200 GeV.

Centrality dependence of $S\sigma$ for $\Delta N_p$ in Au+Au collisions at $\sqrt{s_{NN}}$ = 19.6 GeV.

Centrality dependence of $\kappa\sigma^2$ for $\Delta N_p$ in Au+Au collisions at $\sqrt{s_{NN}}$ = 19.6 GeV.

Centrality dependence of $S\sigma$ for $\Delta N_p$ in Au+Au collisions at $\sqrt{s_{NN}}$ = 62.4 GeV.

Centrality dependence of $\kappa\sigma^2$ for $\Delta N_p$ in Au+Au collisions at $\sqrt{s_{NN}}$ = 62.4 GeV.

Centrality dependence of $S\sigma$ for $\Delta N_p$ in Au+Au collisions at $\sqrt{s_{NN}}$ = 200 GeV.

Centrality dependence of $\kappa\sigma^2$ for $\Delta N_p$ in Au+Au collisions at $\sqrt{s_{NN}}$ = 200 GeV.

Centrality dependence of $S\sigma$ for $\Delta N_p$ in Au+Au collisions from Lattice QCD Calculations.

Centrality dependence of $S\sigma$ for $\Delta N_p$ in Au+Au collisions at $\sqrt{s_{NN}}$ = 19.6 GeV from Transport Model Calculations.

Centrality dependence of $S\sigma$ for $\Delta N_p$ in Au+Au collisions at $\sqrt{s_{NN}}$ = 62.4 GeV from Transport Model Calculations.

Centrality dependence of $S\sigma$ for $\Delta N_p$ in Au+Au collisions at $\sqrt{s_{NN}}$ = 200 GeV from Transport Model Calculations.

Centrality dependence of $\kappa\sigma^2$ for $\Delta N_p$ in Au+Au collisions at $\sqrt{s_{NN}}$ = 200 GeV from Transport Model Calculations.

Centrality dependence of $\kappa\sigma^2$ for $\Delta N_p$ in Au+Au collisions at $\sqrt{s_{NN}}$ = 200 GeV from Transport Model Calculations for net-baryon.

Centrality dependence of $\kappa\sigma^2$ for $\Delta N_p$ in Au+Au collisions at $\sqrt{s_{NN}}$ = 200 GeV from Transport Model Calculations for net-proton (No Decay).

Centrality dependence of $\kappa\sigma^2$ for $\Delta N_p$ in Au+Au collisions at $\sqrt{s_{NN}}$ = 200 GeV from Transport Model Calculations.

Centrality dependence of $\kappa\sigma^2$ for $\Delta N_p$ in Au+Au collisions at $\sqrt{s_{NN}}$ = 200 GeV from Transport Model Calculations.

Centrality dependence of $\kappa\sigma^2$ for $\Delta N_p$ in Au+Au collisions at $\sqrt{s_{NN}}$ = 200 GeV from Transport Model Calculations.

Centrality dependence of $\kappa\sigma^2$ for $\Delta N_p$ in Au+Au collisions at $\sqrt{s_{NN}}$ = 200 GeV from Transport Model Calculations.

Centrality dependence of $\kappa\sigma^2$ for $\Delta N_p$ in Au+Au collisions at $\sqrt{s_{NN}}$ = 200 GeV from Transport Model Calculations.

Centrality dependence of $\kappa\sigma^2$ for $\Delta N_p$ in Au+Au collisions at $\sqrt{s_{NN}}$ = 200 GeV from Transport Model Calculations.

$\sqrt{s_{NN}}$ dependence of $\kappa\sigma^2$ for net proton distributions measured at RHIC. The results are STAR data.

$\sqrt{s_{NN}}$ dependence of $\kappa\sigma^2$ for net proton distributions measured at RHIC. The results are from UrQMD net-proton.

$\sqrt{s_{NN}}$ dependence of $\kappa\sigma^2$ for net proton distributions measured at RHIC. The results are from AMPT default net-proton.

$\sqrt{s_{NN}}$ dependence of $\kappa\sigma^2$ for net proton distributions measured at RHIC. The results are from AMPT String Melting net-proton.

$\sqrt{s_{NN}}$ dependence of $\kappa\sigma^2$ for net proton distributions measured at RHIC. The results are from Therminator net-proton.

$\sqrt{s_{NN}}$ dependence of $\kappa\sigma^2$ for net proton distributions measured at RHIC. The results are from Hijing net-proton.