Projections of the correlation function C.

Projections of the correlation function C.

Projections of the correlation function C.

Projections of the correlation function C.

Pion HBT radii for the 5% most central collisions.

Pion HBT radii for the 5% most central collisions from the STAR experiment at 200 GeV taken from Adams et al. PR C71(2005)044906.

Ratio of Pion HBT radii for the OUT projection to that for the SIDE projection for the 5% most central collisions.

Ratio of Pion HBT radii for the OUT projection to that for the SIDE projection for the 5% most central collisions from the STAR experiment at 200 GeV taken from Adams et al. PR C71(2005)044906.;.

Pion HBT radii at KT=0.3 for the 5% most central collisions as a function of the central charged multiplicity. Data from other experiments is shown for comparison.

Product of the three HBT radii at KT=0.3 for the 5% most central collisions as a function of the central charged multiplicity. Data from other experiments is shown for comparison.

The decoupling time extracted from R_long(KT) as a function of the central charged multiplicity. Data from other experiments is shown for comparison.

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.

Negatively charged kaon spectra from Cu+Cu collisions 62.4 GeV as a function of pT for different centralities.

Negatively charged proton spectra from Cu+Cu collisions 200 GeV as a function of pT for different centralities.

Negatively charged proton spectra from Cu+Cu collisions 62.4 GeV as a function of pT for different centralities.

Positively charged pion spectra from Cu+Cu collisions 200 GeV as a function of pT for different centralities.

Positively charged pion spectra from Cu+Cu collisions 62.4 GeV as a function of pT for different centralities.

Positively charged kaon spectra from Cu+Cu collisions 200 GeV as a function of pT for different centralities.

Positively charged kaon spectra from Cu+Cu collisions 62.4 GeV as a function of pT for different centralities.

Positively charged proton spectra from Cu+Cu collisions 200 GeV as a function of pT for different centralities.

Positively charged proton spectra from Cu+Cu collisions 62.4 GeV as a function of pT for different centralities.

Mean transverse momentum of negatively charged pions, kaons and protons as a function of charged hadron multiplicity.

Mean transverse momentum of positively charged pions, kaons and protons as a function of charged hadron multiplicity [figure not available in published paper].

Integrated yields of negatively charged pions, kaons and protons as a function of charged hadron multiplicity.

Integrated yields of positively charged pions, kaons and protons as a function of charged hadron multiplicity [figure not available in published paper].

particle Ratios -I (pbar/pi^-, k^-/pi^-) versus multiplicity.

particle ratios -II (p/pi^+, k^+/pi^+) versus multiplicity.

particle ratios -III (p+pbar/pi, k/pi) versus multiplicity.

particle Ratios -IV (pi^-/pi^+, k^-/k^+, pbar/p) versus multiplicity [only pbar/p figure available in paper].

Enhancement factors for negatively charged pions, kaons and protons as a function of Npart [Ref. Phys.Rev.C 81, 044902, 2010]. pp dNdy values are from Ref [Phys.Rev.C 79, 034909, 2009].

The kinetic freeze-out temperature (Tkin) and chemical freeze-out temperature (Tch) versus multiplicity.

flow velocity versus multiplicity.

chemical potentials versus multiplicity.

strangeness suppression factor versus multiplicity.

dNch/deta values for different centrality.

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.

Charged hadron azimuthal anisotropy v2 as a function of pT in 50-60% Cu+Cu collisions at 200 GeV using FTPC flow vectors, and that with subtracting the azimuthal correlations in p+p collisions.

Charged hadron azimuthal anisotropy v2 as a function of pT in 40-50% Cu+Cu collisions at 200 GeV using FTPC flow vectors, and that with subtracting the azimuthal correlations in p+p collisions.

Charged hadron azimuthal anisotropy v2 as a function of pT in 30-40% Cu+Cu collisions at 200 GeV using FTPC flow vectors, and that with subtracting the azimuthal correlations in p+p collisions.

Charged hadron azimuthal anisotropy v2 as a function of pT in 20-30% Cu+Cu collisions at 200 GeV using FTPC flow vectors, and that with subtracting the azimuthal correlations in p+p collisions.

Charged hadron azimuthal anisotropy v2 as a function of pT in 10-20% Cu+Cu collisions at 200 GeV using FTPC flow vectors, and that with subtracting the azimuthal correlations in p+p collisions.

Charged hadron azimuthal anisotropy v2 as a function of pT in 0-10% Cu+Cu collisions at 200 GeV using FTPC flow vectors, and that with subtracting the azimuthal correlations in p+p collisions.

Charged hadron azimuthal anisotropy v2 integrated over pT and eta vs centrality in Cu+Cu collisions at 200 GeV for the various methods.

Charged hadron azimuthal anisotropy v2 integrated over pT and eta vs centrality in Cu+Cu collisions at 62.4 GeV for the various methods.

Charged hadron azimuthal anisotropy v2 as a function of pT in 50-60% Cu+Cu collisions at 62.4 GeV. Data tables for 200 GeV are shown in Fig.4.

Charged hadron azimuthal anisotropy v2 as a function of pT in 40-50% Cu+Cu collisions at 62.4 GeV using FTPC flow vectors. Data tables for 200 GeV are shown in Fig.4.

Charged hadron azimuthal anisotropy v2 as a function of pT in 30-40% Cu+Cu collisions at 62.4 GeV using FTPC flow vectors. Data tables for 200 GeV are shown in Fig.4.

Charged hadron azimuthal anisotropy v2 as a function of pT in 20-30% Cu+Cu collisions at 62.4 GeV using FTPC flow vectors. Data tables for 200 GeV are shown in Fig.4.

Charged hadron azimuthal anisotropy v2 as a function of pT in 10-20% Cu+Cu collisions at 62.4 GeV using FTPC flow vectors. Data tables for 200 GeV are shown in Fig.4.

Charged hadron azimuthal anisotropy v2 as a function of pT in 0-10% Cu+Cu collisions at 62.4 GeV using FTPC flow vectors. Data tables for 200 GeV are shown in Fig.4.

Charged hadron azimuthal anisotropy v2 as a function of pT in 0-60% Cu+Cu collisions at 62.4 GeV using FTPC flow vectors. Data tables for 200 GeV are shown in Fig.3.

K0S azimuthal anisotropy v2 vs pT in 0-60% central Cu+Cu collisions at 200 GeV.

K0S azimuthal anisotropy v2 vs pT in 0-20% central Cu+Cu collisions at 200 GeV.

K0S azimuthal anisotropy v2 vs pT in 20-60% central Cu+Cu collisions at 200 GeV.

Lambda azimuthal anisotropy v2 vs pT in 0-60% central Cu+Cu collisions at 200 GeV.

Lambda azimuthal anisotropy v2 vs pT in 0-20% central Cu+Cu collisions at 200 GeV.

Lambda azimuthal anisotropy v2 vs pT in 20-60% central Cu+Cu collisions at 200 GeV.

Xi azimuthal anisotropy v2 vs pT in 0-60% central Cu+Cu collisions at 200 GeV.

phi azimuthal anisotropy v2 vs pT in 0-60% central Cu+Cu collisions at 200 GeV.

$\langle cos(\phi_{\alpha}+\phi_{\beta}−2\phi_{c})\rangle$ as a function of reference multiplicity for different charge combinations, after corrections for acceptance effects. In the legend the signs indicate the charge of particles $\alpha$, $\beta$, and c. The results shown are for Au+Au collisions at 200 GeV obtained in the Full Field.

$\langle cos(\phi_{\alpha}-\phi_{\beta})\rangle$ as a function of centrality for different charge combinations and FF and RFF configurations. The data points corresponding to different charge and field configurations are slightly shifted in the horizontal direction with respect to each other for clarity. The error bars are statistical. Also shown are model predictions described in Section VII.

A comparison of the correlations obtained by selecting the third particle from the main TPC or from the Forward TPCs. The TPC and FTPC points are shifted horizontally relative to one another for clarity purposes. The error bars are statistical.

A comparison of the correlations obtained by selecting the third particle from the main TPC or from the Forward TPCs after scaling by the flow of the third particle. The shaded areas represent the uncertainty from $v_{2,c}$ scaling (see text for details). The TPC and FTPC points are shifted horizontally relative to one another for clarity purposes. The error bars are statistical.

$\langle cos(\phi_{\alpha}+\phi_{\beta}−2\Psi_{RP})\rangle$ in Au+Au and Cu+Cu collisions at $\sqrt{s_{NN}}$ = 200 GeV calculated using Eq. 7. The error-bars show the statistical errors. The shaded area reflects the uncertainty in the elliptic flow values used in calculations, with lower (in magnitude) limit obtained with elliptic flow from two-particle correlations and upper limit from four-particle cumulants. For details, see Section IV. Thick solid (Au+Au) and dashed (Cu+Cu) lines represent possible non-reaction-plane dependent contribution from many-particle clusters as estimated by HIJING, see Section VII A.

$\langle cos(\phi_{\alpha}+\phi_{\beta}−2\Psi_{RP})\rangle$ in Au+Au and Cu+Cu collisions at $\sqrt{s_{NN}}$ = 62 GeV calculated using Eq. 7. The error-bars show the statistical errors. The shaded area reflects the uncertainty in the elliptic flow values used in calculations. For details, see Section IV. Thick solid (Au+Au) and dashed (Cu+Cu) lines represent possible non-reaction-plane dependent contribution from many-particle clusters as estimated by HIJING, see Section VII A.

Au+Au and Cu+Cu collisions at $\sqrt{s_{NN}}$ = 200 GeV. The correlations are scaled with the number of participants and are plotted as function of centrality. The error-bars indicate the statistical errors. The shaded area reflects the uncertainty in the elliptic flow values used in calculations. For details, see Section IV. Thick solid (Au+Au) and dashed (Cu+Cu) lines represent possible non-reaction-plane dependent contribution from many-particle clusters as estimated by HIJING, see Section VII A.

Au+Au and Cu+Cu collisions at $\sqrt{s_{NN}}$ = 200 GeV. The correlations are scaled with the number of participants and are plotted as function of number of participants. The error-bars indicate the statistical errors. The shaded area reflects the uncertainty in the elliptic flow values used in calculations. For details, see Section IV. Thick solid (Au+Au) and dashed (Cu+Cu) lines represent possible non-reaction-plane dependent contribution from many-particle clusters as estimated by HIJING, see Section VII A.

Au+Au at 200 GeV. The correlations dependence on pseudorapidity separation $\Delta\eta$ = |$\eta_{\alpha}$ − $\eta_{\beta}$| for centrality 30-50%. The shaded band indicates uncertainty associated with $v_{2}$ measurements and has been calculated using two- and four-particle cumulant results as the limits.

Au+Au at 200 GeV. The correlations dependence on pseudorapidity separation $\Delta\eta$ = |$\eta_{\alpha}$ − $\eta_{\beta}$| for centrality 10-30%. The shaded band indicates uncertainty associated with $v_{2}$ measurements and has been calculated using two- and four-particle cumulant results as the limits.

Au+Au at 200 GeV. The correlations dependence on (p$_{t,\alpha}$ + p$_{t,\beta}$)/2 for centrality 30-50%. The shaded band indicates uncertainty associated with $v_{2}$ measurements and has been calculated using two- and four-particle cumulant results as the limits.

Au+Au at 200 GeV. The correlations dependence on (p$_{t,\alpha}$ + p$_{t,\beta}$)/2 for centrality 10-30%. The shaded band indicates uncertainty associated with $v_{2}$ measurements and has been calculated using two- and four-particle cumulant results as the limits.

Au+Au at 200 GeV. The correlations dependence on |p$_{t,\alpha}$ - p$_{t,\beta}$| for centrality 30-50%. The shaded band indicates uncertainty associated with $v_{2}$ measurements and has been calculated using two- and four-particle cumulant results as the limits.

Au+Au at 200 GeV. The correlations dependence on |p$_{t,\alpha}$ - p$_{t,\beta}$| for centrality 10-30%. The shaded band indicates uncertainty associated with $v_{2}$ measurements and has been calculated using two- and four-particle cumulant results as the limits.

Three-particle correlator in Au+Au and Cu+Cu collisions compared to HIJING calculations shown as thick lines. All three particles are taken in the main TPC region, |$\eta$| < 1.0. The correlator has been scaled with number of participants and are plotted versus number of participants. In this representation HIJING results for Au+Au and Cu+Cu collisions coincide in the region of overlap.

$\langle cos(\phi_{\alpha} + \phi_{\beta} − 2\Psi_{RP})\rangle$ calculated for 200 GeV Au+Au events with event generators HIJING (with and without an “elliptic flow afterburner”), UrQMD, and MEVSIM. Blue symbols mark opposite-charge correlations, and red are same-charge. Solid stars represent the values from the data to facilitate comparison. Acceptance cuts of 0.15 < p$_{t}$ < 2 GeV/c and |$\eta$| < 1.0 were used in all cases. For MEVSIM, HIJING, and UrQMD points the true reaction plane from the generated event was used for $\Psi_{RP}$ . Thick solid lighter colored lines represent possible non-reactionplane dependent contribution from many-particle clusters as estimated by HIJING and discussed in Section VII A. Corresponding estimates from UrQMD are about factor of two smaller.

Fig. 2. (Color online.) Event-by-event photon multiplicity distributions (solid circles) for $\mathrm{Au}+\mathrm{Au}$ and $\mathrm{Cu}+\mathrm{Cu}$ at $\sqrt{s_{\mathrm{NN}}}=62.4$ and $200 \mathrm{GeV} .$ The distributions for top $0-5 \%$ central $\mathrm{Au}+$ Au collisions and top $0-10 \%$ central $\mathrm{Cu}+\mathrm{Cu}$ collisions are also shown (open circles). The photon multiplicity distributions for central collisions are observed to be Gaussian (solid line). Only statistical errors are shown. NOTE: For points with invisible error bars, the point size was considered as an absolute upper limit for the uncertainty.

Fig. 2. (Color online.) Event-by-event photon multiplicity distributions (solid circles) for $\mathrm{Au}+\mathrm{Au}$ and $\mathrm{Cu}+\mathrm{Cu}$ at $\sqrt{s_{\mathrm{NN}}}=62.4$ and $200 \mathrm{GeV} .$ The distributions for top $0-5 \%$ central $\mathrm{Au}+$ Au collisions and top $0-10 \%$ central $\mathrm{Cu}+\mathrm{Cu}$ collisions are also shown (open circles). The photon multiplicity distributions for central collisions are observed to be Gaussian (solid line). Only statistical errors are shown. NOTE: For points with invisible error bars, the point size was considered as an absolute upper limit for the uncertainty.

Fig. 2. (Color online.) Event-by-event photon multiplicity distributions (solid circles) for $\mathrm{Au}+\mathrm{Au}$ and $\mathrm{Cu}+\mathrm{Cu}$ at $\sqrt{s_{\mathrm{NN}}}=62.4$ and $200 \mathrm{GeV} .$ The distributions for top $0-5 \%$ central $\mathrm{Au}+$ Au collisions and top $0-10 \%$ central $\mathrm{Cu}+\mathrm{Cu}$ collisions are also shown (open circles). The photon multiplicity distributions for central collisions are observed to be Gaussian (solid line). Only statistical errors are shown. NOTE: For points with invisible error bars, the point size was considered as an absolute upper limit for the uncertainty.

Fig. 3. (Color online.) Photon pseudorapidity distributions for $\mathrm{Au}+\mathrm{Au}$ and $\mathrm{Cu}+\mathrm{Cu}$ at $\sqrt{\mathrm{s}_{\mathrm{NN}}}=62.4$ and $200\mathrm{GeV}$. The results for several centrality classes are shown. The solid curves are results of HIJING simulations for central $(0-5 \%$ for $\mathrm{Au}+\mathrm{Au}$ and $0-10 \%$ for $\mathrm{Cu}+\mathrm{Cu})$ and $30-40 \%$ mid-central collisions. The errors shown are systematic, statistical errors are negligible in comparison. NOTE: For points with invisible error bars, the point size was considered as an absolute upper limit for the uncertainty.

Fig. 3. (Color online.) Photon pseudorapidity distributions for $\mathrm{Au}+\mathrm{Au}$ and $\mathrm{Cu}+\mathrm{Cu}$ at $\sqrt{\mathrm{s}_{\mathrm{NN}}}=62.4$ and $200\mathrm{GeV}$. The results for several centrality classes are shown. The solid curves are results of HIJING simulations for central $(0-5 \%$ for $\mathrm{Au}+\mathrm{Au}$ and $0-10 \%$ for $\mathrm{Cu}+\mathrm{Cu})$ and $30-40 \%$ mid-central collisions. The errors shown are systematic, statistical errors are negligible in comparison. NOTE: For points with invisible error bars, the point size was considered as an absolute upper limit for the uncertainty.

Fig. 3. (Color online.) Photon pseudorapidity distributions for $\mathrm{Au}+\mathrm{Au}$ and $\mathrm{Cu}+\mathrm{Cu}$ at $\sqrt{\mathrm{s}_{\mathrm{NN}}}=62.4$ and $200\mathrm{GeV}$. The results for several centrality classes are shown. The solid curves are results of HIJING simulations for central $(0-5 \%$ for $\mathrm{Au}+\mathrm{Au}$ and $0-10 \%$ for $\mathrm{Cu}+\mathrm{Cu})$ and $30-40 \%$ mid-central collisions. The errors shown are systematic, statistical errors are negligible in comparison. NOTE: For points with invisible error bars, the point size was considered as an absolute upper limit for the uncertainty.

Fig. 3. (Color online.) Photon pseudorapidity distributions for $\mathrm{Au}+\mathrm{Au}$ and $\mathrm{Cu}+\mathrm{Cu}$ at $\sqrt{\mathrm{s}_{\mathrm{NN}}}=62.4$ and $200\mathrm{GeV}$. The results for several classes are shown. The solid curves are results of HIJING simulations for central $(0-5 \%$ for $\mathrm{Au}+\mathrm{Au}$ and $0-10 \%$ for $\mathrm{Cu}+\mathrm{Cu})$ and $30-40 \%$ mid-central collisions. The errors shown are systematic, statistical errors are negligible in comparison. NOTE: For points with invisible error bars, the point size was considered as an absolute upper limit for the uncertainty.

Fig. 4. (Color online.) Top panel: The number of photons divided by $\left\langle N_{\text{part}}\right\rangle / 2$ as a function of average number of participating nucleons for $\mathrm{Au}+\mathrm{Au}$ and $\mathrm{Cu}+\mathrm{Cu}$ at $\sqrt{\mathrm{s} _{\mathrm{NN}}}=62.4$ and $200 \mathrm{GeV}$ for $-3.7<\eta<-2.3 .$ Errors shown are systematic only and include those for $\left\langle N_{\text {part }}\right\rangle .$ Results from HIJING are shown as lines (solid for $\mathrm{Au}+\mathrm{Au}$ and dashed for $\mathrm{Cu}+\mathrm{Cu}$ ). Bottom panel: Same as above, for both photons and charged particles from PHOBOS [10]. NOTE: For points with invisible error bars, the point size was considered as an absolute upper limit for the uncertainty.

Fig. 4. (Color online.) Top panel: The number of photons divided by $\left\langle N_{\text{part}}\right\rangle / 2$ as a function of average number of participating nucleons for $\mathrm{Au}+\mathrm{Au}$ and $\mathrm{Cu}+\mathrm{Cu}$ at $\sqrt{\mathrm{s} _{\mathrm{NN}}}=62.4$ and $200 \mathrm{GeV}$ for $-3.7<\eta<-2.3 .$ Errors shown are systematic only and include those for $\left\langle N_{\text {part }}\right\rangle .$ Results from HIJING are shown as lines (solid for $\mathrm{Au}+\mathrm{Au}$ and dashed for $\mathrm{Cu}+\mathrm{Cu}$ ). Bottom panel: Same as above, for both photons and charged particles from PHOBOS [10]. NOTE: For points with invisible error bars, the point size was considered as an absolute upper limit for the uncertainty.

Fig. 5. (Color online.) Photon pseudorapidity distributions normalized by the average number of participating nucleon pairs for different collision centralities are plotted as a function of pseudorapidity shifted by the beam rapidity (-5.36 for $200 \mathrm{GeV}$ and -4.19 for $62.4 \mathrm{GeV}$ ) for $\mathrm{Au}+\mathrm{Au}$ and $\mathrm{Cu}+\mathrm{Cu}$ collisions at $\sqrt{s_{\mathrm{NN}}}=62.4$ and $200 \mathrm{GeV}$. Errors are systematic only, statistical errors are negligible in comparison. For clarity of presentation, results for only four centralities are shown. The $\mathrm{Cu}+$ Cu data are shifted by 0.1 unit in $\eta-y_{\text {beam }} .$ The solid line is a second order polynomial fit to the data (see text for details).

Fig. 5. (Color online.) Photon pseudorapidity distributions normalized by the average number of participating nucleon pairs for different collision centralities are plotted as a function of pseudorapidity shifted by the beam rapidity (-5.36 for $200 \mathrm{GeV}$ and -4.19 for $62.4 \mathrm{GeV}$ ) for $\mathrm{Au}+\mathrm{Au}$ and $\mathrm{Cu}+\mathrm{Cu}$ collisions at $\sqrt{s_{\mathrm{NN}}}=62.4$ and $200 \mathrm{GeV}$. Errors are systematic only, statistical errors are negligible in comparison. For clarity of presentation, results for only four centralities are shown. The $\mathrm{Cu}+$ Cu data are shifted by 0.1 unit in $\eta-y_{\text {beam }} .$ The solid line is a second order polynomial fit to the data (see text for details).

HBT parameters for all centralities in 62.4 GeV Cu+Cu.

HBT parameters for all centralities in 200 GeV Cu+Cu.

HBT parameters for 200 GeV Au+Au, 0-5% and Cu+Cu, 0-10%

HBT parameters for 200 GeV and 62.4 GeV Cu+Cu, 0-10% and 62.4 GeV Au+Au, 0-5%

HBT radii for radii ratios in 200 GeV and 62.4 GeV Cu+Cu, 0-10% and 200 GeV & 62.4 GeV Au+Au, 0-5%. Note that on Figure 8 in the paper the radii ratios are shown while here we list the radii values.

Energy dependence of the pi- HBT Radii for Volume estimate (Rside*Rside*Rlong) in central Au+Au, Pb+Pb, and Pb+Au collisions (AGS,SPS and RHIC) at midrapidity and k_T ~ 0.2-0.3 GeV/c. Note that on Figure 9 in the paper Rside*Rside*Rlong is shown while here we list the Rside and Rlong values.

HBT Radii for Volume estimate(Rside*Rside*Rlong & Rout*Rout*Rlong) comparisons for 200 GeV and 62.4 GeV Au+Au with Npart and dNch/deta. Note that on Figure 10 in the paper specific combinations of the HBT radii are shown while here we list the radii values.

HBT Radii for Volume estimate(Rside*Rside*Rlong & Rout*Rout*Rlong) comparisons for 200 GeV and 62.4 GeV Au+Au & Cu+Cu with dNch/deta. Note that on Figure 11 in the paper specific combinations of the HBT radii are shown while here we list the radii values.

HBT radii comparisons for 200 GeV and 62.4 GeV Au+Au & Cu+Cu with dNch/deta1/3

(Color online) Measured dynamical $K/\pi$ fluctuations in terms of $\nu_{dyn,K\pi}$ for 62.4 and 200 GeV Au+Au compared with $\sigma^{2}_{dyn}$ from central Pb+Pb collisions at 6.3, 7.6, 8.8, 12.3, and 17.3 GeV from NA49 [7]. Statistical errors are shown for the STAR data. Statistical and systematic errors are shown for the NA49 results. The solid line corresponds to a fit to the STAR data of the form $c + d/(dN/d\eta)$.

(Color online) The dN/dη scaled dynamical $K/\pi$ fluctuations for summed charges (stars), same signs (circles), and opposite signs (squares) as a function of $dN/d\eta$. The errors shown are statistical. The open and filled symbols refer to Au+Au collisions at 62.4 GeV and 200 GeV respectively. The dash-dot, dotted, and dashed lines represents HIJING calculations for summed charges, same signs, and opposite signs respectively.

Upper panels. $N_{\scriptsize{\mbox{part}}}$ scaled ($R^{N_{\scriptsize{\mbox{part}}}}_{AA}$) nuclear modification factors as a function of $p_{T}$ of $\phi$ mesons for $0-10\%$ and $20-30\%$ $Cu+Cu$ and $Au+Au$ collisions at $\sqrt{s_{NN}}=200$ GeV. Lower panel. Same as above for $N_{\scriptsize{\mbox{bin}}}$ scaled ($R^{N_{\scriptsize{\mbox{bin}}}}_{AA}$) nuclear modification factor. The error bars represent the statistical and systematic errors added in quadrature. The shaded band in upper panel around 1 at $p_{T}=4.5-5.5$ GeV/$c$ in the right side reflects the uncertainty in $N_{\scriptsize{\mbox{part}}}$ and that on the lower panel for $N_{\scriptsize{\mbox{bin}}}$ calculation for central $Au+Au$ collisions. The respective uncertainties for central $Cu+Cu$ collisions are of similar order.

Upper panels. $N_{\scriptsize{\mbox{part}}}$ scaled ($R^{N_{\scriptsize{\mbox{part}}}}_{AA}$) nuclear modification factors as a function of $p_{T}$ of $\phi$ mesons for $0-10\%$ and $20-30\%$ $Cu+Cu$ and $Au+Au$ collisions at $\sqrt{s_{NN}}=200$ GeV. Lower panel. Same as above for $N_{\scriptsize{\mbox{bin}}}$ scaled ($R^{N_{\scriptsize{\mbox{bin}}}}_{AA}$) nuclear modification factor. The error bars represent the statistical and systematic errors added in quadrature. The shaded band in upper panel around 1 at $p_{T}=4.5-5.5$ GeV/$c$ in the right side reflects the uncertainty in $N_{\scriptsize{\mbox{part}}}$ and that on the lower panel for $N_{\scriptsize{\mbox{bin}}}$ calculation for central $Au+Au$ collisions. The respective uncertainties for central $Cu+Cu$ collisions are of similar order.

Upper panels. $N_{\scriptsize{\mbox{part}}}$ scaled ($R^{N_{\scriptsize{\mbox{part}}}}_{AA}$) nuclear modification factors as a function of $p_{T}$ of $\phi$ mesons for $0-10\%$ and $20-30\%$ $Cu+Cu$ and $Au+Au$ collisions at $\sqrt{s_{NN}}=200$ GeV. Lower panel. Same as above for $N_{\scriptsize{\mbox{bin}}}$ scaled ($R^{N_{\scriptsize{\mbox{bin}}}}_{AA}$) nuclear modification factor. The error bars represent the statistical and systematic errors added in quadrature. The shaded band in upper panel around 1 at $p_{T}=4.5-5.5$ GeV/$c$ in the right side reflects the uncertainty in $N_{\scriptsize{\mbox{part}}}$ and that on the lower panel for $N_{\scriptsize{\mbox{bin}}}$ calculation for central $Au+Au$ collisions. The respective uncertainties for central $Cu+Cu$ collisions are of similar order.

Upper panel. The ratio of the yields of $K^{−}$, $\phi$, $\bar{\Lambda}$ and $\Xi+\bar{\Xi}$ normalized to $N_{\scriptsize{\mbox{part}}}$ in nucleus-nucleus collisions to corresponding yields in inelastic $p+p$ collisions as a function of $N_{\scriptsize{\mbox{part}}}$ at 200 GeV. Lower panel. Same as above for $\phi$ mesons in $Cu+Cu$ collisions at 200 and 62.4 GeV. The $p+p$ collision data at 200 GeV are from STAR [5] and at 62.4 GeV from ISR [29]. The error bars shown represent the statistical and systematic errors added in quadrature.

Upper panel. The ratio of the yields of $K^{−}$, $\phi$, $\bar{\Lambda}$ and $\Xi+\bar{\Xi}$ normalized to $N_{\scriptsize{\mbox{part}}}$ in nucleus-nucleus collisions to corresponding yields in inelastic $p+p$ collisions as a function of $N_{\scriptsize{\mbox{part}}}$ at 200 GeV. Lower panel. Same as above for $\phi$ mesons in $Cu+Cu$ collisions at 200 and 62.4 GeV. The $p+p$ collision data at 200 GeV are from STAR [5] and at 62.4 GeV from ISR [29]. The error bars shown represent the statistical and systematic errors added in quadrature.

Upper panels: same-event (full points) and mixed-event (solid line) $K^{+}K^{-}$ invariant mass distributions at 0.6 < $p_{T}$ < 1.4 GeV/c in p + p 200 GeV collisions (a), 0.8 < $p_{T}$ < 1.2 GeV/c in Au + Au 62.4 GeV collisions (60–80%) (c), and 0.8 < $p_{T}$ < 1.2 GeV/c in Au + Au 200 GeV collisions (0–10%) (e). Lower panels: the corresponding $\phi$ meson mass peaks after subtracting the background. Dashed curves show a Breit-Wigner + linear background function fit in (b), (d). In (f), both linear and quadratic backgrounds are shown as dashed and dot-dashed lines, respectively.

Upper panels: same-event (full points) and mixed-event (solid line) $K^{+}K^{-}$ invariant mass distributions at 0.6 < $p_{T}$ < 1.4 GeV/c in p + p 200 GeV collisions (a), 0.8 < $p_{T}$ < 1.2 GeV/c in Au + Au 62.4 GeV collisions (60–80%) (c), and 0.8 < $p_{T}$ < 1.2 GeV/c in Au + Au 200 GeV collisions (0–10%) (e). Lower panels: the corresponding $\phi$ meson mass peaks after subtracting the background. Dashed curves show a Breit-Wigner + linear background function fit in (b), (d). In (f), both linear and quadratic backgrounds are shown as dashed and dot-dashed lines, respectively.

Upper panels: same-event (full points) and mixed-event (solid line) $K^{+}K^{-}$ invariant mass distributions at 0.6 < $p_{T}$ < 1.4 GeV/c in p + p 200 GeV collisions (a), 0.8 < $p_{T}$ < 1.2 GeV/c in Au + Au 62.4 GeV collisions (60–80%) (c), and 0.8 < $p_{T}$ < 1.2 GeV/c in Au + Au 200 GeV collisions (0–10%) (e). Lower panels: the corresponding $\phi$ meson mass peaks after subtracting the background. Dashed curves show a Breit-Wigner + linear background function fit in (b), (d). In (f), both linear and quadratic backgrounds are shown as dashed and dot-dashed lines, respectively.

Reconstruction efficiency including acceptance of $\phi$ meson as a function of $p_{T}$ in several centrality bins of Au + Au, d + Au, and p + p 200 GeV collisions.

Reconstruction efficiency including acceptance of $\phi$ meson as a function of $p_{T}$ in several centrality bins of Au + Au, d + Au, and p + p 200 GeV collisions.

Reconstruction efficiency including acceptance of $\phi$ meson as a function of $p_{T}$ in several centrality bins of Au + Au, d + Au, and p + p 200 GeV collisions.

Event plane $\Phi_{2}$ resolution as a function of centrality in Au + Au 200 GeV collisions, where the vertical axis starts from 0.3 for clarity.

$\phi−\Phi_{2}$ distribution for $\phi$ meson at 1.5 < $p_{T}$ < 2.0 GeV/c in Au + Au collisions (0–80$\%$) at 200 GeV. The line is the fitting result. Error bars are statistical only.

<cos2($\phi_{K^{+}K^{-}}$-$\Phi_{2}$)> (full red points) and <sin2($\phi_{K^{+}K^{-}}$-$\Phi_{2}$)> (open blue points) as a function of $m_{inv}$ of $K^{+}K^{-}$ pairs at 0.5 < $p_{T}$ < 1.0 GeV/c in Au + Au 200 GeV collisions (0–80\%), where the solid curve is the result of fitting by Eq. (10). The arrow shows the position of the $\phi$ invariant mass peak. The dashed line shows the zero horizontal line.

Masses and widths (FWHMs) of $\phi$ as a function of $p_{T}$ in p + p 200 GeV (NSD), d + Au 200 GeV (0–20$\%$), Au + Au 62.4 GeV (0–20$\%$), and Au + Au 200 GeV (0–5$\%$) collisions, with the corresponding $PDG$ values.

Masses and widths (FWHMs) of $\phi$ as a function of $p_{T}$ in p + p 200 GeV (NSD), d + Au 200 GeV (0–20$\%$), Au + Au 62.4 GeV (0–20$\%$), and Au + Au 200 GeV (0–5$\%$) collisions, with the corresponding $PDG$ values.

Masses and widths (FWHMs) of $\phi$ as a function of $p_{T}$ in p + p 200 GeV (NSD), d + Au 200 GeV (0–20$\%$), Au + Au 62.4 GeV (0–20$\%$), and Au + Au 200 GeV (0–5$\%$) collisions, with the corresponding $PDG$ values.

Masses and widths (FWHMs) of $\phi$ as a function of $p_{T}$ in p + p 200 GeV (NSD), d + Au 200 GeV (0–20$\%$), Au + Au 62.4 GeV (0–20$\%$), and Au + Au 200 GeV (0–5$\%$) collisions, with the corresponding $PDG$ values.

Masses and widths (FWHMs) of $\phi$ as a function of $p_{T}$ in p + p 200 GeV (NSD), d + Au 200 GeV (0–20$\%$), Au + Au 62.4 GeV (0–20$\%$), and Au + Au 200 GeV (0–5$\%$) collisions, with the corresponding $PDG$ values.

Masses and widths (FWHMs) of $\phi$ as a function of $p_{T}$ in p + p 200 GeV (NSD), d + Au 200 GeV (0–20$\%$), Au + Au 62.4 GeV (0–20$\%$), and Au + Au 200 GeV (0–5$\%$) collisions, with the corresponding $PDG$ values.

Masses and widths (FWHMs) of $\phi$ as a function of $p_{T}$ in p + p 200 GeV (NSD), d + Au 200 GeV (0–20$\%$), Au + Au 62.4 GeV (0–20$\%$), and Au + Au 200 GeV (0–5$\%$) collisions, with the corresponding $PDG$ values.

Masses and widths (FWHMs) of $\phi$ as a function of $p_{T}$ in p + p 200 GeV (NSD), d + Au 200 GeV (0–20$\%$), Au + Au 62.4 GeV (0–20$\%$), and Au + Au 200 GeV (0–5$\%$) collisions, with the corresponding $PDG$ values.

Masses and widths (FWHMs) of $\phi$ as a function of $p_{T}$ in p + p 200 GeV (NSD), d + Au 200 GeV (0–20$\%$), Au + Au 62.4 GeV (0–20$\%$), and Au + Au 200 GeV (0–5$\%$) collisions, with the corresponding $PDG$ values.

Masses and widths (FWHMs) of $\phi$ as a function of $p_{T}$ in p + p 200 GeV (NSD), d + Au 200 GeV (0–20$\%$), Au + Au 62.4 GeV (0–20$\%$), and Au + Au 200 GeV (0–5$\%$) collisions, with the corresponding $PDG$ values.

Masses and widths (FWHMs) of $\phi$ as a function of $p_{T}$ in p + p 200 GeV (NSD), d + Au 200 GeV (0–20$\%$), Au + Au 62.4 GeV (0–20$\%$), and Au + Au 200 GeV (0–5$\%$) collisions, with the corresponding $PDG$ values.

Masses and widths (FWHMs) of $\phi$ as a function of $p_{T}$ in p + p 200 GeV (NSD), d + Au 200 GeV (0–20$\%$), Au + Au 62.4 GeV (0–20$\%$), and Au + Au 200 GeV (0–5$\%$) collisions, with the corresponding $PDG$ values.

Masses and widths (FWHMs) of $\phi$ as a function of $p_{T}$ in p + p 200 GeV (NSD), d + Au 200 GeV (0–20$\%$), Au + Au 62.4 GeV (0–20$\%$), and Au + Au 200 GeV (0–5$\%$) collisions, with the corresponding $PDG$ values.

Invariant mass distributions of $\phi$ meson at 0.6 < $p_{T}$ < 1.0 GeV/c in p + p 200 GeV (NSD) and Au + Au 200 GeV (0–5$\%$) collisions. Solid symbols: experimental data. Open symbols: MC simulation. Curves are the results of a Breit-Wigner function fit. Note: Two sets of MC data are shown for Au + Au 200 GeV, and see text for details.

Invariant mass distributions of $\phi$ meson at 0.6 < $p_{T}$ < 1.0 GeV/c in p + p 200 GeV (NSD) and Au + Au 200 GeV (0–5$\%$) collisions. Solid symbols: experimental data. Open symbols: MC simulation. Curves are the results of a Breit-Wigner function fit. Note: Two sets of MC data are shown for Au + Au 200 GeV, and see text for details.

Invariant mass distributions of $\phi$ meson at 0.6 < $p_{T}$ < 1.0 GeV/c in p + p 200 GeV (NSD) and Au + Au 200 GeV (0–5$\%$) collisions. Solid symbols: experimental data. Open symbols: MC simulation. Curves are the results of a Breit-Wigner function fit. Note: Two sets of MC data are shown for Au + Au 200 GeV, and see text for details.

Invariant mass distributions of $\phi$ meson at 0.6 < $p_{T}$ < 1.0 GeV/c in p + p 200 GeV (NSD) and Au + Au 200 GeV (0–5$\%$) collisions. Solid symbols: experimental data. Open symbols: MC simulation. Curves are the results of a Breit-Wigner function fit. Note: Two sets of MC data are shown for Au + Au 200 GeV, and see text for details.

Invariant mass distributions of $\phi$ meson at 0.6 < $p_{T}$ < 1.0 GeV/c in p + p 200 GeV (NSD) and Au + Au 200 GeV (0–5$\%$) collisions. Solid symbols: experimental data. Open symbols: MC simulation. Curves are the results of a Breit-Wigner function fit. Note: Two sets of MC data are shown for Au + Au 200 GeV, and see text for details.

$\phi$ meson transverse mass distributions for different collision systems and different energies. For clarity, distributions for some centrality bins have been scaled by factors indicated in the figure. Curves represent the exponential (solid) and Levy (dashed) function fits to the distributions. Error bars are statistical only. Note that a scale factor of 1.09 is applied to the $\phi$ meson spectra for Au + Au collisions at 200 GeV in Ref. [50] to correct for the kaon identification efficiency effect missed previously.

$\phi$ meson transverse mass distributions for different collision systems and different energies. For clarity, distributions for some centrality bins have been scaled by factors indicated in the figure. Curves represent the exponential (solid) and Levy (dashed) function fits to the distributions. Error bars are statistical only. Note that a scale factor of 1.09 is applied to the $\phi$ meson spectra for Au + Au collisions at 200 GeV in Ref. [50] to correct for the kaon identification efficiency effect missed previously.

$\phi$ meson transverse mass distributions for different collision systems and different energies. For clarity, distributions for some centrality bins have been scaled by factors indicated in the figure. Curves represent the exponential (solid) and Levy (dashed) function fits to the distributions. Error bars are statistical only. Note that a scale factor of 1.09 is applied to the $\phi$ meson spectra for Au + Au collisions at 200 GeV in Ref. [50] to correct for the kaon identification efficiency effect missed previously.

$\phi$ meson transverse mass distributions for different collision systems and different energies. For clarity, distributions for some centrality bins have been scaled by factors indicated in the figure. Curves represent the exponential (solid) and Levy (dashed) function fits to the distributions. Error bars are statistical only. Note that a scale factor of 1.09 is applied to the $\phi$ meson spectra for Au + Au collisions at 200 GeV in Ref. [50] to correct for the kaon identification efficiency effect missed previously.

$\phi$ meson transverse mass distributions for different collision systems and different energies. For clarity, distributions for some centrality bins have been scaled by factors indicated in the figure. Curves represent the exponential (solid) and Levy (dashed) function fits to the distributions. Error bars are statistical only. Note that a scale factor of 1.09 is applied to the $\phi$ meson spectra for Au + Au collisions at 200 GeV in Ref. [50] to correct for the kaon identification efficiency effect missed previously.

$\phi$ meson transverse mass distributions for different collision systems and different energies. For clarity, distributions for some centrality bins have been scaled by factors indicated in the figure. Curves represent the exponential (solid) and Levy (dashed) function fits to the distributions. Error bars are statistical only. Note that a scale factor of 1.09 is applied to the $\phi$ meson spectra for Au + Au collisions at 200 GeV in Ref. [50] to correct for the kaon identification efficiency effect missed previously.

$\phi$ meson transverse mass distributions for different collision systems and different energies. For clarity, distributions for some centrality bins have been scaled by factors indicated in the figure. Curves represent the exponential (solid) and Levy (dashed) function fits to the distributions. Error bars are statistical only. Note that a scale factor of 1.09 is applied to the $\phi$ meson spectra for Au + Au collisions at 200 GeV in Ref. [50] to correct for the kaon identification efficiency effect missed previously.

Comparison of transverse momentum spectra shape among different 200 GeV collision systems: Au + Au (0–5$\%$), d + Au (0–20$\%$), and p + p (inelastic). The spectra are normalized by $N_{bin}$ (top panel) and $N_{part}/2$ (bottom panel).

Comparison of transverse momentum spectra shape among different 200 GeV collision systems: Au + Au (0–5$\%$), d + Au (0–20$\%$), and p + p (inelastic). The spectra are normalized by $N_{bin}$ (top panel) and $N_{part}/2$ (bottom panel).

Comparison of transverse momentum spectra shape among different 200 GeV collision systems: Au + Au (0–5$\%$), d + Au (0–20$\%$), and p + p (inelastic). The spectra are normalized by $N_{bin}$ (top panel) and $N_{part}/2$ (bottom panel).

Comparison of transverse momentum spectra shape among different 200 GeV collision systems: Au + Au (0–5$\%$), d + Au (0–20$\%$), and p + p (inelastic). The spectra are normalized by $N_{bin}$ (top panel) and $N_{part}/2$ (bottom panel).

Comparison of transverse momentum spectra shape among different 200 GeV collision systems: Au + Au (0–5$\%$), d + Au (0–20$\%$), and p + p (inelastic). The spectra are normalized by $N_{bin}$ (top panel) and $N_{part}/2$ (bottom panel).

Comparison of transverse momentum spectra shape among different 200 GeV collision systems: Au + Au (0–5$\%$), d + Au (0–20$\%$), and p + p (inelastic). The spectra are normalized by $N_{bin}$ (top panel) and $N_{part}/2$ (bottom panel).

$N_{part}$ dependence of $(dN/dy)/0.5N_{part}$ in five different collision systems: Au + Au 62.4, 130, 200 GeV; p + p 200 GeV (inelastic); and d + Au 200 GeV. Statistical and systematic errors are included.

$N_{part}$ dependence of $(dN/dy)/0.5N_{part}$ in five different collision systems: Au + Au 62.4, 130, 200 GeV; p + p 200 GeV (inelastic); and d + Au 200 GeV. Statistical and systematic errors are included.

$N_{part}$ dependence of $(dN/dy)/0.5N_{part}$ in five different collision systems: Au + Au 62.4, 130, 200 GeV; p + p 200 GeV (inelastic); and d + Au 200 GeV. Statistical and systematic errors are included.

$N_{part}$ dependence of $(dN/dy)/0.5N_{part}$ in five different collision systems: Au + Au 62.4, 130, 200 GeV; p + p 200 GeV (inelastic); and d + Au 200 GeV. Statistical and systematic errors are included.

$N_{part}$ dependence of $(dN/dy)/0.5N_{part}$ in five different collision systems: Au + Au 62.4, 130, 200 GeV; p + p 200 GeV (inelastic); and d + Au 200 GeV. Statistical and systematic errors are included.

Top panel: $N_{part}$ dependence of $<p_{T}>_{\phi}$ in different collision systems. Bottom panel: Hadron mass dependence of $<p_{T}>$ in central Au + Au collisions at 62.4 and 200 GeV. The band and curve show two hydrodynamic model calculations for central Au + Au collisions at 200 GeV. Note: Hadron masses for the Au + Au 62.4 GeV data are shifted slightly in the $x$-axis direction for clarity, and systematic errors are included for $\phi$.

Top panel: $N_{part}$ dependence of $<p_{T}>_{\phi}$ in different collision systems. Bottom panel: Hadron mass dependence of $<p_{T}>$ in central Au + Au collisions at 62.4 and 200 GeV. The band and curve show two hydrodynamic model calculations for central Au + Au collisions at 200 GeV. Note: Hadron masses for the Au + Au 62.4 GeV data are shifted slightly in the $x$-axis direction for clarity, and systematic errors are included for $\phi$.

Top panel: $N_{part}$ dependence of $<p_{T}>_{\phi}$ in different collision systems. Bottom panel: Hadron mass dependence of $<p_{T}>$ in central Au + Au collisions at 62.4 and 200 GeV. The band and curve show two hydrodynamic model calculations for central Au + Au collisions at 200 GeV. Note: Hadron masses for the Au + Au 62.4 GeV data are shifted slightly in the $x$-axis direction for clarity, and systematic errors are included for $\phi$.

Top panel: $N_{part}$ dependence of $<p_{T}>_{\phi}$ in different collision systems. Bottom panel: Hadron mass dependence of $<p_{T}>$ in central Au + Au collisions at 62.4 and 200 GeV. The band and curve show two hydrodynamic model calculations for central Au + Au collisions at 200 GeV. Note: Hadron masses for the Au + Au 62.4 GeV data are shifted slightly in the $x$-axis direction for clarity, and systematic errors are included for $\phi$.

Top panel: $N_{part}$ dependence of $<p_{T}>_{\phi}$ in different collision systems. Bottom panel: Hadron mass dependence of $<p_{T}>$ in central Au + Au collisions at 62.4 and 200 GeV. The band and curve show two hydrodynamic model calculations for central Au + Au collisions at 200 GeV. Note: Hadron masses for the Au + Au 62.4 GeV data are shifted slightly in the $x$-axis direction for clarity, and systematic errors are included for $\phi$.

Top panel: $N_{part}$ dependence of $<p_{T}>_{\phi}$ in different collision systems. Bottom panel: Hadron mass dependence of $<p_{T}>$ in central Au + Au collisions at 62.4 and 200 GeV. The band and curve show two hydrodynamic model calculations for central Au + Au collisions at 200 GeV. Note: Hadron masses for the Au + Au 62.4 GeV data are shifted slightly in the $x$-axis direction for clarity, and systematic errors are included for $\phi$.

Top panel: Energy dependence of the ratio $\phi/\pi^{-}$ in A + A (full points) and p + p (open points) collisions. Stars are data from the STARexperiment at RHIC.Bottom panel: $N_{part}$ dependence of ratio $\phi/\pi^{-}$ in different collision systems. Systematic errors are included for the STAR data points.

Top panel: Energy dependence of the ratio $\phi/\pi^{-}$ in A + A (full points) and p + p (open points) collisions. Stars are data from the STARexperiment at RHIC.Bottom panel: $N_{part}$ dependence of ratio $\phi/\pi^{-}$ in different collision systems. Systematic errors are included for the STAR data points.

Top panel: Energy dependence of the ratio $\phi/\pi^{-}$ in A + A (full points) and p + p (open points) collisions. Stars are data from the STARexperiment at RHIC.Bottom panel: $N_{part}$ dependence of ratio $\phi/\pi^{-}$ in different collision systems. Systematic errors are included for the STAR data points.

Top panel: Energy dependence of the ratio $\phi/\pi^{-}$ in A + A (full points) and p + p (open points) collisions. Stars are data from the STARexperiment at RHIC.Bottom panel: $N_{part}$ dependence of ratio $\phi/\pi^{-}$ in different collision systems. Systematic errors are included for the STAR data points.

Top panel: Energy dependence of the ratio $\phi/\pi^{-}$ in A + A (full points) and p + p (open points) collisions. Stars are data from the STARexperiment at RHIC.Bottom panel: $N_{part}$ dependence of ratio $\phi/\pi^{-}$ in different collision systems. Systematic errors are included for the STAR data points.

Top panel: Energy dependence of the ratio $\phi/\pi^{-}$ in A + A (full points) and p + p (open points) collisions. Stars are data from the STARexperiment at RHIC.Bottom panel: $N_{part}$ dependence of ratio $\phi/\pi^{-}$ in different collision systems. Systematic errors are included for the STAR data points.

Top panel: Energy dependence of the ratio $\phi/\pi^{-}$ in A + A (full points) and p + p (open points) collisions. Stars are data from the STARexperiment at RHIC.Bottom panel: $N_{part}$ dependence of ratio $\phi/\pi^{-}$ in different collision systems. Systematic errors are included for the STAR data points.

Top panel: Energy dependence of ratio $\phi/K^{-}$ in $A + A$ and elementary ($e + e$ and $p + p$) collisions. Stars are data from STAR experiments at RHIC. Bottom panel: $N_{part}$ dependence of ratio $\phi/K^{-}$ in different collision systems. The dashed line shows results from $UrQMD$ model calculations. Systematic errors are included for the STAR data points.

Top panel: Energy dependence of ratio $\phi/K^{-}$ in $A + A$ and elementary ($e + e$ and $p + p$) collisions. Stars are data from STAR experiments at RHIC. Bottom panel: $N_{part}$ dependence of ratio $\phi/K^{-}$ in different collision systems. The dashed line shows results from $UrQMD$ model calculations. Systematic errors are included for the STAR data points.

Top panel: Energy dependence of ratio $\phi/K^{-}$ in $A + A$ and elementary ($e + e$ and $p + p$) collisions. Stars are data from STAR experiments at RHIC. Bottom panel: $N_{part}$ dependence of ratio $\phi/K^{-}$ in different collision systems. The dashed line shows results from $UrQMD$ model calculations. Systematic errors are included for the STAR data points.

Top panel: Energy dependence of ratio $\phi/K^{-}$ in $A + A$ and elementary ($e + e$ and $p + p$) collisions. Stars are data from STAR experiments at RHIC. Bottom panel: $N_{part}$ dependence of ratio $\phi/K^{-}$ in different collision systems. The dashed line shows results from $UrQMD$ model calculations. Systematic errors are included for the STAR data points.

Top panel: Energy dependence of ratio $\phi/K^{-}$ in $A + A$ and elementary ($e + e$ and $p + p$) collisions. Stars are data from STAR experiments at RHIC. Bottom panel: $N_{part}$ dependence of ratio $\phi/K^{-}$ in different collision systems. The dashed line shows results from $UrQMD$ model calculations. Systematic errors are included for the STAR data points.

Top panel: Energy dependence of ratio $\phi/K^{-}$ in $A + A$ and elementary ($e + e$ and $p + p$) collisions. Stars are data from STAR experiments at RHIC. Bottom panel: $N_{part}$ dependence of ratio $\phi/K^{-}$ in different collision systems. The dashed line shows results from $UrQMD$ model calculations. Systematic errors are included for the STAR data points.

Top panel: Energy dependence of ratio $\phi/K^{-}$ in $A + A$ and elementary ($e + e$ and $p + p$) collisions. Stars are data from STAR experiments at RHIC. Bottom panel: $N_{part}$ dependence of ratio $\phi/K^{-}$ in different collision systems. The dashed line shows results from $UrQMD$ model calculations. Systematic errors are included for the STAR data points.

Top panel: Energy dependence of ratio $\phi/K^{-}$ in $A + A$ and elementary ($e + e$ and $p + p$) collisions. Stars are data from STAR experiments at RHIC. Bottom panel: $N_{part}$ dependence of ratio $\phi/K^{-}$ in different collision systems. The dashed line shows results from $UrQMD$ model calculations. Systematic errors are included for the STAR data points.

The $\Omega/\phi$ ratio vs $p_{T}$ for three centrality bins in $\sqrt{s_{NN}} = 200 GeV$ Au + Au collisions, where the data points for 40–60$\%$ are shifted slightly for clarity. As shown in the legend, the lines represent results from Hwa and Yang [96], Ko et al. [97], and, for SCF, Refs. [98,99].

The $\Omega/\phi$ ratio vs $p_{T}$ for three centrality bins in $\sqrt{s_{NN}} = 200 GeV$ Au + Au collisions, where the data points for 40–60$\%$ are shifted slightly for clarity. As shown in the legend, the lines represent results from Hwa and Yang [96], Ko et al. [97], and, for SCF, Refs. [98,99].

$p_{T}$ dependence of the nuclear modification factor $R_{cp}$ in Au + Au 200 GeV collisions. The top and bottom panels present $R_{cp}$ from midperipheral and most-peripheral collisions, respectively. See legend for symbol and line designations. The rectangular bands showthe uncertainties of binary and participant scalings. Statistical and systematic errors are included.

$p_{T}$ dependence of the nuclear modification factor $R_{cp}$ in Au + Au 200 GeV collisions. The top and bottom panels present $R_{cp}$ from midperipheral and most-peripheral collisions, respectively. See legend for symbol and line designations. The rectangular bands showthe uncertainties of binary and participant scalings. Statistical and systematic errors are included.

$p_{T}$ dependence of the nuclear modification factor $R_{cp}$ in Au + Au 62.4 and 200 GeV and d + Au 200 GeV collisions, where rectangular bands represent the uncertainties of binary and participant scalings (see legend). Statistical and systematic errors are included.

$p_{T}$ dependence of the nuclear modification factor $R_{cp}$ in Au + Au 62.4 and 200 GeV and d + Au 200 GeV collisions, where rectangular bands represent the uncertainties of binary and participant scalings (see legend). Statistical and systematic errors are included.

$p_{T}$ dependence of the nuclear modification factor $R_{cp}$ in Au + Au 62.4 and 200 GeV and d + Au 200 GeV collisions, where rectangular bands represent the uncertainties of binary and participant scalings (see legend). Statistical and systematic errors are included.

$p_{T}$ dependence of the nuclear modification factor $R_{AB}$ for $\phi$ in Au + Au 200 GeV and d + Au 200 GeV collisions. For comparison, data points for $\pi^{+}+\pi^{-}$ in d + Au 200 GeV and $p+\bar{p}$ in d + Au 200 GeV are also shown (see legend). Rectangular bands show the uncertainties of binary (solid line) and participant (dot-dash line) scalings. Systematic errors are included for $\phi$, $\pi^{+}+\pi^{-}$, and $p+\bar{p}$.

$p_{T}$ dependence of the nuclear modification factor $R_{AB}$ for $\phi$ in Au + Au 200 GeV and d + Au 200 GeV collisions. For comparison, data points for $\pi^{+}+\pi^{-}$ in d + Au 200 GeV and $p+\bar{p}$ in d + Au 200 GeV are also shown (see legend). Rectangular bands show the uncertainties of binary (solid line) and participant (dot-dash line) scalings. Systematic errors are included for $\phi$, $\pi^{+}+\pi^{-}$, and $p+\bar{p}$.

$p_{T}$ dependence of the elliptic flow $v_{2}$ of $\phi$, $\Lambda$, and $K_{s}^{0}$ in Au + Au collisions (0–80$\%$) at 200 GeV. Data points for $\phi$ are from the reaction plane method (full up-triangles) and invariant mass method (full circles), where data points from the reaction plane method are shifted slightly along the $x$ axis for clarity. Vertical error bars represent statistical errors, while the square bands represent systematic uncertainties. The magenta curved band represents the $v_{2}$ of $\phi$ meson from the AMPT model with a string melting mechanism [116]. The dash and dot curves represent parametrizations inspired by number-of-quark scaling ideas from Ref. [117] for NQ = 2 and NQ = 3, respectively.

Elliptic flow $v_{2}$ as a function of $p_{T}$ , $v_{2}$($p_{T}$), for the $\phi$ meson from different centralities. The vertical error bars represent the statistical errors, while the square bands represent the systematic uncertainties. For clarity, data points of 10–40$\%$ are shifted in the $p_{T}$ direction slightly.

Elliptic flow $v_{2}$ as a function of $p_{T}$ , $v_{2}$($p_{T}$), for the $\phi$ meson from different centralities. The vertical error bars represent the statistical errors, while the square bands represent the systematic uncertainties. For clarity, data points of 10–40$\%$ are shifted in the $p_{T}$ direction slightly.

(Color online) Dynamical net charge fluctuations, $\nu_{+−,dyn}$, of particles produced with pseudorapidity $|\eta|$ < 0.5 scaled by (a) the multiplicity, $dN_{ch}/d\eta$. The dashed line corresponds to charge conservation effect and the solid line to the prediction for a resonance gas, (b) the number of participants, and (c) the number of binary collisions.

(Color online) Dynamical net charge fluctuations, $\nu_{+−,dyn}$, of particles produced with pseudorapidity $|\eta|$ < 0.5 scaled by (a) the multiplicity, $dN_{ch}/d\eta$. The dashed line corresponds to charge conservation effect and the solid line to the prediction for a resonance gas, (b) the number of participants, and (c) the number of binary collisions.

(Color online) Dynamical fluctuations $\nu_{+−,dyn}$, normalized to their value for |$\eta$| < 1, as function of the integrated pseudorapidity range. (a) data for $Au + Au$ collisions at $\sqrt{s_{NN}}$ = 62.4, 200 GeV (0-5%) along with data for $Cu + Cu$ collisions at $\sqrt{s_{NN}}$ = 62.4, 200 GeV (0-10%), are compared to inclusive $p + p$ data at $\sqrt{s}$ = 200 GeV, and (b) data for $Au + Au$ collisions at $\sqrt{s_{NN}}$ = 62.4, 200 GeV (30-40%) along with data for $Cu + Cu$ collisions at $\sqrt{s_{NN}}$ = 62.4, 200 GeV (0-10%), are compared to inclusive $p + p$ collision data at $\sqrt{}s$ = 200 GeV.

(Color online) Dynamical fluctuations $\nu_{+−,dyn}$, normalized to their value for |$\eta$| < 1, as function of the integrated pseudorapidity range. (a) data for $Au + Au$ collisions at $\sqrt{s_{NN}}$ = 62.4, 200 GeV (0-5%) along with data for $Cu + Cu$ collisions at $\sqrt{s_{NN}}$ = 62.4, 200 GeV (0-10%), are compared to inclusive $p + p$ data at $\sqrt{s}$ = 200 GeV, and (b) data for $Au + Au$ collisions at $\sqrt{s_{NN}}$ = 62.4, 200 GeV (30-40%) along with data for $Cu + Cu$ collisions at $\sqrt{s_{NN}}$ = 62.4, 200 GeV (0-10%), are compared to inclusive $p + p$ collision data at $\sqrt{}s$ = 200 GeV.

(Color online) Dynamical fluctuations $\nu_{+−,dyn}$, as a function of the integrated azimuthal range $\phi$ for selected collision centralities for (a) $Au + Au$ collisions at $\sqrt{s_{NN}}$ = 200 GeV, and (b) $Cu + Cu$ collisions at $\sqrt{s_{NN}}$ = 200 GeV.

(Color online) Dynamical fluctuations $\nu_{+−,dyn}$, as a function of the integrated azimuthal range $\phi$ for selected collision centralities for (a) $Au + Au$ collisions at $\sqrt{s_{NN}}$ = 200 GeV, and (b) $Cu + Cu$ collisions at $\sqrt{s_{NN}}$ = 200 GeV.

(Color online) Dynamical fluctuations $\nu_{+−,dyn}$, as a function of the integrated azimuthal range $\phi$ for similar number of participating nucleons for $Au + Au$ and $Cu + Cu$ $\sqrt{s_{NN}}$ = 200 GeV.

Charged particle $v_1(\eta)$ for 40-80 % centrality in Au+Au collisions at 200 GeV.

Charged particle $v_1(p_t)$ in 200 GeV Au+Au for three centralities and two pseudorapidity ranges.

Charged particle $v_1(\eta)$ for mid-central (30-60%) Au+Au and Cu+Cu collisions at 200 GeV and 62.4 GeV.

Charged particle $v_1$ versus centrality for Au+Au and Cu+Cu collisions at 200 GeV and 62.4 GeV in two pseudorapidity ranges.

Charged particle $v_1$ versus $\eta-y_{beam}$, for 30 − 60% Au+Au collisions at 200 GeV.

Charged particle $v_1$ versus $\eta-y_{beam}$, for 30 − 60% Cu+Cu collisions at 200 GeV.

Charged particle $v_1$ versus $\eta-y_{beam}$, for 30 − 60% Au+Au collisions at 62.4 GeV.

Charged particle $v_1$ versus $\eta-y_{beam}$, for 30 − 60% Cu+Cu collisions at 62.4 GeV.

(Color online) Global polarization of $\Lambda$–hyperons as a function of centrality given as fraction of the total inelastic hadronic cross section. Filled circles show the results for Au+Au collisions at $\sqrt{s_{NN}}$=200 GeV (centrality region 20-70%) and open squares indicate the results for Au+Au collisions at $\sqrt{s_{NN}}$=62.4 GeV (centrality region 0-80%). Only statistical uncertainties are shown.

(Color online) Global polarization of $\overline{\Lambda}$–hyperons as a function of $\overline{\Lambda}$ transverse momentum $p^{\overline{\Lambda}}_{t}$. Filled circles show the results for Au+Au collisions at $\sqrt{s_{NN}}$=200 GeV (centrality region 20-70%) and open squares indicate the results for Au+Au collisions at $\sqrt{s_{NN}}$=62.4 GeV (centrality region 0-80%). Only statistical uncertainties are shown.

(Color online) Global polarization of $\overline{\Lambda}$–hyperons as a function of $\overline{\Lambda}$ pseudorapidity $\eta^{\overline{\Lambda}}$. Filled circles show the results for Au+Au collisions at $\sqrt{s_{NN}}$=200 GeV (centrality region 20-70%). A constant line fit to these data points yields $P_{\Lambda}=(1.8\pm 10.8)\times 10^{-3}$ with $\chi^{2}/ndf=5.5/10$. Open squares show the results for Au+Au collisions as $\sqrt{s_{NN}}$=62.4 GeV (centrality region 0-80%). A constant line fit gives $P_{\Lambda}=(-17.6\pm 11.1)\times 10^{-3}$ with $\chi^{2}/ndf=8.0/10$. Only statistical uncertainties are shown.

(Color online) Global polarization of $\overline{\Lambda}$–hyperons as a function of as a function of centrality. Filled circles show the results for Au+Au collisions at $\sqrt{s_{NN}}$=200 GeV (centrality region 20-70%) and open squares indicate the result for Au+Au collisions at $\sqrt{s_{NN}}$=62.4 GeV (centrality region 0-80%). Only statistical uncertainties are shown.

(Color online) Acceptance correction function $A_{0}(p^{H}_{t}, \eta^{H})$ defined in (10) as a function of $\Lambda$ (filled circles) and $\overline{\Lambda}$ (open squares) transverse momentum (top) and pseudorapidity (bottom). The deviation of this function from unity affects the overall scale of the measured global polarization according to Eq. 9. See the text for details and discussions on $A_{0}$ $p^{H}_{t}$ and $\eta^{H}$ dependence.

(Color online) Acceptance correction function $A_{0}(p^{H}_{t}, \eta^{H})$ defined in (10) as a function of $\Lambda$ (filled circles) and $\overline{\Lambda}$ (open squares) transverse momentum (top) and pseudorapidity (bottom). The deviation of this function from unity affects the overall scale of the measured global polarization according to Eq. 9. See the text for details and discussions on $A_{0}$ $p^{H}_{t}$ and $\eta^{H}$ dependence.

(Color online) Function $A_{2}(p^{H}_{t}, \eta^{H})$ defined in (11) as a function of $\Lambda$ (filled circles) and $\overline{\Lambda}$ (open squares) transverse momentum (top) and pseudorapidity (bottom). The deviation of this function from zero defines the contribution to the observable (3) from $P^{(2)}_{H} (p^{H}_{t}, \eta^{H})$ in the expansion (6).

(Color online) Function $A_{2}(p^{H}_{t}, \eta^{H})$ defined in (11) as a function of $\Lambda$ (filled circles) and $\overline{\Lambda}$ (open squares) transverse momentum (top) and pseudorapidity (bottom). The deviation of this function from zero defines the contribution to the observable (3) from $P^{(2)}_{H} (p^{H}_{t}, \eta^{H})$ in the expansion (6).

Centrality and pT dependence of RCP for pions and protons for Au+Au at sqrt(sNN) = 62.4 GeV. For studying the energy dependence, the corresponding RCP for central 0-12% Au+Au at sqrt(sNN) = 200 GeV are shown. RAA for pions at 62.4 GeV (0-10%) and 200 GeV (0-12%) are also shown.

The proton (anti-proton) over pion ratios versus pT at sqrt(sNN) = 62.4 and 200 GeV for central and peripheral collisions.

Minimum-bias (0–80% of the collision cross section) v2(pT ) for identified hadrons at |η| < 1 from Au+Au collisions at √sNN = 62.4 GeV. To facilitate comparisons between panels, v2 values for inclusive charged hadrons are displayed in each panel. The error bars on the data points represent statistical uncertainties. Systematic uncertainties for the identified particles are shown as shaded bands around v2 = 0.

Minimum-bias (0–80% of the collision cross section) v2(pT ) for identified hadrons at |η| < 1 from Au+Au collisions at √sNN = 62.4 GeV. To facilitate comparisons between panels, v2 values for inclusive charged hadrons are displayed in each panel. The error bars on the data points represent statistical uncertainties. Systematic uncertainties for the identified particles are shown as shaded bands around v2 = 0.

Minimum-bias (0–80% of the collision cross section) v2(pT ) for identified hadrons at |η| < 1 from Au+Au collisions at √sNN = 62.4 GeV. To facilitate comparisons between panels, v2 values for inclusive charged hadrons are displayed in each panel. The error bars on the data points represent statistical uncertainties. Systematic uncertainties for the identified particles are shown as shaded bands around v2 = 0.

Minimum-bias (0–80% of the collision cross section) v2(pT ) for identified hadrons at |η| < 1 from Au+Au collisions at √sNN = 62.4 GeV. To facilitate comparisons between panels, v2 values for inclusive charged hadrons are displayed in each panel. The error bars on the data points represent statistical uncertainties. Systematic uncertainties for the identified particles are shown as shaded bands around v2 = 0.

Minimum-bias (0–80% of the collision cross section) v2(pT ) for identified hadrons at |η| < 1 from Au+Au collisions at √sNN = 62.4 GeV. To facilitate comparisons between panels, v2 values for inclusive charged hadrons are displayed in each panel. The error bars on the data points represent statistical uncertainties. Systematic uncertainties for the identified particles are shown as shaded bands around v2 = 0.

Minimum-bias (0–80% of the collision cross section) v2(pT ) for identified hadrons at |η| < 1 from Au+Au collisions at √sNN = 62.4 GeV. To facilitate comparisons between panels, v2 values for inclusive charged hadrons are displayed in each panel. The error bars on the data points represent statistical uncertainties. Systematic uncertainties for the identified particles are shown as shaded bands around v2 = 0.

Minimum-bias (0–80% of the collision cross section) v2(pT ) for identified hadrons at |η| < 1 from Au+Au collisions at √sNN = 62.4 GeV. To facilitate comparisons between panels, v2 values for inclusive charged hadrons are displayed in each panel. The error bars on the data points represent statistical uncertainties. Systematic uncertainties for the identified particles are shown as shaded bands around v2 = 0.

Minimum-bias (0–80% of the collision cross section) v2(pT ) for identified hadrons at |η| < 1 from Au+Au collisions at √sNN = 62.4 GeV. To facilitate comparisons between panels, v2 values for inclusive charged hadrons are displayed in each panel. The error bars on the data points represent statistical uncertainties. Systematic uncertainties for the identified particles are shown as shaded bands around v2 = 0.

Minimum-bias (0–80% of the collision cross section) v2(pT ) for identified hadrons at |η| < 1 from Au+Au collisions at √sNN = 62.4 GeV. To facilitate comparisons between panels, v2 values for inclusive charged hadrons are displayed in each panel. The error bars on the data points represent statistical uncertainties. Systematic uncertainties for the identified particles are shown as shaded bands around v2 = 0.

Minimum-bias (0–80% of the collision cross section) v2(pT ) for identified hadrons at |η| < 1 from Au+Au collisions at √sNN = 62.4 GeV. To facilitate comparisons between panels, v2 values for inclusive charged hadrons are displayed in each panel. The error bars on the data points represent statistical uncertainties. Systematic uncertainties for the identified particles are shown as shaded bands around v2 = 0.

Minimum-bias (0–80% of the collision cross section) v2(pT ) for identified hadrons at |η| < 1 from Au+Au collisions at √sNN = 62.4 GeV. To facilitate comparisons between panels, v2 values for inclusive charged hadrons are displayed in each panel. The error bars on the data points represent statistical uncertainties. Systematic uncertainties for the identified particles are shown as shaded bands around v2 = 0.

Minimum-bias (0–80% of the collision cross section) v2(pT ) for identified hadrons at |η| < 1 from Au+Au collisions at √sNN = 62.4 GeV. To facilitate comparisons between panels, v2 values for inclusive charged hadrons are displayed in each panel. The error bars on the data points represent statistical uncertainties. Systematic uncertainties for the identified particles are shown as shaded bands around v2 = 0.

Minimum-bias (0–80% of the collision cross section) v2(pT ) for identified hadrons at |η| < 1 from Au+Au collisions at √sNN = 62.4 GeV. To facilitate comparisons between panels, v2 values for inclusive charged hadrons are displayed in each panel. The error bars on the data points represent statistical uncertainties. Systematic uncertainties for the identified particles are shown as shaded bands around v2 = 0.

The unidentified charged hadron, charged pion, $K^0_S$, charged kaon, proton and Λ+Λ v2 as a function of pT for 10%–40%, 0%–10% and 40%–80% of the Au+Au interaction cross section at √sNN = 62.4 GeV. Weak-decay feed-down errors are included in the error bars on the data points while non-flow and tracking error uncertainties are plotted as bands around v2 = 0, which apply to all identified particles. The errors are asymmetric and the portion of the error band above (below) zero represents the negative (positive) error.

The unidentified charged hadron, charged pion, $K^0_S$, charged kaon, proton and Λ+Λ v2 as a function of pT for 10%–40%, 0%–10% and 40%–80% of the Au+Au interaction cross section at √sNN = 62.4 GeV. Weak-decay feed-down errors are included in the error bars on the data points while non-flow and tracking error uncertainties are plotted as bands around v2 = 0, which apply to all identified particles. The errors are asymmetric and the portion of the error band above (below) zero represents the negative (positive) error.

The unidentified charged hadron, charged pion, $K^0_S$, charged kaon, proton and Λ+Λ v2 as a function of pT for 10%–40%, 0%–10% and 40%–80% of the Au+Au interaction cross section at √sNN = 62.4 GeV. Weak-decay feed-down errors are included in the error bars on the data points while non-flow and tracking error uncertainties are plotted as bands around v2 = 0, which apply to all identified particles. The errors are asymmetric and the portion of the error band above (below) zero represents the negative (positive) error.

The unidentified charged hadron, charged pion, $K^0_S$, charged kaon, proton and Λ+Λ v2 as a function of pT for 10%–40%, 0%–10% and 40%–80% of the Au+Au interaction cross section at √sNN = 62.4 GeV. Weak-decay feed-down errors are included in the error bars on the data points while non-flow and tracking error uncertainties are plotted as bands around v2 = 0, which apply to all identified particles. The errors are asymmetric and the portion of the error band above (below) zero represents the negative (positive) error.

The unidentified charged hadron, charged pion, $K^0_S$, charged kaon, proton and Λ+Λ v2 as a function of pT for 10%–40%, 0%–10% and 40%–80% of the Au+Au interaction cross section at √sNN = 62.4 GeV. Weak-decay feed-down errors are included in the error bars on the data points while non-flow and tracking error uncertainties are plotted as bands around v2 = 0, which apply to all identified particles. The errors are asymmetric and the portion of the error band above (below) zero represents the negative (positive) error.

The unidentified charged hadron, charged pion, $K^0_S$, charged kaon, proton and Λ+Λ v2 as a function of pT for 10%–40%, 0%–10% and 40%–80% of the Au+Au interaction cross section at √sNN = 62.4 GeV. Weak-decay feed-down errors are included in the error bars on the data points while non-flow and tracking error uncertainties are plotted as bands around v2 = 0, which apply to all identified particles. The errors are asymmetric and the portion of the error band above (below) zero represents the negative (positive) error.

The unidentified charged hadron, charged pion, $K^0_S$, charged kaon, proton and Λ+Λ v2 as a function of pT for 10%–40%, 0%–10% and 40%–80% of the Au+Au interaction cross section at √sNN = 62.4 GeV. Weak-decay feed-down errors are included in the error bars on the data points while non-flow and tracking error uncertainties are plotted as bands around v2 = 0, which apply to all identified particles. The errors are asymmetric and the portion of the error band above (below) zero represents the negative (positive) error.

The unidentified charged hadron, charged pion, $K^0_S$, charged kaon, proton and Λ+Λ v2 as a function of pT for 10%–40%, 0%–10% and 40%–80% of the Au+Au interaction cross section at √sNN = 62.4 GeV. Weak-decay feed-down errors are included in the error bars on the data points while non-flow and tracking error uncertainties are plotted as bands around v2 = 0, which apply to all identified particles. The errors are asymmetric and the portion of the error band above (below) zero represents the negative (positive) error.

The unidentified charged hadron, charged pion, $K^0_S$, charged kaon, proton and Λ+Λ v2 as a function of pT for 10%–40%, 0%–10% and 40%–80% of the Au+Au interaction cross section at √sNN = 62.4 GeV. Weak-decay feed-down errors are included in the error bars on the data points while non-flow and tracking error uncertainties are plotted as bands around v2 = 0, which apply to all identified particles. The errors are asymmetric and the portion of the error band above (below) zero represents the negative (positive) error.

The unidentified charged hadron, charged pion, $K^0_S$, charged kaon, proton and Λ+Λ v2 as a function of pT for 10%–40%, 0%–10% and 40%–80% of the Au+Au interaction cross section at √sNN = 62.4 GeV. Weak-decay feed-down errors are included in the error bars on the data points while non-flow and tracking error uncertainties are plotted as bands around v2 = 0, which apply to all identified particles. The errors are asymmetric and the portion of the error band above (below) zero represents the negative (positive) error.

The unidentified charged hadron, charged pion, $K^0_S$, charged kaon, proton and Λ+Λ v2 as a function of pT for 10%–40%, 0%–10% and 40%–80% of the Au+Au interaction cross section at √sNN = 62.4 GeV. Weak-decay feed-down errors are included in the error bars on the data points while non-flow and tracking error uncertainties are plotted as bands around v2 = 0, which apply to all identified particles. The errors are asymmetric and the portion of the error band above (below) zero represents the negative (positive) error.

The unidentified charged hadron, charged pion, $K^0_S$, charged kaon, proton and Λ+Λ v2 as a function of pT for 10%–40%, 0%–10% and 40%–80% of the Au+Au interaction cross section at √sNN = 62.4 GeV. Weak-decay feed-down errors are included in the error bars on the data points while non-flow and tracking error uncertainties are plotted as bands around v2 = 0, which apply to all identified particles. The errors are asymmetric and the portion of the error band above (below) zero represents the negative (positive) error.

The unidentified charged hadron, charged pion, $K^0_S$, charged kaon, proton and Λ+Λ v2 as a function of pT for 10%–40%, 0%–10% and 40%–80% of the Au+Au interaction cross section at √sNN = 62.4 GeV. Weak-decay feed-down errors are included in the error bars on the data points while non-flow and tracking error uncertainties are plotted as bands around v2 = 0, which apply to all identified particles. The errors are asymmetric and the portion of the error band above (below) zero represents the negative (positive) error.

The unidentified charged hadron, charged pion, $K^0_S$, charged kaon, proton and Λ+Λ v2 as a function of pT for 10%–40%, 0%–10% and 40%–80% of the Au+Au interaction cross section at √sNN = 62.4 GeV. Weak-decay feed-down errors are included in the error bars on the data points while non-flow and tracking error uncertainties are plotted as bands around v2 = 0, which apply to all identified particles. The errors are asymmetric and the portion of the error band above (below) zero represents the negative (positive) error.

The unidentified charged hadron, charged pion, $K^0_S$, charged kaon, proton and Λ+Λ v2 as a function of pT for 10%–40%, 0%–10% and 40%–80% of the Au+Au interaction cross section at √sNN = 62.4 GeV. Weak-decay feed-down errors are included in the error bars on the data points while non-flow and tracking error uncertainties are plotted as bands around v2 = 0, which apply to all identified particles. The errors are asymmetric and the portion of the error band above (below) zero represents the negative (positive) error.

The unidentified charged hadron, charged pion, $K^0_S$, charged kaon, proton and Λ+Λ v2 as a function of pT for 10%–40%, 0%–10% and 40%–80% of the Au+Au interaction cross section at √sNN = 62.4 GeV. Weak-decay feed-down errors are included in the error bars on the data points while non-flow and tracking error uncertainties are plotted as bands around v2 = 0, which apply to all identified particles. The errors are asymmetric and the portion of the error band above (below) zero represents the negative (positive) error.

The unidentified charged hadron, charged pion, $K^0_S$, charged kaon, proton and Λ+Λ v2 as a function of pT for 10%–40%, 0%–10% and 40%–80% of the Au+Au interaction cross section at √sNN = 62.4 GeV. Weak-decay feed-down errors are included in the error bars on the data points while non-flow and tracking error uncertainties are plotted as bands around v2 = 0, which apply to all identified particles. The errors are asymmetric and the portion of the error band above (below) zero represents the negative (positive) error.

The unidentified charged hadron, charged pion, $K^0_S$, charged kaon, proton and Λ+Λ v2 as a function of pT for 10%–40%, 0%–10% and 40%–80% of the Au+Au interaction cross section at √sNN = 62.4 GeV. Weak-decay feed-down errors are included in the error bars on the data points while non-flow and tracking error uncertainties are plotted as bands around v2 = 0, which apply to all identified particles. The errors are asymmetric and the portion of the error band above (below) zero represents the negative (positive) error.

The unidentified charged hadron, charged pion, $K^0_S$, charged kaon, proton and Λ+Λ v2 as a function of pT for 10%–40%, 0%–10% and 40%–80% of the Au+Au interaction cross section at √sNN = 62.4 GeV. Weak-decay feed-down errors are included in the error bars on the data points while non-flow and tracking error uncertainties are plotted as bands around v2 = 0, which apply to all identified particles. The errors are asymmetric and the portion of the error band above (below) zero represents the negative (positive) error.

The unidentified charged hadron, charged pion, $K^0_S$, charged kaon, proton and Λ+Λ v2 as a function of pT for 10%–40%, 0%–10% and 40%–80% of the Au+Au interaction cross section at √sNN = 62.4 GeV. Weak-decay feed-down errors are included in the error bars on the data points while non-flow and tracking error uncertainties are plotted as bands around v2 = 0, which apply to all identified particles. The errors are asymmetric and the portion of the error band above (below) zero represents the negative (positive) error.

The unidentified charged hadron, charged pion, $K^0_S$, charged kaon, proton and Λ+Λ v2 as a function of pT for 10%–40%, 0%–10% and 40%–80% of the Au+Au interaction cross section at √sNN = 62.4 GeV. Weak-decay feed-down errors are included in the error bars on the data points while non-flow and tracking error uncertainties are plotted as bands around v2 = 0, which apply to all identified particles. The errors are asymmetric and the portion of the error band above (below) zero represents the negative (positive) error.

The unidentified charged hadron, charged pion, $K^0_S$, charged kaon, proton and Λ+Λ v2 as a function of pT for 10%–40%, 0%–10% and 40%–80% of the Au+Au interaction cross section at √sNN = 62.4 GeV. Weak-decay feed-down errors are included in the error bars on the data points while non-flow and tracking error uncertainties are plotted as bands around v2 = 0, which apply to all identified particles. The errors are asymmetric and the portion of the error band above (below) zero represents the negative (positive) error.

The unidentified charged hadron, charged pion, $K^0_S$, charged kaon, proton and Λ+Λ v2 as a function of pT for 10%–40%, 0%–10% and 40%–80% of the Au+Au interaction cross section at √sNN = 62.4 GeV. Weak-decay feed-down errors are included in the error bars on the data points while non-flow and tracking error uncertainties are plotted as bands around v2 = 0, which apply to all identified particles. The errors are asymmetric and the portion of the error band above (below) zero represents the negative (positive) error.

The unidentified charged hadron, charged pion, $K^0_S$, charged kaon, proton and Λ+Λ v2 as a function of pT for 10%–40%, 0%–10% and 40%–80% of the Au+Au interaction cross section at √sNN = 62.4 GeV. Weak-decay feed-down errors are included in the error bars on the data points while non-flow and tracking error uncertainties are plotted as bands around v2 = 0, which apply to all identified particles. The errors are asymmetric and the portion of the error band above (below) zero represents the negative (positive) error.

The ratio of Λ v2 to Λbar v2. The data are from minimum bias Au+Au collisions at √sNN = 62.4 and 200 GeV. The bands show the average values of the ratios within the indicated pT ranges.

The ratio of Λ v2 to Λbar v2. The data are from minimum bias Au+Au collisions at √sNN = 62.4 and 200 GeV. The bands show the average values of the ratios within the indicated pT ranges.

The pT integrated ratio of Λ v2 to Λbar v2 for three centrality intervals 0%–10%, 10%–40%, and 40%–80%. The data are from Au+Au collisions at √sNN = 62.4 and 200 GeV.

The pT integrated ratio of Λ v2 to Λbar v2 for three centrality intervals 0%–10%, 10%–40%, and 40%–80%. The data are from Au+Au collisions at √sNN = 62.4 and 200 GeV.

Identified particle v2 from minimum bias collisions at √sNN = 62.4 GeV scaled by the number of valence quarks in the hadron (nq) and plotted versus pT/nq (a) and (mT − m0)/nq (b). In each case a polynomial curve is fit to all particles except pions. The ratio of v2/nq to the fit function is shown in the bottom panels (c) and (d).

Identified particle v2 from minimum bias collisions at √sNN = 62.4 GeV scaled by the number of valence quarks in the hadron (nq) and plotted versus pT/nq (a) and (mT − m0)/nq (b). In each case a polynomial curve is fit to all particles except pions. The ratio of v2/nq to the fit function is shown in the bottom panels (c) and (d).

Identified particle v2 from minimum bias collisions at √sNN = 62.4 GeV scaled by the number of valence quarks in the hadron (nq) and plotted versus pT/nq (a) and (mT − m0)/nq (b). In each case a polynomial curve is fit to all particles except pions. The ratio of v2/nq to the fit function is shown in the bottom panels (c) and (d).

Identified particle v2 from minimum bias collisions at √sNN = 62.4 GeV scaled by the number of valence quarks in the hadron (nq) and plotted versus pT/nq (a) and (mT − m0)/nq (b). In each case a polynomial curve is fit to all particles except pions. The ratio of v2/nq to the fit function is shown in the bottom panels (c) and (d).

Identified particle v2 from minimum bias collisions at √sNN = 62.4 GeV scaled by the number of valence quarks in the hadron (nq) and plotted versus pT/nq (a) and (mT − m0)/nq (b). In each case a polynomial curve is fit to all particles except pions. The ratio of v2/nq to the fit function is shown in the bottom panels (c) and (d).

Identified particle v2 from minimum bias collisions at √sNN = 62.4 GeV scaled by the number of valence quarks in the hadron (nq) and plotted versus pT/nq (a) and (mT − m0)/nq (b). In each case a polynomial curve is fit to all particles except pions. The ratio of v2/nq to the fit function is shown in the bottom panels (c) and (d).

Identified particle v2 from minimum bias collisions at √sNN = 62.4 GeV scaled by the number of valence quarks in the hadron (nq) and plotted versus pT/nq (a) and (mT − m0)/nq (b). In each case a polynomial curve is fit to all particles except pions. The ratio of v2/nq to the fit function is shown in the bottom panels (c) and (d).

Identified particle v2 from minimum bias collisions at √sNN = 62.4 GeV scaled by the number of valence quarks in the hadron (nq) and plotted versus pT/nq (a) and (mT − m0)/nq (b). In each case a polynomial curve is fit to all particles except pions. The ratio of v2/nq to the fit function is shown in the bottom panels (c) and (d).

Identified particle v2 from minimum bias collisions at √sNN = 62.4 GeV scaled by the number of valence quarks in the hadron (nq) and plotted versus pT/nq (a) and (mT − m0)/nq (b). In each case a polynomial curve is fit to all particles except pions. The ratio of v2/nq to the fit function is shown in the bottom panels (c) and (d).

Identified particle v2 from minimum bias collisions at √sNN = 62.4 GeV scaled by the number of valence quarks in the hadron (nq) and plotted versus pT/nq (a) and (mT − m0)/nq (b). In each case a polynomial curve is fit to all particles except pions. The ratio of v2/nq to the fit function is shown in the bottom panels (c) and (d).

Identified particle v2 from minimum bias collisions at √sNN = 62.4 GeV scaled by the number of valence quarks in the hadron (nq) and plotted versus pT/nq (a) and (mT − m0)/nq (b). In each case a polynomial curve is fit to all particles except pions. The ratio of v2/nq to the fit function is shown in the bottom panels (c) and (d).

Identified particle v2 from minimum bias collisions at √sNN = 62.4 GeV scaled by the number of valence quarks in the hadron (nq) and plotted versus pT/nq (a) and (mT − m0)/nq (b). In each case a polynomial curve is fit to all particles except pions. The ratio of v2/nq to the fit function is shown in the bottom panels (c) and (d).

Identified particle v2 from minimum bias collisions at √sNN = 62.4 GeV scaled by the number of valence quarks in the hadron (nq) and plotted versus pT/nq (a) and (mT − m0)/nq (b). In each case a polynomial curve is fit to all particles except pions. The ratio of v2/nq to the fit function is shown in the bottom panels (c) and (d).

Identified particle v2 from minimum bias collisions at √sNN = 62.4 GeV scaled by the number of valence quarks in the hadron (nq) and plotted versus pT/nq (a) and (mT − m0)/nq (b). In each case a polynomial curve is fit to all particles except pions. The ratio of v2/nq to the fit function is shown in the bottom panels (c) and (d).

Identified particle v2 from minimum bias collisions at √sNN = 62.4 GeV scaled by the number of valence quarks in the hadron (nq) and plotted versus pT/nq (a) and (mT − m0)/nq (b). In each case a polynomial curve is fit to all particles except pions. The ratio of v2/nq to the fit function is shown in the bottom panels (c) and (d).

Identified particle v2 from minimum bias collisions at √sNN = 62.4 GeV scaled by the number of valence quarks in the hadron (nq) and plotted versus pT/nq (a) and (mT − m0)/nq (b). In each case a polynomial curve is fit to all particles except pions. The ratio of v2/nq to the fit function is shown in the bottom panels (c) and (d).

Identified particle v2 from minimum bias collisions at √sNN = 62.4 GeV scaled by the number of valence quarks in the hadron (nq) and plotted versus pT/nq (a) and (mT − m0)/nq (b). In each case a polynomial curve is fit to all particles except pions. The ratio of v2/nq to the fit function is shown in the bottom panels (c) and (d).

Identified particle v2 from minimum bias collisions at √sNN = 62.4 GeV scaled by the number of valence quarks in the hadron (nq) and plotted versus pT/nq (a) and (mT − m0)/nq (b). In each case a polynomial curve is fit to all particles except pions. The ratio of v2/nq to the fit function is shown in the bottom panels (c) and (d).

Identified particle v2 from minimum bias collisions at √sNN = 62.4 GeV scaled by the number of valence quarks in the hadron (nq) and plotted versus pT/nq (a) and (mT − m0)/nq (b). In each case a polynomial curve is fit to all particles except pions. The ratio of v2/nq to the fit function is shown in the bottom panels (c) and (d).

Identified particle v2 from minimum bias collisions at √sNN = 62.4 GeV scaled by the number of valence quarks in the hadron (nq) and plotted versus pT/nq (a) and (mT − m0)/nq (b). In each case a polynomial curve is fit to all particles except pions. The ratio of v2/nq to the fit function is shown in the bottom panels (c) and (d).

Identified particle v2 from minimum bias collisions at √sNN = 62.4 GeV scaled by the number of valence quarks in the hadron (nq) and plotted versus pT/nq (a) and (mT − m0)/nq (b). In each case a polynomial curve is fit to all particles except pions. The ratio of v2/nq to the fit function is shown in the bottom panels (c) and (d).

Identified particle v2 from minimum bias collisions at √sNN = 62.4 GeV scaled by the number of valence quarks in the hadron (nq) and plotted versus pT/nq (a) and (mT − m0)/nq (b). In each case a polynomial curve is fit to all particles except pions. The ratio of v2/nq to the fit function is shown in the bottom panels (c) and (d).

Identified particle v2 from minimum bias collisions at √sNN = 62.4 GeV scaled by the number of valence quarks in the hadron (nq) and plotted versus pT/nq (a) and (mT − m0)/nq (b). In each case a polynomial curve is fit to all particles except pions. The ratio of v2/nq to the fit function is shown in the bottom panels (c) and (d).

Identified particle v2 from minimum bias collisions at √sNN = 62.4 GeV scaled by the number of valence quarks in the hadron (nq) and plotted versus pT/nq (a) and (mT − m0)/nq (b). In each case a polynomial curve is fit to all particles except pions. The ratio of v2/nq to the fit function is shown in the bottom panels (c) and (d).

Identified particle v2 from minimum bias collisions at √sNN = 62.4 GeV scaled by the number of valence quarks in the hadron (nq) and plotted versus pT/nq (a) and (mT − m0)/nq (b). In each case a polynomial curve is fit to all particles except pions. The ratio of v2/nq to the fit function is shown in the bottom panels (c) and (d).

Identified particle v2 from minimum bias collisions at √sNN = 62.4 GeV scaled by the number of valence quarks in the hadron (nq) and plotted versus pT/nq (a) and (mT − m0)/nq (b). In each case a polynomial curve is fit to all particles except pions. The ratio of v2/nq to the fit function is shown in the bottom panels (c) and (d).

Identified particle v2 from minimum bias collisions at √sNN = 62.4 GeV scaled by the number of valence quarks in the hadron (nq) and plotted versus pT/nq (a) and (mT − m0)/nq (b). In each case a polynomial curve is fit to all particles except pions. The ratio of v2/nq to the fit function is shown in the bottom panels (c) and (d).

Identified particle v2 from minimum bias collisions at √sNN = 62.4 GeV scaled by the number of valence quarks in the hadron (nq) and plotted versus pT/nq (a) and (mT − m0)/nq (b). In each case a polynomial curve is fit to all particles except pions. The ratio of v2/nq to the fit function is shown in the bottom panels (c) and (d).

Identified particle v2 from minimum bias collisions at √sNN = 62.4 GeV scaled by the number of valence quarks in the hadron (nq) and plotted versus pT/nq (a) and (mT − m0)/nq (b). In each case a polynomial curve is fit to all particles except pions. The ratio of v2/nq to the fit function is shown in the bottom panels (c) and (d).

v2/nq scaled by the mean eccentricity of the initial overlap region versus (mT − m0)/nq for 0%–10%, 10%–40%, and 40%–80% most central Au+Au collisions at √sNN = 62.4 GeV.

v2/nq scaled by the mean eccentricity of the initial overlap region versus (mT − m0)/nq for 0%–10%, 10%–40%, and 40%–80% most central Au+Au collisions at √sNN = 62.4 GeV.

v2/nq scaled by the mean eccentricity of the initial overlap region versus (mT − m0)/nq for 0%–10%, 10%–40%, and 40%–80% most central Au+Au collisions at √sNN = 62.4 GeV.

v2/nq scaled by the mean eccentricity of the initial overlap region versus (mT − m0)/nq for 0%–10%, 10%–40%, and 40%–80% most central Au+Au collisions at √sNN = 62.4 GeV.

v2/nq scaled by the mean eccentricity of the initial overlap region versus (mT − m0)/nq for 0%–10%, 10%–40%, and 40%–80% most central Au+Au collisions at √sNN = 62.4 GeV.

v2/nq scaled by the mean eccentricity of the initial overlap region versus (mT − m0)/nq for 0%–10%, 10%–40%, and 40%–80% most central Au+Au collisions at √sNN = 62.4 GeV.

v2/nq scaled by the mean eccentricity of the initial overlap region versus (mT − m0)/nq for 0%–10%, 10%–40%, and 40%–80% most central Au+Au collisions at √sNN = 62.4 GeV.

v2/nq scaled by the mean eccentricity of the initial overlap region versus (mT − m0)/nq for 0%–10%, 10%–40%, and 40%–80% most central Au+Au collisions at √sNN = 62.4 GeV.

v2/nq scaled by the mean eccentricity of the initial overlap region versus (mT − m0)/nq for 0%–10%, 10%–40%, and 40%–80% most central Au+Au collisions at √sNN = 62.4 GeV.

v2/nq scaled by the mean eccentricity of the initial overlap region versus (mT − m0)/nq for 0%–10%, 10%–40%, and 40%–80% most central Au+Au collisions at √sNN = 62.4 GeV.

v2/nq scaled by the mean eccentricity of the initial overlap region versus (mT − m0)/nq for 0%–10%, 10%–40%, and 40%–80% most central Au+Au collisions at √sNN = 62.4 GeV.

v2/nq scaled by the mean eccentricity of the initial overlap region versus (mT − m0)/nq for 0%–10%, 10%–40%, and 40%–80% most central Au+Au collisions at √sNN = 62.4 GeV.

v2/nq scaled by the mean eccentricity of the initial overlap region versus (mT − m0)/nq for 0%–10%, 10%–40%, and 40%–80% most central Au+Au collisions at √sNN = 62.4 GeV.

v2/nq scaled by the mean eccentricity of the initial overlap region versus (mT − m0)/nq for 0%–10%, 10%–40%, and 40%–80% most central Au+Au collisions at √sNN = 62.4 GeV.

v2/nq scaled by the mean eccentricity of the initial overlap region versus (mT − m0)/nq for 0%–10%, 10%–40%, and 40%–80% most central Au+Au collisions at √sNN = 62.4 GeV.

v2/nq scaled by the mean eccentricity of the initial overlap region versus (mT − m0)/nq for 0%–10%, 10%–40%, and 40%–80% most central Au+Au collisions at √sNN = 62.4 GeV.

v2/nq scaled by the mean eccentricity of the initial overlap region versus (mT − m0)/nq for 0%–10%, 10%–40%, and 40%–80% most central Au+Au collisions at √sNN = 62.4 GeV.

Top panels - minimum bias $v_4$ for pions, charged kaons, $K^0_S$, anti-protons and Λ + Λbar at √sNN = 62.4 GeV. In the left panel the solid (dashed) line shows thevalue for $v_2^2$ for pions (kaons). Intheright panel the dashedline is $v_2^2$ for Λ + Λbar. Bottom panels - $v_4$ scaled by $v_2^2$ (points where $v_4$ and $v_2$ fluctuate around zero are not plotted). Grey bands correspond to the fit results described in the text and Table II. The systematic errors on the $v_4/v_2^2$ ratio from nonflow are included in the error bars leading to asymmetric errors.

Top panels - minimum bias $v_4$ for pions, charged kaons, $K^0_S$, anti-protons and Λ + Λbar at √sNN = 62.4 GeV. In the left panel the solid (dashed) line shows thevalue for $v_2^2$ for pions (kaons). Intheright panel the dashedline is $v_2^2$ for Λ + Λbar. Bottom panels - $v_4$ scaled by $v_2^2$ (points where $v_4$ and $v_2$ fluctuate around zero are not plotted). Grey bands correspond to the fit results described in the text and Table II. The systematic errors on the $v_4/v_2^2$ ratio from nonflow are included in the error bars leading to asymmetric errors.

Top panels - minimum bias $v_4$ for pions, charged kaons, $K^0_S$, anti-protons and Λ + Λbar at √sNN = 62.4 GeV. In the left panel the solid (dashed) line shows thevalue for $v_2^2$ for pions (kaons). Intheright panel the dashedline is $v_2^2$ for Λ + Λbar. Bottom panels - $v_4$ scaled by $v_2^2$ (points where $v_4$ and $v_2$ fluctuate around zero are not plotted). Grey bands correspond to the fit results described in the text and Table II. The systematic errors on the $v_4/v_2^2$ ratio from nonflow are included in the error bars leading to asymmetric errors.

Top panels - minimum bias $v_4$ for pions, charged kaons, $K^0_S$, anti-protons and Λ + Λbar at √sNN = 62.4 GeV. In the left panel the solid (dashed) line shows thevalue for $v_2^2$ for pions (kaons). Intheright panel the dashedline is $v_2^2$ for Λ + Λbar. Bottom panels - $v_4$ scaled by $v_2^2$ (points where $v_4$ and $v_2$ fluctuate around zero are not plotted). Grey bands correspond to the fit results described in the text and Table II. The systematic errors on the $v_4/v_2^2$ ratio from nonflow are included in the error bars leading to asymmetric errors.

Top panels - minimum bias $v_4$ for pions, charged kaons, $K^0_S$, anti-protons and Λ + Λbar at √sNN = 62.4 GeV. In the left panel the solid (dashed) line shows thevalue for $v_2^2$ for pions (kaons). Intheright panel the dashedline is $v_2^2$ for Λ + Λbar. Bottom panels - $v_4$ scaled by $v_2^2$ (points where $v_4$ and $v_2$ fluctuate around zero are not plotted). Grey bands correspond to the fit results described in the text and Table II. The systematic errors on the $v_4/v_2^2$ ratio from nonflow are included in the error bars leading to asymmetric errors.

Top panels - minimum bias $v_4$ for pions, charged kaons, $K^0_S$, anti-protons and Λ + Λbar at √sNN = 62.4 GeV. In the left panel the solid (dashed) line shows thevalue for $v_2^2$ for pions (kaons). Intheright panel the dashedline is $v_2^2$ for Λ + Λbar. Bottom panels - $v_4$ scaled by $v_2^2$ (points where $v_4$ and $v_2$ fluctuate around zero are not plotted). Grey bands correspond to the fit results described in the text and Table II. The systematic errors on the $v_4/v_2^2$ ratio from nonflow are included in the error bars leading to asymmetric errors.

Top panels - minimum bias $v_4$ for pions, charged kaons, $K^0_S$, anti-protons and Λ + Λbar at √sNN = 62.4 GeV. In the left panel the solid (dashed) line shows thevalue for $v_2^2$ for pions (kaons). Intheright panel the dashedline is $v_2^2$ for Λ + Λbar. Bottom panels - $v_4$ scaled by $v_2^2$ (points where $v_4$ and $v_2$ fluctuate around zero are not plotted). Grey bands correspond to the fit results described in the text and Table II. The systematic errors on the $v_4/v_2^2$ ratio from nonflow are included in the error bars leading to asymmetric errors.

Top panels - minimum bias $v_4$ for pions, charged kaons, $K^0_S$, anti-protons and Λ + Λbar at √sNN = 62.4 GeV. In the left panel the solid (dashed) line shows thevalue for $v_2^2$ for pions (kaons). Intheright panel the dashedline is $v_2^2$ for Λ + Λbar. Bottom panels - $v_4$ scaled by $v_2^2$ (points where $v_4$ and $v_2$ fluctuate around zero are not plotted). Grey bands correspond to the fit results described in the text and Table II. The systematic errors on the $v_4/v_2^2$ ratio from nonflow are included in the error bars leading to asymmetric errors.

Top panel - $v_2$ for pions and protons at √sNN = 62.4 and 17.3 GeV. The 62.4 GeV data are from TOF and dE/dx measurements combined. Middle and bottom panel - ratios of $v_2$ for $\pi^++\pi^-$, $K^0_S$, p+p, Λ+Λbar and at different center-of-mass energies scaled by the values at 62.4 GeV. The grey and yellow bands represent systematic uncertainties in the v2 ratios arising from non-flow effects. The grey bands (above unity) are the uncertainties for the 200 GeV/62.4 GeV data and the yellow bands (below unity) are for the 17.3 GeV/62.4 GeV data.

Top panel - $v_2$ for pions and protons at √sNN = 62.4 and 17.3 GeV. The 62.4 GeV data are from TOF and dE/dx measurements combined. Middle and bottom panel - ratios of $v_2$ for $\pi^++\pi^-$, $K^0_S$, p+p, Λ+Λbar and at different center-of-mass energies scaled by the values at 62.4 GeV. The grey and yellow bands represent systematic uncertainties in the v2 ratios arising from non-flow effects. The grey bands (above unity) are the uncertainties for the 200 GeV/62.4 GeV data and the yellow bands (below unity) are for the 17.3 GeV/62.4 GeV data.

Top panel - $v_2$ for pions and protons at √sNN = 62.4 and 17.3 GeV. The 62.4 GeV data are from TOF and dE/dx measurements combined. Middle and bottom panel - ratios of $v_2$ for $\pi^++\pi^-$, $K^0_S$, p+p, Λ+Λbar and at different center-of-mass energies scaled by the values at 62.4 GeV. The grey and yellow bands represent systematic uncertainties in the v2 ratios arising from non-flow effects. The grey bands (above unity) are the uncertainties for the 200 GeV/62.4 GeV data and the yellow bands (below unity) are for the 17.3 GeV/62.4 GeV data.

Top panel - $v_2$ for pions and protons at √sNN = 62.4 and 17.3 GeV. The 62.4 GeV data are from TOF and dE/dx measurements combined. Middle and bottom panel - ratios of $v_2$ for $\pi^++\pi^-$, $K^0_S$, p+p, Λ+Λbar and at different center-of-mass energies scaled by the values at 62.4 GeV. The grey and yellow bands represent systematic uncertainties in the v2 ratios arising from non-flow effects. The grey bands (above unity) are the uncertainties for the 200 GeV/62.4 GeV data and the yellow bands (below unity) are for the 17.3 GeV/62.4 GeV data.

Top panel - $v_2$ for pions and protons at √sNN = 62.4 and 17.3 GeV. The 62.4 GeV data are from TOF and dE/dx measurements combined. Middle and bottom panel - ratios of $v_2$ for $\pi^++\pi^-$, $K^0_S$, p+p, Λ+Λbar and at different center-of-mass energies scaled by the values at 62.4 GeV. The grey and yellow bands represent systematic uncertainties in the v2 ratios arising from non-flow effects. The grey bands (above unity) are the uncertainties for the 200 GeV/62.4 GeV data and the yellow bands (below unity) are for the 17.3 GeV/62.4 GeV data.

Top panel - $v_2$ for pions and protons at √sNN = 62.4 and 17.3 GeV. The 62.4 GeV data are from TOF and dE/dx measurements combined. Middle and bottom panel - ratios of $v_2$ for $\pi^++\pi^-$, $K^0_S$, p+p, Λ+Λbar and at different center-of-mass energies scaled by the values at 62.4 GeV. The grey and yellow bands represent systematic uncertainties in the v2 ratios arising from non-flow effects. The grey bands (above unity) are the uncertainties for the 200 GeV/62.4 GeV data and the yellow bands (below unity) are for the 17.3 GeV/62.4 GeV data.

Top panel - $v_2$ for pions and protons at √sNN = 62.4 and 17.3 GeV. The 62.4 GeV data are from TOF and dE/dx measurements combined. Middle and bottom panel - ratios of $v_2$ for $\pi^++\pi^-$, $K^0_S$, p+p, Λ+Λbar and at different center-of-mass energies scaled by the values at 62.4 GeV. The grey and yellow bands represent systematic uncertainties in the v2 ratios arising from non-flow effects. The grey bands (above unity) are the uncertainties for the 200 GeV/62.4 GeV data and the yellow bands (below unity) are for the 17.3 GeV/62.4 GeV data.

Top panel - $v_2$ for pions and protons at √sNN = 62.4 and 17.3 GeV. The 62.4 GeV data are from TOF and dE/dx measurements combined. Middle and bottom panel - ratios of $v_2$ for $\pi^++\pi^-$, $K^0_S$, p+p, Λ+Λbar and at different center-of-mass energies scaled by the values at 62.4 GeV. The grey and yellow bands represent systematic uncertainties in the v2 ratios arising from non-flow effects. The grey bands (above unity) are the uncertainties for the 200 GeV/62.4 GeV data and the yellow bands (below unity) are for the 17.3 GeV/62.4 GeV data.

(a) The corrected integrated yield $dN/dy$ at mid-rapidity for $\overline{\Xi}^{+}$ , $\overline{\Lambda}$ and $\overline{p}$ divided by $N_{part}$, normalized to the most peripheral centrality interval (60 − 80%), plotted as a function of $N_{part}$. The gray, black and dashed bands represent the errors on the normalization to the most peripheral bin for the $\overline{\Xi}^{+}$, $\overline{\Lambda}$ and $\overline{p}$. Other errors shown are statistical only. (b) $\gamma_{s}$ as a function of $N_{part}$ calculated from thermal model fits to the measured particle yields ($\pi$,K,p [24], $\Lambda$, $\Xi$, $\Omega$ and their anti-particles) at 200 GeV. Values for $e^{+}+e^{−}$ and $p+\overline{p}$ collisions at $\sqrt{s_{NN}}$ = 91 and 200 GeV respectively and for Pb+Pb SPS collisions at $\sqrt{s_{NN}}$ = 17.2 GeV are shown for comparison [26, 27, 28].

(a) The corrected integrated yield $dN/dy$ at mid-rapidity for $\overline{\Xi}^{+}$ , $\overline{\Lambda}$ and $\overline{p}$ divided by $N_{part}$, normalized to the most peripheral centrality interval (60 − 80%), plotted as a function of $N_{part}$. The gray, black and dashed bands represent the errors on the normalization to the most peripheral bin for the $\overline{\Xi}^{+}$, $\overline{\Lambda}$ and $\overline{p}$. Other errors shown are statistical only. (b) $\gamma_{s}$ as a function of $N_{part}$ calculated from thermal model fits to the measured particle yields ($\pi$,K,p [24], $\Lambda$, $\Xi$, $\Omega$ and their anti-particles) at 200 GeV. Values for $e^{+}+e^{−}$ and $p+\overline{p}$ collisions at $\sqrt{s_{NN}}$ = 91 and 200 GeV respectively and for Pb+Pb SPS collisions at $\sqrt{s_{NN}}$ = 17.2 GeV are shown for comparison [26, 27, 28].

$R_{CP}$ for $\Xi^{-}+\overline{\Xi}^{+}$ and $\Omega^{-}+\overline{\Omega}^{+}$ at mid-rapidity (centrality interval : 0 − 5% vs. 40 − 60%). A dashed line for charged hadrons and gray band for $\Lambda^{-}+\overline{\Lambda}$ are shown as comparison. The gray rectangles represent participant and binary scalings.

$\sqrt{s_{NN}}$ dependence of $<p_t>$ fluctuations for central collisions and STAR full acceptance. CERES fluctuation data (12.3 and 17.3 GeV) were linearly extrapolated to the STAR $\eta$ acceptance

Hijing-1.37 $<p_t>$ fluctuations for central collisions and the STAR acceptance at three energies, for quench-on, quench-off and jets-off collisions.

(Color Online) Variation of $N_{ch}$ normalized to the number of collisions in the FTPC coverage $(2.9 \leq \eta \leq 3.9)$ and $N_{\gamma}$ normalized to number of collisions, in the PMD coverage $(2.3 \leq \eta \leq 3.7)$ as a function of $N_{coll}$. The lower band shows the uncertainty in the ratio due to uncertainties in $N_{coll}$ calculations.

(Color Online) $dN/d\eta$ for charged particles and photons for Au + Au collisions at $\sqrt{s_{NN}}$ = 62.4 GeV for various event centrality classes.

(Color Online) $dN/d\eta$ for charged particles and photons for Au + Au collisions at $\sqrt{s_{NN}}$ = 62.4 GeV for various event centrality classes.

(Color Online) Half width at half maximum of the pseudorapidity distributions ($\eta_{h}$) of charged particles as a function of total charged particle multiplicity ($N_{T}$) normalized to the center of mass energy. The Au + Au collision data are from the PHOBOS [8] experiment and p + p collision data are from the ISR [31] experiments.

(Color Online) Half width at half maximum of the pseudorapidity distributions ($\eta_{h}$) of charged particles as a function of total charged particle multiplicity ($N_{T}$) normalized to the center of mass energy. The Au + Au collision data are from the PHOBOS [8] experiment and p + p collision data are from the ISR [31] experiments.

(Color Online) Variation of $dN_{ch}/d\eta$ normalized to $N_{part}$ with $\eta – y_{beam}$ for central and peripheral collisions for positively charged hadrons ($h^{+}$) and negatively charged hadrons ($h^{−}$).

(Color Online) The top panel shows the variation of pion rapidity density normalized to $N_{part}$ with $y – y_{beam}$ for central collisions at various collision energies. Also shown is the estimated $dN_{\pi^{0}}/dy$ obtained from $dN\_{\gamma}/dy$ normalized to $N_{part}$. The bottom panel shows the variation of net proton rapidity density normalized to $N_{part}$ with $y – y_{beam}$ for central collisions at various collision energies.

(Color Online) The top panel shows the variation of pion rapidity density normalized to $N_{part}$ with $y – y_{beam}$ for central collisions at various collision energies. Also shown is the estimated $dN_{\pi^{0}}/dy$ obtained from $dN\_{\gamma}/dy$ normalized to $N_{part}$. The bottom panel shows the variation of net proton rapidity density normalized to $N_{part}$ with $y – y_{beam}$ for central collisions at various collision energies.

Directed flow of charged particles as a function of pseudorapidity, for centrality 10%-70%.

Directed flow of charged particles as a function of pseudorapidity, for centrality 10%-70%.

Directed flow of charged particles as a function of pseudorapidity, for centrality 10%-70%.

$v_1$ versus rapidity for charged particles.

$v_1$ versus rapidity for pions.

$v_1$ versus rapidity for protons.

Directed flow of charged particles as a function of pseudorapidity for the 60-70% and 70-80% centrality classes.

Directed flow of charged particles as a function of pseudorapidity for the 10-20%, 20-30%, 30-40%, 40-50% and 50-60% centrality classes.

Directed flow of charged particles as a function of pseudorapidity for the 5-10% centrality class.

Directed flow of charged particles as a function of pseudorapidity for the 0-5% centrality class.

$v_1${3} versus $p_t$ measured in the main TPC (|$\eta$| < 1.3), for centrality 10%-70%.

$v_1${3} versus $p_t$ measured in the Forward TPC (2.5 < |$\eta$| < 4.0), for different centralities.

Directed flow of charged particles as a function of impact parameter for the mid-pseudorapidity region (|$\eta$| < 1.3) and the forward pseudorapidity region (2.5 < |$\eta$| < 4.0.

Charged particle $v_1$ for Au+Au collisions (10%-70%) at 200 GeV as a function of $\eta-y_{beam}$.

Charged particle $v_1$ for Au+Au collisions (10%-70%) at 62.4 GeV as a function of $\eta-y_{beam}$.

Charged particle $v_1$ for Pb+Pb collisions (12.5%-33.5%) at 17.2 GeV as a function of $y-y_{beam}$.

Average transverse momentum per event for Au+Au at $\sqrt{s_{NN}}$ = 200 GeV for the 5% most central collisions.

$<\Delta p_{t,i}\Delta p_{t,j}>$ as a function of centrality and incident energy for Au+Au collisions compared with HIJING results.

(d$N/\textrm{d}\eta)<\Delta p_{t,i}\Delta p_{t,j}>$ as a function of centrality and incident energy for Au+Au collisions compared with HIJING results.

$(<\Delta p_{t,i}\Delta p_{t,j}>)^{1/2}/<<p_{t}>>$ as a function of centrality and incident energy for Au+Au collisions compared with HIJING results.

$(<\Delta p_{t,i}\Delta p_{t,j}>)^{1/2}/<<p_{t}>>$ as a function of incident energy for 0-5% most central Au+Au collisions compared with CERES results.

Statistal errors only.

Calculated primary D/S+ cross section. BR2 = BR(D/S+ --> PHI PI+) * BR(PHI --> K+ K-).