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

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.

'Identified transverse momentum spectra of $\pi^{-}$ at midrapidity (|y| < 0.1) in U+U collisions at $\sqrt{s_{\rm NN}}$ = 193 GeV'

'Identified transverse momentum spectra of $K^{-}$ at midrapidity (|y| < 0.1) in U+U collisions at $\sqrt{s_{\rm NN}}$ = 193 GeV'

'Identified transverse momentum spectra of $\bar{p}$ at midrapidity (|y| < 0.1) in U+U collisions at $\sqrt{s_{\rm NN}}$ = 193 GeV'

'Average transverse momentum as a function of $\langle$N$_{part}$$\rangle$ for $\pi^{+}$ at midrapidity (|y| < 0.1) in U+U collisions at $\sqrt{s_{\rm NN}}$ = 193 GeV.'

'Average transverse momentum as a function of $\langle$N$_{part}$$\rangle$ for $K^{+}$ at midrapidity (|y| < 0.1) in U+U collisions at $\sqrt{s_{\rm NN}}$ = 193 GeV.'

'Average transverse momentum as a function of $\langle$N$_{part}$$\rangle$ for p at midrapidity (|y| < 0.1) in U+U collisions at $\sqrt{s_{\rm NN}}$ = 193 GeV.'

'dN/dy scaled by 0.5 × $\langle$N$_{part}$$\rangle$ as a function of $\langle$N$_{part}$$\rangle$ for $\pi^{+}$ at midrapidity (|y| < 0.1) in U+U collisions at $\sqrt{s_{\rm NN}}$ = 193 GeV.'

'dN/dy scaled by 0.5 × $\langle$N$_{part}$$\rangle$ as a function of $\langle$N$_{part}$$\rangle$ for $K^{+}$ at midrapidity (|y| < 0.1) in U+U collisions at $\sqrt{s_{\rm NN}}$ = 193 GeV.'

'dN/dy scaled by 0.5 × $\langle$N$_{part}$$\rangle$ as a function of $\langle$N$_{part}$$\rangle$ for p at midrapidity (|y| < 0.1) in U+U collisions at $\sqrt{s_{\rm NN}}$ = 193 GeV.'

'dN/dy scaled by 0.5 × $\langle$N$_{part}$$\rangle$ as a function of $\langle$N$_{part}$$\rangle$ for $\bar{p}$ at midrapidity (|y| < 0.1) in U+U collisions at $\sqrt{s_{\rm NN}}$ = 193 GeV.'

'$\pi^{−}$/$\pi^{+}$ ratio as a function of $\langle$N$_{part}$$\rangle$ at midrapidity (|y| < 0.1) in U+U collisions at $\sqrt{s_{\rm NN}}$ = 193 GeV.'

'$K^{−}$/$K^{+}$ ratio as a function of $\langle$N$_{part}$$\rangle$ at midrapidity (|y| < 0.1) in U+U collisions at $\sqrt{s_{\rm NN}}$ = 193 GeV.'

'$\bar{p}$/p ratio as a function of $\langle$N$_{part}$$\rangle$ at midrapidity (|y| < 0.1) in U+U collisions at $\sqrt{s_{\rm NN}}$ = 193 GeV.'

'$K^{+}$/$\pi^{+}$ ratio as a function of $\langle$N$_{part}$$\rangle$ at midrapidity (|y| < 0.1) in U+U collisions at $\sqrt{s_{\rm NN}}$ = 193 GeV.'

'$K^{-}$/$\pi^{-}$ ratio as a function of $\langle$N$_{part}$$\rangle$ at midrapidity (|y| < 0.1) in U+U collisions at $\sqrt{s_{\rm NN}}$ = 193 GeV.'

'p/$\pi^{+}$ ratio as a function of $\langle$N$_{part}$$\rangle$ at midrapidity (|y| < 0.1) in U+U collisions at $\sqrt{s_{\rm NN}}$ = 193 GeV.'

'$\bar{p}$/$\pi^{-}$ ratio as a function of $\langle$N$_{part}$$\rangle$ at midrapidity (|y| < 0.1) in U+U collisions at $\sqrt{s_{\rm NN}}$ = 193 GeV.'

The $\phi$-meson integrated nuclear modification factors $R_{AB}$ measured as a function of $N_{part}$

The comparison of $\phi$ meson (a) $v_2$ and (b) $v_2/(\epsilon N_{part}^{1/3})$ measured as a function of $p_T$ in Cu+Au collisions at $\sqrt{s}$ = 200 GeV and U+U collisions at $\sqrt{s}$ = 193 GeVat midrapidity (|η| < 0.35)

The comparison of elliptic flow (a) $v_2$ values measured for $\phi$ mesons as a function of $p_T$ in $0\%–20\%$ Cu+Au collisions

The comparison of elliptic flow (b) $v_2$ values measured for $\phi$ mesons as a function of $KE_T$ in $0\%–20\%$ Cu+Au collisions

The comparison of elliptic flow (c) $v_2/n_q$ values measured for $\phi$ mesons as a function of $KE_T/n_q$ in $0\%–20\%$ Cu+Au collisions

The comparison of elliptic flow (a) $v_2$ values measured for $\phi$ mesons as a function of $p_T$ in $20\%–60\%$ Cu+Au collisions

The comparison of elliptic flow (b) $v_2$ values measured for $\phi$ mesons as a function of $KE_T$ in $20\%–60\%$ Cu+Au collisions

The comparison of elliptic flow (c) $v_2/n_q$ values measured for $\phi$ mesons as a function of $KE_T/n_q$ in $20\%–60\%$ Cu+Au collisions

The comparison of elliptic flow (a) $v_2$ values measured for $\phi$ mesons as a function of $p_T$ in $0\%–50\%$ U+U collisions

The comparison of elliptic flow (b) $v_2$ values measured for $\phi$ mesons as a function of $KE_T$ in $0\%–50\%$ U+U collisions

The comparison of elliptic flow (c) $v_2/n_q$ values measured for $\phi$ mesons as a function of $KE_T/n_q$ in $0\%–50\%$ U+U collisions

Example of combinatorial background subtracted invariant mass distributions and the extracted yields as a function of $\cos \theta^*$ for $\phi$ and $K^{*0}$ mesons. \textbf{a)} example of $\phi \rightarrow K^+ + K^-$ invariant mass distributions, with combinatorial background subtracted, integrated over $\cos \theta^*$; \textbf{b)} example of $K^{*0} (\overline{K^{*0}}) \rightarrow K^{-} \pi^{+} (K^{+} \pi^{-})$ invariant mass distributions, with combinatorial background subtracted, integrated over $\cos \theta^*$; \textbf{c)} extracted yields of $\phi$ as a function of $\cos \theta^*$; \textbf{d)} extracted yields of $K^{*0}$ as a function of $\cos \theta^*$.

Example of combinatorial background subtracted invariant mass distributions and the extracted yields as a function of $\cos \theta^*$ for $\phi$ and $K^{*0}$ mesons. \textbf{a)} example of $\phi \rightarrow K^+ + K^-$ invariant mass distributions, with combinatorial background subtracted, integrated over $\cos \theta^*$; \textbf{b)} example of $K^{*0} (\overline{K^{*0}}) \rightarrow K^{-} \pi^{+} (K^{+} \pi^{-})$ invariant mass distributions, with combinatorial background subtracted, integrated over $\cos \theta^*$; \textbf{c)} extracted yields of $\phi$ as a function of $\cos \theta^*$; \textbf{d)} extracted yields of $K^{*0}$ as a function of $\cos \theta^*$.

Example of combinatorial background subtracted invariant mass distributions and the extracted yields as a function of $\cos \theta^*$ for $\phi$ and $K^{*0}$ mesons. \textbf{a)} example of $\phi \rightarrow K^+ + K^-$ invariant mass distributions, with combinatorial background subtracted, integrated over $\cos \theta^*$; \textbf{b)} example of $K^{*0} (\overline{K^{*0}}) \rightarrow K^{-} \pi^{+} (K^{+} \pi^{-})$ invariant mass distributions, with combinatorial background subtracted, integrated over $\cos \theta^*$; \textbf{c)} extracted yields of $\phi$ as a function of $\cos \theta^*$; \textbf{d)} extracted yields of $K^{*0}$ as a function of $\cos \theta^*$.

Efficiency corrected $\phi$-meson yields as a function of cos$\theta$* and corresponding fits with Eq.1 in the method section.

Efficiency and acceptance corrected $K^{*0}$-meson yields as a function of cos$\theta$* and corresponding fits with Eq.4 in the method section.

$\phi$-meson $\rho_{00}$ obtained from 1st- and 2nd-order event planes. The red stars (gray squares) show the $\phi$-meson $\rho_{00}$ as a function of beam energy, obtained with the 2nd-order (1st-order) EP.

$\phi$-meson $\rho_{00}$ with respect to different quantization axes. $\phi$-meson $\rho_{00}$ as a function of beam energy, for the out-of-plane direction (stars) and the in-plane direction (diamonds). Curves are fits based on theoretical calculations with a $\phi$-meson field. The corresponding $G_s^{(y)}$ values obtained from the fits are shown in the legend.

$\rho_{00}$ as a function of transverse momentum for $\phi$ for different collision energies. The gray squares and red stars are results obtained with the 1st- and 2nd-order EP, respectively.

$\rho_{00}$ as a function of transverse momentum for $\phi$ for different collision energies. The gray squares and red stars are results obtained with the 1st- and 2nd-order EP, respectively.

$\rho_{00}$ as a function of transverse momentum for $\phi$ for different collision energies. The gray squares and red stars are results obtained with the 1st- and 2nd-order EP, respectively.

$\rho_{00}$ as a function of transverse momentum for $\phi$ for different collision energies. The gray squares and red stars are results obtained with the 1st- and 2nd-order EP, respectively.

$\rho_{00}$ as a function of transverse momentum for $\phi$ for different collision energies. The gray squares and red stars are results obtained with the 1st- and 2nd-order EP, respectively.

$\rho_{00}$ as a function of transverse momentum for $\phi$ for different collision energies. The gray squares and red stars are results obtained with the 1st- and 2nd-order EP, respectively.

$\rho_{00}$ as a function of transverse momentum for $\phi$ for different collision energies. The gray squares and red stars are results obtained with the 1st- and 2nd-order EP, respectively.

$\rho_{00}$ as a function of transverse momentum for $\phi$ for different collision energies. The gray squares and red stars are results obtained with the 1st- and 2nd-order EP, respectively.

$\rho_{00}$ as a function of transverse momentum for $\phi$ for different collision energies. The gray squares and red stars are results obtained with the 1st- and 2nd-order EP, respectively.

$\rho_{00}$ as a function of transverse momentum for $\phi$ for different collision energies. The gray squares and red stars are results obtained with the 1st- and 2nd-order EP, respectively.

$\rho_{00}$ as a function of transverse momentum for $\phi$ for different collision energies. The gray squares and red stars are results obtained with the 1st- and 2nd-order EP, respectively.

$\rho_{00}$ as a function of transverse momentum for $\phi$ for different collision energies. The gray squares and red stars are results obtained with the 1st- and 2nd-order EP, respectively.

$\rho_{00}$ as a function of transverse momentum for $K^{*0}$ for different collision energies.

$\rho_{00}$ as a function of transverse momentum for $K^{*0}$ for different collision energies.

$\rho_{00}$ as a function of transverse momentum for $K^{*0}$ for different collision energies.

$\rho_{00}$ as a function of transverse momentum for $K^{*0}$ for different collision energies.

$\rho_{00}$ as a function of transverse momentum for $K^{*0}$ for different collision energies.

$\rho_{00}$ as a function of transverse momentum for $K^{*0}$ for different collision energies.

$\rho_{00}$ as a function of transverse momentum for $K^{*0}$ for different collision energies.

$\rho_{00}$ as a function of centrality for $\phi$ (upper panels) and $K^{*0}$ (lower panels). The solid squares and stars are results for the $\phi$ meson, obtained with the 1st- and 2nd-order EP, respectively. The solid circles are results for the $K^{*0}$ meson, obtained with the 2nd-order EP.}

$\rho_{00}$ as a function of centrality for $\phi$ (upper panels) and $K^{*0}$ (lower panels). The solid squares and stars are results for the $\phi$ meson, obtained with the 1st- and 2nd-order EP, respectively. The solid circles are results for the $K^{*0}$ meson, obtained with the 2nd-order EP.}

$\rho_{00}$ as a function of centrality for $\phi$ (upper panels) and $K^{*0}$ (lower panels). The solid squares and stars are results for the $\phi$ meson, obtained with the 1st- and 2nd-order EP, respectively. The solid circles are results for the $K^{*0}$ meson, obtained with the 2nd-order EP.}

$\rho_{00}$ as a function of centrality for $\phi$ (upper panels) and $K^{*0}$ (lower panels). The solid squares and stars are results for the $\phi$ meson, obtained with the 1st- and 2nd-order EP, respectively. The solid circles are results for the $K^{*0}$ meson, obtained with the 2nd-order EP.}

$\rho_{00}$ as a function of centrality for $\phi$ (upper panels) and $K^{*0}$ (lower panels). The solid squares and stars are results for the $\phi$ meson, obtained with the 1st- and 2nd-order EP, respectively. The solid circles are results for the $K^{*0}$ meson, obtained with the 2nd-order EP.}

$\rho_{00}$ as a function of centrality for $\phi$ (upper panels) and $K^{*0}$ (lower panels). The solid squares and stars are results for the $\phi$ meson, obtained with the 1st- and 2nd-order EP, respectively. The solid circles are results for the $K^{*0}$ meson, obtained with the 2nd-order EP.}

$\rho_{00}$ as a function of centrality for $\phi$ (upper panels) and $K^{*0}$ (lower panels). The solid squares and stars are results for the $\phi$ meson, obtained with the 1st- and 2nd-order EP, respectively. The solid circles are results for the $K^{*0}$ meson, obtained with the 2nd-order EP.}

Global spin alignment measurement of $\phi$ and $K^{*0}$ vector mesons in Au+Au collisions at 0-20\% centrality. The solid squares and stars are results for the $\phi$ meson, obtained with the 1st- and 2nd-order EP, respectively. The solid circles are results for $K^{*0}$-meson, obtained with the 2nd-order EP.

Global spin alignment measurement of $\phi$ and $K^{*0}$ vector mesons in Au+Au collisions at 0-20\% centrality. The solid squares and stars are results for the $\phi$ meson, obtained with the 1st- and 2nd-order EP, respectively. The solid circles are results for $K^{*0}$-meson, obtained with the 2nd-order EP.

Invariant yield of direct photons, every 20% centrality

Invariant yield of nonprompt photons, every 20% centrality

$T_{eff}$ vs charged particle multiplicity

Integrated direct photon yield of 1-5GeV vs charged particle multiplicity

Integrated nonprompt direct photon yield of different $p_T$ integration window vs charged particle multiplicity

Scaling factor $\alpha$ vs charged particle multiplicity

Direct photon spectra for minimum bias (0-86%) Au+Au collision at $\sqrt{s_{NN}} = 62.4$ GeV (a) and $39$ GeV (b). For $62.4$ GeV also centrality bins of 0-20% (c) and 20-40% (d) are shown. Data points are shown with statistical and systematic uncertainties, unless the central value is negative (arrows) or is consistent with zero within the statistical uncertainties (arrows with data point). In these cases upper limit with CL = 95$%$ are given.

Inverse slopes, $T_{eff}$ , obtained from fitting the combined data from central collisions shown in Fig. 11 is compared to the fit results of the individual data sets at $\sqrt{s_{NN}} = 62.4$ GeV, $200$ GeV, and $2760$ GeV. Also included is the value for $\sqrt{s_{NN}} = 39$ GeV obtained from fitting the min. bias data set in the lower $p_T$ range.

Integrated invariant direct-photon yields vs. charged particle multiplicity for $p_{T}$ integrated from (a) 0.4 GeV/c, (b) 1.0 GeV/c, (c) 1.5 GeV/c, and (d) 2.0 to 5.0 GeV/c for all available A+A datasets.

Integrated invariant direct-photon yields vs. charged particle multiplicity for $p_{T}$ integrated from (a) 0.4 GeV/c, (b) 1.0 GeV/c, (c) 1.5 GeV/c, and (d) 2.0 to 5.0 GeV/c for all available A+A datasets.

Integrated invariant direct-photon yields vs. charged particle multiplicity for $p_{T}$ integrated from (a) 0.4 GeV/c, (b) 1.0 GeV/c, (c) 1.5 GeV/c, and (d) 2.0 to 5.0 GeV/c for all available A+A datasets.

Integrated invariant direct-photon yields vs. charged particle multiplicity for $p_{T}$ integrated from (a) 0.4 GeV/c, (b) 1.0 GeV/c, (c) 1.5 GeV/c, and (d) 2.0 to 5.0 GeV/c for all available A+A datasets.

Integrated direct-photon yields from A+A collisions for $p_{T,min}$ of $5$ GeV/c (a) and $8$ GeV/c. The representation is the same as in Fig. 13. Also shown are the results from pQCD calculations scaled by $N_{coll}$

Integrated direct-photon yields from A+A collisions for $p_{T,min}$ of $5$ GeV/c (a) and $8$ GeV/c. The representation is the same as in Fig. 13. Also shown are the results from pQCD calculations scaled by $N_{coll}$

The $\alpha$ values extracted using fits to integrated direct photon yields. The dashed line gives the average $\alpha$ value for the 4 lower $p_{T,min}$ points. Also shown is a model calculation for $\alpha$ discussed in the text.

Azimuthal anisotropy $v_2\{BF\}$ as a function of transverse momentum $p_T$ in $d$+Au collisions at $\sqrt{s_{NN}} =$ 200 GeV.

Azimuthal anisotropy $v_2\{BB\}$ as a function of transverse momentum $p_T$ in $^3$He+Au collisions at $\sqrt{s_{NN}} =$ 200 GeV.

Azimuthal anisotropy $v_2\{BF\}$ as a function of transverse momentum $p_T$ in $^3$He+Au collisions at $\sqrt{s_{NN}} =$ 200 GeV.

Azimuthal anisotropy $v_2$ as a function of transverse momentum $p_T$ in $p$+$p$ collisions at $\sqrt{s_{NN}} =$ 200 GeV.

Azimuthal anisotropy $v_2$ as a function of centrality in $p$+Au collisions at $\sqrt{s_{NN}} =$ 200 GeV.

Azimuthal anisotropy $v_2$ as a function of centrality in $d$+Au collisions at $\sqrt{s_{NN}} =$ 200 GeV.

Azimuthal anisotropy $v_2$ as a function of centrality in $^3$He+Au collisions at $\sqrt{s_{NN}} =$ 200 GeV.

Azimuthal anisotropy $v_2$ as a function of charged particle multiplicity $dN_{ch}/d\eta$ in $p$+Au collisions at $\sqrt{s_{NN}} =$ 200 GeV.

Azimuthal anisotropy $v_2$ as a function of charged particle multiplicity $dN_{ch}/d\eta$ in $d$+Au collisions at $\sqrt{s_{NN}} =$ 200 GeV.

Azimuthal anisotropy $v_2$ as a function of charged particle multiplicity $dN_{ch}/d\eta$ in $^3$He+Au collisions at $\sqrt{s_{NN}} =$ 200 GeV.

Azimuthal anisotropy $v_2$ as a function of charged particle multiplicity $dN_{ch}/d\eta$ in $p$+$p$ collisions at $\sqrt{s_{NN}} =$ 200 GeV.

Azimuthal anisotropy ratio $R=v_2\{BF\}/v_2\{BB\}$ as a function of charged particle multiplicity $dN_{ch}/d\eta$ in $p$+Au collisions at $\sqrt{s_{NN}} =$ 200 GeV.

Azimuthal anisotropy ratio $R=v_2\{BF\}/v_2\{BB\}$ as a function of charged particle multiplicity $dN_{ch}/d\eta$ in $d$+Au collisions at $\sqrt{s_{NN}} =$ 200 GeV.

Azimuthal anisotropy ratio $R=v_2\{BF\}/v_2\{BB\}$ as a function of charged particle multiplicity $dN_{ch}/d\eta$ in $^3$He+Au collisions at $\sqrt{s_{NN}} =$ 200 GeV.

Azimuthal anisotropy ratio $R=v_2\{BF\}/v_2\{BB\}$ as a function of charged particle multiplicity $dN_{ch}/d\eta$ in $p$+$p$ collisions at $\sqrt{s_{NN}} =$ 200 GeV.

J/PSI(1S)-->MU+MU- nuclear modification in p+Al collisions as a function of rapidity. The statistical and systematic uncertainties vary point-to-point and are listed for each measured value. An additional global systematic uncertainty is provided in each column heading, which applies to all data points per column.

PSI(2S)-->MU+MU- nuclear modification in p+Au collisions as a function of rapidity. The statistical and systematic uncertainties vary point-to-point and are listed for each measured value. An additional global systematic uncertainty is provided in each column heading, which applies to all data points per column.

J/PSI(1S)-->MU+MU- nuclear modification in p+Au collisions as a function of rapidity. The statistical and systematic uncertainties vary point-to-point and are listed for each measured value. An additional global systematic uncertainty is provided in each column heading, which applies to all data points per column.

PSI(2S)-->MU+MU- nuclear modification in p+Au collisions as a function of centrality at forward rapidity. The statistical and systematic uncertainties vary point-to-point and are listed for each measured value. An additional global systematic uncertainty is provided in each column heading, which applies to all data points per column.

PSI(2S)-->MU+MU- nuclear modification in p+Au collisions as a function of centrality at backward rapidity. The statistical and systematic uncertainties vary point-to-point and are listed for each measured value. An additional global systematic uncertainty is provided in each column heading, which applies to all data points per column.

J/PSI(1S)-->MU+MU- nuclear modification in p+Au collisions as a function of centrality. The statistical and systematic uncertainties vary point-to-point and are listed for each measured value. An additional global systematic uncertainty is provided in each column heading, which applies to all data points per column.

Invariant yield ratio of PSI(2S)-->MU+MU- with respect to J/PSI(1S)-->MU+MU- in p+p collisions at rapidity 1.2 < |y| < 2.2. The statistical and systematic uncertainties vary point-to-point and are listed for each measured value. The global systematic uncertainty cancels when taking the ratio.

Invariant yield ratio of PSI(2S)-->MU+MU- with respect to J/PSI(1S)-->MU+MU- in p+Au collisions as a function of centrality at forward rapidity. The statistical and systematic uncertainties vary point-to-point and are listed for each measured value. The global systematic uncertainty cancels when taking the ratio.

Invariant yield ratio of PSI(2S)-->MU+MU- with respect to J/PSI(1S)-->MU+MU- in p+Au collisions as a function of centrality at backward rapidity. The statistical and systematic uncertainties vary point-to-point and are listed for each measured value. The global systematic uncertainty cancels when taking the ratio.

Fraction of bottom hadron decayed electrons in 0-80% Au+Au collisions at $\sqrt{s_{\rm NN}}$ = 200 GeV.

Fraction of bottom hadron decayed electrons in 0-20% and 20-40% Au+Au collisions at $\sqrt{s_{\rm NN}}$ = 200 GeV.

Fraction of bottom hadron decayed electrons in 40-80% Au+Au collisions at $\sqrt{s_{\rm NN}}$ = 200 GeV.

Nuclear modification factor $R_{\rm AA}$ of inclusive heavy flavor decayed electrons in 0-80% Au+Au collisions at $\sqrt{s_{\rm NN}}$ = 200 GeV. The 8% global uncertainty from the $p$+$p$ reference and a 8% gloabl uncertainty from $N_{\rm coll.}$ are not included in the table values.

Nuclear modification factor $R_{\rm AA}$ of bottom- and charm-decay electrons in 0-80% Au+Au collisions at $\sqrt{s_{\rm NN}}$ = 200 GeV. The 8% global uncertainty from the $p$+$p$ reference and a 8% gloabl uncertainty from $N_{\rm coll.}$ are not included in the table values.

Nuclear modification factor $R_{\rm AA}$ ratio of bottom- to charm-decay electrons in 0-80% Au+Au collisions at $\sqrt{s_{\rm NN}}$ = 200 GeV.

Null hypothesis for the nuclear modification factor $R_{\rm AA}$ ratio of bottom- to charm-decay electrons in 0-80% Au+Au collisions at $\sqrt{s_{\rm NN}}$ = 200 GeV.

Nuclear modification factor $R_{\rm CP}$ ratio of bottom- and charm-decay electrons in 0-80% Au+Au collisions at $\sqrt{s_{\rm NN}}$ = 200 GeV.

Null hypotheses for the nuclear modification factor $R_{\rm CP}$ ratio of bottom- and charm-decay electrons in 0-80% Au+Au collisions at $\sqrt{s_{\rm NN}}$ = 200 GeV.

Exponent according to the Eq. 5 as a function of transverse momenta extracted from p+Au/p+Al and $^{3}$HeAu/p+Au collision systems. The uncertainties from the $N_{coll}$ calculations and from the overall normalization of p+p cancel in these ratios

Nuclear modification factors in p+Al, p+Au, d+Au, and $^{3}$He+Au in five centrality bins and for inelastic collisions at $\sqrt{s}$ = 200 GeV. The right boxes are the uncertainties of the $N_{coll}$ collisions from the Glauber model, while the left box represents the overall normalization uncertainty from p+p collisions

Average $R_{xA}$ in two different $p_T$ regions versus the number of collisions (panels a,b) and the number of collisions per projectile participant (panels c,d). The tilted error bars represent the anti-correlated uncertainty on the y and x-axis due to the Ncoll calculations.

Average $R_{xA}$ in two different $p_T$ regions versus the number of collisions (panels a,b) and the number of collisions per projectile participant (panels c,d). The tilted error bars represent the anti-correlated uncertainty on the y and x-axis due to the Ncoll calculations. Centrality vs. the number of collisions is tabulated here.

Average $R_{xA}$ in two different $p_T$ regions versus the number of collisions (panels a,b) and the number of collisions per projectile participant (panels c,d). The tilted error bars represent the anti-correlated uncertainty on the y and x-axis due to the Ncoll calculations. Centrality vs. Average $R_{xA}$ is tabulated here.

Average $R_{xA}$ in two different $p_T$ regions versus the number of collisions (panels a,b) and the number of collisions per projectile participant (panels c,d). The tilted error bars represent the anti-correlated uncertainty on the y and x-axis due to the Ncoll calculations.

Average $R_{xA}$ in two different $p_T$ regions versus the number of collisions (panels a,b) and the number of collisions per projectile participant (panels c,d). The tilted error bars represent the anti-correlated uncertainty on the y and x-axis due to the Ncoll calculations. Centrality vs. the number of collisions is tabulated here.

Average $R_{xA}$ in two different $p_T$ regions versus the number of collisions (panels a,b) and the number of collisions per projectile participant (panels c,d). The tilted error bars represent the anti-correlated uncertainty on the y and x-axis due to the Ncoll calculations. Centrality vs. Average $R_{xA}$ is tabulated here.

Average $R_{xA}$ in two different $p_T$ regions versus the number of collisions (panels a,b) and the number of collisions per projectile participant (panels c,d). The tilted error bars represent the anti-correlated uncertainty on the y and x-axis due to the Ncoll calculations.

Average $R_{xA}$ in two different $p_T$ regions versus the number of collisions (panels a,b) and the number of collisions per projectile participant (panels c,d). The tilted error bars represent the anti-correlated uncertainty on the y and x-axis due to the Ncoll calculations. Centrality vs. the number of collisions per projectile participant is tabulated here.

Average $R_{xA}$ in two different $p_T$ regions versus the number of collisions (panels a,b) and the number of collisions per projectile participant (panels c,d). The tilted error bars represent the anti-correlated uncertainty on the y and x-axis due to the Ncoll calculations. Centrality vs. Average $R_{xA}$ is tabulated here.

Average $R_{xA}$ in two different $p_T$ regions versus the number of collisions (panels a,b) and the number of collisions per projectile participant (panels c,d). The tilted error bars represent the anti-correlated uncertainty on the y and x-axis due to the Ncoll calculations.

Average $R_{xA}$ in two different $p_T$ regions versus the number of collisions (panels a,b) and the number of collisions per projectile participant (panels c,d). The tilted error bars represent the anti-correlated uncertainty on the y and x-axis due to the Ncoll calculations. Centrality vs. the number of collisions per projectile participant is tabulated here.

Average $R_{xA}$ in two different $p_T$ regions versus the number of collisions (panels a,b) and the number of collisions per projectile participant (panels c,d). The tilted error bars represent the anti-correlated uncertainty on the y and x-axis due to the Ncoll calculations. Centrality vs. Average $R_{xA}$ is tabulated here.

$R_{xA}$ for inelastic collisions compared to three different nuclear PDF calculations and their uncertainties. The data points include the statistical and systematical uncertainties. Relevant data has been tabulated in Figure 7.

The upper panel (a) shows the <$R_{xA}$> above $p_T$ =8 Gev/c as a function $N_{coll}$/$N_{proj}$. The data are compared to predictions from [49] for the consequences of high-$x$ nucleon size fluctuations. The lower panels show the $R_{xA}$ as a function of $p_T$ for most central on left (b) and most peripheral on right panel (c) collisions. Panel (a) has been tabulated in Figure 10 and Panel (b) and (c) in Figure 9.

The upper panel (a) shows the <$R_{xA}$> above $p_T$ =8 Gev/c as a function $N_{coll}$/$N_{proj}$. The data are compared to predictions from [49] for the consequences of high-$x$ nucleon size fluctuations. The lower panels show the $R_{xA}$ as a function of $p_T$ for most central on left (b) and most peripheral on right panel (c) collisions. Panel (a) has been tabulated in Figure 10 and Panel (b) and (c) in Figure 9. Centrality vs. the number of collisions per projectile participant in panel (a) is tabulated here.

The upper panel (a) shows the <$R_{xA}$> above $p_T$ =8 Gev/c as a function $N_{coll}$/$N_{proj}$. The data are compared to predictions from [49] for the consequences of high-$x$ nucleon size fluctuations. The lower panels show the $R_{xA}$ as a function of $p_T$ for most central on left (b) and most peripheral on right panel (c) collisions. Panel (a) has been tabulated in Figure 10 and Panel (b) and (c) in Figure 9. Centrality vs. Average $R_{xA}$ in panel (a) is tabulated here.

The upper panel (a) shows the <$R_{xA}$> above $p_T$ =8 Gev/c as a function $N_{coll}$/$N_{proj}$. The data are compared to predictions from [49] for the consequences of high-$x$ nucleon size fluctuations. The lower panels show the $R_{xA}$ as a function of $p_T$ for most central on left (b) and most peripheral on right panel (c) collisions. Panel (a) has been tabulated in Figure 10 and Panel (b) and (c) in Figure 9. Panel(b) and panel(c) are tabulated here.

Integrated yields for 1--2 GeV/c in panel (a) and 2--3 GeV/c in panel (b) as a function of charged particle multiplicity density at mid rapidity.

Integrated yields for 1--2 GeV/c in panel (a) and 2--3 GeV/c in panel (b) as a function of charged particle multiplicity density at mid rapidity. Centrality vs. charged particle multiplicity density at mid rapidity is tabulated here.

Integrated yields for 1--2 GeV/c in panel (a) and 2--3 GeV/c in panel (b) as a function of charged particle multiplicity density at mid rapidity. Centrality vs. integrated yield is tabulated here.

Integrated yields for 1--2 GeV/c in panel (a) and 2--3 GeV/c in panel (b) as a function of charged particle multiplicity density at mid rapidity.

Integrated yields for 1--2 GeV/c in panel (a) and 2--3 GeV/c in panel (b) as a function of charged particle multiplicity density at mid rapidity. Centrality vs. charged particle multiplicity density at mid rapidity is tabulated here.

Integrated yields for 1--2 GeV/c in panel (a) and 2--3 GeV/c in panel (b) as a function of charged particle multiplicity density at mid rapidity. Centrality vs. integrated yield is tabulated here.

Forward neutron single spin asymmetries in p+p, p+Al and p+Au collisions as a function of pT in bins of xF correlated with hard collision activity

Forward neutron single spin asymmetries in p+p, p+Al and p+Au collisions as a function of pT in bins of xF without hard collision activity

Forward neutron single spin asymmetries in p+p collisions as a function of xF in bins of pT (anti-)correlated with hard collision activity

Forward neutron single spin asymmetries in p+Al collisions as a function of xF in bins of pT (anti-)correlated with hard collision activity

Forward neutron single spin asymmetries in p+Au collisions as a function of xF in bins of pT (anti-)correlated with hard collision activity

$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

$v_3$ vs $p_T$, p+Au at 200 GeV, 0-5% central, BBCS-FVTXS-CNT detector combination

$v_3$ vs $p_T$, d+Au at 200 GeV, 0-5% central, BBCS-FVTXS-CNT detector combination

$v_3$ vs $p_T$, 3He+Au at 200 GeV, 0-5% central, BBCS-FVTXS-CNT detector combination

$v_2$ vs $p_T$, p+Au at 200 GeV, 0-5% central, FVTXS-CNT-FVTXN detector combination

$v_2$ vs $p_T$, d+Au at 200 GeV, 0-5% central, FVTXS-CNT-FVTXN detector combination

$v_2$ vs $p_T$, 3He+Au at 200 GeV, 0-5% central, FVTXS-CNT-FVTXN detector combination

$v_2$ ratio vs $p_T$, p/d/3He+Au at 200 GeV, 0-5% central

$v_3$ vs $p_T$, p+Au at 200 GeV, 0-5% central, FVTXS-CNT-FVTXN detector combination

$v_3$ vs $p_T$, d+Au at 200 GeV, 0-5% central, FVTXS-CNT-FVTXN detector combination

$v_3$ vs $p_T$, 3He+Au at 200 GeV, 0-5% central, FVTXS-CNT-FVTXN detector combination

$c_2$ vs $\eta$, p+Au at 200 GeV, 0-5% central

$c_2$ vs $\eta$, d+Au at 200 GeV, 0-5% central

$c_2$ vs $\eta$, 3He+Au at 200 GeV, 0-5% central

$c_3$ vs $\eta$, p+Au at 200 GeV, 0-5% central

$c_3$ vs $\eta$, d+Au at 200 GeV, 0-5% central

$c_3$ vs $\eta$, 3He+Au at 200 GeV, 0-5% central

$C(\Delta\phi)$ vs $\Delta\phi$, p+p at 200 GeV, minimum bias, BBCS-FVTXS correlation

$C(\Delta\phi)$ vs $\Delta\phi$, p+Au at 200 GeV, 0-5% central, BBCS-FVTXS correlation

$C(\Delta\phi)$ vs $\Delta\phi$, d+Au at 200 GeV, 0-5% central, BBCS-FVTXS correlation

$C(\Delta\phi)$ vs $\Delta\phi$, 3He+Au at 200 GeV, 0-5% central, BBCS-FVTXS correlation

$C(\Delta\phi)$ vs $\Delta\phi$, p+p at 200 GeV, minimum bias, BBCS-FVTXN correlation

$C(\Delta\phi)$ vs $\Delta\phi$, p+Au at 200 GeV, 0-5% central, BBCS-FVTXN correlation

$C(\Delta\phi)$ vs $\Delta\phi$, d+Au at 200 GeV, 0-5% central, BBCS-FVTXN correlation

$C(\Delta\phi)$ vs $\Delta\phi$, 3He+Au at 200 GeV, 0-5% central, BBCS-FVTXN correlation

$C(\Delta\phi)$ vs $\Delta\phi$, p+p at 200 GeV, minimum bias, FVTXS-FVTXN correlation

$C(\Delta\phi)$ vs $\Delta\phi$, p+Au at 200 GeV, 0-5% central, FVTXS-FVTXN correlation

$C(\Delta\phi)$ vs $\Delta\phi$, d+Au at 200 GeV, 0-5% central, FVTXS-FVTXN correlation

$C(\Delta\phi)$ vs $\Delta\phi$, 3He+Au at 200 GeV, 0-5% central, FVTXS-FVTXN correlation

$C(\Delta\phi)$ vs $\Delta\phi$, p+p at 200 GeV, minimum bias, CNT-BBCS correlation

$C(\Delta\phi)$ vs $\Delta\phi$, p+Au at 200 GeV, 0-5% central, CNT-BBCS correlation

$C(\Delta\phi)$ vs $\Delta\phi$, d+Au at 200 GeV, 0-5% central, CNT-BBCS correlation

$C(\Delta\phi)$ vs $\Delta\phi$, 3He+Au at 200 GeV, 0-5% central, CNT-BBCS correlation

$C(\Delta\phi)$ vs $\Delta\phi$, p+p at 200 GeV, minimum bias, CNT-FVTXS correlation

$C(\Delta\phi)$ vs $\Delta\phi$, p+Au at 200 GeV, 0-5% central, CNT-FVTXS correlation

$C(\Delta\phi)$ vs $\Delta\phi$, d+Au at 200 GeV, 0-5% central, CNT-FVTXS correlation

$C(\Delta\phi)$ vs $\Delta\phi$, 3He+Au at 200 GeV, 0-5% central, CNT-FVTXS correlation

$C(\Delta\phi)$ vs $\Delta\phi$, p+p at 200 GeV, minimum bias, CNT-FVTXN correlation

$C(\Delta\phi)$ vs $\Delta\phi$, p+Au at 200 GeV, 0-5% central, CNT-FVTXN correlation

$C(\Delta\phi)$ vs $\Delta\phi$, d+Au at 200 GeV, 0-5% central, CNT-FVTXN correlation

$C(\Delta\phi)$ vs $\Delta\phi$, 3He+Au at 200 GeV, 0-5% central, CNT-FVTXN correlation

$c_n$ vs detector combination, p+p at 200 GeV, Minimum Bias

$c_n$ vs $p_T$, p+p at 200 GeV, Minimum Bias, BBCS-CNT detector combination

$c_n$ vs $p_T$, p+p at 200 GeV, Minimum Bias, FVTXS-CNT detector combination

$c_n$ vs $p_T$, p+p at 200 GeV, Minimum Bias, FVTXN-CNT detector combination

$c_n$ vs detector combination, p+Au at 200 GeV, 0-5% central

$c_n$ vs $p_T$, p+Au at 200 GeV, 0-5% central, BBCS-CNT detector combination

$c_n$ vs $p_T$, p+Au at 200 GeV, 0-5% central, FVTXS-CNT detector combination

$c_n$ vs $p_T$, p+Au at 200 GeV, 0-5% central, FVTXN-CNT detector combination

$c_n$ vs detector combination, p+Au at 200 GeV, 0-5% central

$c_n$ vs $p_T$, p+Au at 200 GeV, 0-5% central, BBCS-CNT detector combination

$c_n$ vs $p_T$, p+Au at 200 GeV, 0-5% central, FVTXS-CNT detector combination

$c_n$ vs $p_T$, p+Au at 200 GeV, 0-5% central, FVTXN-CNT detector combination

$c_n$ vs detector combination, p+Au at 200 GeV, 0-5% central

$c_n$ vs $p_T$, 3He+Au at 200 GeV, 0-5% central, BBCS-CNT detector combination

$c_n$ vs $p_T$, 3He+Au at 200 GeV, 0-5% central, FVTXS-CNT detector combination

$c_n$ vs $p_T$, 3He+Au at 200 GeV, 0-5% central, FVTXN-CNT detector combination

The centrality dependencies of the $\Delta\gamma\{\psi_\mathrm{ZDC}\}$ for Au+Au collision at $\sqrt{s_{\rm NN}}$=200 GeV.

The centrality dependencies of the A for Au+Au collision at $\sqrt{s_{\rm NN}}$=200 GeV.

The centrality dependencies of the a for Au+Au collision at $\sqrt{s_{\rm NN}}$=200 GeV.

The centrality dependencies of the A/a using full-event method for Au+Au collision at $\sqrt{s_{\rm NN}}$=200 GeV.

The centrality dependencies of the A/a using sub-event method for Au+Au collision at $\sqrt{s_{\rm NN}}$=200 GeV.

The centrality dependencies of the $f_\mathrm{CME}$ using full-event method for Au+Au collision at $\sqrt{s_{\rm NN}}$=200 GeV.

The centrality dependencies of the $f_\mathrm{CME}$ using sub-event method for Au+Au collision at $\sqrt{s_{\rm NN}}$=200 GeV.

The centrality dependencies of the $\Delta\gamma_\mathrm{CME}$ using full-event method for Au+Au collision at $\sqrt{s_{\rm NN}}$=200 GeV.

The centrality dependencies of the $\Delta\gamma_\mathrm{CME}$ using sub-event method for Au+Au collision at $\sqrt{s_{\rm NN}}$=200 GeV.

The $f_\mathrm{CME}$ from different methods for 20-50% centrality Au+Au collision at $\sqrt{s_{\rm NN}}$=200 GeV.

The $f_\mathrm{CME}$ from different methods for 50-80% centrality Au+Au collision at $\sqrt{s_{\rm NN}}$=200 GeV.

The $\Delta\gamma_\mathrm{CME}$ from different methods for 20-50% centrality Au+Au collision at $\sqrt{s_{\rm NN}}$=200 GeV.

The $\Delta\gamma_\mathrm{CME}$ from different methods for 50-80% centrality Au+Au collision at $\sqrt{s_{\rm NN}}$=200 GeV.

Forward neutron single spin asymmetries as a function of PT

$\pi^0$ spectra from figure 2a from 40-60% U+U collisions. Type A uncertainties are uncorrelated point-to-point. Type B uncertainties are correlated point-to-point. Type C uncertainties affect the scale of the data.

$\pi^0$ spectra from figure 2a from 60-80% U+U collisions. Type A uncertainties are uncorrelated point-to-point. Type B uncertainties are correlated point-to-point. Type C uncertainties affect the scale of the data.

$\eta$ spectra from figure 2b from minimum bias U+U collisions. Type A uncertainties are uncorrelated point-to-point. Type B uncertainties are correlated point-to-point. Type C uncertainties affect the scale of the data.

$\eta$ spectra from figure 2b from 0-20% U+U collisions. Type A uncertainties are uncorrelated point-to-point. Type B uncertainties are correlated point-to-point. Type C uncertainties affect the scale of the data.

$\eta$ spectra from figure 2b from 20-40% U+U collisions. Type A uncertainties are uncorrelated point-to-point. Type B uncertainties are correlated point-to-point. Type C uncertainties affect the scale of the data.

$\eta$ spectra from figure 2b from 40-60% U+U collisions. Type A uncertainties are uncorrelated point-to-point. Type B uncertainties are correlated point-to-point. Type C uncertainties affect the scale of the data.

$\eta$ spectra from figure 2b from 60-80% U+U collisions. Type A uncertainties are uncorrelated point-to-point. Type B uncertainties are correlated point-to-point. Type C uncertainties affect the scale of the data.

Ratio of $\pi^0$ yield to the fit from figure 2c from minimum bias U+U collisions. Sys. uncertainties include B and C uncertainties from data points.

Ratio of $\pi^0$ yield to the fit from figure 2c from 0-20% U+U collisions. Sys. uncertainties include B and C uncertainties from data points.

Ratio of $\pi^0$ yield to the fit from figure 2c from 60-80% U+U collisions. Sys. uncertainties include B and C uncertainties from data points.

Ratio of $\eta$ yield to the fit from figure 2d from minimum bias U+U collisions. Sys. uncertainties include B and C uncertainties from data points.

Ratio of $\eta$ yield to the fit from figure 2d from 0-20% U+U collisions. Sys. uncertainties include B and C uncertainties from data points.

Ratio of $\eta$ yield to the fit from figure 2d from 40-60% U+U collisions. Sys. uncertainties include B and C uncertainties from data points.

$\eta/\pi^0$ ratio from figure 3 from minimum bias U+U collisions. Type A uncertainties are uncorrelated point-to-point. Type B uncertainties are correlated point-to-point.

$\eta/\pi^0$ ratio from figure 3 from 0-20% U+U collisions. Type A uncertainties are uncorrelated point-to-point. Type B uncertainties are correlated point-to-point.

$\eta/\pi^0$ ratio from figure 3 from 20-40% U+U collisions. Type A uncertainties are uncorrelated point-to-point. Type B uncertainties are correlated point-to-point.

$\eta/\pi^0$ ratio from figure 3 from 40-60% U+U collisions. Type A uncertainties are uncorrelated point-to-point. Type B uncertainties are correlated point-to-point.

$\eta/\pi^0$ ratio from figure 3 from 60-80% U+U collisions. Type A uncertainties are uncorrelated point-to-point. Type B uncertainties are correlated point-to-point.

$\eta/\pi^0$ ratio from figure 3 from minimum bias Au+Au collisions. Type A uncertainties are uncorrelated point-to-point. Type B uncertainties are correlated point-to-point.

$\pi^0$ $R_{AA}$ from figure 4a from 0-20% U+U collisions. Type A uncertainties are uncorrelated point-to-point. Type B uncertainties are correlated point-to-point. Type C uncertainties affect the scale of the data. Calculation is carried out for Glauber 1 set.

$\eta$ $R_{AA}$ from figure 4a from 0-20% U+U collisions. Type A uncertainties are uncorrelated point-to-point. Type B uncertainties are correlated point-to-point. Type C uncertainties affect the scale of the data. Calculation is carried out for Glauber 1 set.

$\pi^0$ $R_{AA}$ from figure 4b from 20-40% U+U collisions. Type A uncertainties are uncorrelated point-to-point. Type B uncertainties are correlated point-to-point. Type C uncertainties affect the scale of the data. Calculation is carried out for Glauber 1 set.

$\eta$ $R_{AA}$ from figure 4b from 20-40% U+U collisions. Type A uncertainties are uncorrelated point-to-point. Type B uncertainties are correlated point-to-point. Type C uncertainties affect the scale of the data. Calculation is carried out for Glauber 1 set.

$\pi^0$ $R_{AA}$ from figure 4c from 40-60% U+U collisions. Type A uncertainties are uncorrelated point-to-point. Type B uncertainties are correlated point-to-point. Type C uncertainties affect the scale of the data. Calculation is carried out for Glauber 1 set.

$\eta$ $R_{AA}$ from figure 4c from 40-60% U+U collisions. Type A uncertainties are uncorrelated point-to-point. Type B uncertainties are correlated point-to-point. Type C uncertainties affect the scale of the data. Calculation is carried out for Glauber 1 set.

$\pi^0$ $R_{AA}$ from figure 4d from minimum bias U+U collisions. Type A uncertainties are uncorrelated point-to-point. Type B uncertainties are correlated point-to-point. Type C uncertainties affect the scale of the data. Calculation is carried out for Glauber 1 set.

$\eta$ $R_{AA}$ from figure 4d from minimum bias U+U collisions. Type A uncertainties are uncorrelated point-to-point. Type B uncertainties are correlated point-to-point. Type C uncertainties affect the scale of the data. Calculation is carried out for Glauber 1 set.

$\pi^0$ $R_{AA}$ from figure 5a from 0-20% U+U collisions. Type A uncertainties are uncorrelated point-to-point. Type B uncertainties are correlated point-to-point. Type C uncertainties affect the scale of the data. Calculation is carried out for Glauber 2 set.

$\pi^0$ $R_{AA}$ from figure 5a from 0-5% Au+Au collisions. Type A uncertainties are uncorrelated point-to-point. Type B uncertainties are correlated point-to-point. Type C uncertainties affect the scale of the data.

$\pi^0$ $R_{AA}$ from figure 5a from 40-60% U+U collisions. Type A uncertainties are uncorrelated point-to-point. Type B uncertainties are correlated point-to-point. Type C uncertainties affect the scale of the data. Calculation is carried out for Glauber 2 set.

$\pi^0$ $R_{AA}$ from figure 5a from 40-50% Au+Au collisions. Type A uncertainties are uncorrelated point-to-point. Type B uncertainties are correlated point-to-point. Type C uncertainties affect the scale of the data.

$\pi^0$ integrated $R_{AA}$ for $p_T>5$ GeV/$c$ in U+U collisions as a function of centrality from figure 6a. Type A uncertainties are uncorrelated point-to-point. Type B uncertainties are correlated point-to-point. Type C uncertainties affect the scale of the data. Calculation is carried out for Glauber 1 set.

$\eta$ integrated $R_{AA}$ for $p_T>5$ GeV/$c$ in U+U collisions as a function of centrality from figure 6a. Type A uncertainties are uncorrelated point-to-point. Type B uncertainties are correlated point-to-point. Type C uncertainties affect the scale of the data. Calculation is carried out for Glauber 1 set.

$\pi^0$ integrated $R_{AA}$ for $p_T>5$ GeV/$c$ in Au+Au collisions as a function of centrality from figure 6a. Type A uncertainties are uncorrelated point-to-point. Type B uncertainties are correlated point-to-point. Type C uncertainties affect the scale of the data.

$\eta$ integrated $R_{AA}$ for $p_T>5$ GeV/$c$ in U+U collisions as a function of centrality from figure 6b. Type A uncertainties are uncorrelated point-to-point. Type B uncertainties are correlated point-to-point. Type C uncertainties affect the scale of the data. Calculation is carried out for Glauber 2 set.

Number of participant nucleons of in dep U$+$U collisions

Number of participant nucleons of in dep U$+$U collisions

Number of participant nucleons of in dep Au$+$Au collisions