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*+-.

J/psi-->e+e- invariant yield in Cu+Cu collisions as a function of p_T at mid-rapidity for the 60-94 centrality range. The statistical and systematic uncertainties vary point-to-point and are listed for each measured value. An additional global systematic uncertainty is provided in each column heading, which applies to all data points per column.

J/psi-->mu+mu- invariant yield in Cu+Cu collisions as a function of p_T at forward rapidity for the 0-20, 20-40, 40-60 and 60-94 centrality ranges. The statistical and systematic uncertainties vary point-to-point and are listed for each measured value. An additional global systematic uncertainty is provided in each column heading, which applies to all data points per column.

J/psi-->e+e- average momentum squared in Cu+Cu collisions as a function of N_part at forward rapidity (p_T integrated). The statistical and systematic uncertainties vary point-to-point and are listed for each measured value.

J/psi-->mu+mu- average momentum squared in Cu+Cu collisions as a function of N_part at forward rapidity (p_T integrated). The statistical and systematic uncertainties vary point-to-point and are listed for each measured value.

J/psi-->mu+mu- and J/psi-->e+e- nuclear modification in Cu+Cu collisions as a function of p_T at forward and mid-rapidity for the 0-20 centrality range. The statistical and systematic uncertainties vary point-to-point and are listed for each measured value. An additional global systematic uncertainty is provided in each column heading, which applies to all data points per column.

J/psi-->mu+mu- and J/psi-->e+e- nuclear modification in Cu+Cu collisions as a function of y at forward and mid-rapidity for the 0-20 centrality range. The statistical and systematic uncertainties vary point-to-point and are listed for each measured value. An additional global systematic uncertainty is provided in each column heading, which applies to all data points per column.

J/psi-->e+e- nuclear modification in Cu+Cu collisions as a function of N_part at mid-rapidity for the 0-20 centrality range. The statistical and systematic uncertainties vary point-to-point and are listed for each measured value. An additional global systematic uncertainty is provided in each column heading, which applies to all data points per column.

J/psi-->mu+mu- nuclear modification in Cu+Cu collisions as a function of N_part at forward rapidity for the 0-20 centrality range. The statistical and systematic uncertainties vary point-to-point and are listed for each measured value. An additional global systematic uncertainty is provided in each column heading, which applies to all data points per column.

J/psi-->mu+mu- in Cu+Cu collisions as a function of $N_{part}$ at forward rapidity expressed as a ratio with respect to J/psi-->e+e- nuclear modification at mid-rapidity. The statistical and systematic uncertainties vary point-to-point and are listed for each measured value. An additional global systematic uncertainty is provided in each column heading, which applies to all data points per column.

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

The fraction of the direct photon component as a function of $p_T$.

The fraction of the direct photon component as a function of $p_T$.

Invariant yield (Au+Au) of direct photons as a function of $p_T$.

Invariant yield (Au+Au) of direct photons as a function of $p_T$.

Invariant yield (Au+Au) of direct photons as a function of $p_T$.

Invariant cross section ($p$+$p$) of direct photons as a function of $p_T$.

FIG. 1: (a) Raw two-particle correlation signal $Y_2$ (red), background $aB_{inc}F_2$ (solid histogram), and background systematic uncertainty from a (dashed histograms). (b) Background-subtracted two-particle correlation $\hat{Y}_2$ (red), and systematic uncertainties due to a (dashed histograms) and flow (blue histograms). (c) Raw three-particle correlation $Y_3$. (d) $ba^2Y_{inc}^2$ . (e) Sum of trig-corr-bkgd and trigger flow. Data are from 12% central Au+Au collisions. Statistical errors in (a,b) are smaller than the point size.

FIG. 1: (a) Raw two-particle correlation signal $Y_2$ (red), background $aB_{inc}F_2$ (solid histogram), and background systematic uncertainty from a (dashed histograms). (b) Background-subtracted two-particle correlation $\hat{Y}_2$ (red), and systematic uncertainties due to a (dashed histograms) and flow (blue histograms). (c) Raw three-particle correlation $Y_3$. (d) $ba^2Y_{inc}^2$ . (e) Sum of trig-corr-bkgd and trigger flow. Data are from 12% central Au+Au collisions. Statistical errors in (a,b) are smaller than the point size.

FIG. 2: Background subtracted three-particle correlations, $\hat{Y}_3$, for (a) pp, (b) d+Au, (c) 80-50% Au+Au, (d) 50-30% Au+Au, (e) 30-10% Au+Au, and (f) 12% central Au+Au. Statistical errors per bin are approximately $\pm$0.012 in (a) and $\pm$0.006 in (b), both at ($\pi$, $\pi$), and are $\pm$0.022, $\pm$0.049, $\pm$0.099 and $\pm$0.077 from (c) to (f), similar for all bins.

FIG. 2: Background subtracted three-particle correlations, $\hat{Y}_3$, for (a) pp, (b) d+Au, (c) 80-50% Au+Au, (d) 50-30% Au+Au, (e) 30-10% Au+Au, and (f) 12% central Au+Au. Statistical errors per bin are approximately $\pm$0.012 in (a) and $\pm$0.006 in (b), both at ($\pi$, $\pi$), and are $\pm$0.022, $\pm$0.049, $\pm$0.099 and $\pm$0.077 from (c) to (f), similar for all bins.

FIG. 2: Background subtracted three-particle correlations, $\hat{Y}_3$, for (a) pp, (b) d+Au, (c) 80-50% Au+Au, (d) 50-30% Au+Au, (e) 30-10% Au+Au, and (f) 12% central Au+Au. Statistical errors per bin are approximately $\pm$0.012 in (a) and $\pm$0.006 in (b), both at ($\pi$, $\pi$), and are $\pm$0.022, $\pm$0.049, $\pm$0.099 and $\pm$0.077 from (c) to (f), similar for all bins.

FIG. 2: Background subtracted three-particle correlations, $\hat{Y}_3$, for (a) pp, (b) d+Au, (c) 80-50% Au+Au, (d) 50-30% Au+Au, (e) 30-10% Au+Au, and (f) 12% central Au+Au. Statistical errors per bin are approximately $\pm$0.012 in (a) and $\pm$0.006 in (b), both at ($\pi$, $\pi$), and are $\pm$0.022, $\pm$0.049, $\pm$0.099 and $\pm$0.077 from (c) to (f), similar for all bins.

FIG. 2: Background subtracted three-particle correlations, $\hat{Y}_3$, for (a) pp, (b) d+Au, (c) 80-50% Au+Au, (d) 50-30% Au+Au, (e) 30-10% Au+Au, and (f) 12% central Au+Au. Statistical errors per bin are approximately $\pm$0.012 in (a) and $\pm$0.006 in (b), both at ($\pi$, $\pi$), and are $\pm$0.022, $\pm$0.049, $\pm$0.099 and $\pm$0.077 from (c) to (f), similar for all bins.

FIG. 2: Background subtracted three-particle correlations, $\hat{Y}_3$, for (a) pp, (b) d+Au, (c) 80-50% Au+Au, (d) 50-30% Au+Au, (e) 30-10% Au+Au, and (f) 12% central Au+Au. Statistical errors per bin are approximately $\pm$0.012 in (a) and $\pm$0.006 in (b), both at ($\pi$, $\pi$), and are $\pm$0.022, $\pm$0.049, $\pm$0.099 and $\pm$0.077 from (c) to (f), similar for all bins.

Projections of away-side three-particle correlations along the diagonal $\Sigma$ within 0 < $\Delta$ < 0.35 (squares) and along the off-diagonal $\Delta$ within |$\Sigma$| < 0.35 (points) in (a) d+Au and (b) 12% central Au+Au collisions. The shaded areas indicate systematic uncertainties on the off-diagonal projections. The histogram in (a) is the near-side off-diagonal projection. The histogram in (b) is the away-side off-diagonal projection of our result with a = b = 1. NOTE: For points with invisible error bars, the point size was considered as an absolute upper limit for the uncertainty.

Projections of away-side three-particle correlations along the diagonal $\Sigma$ within 0 < $\Delta$ < 0.35 (squares) and along the off-diagonal $\Delta$ within |$\Sigma$| < 0.35 (points) in (a) d+Au and (b) 12% central Au+Au collisions. The shaded areas indicate systematic uncertainties on the off-diagonal projections. The histogram in (a) is the near-side off-diagonal projection. The histogram in (b) is the away-side off-diagonal projection of our result with a = b = 1. NOTE: For points with invisible error bars, the point size was considered as an absolute upper limit for the uncertainty.

Three-particle correlation strength per radian$^2$ versus ($N_{part}$/2)$^{1/3}$ where $N_{part}$ is the number of participants. Some of the data points have been displaced in $N_{part}$ for clarity.