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