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

Experimental correlation moments $R^2_{x2}(q)$ Data. Systematic errors are less than the statistical errors.

Experimental correlation moments $R^2_{x2}(q)$ Fit. Systematic errors are less than the statistical errors.

Experimental correlation moments $R^2_{y2}(q)$ Data. Systematic errors are less than the statistical errors.

Experimental correlation moments $R^2_{y2}(q)$ Fit. Systematic errors are less than the statistical errors.

Experimental correlation moments $R^6_{x6}(q)$ Data. Systematic errors are less than the statistical errors.

Experimental correlation moments $R^6_{x6}(q)$ Fit. Systematic errors are less than the statistical errors.

Experimental correlation moments $R^6_{y6}(q)$ Data. Systematic errors are less than the statistical errors.

Experimental correlation moments $R^6_{y6}(q)$ Fit. Systematic errors are less than the statistical errors.

Experimental correlation moments $R^4_{x2y2}(q)$ Data. Systematic errors are less than the statistical errors.

Experimental correlation moments $R^4_{x2y2}(q)$ Fit. Systematic errors are less than the statistical errors.

Experimental correlation moments $R^4_{x4}(q)$ Data. Systematic errors are less than the statistical errors.

Experimental correlation moments $R^4_{x}(q)$ Fit. Systematic errors are less than the statistical errors.

Experimental correlation moments $R^4_{y4}(q)$ Data. Systematic errors are less than the statistical errors.

Experimental correlation moments $R^4_{y4}(q)$ Fit. Systematic errors are less than the statistical errors.

Experimental correlation moments $R^6_{x2y4}(q)$ Data. Systematic errors are less than the statistical errors.

Experimental correlation moments $R^6_{x2y4}(q)$ Fit. Systematic errors are less than the statistical errors.

Experimental correlation moments $R^6_{x4y2}(q)$ Data. Systematic errors are less than the statistical errors.

Experimental correlation moments $R^6_{x4y2}(q)$ Fit. Systematic errors are less than the statistical errors.

Source function profile $S(r_x)$ Image in the PCMS.

Source function profile $S(r_x)$ Fit in the PCMS. The blue line represents the Hump fit.

Source function profile $S(r_y)$ Image in the PCMS.

Source function profile $S(r_y)$ Fit in the PCMS. The blue line represents the Hump fit.

Source function profile $S(r_z)$ Image in the PCMS.

Source function profile $S(r_z)$ Fit in the PCMS. The blue line represents the Hump fit.

Associated correlation profile $C(q_x)$ Image in the PCMS.

Associated correlation profile $C(q_x)$ Fit in the PCMS.

Associated correlation profile $C(q_y)$ Image in the PCMS.

Associated correlation profile $C(q_y)$ Fit in the PCMS.

Associated correlation profile $C(q_z)$ Image in the PCMS.

Associated correlation profile $C(q_z)$ Fit in the PCMS.

Source function comparison between Therminator calculation and image for $S(r_x)$ in PCMS.

Source function comparison between Therminator calculation and image for $S(r_x)$ in PCMS.

Source function comparison between Therminator calculation and image for $S(r_y)$ in PCMS.

Source function comparison between Therminator calculation and image for $S(r_y)$ in PCMS.

Source function comparison between Therminator calculation and image for $S(r_z)$ in PCMS.

Source function comparison between Therminator calculation and image for $S(r_z)$ in PCMS.

$\Delta$$t_{LCM}$ from Therminator events compared with various assumptions for $\Delta\tau$ and resonance emissions.

$\Delta$$t_{LCM}$ from Therminator events compared with various assumptions for $\Delta\tau$ and resonance emissions.

J/PSI invariant (1/(2PI*PT))*D2(N)/DPT/DYRAP versus PT at mid rapidity (-0.35<y<0.35) in D+AU collisions.

J/PSI invariant (1/(2PI*PT))*D2(N)/DPT/DYRAP versus PT at forward rapidity (1.2<y<2.2) in D+AU collisions.

J/PSI invariant (1/(2PI*PT))*D2(N)/DPT/DYRAP versus centrality at backward rapidity (-2.2<y<-1.2) in D+AU collisions.

J/PSI invariant (1/(2PI*PT))*D2(N)/DPT/DYRAP versus centrality at mid rapidity (-0.35<y<0.35) in D+AU collisions.

J/PSI invariant (1/(2PI*PT))*D2(N)/DPT/DYRAP versus centrality at forward rapidity (1.2<y<2.2) in D+AU collisions.

Nuclear modification factor versus rapidity in D+AU collisions, over 3 bins of rapidity.

Nuclear modification factor versus rapidity in D+AU collisions, over 5 bins of rapidity.

Nuclear modification factor versus PT at backward rapidity (-2.2<y<-1.2) in D+AU collisions.

Nuclear modification factor versus PT at mid rapidity (-0.35<y<0.35) in D+AU collisions.

Nuclear modification factor versus PT at forward rapidity (1.2<y<2.2) in D+AU collisions.

Nuclear modification factor versus centrality at backward rapidity (-2.2<y<-1.2) in D+AU collisions.

Nuclear modification factor versus centrality at mid rapidity (-0.35<y<0.35) in D+AU collisions.

Nuclear modification factor versus centrality at forward rapidity (1.2<y<2.2) in D+AU collisions.

Mean PT^2 versus rapidity, extracted from a fit to the data, for DEUT + Au reactions. The mean value of PT^2 is calculated from the data for PT < 5 Gev/c.

Mean PT^2 versus rapidity, extracted from a fit to the data, for P + P reactions. The mean value of PT^2 is calculated from the data for PT < 5 Gev/c.

Breakup cross section of c-c_bar pairs inside cold nuclear matter,integrated over the whole domain of rapidity.The breakup cross section is calculated with two models of shadowing for nuclear PDFs ; the EKS model and the NDSG model. The uncertainties given, containing statistical and systematical error, are corresponding to one standard deviation.

Breakup cross section of c-c_bar pairs inside cold nuclear matter for different ranges of rapidity.The breakup cross section is calculated with two models of shadowing for nuclear PDFs ; the EKS model and the NDSG model. The uncertainties given, containing statistical and systematical error, are corresponding to one standard deviation.

Breakup cross section of c-c_bar pairs inside cold nuclear matter, integrated over the whole domain of rapidity.The breakup cross section is calculated with two models of shadowing for nuclear PDFs ; the EKS model and the NDSG model. The uncertainties given, containing statistical and systematical error, are corresponding to one standard deviation.

Breakup cross section of c-c_bar pairs inside cold nuclear matter for different ranges of rapidity.The breakup cross section is calculated with two models of shadowing for nuclear PDFs ; the EKS model and the NDSG model. The uncertainties given, containing statistical and systematical error, are corresponding to one standard deviation.

The double transverse spin asymmetry ($A_{TT}$) for neutral pion production at $\sqrt{s}$ = 200 GeV as a function of $p_T$ (GeV/$c$). The correlated for all points scale systematic uncertainty is 9.4% (scales both the values and stat. uncertainties by the same factor).

Mean PT^2 value at forward rapidities : absolute value of y belongs to [1.2;2.2] The mean PT is obtained with a phenomonological fit of the J/PSI distribution in PT of the form (1/(2*PI*PT))*D(SIG)/DPT = A ( 1+(PT/B)^2)^-6 .The systematic error includes the incertainty from the maximum shape deviation permitted by the point-to-point correlated errors and from allowing the exponent of the fit fonctionto be a free parameter.

J/PSI differential cross section times dilepton branching ratio versus rapidity.

Total cross-section for J/PSI production. The value of the total cross section has been determined by fitting the rapidity distribution with many theoretical and phenomenological shapes. The systematic error quotedwas estimated from the maximum variation allowed by shifting the mid and forward rapidity data independantely by their point-to-point correlated systematic errors.There is a +- 18 nanobarn normalization error in addition to the statistical and systematical errors shown.

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