The dielectron $v_2$ as a function of $p_T$ in minimum-bias Au + Au collisions at $\sqrt{s_{NN}}$ = 200 GeV for the $\eta$ mass region.

The dielectron $v_2$ as a function of $p_T$ in minimum-bias Au + Au collisions at $\sqrt{s_{NN}}$ = 200 GeV for the charm + $ ho^0$ mass region.

The dielectron $v_2$ as a function of $p_T$ in minimum-bias Au + Au collisions at $\sqrt{s_{NN}}$ = 200 GeV for the $\omega$ mass region.

The dielectron $v_2$ as a function of $p_T$ in minimum-bias Au + Au collisions at $\sqrt{s_{NN}}$ = 200 GeV for the $\phi$ mass region.

The dielectron $v_2$ as a function of $p_T$ in minimum-bias Au + Au collisions at $\sqrt{s_{NN}}$ = 200 GeV for the charm + thermal radiation mass region.

The dielectron $v_2$ as a function of $p_T$ in minimum-bias Au + Au collisions at $\sqrt{s_{NN}}$ = 200 GeV expected from decays of a mesonic cocktail.

The dielectron $v_2$ as a function of $p_T$ in minimum-bias Au + Au collisions at $\sqrt{s_{NN}}$ = 200 GeV expected from decays of a mesonic cocktail.

The dielectron $v_2$ as a function of $p_T$ in minimum-bias Au + Au collisions at $\sqrt{s_{NN}}$ = 200 GeV expected from decays of a mesonic cocktail.

The dielectron $v_2$ as a function of $p_T$ in minimum-bias Au + Au collisions at $\sqrt{s_{NN}}$ = 200 GeV expected from decays of a mesonic cocktail.

The dielectron $v_2$ as a function of $p_T$ in minimum-bias Au + Au collisions at $\sqrt{s_{NN}}$ = 200 GeV expected from decays of a mesonic cocktail.

The dielectron $v_2$ as a function of $p_T$ in minimum-bias Au + Au collisions at $\sqrt{s_{NN}}$ = 200 GeV expected from the $car{c}$ contribution.

The $p_T$-integrated dielectron $v_2$ as a function of $M_{ee}$ in minimum-bias Au + Au collisions at $\sqrt{s_{NN}}$ = 200 GeV.

The $v_2$ of hadrons $\pi$, $K$, $p$, $\phi$ and $\Lambda$ in minimum-bias Au + Au collisions at $\sqrt{s_{NN}}$ = 200 GeV for the $\pi^0$ mass region.

The $p_T$-integrated dielectron $v_2$ as a function of $M_{ee}$ in minimum-bias Au + Au collisions at $\sqrt{s_{NN}}$ = 200 GeV, expected from the decays of $\pi^0$, $\eta$, $\omega$ and $\phi$ and from the $car{c}$ contribution.

Per-trigger yield versus $\Delta\phi$ for various trigger and partner $p_T$ ($p^a_T \otimes p^b_T$), arranged by increasing pair proxy energy (sum of $p^a_T$ and $p^b_T$), in p + p collisions for 5-10 $\otimes$ 2-3, 4-5 $\otimes$ 4-5, 5-10 $\otimes$ 3-5, and 5-10 $\otimes$ 5-10 GeV/c.

Per-trigger yield versus $\Delta\phi$ for various trigger and partner $p_T$ ($p^a_T \otimes p^b_T$), arranged by increasing pair proxy energy (sum of $p^a_T$ and $p^b_T$), in p + p collisions for 3-4 $\otimes$ 0.4-1, 3-4 $\otimes$ 1-2, 3-4 $\otimes$ 2-3, and 3-4 $\otimes$ 3-4 GeV/c.

Per-trigger yield versus $\Delta\phi$ for various trigger and partner $p_T$ ($p^a_T \otimes p^b_T$), arranged by increasing pair proxy energy (sum of $p^a_T$ and $p^b_T$), in Au + Au collisions for 3-4 $\otimes$ 0.4-1, 3-4 $\otimes$ 1-2, 3-4 $\otimes$ 2-3, and 3-4 $\otimes$ 3-4 GeV/c.

Per-trigger yield versus $\Delta\phi$ for various trigger and partner $p_T$ ($p^a_T \otimes p^b_T$), arranged by increasing pair proxy energy (sum of $p^a_T$ and $p^b_T$), in Au + Au collisions for 5-10 $\otimes$ 2-3, 4-5 $\otimes$ 4-5, 5-10 $\otimes$ 3-5, and 5-10 $\otimes$ 5-10 GeV/c.

RHS versus $p^b_T$ for p + p collisions for four trigger selections.

RHS versus $p^b_T$ for Au + Au collisions for four trigger selections.

RHS versus $N_{part}$ for 2-3 $\otimes$ 2-3 GeV/$c$. Shaded bars (brackets) represent $p_T$ -correlated uncertainties due to elliptic flow (ZYAM procedure). The most left point is from p+p.

Per-trigger yield $\Delta\phi$ distribution and corresponding fits for 2−3 $\otimes$ 2−3 GeV/c in 0-5% Au+Au collisions.

a) D versus $N_{part}$ from FIT 1 and FIT2 for 2 - 3 GeV/$c$ bin. The error bars are the statistical errors; The shaded bars and brackets are the systematic errors due to $v_2$. b) The fraction of the shoulder Gaussian yield relative to the total away-side yield as function of $N_{part}$ determined from FIT2.

Values of D determined from FIT1 as function of partner $p_T$ for 2-3 GeV/c $p_T$ range in 0-20% Au+Au.

Values of D determined from FIT1 as function of partner $p_T$ for 3-4 GeV/c $p_T$ range in 0-20% Au+Au.

Values of D determined from FIT1 as function of partner $p_T$ for 4-5 GeV/c $p_T$ range in 0-20% Au+Au.

Values of D determined from FIT2 as function of partner $p_T$ for 2-3 GeV/c $p_T$ range in 0-20% Au+Au.

Values of D determined from FIT2 as function of partner $p_T$ for 3-4 GeV/c $p_T$ range in 0-20% Au+Au.

Values of D determined from FIT2 as function of partner $p_T$ for 4-5 GeV/c $p_T$ range in 0-20% Au+Au.

Per-trigger yield as function of partner $p_T$ for the HR region for the 0-20% Au + Au collisions. Results for four trigger pT in 2-3, 3-4, 4-5, 5-10 GeV/c are shown.

Per-trigger yield as function of partner $p_T$ for the HR region for the 20-40% Au + Au collisions. Results for four trigger pT in 2-3, 3-4, 4-5, 5-10 GeV/c are shown.

Per-trigger yield as function of partner $p_T$ for the HR region for the 40-60% Au + Au collisions. Results for four trigger pT in 2-3, 3-4, 4-5, 5-10 GeV/c are shown.

Per-trigger yield as function of partner $p_T$ for the HR region for the 60-92% Au + Au collisions. Results for four trigger pT in 2-3, 3-4, 4-5, 5-10 GeV/c are shown.

Per-trigger yield as function of partner $p_T$ for the HR region for the p + p collisions. Results for four trigger pT in 2-3, 3-4, 4-5, 5-10 GeV/c are shown.

Per-trigger yield as function of partner $p_T$ for the SR region for the 0-20% Au + Au collisions. Results for four trigger pT in 2-3, 3-4, 4-5, 5-10 GeV/c are shown.

Per-trigger yield as function of partner $p_T$ for the SR region for the 20-40% Au + Au collisions. Results for four trigger pT in 2-3, 3-4, 4-5, 5-10 GeV/c are shown.

Per-trigger yield as function of partner $p_T$ for the SR region for the 40-60% Au + Au collisions. Results for four trigger pT in 2-3, 3-4, 4-5, 5-10 GeV/c are shown.

Per-trigger yield as function of partner $p_T$ for the SR region for the 60-92% Au + Au collisions. Results for four trigger pT in 2-3, 3-4, 4-5, 5-10 GeV/c are shown.

Per-trigger yield as function of partner $p_T$ for the SR region for the p + p collisions. Results for four trigger pT in 2-3, 3-4, 4-5, 5-10 GeV/c are shown.

Near-side Gaussian widths versus partner $p_T$ for four trigger $p_T$ ranges for 0-20 % Au + Au.

Near-side Gaussian widths versus partner $p_T$ for four trigger $p_T$ ranges for 0-20 % p + p.

Near-side Gaussian widths versus $N_{part}$ for five successively increasing $p^a_T$ $\otimes$ $p^b_T$. The most left point in each panel represents the value from p + p.

Near-side yield in $\lvert \Delta\phi \rvert$ < $\pi$/3 versus partner $p_T$ for four trigger $p_T$ selections for 0-20% p + p.

Near-side yield in $\lvert \Delta\phi \rvert$ < $\pi$/3 versus partner $p_T$ for four trigger $p_T$ selections for 0-20% Au + Au.

Per-trigger yield versus $\Delta\eta$ for 0-20% central Au + Au collisions. Results are shown for four $p^a_T \otimes p^b_T$ selections as indicated.

Per-trigger yield versus $\Delta\eta$ for 0-20% central p + p collisions. Results are shown for four $p^a_T \otimes p^b_T$ selections as indicated.

The projection of 2-D per-trigger yield in $\lvert \Delta\eta \rvert$ < 0.7 $\otimes$ $\lvert \Delta\phi \rvert$ < 0.7 and 0-20% central Au+Au collisions onto $\Delta\eta$ for four $p^a_T \otimes p^b_T$ selections.

The projection of 2-D per-trigger yield in $\lvert \Delta\eta \rvert$ < 0.7 $\otimes$ $\lvert \Delta\phi \rvert$ < 0.7 and 0-20% central Au+Au collisions onto $\Delta\phi$ for four $p^a_T \otimes p^b_T$ selections.

Truncated mean $p_T$, $\langle p_T \rangle$, in 1 <pbT < 5 GeV/c versus $N_{part}$ for the head region for Au + Au and p + p for four trigger $p_T$ bins.

Truncated mean $p_T$, $\langle p_T \rangle$, in 1 <pbT < 5 GeV/c versus $N_{part}$ for the away-side shoulder region for Au + Au and p + p for four trigger $p_T$ bins.

Truncated mean $p_T$, $\langle p_T \rangle$, in 1 <pbT < 5 GeV/c versus $N_{part}$ for the near-side region for Au + Au and p + p for four trigger $p_T$ bins.

Truncated mean $p_T$, $\langle p_T \rangle$, calculated for five partner pT ranges in the HR. Results for triggers in 3 − 4 GeV/c.

Truncated mean $p_T$, $\langle p_T \rangle$, calculated for five partner pT ranges in the HR. Results for triggers in 4 − 5 GeV/c.

The near-side pair suppression factor $J_{AA}$ in 0-20% Au + Au collisions as function of $p^b_T$ for various ranges of $p^a_T$.

The near-side $J_{AA}$ as function of $p^b_T$ and $p^{sum}_T$ for four $p^a_T$ bins and three centrality bins. Data for 20-40%.

The near-side $J_{AA}$ as function of $p^b_T$ and $p^{sum}_T$ for four $p^a_T$ bins and three centrality bins. Data for 40-60%.

The near-side $J_{AA}$ as function of $p^b_T$ and $p^{sum}_T$ for four $p^a_T$ bins and three centrality bins. Data for 60-92%.

The pair suppression factor $J_{AA}$ for the away-side HR in 0-20% Au + Au collisions as function of $p_b^T$ for various ranges of $p_a^T$.

The pair suppression factor $J_{AA}$ for the away-side HR in 20-40% Au + Au collisions as function of $p_b^T$ for various ranges of $p_a^T$.

The pair suppression factor $J_{AA}$ for the away-side HR in 40-60% Au + Au collisions as function of $p_b^T$ for various ranges of $p_a^T$.

The pair suppression factor $J_{AA}$ for the away-side HR in 60-92% Au + Au collisions as function of $p_b^T$ for various ranges of $p_a^T$.

Per-trigger yield versus $\Delta\phi$ for successively increasing trigger and partner $p_T$ ($p^a_T \otimes p^b_T$) in 0-20 % Au + Au collisions.

Per-trigger yield versus $\Delta\phi$ for successively increasing trigger and partner $p_T$ ($p^a_T \otimes p^b_T$) in p + p collisions.

Per-trigger yield versus $\Delta\phi$ for successively increasing trigger and partner $p_T$ ($p^a_T \otimes p^b_T$) in 20-40 % Au + Au collisions.

Per-trigger yield versus $\Delta\phi$ for successively increasing trigger and partner $p_T$ ($p^a_T \otimes p^b_T$) in p + p collisions.

Per-trigger yield versus $\Delta\phi$ for successively increasing trigger and partner $p_T$ ($p^a_T \otimes p^b_T$) in 60-92 % Au + Au collisions.

Per-trigger yield versus $\Delta\phi$ for successively increasing trigger and partner $p_T$ ($p^a_T \otimes p^b_T$) in p + p collisions.

$pi^0$ nuclear modification factors vs reaction plane, $\Delta \phi = 0^{\circ}-15^{\circ}$

$pi^0$ nuclear modification factors vs reaction plane, $\Delta \phi = 15^{\circ}-30^{\circ}$

$pi^0$ nuclear modification factors vs reaction plane, $\Delta \phi = 30^{\circ}-45^{\circ}$

$pi^0$ nuclear modification factors vs reaction plane, $\Delta \phi = 45^{\circ}-60^{\circ}$

$pi^0$ nuclear modification factors vs reaction plane, $\Delta \phi = 60^{\circ}-75^{\circ}$

$pi^0$ nuclear modification factors vs reaction plane, $\Delta \phi = 75^{\circ}-90^{\circ}$

Centrality dependence of $p_T$ ingetrated $R_{AA}(\Delta \phi)$

Total $J/\psi$ efficiency as a function of cos$\theta$ for $2 < p_{T}^{J/\psi} < 3$ GeV/c

Total $J/\psi$ efficiency as a function of cos$\theta$ for $3 < p_{T}^{J/\psi} < 4$ GeV/c

Total $J/\psi$ efficiency as a function of cos$\theta$ for $4 < p_{T}^{J/\psi} < 6$ GeV/c

Different contributions to total $J/\psi$ efficiency as a function of cos$\theta$ for $2 < p_{T}^{J/\psi} < 3$ GeV/c

Different contributions to total $J/\psi$ efficiency as a function of cos$\theta$ for $3 < p_{T}^{J/\psi} < 4$ GeV/c

Different contributions to total $J/\psi$ efficiency as a function of cos$\theta$ for $4 < p_{T}^{J/\psi} < 6$ GeV/c

Corrected cos$\theta$ distribution for $2 < p_{T}^{J/\psi} < 3$ GeV/c

Corrected cos$\theta$ distribution for $3 < p_{T}^{J/\psi} < 4$ GeV/c

Corrected cos$\theta$ distribution for $4 < p_{T}^{J/\psi} < 6$ GeV/c

$J/\psi$ polarization parameter $\lambda_{\theta}$ as a function of $J/\psi$ $p_{T}$ for rapidity |y|<1

The simulated electron pair invariant mass distributions for electrons at $8 < p_{T} < 10$ GeV/c

The distribution of the minimum distance $\Delta$Z (cm) between an electron track projection point at the Barrel Electromagnetic Calorimeter (BEMC) and all BEMC clusters along the beam direction from unlike-sign electron candidate pairs, like-sign electron candidate pairs and unlike-minus-like.

The momentum over energy $p/E_0$ distribution ($p$ - momentum, $E_0$ - the energy of the most energetic tower in a Barrel Electromagnetic Calorimeter cluster) from unlike-sign electron candidate pairs, like-sign electron candidate pairs and unlike-minus-like.

Electron identification efficiency of the cuts on number of TPC points, n$\sigma_e$ and Barrel Electromagnetic Calorimeter cuts in the Run2008 analysis. The total efficiency is the product of all individual ones.

Electron identification efficiency of the cuts on number of TPC points, n$\sigma_e$ and Barrel Electromagnetic Calorimeter cuts in the Run2005 analysis. The total efficiency is the product of all individual ones.

n$\sigma_e$ distribution in the Run2008 analysis for unlike-sign, like-sign and unlike-minus-like pairs at $2.5 < p_T < 3.0$ GeV after applying all the electron identification cuts except the n$\sigma_e$ cut.

n$\sigma_e$ distribution in the Run2008 analysis for unlike-sign, like-sign and unlike-minus-like pairs at $8 < p_T < 10$ GeV after applying all the electron identification cuts except the n$\sigma_e$ cut.

The mean and width of the Gaussian fitting functions for pure photonic electron (unlike - minus-like) n$\sigma_e$ distribution as shown in Fig. 4(a) and Fig 4(b) for each $p_T$ bin.

n$\sigma_e$ distribution for inclusive electrons at (a) $2.5 < p_T < 3.0$ GeV/c in the Run2008 analysis, (b) $2.5< p_T < 3.5$ GeV/c in the Run2005 analysis and (c) $8.0 < p_T < 10.0$ GeV/c in the Run2008 analysis after applying all electron identification cuts except the n$\sigma_e$ cut.

Purity of the inclusive electron sample as a function of $p_T$ in data from Run 2008.

Purity of the inclusive electron sample as a function of $p_T$ in data from Run 2005 with High Tower HT1 trigger.

Purity of the inclusive electron sample as a function of $p_T$ in data from Run 2005 with High Tower HT2 trigger.

Derived $p_T$ spectrum for inclusive photons and the uncertainty represented by the region between the spectra of $\pi^0$ and $\eta$ decay photons and inclusive photon with doubled direct photon yield.

Ratio of photon $p_T$ spectra from Fig. 7(a) to inclusive photons $p_T$ spectrum.

Photonic electron reconstruction efficiency as a function of $p_T$ for $\gamma$ conversion, $\pi^0$ and $\eta$ Dalitz decay for the Run2008 analysis, and for their combination for the Run2008 analysis.

Photonic electron reconstruction efficiency as a function of $p_T$ for $\gamma$ conversion for the Run2005 analysis.

The adc0 distribution for high-tower trigger events. The adc0 is the offline ADC value of a BEMC cluster’s most energetic tower which is one of the high-towers responsible for firing a high-tower trigger.

The $\chi^2$ as a function of the adc0 cut.

Raw inclusive electron $p_T$ spectrum from VPD trigger in Run2008 d+Au collisions before (open squares) and after applying the adc0 > 193 cut.

The adc0 distribution for data at $p_T$ = 4.0 − 5.0 GeV/c

The adc0 distribution for simulations at $p_T$ = 4.0 − 5.0 GeV/c

The $p_T$ dependence of high-tower trigger efficiency from VPD data for Run2008 analysis.

The $p_T$ dependence of high-tower trigger efficiency from simulation for Run2008 analysis.

The $p_T$ dependence of high-tower trigger efficiency for combined VPD data and simulation for Run2008 analysis.

Raw inclusive electron $p_T$ spectrum for minimum-bias (MB) and two high-tower triggers, i.e. HT1 and HT2 for Run2005 analysis.

Raw inclusive electron $p_T$ spectrum for high-tower trigge HT1 for Run2005 analysis.

Raw inclusive electron $p_T$ spectrum for high-tower trigge HT2 for Run2005 analysis.

The $p_T$ dependence of trigger efficiency for the high-tower trigger HT1 data for Run2005.

The $p_T$ dependence of trigger efficiency for the high-tower trigger HT2 data for Run2005.

Distribution of the calculated BBC cross section $\sigma_{BBC}$ before and after removing events at the beginning and the end of Run2008.

Invariant cross section of the electron from decays of $J/\psi$.

Invariant cross section of the electron from decays of $\varUpsilon$.

Invariant cross section of the electron from Drell-Yan decays.

Invariant cross section of the electron from decays of light vector mesons ($\rho,\omega, \phi$).

The uncertainty of the $J/\psi$ feeddown.

Ratio of non-photonic to photonic electron yield from the Run2008 analysis.

Ratio of non-photonic to photonic electron yield from the Run2005 analysis.

Invariant cross section for non-photonic electron production $(e^+ + e^-)/2$ in p+p collisions from the Run2008 analysis.

Invariant cross section for non-photonic electron production $(e^+ + e^-)/2$ in p+p collisions from the Run2005 analysis.

Invariant cross section for non-photonic electron production $(e^+ + e^-)/2$ in p+p collisions from the combined Run2005 and 2008 analyses.

Invariant cross section for non-photonic electron production $(e^+ + e^-)/2$ in p+p collisions (update on PRL 98 (2007) 192301).

Ratio of invariant cross section for non-photonic electron production $(e^+ + e^-)/2$ in p+p collisions from the combined Run2005 and 2008 analyses over FONLL.

Ratio of invariant cross section for non-photonic electron production $(e^+ + e^-)/2$ in p+p collisions (update on PRL 98 (2007) 192301) over FONLL.

Invariant cross section for electrons $(e^+ + e^-)/2$ from bottom meson decay, together with the ratio of the measurements to the FONLL predictions.

Invariant cross section for electrons $(e^+ + e^-)/2$ from charm meson decay, together with the ratio of the measurements to the FONLL predictions.

Amplitude (A) as a function of centrality for Like Sign (LS) and Charge Independent (CI) particles for Au+Au collision at sqrt(snn)=200 GEV.

The third harmonic coefficient as a function of pseudorapidity for 0-5% centrality Au+Au collisions at sqrt(snn)=200 GEV.

The third harmonic coefficient as a function of pseudorapidity for 5-10% centrality Au+Au collisions at sqrt(snn)=200 GEV.

The third harmonic coefficient as a function of pseudorapidity for 10-20% centrality Au+Au collisions at sqrt(snn)=200 GEV.

The third harmonic coefficient as a function of pseudorapidity for 20-30% centrality Au+Au collisions at sqrt(snn)=200 GEV.

The third harmonic coefficient as a function of pseudorapidity for 30-40% centrality Au+Au collisions at sqrt(snn)=200 GEV.

The third harmonic coefficient as a function of pseudorapidity for 40-60% centrality Au+Au collisions at sqrt(snn)=200 GEV.

The third harmonic coefficient as a function of transverse momentum for the wideGaussian method and for fhe event plane in the TPC for 0-5% centrality Au+Au collisions at sqrt(snn)=200 GEV.

The third harmonic coefficient as a function of transverse momentum for the wideGaussian method and for fhe event plane in the TPC for 30-40% centrality Au+Au collisions at sqrt(snn)=200 GEV.

The third harmonic coefficient as a function of transverse momentum for 0-10% centrality Au+Au collisions at sqrt(snn)=200 GEV.

The third harmonic coefficient as a function of transverse momentum for 10-20% centrality Au+Au collisions at sqrt(snn)=200 GEV.

The third harmonic coefficient as a function of transverse momentum for 20-30% centrality Au+Au collisions at sqrt(snn)=200 GEV.

The third harmonic coefficient as a function of transverse momentum for 30-40% centrality Au+Au collisions at sqrt(snn)=200 GEV.

The third harmonic coefficient as a function of transverse momentum for 40-50% centrality Au+Au collisions at sqrt(snn)=200 GEV.

The third harmonic coefficient as a function of transverse momentum for 50-60% centrality Au+Au collisions at sqrt(snn)=200 GEV.

The third harmonic coefficient as a function of centrality from different methods of analysis for Au+Au collisions at sqrt(snn)=200 GEV.

The square of the third harmonic coefficient as a function of pseudorapidity separation (Deta) for Au+Au collisions at sqrt(snn)=200 GEV.

The fourth power of the third harmonic coefficient from four-particle cumulants is plotted as a function of centrality for Au+Au collisions at sqrt(snn)=200 GEV.

The ratio of the fourth power of the third harmonic coefficient from four-particle cumulants and the fourth power of the third harmonic flow is plotted as a function of centrality for Au+Au collisions at sqrt(snn)=200 GEV.

The v_2{TPC} as a function of transverse momentum for 0-5% centrality Au+Au collisions at sqrt(snn)=200 GEV.

The v_2{TPC} as a function of transverse momentum for 20-30% centrality Au+Au collisions at sqrt(snn)=200 GEV.

The v_2{TPC} as a function of transverse momentum for 30-40% centrality Au+Au collisions at sqrt(snn)=200 GEV.

The v_3{TPC} as a function of transverse momentum for 0-5% centrality Au+Au collisions at sqrt(snn)=200 GEV.

The v_3{TPC} as a function of transverse momentum for 20-30% centrality Au+Au collisions at sqrt(snn)=200 GEV.

The v_3{TPC} as a function of transverse momentum for 30-40% centrality Au+Au collisions at sqrt(snn)=200 GEV.

Dphi correlation functions for 2 < pT < 4 GEV/c and 4 < p_T^trig < 6 GEV/c.

Deta correlation functions for 0.15 < pT < 4 GEV/c and 4 < p_T^trig < 6 GEV/c.

Deta correlation functions for 0.15<pT<4 GEV/c and 4 < p_T^trig < 6 GEV/c.

Deta awayside correlation functions for 0.15 < pT < 4 GEV/c and 4 < p_T^trig < 6 GEV/c.

Deta correlation functions for 2<pT<4 GEV/c and 4 < p_T^trig < 6 GEV/c.

Dphi correlation functions for 2<pT<4 GEV/c and 4 < p_T^trig < 6 GEV/c.

Associated multiplicity vs centrality for 4 < p_T^trig < 6 GEV/c.

Associated multiplicity vs centrality for 6 < p_T^trig < 10 GEV/c.

Associated total p_T vs centrality for 4 < p_T^trig < 6 GEV/c.

Associated total p_T vs centrality for 6 < p_T^trig < 10 GEV/c.

Near-side associated charged hadron p_T spectra for 4 < p_T^trig < 6 GEV/c.

Near-side associated charged hadron p_T spectra for 4 < p_T^trig < 6 GEV/c.

Near-side associated charged hadron p_T spectra for 4 < p_T^trig < 6 GEV/c.

Away-side associated charged hadron p_T spectra for 4 < p_T^trig < 6 GEV/c.

Away-side associated charged hadron p_T spectra for 4 < p_T^trig < 6 GEV/c.

Away-side associated charged hadron p_T spectra for 4 < p_T^trig < 6 GEV/c.

Near-side AA/pp ratio for 4 < p_T^trig < 6 GEV/c.

Near-side AA/pp ratio for 4 < p_T^trig < 6 GEV/c.

Away-side AA/pp ratio for 4 < p_T^trig < 6 GEV/c.

Away-side AA/pp ratio for 4 < p_T^trig < 6 GEV/c.

Away-side associated <p_T> vs centrality for 4 < p_T^trig < 6 GEV/c.

Away-side associated <p_T> vs centrality for 6 < p_T^trig < 10 GEV/c.

$N_{J/psi}/<N_{J/\psi}> / N_{ch}/<N_{ch}>$, $|y^{J/\psi}|<0.9$, $p_{T}^{J/\psi}$ 8-15 and 15-40 GeV/c, event selection SPD

$N_{J/psi}/<N_{J/\psi}> / N_{ch}/<N_{ch}>$, $|y^{J/\psi}|<0.9$, $p_{T}^{J/\psi}$ 0-4 and 4-8 GeV/c, event selection V0

$N_{J/psi}/<N_{J/\psi}> / N_{ch}/<N_{ch}>$, $|y^{J/\psi}|<0.9$, $p_{T}^{J/\psi}$ 8-15 and 15-40 GeV/c, event selection V0

Multiplicity in Au+Au collisions at $\sqrt{s_{NN}}$ = 130 GeV

Transverse energy in Au+Au collisions at $\sqrt{s_{NN}}$ = 62.4 GeV

Multiplicity in Au+Au collisions at $\sqrt{s_{NN}}$ = 62.4 GeV

Transverse energy in Au+Au collisions at $\sqrt{s_{NN}}$ = 39 GeV

Multiplicity in Au+Au collisions at $\sqrt{s_{NN}}$ = 39 GeV

Transverse energy in Au+Au collisions at $\sqrt{s_{NN}}$ = 27 GeV

Multiplicity in Au+Au collisions at $\sqrt{s_{NN}}$ = 27 GeV

Transverse energy in Au+Au collisions at $\sqrt{s_{NN}}$ = 19.6 GeV

Multiplicity in Au+Au collisions at $\sqrt{s_{NN}}$ = 19.6 GeV

Transverse energy in Au+Au collisions at $\sqrt{s_{NN}}$ = 14.5 GeV

Multiplicity in Au+Au collisions at $\sqrt{s_{NN}}$ = 14.5 GeV

Transverse energy in Au+Au collisions at $\sqrt{s_{NN}}$ = 7.7 GeV

Multiplicity in Au+Au collisions at $\sqrt{s_{NN}}$ = 7.7 GeV

Transverse energy in Cu+Cu collisions at $\sqrt{s_{NN}}$ = 200 GeV

Multiplicity in Cu+Cu collisions at $\sqrt{s_{NN}}$ = 200 GeV

Transverse energy in Cu+Cu collisions at $\sqrt{s_{NN}}$ = 62.4 GeV

Multiplicity in Cu+Cu collisions at $\sqrt{s_{NN}}$ = 62.4 GeV

Transverse energy in Cu+Au collisions at $\sqrt{s_{NN}}$ = 200 GeV

Multiplicity in Cu+Au collisions at $\sqrt{s_{NN}}$ = 200 GeV

Transverse energy in U+U collisions at $\sqrt{s_{NN}}$ = 193 GeV

Multiplicity in U+U collisions at $\sqrt{s_{NN}}$ = 193 GeV

Transverse energy in d+Au collisions at $\sqrt{s_{NN}}$ = 200 GeV

Multiplicity in d+Au collisions at $\sqrt{s_{NN}}$ = 200 GeV

Transverse energy in He+Au collisions at $\sqrt{s_{NN}}$ = 200 GeV

Multiplicity in He+Au collisions at $\sqrt{s_{NN}}$ = 200 GeV

pion dAu Invariant yields versus $p_T$

proton AuAu Invariant yields versus $p_T$

proton dAu Invariant yields versus $p_T$

AuAu Ratio kaon versus $p_T$

dAu Ratio kaon versus $p_T$

AuAu Ratio pion versus $p_T$

dAu Ratio pion versus $p_T$

AuAu Ratio proton versus $p_T$

dAu Ratio proton versus $p_T$

AuAu Ratio $kaon/\pi$ versus $p_T$

dAu Ratio $kaon/\pi$ versus $p_T$

AuAu Ratio $p/\pi$ versus $p_T$

dAu Ratio $p/\pi$ versus $p_T$

AuAu kaon $R_{CP}$ versus $p_T$

AuAu pion $R_{CP}$ versus $p_T$

AuAu proton $R_{CP}$ versus $p_T$

AuAu kaon $R_{AA}$ versus $p_T$ (Uncertainity type C is not included)

AuAu pion $R_{AA}$ versus $p_T$ (Uncertainity type C is not included)

AuAu proton $R_{AA}$ versus $p_T$ (Uncertainity type C is not included)

dAu kaon $R_{dA}$ versus $p_T$ (Uncertainity type C is not included)

dAu pion $R_{dA}$ versus $p_T$ (Uncertainity type C is not included)

dAu proton $R_{dA}$ versus $p_T$ (Uncertainity type C is not included)

kaon AuAu Ratio of Invariant yields versus $p_T$

pion AuAu Ratio of Invariant yields versus $p_T$

proton AuAu Ratio of Invariant yields versus $p_T$

$\pi^0$ meson $v_4, 20-40\%$ centrality

$\pi^0$ meson $v_4, 40-60\%$ centrality

$\pi^0$ meson $v_4, 20-60\%$ centrality

$\pi^0$ meson $v_4/v_2^2, 0-20\%$ centrality

$\pi^0$ meson $v_4/v_2^2, 20-40\%$ centrality

$\pi^0$ meson $v_4/v_2^2, 40-60\%$ centrality

$\pi^0$ meson $v_4/v_2^2, 20-60\%$ centrality

INVARIANT YIELDS

RAA

RAA

RAA

RAA

pT AVERAGE RAA

Sloss

Sloss

Sloss

n_eff(xT)

n_eff(xT)

n_eff(xT)

n_eff(xT)

$p^{a}_{T} = 4-5$ GeV/$c$

$p^{a}_{T} = 5-7$ GeV/$c$

$p^{a}_{T} = 3-4$ GeV/$c$

$p^{a}_{T} = 4-5$ GeV/$c$

$p^{a}_{T} = 5-7$ GeV/$c$

$I^{out}_{AA}/I^{in}_{AA}$ ratio, central collisions

$I^{out}_{AA}/I^{in}_{AA}$ ratio, central collisions

$I^{out}_{AA}/I^{in}_{AA}$ ratio, mid-central collisions

$I^{out}_{AA}/I^{in}_{AA}$ ratio, mid-central collisions

Per-Trigger Azimuthal Yields, central collisions, 4-7 x 3-4 GeV/c

Per-Trigger Azimuthal Yields, central collisions, 4-7 x 4-5 GeV/c

Per-Trigger Azimuthal Yields, central collisions, 4-7 x 5-7 GeV/c

Per-Trigger Azimuthal Yields, mid-central collisions, 4-7 x 3-4 GeV/c

Per-Trigger Azimuthal Yields, mid-central collisions, 4-7 x 4-5 GeV/c

Per-Trigger Azimuthal Yields, mid-central collisions, 4-7 x 5-7 GeV/c

Per-Trigger Azimuthal Yields, central collisions, 4-7 x 3-4 GeV/c

$p^{a}_{T} = 3-4$ GeV/$c$

$p^{a}_{T} = 4-5$ GeV/$c$

$p^{a}_{T} = 5-7$ GeV/$c$

Inelastic cross section measured in d+Au at $\sqrt{s}$=200 GeV through $\eta \rightarrow \pi^{0} \pi^{+} \pi^{-}$

Invariant yields measured in d+Au for different centrality classes

Invariant yields measured in Au+Au for different centrality classes

Nuclear modification factor $R_{dA}$ for $\eta$ measured in d+Au for different centrality classes

Nuclear modification factor $R_{AA}(p_T)$ for $\eta$ measured in Au+Au for different centrality classes

Ratio of $\eta$ and $\pi^{0}$ measured in p+p at $\sqrt{s}$=200 GeV

Ratio of $\eta$ and $\pi^{0}$ measured in d+Au for different centrality classes

Ratio of $\eta$ and $\pi^{0}$ measured in Au+Au for different centrality classes

Midrapidity (|y| < 0.1) transverse momentum spectra for (d) K− in Au+Au collisions at √sNN = 7.7 GeV for different centralities. The spectra for centralities other than 0–5% are scaled for clarity as shown in the figure. The curves represent the Bose-Einstein, mT -exponential, and double-exponential function fits to 0–5% central data for pions, kaons, and (anti)protons, respectively. The uncertainties are statistical and systematic added in quadrature.

Midrapidity (|y| < 0.1) transverse momentum spectra for (c) K+ in Au+Au collisions at √sNN = 7.7 GeV for different centralities. The spectra for centralities other than 0–5% are scaled for clarity as shown in the figure. The curves represent the Bose-Einstein, mT -exponential, and double-exponential function fits to 0–5% central data for pions, kaons, and (anti)protons, respectively. The uncertainties are statistical and systematic added in quadrature.

Midrapidity (|y| < 0.1) transverse momentum spectra for (f) p¯ in Au+Au collisions at √sNN = 7.7 GeV for different centralities. The spectra for centralities other than 0–5% are scaled for clarity as shown in the figure. The curves represent the Bose-Einstein, mT -exponential, and double-exponential function fits to 0–5% central data for pions, kaons, and (anti)protons, respectively. The uncertainties are statistical and systematic added in quadrature.

Midrapidity (|y| < 0.1) transverse momentum spectra for (e) p in Au+Au collisions at √sNN = 7.7 GeV for different centralities. The spectra for centralities other than 0–5% are scaled for clarity as shown in the figure. The curves represent the Bose-Einstein, mT -exponential, and double-exponential function fits to 0–5% central data for pions, kaons, and (anti)protons, respectively. The uncertainties are statistical and systematic added in quadrature.

Midrapidity (|y| < 0.1) transverse momentum spectra for (b) π− in Au+Au collisions at √sNN = 11.5 GeV for different centralities. The spectra for centralities other than 0–5% are scaled for clarity as shown in the figure. The curves represent the Bose-Einstein, mT -exponential, and double-exponential function fits to 0–5% central data for pions, kaons, and (anti)protons, respectively. The uncertainties are statistical and systematic added in quadrature.

Midrapidity (|y| < 0.1) transverse momentum spectra for (a) π+ in Au+Au collisions at √sNN = 11.5 GeV for different centralities. The spectra for centralities other than 0–5% are scaled for clarity as shown in the figure. The curves represent the Bose-Einstein, mT -exponential, and double-exponential function fits to 0–5% central data for pions, kaons, and (anti)protons, respectively. The uncertainties are statistical and systematic added in quadrature.

Midrapidity (|y| < 0.1) transverse momentum spectra for (d) K− in Au+Au collisions at √sNN = 11.5 GeV for different centralities. The spectra for centralities other than 0–5% are scaled for clarity as shown in the figure. The curves represent the Bose-Einstein, mT -exponential, and double-exponential function fits to 0–5% central data for pions, kaons, and (anti)protons, respectively. The uncertainties are statistical and systematic added in quadrature.

Midrapidity (|y| < 0.1) transverse momentum spectra for (c) K+ in Au+Au collisions at √sNN = 11.5 GeV for different centralities. The spectra for centralities other than 0–5% are scaled for clarity as shown in the figure. The curves represent the Bose-Einstein, mT -exponential, and double-exponential function fits to 0–5% central data for pions, kaons, and (anti)protons, respectively. The uncertainties are statistical and systematic added in quadrature.

Midrapidity (|y| < 0.1) transverse momentum spectra for (f) p¯ in Au+Au collisions at √sNN = 11.5 GeV for different centralities. The spectra for centralities other than 0–5% are scaled for clarity as shown in the figure. The curves represent the Bose-Einstein, mT -exponential, and double-exponential function fits to 0–5% central data for pions, kaons, and (anti)protons, respectively. The uncertainties are statistical and systematic added in quadrature.

Midrapidity (|y| < 0.1) transverse momentum spectra for (e) p in Au+Au collisions at √sNN = 11.5 GeV for different centralities. The spectra for centralities other than 0–5% are scaled for clarity as shown in the figure. The curves represent the Bose-Einstein, mT -exponential, and double-exponential function fits to 0–5% central data for pions, kaons, and (anti)protons, respectively. The uncertainties are statistical and systematic added in quadrature.

Midrapidity (|y| < 0.1) transverse momentum spectra for (b) π− in Au+Au collisions at √sNN = 19.6 GeV for different centralities. The spectra for centralities other than 0–5% are scaled for clarity as shown in the figure. The curves represent the Bose-Einstein, mT -exponential, and double-exponential function fits to 0–5% central data for pions, kaons, and (anti)protons, respectively. The uncertainties are statistical and systematic added in quadrature.

Midrapidity (|y| < 0.1) transverse momentum spectra for (a) π+ in Au+Au collisions at √sNN = 19.6 GeV for different centralities. The spectra for centralities other than 0–5% are scaled for clarity as shown in the figure. The curves represent the Bose-Einstein, mT -exponential, and double-exponential function fits to 0–5% central data for pions, kaons, and (anti)protons, respectively. The uncertainties are statistical and systematic added in quadrature.

Midrapidity (|y| < 0.1) transverse momentum spectra for (d) K− in Au+Au collisions at √sNN = 19.6 GeV for different centralities. The spectra for centralities other than 0–5% are scaled for clarity as shown in the figure. The curves represent the Bose-Einstein, mT -exponential, and double-exponential function fits to 0–5% central data for pions, kaons, and (anti)protons, respectively. The uncertainties are statistical and systematic added in quadrature.

Midrapidity (|y| < 0.1) transverse momentum spectra for (c) K+ in Au+Au collisions at √sNN = 19.6 GeV for different centralities. The spectra for centralities other than 0–5% are scaled for clarity as shown in the figure. The curves represent the Bose-Einstein, mT -exponential, and double-exponential function fits to 0–5% central data for pions, kaons, and (anti)protons, respectively. The uncertainties are statistical and systematic added in quadrature.

Midrapidity (|y| < 0.1) transverse momentum spectra for (f) p¯ in Au+Au collisions at √sNN = 19.6 GeV for different centralities. The spectra for centralities other than 0–5% are scaled for clarity as shown in the figure. The curves represent the Bose-Einstein, mT -exponential, and double-exponential function fits to 0–5% central data for pions, kaons, and (anti)protons, respectively. The uncertainties are statistical and systematic added in quadrature.

Midrapidity (|y| < 0.1) transverse momentum spectra for (e) p in Au+Au collisions at √sNN = 19.6 GeV for different centralities. The spectra for centralities other than 0–5% are scaled for clarity as shown in the figure. The curves represent the Bose-Einstein, mT -exponential, and double-exponential function fits to 0–5% central data for pions, kaons, and (anti)protons, respectively. The uncertainties are statistical and systematic added in quadrature.

Midrapidity (|y| < 0.1) transverse momentum spectra for (b) π− in Au+Au collisions at √sNN = 27 GeV for different centralities. The spectra for centralities other than 0–5% are scaled for clarity as shown in the figure. The curves represent the Bose-Einstein, mT -exponential, and double-exponential function fits to 0–5% central data for pions, kaons, and (anti)protons, respectively. The uncertainties are statistical and systematic added in quadrature.

Midrapidity (|y| < 0.1) transverse momentum spectra for (a) π+ in Au+Au collisions at √sNN = 27 GeV for different centralities. The spectra for centralities other than 0–5% are scaled for clarity as shown in the figure. The curves represent the Bose-Einstein, mT -exponential, and double-exponential function fits to 0–5% central data for pions, kaons, and (anti)protons, respectively. The uncertainties are statistical and systematic added in quadrature.

Midrapidity (|y| < 0.1) transverse momentum spectra for (d) K− in Au+Au collisions at √sNN = 27 GeV for different centralities. The spectra for centralities other than 0–5% are scaled for clarity as shown in the figure. The curves represent the Bose-Einstein, mT -exponential, and double-exponential function fits to 0–5% central data for pions, kaons, and (anti)protons, respectively. The uncertainties are statistical and systematic added in quadrature.

Midrapidity (|y| < 0.1) transverse momentum spectra for (c) K+ in Au+Au collisions at √sNN = 27 GeV for different centralities. The spectra for centralities other than 0–5% are scaled for clarity as shown in the figure. The curves represent the Bose-Einstein, mT -exponential, and double-exponential function fits to 0–5% central data for pions, kaons, and (anti)protons, respectively. The uncertainties are statistical and systematic added in quadrature.

Midrapidity (|y| < 0.1) transverse momentum spectra for (f) p¯ in Au+Au collisions at √sNN = 27 GeV for different centralities. The spectra for centralities other than 0–5% are scaled for clarity as shown in the figure. The curves represent the Bose-Einstein, mT -exponential, and double-exponential function fits to 0–5% central data for pions, kaons, and (anti)protons, respectively. The uncertainties are statistical and systematic added in quadrature.

Midrapidity (|y| < 0.1) transverse momentum spectra for (e) p in Au+Au collisions at √sNN = 27 GeV for different centralities. The spectra for centralities other than 0–5% are scaled for clarity as shown in the figure. The curves represent the Bose-Einstein, mT -exponential, and double-exponential function fits to 0–5% central data for pions, kaons, and (anti)protons, respectively. The uncertainties are statistical and systematic added in quadrature.

Midrapidity (|y| < 0.1) transverse momentum spectra for (b) π− in Au+Au collisions at √sNN = 39 GeV for different centralities. The spectra for centralities other than 0–5% are scaled for clarity as shown in the figure. The curves represent the Bose-Einstein, mT -exponential, and double-exponential function fits to 0–5% central data for pions, kaons, and (anti)protons, respectively. The uncertainties are statistical and systematic added in quadrature.

Midrapidity (|y| < 0.1) transverse momentum spectra for (a) π+ in Au+Au collisions at √sNN = 39 GeV for different centralities. The spectra for centralities other than 0–5% are scaled for clarity as shown in the figure. The curves represent the Bose-Einstein, mT -exponential, and double-exponential function fits to 0–5% central data for pions, kaons, and (anti)protons, respectively. The uncertainties are statistical and systematic added in quadrature.

Midrapidity (|y| < 0.1) transverse momentum spectra for (d) k- in Au+Au collisions at √sNN = 39 GeV for different centralities. The spectra for centralities other than 0–5% are scaled for clarity as shown in the figure. The curves represent the Bose-Einstein, mT -exponential, and double-exponential function fits to 0–5% central data for pions, kaons, and (anti)protons, respectively. The uncertainties are statistical and systematic added in quadrature.

Midrapidity (|y| < 0.1) transverse momentum spectra for (c) k+ in Au+Au collisions at √sNN = 39 GeV for different centralities. The spectra for centralities other than 0–5% are scaled for clarity as shown in the figure. The curves represent the Bose-Einstein, mT -exponential, and double-exponential function fits to 0–5% central data for pions, kaons, and (anti)protons, respectively. The uncertainties are statistical and systematic added in quadrature.

Midrapidity (|y| < 0.1) transverse momentum spectra for (f) pbar in Au+Au collisions at √sNN = 39 GeV for different centralities. The spectra for centralities other than 0–5% are scaled for clarity as shown in the figure. The curves represent the Bose-Einstein, mT -exponential, and double-exponential function fits to 0–5% central data for pions, kaons, and (anti)protons, respectively. The uncertainties are statistical and systematic added in quadrature.

Midrapidity (|y| < 0.1) transverse momentum spectra for (e) p in Au+Au collisions at √sNN = 39 GeV for different centralities. The spectra for centralities other than 0–5% are scaled for clarity as shown in the figure. The curves represent the Bose-Einstein, mT -exponential, and double-exponential function fits to 0–5% central data for pions, kaons, and (anti)protons, respectively. The uncertainties are statistical and systematic added in quadrature.

Centrality dependence of dN/dy normalized by ⟨Npart⟩/2 for π+, π−, K+, K−, p, and p ̄ at midrapidity (|y|<0.1) in Au+Au collisions at √sNN = 7.7 GeV. Errors shown are the quadrature sum of statistical and systematic uncertainties. For clarity, ⟨Npart⟩ uncertainties are not added in quadrature.

Centrality dependence of dN/dy normalized by ⟨Npart⟩/2 for π+, π−, K+, K−, p, and p ̄ at midrapidity (|y|<0.1) in Au+Au collisions at √sNN = 11.5 GeV. Errors shown are the quadrature sum of statistical and systematic uncertainties. For clarity, ⟨Npart⟩ uncertainties are not added in quadrature.

Centrality dependence of dN/dy normalized by ⟨Npart⟩/2 for π+, π−, K+, K−, p, and p ̄ at midrapidity (|y|<0.1) in Au+Au collisions at √sNN = 19.6 GeV. Errors shown are the quadrature sum of statistical and systematic uncertainties. For clarity, ⟨Npart⟩ uncertainties are not added in quadrature.

Centrality dependence of dN/dy normalized by ⟨Npart⟩/2 for π+, π−, K+, K−, p, and p ̄ at midrapidity (|y|<0.1) in Au+Au collisions at √sNN = 27 GeV. Errors shown are the quadrature sum of statistical and systematic uncertainties. For clarity, ⟨Npart⟩ uncertainties are not added in quadrature.

Centrality dependence of dN/dy normalized by ⟨Npart⟩/2 for π+, π−, K+, K−, p, and p ̄ at midrapidity (|y|<0.1) in Au+Au collisions at √sNN = 39 GeV. Errors shown are the quadrature sum of statistical and systematic uncertainties. For clarity, ⟨Npart⟩ uncertainties are not added in quadrature.

Centrality dependences of <pT> for π+, π−, K+, K−, p, and p ̄ at midrapidity (|y|<0.1) in Au+Au collisions at √sNN = 7.7 GeV. Errors shown are quadrature sum of statistical and systematic uncertainties where the latter dominates.

Centrality dependences of <pT> for π+, π−, K+, K−, p, and p ̄ at midrapidity (|y|<0.1) in Au+Au collisions at √sNN = 11.5 GeV. Errors shown are quadrature sum of statistical and systematic uncertainties where the latter dominates.

Centrality dependences of <pT> for π+, π−, K+, K−, p, and p ̄ at midrapidity (|y|<0.1) in Au+Au collisions at √sNN = 19.6 GeV. Errors shown are quadrature sum of statistical and systematic uncertainties where the latter dominates.

Centrality dependences of <pT> for π+, π−, K+, K−, p, and p ̄ at midrapidity (|y|<0.1) in Au+Au collisions at √sNN = 27 GeV. Errors shown are quadrature sum of statistical and systematic uncertainties where the latter dominates.

Centrality dependences of <pT> for π+, π−, K+, K−, p, and p ̄ at midrapidity (|y|<0.1) in Au+Au collisions at √sNN = 39 GeV. Errors shown are quadrature sum of statistical and systematic uncertainties where the latter dominates.

Variation of π−/π+, K−/K+, and p ̄/p ratios as a function of ⟨Npart⟩ at midrapidity (|y| < 0.1) in Au+Au collisions at 7.7 GeV. Errors shown are the quadrature sum of statistical and systematic uncertainties where the latter dominates.

Variation of π−/π+, K−/K+, and p ̄/p ratios as a function of ⟨Npart⟩ at midrapidity (|y| < 0.1) in Au+Au collisions at 11.5 GeV. Errors shown are the quadrature sum of statistical and systematic uncertainties where the latter dominates.

Variation of π−/π+, K−/K+, and p ̄/p ratios as a function of ⟨Npart⟩ at midrapidity (|y| < 0.1) in Au+Au collisions at 19.6 GeV. Errors shown are the quadrature sum of statistical and systematic uncertainties where the latter dominates.

Variation of π−/π+, K−/K+, and p ̄/p ratios as a function of ⟨Npart⟩ at midrapidity (|y| < 0.1) in Au+Au collisions at 27 GeV. Errors shown are the quadrature sum of statistical and systematic uncertainties where the latter dominates.

Variation of π−/π+, K−/K+, and p ̄/p ratios as a function of ⟨Npart⟩ at midrapidity (|y| < 0.1) in Au+Au collisions at 39 GeV. Errors shown are the quadrature sum of statistical and systematic uncertainties where the latter dominates.

Variation of K−/π−, p ̄/π−, K+/π+, and p/π+ ratios as a function of ⟨Npart⟩ at midrapidity (|y| < 0.1) in Au+Au collisions at 7.7 GeV. Errors shown are the quadrature sum of statistical and systematic uncertainties where the latter dominates.

Variation of K−/π−, p ̄/π−, K+/π+, and p/π+ ratios as a function of ⟨Npart⟩ at midrapidity (|y| < 0.1) in Au+Au collisions at 11.5 GeV. Errors shown are the quadrature sum of statistical and systematic uncertainties where the latter dominates.

Variation of K−/π−, p ̄/π−, K+/π+, and p/π+ ratios as a function of ⟨Npart⟩ at midrapidity (|y| < 0.1) in Au+Au collisions at 19.6 GeV. Errors shown are the quadrature sum of statistical and systematic uncertainties where the latter dominates.

Variation of K−/π−, p ̄/π−, K+/π+, and p/π+ ratios as a function of ⟨Npart⟩ at midrapidity (|y| < 0.1) in Au+Au collisions at 27 GeV. Errors shown are the quadrature sum of statistical and systematic uncertainties where the latter dominates.

Variation of K−/π−, p ̄/π−, K+/π+, and p/π+ ratios as a function of ⟨Npart⟩ at midrapidity (|y| < 0.1) in Au+Au collisions at 39 GeV. Errors shown are the quadrature sum of statistical and systematic uncertainties where the latter dominates.

The midrapidity (|y| < 0.1) dN/dy normalized by ⟨Npart⟩/2 as a function of √sNN for π±, K±, and p and p ̄ in 0–5% Au+Au collisions at BES energies. Errors shown are the quadrature sum of statistical and systematic uncertainties where the latter dominates.

⟨mT⟩ − m of π±, K±, and p and p ̄ as a function of √sNN . Midrapidity (|y| < 0.1) results are shown for 0–5% central Au+Au collisions at BES energies. The errors shown are the quadrature sum of statistical and systematic uncertainties where the latter dominates.

π−/π+, K−/K+, and p ̄/p ratios at midrapidity (|y| < 0.1) in central 0–5% Au+Au collisions at √sNN = 7.7, 11.5, 19.6, 27, and 39 GeV. Errors shown are the quadrature sum of statistical and systematic uncertainties where the latter dominates.

K/π ratio at midrapidity (|y| < 0.1) for central 0–5% Au+Au collisions at √sNN = 7.7, 11.5, 19.6, 27, and 39 GeV. Errors shown are the quadrature sum of statistical and systematic uncertainties where the latter dominates.

The GCE model particle yields fits shown along with standard deviations for Au+Au 7.7 and Au+Au 39 GeV in 0–5% central collisions. Uncertainties on experimental data represent statistical and systematic uncertainties added in quadrature.

The GCE model particle ratios fits shown along with standard deviations for Au+Au 7.7 and Au+Au 39 GeV in 0–5% central collisions. Uncertainties on experimental data represent statistical and systematic uncertainties added in quadrature.

The SCE model particle yields fits shown along with standard deviations for Au+Au 7.7 and Au+Au 39 GeV in 0–5% central collisions. Uncertainties on experimental data represent statistical and systematic uncertainties added in quadrature.

The SCE model particle ratios fits shown along with standard deviations for Au+Au 7.7 and Au+Au 39 GeV in 0–5% central collisions. Uncertainties on experimental data represent statistical and systematic uncertainties added in quadrature.

Chemical freeze-out parameter γS plotted vs ⟨Npart⟩ in GCE for particle yields fit. Uncertainties represent systematic errors.

Chemical freeze-out parameter μB plotted vs ⟨Npart⟩ in GCE for particle yields fit. Uncertainties represent systematic errors.

Chemical freeze-out parameter μS plotted vs ⟨Npart⟩ in GCE for particle yields fit. Uncertainties represent systematic errors.

Chemical freeze-out parameter Tch plotted vs ⟨Npart⟩ in GCE for particle yields fit. Uncertainties represent systematic errors.

Chemical freeze-out parameter R plotted vs ⟨Npart⟩ in GCE for particle yields fit. Uncertainties represent systematic errors.

Ratio of chemical freeze-out parameter γS between results from particle yield fits to particle ratio fits in GCE plotted vs ⟨Npart⟩. Uncertainties represent systematic errors.

Ratio of chemical freeze-out parameter μB between results from particle yield fits to particle ratio fits in GCE plotted vs ⟨Npart⟩. Uncertainties represent systematic errors.

Ratio of chemical freeze-out parameter μS between results from particle yield fits to particle ratio fits in GCE plotted vs ⟨Npart⟩. Uncertainties represent systematic errors.

Ratio of chemical freeze-out parameter Tch between results from particle yield fits to particle ratio fits in GCE plotted vs ⟨Npart⟩. Uncertainties represent systematic errors.

Chemical freeze-out parameter γS plotted vs ⟨Npart⟩ in SCE for particle yields fit. Uncertainties represent systematic errors.

Chemical freeze-out parameter μB plotted vs ⟨Npart⟩ in SCE for particle yields fit. Uncertainties represent systematic errors.

Chemical freeze-out parameter Tch plotted vs ⟨Npart⟩ in SCE for particle yields fit. Uncertainties represent systematic errors.

Chemical freeze-out parameter R plotted vs ⟨Npart⟩ in SCE for particle yields fit. Uncertainties represent systematic errors.

Ratio of chemical freeze-out parameter γS between yield and ratio fits in SCE plotted vs ⟨Npart⟩. Uncertainties represent systematic errors.

Ratio of chemical freeze-out parameter μB between yield and ratio fits in SCE plotted vs ⟨Npart⟩. Uncertainties represent systematic errors.

Ratio of chemical freeze-out parameter Tch between yield and ratio fits in SCE plotted vs ⟨Npart⟩. Uncertainties represent systematic errors.

Ratio of chemical freeze-out parameter γS between GCE and SCE results using particle ratios in fits plotted vs ⟨Npart⟩. Uncertainties represent systematic errors.

Ratio of chemical freeze-out parameter μB between GCE and SCE results using particle ratios in fits plotted vs ⟨Npart⟩. Uncertainties represent systematic errors.

Ratio of chemical freeze-out parameter Tch between GCE and SCE results using particle ratios in fits plotted vs ⟨Npart⟩. Uncertainties represent systematic errors.

Ratio of chemical freeze-out parameter γS between GCE and SCE results using particle yields in fits plotted vs ⟨Npart⟩. Uncertainties represent systematic errors.

Ratio of chemical freeze-out parameter μB between GCE and SCE results using particle yields in fits plotted vs ⟨Npart⟩. Uncertainties represent systematic errors.

Ratio of chemical freeze-out parameter Tch between GCE and SCE results using particle yields in fits plotted vs ⟨Npart⟩. Uncertainties represent systematic errors.

Ratio of chemical freeze-out parameter R between GCE and SCE results using particle yields in fits plotted vs ⟨Npart⟩. Uncertainties represent systematic errors.

Extracted chemical freeze-out temperature vs baryon chemical potential for (a) GCE and (b) SCE cases using particle yields as input for fitting. Curves represent two model predictions [81,82]. The gray bands represent the theoretical prediction ranges of the Cleymans et al. model [81]. Uncertainties represent systematic errors.

Extracted chemical freeze-out temperature vs baryon chemical potential for (a) GCE and (b) SCE cases using particle yields as input for fitting. Curves represent two model predictions [81,82]. The gray bands represent the theoretical prediction ranges of the Cleymans et al. model [81]. Uncertainties represent systematic errors.

Extracted chemical freeze-out temperature vs baryon chemical potential for (a) GCE and (b) SCE cases using particle yields as input for fitting. Curves represent two model predictions [81,82]. The gray bands represent the theoretical prediction ranges of the Cleymans et al. model [81]. Uncertainties represent systematic errors.

"Choice on constraints: Extracted chemical freeze-out temperatures shown in panels (a), (c), and (e) and baryon chemical potentials shown in panels (b), (d), and (f) for GCE using particle yields as input for fitting, respectively, for Au+Au collisions at √sNN = 7.7, 19.6, and 39 GeV. Results are compared for three initial conditions: μQ = 0, μQ constrained to B/2Q value, and μQ constrained to B/2Q along with μS constrained to 0. Uncertainties represent systematic errors."

"Choice on constraints: Extracted chemical freeze-out temperatures shown in panels (a), (c), and (e) and baryon chemical potentials shown in panels (b), (d), and (f) for GCE using particle yields as input for fitting, respectively, for Au+Au collisions at √sNN = 7.7, 19.6, and 39 GeV. Results are compared for three initial conditions: μQ = 0, μQ constrained to B/2Q value, and μQ constrained to B/2Q along with μS constrained to 0. Uncertainties represent systematic errors."

"Choice on constraints: Extracted chemical freeze-out temperatures shown in panels (a), (c), and (e) and baryon chemical potentials shown in panels (b), (d), and (f) for GCE using particle yields as input for fitting, respectively, for Au+Au collisions at √sNN = 7.7, 19.6, and 39 GeV. Results are compared for three initial conditions: μQ = 0, μQ constrained to B/2Q value, and μQ constrained to B/2Q along with μS constrained to 0. Uncertainties represent systematic errors."

"Choice on constraints: Extracted chemical freeze-out temperatures shown in panels (a), (c), and (e) and baryon chemical potentials shown in panels (b), (d), and (f) for GCE using particle yields as input for fitting, respectively, for Au+Au collisions at √sNN = 7.7, 19.6, and 39 GeV. Results are compared for three initial conditions: μQ = 0, μQ constrained to B/2Q value, and μQ constrained to B/2Q along with μS constrained to 0. Uncertainties represent systematic errors."

"Choice on constraints: Extracted chemical freeze-out temperatures shown in panels (a), (c), and (e) and baryon chemical potentials shown in panels (b), (d), and (f) for GCE using particle yields as input for fitting, respectively, for Au+Au collisions at √sNN = 7.7, 19.6, and 39 GeV. Results are compared for three initial conditions: μQ = 0, μQ constrained to B/2Q value, and μQ constrained to B/2Q along with μS constrained to 0. Uncertainties represent systematic errors."

"Choice on constraints: Extracted chemical freeze-out temperatures shown in panels (a), (c), and (e) and baryon chemical potentials shown in panels (b), (d), and (f) for GCE using particle yields as input for fitting, respectively, for Au+Au collisions at √sNN = 7.7, 19.6, and 39 GeV. Results are compared for three initial conditions: μQ = 0, μQ constrained to B/2Q value, and μQ constrained to B/2Q along with μS constrained to 0. Uncertainties represent systematic errors."

"Choice on including more particles: Extracted chemical freeze-out parameters (a) Tch, (b) μB, and (c) γS along with (d) χ2/ndf for GCE using particle yields as input for fitting. Results are compared for Au+Au collisions at √sNN = 39 GeV for four different sets of particle yields used in fitting. Uncertainties represent systematic errors."

"Choice on including more particles: Extracted chemical freeze-out parameters (a) Tch, (b) μB, and (c) γS along with (d) χ2/ndf for GCE using particle yields as input for fitting. Results are compared for Au+Au collisions at √sNN = 39 GeV for four different sets of particle yields used in fitting. Uncertainties represent systematic errors."

"Choice on including more particles: Extracted chemical freeze-out parameters (a) Tch, (b) μB, and (c) γS along with (d) χ2/ndf for GCE using particle yields as input for fitting. Results are compared for Au+Au collisions at √sNN = 39 GeV for four different sets of particle yields used in fitting. Uncertainties represent systematic errors."

"Choice on including more particles: Extracted chemical freeze-out parameters (a) Tch, (b) μB, and (c) γS along with (d) χ2/ndf for GCE using particle yields as input for fitting. Results are compared for Au+Au collisions at √sNN = 39 GeV for four different sets of particle yields used in fitting. Uncertainties represent systematic errors."

"Blast wave model fits of π±, K±, p and p p¯ T spectra in 0–5% central Au+Au collisions at √sNN = (a) 7.7, (b) 11.5, (c) 19.6, (d) 27, and (e) 39 GeV. Uncertainties on experimental data represent statistical and systematic uncertainties added in quadrature. Here, the uncertainties are smaller than the symbol size."

"Blast wave model fits of π±, K±, p and p p¯ T spectra in 0–5% central Au+Au collisions at √sNN = (a) 7.7, (b) 11.5, (c) 19.6, (d) 27, and (e) 39 GeV. Uncertainties on experimental data represent statistical and systematic uncertainties added in quadrature. Here, the uncertainties are smaller than the symbol size."

"Blast wave model fits of π±, K±, p and p p¯ T spectra in 0–5% central Au+Au collisions at √sNN = (a) 7.7, (b) 11.5, (c) 19.6, (d) 27, and (e) 39 GeV. Uncertainties on experimental data represent statistical and systematic uncertainties added in quadrature. Here, the uncertainties are smaller than the symbol size."

"Blast wave model fits of π±, K±, p and p p¯ T spectra in 0–5% central Au+Au collisions at √sNN = (a) 7.7, (b) 11.5, (c) 19.6, (d) 27, and (e) 39 GeV. Uncertainties on experimental data represent statistical and systematic uncertainties added in quadrature. Here, the uncertainties are smaller than the symbol size."

"Blast wave model fits of π±, K±, p and p p¯ T spectra in 0–5% central Au+Au collisions at √sNN = (a) 7.7, (b) 11.5, (c) 19.6, (d) 27, and (e) 39 GeV. Uncertainties on experimental data represent statistical and systematic uncertainties added in quadrature. Here, the uncertainties are smaller than the symbol size."

"Variation of Tkin with <β> for different energies and centralities. The centrality increases from left to right for a given energy. The data points other than BES energies are taken from Refs. [43,66]. Uncertainties represent systematic uncertainties."

"Variation of Tkin with <β> for different energies and centralities. The centrality increases from left to right for a given energy. The data points other than BES energies are taken from Refs. [43,66]. Uncertainties represent systematic uncertainties."

"Variation of Tkin with <β> for different energies and centralities. The centrality increases from left to right for a given energy. The data points other than BES energies are taken from Refs. [43,66]. Uncertainties represent systematic uncertainties."

"Variation of Tkin with <β> for different energies and centralities. The centrality increases from left to right for a given energy. The data points other than BES energies are taken from Refs. [43,66]. Uncertainties represent systematic uncertainties."

"Variation of Tkin with <β> for different energies and centralities. The centrality increases from left to right for a given energy. The data points other than BES energies are taken from Refs. [43,66]. Uncertainties represent systematic uncertainties."

" (a) Energy dependence of kinetic and chemical freezeout temperatures for central heavy-ion collisions. The curves represent various theoretical predictions [81,82]. (b) Energy dependence of average transverse radial flow velocity for central heavy-ion collisions. The data points other than BES energies are taken from Refs. [43,53–64,66] and references therein. The BES data points are for 0–5% central collisions, AGS energies are mostly for 0–5%, SPS energies are for mostly 0–7%, and top RHIC and LHC energies are for 0–5% central collisions. Uncertainties represent systematic uncertainties."

Centrality dependence of flow harmonics from $v_2$ to $v_9$.

Centrality dependence of flow harmonics from $v_2$ to $v_9$.

Centrality dependence of flow harmonics from $v_2$ to $v_9$.

Centrality dependence of flow harmonics from $v_2$ to $v_9$.

Centrality dependence of flow harmonics from $v_2$ to $v_9$.

Centrality dependence of linear and non-linear flow modes.

Centrality dependence of linear and non-linear flow modes.

Centrality dependence of linear and non-linear flow modes.

Centrality dependence of linear and non-linear flow modes.

Centrality dependence of linear and non-linear flow modes.

Centrality dependence of linear and non-linear flow modes.

Centrality dependence of symmetry-plane correlations.

Centrality dependence of symmetry-plane correlations.

Centrality dependence of symmetry-plane correlations.

Centrality dependence of symmetry-plane correlations.

Centrality dependence of non-linear flow mode coefficients.

Centrality dependence of non-linear flow mode coefficients.

Centrality dependence of non-linear flow mode coefficients.

Centrality dependence of non-linear flow mode coefficients.

Centrality dependence of non-linear flow mode coefficients.

Ratio of inclusive $\psi(2{\rm S})$ over J/$\psi$ yields times their branching ratios as a function of $y_{\rm cms}$ in p--Pb collisions at $\sqrt{s_{\rm NN}}$ = 8.16 TeV .The first uncertainty is statistical, while the second is the systematic.

Ratio of inclusive $\psi(2{\rm S})$ over J/$\psi$ yields times their branching ratios as a function of $p_{\rm T}$ in p--Pb collisions at $\sqrt{s_{\rm NN}}$ = 8.16 TeV at forward rapidity.The first uncertainty is statistical, while the second is the systematic.

Ratio of inclusive $\psi(2{\rm S})$ over J/$\psi$ yields times their branching ratios as a function of $p_{\rm T}$ in p--Pb collisions at $\sqrt{s_{\rm NN}}$ = 8.16 TeV at backward rapidity.The first uncertainty is statistical, while the second is the systematic.

Inclusive $\psi(2{\rm S})$ $R_{\rm pPb}$ values as a function of $y_{\rm {cms}}$ at $\sqrt{s_{\rm NN}}$ = 8.16 TeV. The first uncertainty is statistical, the second is the uncorrelated systematic, while the third one is a correlated systematic uncertainty.

Inclusive $\psi(2{\rm S})$ $R_{\rm pPb}$ values as a function of $p_{\rm T}$ at $\sqrt{s_{\rm NN}}$ = 8.16 TeV at forward rapidity. The first uncertainty is statistical, the second is the uncorrelated systematic, while the third one is a correlated systematic uncertainty.

Inclusive $\psi(2{\rm S})$ $R_{\rm pPb}$ values as a function of $p_{\rm T}$ at $\sqrt{s_{\rm NN}}$ = 8.16 TeV at backward rapidity. The first uncertainty is statistical, the second is the uncorrelated systematic, while the third one is a correlated systematic uncertainty.

Double ratio of $\psi$(2S) and J/$\psi$ cross sections in p--Pb and pp collisions as a function of rapidity, at $\sqrt{s_{\mathrm{NN}}} = 8.16$ TeV. The first uncertainty is statistical, while the second is the systematic.

Double ratio of $\psi$(2S) and J/$\psi$ cross sections in p--Pb and pp collisions as a function of $p_{\rm T}$, at $\sqrt{s_{\mathrm{NN}}} = 8.16$ TeV at forward rapidity. The first uncertainty is statistical, while the second is the systematic.

Double ratio of $\psi$(2S) and J/$\psi$ cross sections in p--Pb and pp collisions as a function of $p_{\rm T}$, at $\sqrt{s_{\mathrm{NN}}} = 8.16$ TeV at backward rapidity. The first uncertainty is statistical, while the second is the systematic.

Inclusive $\psi(2{\rm S})$ cross section times branching ratio for backward and forward rapidity, in p--Pb collisions at $\sqrt{s_{\rm NN}}$ = 8.16 TeV. The first uncertainty is statistical, the second is the uncorrelated systematic, while the third one is a correlated systematic uncertainty.

Inclusive $\psi(2{\rm S})$ $R_{\rm pPb}$ values for backward and forward rapidity at $\sqrt{s_{\rm NN}}$ = 8.16 TeV. The first uncertainty is statistical, the second is the uncorrelated systematic, while the third one is a correlated systematic uncertainty.

The $pK\pi$ invariant mass distributions (Counts per 10 MeV/c^2 bin) for wrong-sign combinations, scaled by 1/3, in Au+Au collisions at $\sqrt{s_{\rm NN}}$ = 200\,GeV for 10--80\% centrality class.

The $\Lambda_c^{\pm}$ invariant yields measured in the 10-80% centrality class for the different $p_T$ bins, in Au+Au collisions at √sNN = 200 GeV. The $p_T$ bins are 2.5 < $p_T$ < 3.5 GeV/c, 3.5 < $p_T$ < 5.0 GeV/c and 5.0 < $p_T$ < 8.0 GeV/c.

The $\Lambda_c/D^0$ yield ratio measured in the 10-80% centrality class for the different $p_T$ bins, in Au+Au collisions at √sNN = 200 GeV. The $p_T$ bins are 2.5 < $p_T$ < 3.5 GeV/c, 3.5 < $p_T$ < 5.0 GeV/c and 5.0 < $p_T$ < 8.0 GeV/c.

The $\Lambda_c/D^0$ yield ratio measured in 3.0 < $p_T$ < 6.0 GeV/c for different centrality classes, in Au+Au collisions at √sNN = 200 GeV. The centrality bins are 50-80%, 20-50% and 0-20%.

The centrality dependence of e+e− invariant mass spectra within the STAR acceptance from Au+Au collisions and U+U collisions for pair pT < 0.15 GeV/c. The vertical bars on data points depict the statistical uncertainties, while the systematic uncertainties are shown as gray boxes. The hadronic cocktail yields from U+U collisions are ∼5%–12% higher than those from Au+Au collisions in given centrality bins; thus only cocktails for Au+Au collisions are shown here as solid lines, with shaded bands representing the systematic uncertainties for clarity.

The centrality dependence of e+e− invariant mass spectra within the STAR acceptance from Au+Au collisions and U+U collisions for pair pT < 0.15 GeV/c. The vertical bars on data points depict the statistical uncertainties, while the systematic uncertainties are shown as gray boxes. The hadronic cocktail yields from U+U collisions are ∼5%–12% higher than those from Au+Au collisions in given centrality bins; thus only cocktails for Au+Au collisions are shown here as solid lines, with shaded bands representing the systematic uncertainties for clarity.

The centrality dependence of e+e− invariant mass spectra within the STAR acceptance from Au+Au collisions and U+U collisions for pair pT < 0.15 GeV/c. The vertical bars on data points depict the statistical uncertainties, while the systematic uncertainties are shown as gray boxes. The hadronic cocktail yields from U+U collisions are ∼5%–12% higher than those from Au+Au collisions in given centrality bins; thus only cocktails for Au+Au collisions are shown here as solid lines, with shaded bands representing the systematic uncertainties for clarity.

The corresponding ratios of data over cocktail for Figure 1.

The corresponding ratios of data over cocktail for Figure 1.

The corresponding ratios of data over cocktail for Figure 1.

The corresponding ratios of data over cocktail for Figure 1.

The corresponding ratios of data over cocktail for Figure 1.

The corresponding ratios of data over cocktail for Figure 1.

The e+e− pair pT distributions within the STAR acceptance for different mass regions in 60%–80% Au+Au and U+U collisions compared to cocktails. The systematic uncertainties of the data are shown as gray boxes. The gray bands depict the systematic uncertainties of the cocktails.

The e+e− pair pT distributions within the STAR acceptance for different mass regions in 60%–80% Au+Au and U+U collisions compared to cocktails. The systematic uncertainties of the data are shown as gray boxes. The gray bands depict the systematic uncertainties of the cocktails.

The e+e− pair pT distributions within the STAR acceptance for different mass regions in 60%–80% Au+Au and U+U collisions compared to cocktails. The systematic uncertainties of the data are shown as gray boxes. The gray bands depict the systematic uncertainties of the cocktails.

The e+e− pair pT distributions within the STAR acceptance for different mass regions in 60%–80% Au+Au and U+U collisions compared to cocktails. The systematic uncertainties of the data are shown as gray boxes. The gray bands depict the systematic uncertainties of the cocktails.

The e+e− pair pT distributions within the STAR acceptance for different mass regions in 60%–80% Au+Au and U+U collisions compared to cocktails. The systematic uncertainties of the data are shown as gray boxes. The gray bands depict the systematic uncertainties of the cocktails.

The e+e− pair pT distributions within the STAR acceptance for different mass regions in 60%–80% Au+Au and U+U collisions compared to cocktails. The systematic uncertainties of the data are shown as gray boxes. The gray bands depict the systematic uncertainties of the cocktails.

The low-pT (pT < 0.15 GeV/c) e+e− excess mass spectra (data − cocktail) within the STAR acceptance in 60%–80% for Au + Au and U + U collisions, compared with a broadened rho meson model calculation [8]. The contributions of rho and phi from the photonuclear process are shown, as are the contributions of photon-photon process from two models [33,34]. The model calculations are for Au + Au collisions in the corresponding centrality bins.

The low-pT (pT < 0.15 GeV/c) e+e− excess mass spectra (data − cocktail) within the STAR acceptance in 60%–80% for Au + Au and U + U collisions, compared with a broadened rho meson model calculation [8]. The contributions of rho and phi from the photonuclear process are shown, as are the contributions of photon-photon process from two models [33,34]. The model calculations are for Au + Au collisions in the corresponding centrality bins.

The low-pT (pT < 0.15 GeV/c) e+e− excess mass spectra (data − cocktail) within the STAR acceptance in 40%–60% for Au + Au and U + U collisions, compared with a broadened rho meson model calculation [8]. The contributions of rho and phi from the photonuclear process are shown, as are the contributions of photon-photon process from two models [33,34]. The model calculations are for Au + Au collisions in the corresponding centrality bins.

The low-pT (pT < 0.15 GeV/c) e+e− excess mass spectra (data − cocktail) within the STAR acceptance in 40%–60% for Au + Au and U + U collisions, compared with a broadened rho meson model calculation [8]. The contributions of rho and phi from the photonuclear process are shown, as are the contributions of photon-photon process from two models [33,34]. The model calculations are for Au + Au collisions in the corresponding centrality bins.

The centrality dependence of integrated excess yields in the mass regions of 0.4–0.76, 0.76–1.2, and 1.2–2.6 GeV=c2 in Au + Au and U + U collisions. The centrality dependence of hadronic cocktail yields in the mass region of 0.76–1.2 GeV/c2 in both collisions is also shown for comparison. The systematic uncertainties are shown as gray boxes.

The centrality dependence of integrated excess yields in the mass regions of 0.4–0.76, 0.76–1.2, and 1.2–2.6 GeV=c2 in Au + Au and U + U collisions. The centrality dependence of hadronic cocktail yields in the mass region of 0.76–1.2 GeV/c2 in both collisions is also shown for comparison. The systematic uncertainties are shown as gray boxes.

The centrality dependence of integrated excess yields in the mass regions of 0.4–0.76, 0.76–1.2, and 1.2–2.6 GeV=c2 in Au + Au and U + U collisions. The centrality dependence of hadronic cocktail yields in the mass region of 0.76–1.2 GeV/c2 in both collisions is also shown for comparison. The systematic uncertainties are shown as gray boxes.

The centrality dependence of integrated excess yields in the mass regions of 0.4–0.76, 0.76–1.2, and 1.2–2.6 GeV=c2 in Au + Au and U + U collisions. The centrality dependence of hadronic cocktail yields in the mass region of 0.76–1.2 GeV/c2 in both collisions is also shown for comparison. The systematic uncertainties are shown as gray boxes.

The centrality dependence of integrated excess yields in the mass regions of 0.4–0.76, 0.76–1.2, and 1.2–2.6 GeV=c2 in Au + Au and U + U collisions. The centrality dependence of hadronic cocktail yields in the mass region of 0.76–1.2 GeV/c2 in both collisions is also shown for comparison. The systematic uncertainties are shown as gray boxes.

The centrality dependence of integrated excess yields in the mass regions of 0.4–0.76, 0.76–1.2, and 1.2–2.6 GeV/c2 in Au + Au and U + U collisions. The centrality dependence of hadronic cocktail yields in the mass region of 0.76–1.2 GeV/c2 in both collisions is also shown for comparison. The systematic uncertainties are shown as gray boxes.

The p2T distributions of excess yields within the STAR acceptance in the mass regions of 0.4–0.76 GeV/c2 in 60%–80% Au + Au and U + U collisions. The systematic uncertainties are shown as gray boxes. The solid and dotted lines are exponential fits to the data in Au + Au and U + U collisions, respectively. The dot-dashed and dot-dot- dashed lines represent the p2T distributions for the photon-photon process from two models [33,34] within the STAR acceptance in 60%–80% Au + Au collisions. The dashed lines illustrate the corresponding p2T distributions for e+e− pairs from the model [33] traversing 1 fm in a constant magnetic field of 10^14 T perpendicular to the beam line.

The p2T distributions of excess yields within the STAR acceptance in the mass regions of 0.4–0.76 GeV/c2 in 60%–80% Au + Au and U + U collisions. The systematic uncertainties are shown as gray boxes. The solid and dotted lines are exponential fits to the data in Au + Au and U + U collisions, respectively. The dot-dashed and dot-dot- dashed lines represent the p2T distributions for the photon-photon process from two models [33,34] within the STAR acceptance in 60%–80% Au + Au collisions. The dashed lines illustrate the corresponding p2T distributions for e+e− pairs from the model [33] traversing 1 fm in a constant magnetic field of 10^14 T perpendicular to the beam line.

The p2T distributions of excess yields within the STAR acceptance in the mass regions of 0.76–1.2 GeV/c2 in 60%–80% Au + Au and U + U collisions. The systematic uncertainties are shown as gray boxes. The solid and dotted lines are exponential fits to the data in Au + Au and U + U collisions, respectively. The dot-dashed and dot-dot- dashed lines represent the p2T distributions for the photon-photon process from two models [33,34] within the STAR acceptance in 60%–80% Au + Au collisions. The dashed lines illustrate the corresponding p2T distributions for e+e− pairs from the model [33] traversing 1 fm in a constant magnetic field of 10^14 T perpendicular to the beam line.

The p2T distributions of excess yields within the STAR acceptance in the mass regions of 0.76–1.2 GeV/c2 in 60%–80% Au + Au and U + U collisions. The systematic uncertainties are shown as gray boxes. The solid and dotted lines are exponential fits to the data in Au + Au and U + U collisions, respectively. The dot-dashed and dot-dot- dashed lines represent the p2T distributions for the photon-photon process from two models [33,34] within the STAR acceptance in 60%–80% Au + Au collisions. The dashed lines illustrate the corresponding p2T distributions for e+e− pairs from the model [33] traversing 1 fm in a constant magnetic field of 10^14 T perpendicular to the beam line.

The p2T distributions of excess yields within the STAR acceptance in the mass regions of 1.2–2.6 GeV/c2 in 60%–80% Au + Au and U + U collisions. The systematic uncertainties are shown as gray boxes. The solid and dotted lines are exponential fits to the data in Au + Au and U + U collisions, respectively. The dot-dashed and dot-dot- dashed lines represent the p2T distributions for the photon-photon process from two models [33,34] within the STAR acceptance in 60%–80% Au + Au collisions. The dashed lines illustrate the corresponding p2T distributions for e+e− pairs from the model [33] traversing 1 fm in a constant magnetic field of 10^14 T perpendicular to the beam line.

The p2T distributions of excess yields within the STAR acceptance in the mass regions of 1.2–2.6 GeV/c2 in 60%–80% Au + Au and U + U collisions. The systematic uncertainties are shown as gray boxes. The solid and dotted lines are exponential fits to the data in Au + Au and U + U collisions, respectively. The dot-dashed and dot-dot- dashed lines represent the p2T distributions for the photon-photon process from two models [33,34] within the STAR acceptance in 60%–80% Au + Au collisions. The dashed lines illustrate the corresponding p2T distributions for e+e− pairs from the model [33] traversing 1 fm in a constant magnetic field of 10^14 T perpendicular to the beam line.

The corresponding sqrt(p2T) of excess yields. The vertical bars on data points are the combined statistical and systematic uncertainties.

The corresponding sqrt(p2T) of excess yields. The vertical bars on data points are the combined statistical and systematic uncertainties.

Transverse momentum spectra for charged pions at midrapidity.

Transverse momentum spectra for protons at midrapidity.

Transverse momentum spectra for antiprotons at midrapidity.

Transverse momentum spectra for kaons at midrapidity.

dN/dy of pi+ and proton, normalised by N_part/2.

dN/dy of pi+ and proton, normalised by N_part/2.

dN/dy of pi+ and proton, normalised by N_part/2.

dN/dy of K+ and K-, normalised by N_part/2.

dN/dy of K+ and K-, normalised by N_part/2.

dN/dy of K+ and K-, normalised by N_part/2.

Transverse momentum of pi+, K+ and proton.

Transverse momentum of pi+, K+ and proton.

Transverse momentum of pi+, K+ and proton.

K-/K+, K-/pi-, p/pi+ and K+/pi+ vs. N_part.

K-/K+, K-/pi-, p/pi+ and K+/pi+ vs. N_part.

K-/K+, K-/pi-, p/pi+ and K+/pi+ vs. N_part.

dN/d_eta, normalised by N_part/2, as a function of sqrt(S_NN).

dN/dy, normalised by N_part/2; and m_t - m of pi+ and pi- as functions of sqrt(S_NN).

dN/dy, normalised by N_part/2; and m_t - m of K+ and K- as functions of sqrt(S_NN).

pi-/pi+, pbar/p, K-/K+, K+/pi+ and K-/pi- as functions of sqrt(S_NN).

v1 as a function of eta and eta/y-beam from different experiments and methods.

v2 as a function of p_T for different particle species and experiments.

v2 as a function of s_NN from different experiments.

T-mu_B plot for heavy-ion collisions at different s_NN.

$\Omega + \bar{\Omega}$ invariant mass spectra from Cu+Cu $\sqrt{s_{NN}} = 200$ GeV collisions, where $|y| < 0.5$. The uncertainties on the spectra points are statistical and systematic combined.

$K^0_S$ invariant mass spectra from Au+Au $\sqrt{s_{NN}} = 200$ GeV collisions, where $|y| < 0.5$. The uncertainties on the spectra points are statistical and systematic combined.

$\Lambda$ and $\bar{\Lambda}$ invariant mass spectra from Au+Au $\sqrt{s_{NN}} = 200$ GeV collisions, where $|y| < 0.5$. The $\Lambda$ and $\bar{\Lambda}$ yields have not been feed down subtracted from weak decays. The uncertainties on the spectra points are statistical and systematic combined.

$K^0_S$, $\Lambda$, and $\bar{\Lambda}$, spectra divided by $\langle N_{part} \rangle$ for Cu+Cu $0 − 10\%$ ($\langle N_{part} \rangle \sim 99$) and Au+Au $20 − 40\%$ ($\langle N_{part} \rangle \sim 141$) $\sqrt{s_{NN}} = 200$ GeV collisions, where $|y| < 0.5$. The $\Lambda$ and $\bar{\Lambda}$ yields have been feed down subtracted from weak decays. The uncertainties on the spectra points are statistical and systematic; for clarity the uncertainty on $\langle N_{part} \rangle$ has not been included. The curves show the functions described in the text used to extract $dN/dy$.

$\Xi + \bar{\Xi}$ spectra divided by $\langle N_{part} \rangle$ for Cu+Cu $0 − 10\%$ ($\langle N_{part} \rangle \sim 99$) $\sqrt{s_{NN}} = 200$ GeV collisions, where $|y| < 0.5$. The uncertainties on the spectra points are statistical and systematic; for clarity the uncertainty on $\langle N_{part} \rangle$ has not been included. The curves show the functions described in the text used to extract $dN/dy$.

$\Xi + \bar{\Xi}$ spectra divided by $\langle N_{part} \rangle$ for Au+Au $20 − 40\%$ ($\langle N_{part} \rangle \sim 141$) $\sqrt{s_{NN}} = 200$ GeV collisions, where $|y| < 0.5$. The uncertainties on the spectra points are statistical and systematic; for clarity the uncertainty on $\langle N_{part} \rangle$ has not been included. The curves show the functions described in the text used to extract $dN/dy$.

$\Omega + \bar{\Omega}$ spectra divided by $\langle N_{part} \rangle$ for Cu+Cu $0 − 10\%$ ($\langle N_{part} \rangle \sim 99$) $\sqrt{s_{NN}} = 200$ GeV collisions, where $|y| < 0.5$. The uncertainties on the spectra points are statistical and systematic; for clarity the uncertainty on $\langle N_{part} \rangle$ has not been included. The curves show the functions described in the text used to extract $dN/dy$.

$\Omega + \bar{\Omega}$ spectra divided by $\langle N_{part} \rangle$ for Au+Au $20 − 40\%$ ($\langle N_{part} \rangle \sim 141$) $\sqrt{s_{NN}} = 200$ GeV collisions, where $|y| < 0.5$. The uncertainties on the spectra points are statistical and systematic; for clarity the uncertainty on $\langle N_{part} \rangle$ has not been included. The curves show the functions described in the text used to extract $dN/dy$.

The enhancement factor for (multi-) strange particles in Cu+Cu $\sqrt{s_{NN}} = 200$ GeV collisions, where $|y| < 0.5$. The $\Lambda$, and $\bar{\Lambda}$ yields have been feed-down subtracted in all cases. The black bars show the normalization uncertainties, and the uncertainties for the heavy-ion points are the combined statistical and systematic errors. Curves described in the text, where $B_{K} = 2.0$, $B_{\Lambda} = 2.4$, $B_{\Xi} = 5.0$ and $B_{\Omega} = 12.1$.

The enhancement factor for (multi-) strange particles in Au+Au $\sqrt{s_{NN}} = 200$ GeV collisions, where $|y| < 0.5$. The $\Lambda$, and $\bar{\Lambda}$ yields have been feed-down subtracted in all cases. The black bars show the normalization uncertainties, and the uncertainties for the heavy-ion points are the combined statistical and systematic errors. Curves described in the text, where $B_{K} = 2.0$, $B_{\Lambda} = 2.4$, $B_{\Xi} = 5.0$ and $B_{\Omega} = 12.1$.

Ratio of particle yields in central Cu+Cu and midcentral Au+Au collisions when $\langle N_{part} \rangle = 99$ in each case for $|y| < 0.5$. $\pi$ yields are from elsewhere [17]. Boxed uncertainties are from the Glauber calculations and are correlated for every particle. $\langle N_{part>1} \rangle$ refers to the parameterization shown by eq. 1, while the EPOS and AMPT models are described in the text. The default settings are used for each model. The vertical lines show the remaining independent statistical and systematic uncertainties.

The invariant yields $d^2N/(2\pi p_T dp_T dy)$ of $K + p(\bar{p})$ in central Au+Au collisions, and NLO calculations with AKK [9] and DSS [10] FF. The uncertainty of yields due to the scale dependence as evaluated in [10] is about a factor of 2. Bars and boxes (bands) represent statistical and systematic uncertainties, respectively.

The invariant yields $d^2N/(2\pi p_T dp_T dy)$ of $\rho^0$ in central Au+Au collisions, and NLO calculations with AKK [9] and DSS [10] FF. The uncertainty of yields due to the scale dependence as evaluated in [10] is about a factor of 2. Bars and boxes (bands) represent statistical and systematic uncertainties, respectively.

The invariant yields $d^2N/(2\pi p_T dp_T dy)$ of $K^0_S$ in central Au+Au collisions, and NLO calculations with AKK [9] and DSS [10] FF. The uncertainty of yields due to the scale dependence as evaluated in [10] is about a factor of 2. Bars and boxes (bands) represent statistical and systematic uncertainties, respectively.

Yield ratios $\pi^−/\pi^+$, $\bar{p}/p$, $K^−/K^+$, $p/\pi^+$, $\bar{p}/\pi^−$, and $K^{\pm}/\pi^{\pm}$ versus $p_T$ in p+p collisions, and nominal NLO calculations with AKK [9] and DSS [10] FFs without theoretical uncertainties. Bars and boxes (bands) represent statistical and systematic uncertainties, respectively.

Yield ratio $K^0_S$)/$\pi^{\pm}$ versus $p_T$ in p+p collisions, and nominal NLO calculations with AKK [9] and DSS [10] FFs without theoretical uncertainties. Bars and boxes (bands) represent statistical and systematic uncertainties, respectively.

Yield ratios $p/\pi^+$ and $\bar{p}/\pi^−$ versus $p_T$ in central Au+Au collisions [14] with updated uncertainties at high $p_T$. Bars and boxes (bands) represent statistical and systematic uncertainties, respectively.

$R_{AA}$ of $\pi^{\pm} in central Au+Au collisions as a function of $p_T$. Bars and boxes (bands) represent statistical and systematic uncertainties, respectively. The height of the band at unity represents the normalization uncertainty.

$R_{AA}$ of $K^{\pm}+p(\bar{p})$ in central Au+Au collisions as a function of $p_T$. Bars and boxes (bands) represent statistical and systematic uncertainties, respectively. The height of the band at unity represents the normalization uncertainty.

$R_{AA}$ of $K^0_S$ in central Au+Au collisions as a function of $p_T$. The curves are the calculations for $K^0_S$ $R_{AA}$ with and without jet conversion in medium [15]. Bars and boxes (bands) represent statistical and systematic uncertainties, respectively. The height of the band at unity represents the normalization uncertainty.

$R_{AA}$ of $\rho^0$ in central Au+Au collisions as a function of $p_T$. Bars and boxes (bands) represent statistical and systematic uncertainties, respectively. The height of the band at unity represents the normalization uncertainty.

The ratios of $R_{AA}[K^{\pm} + p(\bar{p})]$ to $R_{AA}(\pi^{\pm})$. The boxes and shaded bands represent the systematic uncertainties for $R_{AA}[K^{\pm} +p(\bar{p})]/R_{AA}(\pi^{\pm})$, respectively.

The ratios of $R_{AA}(K^− + \bar{p})$ to $R_{AA}(K^+ p)$. The systematic uncertainties for $R_{AA}(K^− + \bar{p})/R_{AA}(K^+ +p)$ are 2–12% and left off for clarity.

The ratios of $R_{AA}[\rho^0]$ to $R_{AA}(\pi^{\pm})$. The boxes and shaded bands represent the systematic uncertainties for $R_{AA}(\rho^0)/R_{AA}(\pi^{\pm})$, respectively.

I_AA vs xi varying integration range

I_AA vs xi varying integration range