$v_3\{\Psi_1\}$ vs. centrality for $\pi^+$, $\pi^-$, and protons using the event plane method in $\sqrt{s_{\textrm{NN}}}=3$ GeV Au+Au collisions at STAR.

$v_3\{\Psi_1\}$ vs. centrality for $K^+$, and $K^-$ using the event plane method in $\sqrt{s_{\textrm{NN}}}=3$ GeV Au+Au collisions at STAR.

$v_3\{\Psi_1\}$ vs. centrality for $K^+$, and $K^-$ using the event plane method in $\sqrt{s_{\textrm{NN}}}=3$ GeV Au+Au collisions at STAR.

$v_3\{\Psi_1\}$ vs. rapidity for protons in three large centrality bins from a symmetric acceptance across midrapidity in $\sqrt{s_{\textrm{NN}}}=3$ GeV Au+Au collisions at STAR. Protons exhibit an increasingly negative slope as the centrality increases.

$v_3\{\Psi_1\}$ vs. rapidity for protons in three large centrality bins from a symmetric acceptance across midrapidity in $\sqrt{s_{\textrm{NN}}}=3$ GeV Au+Au collisions at STAR. Protons exhibit an increasingly negative slope as the centrality increases.

$v_3\{\Psi_1\}$ vs. rapidity for protons in three large centrality bins from only the backward region in $\sqrt{s_{\textrm{NN}}}=3$ GeV Au+Au collisions at STAR. Note that the $p_{\mathrm{T}}$ acceptance extends to a lower limit than in Fig. 6.

$v_3\{\Psi_1\}$ vs. rapidity for protons in three large centrality bins from only the backward region in $\sqrt{s_{\textrm{NN}}}=3$ GeV Au+Au collisions at STAR. Note that the $p_{\mathrm{T}}$ acceptance extends to a lower limit than in Fig. 6.

$v_3\{\Psi_1\}$ vs. $p_{\mathrm{T}}$ for protons in three large centrality bins in $\sqrt{s_{\textrm{NN}}}=3$ GeV Au+Au collisions at STAR. $v_3\{\Psi_1\}$ becomes increasingly negative as $p_{\mathrm{T}}$ and centrality increase.

$v_3\{\Psi_1\}$ vs. $p_{\mathrm{T}}$ for protons in three large centrality bins in $\sqrt{s_{\textrm{NN}}}=3$ GeV Au+Au collisions at STAR. $v_3\{\Psi_1\}$ becomes increasingly negative as $p_{\mathrm{T}}$ and centrality increase.

$\Delta\phi$ projection of balance function in low $p_{T}$ in 0-10% centrality

$\Delta\phi$ projection of balance function in low $p_{T}$ in 30-40% centrality

1D $\Delta\eta$ projection of balance function in low $p_{T}$ in 70-80% centrality

$\Delta\eta$ projection of balance function in low $p_{T}$ in 0-40 $N_{trk}$

$\Delta\eta$ projection of balance function in low $p_{T}$ in 120-150 $N_{trk}$

$\Delta\eta$ projection of balance function in low $p_{T}$ in 270-300 $N_{trk}$

$\Delta\phi$ projection of balance function in low $p_{T}$ in 0-40 $N_{trk}$

$\Delta\phi$ projection of balance function in low $p_{T}$ in 120-150 $N_{trk}$

$\Delta\phi$ projection of balance function in low $p_{T}$ in 270-300 $N_{trk}$

The width of balance function in $<|\Delta\eta|>$ in low $p_{T}$ for PbPb collisions

Balance function width in low $p_{T}$ PbPb relative $<|\Delta\eta|>/<|\Delta\eta|>_{N_{ch}<65}$

The width of balance function in $<|\Delta\eta|>$ in low $p_{T}$ for pPb collisions

Balance function width in low $p_{T}$ in pPb relative $<|\Delta\eta|>/<|\Delta\eta|>_{N_{ch}<24}$

The width of the balance function in $<|\Delta\phi|>$ in low $p_{T}$ for PbPb collisions

Balance function width in low $p_{T}$ in PbPb for relative $<|\Delta\phi|>/<|\Delta\phi|>_{N_{ch}<65}$

The width of the balance function in $<|\Delta\phi|>$ in low $p_{T}$ for pPb collisions

Balance function width in low $p_{T}$ in pPb for relative $|\Delta\phi|/<|\Delta\phi|>_{N_{ch}<24}$

Balance function projection as a function of $\Delta\eta$ in intermediate $p_{T}$ in PbPb for 0-10% centrality

Balance function projection as a function of $\Delta\eta$ in intermediate $p_{T}$ in PbPb for 30-40% centrality

Balance function projection as a function of $\Delta\eta$ in intermediate $p_{T}$ in PbPb for 70-80% centrality

Balance function projection as a function of $\Delta\phi$ in intermediate $p_{T}$ in PbPb for 0-10% centrality

Balance function projection as a function of $\Delta\phi$ in intermediate $p_{T}$ in PbPb for 30-40% centrality

Balance function projection as a function of $\Delta\phi$ in intermediate $p_{T}$ in PbPb for 70-80% centrality

Balance function projection as a function of $\Delta\eta$ in intermediate $p_{T}$ in pPb for 0-40 $N_{trk}$

Balance function projection as a function of $\Delta\eta$ in intermediate $p_{T}$ in pPb for 120-150 $N_{trk}$

Balance function projection as a function of $\Delta\eta$ in intermediate $p_{T}$ in pPb for 270-300 $N_{trk}$

Balance function projection as a function of $\Delta\phi$ in intermediate $p_{T}$ in pPb for $N_{trk}$ 0-40

Balance function projection as a function of $\Delta\phi$ in intermediate $p_{T}$ in pPb for $N_{trk}$ 120-150

Balance function projection as a function of $\Delta\phi$ in intermediate $p_{T}$ in pPb for $N_{trk}$ 270-300

Balance function projection as a function of $\Delta\eta$ in high $p_{T}$ in PbPb for 0-10% centrality

Balance function projection as a function of $\Delta\eta$ in high $p_{T}$ in PbPb for 30-40% centrality

Balance function projection as a function of $\Delta\eta$ in high $p_{T}$ in PbPb for 70-80% centrality

Balance function projection as a function of $\Delta\phi$ in high $p_{T}$ in PbPb for 0-10% centrality

Balance function projection as a function of $\Delta\phi$ in high $p_{T}$ in PbPb for 30-40% centrality

Balance function projection as a function of $\Delta\phi$ in high $p_{T}$ in PbPb for 70-80% centrality

Balance function projection as a function of $\Delta\eta$ in high $p_{T}$ in pPb for $N_{trk}$ 0-40

Balance function projection as a function of $\Delta\eta$ in high $p_{T}$ in pPb for $N_{trk}$ 120-150

Balance function projection as a function of $\Delta\eta$ in high $p_{T}$ in pPb for $N_{trk}$ 270-300

Balance function projection as a function of $\Delta\phi$ in high $p_{T}$ in PbPb for $N_{trk}$ 0-40

Balance function projection as a function of $\Delta\phi$ in high $p_{T}$ in PbPb for $N_{trk}$ 120-150

Balance function projection as a function of $\Delta\phi$ in high $p_{T}$ in PbPb for $N_{trk}$ 270-300

Average width of the balance function in $\Delta \eta$ with $N_{ch}$ for low $p_{T}$

Average width of the balance function in $\Delta \eta$ with $N_{ch}$ for intermediate $p_{T}$

Average width of the balance function in $\Delta \eta$ with $N_{ch}$ for high $p_{T}$

Average width of the balance function in $\Delta \phi$ with $p_{T}$ for low $p_{T}$

Average width of the balance function in $\Delta \phi$ with $N_{ch}$ for intermediate $p_{T}$

Average width of the balance function in $\Delta \phi$ with $N_{ch}$ for high $p_{T}$

Average width of the balance function in $\Delta \eta$ with $N_{ch}$ for low $p_{T}$

Average width of the balance function in $\Delta \eta$ with $N_{ch}$ for intermediate $p_{T}$

Average width of the balance function in $\Delta \phi$ with $N_{ch}$ for high $p_{T}$

Average width of the balance function in $\Delta \phi$ with $N_{ch}$ for low $p_{T}$

Average width of the balance function in $\Delta \phi$ with $N_{ch}$ for intermediate $p_{T}$

Average width of the balance function in $\Delta \phi$ with $N_{ch}$ for high $p_{T}$

The total coherent $\mathrm{J}/\psi$ photoproduction cross section as a function of photon-nuclear center-of-mass energy per nucleon $W_{\gamma \mathrm{N}}^{\mathrm{Pb}}$, measured in PbPb ultra-peripheral collisions at $\sqrt{s_{\mathrm{NN}}}$ = 5.02 TeV. The $W_{\gamma \mathrm{N}}^{\mathrm{Pb}}$ values used correspond to the center of each rapidity range. The theoretical uncertainties is due to the uncertainties in the photon flux.

The nuclear gluon suppression factor $R_{\mathrm{g}}^{\mathrm{Pb}}$ as a function of Bjorken $x$ extracted from the CMS measurement of the coherent $\mathrm{J}/\psi$ photoproduction in PbPb ultra-peripheral collisions at $\sqrt{s_{\mathrm{NN}}}$ = 5.02 TeV. The $x$ values are evaluated at the centers of their corresponding rapidity ranges. The theoretical uncertainties are due to the uncertainties in the photon flux and the impulse approximation model.

The nuclear gluon suppression factor $R_{\mathrm{g}}^{\mathrm{Pb}}$ as a function of Bjorken $x$ extracted from the CMS measurement of the coherent $\mathrm{J}/\psi$ photoproduction in PbPb ultra-peripheral collisions at $\sqrt{s_{\mathrm{NN}}}$ = 5.02 TeV. The $x$ values are evaluated at the centers of their corresponding rapidity ranges. The theoretical uncertainties are due to the uncertainties in the photon flux and the impulse approximation model.

The total covariance matrix of the total coherent photoproduction cross section as a function of photon-nuclear center-of-mass energy per nucleon $W_{\gamma \mathrm{N}}^{\mathrm{Pb}}$. The covariance matrix includes both the experimental and theoretical (photon flux) uncertainties. The bins are ordered as increasing in $W_{\gamma \mathrm{N}}^{\mathrm{Pb}}$.

The photon flux values from STARLight as a function of rapidity, in different neutron multiplicity classes: 0n0n, 0nXn, and XnXn.

The experimental covariance matrix of the total coherent photoproduction cross section as a function of photon-nuclear center-of-mass energy per nucleon $W_{\gamma \mathrm{N}}^{\mathrm{Pb}}$. The bins are ordered as increasing in $W_{\gamma \mathrm{N}}^{\mathrm{Pb}}$.

The covariance matrix for the flux in the 0n0n neutron multiplicity class.

The theoretical (photon flux) covariance matrix of the total coherent photoproduction cross section as a function of photon-nuclear center-of-mass energy per nucleon $W_{\gamma \mathrm{N}}^{\mathrm{Pb}}$. The bins are ordered as increasing in $W_{\gamma \mathrm{N}}^{\mathrm{Pb}}$.

The covariance matrix for the flux in the 0nXn neutron multiplicity class.

The covariance matrix for the flux in the XnXn neutron multiplicity class.

The total covariance matrix of the total coherent photoproduction cross section as a function of photon-nuclear center-of-mass energy per nucleon $W_{\gamma \mathrm{N}}^{\mathrm{Pb}}$. The covariance matrix includes both the experimental and theoretical (photon flux) uncertainties. The bins are ordered as increasing in $W_{\gamma \mathrm{N}}^{\mathrm{Pb}}$.

The experimental covariance matrix of the total coherent photoproduction cross section as a function of photon-nuclear center-of-mass energy per nucleon $W_{\gamma \mathrm{N}}^{\mathrm{Pb}}$. The bins are ordered as increasing in $W_{\gamma \mathrm{N}}^{\mathrm{Pb}}$.

The theoretical (photon flux) covariance matrix of the total coherent photoproduction cross section as a function of photon-nuclear center-of-mass energy per nucleon $W_{\gamma \mathrm{N}}^{\mathrm{Pb}}$. The bins are ordered as increasing in $W_{\gamma \mathrm{N}}^{\mathrm{Pb}}$.

The elliptic flow after nonflow subtraction, $v_{2}^{sub}/4$, for $f_0(980)$ as a function of $<KE_{T}>/4$ in pPb collision at 8.16 TeV.

The elliptic flow after nonflow subtraction, $v_{2}^{sub}/2$, for $f_0(980)$ as a function of $<p_{T}>/2$ in pPb collision at 8.16 TeV.

The elliptic flow after nonflow subtraction, $v_{2}^{sub}/4$, for $f_0(980)$ as a function of $<p_{T}>/4$ in pPb collision at 8.16 TeV.

Centrality dependence of (a) $v_2${2} and (b) $v_2${4}. (a) The red points indicate no pseudorapidity gap whereas the magenta points indicate a pseudorapidity gap of |$\Delta\eta$| > 2.0. (b) The black points indicate $v_2${4} with no pseudorapidity gap, the blue points indicate a two-subevent method with |$\Delta\eta$| > 2.0 but where some short-range pairs are allowed, and the red points indicate a two-subevent method with |$\Delta\eta$| > 2.0 where no short-range pairs are allowed.

Centrality dependence of (a) $v_2${2} and (b) $v_2${4}. (a) The red points indicate no pseudorapidity gap whereas the magenta points indicate a pseudorapidity gap of |$\Delta\eta$| > 2.0. (b) The black points indicate $v_2${4} with no pseudorapidity gap, the blue points indicate a two-subevent method with |$\Delta\eta$| > 2.0 but where some short-range pairs are allowed, and the red points indicate a two-subevent method with |$\Delta\eta$| > 2.0 where no short-range pairs are allowed.

Multi-particle $v_2$ as a function of centrality in Au+Au collisions at $\sqrt{s_{NN}}$ = 200 GeV. The magenta open diamonds indicate the $v_2${SP }, the blue open squares indicate $v_2${4}, the black open circles indicate $v_2${6}, and the green filled diamonds indicate $v_2${8}.

Multi-particle $v_2$ as a function of centrality in Au+Au collisions at $\sqrt{s_{NN}}$ = 200 GeV. The magenta open diamonds indicate the $v_2${SP }, the blue open squares indicate $v_2${4}, the black open circles indicate $v_2${6}, and the green filled diamonds indicate $v_2${8}.

Multi-particle $v_2$ as a function of centrality in Au+Au collisions at $\sqrt{s_{NN}}$ = 200 GeV. The magenta open diamonds indicate the $v_2${SP }, the blue open squares indicate $v_2${4}, the black open circles indicate $v_2${6}, and the green filled diamonds indicate $v_2${8}.

Multi-particle $v_2$ as a function of centrality in Au+Au collisions at $\sqrt{s_{NN}}$ = 200 GeV. The magenta open diamonds indicate the $v_2${SP }, the blue open squares indicate $v_2${4}, the black open circles indicate $v_2${6}, and the green filled diamonds indicate $v_2${8}.

Cumulant method estimate of $\sigma_{v_{2}}$ / $\langle v_{2} \rangle$ as a function of centrality in Au+Au collisions at $\sqrt{s_{NN}}$ = 200 GeV. The data are shown as black open squares. The same calculation as done in data is done in ampt, shown as a solid green line. Calculations of $\sigma_{\epsilon_{2}}$ /$ \langle\epsilon_{2}\rangle$ performed in the Monte Carlo Glauber model are shown as blue lines. The solid blue line is the Monte Carlo Glauber calculation done using the same estimate as the data, the dashed blue line is the direct calculation of the moments of the MC Glauber $\epsilon_{2}$ distribution.

Centrality dependence of $v_3${2, |$\Delta\eta$| > 2} in Au+Au collisions at $\sqrt{s_{NN}}$ = 200 GeV, shown as magenta diamonds. The systematic uncertainty is indicated as a shaded magenta band. Also shown as black squares are $\sqrt{\langle v^{2}_{3}\rangle}$ as determined from the folding analysis, which is shown in the next section.

Centrality dependence of $v_3${2, |$\Delta\eta$| > 2} in Au+Au collisions at $\sqrt{s_{NN}}$ = 200 GeV, shown as magenta diamonds. The systematic uncertainty is indicated as a shaded magenta band. Also shown as black squares are $\sqrt{\langle v^{2}_{3}\rangle}$ as determined from the folding analysis, which is shown in the next section.

Centrality dependence of $c_3${4} for Au+Au collisions at $\sqrt{s_{NN}}$ = 200 GeV. (a) Calculations using both arms: $c_3${4} (black circles), $c_3${4}$_{ab|ab}$ (blue diamonds), $c_3${4}$_{aa|bb}$ (red squares), and comparison to STAR [36] (black stars). (b) Comparison of $c_3${4} determined using both arms (open symbols) and a single arm (closed symbols). Note that the open black circles are the same in (a) and (b).

Centrality dependence of $c_3${4} for Au+Au collisions at $\sqrt{s_{NN}}$ = 200 GeV. (a) Calculations using both arms: $c_3${4} (black circles), $c_3${4}$_{ab|ab}$ (blue diamonds), $c_3${4}$_{aa|bb}$ (red squares), and comparison to STAR [36] (black stars). (b) Comparison of $c_3${4} determined using both arms (open symbols) and a single arm (closed symbols). Note that the open black circles are the same in (a) and (b).

Centrality dependence of $c_3${4} for Au+Au collisions at $\sqrt{s_{NN}}$ = 200 GeV. (a) Calculations using both arms: $c_3${4} (black circles), $c_3${4}$_{ab|ab}$ (blue diamonds), $c_3${4}$_{aa|bb}$ (red squares), and comparison to STAR [36] (black stars). (b) Comparison of $c_3${4} determined using both arms (open symbols) and a single arm (closed symbols). Note that the open black circles are the same in (a) and (b).

Folding results for (a) $v_2$, (b) $\sigma_{v_{2}}$ , and (c) $\sigma_{v_{2}}$ / $\langle v_{2}\rangle$. The black lines above and below the points indicate the systematic uncertainties. The red (green) boxes indicate the statistical uncertainties at the 68.27% (95.45%) confidence level. In the case of $\langle v_{2}\rangle$, the statistical uncertainties at the 68.27% confidence level are too small to be seen, and the uncertainties at the 95.45% confidence level are visible but noticeably smaller than the marker size. Shown as blue squares are the same quantities as determined using the cumulant based calculation—these points are slightly offset in the x-coordinate to improve visibility.

Folding results for (a) $v_2$, (b) $\sigma_{v_{2}}$ , and (c) $\sigma_{v_{2}}$ / $\langle v_{2}\rangle$. The black lines above and below the points indicate the systematic uncertainties. The red (green) boxes indicate the statistical uncertainties at the 68.27% (95.45%) confidence level. In the case of $\langle v_{2}\rangle$, the statistical uncertainties at the 68.27% confidence level are too small to be seen, and the uncertainties at the 95.45% confidence level are visible but noticeably smaller than the marker size. Shown as blue squares are the same quantities as determined using the cumulant based calculation—these points are slightly offset in the x-coordinate to improve visibility.

Folding results for (a) $v_2$, (b) $\sigma_{v_{2}}$ , and (c) $\sigma_{v_{2}}$ / $\langle v_{2}\rangle$. The black lines above and below the points indicate the systematic uncertainties. The red (green) boxes indicate the statistical uncertainties at the 68.27% (95.45%) confidence level. In the case of $\langle v_{2}\rangle$, the statistical uncertainties at the 68.27% confidence level are too small to be seen, and the uncertainties at the 95.45% confidence level are visible but noticeably smaller than the marker size. Shown as blue squares are the same quantities as determined using the cumulant based calculation—these points are slightly offset in the x-coordinate to improve visibility.

Folding results for (a) $\langle v_{3}\rangle$, (b) $\sigma_{v_{3}}$ , and (c) $\sigma_{v_{3}}$ /$\langle v_{3}\rangle$. The black lines above and below the points indicate the systematic uncertainties. The red (green) boxes indicate the statistical uncertainties at the 68.27% (95.45%) confidence level. The $\sigma_{v_{3}}$ /$\langle v_{3}\rangle$ values are all ≈ 0.52, the apparent limiting value of this quantity for the Bessel-Gaussian distribution.

Folding results for (a) $\langle v_{3}\rangle$, (b) $\sigma_{v_{3}}$ , and (c) $\sigma_{v_{3}}$ /$\langle v_{3}\rangle$. The black lines above and below the points indicate the systematic uncertainties. The red (green) boxes indicate the statistical uncertainties at the 68.27% (95.45%) confidence level. The $\sigma_{v_{3}}$ /$\langle v_{3}\rangle$ values are all ≈ 0.52, the apparent limiting value of this quantity for the Bessel-Gaussian distribution.

Folding results for (a) $\langle v_{3}\rangle$, (b) $\sigma_{v_{3}}$ , and (c) $\sigma_{v_{3}}$ /$\langle v_{3}\rangle$. The black lines above and below the points indicate the systematic uncertainties. The red (green) boxes indicate the statistical uncertainties at the 68.27% (95.45%) confidence level. The $\sigma_{v_{3}}$ /$\langle v_{3}\rangle$ values are all ≈ 0.52, the apparent limiting value of this quantity for the Bessel-Gaussian distribution.

$v_2$ for inclusive charged hadrons in Au+Au at $\sqrt{s_{NN}}$ = 200 GeV.

$v_2$ for inclusive charged hadrons in Au+Au at $\sqrt{s_{NN}}$ = 200 GeV.

$v_2$ for inclusive charged hadrons in Au+Au at $\sqrt{s_{NN}}$ = 200 GeV.

$v_2$ for inclusive charged hadrons in Au+Au at $\sqrt{s_{NN}}$ = 200 GeV.

$v_2$ for inclusive charged hadrons in Au+Au at $\sqrt{s_{NN}}$ = 200 GeV.

$v_2$ for inclusive charged hadrons in Au+Au at $\sqrt{s_{NN}}$ = 62.4 GeV.

$v_2$ for inclusive charged hadrons in Au+Au at $\sqrt{s_{NN}}$ = 62.4 GeV.

$v_2$ for inclusive charged hadrons in Cu+Cu at $\sqrt{s_{NN}}$ = 62.4 GeV compared with 200 GeV.

$v_2$ for inclusive charged hadrons in Cu+Cu at $\sqrt{s_{NN}}$ = 62.4 GeV compared with 200 GeV.

$v_2$ for inclusive charged hadrons in Cu+Cu at $\sqrt{s_{NN}}$ = 62.4 GeV compared with 200 GeV.

Comparison of integrated $v_2$ at $\sqrt{s_{NN}}$ = 62.4 and 200 GeV. Au+Au reaction.

Comparison of integrated $v_2$ at $\sqrt{s_{NN}}$ = 62.4 and 200 GeV. Au+Au reaction.

Comparison of integrated $v_2$ at $\sqrt{s_{NN}}$ = 62.4 and 200 GeV. Cu+Cu reaction.

Comparison of integrated $v_2$ at $\sqrt{s_{NN}}$ = 62.4 and 200 GeV. Cu+Cu reaction.

The comparison of integrated $v_2$ as a function of centrality.

The comparison of the normalized $v_2$/$\epsilon$ vs. centrality.

The comparison of integrated $v_2$ as a function of centrality.

The comparison of the normalized $v_2$/$\epsilon$ vs. centrality.

The comparison of integrated $v_2$ as a function of centrality.

The comparison of the normalized $v_2$/$\epsilon$ vs. centrality.

The comparison of integrated $v_2$ as a function of centrality.

The comparison of the normalized $v_2$/$\epsilon$ vs. centrality.

Comparison of $v_2$ ($p_T$) at 200 GeV for two example systems with different collision size.

Comparison of $v_2$ ($p_T$) at 200 GeV for two example systems with different collision size.

Comparison of $v_2$ ($p_T$) at 200 GeV for two example systems with different collision size.

$v_2$ vs. $p_T$ for $\pi$/$K$/$p$ emitted from Au+Au at 62.4 and 200 GeV and Cu+Cu at 62.4 and 200 GeV for centralities given.

$v_2$ vs. $p_T$ for $\pi$/$K$/$p$ emitted from Au+Au at 62.4 and 200 GeV and Cu+Cu at 62.4 and 200 GeV for centralities given.

$v_2$ vs. $p_T$ for $\pi$/$K$/$p$ emitted from Au+Au at 62.4 and 200 GeV and Cu+Cu at 62.4 and 200 GeV for centralities given.

$v_2$ vs. $p_T$ for $\pi$/$K$/$p$ emitted from Au+Au at 62.4 and 200 GeV and Cu+Cu at 62.4 and 200 GeV for centralities given.

$v_2$ vs. $p_T$ for $\pi$/$K$/$p$ emitted from Au+Au at 62.4 and 200 GeV and Cu+Cu at 62.4 and 200 GeV for centralities given.

$v_2$ vs. $p_T$ for $\pi$/$K$/$p$ emitted from Au+Au at 62.4 and 200 GeV and Cu+Cu at 62.4 and 200 GeV for centralities given.

$v_2$ vs. $p_T$ for $\pi$/$K$/$p$ emitted from Au+Au at 62.4 and 200 GeV and Cu+Cu at 62.4 and 200 GeV for centralities given.

$v_2$ vs. $p_T$ for $\pi$/$K$/$p$ emitted from Au+Au at 62.4 and 200 GeV and Cu+Cu at 62.4 and 200 GeV for centralities given.

$v_2$ vs. $p_T$ for $\pi$/$K$/$p$ emitted from Au+Au at 62.4 and 200 GeV and Cu+Cu at 62.4 and 200 GeV for centralities given.

$v_2$ vs. $p_T$ for $\pi$/$K$/$p$ emitted from Au+Au at 62.4 and 200 GeV and Cu+Cu at 62.4 and 200 GeV for centralities given.

Comparison of $v_2$ between $\sqrt{s_{NN}}$ = 62.4 GeV and 200 GeV for $\pi$/$K$/$p$ emitted from central Au+Au collisions.

Comparison of $v_2$ between $\sqrt{s_{NN}}$ = 62.4 GeV and 200 GeV for $\pi$/$K$/$p$ emitted from central Au+Au collisions.

Comparison of $v_2$ between $\sqrt{s_{NN}}$ = 62.4 GeV and 200 GeV for $\pi$/$K$/$p$ emitted from central Au+Au collisions.

Comparison of $v_2$ between $\sqrt{s_{NN}}$ = 62.4 GeV and 200 GeV for $\pi$/$K$/$p$ emitted from central Au+Au collisions.

Comparison of $v_2$ between $\sqrt{s_{NN}}$ = 62.4 GeV and 200 GeV for $\pi$/$K$/$p$ emitted from central Au+Au collisions.

Comparison of $v_2$ between $\sqrt{s_{NN}}$ = 62.4 GeV and 200 GeV for $\pi$/$K$/$p$ emitted from central Au+Au collisions.

Comparison of $v_2$ between $\sqrt{s_{NN}}$ = 62.4 GeV and 200 GeV for $\pi$/$K$/$p$ emitted from central Au+Au collisions.

Comparison of the $v_2$ of particles, antiparticles, for a minimum bias sample at 200 GeV and central 62.4 GeV Au+Au collisions.

Comparison of the $v_2$ of particles, antiparticles, for a minimum bias sample at 200 GeV and central 62.4 GeV Au+Au collisions.

Comparison of the $v_2$ of particles, antiparticles, for a minimum bias sample at 200 GeV and central 62.4 GeV Au+Au collisions.

Comparison of the $v_2$ of particles, antiparticles, for a minimum bias sample at 200 GeV and central 62.4 GeV Au+Au collisions.

Comparison of the $v_2$ of particles, antiparticles, for a minimum bias sample at 200 GeV and central 62.4 GeV Au+Au collisions.

Comparison of the $v_2$ of particles, antiparticles, for a minimum bias sample at 200 GeV and central 62.4 GeV Au+Au collisions.

The ratio $v_2$/$n_q$ vs. $p_T$/$n_q$ for $\pi$/$K$/$p$ emitted from Au+Au at 62.4 and 200 GeV and Cu+Cu at 200 GeV collisions.

The ratio $v_2$/$n_q$ vs. $p_T$/$n_q$ for $\pi$/$K$/$p$ emitted from Au+Au at 62.4 and 200 GeV and Cu+Cu at 200 GeV collisions.

The ratio $v_2$/$n_q$ vs. $p_T$/$n_q$ for $\pi$/$K$/$p$ emitted from Au+Au at 62.4 and 200 GeV and Cu+Cu at 200 GeV collisions.

The ratio $v_2$/$n_q$ vs. $p_T$/$n_q$ for $\pi$/$K$/$p$ emitted from Au+Au at 62.4 and 200 GeV and Cu+Cu at 200 GeV collisions.

The ratio $v_2$/$n_q$ vs. $p_T$/$n_q$ for $\pi$/$K$/$p$ emitted from Au+Au at 62.4 and 200 GeV and Cu+Cu at 200 GeV collisions.

The ratio $v_2$/$n_q$ vs. $p_T$/$n_q$ for $\pi$/$K$/$p$ emitted from Au+Au at 62.4 and 200 GeV and Cu+Cu at 200 GeV collisions.

The ratio $v_2$/$n_q$ vs. $p_T$/$n_q$ for $\pi$/$K$/$p$ emitted from Au+Au at 62.4 and 200 GeV and Cu+Cu at 200 GeV collisions.

The ratio $v_2$/$n_q$ vs. $p_T$/$n_q$ for $\pi$/$K$/$p$ emitted from Au+Au at 62.4 and 200 GeV and Cu+Cu at 200 GeV collisions.

The ratio $v_2$/$n_q$ vs. $p_T$/$n_q$ for $\pi$/$K$/$p$ emitted from Au+Au at 62.4 and 200 GeV and Cu+Cu at 200 GeV collisions.

The ratio $v_2$/$n_q$ vs. $p_T$/$n_q$ for $\pi$/$K$/$p$ emitted from Au+Au at 62.4 and 200 GeV and Cu+Cu at 200 GeV collisions.

The ratio $v_2$/$n_q$ vs. ${KE}_T$/$n_q$ for $\pi$/$K$/$p$ emitted from Au+Au at 62.4 and 200 GeV and Cu+Cu at 200 GeV collisions.

The ratio $v_2$/$n_q$ vs. ${KE}_T$/$n_q$ for $\pi$/$K$/$p$ emitted from Au+Au at 62.4 and 200 GeV and Cu+Cu at 200 GeV collisions.

The ratio $v_2$/$n_q$ vs. ${KE}_T$/$n_q$ for $\pi$/$K$/$p$ emitted from Au+Au at 62.4 and 200 GeV and Cu+Cu at 200 GeV collisions.

The ratio $v_2$/$n_q$ vs. ${KE}_T$/$n_q$ for $\pi$/$K$/$p$ emitted from Au+Au at 62.4 and 200 GeV and Cu+Cu at 200 GeV collisions.

The ratio $v_2$/$n_q$ vs. ${KE}_T$/$n_q$ for $\pi$/$K$/$p$ emitted from Au+Au at 62.4 and 200 GeV and Cu+Cu at 200 GeV collisions.

The ratio $v_2$/$n_q$ vs. ${KE}_T$/$n_q$ for $\pi$/$K$/$p$ emitted from Au+Au at 62.4 and 200 GeV and Cu+Cu at 200 GeV collisions.

The ratio $v_2$/$n_q$ vs. ${KE}_T$/$n_q$ for $\pi$/$K$/$p$ emitted from Au+Au at 62.4 and 200 GeV and Cu+Cu at 200 GeV collisions.

The ratio $v_2$/$n_q$ vs. ${KE}_T$/$n_q$ for $\pi$/$K$/$p$ emitted from Au+Au at 62.4 and 200 GeV and Cu+Cu at 200 GeV collisions.

The ratio $v_2$/$n_q$ vs. ${KE}_T$/$n_q$ for $\pi$/$K$/$p$ emitted from Au+Au at 62.4 and 200 GeV and Cu+Cu at 200 GeV collisions.

The ratio $v_2$/$n_q$ vs. ${KE}_T$/$n_q$ for $\pi$/$K$/$p$ emitted from Au+Au at 62.4 and 200 GeV and Cu+Cu at 200 GeV collisions.

The ratio $v_2$/$n_q$ vs. ${KE}_T$/$n_q$ for the indicated hadrons emitted from central Au+Au collisions at 62.4 GeV.

The ratio $v_2$/$n_q$ vs. ${KE}_T$/$n_q$ for the indicated hadrons emitted from central Au+Au collisions at 62.4 GeV.

The ratio $v_2$/$n_q$ vs. ${KE}_T$/$n_q$ for the indicated hadrons emitted from central Au+Au collisions at 62.4 GeV.

$v_2$ vs. $p_T$ and $v_2$/($\epsilon * N^{1/3}_{part} * n_q$) vs. ${KE}_T$/$n_q$ for $\pi$/$K$/$p$ in Au+Au at 200 GeV, in Au+Au at 62.4 GeV, and in Cu+Cu at 200 GeV. The values of $v_2$ and $p_T$ in Au+Au at 200 GeV, in Au+Au at 62.4 GeV, and in Cu+Cu at 200 GeV are the same for as figure 14, and the values of $v_2$, $n_q$, and $KE_T$ in Au+Au at 200 GeV, in Au+Au at 62.4 GeV, and in Cu+Cu at 200 GeV are the same for as figure 18.

Transverse momentum spectra for the jet region, obtained as toward-transverse. Events with a leading track with PT>5 GEV at midrapidity are selected. The spectrum is shown in Figure 1 (right panel).

Coalescence parameter B2 as a function of the transverse momentum per nucleon pT/A, in the underlying event (transverse region). Events with a leading track with PT>5 GEV at midrapidity are selected. The B2 is shown in Figure 2 (both panels).

Coalescence parameter B2 as a function of the transverse momentum per nucleon pT/A, in the jet (transverse-toward region). Events with a leading track with PT>5 GEV at midrapidity are selected. The B2 is shown in Figure 2 (both panels).

Net-proton cumulant ratios C4/C2, C5/C1, C6/C2 (0-40$\%$ and 50-60$\%$ centrality) and C3/C1 (0-40$\%$ centrality) in Au+Au collisions from $\sqrt{s_{NN}}$ = 3 - 200 GeV. At $\sqrt{s_{NN}}$ = 3 GeV, measurement is done for protons instead of net-protons, and with halfrapdity coverage (-0.5 $<$ y $<$ 0) while for rest of energies the rapidity coverage is (-0.5 $<$ y $<$ -0.5).

Net-proton cumulant ratios C4/C2, C5/C1, C6/C2 (0-40$\%$ and 50-60$\%$ centrality) and C3/C1 (0-40$\%$ centrality) from UrQMD model in Au+Au collisions from $\sqrt{s_{NN}}$ = 3 -200 GeV. At $\sqrt{s_{NN}}$ = 3 GeV, UrQMD calculation is done for protons instead of net-protons, and with halfrapdity coverage (-0.5 $<$ y $<$ 0) while for rest of energies the rapidity coverage is (-0.5 $<$ y $<$ -0.5).

Fifth-to-first (KS5/KS1) and sixth-to-second (KS6/KS2) order K-statistic ratios of net-proton distributions in Au+Au collisions at $\sqrt{s_{NN}}$ = 7.7 GeV for 0-40$\%$ and 50-60$\%$ centrality.

Net-proton cumulant ratio C6/C2 and corresponding K-statistic ratio (0-40$\%$) in Au+Au collisions at $\sqrt{s_{NN}}$ = 7.7 GeV as a function of sample size. Uncertainties are only statistical.

Cumulant ratios C4/C2, C5/C1, C6/C2 and C3/C1 of net-proton distributions measured with half rapdity coverage (-0.5 $<$ y $<$ 0) in 0-40$\%$ central Au+Au collisions for collider energies from $\sqrt{s_{NN}}$ = 7.7 - 54.4 GeV.

Cumulant ratios C5/C1 of net-proton distributions measured in 40-50$\%$, 60-70$\%$, and 70-80$\%$ Au+Au collisions at $\sqrt{s_{NN}}$ = 200 GeV. This is supplemental data and was not part of the original paper.

The minimum $\chi^2$ and the parameters $T_{fo}$ and $\beta^{max}_T$ for each of the five centrality selections. The best fit parameters are determined by averaging all parameter pairs within the 1$\sigma$ contour. The errors correspond to the standard deviation of the parameter pairs within the 1$\sigma$ $\chi^2$ contour.

Fit parameters for each particle species using Equations $\frac{dN}{dy}$ = $C \cdot$ ($N_{part}$)$^{\alpha_{part}}$ and $\frac{dN}{dy}$ = $C^{\prime} \cdot$ ($N_{coll}$)$^{\alpha_{coll}}$.

Values of the parameters $n_{pp}$ and $x$ from fitting Equation $R \equiv \frac{dN/dy}{N_{part}}$ = (1-$x$) $\cdot n_{pp} \frac{1}{2}$ + $x \cdot n_{pp} \frac{N_{coll}}{N_{part}}$ = $n_{pp}$ [$\frac{1}{2}$ + $x$($\frac{N_{coll}}{N_{part}}$ - $\frac{1}{2})$]to the observed $dN$/$dy$ per $N_{part}$.

Invariant yields for $\pi^{\pm}$ measured in minimum-bias events at midrapidity and normalized to one rapidity unit.

Invariant yields for $K^{\pm}$ measured in minimum-bias events at midrapidity and normalized to one rapidity unit.

Invariant yields for (anti)$p$ measured in minimum-bias events at midrapidity and normalized to one rapidity unit.

Pion invariant yields in each event centrality and $p_T$ bin measured at midrapidity, normalized to one rapidity unit.

Pion invariant yields in each event centrality and $p_T$ bin measured at midrapidity, normalized to one rapidity unit.

Kaon invariant yields in each event centrality and $p_T$ bin measured at midrapidity, normalized to one rapidity unit.

Kaon invariant yields in each event centrality and $p_T$ bin measured at midrapidity, normalized to one rapidity unit.

(Anti)proton invariant yields in each event centrality and $p_T$ bin measured at midrapidity, normalized to one rapidity unit.

(Anti)proton invariant yields in each event centrality and $p_T$ bin measured at midrapidity, normalized to one rapidity unit.

(Anti)proton invariant yields in each event centrality and $p_T$ bin measured at midrapidity, normalized to one rapidity unit.

(Anti)proton invariant yields in each event centrality and $p_T$ bin measured at midrapidity, normalized to one rapidity unit.

(Anti)proton invariant yields in each event centrality and $p_T$ bin measured at midrapidity, normalized to one rapidity unit.

(Anti)proton invariant yields in each event centrality and $p_T$ bin measured at midrapidity, normalized to one rapidity unit.

(Anti)proton invariant yields in each event centrality and $p_T$ bin measured at midrapidity, normalized to one rapidity unit.

Minimum-bias invariant yields for $\pi^{\pm}$ in equal $p_T$ bins.

Minimum-bias invariant yields for $K^{\pm}$ in equal $p_T$ bins.

Minimum-bias invariant yields for (anti)$p$ in equal $p_T$ bins.

Pion invariant yields in each event centrality normalized to one rapidity unit at midrapidity.

Pion invariant yields in each event centrality normalized to one rapidity unit at midrapidity.

Kaon invariant yields in each event centrality normalized to one rapidity unit at midrapidity.

(Anti)proton invariant yields in each event centrality normalized to one rapidity unit at midrapidity.

(Anti)proton invariant yields in each event centrality normalized to one rapidity unit at midrapidity.

The integrated mean $p_T$ for pions, kaons, and (anti)protons produced in the five different classes of event centrality.

The expansion parameters $T_{fo}$ and $\beta^{max}_T$ as a function of the number of participants.

<p>Charged hadron spectra for centrality selected d+Au collisions.</p>

<p>Charged hadron spectra for centrality selected d+Au collisions.</p>

<p>Charged hadron spectra for centrality selected N+Au collisions. All centralities.</p>

<p>Charged hadron spectra for centrality selected N+Au collisions. Centrality A (Central).</p>

<p>Charged hadron spectra for centrality selected N+Au collisions. Centrality B.</p>

<p>Charged hadron spectra for centrality selected N+Au collisions. Centrality C.</p>

<p>Charged hadron spectra for centrality selected N+Au collisions. Centrality D (Peripheral).</p>

<p>$R_{dAu}$ as a function of $p_T$.</p>

<p>$R_{dAu}$ as a function of $p_T$.</p>

<p>$R_{dAu}$ as a function of $p_T$.</p>

<p>$R_{dAu}$ as a function of $p_T$.</p>

<p>$R_{NAu}$ as a function of $p_T$. Centrality A (Peripheral).</p>

<p>$R_{NAu}$ as a function of $p_T$. Centrality B.</p>

<p>$R_{dAu}$ as a function of $p_T$. Centrality C.</p>

<p>$R_{dAu}$ as a function of $p_T$. Centrality D (Central).</p>

<p>$R_{dAu}$ values averaged in three momentum ranges as functions of $\nu−1$. Horizontal bars show the uncertainty in the value of $\nu$ for each centrality class.</p>

<p>$R_{NAu}$ values averaged in three momentum ranges as functions of $\nu−1$. Horizontal bars show the uncertainty in the value of $\nu$ for each centrality class.</p>

<p>Charged hadron spectra for centrality selected d+Au collisions.</p>

The second-order Fourier coefficients $v_2 (p_T)$ for charge-combined identified hadrons $\pi^{\pm}$, $K^{\pm}$, $ p$, and $\bar{p}$ measured at midrapidity in Cu + Au collisions at $\sqrt{S_{NN}}$ = 200 GeV for the centrality classes marked in each panel.

The second-order Fourier coefficients $v_2 (p_T)$ for charge-combined identified hadrons $\pi^{\pm}$, $K^{\pm}$, $ p$, and $\bar{p}$ measured at midrapidity in Cu + Au collisions at $\sqrt{S_{NN}}$ = 200 GeV for the centrality classes marked in each panel.

The second-order Fourier coefficients $v_2 (p_T)$ for charge-combined identified hadrons $\pi^{\pm}$, $K^{\pm}$, $ p$, and $\bar{p}$ measured at midrapidity in Cu + Au collisions at $\sqrt{S_{NN}}$ = 200 GeV for the centrality classes marked in each panel.

The second-order Fourier coefficients $v_2 (p_T)$ for charge-combined identified hadrons $\pi^{\pm}$, $K^{\pm}$, $ p$, and $\bar{p}$ measured at midrapidity in Cu + Au collisions at $\sqrt{S_{NN}}$ = 200 GeV for the centrality classes marked in each panel.

The second-order Fourier coefficients $v_2 (p_T)$ for charge-combined identified hadrons $\pi^{\pm}$, $K^{\pm}$, $ p$, and $\bar{p}$ measured at midrapidity in Cu + Au collisions at $\sqrt{S_{NN}}$ = 200 GeV for the centrality classes marked in each panel.

The third-order Fourier coefficients $v_3(p_T)$ for charge-combined identified hadrons $\pi^{\pm}$, $K^{\pm}$, $p$, and $\bar{p}$ measured at midrapidity in Cu + Au collisions at $\sqrt{S_{NN}}$ = 200 GeV for 0%-30% centrality.

The third-order Fourier coefficients $v_3(p_T)$ for charge-combined identified hadrons $\pi^{\pm}$, $K^{\pm}$, $p$, and $\bar{p}$ measured at midrapidity in Cu + Au collisions at $\sqrt{S_{NN}}$ = 200 GeV for 0%-30% centrality.

The third-order Fourier coefficients $v_3(p_T)$ for charge-combined identified hadrons $\pi^{\pm}$, $K^{\pm}$, $p$, and $\bar{p}$ measured at midrapidity in Cu + Au collisions at $\sqrt{S_{NN}}$ = 200 GeV for 0%-30% centrality.

The first-order Fourier coefficients $v_1(p_T)$ for charge-combined $\pi^{\pm}$, $K^{\pm}$, $ p$, and $\bar{p}$ identified hadrons measured at midrapidity in Cu + Au collisions at $\sqrt{S_{NN}}$ = 200 GeV for 10%-50% centrality.

The first-order Fourier coefficients $v_1(p_T)$ for charge-combined $\pi^{\pm}$, $K^{\pm}$, $ p$, and $\bar{p}$ identified hadrons measured at midrapidity in Cu + Au collisions at $\sqrt{S_{NN}}$ = 200 GeV for 10%-50% centrality.

Scaled second-order Fourier coefficients $v_2(p_T)/\epsilon_2$ for charged hadrons measured at midrapidity in Cu + Cu collisions at $\sqrt{S_{NN}}$ = 200 GeV.

Scaled second-order Fourier coefficients $v_2(p_T)/\epsilon_2$ for charged hadrons measured at midrapidity in Au + Au collisions at $\sqrt{S_{NN}}$ = 200 GeV.

Scaled second-order Fourier coefficients $v_2(p_T)/\epsilon_2$ for charged hadrons measured at midrapidity in Au + Au collisions at $\sqrt{S_{NN}}$ = 200 GeV.

Scaled second-order Fourier coefficients $v_2(p_T)/\epsilon_2$ for charged hadrons measured at midrapidity in Au + Au collisions at $\sqrt{S_{NN}}$ = 200 GeV.

Scaled second-order Fourier coefficients $v_2(p_T)/\epsilon_2$ for charged hadrons measured at midrapidity in Au + Au collisions at $\sqrt{S_{NN}}$ = 200 GeV.

Scaled second-order Fourier coefficients $v_2(p_T)/\epsilon_2$ for charged hadrons measured at midrapidity in Au + Au collisions at $\sqrt{S_{NN}}$ = 200 GeV.

Scaled second-order Fourier coefficients $v_2(p_T)/\epsilon_2$ for charged hadrons measured at midrapidity in Au + Au collisions at $\sqrt{S_{NN}}$ = 200 GeV.

Scaled second-order Fourier coefficients $v_2(p_T)/\epsilon_2$ for charged hadrons measured at midrapidity in Cu + Au collisions at $\sqrt{S_{NN}}$ = 200 GeV.

Scaled second-order Fourier coefficients $v_2(p_T)/(\epsilon_2 N^{1/3}_{part})$ for charged hadrons measured at midrapidity in Cu + Cu collisions at $\sqrt{S_{NN}}$ = 200 GeV.

Scaled second-order Fourier coefficients $v_2(p_T)/(\epsilon_2 N^{1/3}_{part})$ for charged hadrons measured at midrapidity in Au + Au collisions at $\sqrt{S_{NN}}$ = 200 GeV.

Scaled second-order Fourier coefficients $v_2(p_T)/(\epsilon_2 N^{1/3}_{part})$ for charged hadrons measured at midrapidity in Au + Au collisions at $\sqrt{S_{NN}}$ = 200 GeV.

Scaled second-order Fourier coefficients $v_2(p_T)/(\epsilon_2 N^{1/3}_{part})$ for charged hadrons measured at midrapidity in Au + Au collisions at $\sqrt{S_{NN}}$ = 200 GeV.

Scaled second-order Fourier coefficients $v_2(p_T)/(\epsilon_2 N^{1/3}_{part})$ for charged hadrons measured at midrapidity in Au + Au collisions at $\sqrt{S_{NN}}$ = 200 GeV.

Scaled second-order Fourier coefficients $v_2(p_T)/(\epsilon_2 N^{1/3}_{part})$ for charged hadrons measured at midrapidity in Au + Au collisions at $\sqrt{S_{NN}}$ = 200 GeV.

Scaled second-order Fourier coefficients $v_2(p_T)/(\epsilon_2 N^{1/3}_{part})$ for charged hadrons measured at midrapidity in Au + Au collisions at $\sqrt{S_{NN}}$ = 200 GeV.

Scaled second-order Fourier coefficients $v_2(p_T)/(\epsilon_2 N^{1/3}_{part})$ for charged hadrons measured at midrapidity in Cu + Au collisions at $\sqrt{S_{NN}}$ = 200 GeV.

Scaled third-order Fourier coefficients $v_3(p_T)/\epsilon_3$ for charged hadrons measured at midrapidity in Au + Au collisions at $\sqrt{S_{NN}}$ = 200 GeV.

Scaled third-order Fourier coefficients $v_3(p_T)/\epsilon_3$ for charged hadrons measured at midrapidity in Au + Au collisions at $\sqrt{S_{NN}}$ = 200 GeV.

Scaled third-order Fourier coefficients $v_3(p_T)/\epsilon_3$ for charged hadrons measured at midrapidity in Cu + Au collisions at $\sqrt{S_{NN}}$ = 200 GeV.

Scaled third-order Fourier coefficients $\frac{v_3}{($_3 N_{part}^{1/3})}$ for charged hadrons measured at midrapidity in Au + Au collisions at $\sqrt{S_{NN}}$ = 200 GeV.

Scaled third-order Fourier coefficients $\frac{v_3}{(e_3 N_{part}^{1/3})}$ for charged hadrons measured at midrapidity in Au + Au collisions at $\sqrt{S_{NN}}$ = 200 GeV.

Scaled third-order Fourier coefficients $\frac{v_3}{(e_3 N_{part}^{1/3})}$ for charged hadrons measured at midrapidity in Cu + Au collisions at $\sqrt{S_{NN}}$ = 200 GeV.

The second-order Fourier coefficients $v_2$($p_T$ ) for charged hadrons measured at midrapidity in Cu + Au collisions at $\sqrt{S_{NN}}$ = 200 GeV in comparison with hydrodynamics calculations for the centrality classes marked in each panel.

The third-order Fourier coefficients $v_3(p_T)$ for charged hadrons measured at midrapidity in Cu + Au collisions at $\sqrt{S_{NN}}$ = 200 GeV in comparison with hydrodynamics calculations for the centrality classes marked in each panel.

The ratio of $\psi(2S)$ to $J/\psi(1S)$ cross sections, as measured in the dimuon decay channel in pp collisions at 200 GeV over the rapidity interval $1.2<|y|<2.2$, as a function of transverse momentum $p_{T}$.

The double ratios $[\sigma_{\psi(2S)}/\sigma_{\psi(1S)}]_{pA}/[\sigma_{\psi(2S)}/\sigma_{\psi(1S)}]_{pp}$. The global systematic uncertainty is $\pm 14\%$.

The charged-particle multiplicity $dN_{ch}/d\eta$ determined from AMPT simulations and the nuclear overlap $S_T$ from Glauber simulations for the data shown in Fig. 5. Here we assume $dN/d\eta$ = $\frac{3}{2}dN_{ch}/d\eta$.

The charged-particle multiplicity $dN_{ch}/d\eta$ determined from AMPT simulations and the nuclear overlap $S_T$ from Glauber simulations for the data shown in Fig. 5. Here we assume $dN/d\eta$ = $\frac{3}{2}dN_{ch}/d\eta$.

The charged-particle multiplicity $dN_{ch}/d\eta$ determined from AMPT simulations and the nuclear overlap $S_T$ from Glauber simulations for the data shown in Fig. 5. Here we assume $dN/d\eta$ = $\frac{3}{2}dN_{ch}/d\eta$.

Data table for $\sigma_{\psi(2S)}/\sigma_{\psi(1S)}$ ratios, measured in the dimuon decay channel. Here "fwd" refers to the forward rapidity interval $1.2<y<2.2$, "bwd" refers to the backward rapidity interval $-2.2<y<-1.2$, and "comb" refers to the rapidity interval $1.2<|y|<2.2$.

The ratios between the hydrodynamic probes and their Taylor expansions. The ratios are plotted as functions of centrality in PbPb collisions at $\sqrt{s_{NN}}=5.02~\textrm{Te}\textrm{V}$. The measurement is performed with charged particles within the acceptance region.

The ratios between the hydrodynamic probes and their Taylor expansions. The ratios are plotted as functions of centrality in PbPb collisions at $\sqrt{s_{NN}}=5.02~\textrm{Te}\textrm{V}$. The measurement is performed with charged particles within the acceptance region.

The magnitudes of the measured skewness $\gamma_1^\textrm{exp}$ and "corrected" skewness $\gamma_{1,\textrm{corr}}^\textrm{exp}$ as functions of centrality in PbPb collisions at $\sqrt{s_{NN}}=5.02~\textrm{Te}\textrm{V}$. The measurement is performed with charged particles within the acceptance region.

The magnitudes of the measured kurtosis $\gamma_2^\textrm{exp}$ and "corrected" kurtosis $\gamma_{2,\textrm{corr}}^\textrm{exp}$ as functions of centrality in PbPb collisions at $\sqrt{s_{NN}}=5.02~\textrm{Te}\textrm{V}$. The measurement is performed with charged particles within the acceptance region.

The magnitudes of the measured superskewness $\gamma_3^\textrm{exp}$ and "corrected" superskewness $\gamma_{3,\textrm{corr}}^\textrm{exp}$ as functions of centrality in PbPb collisions at $\sqrt{s_{NN}}=5.02~\textrm{Te}\textrm{V}$. The measurement is performed with charged particles within the acceptance region.

Efficiency corrected cumulants of net-charge distributions as a function of $\langle N_{part} \rangle$ from Au+Au collisions at different collision energies.

$\langle N_{part} \rangle$ dependence of efficiency corrected $\mu / \sigma^2$ of net-charge distributions for Au+Au collisions at different collision energies.

$\langle N_{part} \rangle$ dependence of efficiency corrected $S \sigma$ of net-charge distributions for Au+Au collisions at different collision energies.

$\langle N_{part} \rangle$ dependence of efficiency corrected $\kappa \sigma^2$ of net-charge distributions for Au+Au collisions at different collision energies.

$\langle N_{part} \rangle$ dependence of efficiency corrected $S \sigma^3 / \mu$ of net-charge distributions for Au+Au collisions at different collision energies.

The energy dependence of efficiency corrected $\mu / \sigma^2$, $S \sigma$, $\kappa \sigma^2$, and $S \sigma^3 / \mu$ of netcharge distributions for central (0%–5%) Au+Au collisions.

The energy dependence of the chemical freeze-out parameter $\mu_B$.

Away-side jet widths from a Gaussian fit by $h^{\pm}$ partner momentum for various $\pi^0$ trigger momenta in Au+Au collisions.

Away-side $I_{AA}$ for the entire away-side $|\Delta \phi - \pi| < \pi /2$ selection vs. $h^{\pm}$ partner momentum for various $\pi^0$ trigger momenta. A point-to-point uncorrelated 6% normalization uncertainty (mainly due to efficiency corrections) applies to all measurements.

Away-side $I_{AA}$ for the entire away-side $|\Delta \phi - \pi| < \pi /2$ selection vs. $h^{\pm}$ partner momentum for various $\pi^0$ trigger momenta. A point-to-point uncorrelated 6% normalization uncertainty (mainly due to efficiency corrections) applies to all measurements.

Away-side $I_{AA}$ for a narrow “head” $|\Delta \phi - \pi| < \pi /6$ selection vs. $h^{\pm}$ partner momentum for various $\pi^0$ trigger momenta. A point-to-point uncorrelated 6% normalization uncertainty (mainly due to efficiency corrections) applies to all measurements.

Away-side $I_{AA}$ for a narrow “head” $|\Delta \phi - \pi| < \pi /6$ selection vs. $h^{\pm}$ partner momentum for various $\pi^0$ trigger momenta. A point-to-point uncorrelated 6% normalization uncertainty (mainly due to efficiency corrections) applies to all measurements.

Supplemental near-side jet widths from a Gaussian fit by $h^{\pm}$ partner momentum for various $\pi^0$ trigger momenta in $p+p$ collisions.

Supplemental near-side jet widths from a Gaussian fit by $h^{\pm}$ partner momentum for various $\pi^0$ trigger momenta in Au+Au collisions.

Supplemental near-side jet widths from a Gaussian fit by $h^{\pm}$ partner momentum for various $\pi^0$ trigger momenta in Au+Au collisions.

Supplemental away-side $I_{AA}$ for $|\Delta \phi - \pi| < \pi /3$ selection vs. $h^{\pm}$ partner momentum for various $\pi^0$ trigger momenta. A point-to-point uncorrelated 6% normalization uncertainty (mainly due to efficiency corrections) applies to all measurements.

Supplemental away-side $I_{AA}$ for $|\Delta \phi - \pi| < \pi /3$ selection vs. $h^{\pm}$ partner momentum for various $\pi^0$ trigger momenta. A point-to-point uncorrelated 6% normalization uncertainty (mainly due to efficiency corrections) applies to all measurements.

$p+p$ yield vs $p_T^a$ supplemental tables.

Au+Au yield vs $p_T^a$ supplemental tables.

Au+Au yield vs $p_T^a$ supplemental tables.

$z_{g}$ for Matched Trigger jets in pp$+$AuAu Data anti-kT R$=$0.4

$z_{g}$ for HardCore Recoil jets in AuAu Data anti-kT R$=$0.4

$z_{g}$ for HardCore Recoil jets in pp$+$AuAu Data anti-kT R$=$0.4

$z_{g}$ for Matched Recoil jets in AuAu Data anti-kT R$=$0.4

$z_{g}$ for Matched Recoil jets in pp$+$AuAu Data anti-kT R$=$0.4

$R_{g}$ for HardCore Trigger jets in AuAu Data anti-kT R$=$0.4

$R_{g}$ for HardCore Trigger jets in pp$+$AuAu Data anti-kT R$=$0.4

$R_{g}$ for Matched Trigger jets in AuAu Data anti-kT R$=$0.4

$R_{g}$ for Matched Trigger jets in pp$+$AuAu Data anti-kT R$=$0.4

$R_{g}$ for HardCore Recoil jets in AuAu Data anti-kT R$=$0.4

$R_{g}$ for HardCore Recoil jets in pp$+$AuAu Data anti-kT R$=$0.4

$R_{g}$ for Matched Recoil jets in AuAu Data anti-kT R$=$0.4

$R_{g}$ for Matched Recoil jets in pp$+$AuAu Data anti-kT R$=$0.4

Fully corrected two subjet $z_{sj}$ for $15 < p_{\rm{T, jet}} < 20$ GeV/c, R$=$0.4 anti-kT and R$=$0.1 subjets

Fully corrected two subjet $z_{sj}$ for $20 < p_{\rm{T, jet}} < 25$ GeV/c, R$=$0.4 anti-kT and R$=$0.1 subjets

Fully corrected two subjet $z_{sj}$ for $25 < p_{\rm{T, jet}} < 30$ GeV/c, R$=$0.4 anti-kT and R$=$0.1 subjets

Fully corrected two subjet $z_{sj}$ for $30 < p_{\rm{T, jet}} < 40$ GeV/c, R$=$0.4 anti-kT and R$=$0.1 subjets

Fully corrected two subjet $\theta_{sj}$ for $15 < p_{\rm{T, jet}} < 20$ GeV/c, R$=$0.4 anti-kT and R$=$0.1 subjets

Fully corrected two subjet $\theta_{sj}$ for $20 < p_{\rm{T, jet}} < 25$ GeV/c, R$=$0.4 anti-kT and R$=$0.1 subjets

Fully corrected two subjet $\theta_{sj}$ for $25 < p_{\rm{T, jet}} < 30$ GeV/c, R$=$0.4 anti-kT and R$=$0.1 subjets

Fully corrected two subjet $\theta_{sj}$ for $30 < p_{\rm{T, jet}} < 40$ GeV/c, R$=$0.4 anti-kT and R$=$0.1 subjets

$z_{sj}$ for HardCore Trigger jets in AuAu Data anti-kT R$=$0.4

$z_{sj}$ for HardCore Trigger jets in pp$+$AuAu Data anti-kT R$=$0.4

$z_{sj}$ for Matched Trigger jets in AuAu Data anti-kT R$=$0.4

$z_{sj}$ for Matched Trigger jets in pp$+$AuAu Data anti-kT R$=$0.4

$z_{sj}$ for HardCore Recoil jets in AuAu Data anti-kT R$=$0.4

$z_{sj}$ for HardCore Recoil jets in pp$+$AuAu Data anti-kT R$=$0.4

$z_{sj}$ for Matched Recoil jets in AuAu Data anti-kT R$=$0.4

$z_{sj}$ for Matched Recoil jets in pp$+$AuAu Data anti-kT R$=$0.4

$\theta_{sj}$ for HardCore Trigger jets in AuAu Data anti-kT R$=$0.4

$\theta_{sj}$ for HardCore Trigger jets in pp$+$AuAu Data anti-kT R$=$0.4

$\theta_{sj}$ for Matched Trigger jets in AuAu Data anti-kT R$=$0.4

$\theta_{sj}$ for Matched Trigger jets in pp$+$AuAu Data anti-kT R$=$0.4

$\theta_{sj}$ for HardCore Recoil jets in AuAu Data anti-kT R$=$0.4

$\theta_{sj}$ for HardCore Recoil jets in pp$+$AuAu Data anti-kT R$=$0.4

$\theta_{sj}$ for Matched Recoil jets in AuAu Data anti-kT R$=$0.4

$\theta_{sj}$ for Matched Recoil jets in pp$+$AuAu Data anti-kT R$=$0.4

$A_{J}$ for HardCore R=0.4 anti-kT jets with $\theta_{sj} > 0.1$ in Au$+$Au collisions

$A_{J}$ for HardCore R=0.4 anti-kT jets with $\theta_{sj} > 0.1$ in pp$+$AuAu reference

$A_{J}$ for HardCore R=0.4 anti-kT jets with $0.1 < \theta_{sj} < 0.2$ in Au$+$Au collisions

$A_{J}$ for HardCore R=0.4 anti-kT jets with $0.1 < \theta_{sj} < 0.2$ in pp$+$AuAu reference

$A_{J}$ for HardCore R=0.4 anti-kT jets with $0.2 < \theta_{sj} < 0.3$ in Au$+$Au collisions

$A_{J}$ for HardCore R=0.4 anti-kT jets with $0.2 < \theta_{sj} < 0.3$ in pp$+$AuAu reference

$A_{J}$ for Matched R=0.4 anti-kT jets with $\theta_{sj} > 0.1$ in Au$+$Au collisions

$A_{J}$ for Matched R=0.4 anti-kT jets with $\theta_{sj} > 0.1$ in pp$+$AuAu reference

$A_{J}$ for Matched R=0.4 anti-kT jets with $0.1 < \theta_{sj} < 0.2$ in Au$+$Au collisions

$A_{J}$ for Matched R=0.4 anti-kT jets with $0.1 < \theta_{sj} < 0.2$ in pp$+$AuAu reference

$A_{J}$ for Matched R=0.4 anti-kT jets with $0.2 < \theta_{sj} < 0.3$ in Au$+$Au collisions

$A_{J}$ for Matched R=0.4 anti-kT jets with $0.2 < \theta_{sj} < 0.3$ in pp$+$AuAu reference

tracking efficiency in pp collisions as a function of track momentum

tracking efficiency in AuAu collisions as a function of track momentum

absolutely normalized dijet yields scaled by 1/<TAA> in 20-40% central PbPb collisions

absolutely normalized dijet yields scaled by 1/<TAA> in 40-60% central PbPb collisions

absolutely normalized dijet yields scaled by 1/<TAA> in 60-80% central PbPb collisions

self normalized dijets from pp collisions

self normalized dijet distributions in 0-10% central PbPb collisions

self normalized dijet distributions in 10-20% central PbPb collisions

self normalized dijet distributions in 20-40% central PbPb collisions

self normalized dijet distributions in 40-60% central PbPb collisions

self normalized dijet distributions in 60-80% central PbPb collisions

leading jet RAA^pair in 0-10% central PbPb collisions

subleading jet RAA^pair in 0-10% central PbPb collisions

leading jet RAA^pair in 10-20% central PbPb collisions

subleading jet RAA^pair in 10-20% central PbPb collisions

leading jet RAA^pair in 20-40% central PbPb collisions

subleading jet RAA^pair in 20-40% central PbPb collisions

leading jet RAA^pair in 40-60% central PbPb collisions

subleading jet RAA^pair in 40-60% central PbPb collisions

leading jet RAA^pair in 60-80% central PbPb collisions

subleading jet RAA^pair in 60-80% central PbPb collisions

ratio of subleading jet RAA^pair to leading jet RAA^pair in PbPb collisions

\Delta v_1 between proton and anti-proton vs rapidity in Au+Au 200 GeV 50-80% centrality.

\Delta v_1 between proton and anti-proton vs rapidity in Zr+Zr and Ru+Ru (combined) 200 GeV 50-80% centrality.

\Delta v_1 between proton and anti-proton vs rapidity in Au+Au 27 GeV 50-80% centrality.

\Delta dv1/dy as a function of centrality in Au+Au 200 GeV.

\Delta dv1/dy as a function of centrality in Ru+Ru and Zr+Zr 200 GeV

\Delta dv1/dy as a function of centrality in Au+Au 27 GeV.

Forward-rapidity $J/\psi$ —> $\mu^+\mu^-$ $d$+Au Nuclear Dependence at $\sqrt{s}$ = 200 GeV. (sys. A, B systematics are relative, i.e. they multiply the $R_{dAu}$ value. The sysA uncertainty includes statistical uncertainties as well as point-to-point uncorrelated systematics, sys. B represents uncertainties that are correlated from point to point, and sys. C represents uncertainties in the overall normalization.)

Mid rapidity $d$+Au —> $e^+e^-$ $J/\psi$ $R_{dAu}$ at $\sqrt{s}$=200 GeV. (Sys. A, B systematics are absolute, i.e. they add/subtract directly from $R_{dAu}$. The sys. A uncertainty includes statistical uncertainties as well as point-to-point uncorrelated systematics, sys. B represents uncertainties that are correlated from point to point, and sys. C represents uncertainties in the overall normalization.)

$d$+Au vs centrality vs $y$. (All uncertainties are absolute. The sys. A uncertainty includes both the statistical uncertainty and the point-to-point uncorrelated systematic, sys. B represents uncertainties that are correlated from point to point, and sys. C represents uncertainties in the overall normalization.)

Forward-rapidity $J/\psi$ —> $\mu^+\mu^-$ $d$+Au Nuclear Dependence at $\sqrt{s}$ = 200 GeV. (sys. A, B systematics are relative, i.e. they multiply the $R_{dAu}$ value. The sysA uncertainty includes statistical uncertainties as well as point-to-point uncorrelated systematics, sys. B represents uncertainties that are correlated from point to point, and sys. C represents uncertainties in the overall normalization.)

Mid rapidity $d$+Au —> $e^+e^-$ $J/\psi$ $R_{dAu}$ at $\sqrt{s}$=200 GeV. (Sys. A, B systematics are absolute, i.e. they add/subtract directly from $R_{dAu}$. The sys. A uncertainty includes statistical uncertainties as well as point-to-point uncorrelated systematics, sys. B represents uncertainties that are correlated from point to point, and sys. C represents uncertainties in the overall normalization.)

Forward-rapidity $J/\psi$ —> $\mu^+\mu^-$ $d$+Au Nuclear Dependence at $\sqrt{s}$ = 200 GeV. (sys. A, B systematics are relative, i.e. they multiply the $R_{CP}$ value. The sys. A uncertainty includes statistical uncertainties as well as point-to-point uncorrelated systematics, sys. B represents uncertainties that are correlated from point to point, and sys. C represents uncertainties in the overall normalization.)

Mid rapidity $d$+Au —> $e^+e^-$ $J/\psi$ $R_{CP}$ at $\sqrt{s}$=200 GeV. (Sys. A, B systematics are absolute, i.e. they add/subtract directly from $R_{CP}$. The sys. A uncertainty includes statistical uncertainties as well as point-to-point uncorrelated systematics, sys. B represents uncertainties that are correlated from point to point, and sys. C represents uncertainties in the overall normalization.)

Nuclear modification factor of the $\psi$(2S) shown as a function of transverse momentum

Ratio of the $\psi$(2S) over J/$\psi$ cross sections, not corrected for the branching ratio, shown as a function of transverse momentum

Double ratio of the $\psi$(2S) over J/$\psi$ cross sections in Pb--Pb and pp collisions shown as a function of transverse momentum

$P_z$ sine coefficients of $\Lambda+\bar{\Lambda}$ as a function of transverse momentum in isobar collisions at 200 GeV.

$P_z$ sine coefficients of $\Lambda+\bar{\Lambda}$ as a function of the number of participants in isobar collisions at 200 GeV. The number of participants for isobar collisions were taken from Table III of Phys.Rev.C105, 014901 (2022) and averaged. Same data as the centrality dependence in Fig. 3.

The $\Delta p_{\mathrm{T}}$ of photon-tagged jets as a function of jet $p_{\mathrm{T}}$.

The double $R_{\mathrm{AA}}$ ratio of $R_{\mathrm{AA}}^{\gamma\text{-jet}}/R_{\mathrm{AA}}^{\text{inclusive jet}}$ as a function of jet $p_{\mathrm{T}}$.

The differential cross-section of photon-tagged jets as a function of jet $p_{\mathrm{T}}$ in pp collisions.

The yields of photon-tagged jets as a function of jet $p_{\mathrm{T}}$ in Pb+Pb collisions for different centrality intervals.

The nuclear modification factor of photon-tagged jets as a function of jet $p_{\mathrm{T}}$ for different centrality intervals.

The $\Delta p_{\mathrm{T}}$ of photon-tagged jets as a function of jet $p_{\mathrm{T}}$.

The $R_{\mathrm{CP}}$ of photon-tagged jets as a function of jet $p_{\mathrm{T}}$.

(Color online) Invariant mass spectrum of $e^+e^-$ pairs inclusive in $p_T$ compared to expectations from the model of hadron decays for $p$+$p$ and for different Au+Au centrality classes. The charmed meson decay contribution based on PYTHIA [55] is included in the sum of sources (solid black line). The dotted line shows the contribution from charm calculated assuming an isotropic angular distribution. Statistical (bars) and systematic (boxes) uncertainties are shown separately. The systematic uncertainty on the expected hadronic sources is not shown: it ranges from ~10% in the $\pi^o$ region to ~30% in the region of the vector mesons. The uncertainty on the charm cross section, which dominates the IMR, is ~30% in both $p$+$p$ and in Au+Au collisions.

(Color online) Invariant mass spectrum of $e^+e^-$ pairs inclusive in $p_T$ compared to expectations from the model of hadron decays for $p$+$p$ and for different Au+Au centrality classes. The charmed meson decay contribution based on PYTHIA [55] is included in the sum of sources (solid black line). The dotted line shows the contribution from charm calculated assuming an isotropic angular distribution. Statistical (bars) and systematic (boxes) uncertainties are shown separately. The systematic uncertainty on the expected hadronic sources is not shown: it ranges from ~10% in the $\pi^o$ region to ~30% in the region of the vector mesons. The uncertainty on the charm cross section, which dominates the IMR, is ~30% in both $p$+$p$ and in Au+Au collisions.

(Color online) Invariant mass spectrum of $e^+e^-$ pairs inclusive in $p_T$ compared to expectations from the model of hadron decays for $p$+$p$ and for different Au+Au centrality classes. The charmed meson decay contribution based on PYTHIA [55] is included in the sum of sources (solid black line). The dotted line shows the contribution from charm calculated assuming an isotropic angular distribution. Statistical (bars) and systematic (boxes) uncertainties are shown separately. The systematic uncertainty on the expected hadronic sources is not shown: it ranges from ~10% in the $\pi^o$ region to ~30% in the region of the vector mesons. The uncertainty on the charm cross section, which dominates the IMR, is ~30% in both $p$+$p$ and in Au+Au collisions.

(Color online) Invariant mass spectrum of $e^+e^-$ pairs inclusive in $p_T$ compared to expectations from the model of hadron decays for $p$+$p$ and for different Au+Au centrality classes. The charmed meson decay contribution based on PYTHIA [55] is included in the sum of sources (solid black line). The dotted line shows the contribution from charm calculated assuming an isotropic angular distribution. Statistical (bars) and systematic (boxes) uncertainties are shown separately. The systematic uncertainty on the expected hadronic sources is not shown: it ranges from ~10% in the $\pi^o$ region to ~30% in the region of the vector mesons. The uncertainty on the charm cross section, which dominates the IMR, is ~30% in both $p$+$p$ and in Au+Au collisions.

(Color online) Invariant mass spectrum of $e^+e^-$ pairs inclusive in $p_T$ compared to expectations from the model of hadron decays for $p$+$p$ and for different Au+Au centrality classes. The charmed meson decay contribution based on PYTHIA [55] is included in the sum of sources (solid black line). The dotted line shows the contribution from charm calculated assuming an isotropic angular distribution. Statistical (bars) and systematic (boxes) uncertainties are shown separately. The systematic uncertainty on the expected hadronic sources is not shown: it ranges from ~10% in the $\pi^o$ region to ~30% in the region of the vector mesons. The uncertainty on the charm cross section, which dominates the IMR, is ~30% in both $p$+$p$ and in Au+Au collisions.

(Color online) Invariant mass spectrum of $e^+e^-$ pairs inclusive in $p_T$ compared to expectations from the model of hadron decays for $p$+$p$ and for different Au+Au centrality classes. The charmed meson decay contribution based on PYTHIA [55] is included in the sum of sources (solid black line). The dotted line shows the contribution from charm calculated assuming an isotropic angular distribution. Statistical (bars) and systematic (boxes) uncertainties are shown separately. The systematic uncertainty on the expected hadronic sources is not shown: it ranges from ~10% in the $\pi^o$ region to ~30% in the region of the vector mesons. The uncertainty on the charm cross section, which dominates the IMR, is ~30% in both $p$+$p$ and in Au+Au collisions.

(Color online) Invariant mass spectrum of $e^+e^-$ pairs inclusive in $p_T$ compared to expectations from the model of hadron decays for $p$+$p$ and for different Au+Au centrality classes. The charmed meson decay contribution based on PYTHIA [55] is included in the sum of sources (solid black line). The dotted line shows the contribution from charm calculated assuming an isotropic angular distribution. Statistical (bars) and systematic (boxes) uncertainties are shown separately. The systematic uncertainty on the expected hadronic sources is not shown: it ranges from ~10% in the $\pi^o$ region to ~30% in the region of the vector mesons. The uncertainty on the charm cross section, which dominates the IMR, is ~30% in both $p$+$p$ and in Au+Au collisions.

(Color online) Invariant mass spectrum of $e^+e^-$ pairs inclusive in $p_T$ compared to expectations from the model of hadron decays for $p$+$p$ and for different Au+Au centrality classes. The charmed meson decay contribution based on PYTHIA [55] is included in the sum of sources (solid black line). The dotted line shows the contribution from charm calculated assuming an isotropic angular distribution. Statistical (bars) and systematic (boxes) uncertainties are shown separately. The systematic uncertainty on the expected hadronic sources is not shown: it ranges from ~10% in the $\pi^o$ region to ~30% in the region of the vector mesons. The uncertainty on the charm cross section, which dominates the IMR, is ~30% in both $p$+$p$ and in Au+Au collisions.

(Color online) Invariant mass spectrum of $e^+e^-$ pairs inclusive in $p_T$ compared to expectations from the model of hadron decays for $p$+$p$ and for different Au+Au centrality classes. The charmed meson decay contribution based on PYTHIA [55] is included in the sum of sources (solid black line). The dotted line shows the contribution from charm calculated assuming an isotropic angular distribution. Statistical (bars) and systematic (boxes) uncertainties are shown separately. The systematic uncertainty on the expected hadronic sources is not shown: it ranges from ~10% in the $\pi^o$ region to ~30% in the region of the vector mesons. The uncertainty on the charm cross section, which dominates the IMR, is ~30% in both $p$+$p$ and in Au+Au collisions.

(Color online) Invariant mass spectrum of $e^+e^-$ pairs inclusive in $p_T$ compared to expectations from the model of hadron decays for $p$+$p$ and for different Au+Au centrality classes. The charmed meson decay contribution based on PYTHIA [55] is included in the sum of sources (solid black line). The dotted line shows the contribution from charm calculated assuming an isotropic angular distribution. Statistical (bars) and systematic (boxes) uncertainties are shown separately. The systematic uncertainty on the expected hadronic sources is not shown: it ranges from ~10% in the $\pi^o$ region to ~30% in the region of the vector mesons. The uncertainty on the charm cross section, which dominates the IMR, is ~30% in both $p$+$p$ and in Au+Au collisions.

(Color online) Invariant mass spectrum of $e^+e^-$ pairs inclusive in $p_T$ compared to expectations from the model of hadron decays for $p$+$p$ and for different Au+Au centrality classes. The charmed meson decay contribution based on PYTHIA [55] is included in the sum of sources (solid black line). The dotted line shows the contribution from charm calculated assuming an isotropic angular distribution. Statistical (bars) and systematic (boxes) uncertainties are shown separately. The systematic uncertainty on the expected hadronic sources is not shown: it ranges from ~10% in the $\pi^o$ region to ~30% in the region of the vector mesons. The uncertainty on the charm cross section, which dominates the IMR, is ~30% in both $p$+$p$ and in Au+Au collisions.

(Color online) Dielectron yield per binary collision in the mass range 1.2 to 2.8 GeV/$c^2$ as a function of $N_{part}$. Statistical and systematic uncertainties are shown separately. Also shown are two bands corresponding to different estimates of the contribution from charmed meson decays. The width of the bands reflects the uncertainty of the charm cross section only.

(Color online) Dielectron yield per binary collision in the mass range 1.2 to 2.8 GeV/$c^2$ as a function of $N_{part}$. Statistical and systematic uncertainties are shown separately. Also shown are two bands corresponding to different estimates of the contribution from charmed meson decays. The width of the bands reflects the uncertainty of the charm cross section only.

(Color online) Dielectron yield per participating nucleon pair ($N_{part}/2$) as function of $N_{part}$ for two different mass ranges (a: $0.15<m_{ee}<0.75$ GeV/$c^2$, b: $0<m_{ee}<0.1$ GeV/$c^2$) compared to the expected yield from the hadron decay model. The two lines give the systematic uncertainty of the yield from cocktail and charmed hadron decays. For the data statistical and systematic uncertainties are shown separately.

(Color online) Dielectron yield per participating nucleon pair ($N_{part}/2$) as function of $N_{part}$ for two different mass ranges (a: $0.15<m_{ee}<0.75$ GeV/$c^2$, b: $0<m_{ee}<0.1$ GeV/$c^2$) compared to the expected yield from the hadron decay model. The two lines give the systematic uncertainty of the yield from cocktail and charmed hadron decays. For the data statistical and systematic uncertainties are shown separately.

(Color online) $e^+e^-$ pair invariant mass distributions in $p$+$p$ (left) and minimum bias Au+Au collisions (right). The $p_T$ ranges are shown in the legend. The solid curves represent the cocktail of hadronic sources (see Sec. IV) and include contribution from charm calculated by PYTHIA using the cross section from Ref. [48] scaled by $N_{coll}$.

(Color online) $e^+e^-$ pair invariant mass distributions in $p$+$p$ (left) and minimum bias Au+Au collisions (right). The $p_T$ ranges are shown in the legend. The solid curves represent the cocktail of hadronic sources (see Sec. IV) and include contribution from charm calculated by PYTHIA using the cross section from Ref. [48] scaled by $N_{coll}$.

(Color online) $e^+e^-$ pair invariant mass distributions in $p$+$p$ (left) and minimum bias Au+Au collisions (right). The $p_T$ ranges are shown in the legend. The solid curves represent the cocktail of hadronic sources (see Sec. IV) and include contribution from charm calculated by PYTHIA using the cross section from Ref. [48] scaled by $N_{coll}$.

(Color online) $e^+e^-$ pair invariant mass distributions in $p$+$p$ (left) and minimum bias Au+Au collisions (right). The $p_T$ ranges are shown in the legend. The solid curves represent the cocktail of hadronic sources (see Sec. IV) and include contribution from charm calculated by PYTHIA using the cross section from Ref. [48] scaled by $N_{coll}$.

(Color online) $e^+e^-$ pair invariant mass distributions in $p$+$p$ (left) and minimum bias Au+Au collisions (right). The $p_T$ ranges are shown in the legend. The solid curves represent the cocktail of hadronic sources (see Sec. IV) and include contribution from charm calculated by PYTHIA using the cross section from Ref. [48] scaled by $N_{coll}$.

(Color online) $e^+e^-$ pair invariant mass distributions in $p$+$p$ (left) and minimum bias Au+Au collisions (right). The $p_T$ ranges are shown in the legend. The solid curves represent the cocktail of hadronic sources (see Sec. IV) and include contribution from charm calculated by PYTHIA using the cross section from Ref. [48] scaled by $N_{coll}$.

(Color online) $e^+e^-$ pair invariant mass distributions in $p$+$p$ (left) and minimum bias Au+Au collisions (right). The $p_T$ ranges are shown in the legend. The solid curves represent the cocktail of hadronic sources (see Sec. IV) and include contribution from charm calculated by PYTHIA using the cross section from Ref. [48] scaled by $N_{coll}$.

(Color online) $e^+e^-$ pair invariant mass distributions in $p$+$p$ (left) and minimum bias Au+Au collisions (right). The $p_T$ ranges are shown in the legend. The solid curves represent the cocktail of hadronic sources (see Sec. IV) and include contribution from charm calculated by PYTHIA using the cross section from Ref. [48] scaled by $N_{coll}$.

(Color online) $e^+e^-$ pair invariant mass distributions in $p$+$p$ (left) and minimum bias Au+Au collisions (right). The $p_T$ ranges are shown in the legend. The solid curves represent the cocktail of hadronic sources (see Sec. IV) and include contribution from charm calculated by PYTHIA using the cross section from Ref. [48] scaled by $N_{coll}$.

(Color online) $e^+e^-$ pair invariant mass distributions in $p$+$p$ (left) and minimum bias Au+Au collisions (right). The $p_T$ ranges are shown in the legend. The solid curves represent the cocktail of hadronic sources (see Sec. IV) and include contribution from charm calculated by PYTHIA using the cross section from Ref. [48] scaled by $N_{coll}$.

(Color online) $e^+e^-$ pair invariant mass distributions in $p$+$p$ (left) and minimum bias Au+Au collisions (right). The $p_T$ ranges are shown in the legend. The solid curves represent the cocktail of hadronic sources (see Sec. IV) and include contribution from charm calculated by PYTHIA using the cross section from Ref. [48] scaled by $N_{coll}$.

(Color online) $e^+e^-$ pair invariant mass distributions in $p$+$p$ (left) and minimum bias Au+Au collisions (right). The $p_T$ ranges are shown in the legend. The solid curves represent the cocktail of hadronic sources (see Sec. IV) and include contribution from charm calculated by PYTHIA using the cross section from Ref. [48] scaled by $N_{coll}$.

(Color online) $e^+e^-$ pair invariant mass distributions in $p$+$p$ (left) and minimum bias Au+Au collisions (right). The $p_T$ ranges are shown in the legend. The solid curves represent the cocktail of hadronic sources (see Sec. IV) and include contribution from charm calculated by PYTHIA using the cross section from Ref. [48] scaled by $N_{coll}$.

(Color online) $e^+e^-$ pair invariant mass distributions in $p$+$p$ (left) and minimum bias Au+Au collisions (right). The $p_T$ ranges are shown in the legend. The solid curves represent the cocktail of hadronic sources (see Sec. IV) and include contribution from charm calculated by PYTHIA using the cross section from Ref. [48] scaled by $N_{coll}$.

(Color online) $e^+e^-$ pair invariant mass distributions in $p$+$p$ (left) and minimum bias Au+Au collisions (right). The $p_T$ ranges are shown in the legend. The solid curves represent the cocktail of hadronic sources (see Sec. IV) and include contribution from charm calculated by PYTHIA using the cross section from Ref. [48] scaled by $N_{coll}$.

(Color online) $e^+e^-$ pair invariant mass distributions in $p$+$p$ (left) and minimum bias Au+Au collisions (right). The $p_T$ ranges are shown in the legend. The solid curves represent the cocktail of hadronic sources (see Sec. IV) and include contribution from charm calculated by PYTHIA using the cross section from Ref. [48] scaled by $N_{coll}$.

(Color online) $e^+e^-$ pair invariant mass distributions in $p$+$p$ (left) and minimum bias Au+Au collisions (right). The $p_T$ ranges are shown in the legend. The solid curves represent the cocktail of hadronic sources (see Sec. IV) and include contribution from charm calculated by PYTHIA using the cross section from Ref. [48] scaled by $N_{coll}$.

(Color online) $e^+e^-$ pair invariant mass distributions in $p$+$p$ (left) and minimum bias Au+Au collisions (right). The $p_T$ ranges are shown in the legend. The solid curves represent the cocktail of hadronic sources (see Sec. IV) and include contribution from charm calculated by PYTHIA using the cross section from Ref. [48] scaled by $N_{coll}$.

(Color online) $e^+e^-$ pair invariant mass distributions in $p$+$p$ (left) and minimum bias Au+Au collisions (right). The $p_T$ ranges are shown in the legend. The solid curves represent the cocktail of hadronic sources (see Sec. IV) and include contribution from charm calculated by PYTHIA using the cross section from Ref. [48] scaled by $N_{coll}$.

(Color online) $e^+e^-$ pair invariant mass distributions in $p$+$p$ (left) and minimum bias Au+Au collisions (right). The $p_T$ ranges are shown in the legend. The solid curves represent the cocktail of hadronic sources (see Sec. IV) and include contribution from charm calculated by PYTHIA using the cross section from Ref. [48] scaled by $N_{coll}$.

(Color online) Electron pair mass distribution for Au+Au (Min.Bias) for $1.0<p_T<1.5$ GeV/$c$. The two component fit is explained in the text. The fit range is $0.12<m_{ee}<0.3$ GeV/$c^2$. The dashed (black) curve at greater $m_{ee}$ shows $f(m_{ee})$ outside of the fit range.

(Color online) Electron pair mass distribution for Au+Au (Min.Bias) for $1.0<p_T<1.5$ GeV/$c$. The two component fit is explained in the text. The fit range is $0.12<m_{ee}<0.3$ GeV/$c^2$. The dashed (black) curve at greater $m_{ee}$ shows $f(m_{ee})$ outside of the fit range.

(Color online) Electron pair mass distribution for Au+Au (Min.Bias) for $1.0<p_T<1.5$ GeV/$c$. The two component fit is explained in the text. The fit range is $0.12<m_{ee}<0.3$ GeV/$c^2$. The dashed (black) curve at greater $m_{ee}$ shows $f(m_{ee})$ outside of the fit range.

(Color online) Electron pair mass distribution for Au+Au (Min.Bias) for $1.0<p_T<1.5$ GeV/$c$. The two component fit is explained in the text. The fit range is $0.12<m_{ee}<0.3$ GeV/$c^2$. The dashed (black) curve at greater $m_{ee}$ shows $f(m_{ee})$ outside of the fit range.

(Color online) Ratio R=(data-cocktail)/$f_{dir}(m_{ee})$ of electron pairs for different $p_T$ bins in Min.Bias Au+Au collisions. The $p_T$ range of each panel is indicated in the figure.

(Color online) The fraction of the direct photon component as a function of $p_T$. The error bars and the error band represent statistical and systematic uncertainties, respectively. The curves are from a NLO pQCD calculation (see text).

(Color online) The fraction of the direct photon component as a function of $p_T$. The error bars and the error band represent statistical and systematic uncertainties, respectively. The curves are from a NLO pQCD calculation (see text).

(Color online) Invariant cross section ($p$+$p$) and invariant yield (Au+Au) of direct photons as a function of $p_T$. The filled points are from this analysis and open points are from [81,82]. The three curves on the $p$+$p$ data represent NLO pQCD calculations, and the dashed curves show a modified power-law fit to the $p$+$p$ data, scaled by $T_{AA}$. The dashed (black) curves are exponential plus the $T_{AA}$ scaled $p$+$p$ fit.

(Color online) Invariant cross section ($p$+$p$) and invariant yield (Au+Au) of direct photons as a function of $p_T$. The filled points are from this analysis and open points are from [81,82]. The three curves on the $p$+$p$ data represent NLO pQCD calculations, and the dashed curves show a modified power-law fit to the $p$+$p$ data, scaled by $T_{AA}$. The dashed (black) curves are exponential plus the $T_{AA}$ scaled $p$+$p$ fit.

(Color online) Invariant cross section ($p$+$p$) and invariant yield (Au+Au) of direct photons as a function of $p_T$. The filled points are from this analysis and open points are from [81,82]. The three curves on the $p$+$p$ data represent NLO pQCD calculations, and the dashed curves show a modified power-law fit to the $p$+$p$ data, scaled by $T_{AA}$. The dashed (black) curves are exponential plus the $T_{AA}$ scaled $p$+$p$ fit.

(Color online) Invariant cross section ($p$+$p$) and invariant yield (Au+Au) of direct photons as a function of $p_T$. The filled points are from this analysis and open points are from [81,82]. The three curves on the $p$+$p$ data represent NLO pQCD calculations, and the dashed curves show a modified power-law fit to the $p$+$p$ data, scaled by $T_{AA}$. The dashed (black) curves are exponential plus the $T_{AA}$ scaled $p$+$p$ fit.

(Color online) The $e^+e^-$ pair invariant mass distributions in minimum bias Au+Au collisions for the low-$p_T$ range. The solid curves represent the cocktail of hadronic sources (see Sec. IV) and include contribution from charm calculated by PYTHIA using the cross section from Ref. [48] scaled by $N_{coll}$.

(Color online) The $e^+e^-$ pair invariant mass distributions in minimum bias Au+Au collisions for the low-$p_T$ range. The solid curves represent the cocktail of hadronic sources (see Sec. IV) and include contribution from charm calculated by PYTHIA using the cross section from Ref. [48] scaled by $N_{coll}$.

(Color online) The $e^+e^-$ pair invariant mass distributions in minimum bias Au+Au collisions for the low-$p_T$ range. The solid curves represent the cocktail of hadronic sources (see Sec. IV) and include contribution from charm calculated by PYTHIA using the cross section from Ref. [48] scaled by $N_{coll}$.

(Color online) The $e^+e^-$ pair invariant mass distributions in minimum bias Au+Au collisions for the low-$p_T$ range. The solid curves represent the cocktail of hadronic sources (see Sec. IV) and include contribution from charm calculated by PYTHIA using the cross section from Ref. [48] scaled by $N_{coll}$.

(Color online) The $e^+e^-$ pair invariant mass distributions in minimum bias Au+Au collisions for the low-$p_T$ range. The solid curves represent the cocktail of hadronic sources (see Sec. IV) and include contribution from charm calculated by PYTHIA using the cross section from Ref. [48] scaled by $N_{coll}$.

(Color online) The $e^+e^-$ pair invariant mass distributions in minimum bias Au+Au collisions for the low-$p_T$ range. The solid curves represent the cocktail of hadronic sources (see Sec. IV) and include contribution from charm calculated by PYTHIA using the cross section from Ref. [48] scaled by $N_{coll}$.

(Color online) The $e^+e^-$ pair invariant mass distributions in minimum bias Au+Au collisions for the low-$p_T$ range. The solid curves represent the cocktail of hadronic sources (see Sec. IV) and include contribution from charm calculated by PYTHIA using the cross section from Ref. [48] scaled by $N_{coll}$.

Ratio of R = (data − cocktail)/$f_{dir}(m_{ee})$ for 0.8 $< p_T <$ 1.0 GeV/c in minimum bias Au$+$Au collisions. The yellow band in each panel shows $\pm1\sigma$ band of a constant fit value to the data points.

Ratio of R = (data − cocktail)/$f_{dir}(m_{ee})$ for 0.6 $< p_T <$ 0.8 GeV/c in minimum bias Au$+$Au collisions. The yellow band in each panel shows $\pm1\sigma$ band of a constant fit value to the data points.

Ratio of R = (data − cocktail)/$f_{dir}(m_{ee})$ for 0.4 $ < p_T <$ 0.6 GeV/c in minimum bias Au$+$Au collisions. The yellow band in each panel shows $\pm1\sigma$ band of a constant fit value to the data points.

(Color online) $p_T$ spectra of $e^+e^-$ pairs in $p$+$p$ (left) and Au+Au (right) collisions for different mass bins, which are fully acceptance corrected. Au+Au spectra are divided by $N_{part}/2$. The solid curves show the expectations from the sum of the hadronic decay cocktail and the contribution from charmed mesons. The dashed curves show the sum of the cocktail and charmed meson contributions plus the contribution from direct photons calculated by converting the photon yield from Fig. 34 to the $e^+e^-$ pair yield using Eqs. (31) and (B14). $p$$+$$p$ collision data shown.

(Color online) $p_T$ spectra of $e^+e^-$ pairs in $p$+$p$ (left) and Au+Au (right) collisions for different mass bins, which are fully acceptance corrected. Au+Au spectra are divided by $N_{part}/2$. The solid curves show the expectations from the sum of the hadronic decay cocktail and the contribution from charmed mesons. The dashed curves show the sum of the cocktail and charmed meson contributions plus the contribution from direct photons calculated by converting the photon yield from Fig. 34 to the $e^+e^-$ pair yield using Eqs. (31) and (B14). $p$$+$$p$ collision data shown.

(Color online) $p_T$ spectra of $e^+e^-$ pairs in $p$+$p$ (left) and Au+Au (right) collisions for different mass bins, which are fully acceptance corrected. Au+Au spectra are divided by $N_{part}/2$. The solid curves show the expectations from the sum of the hadronic decay cocktail and the contribution from charmed mesons. The dashed curves show the sum of the cocktail and charmed meson contributions plus the contribution from direct photons calculated by converting the photon yield from Fig. 34 to the $e^+e^-$ pair yield using Eqs. (31) and (B14). $p$$+$$p$ collision data shown.

(Color online) $p_T$ spectra of $e^+e^-$ pairs in $p$+$p$ (left) and Au+Au (right) collisions for different mass bins, which are fully acceptance corrected. Au+Au spectra are divided by $N_{part}/2$. The solid curves show the expectations from the sum of the hadronic decay cocktail and the contribution from charmed mesons. The dashed curves show the sum of the cocktail and charmed meson contributions plus the contribution from direct photons calculated by converting the photon yield from Fig. 34 to the $e^+e^-$ pair yield using Eqs. (31) and (B14). $p$$+$$p$ collision data shown.

(Color online) $p_T$ spectra of $e^+e^-$ pairs in $p$+$p$ (left) and Au+Au (right) collisions for different mass bins, which are fully acceptance corrected. Au+Au spectra are divided by $N_{part}/2$. The solid curves show the expectations from the sum of the hadronic decay cocktail and the contribution from charmed mesons. The dashed curves show the sum of the cocktail and charmed meson contributions plus the contribution from direct photons calculated by converting the photon yield from Fig. 34 to the $e^+e^-$ pair yield using Eqs. (31) and (B14). $p$$+$$p$ collision data shown.

(Color online) $p_T$ spectra of $e^+e^-$ pairs in $p$+$p$ (left) and Au+Au (right) collisions for different mass bins, which are fully acceptance corrected. Au+Au spectra are divided by $N_{part}/2$. The solid curves show the expectations from the sum of the hadronic decay cocktail and the contribution from charmed mesons. The dashed curves show the sum of the cocktail and charmed meson contributions plus the contribution from direct photons calculated by converting the photon yield from Fig. 34 to the $e^+e^-$ pair yield using Eqs. (31) and (B14). $p$$+$$p$ collision data shown.

(Color online) $p_T$ spectra of $e^+e^-$ pairs in $p$+$p$ (left) and Au+Au (right) collisions for different mass bins, which are fully acceptance corrected. Au+Au spectra are divided by $N_{part}/2$. The solid curves show the expectations from the sum of the hadronic decay cocktail and the contribution from charmed mesons. The dashed curves show the sum of the cocktail and charmed meson contributions plus the contribution from direct photons calculated by converting the photon yield from Fig. 34 to the $e^+e^-$ pair yield using Eqs. (31) and (B14). Au$+$Au collision data shown.

(Color online) $p_T$ spectra of $e^+e^-$ pairs in $p$+$p$ (left) and Au+Au (right) collisions for different mass bins, which are fully acceptance corrected. Au+Au spectra are divided by $N_{part}/2$. The solid curves show the expectations from the sum of the hadronic decay cocktail and the contribution from charmed mesons. The dashed curves show the sum of the cocktail and charmed meson contributions plus the contribution from direct photons calculated by converting the photon yield from Fig. 34 to the $e^+e^-$ pair yield using Eqs. (31) and (B14). Au$+$Au collision data shown.

(Color online) $p_T$ spectra of $e^+e^-$ pairs in $p$+$p$ (left) and Au+Au (right) collisions for different mass bins, which are fully acceptance corrected. Au+Au spectra are divided by $N_{part}/2$. The solid curves show the expectations from the sum of the hadronic decay cocktail and the contribution from charmed mesons. The dashed curves show the sum of the cocktail and charmed meson contributions plus the contribution from direct photons calculated by converting the photon yield from Fig. 34 to the $e^+e^-$ pair yield using Eqs. (31) and (B14). Au$+$Au collision data shown.

(Color online) $p_T$ spectra of $e^+e^-$ pairs in $p$+$p$ (left) and Au+Au (right) collisions for different mass bins, which are fully acceptance corrected. Au+Au spectra are divided by $N_{part}/2$. The solid curves show the expectations from the sum of the hadronic decay cocktail and the contribution from charmed mesons. The dashed curves show the sum of the cocktail and charmed meson contributions plus the contribution from direct photons calculated by converting the photon yield from Fig. 34 to the $e^+e^-$ pair yield using Eqs. (31) and (B14). Au$+$Au collision data shown.

(Color online) $p_T$ spectra of $e^+e^-$ pairs in $p$+$p$ (left) and Au+Au (right) collisions for different mass bins, which are fully acceptance corrected. Au+Au spectra are divided by $N_{part}/2$. The solid curves show the expectations from the sum of the hadronic decay cocktail and the contribution from charmed mesons. The dashed curves show the sum of the cocktail and charmed meson contributions plus the contribution from direct photons calculated by converting the photon yield from Fig. 34 to the $e^+e^-$ pair yield using Eqs. (31) and (B14). Au$+$Au collision data shown.

(Color online) $p_T$ spectra of $e^+e^-$ pairs in $p$+$p$ (left) and Au+Au (right) collisions for different mass bins, which are fully acceptance corrected. Au+Au spectra are divided by $N_{part}/2$. The solid curves show the expectations from the sum of the hadronic decay cocktail and the contribution from charmed mesons. The dashed curves show the sum of the cocktail and charmed meson contributions plus the contribution from direct photons calculated by converting the photon yield from Fig. 34 to the $e^+e^-$ pair yield using Eqs. (31) and (B14). Au$+$Au collision data shown.

(Color online) $p_T$ spectra of $e^+e^-$ pairs in $p$+$p$ (left) and Au+Au (right) collisions for different mass bins, which are fully acceptance corrected. Au+Au spectra are divided by $N_{part}/2$. The solid curves show the expectations from the sum of the hadronic decay cocktail and the contribution from charmed mesons. The dashed curves show the sum of the cocktail and charmed meson contributions plus the contribution from direct photons calculated by converting the photon yield from Fig. 34 to the $e^+e^-$ pair yield using Eqs. (31) and (B14). Au$+$Au collision data shown.

(Color online) $p_T$ spectra of $e^+e^-$ pairs in $p$+$p$ (left) and Au+Au (right) collisions for different mass bins, which are fully acceptance corrected. Au+Au spectra are divided by $N_{part}/2$. The solid curves show the expectations from the sum of the hadronic decay cocktail and the contribution from charmed mesons. The dashed curves show the sum of the cocktail and charmed meson contributions plus the contribution from direct photons calculated by converting the photon yield from Fig. 34 to the $e^+e^-$ pair yield using Eqs. (31) and (B14). Au$+$Au collision data shown.

(Color online) $p_T$ spectra of $e^+e^-$ pairs in $p$+$p$ (left) and Au+Au (right) collisions for different mass bins, which are fully acceptance corrected. Au+Au spectra are divided by $N_{part}/2$. The solid curves show the expectations from the sum of the hadronic decay cocktail and the contribution from charmed mesons. The dashed curves show the sum of the cocktail and charmed meson contributions plus the contribution from direct photons calculated by converting the photon yield from Fig. 34 to the $e^+e^-$ pair yield using Eqs. (31) and (B14). Au$+$Au collision data shown.

(Color online) $p_T$ spectra of $e^+e^-$ pairs in $p$+$p$ (left) and Au+Au (right) collisions for different mass bins, which are fully acceptance corrected. Au+Au spectra are divided by $N_{part}/2$. The solid curves show the expectations from the sum of the hadronic decay cocktail and the contribution from charmed mesons. The dashed curves show the sum of the cocktail and charmed meson contributions plus the contribution from direct photons calculated by converting the photon yield from Fig. 34 to the $e^+e^-$ pair yield using Eqs. (31) and (B14). Au$+$Au collision data shown.

The $m_{T} - m_{0}$ spectrum for the mass range 0.3 < $m_{ee}$ < 0.75 GeV/$c^{2}$ after subtracting contributions from cocktail and charm. The spectrum is fully acceptance corrected. The systematic error band includes the difference in charm yields in this mass range. The spectrum is fit to the sum of two exponential functions which are also shown separately as the dashed and dotted lines. The solid line is the sum.

Local inverse slope of the $m_{T}$ spectra of electron pairs, after subtracting the cocktail and the charm contribution, for different mass bins. The local slope is calculated in different mass ranges, 0 < $m_{T} - m_{0}$ < 0.6 GeV/$c^{2}$ and 0.6 < $m_{T} - m_{0}$ < 2.5 GeV/$c^{2}$. The solid and dashed lines show the local slope of the cocktail for the corresponding mass ranges.

Local inverse slope of the $m_{T}$ spectra of electron pairs, after subtracting the cocktail and the charm contribution, for different mass bins. The local slope is calculated in different mass ranges, 0 < $m_{T} - m_{0}$ < 0.6 GeV/$c^{2}$ and 0.6 < $m_{T} - m_{0}$ < 2.5 GeV/$c^{2}$. The solid and dashed lines show the local slope of the cocktail for the corresponding mass ranges.

Invariant mass spectrum of $e^{+}e^{-}$ pairs in MB Au+Au compared to predictions from Ralf Rapp. Invariant mass spectra of $e^{+}e^{-}$ pairs in Min. Bias Au + Au collisions in the IMR. (top left) The data are compared to the sum of cocktail+charm. The data are also compared to the sum of cocktail+charm and partonic contributions from different models. The calculations are from (center) Rapp and van Hees [15, 18, 83] and (right) Dusling and Zahed [19, 84, 85]. The partonic yields (PY) have been added to the two scenarios for charmed mesons decays, i.e. (i) $PYTHIA$ and (ii) random $c\bar{c}$ correlation.

Invariant mass spectrum of $e^{+}e^{-}$ pairs in MB Au+Au compared to predictions from Kevin Dusling. Invariant mass spectra of $e^{+}e^{-}$ pairs in Min. Bias Au + Au collisions in the IMR. (top left) The data are compared to the sum of cocktail+charm. The data are also compared to the sum of cocktail+charm and partonic contributions from different models. The calculations are from (center) Rapp and van Hees [15, 18, 83] and (right) Dusling and Zahed [19, 84, 85]. The partonic yields (PY) have been added to the two scenarios for charmed mesons decays, i.e. (i) $PYTHIA$ and (ii) random $c\bar{c}$ correlation.

Invariant mass spectrum of $e^{+}e^{-}$ pairs in MB Au+Au compared to predictions from Elena Bratkovskaya. Invariant mass spectra of $e^{+}e^{-}$ pairs in Min. Bias Au + Au collisions in the IMR. (top left) The data are compared to the sum of cocktail+charm. The data are also compared to the sum of cocktail+charm and partonic contributions from different models. The calculations are from (center) Rapp and van Hees [15, 18, 83] and (right) Dusling and Zahed [19, 84, 85]. The partonic yields (PY) have been added to the two scenarios for charmed mesons decays, i.e. (i) $PYTHIA$ and (ii) random $c\bar{c}$ correlation.

invariant mass spectrum of e+e- pairs in MB Au+Au compared to predictions from Ralf Rapp.

invariant mass spectrum of e+e- pairs in MB Au+Au compared to predictions from Ralf Rapp.

Invariant mass spectrum of $e^{+}e^{-}$ pairs in MB Au+Au for different $p_{T}$ ranges compared to predictions from Kevin Dusling. Invariant mass spectra of $e^{+}e^{-}$ pairs in Au + Au collisions in the LMR. The data are compared to the sum of cocktail+charm (top left). The data are also compared to the sum of cocktail+charm and hadronic+partonic contributions from different models. The calculations are from The calculations are from (top right) Rapp and van Hees [15, 18, 83], (bottom right) Dusling and Zahed [19, 84, 85], and Cassing and Bratkovskaya [20, 27, 86, 87].

Invariant mass spectrum of $e^{+}e^{-}$ pairs in MB Au+Au for different $p_{T}$ ranges compared to predictions from Kevin Dusling. Invariant mass spectra of $e^{+}e^{-}$ pairs in Au + Au collisions in the LMR. The data are compared to the sum of cocktail+charm (top left). The data are also compared to the sum of cocktail+charm and hadronic+partonic contributions from different models. The calculations are from The calculations are from (top right) Rapp and van Hees [15, 18, 83], (bottom right) Dusling and Zahed [19, 84, 85], and Cassing and Bratkovskaya [20, 27, 86, 87].

Invariant mass spectrum of $e^{+}e^{-}$ pairs in MB Au+Au for different $p_{T}$ ranges compared to predictions from Elena Bratkovskaya. Invariant mass spectra of $e^{+}e^{-}$ pairs in Au + Au collisions in the LMR. The data are compared to the sum of cocktail+charm (top left). The data are also compared to the sum of cocktail+charm and hadronic+partonic contributions from different models. The calculations are from The calculations are from (top right) Rapp and van Hees [15, 18, 83], (bottom right) Dusling and Zahed [19, 84, 85], and Cassing and Bratkovskaya [20, 27, 86, 87].

Invariant mass spectrum of $e^{+}e^{-}$ pairs in MB Au+Au for different $p_{T}$ ranges compared to predictions from Elena Bratkovskaya. Invariant mass spectra of $e^{+}e^{-}$ pairs in Au + Au collisions in the LMR. The data are compared to the sum of cocktail+charm (top left). The data are also compared to the sum of cocktail+charm and hadronic+partonic contributions from different models. The calculations are from The calculations are from (top right) Rapp and van Hees [15, 18, 83], (bottom right) Dusling and Zahed [19, 84, 85], and Cassing and Bratkovskaya [20, 27, 86, 87].

Invariant mass spectrum of $e^{+}e^{-}$ pairs in MB Au+Au for different $p_{T}$ ranges compared to predictions from Ralf Rapp (0<$p_{T}$<0.5 GeV/$c$). Invariant mass spectra of $e^{+}e^{-}$ pairs in Min. Bias Au + Au collisions for different $p_{T}$ windows compared to the expectations from the calculations of Rapp and van Hees [15, 18, 83], separately showing the partonic and the hadronic yields and the different scenarios for the $\rho$ spectral function, namely “Hadron Many Body Theory” (HMBT) and “Dropping Mass” (DM). The calculations have been added to the cocktail of hadronic decays (where the contribution of the freeze-out $\rho$ meson is subtracted) and charmed meson decays products.

Invariant mass spectrum of $e^{+}e^{-}$ pairs in MB Au+Au for different $p_{T}$ ranges compared to predictions from Ralf Rapp (0<$p_{T}$<0.5 GeV/$c$). Invariant mass spectra of $e^{+}e^{-}$ pairs in Min. Bias Au + Au collisions for different $p_{T}$ windows compared to the expectations from the calculations of Rapp and van Hees [15, 18, 83], separately showing the partonic and the hadronic yields and the different scenarios for the $\rho$ spectral function, namely “Hadron Many Body Theory” (HMBT) and “Dropping Mass” (DM). The calculations have been added to the cocktail of hadronic decays (where the contribution of the freeze-out $\rho$ meson is subtracted) and charmed meson decays products.

Invariant mass spectrum of $e^{+}e^{-}$ pairs in MB Au+Au for different $p_{T}$ ranges compared to predictions from Ralf Rapp (0<$p_{T}$<0.5 GeV/$c$). Invariant mass spectra of $e^{+}e^{-}$ pairs in Min. Bias Au + Au collisions for different $p_{T}$ windows compared to the expectations from the calculations of Rapp and van Hees [15, 18, 83], separately showing the partonic and the hadronic yields and the different scenarios for the $\rho$ spectral function, namely “Hadron Many Body Theory” (HMBT) and “Dropping Mass” (DM). The calculations have been added to the cocktail of hadronic decays (where the contribution of the freeze-out $\rho$ meson is subtracted) and charmed meson decays products.

Invariant mass spectrum of $e^{+}e^{-}$ pairs in MB Au+Au for different $p_{T}$ ranges compared to predictions from Ralf Rapp (0<$p_{T}$<0.5 GeV/$c$). Invariant mass spectra of $e^{+}e^{-}$ pairs in Min. Bias Au + Au collisions for different $p_{T}$ windows compared to the expectations from the calculations of Rapp and van Hees [15, 18, 83], separately showing the partonic and the hadronic yields and the different scenarios for the $\rho$ spectral function, namely “Hadron Many Body Theory” (HMBT) and “Dropping Mass” (DM). The calculations have been added to the cocktail of hadronic decays (where the contribution of the freeze-out $\rho$ meson is subtracted) and charmed meson decays products.

Invariant mass spectrum of $e^{+}e^{-}$ pairs in MB Au+Au for different $p_{T}$ ranges compared to predictions from Kevin Dusling (0<$p_{T}$<0.5 GeV/$c$). Invariant mass spectra of $e^{+}e^{-}$ pairs in Min. Bias Au + Au collisions for different $p_{T}$ windows compared to the expectations from the calculations of Dusling and Zahed [19, 84, 85], separately showing the partonic and the hadronic yields. The calculations have been added to the cocktail of hadronic decays (where the contribution of the freeze-out $\rho$ meson is subtracted) and charmed meson decays products.

Invariant mass spectrum of $e^{+}e^{-}$ pairs in MB Au+Au for different $p_{T}$ ranges compared to predictions from Kevin Dusling (0<$p_{T}$<0.5 GeV/$c$). Invariant mass spectra of $e^{+}e^{-}$ pairs in Min. Bias Au + Au collisions for different $p_{T}$ windows compared to the expectations from the calculations of Dusling and Zahed [19, 84, 85], separately showing the partonic and the hadronic yields. The calculations have been added to the cocktail of hadronic decays (where the contribution of the freeze-out $\rho$ meson is subtracted) and charmed meson decays products.

Invariant mass spectrum of $e^{+}e^{-}$ pairs in MB Au+Au for different $p_{T}$ ranges compared to predictions from Kevin Dusling (0.5<$p_{T}$<1.0 GeV/$c$). Invariant mass spectra of $e^{+}e^{−}$ pairs in Min. Bias Au + Au collisions for different $p_{T}$ windows compared to the expectations from the calculations of Dusling and Zahed [19, 84, 85], separately showing the partonic and the hadronic yields. The calculations have been added to the cocktail of hadronic decays (where the contribution of the freeze-out $\rho$ meson is subtracted) and charmed meson decays products.

Invariant mass spectrum of $e^{+}e^{-}$ pairs in MB Au+Au for different $p_{T}$ ranges compared to predictions from Kevin Dusling (0.5<$p_{T}$<1.0 GeV/$c$). Invariant mass spectra of $e^{+}e^{−}$ pairs in Min. Bias Au + Au collisions for different $p_{T}$ windows compared to the expectations from the calculations of Dusling and Zahed [19, 84, 85], separately showing the partonic and the hadronic yields. The calculations have been added to the cocktail of hadronic decays (where the contribution of the freeze-out $\rho$ meson is subtracted) and charmed meson decays products.

Invariant mass spectrum of $e^{+}e^{-}$ pairs in MB Au+Au for different $p_{T}$ ranges compared to predictions from Kevin Dusling 1.0<$p_{T}$<1.5 GeV/$c$. Invariant mass spectra of $e^{+}e^{−}$ pairs in Min. Bias Au + Au collisions for different $p_{T}$ windows compared to the expectations from the calculations of Dusling and Zahed [19, 84, 85], separately showing the partonic and the hadronic yields. The calculations have been added to the cocktail of hadronic decays (where the contribution of the freeze-out $\rho$ meson is subtracted) and charmed meson decays products.

Invariant mass spectrum of $e^{+}e^{-}$ pairs in MB Au+Au for different $p_{T}$ ranges compared to predictions from Kevin Dusling 1.0<$p_{T}$<1.5 GeV/$c$. Invariant mass spectra of $e^{+}e^{−}$ pairs in Min. Bias Au + Au collisions for different $p_{T}$ windows compared to the expectations from the calculations of Dusling and Zahed [19, 84, 85], separately showing the partonic and the hadronic yields. The calculations have been added to the cocktail of hadronic decays (where the contribution of the freeze-out $\rho$ meson is subtracted) and charmed meson decays products.

Invariant mass spectrum of $e^{+}e^{-}$ pairs in MB Au+Au for different $p_{T}$ ranges compared to predictions from Kevin Dusling (1.5<$p_{T}$<2.0 GeV/$c$). Invariant mass spectra of $e^{+}e^{−}$ pairs in Min. Bias Au + Au collisions for different $p_{T}$ windows compared to the expectations from the calculations of Dusling and Zahed [19, 84, 85], separately showing the partonic and the hadronic yields. The calculations have been added to the cocktail of hadronic decays (where the contribution of the freeze-out $\rho$ meson is subtracted) and charmed meson decays products.

Invariant mass spectrum of $e^{+}e^{-}$ pairs in MB Au+Au for different $p_{T}$ ranges compared to predictions from Kevin Dusling (1.5<$p_{T}$<2.0 GeV/$c$). Invariant mass spectra of $e^{+}e^{−}$ pairs in Min. Bias Au + Au collisions for different $p_{T}$ windows compared to the expectations from the calculations of Dusling and Zahed [19, 84, 85], separately showing the partonic and the hadronic yields. The calculations have been added to the cocktail of hadronic decays (where the contribution of the freeze-out $\rho$ meson is subtracted) and charmed meson decays products.

Invariant mass spectrum of $e^{+}e^{-}$ pairs in MB Au+Au for different $p_{T}$ ranges compared to predictions from Elena Bratkovskaya (0<$p_{T}$<0.5 GeV/$c$). Invariant mass spectra of $e^{+}e^{−}$ pairs in Min. Bias Au + Au collisions for different $p_{T}$ windows collisions compared to the expectations from the calculations of Cassing and Bratkovskaya [20, 27, 86, 87], separately showing the partonic and the hadronic yields calculated with different implementations of the $\rho$ spectral function, namely according to collisional broadening, with or without a dropping mass scenario. The calculations which include the dropping mass scenario have been added to the cocktail of hadronic decays (which is calculated by the HSD model itself) and charmed meson decays products.

Invariant mass spectrum of $e^{+}e^{-}$ pairs in MB Au+Au for different $p_{T}$ ranges compared to predictions from Elena Bratkovskaya (0<$p_{T}$<0.5 GeV/$c$). Invariant mass spectra of $e^{+}e^{−}$ pairs in Min. Bias Au + Au collisions for different $p_{T}$ windows collisions compared to the expectations from the calculations of Cassing and Bratkovskaya [20, 27, 86, 87], separately showing the partonic and the hadronic yields calculated with different implementations of the $\rho$ spectral function, namely according to collisional broadening, with or without a dropping mass scenario. The calculations which include the dropping mass scenario have been added to the cocktail of hadronic decays (which is calculated by the HSD model itself) and charmed meson decays products.

Invariant mass spectrum of $e^{+}e^{-}$ pairs in MB Au+Au for different $p_{T}$ ranges compared to predictions from Elena Bratkovskaya (0<$p_{T}$<0.5 GeV/$c$). Invariant mass spectra of $e^{+}e^{−}$ pairs in Min. Bias Au + Au collisions for different $p_{T}$ windows collisions compared to the expectations from the calculations of Cassing and Bratkovskaya [20, 27, 86, 87], separately showing the partonic and the hadronic yields calculated with different implementations of the $\rho$ spectral function, namely according to collisional broadening, with or without a dropping mass scenario. The calculations which include the dropping mass scenario have been added to the cocktail of hadronic decays (which is calculated by the HSD model itself) and charmed meson decays products.

Invariant mass spectrum of $e^{+}e^{-}$ pairs in MB Au+Au for different $p_{T}$ ranges compared to predictions from Elena Bratkovskaya (0.5<$p_{T}$<1.0 GeV/$c$). Invariant mass spectra of $e^{+}e^{−}$ pairs in Min. Bias Au + Au collisions for different $p_{T}$ windows collisions compared to the expectations from the calculations of Cassing and Bratkovskaya [20, 27, 86, 87], separately showing the partonic and the hadronic yields calculated with different implementations of the $\rho$ spectral function, namely according to collisional broadening, with or without a dropping mass scenario. The calculations which include the dropping mass scenario have been added to the cocktail of hadronic decays (which is calculated by the HSD model itself) and charmed meson decays products.

Invariant mass spectrum of $e^{+}e^{-}$ pairs in MB Au+Au for different $p_{T}$ ranges compared to predictions from Elena Bratkovskaya (0.5<$p_{T}$<1.0 GeV/$c$). Invariant mass spectra of $e^{+}e^{−}$ pairs in Min. Bias Au + Au collisions for different $p_{T}$ windows collisions compared to the expectations from the calculations of Cassing and Bratkovskaya [20, 27, 86, 87], separately showing the partonic and the hadronic yields calculated with different implementations of the $\rho$ spectral function, namely according to collisional broadening, with or without a dropping mass scenario. The calculations which include the dropping mass scenario have been added to the cocktail of hadronic decays (which is calculated by the HSD model itself) and charmed meson decays products.

Invariant mass spectrum of $e^{+}e^{-}$ pairs in MB Au+Au for different $p_{T}$ ranges compared to predictions from Elena Bratkovskaya (0.5<$p_{T}$<1.0 GeV/$c$). Invariant mass spectra of $e^{+}e^{−}$ pairs in Min. Bias Au + Au collisions for different $p_{T}$ windows collisions compared to the expectations from the calculations of Cassing and Bratkovskaya [20, 27, 86, 87], separately showing the partonic and the hadronic yields calculated with different implementations of the $\rho$ spectral function, namely according to collisional broadening, with or without a dropping mass scenario. The calculations which include the dropping mass scenario have been added to the cocktail of hadronic decays (which is calculated by the HSD model itself) and charmed meson decays products.

Invariant mass spectrum of $e^{+}e^{-}$ pairs in MB Au+Au for different $p_{T}$ ranges compared to predictions from Elena Bratkovskaya (1.0<$p_{T}$<1.5 GeV/$c$). Invariant mass spectra of $e^{+}e^{−}$ pairs in Min. Bias Au + Au collisions for different $p_{T}$ windows collisions compared to the expectations from the calculations of Cassing and Bratkovskaya [20, 27, 86, 87], separately showing the partonic and the hadronic yields calculated with different implementations of the $\rho$ spectral function, namely according to collisional broadening, with or without a dropping mass scenario. The calculations which include the dropping mass scenario have been added to the cocktail of hadronic decays (which is calculated by the HSD model itself) and charmed meson decays products.

Invariant mass spectrum of $e^{+}e^{-}$ pairs in MB Au+Au for different $p_{T}$ ranges compared to predictions from Elena Bratkovskaya (1.0<$p_{T}$<1.5 GeV/$c$). Invariant mass spectra of $e^{+}e^{−}$ pairs in Min. Bias Au + Au collisions for different $p_{T}$ windows collisions compared to the expectations from the calculations of Cassing and Bratkovskaya [20, 27, 86, 87], separately showing the partonic and the hadronic yields calculated with different implementations of the $\rho$ spectral function, namely according to collisional broadening, with or without a dropping mass scenario. The calculations which include the dropping mass scenario have been added to the cocktail of hadronic decays (which is calculated by the HSD model itself) and charmed meson decays products.

Invariant mass spectrum of $e^{+}e^{-}$ pairs in MB Au+Au for different $p_{T}$ ranges compared to predictions from Elena Bratkovskaya (1.5<$p_{T}$<2.0 GeV/$c$). Invariant mass spectra of $e^{+}e^{−}$ pairs in Min. Bias Au + Au collisions for different $p_{T}$ windows collisions compared to the expectations from the calculations of Cassing and Bratkovskaya [20, 27, 86, 87], separately showing the partonic and the hadronic yields calculated with different implementations of the $\rho$ spectral function, namely according to collisional broadening, with or without a dropping mass scenario. The calculations which include the dropping mass scenario have been added to the cocktail of hadronic decays (which is calculated by the HSD model itself) and charmed meson decays products.

Invariant mass spectrum of $e^{+}e^{-}$ pairs in MB Au+Au for different $p_{T}$ ranges compared to predictions from Elena Bratkovskaya (1.5<$p_{T}$<2.0 GeV/$c$). Invariant mass spectra of $e^{+}e^{−}$ pairs in Min. Bias Au + Au collisions for different $p_{T}$ windows collisions compared to the expectations from the calculations of Cassing and Bratkovskaya [20, 27, 86, 87], separately showing the partonic and the hadronic yields calculated with different implementations of the $\rho$ spectral function, namely according to collisional broadening, with or without a dropping mass scenario. The calculations which include the dropping mass scenario have been added to the cocktail of hadronic decays (which is calculated by the HSD model itself) and charmed meson decays products.

Subtracted $p_{T}$ spectrum in 300-750 compared to calculations from Ralf Rapp. $p_{T}$ spectra of $e^{+}e^{−}$ pairs for 0.3 < $m_{ee}$ < 0.75 GeV/$c^{2}$ in Min. Bias Au + Au collisions compared to the expectations from the calculations of respectively R. Rapp and van Hees [15, 18, 83], Dusling and Zahed [19, 84, 85], Cassing and Bratkovskaya [20, 27, 86, 87]. The spectra are fully acceptance corrected. The curves show separately partonic and hadronic yields. For the curves of Rapp and van Hees [15, 18, 83] the two scenarios: Hadron Many Body Theory (HMBT) and Dropping Mass (DM) are shown. The sum is calculated with HMBT. The calculations are compared to the data from which the contributions of the cocktail of hadronic decays and charmed meson decays have been subtracted.

Subtracted $p_{T}$ spectrum in 300-750 compared to calculations from Kevin Dusling. $p_{T}$ spectra of $e^{+}e^{−}$ pairs for 0.3 < $m_{ee}$ < 0.75 GeV/$c^{2}$ in Min. Bias Au + Au collisions compared to the expectations from the calculations of respectively R. Rapp and van Hees [15, 18, 83], Dusling and Zahed [19, 84, 85], Cassing and Bratkovskaya [20, 27, 86, 87]. The spectra are fully acceptance corrected. The curves show separately partonic and hadronic yields. For the curves of Rapp and van Hees [15, 18, 83] the two scenarios: Hadron Many Body Theory (HMBT) and Dropping Mass (DM) are shown. The sum is calculated with HMBT. The calculations are compared to the data from which the contributions of the cocktail of hadronic decays and charmed meson decays have been subtracted.

Subtracted $p_{T}$ spectrum in 300-750 compared to calculations from Elena Bratkovskaya. $p_{T}$ spectra of $e^{+}e^{−}$ pairs for 0.3 < $m_{ee}$ < 0.75 GeV/$c^{2}$ in Min. Bias Au + Au collisions compared to the expectations from the calculations of respectively R. Rapp and van Hees [15, 18, 83], Dusling and Zahed [19, 84, 85], Cassing and Bratkovskaya [20, 27, 86, 87]. The spectra are fully acceptance corrected. The curves show separately partonic and hadronic yields. For the curves of Rapp and van Hees [15, 18, 83] the two scenarios: Hadron Many Body Theory (HMBT) and Dropping Mass (DM) are shown. The sum is calculated with HMBT. The calculations are compared to the data from which the contributions of the cocktail of hadronic decays and charmed meson decays have been subtracted.

Subtracted $p_{T}$ spectrum in 300-750 compared to calculations from Elena Bratkovskaya. $p_{T}$ spectra of $e^{+}e^{−}$ pairs for 0.3 < $m_{ee}$ < 0.75 GeV/$c^{2}$ in Min. Bias Au + Au collisions compared to the expectations from the calculations of respectively R. Rapp and van Hees [15, 18, 83], Dusling and Zahed [19, 84, 85], Cassing and Bratkovskaya [20, 27, 86, 87]. The spectra are fully acceptance corrected. The curves show separately partonic and hadronic yields. For the curves of Rapp and van Hees [15, 18, 83] the two scenarios: Hadron Many Body Theory (HMBT) and Dropping Mass (DM) are shown. The sum is calculated with HMBT. The calculations are compared to the data from which the contributions of the cocktail of hadronic decays and charmed meson decays have been subtracted.

Measured coherent J/psi cross section for the XNXN forward class. Note that for each rapidity range the 0n0n uncertainty related to migrations is preceded by a ∓, while the other neutron classes have a ±; this means that these uncertainties are anti-correlated.

Photonuclear cross sections extracted from the UPC measurements.

Nuclear suppression factor extracted from the UPC measurements.

High $p_T$ ($p_T$ > 6 GeV/$c$) integrated $v_2$ vs $N_{part}$ for (a) $\pi^0$ (b) inclusive photon, and (c) direct photon. Results are shown with both reaction plane detectors$:$ (solid black circles) BBC and (solid red squares) RXN. Each point represents a 10% wide centrality bin from 60%-0%.

The $p_T$ spectra of electrons from the decays of open heavy flavor hadrons produced in Cu+Cu collisions, separated by centrality.

The $p_T$ spectra of electrons from the decays of open heavy flavor hadrons produced in Cu+Cu collisions, separated by centrality.

The $p_T$ spectra of electrons from the decays of open heavy flavor hadrons produced in Cu+Cu collisions, separated by centrality.

The nuclear modification factor for MB (0%–94%).

The nuclear modification factor for 0%–10%.

The nuclear modification factor for 0%–20%.

The nuclear modification factor for 20%–40%.

The nuclear modification factor for 40%–60%.

The nuclear modification factor for 60%–94%.

The nuclear modification factors for $e_{HF}$ at midrapidity in Cu+Cu collisions plotted as a function of $\langle N_{coll} \rangle$.

The nuclear modification factors for $e_{HF}$ at midrapidity in Cu+Cu collisions plotted as a function of $\langle N_{coll} \rangle$.

The nuclear modification factors for $e_{HF}$ at midrapidity in Cu+Cu collisions plotted as a function of $\langle N_{part} \rangle$.

The nuclear modification factors for $e_{HF}$ at midrapidity in Cu+Cu collisions plotted as a function of $\langle N_{part} \rangle$.

The relative nuclear-modification factor for U+U and Au+Au as a function of collision centrality.

The ratio of the invariant yield for U+U to that for Au+Au, as a function of collision centrality.

Figure 4a

Figure 4b

Figure 4c

Figure 5a

Figure 5b

Figure 5c

Figure 7

Figure 8a

Figure 8b

Figure 8c

Figure 9a

Figure 9b

Figure 9c

Figure 10a

Figure 10b

Figure 10c

$P_{\rm H}$ measurements are shown as a function of hyperon $p_{\rm T}$ at $\sqrt{s_{\rm NN}}$=19.6 and 27 GeV. Statistical uncertainties are represented as lines while systematic uncertainties are represented as boxes. There is no observed dependence of $P_{\rm H}$ on $p_{\rm T}$ at $\sqrt{s_{\rm NN}}$=19.6 or 27 GeV, consistent with previous observations.

$P_{\rm H}$ measurements are shown as a function of hyperon $y$ at $\sqrt{s_{\rm NN}}$=19.6 and 27 GeV. Statistical uncertainties are represented as lines while systematic uncertainties are represented as boxes. The data set at $\sqrt{s_{\rm NN}}$=19.6 GeV takes advantage of STAR upgrades to reach larger $|y|$.

HFE (electrons from semileptonic decays of heavy-flavor hadrons) invariant yields in different centrality intervals of Au+Au collisions at $\sqrt{s_{\rm NN}}=200$ GeV. The vertical bars and the boxes represent statistical and systematic uncertainties, respectively. The horizontal bars indicate the bin width.

HFE (electrons from semileptonic decays of heavy-flavor hadrons) $R_{\rm AA}$ (red circles) as a function of $p_{\rm T}$ in different centrality intervals of Au+Au collisions at $\sqrt{s_{\rm NN}}=200$ GeV, compared with STAR (yellow stars) and PHENIX (green squares) published results, and Duke ((modified Langevin transport model, blue line) and PHSD (parton-hadron-string dynamics model, orange line) model calculations. Vertical bars and boxes around data points represent combined statistical and systematic uncertainties from both Au+Au and $p$+$p$ measurements, respectively. Boxes at unity show the global uncertainties, which for this analysis include the 8% global uncertainty on $p$+$p$ reference and the $N_{\rm coll}$ uncertainties. The left box is for PHENIX and the right one for STAR.

HFE (electrons from semileptonic decays of heavy-flavor hadrons) $R_{\rm AA}$ (red circles) as a function of $N_{\rm part}$ in Au+Au collisions at $\sqrt{s_{\rm NN}}=200$ GeV, compared with PHENIX measurements (green squares), and Duke (blue line) and PHSD (orange line) model calculations. Vertical bars and boxes around data points represent statistical and systematic uncertainties from Au+Au measurements, respectively. The gray band represents the $N_{\rm coll}$ uncertainties. The boxes at unity show the global uncertainties including the total uncertainties of the $p$+$p$ reference.

Lévy scale parameter $R$ versus average $m_T$ of the pair. The graphical representation of statistical and systematic uncertainties is the same as in Fig. 4.

Lévy index parameter $\alpha$ versus average $m_T$ of the pair. Statistical and systematic uncertainties are indicated similarly to Fig. 4. The horizontal line, %\aplha$ = 1.207, represents the 0%-30% centrality average value of %\aplha$.

Inverse square of the Lévy scale parameter 1/$R^2$ versus average $m_T$ of the pair. Statistical and systematic uncertainties shown as bars and boxes, respectively.

New scale parameter $\hat{R}$ versus average $m_T$ of the pair, with a linear fit. Statistical and systematic uncertainties shown as bars and boxes, respectively.

Normalized correlation strength parameter $\lambda/\lambda_{max}$ versus average $m_T$ of the pair. The data are compared with parameter scans from Refs. [58, 59] using different values of in-medium $\eta^{\prime}$ mass $m^{*}_{\eta^{\prime}}$ and slope parameter $B^{-1}_{\eta^{\prime}}$ . A best fit with Eq. (63) and the resulting $H$ and $\sigma$ parameters are also shown.

Normalized correlation strength parameter $\lambda/\lambda_{max}$ versus average $m_T$ of the pair. The data are compared with parameter scans from Refs. [58, 59] using different values of in-medium $\eta^{\prime}$ mass $m^{*}_{\eta^{\prime}}$ and slope parameter $B^{-1}_{\eta^{\prime}}$ . A best fit with Eq. (63) and the resulting $H$ and $\sigma$ parameters are also shown.

Normalized correlation strength parameter $\lambda/\lambda_{max}$ versus average $m_T$ of the pair. The data are compared with parameter scans from Refs. [58, 59] using different values of in-medium $\eta^{\prime}$ mass $m^{*}_{\eta^{\prime}}$ and slope parameter $B^{-1}_{\eta^{\prime}}$ . A best fit with Eq. (63) and the resulting $H$ and $\sigma$ parameters are also shown.

Normalized correlation strength parameter $\lambda/\lambda_{max}$ versus average $m_T$ of the pair. The data are compared with parameter scans from Refs. [58, 59] using different values of in-medium $\eta^{\prime}$ mass $m^{*}_{\eta^{\prime}}$ and slope parameter $B^{-1}_{\eta^{\prime}}$ . A best fit with Eq. (63) and the resulting $H$ and $\sigma$ parameters are also shown.

Normalized correlation strength parameter $\lambda/\lambda_{max}$ versus average $m_T$ of the pair. The data are compared with parameter scans from Refs. [58, 59] using different values of in-medium $\eta^{\prime}$ mass $m^{*}_{\eta^{\prime}}$ and slope parameter $B^{-1}_{\eta^{\prime}}$ . A best fit with Eq. (63) and the resulting $H$ and $\sigma$ parameters are also shown.

Dependence of $\overline{p}$-$\overline{d}$ correlation on pseudorapidity acceptance of $\overline{p}$ and $\overline{d}$ selection in Pb-Pb collisions at $\sqrt{s_\mathrm{NN}}$ = 5.02 TeV. Results are for 20.0--30.0$\%$ collision centrality.

Dependence of $\overline{p}$-$\overline{d}$ correlation on pseudorapidity acceptance of $\overline{p}$ and $\overline{d}$ selection in Pb-Pb collisions at $\sqrt{s_\mathrm{NN}}$ = 5.02 TeV. Results are for 50.0--60.0$\%$ collision centrality.

Nuclear modification factor $R_{AA}$ in 10% most central $Au$+$Au$ collisions.

Nuclear modification factor $R_{dA}$ for ($h^+$+$h^-$)/2 in minimum bias $d$+$Au$.

$R_{AA}$ in the 10% most central $Au$+$Au$ collisions.

$R_{dA}$ for ($h^+$+$h^-$)/2 measurement in $d$+$Au$.

$R_{dA}$ for $\pi^0$ measurement in $d$+$Au$.

Ratio of $R_{dA}$ in minimum bias $d$+$Au$ and $R_{pA}$ from neutron tagged $d$+$Au$ collisions for ($h^+$+$h^-$)/2.

Ratio of $R_{dA}$ in minimum bias $d$+$Au$ and $R_{pA}$ from neutron tagged $d$+$Au$ collisions for $\pi^0$

HBT parameters of charge-combined kaon pairs, shown as value $\pm$ statistical uncertainty [absolute value] $\pm$ systematic uncertainty [%] for the centrality bins shown in Fig. 3.

$R_o/R_s$ of positive pion pairs shown as value $\pm$ statistical uncertainty [absolute value] $\pm$ systematic uncertainty [%] for the centrality bins shown in Fig. 4.

$R_o/R_s$ of negative pion pairs shown as value $\pm$ statistical uncertainty [absolute value] $\pm$ systematic uncertainty [%] for the centrality bins shown in Fig. 4.

$R_o/R_s$ of charge-combined kaon pairs shown as value $\pm$ statistical uncertainty [absolute value] $\pm$ systematic uncertainty [%] for the centrality bins shown in Fig. 4.

Azimuthal angle dependence of HBT radii of charged kaons, shown as value $\pm$ statistical uncertainty $\pm$ systematic uncertainty for the centrality bins shown in Fig. 8.

Azimuthal angle dependence of HBT radii of charged pions, shown as value $\pm$ statistical uncertainty $\pm$ systematic uncertainty for the 0%–20% and 20%–60% centrality bins.

Oscillation amplitudes relative to the event plane for charged pions, shown as value $\pm$ statistical uncertainty $\pm$ systematic uncertainty for the 0%–20% and 20%–60% centrality bins.

Oscillation amplitudes relative to the event plane for charged kaons, shown as value $\pm$ statistical uncertainty $\pm$ systematic uncertainty for the 0%–20% and 20%–60% centrality bins.

Summary of extracted parameters in the blast-wave model fit for two fitting conditions. The $T_f$ and $\rho_0$ parameters were obtained by fits to $p_T$ spectra, and the other parameters were obtained by a simultaneous fit to $v_2$ and the HBT radii.

Charged hadron $v_2$($p_T$) in Au+Au collisions at $\sqrt{s_{NN}}$ = 200 GeV from the two-particle cumulant method, the BBC event plane, and the ZDC-SMD event plane for the indicated centralities.

Charged hadron $v_2$($p_T$) in Au+Au collisions at $\sqrt{s_{NN}}$ = 200 GeV from the two-particle cumulant method, the BBC event plane, and the ZDC-SMD event plane for the indicated centralities.

Charged hadron $v_2$($p_T$) in Au+Au collisions at $\sqrt{s_{NN}}$ = 200 GeV from the two-particle cumulant method, the BBC event plane, and the ZDC-SMD event plane for the indicated centralities.

Charged hadron $v_2$($p_T$) in Au+Au collisions at $\sqrt{s_{NN}}$ = 200 GeV from the two-particle cumulant method, the BBC event plane, and the ZDC-SMD event plane for the indicated centralities.

Charged hadron $v_2$($p_T$) in Au+Au collisions at $\sqrt{s_{NN}}$ = 200 GeV from the two-particle cumulant method, the BBC event plane, and the ZDC-SMD event plane for the indicated centralities.

Charged hadron $v_2$($p_T$) in Au+Au collisions at $\sqrt{s_{NN}}$ = 200 GeV from the two-particle cumulant method, the BBC event plane, and the ZDC-SMD event plane for the indicated centralities.

Charged hadron $v_2$($p_T$) in Au+Au collisions at $\sqrt{s_{NN}}$ = 200 GeV from the two-particle cumulant method, the BBC event plane, and the ZDC-SMD event plane for the indicated centralities.

Charged hadron $v_2$($p_T$) in Au+Au collisions at $\sqrt{s_{NN}}$ = 200 GeV from the two-particle cumulant method, the BBC event plane, and the ZDC-SMD event plane for the indicated centralities.

Charged hadron $v_2$($p_T$) in Au+Au collisions at $\sqrt{s_{NN}}$ = 200 GeV from the two-particle cumulant method, the BBC event plane, and the ZDC-SMD event plane for the indicated centralities.

Charged hadron $v_2$($p_T$) in Au+Au collisions at $\sqrt{s_{NN}}$ = 200 GeV from the two-particle cumulant method, the BBC event plane, and the ZDC-SMD event plane for the indicated centralities.

Charged hadron $v_2$($p_T$) in Au+Au collisions at $\sqrt{s_{NN}}$ = 200 GeV from the two-particle cumulant method, the BBC event plane, and the ZDC-SMD event plane for the indicated centralities.

Anisotropy parameter $v_2$ as a function of pseudorapidity within the PHENIX central arms using event planes from the BBC and ZDC-SMD, and from the two-particle cumulant method for centrality 20–40%. The results are shown for three $p_T$ bins, which are from top to bottom: 2.0–3.0, 1.2–1.4 and 0.6–0.8 GeV/$c$.

Comparison of the $v_2${BBC} and $v_2${ZDC-SMD} obtained from the S-N and ZDC-BBC-CNT subevents as a function of pT in the 20–60% centrality range.

Charged hadron $v_2${2} for the PHENIX experiment as a function of $p_T$ in centrality 20–60%.

The PHENIX $v_2${BBC} and $v_2${ZDC-SMD} from the modified event-plane method for charged hadrons as a function of $p_T$ in different centralities.

The PHENIX $v_2${BBC} and $v_2${ZDC-SMD} from the modified event-plane method for charged hadrons as a function of $p_T$ in different centralities.

The PHENIX $v_2${BBC} and $v_2${ZDC-SMD} from the modified event-plane method for charged hadrons as a function of $p_T$ in different centralities.

The PHENIX $v_2${BBC} and $v_2${ZDC-SMD} from the modified event-plane method for charged hadrons as a function of $p_T$ in different centralities.

The PHENIX $v_2${BBC} and $v_2${ZDC-SMD} from the modified event-plane method for charged hadrons as a function of $p_T$ in different centralities.

The PHENIX $v_2${BBC} and $v_2${ZDC-SMD} from the modified event-plane method for charged hadrons as a function of $p_T$ in different centralities.

$I_{AA}$ versus $p_T^B$ for four trigger $p_T$ bins in HR+SR ($|\Delta\phi - \pi|$ < $\pi/2$) and HR ($|\Delta\phi - \pi|$ < $\pi/6$).

$I_{AA}$ versus $p_T^B$ for four trigger $p_T$ bins in HR+SR ($|\Delta\phi - \pi|$ < $\pi/2$) and HR ($|\Delta\phi - \pi|$ < $\pi/6$).

Truncated mean $<p_T^{\prime}>$ in 1 < $p_T^B$ < 5 GeV/$c$ versus $N_{part}$ for the near-side, away-side shoulder, and head regions for Au+Au and $p$+$p$ for three trigger $p_T$ bins.

Additional information containing number of events which were used to reconstruct the numvers matching to Figure 1 and 2.

Additional information containing number of events which were used to reconstruct the numvers matching to Figure 1 and 2.

Additional information containing number of events which were used to reconstruct the numvers matching to Figure 1 and 2.

The uncorrected multiplicity distributions of charged hadrons with 0.2 < $p_T$ < 2.0 GeV/$c$ for 200 GeV Au+Au collisions.

The uncorrected multiplicity distributions of charged hadrons with 0.2 < $p_T$ < 2.0 GeV/$c$ for 200 GeV Au+Au collisions.

The uncorrected multiplicity distributions of charged hadrons with 0.2 < $p_T$ < 2.0 GeV/$c$ for 200 GeV Au+Au collisions.

The uncorrected multiplicity distributions of charged hadrons with 0.2 < $p_T$ < 2.0 GeV/$c$ for 200 GeV Au+Au collisions.

The uncorrected multiplicity distributions of charged hadrons with 0.2 < $p_T$ < 2.0 GeV/$c$ for 200 GeV Au+Au collisions.

The uncorrected multiplicity distributions of charged hadrons with 0.2 < $p_T$ < 2.0 GeV/$c$ for 200 GeV Au+Au collisions.

The uncorrected multiplicity distributions of charged hadrons with 0.2 < $p_T$ < 2.0 GeV/$c$ for 200 GeV Au+Au collisions.

The uncorrected multiplicity distributions of charged hadrons with 0.2 < $p_T$ < 2.0 GeV/$c$ for 200 GeV Au+Au collisions.

The uncorrected multiplicity distributions of charged hadrons with 0.2 < $p_T$ < 2.0 GeV/$c$ for 200 GeV Au+Au collisions.

The uncorrected multiplicity distributions of charged hadrons with 0.2 < $p_T$ < 2.0 GeV/$c$ for 200 GeV Au+Au collisions.

The uncorrected multiplicity distributions of charged hadrons with 0.2 < $p_T$ < 2.0 GeV/$c$ for 200 GeV Au+Au collisions.

The uncorrected multiplicity distributions of charged hadrons with 0.2 < $p_T$ < 2.0 GeV/$c$ for 200 GeV Au+Au collisions.

The uncorrected multiplicity distributions of charged hadrons with 0.2 < $p_T$ < 2.0 GeV/$c$ for 62.4 GeV Au+Au collisions.

The uncorrected multiplicity distributions of charged hadrons with 0.2 < $p_T$ < 2.0 GeV/$c$ for 62.4 GeV Au+Au collisions.

The uncorrected multiplicity distributions of charged hadrons with 0.2 < $p_T$ < 2.0 GeV/$c$ for 62.4 GeV Au+Au collisions.

The uncorrected multiplicity distributions of charged hadrons with 0.2 < $p_T$ < 2.0 GeV/$c$ for 62.4 GeV Au+Au collisions.

The uncorrected multiplicity distributions of charged hadrons with 0.2 < $p_T$ < 2.0 GeV/$c$ for 62.4 GeV Au+Au collisions.

The uncorrected multiplicity distributions of charged hadrons with 0.2 < $p_T$ < 2.0 GeV/$c$ for 62.4 GeV Au+Au collisions.

The uncorrected multiplicity distributions of charged hadrons with 0.2 < $p_T$ < 2.0 GeV/$c$ for 62.4 GeV Au+Au collisions.

The uncorrected multiplicity distributions of charged hadrons with 0.2 < $p_T$ < 2.0 GeV/$c$ for 62.4 GeV Au+Au collisions.

The uncorrected multiplicity distributions of charged hadrons with 0.2 < $p_T$ < 2.0 GeV/$c$ for 62.4 GeV Au+Au collisions.

The uncorrected multiplicity distributions of charged hadrons with 0.2 < $p_T$ < 2.0 GeV/$c$ for 62.4 GeV Au+Au collisions.

The uncorrected multiplicity distributions of charged hadrons with 0.2 < $p_T$ < 2.0 GeV/$c$ for 200 GeV Cu+Cu collisions.

The uncorrected multiplicity distributions of charged hadrons with 0.2 < $p_T$ < 2.0 GeV/$c$ for 200 GeV Cu+Cu collisions.

The uncorrected multiplicity distributions of charged hadrons with 0.2 < $p_T$ < 2.0 GeV/$c$ for 200 GeV Cu+Cu collisions.

The uncorrected multiplicity distributions of charged hadrons with 0.2 < $p_T$ < 2.0 GeV/$c$ for 200 GeV Cu+Cu collisions.

The uncorrected multiplicity distributions of charged hadrons with 0.2 < $p_T$ < 2.0 GeV/$c$ for 200 GeV Cu+Cu collisions.

The uncorrected multiplicity distributions of charged hadrons with 0.2 < $p_T$ < 2.0 GeV/$c$ for 200 GeV Cu+Cu collisions.

The uncorrected multiplicity distributions of charged hadrons with 0.2 < $p_T$ < 2.0 GeV/$c$ for 200 GeV Cu+Cu collisions.

The uncorrected multiplicity distributions of charged hadrons with 0.2 < $p_T$ < 2.0 GeV/$c$ for 62.4 GeV Cu+Cu collisions.

The uncorrected multiplicity distributions of charged hadrons with 0.2 < $p_T$ < 2.0 GeV/$c$ for 62.4 GeV Cu+Cu collisions.

The uncorrected multiplicity distributions of charged hadrons with 0.2 < $p_T$ < 2.0 GeV/$c$ for 62.4 GeV Cu+Cu collisions.

The uncorrected multiplicity distributions of charged hadrons with 0.2 < $p_T$ < 2.0 GeV/$c$ for 62.4 GeV Cu+Cu collisions.

The uncorrected multiplicity distributions of charged hadrons with 0.2 < $p_T$ < 2.0 GeV/$c$ for 62.4 GeV Cu+Cu collisions.

The uncorrected multiplicity distributions of charged hadrons with 0.2 < $p_T$ < 2.0 GeV/$c$ for 22.5 GeV Cu+Cu collisions.

The uncorrected multiplicity distributions of charged hadrons with 0.2 < $p_T$ < 2.0 GeV/$c$ for 22.5 GeV Cu+Cu collisions.

The uncorrected multiplicity distributions of charged hadrons with 0.2 < $p_T$ < 2.0 GeV/$c$ for 200 GeV $p$+$p$ collisions.

The mean from the NBD fit as a function of $N_{part}$ for 200 GeV Au+Au collisions over the range 0.2 < $p_T$ < 2.0 GeV/$c$.

The mean from the NBD fit as a function of $N_{part}$ for 62.4 GeV Au+Au collisions over the range 0.2 < $p_T$ < 2.0 GeV/$c$.

The mean from the NBD fit as a function of $N_{part}$ for 200 GeV Cu+Cu collisions over the range 0.2 < $p_T$ < 2.0 GeV/$c$.

The mean from the NBD fit as a function of $N_{part}$ for 62.4 GeV Cu+Cu collisions over the range 0.2 < $p_T$ < 2.0 GeV/$c$.

The mean from the NBD fit as a function of $N_{part}$ for 22.5 GeV Cu+Cu collisions over the range 0.2 < $p_T$ < 2.0 GeV/$c$.

Fluctuations expressed as the scaled variance as a function of centrality for 200 GeV Au+Au collisions in the range 0.2 < $p_T$ < 2.0 GeV/$c$.

Fluctuations expressed as the scaled variance as a function of $N_{part}$ for 200 GeV Au+Au collisions for 0.2 < $p_T$ < 2.0 GeV/$c$.

Fluctuations expressed as the scaled variance as a function of $N_{part}$ for 62.4 GeV Au+Au collisions for 0.2 < $p_T$ < 2.0 GeV/$c$.

Fluctuations expressed as the scaled variance as a function of $N_{part}$ for 200 GeV Cu+Cu collisions for 0.2 < $p_T$ < 2.0 GeV/$c$.

Fluctuations expressed as the scaled variance as a function of $N_{part}$ for 62.4 GeV Cu+Cu collisions for 0.2 < $p_T$ < 2.0 GeV/$c$.

Fluctuations expressed as the scaled variance as a function of $N_{part}$ for 22.5 GeV Cu+Cu collisions for 0.2 < $p_T$ < 2.0 GeV/$c$.

Scaled variance for 200 GeV Au+Au collisions plotted as a function of $N_{part}$ for several $p_T$ ranges.

Scaled variance for 200 GeV Au+Au collisions plotted as a function of $N_{part}$ for several $p_T$ ranges.

Scaled variance for 62.4 GeV Au+Au collisions plotted as a function of $N_{part}$ for several $p_T$ ranges.

Scaled variance for 62.4 GeV Au+Au collisions plotted as a function of $N_{part}$ for several $p_T$ ranges.

Scaled variance for 62.4 GeV Au+Au collisions plotted as a function of $N_{part}$ for several $p_T$ ranges.

Scaled variance for 200 GeV Cu+Cu collisions plotted as a function of $N_{part}$ for several $p_T$ ranges.

Scaled variance for 200 GeV Cu+Cu collisions plotted as a function of $N_{part}$ for several $p_T$ ranges.

Scaled variance for 200 GeV Cu+Cu collisions plotted as a function of $N_{part}$ for several $p_T$ ranges.

Scaled variance for 62.4 GeV Cu+Cu collisions plotted as a function of $N_{part}$ for several $p_T$ ranges.

Scaled variance for 62.4 GeV Cu+Cu collisions plotted as a function of $N_{part}$ for several $p_T$ ranges.

The inverse of the parameter $k_{NBD}$ from the Negative Binomial Distribution fits for 200 GeV Au+Au collisions. The fluctuations are plotted as a function of $N_{part}$ for several $p_T$ ranges.

The inverse of the parameter $k_{NBD}$ from the Negative Binomial Distribution fits for 200 GeV Au+Au collisions. The fluctuations are plotted as a function of $N_{part}$ for several $p_T$ ranges.

The inverse of the parameter $k_{NBD}$ from the Negative Binomial Distribution fits for 62.4 GeV Au+Au collisions. The fluctuations are plotted as a function of $N_{part}$ for several $p_T$ ranges.

The inverse of the parameter $k_{NBD}$ from the Negative Binomial Distribution fits for 62.4 GeV Au+Au collisions. The fluctuations are plotted as a function of $N_{part}$ for several $p_T$ ranges.

The inverse of the parameter $k_{NBD}$ from the Negative Binomial Distribution fits for 62.4 GeV Au+Au collisions. The fluctuations are plotted as a function of $N_{part}$ for several $p_T$ ranges.

The inverse of the parameter $k_{NBD}$ from the Negative Binomial Distribution fits for 200 GeV Cu+Cu collisions. The fluctuations are plotted as a function of $N_{part}$ for several $p_T$ ranges.

The inverse of the parameter $k_{NBD}$ from the Negative Binomial Distribution fits for 200 GeV Cu+Cu collisions. The fluctuations are plotted as a function of $N_{part}$ for several $p_T$ ranges.

The inverse of the parameter $k_{NBD}$ from the Negative Binomial Distribution fits for 200 GeV Cu+Cu collisions. The fluctuations are plotted as a function of $N_{part}$ for several $p_T$ ranges.

The inverse of the parameter $k_{NBD}$ from the Negative Binomial Distribution fits for 62.4 GeV Cu+Cu collisions. The fluctuations are plotted as a function of $N_{part}$ for several $p_T$ ranges.

The inverse of the parameter $k_{NBD}$ from the Negative Binomial Distribution fits for 62.4 GeV Cu+Cu collisions. The fluctuations are plotted as a function of $N_{part}$ for several $p_T$ ranges.

The scaled variance as a funciton of $N_{part}$ for 200 GeV Au+Au collisions in the range 0.2 < $p_T$ < 2.0 GeV/$c$.

The correlation of the clan model parameters $\bar{n_c}$ and $\bar{N_c}$ for all of the collision species measured as a function of centrality.

The correlation of the clan model parameters $\bar{n_c}$ and $\bar{N_c}$ for all of the collision species measured as a function of centrality.

The correlation of the clan model parameters $\bar{n_c}$ and $\bar{N_c}$ for all of the collision species measured as a function of centrality.

The correlation of the clan model parameters $\bar{n_c}$ and $\bar{N_c}$ for all of the collision species measured as a function of centrality.

The correlation of the clan model parameters $\bar{n_c}$ and $\bar{N_c}$ for all of the collision species measured as a function of centrality.

The correlation of the clan model parameters $\bar{n_c}$ and $\bar{N_c}$ for all of the collision species measured as a function of centrality.

The correlation of the clan model parameters $\bar{n_c}$ and $\bar{N_c}$ for all of the collision species measured as a function of centrality.

The correlation of the clan model parameters $\bar{n_c}$ and $\bar{N_c}$ for all of the collision species measured as a function of centrality.

The correlation of the clan model parameters $\bar{n_c}$ and $\bar{N_c}$ for all of the collision species measured as a function of centrality.

The correlation of the clan model parameters $\bar{n_c}$ and $\bar{N_c}$ for all of the collision species measured as a function of centrality.

The correlation of the clan model parameters $\bar{n_c}$ and $\bar{N_c}$ for all of the collision species measured as a function of centrality.

$1/{N_{trig}}$ ${dN}/{d\Delta\phi}$ distributions for charge selected $\bar{p}$ and $p$ triggers both with associated $p$ for six centrality bins. Triggers have 2.5 < $p_T$ < 4.0 GeV/$c$ and associated particles have 1.8 < $p_T$ < 2.5 GeV/$c$.

$1/{N_{trig}}$ ${dN}/{d\Delta\phi}$ distributions for charge selected $\bar{p}$ and $p$ triggers both with associated $p$ for six centrality bins. Triggers have 2.5 < $p_T$ < 4.0 GeV/$c$ and associated particles have 1.8 < $p_T$ < 2.5 GeV/$c$.

$1/{N_{trig}}$ ${dN}/{d\Delta\phi}$ distributions for charge-inclusive baryon and meson triggers and associated mesons for six centrality bins. Triggers have 2.5 < $p_T$ < 4.0 GeV/$c$ and associated particles have 1.8 < $p_T$ < 2.5 GeV/$c$.

Conditional yields per trigger on the near side for charge selected and charge selected $p$ and $\bar{p}$ correlations. Triggers have 2.5 < $p_T$ < 4.0 GeV/$c$ and associated particles have 1.8 < $p_T$ < 2.5 GeV/$c$. There is an 11.4% additional normalization error on baryon associated particle points and 8.9% each on the $p$ and $\bar{p}$ associated particle points.

Conditional yields per trigger on the away side for charge selected and charge selected $p$ and $\bar{p}$ correlations. Triggers have 2.5 < $p_T$ < 4.0 GeV/$c$ and associated particles have 1.8 < $p_T$ < 2.5 GeV/$c$.

Conditional yields per trigger for baryon and meson triggers with associated mesons. Triggers have 2.5 < $p_T$ < 4.0 GeV/$c$ and associated particles have 1.8 < $p_T$ < 2.5 GeV/$c$. There is an additional 13.6% normalization error.

Fit results based on $k(\delta_{\eta})$=$1/{{2\alpha\xi}/{\delta_{\eta}}}$ ($\xi << \delta_{\eta}$).

$F_{p_T}$ (in percent) of non-random fluctuations as a function of the $p_T$ range over which $M_{p_T}$ is calculated, 0.2 GeV/$c$ < $p_T$ < $p^{max}_T$, for the 20-25% centrality class ($N_{part}$ = 181.6).

$v_2$ vs $p_T$ for identified particle species obtained in minimum bias Au+Au collisions.

$v_2$ vs $p_T$ for identified particle species obtained in minimum bias Au+Au collisions.

$v_2$ vs $p_T$ for identified particle species obtained in minimum bias Au+Au collisions.

$v_2$ vs $KE_T$ for identified particle species obtained in minimum bias Au+Au collisions.

$v_2$ vs $KE_T$ for identified particle species obtained in minimum bias Au+Au collisions.

$v_2$ vs $KE_T$ for identified particle species obtained in minimum bias Au+Au collisions.

${v_2}/{n_q}$ vs ${p_T}/{n_q}$ for identified particle species obtained in minimum bias Au+Au collisions.

${v_2}/{n_q}$ vs ${p_T}/{n_q}$ for identified particle species obtained in minimum bias Au+Au collisions.

${v_2}/{n_q}$ vs ${p_T}/{n_q}$ for identified particle species obtained in minimum bias Au+Au collisions.

${v_2}/{n_q}$ vs ${KE_T}/{n_q}$ for identified particle species obtained in minimum bias Au+Au collisions.

${v_2}/{n_q}$ vs ${KE_T}/{n_q}$ for identified particle species obtained in minimum bias Au+Au collisions.

${v_2}/{n_q}$ vs ${KE_T}/{n_q}$ for identified particle species obtained in minimum bias Au+Au collisions.

$v(Q)$, for central events, as a function of the azimuthal coverage of the detector for data. (Angles in degrees.)

The effect on $v(Q)$ of varying the acceptance in $\varphi_r$. The solid curve shows the expectation from global charge conservation. The 10% most central events are used for data and RQMD. (Angles in degrees.)