FIG. 2: a) left side: The correlation data for the ratio of the histograms of same-event-pairs to mixed-event-pairs for unlike-sign charged pairs, shown in a two-dimensional (2-D) perspective plot $\Delta\phi$ - $\Delta\eta$. The plot was normalized to a mean of 1. b) right side: The similar correlation data for like-sign charge pairs.

FIG. 3: a) left side: Folded correlation data for unlike-sign charge pairs on $\Delta\phi$ and $\Delta\eta$ based on demonstrated symmetries in the data. This increases the statistics in a typical bin by a factor of four. Henceforth we will be dealing with folded data only. b) right side: Similar folded data for like-sign charge pairs.

FIG. 3: a) left side: Folded correlation data for unlike-sign charge pairs on $\Delta\phi$ and $\Delta\eta$ based on demonstrated symmetries in the data. This increases the statistics in a typical bin by a factor of four. Henceforth we will be dealing with folded data only. b) right side: Similar folded data for like-sign charge pairs.

FIG. 4: “(Color online)”a) left side: The 2 dimensional (2-D) approximately gaussian signal shape equation (4) from the fit to the normalized and folded unlike-sign charge pairs signal data. b) right side: The unlike-sign charge pairs signal data corresponding to the adjoining signal function on the left side.

FIG. 4: “(Color online)”a) left side: The 2 dimensional (2-D) approximately gaussian signal shape equation (4) from the fit to the normalized and folded unlike-sign charge pairs signal data. b) right side: The unlike-sign charge pairs signal data corresponding to the adjoining signal function on the left side.

FIG. 5 :“(Color online)”a) left side: A perspective plot of the fit to the 2-D like-sign charge pairs signal shape of equation (5). b) right side: Like-sign charge pairs signal data corresponding to the adjacent signal fit on the left side.

FIG. 5: “(Color online)”a) left side: A perspective plot of the fit to the 2-D like-sign charge pairs signal shape of equation (5). b) right side: Like-sign charge pairs signal data corresponding to the adjacent signal fit on the left side.

FIG. 6: “(Color online)”a) left side: The charge dependent (CD) signal shape is a 2 dimensional (2-D) approximate gaussian which is symmetrical. The average gaussian width is 27.5$^{\circ}$ $\pm$ 3.2$^{\circ}$ (0.48 $\pm$ 0.056 rad) both in $\Delta\phi$ and $\Delta\eta$. Thus we are observing a greater probability for unlike-sign charge pairs of particles than for like-sign charge pairs emitted randomly on 2-D $\eta\phi$ directions. b) right side: The CD signal data corresponding to the adjacent signal fit on the left.

FIG. 6: “(Color online)”a) left side: The charge dependent (CD) signal shape is a 2 dimensional (2-D) approximate gaussian which is symmetrical. The average gaussian width is 27.5$^{\circ}$ $\pm$ 3.2$^{\circ}$ (0.48 $\pm$ 0.056 rad) both in $\Delta\phi$ and $\Delta\eta$. Thus we are observing a greater probability for unlike-sign charge pairs of particles than for like-sign charge pairs emitted randomly on 2-D $\eta\phi$ directions. b) right side: The CD signal data corresponding to the adjacent signal fit on the left.

FIG. 7: “(Color online)”a) left side: The fit to the charge independent (CI) signal shape which is the sum of the fits of like-sign plus unlike-sign charge pairs signals. The 2-D gaussian equivalent RMS signal has a $\Delta\phi$ width of about 32$^{\circ}$ (0.52 rad), and $\Delta\eta$ width of about 66$^{\circ}$ (1.15 rad) which is about double the $\Delta\phi$ width. b) right side: The CI signal data corresponding to the adjacent fit on the left.

FIG. 7: “(Color online)”a) left side: The fit to the charge independent (CI) signal shape which is the sum of the fits of like-sign plus unlike-sign charge pairs signals. The 2-D gaussian equivalent RMS signal has a $\Delta\phi$ width of about 32$^{\circ}$ (0.52 rad), and $\Delta\eta$ width of about 66$^{\circ}$ (1.15 rad) which is about double the $\Delta\phi$ width. b) right side: The CI signal data corresponding to the adjacent fit on the left.

FIG. 8: a) left side: The charge independent (CI) signal data with the lower range of elliptic flow amplitude and the best $\chi^2$ (the same as Figure 7b) plotted as a 2-D perspective plot on $\Delta\phi$ vs. $\Delta\eta$. b) right side: The charge independent (CI) signal data with our maximum elliptic flow amplitude used in our analysis ($\chi^2$ was worse by 1$\sigma$).

FIG. 8: a) left side: The charge independent (CI) signal data with the lower range of elliptic flow amplitude and the best $\chi^2$ (the same as Figure 7b) plotted as a 2-D perspective plot on $\Delta\phi$ vs. $\Delta\eta$. b) right side: The charge independent (CI) signal data with our maximum elliptic flow amplitude used in our analysis ($\chi^2$ was worse by 1$\sigma$).

Non photonic electrons yield versus transverse momentum in D+AU collisions.

Non photonic electrons yield versus transverse momentum in P+P collisions.

D0 differential yield at mid rapidity and corresponding calculated differential charm cross section at midrapidity.

Sum of the non photonic electrons/positrons and D0 differential yield at mid rapidity and corresponding calculated differential charm cross section at midrapidity.

Total charm cross section per nucleon nucleon collision IN D+AU collisions.

xT scaling (with xT.

J/psi nuclear modification factor RAA from sqrt(s).

J/psi nuclear modification factor RAA from sqrt(s).

J/psi-hadron azimuthal correlations in sqrt(s).

Contribution from B-hadron feed-down to inclusive J/psi production in sqrt(s).

The jet-like yield in $\mid\Delta\eta\mid<$0.78 as a function of $p_T^{trigger}$ for $K_S^0$-h and $\Lambda$-h correlations for 1.5 GeV/$c$ < $p_T^{associated}$ < $p_T^{trigger}$ in minimum bias $d$+Au collisions at $\sqrt{s_{NN}}$ = 200 GeV. Statistical errors listed in table. 5.4% systematic error on all points.

The jet-like yield in $\mid\Delta\eta\mid<$0.78 as a function of $p_T^{trigger}$ for $K_S^0$-h and $\Lambda$-h correlations for 1.5 GeV/$c$ < $p_T^{associated}$ < $p_T^{trigger}$ in 0-60% Cu+Cu collisions at $\sqrt{s_{NN}}$ = 200 GeV. Statistical errors listed in table. 5.4% systematic error on all points.

The jet-like yield in $\mid\Delta\eta\mid<$0.78 as a function of $p_T^{trigger}$ for $K_S^0$-h and $\Lambda$-h correlations for 1.5 GeV/$c$ < $p_T^{associated}$ < $p_T^{trigger}$ in 40-80% Au+Au collisions at $\sqrt{s_{NN}}$ = 200 GeV. Statistical errors listed in table. 5.4% systematic error on all points.

The jet-like yield in $\mid\Delta\eta\mid<$0.78 as a function of $p_T^{trigger}$ for $K_S^0$-h and $\Lambda$-h correlations for 1.5 GeV/$c$ < $p_T^{associated}$ < $p_T^{trigger}$ in 0-12% Au+Au collisions at $\sqrt{s_{NN}}$ = 200 GeV. Statistical errors listed in table. 5.4% systematic error on all points.

Centrality dependence of the jet-like yield of $K_S^0$-h and $\Lambda$-h correlations for 3 < $p_T^{trigger}$ < 6 GeV/$c$ and 1.5 GeV/$c$ < $p_T^{associated}$ < $p_T^{trigger}$ in Cu+Cu collisions at $\sqrt{s_{NN}}$ = 200 GeV. 5.4% systematic error on all points.

Centrality dependence of the jet-like yield of $K_S^0$-h and $\Lambda$-h correlations for 3 < $p_T^{trigger}$ < 6 GeV/$c$ and 1.5 GeV/$c$ < $p_T^{associated}$ < $p_T^{trigger}$ in Au+Au collisions at $\sqrt{s_{NN}}$ = 200 GeV. 5.4% systematic error on all points.

Centrality dependence of the jet-like yield of $K_S^0$-h and $\Lambda$-h correlations for 3 < $p_T^{trigger}$ < 6 GeV/$c$ and 1.5 GeV/$c$ < $p_T^{associated}$ < $p_T^{trigger}$ in $d$+Au collisions at $\sqrt{s_{NN}}$ = 200 GeV. 5.4% systematic error on all points.

Centrality dependence of the jet-like yield of $K_S^0$-h and $\Lambda$-h correlations for 3 < $p_T^{trigger}$ < 6 GeV/$c$ and 1.5 GeV/$c$ < $p_T^{associated}$ < $p_T^{trigger}$ in minimum bias $d$+Au and 0-60% central Cu+Cu collisions at $\sqrt{s_{NN}}$ = 200 GeV. 5.4% systematic error on all points.

Centrality dependence of moments of $\Delta N_p$ distributions for Au+Au collisions at $\sqrt{s_{NN}}$ = 19.6 GeV.

Centrality dependence of moments of $\Delta N_p$ distributions for Au+Au collisions at $\sqrt{s_{NN}}$ = 19.6 GeV.

Centrality dependence of moments of $\Delta N_p$ distributions for Au+Au collisions at $\sqrt{s_{NN}}$ = 19.6 GeV.

Centrality dependence of moments of $\Delta N_p$ distributions for Au+Au collisions at $\sqrt{s_{NN}}$ = 19.6 GeV.

Centrality dependence of moments of $\Delta N_p$ distributions for Au+Au collisions at $\sqrt{s_{NN}}$ = 62.4 GeV.

Centrality dependence of moments of $\Delta N_p$ distributions for Au+Au collisions at $\sqrt{s_{NN}}$ = 62.4 GeV.

Centrality dependence of moments of $\Delta N_p$ distributions for Au+Au collisions at $\sqrt{s_{NN}}$ = 62.4 GeV.

Centrality dependence of moments of $\Delta N_p$ distributions for Au+Au collisions at $\sqrt{s_{NN}}$ = 62.4 GeV.

Centrality dependence of moments of $\Delta N_p$ distributions for Au+Au collisions at $\sqrt{s_{NN}}$ = 200 GeV.

Centrality dependence of moments of $\Delta N_p$ distributions for Au+Au collisions at $\sqrt{s_{NN}}$ = 200 GeV.

Centrality dependence of moments of $\Delta N_p$ distributions for Au+Au collisions at $\sqrt{s_{NN}}$ = 200 GeV.

Centrality dependence of moments of $\Delta N_p$ distributions for Au+Au collisions at $\sqrt{s_{NN}}$ = 200 GeV.

Centrality dependence of $S\sigma$ for $\Delta N_p$ in Au+Au collisions at $\sqrt{s_{NN}}$ = 19.6 GeV.

Centrality dependence of $\kappa\sigma^2$ for $\Delta N_p$ in Au+Au collisions at $\sqrt{s_{NN}}$ = 19.6 GeV.

Centrality dependence of $S\sigma$ for $\Delta N_p$ in Au+Au collisions at $\sqrt{s_{NN}}$ = 62.4 GeV.

Centrality dependence of $\kappa\sigma^2$ for $\Delta N_p$ in Au+Au collisions at $\sqrt{s_{NN}}$ = 62.4 GeV.

Centrality dependence of $S\sigma$ for $\Delta N_p$ in Au+Au collisions at $\sqrt{s_{NN}}$ = 200 GeV.

Centrality dependence of $\kappa\sigma^2$ for $\Delta N_p$ in Au+Au collisions at $\sqrt{s_{NN}}$ = 200 GeV.

Centrality dependence of $S\sigma$ for $\Delta N_p$ in Au+Au collisions from Lattice QCD Calculations.

Centrality dependence of $S\sigma$ for $\Delta N_p$ in Au+Au collisions at $\sqrt{s_{NN}}$ = 19.6 GeV from Transport Model Calculations.

Centrality dependence of $S\sigma$ for $\Delta N_p$ in Au+Au collisions at $\sqrt{s_{NN}}$ = 62.4 GeV from Transport Model Calculations.

Centrality dependence of $S\sigma$ for $\Delta N_p$ in Au+Au collisions at $\sqrt{s_{NN}}$ = 200 GeV from Transport Model Calculations.

Centrality dependence of $\kappa\sigma^2$ for $\Delta N_p$ in Au+Au collisions at $\sqrt{s_{NN}}$ = 200 GeV from Transport Model Calculations.

Centrality dependence of $\kappa\sigma^2$ for $\Delta N_p$ in Au+Au collisions at $\sqrt{s_{NN}}$ = 200 GeV from Transport Model Calculations for net-baryon.

Centrality dependence of $\kappa\sigma^2$ for $\Delta N_p$ in Au+Au collisions at $\sqrt{s_{NN}}$ = 200 GeV from Transport Model Calculations for net-proton (No Decay).

Centrality dependence of $\kappa\sigma^2$ for $\Delta N_p$ in Au+Au collisions at $\sqrt{s_{NN}}$ = 200 GeV from Transport Model Calculations.

Centrality dependence of $\kappa\sigma^2$ for $\Delta N_p$ in Au+Au collisions at $\sqrt{s_{NN}}$ = 200 GeV from Transport Model Calculations.

Centrality dependence of $\kappa\sigma^2$ for $\Delta N_p$ in Au+Au collisions at $\sqrt{s_{NN}}$ = 200 GeV from Transport Model Calculations.

Centrality dependence of $\kappa\sigma^2$ for $\Delta N_p$ in Au+Au collisions at $\sqrt{s_{NN}}$ = 200 GeV from Transport Model Calculations.

Centrality dependence of $\kappa\sigma^2$ for $\Delta N_p$ in Au+Au collisions at $\sqrt{s_{NN}}$ = 200 GeV from Transport Model Calculations.

Centrality dependence of $\kappa\sigma^2$ for $\Delta N_p$ in Au+Au collisions at $\sqrt{s_{NN}}$ = 200 GeV from Transport Model Calculations.

$\sqrt{s_{NN}}$ dependence of $\kappa\sigma^2$ for net proton distributions measured at RHIC. The results are STAR data.

$\sqrt{s_{NN}}$ dependence of $\kappa\sigma^2$ for net proton distributions measured at RHIC. The results are from UrQMD net-proton.

$\sqrt{s_{NN}}$ dependence of $\kappa\sigma^2$ for net proton distributions measured at RHIC. The results are from AMPT default net-proton.

$\sqrt{s_{NN}}$ dependence of $\kappa\sigma^2$ for net proton distributions measured at RHIC. The results are from AMPT String Melting net-proton.

$\sqrt{s_{NN}}$ dependence of $\kappa\sigma^2$ for net proton distributions measured at RHIC. The results are from Therminator net-proton.

$\sqrt{s_{NN}}$ dependence of $\kappa\sigma^2$ for net proton distributions measured at RHIC. The results are from Hijing net-proton.

Sample correlations in $\Delta\phi$ ($|\Delta\eta|<$ 1.78) for 3 $<$ $p_T^{trigger}$ $<$ 6 GeV/$c$ and 1.5 GeV/$c$ $<$ $p_T^{associated}$ $<$ $p_T^{trigger}$ for 0-60% Cu+Cu at $\sqrt{s_{NN}}$ = 62.4 GeV, 0-80% Au+Au at $\sqrt{s_{NN}}$ = 62.4 GeV, 0-95% $d$+Au at $\sqrt{s_{NN}}$ = 200 GeV, 0-60% Cu+Cu at $\sqrt{s_{NN}}$ = 200 GeV, 40-80% Au+Au at $\sqrt{s_{NN}}$ = 200 GeV, and 0-12% central Au+Au at $\sqrt{s_{NN}}$ = 200 GeV. The data are averaged between positive and negative $\Delta\phi$ and reflected in the plot.

Background subtracted sample correlations for 3 $<$ $p_T^{trigger}$ $<$ 6 GeV/$c$ and 1.5 GeV/$c$ $<$ $p_T^{associated}$ $<$ $p_T^{trigger}$ on the near-side for 0-60% Cu+Cu at $\sqrt{s_{NN}}$ = 62.4 GeV, 0-80% Au+Au at $\sqrt{s_{NN}}$ = 62.4 GeV, 0-95% $d$+Au collisions at $\sqrt{s_{NN}}$ = 200 GeV, 0-60% Cu+Cu at $\sqrt{s_{NN}}$ = 200 GeV, 40-80% Au+Au at $\sqrt{s_{NN}}$ = 200 GeV and 0-12% central Au+Au at $\sqrt{s_{NN}}$ = 200 GeV. The dependence of the jet-like correlation is shown as a function of $\Delta\eta$ ($|\Delta\phi|<$ 0.78). The data are averaged between positive and negative $\Delta\eta$.

Background subtracted sample correlations for 3 $<$ $p_T^{trigger}$ $<$ 6 GeV/$c$ and 1.5 GeV/$c$ $<$ $p_T^{associated}$ $<$ $p_T^{trigger}$ on the near-side for 0-60% Cu+Cu at $\sqrt{s_{NN}}$ = 62.4 GeV, 0-80% Au+Au at $\sqrt{s_{NN}}$ = 62.4 GeV, 0-95% $d$+Au collisions at $\sqrt{s_{NN}}$ = 200 GeV, 0-60% Cu+Cu at $\sqrt{s_{NN}}$ = 200 GeV, 40-80% Au+Au at $\sqrt{s_{NN}}$ = 200 GeV and 0-12% central Au+Au at $\sqrt{s_{NN}}$ = 200 GeV. The dependence of the jet-like correlation is shown as a function of $\Delta\phi$ ($|\Delta\eta|<$ 1.78). The data are averaged between positive and negative $\Delta\phi$.

Dependence of jet-like yield on $N_{part}$ for 3 $<$ $p_T^{trigger}$ $<$ 6 GeV/$c$ and 1.5 GeV/$c$ $<$ $p_T^{associated}$ $<$ $p_T^{trigger}$ for Cu+Cu at $\sqrt{s_{NN}}$ = 62.4 GeV. 5% systematic error due to the uncertainty on the associated particle's efficiency is not shown.

Dependence of jet-like yield on $N_{part}$ for 3 $<$ $p_T^{trigger}$ $<$ 6 GeV/$c$ and 1.5 GeV/$c$ $<$ $p_T^{associated}$ $<$ $p_T^{trigger}$ for Au+Au at $\sqrt{s_{NN}}$ = 62.4 GeV. 5% systematic error due to the uncertainty on the associated particle's efficiency is not shown.

Dependence of jet-like yield on $N_{part}$ for 3 $<$ $p_T^{trigger}$ $<$ 6 GeV/$c$ and 1.5 GeV/$c$ $<$ $p_T^{associated}$ $<$ $p_T^{trigger}$ for $d$+Au at $\sqrt{s_{NN}}$ = 200 GeV. 5% systematic error due to the uncertainty on the associated particle's efficiency is not shown.

Dependence of jet-like yield on $N_{part}$ for 3 $<$ $p_T^{trigger}$ $<$ 6 GeV/$c$ and 1.5 GeV/$c$ $<$ $p_T^{associated}$ $<$ $p_T^{trigger}$ for Cu+Cu at $\sqrt{s_{NN}}$ = 200 GeV. 5% systematic error due to the uncertainty on the associated particle's efficiency is not shown.

Dependence of jet-like yield on $N_{part}$ for 3 $<$ $p_T^{trigger}$ $<$ 6 GeV/$c$ and 1.5 GeV/$c$ $<$ $p_T^{associated}$ $<$ $p_T^{trigger}$ for Au+Au at $\sqrt{s_{NN}}$ = 200 GeV. 5% systematic error due to the uncertainty on the associated particle's efficiency is not shown.

Dependence of jet-like yield on $p_T^{\mathrm{trigger}}$ for 0-60% Cu+Cu at $\sqrt{s_{NN}}$ = 62.4 GeV, 0-80% Au+Au at $\sqrt{s_{NN}}$ = 62.4 GeV, 0-95% $d$+Au collisions at $\sqrt{s_{NN}}$ = 200 GeV, 0-60% Cu+Cu at $\sqrt{s_{NN}}$ = 200 GeV, 40-80% Au+Au at $\sqrt{s_{NN}}$ = 200 GeV and 0-12% central Au+Au at $\sqrt{s_{NN}}$ = 200 GeV. The 5% systematic error due to the uncertainty on the associated particle's efficiency is not shown.

Dependence of jet-like yield on $p_T^{\mathrm{associated}}$ for 3 $<$ $p_T^{trigger}$ $<$ 6 GeV/$c$ for 0-60% Cu+Cu and 0-80% Au+Au collisions at $\sqrt{s_{NN}}$ = 62.4 GeV and 0-95% $d$+Au, 0-60% Cu+Cu, 0-12% Au+Au and 40-80% Au+Au collisions at $\sqrt{s_{NN}}$ = 200 GeV. The 5% systematic error due to the uncertainty on the associated particle's efficiency is not shown.

Dependence of the widths in $\Delta\phi$ on $p_T^{\mathrm{trigger}}$ for 1.5 GeV/$c$ $<$ $p_T^{associated}$ $<$ $p_T^{trigger}$ for 0-95% $d$+Au, 0-60% Cu+Cu at $\sqrt{s_{NN}}$ = 62.4 GeV and $\sqrt{s_{NN}}$ = 200 GeV, 0-80% Au+Au at $\sqrt{s_{NN}}$ = 62.4 GeV, and 0-12% and 40-80% Au+Au at $\sqrt{s_{NN}}$ = 200 GeV.

Dependence of the widths in $\Delta\phi$ on $p_T^{\mathrm{associated}}$ for 3 $<$ $p_T^{trigger}$ $<$ 6 GeV/$c$ for 0-95% $d$+Au, 0-60% Cu+Cu at $\sqrt{s_{NN}}$ = 62.4 GeV and $\sqrt{s_{NN}}$ = 200 GeV, 0-80% Au+Au at $\sqrt{s_{NN}}$ = 62.4 GeV, and 0-12% and 40-80% Au+Au at $\sqrt{s_{NN}}$ = 200 GeV.

Dependence of the widths in $\Delta\phi$ on $N_{part}$ for 3 $<$ $p_T^{trigger}$ $<$ 6 GeV/$c$ and 1.5 GeV/$c$ $<$ $p_T^{associated}$ $<$ $p_T^{trigger}$ for 0-60% Cu+Cu at $\sqrt{s_{NN}}$ = 62.4 GeV.

Dependence of the widths in $\Delta\phi$ on $N_{part}$ for 3 $<$ $p_T^{trigger}$ $<$ 6 GeV/$c$ and 1.5 GeV/$c$ $<$ $p_T^{associated}$ $<$ $p_T^{trigger}$ for 0-80% Au+Au at $\sqrt{s_{NN}}$ = 62.4 GeV.

Dependence of the widths in $\Delta\phi$ on $N_{part}$ for 3 $<$ $p_T^{trigger}$ $<$ 6 GeV/$c$ and 1.5 GeV/$c$ $<$ $p_T^{associated}$ $<$ $p_T^{trigger}$ for 0-95% $d$+Au at $\sqrt{s_{NN}}$ = 200 GeV.

Dependence of the widths in $\Delta\phi$ on $N_{part}$ for 3 $<$ $p_T^{trigger}$ $<$ 6 GeV/$c$ and 1.5 GeV/$c$ $<$ $p_T^{associated}$ $<$ $p_T^{trigger}$ for 0-60% Cu+Cu at $\sqrt{s_{NN}}$ = 200 GeV.

Dependence of the widths in $\Delta\phi$ on $N_{part}$ for 3 $<$ $p_T^{trigger}$ $<$ 6 GeV/$c$ and 1.5 GeV/$c$ $<$ $p_T^{associated}$ $<$ $p_T^{trigger}$ for Au+Au at $\sqrt{s_{NN}}$ = 200 GeV.

Dependence of the widths in $\Delta\eta$ on $p_T^{\mathrm{trigger}}$ for 1.5 GeV/$c$ $<$ $p_T^{associated}$ $<$ $p_T^{trigger}$ for 0-95% $d$+Au, 0-60% Cu+Cu at $\sqrt{s_{NN}}$ = 62.4 GeV and $\sqrt{s_{NN}}$ = 200 GeV, 0-80% Au+Au at $\sqrt{s_{NN}}$ = 62.4 GeV, and 0-12% and 40-80% Au+Au at $\sqrt{s_{NN}}$ = 200 GeV.

Dependence of the widths in $\Delta\eta$ on $p_T^{\mathrm{associated}}$ for 3 $<$ $p_T^{trigger}$ $<$ 6 GeV/$c$ for 0-95% $d$+Au, 0-60% Cu+Cu at $\sqrt{s_{NN}}$ = 62.4 GeV and $\sqrt{s_{NN}}$ = 200 GeV, 0-80% Au+Au at $\sqrt{s_{NN}}$ = 62.4 GeV, and 0-12% and 40-80% Au+Au at $\sqrt{s_{NN}}$ = 200 GeV.

Dependence of the widths in $\Delta\eta$ on $N_{part}$ for 3 $<$ $p_T^{trigger}$ $<$ 6 GeV/$c$ and 1.5 GeV/$c$ $<$ $p_T^{associated}$ $<$ $p_T^{trigger}$ for 0-60% Cu+Cu at $\sqrt{s_{NN}}$ = 62.4 GeV.

Dependence of the widths in $\Delta\eta$ on $N_{part}$ for 3 $<$ $p_T^{trigger}$ $<$ 6 GeV/$c$ and 1.5 GeV/$c$ $<$ $p_T^{associated}$ $<$ $p_T^{trigger}$ for 0-80% Au+Au at $\sqrt{s_{NN}}$ = 62.4 GeV.

Dependence of the widths in $\Delta\eta$ on $N_{part}$ for 3 $<$ $p_T^{trigger}$ $<$ 6 GeV/$c$ and 1.5 GeV/$c$ $<$ $p_T^{associated}$ $<$ $p_T^{trigger}$ for 0-95% $d$+Au at $\sqrt{s_{NN}}$ = 200 GeV.

Dependence of the widths in $\Delta\eta$ on $N_{part}$ for 3 $<$ $p_T^{trigger}$ $<$ 6 GeV/$c$ and 1.5 GeV/$c$ $<$ $p_T^{associated}$ $<$ $p_T^{trigger}$ for 0-60% Cu+Cu at $\sqrt{s_{NN}}$ = 200 GeV.

Dependence of the widths in $\Delta\eta$ on $N_{part}$ for 3 $<$ $p_T^{trigger}$ $<$ 6 GeV/$c$ and 1.5 GeV/$c$ $<$ $p_T^{associated}$ $<$ $p_T^{trigger}$ for Au+Au at $\sqrt{s_{NN}}$ = 200 GeV.

Dependence of ridge yield on $N_{part}$ for 3 $<$ $p_T^{trigger}$ $<$ 6 GeV/$c$ and 1.5 GeV/$c$ $<$ $p_T^{associated}$ $<$ $p_T^{trigger}$ for Cu+Cu at $\sqrt{s_{NN}}$ = 62.4 GeV.

Dependence of ridge yield on $N_{part}$ for 3 $<$ $p_T^{trigger}$ $<$ 6 GeV/$c$ and 1.5 GeV/$c$ $<$ $p_T^{associated}$ $<$ $p_T^{trigger}$ for Au+Au at $\sqrt{s_{NN}}$ = 62.4 GeV.

Dependence of ridge yield on $N_{part}$ for 3 $<$ $p_T^{trigger}$ $<$ 6 GeV/$c$ and 1.5 GeV/$c$ $<$ $p_T^{associated}$ $<$ $p_T^{trigger}$ for Cu+Cu at $\sqrt{s_{NN}}$ = 200 GeV.

Dependence of ridge yield on $N_{part}$ for 3 $<$ $p_T^{trigger}$ $<$ 6 GeV/$c$ and 1.5 GeV/$c$ $<$ $p_T^{associated}$ $<$ $p_T^{trigger}$ for Au+Au at $\sqrt{s_{NN}}$ = 200 GeV.

Ratio of the ridge and jet-like yields on $N_{part}$ for 3 $<$ $p_T^{trigger}$ $<$ 6 GeV/$c$ and 1.5 GeV/$c$ $<$ $p_T^{associated}$ $<$ $p_T^{trigger}$ for Cu+Cu at $\sqrt{s_{NN}}$ = 62.4 GeV.

Ratio of the ridge and jet-like yields on $N_{part}$ for 3 $<$ $p_T^{trigger}$ $<$ 6 GeV/$c$ and 1.5 GeV/$c$ $<$ $p_T^{associated}$ $<$ $p_T^{trigger}$ for Au+Au at $\sqrt{s_{NN}}$ = 62.4 GeV.

Ratio of the ridge and jet-like yields on $N_{part}$ for 3 $<$ $p_T^{trigger}$ $<$ 6 GeV/$c$ and 1.5 GeV/$c$ $<$ $p_T^{associated}$ $<$ $p_T^{trigger}$ for Cu+Cu at $\sqrt{s_{NN}}$ = 200 GeV.

Ratio of the ridge and jet-like yields on $N_{part}$ for 3 $<$ $p_T^{trigger}$ $<$ 6 GeV/$c$ and 1.5 GeV/$c$ $<$ $p_T^{associated}$ $<$ $p_T^{trigger}$ for Au+Au at $\sqrt{s_{NN}}$ = 200 GeV.

Dependence of $V_{3\Delta}/V_{2\Delta}$ on $N_{part}$ for 3 $<$ $p_T^{trigger}$ $<$ 6 GeV/$c$ and 1.5 GeV/$c$ $<$ $p_T^{associated}$ $<$ $p_T^{trigger}$ for Cu+Cu at $\sqrt{s_{NN}}$ = 62.4 GeV.

Dependence of $V_{3\Delta}/V_{2\Delta}$ on $N_{part}$ for 3 $<$ $p_T^{trigger}$ $<$ 6 GeV/$c$ and 1.5 GeV/$c$ $<$ $p_T^{associated}$ $<$ $p_T^{trigger}$ for Au+Au at $\sqrt{s_{NN}}$ = 62.4 GeV.

Dependence of $V_{3\Delta}/V_{2\Delta}$ on $N_{part}$ for 3 $<$ $p_T^{trigger}$ $<$ 6 GeV/$c$ and 1.5 GeV/$c$ $<$ $p_T^{associated}$ $<$ $p_T^{trigger}$ for $d$+Au at $\sqrt{s_{NN}}$ = 200 GeV.

Dependence of $V_{3\Delta}/V_{2\Delta}$ on $N_{part}$ for 3 $<$ $p_T^{trigger}$ $<$ 6 GeV/$c$ and 1.5 GeV/$c$ $<$ $p_T^{associated}$ $<$ $p_T^{trigger}$ for Cu+Cu at $\sqrt{s_{NN}}$ = 200 GeV.

Dependence of $V_{3\Delta}/V_{2\Delta}$ on $N_{part}$ for 3 $<$ $p_T^{trigger}$ $<$ 6 GeV/$c$ and 1.5 GeV/$c$ $<$ $p_T^{associated}$ $<$ $p_T^{trigger}$ for Au+Au at $\sqrt{s_{NN}}$ = 200 GeV.

Charged hadron v2 for the centrality bin 0-50%.

v2(4) vs. p_T for pions for 20 - 60% centrality

v2(4) vs. p_T for kaons for 20 - 60% centrality

v2(4) vs. p_T for anti-protons for 20 - 60% centrality

v2(2) vs. centrality for pions

v2(2) vs. centrality for kaons

v2(2) vs. centrality for anti-protons

v2 vs. p_T for particles identified in the RICH (min. bias) for pions+kaons and protons/antiprotons

v2 for Kshort for min. bias

v2 for kaons from kinks for min. bias

v2 for kaons from dE/dx for min. bias

v2 for Lambda+Lambdabar for min. bias [PRL/92/052302/2004]. The PHENIX data are from [PRL/91/182301/2003]. The hydro calculations are from [P. Huovinen/private communication].

v2 for charged pions for min. bias. kaon data points are from Figure 9.

v2 for anti-protons for min. bias. kaon data points are from Figure 9.

Azimuthal correlations (\sum_i cos(2(\phi_pT -\phi_i))) in Au+Au collisions as a function of p_T along with results from p+p collisions.

Azimuthal correlations (\sum_i cos(2(\phi_pT -\phi_i))) from d+Au as a function of p_T for min bias collisions.

Azimuthal correlations (\sum_i cos(2(\phi_pT -\phi_i))) from d+Au as a function of p_T for three centralities.

Charged hadron v2 and v4{EP_2} as a function of pseudorapidity for minimum bias.

v2(2) and v2(4) vs. pseudorapidity

v2(2) and v2(4) vs. pseudorapidity,FTPC only

v2(4) 20-70% in TPC+FTPC and v2(2) 20-70% in FTPC.

MC v2 with FTPC momentum resolution.

Charged hadron v2 for the centrality bin 0-5% in standard and modified event plane methods.

Charged hadron v2 for the centrality bin 5-10% in standard and modified event plane methods.

Charged hadron v2 for the centrality bin 10-20% in standard and modified event plane methods.

Charged hadron v2 for the centrality bin 20-30% in standard and modified event plane methods.

Charged hadron v2 for the centrality bin 30-40% in standard and modified event plane methods.

Charged hadron v2 for the centrality bin 40-60% in standard and modified event plane methods.

Charged hadron v2{2} from modified reaction plane method and two-particle cumulant results (after subtracting the correlations measured in p+p collisions) for centrality bin 20-60% .

v4(EP2) for charged hadrons min. bias

v4(EP2) vs. centrality for charged hadrons

v4(EP2) for pions in min. bias

v4(EP2) for protons in min. bias

v4(EP2) for kaons in min. bias

v4(EP2) for Lambda in min. bias

Charged hadron v4{EP_2} as a function of pseudorapidity for minimum bias (TPC only).

charged hadrons high p_T (3 - 6 GeV/c) v4(3)

charged hadron integrated v2 as a function of centrality for the different methods

charged hadron v2(6) vs. centrality bin. v2(2) and v2(4) are from Figure 29(a).

g2 as a function of centrality from 200 GeV, 130 GeV,NA49 and q-dist methods. The solid points are from the cumulant method, while open circles are from the q distribution method. Please note that centrality percentage differs from different method. For 200 GeV-cumulant (0, 64.1, 53.9,44.7,35.2,25.4,15.1,7.7,0), for 200 GeV-q-dist (73.8,64.1,53.9,44.7,35.2,25.4,15.1,7.7,2.3), for 130 GeV (65,47,36,27.5,20,13,7.5,2.5,0) and for 17 GeV (0,0,0,61.2,38.2,29.1,18.2,8.5,3.8).

wounded nucleons as a function of centrality from 200 GeV, 130 GeV,NA49 and q-dist methods.Please note that centrality percentage differs from different method. For 200 GeV-cumulant (0, 64.1, 53.9,44.7,35.2,25.4,15.1,7.7,0), for 200 GeV-q-dist (73.8,64.1,53.9,44.7,35.2,25.4,15.1,7.7,2.3), for 130 GeV (65,47,36,27.5,20,13,7.5,2.5,0) and for 17 GeV (0,0,0,61.2,38.2,29.1,18.2,8.5,3.8).

v2/n for pions for three centrality bins

v2/n for protons for three centrality bins

v2/n for Kshort for three centrality bins

v2/n for Lambda + Lambda bar for three centrality bins

Blast Wave parameters rho_2 from pion v2(2),pion v2 and hadron v2

Blast Wave parameters rho_4 from hadron v2

Blast Wave parameters epsilon from pion v2(2),pion v2, hadron v2 and initial

Blast Wave parameters rho_0 from pion v2(2),pion v2 and hadron v2

Blast Wave parameters T from pion v2(2),pion v2 and hadron v2

HBT parameters for the three possible fitting procedures to the correlation functions described in this paper depending on how Coulomb interaction is taken into account from the 0-5% most central events.

HBT parameters for the three possible fitting procedures to the correlation functions described in this paper depending on how Coulomb interaction is taken into account from the 0-5% most central events.

HBT parameters for the three possible fitting procedures to the correlation functions described in this paper depending on how Coulomb interaction is taken into account from the 0-5% most central events.

HBT parameters for the three possible fitting procedures to the correlation functions described in this paper depending on how Coulomb interaction is taken into account from the 0-5% most central events.

1D correlation function for pi+pi- compared to standard, Bowler-Sinyukov, and a thoretical calculation that includes Coulomb and strong interactions.

Momentum resolution for pions at midrapidity expressed by the widths dp_T/p_T, dphi and dtheta as a function of p.

Squared HBT radii relative to the reaction plane angle, without and with the reaction plane resolution applied, for two different centrality-kT ranges. 20-30%, 0.15 < kT < 0.25, RP uncorrected.

Squared HBT radii relative to the reaction plane angle, without and with the reaction plane resolution applied, for two different centrality-kT ranges. 20-30%, 0.15 < kT < 0.25, RP res. corrected.

Squared HBT radii relative to the reaction plane angle, without and with the reaction plane resolution applied, for two different centrality-kT ranges. 10-20%, 0.35 < kT < 0.45, RP uncorrected.

Squared HBT radii relative to the reaction plane angle, without and with the reaction plane resolution applied, for two different centrality-kT ranges. 10-20%, 0.35 < kT < 0.45, RP res. uncorrected.

Squared HBT radii relative to the reaction plane angle, for the case where the lambda parameter is averaged and fixed in the fit, and the case where lambda is a free fit parameter for each Phi. 20-30%, 0.15 < kT < 0.25, fixed lambda.

Squared HBT radii relative to the reaction plane angle, for the case where the lambda parameter is averaged and fixed in the fit, and the case where lambda is a free fit parameter for each Phi. 20-30%, 0.15 < kT < 0.25, lambda varied.

Squared HBT radii relative to the reaction plane angle, for the case where the lambda parameter is averaged and fixed in the fit, and the case where lambda is a free fit parameter for each Phi. 10-20%, 0.35 < kT < 0.45, lambda fixed.

Squared HBT radii relative to the reaction plane angle, for the case where the lambda parameter is averaged and fixed in the fit, and the case where lambda is a free fit parameter for each Phi. 10-20%, 0.35 < kT < 0.45, lambda varied.

Fourier coefficients as a function of the maximum fraction of merged hits for the 5% most central events and kT between 150 and 250 MeV/c.

Projections of the 3 dimensional correlation function and fits to Eq. (10) and with the Edgeworth expansion to Eq. (17) to 6th order.

HBT parameters for the 0-5% most central events for fits to Eq. (10) and to Eq. (17) with no expansion.

HBT parameters for the 0-5% most central events for fits to Eq. (10) and to Eq. (17) upto 4th order expansion.

HBT parameters for the 0-5% most central events for fits to Eq. (10) and to Eq. (17) upto 6th order expansion.

HBT parameters for the 0-5% most central events for pi+pi+ and pi-pi- correlation functions.

HBT parameters from STAR and PHENIX at the same beam energy for the 0-30% most central events.

Energy dependence of the pi- HBT parameters for central Au+Au, Pb+Pb, and Pb+Au collisions at midrapidity and k_T ~ 0.2 GeV/c.

HBT parameters vs. m_T for 6 different centralities.

Fourier coefficients of azimuthal oscillations of HBT radii vs. numer of participating nucleons, for pi+ and pi- pairs separately (0.25 < kT < 0.35 GeV/c).

Comparison between the HBT radii obtained from the azimuthally integrated (traditional) HBT analysis and the 0th-order Fourier coefficients from the azimuthally-sensitive HBT analysis.

HBT radius parametres for 6 different centralities.

Extracted parameters R' and alpha, from the power-law fits to the HBT radius parameters.

Extracted freeze-out source radius extracted from a blast wave fit; source radius Rgeom from fits to Rside; and 2*Rside for the lowest k_T bin as a function of the number of participants.

HBT parameter Rside.

R - Rinitial (top panel) and R/Rinitial (bottom panel) for the azimuthally integrated analysis for the azimuthally-sensitive case vs. number of participants.

R - Rinitial (top panel) and R/Rinitial (bottom panel) for in the x (in-plane) and y (out-of-plane) directions for the azimuthally-sensitive case vs. number of participants.

R - Rinitial for the azimuthally integrated analysis for the azimuthally-sensitive case vs. (dN/dy)/Rinitial.

R - Rinitial in the x (in-plane) direction for the azimuthally-sensitive case vs. (dN/dy)/Rinitial.

R - Rinitial in the y (out-of-plane) direction for the azimuthally-sensitive case vs. (dN/dy)/Rinitial.

Longitudinal HBT radius Rl.

Evolution time tau vs. number of participants as extracted from a fit to Rl, lines in Fig. 26 (triangles), and from a blast wave fit to HBT parameters and spectra (circles).

R - Rinital for the azimuthally integrated analysis and for the in-plane and out-of-plane directions vs. beta_T,max*tau.

R - Rinital for the azimuthally integrated analysis and for the in-plane and out-of-plane directions vs. beta_T,max*tau.

R - Rinital for the azimuthally integrated analysis and for the in-plane and out-of-plane directions vs. beta_T,max*tau.

Emission duration time dTau vs. number of participants as extracted using a blast wave fit to HBT parameters and spectra.

Final source eccentricity as calculated from the Fourier coefficients (2R_{s,2}^2/R_{s,0}^2) and from the final in-plane and out-of-plane radii vs. initial eccentricity.

Transverse momentum distribution for $\pi^+$ production in d+Au collisions with centrality 20-40% in the mid rapidity region, $|y|<0.5$.

Transverse momentum distribution for $\pi^+$ production in d+Au collisions with centrality 40-100% in the mid rapidity region, $|y|<0.5$.

Transverse momentum distribution for $\pi^-$ production in d+Au minbias events in the mid rapidity region, $|y|<0.5$.

Transverse momentum distribution for $\pi^-$ production in p+p NSD events in the mid rapidity region, $|y|<0.5$.

Transverse momentum distribution for $\pi^-$ production in d+Au collisions with centrality 0-20% in the mid rapidity region, $|y|<0.5$.

Transverse momentum distribution for $\pi^-$ production in d+Au collisions with centrality 20-40% in the mid rapidity region, $|y|<0.5$.

Transverse momentum distribution for $\pi^-$ production in d+Au collisions with centrality 40-100% in the mid rapidity region, $|y|<0.5$.

Transverse momentum distribution for proton production in d+Au minbias events in the mid rapidity region, $|y|<0.5$.

Transverse momentum distribution for proton production in p+p NSD events in the mid rapidity region, $|y|<0.5$.

Transverse momentum distribution for proton production in d+Au collisions with centrality 0-20% in the mid rapidity region, $|y|<0.5$.

Transverse momentum distribution for proton production in d+Au collisions with centrality 20-40% in the mid rapidity region, $|y|<0.5$.

Transverse momentum distribution for proton production in d+Au collisions with centrality 40-100% in the mid rapidity region, $|y|<0.5$.

Transverse momentum distribution for pbar production in d+Au minbias events in the mid rapidity region, $|y|<0.5$.

Transverse momentum distribution for pbar production in p+p NSD events in the mid rapidity region, $|y|<0.5$.

Transverse momentum distribution for pbar production in d+Au collisions with centrality 0-20% in the mid rapidity region, $|y|<0.5$.

Transverse momentum distribution for pbar production in d+Au collisions with centrality 20-40% in the mid rapidity region, $|y|<0.5$.

Transverse momentum distribution for pbar production in d+Au collisions with centrality 40-100% in the mid rapidity region, $|y|<0.5$.

Nuclear modification factors $R_{dAu}$ for $\pi^+$ + $\pi^-$ production in the mid rapidity region, $|y|<0.5$. There is an overall additional normalization uncertainty of the order of 17%. For the measurement of $R_{dAu}$ there is an additional uncertainty of 5% due to uncertainty on $N_{bin}$ for minbias collisions.

Nuclear modification factors $R_{dAu}$ for $p + \bar{p}$ production in the mid rapidity region, $|y|<0.5$. There is an overall additional normalization uncertainty of the order of 17%. For the measurement of $R_{dAu}$ there is an additional uncertainty of 5% due to uncertainty on $N_{bin}$ for minbias collisions.

Ratios of $ \pi^-/\pi^+ $ in DEUT AU and P P minbias collisions in the mid rapidity region, $|y|<0.5$. In addition there are approximately 7% correlated errors for the ratios in P P collisions and about 10% in the DEUT AU collisions.

Ratios of $ \bar{p}/p $ in DEUT AU and P P minbias collisions in the mid rapidity region, $|y|<0.5$. In addition there are approximately 7% correlated errors for the ratios in P P collisions and about 10% in the DEUT AU collisions.

Ratios of $ p/\pi^+ $ in DEUT AU and P P minbias collisions in the mid rapidity region, $|y|<0.5$. In addition there are approximately 7% correlated errors for the ratios in P P collisions and about 10% in the DEUT AU collisions.

Ratios of $ \bar{p}/\pi^- $ in DEUT AU and P P minbias collisions in the mid rapidity region, $|y|<0.5$. In addition there are approximately 7% correlated errors for the ratios in P P collisions and about 10% in the DEUT AU collisions.