Cross sections for final states consisting of two charm mesons and an extra pion.

Total charm cross sections for exclusive and inclusive final states.

The measured total hadronic cross section and R value with correction for initial state radiation.

Analyzing powers A_N(-x_F) in x_F bins at < eta > =3.3 and x_F < 0.

Analyzing powers versus pi^0 transverse momentum (p_T) for events with scaled pi^0 longitudinal momentum x_F > 0.4.

Analyzing powers versus pi^0 transverse momentum (p_T) for events with scaled pi^0 longitudinal momentum x_F < -0.4.

Analyzing powers versus pi^0 transverse momentum (p_T) in 0.25 < x_F < 0.30 bin.

Analyzing powers versus pi^0 transverse momentum (p_T) in 0.30 < x_F < 0.35 bin.

Analyzing powers versus pi^0 transverse momentum (p_T) in 0.35 < x_F < 0.40 bin.

Analyzing powers versus pi^0 transverse momentum (p_T) in 0.40 < x_F < 0.47 bin.

Analyzing powers versus pi^0 transverse momentum (p_T) in 0.47 < x_F < 0.56 bin.

Analyzing powers versus pi^0 transverse momentum (p_T) in 0.56 < x_F < 0.80 bin.

$v_{2}$ as a function of $p_{T}$ for charged hadrons with $|\eta| < 1.0$ in 10–40$%$ $Au+Au$ collisions, at $\sqrt{s_{NN}} = 200$ GeV, from the Lee-Yang Zero method (solid circles) with Sum Generating Function.

$v_{2}$ as a function of $p_{T}$ for charged hadrons with $|\eta| < 1.0$ in 0–80$%$ $Au+Au$ collisions, at $\sqrt{s_{NN}} = 200$ GeV, from the Event Plane method.

$p_{T}$ integrated charged hadron $v_{2}$ in the TPC as a function of geometrical cross section. Shown are the Event Plane method ($v_{2}${EP}) (open circles). All data are from $Au+Au$ collisions, at $\sqrt{s_{NN}} = 200$ GeV. For the TPC, $|\eta|< 1.0$ was used.

$p_{T}$ integrated charged hadron $v_{2}$ in the TPC as a function of geometrical cross section. Shown are the 4 Cumulant method ($v_{2}${4}) (solid squares). All data are from $Au+Au$ collisions, at $\sqrt{s_{NN}} = 200$ GeV. For the TPC, $|\eta|< 1.0$ was used.

$p_{T}$ integrated charged hadron $v_{2}$ in the TPC as a function of geometrical cross section. Shown are the Lee-Yang Zero method with Sum Generating Function ($v_{2}$) (solid circles). All data are from $Au+Au$ collisions, at $\sqrt{s_{NN}} = 200$ GeV. For the TPC, $|\eta|< 1.3$ was used.

$p_{T}$ integrated charged hadron $v_{2}$ in the TPC as a function of geometrical cross section. Shown are the Lee-Yang Zero method with Product Generating Function ($v_{2}$) (open stars). All data are from $Au+Au$ collisions, at $\sqrt{s_{NN}} = 200$ GeV. For the TPC, $|\eta|< 1.3$ was used.

$v_{2}$ as a function of $p_{T}$ for $10–40\%$ centrality using $\eta$-subevent method (open crosses), are shown in (a) for $\Lambda$. All data are from $\sqrt{s_{NN}} = 200$ GeV $Au+Au$ collisions.

$v_{2}$ as a function of $p_{T}$ for $10–40\%$ centrality using Event Plane method (open circles), are shown in (a) for $\Lambda$. All data are from $\sqrt{s_{NN}} = 200$ GeV $Au+Au$ collisions.

$v_{2}$ as a function of $p_{T}$ for $10–40\%$ centrality using $\eta$-subevent method (open crosses), are shown in (b) for $K_{S}^{0}$. All data are from $\sqrt{s_{NN}} = 200$ GeV $Au+Au$ collisions.

$v_{2}$ as a function of $p_{T}$ for $10–40\%$ centrality using Event Plane method (open circles), are shown in (b) for $K_{S}^{0}$. All data are from $\sqrt{s_{NN}} = 200$ GeV $Au+Au$ collisions.

$v_{2}$ as a function of $p_{T}$ for $10–40\%$ centrality using Lee-Yang Zero method with Sum Generating Function (solid circles), are shown in (a) and (b) for $\Lambda$ and $K_{S}^{0}$ respectively. All data are from $\sqrt{s_{NN}} = 200$ GeV $Au+Au$ collisions.

$v_{2}$ of $K_{S}^{0}$ (open circles) as a function of $p_{T}$ for (a) $0–80\%$, (b) $40–80\%$, (c) $10–40\%$ and (d) $0–10\%$ centralities in $Au+Au$ collisions at $\sqrt{s_{NN}} = 200$ GeV.

$v_{2}$ of $\Lambda+\bar{\Lambda}$ (open squares) as a function of $p_{T}$ for (a) $0–80\%$, (b) $40–80\%$, (c) $10–40\%$ and (d) $0–10\%$ centralities in $Au+Au$ collisions at $\sqrt{s_{NN}} = 200$ GeV.

$v_{2}$ of $\Xi^{-}+\bar{\Xi}^{+}$ (filled triangles) as a function of $p_{T}$ for (a) $0–80\%$ centrality in $Au+Au$ collisions at $\sqrt{s_{NN}} = 200$ GeV.

$v_{2}$ of $\Xi^{-}+\bar{\Xi}^{+}$ (filled triangles) as a function of $p_{T}$ for (b) $40–80\%$ centrality in $Au+Au$ collisions at $\sqrt{s_{NN}} = 200$ GeV.

$v_{2}$ of $\Xi^{-}+\bar{\Xi}^{+}$ (filled triangles) as a function of $p_{T}$ for (c) $10–40\%$ centrality in $Au+Au$ collisions at $\sqrt{s_{NN}} = 200$ GeV.

$v_{2}$ of $\Xi^{-}+\bar{\Xi}^{+}$ (filled triangles) as a function of $p_{T}$ for (d) $0–10\%$ centrality in $Au+Au$ collisions at $\sqrt{s_{NN}} = 200$ GeV.

$v_{2}$ of $\Omega^{-}+\bar{\Omega}^{+}$ (filled circles) as a function of $p_{T}$ for (a) $0–80\%$ centrality in $Au+Au$ collisions at $\sqrt{s_{NN}} = 200$ GeV.

$v_{2}$ of $\Omega^{-}+\bar{\Omega}^{+}$ (filled circles) as a function of $p_{T}$ for (d) $0–10\%$ centrality in $Au+Au$ collisions at $\sqrt{s_{NN}} = 200$ GeV.

$v_{2}$ of $p+\bar{p}$ (filled squares) and $\pi^{+}+\pi^{-}$ (filled stars) as a function of $p_{T}$ for (a) $0–80\%$ centrality in $Au+Au$ collisions at $\sqrt{s_{NN}} = 200$ GeV.

Fig12d inset $\Lambda$ (open squares) expansion at low $p_{T}$

Fig12d inset for $K_{S}^{0}$ (filled circles) expansion at low $p_{T}$

Fig13 Centrality dependence of $v_{2}$ versus number of participants $N_{part}$ for charged hadrons (crosses) in $Au+Au$ collisions at $\sqrt{s_{NN}} = 200$ GeV. For identified particle's value, $v_{2}$ in Tab. II scaled by participant eccentricity in Tab. I versus number of participant in Tab. I.

Fig14 $p_{T}$ dependence of $v_{2}$ of $K_{S}^{0}$ at 200 GeV.

Fig14 $p_{T}$ dependence of $v_{2}$ of $\Lambda$ at 200 GeV.

Delta++ reconstruction efficiency times detector acceptance as a function of the invariant mass for minimum bias d+Au collisions. The error shown is due to the available statistics in the simulation.

Delta++ reconstruction efficiency times detector acceptance as a function of the invariant mass for p+p collisions. The error shown is due to the available statistics in the simulation.

Lambda* reconstruction efficiency times detector acceptance as a function of pT for minimum bias d+Au collisions. The error shown is due to the available statistics in the simulation.

Sigma*+- reconstruction efficiency times detector acceptance as a function of pT for minimum bias d+Au collisions. The error shown is due to the available statistics in the simulation.

Sigma*+- reconstruction efficiency times detector acceptance as a function of pT for d+Au at 40-100% centrality. The error shown is due to the available statistics in the simulation.

Sigma*+- reconstruction efficiency times detector acceptance as a function of pT for d+Au at 20-40% centrality . The error shown is due to the available statistics in the simulation.

Sigma*+- reconstruction efficiency times detector acceptance as a function of pT for d+Au at 0-20% centrality . The error shown is due to the available statistics in the simulation.

Raw pi+pi- invariant mass distribution at midrapidity after subtraction of the like-sign reference distribution for minimum bias d+Au collisions with the hadronic cocktail fit.

rho0 mass as a function of pT at |y| < 0.5 for d+A at 0-20 % centrality.

rho0 mass as a function of pT at |y| < 0.5 for d+A at 20-40 % centrality.

rho0 mass as a function of pT at |y| < 0.5 for d+A at 40-100 % centrality.

rho0 mass as a function of pT at |y| < 0.5 for minimum bias d+A collisions.

rho0 mass as a function of pT at |y| < 0.5 for p+p collisions.

rho0 mass as a function of pT at |y| < 0.5 for high multiplicity p+p collisions.

rho0 mass as a function of pT at |y| < 0.5 for minimum bias d+Au collisions.

Invariant pi+pi- mass distribution after background subtraction for minimum bias p+p collisions measured by NA27.

Mixed-event background subtracted K pi raw invariant mass distribution for a particular K*0 pT bin at |y| < 0.5 for minimum bias d+Au interactions.

Mixed-event background subtracted KS*0 pi raw invariant mass distribution integrated over nvariant mass distribution for a particular K*+- pT bin at |y| < 0.5 for minimum bias d+Au interactions.

K*0 mass as a function of pT at |y| < 0.5 for minimum bias d+Au collisions.

K*+- mass as a function of pT at |y| < 0.5 for minimum bias d+Au collisions.

K*0 width as a function of pT at |y| < 0.5 for minimum bias d+Au collisions.

K*+- width as a function of pT at |y| < 0.5 for minimum bias d+Au collisions.

Mixed-event background subtracted p pi raw invariant mass distribution for a particular pT bin at |y| < 0.5 for minimum bias d+Au interactions.

Mixed-event background subtracted p pi raw invariant mass distribution for a particular pT bin at |y| < 0.5 for minimum bias p+p collisions.

Delta++ mass as a function of pT at |y| < 0.5 for minimum bias d+Au collisions.

Delta++ width as a function of pT at |y| < 0.5 for minimum bias d+Au collisions.

Delta++ mass as a function of pT at |y| < 0.5 for minimum bias p+p collisions.

Delta++ width as a function of pT at |y| < 0.5 for minimum bias p+p collisions.

Mixed-event background subtracted Lambda pi raw invariant mass distribution integrated over the Sigma* pT at |y| < 0.5 for minimum bias d+Au interactions.

Mixed-event background subtracted p K raw invariant mass distribution integrated over the Lambda* pT at |y| < 0.5 for minimum bias d+Au interactions.

Sigma* mass as a function of pT at |y| < 0.75 for minimum bias d+Au collisions.

Lambda* mass as a function of pT at |y| < 0.5 for minimum bias d+Au collisions.

Lambda* width as a function of pT at |y| < 0.5 for minimum bias d+Au collisions.

rho0 invariant yields as a function of pT at |y| < 0.5 for d+Au collisions at 40-100% centrality.

rho0 invariant yields as a function of pT at |y| < 0.5 for d+Au collisions at 20-40% centrality.

rho0 invariant yields as a function of pT at |y| < 0.5 for d+Au collisions at 0-20% centrality.

rho0 invariant yields as a function of pT at |y| < 0.5 for minimum bias d+Au collisions.

(K*0+K*0bar)/2 invariant yields as a function of pT at |y| < 0.5 for minimum bias d+Au collisions.

(K*0+K*0bar)/2 invariant yields as a function of pT at |y| < 0.5 for d+Au collisions at 0-20% centrality.

(K*0 + K*0bar)/2 invariant yields as a function of pT at |y| < 0.5 for d+Au collisions at 20-40% centrality.

(K*0 + K*0bar)/2 invariant yields as a function of pT at |y| < 0.5 for d+Au collisions at 40-100% centrality.

(K*+ + K*-)/2 invariant yields as a function of pT at |y| < 0.5 for minimum bias d+Au collisions.

(K*+ +K*-)/2 invariant yields as a function of pT at |y| < 0.5 for d+Au collisions at 0-20% centrality.

(K*+ +K*-)/2 invariant yields as a function of pT at |y| < 0.5 for d+Au collisions at 20-40% centrality.

(K*+ +K*-)/2 invariant yields as a function of pT at |y| < 0.5 for d+Au collisions at 40-100% centrality.

(Delta++ + Delta--bar)/2 invariant yields as a function of pT at |y| < 0.5 for minimum bias d+Au collisions.

(Delta++ + Delta--bar)/2 invariant yields as a function of pT at |y| < 0.5 for d+Au collisions at 0-20% centrality.

(Delta++ + Delta--bar)/2 invariant yields as a function of pT at |y| < 0.5 for d+Au collisions at 20-40% centrality.

(Delta++ + Delta--bar)/2 invariant yields as a function of pT at |y| < 0.5 for d+Au collisions at 40-100% centrality.

(Delta++ + Delta--bar)/2 invariant yields as a function of pT at |y| < 0.5 for minimum bias p+p collisions. The errors shown also include the systematic uncertainties.

Sigma* invariant yields as a function of pT at |y| < 0.5 for minimum bias d+Au collisions.

Sigma* invariant yields as a function of pT at |y| < 0.5 for d+Au collisions at 0-20% centrality.

Sigma* invariant yields as a function of pT at |y| < 0.5 for d+Au collisions at 20-40% centrality.

Sigma* invariant yields as a function of pT at |y| < 0.5 for d+Au collisions at 40-100% centrality.

Lambda* invariant yields as a function of pT at |y| < 0.5 for minimum bias d+Au collisions.

Rescaled spectra of pion minimum bias d+Au.

Rescaled spectra of p minimum bias d+Au.

Rescaled spectra of rho0 minimum bias d+Au.

Rescaled spectra of K* minimum bias d+Au.

Rescaled spectra of Delta minimum bias d+Au.

Rescaled spectra of Sigma* minimum bias d+Au.

Rescaled spectra of Lambda* minimum bias d+Au.

pion as a function of < dNch > /deta for minimum bias p+p, minimum bias d+Au, and 0-20%, 20-40%, 40-100% of the total d+Au cross section.

proton as a function of < dNch > /deta for minimum bias p+p, minimum bias d+Au, and 0-20%, 20-40%, 40-100% of the total d+Au cross section.

K as a function of < dNch > /deta for minimum bias p+p, minimum bias d+Au, and 0-20%, 20-40%, 40-100% of the total d+Au cross section.

rho0 as a function of < dNch > /deta for minimum bias p+p, minimum bias d+Au, and 0-20%, 20-40%, 40-100% of the total d+Au cross section.

K* as a function of < dNch > /deta for minimum bias p+p, minimum bias d+Au, and 0-20%, 20-40%, 40-100% of the total d+Au cross section.

pion as a function of < dNch > /deta for minimum bias p+p, minimum bias d+Au, and 0-20%, 20-40%, 40-100% of the total d+Au cross section.

proton as a function of < dNch > /deta for minimum bias p+p, minimum bias d+Au, and 0-20%, 20-40%, 40-100% of the total d+Au cross section.

K as a function of < dNch > /deta for minimum bias p+p, minimum bias d+Au, and 0-20%, 20-40%, 40-100% of the total d+Au cross section.

Delta++ as a function of < dNch > /deta for minimum bias p+p, minimum bias d+Au, and 0-20%, 20-40%, 40-100% of the total d+Au cross section.

Sigma* as a function of < dNch > /deta for minimum bias p+p, minimum bias d+Au, and 0-20%, 20-40%, 40-100% of the total d+Au cross section.

Lambda* as a function of < dNch > /deta for minimum bias p+p, minimum bias d+Au, and 0-20%, 20-40%, 40-100% of the total d+Au cross section.

rho0/pi ratios in p+p, various centralities in d+Au, and in peripheral Au+Au collisions as a function of dNch/deta. The errors shown are the quadratic sum of the statistical and systematic uncertainties.

K*/k- ratios in p+p, various centralities in d+Au and Au+Au collisions as a function of dNch/deta. The errors shown are the quadratic sum of the statistical and systematic uncertainties.

Delta++/p ratios in p+p, various centralities in d+Au as a function of dNch/deta. The errors shown are the quadratic sum of the statistical and systematic uncertainties.

Sigma*/Lambda ratios in p+p, various centralities in d+Au, and in central Au+Au collisions as a function of dNch/deta. The errors shown are the quadratic sum of the statistical and systematic uncertainties.

Lambda*/Lambda ratios in p+p, minimum bias d+Au, and two different centralities in Au+Au collisions as a function of dNch/deta. The errors shown are the quadratic sum of the statistical and systematic uncertainties.

Delta++/p ratios normalized by their ratio measured in p+p collisions at the same beam energy as a function of dNch/deta. The errors shown are the quadratic sum of the statistical and systematic uncertainties.

K*/K- ratios normalized by their ratio measured in p+p collisions at the same beam energy as a function of dNch/deta. The errors shown are the quadratic sum of the statistical and systematic uncertainties.

rho0/pi- ratios normalized by their ratio measured in p+p collisions at the same beam energy as a function of dNch/deta. The errors shown are the quadratic sum of the statistical and systematic uncertainties.

Lambda*/Lambda ratios normalized by their ratio measured in p+p collisions at the same beam energy as a function of dNch/deta. The errors shown are the quadratic sum of the statistical and systematic uncertainties.

Sigma*/Lambda ratios normalized by their ratio measured in p+p collisions at the same beam energy as a function of dNch/deta. The errors shown are the quadratic sum of the statistical and systematic uncertainties.

RdAu for pion.

RdAu for charged hadrons.

RdAu for charged rho0.

RdAu for pion.

RdAu for charged hadrons.

RdAu for charged K*0.

RdAu for charged K*+-.

Per-trigger yield versus $\Delta\phi$ for various trigger and partner $p_T$ ($p^a_T \otimes p^b_T$), arranged by increasing pair proxy energy (sum of $p^a_T$ and $p^b_T$), in Au + Au collisions for 5-10 $\otimes$ 2-3, 4-5 $\otimes$ 4-5, 5-10 $\otimes$ 3-5, and 5-10 $\otimes$ 5-10 GeV/c.

RHS versus $p^b_T$ for p + p collisions for four trigger selections.

RHS versus $p^b_T$ for Au + Au collisions for four trigger selections.

RHS versus $N_{part}$ for 2-3 $\otimes$ 2-3 GeV/$c$. Shaded bars (brackets) represent $p_T$ -correlated uncertainties due to elliptic flow (ZYAM procedure). The most left point is from p+p.

Per-trigger yield $\Delta\phi$ distribution and corresponding fits for 2-3 $\otimes$ 2-3 GeV/c in 0-5% Au+Au collisions.

a) D versus centrality from FIT 1 and FIT2 for 2 - 3 GeV/$c$ bin. The error bars are the statistical errors; The shaded bars and brackets are the systematic errors due to $v_2$. b) The fraction of the shoulder Gaussian yield relative to the total away-side yield as function of centrality determined from FIT2.

Values of D determined from FIT1 as function of partner $p_T$ for 2-3 GeV/c $p_T$ range in 0-20% Au+Au.

Values of D determined from FIT1 as function of partner $p_T$ for 3-4 GeV/c $p_T$ range in 0-20% Au+Au.

Values of D determined from FIT1 as function of partner $p_T$ for 4-5 GeV/c $p_T$ range in 0-20% Au+Au.

Values of D determined from FIT2 as function of partner $p_T$ for 2-3 GeV/c $p_T$ range in 0-20% Au+Au.

Values of D determined from FIT2 as function of partner $p_T$ for 3-4 GeV/c $p_T$ range in 0-20% Au+Au.

Values of D determined from FIT2 as function of partner $p_T$ for 4-5 GeV/c $p_T$ range in 0-20% Au+Au.

Per-trigger yield as function of partner $p_T$ for the HR region for the 0-20% Au + Au collisions. Results for four trigger pT in 2-3, 3-4, 4-5, 5-10 GeV/c are shown.

Per-trigger yield as function of partner $p_T$ for the HR region for the 20-40% Au + Au collisions. Results for four trigger pT in 2-3, 3-4, 4-5, 5-10 GeV/c are shown.

Per-trigger yield as function of partner $p_T$ for the HR region for the 40-60% Au + Au collisions. Results for four trigger pT in 2-3, 3-4, 4-5, 5-10 GeV/c are shown.

Per-trigger yield as function of partner $p_T$ for the HR region for the 60-92% Au + Au collisions. Results for four trigger pT in 2-3, 3-4, 4-5, 5-10 GeV/c are shown.

Per-trigger yield as function of partner $p_T$ for the HR region for the p + p collisions. Results for four trigger pT in 2-3, 3-4, 4-5, 5-10 GeV/c are shown.

Per-trigger yield as function of partner $p_T$ for the SR region for the 0-20% Au + Au collisions. Results for four trigger pT in 2-3, 3-4, 4-5, 5-10 GeV/c are shown.

Per-trigger yield as function of partner $p_T$ for the SR region for the 20-40% Au + Au collisions. Results for four trigger pT in 2-3, 3-4, 4-5, 5-10 GeV/c are shown.

Per-trigger yield as function of partner $p_T$ for the SR region for the 40-60% Au + Au collisions. Results for four trigger pT in 2-3, 3-4, 4-5, 5-10 GeV/c are shown.

Per-trigger yield as function of partner $p_T$ for the SR region for the 60-92% Au + Au collisions. Results for four trigger pT in 2-3, 3-4, 4-5, 5-10 GeV/c are shown.

Per-trigger yield as function of partner $p_T$ for the SR region for the p + p collisions. Results for four trigger pT in 2-3, 3-4, 4-5, 5-10 GeV/c are shown.

Near-side Gaussian widths versus partner $p_T$ for four trigger $p_T$ ranges for 0-20 % Au + Au.

Near-side Gaussian widths versus partner $p_T$ for four trigger $p_T$ ranges for 0-20 % p + p.

Near-side Gaussian widths versus $N_{part}$ for five successively increasing $p^a_T$ $\otimes$ $p^b_T$. The most left point in each panel represents the value from p + p.

Near-side yield in $\lvert \Delta\phi \rvert$ < $\pi$/3 versus partner $p_T$ for four trigger $p_T$ selections for 0-20% p + p.

Near-side yield in $\lvert \Delta\phi \rvert$ < $\pi$/3 versus partner $p_T$ for four trigger $p_T$ selections for 0-20% Au + Au.

Near-side yield in $\lvert \Delta\phi \rvert$ < $\pi$/3 versus partner $p_T$ for four trigger $p_T$ selections for 0-20% Au + Au.

Near-side yield in $\lvert \Delta\phi \rvert$ < $\pi$/3 versus partner $p_T$ for four trigger $p_T$ selections for 0-20% Au + Au.

Near-side yield in $\lvert \Delta\phi \rvert$ < $\pi$/3 versus partner $p_T$ for four trigger $p_T$ selections for 0-20% Au + Au.

Per-trigger yield versus $\Delta\eta$ for 0-20% central Au + Au collisions. Results are shown for four $p^a_T \otimes p^b_T$ selections as indicated.

Per-trigger yield versus $\Delta\eta$ for 0-20% central p + p collisions. Results are shown for four $p^a_T \otimes p^b_T$ selections as indicated.

The projection of 2-D per-trigger yield in $\lvert \Delta\eta \rvert$ < 0.7 $\otimes$ $\lvert \Delta\phi \rvert$ < 0.7 and 0-20% central Au+Au collisions onto $\Delta\eta$ for four $p^a_T \otimes p^b_T$ selections.

The projection of 2-D per-trigger yield in $\lvert \Delta\eta \rvert$ < 0.7 $\otimes$ $\lvert \Delta\phi \rvert$ < 0.7 and 0-20% central Au+Au collisions onto $\Delta\phi$ for four $p^a_T \otimes p^b_T$ selections.

Truncated mean $p_T$, $\langle p_T \rangle$, in 1 <pbT < 5 GeV/c versus centrality for the head region for Au + Au for four trigger $p_T$ bins.

Truncated mean $p_T$, $\langle p_T \rangle$, in 1 <pbT < 5 GeV/c versus centrality for the away-side shoulder region for Au + Au for four trigger $p_T$ bins.

Truncated mean $p_T$, $\langle p_T \rangle$, in 1 <pbT < 5 GeV/c versus centrality for the near-side region for Au + Au for four trigger $p_T$ bins.

Truncated mean $p_T$, $\langle p_T \rangle$, in 1 <pbT < 5 GeV/c versus centrality for the head region for p + p for four trigger $p_T$ bins.

Truncated mean $p_T$, $\langle p_T \rangle$, in 1 <pbT < 5 GeV/c versus centrality for the away-side shoulder region for p + p for four trigger $p_T$ bins.

Truncated mean $p_T$, $\langle p_T \rangle$, in 1 <pbT < 5 GeV/c versus centrality for the near-side region for p + p for four trigger $p_T$ bins.

Truncated mean $p_T$, $\langle p_T \rangle$, calculated for five partner pT ranges in the HR. Results for triggers in 3-4 GeV/c.

Truncated mean $p_T$, $\langle p_T \rangle$, calculated for five partner pT ranges in the HR. Results for triggers in 4-5 GeV/c.

Truncated mean $p_T$, $\langle p_T \rangle$, calculated for five partner pT ranges in the HR. Results for triggers in 3-4 GeV/c.

Truncated mean $p_T$, $\langle p_T \rangle$, calculated for five partner pT ranges in the HR. Results for triggers in 4-5 GeV/c.

The near-side pair suppression factor $J_{AA}$ in 0-20% Au + Au collisions as function of $p^b_T$ for various ranges of $p^a_T$.

The near-side $J_{AA}$ as function of $p^b_T$ and $p^{sum}_T$ for four $p^a_T$ bins and three centrality bins. Data for 20-40%.

The near-side $J_{AA}$ as function of $p^b_T$ and $p^{sum}_T$ for four $p^a_T$ bins and three centrality bins. Data for 40-60%.

The near-side $J_{AA}$ as function of $p^b_T$ and $p^{sum}_T$ for four $p^a_T$ bins and three centrality bins. Data for 60-92%.

The pair suppression factor $J_{AA}$ for the away-side HR in 0-20% Au + Au collisions as function of $p_b^T$ for various ranges of $p_a^T$.

The pair suppression factor $J_{AA}$ for the away-side HR in 20-40% Au + Au collisions as function of $p_b^T$ for various ranges of $p_a^T$.

The pair suppression factor $J_{AA}$ for the away-side HR in 40-60% Au + Au collisions as function of $p_b^T$ for various ranges of $p_a^T$.

The pair suppression factor $J_{AA}$ for the away-side HR in 60-92% Au + Au collisions as function of $p_b^T$ for various ranges of $p_a^T$.

Per-trigger yield versus $\Delta\phi$ for successively increasing trigger and partner $p_T$ ($p^a_T \otimes p^b_T$) in 0-20 % Au + Au collisions.

Per-trigger yield versus $\Delta\phi$ for successively increasing trigger and partner $p_T$ ($p^a_T \otimes p^b_T$) in 20-40 % Au + Au collisions.

Per-trigger yield versus $\Delta\phi$ for successively increasing trigger and partner $p_T$ ($p^a_T \otimes p^b_T$) in 40-60 % Au + Au collisions.

Per-trigger yield versus $\Delta\phi$ for successively increasing trigger and partner $p_T$ ($p^a_T \otimes p^b_T$) in 60-92 % Au + Au collisions.

Per-trigger yield versus $\Delta\phi$ for successively increasing trigger and partner $p_T$ ($p^a_T \otimes p^b_T$) in p + p collisions.

$\pi^0$ invariant yields for different centralities. The bin range is not an uncertainty in the x-axis because the actual uncertainty by having the finite bin width is corrected for by the bin-shift correction. These bins were constructed using the corrected finite values as centers.

$\pi^0$ invariant yields for different centralities. The bin range is not an uncertainty in the x-axis because the actual uncertainty by having the finite bin width is corrected for by the bin-shift correction. These bins were constructed using the corrected finite values as centers.

Nuclear modification factor ($R_{AA}$) for $\pi^0$s. The bin range is not an uncertainty in the x-axis because the actual uncertainty by having the finite bin width is corrected for by the bin-shift correction. These bins were constructed using the corrected finite values as centers.

Nuclear modification factor ($R_{AA}$) for $\pi^0$s. The bin range is not an uncertainty in the x-axis because the actual uncertainty by having the finite bin width is corrected for by the bin-shift correction. These bins were constructed using the corrected finite values as centers.

Nuclear modification factor ($R_{AA}$) for $\pi^0$s. The bin range is not an uncertainty in the x-axis because the actual uncertainty by having the finite bin width is corrected for by the bin-shift correction. These bins were constructed using the corrected finite values as centers.

Nuclear modification factor ($R_{AA}$) for $\pi^0$s. The bin range is not an uncertainty in the x-axis because the actual uncertainty by having the finite bin width is corrected for by the bin-shift correction. These bins were constructed using the corrected finite values as centers.

Nuclear modification factor ($R_{AA}$) for $\pi^0$s. The bin range is not an uncertainty in the x-axis because the actual uncertainty by having the finite bin width is corrected for by the bin-shift correction. These bins were constructed using the corrected finite values as centers.

$\pi^0$ $R_{AA}$ for the most central (0-5%) Au+Au collisions and PQM model calculations for different values of $<\bar{q}>$. The bin range is not an uncertainty in the x-axis because the actual uncertainty by having the finite bin width is corrected for by the bin-shift correction. These bins were constructed using the corrected finite values as centers.

Integrated nuclear modification factor ($R_{AA}$) for $\pi^0$ as a function of collision centrality.

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