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

Spectrum in transverse momentum of electrons created in open heavy flavor decays, for 0-10% centrality.

Spectrum in transverse momentum of electrons created in open heavy flavor decays, for 10-20% centrality.

Spectrum in transverse momentum of electrons created in open heavy flavor decays, for 20-40% centrality.

Spectrum in transverse momentum of electrons created in open heavy flavor decays, for 40-60% centrality.

Spectra in transverse momentum of electrons created in open heavy flavor decays, for 60-92% centrality.

Non photonic electrons differential yield scaled by the number of collisions, as a function of centrality. PT belongs to 0.8-4.0 GeV/c.

Non photonic electrons differential yield scaled by the number of collisions, for minimum bias collisions. PT belongs to 0.8-4.0 GeV/c.

Number of collisions as a fonction of centrality There are only systematic errors due to the calculation of Ncoll.

Number of collisions for minimum bias events There are only systematic errors due to the calculation of Ncoll.

Nuclear overlap function TAA as a function of centrality The axis is labelled with 1/SIG, because the unit of the nuclear overlap function is barn-1 There are only systematic errors due to the calculation of TAA.

Nuclear overlap function TAA for minimum bias events The axis is labelled with 1/SIG, as the unit of the nuclear overlap function is barn-1 There are only systematic errors due to the calculation of TAA.

Differential charm cross section scaled by the number of collisions, at mid rapidity, as a function of centrality The differential cross section per number of collisions D(SIG)/Ncoll is calculated as following : D(SIG)/Ncoll = D(N(cc_bar))/TAA.

Differential charm cross section scaled by the number of collisions, at mid rapidity, for minimum bias events The differential cross section per number of collisions D(SIG)/Ncoll is calculated as following : D(SIG)/Ncoll = D(N(cc_bar))/TAA.

Total charm cross section scaled by the number of collisions, as a function of centrality The total cross section per number of collisions SIG/Ncoll is calculated as following : SIG/Ncoll = N(cc_bar)/TAA.

Total charm cross section scaled by the number of collisions, for minimum bias events The total cross section per number of collisions SIG/Ncoll is calculated as following : SIG/Ncoll = N(cc_bar)/TAA.

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.

Mean PT value for J/PSI production. The mean PT is obtained with a phenomonological fit of the J/PSI distribution in PT of the form (1/(2*PI*PT))*D(SIG)/DPT = A ( 1+(PT/B)^2)^-6 .The value given is a combined value of the mean PT for the electron channel results and the muon channel results. The systematic uncertainties were estimated from the spread in mean pt from a weighted mean of the binned data, the phenomenological fit and an exponential fit.An additionnal 3% incertainty is added in the systenatic error for the muon PT due to the incertainty in momentum scale.

Value of the total cross-section for J/PSI production. The total cross section is estimated using a pythia calculation, which reproduces the sigma distribution in rapidity the best, normalized to our data.The systematic error given is obtained by setting the measured cross sections to their upper and lower systematic error limits and noting how the total cross section changed.The choice of the PDF and model was estimated to have little impact in the value of the total cross section (~3%).In addition to the systematic error given, there is a +-40 microbarn error due to normalization.

Centrality dependence of the $p_{T}$ distribution for $K^{-}$ in Au+Au collisions at $\sqrt{s_{NN}}$ = 200 GeV. Errors are statistical only.

Centrality dependence of the $p_{T}$ distribution for protons in Au+Au collisions at $\sqrt{s_{NN}}$ = 200 GeV. Errors are statistical only.

Centrality dependence of the $p_{T}$ distribution for antiprotons in Au+Au collisions at $\sqrt{s_{NN}}$ = 200 GeV. Errors are statistical only.

Centrality dependence of the $p_{T}$ distribution for antiprotons in Au+Au collisions at $\sqrt{s_{NN}}$ = 200 GeV. Errors are statistical only.

Inverse slope parameters for $\pi^{+}$, $\pi^{-}$, $K^{+}$, $K^{-}$, proton, and antiproton for the 0-5%, 40-50%, and 60-92% centrality bins. Errors are statistical only.

Fit paramaters $T_{0}^{+}$ and $T_{0}^{-}$. Errors are statistical only.

Fit paramaters $u_{t}^{+}$ and $u_{t}^{-}$. Errors are statistical only.

Relationship between $N_{part}$ and centrality.

Mean transverse momentum as a function of centrality for pions, kaons, protons,and antiprotons in Au+Au collisions at $\sqrt{s_{NN}}$ = 200 GeV. Errors are systematic only.

dN/dy as a function of centrality for pions, kaons, protons,and antiprotons in Au+Au collisions at $\sqrt{s_{NN}}$ = 200 GeV. Errors are systematic only.

Particle ratios of $\pi^{-}$/$\pi^{+}$ for centrality 0-5% and 60-92%. Errors are statistical and systematic.

Particle ratios of $K^{-}$/$K^{+}$ for centrality 0-5% and 60-92%. Errors are statistical and systematic.

Particle ratios of antiprotons/protons for centrality 0-5%, 60-92%, and minimum bias. Errors are statistical and systematic.

Particle ratios of kaons/pions for centrality 0-5% and 60-92%. Errors are statistical and systematic.

Particle ratios of proton/$\pi^{+}$ and antiproton/$\pi^{-}$ for centrality 0-10%, 20-30%, and 60-92%. Errors are statistical and systematic.

Particle ratios of proton/$\pi^{0}$ and antiproton/$\pi^{0}$ for centrality 0-10%, 20-30%, and 60-92%. Errors are statistical and systematic.

Centrality dependence of particle ratios for $\pi^{-}$/$\pi^{+}$, $K^{-}$/$K^{+}$, and antiprotons/protons. Errors are statistical and systematic.

Rcp as a function of centrality 0-10% and 60-92% for ($\pi^{+}$ + $\pi^{-}$)/2. Errors are statistical only.

Rcp as a function of centrality 0-10% and 60-92% for ($K^{+}$ + $K^{-}$)/2. Errors are statistical only.

Rcp as a function of centrality 0-10% and 60-92% for charged (p+pbar)/2. Errors are statistical only.

Rcp as a function of centrality 0-10% and 60-92% for neutral pions. Errors are statistical only.

Integrated Rcp above $p_{T}$ = 1.5 GeV/c as a function of centrality 60-92% for pions, kaons, and (p+pbar)/2. Errors are statistical only.

The minimum $\chi^2$ and the parameters $T_{fo}$ and $\beta^{max}_T$ for each of the five centrality selections. The best fit parameters are determined by averaging all parameter pairs within the 1$\sigma$ contour. The errors correspond to the standard deviation of the parameter pairs within the 1$\sigma$ $\chi^2$ contour.

Fit parameters for each particle species using Equations $\frac{dN}{dy}$ = $C \cdot$ ($N_{part}$)$^{\alpha_{part}}$ and $\frac{dN}{dy}$ = $C^{\prime} \cdot$ ($N_{coll}$)$^{\alpha_{coll}}$.

Values of the parameters $n_{pp}$ and $x$ from fitting Equation $R \equiv \frac{dN/dy}{N_{part}}$ = (1-$x$) $\cdot n_{pp} \frac{1}{2}$ + $x \cdot n_{pp} \frac{N_{coll}}{N_{part}}$ = $n_{pp}$ [$\frac{1}{2}$ + $x$($\frac{N_{coll}}{N_{part}}$ - $\frac{1}{2})$]to the observed $dN$/$dy$ per $N_{part}$.

Invariant yields for $\pi^{\pm}$ measured in minimum-bias events at midrapidity and normalized to one rapidity unit.

Invariant yields for $K^{\pm}$ measured in minimum-bias events at midrapidity and normalized to one rapidity unit.

Invariant yields for (anti)$p$ measured in minimum-bias events at midrapidity and normalized to one rapidity unit.

Pion invariant yields in each event centrality and $p_T$ bin measured at midrapidity, normalized to one rapidity unit.

Pion invariant yields in each event centrality and $p_T$ bin measured at midrapidity, normalized to one rapidity unit.

Kaon invariant yields in each event centrality and $p_T$ bin measured at midrapidity, normalized to one rapidity unit.

Kaon invariant yields in each event centrality and $p_T$ bin measured at midrapidity, normalized to one rapidity unit.

(Anti)proton invariant yields in each event centrality and $p_T$ bin measured at midrapidity, normalized to one rapidity unit.

(Anti)proton invariant yields in each event centrality and $p_T$ bin measured at midrapidity, normalized to one rapidity unit.

(Anti)proton invariant yields in each event centrality and $p_T$ bin measured at midrapidity, normalized to one rapidity unit.

(Anti)proton invariant yields in each event centrality and $p_T$ bin measured at midrapidity, normalized to one rapidity unit.

(Anti)proton invariant yields in each event centrality and $p_T$ bin measured at midrapidity, normalized to one rapidity unit.

(Anti)proton invariant yields in each event centrality and $p_T$ bin measured at midrapidity, normalized to one rapidity unit.

(Anti)proton invariant yields in each event centrality and $p_T$ bin measured at midrapidity, normalized to one rapidity unit.

Minimum-bias invariant yields for $\pi^{\pm}$ in equal $p_T$ bins.

Minimum-bias invariant yields for $K^{\pm}$ in equal $p_T$ bins.

Minimum-bias invariant yields for (anti)$p$ in equal $p_T$ bins.

Pion invariant yields in each event centrality normalized to one rapidity unit at midrapidity.

Pion invariant yields in each event centrality normalized to one rapidity unit at midrapidity.

Kaon invariant yields in each event centrality normalized to one rapidity unit at midrapidity.

(Anti)proton invariant yields in each event centrality normalized to one rapidity unit at midrapidity.

(Anti)proton invariant yields in each event centrality normalized to one rapidity unit at midrapidity.

The integrated mean $p_T$ for pions, kaons, and (anti)protons produced in the five different classes of event centrality.

The expansion parameters $T_{fo}$ and $\beta^{max}_T$ as a function of the number of participants.

Transverse momentum dependepnce of $v_2$ for identified particle $\pi^-$.

Transverse momentum dependepnce of $v_2$ for identified particle $K^-$.

Transverse momentum dependepnce of $v_2$ for identified particle $\bar{p}$.

Transverse momentum dependepnce of $v_2$ for identified particle $\pi^+$.

Transverse momentum dependepnce of $v_2$ for identified particle $K$.

Transverse momentum dependepnce of $v_2$ for identified particle $p$.

Transverse momentum dependepnce of $v_2$ for combined particles $\pi^+$ and $\pi^-$.

Transverse momentum dependepnce of $v_2$ for combined particles $K^+$ and $K^-$.

Transverse momentum dependepnce of $v_2$ for combined particles $p$ and $\bar{p}$.

Transverse momentum dependepnce of $v_2$ for combined particles $\pi^+$ and $\pi^-$ scaled by quarks for each particle as motivated by a quark coalescence model.

Transverse momentum dependepnce of $v_2$ for combined particles $K^+$ and $K^-$ scaled by quarks for each particle as motivated by a quark coalescence model.

Transverse momentum dependepnce of $v_2$ for combined particles $p$ and $\bar{p}$ scaled by quarks for each particle as motivated by a quark coalescence model.

Transverse momentum dependence of $v_2$ for negative charged particle distributions at 0-20% centrality.

Transverse momentum dependence of $v_2$ for combined $\pi^-$ and $K^-$ at 0-20% centrality.

Transverse momentum dependence of $v_2$ for $\bar{p}$ at 0-20% centrality.

Transverse momentum dependence of $v_2$ for positive charged particle distributions at 0-20% centrality.

Transverse momentum dependence of $v_2$ for combined $\pi^+$ and $K^+$ at 0-20% centrality.

Transverse momentum dependence of $v_2$ for $p$ at 0-20% centrality.

Transverse momentum dependence of $v_2$ for negative charged particle distributions at 20-40% centrality.

Transverse momentum dependence of $v_2$ for combined $\pi^-$ and $K^-$ at 20-40% centrality.

Transverse momentum dependence of $v_2$ for $\bar{p}$ at 20-40% centrality.

Transverse momentum dependence of $v_2$ for positive charged particle distributions at 20-40% centrality.

Transverse momentum dependence of $v_2$ for combined $\pi^+$ and $K^+$ at 20-40% centrality.

Transverse momentum dependence of $v_2$ for $p$ at 20-40% centrality.

Transverse momentum dependence of $v_2$ for negative charged particle distributions at 40-60% centrality.

Transverse momentum dependence of $v_2$ for combined $\pi^-$ and $K^-$ at 40-60% centrality.

Transverse momentum dependence of $v_2$ for $\bar{p}$ at 40-60% centrality.

Transverse momentum dependence of $v_2$ for positive charged particle distributions at 40-60% centrality.

Transverse momentum dependence of $v_2$ for combined $\pi^+$ and $K^+$ at 40-60% centrality.

Transverse momentum dependence of $v_2$ for $p$ at 40-60% centrality.

Nucelar modification factor $R_{CP}$ for ($p+\bar{p}$)/2.

Nucelar modification factor $R_{CP}$ for $\pi^0$.

Charged hadron to $\pi^0$ ratio in central (0-10%) Au+Au collisions.

Charged hadron to $\pi^0$ ratio in peripheral (60-92%) Au+Au collisions. The peripheral data points are offset by +130 MeV/$c$ for clarity.

Systematic error on normalization a) common to both centralities b) that can move the $h$/$\pi^0$ ratios for each centrality independent of the other.

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.

Centrality dependence of $⟨p^{trunc}_T⟩$, the average $p_T$ of charged particles with $p_T$ above a threshold $p^{min}_T$ minus the threshold $p^{min}_T$. Shown are values for two $p^{min}_T$ cuts, one at $p_T >$ 0.5 GeV/$c$ representing all data presented in Fig. 2 and the other one at $p_T >$ 2 GeV/$c$.

Nuclear modification factor ($R_{AA}$) for the 60-80%, 30-60%, 15-30%, and 5-15% centrality selections compared to the one for the most central sample (0-5%). Due to insufficient statistics $R_{AA}$ is not shown for the 80-92% sample.

The top panel gives the nuclear modification factor $R_{AA}$ for three exclusive $p_T$ regions as a function of the centrality of the collision. The lower panel shows essentially the same quantity but normalized to the number of participant pairs rather than to the number of binary collisions. The dotted line indicates the expectation for scaling with the number of binary collisions (top) or with the number of participants (bottom). Only statistical errors are shown. The systematic error on the scale and the centrality dependence are identical to the errors shown in Fig. 5. These errors are correlated, i.e. take their maximum or minimum value simultaneously for all centrality and $p_T$ selections. In addition, there are also $p_T$ dependent systematic errors, which are given in Table 1. The systematic errors do not alter the trends in the data.

Spectra of inclusive $\Lambda$ ($\bar{\Lambda}$) for minimum-bias events.

Spectra of inclusive feed-down corrected protons (antiprotons) for minimum-bias events.

Inclusive $\Lambda$ and feed-down corrected proton yields for minimum bias and for the 5% most central events. Systematic errors are 17% for protons and 16% for $\Lambda$

Inclusive $\Lambda$ and feed-down corrected average transverse momentum for minimum bias and for the 5% most central events. Systematic errors are 17% for protons and 16% for $\Lambda$

Inclusive $\Lambda$ and feed-down corrected slope parameters for minimum bias and for the 5% most central events. Systematic errors are 17% for protons and 16% for $\Lambda$

Total measured and predicted baryon yields for minimum-bias and for the 5% most central events. The proton yields are corrected for feed-down. The systematic errors are 24% for the protons and 23% for the $\Lambda$

$v_2$ vs Fixed $p_T$ for several centrality selections. [F] and [A] follow the notation Fig. 2. Systematic errors are estimated to be $\sim 5$%.

$v_2$ vs Assorted $p_T$ for several centrality selections. [F] and [A] follow the notation Fig. 2. Systematic errors are estimated to be $\sim 5$%.

$v_2$ vs Assorted $p_T$ for several centrality selections. [F] and [A] follow the notation Fig. 2. Systematic errors are estimated to be $\sim 5$%.

The centrality dependence of $A_2$ for two different $p_T$ selections. [F] and [A] follow the notation in Fig. 2. Systematic errors are estimated to be $\sim 10$%, dominated by the normalization of the correction function and the model determination of$ \varepsilon$.

The centrality dependence of $A_2$ for two different $p_T$ selections. [F] and [A] follow the notation in Fig. 2. Systematic errors are estimated to be $\sim 10$%, dominated by the normalization of the correction function and the model determination of$ \varepsilon$.

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

Ratio of the electron data to the calculated background as a function of $p_T$ in minimum bias collisions (upper panel) and relative contributions to the background from various sources (lower panel). The curves for $\omega$ and $\phi$ show the sum of the Dalitz and the di-electron decay modes.

The background-subtracted electron spectra for central (0-10%) collisions compared with the expected contributions from open charm decays. Also shown, for central collisions only, are the expected contribution from bottom decays (dashed) and the conversion electron spectrum from a direct phonon prediction (dotted).

The background-subtracted electron spectra for minimum bias (0-92%) collisions compared with the expected contributions from open charm decays. Also shown, for central collisions only, are the expected contribution from bottom decays (dashed) and the conversion electron spectrum from a direct phonon prediction (dotted).

Single electron cross sections $\frac{d\sigma_e}{dy}|_{y=0}$ of this measurement and ISR experiments [9,11,31] are displayed (bottom of figure, right-hand scale) with charm decay contributions calculated with PYTHIA. Open and filled symbols for 1.0 < $p_T$ < 1.4 GeV/$c$ and $p_T$ > 1.4 GeV/$c$, respectively. The derived charm cross section of this measurement is compared with charm cross sections from SPS/FNAL experiments (top of figure, left-hand scale). The thick curve and the shaded band represent the charm cross section in the PYTHIA model and in a NLO pQCD calculation [30], respectively.

Single electron cross sections $\frac{d\sigma_e}{dy}|_{y=0}$ of this measurement and ISR experiments [9,11,31] are displayed (bottom of figure, right-hand scale) with charm decay contributions calculated with PYTHIA. Open and filled symbols for 1.0 < $p_T$ < 1.4 GeV/$c$ and $p_T$ > 1.4 GeV/$c$, respectively. The derived charm cross section of this measurement is compared with charm cross sections from SPS/FNAL experiments (top of figure, left-hand scale). The thick curve and the shaded band represent the charm cross section in the PYTHIA model and in a NLO pQCD calculation [30], respectively.

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

Transverse momentum spectra for K- in the midrapidity range for the centrality region 0 to 5 PCT. Errors are combined statistical and systematics.

Transverse momentum spectra for P in the midrapidity range for the centrality region 0 to 5 PCT. Errors are combined statistical and systematics.

Transverse momentum spectra for PBAR in the midrapidity range for the centrality region 0 to 5 PCT. Errors are combined statistical and systematics.

Transverse momentum spectra for PI+ in the midrapidity range for the centrality region 15 to 30 PCT. Errors are combined statistical and systematics.

Transverse momentum spectra for PI- in the midrapidity range for the centrality region 15 to 30 PCT. Errors are combined statistical and systematics.

Transverse momentum spectra for K+ in the midrapidity range for the centrality region 15 to 30 PCT. Errors are combined statistical and systematics.

Transverse momentum spectra for K- in the midrapidity range for the centrality region 15 to 30 PCT. Errors are combined statistical and systematics.

Transverse momentum spectra for P in the midrapidity range for the centrality region 15 to 30 PCT. Errors are combined statistical and systematics.

Transverse momentum spectra for PBAR in the midrapidity range for the centrality region 15 to 30 PCT. Errors are combined statistical and systematics.

Transverse momentum spectra for PI+ in the midrapidity range for the centrality region 60 to 92 PCT. Errors are combined statistical and systematics.

Transverse momentum spectra for PI- in the midrapidity range for the centrality region 60 to 92 PCT. Errors are combined statistical and systematics.

Transverse momentum spectra for K+ in the midrapidity range for the centrality region 60 to 92 PCT. Errors are combined statistical and systematics.

Transverse momentum spectra for K- in the midrapidity range for the centrality region 60 to 92 PCT. Errors are combined statistical and systematics.

Transverse momentum spectra for P in the midrapidity range for the centrality region 60 to 92 PCT. Errors are combined statistical and systematics.

Transverse momentum spectra for PBAR in the midrapidity range for the centrality region 60 to 92 PCT. Errors are combined statistical and systematics.

Average transverse momentum for PI+ and PI- as a function of the number of participation collisions (N(C=PART)).

Average transverse momentum for K+ and K- as a function of the number of participation collisions (N(C=PART)).

Average transverse momentum for P and PBAR as a function of the number of participation collisions (N(C=PART)).

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.