$\pi^0$ spectra from figure 2a from 40-60% U+U 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 2a from 60-80% U+U 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 2b from minimum bias U+U 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 2b from 0-20% U+U 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 2b from 20-40% U+U 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 2b from 40-60% U+U 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 2b from 60-80% U+U 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$ yield to the fit from figure 2c from minimum bias U+U collisions. Sys. uncertainties include B and C uncertainties from data points.

Ratio of $\pi^0$ yield to the fit from figure 2c from 0-20% U+U collisions. Sys. uncertainties include B and C uncertainties from data points.

Ratio of $\pi^0$ yield to the fit from figure 2c from 60-80% U+U collisions. Sys. uncertainties include B and C uncertainties from data points.

Ratio of $\eta$ yield to the fit from figure 2d from minimum bias U+U collisions. Sys. uncertainties include B and C uncertainties from data points.

Ratio of $\eta$ yield to the fit from figure 2d from 0-20% U+U collisions. Sys. uncertainties include B and C uncertainties from data points.

Ratio of $\eta$ yield to the fit from figure 2d from 40-60% U+U collisions. Sys. uncertainties include B and C uncertainties from data points.

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

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

$\eta/\pi^0$ ratio from figure 3 from 20-40% U+U collisions. Type A uncertainties are uncorrelated point-to-point. Type B uncertainties are correlated point-to-point.

$\eta/\pi^0$ ratio from figure 3 from 40-60% U+U collisions. Type A uncertainties are uncorrelated point-to-point. Type B uncertainties are correlated point-to-point.

$\eta/\pi^0$ ratio from figure 3 from 60-80% U+U collisions. Type A uncertainties are uncorrelated point-to-point. Type B uncertainties are correlated point-to-point.

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

$\pi^0$ $R_{AA}$ from figure 4a from 0-20% U+U 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. Calculation is carried out for Glauber 1 set.

$\eta$ $R_{AA}$ from figure 4a from 0-20% U+U 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. Calculation is carried out for Glauber 1 set.

$\pi^0$ $R_{AA}$ from figure 4b from 20-40% U+U 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. Calculation is carried out for Glauber 1 set.

$\eta$ $R_{AA}$ from figure 4b from 20-40% U+U 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. Calculation is carried out for Glauber 1 set.

$\pi^0$ $R_{AA}$ from figure 4c from 40-60% U+U 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. Calculation is carried out for Glauber 1 set.

$\eta$ $R_{AA}$ from figure 4c from 40-60% U+U 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. Calculation is carried out for Glauber 1 set.

$\pi^0$ $R_{AA}$ from figure 4d from minimum bias U+U 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. Calculation is carried out for Glauber 1 set.

$\eta$ $R_{AA}$ from figure 4d from minimum bias U+U 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. Calculation is carried out for Glauber 1 set.

$\pi^0$ $R_{AA}$ from figure 5a from 0-20% U+U 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. Calculation is carried out for Glauber 2 set.

$\pi^0$ $R_{AA}$ from figure 5a from 0-5% 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_{AA}$ from figure 5a from 40-60% U+U 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. Calculation is carried out for Glauber 2 set.

$\pi^0$ $R_{AA}$ from figure 5a from 40-50% 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$ integrated $R_{AA}$ for $p_T>5$ GeV/$c$ in U+U collisions as a function of centrality from figure 6a. 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. Calculation is carried out for Glauber 1 set.

$\eta$ integrated $R_{AA}$ for $p_T>5$ GeV/$c$ in U+U collisions as a function of centrality from figure 6a. 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. Calculation is carried out for Glauber 1 set.

$\pi^0$ integrated $R_{AA}$ for $p_T>5$ GeV/$c$ in Au+Au collisions as a function of centrality from figure 6a. 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$ integrated $R_{AA}$ for $p_T>5$ GeV/$c$ in U+U collisions as a function of centrality from figure 6b. 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. Calculation is carried out for Glauber 2 set.

Number of participant nucleons of in dep U$+$U collisions

Number of participant nucleons of in dep U$+$U collisions

Number of participant nucleons of in dep Au$+$Au collisions

Integrated away-side $\gamma_{dir}$-h per-trigger $I_{AA}$ and $I_{dA}$ of Au+Au, d+Au, and p+p, as a function of $\xi$.

$I_{AA}$ vs $\xi$ for direct photon $p_{T}^{\gamma}$ of 5-7 GeV/$c$, 7-9 GeV/$c$, and 9-12 GeV/$c$.

$I_{AA}$ vs $\xi$ for direct photon $p_{T}^{\gamma}$ of 5-7 GeV/$c$, 7-9 GeV/$c$, and 9-12 GeV/$c$.

$I_{AA}$ vs $\xi$ for direct photon $p_{T}^{\gamma}$ of 5-7 GeV/$c$, 7-9 GeV/$c$, and 9-12 GeV/$c$.

$I_{AA}$ as a function of $\xi$ for direct photon $p_{T}^{\gamma}$ of 5-7, 7-9, and 9-12 GeV/$c$. Three away-side integration ranges are chosen to calculate the per-trigger yield and the corresponding $I_{AA}$, $|\Delta\phi-\pi|<\pi/2$, $|\Delta\phi-\pi|<\pi/3$, and $|\Delta\phi-\pi|<\pi/6$.

$I_{AA}$ as a function of $\xi$ for direct photon $p_{T}^{\gamma}$ of 5-7, 7-9, and 9-12 GeV/$c$. Three away-side integration ranges are chosen to calculate the per-trigger yield and the corresponding $I_{AA}$, $|\Delta\phi-\pi|<\pi/2$, $|\Delta\phi-\pi|<\pi/3$, and $|\Delta\phi-\pi|<\pi/6$.

$I_{AA}$ as a function of $\xi$ for direct photon $p_{T}^{\gamma}$ of 5-7, 7-9, and 9-12 GeV/$c$. Three away-side integration ranges are chosen to calculate the per-trigger yield and the corresponding $I_{AA}$, $|\Delta\phi-\pi|<\pi/2$, $|\Delta\phi-\pi|<\pi/3$, and $|\Delta\phi-\pi|<\pi/6$.

Ratios of $I_{AA}$ as a function of direct photon $p_T$ for three different away-side integration ranges.

Measured $I_{AA}$ for direct photon $p_T$ of 5-7, 7-9, and 9-12 GeV/$c$, as a function of $\xi$, are compared with theoretical model calculations.

Measured $I_{AA}$ for direct photon $p_T$ of 5-7, 7-9, and 9-12 GeV/$c$, as a function of $\xi$, are compared with theoretical model calculations.

Measured $I_{AA}$ for direct photon $p_T$ of 5-7, 7-9, and 9-12 GeV/$c$, as a function of $\xi$, are compared with theoretical model calculations.

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

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.

$\Psi_{2}$-dependent $C(\Delta\phi)$ for $2<p_{T}^{trig}<4$ GeV/c and $1<p_{T}^{assoc}<2$ GeV/c.

$\Psi_{2}$-dependent $C(\Delta\phi)$ for $2<p_{T}^{trig}<4$ GeV/c and $2<p_{T}^{assoc}<4$ GeV/c.

$\Psi_{2}$-dependent $C(\Delta\phi)$ for $4<p_{T}^{trig}<10$ GeV/c and $2<p_{T}^{assoc}<4$ GeV/c.

$\Psi_{3}$-dependent $C(\Delta\phi)$ for $2<p_{T}^{trig}<4$ GeV/c and $1<p_{T}^{assoc}<2$ GeV/c.

$\Psi_{3}$-dependent $C(\Delta\phi)$ for $2<p_{T}^{trig}<4$ GeV/c and $2<p_{T}^{assoc}<4$ GeV/c.

$\Psi_{3}$-dependent $C(\Delta\phi)$ for $4<p_{T}^{trig}<10$ GeV/c and $1<p_{T}^{assoc}<2$ GeV/c.

$\Psi_{2}$-dependent $Y(\Delta\phi)$ for $2<p_{T}^{trig}<4$ GeV/c and $1<p_{T}^{assoc}<2$ GeV/c.

$\Psi_{2}$-dependent $C(\Delta\phi)$ for $2<p_{T}^{trig}<4$ GeV/c and $2<p_{T}^{assoc}<4$ GeV/c.

$\Psi_{2}$-dependent $C(\Delta\phi)$ for $4<p_{T}^{trig}<10$ GeV/c and $2<p_{T}^{assoc}<4$ GeV/c.

$\Psi_{3}$-dependent $C(\Delta\phi)$ for $2<p_{T}^{trig}<4$ GeV/c and $1<p_{T}^{assoc}<2$ GeV/c.

$\Psi_{3}$-dependent $C(\Delta\phi)$ for $2<p_{T}^{trig}<4$ GeV/c and $2<p_{T}^{assoc}<4$ GeV/c.

$\Psi_{3}$-dependent $C(\Delta\phi)$ for $4<p_{T}^{trig}<10$ GeV/c and $2<p_{T}^{assoc}<4$ GeV/c.

$\Psi_2$-dependent C for 2 $< p_T^t <$ 4 and 1 $< p_T^a <$ 2 GeV/c.

$\Psi_2$-dependent C for 2 $< p_T^t <$ 4 and 2 $< p_T^a <$ 4 GeV/c.

$\Psi_2$-dependent C for 4 $< p_T^t <$ 10 and 2 $< p_T^a <$ 4 GeV/c.

$\Psi_2$-dependent Y for 2 $< p_T^t <$ 4 and 1 $< p_T^a <$ 2 GeV/c.

$\Psi_2$-dependent Y for 2 $< p_T^t <$ 4 and 2 $< p_T^a <$ 4 GeV/c.

$\Psi_2$-dependent Y for 4 $< p_T^t <$ 10 and 2 $< p_T^a <$ 4 GeV/c.

$\Psi_3$-dependent C for 2 $< p_T^t <$ 4 and 1 $< p_T^a <$ 2 GeV/c.

$\Psi_3$-dependent C for 2 $< p_T^t <$ 4 and 2 $< p_T^a <$ 4 GeV/c.

$\Psi_3$-dependent C for 4 $< p_T^t <$ 10 and 2 $< p_T^a <$ 4 GeV/c.

$\Psi_3$-dependent Y for 2 $< p_T^t <$ 4 and 1 $< p_T^a <$ 2 GeV/c.

$\Psi_3$-dependent Y for 2 $< p_T^t <$ 4 and 2 $< p_T^a <$ 4 GeV/c.

$\Psi_3$-dependent Y for 4 $< p_T^t <$ 10 and 2 $< p_T^a <$ 4 GeV/c.

$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

Differential invariant yield of electrons (($N_{e^+}$+$N_{e^-}$)/2) from heavy-flavor decays for different centrality classes.

Differential invariant yield of electrons (($N_{e^+}$+$N_{e^-}$)/2) from heavy-flavor decays for different centrality classes.

Differential invariant yield of electrons (($N_{e^+}$+$N_{e^-}$)/2) from heavy-flavor decays for different centrality classes.

Differential invariant yield of electrons (($N_{e^+}$+$N_{e^-}$)/2) from heavy-flavor decays for different centrality classes.

$R_{AuAu}(p^e_T)$ for different centrality classes.

$R_{AuAu}(p^e_T)$ for different centrality classes.

$R_{AuAu}(p^e_T)$ for different centrality classes.

$p^e_T$-integrated nuclear modification factors $R_{AuAu}$.

Nonphotonic electron $v_2$ as a function of $p_T$ in Au+Au collisions at $\sqrt{s_{NN}}$ = 200 GeV.

Experimental correlation moments $R^2_{x2}(q)$ Data. Systematic errors are less than the statistical errors.

Experimental correlation moments $R^2_{x2}(q)$ Fit. Systematic errors are less than the statistical errors.

Experimental correlation moments $R^2_{y2}(q)$ Data. Systematic errors are less than the statistical errors.

Experimental correlation moments $R^2_{y2}(q)$ Fit. Systematic errors are less than the statistical errors.

Experimental correlation moments $R^6_{x6}(q)$ Data. Systematic errors are less than the statistical errors.

Experimental correlation moments $R^6_{x6}(q)$ Fit. Systematic errors are less than the statistical errors.

Experimental correlation moments $R^6_{y6}(q)$ Data. Systematic errors are less than the statistical errors.

Experimental correlation moments $R^6_{y6}(q)$ Fit. Systematic errors are less than the statistical errors.

Experimental correlation moments $R^4_{x2y2}(q)$ Data. Systematic errors are less than the statistical errors.

Experimental correlation moments $R^4_{x2y2}(q)$ Fit. Systematic errors are less than the statistical errors.

Experimental correlation moments $R^4_{x4}(q)$ Data. Systematic errors are less than the statistical errors.

Experimental correlation moments $R^4_{x}(q)$ Fit. Systematic errors are less than the statistical errors.

Experimental correlation moments $R^4_{y4}(q)$ Data. Systematic errors are less than the statistical errors.

Experimental correlation moments $R^4_{y4}(q)$ Fit. Systematic errors are less than the statistical errors.

Experimental correlation moments $R^6_{x2y4}(q)$ Data. Systematic errors are less than the statistical errors.

Experimental correlation moments $R^6_{x2y4}(q)$ Fit. Systematic errors are less than the statistical errors.

Experimental correlation moments $R^6_{x4y2}(q)$ Data. Systematic errors are less than the statistical errors.

Experimental correlation moments $R^6_{x4y2}(q)$ Fit. Systematic errors are less than the statistical errors.

Source function profile $S(r_x)$ Image in the PCMS.

Source function profile $S(r_x)$ Fit in the PCMS. The blue line represents the Hump fit.

Source function profile $S(r_y)$ Image in the PCMS.

Source function profile $S(r_y)$ Fit in the PCMS. The blue line represents the Hump fit.

Source function profile $S(r_z)$ Image in the PCMS.

Source function profile $S(r_z)$ Fit in the PCMS. The blue line represents the Hump fit.

Associated correlation profile $C(q_x)$ Image in the PCMS.

Associated correlation profile $C(q_x)$ Fit in the PCMS.

Associated correlation profile $C(q_y)$ Image in the PCMS.

Associated correlation profile $C(q_y)$ Fit in the PCMS.

Associated correlation profile $C(q_z)$ Image in the PCMS.

Associated correlation profile $C(q_z)$ Fit in the PCMS.

Source function comparison between Therminator calculation and image for $S(r_x)$ in PCMS.

Source function comparison between Therminator calculation and image for $S(r_x)$ in PCMS.

Source function comparison between Therminator calculation and image for $S(r_y)$ in PCMS.

Source function comparison between Therminator calculation and image for $S(r_y)$ in PCMS.

Source function comparison between Therminator calculation and image for $S(r_z)$ in PCMS.

Source function comparison between Therminator calculation and image for $S(r_z)$ in PCMS.

$\Delta$$t_{LCM}$ from Therminator events compared with various assumptions for $\Delta\tau$ and resonance emissions.

$\Delta$$t_{LCM}$ from Therminator events compared with various assumptions for $\Delta\tau$ and resonance emissions.

<p>Jet-pair distributions for associated baryons for $1 < p_{T,assoc} < 1.3\ \mathrm{GeV}/c$ and $1.6 < p_{T,assoc} < 2.0\ \mathrm{GeV}/c$. Results are for a hadron trigger $2.5 < p_T < 4.0\ \mathrm{GeV}/c$ and centrality selections of 0-20% and 20-40%.</p> <p><i>Note that only statistical uncertainties are available.</i></p>

<p>Conditional jet yeilds for associated mesons for near-side jets.</p> <p><i>Note that only statistical uncertainties are available.</i></p>

<p>Conditional jet yeilds for associated mesons for away-side jets.</p> <p><i>Note that only statistical uncertainties are available.</i></p>

<p>Conditional jet yeilds for associated baryons for near-side jets.</p> <p><i>Note that only statistical uncertainties are available.</i></p>

<p>Conditional jet yeilds for associated baryons for away-side jets.</p> <p><i>Note that only statistical uncertainties are available.</i></p>

<p>Ratio of jet-associated baryons to jet-associated mesons for near-side jets.</p> <p><i>Note that only statistical uncertainties are available.</i></p>

<p>Ratio of jet-associated baryons to jet-associated mesons for away-side jets.</p> <p><i>Note that only statistical uncertainties are available.</i></p>

<p>Charged hadron spectra for centrality selected d+Au collisions.</p>

<p>Charged hadron spectra for centrality selected d+Au collisions.</p>

<p>Charged hadron spectra for centrality selected N+Au collisions. All centralities.</p>

<p>Charged hadron spectra for centrality selected N+Au collisions. Centrality A (Central).</p>

<p>Charged hadron spectra for centrality selected N+Au collisions. Centrality B.</p>

<p>Charged hadron spectra for centrality selected N+Au collisions. Centrality C.</p>

<p>Charged hadron spectra for centrality selected N+Au collisions. Centrality D (Peripheral).</p>

<p>$R_{dAu}$ as a function of $p_T$.</p>

<p>$R_{dAu}$ as a function of $p_T$.</p>

<p>$R_{dAu}$ as a function of $p_T$.</p>

<p>$R_{dAu}$ as a function of $p_T$.</p>

<p>$R_{NAu}$ as a function of $p_T$. Centrality A (Peripheral).</p>

<p>$R_{NAu}$ as a function of $p_T$. Centrality B.</p>

<p>$R_{dAu}$ as a function of $p_T$. Centrality C.</p>

<p>$R_{dAu}$ as a function of $p_T$. Centrality D (Central).</p>

<p>$R_{dAu}$ values averaged in three momentum ranges as functions of $\nu−1$. Horizontal bars show the uncertainty in the value of $\nu$ for each centrality class.</p>

<p>$R_{NAu}$ values averaged in three momentum ranges as functions of $\nu−1$. Horizontal bars show the uncertainty in the value of $\nu$ for each centrality class.</p>

<p>Charged hadron spectra for centrality selected d+Au collisions.</p>

The double transverse spin asymmetry ($A_{TT}$) for neutral pion production at $\sqrt{s}$ = 200 GeV as a function of $p_T$ (GeV/$c$). The correlated for all points scale systematic uncertainty is 9.4% (scales both the values and stat. uncertainties by the same factor).

Invariant cross section of $\omega$ production in $p$+$p$ and $d$+Au collisions at $\sqrt{s_{NN}}$ = 200 GeV measured in $\omega \rightarrow \pi^0\pi^+\pi^-$ and $\omega \rightarrow \pi^0\gamma$ decay channels.

Invariant cross section of $\omega$ production in $p$+$p$ and $d$+Au collisions at $\sqrt{s_{NN}}$ = 200 GeV measured in $\omega \rightarrow \pi^0\pi^+\pi^-$ and $\omega \rightarrow \pi^0\gamma$ decay channels.

Invariant cross section of $\omega$ production in $p$+$p$ and $d$+Au collisions at $\sqrt{s_{NN}}$ = 200 GeV measured in $\omega \rightarrow \pi^0\pi^+\pi^-$ and $\omega \rightarrow \pi^0\gamma$ decay channels.

Measured $\omega/{\pi^0}$ vs $p_T$ in $p$+$p$ collisions at $\sqrt{s_{NN}}$ = 200 GeV.

Measured $\omega/{\pi^0}$ vs $p_T$ in $p$+$p$ collisions at $\sqrt{s_{NN}}$ = 200 GeV.

Measured $\omega/{\pi^0}$ vs $p_T$ in $d$+Au collisions at $\sqrt{s_{NN}}$ = 200 GeV.

Measured $\omega/{\pi^0}$ vs $p_T$ in $d$+Au collisions at $\sqrt{s_{NN}}$ = 200 GeV.

Measured $R_{dA}$ vs $p_T$ for neutral mesons in $d$+Au collisions at $\sqrt{s_{NN}}$ = 200 GeV for minimum bias.

Measured $R_{dA}$ vs $p_T$ for neutral mesons in $d$+Au collisions at $\sqrt{s_{NN}}$ = 200 GeV for minimum bias.

Measured $R_{dA}$ vs $p_T$ for neutral mesons in $d$+Au collisions at $\sqrt{s_{NN}}$ = 200 GeV for 0-20% centrality.

Measured $R_{dA}$ vs $p_T$ for neutral mesons in $d$+Au collisions at $\sqrt{s_{NN}}$ = 200 GeV for 0-20% centrality.

Reconstructed $\omega$ mass measured in the $\pi^0\pi^+\pi^-$ decay channel vs $p_T$ in $p$+$p$ collisions at $\sqrt{s_{NN}}$ = 200 GeV.

Reconstructed $\omega$ mass measured in the $\pi^0\pi^+\pi^-$ decay channel vs $p_T$ in $d$+Au collisions at $\sqrt{s_{NN}}$ = 200 GeV.

Neutral pion invariant yields as a function of $p_T$ measured in minimum bias and 9 centrality classes in Au+Au collisions at $\sqrt{s_{NN}}$ = 200 GeV.

Neutral pion invariant yields as a function of $p_T$ measured in minimum bias and 9 centrality classes in Au+Au collisions at $\sqrt{s_{NN}}$ = 200 GeV.

$R_{AA}$ vs. $\Delta\phi$ for $\pi^0$ yields integrated over 3 < $p_T$ < 5 GeV/$c$.

$<R_{AA}>$ for $\pi^0$ yields integrated over 3 < $p_T$ < 5 GeV/$c$.

$R_{AA}$ vs. $\Delta\phi$ for $\pi^0$ yields integrated over 5 < $p_T$ < 8 GeV/$c$.

$<R_{AA}>$ for $\pi^0$ yields integrated over 5 < $p_T$ < 8 GeV/$c$.

Cross sections at mid-rapidity in $p$+$p$ and $d$+Au collisions as a function of $p_T$.

Cross sections at mid-rapidity in $p$+$p$ and $d$+Au collisions as a function of $p_T$.

Cross sections at mid-rapidity in $p$+$p$ and $d$+Au collisions as a function of $p_T$.

Cross sections at mid-rapidity in $p$+$p$ and $d$+Au collisions as a function of $p_T$.

Nuclear modification factor $R_{dA}$ for $\pi^0$ and $\eta$ in differenet centrality selections and min. bias data.

Nuclear modification factor $R_{dA}$ for $\pi^0$ and $\eta$ in differenet centrality selections and min. bias data.

Nuclear modification factor $R_{dA}$ for $\pi^0$ and $\eta$ in differenet centrality selections and min. bias data.

Correlation function, C(q) for $\pi^+\pi^+$ and $\pi^-\pi^-$ pairs.

Correlation function, C(q) for $\pi^+\pi^+$ and $\pi^-\pi^-$ pairs.

Correlation function, C(q) for $\pi^+\pi^+$ and $\pi^-\pi^-$ pairs.

Correlation functions, C(q), for $\pi^+\pi^+$ and $\pi^-\pi^-$ pairs.

Correlation functions, C(q), for $\pi^+\pi^+$ and $\pi^-\pi^-$ pairs.

Correlation functions, C(q), for $\pi^+\pi^+$ and $\pi^-\pi^-$ pairs.

Correlation functions, C(q), for $\pi^+\pi^+$ and $\pi^-\pi^-$ pairs.

Correlation functions, C(q), for $\pi^+\pi^+$ and $\pi^-\pi^-$ pairs.

Correlation functions, C(q), for $\pi^+\pi^+$ and $\pi^-\pi^-$ pairs.

Correlation functions, C(q), for $\pi^+\pi^+$ and $\pi^-\pi^-$ pairs.

Correlation functions, C(q), for $\pi^+\pi^+$ and $\pi^-\pi^-$ pairs.

$k_T$ dependence of extracted spheroidal source parameters for pion source functions (a) $R_T$ , (b) $a$, and (c) $\lambda$ for most central collisions & perhipheral collisions.

$k_T$ dependence of extracted spheroidal source parameters for pion source functions (a) $R_T$ , (b) $a$, and (c) $\lambda$ for most central collisions & perhipheral collisions.

$k_T$ dependence of extracted spheroidal source parameters for pion source functions (a) $R_T$ , (b) $a$, and (c) $\lambda$ for most central collisions & perhipheral collisions.

Transverse momentum in GeV/$c$ for $K^{\pm}$.

Transverse momentum in GeV/$c$ for $K^{\pm}$.

Transverse momentum in GeV/$c$ for $p$ and $\bar{p}$.

$dN/dy$ and $T$ for different centrality bins.

Blast wave model parameters as a function of centrality from fitting $\pi^±$, $K^±$, $p$, and $\bar{p}$ spectra.

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

Centrality dependence of the $\phi$ mass centroid and $\phi$ intrinsic width.

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

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

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

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.

Inverse slope parameteres $T$ in $m_T$ spectra.

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

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.

$\Delta\phi$ distributions for meson triggers with 2.5 < $p_T$ < 4.0 GeV/$c$ and associated charged hadrons with 1.7 < $p_T$ < 2.5 GeV/$c$ for 60-70% centrality in Au+Au collisions.

$\Delta\phi$ distributions for baryon triggers with 2.5 < $p_T$ < 4.0 GeV/$c$ and associated charged hadrons with 1.7 < $p_T$ < 2.5 GeV/$c$ for 60-70% centrality in Au+Au collisions.

Combinatorial Background Levels from Figure 1 Top Panels.

Combinatorial Background Levels from Figure 1 Bottom Panels.

Trigger Particle $v_2$ values.

Partner Particle $v_2$ values.

Yield per trigger for associated charged hadrons between 1.7 < $p_T$ < 2.5 GeV/$c$ for the near- (top) and away- (bottom) side jets.

Yield per trigger for associated charged hadrons between 1.7 < $p_T$ < 2.5 GeV/$c$ for the near- (top) and away- (bottom) side jets.

Uncertainties in Figure 2 as a function of $N_{part}$ and trigger type for AuAu and $d$+Au.

Uncertainties in Figure 2 as a function of $N_{part}$ and trigger type for AuAu and $d$+Au.

$p_T$ spectra of the near side associated charged hadrons corrected to the full jet yield for meson triggers at 2.5 < $p_T$ < 4.0 GeV/$c$ and $|\eta|$ < 0.35 for six centralities in Au+Au and $d$+Au collisions.

$p_T$ spectra of the near side associated charged hadrons corrected to the full jet yield for meson triggers at 2.5 < $p_T$ < 4.0 GeV/$c$ and $|\eta|$ < 0.35 for six centralities in Au+Au and $d$+Au collisions.

$p_T$ spectra of the near side associated charged hadrons corrected to the full jet yield for baryon triggers at 2.5 < $p_T$ < 4.0 GeV/$c$ and $|\eta|$ < 0.35 for six centralities in Au+Au, $d$+Au and $p$+$p$ collisions.

Inverse slopes from the fits in Figure 3.

Inverse slopes from the fits in Figure 3.

Inverse slopes from the fits in Figure 3.

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.

$F_{p_T}$ (in percent) of non-random fluctuations as a function of the $p_T$ range over which $M_{p_T}$ is calculated, 0.2 GeV/$c$ < $p_T$ < $p^{max}_T$, for the 20-25% centrality class ($N_{part}$ = 181.6).

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

Invariant yields from the $p$+$p$ reference data. The normalization uncertainty is labeled norm, the uncertainty from the fit is labeled fit, and the uncertainty $\pi^0$ is from the $\pi^0$ ratio.

Nuclear modification factor $R_{AA}$ in 10% most central $Au$+$Au$ collisions.

Nuclear modification factor $R_{dA}$ for ($h^+$+$h^-$)/2 in minimum bias $d$+$Au$.

$R_{AA}$ in the 10% most central $Au$+$Au$ collisions.

$R_{dA}$ for ($h^+$+$h^-$)/2 measurement in $d$+$Au$.

$R_{dA}$ for $\pi^0$ measurement in $d$+$Au$.

Ratio of $R_{dA}$ in minimum bias $d$+$Au$ and $R_{pA}$ from neutron tagged $d$+$Au$ collisions for ($h^+$+$h^-$)/2.

Ratio of $R_{dA}$ in minimum bias $d$+$Au$ and $R_{pA}$ from neutron tagged $d$+$Au$ collisions for $\pi^0$

Invariant $\pi^0$ yields at midrapidity as a function of $p_T$ for minimum bias and nine centralities in $Au\ +\ Au$ at $\sqrt{s_{NN}} = 200\ GeV$ [0%–10% (80%–92%) is most central (peripheral)]. The labels "uncorr." and "corr." include systematic errors that are uncorrelated and correlated point-to-point, respectively.

Invariant $\pi^0$ yields at midrapidity as a function of $p_T$ for minimum bias and nine centralities in $Au\ +\ Au$ at $\sqrt{s_{NN}} = 200\ GeV$ [0%–10% (80%–92%) is most central (peripheral)]. The labels "uncorr." and "corr." include systematic errors that are uncorrelated and correlated point-to-point, respectively.

Invariant $\pi^0$ yields at midrapidity as a function of $p_T$ for minimum bias and nine centralities in $Au\ +\ Au$ at $\sqrt{s_{NN}} = 200\ GeV$ [0%–10% (80%–92%) is most central (peripheral)]. The labels "uncorr." and "corr." include systematic errors that are uncorrelated and correlated point-to-point, respectively.

Invariant $\pi^0$ yields at midrapidity as a function of $p_T$ for minimum bias and nine centralities in $Au\ +\ Au$ at $\sqrt{s_{NN}} = 200\ GeV$ [0%–10% (80%–92%) is most central (peripheral)]. The labels "uncorr." and "corr." include systematic errors that are uncorrelated and correlated point-to-point, respectively.

Invariant $\pi^0$ yields at midrapidity as a function of $p_T$ for minimum bias and nine centralities in $Au\ +\ Au$ at $\sqrt{s_{NN}} = 200\ GeV$ [0%–10% (80%–92%) is most central (peripheral)]. The labels "uncorr." and "corr." include systematic errors that are uncorrelated and correlated point-to-point, respectively.

Invariant $\pi^0$ yields at midrapidity as a function of $p_T$ for minimum bias and nine centralities in $Au\ +\ Au$ at $\sqrt{s_{NN}} = 200\ GeV$ [0%–10% (80%–92%) is most central (peripheral)]. The labels "uncorr." and "corr." include systematic errors that are uncorrelated and correlated point-to-point, respectively.

Invariant $\pi^0$ yields at midrapidity as a function of $p_T$ for minimum bias and nine centralities in $Au\ +\ Au$ at $\sqrt{s_{NN}} = 200\ GeV$ [0%–10% (80%–92%) is most central (peripheral)]. The labels "uncorr." and "corr." include systematic errors that are uncorrelated and correlated point-to-point, respectively.

Invariant $\pi^0$ yields at midrapidity as a function of $p_T$ for minimum bias and nine centralities in $Au\ +\ Au$ at $\sqrt{s_{NN}} = 200\ GeV$ [0%–10% (80%–92%) is most central (peripheral)]. The labels "uncorr." and "corr." include systematic errors that are uncorrelated and correlated point-to-point, respectively.

Invariant $\pi^0$ yields at midrapidity as a function of $p_T$ for minimum bias and nine centralities in $Au\ +\ Au$ at $\sqrt{s_{NN}} = 200\ GeV$ [0%–10% (80%–92%) is most central (peripheral)]. The labels "uncorr." and "corr." include systematic errors that are uncorrelated and correlated point-to-point, respectively.

Invariant $\pi^0$ yields at midrapidity as a function of $p_T$ for minimum bias and nine centralities in $Au\ +\ Au$ at $\sqrt{s_{NN}} = 200\ GeV$ [0%–10% (80%–92%) is most central (peripheral)]. The labels "uncorr." and "corr." include systematic errors that are uncorrelated and correlated point-to-point, respectively.

Nuclear modification factor $R_{AA}(p_T)$ for $\pi^0$ in Au+Au collisions at $\sqrt{s_{NN}} = 200$ GeV. The error bars include all point-to-point experimental ($p+p$, Au+Au) errors. The label "sys(uncorr)" refers to the uncorrelated point-to-point uncertainty (which includes statistical errors), while "sys(norm)" is the normalization uncertainty.

Nuclear modification factor $R_{AA}(p_T)$ for $\pi^0$ in Au+Au collisions at $\sqrt{s_{NN}} = 200$ GeV. The error bars include all point-to-point experimental ($p+p$, Au+Au) errors. The label "sys(uncorr)" refers to the uncorrelated point-to-point uncertainty (which includes statistical errors), while "sys(norm)" is the normalization uncertainty.

Nuclear modification factor $R_{AA}(p_T)$ for $\pi^0$ in Au+Au collisions at $\sqrt{s_{NN}} = 200$ GeV. The error bars include all point-to-point experimental ($p+p$, Au+Au) errors. The label "sys(uncorr)" refers to the uncorrelated point-to-point uncertainty (which includes statistical errors), while "sys(norm)" is the normalization uncertainty.

Nuclear modification factor $R_{AA}(p_T)$ for $\pi^0$ in Au+Au collisions at $\sqrt{s_{NN}} = 200$ GeV. The error bars include all point-to-point experimental ($p+p$, Au+Au) errors. The label "sys(uncorr)" refers to the uncorrelated point-to-point uncertainty (which includes statistical errors), while "sys(norm)" is the normalization uncertainty.

Nuclear modification factor $R_{AA}(p_T)$ for $\pi^0$ in Au+Au collisions at $\sqrt{s_{NN}} = 200$ GeV. The error bars include all point-to-point experimental ($p+p$, Au+Au) errors. The label "sys(uncorr)" refers to the uncorrelated point-to-point uncertainty (which includes statistical errors), while "sys(norm)" is the normalization uncertainty.

Nuclear modification factor $R_{AA}(p_T)$ for $\pi^0$ in Au+Au collisions at $\sqrt{s_{NN}} = 200$ GeV. The error bars include all point-to-point experimental ($p+p$, Au+Au) errors. The label "sys(uncorr)" refers to the uncorrelated point-to-point uncertainty (which includes statistical errors), while "sys(norm)" is the normalization uncertainty.

Nuclear modification factor $R_{AA}(p_T)$ for $\pi^0$ in Au+Au collisions at $\sqrt{s_{NN}} = 200$ GeV. The error bars include all point-to-point experimental ($p+p$, Au+Au) errors. The label "sys(uncorr)" refers to the uncorrelated point-to-point uncertainty (which includes statistical errors), while "sys(norm)" is the normalization uncertainty.

Nuclear modification factor $R_{AA}(p_T)$ for $\pi^0$ in Au+Au collisions at $\sqrt{s_{NN}} = 200$ GeV. The error bars include all point-to-point experimental ($p+p$, Au+Au) errors. The label "sys(uncorr)" refers to the uncorrelated point-to-point uncertainty (which includes statistical errors), while "sys(norm)" is the normalization uncertainty.

Nuclear modification factor $R_{AA}(p_T)$ for $\pi^0$ in Au+Au collisions at $\sqrt{s_{NN}} = 200$ GeV. The error bars include all point-to-point experimental ($p+p$, Au+Au) errors. The label "sys(uncorr)" refers to the uncorrelated point-to-point uncertainty (which includes statistical errors), while "sys(norm)" is the normalization uncertainty.

Nuclear modification factor $R_{AA}(p_T)$ for $\pi^0$ in Au+Au collisions at $\sqrt{s_{NN}} = 200$ GeV. The error bars include all point-to-point experimental ($p+p$, Au+Au) errors. The label "sys(uncorr)" refers to the uncorrelated point-to-point uncertainty (which includes statistical errors), while "sys(norm)" is the normalization uncertainty.

Nuclear modification factor $R_{AA}(p_T)$ for $\pi^0$ in Au+Au collisions at $\sqrt{s_{NN}} = 200$ GeV. The error bars include all point-to-point experimental ($p+p$, Au+Au) errors. The label "sys(uncorr)" refers to the uncorrelated point-to-point uncertainty (which includes statistical errors), while "sys(norm)" is the normalization uncertainty.

Nuclear modification factor $R_{AA}(p_T)$ for $\pi^0$ in Au+Au collisions at $\sqrt{s_{NN}} = 200$ GeV. The error bars include all point-to-point experimental ($p+p$, Au+Au) errors. The label "sys(uncorr)" refers to the uncorrelated point-to-point uncertainty (which includes statistical errors), while "sys(norm)" is the normalization uncertainty.

Nuclear modification factor $R_{AA}(p_T)$ for $\pi^0$ in Au+Au collisions at $\sqrt{s_{NN}} = 200$ GeV. The error bars include all point-to-point experimental ($p+p$, Au+Au) errors. The label "sys(uncorr)" refers to the uncorrelated point-to-point uncertainty (which includes statistical errors), while "sys(norm)" is the normalization uncertainty.

The $M_{e_T}$ distributions for four different centrality classes. The curves are the random baseline mixed event distributions. The source of differences in the data and mixed event distributions are addressed in the text.

The residual distribution between the data and mixed event $M_{e_T}$ spectra as a function of $M_{e_T}$ in units of standard deviation. The total ${\chi}^2$ and the number of degrees of freedom for the $0-5\%$, $0-10\%$, $10-20\%$, $20-30\%$ centrality classes are 310.0/32, 896.4/36, 678.7/47, and 553.9/53, respectively. A large fraction of the residual contributions are due to the effects of cluster merging.

$v(Q)$, for central events, as a function of the azimuthal coverage of the detector for data. (Angles in degrees.)

The effect on $v(Q)$ of varying the acceptance in $\varphi_r$. The solid curve shows the expectation from global charge conservation. The 10% most central events are used for data and RQMD. (Angles in degrees.)

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.