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

$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

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

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.

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

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

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.

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.