Transverse momentum spectra of antideuterons measured in pp collisions at centre-of-mass per nucleon-nucleon energy of 13 TeV, as shown in Fig. 1 (right panel). Rapidity interval 0.3 to 0.4.

Transverse momentum spectra of antideuterons measured in pp collisions at centre-of-mass per nucleon-nucleon energy of 13 TeV, as shown in Fig. 1 (right panel). Rapidity interval 0.4 to 0.5.

Transverse momentum spectra of antideuterons measured in pp collisions at centre-of-mass per nucleon-nucleon energy of 13 TeV, as shown in Fig. 1 (right panel). Rapidity interval 0.5 to 0.6.

Transverse momentum spectra of antideuterons measured in pp collisions at centre-of-mass per nucleon-nucleon energy of 13 TeV, as shown in Fig. 1 (right panel). Rapidity interval 0.6 to 0.7.

Transverse momentum spectra of antiprotons measured in pp collisions at centre-of-mass per nucleon-nucleon energy of 13 TeV, as shown in Fig. 1 (left panel). Rapidity interval 0 to 0.1.

Transverse momentum spectra of antiprotons measured in pp collisions at centre-of-mass per nucleon-nucleon energy of 13 TeV, as shown in Fig. 1 (left panel). Rapidity interval 0.1 to 0.2.

Transverse momentum spectra of antiprotons measured in pp collisions at centre-of-mass per nucleon-nucleon energy of 13 TeV, as shown in Fig. 1 (left panel). Rapidity interval 0.2 to 0.3.

Transverse momentum spectra of antiprotons measured in pp collisions at centre-of-mass per nucleon-nucleon energy of 13 TeV, as shown in Fig. 1 (left panel). Rapidity interval 0.3 to 0.4.

Transverse momentum spectra of antiprotons measured in pp collisions at centre-of-mass per nucleon-nucleon energy of 13 TeV, as shown in Fig. 1 (left panel). Rapidity interval 0.4 to 0.5.

Transverse momentum spectra of antiprotons measured in pp collisions at centre-of-mass per nucleon-nucleon energy of 13 TeV, as shown in Fig. 1 (left panel). Rapidity interval 0.5 to 0.6.

Transverse momentum spectra of antiprotons measured in pp collisions at centre-of-mass per nucleon-nucleon energy of 13 TeV, as shown in Fig. 1 (left panel). Rapidity interval 0.6 to 0.7.

Integrated yields of antideuterons measured in pp collisions at centre-of-mass per nucleon-nucleon energy of 13 TeV, as a function of absolute value of rapidity. The figure is shown in Fig. 2 (right panel).

Integrated yields of antiprotons measured in pp collisions at centre-of-mass per nucleon-nucleon energy of 13 TeV, as a function of absolute value of rapidity. The figure is shown in Fig. 2 (left panel).

coalescence parameter B2 measured in pp collisions at centre-of-mass per nucleon-nucleon energy of 13 TeV, as a function of transverse momentum per nucleon. The figure is shown in Fig. 3, rapidity interval 0.0 to 0.1, in absolute value.

coalescence parameter B2 measured in pp collisions at centre-of-mass per nucleon-nucleon energy of 13 TeV, as a function of transverse momentum per nucleon. The figure is shown in Fig. 3, rapidity interval 0.1 to 0.2, in absolute value.

coalescence parameter B2 measured in pp collisions at centre-of-mass per nucleon-nucleon energy of 13 TeV, as a function of transverse momentum per nucleon. The figure is shown in Fig. 3, rapidity interval 0.2 to 0.3, in absolute value.

coalescence parameter B2 measured in pp collisions at centre-of-mass per nucleon-nucleon energy of 13 TeV, as a function of transverse momentum per nucleon. The figure is shown in Fig. 3, rapidity interval 0.3 to 0.4, in absolute value.

coalescence parameter B2 measured in pp collisions at centre-of-mass per nucleon-nucleon energy of 13 TeV, as a function of transverse momentum per nucleon. The figure is shown in Fig. 3, rapidity interval 0.4 to 0.5, in absolute value.

coalescence parameter B2 measured in pp collisions at centre-of-mass per nucleon-nucleon energy of 13 TeV, as a function of transverse momentum per nucleon. The figure is shown in Fig. 3, rapidity interval 0.5 to 0.6, in absolute value.

coalescence parameter B2 measured in pp collisions at centre-of-mass per nucleon-nucleon energy of 13 TeV, as a function of transverse momentum per nucleon. The figure is shown in Fig. 3, rapidity interval 0.6 to 0.7, in absolute value.

coalescence parameter B2 measured in pp collisions at centre-of-mass per nucleon-nucleon energy of 13 TeV, as a function of rapidity, for fixed intervals of transverse momentum per nucleon. The data is shown in Fig. 4, in the PT per nucleon interval from 0.9 to 1.0.

coalescence parameter B2 measured in pp collisions at centre-of-mass per nucleon-nucleon energy of 13 TeV, as a function of rapidity, for fixed intervals of transverse momentum per nucleon. The data is shown in Fig. 4, in the PT per nucleon interval from 1.0 to 1.1.

coalescence parameter B2 measured in pp collisions at centre-of-mass per nucleon-nucleon energy of 13 TeV, as a function of rapidity, for fixed intervals of transverse momentum per nucleon. The data is shown in Fig. 4, in the PT per nucleon interval from 1.1 to 1.3.

coalescence parameter B2 measured in pp collisions at centre-of-mass per nucleon-nucleon energy of 13 TeV, as a function of rapidity, for fixed intervals of transverse momentum per nucleon. The data is shown in Fig. 4, in the PT per nucleon interval from 1.3 to 1.5.

coalescence parameter B2 measured in pp collisions at centre-of-mass per nucleon-nucleon energy of 13 TeV, as a function of rapidity, for fixed intervals of transverse momentum per nucleon. The data is shown in Fig. 4, in the PT per nucleon interval from 1.5 to 1.7.

Transverse momentum spectra of deuterons measured in Xe--Xe collisions at centre-of-mass per nucleon-nucleon energy of 5.44 TeV, as shown in Fig. 3 (left panel). Centrality class 20 to 40 percent.

Transverse momentum spectra of deuterons measured in Xe--Xe collisions at centre-of-mass per nucleon-nucleon energy of 5.44 TeV, as shown in Fig. 3 (left panel). Centrality class 40 to 60 percent.

Transverse momentum spectra of deuterons measured in Xe--Xe collisions at centre-of-mass per nucleon-nucleon energy of 5.44 TeV, as shown in Fig. 3 (left panel). Centrality class 60 to 90 percent.

Coalescence parameter B2 as a function of the transverse momentum per nucleon pT/A, in the centrality class 0 to 10 percent. The B2 is shown in Figure 7 (left panel).

Coalescence parameter B2 as a function of the transverse momentum per nucleon pT/A, in the centrality class 10 to 20 percent. The B2 is shown in Figure 7 (left panel).

Coalescence parameter B2 as a function of the transverse momentum per nucleon pT/A, in the centrality class 20 to 40 percent. The B2 is shown in Figure 7 (left panel).

Coalescence parameter B2 as a function of the transverse momentum per nucleon pT/A, in the centrality class 40 to 60 percent. The B2 is shown in Figure 7 (left panel).

Coalescence parameter B2 as a function of the transverse momentum per nucleon pT/A, in the centrality class 60 to 90 percent. The B2 is shown in Figure 7 (left panel).

Coalescence parameter B3 as a function of the transverse momentum per nucleon pT/A, in the centrality class 0 to 90 percent. The B3 is shown in Figure 7 (right panel).

Ratio of integrated yields of deuterons and of protons, as a function of the average charged particle multiplicity density. The d-to-p ratio is shown in Figure 5 (top panel).

Ratio of integrated yields of He3 and of protons, as a function of the average charged particle multiplicity density. The He3-to-p ratio is shown in Figure 5 (bottom panel).

Ratio of integrated yields of deuterons and of pions, as a function of the average charged particle multiplicity density. The d-to-pi ratio is shown in Figure 6 (top panel).

Ratio of integrated yields of He3 and of pions, as a function of the average charged particle multiplicity density. The He3-to-pi ratio is shown in Figure 6 (bottom panel).

B2 as a function of the average charged particle multiplicity density for pT-over-A equal to 0.75 GEV. The B2 vs. multiplicity is not shown in the paper.

(Anti)deuteron $v_2$ measured at midrapidity in the 0-20 centrality class, as shown in Fig. 8 (left panel).

(Anti)deuteron $v_2$ measured at midrapidity in the 20-40 centrality class, as shown in Fig. 8 (right panel).

The double ratio of particle yield of f0((980) to K*(892)0

Transverse-momentum-differential particle yield ratio of f$_{0}$(980) to charged pions and its double ratio

Transverse-momentum-differential particle yield ratio of f$_{0}$(980) to K$^{*}(892)$^{0}$ and its double ratio

Nuclear modification factor of f0(980) in p-Pb collisions at 5.02 TeV

Fully corrected assorted $\pi^{0}$-hadron conditional pair distributions for $d$+Au collisions centrality 0-88% and $p$+$p$ collisions. The trigger $\pi^0$s are within 5 < $p_{T,trig}$ < 10 GeV/$c$ and are correlated with hadrons with $p_{T,assoc}$ of 1.5-2 GeV/$c$, 2-3GeV/$c$, 3-4 GeV/$c$, and 4-5 GeV/$c$.

Near- and far-side widths as a function of $p_{T, assoc}$ for charged hadron azimuthal correlations from minimumbias $d$+Au collisions.

Near-side width, far-side width as a function of $p_{T,assoc}$ for a charged pion and neutral pion triggers from the $p_{T,trig}$ range of 5-10 GeV/$c$ in minimum bias $d$+Au collisions.

Near-side width, far-side width as a function of $p_{T,assoc}$ for a charged pion and neutral pion triggers from the $p_{T,trig}$ range of 5-10 GeV/$c$ in minimum bias $d$+Au collisions.

Near- and far-side conditional yields as a function of $p_{T, assoc}$ for charged hadron triggers (2.5−4 GeV/$c$) and associated charged hadrons from $d$+Au collisions.

Near- and far-side conditional yields as a function of $p_{T, assoc}$ for charged hadron triggers (4−6 GeV/$c$) and associated charged hadrons from $d$+Au collisions.

Near- and far-side conditional yields as a function of $p_{T, assoc}$ for charged pion triggers (5−10 GeV/$c$) and associated charged hadrons from minimum-bias $d$+Au collisions.

Near- and far-side conditional yields as a function of $p_{T, assoc}$ for neutral pion triggers (5−10 GeV/$c$) and associated charged hadrons from minimum-bias $d$+Au collisions.

Near- and far-side widths and conditional yields as a function of $N_{coll}$ for charged hadron triggers (2.5−4 GeV/$c$) and associated charged hadrons (1–2.5 GeV/$c$) from $d$+Au collisions.

Near- and far-side widths and conditional yields as a function of centrality for neutral pion triggers (5−10 GeV/$c$) and associated charged hadrons (2–3 GeV/$c$) from $d$+Au collisions.

The extracted $\sqrt{<j^2_T>}$ for minimum-bias $d$+Au collisions from all four dihadron correlations. The trigger $p_T$ ranges are 3 - 5 GeV/$c$, 5 - 10 GeV/$c$ and 5 - 10 GeV/$c$ for $h^{\pm}$ - $h^{\pm}$, $\pi^0$ - $h^{\pm}$ and $\pi^{\pm}$ - $h^{\pm}$ assorted-$p_T$ correlations.

The extracted $\sqrt{<j^2_T>}$ for minimum-bias $d$+Au collisions from all four dihadron correlations. The trigger $p_T$ ranges are 3 - 5 GeV/$c$, 5 - 10 GeV/$c$ and 5 - 10 GeV/$c$ for $h^{\pm}$ - $h^{\pm}$, $\pi^0$ - $h^{\pm}$ and $\pi^{\pm}$ - $h^{\pm}$ assorted-$p_T$ correlations.

The extracted $\sqrt{<j^2_T>}$ for minimum-bias $d$+Au collisions from all four dihadron correlations. The trigger $p_T$ ranges are 3 - 5 GeV/$c$, 5 - 10 GeV/$c$ and 5 - 10 GeV/$c$ for $h^{\pm}$ - $h^{\pm}$, $\pi^0$ - $h^{\pm}$ and $\pi^{\pm}$ - $h^{\pm}$ assorted-$p_T$ correlations.

The extracted $\sqrt{<j^2_T>}$ for minimum-bias $d$+Au collisions from all four dihadron correlations. The trigger $p_T$ ranges are 3 - 5 GeV/$c$, 5 - 10 GeV/$c$ and 5 - 10 GeV/$c$ for $h^{\pm}$ - $h^{\pm}$, $\pi^0$ - $h^{\pm}$ and $\pi^{\pm}$ - $h^{\pm}$ assorted-$p_T$ correlations.

$<sin^2(\phi_{jj})>$ for minimum bias $d$+Au collisions as function of associated particle $p_T$ for $\pi^{\pm}$ - $h^{\pm}$ and $\pi^0$ - $h^{\pm}$ correlations where the trigger particles have a $p_T$ between 5 - 10 GeV/$c$.

$<sin^2(\phi_{jj})>$ for minimum bias $d$+Au collisions as function of associated particle $p_T$ for $\pi^{\pm}$ - $h^{\pm}$ and $\pi^0$ - $h^{\pm}$ correlations where the trigger particles have a $p_T$ between 5 - 10 GeV/$c$.

$<sin^2(\phi_{jj})>$ for minimum bias $d$+Au collisions as function of trigger particle $p_T$ for $\pi^{\pm}$ - $h^{\pm}$ and $\pi^0$ - $h^{\pm}$ correlations.

$<sin^2(\phi_{jj})>$ for minimum bias $d$+Au collisions as function of trigger particle $p_T$ for $\pi^{\pm}$ - $h^{\pm}$ and $\pi^0$ - $h^{\pm}$ correlations.

Near-side and far-side $p_{out}$ distributions for minimum bias $d$+Au collisions obtained from $\pi^{\pm}$ - $h^{\pm}$ correlation. The trigger is 5 < $p_{T,trig}$ < 10 GeV/$c$, the associated particle is 0.5 < $p_{T,assoc}$ < 5.0 GeV/$c$.

Conditional yield as a function of $x_E$ for near-side and far-side for $\pi^{\pm}$ - $h^{\pm}$ and $\pi^0$ - $h^{\pm}$ from minimum-bias $d$+Au collisions.

Conditional yield as a function of $x_E$ for near-side and far-side for $\pi^{\pm}$ - $h^{\pm}$ and $\pi^0$ - $h^{\pm}$ from minimum-bias $d$+Au collisions.

Conditional yield as a function of $x_E$ for near-side and far-side correlations for $\pi^{\pm}$ - $h^{\pm}$ correlations from minimum-bias $d$+Au collisions. The trigger pions are 5 < $p_{T,trig}$ < 6 GeV/$c$.

Conditional yield as a function of $x_E$ for near-side and far-side correlation for $\pi^{\pm}$ - $h^{\pm}$ correlation for several different trigger $p_T$s for minimum-bias $d$+Au collisions.

Far-side conditional yield as a function of $p_{T,trig}$ for different ranges of $x_E$ for $\pi^{\pm}$ - $h^{\pm}$ correlation for minimum-bias $d$+Au collisions.

Far-side conditional yield as a function of $p_{T,trig}$ for different ranges of $x_E$ for $\pi^{\pm}$ - $h^{\pm}$ correlation for minimum-bias $d$+Au collisions.

Far-side conditional yield as a function of $p_{T,trig}$ for different ranges of $x_E$ for $\pi^{\pm}$ - $h^{\pm}$ correlation for minimum-bias $d$+Au collisions.

Far-side conditional yield as a function of $p_{T,trig}$ for different ranges of $x_E$ for $\pi^{\pm}$ - $h^{\pm}$ correlation for minimum-bias $d$+Au collisions.

The comparison of the $\sqrt{<j^2_T>}$ values between $d$+Au and $p$+$p$ for the $\pi^{\pm}$ - $h^{\pm}$ correlations and $\pi^0$ - $h^{\pm}$ correlations. The trigger pion range is 5 - 10 GeV/$c$.

The comparison of the $\sqrt{<j^2_T>}$ values between $d$+Au and $p$+$p$ for the $\pi^{\pm}$ - $h^{\pm}$ correlations and $\pi^0$ - $h^{\pm}$ correlations. The trigger pion range is 5 - 10 GeV/$c$.

Quadrature difference between minimum-bias $d$+Au and $p$+$p$ $<sin^2(\phi_{jj})>$ values.

Quadrature difference between minimum-bias $d$+Au and $p$+$p$ $<sin^2(\phi_{jj})>$ values.

The comparison of the $p_{out}$ distribution at the near-side and far-side between central $d$+Au collisions and $p$+$p$ collisions. Results are obtained for $\pi^{\pm}$ - $h^{\pm}$ correlations with the associated hadron range 0.5 < $p_{T,assoc}$ < 5 GeV/$c$ and trigger pion range of 5 < $p_{T,trig}$ < 10 GeV/$c$.

The comparison of the $x_E$ distribution from $\pi^{\pm}$ - $h^{\pm}$ correlation at the near-side and far-side between minimum bias $d$+Au collisions and $p$+$p$ collisions. The trigger $\pi^{\pm}$ are from 5 - 10 GeV/$c$.

Far-side conditional yield as a function of $p_{T,trig}$ for different ranges of $x_E$ from $p$+$p$ collisions, triggers are $\pi^{\pm}$ from 5 - 10 GeV/$c$.

Far-side conditional yield as a function of $p_{T,trig}$ for different ranges of $x_E$ from $p$+$p$ collisions, triggers are $\pi^{\pm}$ from 5 - 10 GeV/$c$.

Far-side conditional yield as a function of $p_{T,trig}$ for different ranges of $x_E$ from $p$+$p$ collisions, triggers are $\pi^{\pm}$ from 5 - 10 GeV/$c$.

The fitted fractional change in $dN/d_{x_E}$ per unit $p_{T,trig}$ of the far-side conditional yield for different ranges of $x_E$.

Centrality dependent near and far $CY(p_T)$ for $\pi^{\pm}$ - $h^{\pm}$ correlations and $\pi^0$ - $h^{\pm}$ correlations.

Centrality dependent near and far $CY(p_T)$ for $\pi^{\pm}$ - $h^{\pm}$ correlations and $\pi^0$ - $h^{\pm}$ correlations.

Centrality dependence to the ratio of near and far $CY(p_T)$ for $d$+Au to $p$+$p$.

Centrality dependence to the ratio of near and far $CY(p_T)$ for $d$+Au to $p$+$p$.

Near- and far-side widths as a function of $p_{T, trig}$ for charged pion triggers and associated charged hadrons (2–4.5 Gev/$c$) from minimum-bias $d$+Au collisions.

Near- and far-side widths as a function of $p_{T, trig}$ for neutral pion triggers and associated charged hadrons (2–4.5 Gev/$c$) from minimum-bias $d$+Au collisions.

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

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

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

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

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

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

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

Charged hadron $v_2$($p_T$) in Au+Au collisions at $\sqrt{s_{NN}}$ = 200 GeV from the two-particle cumulant method, the BBC event plane, and the ZDC-SMD event plane for the indicated centralities.

Charged hadron $v_2$($p_T$) in Au+Au collisions at $\sqrt{s_{NN}}$ = 200 GeV from the two-particle cumulant method, the BBC event plane, and the ZDC-SMD event plane for the indicated centralities.

Charged hadron $v_2$($p_T$) in Au+Au collisions at $\sqrt{s_{NN}}$ = 200 GeV from the two-particle cumulant method, the BBC event plane, and the ZDC-SMD event plane for the indicated centralities.

Charged hadron $v_2$($p_T$) in Au+Au collisions at $\sqrt{s_{NN}}$ = 200 GeV from the two-particle cumulant method, the BBC event plane, and the ZDC-SMD event plane for the indicated centralities.

Charged hadron $v_2$($p_T$) in Au+Au collisions at $\sqrt{s_{NN}}$ = 200 GeV from the two-particle cumulant method, the BBC event plane, and the ZDC-SMD event plane for the indicated centralities.

Charged hadron $v_2$($p_T$) in Au+Au collisions at $\sqrt{s_{NN}}$ = 200 GeV from the two-particle cumulant method, the BBC event plane, and the ZDC-SMD event plane for the indicated centralities.

Charged hadron $v_2$($p_T$) in Au+Au collisions at $\sqrt{s_{NN}}$ = 200 GeV from the two-particle cumulant method, the BBC event plane, and the ZDC-SMD event plane for the indicated centralities.

Charged hadron $v_2$($p_T$) in Au+Au collisions at $\sqrt{s_{NN}}$ = 200 GeV from the two-particle cumulant method, the BBC event plane, and the ZDC-SMD event plane for the indicated centralities.

Charged hadron $v_2$($p_T$) in Au+Au collisions at $\sqrt{s_{NN}}$ = 200 GeV from the two-particle cumulant method, the BBC event plane, and the ZDC-SMD event plane for the indicated centralities.

Charged hadron $v_2$($p_T$) in Au+Au collisions at $\sqrt{s_{NN}}$ = 200 GeV from the two-particle cumulant method, the BBC event plane, and the ZDC-SMD event plane for the indicated centralities.

Charged hadron $v_2$($p_T$) in Au+Au collisions at $\sqrt{s_{NN}}$ = 200 GeV from the two-particle cumulant method, the BBC event plane, and the ZDC-SMD event plane for the indicated centralities.

Anisotropy parameter $v_2$ as a function of pseudorapidity within the PHENIX central arms using event planes from the BBC and ZDC-SMD, and from the two-particle cumulant method for centrality 20–40%. The results are shown for three $p_T$ bins, which are from top to bottom: 2.0–3.0, 1.2–1.4 and 0.6–0.8 GeV/$c$.

Comparison of the $v_2${BBC} and $v_2${ZDC-SMD} obtained from the S-N and ZDC-BBC-CNT subevents as a function of pT in the 20–60% centrality range.

Charged hadron $v_2${2} for the PHENIX experiment as a function of $p_T$ in centrality 20–60%.

The PHENIX $v_2${BBC} and $v_2${ZDC-SMD} from the modified event-plane method for charged hadrons as a function of $p_T$ in different centralities.

The PHENIX $v_2${BBC} and $v_2${ZDC-SMD} from the modified event-plane method for charged hadrons as a function of $p_T$ in different centralities.

The PHENIX $v_2${BBC} and $v_2${ZDC-SMD} from the modified event-plane method for charged hadrons as a function of $p_T$ in different centralities.

The PHENIX $v_2${BBC} and $v_2${ZDC-SMD} from the modified event-plane method for charged hadrons as a function of $p_T$ in different centralities.

The PHENIX $v_2${BBC} and $v_2${ZDC-SMD} from the modified event-plane method for charged hadrons as a function of $p_T$ in different centralities.

The PHENIX $v_2${BBC} and $v_2${ZDC-SMD} from the modified event-plane method for charged hadrons as a function of $p_T$ in different centralities.

Additional information containing number of events which were used to reconstruct the numvers matching to Figure 1 and 2.

Additional information containing number of events which were used to reconstruct the numvers matching to Figure 1 and 2.

Additional information containing number of events which were used to reconstruct the numvers matching to Figure 1 and 2.

The uncorrected multiplicity distributions of charged hadrons with 0.2 < $p_T$ < 2.0 GeV/$c$ for 200 GeV Au+Au collisions.

The uncorrected multiplicity distributions of charged hadrons with 0.2 < $p_T$ < 2.0 GeV/$c$ for 200 GeV Au+Au collisions.

The uncorrected multiplicity distributions of charged hadrons with 0.2 < $p_T$ < 2.0 GeV/$c$ for 200 GeV Au+Au collisions.

The uncorrected multiplicity distributions of charged hadrons with 0.2 < $p_T$ < 2.0 GeV/$c$ for 200 GeV Au+Au collisions.

The uncorrected multiplicity distributions of charged hadrons with 0.2 < $p_T$ < 2.0 GeV/$c$ for 200 GeV Au+Au collisions.

The uncorrected multiplicity distributions of charged hadrons with 0.2 < $p_T$ < 2.0 GeV/$c$ for 200 GeV Au+Au collisions.

The uncorrected multiplicity distributions of charged hadrons with 0.2 < $p_T$ < 2.0 GeV/$c$ for 200 GeV Au+Au collisions.

The uncorrected multiplicity distributions of charged hadrons with 0.2 < $p_T$ < 2.0 GeV/$c$ for 200 GeV Au+Au collisions.

The uncorrected multiplicity distributions of charged hadrons with 0.2 < $p_T$ < 2.0 GeV/$c$ for 200 GeV Au+Au collisions.

The uncorrected multiplicity distributions of charged hadrons with 0.2 < $p_T$ < 2.0 GeV/$c$ for 200 GeV Au+Au collisions.

The uncorrected multiplicity distributions of charged hadrons with 0.2 < $p_T$ < 2.0 GeV/$c$ for 200 GeV Au+Au collisions.

The uncorrected multiplicity distributions of charged hadrons with 0.2 < $p_T$ < 2.0 GeV/$c$ for 200 GeV Au+Au collisions.

The uncorrected multiplicity distributions of charged hadrons with 0.2 < $p_T$ < 2.0 GeV/$c$ for 62.4 GeV Au+Au collisions.

The uncorrected multiplicity distributions of charged hadrons with 0.2 < $p_T$ < 2.0 GeV/$c$ for 62.4 GeV Au+Au collisions.

The uncorrected multiplicity distributions of charged hadrons with 0.2 < $p_T$ < 2.0 GeV/$c$ for 62.4 GeV Au+Au collisions.

The uncorrected multiplicity distributions of charged hadrons with 0.2 < $p_T$ < 2.0 GeV/$c$ for 62.4 GeV Au+Au collisions.

The uncorrected multiplicity distributions of charged hadrons with 0.2 < $p_T$ < 2.0 GeV/$c$ for 62.4 GeV Au+Au collisions.

The uncorrected multiplicity distributions of charged hadrons with 0.2 < $p_T$ < 2.0 GeV/$c$ for 62.4 GeV Au+Au collisions.

The uncorrected multiplicity distributions of charged hadrons with 0.2 < $p_T$ < 2.0 GeV/$c$ for 62.4 GeV Au+Au collisions.

The uncorrected multiplicity distributions of charged hadrons with 0.2 < $p_T$ < 2.0 GeV/$c$ for 62.4 GeV Au+Au collisions.

The uncorrected multiplicity distributions of charged hadrons with 0.2 < $p_T$ < 2.0 GeV/$c$ for 62.4 GeV Au+Au collisions.

The uncorrected multiplicity distributions of charged hadrons with 0.2 < $p_T$ < 2.0 GeV/$c$ for 62.4 GeV Au+Au collisions.

The uncorrected multiplicity distributions of charged hadrons with 0.2 < $p_T$ < 2.0 GeV/$c$ for 200 GeV Cu+Cu collisions.

The uncorrected multiplicity distributions of charged hadrons with 0.2 < $p_T$ < 2.0 GeV/$c$ for 200 GeV Cu+Cu collisions.

The uncorrected multiplicity distributions of charged hadrons with 0.2 < $p_T$ < 2.0 GeV/$c$ for 200 GeV Cu+Cu collisions.

The uncorrected multiplicity distributions of charged hadrons with 0.2 < $p_T$ < 2.0 GeV/$c$ for 200 GeV Cu+Cu collisions.

The uncorrected multiplicity distributions of charged hadrons with 0.2 < $p_T$ < 2.0 GeV/$c$ for 200 GeV Cu+Cu collisions.

The uncorrected multiplicity distributions of charged hadrons with 0.2 < $p_T$ < 2.0 GeV/$c$ for 200 GeV Cu+Cu collisions.

The uncorrected multiplicity distributions of charged hadrons with 0.2 < $p_T$ < 2.0 GeV/$c$ for 200 GeV Cu+Cu collisions.

The uncorrected multiplicity distributions of charged hadrons with 0.2 < $p_T$ < 2.0 GeV/$c$ for 62.4 GeV Cu+Cu collisions.

The uncorrected multiplicity distributions of charged hadrons with 0.2 < $p_T$ < 2.0 GeV/$c$ for 62.4 GeV Cu+Cu collisions.

The uncorrected multiplicity distributions of charged hadrons with 0.2 < $p_T$ < 2.0 GeV/$c$ for 62.4 GeV Cu+Cu collisions.

The uncorrected multiplicity distributions of charged hadrons with 0.2 < $p_T$ < 2.0 GeV/$c$ for 62.4 GeV Cu+Cu collisions.

The uncorrected multiplicity distributions of charged hadrons with 0.2 < $p_T$ < 2.0 GeV/$c$ for 62.4 GeV Cu+Cu collisions.

The uncorrected multiplicity distributions of charged hadrons with 0.2 < $p_T$ < 2.0 GeV/$c$ for 22.5 GeV Cu+Cu collisions.

The uncorrected multiplicity distributions of charged hadrons with 0.2 < $p_T$ < 2.0 GeV/$c$ for 22.5 GeV Cu+Cu collisions.

The uncorrected multiplicity distributions of charged hadrons with 0.2 < $p_T$ < 2.0 GeV/$c$ for 200 GeV $p$+$p$ collisions.

The mean from the NBD fit as a function of $N_{part}$ for 200 GeV Au+Au collisions over the range 0.2 < $p_T$ < 2.0 GeV/$c$.

The mean from the NBD fit as a function of $N_{part}$ for 62.4 GeV Au+Au collisions over the range 0.2 < $p_T$ < 2.0 GeV/$c$.

The mean from the NBD fit as a function of $N_{part}$ for 200 GeV Cu+Cu collisions over the range 0.2 < $p_T$ < 2.0 GeV/$c$.

The mean from the NBD fit as a function of $N_{part}$ for 62.4 GeV Cu+Cu collisions over the range 0.2 < $p_T$ < 2.0 GeV/$c$.

The mean from the NBD fit as a function of $N_{part}$ for 22.5 GeV Cu+Cu collisions over the range 0.2 < $p_T$ < 2.0 GeV/$c$.

Fluctuations expressed as the scaled variance as a function of centrality for 200 GeV Au+Au collisions in the range 0.2 < $p_T$ < 2.0 GeV/$c$.

Fluctuations expressed as the scaled variance as a function of $N_{part}$ for 200 GeV Au+Au collisions for 0.2 < $p_T$ < 2.0 GeV/$c$.

Fluctuations expressed as the scaled variance as a function of $N_{part}$ for 62.4 GeV Au+Au collisions for 0.2 < $p_T$ < 2.0 GeV/$c$.

Fluctuations expressed as the scaled variance as a function of $N_{part}$ for 200 GeV Cu+Cu collisions for 0.2 < $p_T$ < 2.0 GeV/$c$.

Fluctuations expressed as the scaled variance as a function of $N_{part}$ for 62.4 GeV Cu+Cu collisions for 0.2 < $p_T$ < 2.0 GeV/$c$.

Fluctuations expressed as the scaled variance as a function of $N_{part}$ for 22.5 GeV Cu+Cu collisions for 0.2 < $p_T$ < 2.0 GeV/$c$.

Scaled variance for 200 GeV Au+Au collisions plotted as a function of $N_{part}$ for several $p_T$ ranges.

Scaled variance for 200 GeV Au+Au collisions plotted as a function of $N_{part}$ for several $p_T$ ranges.

Scaled variance for 62.4 GeV Au+Au collisions plotted as a function of $N_{part}$ for several $p_T$ ranges.

Scaled variance for 62.4 GeV Au+Au collisions plotted as a function of $N_{part}$ for several $p_T$ ranges.

Scaled variance for 62.4 GeV Au+Au collisions plotted as a function of $N_{part}$ for several $p_T$ ranges.

Scaled variance for 200 GeV Cu+Cu collisions plotted as a function of $N_{part}$ for several $p_T$ ranges.

Scaled variance for 200 GeV Cu+Cu collisions plotted as a function of $N_{part}$ for several $p_T$ ranges.

Scaled variance for 200 GeV Cu+Cu collisions plotted as a function of $N_{part}$ for several $p_T$ ranges.

Scaled variance for 62.4 GeV Cu+Cu collisions plotted as a function of $N_{part}$ for several $p_T$ ranges.

Scaled variance for 62.4 GeV Cu+Cu collisions plotted as a function of $N_{part}$ for several $p_T$ ranges.

The inverse of the parameter $k_{NBD}$ from the Negative Binomial Distribution fits for 200 GeV Au+Au collisions. The fluctuations are plotted as a function of $N_{part}$ for several $p_T$ ranges.

The inverse of the parameter $k_{NBD}$ from the Negative Binomial Distribution fits for 200 GeV Au+Au collisions. The fluctuations are plotted as a function of $N_{part}$ for several $p_T$ ranges.

The inverse of the parameter $k_{NBD}$ from the Negative Binomial Distribution fits for 62.4 GeV Au+Au collisions. The fluctuations are plotted as a function of $N_{part}$ for several $p_T$ ranges.

The inverse of the parameter $k_{NBD}$ from the Negative Binomial Distribution fits for 62.4 GeV Au+Au collisions. The fluctuations are plotted as a function of $N_{part}$ for several $p_T$ ranges.

The inverse of the parameter $k_{NBD}$ from the Negative Binomial Distribution fits for 62.4 GeV Au+Au collisions. The fluctuations are plotted as a function of $N_{part}$ for several $p_T$ ranges.

The inverse of the parameter $k_{NBD}$ from the Negative Binomial Distribution fits for 200 GeV Cu+Cu collisions. The fluctuations are plotted as a function of $N_{part}$ for several $p_T$ ranges.

The inverse of the parameter $k_{NBD}$ from the Negative Binomial Distribution fits for 200 GeV Cu+Cu collisions. The fluctuations are plotted as a function of $N_{part}$ for several $p_T$ ranges.

The inverse of the parameter $k_{NBD}$ from the Negative Binomial Distribution fits for 200 GeV Cu+Cu collisions. The fluctuations are plotted as a function of $N_{part}$ for several $p_T$ ranges.

The inverse of the parameter $k_{NBD}$ from the Negative Binomial Distribution fits for 62.4 GeV Cu+Cu collisions. The fluctuations are plotted as a function of $N_{part}$ for several $p_T$ ranges.

The inverse of the parameter $k_{NBD}$ from the Negative Binomial Distribution fits for 62.4 GeV Cu+Cu collisions. The fluctuations are plotted as a function of $N_{part}$ for several $p_T$ ranges.

The scaled variance as a funciton of $N_{part}$ for 200 GeV Au+Au collisions in the range 0.2 < $p_T$ < 2.0 GeV/$c$.

The correlation of the clan model parameters $\bar{n_c}$ and $\bar{N_c}$ for all of the collision species measured as a function of centrality.

The correlation of the clan model parameters $\bar{n_c}$ and $\bar{N_c}$ for all of the collision species measured as a function of centrality.

The correlation of the clan model parameters $\bar{n_c}$ and $\bar{N_c}$ for all of the collision species measured as a function of centrality.

The correlation of the clan model parameters $\bar{n_c}$ and $\bar{N_c}$ for all of the collision species measured as a function of centrality.

The correlation of the clan model parameters $\bar{n_c}$ and $\bar{N_c}$ for all of the collision species measured as a function of centrality.

The correlation of the clan model parameters $\bar{n_c}$ and $\bar{N_c}$ for all of the collision species measured as a function of centrality.

The correlation of the clan model parameters $\bar{n_c}$ and $\bar{N_c}$ for all of the collision species measured as a function of centrality.

The correlation of the clan model parameters $\bar{n_c}$ and $\bar{N_c}$ for all of the collision species measured as a function of centrality.

The correlation of the clan model parameters $\bar{n_c}$ and $\bar{N_c}$ for all of the collision species measured as a function of centrality.

The correlation of the clan model parameters $\bar{n_c}$ and $\bar{N_c}$ for all of the collision species measured as a function of centrality.

The correlation of the clan model parameters $\bar{n_c}$ and $\bar{N_c}$ for all of the collision species measured as a function of centrality.

$1/{N_{trig}}$ ${dN}/{d\Delta\phi}$ distributions for charge selected $\bar{p}$ and $p$ triggers both with associated $p$ for six centrality bins. Triggers have 2.5 < $p_T$ < 4.0 GeV/$c$ and associated particles have 1.8 < $p_T$ < 2.5 GeV/$c$.

$1/{N_{trig}}$ ${dN}/{d\Delta\phi}$ distributions for charge selected $\bar{p}$ and $p$ triggers both with associated $p$ for six centrality bins. Triggers have 2.5 < $p_T$ < 4.0 GeV/$c$ and associated particles have 1.8 < $p_T$ < 2.5 GeV/$c$.

$1/{N_{trig}}$ ${dN}/{d\Delta\phi}$ distributions for charge-inclusive baryon and meson triggers and associated mesons for six centrality bins. Triggers have 2.5 < $p_T$ < 4.0 GeV/$c$ and associated particles have 1.8 < $p_T$ < 2.5 GeV/$c$.

Conditional yields per trigger on the near side for charge selected and charge selected $p$ and $\bar{p}$ correlations. Triggers have 2.5 < $p_T$ < 4.0 GeV/$c$ and associated particles have 1.8 < $p_T$ < 2.5 GeV/$c$. There is an 11.4% additional normalization error on baryon associated particle points and 8.9% each on the $p$ and $\bar{p}$ associated particle points.

Conditional yields per trigger on the away side for charge selected and charge selected $p$ and $\bar{p}$ correlations. Triggers have 2.5 < $p_T$ < 4.0 GeV/$c$ and associated particles have 1.8 < $p_T$ < 2.5 GeV/$c$.

Conditional yields per trigger for baryon and meson triggers with associated mesons. Triggers have 2.5 < $p_T$ < 4.0 GeV/$c$ and associated particles have 1.8 < $p_T$ < 2.5 GeV/$c$. There is an additional 13.6% normalization error.

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

$v_2$ vs $p_T$ for identified particle species obtained in minimum bias Au+Au collisions.

$v_2$ vs $p_T$ for identified particle species obtained in minimum bias Au+Au collisions.

$v_2$ vs $p_T$ for identified particle species obtained in minimum bias Au+Au collisions.

$v_2$ vs $KE_T$ for identified particle species obtained in minimum bias Au+Au collisions.

$v_2$ vs $KE_T$ for identified particle species obtained in minimum bias Au+Au collisions.

$v_2$ vs $KE_T$ for identified particle species obtained in minimum bias Au+Au collisions.

${v_2}/{n_q}$ vs ${p_T}/{n_q}$ for identified particle species obtained in minimum bias Au+Au collisions.

${v_2}/{n_q}$ vs ${p_T}/{n_q}$ for identified particle species obtained in minimum bias Au+Au collisions.

${v_2}/{n_q}$ vs ${p_T}/{n_q}$ for identified particle species obtained in minimum bias Au+Au collisions.

${v_2}/{n_q}$ vs ${KE_T}/{n_q}$ for identified particle species obtained in minimum bias Au+Au collisions.

${v_2}/{n_q}$ vs ${KE_T}/{n_q}$ for identified particle species obtained in minimum bias Au+Au collisions.

${v_2}/{n_q}$ vs ${KE_T}/{n_q}$ for identified particle species obtained in minimum bias Au+Au collisions.

$J/\psi$ $R_{CP}$ for 0-20%, 20-40%, and 40-60% (central) relative to 60-86% (peripheral) Au+Au collisions at 62.4 GeV. Type A errors are point-to-point uncorrelated, Type B errors are point-to-point correlated, and Type C are global systematic from the uncertainty in the peripheral $<N_{coll}>$ value.

$J/\psi$ $R_{CP}$ for 0-40% (central) relative to 40-86% (peripheral) Au+Au collisions at 39 GeV. Type A errors are point-to-point uncorrelated, Type B errors are point-to-point correlated, and Type C are global systematic from the uncertainty in the peripheral $<N_{coll}>$ value.

$J/\psi$ $R_{AA}$ at $\sqrt{s_{NN}}$ = 39 and 62.4 GeV. Type A errors are point-to-point uncorrelated, Type B errors are point-to-point correlated, and Type C are global systematic from the uncertainty in the peripheral $<N_{coll}>$ value.

$J/\psi$ $R_{AA}$ at $\sqrt{s_{NN}}$ = 39 and 62.4 GeV.

Parameters of the Tsallis fit with parameter $n$ constrained to a fixed linear dependence on mass (for mesons). Cross sections are in $\mu$b for $J/\psi$ and $\psi^{\prime}$ and in mb for all other particles.

Constant and linear fits to the Tsallis parameter $T$ of mesons with fixed parameter $n$. The last column (Prob.) gives the probability estimated by the $\chi^2$/$n.d.f.$ of the fit.

Parameters of the Tsallis fit with Eq. 8 in the paper with parameters $n$ and $T$ constrained to have a fixed linear dependence on mass (for mesons). Cross sections are in $\mu$b for $J/\psi$ and $\psi^{\prime}$ and in mb for all other particles.

Cross sections in mb and $<p_T>$ in GeV/$c$ of different particles in $p$+$p$ collisions at $\sqrt{s}$ = 200 GeV. PHENIX columns show the values obtained by fits to the experiment spectra with the Tsallis functional form as described in the text. One should state explicitly that these values do not supersede values given in the ”Published” column by the experiments, in their publications listed in the last column, or elsewhere. An additional 9.7% systematic uncertainty should be added to all $d\sigma$/$dy$ values listed in the column “PHENIX” to account for the trigger uncertainties. Values in the column “Published” are also given without these systematic uncertainties. The column “S.M.” is the prediction of the statistical model discussed in the text. The characteristic widths of the particle spectra are $\langle p_T \rangle$ for all species except for $J/\psi$ and $\psi^{\prime}$ for which the values given in the Table are $\langle p_T^2 \rangle$. For $\psi^{\prime}$ the integration is done in the $p_T$ regionbelow 5 GeV/$c$. All errors are the combined statistical and systematic uncertainties.

Cross sections in mb and $<p_T>$ in GeV/$c$ of different particles in $p$+$p$ collisions at $\sqrt{s}$ = 200 GeV. PHENIX columns show the values obtained by fits to the experiment spectra with the Tsallis functional form as described in the text. One should state explicitly that these values do not supersede values given in the ”Published” column by the experiments, in their publications listed in the last column, or elsewhere. An additional 9.7% systematic uncertainty should be added to all $d\sigma$/$dy$ values listed in the column “PHENIX” to account for the trigger uncertainties. Values in the column “Published” are also given without these systematic uncertainties. The column “S.M.” is the prediction of the statistical model discussed in the text. The characteristic widths of the particle spectra are $\langle p_T \rangle$ for all species except for $J/\psi$ and $\psi^{\prime}$ for which the values given in the Table are $\langle p_T^2 \rangle$. For $\psi^{\prime}$ the integration is done in the $p_T$ regionbelow 5 GeV/$c$. All errors are the combined statistical and systematic uncertainties.

The invariant differential cross section ${d^2 \sigma}$/${2\pi p_T dy p_T}$ of $\omega$ meson production measured in the indicated decay channel. Type A uncertainties are $p_T$ independent. Type C & B systematic uncertainties are residual with Type B having an unknown $p_T$ dependence. For each $y$-coordinate in each $p_T$ bin was varied by the same amount according to the uncertainty Type C.

The invariant differential cross section ${d^2 \sigma}$/${2\pi p_T dy p_T}$ of $\omega$ meson production measured in the indicated decay channel. Type A uncertainties are $p_T$ independent. Type C & B systematic uncertainties are residual with Type B having an unknown $p_T$ dependence. For each $y$-coordinate in each $p_T$ bin was varied by the same amount according to the uncertainty Type C.

The invariant differential cross section ${d^2 \sigma}$/${2\pi p_T dy p_T}$ of $\omega$ meson production measured in the indicated decay channel. Type A uncertainties are $p_T$ independent. Type C & B systematic uncertainties are residual with Type B having an unknown $p_T$ dependence. For each $y$-coordinate in each $p_T$ bin was varied by the same amount according to the uncertainty Type C.

The invariant differential cross section ${d^2 \sigma}$/${2\pi p_T dy p_T}$ of $\omega$ meson production measured in the indicated decay channel. Type A uncertainties are $p_T$ independent. Type C & B systematic uncertainties are residual with Type B having an unknown $p_T$ dependence. For each $y$-coordinate in each $p_T$ bin was varied by the same amount according to the uncertainty Type C.

The invariant differential cross section ${d^2 \sigma}$/${2\pi p_T dy p_T}$ of $\omega$ meson production measured in the indicated decay channel. Type A uncertainties are $p_T$ independent. Type C & B systematic uncertainties are residual with Type B having an unknown $p_T$ dependence. For each $y$-coordinate in each $p_T$ bin was varied by the same amount according to the uncertainty Type C.

The invariant differential cross section ${d^2 \sigma}$/${2\pi p_T dy p_T}$ of $\omega$ meson production measured in the indicated decay channel. Type A uncertainties are $p_T$ independent. Type C & B systematic uncertainties are residual with Type B having an unknown $p_T$ dependence. For each $y$-coordinate in each $p_T$ bin was varied by the same amount according to the uncertainty Type C.

The invariant differential cross section ${d^2 \sigma}$/${2\pi p_T dy p_T}$ of $\omega$ meson production measured in the indicated decay channel. Type A uncertainties are $p_T$ independent. Type C & B systematic uncertainties are residual with Type B having an unknown $p_T$ dependence. For each $y$-coordinate in each $p_T$ bin was varied by the same amount according to the uncertainty Type C.

$p_T$ dependence of $v_2$ for charged hadrons for several centrality selections as indicated.

$p_T$ dependence of $v_4$ for charged hadrons for several centrality selections as indicated.

$v_2$ vs. $N_{part}$ and $v_4$ vs. $N_{part}$ for charged hadrons for several $p_T$ selections as indicated.

The ratio ${v_4}/{(v_2)^2}$ vs. $N_{part}$ for the same $p_T$ selections.

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

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

Collision centrality, energy, and system size dependence of shape parameters.

Collision centrality, energy, and system size dependence of shape parameters.

Collision centrality, energy, and system size dependence of shape parameters.

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

Identified hadron $v_2$ in midcentral (20–60%, right panels) Au + Au collisions at $\sqrt{s_{NN}}$ = 200 GeV. Panels (a) and (b) show $v_2$ as a function of transverse momentum $p_T$. The $v_2$ of all species for centrality 0–20% has been scaled up by a factor of 1.6 for better comparison with results of 20–60% centrality. The error bars (shaded boxes) represent the statistical (systematic) uncertainties. The systematic uncertainties shown are type A and B only.

Identified hadron $v_2$ in midcentral (20–60%, right panels) Au + Au collisions at $\sqrt{s_{NN}}$ = 200 GeV. Panels (a) and (b) show $v_2$ as a function of transverse momentum $p_T$. The $v_2$ of all species for centrality 0–20% has been scaled up by a factor of 1.6 for better comparison with results of 20–60% centrality. The error bars (shaded boxes) represent the statistical (systematic) uncertainties. The systematic uncertainties shown are type A and B only.

Identified hadron $v_2$ in midcentral (20–60%, right panels) Au + Au collisions at $\sqrt{s_{NN}}$ = 200 GeV. Panels (a) and (b) show $v_2$ as a function of transverse momentum $p_T$. The $v_2$ of all species for centrality 0–20% has been scaled up by a factor of 1.6 for better comparison with results of 20–60% centrality. The error bars (shaded boxes) represent the statistical (systematic) uncertainties. The systematic uncertainties shown are type A and B only.

The quark-number-scaled $v_2$ ($v_2/n_q$) of identified hadrons are shown as a function of the kinetic energy per quark, KE$_T/n_q$ in 0–10% centrality [panel (a)] in Au + Au collisions at $\sqrt{s_{NN}}$ = 200 GeV. The error bars (shaded boxes) represent the statistical (systematic) uncertainties. The systematic uncertainties shown are type A and B only.

The quark-number-scaled $v_2$ ($v_2/n_q$) of identified hadrons are shown as a function of the kinetic energy per quark, KE$_T/n_q$ in 0–10% centrality [panel (a)] in Au + Au collisions at $\sqrt{s_{NN}}$ = 200 GeV. The error bars (shaded boxes) represent the statistical (systematic) uncertainties. The systematic uncertainties shown are type A and B only.

The quark-number-scaled $v_2$ ($v_2/n_q$) of identified hadrons are shown as a function of the kinetic energy per quark, KE$_T/n_q$ in 0–10% centrality [panel (a)] in Au + Au collisions at $\sqrt{s_{NN}}$ = 200 GeV. The error bars (shaded boxes) represent the statistical (systematic) uncertainties. The systematic uncertainties shown are type A and B only.

The quark-number-scaled $v_2$ ($v_2/n_q$) of identified hadrons are shown as a function of the kinetic energy per quark, KE$_T/n_q$ in 10–20% centrality [panel (b)] in Au + Au collisions at $\sqrt{s_{NN}}$ = 200 GeV. The error bars (shaded boxes) represent the statistical (systematic) uncertainties. The systematic uncertainties shown are type A and B only.

The quark-number-scaled $v_2$ ($v_2/n_q$) of identified hadrons are shown as a function of the kinetic energy per quark, KE$_T/n_q$ in 10–20% centrality [panel (b)] in Au + Au collisions at $\sqrt{s_{NN}}$ = 200 GeV. The error bars (shaded boxes) represent the statistical (systematic) uncertainties. The systematic uncertainties shown are type A and B only.

The quark-number-scaled $v_2$ ($v_2/n_q$) of identified hadrons are shown as a function of the kinetic energy per quark, KE$_T/n_q$ in 10–20% centrality [panel (b)] in Au + Au collisions at $\sqrt{s_{NN}}$ = 200 GeV. The error bars (shaded boxes) represent the statistical (systematic) uncertainties. The systematic uncertainties shown are type A and B only.

The quark-number-scaled $v_2$ ($v_2/n_q$) of identified hadrons are shown as a function of the kinetic energy per quark, KE$_T/n_q$ in 20–40% centrality [panel (c)] in Au + Au collisions at $\sqrt{s_{NN}}$ = 200 GeV. The error bars (shaded boxes) represent the statistical (systematic) uncertainties. The systematic uncertainties shown are type A and B only.

The quark-number-scaled $v_2$ ($v_2/n_q$) of identified hadrons are shown as a function of the kinetic energy per quark, KE$_T/n_q$ in 20–40% centrality [panel (c)] in Au + Au collisions at $\sqrt{s_{NN}}$ = 200 GeV. The error bars (shaded boxes) represent the statistical (systematic) uncertainties. The systematic uncertainties shown are type A and B only.

The quark-number-scaled $v_2$ ($v_2/n_q$) of identified hadrons are shown as a function of the kinetic energy per quark, KE$_T/n_q$ in 20–40% centrality [panel (c)] in Au + Au collisions at $\sqrt{s_{NN}}$ = 200 GeV. The error bars (shaded boxes) represent the statistical (systematic) uncertainties. The systematic uncertainties shown are type A and B only.

The quark-number-scaled $v_2$ ($v_2/n_q$) of identified hadrons are shown as a function of the kinetic energy per quark, KE$_T/n_q$ in 40–60% centrality [panel (d)] in Au + Au collisions at $\sqrt{s_{NN}}$ = 200 GeV. The error bars (shaded boxes) represent the statistical (systematic) uncertainties. The systematic uncertainties shown are type A and B only.

The quark-number-scaled $v_2$ ($v_2/n_q$) of identified hadrons are shown as a function of the kinetic energy per quark, KE$_T/n_q$ in 40–60% centrality [panel (d)] in Au + Au collisions at $\sqrt{s_{NN}}$ = 200 GeV. The error bars (shaded boxes) represent the statistical (systematic) uncertainties. The systematic uncertainties shown are type A and B only.

The quark-number-scaled $v_2$ ($v_2/n_q$) of identified hadrons are shown as a function of the kinetic energy per quark, KE$_T/n_q$ in 40–60% centrality [panel (d)] in Au + Au collisions at $\sqrt{s_{NN}}$ = 200 GeV. The error bars (shaded boxes) represent the statistical (systematic) uncertainties. The systematic uncertainties shown are type A and B only.

The quark-number-scaled $v_2$ ($v_2/n_q$) of identified hadrons are shown as a function of the kinetic energy per quark, KE$_T/n_q$ in 10–40% centrality [panel (b)] in Au + Au collisions at $\sqrt{s_{NN}}$ = 200 GeV. The $v_2$ of $\Lambda$ and K$^0_S$ are measured by STAR collaboration [21]. The error bars (open boxes) represent the statistical (systematic) uncertainties. The systematic uncertainties shown on the results from this study are type A and B only.

The quark-number-scaled $v_2$ ($v_2/n_q$) of identified hadrons are shown as a function of the kinetic energy per quark, KE$_T/n_q$ in 10–40% centrality [panel (b)] in Au + Au collisions at $\sqrt{s_{NN}}$ = 200 GeV. The $v_2$ of $\Lambda$ and K$^0_S$ are measured by STAR collaboration [21]. The error bars (open boxes) represent the statistical (systematic) uncertainties. The systematic uncertainties shown on the results from this study are type A and B only.

The quark-number-scaled $v_2$ ($v_2/n_q$) of identified hadrons are shown as a function of the kinetic energy per quark, KE$_T/n_q$ in 10–40% centrality [panel (b)] in Au + Au collisions at $\sqrt{s_{NN}}$ = 200 GeV. The $v_2$ of $\Lambda$ and K$^0_S$ are measured by STAR collaboration [21]. The error bars (open boxes) represent the statistical (systematic) uncertainties. The systematic uncertainties shown on the results from this study are type A and B only.

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

Transverse momentum dependence of $v_2$ for identified pions and protons within $|\eta|$ < 0.35 in 0-5% central $^3$He+Au collisions.

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.

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