Parameter $\alpha$ from the fit to the data.

Results of the measurements by PHENIX at $\sqrt{s_{NN}}$ = 200 GeV.

Results of the measurements by PHENIX at $\sqrt{s_{NN}}$ = 130 GeV.

Results of the measurements by PHENIX at $\sqrt{s_{NN}}$ = 19.6 GeV.

Ratios of measured quantities at 200 GeV/130 GeV and 200 GeV/19.6 GeV. The number of $N_P$ is the average between two energies.

Ratios of measured quantities at 200 GeV/130 GeV and 200 GeV/19.6 GeV. The number of $N_P$ is the average between two energies.

200, 130, and 19.6 GeV are RHIC average values, and 17.2, 8.7, and 4.8 GeV are SPS average values of $dN_{ch}$/$d\eta$/($0.5N_p$) at different $\sqrt{s_{NN}}$. An additional 5% error (recalc.) is added for collision energies 4.8 - 17.6 GeV for the uncertainty related to recalculation to the Center of Mass system.

200, 130, and 19.6 GeV are RHIC average values, and 17.2, 8.7, and 4.8 GeV are SPS average values of $dN_{ch}$/$d\eta$/($0.5N_p$) at different $\sqrt{s_{NN}}$. An additional 5% error (recalc.) is added for collision energies 4.8 - 17.6 GeV for the uncertainty related to recalculation to the Center of Mass system.

200, 130, and 19.6 GeV are RHIC average values, and 17.2, 8.7, and 4.8 GeV are SPS average values of $dN_{ch}$/$d\eta$/($0.5N_p$) at different $\sqrt{s_{NN}}$. An additional 5% error (recalc.) is added for collision energies 4.8 - 17.6 GeV for the uncertainty related to recalculation to the Center of Mass system.

200, 130, and 19.6 GeV are RHIC average values, and 17.2, 8.7, and 4.8 GeV are SPS average values of $dN_{ch}$/$d\eta$/($0.5N_p$) at different $\sqrt{s_{NN}}$. An additional 5% error (recalc.) is added for collision energies 4.8 - 17.6 GeV for the uncertainty related to recalculation to the Center of Mass system.

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.

Fit of elastic ($m_Y=m_p$) and proton-dissociative ($1<m_Y<10$ GeV) $\pi^{+}\pi^{-}$ photoproduction off protons with a Soeding-inspired analytic function --- statistical correlations only. Only one half of the (symmetric) matrix is stored.

Elastic ($m_Y=m_p$) and proton-dissociative ($1<m_Y<10$ GeV) $\pi^{+}\pi^{-}$ photoproduction off protons, differential in the dipion mass and in bins of $W$. The tabulated cross sections are $\gamma p$ cross sections but can be converted to $ep$ cross sections using the effective photon flux $\Phi_{\gamma/e}$.

Elastic ($m_Y=m_p$) and proton-dissociative ($1<m_Y<10$ GeV) $\pi^{+}\pi^{-}$ photoproduction off protons, differential in the dipion mass and in bins of $W$ --- statistical correlations coefficients $\rho_{ij}$ only. Only one half of the (symmetric) matrix is stored. Bins are identified by their global bin number.

Elastic ($m_Y=m_p$) and proton-dissociative ($1<m_Y<10$ GeV) $\rho^0(770)$ photoproduction off protons, in bins of the photon-proton energy $W$. The cross section is defined as the integral of the relativistic Breit Wigner resonance in the dipion mass over the range $2m_\pi<m_{\pi\pi}<1.53$ GeV. The tabulated cross sections are $\gamma p$ cross sections but can be converted to $ep$ cross sections using the effective photon flux $\Phi_{\gamma/e}$.

Elastic ($m_Y=m_p$) and proton-dissociative ($1<m_Y<10$ GeV) $\rho^0(770)$ photoproduction off protons in bins of the photon-proton energy $W$ --- statistical correlations coefficients $\rho_{ij}$ only. Only one half of the (symmetric) matrix is stored. Bins are identified by their global bin number.

Fit of elastic ($m_Y=m_p$) and proton-dissociative ($1<m_Y<10$ GeV) $\rho^0(770)$ photoproduction cross sections off protons as a function of energy. Parameters with subscript "el" and "pd" correspond to elastic and proton-dissociative cross sections, respectively.

Fit of elastic ($m_Y=m_p$) and proton-dissociative ($1<m_Y<10$ GeV) $\rho^0(770)$ photoproduction cross sections off protons as a function of energy --- statistical correlations coefficients $\rho_{ij}$ only. Only one half of the (symmetric) matrix is stored. Bins are identified by their global bin number.

Fit of elastic $\rho^0(770)$ photoproduction cross section off protons as a function of energy (various experiments)

Fit of elastic $\rho^0(770)$ photoproduction cross sections off protons as a function of energy (various experiments) --- statistical correlations coefficients $\rho_{ij}$ only. Only one half of the (symmetric) matrix is stored. Bins are identified by their global bin number.

Elastic ($m_Y=m_p$) and proton-dissociative ($1<m_Y<10$ GeV) $\pi^{+}\pi^{-}$ photoproduction cross section off protons, double-differential in the dipion mass and the momentum transfer $\vert t\vert$. The tabulated cross sections are $\gamma p$ cross sections but can be converted to $ep$ cross sections using the effective photon flux $\Phi_{\gamma/e}$.

Elastic ($m_Y=m_p$) and proton-dissociative ($1<m_Y<10$ GeV) $\pi^{+}\pi^{-}$ photoproduction cross section off protons, double-differential in the dipion mass and the momentum transfer $\vert t\vert$ --- statistical correlations coefficients $\rho_{ij}$ only. Only one half of the (symmetric) matrix is stored. Bins are identified by their global bin number.

Elastic ($m_Y=m_p$) and proton-dissociative ($1<m_Y<10$ GeV) $\rho^0(770)$ photoproduction off protons, single-differential in the momentum transfer $\vert t\vert$. The cross section is defined as the integral of the relativistic Breit Wigner resonance in the dipion mass over the range $2m_\pi<m_{\pi\pi}<1.53$ GeV. The tabulated cross sections are $\gamma p$ cross sections but can be converted to $ep$ cross sections using the effective photon flux $\Phi_{\gamma/e}$.

Elastic ($m_Y=m_p$) and proton-dissociative ($1<m_Y<10$ GeV) $\rho^0(770)$ photoproduction off protons, single differential in the momentum transfer $\vert t\vert$ --- statistical correlations coefficients $\rho_{ij}$ only. Only one half of the (symmetric) matrix is stored. Bins are identified by their global bin number.

Fit of elastic ($m_Y=m_p$) and proton-dissociative ($1<m_Y<10$ GeV) $\rho^0(770)$ photoproduction single-differential cross sections off protons as a function of momentum transfer $t$. Parameters with subscript "el" and "pd" correspond to elastic and proton-dissociative cross sections, respectively.

Fit of elastic ($m_Y=m_p$) and proton-dissociative ($1<m_Y<10$ GeV) $\rho^0(770)$ photoproduction single-differential cross sections off protons as a function of momentum transfer $t$ --- statistical correlations coefficients $\rho_{ij}$ only. Only one half of the (symmetric) matrix is stored. Bins are identified by their global bin number.

Elastic ($m_Y=m_p$) and proton-dissociative ($1<m_Y<10$ GeV) $\pi^{+}\pi^{-}$ photoproduction double-differenial cross section off protons, as a function of the dipion mass and the momentum transfer $\vert t\vert$, in bins of the photon-proton energy $W$ The tabulated cross sections are $\gamma p$ cross sections but can be converted to $ep$ cross sections using the effective photon flux $\Phi_{\gamma/e}$.

Elastic ($m_Y=m_p$) and proton-dissociative ($1<m_Y<10$ GeV) $\pi^{+}\pi^{-}$ photoproduction double-differenial cross section off protons, as a function of the dipion mass and the momentum transfer $\vert t\vert$, in bins of the photon-proton energy $W$ --- statistical correlations coefficients $\rho_{ij}$ only. Only one half of the (symmetric) matrix is stored. Bins are identified by their global bin number.

Elastic ($m_Y=m_p$) and proton-dissociative ($1<m_Y<10$ GeV) $\rho^0(770)$ photoproduction off protons, single-differential in the momentum transfer $\vert t\vert$, in bins of the photon-proton energy $W$. The cross section is defined as the integral of the relativistic Breit Wigner resonance in the dipion mass over the range $2m_\pi<m_{\pi\pi}<1.53$ GeV. The tabulated cross sections are $\gamma p$ cross sections but can be converted to $ep$ cross sections using the effective photon flux $\Phi_{\gamma/e}$.

Elastic ($m_Y=m_p$) and proton-dissociative ($1<m_Y<10$ GeV) $\rho^0(770)$ photoproduction cross section off protons, single differential in the momentum transfer $\vert t\vert$ and in bins of the photon-proton energy $W$ --- statistical correlations coefficients $\rho_{ij}$ only. Only one half of the (symmetric) matrix is stored. Bins are identified by their global bin number.

Regge fit of elastic ($m_Y=m_p$) and proton-dissociative ($1<m_Y<10$ GeV) $\rho^0(770)$ photoproduction single-differential cross sections off protons as a function of momentum transfer $t$ and photon-proton energy $W$. Parameters with subscript "el" and "pd" correspond to elastic and proton-dissociative cross sections, respectively.

Regge fit of elastic ($m_Y=m_p$) and proton-dissociative ($1<m_Y<10$ GeV) $\rho^0(770)$ photoproduction single-differential cross sections off protons as a function of momentum transfer $t$ and photon-proton energy $W$ --- statistical correlations coefficients $\rho_{ij}$ only. Only one half of the (symmetric) matrix is stored. Bins are identified by their global bin number.

Fit of $b$-slopes in elastic ($m_Y=m_p$) $\rho^0(770)$ photoproduction single-differential cross sections off protons as a function of momentum transfer $t$ in bins of photon-proton energy $W$.

Fit of $b$-slopes in elastic ($m_Y=m_p$) $\rho^0(770)$ photoproduction single-differential cross sections off protons as a function of momentum transfer $t$ in bins of photon-proton energy $W$ --- statistical correlations coefficients $\rho_{ij}$ only. Only one half of the (symmetric) matrix is stored. Bins are identified by their global bin number.

Fit of Pomeron trajectories in elastic ($m_Y=m_p$) $\rho^0(770)$ photoproduction single-differential cross sections off protons as a function of momentum transfer $t$ in bins of photon-proton energy $W$.

Fit of Pomeron trajectories in elastic ($m_Y=m_p$) $\rho^0(770)$ photoproduction single-differential cross sections off protons as a function of momentum transfer $t$ in bins of photon-proton energy $W$ --- statistical correlations coefficients $\rho_{ij}$ only. Only one half of the (symmetric) matrix is stored. Bins are identified by their global bin number.

The TPC mid-rapidity multiplicity distributions ($|\eta|$ < 0.5) for the corresponding E-FTPC selected centrality bins.

Uncorrected charged particle multiplicity distribution measured in the TPC within $|\eta|$ < 0.5 in pp collisions at 200 GeV.

The ratio of the transverse overlap area $S_{⊥}$ to $(N_{part}/2)^{2/3}$ versus $(N_{part}/2)^{2/3}$.

The ratio of the transverse overlap area $S_{⊥}$ to $(N_{part}/2)^{2/3}$ versus $(N_{part}/2)^{2/3}$.

The ratio of the transverse overlap area $S_{⊥}$ to $(N_{part}/2)^{2/3}$ versus $(N_{part}/2)^{2/3}$.

The ratio of the charged pion multiplicity to the transverse overlap area $dN_{\pi}/dy/S_{⊥}$ versus $N_{coll}/N_{part}$. Errors shown are total errors, dominated by systematic uncertainties. The systematic uncertainties are correlated between Npart, Ncoll, and S⊥, and are largely canceled in the plotted ratio quantities.

The ratio of the charged pion multiplicity to the transverse overlap area $dN_{\pi}/dy/S_{⊥}$ versus $N_{coll}/N_{part}$. Errors shown are total errors, dominated by systematic uncertainties. The systematic uncertainties are correlated between Npart, Ncoll, and S⊥, and are largely canceled in the plotted ratio quantities.

The ratio of the charged pion multiplicity to the transverse overlap area $dN_{\pi}/dy/S_{⊥}$ versus $N_{coll}/N_{part}$. Errors shown are total errors, dominated by systematic uncertainties. The systematic uncertainties are correlated between Npart, Ncoll, and S⊥, and are largely canceled in the plotted ratio quantities.

Energy loss effect for $\pi^{+-}$ (a), $K^{+-}$ (b), and p and pbar (c) at mid-rapidity (|y| < 0.1) as a function of particle momentum magnitude in 200 GeV pp and 62.4 GeV central 0-5% Au+Au collisions. Only negative particles are shown; energy loss for particles and antiparticles are the same. Errors shown are statistical only. The pion energy loss is already corrected by the track reconstruction algorithm.

Energy loss effect for $\pi^{+-}$ (a), $K^{+-}$ (b), and p and pbar (c) at mid-rapidity (|y| < 0.1) as a function of particle momentum magnitude in 200 GeV pp and 62.4 GeV central 0-5% Au+Au collisions. Only negative particles are shown; energy loss for particles and antiparticles are the same. Errors shown are statistical only. The pion energy loss is already corrected by the track reconstruction algorithm.

Energy loss effect for $\pi^{+-}$ (a), $K^{+-}$ (b), and p and pbar (c) at mid-rapidity (|y| < 0.1) as a function of particle momentum magnitude in 200 GeV pp and 62.4 GeV central 0-5% Au+Au collisions. Only negative particles are shown; energy loss for particles and antiparticles are the same. Errors shown are statistical only. The pion energy loss is already corrected by the track reconstruction algorithm.

Vertex-finding efficiency (ǫvtx) and fake vertex rate (δfake) as a function of the number of good global tracks (a) and the number of good primary tracks (b). Errors shown or smaller than the point size are statistical only.

Vertex-finding efficiency (ǫvtx) and fake vertex rate (δfake) as a function of the number of good global tracks (a) and the number of good primary tracks (b). Errors shown or smaller than the point size are statistical only.

The $p_{T}$ dependent correction to particle spectra due to fake vertex events, $\epsilon_{fake}(p_{T})$, in 200 GeV minimum bias pp and d+Au collisions. Errors shown are statistical only.

The $p_{T}$ dependent correction to particle spectra due to fake vertex events, $\epsilon_{fake}(p_{T})$, in 200 GeV minimum bias pp and d+Au collisions. Errors shown are statistical only.

Efficiency (product of tracking efficiency and detector acceptance) of π−, K−, and p in pp (a) and d+Au collisions (b) at 200 GeV as a function of input MC p⊥. Errors shown are binomial errors. The curves are parameterizations to the efficiency data and are used for corrections in the analysis.

Efficiency (product of tracking efficiency and detector acceptance) of π−, K−, and p in pp (a) and d+Au collisions (b) at 200 GeV as a function of input MC p⊥. Errors shown are binomial errors. The curves are parameterizations to the efficiency data and are used for corrections in the analysis.

Efficiency (product of tracking efficiency and detector acceptance) of π−, K−, and p in 70-80% peripheral Au+Au (a) and 0-5% central Au+Au collisions (b) at 62.4 GeV as a function of input MC p⊥. Errors shown are binomial errors. The curves are parameterizations to the efficiency data and are used for corrections in the analysis.

Efficiency (product of tracking efficiency and detector acceptance) of π−, K−, and p in 70-80% peripheral Au+Au (a) and 0-5% central Au+Au collisions (b) at 62.4 GeV as a function of input MC p⊥. Errors shown are binomial errors. The curves are parameterizations to the efficiency data and are used for corrections in the analysis.

Efficiency (product of tracking efficiency and detector acceptance) of π−, K−, and p in 70-80% peripheral Au+Au (a) and 0-5% central Au+Au collisions (b) at 62.4 GeV as a function of input MC p⊥. Errors shown are binomial errors. The curves are parameterizations to the efficiency data and are used for corrections in the analysis.

Efficiency (product of tracking efficiency and detector acceptance) of π−, K−, and p in 70-80% peripheral Au+Au (a) and 0-5% central Au+Au collisions (b) at 62.4 GeV as a function of input MC p⊥. Errors shown are binomial errors. The curves are parameterizations to the efficiency data and are used for corrections in the analysis.

The dca distributions of protons and antiprotons for 0.40 < p⊥ < 0.45 GeV/c in 200 GeV minimum bias d+Au. Errors shown are statistical only.

The dca distributions of protons and antiprotons for 0.70 < p⊥ < 0.75 GeV/c in 200 GeV minimum bias d+Au. Errors shown are statistical only.

The dca distributions of protons and antiprotons for 0.40 < p⊥ < 0.45 GeV/c in 62.4 GeV 0-5% central Au+Au collisions. Errors shown are statistical only.

The dca distributions of protons and antiprotons for 0.70 < p⊥ < 0.75 GeV/c in 62.4 GeV 0-5% central Au+Au collisions. Errors shown are statistical only.

Pion background fraction from weak decays ($\Lambda, K_S^0$) and $\mu^{+-}$ contamination as a function of p⊥ in minimum bias d+Au collisions at 200 GeV. Errors shown are statistical only.

Pion background fraction from weak decays ($\Lambda, K_S^0$) and $\mu^{+-}$ contamination as a function of p⊥ in minimum bias d+Au collisions at 200 GeV. Errors shown are statistical only.

Mid-rapidity identified antiproton spectra in 200 GeV minimum bias d+Au (a) and 62.4 GeV central Au+Au collisions (b) measured by dE/dx together with those by TOF [46, 47]. The dE/dx data are from |y| < 0.1 and the TOF data are from |y| < 0.5. The curves are various fits to the dE/dx data for extrapolation. The quadratic sum of statistical errors and point-to-point systematic errors are plotted, but are smaller than the point size.

Mid-rapidity identified antiproton spectra in 200 GeV minimum bias d+Au (a) and 62.4 GeV central Au+Au collisions (b) measured by dE/dx together with those by TOF [46, 47]. The dE/dx data are from |y| < 0.1 and the TOF data are from |y| < 0.5. The curves are various fits to the dE/dx data for extrapolation. The quadratic sum of statistical errors and point-to-point systematic errors are plotted, but are smaller than the point size.

Mid-rapidity identified antiproton spectra in 200 GeV minimum bias d+Au (a) and 62.4 GeV central Au+Au collisions (b) measured by dE/dx together with those by TOF [46, 47]. The dE/dx data are from |y| < 0.1 and the TOF data are from |y| < 0.5. The curves are various fits to the dE/dx data for extrapolation. The quadratic sum of statistical errors and point-to-point systematic errors are plotted, but are smaller than the point size.

Mid-rapidity identified antiproton spectra in 200 GeV minimum bias d+Au (a) and 62.4 GeV central Au+Au collisions (b) measured by dE/dx together with those by TOF [46, 47]. The dE/dx data are from |y| < 0.1 and the TOF data are from |y| < 0.5. The curves are various fits to the dE/dx data for extrapolation. The quadratic sum of statistical errors and point-to-point systematic errors are plotted, but are smaller than the point size.

Mid-rapidity identified antiproton spectra in 200 GeV minimum bias d+Au (a) and 62.4 GeV central Au+Au collisions (b) measured by dE/dx together with those by TOF [46, 47]. The dE/dx data are from |y| < 0.1 and the TOF data are from |y| < 0.5. The curves are various fits to the dE/dx data for extrapolation. The quadratic sum of statistical errors and point-to-point systematic errors are plotted, but are smaller than the point size.

Mid-rapidity identified antiproton spectra in 200 GeV minimum bias d+Au (a) and 62.4 GeV central Au+Au collisions (b) measured by dE/dx together with those by TOF [46, 47]. The dE/dx data are from |y| < 0.1 and the TOF data are from |y| < 0.5. The curves are various fits to the dE/dx data for extrapolation. The quadratic sum of statistical errors and point-to-point systematic errors are plotted, but are smaller than the point size.

Mid-rapidity identified antiproton spectra in 200 GeV minimum bias d+Au (a) and 62.4 GeV central Au+Au collisions (b) measured by dE/dx together with those by TOF [46, 47]. The dE/dx data are from |y| < 0.1 and the TOF data are from |y| < 0.5. The curves are various fits to the dE/dx data for extrapolation. The quadratic sum of statistical errors and point-to-point systematic errors are plotted, but are smaller than the point size.

Mid-rapidity identified antiproton spectra in 200 GeV minimum bias d+Au (a) and 62.4 GeV central Au+Au collisions (b) measured by dE/dx together with those by TOF [46, 47]. The dE/dx data are from |y| < 0.1 and the TOF data are from |y| < 0.5. The curves are various fits to the dE/dx data for extrapolation. The quadratic sum of statistical errors and point-to-point systematic errors are plotted, but are smaller than the point size.

Mid-rapidity identified antiproton spectra in 200 GeV minimum bias d+Au (a) and 62.4 GeV central Au+Au collisions (b) measured by dE/dx together with those by TOF [46, 47]. The dE/dx data are from |y| < 0.1 and the TOF data are from |y| < 0.5. The curves are various fits to the dE/dx data for extrapolation. The quadratic sum of statistical errors and point-to-point systematic errors are plotted, but are smaller than the point size.

Mid-rapidity identified antiproton spectra in 200 GeV minimum bias d+Au (a) and 62.4 GeV central Au+Au collisions (b) measured by dE/dx together with those by TOF [46, 47]. The dE/dx data are from |y| < 0.1 and the TOF data are from |y| < 0.5. The curves are various fits to the dE/dx data for extrapolation. The quadratic sum of statistical errors and point-to-point systematic errors are plotted, but are smaller than the point size.

Mid-rapidity (|y| < 0.1) identified particle spectra in d+Au collisions at 200 GeV. The p and pbar spectra are inclusive, including weak decay products. Spectra are plotted for three centrality bins and for minimum bias events. Spectra from top to bottom are for 0-20% scaled by 4, 20-40% scaled by 2, minimum bias not scaled, and 40-100% scaled by 1/2. Errors plotted are statistical and point-to-point systematic errors added in quadrature, but are smaller than the point size. The curves are the blast-wave model fits to the minimum bias data; the normalizations of the curves are fixed by the corresponding negative particle spectra.

Mid-rapidity (|y| < 0.1) identified particle spectra in d+Au collisions at 200 GeV. The p and pbar spectra are inclusive, including weak decay products. Spectra are plotted for three centrality bins and for minimum bias events. Spectra from top to bottom are for 0-20% scaled by 4, 20-40% scaled by 2, minimum bias not scaled, and 40-100% scaled by 1/2. Errors plotted are statistical and point-to-point systematic errors added in quadrature, but are smaller than the point size. The curves are the blast-wave model fits to the minimum bias data; the normalizations of the curves are fixed by the corresponding negative particle spectra.

Mid-rapidity (|y| < 0.1) identified particle spectra in d+Au collisions at 200 GeV. The p and pbar spectra are inclusive, including weak decay products. Spectra are plotted for three centrality bins and for minimum bias events. Spectra from top to bottom are for 0-20% scaled by 4, 20-40% scaled by 2, minimum bias not scaled, and 40-100% scaled by 1/2. Errors plotted are statistical and point-to-point systematic errors added in quadrature, but are smaller than the point size. The curves are the blast-wave model fits to the minimum bias data; the normalizations of the curves are fixed by the corresponding negative particle spectra.

Mid-rapidity (|y| < 0.1) identified particle spectra in Au+Au collisions at 62.4 GeV. The p and p spectra are inclusive, including weak decay products. Spectra are plotted for nine centrality bins, from top to bottom, 0-5%, 5-10%, 10-20%, 20-30%, 30-40%, 40-50%, 50-60%, 60-70%, and 70-80%. Errors plotted are statistical and point-to-point systematic errors added in quadrature, but are smaller than the point size. The curves are the blast-wave model fits to the spectra; the normalizations of the curves in (a,b) are fixed by the corresponding negative particle spectra.

Mid-rapidity (|y| < 0.1) identified particle spectra in Au+Au collisions at 62.4 GeV. The p and p spectra are inclusive, including weak decay products. Spectra are plotted for nine centrality bins, from top to bottom, 0-5%, 5-10%, 10-20%, 20-30%, 30-40%, 40-50%, 50-60%, 60-70%, and 70-80%. Errors plotted are statistical and point-to-point systematic errors added in quadrature, but are smaller than the point size. The curves are the blast-wave model fits to the spectra; the normalizations of the curves in (a,b) are fixed by the corresponding negative particle spectra.

Mid-rapidity (|y| < 0.1) identified particle spectra in Au+Au collisions at 62.4 GeV. The p and p spectra are inclusive, including weak decay products. Spectra are plotted for nine centrality bins, from top to bottom, 0-5%, 5-10%, 10-20%, 20-30%, 30-40%, 40-50%, 50-60%, 60-70%, and 70-80%. Errors plotted are statistical and point-to-point systematic errors added in quadrature, but are smaller than the point size. The curves are the blast-wave model fits to the spectra; the normalizations of the curves in (a,b) are fixed by the corresponding negative particle spectra.

Mid-rapidity (|y| < 0.1) identified particle spectra in Au+Au collisions at 62.4 GeV. The p and p spectra are inclusive, including weak decay products. Spectra are plotted for nine centrality bins, from top to bottom, 0-5%, 5-10%, 10-20%, 20-30%, 30-40%, 40-50%, 50-60%, 60-70%, and 70-80%. Errors plotted are statistical and point-to-point systematic errors added in quadrature, but are smaller than the point size. The curves are the blast-wave model fits to the spectra; the normalizations of the curves in (a,b) are fixed by the corresponding negative particle spectra.

Mid-rapidity (|y| < 0.1) identified pion spectra in Au+Au collisions at 130 GeV. Spectra are plotted for eight centrality bins, from top to bottom, 0-6%, 6-11%, 11-18%, 18- 26%, 26-34%, 34-45%, 45-58%, and 58-85%. Errors plotted are statistical and point-to-point systematic errors added in quadrature, but they are smaller than the data point size. The curves are the Bose-Einstein fits to the spectra; the normalizations of the curves are fixed by the corresponding negative particle spectra.

Estimate of the product of the Bjorken energy density and the formation time ($\epsilon_{B_j}\tau$) as a function of centrality Npart. Errors shown are the quadratic sum of statistical and systematic uncertainties.

Estimate of the product of the Bjorken energy density and the formation time ($\epsilon_{B_j}\tau$) as a function of centrality Npart. Errors shown are the quadratic sum of statistical and systematic uncertainties.

Estimate of the product of the Bjorken energy density and the formation time ($\epsilon_{B_j}\tau$) as a function of centrality Npart. Errors shown are the quadratic sum of statistical and systematic uncertainties.

The K+/π+ and K−/π− ratios as a function of the collision energy in pp and central heavy-ion collisions.

The K+/K− ratio as a function of the collision energy in central heavy-ion collisions.

The K−/π− ratio as a function of the number of participants Npart in heavy-ion collisions at RHIC. Errors shown are the quadratic sum of statistical and systematic uncertainties.

The K−/π− ratio as a function of the number of participants Npart in heavy-ion collisions at RHIC. Errors shown are the quadratic sum of statistical and systematic uncertainties.

The K−/π− ratio as a function of the number of participants Npart in heavy-ion collisions at RHIC. Errors shown are the quadratic sum of statistical and systematic uncertainties.

The K−/π− ratio as a function of dNπ/dy/S⊥ in heavy-ion collisions at RHIC. Errors shown are the quadratic sum of statistical and systematic uncertainties.

The K−/π− ratio as a function of dNπ/dy/S⊥ in heavy-ion collisions at RHIC. Errors shown are the quadratic sum of statistical and systematic uncertainties.

The K−/π− ratio as a function of dNπ/dy/S⊥ in heavy-ion collisions at RHIC. Errors shown are the quadratic sum of statistical and systematic uncertainties.

Baryon chemical potential extracted for central heavy-ion collisions as a function of the collision energy. Errors shown are the total statistical and systematic errors.

The extracted chemical freeze-out temperatures for central heavy-ion collisions as a function of the collision energy. Errors shown are the total statistical and systematic errors.

The extracted chemical freeze-out temperatures for central heavy-ion collisions as a function of the collision energy. Errors shown are the total statistical and systematic errors.

The extracted kinetic freeze-out temperatures for central heavy-ion collisions as a function of the collision energy. Errors shown are the total statistical and systematic errors.

Average transverse radial flow velocity extracted from the blast-wave model for central heavy-ion collisions as a function of the collision energy. Errors shown are the total statistical and systematic errors.

Phase diagram plot of chemical freezeout temperature versus baryon chemical potential extracted from chemical equilibrium models. Errors shown are the total statistical and systematic errors.

Phase diagram plot of chemical freezeout temperature versus baryon chemical potential extracted from chemical equilibrium models. Errors shown are the total statistical and systematic errors.

Differential cross-sections obtained from the optical and MC Glauber calculations for Au+Au collisions at 200 GeV. Statistical errors are smaller than the point size.

Differential cross-sections obtained from the optical and MC Glauber calculations for Au+Au collisions at 200 GeV. Statistical errors are smaller than the point size.

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.

$\mathrm{p}/\mathrm{\pi}^{+}$ versus $p_{\mathrm{T}}$ for $\mathrm{\pi}^{+}$, $\mathrm{p}$ in $\mathrm{Au}-\mathrm{Au}$ at $\sqrt{s_{\mathrm{NN}}}=200\,\mathrm{Ge\!V}$

$\mathrm{p}/\mathrm{\pi}^{+}$ versus $p_{\mathrm{T}}$ for $\mathrm{\pi}^{+}$, $\mathrm{p}$ in $\mathrm{Au}-\mathrm{Au}$ at $\sqrt{s_{\mathrm{NN}}}=200\,\mathrm{Ge\!V}$

$\mathrm{p}/\mathrm{\pi}^{+}$ versus $p_{\mathrm{T}}$ for $\mathrm{\pi}^{+}$, $\mathrm{p}$ in $\mathrm{Au}-\mathrm{Au}$ at $\sqrt{s_{\mathrm{NN}}}=200\,\mathrm{Ge\!V}$

$\mathrm{p}/\mathrm{\pi}^{+}$ versus $p_{\mathrm{T}}$ for $\mathrm{\pi}^{+}$, $\mathrm{p}$ in $\mathrm{p}\mathrm{p}$ at $\sqrt{s}=200\,\mathrm{Ge\!V}$

$\mathrm{p}/\mathrm{\pi}^{+}$ versus $p_{\mathrm{T}}$ for $\mathrm{\pi}^{+}$, $\mathrm{p}$ in $\mathrm{p}\mathrm{p}$ at $\sqrt{s}=200\,\mathrm{Ge\!V}$

$\mathrm{p}/\mathrm{\pi}^{+}$ versus $p_{\mathrm{T}}$ for $\mathrm{\pi}^{+}$, $\mathrm{p}$ in $\mathrm{p}\mathrm{p}$ at $\sqrt{s}=200\,\mathrm{Ge\!V}$

$\mathrm{p}/\mathrm{\pi}^{+}$ versus $p_{\mathrm{T}}$ for $\mathrm{\pi}^{+}$, $\mathrm{p}$ in $\mathrm{p}\mathrm{p}$ at $\sqrt{s}=200\,\mathrm{Ge\!V}$

$\mathrm{p}/\mathrm{\pi}^{+}$ versus $p_{\mathrm{T}}$ for $\mathrm{\pi}^{+}$, $\mathrm{p}$ in $\mathrm{p}\mathrm{p}$ at $\sqrt{s}=200\,\mathrm{Ge\!V}$

$\overline{\mathrm{p}}/\mathrm{\pi}^{-}$ versus $p_{\mathrm{T}}$ for $\mathrm{\pi}^{-}$, $\overline{\mathrm{p}}$ in $\mathrm{Au}-\mathrm{Au}$ at $\sqrt{s_{\mathrm{NN}}}=200\,\mathrm{Ge\!V}$

$\overline{\mathrm{p}}/\mathrm{\pi}^{-}$ versus $p_{\mathrm{T}}$ for $\mathrm{\pi}^{-}$, $\overline{\mathrm{p}}$ in $\mathrm{Au}-\mathrm{Au}$ at $\sqrt{s_{\mathrm{NN}}}=200\,\mathrm{Ge\!V}$

$\overline{\mathrm{p}}/\mathrm{\pi}^{-}$ versus $p_{\mathrm{T}}$ for $\mathrm{\pi}^{-}$, $\overline{\mathrm{p}}$ in $\mathrm{Au}-\mathrm{Au}$ at $\sqrt{s_{\mathrm{NN}}}=200\,\mathrm{Ge\!V}$

$\overline{\mathrm{p}}/\mathrm{\pi}^{-}$ versus $p_{\mathrm{T}}$ for $\mathrm{\pi}^{-}$, $\overline{\mathrm{p}}$ in $\mathrm{Au}-\mathrm{Au}$ at $\sqrt{s_{\mathrm{NN}}}=200\,\mathrm{Ge\!V}$

$\overline{\mathrm{p}}/\mathrm{\pi}^{-}$ versus $p_{\mathrm{T}}$ for $\mathrm{\pi}^{-}$, $\overline{\mathrm{p}}$ in $\mathrm{Au}-\mathrm{Au}$ at $\sqrt{s_{\mathrm{NN}}}=200\,\mathrm{Ge\!V}$

$\overline{\mathrm{p}}/\mathrm{\pi}^{-}$ versus $p_{\mathrm{T}}$ for $\mathrm{\pi}^{-}$, $\overline{\mathrm{p}}$ in $\mathrm{Au}-\mathrm{Au}$ at $\sqrt{s_{\mathrm{NN}}}=200\,\mathrm{Ge\!V}$

$\overline{\mathrm{p}}/\mathrm{\pi}^{-}$ versus $p_{\mathrm{T}}$ for $\mathrm{\pi}^{-}$, $\overline{\mathrm{p}}$ in $\mathrm{p}\mathrm{p}$ at $\sqrt{s}=200\,\mathrm{Ge\!V}$

$\overline{\mathrm{p}}/\mathrm{\pi}^{-}$ versus $p_{\mathrm{T}}$ for $\mathrm{\pi}^{-}$, $\overline{\mathrm{p}}$ in $\mathrm{p}\mathrm{p}$ at $\sqrt{s}=200\,\mathrm{Ge\!V}$

$\overline{\mathrm{p}}/\mathrm{\pi}^{-}$ versus $p_{\mathrm{T}}$ for $\mathrm{\pi}^{-}$, $\overline{\mathrm{p}}$ in $\mathrm{p}\mathrm{p}$ at $\sqrt{s}=200\,\mathrm{Ge\!V}$

$\overline{\mathrm{p}}/\mathrm{\pi}^{-}$ versus $p_{\mathrm{T}}$ for $\mathrm{\pi}^{-}$, $\overline{\mathrm{p}}$ in $\mathrm{p}\mathrm{p}$ at $\sqrt{s}=200\,\mathrm{Ge\!V}$

$\overline{\mathrm{p}}/\mathrm{\pi}^{-}$ versus $p_{\mathrm{T}}$ for $\mathrm{\pi}^{-}$, $\overline{\mathrm{p}}$ in $\mathrm{p}\mathrm{p}$ at $\sqrt{s}=200\,\mathrm{Ge\!V}$

$\mathrm{p}/\mathrm{\pi}^{+}$ versus $p_{\mathrm{T}}$ for $\mathrm{\pi}^{+}$, $\mathrm{p}$ in $\mathrm{Au}-\mathrm{Au}$ at $\sqrt{s_{\mathrm{NN}}}=200\,\mathrm{Ge\!V}$

$\mathrm{p}/\mathrm{\pi}^{+}$ versus $p_{\mathrm{T}}$ for $\mathrm{\pi}^{+}$, $\mathrm{p}$ in $\mathrm{Au}-\mathrm{Au}$ at $\sqrt{s_{\mathrm{NN}}}=62.4\,\mathrm{Ge\!V}$

$\mathrm{p}/\mathrm{\pi}^{+}$ versus $p_{\mathrm{T}}$ for $\mathrm{\pi}^{+}$, $\mathrm{p}$ in $\mathrm{p}\mathrm{p}$ at $\sqrt{s}=62.4\,\mathrm{Ge\!V}$

$\mathrm{p}/\mathrm{\pi}^{+}$ versus $p_{\mathrm{T}}$ for $\mathrm{\pi}^{+}$, $\mathrm{p}$ in $\mathrm{Au}-\mathrm{Au}$ at $\sqrt{s_{\mathrm{NN}}}=62.4\,\mathrm{Ge\!V}$

$\mathrm{p}/\mathrm{\pi}^{+}$ versus $p_{\mathrm{T}}$ for $\mathrm{\pi}^{+}$, $\mathrm{p}$ in $\mathrm{Au}-\mathrm{Au}$ at $\sqrt{s_{\mathrm{NN}}}=62.4\,\mathrm{Ge\!V}$

$\mathrm{p}/\mathrm{\pi}^{+}$ versus $p_{\mathrm{T}}$ for $\mathrm{\pi}^{+}$, $\mathrm{p}$ in $\mathrm{Au}-\mathrm{Au}$ at $\sqrt{s_{\mathrm{NN}}}=62.4\,\mathrm{Ge\!V}$

$\overline{\mathrm{p}}/\mathrm{\pi}^{-}$, $\mathrm{p}/\mathrm{\pi}^{+}$ versus $p_{\mathrm{T}}$ for $\mathrm{\pi}^{-}$, $\mathrm{\pi}^{+}$, $\mathrm{p}$, $\overline{\mathrm{p}}$ in $\mathrm{Au}-\mathrm{Au}$ at $\sqrt{s_{\mathrm{NN}}}=200\,\mathrm{Ge\!V}$

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}$.

$\frac{1}{2\pi p_{\mathrm{T}}}\frac{\mathrm{d}^2N}{\mathrm{d}p_{\mathrm{T}}\mathrm{d}y}$ versus $p_{\mathrm{T}}$ for $\mathrm{\pi}^{+}$ in $\mathrm{Au}-\mathrm{Au}$ at $\sqrt{s_{\mathrm{NN}}}=62.4\,\mathrm{Ge\!V}$ near $y=0.9-1.1$ for $0-10$% central

$\frac{1}{2\pi p_{\mathrm{T}}}\frac{\mathrm{d}^2N}{\mathrm{d}p_{\mathrm{T}}\mathrm{d}y}$ versus $p_{\mathrm{T}}$ for $\mathrm{\pi}^{+}$ in $\mathrm{Au}-\mathrm{Au}$ at $\sqrt{s_{\mathrm{NN}}}=62.4\,\mathrm{Ge\!V}$ near $y=3.05-3.15$ for $0-10$% central

$\frac{1}{2\pi p_{\mathrm{T}}}\frac{\mathrm{d}^2N}{\mathrm{d}p_{\mathrm{T}}\mathrm{d}y}$ versus $p_{\mathrm{T}}$ for $\mathrm{\pi}^{+}$ in $\mathrm{Au}-\mathrm{Au}$ at $\sqrt{s_{\mathrm{NN}}}=62.4\,\mathrm{Ge\!V}$ near $y=3.15-3.25$ for $0-10$% central

$\frac{1}{2\pi p_{\mathrm{T}}}\frac{\mathrm{d}^2N}{\mathrm{d}p_{\mathrm{T}}\mathrm{d}y}$ versus $p_{\mathrm{T}}$ for $\mathrm{\pi}^{+}$ in $\mathrm{Au}-\mathrm{Au}$ at $\sqrt{s_{\mathrm{NN}}}=62.4\,\mathrm{Ge\!V}$ near $y=3.4-3.6$ for $0-10$% central

$\frac{1}{2\pi p_{\mathrm{T}}}\frac{\mathrm{d}^2N}{\mathrm{d}p_{\mathrm{T}}\mathrm{d}y}$ versus $p_{\mathrm{T}}$ for $\mathrm{\pi}^{-}$ in $\mathrm{Au}-\mathrm{Au}$ at $\sqrt{s_{\mathrm{NN}}}=62.4\,\mathrm{Ge\!V}$ near $y=-0.2-0.0$ for $0-10$% central

$\frac{1}{2\pi p_{\mathrm{T}}}\frac{\mathrm{d}^2N}{\mathrm{d}p_{\mathrm{T}}\mathrm{d}y}$ versus $p_{\mathrm{T}}$ for $\mathrm{\pi}^{-}$ in $\mathrm{Au}-\mathrm{Au}$ at $\sqrt{s_{\mathrm{NN}}}=62.4\,\mathrm{Ge\!V}$ near $y=0.0-0.2$ for $0-10$% central

$\frac{1}{2\pi p_{\mathrm{T}}}\frac{\mathrm{d}^2N}{\mathrm{d}p_{\mathrm{T}}\mathrm{d}y}$ versus $p_{\mathrm{T}}$ for $\mathrm{\pi}^{-}$ in $\mathrm{Au}-\mathrm{Au}$ at $\sqrt{s_{\mathrm{NN}}}=62.4\,\mathrm{Ge\!V}$ near $y=0.7-0.9$ for $0-10$% central

$\frac{1}{2\pi p_{\mathrm{T}}}\frac{\mathrm{d}^2N}{\mathrm{d}p_{\mathrm{T}}\mathrm{d}y}$ versus $p_{\mathrm{T}}$ for $\mathrm{\pi}^{-}$ in $\mathrm{Au}-\mathrm{Au}$ at $\sqrt{s_{\mathrm{NN}}}=62.4\,\mathrm{Ge\!V}$ near $y=0.9-1.1$ for $0-10$% central

$\frac{1}{2\pi p_{\mathrm{T}}}\frac{\mathrm{d}^2N}{\mathrm{d}p_{\mathrm{T}}\mathrm{d}y}$ versus $p_{\mathrm{T}}$ for $\mathrm{\pi}^{-}$ in $\mathrm{Au}-\mathrm{Au}$ at $\sqrt{s_{\mathrm{NN}}}=62.4\,\mathrm{Ge\!V}$ near $y=3.05-3.15$ for $0-10$% central

$\frac{1}{2\pi p_{\mathrm{T}}}\frac{\mathrm{d}^2N}{\mathrm{d}p_{\mathrm{T}}\mathrm{d}y}$ versus $p_{\mathrm{T}}$ for $\mathrm{\pi}^{-}$ in $\mathrm{Au}-\mathrm{Au}$ at $\sqrt{s_{\mathrm{NN}}}=62.4\,\mathrm{Ge\!V}$ near $y=3.15-3.25$ for $0-10$% central

$\frac{1}{2\pi p_{\mathrm{T}}}\frac{\mathrm{d}^2N}{\mathrm{d}p_{\mathrm{T}}\mathrm{d}y}$ versus $p_{\mathrm{T}}$ for $\mathrm{\pi}^{-}$ in $\mathrm{Au}-\mathrm{Au}$ at $\sqrt{s_{\mathrm{NN}}}=62.4\,\mathrm{Ge\!V}$ near $y=3.4-3.6$ for $0-10$% central

$\frac{1}{2\pi p_{\mathrm{T}}}\frac{\mathrm{d}^2N}{\mathrm{d}p_{\mathrm{T}}\mathrm{d}y}$ versus $p_{\mathrm{T}}$ for $\mathrm{K}^{+}$ in $\mathrm{Au}-\mathrm{Au}$ at $\sqrt{s_{\mathrm{NN}}}=62.4\,\mathrm{Ge\!V}$ near $y=-0.15-0.15$ for $0-10$% central

$\frac{1}{2\pi p_{\mathrm{T}}}\frac{\mathrm{d}^2N}{\mathrm{d}p_{\mathrm{T}}\mathrm{d}y}$ versus $p_{\mathrm{T}}$ for $\mathrm{K}^{+}$ in $\mathrm{Au}-\mathrm{Au}$ at $\sqrt{s_{\mathrm{NN}}}=62.4\,\mathrm{Ge\!V}$ near $y=0.6-0.8$ for $0-10$% central

$\frac{1}{2\pi p_{\mathrm{T}}}\frac{\mathrm{d}^2N}{\mathrm{d}p_{\mathrm{T}}\mathrm{d}y}$ versus $p_{\mathrm{T}}$ for $\mathrm{K}^{+}$ in $\mathrm{Au}-\mathrm{Au}$ at $\sqrt{s_{\mathrm{NN}}}=62.4\,\mathrm{Ge\!V}$ near $y=0.8-1.0$ for $0-10$% central

$\frac{1}{2\pi p_{\mathrm{T}}}\frac{\mathrm{d}^2N}{\mathrm{d}p_{\mathrm{T}}\mathrm{d}y}$ versus $p_{\mathrm{T}}$ for $\mathrm{K}^{+}$ in $\mathrm{Au}-\mathrm{Au}$ at $\sqrt{s_{\mathrm{NN}}}=62.4\,\mathrm{Ge\!V}$ near $y=2.55-2.65$ for $0-10$% central

$\frac{1}{2\pi p_{\mathrm{T}}}\frac{\mathrm{d}^2N}{\mathrm{d}p_{\mathrm{T}}\mathrm{d}y}$ versus $p_{\mathrm{T}}$ for $\mathrm{K}^{+}$ in $\mathrm{Au}-\mathrm{Au}$ at $\sqrt{s_{\mathrm{NN}}}=62.4\,\mathrm{Ge\!V}$ near $y=2.65-2.75$ for $0-10$% central

$\frac{1}{2\pi p_{\mathrm{T}}}\frac{\mathrm{d}^2N}{\mathrm{d}p_{\mathrm{T}}\mathrm{d}y}$ versus $p_{\mathrm{T}}$ for $\mathrm{K}^{+}$ in $\mathrm{Au}-\mathrm{Au}$ at $\sqrt{s_{\mathrm{NN}}}=62.4\,\mathrm{Ge\!V}$ near $y=3.15-3.25$ for $0-10$% central

$\frac{1}{2\pi p_{\mathrm{T}}}\frac{\mathrm{d}^2N}{\mathrm{d}p_{\mathrm{T}}\mathrm{d}y}$ versus $p_{\mathrm{T}}$ for $\mathrm{K}^{+}$ in $\mathrm{Au}-\mathrm{Au}$ at $\sqrt{s_{\mathrm{NN}}}=62.4\,\mathrm{Ge\!V}$ near $y=3.25-3.35$ for $0-10$% central

$\frac{1}{2\pi p_{\mathrm{T}}}\frac{\mathrm{d}^2N}{\mathrm{d}p_{\mathrm{T}}\mathrm{d}y}$ versus $p_{\mathrm{T}}$ for $\mathrm{K}^{-}$ in $\mathrm{Au}-\mathrm{Au}$ at $\sqrt{s_{\mathrm{NN}}}=62.4\,\mathrm{Ge\!V}$ near $y=-0.15-0.15$ for $0-10$% central

$\frac{1}{2\pi p_{\mathrm{T}}}\frac{\mathrm{d}^2N}{\mathrm{d}p_{\mathrm{T}}\mathrm{d}y}$ versus $p_{\mathrm{T}}$ for $\mathrm{K}^{-}$ in $\mathrm{Au}-\mathrm{Au}$ at $\sqrt{s_{\mathrm{NN}}}=62.4\,\mathrm{Ge\!V}$ near $y=0.6-0.8$ for $0-10$% central

$\frac{1}{2\pi p_{\mathrm{T}}}\frac{\mathrm{d}^2N}{\mathrm{d}p_{\mathrm{T}}\mathrm{d}y}$ versus $p_{\mathrm{T}}$ for $\mathrm{K}^{-}$ in $\mathrm{Au}-\mathrm{Au}$ at $\sqrt{s_{\mathrm{NN}}}=62.4\,\mathrm{Ge\!V}$ near $y=0.8-1.0$ for $0-10$% central

$\frac{1}{2\pi p_{\mathrm{T}}}\frac{\mathrm{d}^2N}{\mathrm{d}p_{\mathrm{T}}\mathrm{d}y}$ versus $p_{\mathrm{T}}$ for $\mathrm{K}^{-}$ in $\mathrm{Au}-\mathrm{Au}$ at $\sqrt{s_{\mathrm{NN}}}=62.4\,\mathrm{Ge\!V}$ near $y=2.55-2.65$ for $0-10$% central

$\frac{1}{2\pi p_{\mathrm{T}}}\frac{\mathrm{d}^2N}{\mathrm{d}p_{\mathrm{T}}\mathrm{d}y}$ versus $p_{\mathrm{T}}$ for $\mathrm{K}^{-}$ in $\mathrm{Au}-\mathrm{Au}$ at $\sqrt{s_{\mathrm{NN}}}=62.4\,\mathrm{Ge\!V}$ near $y=2.65-2.75$ for $0-10$% central

$\frac{1}{2\pi p_{\mathrm{T}}}\frac{\mathrm{d}^2N}{\mathrm{d}p_{\mathrm{T}}\mathrm{d}y}$ versus $p_{\mathrm{T}}$ for $\mathrm{K}^{-}$ in $\mathrm{Au}-\mathrm{Au}$ at $\sqrt{s_{\mathrm{NN}}}=62.4\,\mathrm{Ge\!V}$ near $y=3.15-3.25$ for $0-10$% central

$\frac{1}{2\pi p_{\mathrm{T}}}\frac{\mathrm{d}^2N}{\mathrm{d}p_{\mathrm{T}}\mathrm{d}y}$ versus $p_{\mathrm{T}}$ for $\mathrm{K}^{-}$ in $\mathrm{Au}-\mathrm{Au}$ at $\sqrt{s_{\mathrm{NN}}}=62.4\,\mathrm{Ge\!V}$ near $y=3.25-3.35$ for $0-10$% central

$\frac{\mathrm{d}N}{\mathrm{d}y}$ versus $y$ for $\mathrm{\pi}^{+}$ in $\mathrm{Au}-\mathrm{Au}$ at $\sqrt{s_{\mathrm{NN}}}=62.4\,\mathrm{Ge\!V}$ for $0-10$% central

$\frac{\mathrm{d}N}{\mathrm{d}y}$ versus $y$ for $\mathrm{\pi}^{-}$ in $\mathrm{Au}-\mathrm{Au}$ at $\sqrt{s_{\mathrm{NN}}}=62.4\,\mathrm{Ge\!V}$ for $0-10$% central

$\frac{\mathrm{d}N}{\mathrm{d}y}$ versus $y$ for $\mathrm{K}^{+}$ in $\mathrm{Au}-\mathrm{Au}$ at $\sqrt{s_{\mathrm{NN}}}=62.4\,\mathrm{Ge\!V}$ for $0-10$% central

$\frac{\mathrm{d}N}{\mathrm{d}y}$ versus $y$ for $\mathrm{K}^{-}$ in $\mathrm{Au}-\mathrm{Au}$ at $\sqrt{s_{\mathrm{NN}}}=62.4\,\mathrm{Ge\!V}$ for $0-10$% central

$\mathrm{\pi}^{-}/\mathrm{\pi}^{+}$ versus $y$ for $\mathrm{\pi}^{-}$,$\mathrm{\pi}^{+}$ in $\mathrm{Au}-\mathrm{Au}$ at $\sqrt{s_{\mathrm{NN}}}=62.4\,\mathrm{Ge\!V}$ for $0-10$% central

$\mathrm{K}^{-}/\mathrm{K}^{+}$ versus $y$ for $\mathrm{K}^{-}$,$\mathrm{K}^{+}$ in $\mathrm{Au}-\mathrm{Au}$ at $\sqrt{s_{\mathrm{NN}}}=62.4\,\mathrm{Ge\!V}$ for $0-10$% central

$\overline{\mathrm{p}}/\mathrm{p}$ versus $y$ for $\overline{\mathrm{p}}$,$\mathrm{p}$ in $\mathrm{Au}-\mathrm{Au}$ at $\sqrt{s_{\mathrm{NN}}}=62.4\,\mathrm{Ge\!V}$ for $0-10$% central

$\mathrm{K}^{+}/\mathrm{\pi}^{+}$ versus $y$ for $\mathrm{K}^{+}$,$\mathrm{\pi}^{+}$ in $\mathrm{Au}-\mathrm{Au}$ at $\sqrt{s_{\mathrm{NN}}}=62.4\,\mathrm{Ge\!V}$ for $0-10$% central

$\mathrm{K}^{-}/\mathrm{\pi}^{-}$ versus $y$ for $\mathrm{K}^{-}$,$\mathrm{\pi}^{-}$ in $\mathrm{Au}-\mathrm{Au}$ at $\sqrt{s_{\mathrm{NN}}}=62.4\,\mathrm{Ge\!V}$ for $0-10$% central

$\mathrm{K}^{+}/\mathrm{K}^{-}$ versus $\overline{\mathrm{p}}/\mathrm{p}$ for $\mathrm{K}^{+}$,$\mathrm{K}^{-}$,$\overline{\mathrm{p}}$,$\mathrm{p}$ in $\mathrm{Au}-\mathrm{Au}$ at $\sqrt{s_{\mathrm{NN}}}=62.4\,\mathrm{Ge\!V}$ for $0-10$% central

$\mathrm{K}^{+}/\mathrm{\pi}^{+}$,$\mathrm{K}^{-}/\mathrm{\pi}^{-}$ versus $\overline{\mathrm{p}}/\mathrm{p}$ for $\mathrm{K}^{+}$,$\mathrm{K}^{-}$,$\mathrm{\pi}^{+}$,$\mathrm{\pi}^{-}$,$\overline{\mathrm{p}}$,$\mathrm{p}$ in $\mathrm{Au}-\mathrm{Au}$ at $\sqrt{s_{\mathrm{NN}}}=62.4\,\mathrm{Ge\!V}$ for $0-10$% central

$T_{\mathrm{K^{+}}}$,$T_{\mathrm{K^{-}}}$ versus $\overline{\mathrm{p}}/\mathrm{p}$ for $\mathrm{K}^{+}$,$\mathrm{K}^{-}$,$\mathrm{p}$,$\overline{\mathrm{p}}$ in $\mathrm{Au}-\mathrm{Au}$ at $\sqrt{s_{\mathrm{NN}}}=62.4\,\mathrm{Ge\!V}$ for $0-10$% central

(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).

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.