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.

Fig. 2. (Color online.) Event-by-event photon multiplicity distributions (solid circles) for $\mathrm{Au}+\mathrm{Au}$ and $\mathrm{Cu}+\mathrm{Cu}$ at $\sqrt{s_{\mathrm{NN}}}=62.4$ and $200 \mathrm{GeV} .$ The distributions for top $0-5 \%$ central $\mathrm{Au}+$ Au collisions and top $0-10 \%$ central $\mathrm{Cu}+\mathrm{Cu}$ collisions are also shown (open circles). The photon multiplicity distributions for central collisions are observed to be Gaussian (solid line). Only statistical errors are shown. NOTE: For points with invisible error bars, the point size was considered as an absolute upper limit for the uncertainty.

Fig. 2. (Color online.) Event-by-event photon multiplicity distributions (solid circles) for $\mathrm{Au}+\mathrm{Au}$ and $\mathrm{Cu}+\mathrm{Cu}$ at $\sqrt{s_{\mathrm{NN}}}=62.4$ and $200 \mathrm{GeV} .$ The distributions for top $0-5 \%$ central $\mathrm{Au}+$ Au collisions and top $0-10 \%$ central $\mathrm{Cu}+\mathrm{Cu}$ collisions are also shown (open circles). The photon multiplicity distributions for central collisions are observed to be Gaussian (solid line). Only statistical errors are shown. NOTE: For points with invisible error bars, the point size was considered as an absolute upper limit for the uncertainty.

Fig. 2. (Color online.) Event-by-event photon multiplicity distributions (solid circles) for $\mathrm{Au}+\mathrm{Au}$ and $\mathrm{Cu}+\mathrm{Cu}$ at $\sqrt{s_{\mathrm{NN}}}=62.4$ and $200 \mathrm{GeV} .$ The distributions for top $0-5 \%$ central $\mathrm{Au}+$ Au collisions and top $0-10 \%$ central $\mathrm{Cu}+\mathrm{Cu}$ collisions are also shown (open circles). The photon multiplicity distributions for central collisions are observed to be Gaussian (solid line). Only statistical errors are shown. NOTE: For points with invisible error bars, the point size was considered as an absolute upper limit for the uncertainty.

Fig. 3. (Color online.) Photon pseudorapidity distributions for $\mathrm{Au}+\mathrm{Au}$ and $\mathrm{Cu}+\mathrm{Cu}$ at $\sqrt{\mathrm{s}_{\mathrm{NN}}}=62.4$ and $200\mathrm{GeV}$. The results for several centrality classes are shown. The solid curves are results of HIJING simulations for central $(0-5 \%$ for $\mathrm{Au}+\mathrm{Au}$ and $0-10 \%$ for $\mathrm{Cu}+\mathrm{Cu})$ and $30-40 \%$ mid-central collisions. The errors shown are systematic, statistical errors are negligible in comparison. NOTE: For points with invisible error bars, the point size was considered as an absolute upper limit for the uncertainty.

Fig. 3. (Color online.) Photon pseudorapidity distributions for $\mathrm{Au}+\mathrm{Au}$ and $\mathrm{Cu}+\mathrm{Cu}$ at $\sqrt{\mathrm{s}_{\mathrm{NN}}}=62.4$ and $200\mathrm{GeV}$. The results for several centrality classes are shown. The solid curves are results of HIJING simulations for central $(0-5 \%$ for $\mathrm{Au}+\mathrm{Au}$ and $0-10 \%$ for $\mathrm{Cu}+\mathrm{Cu})$ and $30-40 \%$ mid-central collisions. The errors shown are systematic, statistical errors are negligible in comparison. NOTE: For points with invisible error bars, the point size was considered as an absolute upper limit for the uncertainty.

Fig. 3. (Color online.) Photon pseudorapidity distributions for $\mathrm{Au}+\mathrm{Au}$ and $\mathrm{Cu}+\mathrm{Cu}$ at $\sqrt{\mathrm{s}_{\mathrm{NN}}}=62.4$ and $200\mathrm{GeV}$. The results for several centrality classes are shown. The solid curves are results of HIJING simulations for central $(0-5 \%$ for $\mathrm{Au}+\mathrm{Au}$ and $0-10 \%$ for $\mathrm{Cu}+\mathrm{Cu})$ and $30-40 \%$ mid-central collisions. The errors shown are systematic, statistical errors are negligible in comparison. NOTE: For points with invisible error bars, the point size was considered as an absolute upper limit for the uncertainty.

Fig. 3. (Color online.) Photon pseudorapidity distributions for $\mathrm{Au}+\mathrm{Au}$ and $\mathrm{Cu}+\mathrm{Cu}$ at $\sqrt{\mathrm{s}_{\mathrm{NN}}}=62.4$ and $200\mathrm{GeV}$. The results for several classes are shown. The solid curves are results of HIJING simulations for central $(0-5 \%$ for $\mathrm{Au}+\mathrm{Au}$ and $0-10 \%$ for $\mathrm{Cu}+\mathrm{Cu})$ and $30-40 \%$ mid-central collisions. The errors shown are systematic, statistical errors are negligible in comparison. NOTE: For points with invisible error bars, the point size was considered as an absolute upper limit for the uncertainty.

Fig. 4. (Color online.) Top panel: The number of photons divided by $\left\langle N_{\text{part}}\right\rangle / 2$ as a function of average number of participating nucleons for $\mathrm{Au}+\mathrm{Au}$ and $\mathrm{Cu}+\mathrm{Cu}$ at $\sqrt{\mathrm{s} _{\mathrm{NN}}}=62.4$ and $200 \mathrm{GeV}$ for $-3.7<\eta<-2.3 .$ Errors shown are systematic only and include those for $\left\langle N_{\text {part }}\right\rangle .$ Results from HIJING are shown as lines (solid for $\mathrm{Au}+\mathrm{Au}$ and dashed for $\mathrm{Cu}+\mathrm{Cu}$ ). Bottom panel: Same as above, for both photons and charged particles from PHOBOS [10]. NOTE: For points with invisible error bars, the point size was considered as an absolute upper limit for the uncertainty.

Fig. 4. (Color online.) Top panel: The number of photons divided by $\left\langle N_{\text{part}}\right\rangle / 2$ as a function of average number of participating nucleons for $\mathrm{Au}+\mathrm{Au}$ and $\mathrm{Cu}+\mathrm{Cu}$ at $\sqrt{\mathrm{s} _{\mathrm{NN}}}=62.4$ and $200 \mathrm{GeV}$ for $-3.7<\eta<-2.3 .$ Errors shown are systematic only and include those for $\left\langle N_{\text {part }}\right\rangle .$ Results from HIJING are shown as lines (solid for $\mathrm{Au}+\mathrm{Au}$ and dashed for $\mathrm{Cu}+\mathrm{Cu}$ ). Bottom panel: Same as above, for both photons and charged particles from PHOBOS [10]. NOTE: For points with invisible error bars, the point size was considered as an absolute upper limit for the uncertainty.

Fig. 5. (Color online.) Photon pseudorapidity distributions normalized by the average number of participating nucleon pairs for different collision centralities are plotted as a function of pseudorapidity shifted by the beam rapidity (-5.36 for $200 \mathrm{GeV}$ and -4.19 for $62.4 \mathrm{GeV}$ ) for $\mathrm{Au}+\mathrm{Au}$ and $\mathrm{Cu}+\mathrm{Cu}$ collisions at $\sqrt{s_{\mathrm{NN}}}=62.4$ and $200 \mathrm{GeV}$. Errors are systematic only, statistical errors are negligible in comparison. For clarity of presentation, results for only four centralities are shown. The $\mathrm{Cu}+$ Cu data are shifted by 0.1 unit in $\eta-y_{\text {beam }} .$ The solid line is a second order polynomial fit to the data (see text for details).

Fig. 5. (Color online.) Photon pseudorapidity distributions normalized by the average number of participating nucleon pairs for different collision centralities are plotted as a function of pseudorapidity shifted by the beam rapidity (-5.36 for $200 \mathrm{GeV}$ and -4.19 for $62.4 \mathrm{GeV}$ ) for $\mathrm{Au}+\mathrm{Au}$ and $\mathrm{Cu}+\mathrm{Cu}$ collisions at $\sqrt{s_{\mathrm{NN}}}=62.4$ and $200 \mathrm{GeV}$. Errors are systematic only, statistical errors are negligible in comparison. For clarity of presentation, results for only four centralities are shown. The $\mathrm{Cu}+$ Cu data are shifted by 0.1 unit in $\eta-y_{\text {beam }} .$ The solid line is a second order polynomial fit to the data (see text for details).

FIG. 1. (Color online) Perspective views of $2 \mathrm{D}$ charge-independent angular correlations $\Delta \rho / \sqrt{\rho_{\mathrm{ref}}}$ on $\left(\eta_{\Delta}, \phi_{\Delta}\right)$ for Au-Au collisions at $\sqrt{s_{N N}}=200$ and $62 \mathrm{GeV}$ (top and bottom rows, respectively). Centrality increases left to right from most peripheral to most central. Corrected total cross-section fractions are (left to right) $84 \%-93 \%, 55 \%-64 \%, 18 \%-28 \%,$ and $0 \%-5 \%$ for the $200-\mathrm{GeV}$ data and $84 \%-95 \%, 56 \%-65 \%$ $18 \%-28 \%,$ and $0 \%-5 \%$ for the $62 \mathrm{GeV}$ data (see Tables III and IV).

FIG. 1. (Color online) Perspective views of $2 \mathrm{D}$ charge-independent angular correlations $\Delta \rho / \sqrt{\rho_{\mathrm{ref}}}$ on $\left(\eta_{\Delta}, \phi_{\Delta}\right)$ for Au-Au collisions at $\sqrt{s_{N N}}=200$ and $62 \mathrm{GeV}$ (top and bottom rows, respectively). Centrality increases left to right from most peripheral to most central. Corrected total cross-section fractions are (left to right) $84 \%-93 \%, 55 \%-64 \%, 18 \%-28 \%,$ and $0 \%-5 \%$ for the $200-\mathrm{GeV}$ data and $84 \%-95 \%, 56 \%-65 \%$ $18 \%-28 \%,$ and $0 \%-5 \%$ for the $62 \mathrm{GeV}$ data (see Tables III and IV).

FIG. 1. (Color online) Perspective views of $2 \mathrm{D}$ charge-independent angular correlations $\Delta \rho / \sqrt{\rho_{\mathrm{ref}}}$ on $\left(\eta_{\Delta}, \phi_{\Delta}\right)$ for Au-Au collisions at $\sqrt{s_{N N}}=200$ and $62 \mathrm{GeV}$ (top and bottom rows, respectively). Centrality increases left to right from most peripheral to most central. Corrected total cross-section fractions are (left to right) $84 \%-93 \%, 55 \%-64 \%, 18 \%-28 \%,$ and $0 \%-5 \%$ for the $200-\mathrm{GeV}$ data and $84 \%-95 \%, 56 \%-65 \%$ $18 \%-28 \%,$ and $0 \%-5 \%$ for the $62 \mathrm{GeV}$ data (see Tables III and IV).

FIG. 1. (Color online) Perspective views of $2 \mathrm{D}$ charge-independent angular correlations $\Delta \rho / \sqrt{\rho_{\mathrm{ref}}}$ on $\left(\eta_{\Delta}, \phi_{\Delta}\right)$ for Au-Au collisions at $\sqrt{s_{N N}}=200$ and $62 \mathrm{GeV}$ (top and bottom rows, respectively). Centrality increases left to right from most peripheral to most central. Corrected total cross-section fractions are (left to right) $84 \%-93 \%, 55 \%-64 \%, 18 \%-28 \%,$ and $0 \%-5 \%$ for the $200-\mathrm{GeV}$ data and $84 \%-95 \%, 56 \%-65 \%$ $18 \%-28 \%,$ and $0 \%-5 \%$ for the $62 \mathrm{GeV}$ data (see Tables III and IV).

FIG. 1. (Color online) Perspective views of $2 \mathrm{D}$ charge-independent angular correlations $\Delta \rho / \sqrt{\rho_{\mathrm{ref}}}$ on $\left(\eta_{\Delta}, \phi_{\Delta}\right)$ for Au-Au collisions at $\sqrt{s_{N N}}=200$ and $62 \mathrm{GeV}$ (top and bottom rows, respectively). Centrality increases left to right from most peripheral to most central. Corrected total cross-section fractions are (left to right) $84 \%-93 \%, 55 \%-64 \%, 18 \%-28 \%,$ and $0 \%-5 \%$ for the $200-\mathrm{GeV}$ data and $84 \%-95 \%, 56 \%-65 \%$ $18 \%-28 \%,$ and $0 \%-5 \%$ for the $62 \mathrm{GeV}$ data (see Tables III and IV).

FIG. 2. (Color online) Fit decomposition of the $46 \%-56 \%$ centrality data for 62 -GeV Au-Au collisions. The top panels show from left to right the corrected data, model fit, fit residuals (data - model), and same-side 2D Gaussian. The bottom panels similarly show the away-side azimuth dipole, the nonjet azimuth quadrupole, the 1D $\eta_{\Delta}$ Gaussian, and the 2D exponential. This centrality is just below the sharp transition at $v_{\text {trans }}=3$. Fit residuals (c) are scaled up eightfold relative to the data.

FIG. 3. Fit parameters for $\left(\eta_{\Delta}, \phi_{\Delta}\right)$ correlation data from Au-Au collisions at $\sqrt{s_{N N}}=62$ (open symbols) and 200 GeV (solid symbols) versus centrality measure $\nu$ computed at fixed energy $(200 \mathrm{GeV})$. The SS $2 \mathrm{D}$ Gaussian amplitudes, $\eta_{\Delta}$ widths, and $\phi_{\Delta}$ widths are shown in the left, center, and right panels, respectively of the top row. The bottom row shows from left to right the amplitudes for the dipole, quadrupole, and SS peak width aspect ratio $\sigma_{\eta_{\Delta}} / \sigma_{\phi_{\Delta}} .$ Fitting errors are indicated by error bars where larger than the symbols. Solid lines connect the points for clarity. The dotted and dashed curves indicate Glauber linear superposition estimates for 62 - and 200 -GeV peak amplitudes respectively, as discussed in the text. The quadrupole data are consistent with Ref. [60]. The hatched regions indicate the full range of systematic uncertainties listed in Appendix F. The vertical dark bands indicate estimated $v$ equivalents for $N-N$ collisions and $b=0$ Au-Au collisions.

FIG. 3. Fit parameters for $\left(\eta_{\Delta}, \phi_{\Delta}\right)$ correlation data from Au-Au collisions at $\sqrt{s_{N N}}=62$ (open symbols) and 200 GeV (solid symbols) versus centrality measure $\nu$ computed at fixed energy $(200 \mathrm{GeV})$. The SS $2 \mathrm{D}$ Gaussian amplitudes, $\eta_{\Delta}$ widths, and $\phi_{\Delta}$ widths are shown in the left, center, and right panels, respectively of the top row. The bottom row shows from left to right the amplitudes for the dipole, quadrupole, and SS peak width aspect ratio $\sigma_{\eta_{\Delta}} / \sigma_{\phi_{\Delta}} .$ Fitting errors are indicated by error bars where larger than the symbols. Solid lines connect the points for clarity. The dotted and dashed curves indicate Glauber linear superposition estimates for 62 - and 200 -GeV peak amplitudes respectively, as discussed in the text. The quadrupole data are consistent with Ref. [60]. The hatched regions indicate the full range of systematic uncertainties listed in Appendix F. The vertical dark bands indicate estimated $v$ equivalents for $N-N$ collisions and $b=0$ Au-Au collisions.

FIG. 3. Fit parameters for $\left(\eta_{\Delta}, \phi_{\Delta}\right)$ correlation data from Au-Au collisions at $\sqrt{s_{N N}}=62$ (open symbols) and 200 GeV (solid symbols) versus centrality measure $\nu$ computed at fixed energy $(200 \mathrm{GeV})$. The SS $2 \mathrm{D}$ Gaussian amplitudes, $\eta_{\Delta}$ widths, and $\phi_{\Delta}$ widths are shown in the left, center, and right panels, respectively of the top row. The bottom row shows from left to right the amplitudes for the dipole, quadrupole, and SS peak width aspect ratio $\sigma_{\eta_{\Delta}} / \sigma_{\phi_{\Delta}} .$ Fitting errors are indicated by error bars where larger than the symbols. Solid lines connect the points for clarity. The dotted and dashed curves indicate Glauber linear superposition estimates for 62 - and 200 -GeV peak amplitudes respectively, as discussed in the text. The quadrupole data are consistent with Ref. [60]. The hatched regions indicate the full range of systematic uncertainties listed in Appendix F. The vertical dark bands indicate estimated $v$ equivalents for $N-N$ collisions and $b=0$ Au-Au collisions.

FIG. 3. Fit parameters for $\left(\eta_{\Delta}, \phi_{\Delta}\right)$ correlation data from Au-Au collisions at $\sqrt{s_{N N}}=62$ (open symbols) and 200 GeV (solid symbols) versus centrality measure $\nu$ computed at fixed energy $(200 \mathrm{GeV})$. The SS $2 \mathrm{D}$ Gaussian amplitudes, $\eta_{\Delta}$ widths, and $\phi_{\Delta}$ widths are shown in the left, center, and right panels, respectively of the top row. The bottom row shows from left to right the amplitudes for the dipole, quadrupole, and SS peak width aspect ratio $\sigma_{\eta_{\Delta}} / \sigma_{\phi_{\Delta}} .$ Fitting errors are indicated by error bars where larger than the symbols. Solid lines connect the points for clarity. The dotted and dashed curves indicate Glauber linear superposition estimates for 62 - and 200 -GeV peak amplitudes respectively, as discussed in the text. The quadrupole data are consistent with Ref. [60]. The hatched regions indicate the full range of systematic uncertainties listed in Appendix F. The vertical dark bands indicate estimated $v$ equivalents for $N-N$ collisions and $b=0$ Au-Au collisions.

FIG. 3. Fit parameters for $\left(\eta_{\Delta}, \phi_{\Delta}\right)$ correlation data from Au-Au collisions at $\sqrt{s_{N N}}=62$ (open symbols) and 200 GeV (solid symbols) versus centrality measure $\nu$ computed at fixed energy $(200 \mathrm{GeV})$. The SS $2 \mathrm{D}$ Gaussian amplitudes, $\eta_{\Delta}$ widths, and $\phi_{\Delta}$ widths are shown in the left, center, and right panels, respectively of the top row. The bottom row shows from left to right the amplitudes for the dipole, quadrupole, and SS peak width aspect ratio $\sigma_{\eta_{\Delta}} / \sigma_{\phi_{\Delta}} .$ Fitting errors are indicated by error bars where larger than the symbols. Solid lines connect the points for clarity. The dotted and dashed curves indicate Glauber linear superposition estimates for 62 - and 200 -GeV peak amplitudes respectively, as discussed in the text. The quadrupole data are consistent with Ref. [60]. The hatched regions indicate the full range of systematic uncertainties listed in Appendix F. The vertical dark bands indicate estimated $v$ equivalents for $N-N$ collisions and $b=0$ Au-Au collisions.

FIG. 3. Fit parameters for $\left(\eta_{\Delta}, \phi_{\Delta}\right)$ correlation data from Au-Au collisions at $\sqrt{s_{N N}}=62$ (open symbols) and 200 GeV (solid symbols) versus centrality measure $\nu$ computed at fixed energy $(200 \mathrm{GeV})$. The SS $2 \mathrm{D}$ Gaussian amplitudes, $\eta_{\Delta}$ widths, and $\phi_{\Delta}$ widths are shown in the left, center, and right panels, respectively of the top row. The bottom row shows from left to right the amplitudes for the dipole, quadrupole, and SS peak width aspect ratio $\sigma_{\eta_{\Delta}} / \sigma_{\phi_{\Delta}} .$ Fitting errors are indicated by error bars where larger than the symbols. Solid lines connect the points for clarity. The dotted and dashed curves indicate Glauber linear superposition estimates for 62 - and 200 -GeV peak amplitudes respectively, as discussed in the text. The quadrupole data are consistent with Ref. [60]. The hatched regions indicate the full range of systematic uncertainties listed in Appendix F. The vertical dark bands indicate estimated $v$ equivalents for $N-N$ collisions and $b=0$ Au-Au collisions.

FIG. 4. (Color online) (a) 2D data histogram from $0 \%-5 \%$ central 200 GeV Au-Au collisions $[(0,0)$ bin suppressed]. (b) Additional model component fitted to the $2 \mathrm{D}$ histogram and consisting of an AS azimuth dipole modulated by function $F\left(\eta_{\Delta}\right)$ described in the text. (c) Fit residuals including the new model component. (d) Data histogram in the first panel minus the additional model component and fitted offset. Remaining structure is described accurately by an AS dipole and SS 2D Gaussian.

FIG. 9. (a) Amplitude of SS 2D peak (jet-correlated pairs per final-state hadron) for data (solid points), GLS extrapolation from $p$ - $p$ data (dashed curve) and from the HIJING Monte Carlo (open points, dotted curve). (b) Single-particle production extrapolated from $p-p$ data (dashed curve, GLS), and from more-central Au-Au data (solid line) and HIJING (dash-dotted line, open points). (c) Single-particle hard-component production per $N-N$ binary collision from $p-p$ data, from more-central Au-Au data and from HIJING (open points). (d) Jet-correlated pair production (SS 2D peak) per $N-N$ binary collision extrapolated from $p-p$ data (GLS), and from Au-Au data (solid points) and HIJING (open points).

FIG. 9. (a) Amplitude of SS 2D peak (jet-correlated pairs per final-state hadron) for data (solid points), GLS extrapolation from $p$ - $p$ data (dashed curve) and from the HIJING Monte Carlo (open points, dotted curve). (b) Single-particle production extrapolated from $p-p$ data (dashed curve, GLS), and from more-central Au-Au data (solid line) and HIJING (dash-dotted line, open points). (c) Single-particle hard-component production per $N-N$ binary collision from $p-p$ data, from more-central Au-Au data and from HIJING (open points). (d) Jet-correlated pair production (SS 2D peak) per $N-N$ binary collision extrapolated from $p-p$ data (GLS), and from Au-Au data (solid points) and HIJING (open points).

FIG. 9. (a) Amplitude of SS 2D peak (jet-correlated pairs per final-state hadron) for data (solid points), GLS extrapolation from $p$ - $p$ data (dashed curve) and from the HIJING Monte Carlo (open points, dotted curve). (b) Single-particle production extrapolated from $p-p$ data (dashed curve, GLS), and from more-central Au-Au data (solid line) and HIJING (dash-dotted line, open points). (c) Single-particle hard-component production per $N-N$ binary collision from $p-p$ data, from more-central Au-Au data and from HIJING (open points). (d) Jet-correlated pair production (SS 2D peak) per $N-N$ binary collision extrapolated from $p-p$ data (GLS), and from Au-Au data (solid points) and HIJING (open points).

FIG. 9. (a) Amplitude of SS 2D peak (jet-correlated pairs per final-state hadron) for data (solid points), GLS extrapolation from $p$ - $p$ data (dashed curve) and from the HIJING Monte Carlo (open points, dotted curve). (b) Single-particle production extrapolated from $p-p$ data (dashed curve, GLS), and from more-central Au-Au data (solid line) and HIJING (dash-dotted line, open points). (c) Single-particle hard-component production per $N-N$ binary collision from $p-p$ data, from more-central Au-Au data and from HIJING (open points). (d) Jet-correlated pair production (SS 2D peak) per $N-N$ binary collision extrapolated from $p-p$ data (GLS), and from Au-Au data (solid points) and HIJING (open points).

FIG. 11. (Color online) (Left) Uncorrected angular correlations from 62-GeV $37 \%-46 \%$ central Au-Au collisions showing pileup distortions, especially evident as the W-shaped nonuniformity of the AS ridge on $\eta_{\Delta}$. (Right) The same data with pileup correction applied.

Elliptic flow versus $p_{T}$ for charged hadrons from Cu+Cu collisions at 0–60% centrality at $\sqrt{s_{NN}}$ = 22.4 GeV measured with the subevent method with a pseudorapidity gap of 0.3 units compared with previously published STAR results for 200 and 62.4 GeV Cu+Cu [27] measured with full TPC event plane method $v_{2}${TPC} and full FTPC event plane method $v_{2}${FTPC}. The error bars are statistical. Results are also compared to $v_{2}$($p_{T}$) model calculations.

Elliptic flow $v_{2}(\eta)$ for charged hadrons fromCu+Cu collisions at 0–60% centrality at $\sqrt{s_{NN}}$ = 22.4 GeV. The present STAR results are compared to the measurement from the PHOBOS [29] collaboration for Cu+Cu at 22.4 GeV. The PHOBOS results include statistical and systematic errors whereas the STAR results are plotted with statistical uncertainties only, and systematic errors are discussed in Section III C. Results are also compared to $v_{2}(\eta)$ calculations from the indicated models.

Elliptic flow $v_{2}${TPC} and $v_{2}${4} as a function of centrality for charged hadrons from Cu+Cu collisions at 0 − 60% centrality at $\sqrt{s_{NN}}$ = 22.4 GeV, compared with previously published results from the STAR collaboration at 200 GeV and 62.4 GeV.

Charged hadron azimuthal anisotropy v2 as a function of pT in 50-60% Cu+Cu collisions at 200 GeV using FTPC flow vectors, and that with subtracting the azimuthal correlations in p+p collisions.

Charged hadron azimuthal anisotropy v2 as a function of pT in 40-50% Cu+Cu collisions at 200 GeV using FTPC flow vectors, and that with subtracting the azimuthal correlations in p+p collisions.

Charged hadron azimuthal anisotropy v2 as a function of pT in 30-40% Cu+Cu collisions at 200 GeV using FTPC flow vectors, and that with subtracting the azimuthal correlations in p+p collisions.

Charged hadron azimuthal anisotropy v2 as a function of pT in 20-30% Cu+Cu collisions at 200 GeV using FTPC flow vectors, and that with subtracting the azimuthal correlations in p+p collisions.

Charged hadron azimuthal anisotropy v2 as a function of pT in 10-20% Cu+Cu collisions at 200 GeV using FTPC flow vectors, and that with subtracting the azimuthal correlations in p+p collisions.

Charged hadron azimuthal anisotropy v2 as a function of pT in 0-10% Cu+Cu collisions at 200 GeV using FTPC flow vectors, and that with subtracting the azimuthal correlations in p+p collisions.

Charged hadron azimuthal anisotropy v2 integrated over pT and eta vs centrality in Cu+Cu collisions at 200 GeV for the various methods.

Charged hadron azimuthal anisotropy v2 integrated over pT and eta vs centrality in Cu+Cu collisions at 62.4 GeV for the various methods.

Charged hadron azimuthal anisotropy v2 as a function of pT in 50-60% Cu+Cu collisions at 62.4 GeV. Data tables for 200 GeV are shown in Fig.4.

Charged hadron azimuthal anisotropy v2 as a function of pT in 40-50% Cu+Cu collisions at 62.4 GeV using FTPC flow vectors. Data tables for 200 GeV are shown in Fig.4.

Charged hadron azimuthal anisotropy v2 as a function of pT in 30-40% Cu+Cu collisions at 62.4 GeV using FTPC flow vectors. Data tables for 200 GeV are shown in Fig.4.

Charged hadron azimuthal anisotropy v2 as a function of pT in 20-30% Cu+Cu collisions at 62.4 GeV using FTPC flow vectors. Data tables for 200 GeV are shown in Fig.4.

Charged hadron azimuthal anisotropy v2 as a function of pT in 10-20% Cu+Cu collisions at 62.4 GeV using FTPC flow vectors. Data tables for 200 GeV are shown in Fig.4.

Charged hadron azimuthal anisotropy v2 as a function of pT in 0-10% Cu+Cu collisions at 62.4 GeV using FTPC flow vectors. Data tables for 200 GeV are shown in Fig.4.

Charged hadron azimuthal anisotropy v2 as a function of pT in 0-60% Cu+Cu collisions at 62.4 GeV using FTPC flow vectors. Data tables for 200 GeV are shown in Fig.3.

K0S azimuthal anisotropy v2 vs pT in 0-60% central Cu+Cu collisions at 200 GeV.

K0S azimuthal anisotropy v2 vs pT in 0-20% central Cu+Cu collisions at 200 GeV.

K0S azimuthal anisotropy v2 vs pT in 20-60% central Cu+Cu collisions at 200 GeV.

Lambda azimuthal anisotropy v2 vs pT in 0-60% central Cu+Cu collisions at 200 GeV.

Lambda azimuthal anisotropy v2 vs pT in 0-20% central Cu+Cu collisions at 200 GeV.

Lambda azimuthal anisotropy v2 vs pT in 20-60% central Cu+Cu collisions at 200 GeV.

Xi azimuthal anisotropy v2 vs pT in 0-60% central Cu+Cu collisions at 200 GeV.

phi azimuthal anisotropy v2 vs pT in 0-60% central Cu+Cu collisions at 200 GeV.

$\langle cos(\phi_{\alpha}+\phi_{\beta}−2\phi_{c})\rangle$ as a function of reference multiplicity for different charge combinations, after corrections for acceptance effects. In the legend the signs indicate the charge of particles $\alpha$, $\beta$, and c. The results shown are for Au+Au collisions at 200 GeV obtained in the Full Field.

$\langle cos(\phi_{\alpha}-\phi_{\beta})\rangle$ as a function of centrality for different charge combinations and FF and RFF configurations. The data points corresponding to different charge and field configurations are slightly shifted in the horizontal direction with respect to each other for clarity. The error bars are statistical. Also shown are model predictions described in Section VII.

A comparison of the correlations obtained by selecting the third particle from the main TPC or from the Forward TPCs. The TPC and FTPC points are shifted horizontally relative to one another for clarity purposes. The error bars are statistical.

A comparison of the correlations obtained by selecting the third particle from the main TPC or from the Forward TPCs after scaling by the flow of the third particle. The shaded areas represent the uncertainty from $v_{2,c}$ scaling (see text for details). The TPC and FTPC points are shifted horizontally relative to one another for clarity purposes. The error bars are statistical.

$\langle cos(\phi_{\alpha}+\phi_{\beta}−2\Psi_{RP})\rangle$ in Au+Au and Cu+Cu collisions at $\sqrt{s_{NN}}$ = 200 GeV calculated using Eq. 7. The error-bars show the statistical errors. The shaded area reflects the uncertainty in the elliptic flow values used in calculations, with lower (in magnitude) limit obtained with elliptic flow from two-particle correlations and upper limit from four-particle cumulants. For details, see Section IV. Thick solid (Au+Au) and dashed (Cu+Cu) lines represent possible non-reaction-plane dependent contribution from many-particle clusters as estimated by HIJING, see Section VII A.

$\langle cos(\phi_{\alpha}+\phi_{\beta}−2\Psi_{RP})\rangle$ in Au+Au and Cu+Cu collisions at $\sqrt{s_{NN}}$ = 62 GeV calculated using Eq. 7. The error-bars show the statistical errors. The shaded area reflects the uncertainty in the elliptic flow values used in calculations. For details, see Section IV. Thick solid (Au+Au) and dashed (Cu+Cu) lines represent possible non-reaction-plane dependent contribution from many-particle clusters as estimated by HIJING, see Section VII A.

Au+Au and Cu+Cu collisions at $\sqrt{s_{NN}}$ = 200 GeV. The correlations are scaled with the number of participants and are plotted as function of centrality. The error-bars indicate the statistical errors. The shaded area reflects the uncertainty in the elliptic flow values used in calculations. For details, see Section IV. Thick solid (Au+Au) and dashed (Cu+Cu) lines represent possible non-reaction-plane dependent contribution from many-particle clusters as estimated by HIJING, see Section VII A.

Au+Au and Cu+Cu collisions at $\sqrt{s_{NN}}$ = 200 GeV. The correlations are scaled with the number of participants and are plotted as function of number of participants. The error-bars indicate the statistical errors. The shaded area reflects the uncertainty in the elliptic flow values used in calculations. For details, see Section IV. Thick solid (Au+Au) and dashed (Cu+Cu) lines represent possible non-reaction-plane dependent contribution from many-particle clusters as estimated by HIJING, see Section VII A.

Au+Au at 200 GeV. The correlations dependence on pseudorapidity separation $\Delta\eta$ = |$\eta_{\alpha}$ − $\eta_{\beta}$| for centrality 30-50%. The shaded band indicates uncertainty associated with $v_{2}$ measurements and has been calculated using two- and four-particle cumulant results as the limits.

Au+Au at 200 GeV. The correlations dependence on pseudorapidity separation $\Delta\eta$ = |$\eta_{\alpha}$ − $\eta_{\beta}$| for centrality 10-30%. The shaded band indicates uncertainty associated with $v_{2}$ measurements and has been calculated using two- and four-particle cumulant results as the limits.

Au+Au at 200 GeV. The correlations dependence on (p$_{t,\alpha}$ + p$_{t,\beta}$)/2 for centrality 30-50%. The shaded band indicates uncertainty associated with $v_{2}$ measurements and has been calculated using two- and four-particle cumulant results as the limits.

Au+Au at 200 GeV. The correlations dependence on (p$_{t,\alpha}$ + p$_{t,\beta}$)/2 for centrality 10-30%. The shaded band indicates uncertainty associated with $v_{2}$ measurements and has been calculated using two- and four-particle cumulant results as the limits.

Au+Au at 200 GeV. The correlations dependence on |p$_{t,\alpha}$ - p$_{t,\beta}$| for centrality 30-50%. The shaded band indicates uncertainty associated with $v_{2}$ measurements and has been calculated using two- and four-particle cumulant results as the limits.

Au+Au at 200 GeV. The correlations dependence on |p$_{t,\alpha}$ - p$_{t,\beta}$| for centrality 10-30%. The shaded band indicates uncertainty associated with $v_{2}$ measurements and has been calculated using two- and four-particle cumulant results as the limits.

Three-particle correlator in Au+Au and Cu+Cu collisions compared to HIJING calculations shown as thick lines. All three particles are taken in the main TPC region, |$\eta$| < 1.0. The correlator has been scaled with number of participants and are plotted versus number of participants. In this representation HIJING results for Au+Au and Cu+Cu collisions coincide in the region of overlap.

$\langle cos(\phi_{\alpha} + \phi_{\beta} − 2\Psi_{RP})\rangle$ calculated for 200 GeV Au+Au events with event generators HIJING (with and without an “elliptic flow afterburner”), UrQMD, and MEVSIM. Blue symbols mark opposite-charge correlations, and red are same-charge. Solid stars represent the values from the data to facilitate comparison. Acceptance cuts of 0.15 < p$_{t}$ < 2 GeV/c and |$\eta$| < 1.0 were used in all cases. For MEVSIM, HIJING, and UrQMD points the true reaction plane from the generated event was used for $\Psi_{RP}$ . Thick solid lighter colored lines represent possible non-reactionplane dependent contribution from many-particle clusters as estimated by HIJING and discussed in Section VII A. Corresponding estimates from UrQMD are about factor of two smaller.

Negatively charged kaon spectra from Cu+Cu collisions 62.4 GeV as a function of pT for different centralities.

Negatively charged proton spectra from Cu+Cu collisions 200 GeV as a function of pT for different centralities.

Negatively charged proton spectra from Cu+Cu collisions 62.4 GeV as a function of pT for different centralities.

Positively charged pion spectra from Cu+Cu collisions 200 GeV as a function of pT for different centralities.

Positively charged pion spectra from Cu+Cu collisions 62.4 GeV as a function of pT for different centralities.

Positively charged kaon spectra from Cu+Cu collisions 200 GeV as a function of pT for different centralities.

Positively charged kaon spectra from Cu+Cu collisions 62.4 GeV as a function of pT for different centralities.

Positively charged proton spectra from Cu+Cu collisions 200 GeV as a function of pT for different centralities.

Positively charged proton spectra from Cu+Cu collisions 62.4 GeV as a function of pT for different centralities.

Mean transverse momentum of negatively charged pions, kaons and protons as a function of charged hadron multiplicity.

Mean transverse momentum of positively charged pions, kaons and protons as a function of charged hadron multiplicity [figure not available in published paper].

Integrated yields of negatively charged pions, kaons and protons as a function of charged hadron multiplicity.

Integrated yields of positively charged pions, kaons and protons as a function of charged hadron multiplicity [figure not available in published paper].

particle Ratios -I (pbar/pi^-, k^-/pi^-) versus multiplicity.

particle ratios -II (p/pi^+, k^+/pi^+) versus multiplicity.

particle ratios -III (p+pbar/pi, k/pi) versus multiplicity.

particle Ratios -IV (pi^-/pi^+, k^-/k^+, pbar/p) versus multiplicity [only pbar/p figure available in paper].

Enhancement factors for negatively charged pions, kaons and protons as a function of Npart [Ref. Phys.Rev.C 81, 044902, 2010]. pp dNdy values are from Ref [Phys.Rev.C 79, 034909, 2009].

The kinetic freeze-out temperature (Tkin) and chemical freeze-out temperature (Tch) versus multiplicity.

flow velocity versus multiplicity.

chemical potentials versus multiplicity.

strangeness suppression factor versus multiplicity.

dNch/deta values for different centrality.

The $v_2$ as a function of $p_T$ for 20-30% central Au + Au collisions at midrapidity for $\sqrt{s_{NN}}$ = 7.7 GeV, 11.5 GeV, 19.6 GeV, 27 GeV and 39 GeV.

$\varepsilon$ (Glauber) as a function of $p_T$ for various collision centralities (10-20%, 30-40% and 50-60%) in Au + Au collisions at midrapidity for $\sqrt{s_{NN}}$ = 7.7 GeV, 11.5 GeV, 19.6 GeV, 27 GeV and 39 GeV.

$\varepsilon$ (CGC) as a function of $p_T$ for various collision centralities (10-20%, 30-40% and 50-60%) in Au + Au collisions at midrapidity for $\sqrt{s_{NN}}$ = 7.7 GeV, 11.5 GeV, 19.6 GeV, 27 GeV and 39 GeV.

$v_2${EtaSubs} as a function of $p_T$ for various collision centralities (10-20%, 30-40% and 50-60%) in Au + Au collisions at midrapidity for $\sqrt{s_{NN}}$ = 7.7 GeV, 11.5 GeV, 19.6 GeV, 27 GeV and 39 GeV.

The $v_2${EP} vs. $\eta$ for 10-40% centrality in Au + Au collisions at $\sqrt{s_{NN}}$ = 7.7 GeV, 11.5 GeV, 19.6 GeV, 27 GeV and 39 GeV.

The $v_2${EP} vs. $\eta$ for 10-40% centrality in Au + Au collisions at $\sqrt{s_{NN}}$ = 62.4 GeV.

The $v_2${4} vs. $p_T$ at midrapidity for various collision energies ($\sqrt{s_{NN}}$ = 7.7 GeV, 11.5 GeV, 19.6 GeV, 27 GeV and 39 GeV).

The $v_2${4} vs. $p_T$ at midrapidity for $\sqrt{s_{NN}}$ = 62.4 GeV.

The $v_2${4} vs. $p_T$ at midrapidity for $\sqrt{s_{NN}}$ = 200 GeV.

Efficiency corrected pT spectra for the different centrality bins and for the various particles. Note that 7 centrality bins have been used for the KS0 and the Λ while only 6 and 3 have been used for the Ξ and Ω, respectively. Errors are statistical only. The Λ spectra are corrected for the feed-down of the Ξ decay.

Efficiency corrected pT spectra for the different centrality bins and for the various particles. Note that 7 centrality bins have been used for the KS0 and the Λ while only 6 and 3 have been used for the Ξ and Ω, respectively. Errors are statistical only. The Λ spectra are corrected for the feed-down of the Ξ decay.

Efficiency corrected pT spectra for the different centrality bins and for the various particles. Note that 7 centrality bins have been used for the KS0 and the Λ while only 6 and 3 have been used for the Ξ and Ω, respectively. Errors are statistical only. The Λ spectra are corrected for the feed-down of the Ξ decay.

Efficiency corrected pT spectra for the different centrality bins and for the various particles. Note that 7 centrality bins have been used for the KS0 and the Λ while only 6 and 3 have been used for the Ξ and Ω, respectively. Errors are statistical only. The Λ spectra are corrected for the feed-down of the Ξ decay.

Efficiency corrected pT spectra for the different centrality bins and for the various particles. Note that 7 centrality bins have been used for the KS0 and the Λ while only 6 and 3 have been used for the Ξ and Ω, respectively. Errors are statistical only. The Λ spectra are corrected for the feed-down of the Ξ decay.

Efficiency corrected pT spectra for the different centrality bins and for the various particles. Note that 7 centrality bins have been used for the KS0 and the Λ while only 6 and 3 have been used for the Ξ and Ω, respectively. Errors are statistical only. The Λ spectra are corrected for the feed-down of the Ξ decay.

Efficiency corrected pT spectra for the different centrality bins and for the various particles. Note that 7 centrality bins have been used for the KS0 and the Λ while only 6 and 3 have been used for the Ξ and Ω, respectively. Errors are statistical only. The Λ spectra are corrected for the feed-down of the Ξ decay.

Efficiency corrected pT spectra for the different centrality bins and for the various particles. Note that 7 centrality bins have been used for the KS0 and the Λ while only 6 and 3 have been used for the Ξ and Ω, respectively. Errors are statistical only. The Λ spectra are corrected for the feed-down of the Ξ decay.

Efficiency corrected pT spectra for the different centrality bins and for the various particles. Note that 7 centrality bins have been used for the KS0 and the Λ while only 6 and 3 have been used for the Ξ and Ω, respectively. Errors are statistical only. The Λ spectra are corrected for the feed-down of the Ξ decay.

Efficiency corrected pT spectra for the different centrality bins and for the various particles. Note that 7 centrality bins have been used for the KS0 and the Λ while only 6 and 3 have been used for the Ξ and Ω, respectively. Errors are statistical only. The Λ spectra are corrected for the feed-down of the Ξ decay.

Efficiency corrected pT spectra for the different centrality bins and for the various particles. Note that 7 centrality bins have been used for the KS0 and the Λ while only 6 and 3 have been used for the Ξ and Ω, respectively. Errors are statistical only. The Λ spectra are corrected for the feed-down of the Ξ decay.

Efficiency corrected pT spectra for the different centrality bins and for the various particles. Note that 7 centrality bins have been used for the KS0 and the Λ while only 6 and 3 have been used for the Ξ and Ω, respectively. Errors are statistical only. The Λ spectra are corrected for the feed-down of the Ξ decay.

Efficiency corrected pT spectra for the different centrality bins and for the various particles. Note that 7 centrality bins have been used for the KS0 and the Λ while only 6 and 3 have been used for the Ξ and Ω, respectively. Errors are statistical only. The Λ spectra are corrected for the feed-down of the Ξ decay.

Extrapolated average transverse momenta ⟨pT ⟩ as a function of dNch/dy for different particle species in Au+Au collisions at 62.4 GeV. Statistical uncertainties are represented by the error bars at the points while the systematic uncertainties are represented by the gray bars. The π, charged K and p data were extracted from Ref. [14].

Extrapolated average transverse momenta ⟨pT ⟩ as a function of dNch/dy for different particle species in Au+Au collisions at 62.4 GeV. Statistical uncertainties are represented by the error bars at the points while the systematic uncertainties are represented by the gray bars. The π, charged K and p data were extracted from Ref. [14].

KS0 dN/dpT spectra compared to the charged Kaon spectra for the event centrality of 0-5% and 30-40%. The charged Kaons data points are for rapidity range of |y| < 0.1 and were extracted from Ref. [14].

KS0 dN/dpT spectra compared to the charged Kaon spectra for the event centrality of 0-5% and 30-40%. The charged Kaons data points are for rapidity range of |y| < 0.1 and were extracted from Ref. [14].

KS0 dN/dpT spectra compared to the charged Kaon spectra for the event centrality of 0-5% and 30-40%. The charged Kaons data points are for rapidity range of |y| < 0.1 and were extracted from Ref. [14].

KS0 dN/dpT spectra compared to the charged Kaon spectra for the event centrality of 0-5% and 30-40%. The charged Kaons data points are for rapidity range of |y| < 0.1 and were extracted from Ref. [14].

Strange particle production yields at mid-rapidity in central Au+Au and Pb+Pb collisions versus the center of mass energy √sNN. The top panel shows results for K0S and Λ. The AGS values are from E896 [1] (centrality 0 − 5 %). The SPS values are from NA49 [20] (centrality 0 − 7 %) and the RHIC values are from STAR [4, 15] (centrality 0 − 5 %). For the multi-strange baryons Ξ and Ω (bottom panel), the SPS results are from NA57 [2] (centrality 0 − 11 %) and the RHIC values are from STAR [15, 21] (centrality 0 − 20 %).

Strange particle production yields at mid-rapidity in central Au+Au and Pb+Pb collisions versus the center of mass energy √sNN. The top panel shows results for K0S and Λ. The AGS values are from E896 [1] (centrality 0 − 5 %). The SPS values are from NA49 [20] (centrality 0 − 7 %) and the RHIC values are from STAR [4, 15] (centrality 0 − 5 %). For the multi-strange baryons Ξ and Ω (bottom panel), the SPS results are from NA57 [2] (centrality 0 − 11 %) and the RHIC values are from STAR [15, 21] (centrality 0 − 20 %).

Strange particle production yields at mid-rapidity in central Au+Au and Pb+Pb collisions versus the center of mass energy √sNN. The top panel shows results for K0S and Λ. The AGS values are from E896 [1] (centrality 0 − 5 %). The SPS values are from NA49 [20] (centrality 0 − 7 %) and the RHIC values are from STAR [4, 15] (centrality 0 − 5 %). For the multi-strange baryons Ξ and Ω (bottom panel), the SPS results are from NA57 [2] (centrality 0 − 11 %) and the RHIC values are from STAR [15, 21] (centrality 0 − 20 %).

Strange particle production yields at mid-rapidity in central Au+Au and Pb+Pb collisions versus the center of mass energy √sNN. The top panel shows results for K0S and Λ. The AGS values are from E896 [1] (centrality 0 − 5 %). The SPS values are from NA49 [20] (centrality 0 − 7 %) and the RHIC values are from STAR [4, 15] (centrality 0 − 5 %). For the multi-strange baryons Ξ and Ω (bottom panel), the SPS results are from NA57 [2] (centrality 0 − 11 %) and the RHIC values are from STAR [15, 21] (centrality 0 − 20 %).

Anti-baryon to baryon yield ratios for strange baryons versus the center of mass energy √sNN. Λ/Λ is shown in the top panel while the multi-strange baryons are on the bottom panel. The data from AGS are not corrected for the weak decay feed-down from the multistrange baryons while the data from SPS and RHIC are corrected. The lines are the results of a thermal model calculation (see text section IV A). The AGS values are from E896 [1] (centrality 0 − 5 %). The SPS values are from NA49 [20] (centrality 0 − 7 %) and the RHIC values are from STAR [4, 15] (centrality 0 − 5 %). For the multi- strange baryons Ξ and Ω (bottom panel), the SPS results are from NA57 [2] (centrality 0 − 11 %) and the RHIC values are from STAR [15, 21] (centrality 0 − 20 %).

Anti-baryon to baryon yield ratios for strange baryons versus the center of mass energy √sNN. Λ/Λ is shown in the top panel while the multi-strange baryons are on the bottom panel. The data from AGS are not corrected for the weak decay feed-down from the multistrange baryons while the data from SPS and RHIC are corrected. The lines are the results of a thermal model calculation (see text section IV A). The AGS values are from E896 [1] (centrality 0 − 5 %). The SPS values are from NA49 [20] (centrality 0 − 7 %) and the RHIC values are from STAR [4, 15] (centrality 0 − 5 %). For the multi- strange baryons Ξ and Ω (bottom panel), the SPS results are from NA57 [2] (centrality 0 − 11 %) and the RHIC values are from STAR [15, 21] (centrality 0 − 20 %).

Antibaryon-to-baryon yield ratios for strange particles and protons as a function of dNch/dy at √sNN=62.4 and 200 GeV. The p data were extracted from Ref. [14]. The √sNN=200 GeV strange hadron data were extracted from Ref. [15].

Antibaryon-to-baryon yield ratios for strange particles and protons as a function of dNch/dy at √sNN=62.4 and 200 GeV. The p data were extracted from Ref. [14]. The √sNN=200 GeV strange hadron data were extracted from Ref. [15].

Particle-yield ratios as obtained by measurements (black dots) for the most central (0–5%) Au+Au collisions at 62.4 GeV and statistical model predictions (lines). The ratios indicated by the dashed lines (blue) were obtained by using only π, K, and protons, whereas the ratios indicated by the full lines (green) were obtained by also using the hyperons in the fit.

Particle-yield ratios as obtained by measurements (black dots) for the most central (0–5%) Au+Au collisions at 62.4 GeV and statistical model predictions (lines). The ratios indicated by the dashed lines (blue) were obtained by using only π, K, and protons, whereas the ratios indicated by the full lines (green) were obtained by also using the hyperons in the fit.

Chemical freeze-out temperature Tch (a) and strangeness saturation factor γs (b) as a function of the mean number of participants.

Chemical freeze-out temperature Tch (a) and strangeness saturation factor γs (b) as a function of the mean number of participants.

Chemical freeze-out temperature Tch (a) and strangeness saturation factor γs (b) as a function of the mean number of participants.

Chemical freeze-out temperature Tch (a) and strangeness saturation factor γs (b) as a function of the mean number of participants.

Temperature and baryon chemical potential obtained from thermal model fits as a function of √sNN (see Ref. [22]). The dashed lines correspond to the parametrizations given in Ref. [22]. The solid stars show the result for √sNN=62.4 and 200 GeV.

Temperature and baryon chemical potential obtained from thermal model fits as a function of √sNN (see Ref. [22]). The dashed lines correspond to the parametrizations given in Ref. [22]. The solid stars show the result for √sNN=62.4 and 200 GeV.

Temperature and baryon chemical potential obtained from thermal model fits as a function of √sNN (see Ref. [22]). The dashed lines correspond to the parametrizations given in Ref. [22]. The solid stars show the result for √sNN=62.4 and 200 GeV.

Temperature and baryon chemical potential obtained from thermal model fits as a function of √sNN (see Ref. [22]). The dashed lines correspond to the parametrizations given in Ref. [22]. The solid stars show the result for √sNN=62.4 and 200 GeV.

Ratio of baryon (solid symbols) and antibaryon (open symbols) to π+ as a function of dNch/dy for √sNN=62.4 GeV (left) and √sNN=200 GeV (right). The π and p data were extracted from Ref. [14].

Ratio of baryon (solid symbols) and antibaryon (open symbols) to π- as a function of dNch/dy for √sNN=62.4 GeV (left) and √sNN=200 GeV (right). The π and p data were extracted from Ref. [14].

Ratio of baryon (solid symbols) and antibaryon (open symbols) to π+ as a function of dNch/dy for √sNN=62.4 GeV (left) and √sNN=200 GeV (right). The π and p data were extracted from Ref. [14].

Ratio of baryon (solid symbols) and antibaryon (open symbols) to π- as a function of dNch/dy for √sNN=62.4 GeV (left) and √sNN=200 GeV (right). The π and p data were extracted from Ref. [14].

Ratio of baryon (solid symbols) and antibaryon (open symbols) to π+ as a function of dNch/dy for √sNN=62.4 GeV (left) and √sNN=200 GeV (right). The π and p data were extracted from Ref. [14].

Ratio of baryon (solid symbols) and antibaryon (open symbols) to π- as a function of dNch/dy for √sNN=62.4 GeV (left) and √sNN=200 GeV (right). The π and p data were extracted from Ref. [14].

Ratio of baryon (solid symbols) and antibaryon (open symbols) to π+ as a function of dNch/dy for √sNN=62.4 GeV (left) and √sNN=200 GeV (right). The π and p data were extracted from Ref. [14].

Ratio of baryon (solid symbols) and antibaryon (open symbols) to π- as a function of dNch/dy for √sNN=62.4 GeV (left) and √sNN=200 GeV (right). The π and p data were extracted from Ref. [14].

Ratio of baryon (solid symbols) and antibaryon (open symbols) to π− as a function of √sNN. The lines are the results of the thermal model calculation (see text Sec. 4a). The SPS values are from NA49 [20] (centrality 0–7%) and the RHIC values are from STAR [4, 15] (centrality 0–5%). For the multistrange baryons Ξ and Ω (bottom), the SPS results are from NA57 [2] (centrality 0–11%) and the RHIC values are from STAR [15, 21] (centrality 0–20%).

Ratio of baryon (solid symbols) and antibaryon (open symbols) to π− as a function of √sNN. The lines are the results of the thermal model calculation (see text Sec. 4a). The SPS values are from NA49 [20] (centrality 0–7%) and the RHIC values are from STAR [4, 15] (centrality 0–5%). For the multistrange baryons Ξ and Ω (bottom), the SPS results are from NA57 [2] (centrality 0–11%) and the RHIC values are from STAR [15, 21] (centrality 0–20%).

Nuclear modification factor RCP, calculated as the ratio between 0–10% central spectra and 40–80% peripheral spectra, for π, K0S, Λ, and Ξ particles in Au+Au collisions at 62.4 GeV. The π RCP values were extracted from Ref. [10]. The gray band on the right side of the plot shows the uncertainties on the estimation of the number of binary collisions and the gray band on the lower left side indicates the uncertainties on the number of participants.

Nuclear modification factor RCP, calculated as the ratio between 0–10% central spectra and 40–80% peripheral spectra, for π, K0S, Λ, and Ξ particles in Au+Au collisions at 62.4 GeV. The π RCP values were extracted from Ref. [10]. The gray band on the right side of the plot shows the uncertainties on the estimation of the number of binary collisions and the gray band on the lower left side indicates the uncertainties on the number of participants.

Nuclear modification factor RCP, calculated as the ratio between 0–5% central spectra and 40–60% peripheral spectra, for Λ and Ξ particles measured in Au + Au collisions at 62.4 GeV. The gray band corresponds to the equivalent RCP curve for the Λ particles measured in Au+Au collisions at 200 GeV [15].

Nuclear modification factor RCP, calculated as the ratio between 0–5% central spectra and 40–60% peripheral spectra, for Λ and Ξ particles measured in Au + Au collisions at 62.4 GeV. The gray band corresponds to the equivalent RCP curve for the Λ particles measured in Au+Au collisions at 200 GeV [15].

Λ/K0S ratio as a function of transverse momentum for different centrality classes. 0–5% (solid circles), 40–60% (open squares), and 60–80% (solid triangles) in Au+Au collisions at 62.4 GeV.

Λ/K0S ratio as a function of transverse momentum for different centrality classes. 0–5% (solid circles), 40–60% (open squares), and 60–80% (solid triangles) in Au+Au collisions at 62.4 GeV.

Maximum value of the Λ/K0S ratio from Au+Au collisions at 62.4 GeV (solid circles) and 200 GeV (open circles) [11] as a function of ⟨Npart⟩ for different centrality classes. The lowest ⟨Npart⟩ point corresponds to p+p collisions at 200 GeV [44]. The maximum of the Λ––/K0S from Au+Au collisions at 62.4 GeV is shown as solid triangles.

Maximum value of the Λ/K0S ratio from Au+Au collisions at 62.4 GeV (solid circles) and 200 GeV (open circles) [11] as a function of ⟨Npart⟩ for different centrality classes. The lowest ⟨Npart⟩ point corresponds to p+p collisions at 200 GeV [44]. The maximum of the Λ––/K0S from Au+Au collisions at 62.4 GeV is shown as solid triangles.

The K$\pi$ pair invariant mass distribution integrated over the $K^{*0}$ $p_T$ for minimum bias Cu+Cu collisions at $\sqrt{s_{NN}}$ =62.4 GeV after mixed-event background subtraction.

The Kπ pair invariant mass distribution for various pT bins (top left) pT = 0.4–0.6 GeV/c in Au+Au collisions at √sNN = 200 GeV after the mixed-event background subtraction.

The Kπ pair invariant mass distribution for various pT bins (top right) pT = 0.6–0.8 GeV/c in Au+Au collisions at √sNN = 62.4 GeV after the mixed-event background subtraction.

The Kπ pair invariant mass distribution for various pT bins (bottom left) pT = 0.8–1.0 GeV/c in Au+Au collisions at √sNN = 200 GeV after the mixed-event background subtraction.

The Kπ pair invariant mass distribution for various pT bins (bottom right) pT = 1.0–1.2 GeV/c in Au+Au collisions at √sNN = 62.4 GeV after the mixed-event background subtraction.

The signal-to-background ratio for $K^{*0}$ measurements as a function of $p_T$ for different collision centrality bins (0-10%, 10-40%, 40-60%, 60-80%) in Au+Au collisions at 200 GeV.

$K^{*0}$ mass as a function of $p_T$ for minimum bias Au+Au collisions at $\sqrt{s_{NN}}$ = 62.4 GeV.

$K^{*0}$ mass as a function of $p_T$ for minimum bias Au+Au collisions at $\sqrt{s_{NN}}$ = 200 GeV.

$K^{*0}$ mass as a function of $p_T$ for minimum bias Cu+Cu collisions at $\sqrt{s_{NN}}$ = 62.4 GeV

$K^{*0}$ mass as a function of $p_T$ for minimum bias Cu+Cu collisions at $\sqrt{s_{NN}}$ = 200 GeV

$K^{*0}$ width as a function of $p_T$ for minimum bias Au+Au collisions at $\sqrt{s_{NN}}$ = 62.4 GeV

$K^{*0}$ width as a function of $p_T$ for minimum bias Au+Au collisions at $\sqrt{s_{NN}}$ = 200 GeV

$K^{*0}$ width as a function of $p_T$ for minimum bias Cu+Cu collisions at $\sqrt{s_{NN}}$ = 62.4 GeV

$K^{*0}$ width as a function of $p_T$ for minimum bias Cu+Cu collisions at $\sqrt{s_{NN}}$ = 200 GeV

The $K^{*0}$ reconstruction efficiency multiplied by the detector acceptance as a function of $p_T$ in Au+Au (|$\eta$| < 0.8) collisions at 200 GeV for different collision centrality bins (0-20% ,20-40% , 40-60%)

The $K^{*0}$ reconstruction efficiency multiplied by the detector acceptance as a function of $p_T$ in Cu+Cu (|$\eta$| < 1.0) collisions at 200 GeV for different collision centrality bins (0-20% ,20-40% , 40-60%)

Mid-rapidity $K^{*0}$ $p_T$ spectra for various collision centrality bins (0-20%, 20-40%, 40-60%, 60-80%) in Au+Au collisions at $\sqrt{s_{NN}}$ = 62.4 GeV

Mid-rapidity $K^{*0}$ $p_T$ spectra for various collision centrality bins (0-20%, 20-40%, 40-60%) in Cu+Cu collisions at $\sqrt{s_{NN}}$ = 62.4 GeV

Mid-rapidity $K^{*0}$ $p_T$ spectra for various collision centrality bins (0-20%, 20-40%, 40-60%, 60-80%) in Au+Au collisions at $\sqrt{s_{NN}}$ = 200 GeV

Mid-rapidity $K^{*0}$ $p_T$ spectra for various collision centrality bins (0-20%, 20-40%, 40-60%) in Cu+Cu collisions at $\sqrt{s_{NN}}$ = 200 GeV

The mid-rapidity yields dN/dy of $K^{*0}$ as a function of the average number of participating nucleons, $⟨N_{part}⟩$, for Au+Au collisions at $\sqrt{s_{NN}}$ = 62.4 GeV

The mid-rapidity yields dN/dy of $K^{*0}$ as a function of the average number of participating nucleons, $⟨N_{part}⟩$, for Cu+Cu collisions at $\sqrt{s_{NN}}$ = 62.4 GeV

The mid-rapidity yields dN/dy of $K^{*0}$ as a function of the average number of participating nucleons, $⟨N_{part}⟩$, for Au+Au collisions at $\sqrt{s_{NN}}$ = 200 GeV

The mid-rapidity yields dN/dy of $K^{*0}$ as a function of the average number of participating nucleons, $⟨N_{part}⟩$, for Cu+Cu collisions at $\sqrt{s_{NN}}$ = 200 GeV

The mid-rapidity $K^{*0}$ $⟨p_T⟩$ as a function $⟨N_{part}⟩$ for Au+Au collisions at $\sqrt{s_{NN}}$ = 62.4 GeV

The mid-rapidity $K^{*0}$ $⟨p_T⟩$ as a function $⟨N_{part}⟩$ for Cu+Cu collisions at $\sqrt{s_{NN}}$ = 62.4 GeV

The mid-rapidity $K^{*0}$ $⟨p_T⟩$ as a function $⟨N_{part}⟩$ for Au+Au collisions at $\sqrt{s_{NN}}$ = 200 GeV

The mid-rapidity $K^{*0}$ $⟨p_T⟩$ as a function $⟨N_{part}⟩$ for Cu+Cu collisions at $\sqrt{s_{NN}}$ = 200 GeV

The mid-rapidity $⟨p_T⟩$ of $\pi$, K, p and $K^{*0}$ as a function of $⟨N_{part}⟩$ for Au+Au collisions at $\sqrt{s_{NN}}$ = 62.4 GeV.

Mid-rapidity $N(K^{*0})/N(K^-)$ ratio for Au+Au at $\sqrt{s_{NN}}$ = 62.4 GeV as a function of $⟨N_{part}⟩$

Mid-rapidity $N(K^{*0})/N(K^-)$ ratio for Cu+Cu at $\sqrt{s_{NN}}$ = 62.4 GeV as a function of $⟨N_{part}⟩$

Mid-rapidity $N(K^{*0})/N(K^-)$ ratio for Au+Au at $\sqrt{s_{NN}}$ = 200 GeV as a function of $⟨N_{part}⟩$

Mid-rapidity $N(K^{*0})/N(K^-)$ ratio for Cu+Cu at $\sqrt{s_{NN}}$ = 200 GeV as a function of $⟨N_{part}⟩$

Mid-rapidity $N(K^{*0})N(K^-)$ in Au+Au collisions divided by $N(K^{*0})N(K^-)$ ratio in p+p collisions at $\sqrt{s_{NN}}$=200 GeV as a function of $⟨N_{part}⟩$.

Mid-rapidity $N(K^{*0})N(K^-)$ in Cu+Cu collisions divided by $N(K^{*0})N(K^-)$ ratio in p+p collisions at $\sqrt{s_{NN}}$=200 GeV as a function of $⟨N_{part}⟩$

Mid-rapidity $N(K^{*0})N(K^-)$ in d+Au collisions divided by $N(K^{*0})N(K^-)$ ratio in d+Au collisions at $\sqrt{s_{NN}}$=200 GeV as a function of $⟨N_{part}⟩$

Mid-rapidity $N(K^{*0})/N(K^-)$ ratio in minimum bias Au+Au collisions as a function of $\sqrt{s_{NN}}.

Mid-rapidity $N(K^{*0})/N(K^-)$ ratio in minimum bias Cu+Cu collisions as a function of $\sqrt{s_{NN}}.

Mid-rapidity $N(K^{*0})/N(K^-)$ ratio in minimum bias p+p collisions as a function of $\sqrt{s_{NN}}.

Mid-rapidity $N(K^{*0})/N(K^-)$ ratio in minimum bias Au+Au collisions as a function of $\sqrt{s_{NN}}.

Mid-rapidity $N(K^{*0})/N(K^-)$ ratio in minimum bias Cu+Cu collisions as a function of $\sqrt{s_{NN}}.

Mid-rapidity $N(K^{*0})/N(K^-)$ ratio in minimum bias p+p collisions as a function of $\sqrt{s_{NN}}.

Mid-rapidity $N(\phi)/N(K^{*0})$ ratio for Au+Au at $\sqrt{s_{NN}}$ = 62.4 GeV as a function of $⟨N_{part}⟩$

Mid-rapidity $N(\phi)/N(K^{*0})$ ratio for Cu+Cu at $\sqrt{s_{NN}}$ = 62.4 GeV as a function of $⟨N_{part}⟩$

Mid-rapidity $N(\phi)/N(K^{*0})$ ratio for Au+Au at $\sqrt{s_{NN}}$ = 200 GeV as a function of $⟨N_{part}⟩$

Mid-rapidity $N(\phi)/N(K^{*0})$ ratio for Cu+Cu at $\sqrt{s_{NN}}$ = 200 GeV as a function of $⟨N_{part}⟩$

Mid-rapidity $[N(\phi)/N(K^{*0})]$ in Au+Au collisions divided by $[N(\phi)/N(K^{*0})]$ ratio in p+p collisions at $\sqrt{s_{NN}}$=200 GeV as a function of $⟨N_{part}⟩$

Mid-rapidity $[N(\phi)/N(K^{*0})]$ in Cu+Cu collisions divided by $[N(\phi)/N(K^{*0})]$ ratio in p+p collisions at $\sqrt{s_{NN}}$=200 GeV as a function of $⟨N_{part}⟩$

Mid-rapidity $[N(\phi)/N(K^{*0})]$ in d+Au collisions divided by $[N(\phi)/N(K^{*0})]$ ratio in p+p collisions at $\sqrt{s_{NN}}$=200 GeV as a function of $⟨N_{part}⟩$

Mid-rapidity $N(\phi)/N(K^{*0})$ ratio in minimum bias Au+Au collisions as a function of $\sqrt{s_{NN}}$.

Mid-rapidity $N(\phi)/N(K^{*0})$ ratio in minimum bias Cu+Cu collisions as a function of $\sqrt{s_{NN}}$.

Mid-rapidity $N(\phi)/N(K^{*0})$ ratio in minimum bias p+p collisions as a function of $\sqrt{s_{NN}}$.

Mid-rapidity $N(\phi)/N(K^{*0})$ ratio in minimum bias Au+Au collisions as a function of $\sqrt{s_{NN}}$.

Mid-rapidity $N(\phi)/N(K^{*0})$ ratio in minimum bias Cu+Cu collisions as a function of $\sqrt{s_{NN}}$.

Mid-rapidity $N(\phi)/N(K^{*0})$ ratio in minimum bias p+p collisions as a function of $\sqrt{s_{NN}}$.

The $K^{*0}$ $v_2$ (Run IV) as a function of $p_T$ in minimum bias Au+Au collisions at $\sqrt{s_{NN}}$ = 200 GeV.

The $K^{*0}$ $v_2$ (Run II) as a function of $p_T$ in minimum bias Au+Au collisions at $\sqrt{s_{NN}}$ = 200 GeV.

The $K^{*0}$ $R_{CP}$ as a function of $p_T$ in Au+Au collisions at 62.4 and 200 GeV compared to the $R_{CP}$ of $K^0_S$ and $\Lambda$ at 200 GeV.

The $K^{*0}$ $R_{CP}$ as a function of $p_T$ in Au+Au collisions at 62.4 and 200 GeV compared to the $R_{CP}$ of $K^0_S$ and $\Lambda$ at 200 GeV.

The $K^{*0}$ $R_{CP}$ as a function of $p_T$ in Au+Au collisions at 62.4 and 200 GeV compared to the $R_{CP}$ of $K^0_S$ and $\Lambda$ at 200 GeV.

The $K^{*0}$ ~$R_{CP}$~ as a function of $p_T$ in Au+Au collisions at 62.4 and 200 GeV compared to the $R_{CP}$ of $K^0_S$ and $\Lambda$ at 200 GeV.