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

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

The performance of jet energy resolution (JER) for jets with |y| < 2.1 evaluated as a function of pT_truth in different centrality bins. Simulated hard scatter events were overlaid onto events from a dedicated sample of minimum-bias Xe+Xe data. The fit parameters are listed in a sperate table (Extras 1)

The relative magnitude of systematic uncertainties for per-pair normalized xJ distributions in 0-10% Xe+Xe centrality

The relative magnitude of systematic uncertainties for absolutely normalized xJ distributions in 0-10% Xe+Xe centrality

The relative magnitude of systematic uncertainties for rho distributions for leading jets in 0-10% Xe+Xe centrality

Per-pair normalized xJ distribution evaluated in four centrality intervals and given pT1 interval.

Per-pair normalized xJ distribution evaluated in four centrality intervals and given pT1 interval.

Per-pair normalized xJ distribution evaluated in four centrality intervals and given pT1 interval.

Absolutely normalized xJ distribution evaluated in four centrality intervals and given pT1 interval.

Absolutely normalized xJ distribution evaluated in four centrality intervals and given pT1 interval.

Absolutely normalized xJ distribution evaluated in four centrality intervals and given pT1 interval.

Per-pair normalized xJ distribution in Xe+Xe collisions.

Per-pair normalized xJ distribution in Pb+Pb collisions.

Per-pair normalized xJ distribution in Xe+Xe collisions.

Per-pair normalized xJ distribution in Pb+Pb collisions.

Per-pair normalized xJ distribution in Xe+Xe collisions.

Per-pair normalized xJ distribution in Pb+Pb collisions.

Absolutely normalized xJ distribution in Xe+Xe collisions.

Absolutely normalized xJ distribution in Xe+Xe collisions.

Absolutely normalized xJ distribution in Xe+Xe collisions.

Absolutely normalized xJ distribution in Xe+Xe collisions.

Absolutely normalized xJ distribution in Xe+Xe collisions.

Absolutely normalized xJ distribution in Xe+Xe collisions.

The ratios of Xe+Xe and Pb+Pb pair nuclear-modification factors, rho, evaluated as a function of leading jet pT in the same centrality intervals.

The ratios of Xe+Xe and Pb+Pb pair nuclear-modification factors, rho, evaluated as a function of leading jet pT in the same centrality intervals.

The ratios of Xe+Xe and Pb+Pb pair nuclear-modification factors, rho, evaluated as a function of leading jet pT in the same centrality intervals.

The ratios of Xe+Xe and Pb+Pb pair nuclear-modification factors, rho, evaluated as a function of subleading jet pT in the same centrality intervals.

The ratios of Xe+Xe and Pb+Pb pair nuclear-modification factors, rho, evaluated as a function of subleading jet pT in the same centrality intervals.

The ratios of Xe+Xe and Pb+Pb pair nuclear-modification factors, rho, evaluated as a function of subleading jet pT in the same centrality intervals.

The ratios of Xe+Xe and Pb+Pb pair nuclear-modification factors, rho, evaluated as a function of leading jet pT in the same SUM ETFCal intervals (selecting equivalent event activity)

The ratios of Xe+Xe and Pb+Pb pair nuclear-modification factors, rho, evaluated as a function of leading jet pT in the same SUM ETFCal intervals (selecting equivalent event activity)

The ratios of Xe+Xe and Pb+Pb pair nuclear-modification factors, rho, evaluated as a function of leading jet pT in the same SUM ETFCal intervals (selecting equivalent event activity)

The ratios of Xe+Xe and Pb+Pb pair nuclear-modification factors, rho, evaluated as a function of subleading jet pT in the same SUM ETFCal intervals (selecting equivalent event activity)

The ratios of Xe+Xe and Pb+Pb pair nuclear-modification factors, rho, evaluated as a function of subleading jet pT in the same SUM ETFCal intervals (selecting equivalent event activity)

The ratios of Xe+Xe and Pb+Pb pair nuclear-modification factors, rho, evaluated as a function of subleading jet pT in the same SUM ETFCal intervals (selecting equivalent event activity)

Parameter a,b, and c from JER fits in Figure 1b.

Per-pair normalized xJ distribution in Xe+Xe collisions for selected Pb+Pb centrality and pT1 bin.

Per-pair normalized xJ distribution in Xe+Xe collisions for selected Pb+Pb centrality and pT1 bin.

Per-pair normalized xJ distribution in Xe+Xe collisions for selected Pb+Pb centrality and pT1 bin.

Per-pair normalized xJ distribution in Xe+Xe collisions for selected Pb+Pb centrality and pT1 bin.

Per-pair normalized xJ distribution in Xe+Xe collisions for selected Pb+Pb centrality and pT1 bin.

Per-pair normalized xJ distribution in Xe+Xe collisions for selected Pb+Pb centrality and pT1 bin.

Per-pair normalized xJ distribution in Xe+Xe collisions for selected Pb+Pb centrality and pT1 bin.

Per-pair normalized xJ distribution in Xe+Xe collisions for selected Pb+Pb centrality and pT1 bin.

Per-pair normalized xJ distribution in Xe+Xe collisions for selected Pb+Pb centrality and pT1 bin.

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.

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.

The $p_{T}$ dependencies of $v_{1}$ within $-0.4<y<-0.1$ for $t$ in 10-40% mid-central Au+Au collisions at 3 GeV.

The rapidity and $p_{T}$ dependencies of $v_{1}$ for $^{3}$He in 10-40% mid-central Au+Au collisions at 3 GeV.

The $p_{T}$ dependencies of $v_{1}$ within $-0.1<y<0$ for $^4$He in 10-40% mid-central Au+Au collisions at 3 GeV.

The $p_{T}$ dependencies of $v_{1}$ within $-0.4<y<-0.1$ for $^{4}$He in 10-40% mid-central Au+Au collisions at 3 GeV.

The rapidity and $p_{T}$ dependencies of $v_{2}$ for $p$ in 10-40% mid-central Au+Au collisions at 3 GeV.

The rapidity and $p_{T}$ dependencies of $v_{2}$ for $d$ in 10-40% mid-central Au+Au collisions at 3 GeV.

The $p_{T}$ dependencies of $v_{2}$ within $-0.1<y<0$ for $t$ in 10-40% mid-central Au+Au collisions at 3 GeV.

The $p_{T}$ dependencies of $v_{2}$ within $-0.4<y<-0.1$ for $t$ in 10-40% mid-central Au+Au collisions at 3 GeV.

The rapidity and $p_{T}$ dependencies of $v_{2}$ for $^{3}$He in 10-40% mid-central Au+Au collisions at 3 GeV.

The $p_{T}$ dependencies of $v_{2}$ within $-0.1<y<0$ for $^4$He in 10-40% mid-central Au+Au collisions at 3 GeV.

The $p_{T}$ dependencies of $v_{2}$ within $-0.4<y<-0.1$ for $^{4}$He in 10-40% mid-central Au+Au collisions at 3 GeV.

The rapidity dependencies of $v_{1}$ for $p$, $d$, and $^{3}$He in 10-40% mid-central Au+Au collisions at 3 GeV.

The rapidity dependencies of $v_{1}$ for $t$ and $^{4}$He in 10-40% mid-central Au+Au collisions at 3 GeV.

The rapidity dependencies of $v_{1}$ from JAM plus coalescence calculations in 10-40% mid-central Au+Au collisions at 3 GeV.

The rapidity dependencies of $v_{2}$ for $p$, $d$, and $^{3}$He in 10-40% mid-central Au+Au collisions at 3 GeV.

The rapidity dependencies of $v_{2}$ for $t$ and $^{4}$He in 10-40% mid-central Au+Au collisions at 3 GeV.

The rapidity dependencies of $v_{2}$ from JAM plus coalescence calculations in 10-40% mid-central Au+Au collisions at 3 GeV.

Light nucleus scaled $v_{1}$ slopes as a function os collision energy in 10-40 mid-cantral Au+Au collisions.

Light nucleus scaled $v_{1}$ slopes from JAM plus coalescence in 10-40 mid-cantral Au+Au collisions at 3 GeV.

Centrality dependence of the proton cumulant ratios for Au+Au collisions at $\sqrt{s_{NN}}$ = 3.0 GeV. Protons are from $-0.5 < y < 0$ and $0.4 < p_{T} < 2.0$ GeV/$c$. Systematic uncertainties are represented by gray bars. Statistical uncertainties are smaller than marker size. CBWC is applied to all cumulant ratios. While open squares represent the data without the VFC correction, blue triangles and red circles are the results with VFC using the $\langle N_{\rm{part}} \rangle$ distributions from the UrQMD and Glauber models, respectively. UrQMD model results are represented as gold dashed line.

Centrality dependence of the proton cumulant ratios for Au+Au collisions at $\sqrt{s_{NN}}$ = 3.0 GeV. Protons are from $-0.5 < y < 0$ and $0.4 < p_{T} < 2.0$ GeV/$c$. Systematic uncertainties are represented by gray bars. Statistical uncertainties are smaller than marker size. CBWC is applied to all cumulant ratios. While open squares represent the data without the VFC correction, blue triangles and red circles are the results with VFC using the $\langle N_{\rm{part}} \rangle$ distributions from the UrQMD and Glauber models, respectively. UrQMD model results are represented as gold dashed line.

Rapidity and transverse momentum dependence of the proton cumulant ratios for 0–5% central collisions. Black-squares, red-dots and blue-triangles stand for data without and with the VFC using Glauber and UrQMD, respectively.

Rapidity and transverse momentum dependence of the proton cumulant ratios for 0–5% central collisions. Black-squares, red-dots and blue-triangles stand for data without and with the VFC using Glauber and UrQMD, respectively.

Rapidity and transverse momentum dependence of the proton cumulant ratios for 0–5% central collisions. Black-squares, red-dots and blue-triangles stand for data without and with the VFC using Glauber and UrQMD, respectively.

Rapidity and transverse momentum dependence of the proton cumulant ratios for 0–5% central collisions. Black-squares, red-dots and blue-triangles stand for data without and with the VFC using Glauber and UrQMD, respectively.

Collision energy dependence of the ratios of cumulants, $C_4/C_2$, for proton (squares) and net-proton (red circles) from top 0-5% Au+Au collisions at RHIC. The points for protons are shifted horizontally for clarity. The new result for proton from $\sqrt{s_{NN}}$ = 3.0 GeV collisions is shown as a filled square. HADES data of $\sqrt{s_{NN}}$ = 2.4 GeV 0-10% collisions is also shown. The vertical black and gray bars are the statistical and systematic uncertainties, respectively. In addition, results from the HRG model, based on both Canonical Ensemble (CE) and Grand-Canonical Ensemble (GCE), and transport model UrQMD are presented.

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.

The correlation functions (corrected for nonuniform detector efficiency in $\phi$; not corrected for the absolute detection efficiency) vs. azimuthal angle difference between forward ($2.6<\eta<4.0$) $\pi^{0}$s in $p$+$p$ collisions at $\sqrt{s_{\mathrm{_{NN}}}}=200$ GeV at high $p_{T}$ ($p^{trig}_{T}$=2.5-3 GeV/c, $p^{asso}_{T}$=2-2.5 GeV/c)

The correlation functions (corrected for nonuniform detector efficiency in $\phi$; not corrected for the absolute detection efficiency) vs. azimuthal angle difference between forward ($2.6<\eta<4.0$) $\pi^{0}$s in $p$+Al collisions at $\sqrt{s_{\mathrm{_{NN}}}}=200$ GeV at high $p_{T}$ ($p^{trig}_{T}$=2.5-3 GeV/c, $p^{asso}_{T}$=2-2.5 GeV/c)

The correlation functions (corrected for nonuniform detector efficiency in $\phi$; not corrected for the absolute detection efficiency) vs. azimuthal angle difference between forward ($2.6<\eta<4.0$) $\pi^{0}$s in $p$+$Au collisions at $\sqrt{s_{\mathrm{_{NN}}}}=200$ GeV at high $p_{T}$ ($p^{trig}_{T}$=2.5-3 GeV/c, $p^{asso}_{T}$=2-2.5 GeV/c)

Relative area of MinBias $pAu/pp$ of back-to-back di-$\pi^0$ correlations at forward pseudorapidities ($2.6<\eta<4.0$) from data

Relative area of MinBias $pAl/pp$ of back-to-back di-$\pi^0$ correlations at forward pseudorapidities ($2.6<\eta<4.0$) from data

Relative area of $pAu/pp$ at b=0 of back-to-back di-$\pi^0$ correlations at forward pseudorapidities ($2.6<\eta<4.0$) from prediction

Relative width of MinBias $pAu/pp$ of back-to-back di-$\pi^0$ correlations at forward pseudorapidities ($2.6<\eta<4.0$) from data

Relative width of MinBias $pAl/pp$ of back-to-back di-$\pi^0$ correlations at forward pseudorapidities ($2.6<\eta<4.0$) from data

Relative pedestal of MinBias $pAu/pp$ of back-to-back di-$\pi^0$ correlations at forward pseudorapidities ($2.6<\eta<4.0$) from data

Relative pedestal of MinBias $pAl/pp$ of back-to-back di-$\pi^0$ correlations at forward pseudorapidities ($2.6<\eta<4.0$) from data

Relative area of MinBias $pp/pp$ of back-to-back di-$\pi^0$ correlations at forward pseudorapidities ($2.6<\eta<4.0$), the ratio is unity, error is zero

Relative area of MinBias $pA/pp$ of back-to-back di-$\pi^0$ correlations at forward pseudorapidities ($2.6<\eta<4.0$) from data

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.