The top panel shows the measured $R_{side}$ from identical pions for PHENIX. The bottom panel shows the ratio $R_{out}/R_{side}$ as a function of $k_T$. Longitudinal Co-Moving System (LCMS) frame for $\pi^-$

The top panel shows the measured $R_{side}$ from identical pions for PHENIX. The bottom panel shows the ratio $R_{out}/R_{side}$ as a function of $k_T$. Pair Center-of-Mass System (PCMS) frame for $\pi^+$

The top panel shows the measured $R_{side}$ from identical pions for PHENIX. The bottom panel shows the ratio $R_{out}/R_{side}$ as a function of $k_T$. Pair Center-of-Mass System (PCMS) frame for $\pi^-$

The yields per event at mid-rapidity for neutral pions as a function of $p_T$ for 60-80% from the PbGl detector.

The yields per event at mid-rapidity for charged hadrons as a function of $p_T$ for 0-10%.

The yields per event at mid-rapidity for charged hadrons as a functon of $p_T$ for 60-80%.

The ratio $R_{AA}$ for neutral pions in central Au+Au collisions from the PbSc detector.

The ratio $R_{AA}$ for neutral pions in central Au+Au collisions from the PbGl detector.

The ratio $R_{AA}$ for charged hadrons in central Au+Au collisions.

Ratio of yield per event in central vs peripheral Au+Au collisions, with each divided by $N_{binary}$ for that class.

Ratio of yield per event in central vs peripheral Au+Au collisions, with each divided by N_binary for that class.

Ratio of yield her event in central vs peripheral Au+Au collisions, with each divided by N_binary for that class.

Integrated per-trigger yields as a function of associated particle angle relative to $\psi_2$ event plane $\phi^a - \psi_2$ integrated over the near side $|\Delta\phi|<\pi/4$.

Integrated per-trigger yields as a function of associated particle angle relative to $\psi_2$ event plane $\phi^a - \psi_2$ integrated over the away side $|\Delta\phi|<\pi/4$.

Integrated per-trigger yields as a function of associated particle angle relative to $\psi_2$ event plane for near-side $|\Delta\phi|<\pi/4$ and for away-side $|\Delta\phi|<\pi/4$.

Integrated per-trigger yields as a function of associated particle angle relative to $\psi_2$ event plane for near-side $|\Delta\phi|<\pi/4$ and for away-side $|\Delta\phi|<\pi/4$.

Per-trigger yields $Y(\Delta\phi)$ of dihadrons pairs measured in Au$+$Au collisions at $\sqrt{s_{_{NN}}}=200~{\rm GeV}$ after subtracting the underlying event model with several $p_T$ selections. Systematic uncertainties due to track matching are given.

Per-trigger yields $Y(\Delta\phi)$ of dihadrons pairs measured in Au$+$Au collisions at $\sqrt{s_{_{NN}}}=200~{\rm GeV}$ after subtracting the underlying event model with several $p_T$ selections. Systematic uncertainties due to track matching are given.

$E\frac{dN^3}{dp^3}$ vs. $p_T$, 0% to 20% centrality $Au+Au$. 90% Limit on 18-20 and 20-22 GeV/c bins.

$E\frac{dN^3}{dp^3}$ vs. $p_T$, 20% to 40% centrality $Au+Au$. 90% Limit for 20-22 GeV/c bin.

$E\frac{dN^3}{dp^3}$ vs. $p_T$, 40% to 60% centrality $Au+Au$. 90% Limit for 18-20 GeV/c bin.

$E\frac{dN^3}{dp^3}$ vs. $p_T$, 20% to 60% centrality $Au+Au$. 90% Limit for 18-20 GeV/c bin.

$E\frac{dN^3}{dp^3}$ vs. $p_T$, 60% to 92% centrality $Au+Au$. 90% Limit for 18-20 GeV/c bin.

$E\frac{dN^3}{dp^3}$ vs. $p_T$, 0% to 92% centrality $Au+Au$.

$E\frac{dN^3}{dp^3}$ vs. $p_T$, $p+p$.

${\eta}$ $R_{AA}$, 0% to 5% centrality $Au+Au$.

${\eta}$ $R_{AA}$, 0% to 10% centrality $Au+Au$.

${\eta}$ $R_{AA}$, 10% to 20% centrality $Au+Au$.

${\eta}$ $R_{AA}$, 0% to 20% centrality $Au+Au$.

${\eta}$ $R_{AA}$, 20% to 40% centrality $Au+Au$.

${\eta}$ $R_{AA}$, 40% to 60% centrality $Au+Au$.

${\eta}$ $R_{AA}$, 20% to 60% centrality $Au+Au$.

${\eta}$ $R_{AA}$, 60% to 92% centrality $Au+Au$.

${\eta}$ $R_{AA}$, 0% to 92% centrality $Au+Au$.

${\pi^{0}}$ $R_{AA}$, 0% to 92% centrality $Au+Au$.

${\eta}$ ${R_{AA}}$, 0% to 92% centrality $Au+Au$.

${\eta}$ $R_{AA}$ vs centrality ($p_T>$20 GeV/c).

${\eta}$ $R_{AA}$ vs centrality ($p_T>$5 GeV/c).

Slope ${\eta}$ $R_{AA}$ ($5<p_T<$ 18 GeV/c) vs centrality.

$\chi^{2}/{NDF}$, Prob%, ${\eta}$ $R_{AA}$ ($5<p_T<$ 18 GeV/c) Prob, $R_{AA}$ ($p_T>$ 5 GeV/c) Prob vs centrality

Bertsch-Pratt radius parameters versus the cube root of the number of participants ($N^{1/3}_{part}$) for $\pi^+\pi^+$ (blue squares) and $\pi^-\pi^-$ (red triangle) for the fits to a core-halo structure by Eq. (2) with statistical error bars and systematic error bands, in the 0.2 GeV/$c$ < $k_T$ < 2.0 GeV/$c$ range, with $<k_T>$ ~ 0.45 GeV/$c$. The solid line shows to $p_0$+$p_1$*$N^{1/3}_{part}$ for both signs.

Bertsch-Pratt radius parameters versus the cube root of the number of participants ($N^{1/3}_{part}$) for $\pi^+\pi^+$ (blue squares) and $\pi^-\pi^-$ (red triangle) for the fits to a core-halo structure by Eq. (2) with statistical error bars and systematic error bands, in the 0.2 GeV/$c$ < $k_T$ < 2.0 GeV/$c$ range, with $<k_T>$ ~ 0.45 GeV/$c$. The solid line shows to $p_0$+$p_1$*$N^{1/3}_{part}$ for both signs.

$\pi^0$ spectra from figure 3a from 0-20% central Cu+Au collisions. Type A uncertainties are uncorrelated point-to-point. Type B uncertainties are correlated point-to-point. Type C uncertainties affect the scale of the data.

$\pi^0$ spectra from figure 3a from 20-40% semi-central Cu+Au collisions. Type A uncertainties are uncorrelated point-to-point. Type B uncertainties are correlated point-to-point. Type C uncertainties affect the scale of the data.

$\pi^0$ spectra from figure 3a from 40-60% semi-peripheral Cu+Au collisions. Type A uncertainties are uncorrelated point-to-point. Type B uncertainties are correlated point-to-point. Type C uncertainties affect the scale of the data.

$\pi^0$ spectra supplemental to figure 3a from 60-90% peripheral Cu+Au collisions. Type A uncertainties are uncorrelated point-to-point. Type B uncertainties are correlated point-to-point. Type C uncertainties affect the scale of the data.

Ratio of $\pi^0$ spectra from the PbSc to the combined spectra from figure 3b from 0-20% central Cu+Au collisions. Type A uncertainties are uncorrelated point-to-point. Systematic noncorr. uncertainties are uncorrelated between PbSc and PbGl subsystems.

Ratio of $\pi^0$ spectra from the PbGl to the combined spectra from figure 3b from 0-20% central Cu+Au collisions. Type A uncertainties are uncorrelated point-to-point. Systematic noncorr. uncertainties are uncorrelated between PbSc and PbGl subsystems.

Ratio of $\pi^0$ spectra from the PbSc to the combined spectra from figure 3c from 20-40% semi-central Cu+Au collisions. Type A uncertainties are uncorrelated point-to-point. Systematic noncorr. uncertainties are uncorrelated between PbSc and PbGl subsystems.

Ratio of $\pi^0$ spectra from the PbGl to the combined spectra from figure 3c from 20-40% semi-central Cu+Au collisions. Type A uncertainties are uncorrelated point-to-point. Systematic noncorr. uncertainties are uncorrelated between PbSc and PbGl subsystems.

Ratio of $\pi^0$ spectra from the PbSc to the combined spectra from figure 3d from 40-60% semi-peripheral Cu+Au collisions. Type A uncertainties are uncorrelated point-to-point. Systematic noncorr. uncertainties are uncorrelated between PbSc and PbGl subsystems.

Ratio of $\pi^0$ spectra from the PbGl to the combined spectra from figure 3d from 40-60% semi-peripheral Cu+Au collisions. Type A uncertainties are uncorrelated point-to-point. Systematic noncorr. uncertainties are uncorrelated between PbSc and PbGl subsystems.

$\eta$ spectra from figure 3e from minimum bias Cu+Au collisions. Type A uncertainties are uncorrelated point-to-point. Type B uncertainties are correlated point-to-point. Type C uncertainties affect the scale of the data.

$\eta$ spectra from figure 3e from 0-20% central Cu+Au collisions. Type A uncertainties are uncorrelated point-to-point. Type B uncertainties are correlated point-to-point. Type C uncertainties affect the scale of the data.

$\eta$ spectra from figure 3e from 20-40% semi-central Cu+Au collisions. Type A uncertainties are uncorrelated point-to-point. Type B uncertainties are correlated point-to-point. Type C uncertainties affect the scale of the data.

$\eta$ spectra from figure 3e from 40-60% semi-peripheral Cu+Au collisions. Type A uncertainties are uncorrelated point-to-point. Type B uncertainties are correlated point-to-point. Type C uncertainties affect the scale of the data.

$\eta$ spectra supplemental to figure 3e from 60-90% peripheral Cu+Au collisions. Type A uncertainties are uncorrelated point-to-point. Type B uncertainties are correlated point-to-point. Type C uncertainties affect the scale of the data.

Ratio of $\eta$ spectra from the PbSc to the combined spectra from figure 3f from 0-20% central Cu+Au collisions. Type A uncertainties are uncorrelated point-to-point. Systematic noncorr. uncertainties are uncorrelated between PbSc and PbGl subsystems.

Ratio of $\eta$ spectra from the PbGl to the combined spectra from figure 3f from 0-20% central Cu+Au collisions. Type A uncertainties are uncorrelated point-to-point. Systematic noncorr. uncertainties are uncorrelated between PbSc and PbGl subsystems.

Ratio of $\eta$ spectra from the PbSc to the combined spectra from figure 3g from 20-40% semi-central Cu+Au collisions. Type A uncertainties are uncorrelated point-to-point. Systematic noncorr. uncertainties are uncorrelated between PbSc and PbGl subsystems.

Ratio of $\eta$ spectra from the PbGl to the combined spectra from figure 3g from 20-40% semi-central Cu+Au collisions. Type A uncertainties are uncorrelated point-to-point. Systematic noncorr. uncertainties are uncorrelated between PbSc and PbGl subsystems.

Ratio of $\eta$ spectra from the PbSc to the combined spectra from figure 3h from 40-60% semi-peripheral Cu+Au collisions. Type A uncertainties are uncorrelated point-to-point. Systematic noncorr. uncertainties are uncorrelated between PbSc and PbGl subsystems.

Ratio of $\eta$ spectra from the PbGl to the combined spectra from figure 3h from 40-60% semi-peripheral Au+Cu collisions. Type A uncertainties are uncorrelated point-to-point. Systematic noncorr. uncertainties are uncorrelated between PbSc and PbGl subsystems.

$\eta/\pi^0$ ratio from figure 4 from minimum bias Cu+Au collisions. Type A uncertainties are uncorrelated point-to-point. Type B uncertainties are correlated point-to-point.

$\eta/\pi^0$ ratio from figure 4 from 0-20% central Cu+Au collisions. Type A uncertainties are uncorrelated point-to-point. Type B uncertainties are correlated point-to-point.

$\eta/\pi^0$ ratio from figure 4 from 20-40% semi-central Cu+Au collisions. Type A uncertainties are uncorrelated point-to-point. Type B uncertainties are correlated point-to-point.

$\eta/\pi^0$ ratio from figure 4 from 40-60% semi-central Cu+Au collisions. Type A uncertainties are uncorrelated point-to-point. Type B uncertainties are correlated point-to-point.

$\eta/\pi^0$ ratio supplemental to figure 4 from 60-90% peripheral Cu+Au collisions. Type A uncertainties are uncorrelated point-to-point. Type B uncertainties are correlated point-to-point.

$\pi^0$ $R_{AB}$ from figure 5a from 0-10% central Cu+Au collisions. Type A uncertainties are uncorrelated point-to-point. Type B uncertainties are correlated point-to-point. Type C uncertainties affect the scale of the data.

$\pi^0$ $R_{AB}$ from figure 5b from 10-20% central Cu+Au collisions. Type A uncertainties are uncorrelated point-to-point. Type B uncertainties are correlated point-to-point. Type C uncertainties affect the scale of the data.

$\pi^0$ $R_{AB}$ from figure 5c from 0-20% central Cu+Au collisions. Type A uncertainties are uncorrelated point-to-point. Type B uncertainties are correlated point-to-point. Type C uncertainties affect the scale of the data.

$\eta$ $R_{AB}$ from figure 5a from 0-20% central Cu+Au collisions. Type A uncertainties are uncorrelated point-to-point. Type B uncertainties are correlated point-to-point. Type C uncertainties affect the scale of the data.

$\pi^0$ $R_{AB}$ from figure 5d from 20-40% semi-central Cu+Au collisions. Type A uncertainties are uncorrelated point-to-point. Type B uncertainties are correlated point-to-point. Type C uncertainties affect the scale of the data.

$\eta$ $R_{AB}$ from figure 5b from 20-40% semi-central Cu+Au collisions. Type A uncertainties are uncorrelated point-to-point. Type B uncertainties are correlated point-to-point. Type C uncertainties affect the scale of the data.

$\pi^0$ $R_{AB}$ from figure 5e from 40-60% semi-peripheral Cu+Au collisions. Type A uncertainties are uncorrelated point-to-point. Type B uncertainties are correlated point-to-point. Type C uncertainties affect the scale of the data.

$\eta$ $R_{AB}$ from figure 5c from 40-60% semi-peripheral Cu+Au collisions. Type A uncertainties are uncorrelated point-to-point. Type B uncertainties are correlated point-to-point. Type C uncertainties affect the scale of the data.

$\pi^0$ $R_{AB}$ from figure 5f from minimum bias Cu+Au collisions. Type A uncertainties are uncorrelated point-to-point. Type B uncertainties are correlated point-to-point. Type C uncertainties affect the scale of the data.

$\eta$ $R_{AB}$ from figure 5f from minimum bias Cu+Au collisions. Type A uncertainties are uncorrelated point-to-point. Type B uncertainties are correlated point-to-point. Type C uncertainties affect the scale of the data.

$\pi^0$ $R_{AB}$ supplemental to figure 5 from 60-90% peripheral Cu+Au collisions. Type A uncertainties are uncorrelated point-to-point. Type B uncertainties are correlated point-to-point. Type C uncertainties affect the scale of the data.

$\eta$ $R_{AB}$ supplemental to figure 5 from 60-90% peripheral Cu+Au collisions. Type A uncertainties are uncorrelated point-to-point. Type B uncertainties are correlated point-to-point. Type C uncertainties affect the scale of the data.

$\pi^0$ $R_{AB}$ supplemental to figure 6a from 20-30% central Au+Au collisions. Type A uncertainties are uncorrelated point-to-point. Type B uncertainties are correlated point-to-point. Type C uncertainties affect the scale of the data.

$\pi^0$ $R_{AB}$ supplemental to figure 6b from 0-10% central Cu+Cu collisions. Type A uncertainties are uncorrelated point-to-point. Type B uncertainties are correlated point-to-point. Type C uncertainties affect the scale of the data.

$\pi^0$ $R_{AB}$ supplemental to figure 6b from 40-50% semi-central Au+Au collisions. Type A uncertainties are uncorrelated point-to-point. Type B uncertainties are correlated point-to-point. Type C uncertainties affect the scale of the data.

$\pi^0$ for $p_T>5$ GeV/$c$ in Cu+Au collisions as a function of centrality from figure 7a. Type A uncertainties are uncorrelated point-to-point. Type B uncertainties are correlated point-to-point. Type C uncertainties affect the scale of the data.

$\pi^0$ for $p_T>5$ GeV/$c$ in Au+Au collisions as a function of centrality from figure 7a. Type A uncertainties are uncorrelated point-to-point. Type B uncertainties are correlated point-to-point. Type C uncertainties affect the scale of the data.

$\pi^0$ for $p_T>10$ GeV/$c$ in Cu+Au collisions as a function of centrality from figure 7b. Type A uncertainties are uncorrelated point-to-point. Type B uncertainties are correlated point-to-point. Type C uncertainties affect the scale of the data.

$\pi^0$ for $p_T>10$ GeV/$c$ in Au+Au collisions as a function of centrality from figure 7b. Type A uncertainties are uncorrelated point-to-point. Type B uncertainties are correlated point-to-point. Type C uncertainties affect the scale of the data.

Minimum-bias $m_T$ spectra of $\phi$ for different subsystem combinations.

$\phi$ $m_T$ spectra in different centrality bins.

Centrality dependence of $\phi$ yield at midrapidity.

Centrality dependence of $\phi$ yield at midrapidity.

Centrality dependence of the inverse slope, $T$.

Centrality dependence of particle ratios for $K^+$/$\pi^+$ and $K^-$/$\pi^-$.

Centrality dependence of particle ratios for $\phi$/$\pi$ and $\phi$/$K$.

Centrality dependence of $\phi$ yields at different collision energies.

Centrality dependence of $\phi$ yields at different collision energies.

$\phi$/$\pi$ and $\phi$/$K^-$ ratios as a function of collision energies.

$N_{coll}$ scaled central to peripheral ratio $R_{CP}$ for $\phi$. The total systematic error is 19%, but 7% is uncorrelated between the $\phi$ and $p$ $R_{CP}$.

$p_T$ spectra of charged hadrons for minimum bias collisions along with spectra for 9 centrality classes derived from the pseudo-rapidity region $|\eta|$ < 0.18. Stat. stands for statistical error, syst. stands for the systematic errors and occ. stands for occupancy error.

$p_T$ spectra of charged hadrons for minimum bias collisions along with spectra for 9 centrality classes derived from the pseudo-rapidity region $|\eta|$ < 0.18. Stat. stands for statistical error, syst. stands for the systematic errors and occ. stands for occupancy error.

Centrality dependence of $<p^{trunc}_T>$. Data points for truncated mean $p_T$ calculation for three $p_T$ ranges. The integral region is from 8 GeV/$c$ > $p_T$ > $p_{T_{min}}$ for the most peripheral bin. Since the data stop around 6.5 GeV/$c$, truncated mean $p_T$ is lower.

Ratio of charged hadron yields per nucleon-nucleon collision between central (0-10%) and peripheral (60-92%) Au + Au collisions. Stat. stands for statistical error, syst. stands for the systematic errors and occ. stands for occupancy error.

$R_{AA}$ for $(h^++h^-)/2$ and $\pi^0$ as function of $p_{T}$ for minimum bias and 9 centrality classes. Stat. stands for statistical error from Au+Au. The asymmetric uncertainty from the fit of the $\pi^0$ mass peak is denoted by fit. $\sigma_{ENBC}$ stands for error on number of binary collisions. Error from $pp$ parametrization that is not $p_T$ dependent. This error is not included in the calculation of the Upper bound of the systematic error and the Lower bound of the systematic error.

$R_{AA}$ for $(h^++h^-)/2$ and $\pi^0$ as function of $p_{T}$ for minimum bias and 9 centrality classes. Stat. stands for statistical error from Au+Au. The asymmetric uncertainty from the fit of the $\pi^0$ mass peak is denoted by fit. $\sigma_{ENBC}$ stands for error on number of binary collisions. Error from $pp$ parametrization that is not $p_T$ dependent. This error is not included in the calculation of the Upper bound of the systematic error and the Lower bound of the systematic error.

$R_{AA}$ for $(h^++h^-)/2$ and $\pi^0$ as function of $p_{T}$ for minimum bias and 9 centrality classes. Stat. stands for statistical error from Au+Au. The asymmetric uncertainty from the fit of the $\pi^0$ mass peak is denoted by fit. $\sigma_{ENBC}$ stands for error on number of binary collisions. Error from $pp$ parametrization that is not $p_T$ dependent. This error is not included in the calculation of the Upper bound of the systematic error and the Lower bound of the systematic error.

Charged hadron to $\pi^0$ ratios for minimum bias events and 9 centrality classes.

Charged hadron to $\pi^0$ ratios for minimum bias events and 9 centrality classes.

Charged hadron to $\pi^0$ ratios for minimum bias events and 9 centrality classes.

Charged hadron to $\pi^0$ ratios for minimum bias events and 9 centrality classes.

Charged hadron to $\pi^0$ ratios for minimum bias events and 9 centrality classes.

Charged hadron to $\pi^0$ ratios for minimum bias events and 9 centrality classes.

Au + Au yield integrated for $p_T$ > 4.5 (GeV/$c$) over the $N$ + $N$ yield, normalized using either $N_{coll}$ or $N_{part}$, plot as a function of $<N_{part}>$. Statistical and systematic error from Au+Au and $N$+$N$ reference includes: a) $NN$ reference error with 10.4% normalization error subtracted, this error is 11.2%, b) systematic error on Au+Au is 8.5% (correction), and c) 5% occupancy correction (it is centrality dependent). Error on $N_{coll}$+10.4% error on $N$+$N$ reference.

Au + Au yield integrated for $p_T$ > 4.5 (GeV/$c$) over the $N$ + $N$ yield, normalized using either $N_{coll}$ or $N_{part}$, plot as a function of $<N_{part}>$. Statistical and systematic error from Au+Au and $N$+$N$ reference includes: a) $NN$ reference error with 10.4% normalization error subtracted, this error is 11.2%, b) systematic error on Au+Au is 8.5% (correction), and c) 5% occupancy correction (it is centrality dependent). Error on $N_{coll}$+10.4% error on $N$+$N$ reference.

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.