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

Signal-to-background (S/B) ratio as a function of transverse momentum, Au+Au 62.4 GeV, 0-60% central events with minimum bias trigger

Inclusive electron elliptic flow v2\{2\} as a function of transverse momentum, Au+Au 200 GeV, 0-60% central events with minimum bias trigger

Inclusive electron elliptic flow v2\{4\} as a function of transverse momentum, Au+Au 200 GeV, 0-60% central events with minimum bias trigger

Inclusive electron elliptic flow v2\{2\} as a function of transverse momentum, Au+Au 39 GeV, 0-60% central events with minimum bias trigger

Inclusive electron elliptic flow v2\{2\} as a function of transverse momentum, Au+Au 62.4 GeV, 0-60% central events with minimum bias trigger

Inclusive electron elliptic flow v2\{2\} as a function of transverse momentum, Au+Au 200 GeV, 0-60% central events with High Tower (high pT) trigger

Photonic electron elliptic flow v2 as a function of transverse momentum, Au+Au 39 GeV, 0-60% central events with minimum bias trigger

Photonic electron elliptic flow v2 as a function of transverse momentum, Au+Au 62.4 GeV, 0-60% central events with minimum bias trigger

Photonic electron elliptic flow v2 as a function of transverse momentum, Au+Au 200 GeV, 0-60% central events with minimum bias trigger

Elliptic flow v2\{2\} of electrons from heavy-flavor hadron decays, Au+Au 200 GeV, 0-60% central events with minimum bias trigger

Elliptic flow v2\{4\} of electrons from heavy-flavor hadron decays, Au+Au 200 GeV, 0-60% central events with minimum bias trigger

Elliptic flow v2\{2\} of electrons from heavy-flavor hadron decays, Au+Au 200 GeV, 0-60% central events with High Tower (high pT) trigger

Elliptic flow v2\{EP\} of electrons from heavy-flavor hadron decays, Au+Au 200 GeV, 0-60% central events with High Tower (high pT) trigger

Elliptic flow v2\{2\} of electrons from heavy-flavor hadron decays, Au+Au 39 GeV, 0-60% central events with minimum bias trigger

Elliptic flow v2\{2}\ of electrons from heavy-flavor hadron decays, Au+Au 62.4 GeV, 0-60% central events with minimum bias trigger

Nuclear modification factors Rcp for $\pi^{+}$ + $\pi^{-}$ and p + $\bar{p}$ in 200 GeV AuAu collisions.

Nuclear modification factors Rcp for $\pi^{+}$ + $\pi^{-}$ and p + $\bar{p}$ in 200 GeV AuAu collisions.

Nuclear modification factors Rcp for $\pi^{+}$ + $\pi^{-}$ and p + $\bar{p}$ in 200 GeV AuAu collisions.

The $\pi^{-}$/$\pi{+}$ and $\bar{p}$/p ratios in 12 $\%$ central, MB Au+Au and d+Au collisions at $\sqrt{s}_{NN}$ = 200 GeV. Errors are statistical and point-to-point systematic.

The p/$\pi^{+}$ and $\bar{p}$/$\pi^{-}$ ratios from d+Au and Au+Au collisions at $\sqrt{s}_{NN}$ = 200 GeV. The systematic uncertainties for 60-80$\%$ Au+Au collisions are similar.

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

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

Measured $\pi^0$ $R_{AA}$ as a function of $p_T$ for the 0-10% most central collisions at $\sqrt{s_{NN}}$ = 22.4 GeV in comparison to a jet quenching calculation. The error (tot.) includes the quadratic sum of the statistical uncertainties and the point-to-point uncorrelated and correlated systematic uncertainties.

Measured $\pi^0$ $R_{AA}$ as a function of $p_T$ for the 0-10% most central collisions at $\sqrt{s_{NN}}$ = 62.4 GeV in comparison to a jet quenching calculation. The error (tot.) includes the quadratic sum of the statistical uncertainties and the point-to-point uncorrelated and correlated systematic uncertainties.

Measured $\pi^0$ $R_{AA}$ as a function of $p_T$ for the 0-10% most central collisions at $\sqrt{s_{NN}}$ = 200.0 GeV in comparison to a jet quenching calculation. The error (tot.) includes the quadratic sum of the statistical uncertainties and the point-to-point uncorrelated and correlated systematic uncertainties.

The average $R_{AA}$ in the interval 2.5 < $p_T$ < 3.5 GeV/$c$ as a function of centrality for Cu+Cu collisions at $\sqrt{s_{NN}}$ = 22.4 GeV. The error (sys.) includes the normalization and $<N_{coll}>$ uncertainties for a typical $N_{coll}$ uncertainty of 12%.

The average $R_{AA}$ in the interval 2.5 < $p_T$ < 3.5 GeV/$c$ as a function of centrality for Cu+Cu collisions at $\sqrt{s_{NN}}$ = 62.4 GeV. The error (sys.) includes the normalization and $<N_{coll}>$ uncertainties for a typical $N_{coll}$ uncertainty of 12%.

The average $R_{AA}$ in the interval 2.5 < $p_T$ < 3.5 GeV/$c$ as a function of centrality for Cu+Cu collisions at $\sqrt{s_{NN}}$ = 200 GeV. The error (sys.) includes the normalization and $<N_{coll}>$ uncertainties for a typical $N_{coll}$ uncertainty of 12%.

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