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

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.

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.

$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

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.

Event by event net-proton multiplicity, $\Delta N_{p}$, distributions for Au+Au collisions at √sNN = 27, 54.4, and 200 GeV in 0-10% and 30-40% centralities at midrapidity (|y| < 0.5) for the transverse momentum range of 0.4 < $p_{T}$ (GeV/c) < 2.0. These distributions are normalized by the corresponding numbers of events and are not corrected for detector efficiencies. Statistical uncertainties are shown as vertical lines. The dashed lines show the Skellam distributions for each collision energy and centrality. The bottom panel shows the ratio of the data to the Skellam expectations.

Event by event net-proton multiplicity, $\Delta N_{p}$, distributions for Au+Au collisions at √sNN = 27, 54.4, and 200 GeV in 0-10% and 30-40% centralities at midrapidity (|y| < 0.5) for the transverse momentum range of 0.4 < $p_{T}$ (GeV/c) < 2.0. These distributions are normalized by the corresponding numbers of events and are not corrected for detector efficiencies. Statistical uncertainties are shown as vertical lines. The dashed lines show the Skellam distributions for each collision energy and centrality. The bottom panel shows the ratio of the data to the Skellam expectations.

Event by event net-proton multiplicity, $\Delta N_{p}$, distributions for Au+Au collisions at √sNN = 27, 54.4, and 200 GeV in 0-10% and 30-40% centralities at midrapidity (|y| < 0.5) for the transverse momentum range of 0.4 < $p_{T}$ (GeV/c) < 2.0. These distributions are normalized by the corresponding numbers of events and are not corrected for detector efficiencies. Statistical uncertainties are shown as vertical lines. The dashed lines show the Skellam distributions for each collision energy and centrality. The bottom panel shows the ratio of the data to the Skellam expectations.

Net-proton $C_{6}/C_{2}$ as a function of rapidity (left) and transverse momentum acceptance (right) from $\sqrt{s_{NN}}$ = 27 GeV (crosses), 54.4 (open squares), and 200 GeV (filled circles) Au+Au collisions. The upper and lower plots are for 0-10% and 30-40% centralities, respectively. The error bars are statistical and caps are systematic errors. Points for different beam energies are staggered horizontally to improve clarity. UrQMD transport model results are shown as shaded and hatched bands. The Skellam expectation ($C_{6}/C_{2} = 1) is shown as long-dashed lines.

Net-proton $C_{6}/C_{2}$ as a function of rapidity (left) and transverse momentum acceptance (right) from $\sqrt{s_{NN}}$ = 27 GeV (crosses), 54.4 (open squares), and 200 GeV (filled circles) Au+Au collisions. The upper and lower plots are for 0-10% and 30-40% centralities, respectively. The error bars are statistical and caps are systematic errors. Points for different beam energies are staggered horizontally to improve clarity. UrQMD transport model results are shown as shaded and hatched bands. The Skellam expectation ($C_{6}/C_{2} = 1) is shown as long-dashed lines.

Net-proton $C_{6}/C_{2}$ as a function of rapidity (left) and transverse momentum acceptance (right) from $\sqrt{s_{NN}}$ = 27 GeV (crosses), 54.4 (open squares), and 200 GeV (filled circles) Au+Au collisions. The upper and lower plots are for 0-10% and 30-40% centralities, respectively. The error bars are statistical and caps are systematic errors. Points for different beam energies are staggered horizontally to improve clarity. UrQMD transport model results are shown as shaded and hatched bands. The Skellam expectation ($C_{6}/C_{2} = 1) is shown as long-dashed lines.

Net-proton $C_{6}/C_{2}$ as a function of rapidity (left) and transverse momentum acceptance (right) from $\sqrt{s_{NN}}$ = 27 GeV (crosses), 54.4 (open squares), and 200 GeV (filled circles) Au+Au collisions. The upper and lower plots are for 0-10% and 30-40% centralities, respectively. The error bars are statistical and caps are systematic errors. Points for different beam energies are staggered horizontally to improve clarity. UrQMD transport model results are shown as shaded and hatched bands. The Skellam expectation ($C_{6}/C_{2} = 1) is shown as long-dashed lines.

Net-proton $C_{6}/C_{2}$ as a function of rapidity (left) and transverse momentum acceptance (right) from $\sqrt{s_{NN}}$ = 27 GeV (crosses), 54.4 (open squares), and 200 GeV (filled circles) Au+Au collisions. The upper and lower plots are for 0-10% and 30-40% centralities, respectively. The error bars are statistical and caps are systematic errors. Points for different beam energies are staggered horizontally to improve clarity. UrQMD transport model results are shown as shaded and hatched bands. The Skellam expectation ($C_{6}/C_{2} = 1) is shown as long-dashed lines.

Net-proton $C_{6}/C_{2}$ as a function of rapidity (left) and transverse momentum acceptance (right) from $\sqrt{s_{NN}}$ = 27 GeV (crosses), 54.4 (open squares), and 200 GeV (filled circles) Au+Au collisions. The upper and lower plots are for 0-10% and 30-40% centralities, respectively. The error bars are statistical and caps are systematic errors. Points for different beam energies are staggered horizontally to improve clarity. UrQMD transport model results are shown as shaded and hatched bands. The Skellam expectation ($C_{6}/C_{2} = 1) is shown as long-dashed lines.

Collisions centrality dependence of net-proton $C_{6}/C_{2}$ in Au+Au collisions for |$y$| < 0.5 and 0.4 < $p_{T}$ (GeV/c) < 2.0. The error bars are statistical and caps are systematic errors. Points for different beam energies are staggered horizontally to improve clarity. A shaded band shows the results from UrQMD model calculations. UrQMD calculations from the above three collision energies are consistent among them so they are merged in order to reduce statistical fluctuations. Details on these calculations can be found in the Supplemental Material at [URL will be inserted by publisher]. The lattice QCD calculations [16, 17] for T = 160 MeV and $\mu_{B}$ < 110 MeV. are shown as a blue band at $\langle N_{part}\rangle$ $\approx$ 340.

Collisions centrality dependence of net-proton $C_{6}/C_{2}$ in Au+Au collisions for |$y$| < 0.5 and 0.4 < $p_{T}$ (GeV/c) < 2.0. The error bars are statistical and caps are systematic errors. Points for different beam energies are staggered horizontally to improve clarity. A shaded band shows the results from UrQMD model calculations. UrQMD calculations from the above three collision energies are consistent among them so they are merged in order to reduce statistical fluctuations. Details on these calculations can be found in the Supplemental Material at [URL will be inserted by publisher]. The lattice QCD calculations [16, 17] for T = 160 MeV and $\mu_{B}$ < 110 MeV. are shown as a blue band at $\langle N_{part}\rangle$ $\approx$ 340.

Collisions centrality dependence of net-proton $C_{6}/C_{2}$ in Au+Au collisions for |$y$| < 0.5 and 0.4 < $p_{T}$ (GeV/c) < 2.0. The error bars are statistical and caps are systematic errors. Points for different beam energies are staggered horizontally to improve clarity. A shaded band shows the results from UrQMD model calculations. UrQMD calculations from the above three collision energies are consistent among them so they are merged in order to reduce statistical fluctuations. Details on these calculations can be found in the Supplemental Material at [URL will be inserted by publisher]. The lattice QCD calculations [16, 17] for T = 160 MeV and $\mu_{B}$ < 110 MeV. are shown as a blue band at $\langle N_{part}\rangle$ $\approx$ 340.

Distributions of reconstructed protons (black circles) from embedding simulations in 200 GeV Au+Au collisions at 0-2.5% centrality. Red lines are fits with the binomial distribution, and green dotted lines represent the fit with the beta-binomial distributions using α that gives the minimal chi2/ndf. Each panel shows the result from the given combinations of embedded protons and antiprotons. The ratio of fits to the embedding data is shown for each panel.

Distributions of reconstructed protons (black circles) from embedding simulations in 200 GeV Au+Au collisions at 0-2.5% centrality. Red lines are fits with the binomial distribution, and green dotted lines represent the fit with the beta-binomial distributions using α that gives the minimal chi2/ndf. Each panel shows the result from the given combinations of embedded protons and antiprotons. The ratio of fits to the embedding data is shown for each panel.

Distributions of reconstructed protons (black circles) from embedding simulations in 200 GeV Au+Au collisions at 0-2.5% centrality. Red lines are fits with the binomial distribution, and green dotted lines represent the fit with the beta-binomial distributions using α that gives the minimal chi2/ndf. Each panel shows the result from the given combinations of embedded protons and antiprotons. The ratio of fits to the embedding data is shown for each panel.

Distributions of reconstructed protons (black circles) from embedding simulations in 200 GeV Au+Au collisions at 0-2.5% centrality. Red lines are fits with the binomial distribution, and green dotted lines represent the fit with the beta-binomial distributions using α that gives the minimal chi2/ndf. Each panel shows the result from the given combinations of embedded protons and antiprotons. The ratio of fits to the embedding data is shown for each panel.

Distributions of reconstructed protons (black circles) from embedding simulations in 200 GeV Au+Au collisions at 0-2.5% centrality. Red lines are fits with the binomial distribution, and green dotted lines represent the fit with the beta-binomial distributions using α that gives the minimal chi2/ndf. Each panel shows the result from the given combinations of embedded protons and antiprotons. The ratio of fits to the embedding data is shown for each panel.

Distributions of reconstructed protons (black circles) from embedding simulations in 200 GeV Au+Au collisions at 0-2.5% centrality. Red lines are fits with the binomial distribution, and green dotted lines represent the fit with the beta-binomial distributions using α that gives the minimal chi2/ndf. Each panel shows the result from the given combinations of embedded protons and antiprotons. The ratio of fits to the embedding data is shown for each panel.

Distributions of reconstructed protons (black circles) from embedding simulations in 200 GeV Au+Au collisions at 0-2.5% centrality. Red lines are fits with the binomial distribution, and green dotted lines represent the fit with the beta-binomial distributions using α that gives the minimal chi2/ndf. Each panel shows the result from the given combinations of embedded protons and antiprotons. The ratio of fits to the embedding data is shown for each panel.

Distributions of reconstructed protons (black circles) from embedding simulations in 200 GeV Au+Au collisions at 0-2.5% centrality. Red lines are fits with the binomial distribution, and green dotted lines represent the fit with the beta-binomial distributions using α that gives the minimal chi2/ndf. Each panel shows the result from the given combinations of embedded protons and antiprotons. The ratio of fits to the embedding data is shown for each panel.

Distributions of reconstructed protons (black circles) from embedding simulations in 200 GeV Au+Au collisions at 0-2.5% centrality. Red lines are fits with the binomial distribution, and green dotted lines represent the fit with the beta-binomial distributions using α that gives the minimal chi2/ndf. Each panel shows the result from the given combinations of embedded protons and antiprotons. The ratio of fits to the embedding data is shown for each panel.

Distributions of reconstructed protons (black circles) from embedding simulations in 200 GeV Au+Au collisions at 0-2.5% centrality. Red lines are fits with the binomial distribution, and green dotted lines represent the fit with the beta-binomial distributions using α that gives the minimal chi2/ndf. Each panel shows the result from the given combinations of embedded protons and antiprotons. The ratio of fits to the embedding data is shown for each panel.

Distributions of reconstructed protons (black circles) from embedding simulations in 200 GeV Au+Au collisions at 0-2.5% centrality. Red lines are fits with the binomial distribution, and green dotted lines represent the fit with the beta-binomial distributions using α that gives the minimal chi2/ndf. Each panel shows the result from the given combinations of embedded protons and antiprotons. The ratio of fits to the embedding data is shown for each panel.

Distributions of reconstructed protons (black circles) from embedding simulations in 200 GeV Au+Au collisions at 0-2.5% centrality. Red lines are fits with the binomial distribution, and green dotted lines represent the fit with the beta-binomial distributions using α that gives the minimal chi2/ndf. Each panel shows the result from the given combinations of embedded protons and antiprotons. The ratio of fits to the embedding data is shown for each panel.

Distributions of reconstructed protons (black circles) from embedding simulations in 200 GeV Au+Au collisions at 0-2.5% centrality. Red lines are fits with the binomial distribution, and green dotted lines represent the fit with the beta-binomial distributions using α that gives the minimal chi2/ndf. Each panel shows the result from the given combinations of embedded protons and antiprotons. The ratio of fits to the embedding data is shown for each panel.

Distributions of reconstructed protons (black circles) from embedding simulations in 200 GeV Au+Au collisions at 0-2.5% centrality. Red lines are fits with the binomial distribution, and green dotted lines represent the fit with the beta-binomial distributions using α that gives the minimal chi2/ndf. Each panel shows the result from the given combinations of embedded protons and antiprotons. The ratio of fits to the embedding data is shown for each panel.

Distributions of reconstructed protons (black circles) from embedding simulations in 200 GeV Au+Au collisions at 0-2.5% centrality. Red lines are fits with the binomial distribution, and green dotted lines represent the fit with the beta-binomial distributions using α that gives the minimal chi2/ndf. Each panel shows the result from the given combinations of embedded protons and antiprotons. The ratio of fits to the embedding data is shown for each panel.

Distributions of reconstructed protons (black circles) from embedding simulations in 200 GeV Au+Au collisions at 0-2.5% centrality. Red lines are fits with the binomial distribution, and green dotted lines represent the fit with the beta-binomial distributions using α that gives the minimal chi2/ndf. Each panel shows the result from the given combinations of embedded protons and antiprotons. The ratio of fits to the embedding data is shown for each panel.

Distributions of reconstructed protons (black circles) from embedding simulations in 200 GeV Au+Au collisions at 0-2.5% centrality. Red lines are fits with the binomial distribution, and green dotted lines represent the fit with the beta-binomial distributions using α that gives the minimal chi2/ndf. Each panel shows the result from the given combinations of embedded protons and antiprotons. The ratio of fits to the embedding data is shown for each panel.

Distributions of reconstructed protons (black circles) from embedding simulations in 200 GeV Au+Au collisions at 0-2.5% centrality. Red lines are fits with the binomial distribution, and green dotted lines represent the fit with the beta-binomial distributions using α that gives the minimal chi2/ndf. Each panel shows the result from the given combinations of embedded protons and antiprotons. The ratio of fits to the embedding data is shown for each panel.

Distributions of reconstructed protons (black circles) from embedding simulations in 200 GeV Au+Au collisions at 0-2.5% centrality. Red lines are fits with the binomial distribution, and green dotted lines represent the fit with the beta-binomial distributions using α that gives the minimal chi2/ndf. Each panel shows the result from the given combinations of embedded protons and antiprotons. The ratio of fits to the embedding data is shown for each panel.

Distributions of reconstructed protons (black circles) from embedding simulations in 200 GeV Au+Au collisions at 0-2.5% centrality. Red lines are fits with the binomial distribution, and green dotted lines represent the fit with the beta-binomial distributions using α that gives the minimal chi2/ndf. Each panel shows the result from the given combinations of embedded protons and antiprotons. The ratio of fits to the embedding data is shown for each panel.

Distributions of reconstructed protons (black circles) from embedding simulations in 200 GeV Au+Au collisions at 0-2.5% centrality. Red lines are fits with the binomial distribution, and green dotted lines represent the fit with the beta-binomial distributions using α that gives the minimal chi2/ndf. Each panel shows the result from the given combinations of embedded protons and antiprotons. The ratio of fits to the embedding data is shown for each panel.

Distributions of reconstructed protons (black circles) from embedding simulations in 200 GeV Au+Au collisions at 0-2.5% centrality. Red lines are fits with the binomial distribution, and green dotted lines represent the fit with the beta-binomial distributions using α that gives the minimal chi2/ndf. Each panel shows the result from the given combinations of embedded protons and antiprotons. The ratio of fits to the embedding data is shown for each panel.

Distributions of reconstructed protons (black circles) from embedding simulations in 200 GeV Au+Au collisions at 0-2.5% centrality. Red lines are fits with the binomial distribution, and green dotted lines represent the fit with the beta-binomial distributions using α that gives the minimal chi2/ndf. Each panel shows the result from the given combinations of embedded protons and antiprotons. The ratio of fits to the embedding data is shown for each panel.

Distributions of reconstructed protons (black circles) from embedding simulations in 200 GeV Au+Au collisions at 0-2.5% centrality. Red lines are fits with the binomial distribution, and green dotted lines represent the fit with the beta-binomial distributions using α that gives the minimal chi2/ndf. Each panel shows the result from the given combinations of embedded protons and antiprotons. The ratio of fits to the embedding data is shown for each panel.

Distributions of reconstructed protons (black circles) from embedding simulations in 200 GeV Au+Au collisions at 0-2.5% centrality. Red lines are fits with the binomial distribution, and green dotted lines represent the fit with the beta-binomial distributions using α that gives the minimal chi2/ndf. Each panel shows the result from the given combinations of embedded protons and antiprotons. The ratio of fits to the embedding data is shown for each panel.

Distributions of reconstructed protons (black circles) from embedding simulations in 200 GeV Au+Au collisions at 0-2.5% centrality. Red lines are fits with the binomial distribution, and green dotted lines represent the fit with the beta-binomial distributions using α that gives the minimal chi2/ndf. Each panel shows the result from the given combinations of embedded protons and antiprotons. The ratio of fits to the embedding data is shown for each panel.

Distributions of reconstructed protons (black circles) from embedding simulations in 200 GeV Au+Au collisions at 0-2.5% centrality. Red lines are fits with the binomial distribution, and green dotted lines represent the fit with the beta-binomial distributions using α that gives the minimal chi2/ndf. Each panel shows the result from the given combinations of embedded protons and antiprotons. The ratio of fits to the embedding data is shown for each panel.

Distributions of reconstructed protons (black circles) from embedding simulations in 200 GeV Au+Au collisions at 0-2.5% centrality. Red lines are fits with the binomial distribution, and green dotted lines represent the fit with the beta-binomial distributions using α that gives the minimal chi2/ndf. Each panel shows the result from the given combinations of embedded protons and antiprotons. The ratio of fits to the embedding data is shown for each panel.

Distributions of reconstructed protons (black circles) from embedding simulations in 200 GeV Au+Au collisions at 0-2.5% centrality. Red lines are fits with the binomial distribution, and green dotted lines represent the fit with the beta-binomial distributions using α that gives the minimal chi2/ndf. Each panel shows the result from the given combinations of embedded protons and antiprotons. The ratio of fits to the embedding data is shown for each panel.

Distributions of reconstructed protons (black circles) from embedding simulations in 200 GeV Au+Au collisions at 0-2.5% centrality. Red lines are fits with the binomial distribution, and green dotted lines represent the fit with the beta-binomial distributions using α that gives the minimal chi2/ndf. Each panel shows the result from the given combinations of embedded protons and antiprotons. The ratio of fits to the embedding data is shown for each panel.

Distributions of reconstructed protons (black circles) from embedding simulations in 200 GeV Au+Au collisions at 0-2.5% centrality. Red lines are fits with the binomial distribution, and green dotted lines represent the fit with the beta-binomial distributions using α that gives the minimal chi2/ndf. Each panel shows the result from the given combinations of embedded protons and antiprotons. The ratio of fits to the embedding data is shown for each panel.

Distributions of reconstructed protons (black circles) from embedding simulations in 200 GeV Au+Au collisions at 0-2.5% centrality. Red lines are fits with the binomial distribution, and green dotted lines represent the fit with the beta-binomial distributions using α that gives the minimal chi2/ndf. Each panel shows the result from the given combinations of embedded protons and antiprotons. The ratio of fits to the embedding data is shown for each panel.

Distributions of reconstructed protons (black circles) from embedding simulations in 200 GeV Au+Au collisions at 0-2.5% centrality. Red lines are fits with the binomial distribution, and green dotted lines represent the fit with the beta-binomial distributions using α that gives the minimal chi2/ndf. Each panel shows the result from the given combinations of embedded protons and antiprotons. The ratio of fits to the embedding data is shown for each panel.

Cumulants and their ratios up to the sixth order corrected for non-binomial efficiencies for 200 GeV Au+Au collisions at 0-5% centrality. The CBWC is applied for 2.5% centrality bin width. Results from the conventional efficiency correction are shown as black filled circles, results from the unfolding with the binomial detector response are shown as black open circles, and results from beta-binomial detector response with $\alpha+\sigma$, $\alpha$ and $\alpha-\sigma$ are shown in green triangles, red squares and blue triangles, respectively. C5, C6, C2/C1, C5/C1 and C6/C2 are scaled by constant shown in each column.

Collision centrality dependence of net-proton C6/C2 in Au+Au collisions for $\sqrt{s_{NN}}$ = 200 GeV within |y| < 0.5 and 0.4 < pT (GeV/c) < 2.0. Results with and without the CBWC are overlaid. The results are corrected for detector efficiencies. Points for different calculation methods are staggered horizontally to improve clarity.

Collision centrality dependence of net-proton C6/C2 in Au+Au collisions for $\sqrt{s_{NN}}$ = 27, 54.4, and 200 GeV within |y| < 0.5 and 0.4 < pT (GeV/c) < 2.0. Points for different beam energies are staggered horizontally to improve clarity. Shaded and hatched bands show the results from UrQMD model calculations The lattice QCD calculations [13, 14] for T = 160 MeV and $\mu_{B}$ < 110 MeV. are shown as a blue band at $\langle N_{part}\rangle$ $\approx$ 340.

The uncorrected SoftDrop groomed jet mass distribution for $R = 0.4$ anti-$k_{\rm{T}}$ jets with $20 < p_{\rm{T,jet}} < 25$ GeV$/c$. Updated to correct a small bug that had shifted the jet mass to slightly smaller values.

The fully corrected jet mass distribution for $R = 0.4$ anti-$k_{\rm{T}}$ jets with $20 < p_{\rm{T,jet}} < 25$ GeV$/c$. Data are reported beyond the x-axis upper limit of the figure, for future reference.

The fully corrected jet mass distribution for $R = 0.4$ anti-$k_{\rm{T}}$ jets with $20 < p_{\rm{T,jet}} < 25$ GeV$/c$. Data are reported beyond the x-axis upper limit of the figure, for future reference. Updated to correct a small bug that had shifted the jet mass to slightly smaller values.

The fully corrected jet mass distribution for $R = 0.4$ anti-$k_{\rm{T}}$ jets with $25 < p_{\rm{T,jet}} < 30$ GeV$/c$. Data are reported beyond the x-axis upper limit of the figure, for future reference.

The fully corrected jet mass distribution for $R = 0.4$ anti-$k_{\rm{T}}$ jets with $25 < p_{\rm{T,jet}} < 30$ GeV$/c$. Data are reported beyond the x-axis upper limit of the figure, for future reference. Updated to correct a small bug that had shifted the jet mass to slightly smaller values.

The fully corrected jet mass distribution for $R = 0.4$ anti-$k_{\rm{T}}$ jets with $30 < p_{\rm{T,jet}} < 40$ GeV$/c$. Data are reported beyond the x-axis upper limit of the figure, for future reference.

The fully corrected jet mass distribution for $R = 0.4$ anti-$k_{\rm{T}}$ jets with $30 < p_{\rm{T,jet}} < 40$ GeV$/c$. Data are reported beyond the x-axis upper limit of the figure, for future reference. Updated to correct a small bug that had shifted the jet mass to slightly smaller values.

The fully corrected SoftDrop groomed jet mass distribution for $R = 0.4$ anti-$k_{\rm{T}}$ jets with $20 < p_{\rm{T,jet}} < 25$ GeV$/c$. Data are reported beyond the x-axis upper limit of the figure, for future reference.

The fully corrected SoftDrop groomed jet mass distribution for $R = 0.4$ anti-$k_{\rm{T}}$ jets with $20 < p_{\rm{T,jet}} < 25$ GeV$/c$. Data are reported beyond the x-axis upper limit of the figure, for future reference. Updated to correct a small bug that had shifted the jet mass to slightly smaller values.

The fully corrected SoftDrop groomed jet mass distribution for $R = 0.4$ anti-$k_{\rm{T}}$ jets with $25 < p_{\rm{T,jet}} < 30$ GeV$/c$. Data are reported beyond the x-axis upper limit of the figure, for future reference.

The fully corrected SoftDrop groomed jet mass distribution for $R = 0.4$ anti-$k_{\rm{T}}$ jets with $25 < p_{\rm{T,jet}} < 30$ GeV$/c$. Data are reported beyond the x-axis upper limit of the figure, for future reference. Updated to correct a small bug that had shifted the jet mass to slightly smaller values.

The fully corrected SoftDrop groomed jet mass distribution for $R = 0.4$ anti-$k_{\rm{T}}$ jets with $30 < p_{\rm{T,jet}} < 40$ GeV$/c$. Data are reported beyond the x-axis upper limit of the figure, for future reference.

The fully corrected SoftDrop groomed jet mass distribution for $R = 0.4$ anti-$k_{\rm{T}}$ jets with $30 < p_{\rm{T,jet}} < 40$ GeV$/c$. Data are reported beyond the x-axis upper limit of the figure, for future reference. Updated to correct a small bug that had shifted the jet mass to slightly smaller values.

The fully corrected jet mass distribution for $R = 0.2$ anti-$k_{\rm{T}}$ jets with $30 < p_{\rm{T,jet}} < 40$ GeV$/c$.

The fully corrected jet mass distribution for $R = 0.2$ anti-$k_{\rm{T}}$ jets with $30 < p_{\rm{T,jet}} < 40$ GeV$/c$. Updated to correct a small bug that had shifted the jet mass to slightly smaller values.

The fully corrected jet mass distribution for $R = 0.4$ anti-$k_{\rm{T}}$ jets with $30 < p_{\rm{T,jet}} < 40$ GeV$/c$.

The fully corrected jet mass distribution for $R = 0.4$ anti-$k_{\rm{T}}$ jets with $30 < p_{\rm{T,jet}} < 40$ GeV$/c$. Updated to correct a small bug that had shifted the jet mass to slightly smaller values.

The fully corrected jet mass distribution for $R = 0.6$ anti-$k_{\rm{T}}$ jets with $30 < p_{\rm{T,jet}} < 40$ GeV$/c$.

The fully corrected jet mass distribution for $R = 0.6$ anti-$k_{\rm{T}}$ jets with $30 < p_{\rm{T,jet}} < 40$ GeV$/c$. Updated to correct a small bug that had shifted the jet mass to slightly smaller values.

The fully corrected SoftDrop groomed jet mass distribution for $R = 0.2$ anti-$k_{\rm{T}}$ jets with $30 < p_{\rm{T,jet}} < 40$ GeV$/c$.

The fully corrected SoftDrop groomed jet mass distribution for $R = 0.2$ anti-$k_{\rm{T}}$ jets with $30 < p_{\rm{T,jet}} < 40$ GeV$/c$. Updated to correct a small bug that had shifted the jet mass to slightly smaller values.

The fully corrected SoftDrop groomed jet mass distribution for $R = 0.4$ anti-$k_{\rm{T}}$ jets with $30 < p_{\rm{T,jet}} < 40$ GeV$/c$.

The fully corrected SoftDrop groomed jet mass distribution for $R = 0.4$ anti-$k_{\rm{T}}$ jets with $30 < p_{\rm{T,jet}} < 40$ GeV$/c$. Updated to correct a small bug that had shifted the jet mass to slightly smaller values.

The fully corrected SoftDrop groomed jet mass distribution for $R = 0.6$ anti-$k_{\rm{T}}$ jets with $30 < p_{\rm{T,jet}} < 40$ GeV$/c$.

The fully corrected SoftDrop groomed jet mass distribution for $R = 0.6$ anti-$k_{\rm{T}}$ jets with $30 < p_{\rm{T,jet}} < 40$ GeV$/c$. Updated to correct a small bug that had shifted the jet mass to slightly smaller values.

The mean of the fully corrected jet mass distribution for $R=0.2$ anti-$k_{\rm{T}}$ jets as a function of $p_{\rm{T,jet}}$.

The mean of the fully corrected jet mass distribution for $R=0.2$ anti-$k_{\rm{T}}$ jets as a function of $p_{\rm{T,jet}}$. Updated to correct a small bug that had shifted the jet mass to slightly smaller values.

The mean of the fully corrected jet mass distribution for $R=0.4$ anti-$k_{\rm{T}}$ jets as a function of $p_{\rm{T,jet}}$.

The mean of the fully corrected jet mass distribution for $R=0.4$ anti-$k_{\rm{T}}$ jets as a function of $p_{\rm{T,jet}}$. Updated to correct a small bug that had shifted the jet mass to slightly smaller values.

The mean of the fully corrected jet mass distribution for $R=0.6$ anti-$k_{\rm{T}}$ jets as a function of $p_{\rm{T,jet}}$.

The mean of the fully corrected jet mass distribution for $R=0.6$ anti-$k_{\rm{T}}$ jets as a function of $p_{\rm{T,jet}}$. Updated to correct a small bug that had shifted the jet mass to slightly smaller values.

The mean of the fully corrected SoftDrop groomed jet mass distribution for $R=0.2$ anti-$k_{\rm{T}}$ jets as a function of $p_{\rm{T,jet}}$.

The mean of the fully corrected SoftDrop groomed jet mass distribution for $R=0.2$ anti-$k_{\rm{T}}$ jets as a function of $p_{\rm{T,jet}}$. Updated to correct a small bug that had shifted the jet mass to slightly smaller values.

The mean of the fully corrected SoftDrop groomed jet mass distribution for $R=0.4$ anti-$k_{\rm{T}}$ jets as a function of $p_{\rm{T,jet}}$.

The mean of the fully corrected SoftDrop groomed jet mass distribution for $R=0.4$ anti-$k_{\rm{T}}$ jets as a function of $p_{\rm{T,jet}}$. Updated to correct a small bug that had shifted the jet mass to slightly smaller values.

The mean of the fully corrected SoftDrop groomed jet mass distribution for $R=0.6$ anti-$k_{\rm{T}}$ jets as a function of $p_{\rm{T,jet}}$.

The mean of the fully corrected SoftDrop groomed jet mass distribution for $R=0.6$ anti-$k_{\rm{T}}$ jets as a function of $p_{\rm{T,jet}}$. Updated to correct a small bug that had shifted the jet mass to slightly smaller values.

The fully corrected jet mass distribution for $R = 0.6$ anti-$k_{\rm{T}}$ jets with $20 < p_{\rm{T,jet}} < 25$ GeV$/c$. Data are reported beyond the x-axis upper limit of the figure, for future reference.

The fully corrected jet mass distribution for $R = 0.6$ anti-$k_{\rm{T}}$ jets with $20 < p_{\rm{T,jet}} < 25$ GeV$/c$. Data are reported beyond the x-axis upper limit of the figure, for future reference. Updated to correct a small bug that had shifted the jet mass to slightly smaller values.

The fully corrected SoftDrop groomed jet mass distribution for $R = 0.6$ anti-$k_{\rm{T}}$ jets with $20 < p_{\rm{T,jet}} < 25$ GeV$/c$. Data are reported beyond the x-axis upper limit of the figure, for future reference.

The fully corrected SoftDrop groomed jet mass distribution for $R = 0.6$ anti-$k_{\rm{T}}$ jets with $20 < p_{\rm{T,jet}} < 25$ GeV$/c$. Data are reported beyond the x-axis upper limit of the figure, for future reference. Updated to correct a small bug that had shifted the jet mass to slightly smaller values.