Rapidity($y$) dependence of $v_1$ (top panels) and $v_2$ (bottom panels) of proton and $\Lambda$ baryons (left panels), pions (middle panels) and kaons (right panels) in 10-40% centrality for the $\sqrt{s_{NN}}$ = 3GeV Au+Au collisions. Statistical and systematic uncertainties are shown as bars and gray bands, respectively. Some uncertainties are smaller than the data points. The UrQMD and JAM results are shown as bands:golden, red and blue bands stand for JAM mean-field, UrQMD mean-field and UrQMD cascade mode, respectively. The value of the incompressibility $\kappa$ = 380 MeV is used in the mean-field option. More detailed model descriptions and data comparisons can be found in Supplemental Material.

Rapidity($y$) dependence of $v_1$ (top panels) and $v_2$ (bottom panels) of proton and $\Lambda$ baryons (left panels), pions (middle panels) and kaons (right panels) in 10-40% centrality for the $\sqrt{s_{NN}}$ = 3GeV Au+Au collisions. Statistical and systematic uncertainties are shown as bars and gray bands, respectively. Some uncertainties are smaller than the data points. The UrQMD and JAM results are shown as bands:golden, red and blue bands stand for JAM mean-field, UrQMD mean-field and UrQMD cascade mode, respectively. The value of the incompressibility $\kappa$ = 380 MeV is used in the mean-field option. More detailed model descriptions and data comparisons can be found in Supplemental Material.

Rapidity($y$) dependence of $v_1$ (top panels) and $v_2$ (bottom panels) of proton and $\Lambda$ baryons (left panels), pions (middle panels) and kaons (right panels) in 10-40% centrality for the $\sqrt{s_{NN}}$ = 3GeV Au+Au collisions. Statistical and systematic uncertainties are shown as bars and gray bands, respectively. Some uncertainties are smaller than the data points. The UrQMD and JAM results are shown as bands:golden, red and blue bands stand for JAM mean-field, UrQMD mean-field and UrQMD cascade mode, respectively. The value of the incompressibility $\kappa$ = 380 MeV is used in the mean-field option. More detailed model descriptions and data comparisons can be found in Supplemental Material.

$v_2$ scaled by the number of constituent quarks, $v_2/n_q$ , as a function of scaled transverse kinetic energy ($(m_T − m_0)/n_q$) for pions, kaons and protons from Au+Au collisions in 10-40% centrality at $\sqrt{s_{NN}}$ = 3, 27, and 54.4 GeV for positive charged particles (left panel) and negative charged particles (right panel). Colored dashed lines represent the scaling fit to data in 7.7, 14.5, 27, 54.4, and 200 GeV Au+Au collisions from STAR experiment at RHIC [43–45]. Statistical and systematic uncertainties are shown as bars and gray bands, respectively. Some uncertainties are smaller than the data points.

$v_2$ scaled by the number of constituent quarks, $v_2/n_q$ , as a function of scaled transverse kinetic energy ($(m_T − m_0)/n_q$) for pions, kaons and protons from Au+Au collisions in 10-40% centrality at $\sqrt{s_{NN}}$ = 3, 27, and 54.4 GeV for positive charged particles (left panel) and negative charged particles (right panel). Colored dashed lines represent the scaling fit to data in 7.7, 14.5, 27, 54.4, and 200 GeV Au+Au collisions from STAR experiment at RHIC [43–45]. Statistical and systematic uncertainties are shown as bars and gray bands, respectively. Some uncertainties are smaller than the data points.

$v_2$ scaled by the number of constituent quarks, $v_2/n_q$ , as a function of scaled transverse kinetic energy ($(m_T − m_0)/n_q$) for pions, kaons and protons from Au+Au collisions in 10-40% centrality at $\sqrt{s_{NN}}$ = 3, 27, and 54.4 GeV for positive charged particles (left panel) and negative charged particles (right panel). Colored dashed lines represent the scaling fit to data in 7.7, 14.5, 27, 54.4, and 200 GeV Au+Au collisions from STAR experiment at RHIC [43–45]. Statistical and systematic uncertainties are shown as bars and gray bands, respectively. Some uncertainties are smaller than the data points.

$v_2$ scaled by the number of constituent quarks, $v_2/n_q$ , as a function of scaled transverse kinetic energy ($(m_T − m_0)/n_q$) for pions, kaons and protons from Au+Au collisions in 10-40% centrality at $\sqrt{s_{NN}}$ = 3, 27, and 54.4 GeV for positive charged particles (left panel) and negative charged particles (right panel). Colored dashed lines represent the scaling fit to data in 7.7, 14.5, 27, 54.4, and 200 GeV Au+Au collisions from STAR experiment at RHIC [43–45]. Statistical and systematic uncertainties are shown as bars and gray bands, respectively. Some uncertainties are smaller than the data points.

$v_2$ scaled by the number of constituent quarks, $v_2/n_q$ , as a function of scaled transverse kinetic energy ($(m_T − m_0)/n_q$) for pions, kaons and protons from Au+Au collisions in 10-40% centrality at $\sqrt{s_{NN}}$ = 3, 27, and 54.4 GeV for positive charged particles (left panel) and negative charged particles (right panel). Colored dashed lines represent the scaling fit to data in 7.7, 14.5, 27, 54.4, and 200 GeV Au+Au collisions from STAR experiment at RHIC [43–45]. Statistical and systematic uncertainties are shown as bars and gray bands, respectively. Some uncertainties are smaller than the data points.

$v_2$ scaled by the number of constituent quarks, $v_2/n_q$ , as a function of scaled transverse kinetic energy ($(m_T − m_0)/n_q$) for pions, kaons and protons from Au+Au collisions in 10-40% centrality at $\sqrt{s_{NN}}$ = 3, 27, and 54.4 GeV for positive charged particles (left panel) and negative charged particles (right panel). Colored dashed lines represent the scaling fit to data in 7.7, 14.5, 27, 54.4, and 200 GeV Au+Au collisions from STAR experiment at RHIC [43–45]. Statistical and systematic uncertainties are shown as bars and gray bands, respectively. Some uncertainties are smaller than the data points.

$v_2$ scaled by the number of constituent quarks, $v_2/n_q$ , as a function of scaled transverse kinetic energy ($(m_T − m_0)/n_q$) for pions, kaons and protons from Au+Au collisions in 10-40% centrality at $\sqrt{s_{NN}}$ = 3, 27, and 54.4 GeV for positive charged particles (left panel) and negative charged particles (right panel). Colored dashed lines represent the scaling fit to data in 7.7, 14.5, 27, 54.4, and 200 GeV Au+Au collisions from STAR experiment at RHIC [43–45]. Statistical and systematic uncertainties are shown as bars and gray bands, respectively. Some uncertainties are smaller than the data points.

$v_2$ scaled by the number of constituent quarks, $v_2/n_q$ , as a function of scaled transverse kinetic energy ($(m_T − m_0)/n_q$) for pions, kaons and protons from Au+Au collisions in 10-40% centrality at $\sqrt{s_{NN}}$ = 3, 27, and 54.4 GeV for positive charged particles (left panel) and negative charged particles (right panel). Colored dashed lines represent the scaling fit to data in 7.7, 14.5, 27, 54.4, and 200 GeV Au+Au collisions from STAR experiment at RHIC [43–45]. Statistical and systematic uncertainties are shown as bars and gray bands, respectively. Some uncertainties are smaller than the data points.

$v_2$ scaled by the number of constituent quarks, $v_2/n_q$ , as a function of scaled transverse kinetic energy ($(m_T − m_0)/n_q$) for pions, kaons and protons from Au+Au collisions in 10-40% centrality at $\sqrt{s_{NN}}$ = 3, 27, and 54.4 GeV for positive charged particles (left panel) and negative charged particles (right panel). Colored dashed lines represent the scaling fit to data in 7.7, 14.5, 27, 54.4, and 200 GeV Au+Au collisions from STAR experiment at RHIC [43–45]. Statistical and systematic uncertainties are shown as bars and gray bands, respectively. Some uncertainties are smaller than the data points.

$v_2$ scaled by the number of constituent quarks, $v_2/n_q$ , as a function of scaled transverse kinetic energy ($(m_T − m_0)/n_q$) for pions, kaons and protons from Au+Au collisions in 10-40% centrality at $\sqrt{s_{NN}}$ = 3, 27, and 54.4 GeV for positive charged particles (left panel) and negative charged particles (right panel). Colored dashed lines represent the scaling fit to data in 7.7, 14.5, 27, 54.4, and 200 GeV Au+Au collisions from STAR experiment at RHIC [43–45]. Statistical and systematic uncertainties are shown as bars and gray bands, respectively. Some uncertainties are smaller than the data points.

$v_2$ scaled by the number of constituent quarks, $v_2/n_q$ , as a function of scaled transverse kinetic energy ($(m_T − m_0)/n_q$) for pions, kaons and protons from Au+Au collisions in 10-40% centrality at $\sqrt{s_{NN}}$ = 3, 27, and 54.4 GeV for positive charged particles (left panel) and negative charged particles (right panel). Colored dashed lines represent the scaling fit to data in 7.7, 14.5, 27, 54.4, and 200 GeV Au+Au collisions from STAR experiment at RHIC [43–45]. Statistical and systematic uncertainties are shown as bars and gray bands, respectively. Some uncertainties are smaller than the data points.

$v_2$ scaled by the number of constituent quarks, $v_2/n_q$ , as a function of scaled transverse kinetic energy ($(m_T − m_0)/n_q$) for pions, kaons and protons from Au+Au collisions in 10-40% centrality at $\sqrt{s_{NN}}$ = 3, 27, and 54.4 GeV for positive charged particles (left panel) and negative charged particles (right panel). Colored dashed lines represent the scaling fit to data in 7.7, 14.5, 27, 54.4, and 200 GeV Au+Au collisions from STAR experiment at RHIC [43–45]. Statistical and systematic uncertainties are shown as bars and gray bands, respectively. Some uncertainties are smaller than the data points.

$v_2$ scaled by the number of constituent quarks, $v_2/n_q$ , as a function of scaled transverse kinetic energy ($(m_T − m_0)/n_q$) for pions, kaons and protons from Au+Au collisions in 10-40% centrality at $\sqrt{s_{NN}}$ = 3, 27, and 54.4 GeV for positive charged particles (left panel) and negative charged particles (right panel). Colored dashed lines represent the scaling fit to data in 7.7, 14.5, 27, 54.4, and 200 GeV Au+Au collisions from STAR experiment at RHIC [43–45]. Statistical and systematic uncertainties are shown as bars and gray bands, respectively. Some uncertainties are smaller than the data points.

$v_2$ scaled by the number of constituent quarks, $v_2/n_q$ , as a function of scaled transverse kinetic energy ($(m_T − m_0)/n_q$) for pions, kaons and protons from Au+Au collisions in 10-40% centrality at $\sqrt{s_{NN}}$ = 3, 27, and 54.4 GeV for positive charged particles (left panel) and negative charged particles (right panel). Colored dashed lines represent the scaling fit to data in 7.7, 14.5, 27, 54.4, and 200 GeV Au+Au collisions from STAR experiment at RHIC [43–45]. Statistical and systematic uncertainties are shown as bars and gray bands, respectively. Some uncertainties are smaller than the data points.

$v_2$ scaled by the number of constituent quarks, $v_2/n_q$ , as a function of scaled transverse kinetic energy ($(m_T − m_0)/n_q$) for pions, kaons and protons from Au+Au collisions in 10-40% centrality at $\sqrt{s_{NN}}$ = 3, 27, and 54.4 GeV for positive charged particles (left panel) and negative charged particles (right panel). Colored dashed lines represent the scaling fit to data in 7.7, 14.5, 27, 54.4, and 200 GeV Au+Au collisions from STAR experiment at RHIC [43–45]. Statistical and systematic uncertainties are shown as bars and gray bands, respectively. Some uncertainties are smaller than the data points.

$v_2$ scaled by the number of constituent quarks, $v_2/n_q$ , as a function of scaled transverse kinetic energy ($(m_T − m_0)/n_q$) for pions, kaons and protons from Au+Au collisions in 10-40% centrality at $\sqrt{s_{NN}}$ = 3, 27, and 54.4 GeV for positive charged particles (left panel) and negative charged particles (right panel). Colored dashed lines represent the scaling fit to data in 7.7, 14.5, 27, 54.4, and 200 GeV Au+Au collisions from STAR experiment at RHIC [43–45]. Statistical and systematic uncertainties are shown as bars and gray bands, respectively. Some uncertainties are smaller than the data points.

$v_2$ scaled by the number of constituent quarks, $v_2/n_q$ , as a function of scaled transverse kinetic energy ($(m_T − m_0)/n_q$) for pions, kaons and protons from Au+Au collisions in 10-40% centrality at $\sqrt{s_{NN}}$ = 3, 27, and 54.4 GeV for positive charged particles (left panel) and negative charged particles (right panel). Colored dashed lines represent the scaling fit to data in 7.7, 14.5, 27, 54.4, and 200 GeV Au+Au collisions from STAR experiment at RHIC [43–45]. Statistical and systematic uncertainties are shown as bars and gray bands, respectively. Some uncertainties are smaller than the data points.

Collision energy dependence of (top panel) directed flow slope $dv_1/dy|_{y=0}$ for p, $\Lambda, \pi$(combined from $\pi^{\pm}$), K(combined from $K^{\pm}$ and $K_S^0$), $\phi$ and $\Xi^−$, and (bottom panel) elliptic flow $v_2$ for p, $\pi$(combined from $\pi^{\pm}$) in heavy- ion collisions. Collision centrality for all data from RHIC is 10-40%, except for 4.5 GeV where 10-30% is for $dv_1/dy|_{y=0}$ and 0-30% is for $v_2$. Statistical and systematic uncertainties are shown as bars and gray bands, respectively. Some uncertainties are smaller than the data points. The JAM and UrQMD results are shown as colored bands:golden, red and blue bands stand for JAM mean-field, UrQMD mean-field and UrQMD cascade mode, respectively. For clarity the x-axis value of the data points have been shifted.

Collision energy dependence of (top panel) directed flow slope $dv_1/dy|_{y=0}$ for p, $\Lambda, \pi$(combined from $\pi^{\pm}$), K(combined from $K^{\pm}$ and $K_S^0$), $\phi$ and $\Xi^−$, and (bottom panel) elliptic flow $v_2$ for p, $\pi$(combined from $\pi^{\pm}$) in heavy- ion collisions. Collision centrality for all data from RHIC is 10-40%, except for 4.5 GeV where 10-30% is for $dv_1/dy|_{y=0}$ and 0-30% is for $v_2$. Statistical and systematic uncertainties are shown as bars and gray bands, respectively. Some uncertainties are smaller than the data points. The JAM and UrQMD results are shown as colored bands:golden, red and blue bands stand for JAM mean-field, UrQMD mean-field and UrQMD cascade mode, respectively. For clarity the x-axis value of the data points have been shifted.

Collision energy dependence of (top panel) directed flow slope $dv_1/dy|_{y=0}$ for p, $\Lambda, \pi$(combined from $\pi^{\pm}$), K(combined from $K^{\pm}$ and $K_S^0$), $\phi$ and $\Xi^−$, and (bottom panel) elliptic flow $v_2$ for p, $\pi$(combined from $\pi^{\pm}$) in heavy- ion collisions. Collision centrality for all data from RHIC is 10-40%, except for 4.5 GeV where 10-30% is for $dv_1/dy|_{y=0}$ and 0-30% is for $v_2$. Statistical and systematic uncertainties are shown as bars and gray bands, respectively. Some uncertainties are smaller than the data points. The JAM and UrQMD results are shown as colored bands:golden, red and blue bands stand for JAM mean-field, UrQMD mean-field and UrQMD cascade mode, respectively. For clarity the x-axis value of the data points have been shifted.

Collision energy dependence of (top panel) directed flow slope $dv_1/dy|_{y=0}$ for p, $\Lambda, \pi$(combined from $\pi^{\pm}$), K(combined from $K^{\pm}$ and $K_S^0$), $\phi$ and $\Xi^−$, and (bottom panel) elliptic flow $v_2$ for p, $\pi$(combined from $\pi^{\pm}$) in heavy- ion collisions. Collision centrality for all data from RHIC is 10-40%, except for 4.5 GeV where 10-30% is for $dv_1/dy|_{y=0}$ and 0-30% is for $v_2$. Statistical and systematic uncertainties are shown as bars and gray bands, respectively. Some uncertainties are smaller than the data points. The JAM and UrQMD results are shown as colored bands:golden, red and blue bands stand for JAM mean-field, UrQMD mean-field and UrQMD cascade mode, respectively. For clarity the x-axis value of the data points have been shifted.

Collision energy dependence of (top panel) directed flow slope $dv_1/dy|_{y=0}$ for p, $\Lambda, \pi$(combined from $\pi^{\pm}$), K(combined from $K^{\pm}$ and $K_S^0$), $\phi$ and $\Xi^−$, and (bottom panel) elliptic flow $v_2$ for p, $\pi$(combined from $\pi^{\pm}$) in heavy- ion collisions. Collision centrality for all data from RHIC is 10-40%, except for 4.5 GeV where 10-30% is for $dv_1/dy|_{y=0}$ and 0-30% is for $v_2$. Statistical and systematic uncertainties are shown as bars and gray bands, respectively. Some uncertainties are smaller than the data points. The JAM and UrQMD results are shown as colored bands:golden, red and blue bands stand for JAM mean-field, UrQMD mean-field and UrQMD cascade mode, respectively. For clarity the x-axis value of the data points have been shifted.

Collision energy dependence of (top panel) directed flow slope $dv_1/dy|_{y=0}$ for p, $\Lambda, \pi$(combined from $\pi^{\pm}$), K(combined from $K^{\pm}$ and $K_S^0$), $\phi$ and $\Xi^−$, and (bottom panel) elliptic flow $v_2$ for p, $\pi$(combined from $\pi^{\pm}$) in heavy- ion collisions. Collision centrality for all data from RHIC is 10-40%, except for 4.5 GeV where 10-30% is for $dv_1/dy|_{y=0}$ and 0-30% is for $v_2$. Statistical and systematic uncertainties are shown as bars and gray bands, respectively. Some uncertainties are smaller than the data points. The JAM and UrQMD results are shown as colored bands:golden, red and blue bands stand for JAM mean-field, UrQMD mean-field and UrQMD cascade mode, respectively. For clarity the x-axis value of the data points have been shifted.

Collision energy dependence of (top panel) directed flow slope $dv_1/dy|_{y=0}$ for p, $\Lambda, \pi$(combined from $\pi^{\pm}$), K(combined from $K^{\pm}$ and $K_S^0$), $\phi$ and $\Xi^−$, and (bottom panel) elliptic flow $v_2$ for p, $\pi$(combined from $\pi^{\pm}$) in heavy- ion collisions. Collision centrality for all data from RHIC is 10-40%, except for 4.5 GeV where 10-30% is for $dv_1/dy|_{y=0}$ and 0-30% is for $v_2$. Statistical and systematic uncertainties are shown as bars and gray bands, respectively. Some uncertainties are smaller than the data points. The JAM and UrQMD results are shown as colored bands:golden, red and blue bands stand for JAM mean-field, UrQMD mean-field and UrQMD cascade mode, respectively. For clarity the x-axis value of the data points have been shifted.

Collision energy dependence of (top panel) directed flow slope $dv_1/dy|_{y=0}$ for p, $\Lambda, \pi$(combined from $\pi^{\pm}$), K(combined from $K^{\pm}$ and $K_S^0$), $\phi$ and $\Xi^−$, and (bottom panel) elliptic flow $v_2$ for p, $\pi$(combined from $\pi^{\pm}$) in heavy- ion collisions. Collision centrality for all data from RHIC is 10-40%, except for 4.5 GeV where 10-30% is for $dv_1/dy|_{y=0}$ and 0-30% is for $v_2$. Statistical and systematic uncertainties are shown as bars and gray bands, respectively. Some uncertainties are smaller than the data points. The JAM and UrQMD results are shown as colored bands:golden, red and blue bands stand for JAM mean-field, UrQMD mean-field and UrQMD cascade mode, respectively. For clarity the x-axis value of the data points have been shifted.

Collision energy dependence of (top panel) directed flow slope $dv_1/dy|_{y=0}$ for p, $\Lambda, \pi$(combined from $\pi^{\pm}$), K(combined from $K^{\pm}$ and $K_S^0$), $\phi$ and $\Xi^−$, and (bottom panel) elliptic flow $v_2$ for p, $\pi$(combined from $\pi^{\pm}$) in heavy- ion collisions. Collision centrality for all data from RHIC is 10-40%, except for 4.5 GeV where 10-30% is for $dv_1/dy|_{y=0}$ and 0-30% is for $v_2$. Statistical and systematic uncertainties are shown as bars and gray bands, respectively. Some uncertainties are smaller than the data points. The JAM and UrQMD results are shown as colored bands:golden, red and blue bands stand for JAM mean-field, UrQMD mean-field and UrQMD cascade mode, respectively. For clarity the x-axis value of the data points have been shifted.

The centrality dependence of Global $\Lambda$-hyperon Polarization at $\sqrt{s_{\rm NN}}=3$ GeV.

The $p_{\rm T}$ dependence of Global $\Lambda$-hyperon Polarization at $\sqrt{s_{\rm NN}}=3$ GeV.

The rapidity dependence of Global $\Lambda$-hyperon Polarization at $\sqrt{s_{\rm NN}}=3$ GeV.

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.

Dijet signal efficiencies as a function of the signal mass and lifetime for events satisfying all event and vertex requirements, with corrections based on systematic differences in the vertex reconstruction efficiency between data and simulation.

The distribution of $d_{\mathrm{BV}}$ for $\geq$5-track one-vertex events in data and three simulated multijet signal samples each with a mass of 1600 GeV. The production cross section for each signal model is assumed to be the lower limit excluded by CMS-EXO-17-018, corresponding to values of 0.8, 0.25, and 0.15 fb for the samples with $c\tau =$ 0.3, 1.0, and 10 mm, respectively. The last bin includes the overflow events. This bin includes one event in data with a vertex with large $d_{\mathrm{BV}}$ that appears to arise from tracks originating from separate pp interaction vertices, consistent with background.

Distribution of the $x$-$y$ distances between vertices, $d_{\mathrm{VV}}$, for 2017 and 2018 data. The background distribution $d_{\mathrm{VV}}^{\kern 0.15em\mathrm{C}}$ (blue continuous line) is constructed from one-vertex events in data, and is normalized to the number of two-vertex events in data with two 3-track vertices. The two vertical red dashed lines separate the regions used in the fit.

Distribution of the $x$-$y$ distances between vertices, $d_{\mathrm{VV}}$, for 2017 and 2018 data. The background distribution $d_{\mathrm{VV}}^{\kern 0.15em\mathrm{C}}$ (blue continuous line) is constructed from one-vertex events in data, and is normalized to the number of two-vertex events in data which have exactly one 4-track vertex and one 3-track vertex. The two vertical red dashed lines separate the regions used in the fit.

Distribution of the $x$-$y$ distances between vertices, $d_{\mathrm{VV}}$, for 2017 and 2018 data. The background distribution $d_{\mathrm{VV}}^{\kern 0.15em\mathrm{C}}$ (blue continuous line) is constructed from one-vertex events in data, and is normalized to the number of two-vertex events in data with two 4-track vertices. The two vertical red dashed lines separate the regions used in the fit.

Distribution of the $x$-$y$ distances between vertices, $d_{\mathrm{VV}}$, for 2017 and 2018 data. The background distribution $d_{\mathrm{VV}}^{\kern 0.15em\mathrm{C}}$ (blue continuous line) is constructed from one-vertex events in data, and is normalized using $\geq$5-track one-vertex event information. The two vertical red dashed lines separate the regions used in the fit.

Predicted yields for the background-only normalized template, predicted yields for three simulated multijet signals each with a mass of 1600 GeV, and the observed yield in each $d_{\mathrm{VV}}$ bin. The production cross section for each signal model is assumed to be the lower limit excluded by CMS-EXO-17-018, corresponding to values of 0.8, 0.25, and 0.15 fb for samples with $c\tau =$ 0.3, 1.0, and 10 mm, respectively. The uncertainties in the signal yields and the systematic uncertainties in the background prediction reflect the systematic uncertainties given in the text.

Observed 95% CL upper limits on the product of cross section and branching fraction squared for the multijet signals, as a function of mass and $c\tau$. The overlaid mass-lifetime exclusion curves assume pair-production cross sections for the neutralino (red) and gluino (purple) with 100% branching fraction to each model's respective decay mode specified. The solid black (dashed colored) lines represent the observed (median expected) limits at 95% CL. The thin black lines represent the variation of the observed limit within theoretical uncertainties of the signal cross section. The thin dashed colored lines represent the region containing 68% of the expected limit distribution under the background-only hypothesis. The observed limits from the CMS displaced jets search (CMS-EXO-19-021) are also shown in teal for comparison.

Observed 95% CL upper limits on the product of cross section and branching fraction squared for the multijet signals, as a function of mass and $c\tau$. The overlaid mass-lifetime exclusion curves assume pair-production cross sections for the neutralino (red) and gluino (purple) with 100% branching fraction to each model's respective decay mode specified. The solid black (dashed colored) lines represent the observed (median expected) limits at 95% CL. The thin black lines represent the variation of the observed limit within theoretical uncertainties of the signal cross section. The thin dashed colored lines represent the region containing 68% of the expected limit distribution under the background-only hypothesis. The observed limits from the CMS displaced jets search (CMS-EXO-19-021) are also shown in teal for comparison.

Observed 95% CL upper limits on the product of cross section and branching fraction squared for the multijet signals, as a function of mass and $c\tau$. The overlaid mass-lifetime exclusion curves assume pair-production cross sections for the neutralino (red) and gluino (purple) with 100% branching fraction to each model's respective decay mode specified. The solid black (dashed colored) lines represent the observed (median expected) limits at 95% CL. The thin black lines represent the variation of the observed limit within theoretical uncertainties of the signal cross section. The thin dashed colored lines represent the region containing 68% of the expected limit distribution under the background-only hypothesis. The observed limits from the CMS displaced jets search (CMS-EXO-19-021) are also shown in teal for comparison.

Observed 95% CL upper limits on the product of cross section and branching fraction squared for the dijet signals, as a function of mass and $c\tau$. The overlaid mass-lifetime exclusion curves assume pair-production cross sections for the top squark with 100% branching fraction to each model's respective decay mode specified. The solid black (dashed colored) lines represent the observed (median expected) limits at 95% CL. The thin black lines represent the variation of the observed limit within theoretical uncertainties of the signal cross section. The thin dashed colored lines represent the region containing 68% of the expected limit distribution under the background-only hypothesis. The observed limits from the CMS displaced jets search (CMS-EXO-19-021) are also shown in teal for comparison.

Observed 95% CL upper limits on the product of cross section and branching fraction squared for the dijet signals, as a function of mass and $c\tau$. The overlaid mass-lifetime exclusion curves assume pair-production cross sections for the top squark with 100% branching fraction to each model's respective decay mode specified. The solid black (dashed colored) lines represent the observed (median expected) limits at 95% CL. The thin black lines represent the variation of the observed limit within theoretical uncertainties of the signal cross section. The thin dashed colored lines represent the region containing 68% of the expected limit distribution under the background-only hypothesis. The observed limits from the CMS displaced jets search (CMS-EXO-19-021) are also shown in teal for comparison.

Observed and expected 95% CL upper limits on the product of cross section and branching fraction squared, as a function of mass for multijet signals, for a fixed $c\tau$ of 300um in the full Run-2 data set. The neutralino and gluino pair production cross sections are shown for the multijet signals.

Observed and expected 95% CL upper limits on the product of cross section and branching fraction squared, as a function of mass for dijet signals, for a fixed $c\tau$ of 300um in the full Run-2 data set. The top squark pair-production cross section is shown for the dijet signals.

Observed and expected 95% CL upper limits on the product of cross section and branching fraction squared, as a function of mass for multijet signals, for a fixed $c\tau$ of 1 mm in the full Run-2 data set. The neutralino and gluino pair production cross sections are shown for the multijet signals.

Observed and expected 95% CL upper limits on the product of cross section and branching fraction squared, as a function of mass for dijet signals, for a fixed $c\tau$ of 1 mm in the full Run-2 data set. The top squark pair-production cross section is shown for the dijet signals.

Observed and expected 95% CL upper limits on the product of cross section and branching fraction squared, as a function of mass for multijet signals, for a fixed $c\tau$ of 10 mm in the full Run-2 data set. The neutralino and gluino pair production cross sections are shown for the multijet signals.

Observed and expected 95% CL upper limits on the product of cross section and branching fraction squared, as a function of mass for dijet signals, for a fixed $c\tau$ of 10 mm in the full Run-2 data set. The top squark pair-production cross section is shown for the dijet signals.

Observed and expected 95% CL upper limits on the product of cross section and branching fraction squared, as a function of $c\tau$ for multijet signals, for a fixed mass of 800 GeV in the full Run-2 data set. The neutralino and gluino pair production cross sections are shown for the multijet signals.

Observed and expected 95% CL upper limits on the product of cross section and branching fraction squared, as a function of $c\tau$ for dijet signals, for a fixed mass of 800 GeV in the full Run-2 data set. The top squark pair-production cross section is shown for the dijet signals.

Observed and expected 95% CL upper limits on the product of cross section and branching fraction squared, as a function of $c\tau$ for multijet signals, for a fixed mass of 1600 GeV in the full Run-2 data set. The neutralino and gluino pair production cross sections are shown for the multijet signals.

Observed and expected 95% CL upper limits on the product of cross section and branching fraction squared, as a function of $c\tau$ for dijet signals, for a fixed mass of 1600 GeV in the full Run-2 data set. The top squark pair-production cross section is shown for the dijet signals.

Observed and expected 95% CL upper limits on the product of cross section and branching fraction squared, as a function of $c\tau$ for multijet signals, for a fixed mass of 2400 GeV in the full Run-2 data set. The neutralino and gluino pair production cross sections are shown for the multijet signals. For $m$ = 2400 GeV, the expected neutralino cross section is $\approx 8\times 10^{-5}$ fb and is not shown.

Observed and expected 95% CL upper limits on the product of cross section and branching fraction squared, as a function of $c\tau$ for dijet signals, for a fixed mass of 2400 GeV in the full Run-2 data set. The top squark pair-production cross section is shown for the dijet signals.

Data-to-simulation efficiency correction factors, shown for multijet and dijet signal topologies in several ranges of $c\tau$. Note that these correction factors account for the two long-lived particles in the simulated events, and are therefore the total correction factors used to scale event yields rather than the correction factors one would apply to individual vertices.

Distribution of the azimuthal angle between vertices, $\Delta\phi_{\mathrm{VV}}$, for 2017 and 2018 data. The background distribution (blue continuous line) is constructed from 3-track one-vertex events in data, and is normalized to the number of 3-track two-vertex events in data.

Distribution of the azimuthal angle between vertices, $\Delta\phi_{\mathrm{VV}}$, for 2017 and 2018 data. The background distribution (blue continuous line) is constructed from 4-track and 3-track one-vertex events in data, and is normalized to the number of two-vertex events in data which have exactly one 4-track vertex and one 3-track vertex.

Distribution of the azimuthal angle between vertices, $\Delta\phi_{\mathrm{VV}}$, for 2017 and 2018 data. The background distribution (blue continuous line) is constructed from 4-track one-vertex events in data, and is normalized to the number of 4-track two-vertex events in data.

Distribution of the azimuthal angle between vertices, $\Delta\phi_{\mathrm{VV}}$, for 2017 and 2018 data. The background distribution (blue continuous line) is constructed from $\geq$5-track one-vertex events in data, and is normalized using one-vertex event information. No $\geq$5-track two-vertex data events pass the selection.

The Z boson differential cross section as a function of |yZ|.

The Z boson differential cross section as a function of |yZ|.

The Z boson differential cross section as a function of |yZ|.

The Z boson differential cross section as a function of pTZ.

The Z boson differential cross section as a function of pTZ.

The Z boson differential cross section as a function of pTZ.

The TAA-normalized yields of Z bosons as a function of centrality.

The TAA-normalized yields of Z bosons as a function of centrality.

The TAA-normalized yields of Z bosons as a function of centrality.

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. 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. 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. 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. 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$. 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$. 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$. 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$. 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$. 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$. 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}}$. 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}}$. 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}}$. 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}}$. 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}}$. 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}}$. 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. 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. Updated to correct a small bug that had shifted the jet mass to slightly smaller values.

The $p_{\text{T}}$ dependence of $v_{n}$ coefficients at mid-rapidity ($|y| <$ 1) in minimum bias U+U collisions at $\sqrt{s_{NN}}$ = 193 GeV. The error bars represent statistical uncertainties. The bands represent point-by-point systematic uncertainties.

The $p_{\text{T}}$ dependence of $v_{n}$ coefficients at mid-rapidity ($|y| <$ 1) in minimum bias U+U collisions at $\sqrt{s_{NN}}$ = 193 GeV. The error bars represent statistical uncertainties. The bands represent point-by-point systematic uncertainties.

The $p_{\text{T}}$ dependence of $v_{n}$ coefficients at mid-rapidity ($|y| <$ 1) in minimum bias U+U collisions at $\sqrt{s_{NN}}$ = 193 GeV. The error bars represent statistical uncertainties. The bands represent point-by-point systematic uncertainties.

The $p_{\text{T}}$ dependence of $v_{n}$ coefficients at mid-rapidity ($|y| <$ 1) in minimum bias U+U collisions at $\sqrt{s_{NN}}$ = 193 GeV. The error bars represent statistical uncertainties. The bands represent point-by-point systematic uncertainties.

The $p_{\text{T}}$ dependence of $v_{n}$ coefficients at mid-rapidity ($|y| <$ 1) in minimum bias U+U collisions at $\sqrt{s_{NN}}$ = 193 GeV. The error bars represent statistical uncertainties. The bands represent point-by-point systematic uncertainties.

The $p_{\text{T}}$ dependence of $v_{n}$ coefficients at mid-rapidity ($|y| <$ 1) in minimum bias U+U collisions at $\sqrt{s_{NN}}$ = 193 GeV. The error bars represent statistical uncertainties. The bands represent point-by-point systematic uncertainties.

The $p_{\text{T}}$ dependence of $v_{n}$ coefficients at mid-rapidity ($|y| <$ 1) in minimum bias U+U collisions at $\sqrt{s_{NN}}$ = 193 GeV. The error bars represent statistical uncertainties. The bands represent point-by-point systematic uncertainties.

The $p_{\text{T}}$ dependence of $v_{n}$ coefficients at mid-rapidity ($|y| <$ 1) in minimum bias U+U collisions at $\sqrt{s_{NN}}$ = 193 GeV. The error bars represent statistical uncertainties. The bands represent point-by-point systematic uncertainties.

The $p_{\text{T}}$ dependence of $v_{n}$ coefficients at mid-rapidity ($|y| <$ 1) in minimum bias U+U collisions at $\sqrt{s_{NN}}$ = 193 GeV. The error bars represent statistical uncertainties. The bands represent point-by-point systematic uncertainties.

The $p_{\text{T}}$ dependence of $v_{n}$ coefficients at mid-rapidity ($|y| <$ 1) in minimum bias U+U collisions at $\sqrt{s_{NN}}$ = 193 GeV. The error bars represent statistical uncertainties. The bands represent point-by-point systematic uncertainties.

The $p_{\text{T}}$ dependence of $v_{n}$ coefficients at mid-rapidity ($|y| <$ 1) in minimum bias U+U collisions at $\sqrt{s_{NN}}$ = 193 GeV. The error bars represent statistical uncertainties. The bands represent point-by-point systematic uncertainties.

The $p_{\text{T}}$ dependence of $v_{n}$ coefficients at mid-rapidity ($|y| <$ 1) in minimum bias U+U collisions at $\sqrt{s_{NN}}$ = 193 GeV. The error bars represent statistical uncertainties. The bands represent point-by-point systematic uncertainties.

The $p_{\text{T}}$ dependence of $v_{n}$ coefficients at mid-rapidity ($|y| <$ 1) in minimum bias U+U collisions at $\sqrt{s_{NN}}$ = 193 GeV. The error bars represent statistical uncertainties. The bands represent point-by-point systematic uncertainties.

The flow coefficients $v_{n}$ as a function of $p_{\text{T}}$ at mid-rapidity ($|y| <$ 1) in U+U collisions at $\sqrt{s_{NN}}$ = 193 GeV for different centrality classes. The error bars represent statistical uncertainties. The bands represent point-by-point systematic uncertainties.

The flow coefficients $v_{n}$ as a function of $p_{\text{T}}$ at mid-rapidity ($|y| <$ 1) in U+U collisions at $\sqrt{s_{NN}}$ = 193 GeV for different centrality classes. The error bars represent statistical uncertainties. The bands represent point-by-point systematic uncertainties.

The flow coefficients $v_{n}$ as a function of $p_{\text{T}}$ at mid-rapidity ($|y| <$ 1) in U+U collisions at $\sqrt{s_{NN}}$ = 193 GeV for different centrality classes. The error bars represent statistical uncertainties. The bands represent point-by-point systematic uncertainties.

The flow coefficients $v_{n}$ as a function of $p_{\text{T}}$ at mid-rapidity ($|y| <$ 1) in U+U collisions at $\sqrt{s_{NN}}$ = 193 GeV for different centrality classes. The error bars represent statistical uncertainties. The bands represent point-by-point systematic uncertainties.

The flow coefficients $v_{n}$ as a function of $p_{\text{T}}$ at mid-rapidity ($|y| <$ 1) in U+U collisions at $\sqrt{s_{NN}}$ = 193 GeV for different centrality classes. The error bars represent statistical uncertainties. The bands represent point-by-point systematic uncertainties.

The flow coefficients $v_{n}$ as a function of $p_{\text{T}}$ at mid-rapidity ($|y| <$ 1) in U+U collisions at $\sqrt{s_{NN}}$ = 193 GeV for different centrality classes. The error bars represent statistical uncertainties. The bands represent point-by-point systematic uncertainties.

The flow coefficients $v_{n}$ as a function of $p_{\text{T}}$ at mid-rapidity ($|y| <$ 1) in U+U collisions at $\sqrt{s_{NN}}$ = 193 GeV for different centrality classes. The error bars represent statistical uncertainties. The bands represent point-by-point systematic uncertainties.

The flow coefficients $v_{n}$ as a function of $p_{\text{T}}$ at mid-rapidity ($|y| <$ 1) in U+U collisions at $\sqrt{s_{NN}}$ = 193 GeV for different centrality classes. The error bars represent statistical uncertainties. The bands represent point-by-point systematic uncertainties.

The flow coefficients $v_{n}$ as a function of $p_{\text{T}}$ at mid-rapidity ($|y| <$ 1) in U+U collisions at $\sqrt{s_{NN}}$ = 193 GeV for different centrality classes. The error bars represent statistical uncertainties. The bands represent point-by-point systematic uncertainties.

The flow coefficients $v_{n}$ as a function of $p_{\text{T}}$ at mid-rapidity ($|y| <$ 1) in U+U collisions at $\sqrt{s_{NN}}$ = 193 GeV for different centrality classes. The error bars represent statistical uncertainties. The bands represent point-by-point systematic uncertainties.

The flow coefficients $v_{n}$ as a function of $p_{\text{T}}$ at mid-rapidity ($|y| <$ 1) in U+U collisions at $\sqrt{s_{NN}}$ = 193 GeV for different centrality classes. The error bars represent statistical uncertainties. The bands represent point-by-point systematic uncertainties.

The flow coefficients $v_{n}$ as a function of $p_{\text{T}}$ at mid-rapidity ($|y| <$ 1) in U+U collisions at $\sqrt{s_{NN}}$ = 193 GeV for different centrality classes. The error bars represent statistical uncertainties. The bands represent point-by-point systematic uncertainties.

The flow coefficients $v_{n}$ as a function of $p_{\text{T}}$ at mid-rapidity ($|y| <$ 1) in U+U collisions at $\sqrt{s_{NN}}$ = 193 GeV for different centrality classes. The error bars represent statistical uncertainties. The bands represent point-by-point systematic uncertainties.

The flow coefficients $v_{n}$ as a function of $p_{\text{T}}$ at mid-rapidity ($|y| <$ 1) in U+U collisions at $\sqrt{s_{NN}}$ = 193 GeV for different centrality classes. The error bars represent statistical uncertainties. The bands represent point-by-point systematic uncertainties.

The flow coefficients $v_{n}$ as a function of $p_{\text{T}}$ at mid-rapidity ($|y| <$ 1) in U+U collisions at $\sqrt{s_{NN}}$ = 193 GeV for different centrality classes. The error bars represent statistical uncertainties. The bands represent point-by-point systematic uncertainties.

The flow coefficients $v_{n}$ as a function of $p_{\text{T}}$ at mid-rapidity ($|y| <$ 1) in U+U collisions at $\sqrt{s_{NN}}$ = 193 GeV for different centrality classes. The error bars represent statistical uncertainties. The bands represent point-by-point systematic uncertainties.

The flow coefficients $v_{n}$ as a function of $p_{\text{T}}$ at mid-rapidity ($|y| <$ 1) in U+U collisions at $\sqrt{s_{NN}}$ = 193 GeV for different centrality classes. The error bars represent statistical uncertainties. The bands represent point-by-point systematic uncertainties.

The flow coefficients $v_{n}$ as a function of $p_{\text{T}}$ at mid-rapidity ($|y| <$ 1) in U+U collisions at $\sqrt{s_{NN}}$ = 193 GeV for different centrality classes. The error bars represent statistical uncertainties. The bands represent point-by-point systematic uncertainties.

The flow coefficients $v_{n}$ as a function of $p_{\text{T}}$ at mid-rapidity ($|y| <$ 1) in U+U collisions at $\sqrt{s_{NN}}$ = 193 GeV for different centrality classes. The error bars represent statistical uncertainties. The bands represent point-by-point systematic uncertainties.

The flow coefficients $v_{n}$ as a function of $p_{\text{T}}$ at mid-rapidity ($|y| <$ 1) in U+U collisions at $\sqrt{s_{NN}}$ = 193 GeV for different centrality classes. The error bars represent statistical uncertainties. The bands represent point-by-point systematic uncertainties.

The flow coefficients $v_{n}$ as a function of $p_{\text{T}}$ at mid-rapidity ($|y| <$ 1) in U+U collisions at $\sqrt{s_{NN}}$ = 193 GeV for different centrality classes. The error bars represent statistical uncertainties. The bands represent point-by-point systematic uncertainties.

The flow coefficients $v_{n}$ as a function of $p_{\text{T}}$ at mid-rapidity ($|y| <$ 1) in U+U collisions at $\sqrt{s_{NN}}$ = 193 GeV for different centrality classes. The error bars represent statistical uncertainties. The bands represent point-by-point systematic uncertainties.

The flow coefficients $v_{n}$ as a function of $p_{\text{T}}$ at mid-rapidity ($|y| <$ 1) in U+U collisions at $\sqrt{s_{NN}}$ = 193 GeV for different centrality classes. The error bars represent statistical uncertainties. The bands represent point-by-point systematic uncertainties.

The flow coefficients $v_{n}$ as a function of $p_{\text{T}}$ at mid-rapidity ($|y| <$ 1) in U+U collisions at $\sqrt{s_{NN}}$ = 193 GeV for different centrality classes. The error bars represent statistical uncertainties. The bands represent point-by-point systematic uncertainties.

The flow coefficients $v_{n}$ as a function of $p_{\text{T}}$ at mid-rapidity ($|y| <$ 1) in U+U collisions at $\sqrt{s_{NN}}$ = 193 GeV for different centrality classes. The error bars represent statistical uncertainties. The bands represent point-by-point systematic uncertainties.

The flow coefficients $v_{n}$ as a function of $p_{\text{T}}$ at mid-rapidity ($|y| <$ 1) in U+U collisions at $\sqrt{s_{NN}}$ = 193 GeV for different centrality classes. The error bars represent statistical uncertainties. The bands represent point-by-point systematic uncertainties.

The flow coefficients $v_{n}$ as a function of $p_{\text{T}}$ at mid-rapidity ($|y| <$ 1) in U+U collisions at $\sqrt{s_{NN}}$ = 193 GeV for different centrality classes. The error bars represent statistical uncertainties. The bands represent point-by-point systematic uncertainties.

The flow coefficients $v_{n}$ as a function of $p_{\text{T}}$ at mid-rapidity ($|y| <$ 1) in U+U collisions at $\sqrt{s_{NN}}$ = 193 GeV for different centrality classes. The error bars represent statistical uncertainties. The bands represent point-by-point systematic uncertainties.

The flow coefficients $v_{n}$ as a function of $p_{\text{T}}$ at mid-rapidity ($|y| <$ 1) in U+U collisions at $\sqrt{s_{NN}}$ = 193 GeV for different centrality classes. The error bars represent statistical uncertainties. The bands represent point-by-point systematic uncertainties.

The flow coefficients $v_{n}$ as a function of $p_{\text{T}}$ at mid-rapidity ($|y| <$ 1) in U+U collisions at $\sqrt{s_{NN}}$ = 193 GeV for different centrality classes. The error bars represent statistical uncertainties. The bands represent point-by-point systematic uncertainties.

The flow coefficients $v_{n}$ as a function of $p_{\text{T}}$ at mid-rapidity ($|y| <$ 1) in U+U collisions at $\sqrt{s_{NN}}$ = 193 GeV for different centrality classes. The error bars represent statistical uncertainties. The bands represent point-by-point systematic uncertainties.

The flow coefficients $v_{n}$ as a function of $p_{\text{T}}$ at mid-rapidity ($|y| <$ 1) in U+U collisions at $\sqrt{s_{NN}}$ = 193 GeV for different centrality classes. The error bars represent statistical uncertainties. The bands represent point-by-point systematic uncertainties.

The flow coefficients $v_{n}$ as a function of $p_{\text{T}}$ at mid-rapidity ($|y| <$ 1) in U+U collisions at $\sqrt{s_{NN}}$ = 193 GeV for different centrality classes. The error bars represent statistical uncertainties. The bands represent point-by-point systematic uncertainties.

The flow coefficients $v_{n}$ as a function of $p_{\text{T}}$ at mid-rapidity ($|y| <$ 1) in U+U collisions at $\sqrt{s_{NN}}$ = 193 GeV for different centrality classes. The error bars represent statistical uncertainties. The bands represent point-by-point systematic uncertainties.

The flow coefficients $v_{n}$ as a function of $p_{\text{T}}$ at mid-rapidity ($|y| <$ 1) in U+U collisions at $\sqrt{s_{NN}}$ = 193 GeV for different centrality classes. The error bars represent statistical uncertainties. The bands represent point-by-point systematic uncertainties.

The flow coefficients $v_{n}$ as a function of $p_{\text{T}}$ at mid-rapidity ($|y| <$ 1) in U+U collisions at $\sqrt{s_{NN}}$ = 193 GeV for different centrality classes. The error bars represent statistical uncertainties. The bands represent point-by-point systematic uncertainties.

The flow coefficients $v_{n}$ as a function of $p_{\text{T}}$ at mid-rapidity ($|y| <$ 1) in U+U collisions at $\sqrt{s_{NN}}$ = 193 GeV for different centrality classes. The error bars represent statistical uncertainties. The bands represent point-by-point systematic uncertainties.

The flow coefficients $v_{n}$ as a function of $p_{\text{T}}$ at mid-rapidity ($|y| <$ 1) in U+U collisions at $\sqrt{s_{NN}}$ = 193 GeV for different centrality classes. The error bars represent statistical uncertainties. The bands represent point-by-point systematic uncertainties.

The flow coefficients $v_{n}$ as a function of $p_{\text{T}}$ at mid-rapidity ($|y| <$ 1) in U+U collisions at $\sqrt{s_{NN}}$ = 193 GeV for different centrality classes. The error bars represent statistical uncertainties. The bands represent point-by-point systematic uncertainties.

The flow coefficients $v_{n}$ as a function of $p_{\text{T}}$ at mid-rapidity ($|y| <$ 1) in U+U collisions at $\sqrt{s_{NN}}$ = 193 GeV for different centrality classes. The error bars represent statistical uncertainties. The bands represent point-by-point systematic uncertainties.

Flow coefficients $v_{n}$ as a function of transverse kinetic energy $KE_{\text{T}}/n_{q}$ for various particles at mid-rapidity ($|y| <$ 1) in U+U collisions at $\sqrt{s_{NN}}$ = 193 GeV, scaled by the number of constituent quarks $(n_{q})$ to the power $n/2$. The error bars represent statistical uncertainties. The bands represent point-by-point systematic uncertainties.

Flow coefficients $v_{n}$ as a function of transverse kinetic energy $KE_{\text{T}}/n_{q}$ for various particles at mid-rapidity ($|y| <$ 1) in U+U collisions at $\sqrt{s_{NN}}$ = 193 GeV, scaled by the number of constituent quarks $(n_{q})$ to the power $n/2$. The error bars represent statistical uncertainties. The bands represent point-by-point systematic uncertainties.

Flow coefficients $v_{n}$ as a function of transverse kinetic energy $KE_{\text{T}}/n_{q}$ for various particles at mid-rapidity ($|y| <$ 1) in U+U collisions at $\sqrt{s_{NN}}$ = 193 GeV, scaled by the number of constituent quarks $(n_{q})$ to the power $n/2$. The error bars represent statistical uncertainties. The bands represent point-by-point systematic uncertainties.

Flow coefficients $v_{n}$ as a function of transverse kinetic energy $KE_{\text{T}}/n_{q}$ for various particles at mid-rapidity ($|y| <$ 1) in U+U collisions at $\sqrt{s_{NN}}$ = 193 GeV, scaled by the number of constituent quarks $(n_{q})$ to the power $n/2$. The error bars represent statistical uncertainties. The bands represent point-by-point systematic uncertainties.

Flow coefficients $v_{n}$ as a function of transverse kinetic energy $KE_{\text{T}}/n_{q}$ for various particles at mid-rapidity ($|y| <$ 1) in U+U collisions at $\sqrt{s_{NN}}$ = 193 GeV, scaled by the number of constituent quarks $(n_{q})$ to the power $n/2$. The error bars represent statistical uncertainties. The bands represent point-by-point systematic uncertainties.

Flow coefficients $v_{n}$ as a function of transverse kinetic energy $KE_{\text{T}}/n_{q}$ for various particles at mid-rapidity ($|y| <$ 1) in U+U collisions at $\sqrt{s_{NN}}$ = 193 GeV, scaled by the number of constituent quarks $(n_{q})$ to the power $n/2$. The error bars represent statistical uncertainties. The bands represent point-by-point systematic uncertainties.

Flow coefficients $v_{n}$ as a function of transverse kinetic energy $KE_{\text{T}}/n_{q}$ for various particles at mid-rapidity ($|y| <$ 1) in U+U collisions at $\sqrt{s_{NN}}$ = 193 GeV, scaled by the number of constituent quarks $(n_{q})$ to the power $n/2$. The error bars represent statistical uncertainties. The bands represent point-by-point systematic uncertainties.

Flow coefficients $v_{n}$ as a function of transverse kinetic energy $KE_{\text{T}}/n_{q}$ for various particles at mid-rapidity ($|y| <$ 1) in U+U collisions at $\sqrt{s_{NN}}$ = 193 GeV, scaled by the number of constituent quarks $(n_{q})$ to the power $n/2$. The error bars represent statistical uncertainties. The bands represent point-by-point systematic uncertainties.

Flow coefficients $v_{n}$ as a function of transverse kinetic energy $KE_{\text{T}}/n_{q}$ for various particles at mid-rapidity ($|y| <$ 1) in U+U collisions at $\sqrt{s_{NN}}$ = 193 GeV, scaled by the number of constituent quarks $(n_{q})$ to the power $n/2$. The error bars represent statistical uncertainties. The bands represent point-by-point systematic uncertainties.

Flow coefficients $v_{n}$ as a function of transverse kinetic energy $KE_{\text{T}}/n_{q}$ for various particles at mid-rapidity ($|y| <$ 1) in U+U collisions at $\sqrt{s_{NN}}$ = 193 GeV, scaled by the number of constituent quarks $(n_{q})$ to the power $n/2$. The error bars represent statistical uncertainties. The bands represent point-by-point systematic uncertainties.

Flow coefficients $v_{n}$ as a function of transverse kinetic energy $KE_{\text{T}}/n_{q}$ for various particles at mid-rapidity ($|y| <$ 1) in U+U collisions at $\sqrt{s_{NN}}$ = 193 GeV, scaled by the number of constituent quarks $(n_{q})$ to the power $n/2$. The error bars represent statistical uncertainties. The bands represent point-by-point systematic uncertainties.

Flow coefficients $v_{n}$ as a function of transverse kinetic energy $KE_{\text{T}}/n_{q}$ for various particles at mid-rapidity ($|y| <$ 1) in U+U collisions at $\sqrt{s_{NN}}$ = 193 GeV, scaled by the number of constituent quarks $(n_{q})$ to the power $n/2$. The error bars represent statistical uncertainties. The bands represent point-by-point systematic uncertainties.

Flow coefficients $v_{n}$ as a function of transverse kinetic energy $KE_{\text{T}}/n_{q}$ for various particles at mid-rapidity ($|y| <$ 1) in U+U collisions at $\sqrt{s_{NN}}$ = 193 GeV, scaled by the number of constituent quarks $(n_{q})$ to the power $n/2$. The error bars represent statistical uncertainties. The bands represent point-by-point systematic uncertainties.

Flow coefficients $v_{n}$ as a function of transverse kinetic energy $KE_{\text{T}}/n_{q}$ for various particles at mid-rapidity ($|y| <$ 1) in U+U collisions at $\sqrt{s_{NN}}$ = 193 GeV, scaled by the number of constituent quarks $(n_{q})$ to the power $n/2$. The error bars represent statistical uncertainties. The bands represent point-by-point systematic uncertainties.

Flow coefficients $v_{n}$ as a function of transverse kinetic energy $KE_{\text{T}}/n_{q}$ for various particles at mid-rapidity ($|y| <$ 1) in U+U collisions at $\sqrt{s_{NN}}$ = 193 GeV, scaled by the number of constituent quarks $(n_{q})$ to the power $n/2$. The error bars represent statistical uncertainties. The bands represent point-by-point systematic uncertainties.

Flow coefficients $v_{n}$ as a function of transverse kinetic energy $KE_{\text{T}}/n_{q}$ for various particles at mid-rapidity ($|y| <$ 1) in U+U collisions at $\sqrt{s_{NN}}$ = 193 GeV, scaled by the number of constituent quarks $(n_{q})$ to the power $n/2$. The error bars represent statistical uncertainties. The bands represent point-by-point systematic uncertainties.

Flow coefficients $v_{n}$ as a function of transverse kinetic energy $KE_{\text{T}}/n_{q}$ for various particles at mid-rapidity ($|y| <$ 1) in U+U collisions at $\sqrt{s_{NN}}$ = 193 GeV, scaled by the number of constituent quarks $(n_{q})$ to the power $n/2$. The error bars represent statistical uncertainties. The bands represent point-by-point systematic uncertainties.

Flow coefficients $v_{n}$ as a function of transverse kinetic energy $KE_{\text{T}}/n_{q}$ for various particles at mid-rapidity ($|y| <$ 1) in U+U collisions at $\sqrt{s_{NN}}$ = 193 GeV, scaled by the number of constituent quarks $(n_{q})$ to the power $n/2$. The error bars represent statistical uncertainties. The bands represent point-by-point systematic uncertainties.

Flow coefficients $v_{n}$ as a function of transverse kinetic energy $KE_{\text{T}}/n_{q}$ for various particles at mid-rapidity ($|y| <$ 1) in U+U collisions at $\sqrt{s_{NN}}$ = 193 GeV, scaled by the number of constituent quarks $(n_{q})$ to the power $n/2$. The error bars represent statistical uncertainties. The bands represent point-by-point systematic uncertainties.

Flow coefficients $v_{n}$ as a function of transverse kinetic energy $KE_{\text{T}}/n_{q}$ for various particles at mid-rapidity ($|y| <$ 1) in U+U collisions at $\sqrt{s_{NN}}$ = 193 GeV, scaled by the number of constituent quarks $(n_{q})$ to the power $n/2$. The error bars represent statistical uncertainties. The bands represent point-by-point systematic uncertainties.

Flow coefficients $v_{n}$ as a function of transverse kinetic energy $KE_{\text{T}}/n_{q}$ for various particles at mid-rapidity ($|y| <$ 1) in U+U collisions at $\sqrt{s_{NN}}$ = 193 GeV, scaled by the number of constituent quarks $(n_{q})$ to the power $n/2$. The error bars represent statistical uncertainties. The bands represent point-by-point systematic uncertainties.

Flow coefficients $v_{n}$ as a function of transverse kinetic energy $KE_{\text{T}}/n_{q}$ for various particles at mid-rapidity ($|y| <$ 1) in U+U collisions at $\sqrt{s_{NN}}$ = 193 GeV, scaled by the number of constituent quarks $(n_{q})$ to the power $n/2$. The error bars represent statistical uncertainties. The bands represent point-by-point systematic uncertainties.

Flow coefficients $v_{n}$ as a function of transverse kinetic energy $KE_{\text{T}}/n_{q}$ for various particles at mid-rapidity ($|y| <$ 1) in U+U collisions at $\sqrt{s_{NN}}$ = 193 GeV, scaled by the number of constituent quarks $(n_{q})$ to the power $n/2$. The error bars represent statistical uncertainties. The bands represent point-by-point systematic uncertainties.

Flow coefficients $v_{n}$ as a function of transverse kinetic energy $KE_{\text{T}}/n_{q}$ for various particles at mid-rapidity ($|y| <$ 1) in U+U collisions at $\sqrt{s_{NN}}$ = 193 GeV, scaled by the number of constituent quarks $(n_{q})$ to the power $n/2$. The error bars represent statistical uncertainties. The bands represent point-by-point systematic uncertainties.

Flow coefficients $v_{n}$ as a function of transverse kinetic energy $KE_{\text{T}}/n_{q}$ for various particles at mid-rapidity ($|y| <$ 1) in U+U collisions at $\sqrt{s_{NN}}$ = 193 GeV, scaled by the number of constituent quarks $(n_{q})$ to the power $n/2$. The error bars represent statistical uncertainties. The bands represent point-by-point systematic uncertainties.

Flow coefficients $v_{n}$ as a function of transverse kinetic energy $KE_{\text{T}}/n_{q}$ for various particles at mid-rapidity ($|y| <$ 1) in U+U collisions at $\sqrt{s_{NN}}$ = 193 GeV, scaled by the number of constituent quarks $(n_{q})$ to the power $n/2$. The error bars represent statistical uncertainties. The bands represent point-by-point systematic uncertainties.

Flow coefficients $v_{n}$ as a function of transverse kinetic energy $KE_{\text{T}}/n_{q}$ for various particles at mid-rapidity ($|y| <$ 1) in U+U collisions at $\sqrt{s_{NN}}$ = 193 GeV, scaled by the number of constituent quarks $(n_{q})$ to the power $n/2$. The error bars represent statistical uncertainties. The bands represent point-by-point systematic uncertainties.

Flow coefficients $v_{n}$ as a function of transverse kinetic energy $KE_{\text{T}}/n_{q}$ for various particles at mid-rapidity ($|y| <$ 1) in U+U collisions at $\sqrt{s_{NN}}$ = 193 GeV, scaled by the number of constituent quarks $(n_{q})$ to the power $n/2$. The error bars represent statistical uncertainties. The bands represent point-by-point systematic uncertainties.

Flow coefficients $v_{n}$ as a function of transverse kinetic energy $KE_{\text{T}}/n_{q}$ for various particles at mid-rapidity ($|y| <$ 1) in U+U collisions at $\sqrt{s_{NN}}$ = 193 GeV, scaled by the number of constituent quarks $(n_{q})$ to the power $n/2$. The error bars represent statistical uncertainties. The bands represent point-by-point systematic uncertainties.

Flow coefficients $v_{n}$ as a function of transverse kinetic energy $KE_{\text{T}}/n_{q}$ for various particles at mid-rapidity ($|y| <$ 1) in U+U collisions at $\sqrt{s_{NN}}$ = 193 GeV, scaled by the number of constituent quarks $(n_{q})$ to the power $n/2$. The error bars represent statistical uncertainties. The bands represent point-by-point systematic uncertainties.

The flow coefficients $v_{n}$ scaled by $\varepsilon_{n}\left\lbrace 2\right\rbrace$ as a function of $p_{\text{T}}$ at mid-rapidity ($|y| <$ 1) in U+U collisions at $\sqrt{s_{NN}}$ = 193 GeV. The error bars represent statistical uncertainties. The bands represent point-by-point systematic uncertainties.

The flow coefficients $v_{n}$ scaled by $\varepsilon_{n}\left\lbrace 2\right\rbrace$ as a function of $p_{\text{T}}$ at mid-rapidity ($|y| <$ 1) in U+U collisions at $\sqrt{s_{NN}}$ = 193 GeV. The error bars represent statistical uncertainties. The bands represent point-by-point systematic uncertainties.

The flow coefficients $v_{n}$ scaled by $\varepsilon_{n}\left\lbrace 2\right\rbrace$ as a function of $p_{\text{T}}$ at mid-rapidity ($|y| <$ 1) in U+U collisions at $\sqrt{s_{NN}}$ = 193 GeV. The error bars represent statistical uncertainties. The bands represent point-by-point systematic uncertainties.

The flow coefficients $v_{n}$ scaled by $\varepsilon_{n}\left\lbrace 2\right\rbrace$ as a function of $p_{\text{T}}$ at mid-rapidity ($|y| <$ 1) in U+U collisions at $\sqrt{s_{NN}}$ = 193 GeV. The error bars represent statistical uncertainties. The bands represent point-by-point systematic uncertainties.

The flow coefficients $v_{n}$ scaled by $\varepsilon_{n}\left\lbrace 2\right\rbrace$ as a function of $p_{\text{T}}$ at mid-rapidity ($|y| <$ 1) in U+U collisions at $\sqrt{s_{NN}}$ = 193 GeV. The error bars represent statistical uncertainties. The bands represent point-by-point systematic uncertainties.

The flow coefficients $v_{n}$ scaled by $\varepsilon_{n}\left\lbrace 2\right\rbrace$ as a function of $p_{\text{T}}$ at mid-rapidity ($|y| <$ 1) in U+U collisions at $\sqrt{s_{NN}}$ = 193 GeV. The error bars represent statistical uncertainties. The bands represent point-by-point systematic uncertainties.

The flow coefficients $v_{n}$ scaled by $\varepsilon_{n}\left\lbrace 2\right\rbrace$ as a function of $p_{\text{T}}$ at mid-rapidity ($|y| <$ 1) in U+U collisions at $\sqrt{s_{NN}}$ = 193 GeV. The error bars represent statistical uncertainties. The bands represent point-by-point systematic uncertainties.

The flow coefficients $v_{n}$ scaled by $\varepsilon_{n}\left\lbrace 2\right\rbrace$ as a function of $p_{\text{T}}$ at mid-rapidity ($|y| <$ 1) in U+U collisions at $\sqrt{s_{NN}}$ = 193 GeV. The error bars represent statistical uncertainties. The bands represent point-by-point systematic uncertainties.

The flow coefficients $v_{n}$ scaled by $\varepsilon_{n}\left\lbrace 2\right\rbrace$ as a function of $p_{\text{T}}$ at mid-rapidity ($|y| <$ 1) in U+U collisions at $\sqrt{s_{NN}}$ = 193 GeV. The error bars represent statistical uncertainties. The bands represent point-by-point systematic uncertainties.

The flow coefficients $v_{n}$ scaled by $\varepsilon_{n}\left\lbrace 2\right\rbrace$ as a function of $p_{\text{T}}$ at mid-rapidity ($|y| <$ 1) in U+U collisions at $\sqrt{s_{NN}}$ = 193 GeV. The error bars represent statistical uncertainties. The bands represent point-by-point systematic uncertainties.

The flow coefficients $v_{n}$ scaled by $\varepsilon_{n}\left\lbrace 2\right\rbrace$ as a function of $p_{\text{T}}$ at mid-rapidity ($|y| <$ 1) in U+U collisions at $\sqrt{s_{NN}}$ = 193 GeV. The error bars represent statistical uncertainties. The bands represent point-by-point systematic uncertainties.

The flow coefficients $v_{n}$ scaled by $\varepsilon_{n}\left\lbrace 2\right\rbrace$ as a function of $p_{\text{T}}$ at mid-rapidity ($|y| <$ 1) in U+U collisions at $\sqrt{s_{NN}}$ = 193 GeV. The error bars represent statistical uncertainties. The bands represent point-by-point systematic uncertainties.

The flow coefficients $v_{n}$ scaled by $\varepsilon_{n}\left\lbrace 2\right\rbrace$ as a function of $p_{\text{T}}$ at mid-rapidity ($|y| <$ 1) in U+U collisions at $\sqrt{s_{NN}}$ = 193 GeV. The error bars represent statistical uncertainties. The bands represent point-by-point systematic uncertainties.

The flow coefficients $v_{n}$ scaled by $\varepsilon_{n}\left\lbrace 2\right\rbrace$ as a function of $p_{\text{T}}$ at mid-rapidity ($|y| <$ 1) in U+U collisions at $\sqrt{s_{NN}}$ = 193 GeV. The error bars represent statistical uncertainties. The bands represent point-by-point systematic uncertainties.

The flow coefficients $v_{n}$ scaled by $\varepsilon_{n}\left\lbrace 2\right\rbrace$ as a function of $p_{\text{T}}$ at mid-rapidity ($|y| <$ 1) in U+U collisions at $\sqrt{s_{NN}}$ = 193 GeV. The error bars represent statistical uncertainties. The bands represent point-by-point systematic uncertainties.

The flow coefficients $v_{n}$ scaled by $\varepsilon_{n}\left\lbrace 2\right\rbrace$ as a function of $p_{\text{T}}$ at mid-rapidity ($|y| <$ 1) in U+U collisions at $\sqrt{s_{NN}}$ = 193 GeV. The error bars represent statistical uncertainties. The bands represent point-by-point systematic uncertainties.

The flow coefficients $v_{n}$ scaled by $\varepsilon_{n}\left\lbrace 2\right\rbrace$ as a function of $p_{\text{T}}$ at mid-rapidity ($|y| <$ 1) in U+U collisions at $\sqrt{s_{NN}}$ = 193 GeV. The error bars represent statistical uncertainties. The bands represent point-by-point systematic uncertainties.

The flow coefficients $v_{n}$ scaled by $\varepsilon_{n}\left\lbrace 2\right\rbrace$ as a function of $p_{\text{T}}$ at mid-rapidity ($|y| <$ 1) in U+U collisions at $\sqrt{s_{NN}}$ = 193 GeV. The error bars represent statistical uncertainties. The bands represent point-by-point systematic uncertainties.

The flow coefficients $v_{n}$ scaled by $\varepsilon_{n}\left\lbrace 2\right\rbrace$ as a function of $p_{\text{T}}$ at mid-rapidity ($|y| <$ 1) in U+U collisions at $\sqrt{s_{NN}}$ = 193 GeV. The error bars represent statistical uncertainties. The bands represent point-by-point systematic uncertainties.

The flow coefficients $v_{n}$ scaled by $\varepsilon_{n}\left\lbrace 2\right\rbrace$ as a function of $p_{\text{T}}$ at mid-rapidity ($|y| <$ 1) in U+U collisions at $\sqrt{s_{NN}}$ = 193 GeV. The error bars represent statistical uncertainties. The bands represent point-by-point systematic uncertainties.

The flow coefficients $v_{n}$ scaled by $\varepsilon_{n}\left\lbrace 2\right\rbrace$ as a function of $p_{\text{T}}$ at mid-rapidity ($|y| <$ 1) in U+U collisions at $\sqrt{s_{NN}}$ = 193 GeV. The error bars represent statistical uncertainties. The bands represent point-by-point systematic uncertainties.

The flow coefficients $v_{n}$ scaled by $\varepsilon_{n}\left\lbrace 2\right\rbrace$ as a function of $p_{\text{T}}$ at mid-rapidity ($|y| <$ 1) in U+U collisions at $\sqrt{s_{NN}}$ = 193 GeV. The error bars represent statistical uncertainties. The bands represent point-by-point systematic uncertainties.

The flow coefficients $v_{n}$ scaled by $\varepsilon_{n}\left\lbrace 2\right\rbrace$ as a function of $p_{\text{T}}$ at mid-rapidity ($|y| <$ 1) in U+U collisions at $\sqrt{s_{NN}}$ = 193 GeV. The error bars represent statistical uncertainties. The bands represent point-by-point systematic uncertainties.

The flow coefficients $v_{n}$ scaled by $\varepsilon_{n}\left\lbrace 2\right\rbrace$ as a function of $p_{\text{T}}$ at mid-rapidity ($|y| <$ 1) in U+U collisions at $\sqrt{s_{NN}}$ = 193 GeV. The error bars represent statistical uncertainties. The bands represent point-by-point systematic uncertainties.

$v_{n}$ coefficients, scaled by the number of constituent quarks $(n_{q})$ to the power $n/2$ and participant eccentricity $\varepsilon_{n}$, of identified particles versus $(m_{T}-m_{0})/n_{q}$ for three centrality bins in U+U collisions at $\sqrt{s_{NN}}$ = 193 GeV. The error bars represent statistical uncertainties. The bands represent point-by-point systematic uncertainties.

$v_{n}$ coefficients, scaled by the number of constituent quarks $(n_{q})$ to the power $n/2$ and participant eccentricity $\varepsilon_{n}$, of identified particles versus $(m_{T}-m_{0})/n_{q}$ for three centrality bins in U+U collisions at $\sqrt{s_{NN}}$ = 193 GeV. The error bars represent statistical uncertainties. The bands represent point-by-point systematic uncertainties.

$v_{n}$ coefficients, scaled by the number of constituent quarks $(n_{q})$ to the power $n/2$ and participant eccentricity $\varepsilon_{n}$, of identified particles versus $(m_{T}-m_{0})/n_{q}$ for three centrality bins in U+U collisions at $\sqrt{s_{NN}}$ = 193 GeV. The error bars represent statistical uncertainties. The bands represent point-by-point systematic uncertainties.

$v_{n}$ coefficients, scaled by the number of constituent quarks $(n_{q})$ to the power $n/2$ and participant eccentricity $\varepsilon_{n}$, of identified particles versus $(m_{T}-m_{0})/n_{q}$ for three centrality bins in U+U collisions at $\sqrt{s_{NN}}$ = 193 GeV. The error bars represent statistical uncertainties. The bands represent point-by-point systematic uncertainties.

$v_{n}$ coefficients, scaled by the number of constituent quarks $(n_{q})$ to the power $n/2$ and participant eccentricity $\varepsilon_{n}$, of identified particles versus $(m_{T}-m_{0})/n_{q}$ for three centrality bins in U+U collisions at $\sqrt{s_{NN}}$ = 193 GeV. The error bars represent statistical uncertainties. The bands represent point-by-point systematic uncertainties.

$v_{n}$ coefficients, scaled by the number of constituent quarks $(n_{q})$ to the power $n/2$ and participant eccentricity $\varepsilon_{n}$, of identified particles versus $(m_{T}-m_{0})/n_{q}$ for three centrality bins in U+U collisions at $\sqrt{s_{NN}}$ = 193 GeV. The error bars represent statistical uncertainties. The bands represent point-by-point systematic uncertainties.

$v_{n}$ coefficients, scaled by the number of constituent quarks $(n_{q})$ to the power $n/2$ and participant eccentricity $\varepsilon_{n}$, of identified particles versus $(m_{T}-m_{0})/n_{q}$ for three centrality bins in U+U collisions at $\sqrt{s_{NN}}$ = 193 GeV. The error bars represent statistical uncertainties. The bands represent point-by-point systematic uncertainties.

$v_{n}$ coefficients, scaled by the number of constituent quarks $(n_{q})$ to the power $n/2$ and participant eccentricity $\varepsilon_{n}$, of identified particles versus $(m_{T}-m_{0})/n_{q}$ for three centrality bins in U+U collisions at $\sqrt{s_{NN}}$ = 193 GeV. The error bars represent statistical uncertainties. The bands represent point-by-point systematic uncertainties.

$v_{n}$ coefficients, scaled by the number of constituent quarks $(n_{q})$ to the power $n/2$ and participant eccentricity $\varepsilon_{n}$, of identified particles versus $(m_{T}-m_{0})/n_{q}$ for three centrality bins in U+U collisions at $\sqrt{s_{NN}}$ = 193 GeV. The error bars represent statistical uncertainties. The bands represent point-by-point systematic uncertainties.

$v_{n}$ coefficients, scaled by the number of constituent quarks $(n_{q})$ to the power $n/2$ and participant eccentricity $\varepsilon_{n}$, of identified particles versus $(m_{T}-m_{0})/n_{q}$ for three centrality bins in U+U collisions at $\sqrt{s_{NN}}$ = 193 GeV. The error bars represent statistical uncertainties. The bands represent point-by-point systematic uncertainties.

$v_{n}$ coefficients, scaled by the number of constituent quarks $(n_{q})$ to the power $n/2$ and participant eccentricity $\varepsilon_{n}$, of identified particles versus $(m_{T}-m_{0})/n_{q}$ for three centrality bins in U+U collisions at $\sqrt{s_{NN}}$ = 193 GeV. The error bars represent statistical uncertainties. The bands represent point-by-point systematic uncertainties.

$v_{n}$ coefficients, scaled by the number of constituent quarks $(n_{q})$ to the power $n/2$ and participant eccentricity $\varepsilon_{n}$, of identified particles versus $(m_{T}-m_{0})/n_{q}$ for three centrality bins in U+U collisions at $\sqrt{s_{NN}}$ = 193 GeV. The error bars represent statistical uncertainties. The bands represent point-by-point systematic uncertainties.

$v_{n}$ coefficients, scaled by the number of constituent quarks $(n_{q})$ to the power $n/2$ and participant eccentricity $\varepsilon_{n}$, of identified particles versus $(m_{T}-m_{0})/n_{q}$ for three centrality bins in U+U collisions at $\sqrt{s_{NN}}$ = 193 GeV. The error bars represent statistical uncertainties. The bands represent point-by-point systematic uncertainties.

$v_{n}$ coefficients, scaled by the number of constituent quarks $(n_{q})$ to the power $n/2$ and participant eccentricity $\varepsilon_{n}$, of identified particles versus $(m_{T}-m_{0})/n_{q}$ for three centrality bins in U+U collisions at $\sqrt{s_{NN}}$ = 193 GeV. The error bars represent statistical uncertainties. The bands represent point-by-point systematic uncertainties.

$v_{n}$ coefficients, scaled by the number of constituent quarks $(n_{q})$ to the power $n/2$ and participant eccentricity $\varepsilon_{n}$, of identified particles versus $(m_{T}-m_{0})/n_{q}$ for three centrality bins in U+U collisions at $\sqrt{s_{NN}}$ = 193 GeV. The error bars represent statistical uncertainties. The bands represent point-by-point systematic uncertainties.

$v_{n}$ coefficients, scaled by the number of constituent quarks $(n_{q})$ to the power $n/2$ and participant eccentricity $\varepsilon_{n}$, of identified particles versus $(m_{T}-m_{0})/n_{q}$ for three centrality bins in U+U collisions at $\sqrt{s_{NN}}$ = 193 GeV. The error bars represent statistical uncertainties. The bands represent point-by-point systematic uncertainties.

$v_{n}$ coefficients, scaled by the number of constituent quarks $(n_{q})$ to the power $n/2$ and participant eccentricity $\varepsilon_{n}$, of identified particles versus $(m_{T}-m_{0})/n_{q}$ for three centrality bins in U+U collisions at $\sqrt{s_{NN}}$ = 193 GeV. The error bars represent statistical uncertainties. The bands represent point-by-point systematic uncertainties.

$v_{n}$ coefficients, scaled by the number of constituent quarks $(n_{q})$ to the power $n/2$ and participant eccentricity $\varepsilon_{n}$, of identified particles versus $(m_{T}-m_{0})/n_{q}$ for three centrality bins in U+U collisions at $\sqrt{s_{NN}}$ = 193 GeV. The error bars represent statistical uncertainties. The bands represent point-by-point systematic uncertainties.

Ratios of $v_{n}$ coefficients at mid-rapidity ($|y| <$ 1.0) in minimum bias U+U collisions at $\sqrt{s_{NN}}$ = 193 GeV. Error bars represent statistical uncertainties.

Ratios of $v_{n}$ coefficients at mid-rapidity ($|y| <$ 1.0) in minimum bias U+U collisions at $\sqrt{s_{NN}}$ = 193 GeV. Error bars represent statistical uncertainties.

Ratios of $v_{n}$ coefficients at mid-rapidity ($|y| <$ 1.0) in minimum bias U+U collisions at $\sqrt{s_{NN}}$ = 193 GeV. Error bars represent statistical uncertainties.

Ratios of $v_{n}$ coefficients at mid-rapidity ($|y| <$ 1.0) in minimum bias U+U collisions at $\sqrt{s_{NN}}$ = 193 GeV. Error bars represent statistical uncertainties.

Ratios of $v_{n}$ coefficients at mid-rapidity ($|y| <$ 1.0) in minimum bias U+U collisions at $\sqrt{s_{NN}}$ = 193 GeV. Error bars represent statistical uncertainties.

Ratios of $v_{n}$ coefficients at mid-rapidity ($|y| <$ 1.0) in minimum bias U+U collisions at $\sqrt{s_{NN}}$ = 193 GeV. Error bars represent statistical uncertainties.

Ratios of $v_{n}$ coefficients at mid-rapidity ($|y| <$ 1.0) in minimum bias U+U collisions at $\sqrt{s_{NN}}$ = 193 GeV. Error bars represent statistical uncertainties.

Ratios of $v_{n}$ coefficients at mid-rapidity ($|y| <$ 1.0) in minimum bias U+U collisions at $\sqrt{s_{NN}}$ = 193 GeV. Error bars represent statistical uncertainties.

Ratios of $v_{n}$ coefficients at mid-rapidity ($|y| <$ 1.0) in minimum bias U+U collisions at $\sqrt{s_{NN}}$ = 193 GeV. Error bars represent statistical uncertainties.