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

The $\pi^0$ invariant mass peak positions and widths as a function of $p_{T}$ in 0-80% Au+Au collisions at $\sqrt{s_{NN}}=200$ GeV.

The $\pi^0$ invariant mass peak positions and widths as a function of $p_{T}$ in 0-80% Au+Au collisions at $\sqrt{s_{NN}}=200$ GeV.

The $\pi^0$ invariant mass peak positions and widths as a function of $p_{T}$ in 0-80% Au+Au collisions at $\sqrt{s_{NN}}=200$ GeV.

The $\pi^0$ invariant mass peak positions and widths as a function of $p_{T}$ in 0-80% Au+Au collisions at $\sqrt{s_{NN}}=200$ GeV.

The $\pi^0$ invariant mass peak positions and widths as a function of $p_{T}$ in 0-80% Au+Au collisions at $\sqrt{s_{NN}}=200$ GeV.

The $\pi^0$ invariant mass peak positions and widths as a function of $p_{T}$ in 0-80% Au+Au collisions at $\sqrt{s_{NN}}=200$ GeV.

Photon conversion point radius distribution from real data and MC simulation in Au+Au collisions at $\sqrt{s_{NN}}=200$ GeV.

Conversion probability correction factor for $\pi^0$ as a function of $p_{T}$

The overall detection efficiency of $\pi^0$ from embedding studt in 0-80% Au+Au collisions at $\sqrt{s_{NN}}=200$ GeV.

The overall detection efficiency of $\pi^0$ from embedding studt in 0-80% Au+Au collisions at $\sqrt{s_{NN}}=200$ GeV.

The overall detection efficiency of $\pi^0$ from embedding studt in 0-80% Au+Au collisions at $\sqrt{s_{NN}}=200$ GeV.

Invariant yield of STAR $\pi^0$ as a function of $p_{T}$ at midrapidity for different collision centrality bins in Au+Au collisions at $\sqrt{s_{NN}}=200$ GeV.

Invariant yield of STAR $\pi^0$ as a function of $p_{T}$ at midrapidity for different collision centrality bins in Au+Au collisions at $\sqrt{s_{NN}}=200$ GeV.

Invariant yield of STAR $\pi^0$ as a function of $p_{T}$ at midrapidity for different collision centrality bins in Au+Au collisions at $\sqrt{s_{NN}}=200$ GeV.

The ratios of STAR $\pi^0$ spectra over $\pi^{\pm}$ from STAR and $\pi^0$ from PHENIX for different collision centrality bins in Au+Au collisions at $\sqrt{s_{NN}}=200$ GeV.

The nuclear modification factor $R_{CP}$ as a functon of $p_{T}$ of STAR $\pi^0$ in Au+Au collisions at $\sqrt{s_{NN}}=200$ GeV.

The nuclear modification factor $R_{AA}$ as a functon of $p_{T}$ of STAR $\pi^0$ in Au+Au collisions at $\sqrt{s_{NN}}=200$ GeV.

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.

$e^+$ $e^-$ invariant mass pair distributions of signal pairs compared to the inclusive unlike sign and reconstructed background pairs in 200 GeV Au+Au minimum-bias collisions.

Signal-to-background ratios in minimum bias p+p and Au+Au collisions.

Signal-to-background ratios in central p+p and Au+Au collisions.

Invariant mass spectrum in the STAR acceptance ($p_T^e$ > 0.2 GeV/c, |$\eta^e$| < 1, and |$y_{ee}$|) from SQRT($s_{NN}$) = 200 GeV Au+Au minimum-bias collisions.

Invariant mass spectrum in the STAR acceptance ($p_T^e$ > 0.2 GeV/c, |$\eta^e$| < 1, and |$y_{ee}$|) from SQRT($s_{NN}$) = 200 GeV Au+Au minimum-bias collisions.

Invariant mass spectra from SQRT($s_{NN}$ = 200 GeV Au+Au minimum-bias collisions in pT range 0 - 5 GeV/c. The ratio of dielectron yield over cocktail for different pT bins, and the comparison with model calculations.

Invariant mass spectra from SQRT($s_{NN}$ = 200 GeV Au+Au minimum-bias collisions in pT range 5 - 10 GeV/c. The ratio of dielectron yield over cocktail for different pT bins, and the comparison with model calculations.

Invariant mass spectra from SQRT($s_{NN}$ = 200 GeV Au+Au minimum-bias collisions in pT range 1 - 1.5 GeV/c. The ratio of dielectron yield over cocktail for different pT bins, and the comparison with model calculations.

Invariant mass spectra from SQRT($s_{NN}$ = 200 GeV Au+Au minimum-bias collisions in pT range 1.5 - 2 GeV/c. The ratio of dielectron yield over cocktail for different pT bins, and the comparison with model calculations.

Invariant mass spectra from SQRT($s_{NN}$ = 200 GeV Au+Au minimum-bias collisions integrated. The ratio of dielectron yield over cocktail for different pT bins, and the comparison with model calculations.

The integrated dielectron yield as a function of pair pT in invariant mass range 0 - 0.15 GeV/c^2 compared with cocktail. The ratio of dielectron yield over cocktail for different mass ranges as a function of pair pT.

The integrated dielectron yield as a function of pair pT in invariant mass range 0.15 - 0.3 GeV/c^2 compared with cocktail. The ratio of dielectron yield over cocktail for different mass ranges as a function of pair pT.

The integrated dielectron yield as a function of pair pT in invariant mass range 0.3 - 0.76 GeV/c^2 compared with cocktail. The ratio of dielectron yield over cocktail for different mass ranges as a function of pair pT.

The integrated dielectron yield as a function of pair pT in invariant mass range 0.76 - 1.05 GeV/c^2 compared with cocktail. The ratio of dielectron yield over cocktail for different mass ranges as a function of pair pT.

The integrated dielectron yield as a function of pair pT in invariant mass range 1.05 - 1.8 GeV/c^2 compared with cocktail. The ratio of dielectron yield over cocktail for different mass ranges as a function of pair pT.

The integrated dielectron yield as a function of pair pT in invariant mass range 1.8 - 2.8 GeV/c^2 compared with cocktail. The ratio of dielectron yield over cocktail for different mass ranges as a function of pair pT.

The integrated dielectron yield as a function of pair pT in invariant mass range 2.8 to 3.5 GeV/c^2 compared with cocktail. The ratio of dielectron yield over cocktail for different mass ranges as a function of pair pT.

Invariant mass spectra from SQRT($s_{NN}$) = 200 GeV Au+Au collisions in different centralities. The ratio of dielectron yield over cocktail for a centrality of 0-10%.

Invariant mass spectra from SQRT($s_{NN}$) = 200 GeV Au+Au collisions in different centralities. The ratio of dielectron yield over cocktail for a centrality of 10-40%.

Invariant mass spectra from SQRT($s_{NN}$) = 200 GeV Au+Au collisions in different centralities. The ratio of dielectron yield over cocktail for a centrality of 40-80%.

Invariant mass spectra from SQRT($s_{NN}$) = 200 GeV Au+Au collisions in different centralities. The ratio of dielectron yield over cocktail for minimum bias.

The integrated dielectron yield and ratio of dielectron yield over cocktail in different centralities for mass window 0 - 0.15 GeV/c^2, and the comparison with model calculations.

The integrated dielectron yield and ratio of dielectron yield over cocktail in different centralities for mass window 0.15 - 0.3 GeV/c^2, and the comparison with model calculations.

The integrated dielectron yield and ratio of dielectron yield over cocktail in different centralities for mass window 0.3 - 0.76 GeV/c^2, and the comparison with model calculations.

The integrated dielectron yield and ratio of dielectron yield over cocktail in different centralities for mass window 0.76 - 1.05 GeV/c^2, and the comparison with model calculations.

The integrated dielectron yield and ratio of dielectron yield over cocktail in different centralities for mass window 1.05 - 1.8 GeV/c^2, and the comparison with model calculations.

The integrated dielectron yield and ratio of dielectron yield over cocktail in different centralities for mass window 1.8 - 2.8 GeV/c^2, and the comparison with model calculations.

The integrated dielectron yield and ratio of dielectron yield over cocktail in different centralities for mass window 2.8 - 3.5 GeV/c^2, and the comparison with model calculations.

Dielectron invariant mass spectra from minimum-bias.

Dielectron invariant mass spectra from the most central (0-10%) collisions that we are able to achieve most statistics at present.

The ratio of the Npart-scaled dielectron yield between minimum-bias and the most central collisions.

Mass spectrum of the excess (data - cocktail) in the low-mass region in Au+Au minimum-bias collisions compared to model calculations.

The yields scaled by Npart for the $\rho$-like region with the cocktail subtracted.

The Omega-like region without cocktail subtraction as a function of Npart.

The Phi-like region without cocktail subtraction as a function of Npart.

Correlated charm contributions to the dielectron mass spectra for different assumptions of the correlation strength.

Slope parameters Teff versus invariant mass for dielectrons from charm hadron decays.

Omega meson invariant mass distribution from SQRT($s_{NN}$) = 200 GeV Au+Au minimum-bias collisions after subtraction of the combinatorial background using the mixed-event method.

Phi meson invariant mass distribution from SQRT($s_{NN}$) = 200 GeV Au+Au minimum-bias collisions after subtraction of the combinatorial background using the mixed-event method.

The widths and mass positions of the omega signal from data compared to the values from the PDG and the full Geant simulation.

The widths and mass positions of the phi signal from data compared to the values from the PDG and the full Geant simulation.

The efficiency and acceptance correction factor as function of pT for mid-rapidity omega and phi mesons.

The pT distributions of the omega meson invariant yields from SQRT($s_{NN}$) = 200 GeV Au+Au minimum-bias collisions.

The pT distributions of the phi meson invariant yields from SQRT($s_{NN}$) = 200 GeV Au+Au minimum-bias collisions.

Unlike-sign/like-sign pair acceptance difference correction factor with the PHENIX phi acceptance compared with the full acceptance.

Efficiency corrected invariant mass spectra calculated using the STAR data filtered with the PHENIX azimuthal angle acceptance. The data points are compared to cocktail simulations.

The same data points compared to theoretical model calculations of medium vector meson and QGP contributions.