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.

Parton inclusive-jet $p_T$ and $A_{LL}$ values with associated uncertainties for jet-$\eta$ region $0.5<|\eta|<1$.

Parton inclusive-jet $p_T$ and $A_{LL}$ values with associated uncertainties for jet-$\eta$ region $|\eta|<0.5$.

Parton dijet invariant mass $M_{inv}$ and $A_{LL}$ values with associated uncertainties for the $\textrm{sign}(\eta_1) = \textrm{sign}(\eta_2)$ topology.

Parton dijet invariant mass $M_{inv}$ and $A_{LL}$ values with associated uncertainties for the $\textrm{sign}(\eta_1) \neq \textrm{sign}(\eta_2)$ topology.

Parton inclusive-jet $p_T$, $x_T$, and $A_{LL}$ values with associated uncertainties for jet-$\eta$ region $|\eta|<1$.

The correlation matrix for the point-to-point uncertainties for inclusive jet measurements with jets in the $0.5<|\eta|<1$ region. The $A_{LL}$ uncertainty contribution of $0.0007$ from uncertainty in the relative luminosity measurement and $6.1\%$ from the beam polarization uncertainty, which are common to all the data points, are separated from the listed values.

The correlation matrix for the point-to-point uncertainties for inclusive jet measurements with jets in the $|\eta|<0.5$ region. The $A_{LL}$ uncertainty contribution of $0.0007$ from uncertainty in the relative luminosity measurement and $6.1\%$ from the beam polarization uncertainty, which are common to all the data points, are separated from the listed values.

The correlation matrix for the point-to-point uncertainties for inclusive jet measurements with jets in the $|\eta|<1$ region. The $A_{LL}$ uncertainty contribution of $0.0007$ from uncertainty in the relative luminosity measurement and $6.1\%$ from the beam polarization uncertainty, which are common to all the data points, are separated from the listed values.

The correlation matrix for the point-to-point uncertainties coupling inclusive jet measurements with jets in the $|\eta|<0.5$ region and jets in the $0.5<|\eta|<1$ region. The $A_{LL}$ uncertainty contribution of $0.0007$ from uncertainty in the relative luminosity measurement and $6.1\%$ from the beam polarization uncertainty, which are common to all the data points, are separated from the listed values.

The correlation matrix for the point-to-point uncertainties coupling inclusive jet measurements with jets in the $0.5<|\eta|<1$ region with dijet measurements with the the $\textrm{sign}(\eta_1) = \textrm{sign}(\eta_2)$ topology. The $A_{LL}$ uncertainty contribution of $0.0007$ from uncertainty in the relative luminosity measurement and $6.1\%$ from the beam polarization uncertainty, which are common to all the data points, are separated from the listed values.

The correlation matrix for the point-to-point uncertainties coupling inclusive jet measurements with jets in the $0.5<|\eta|<1$ region with dijet measurements with the the $\textrm{sign}(\eta_1) \neq \textrm{sign}(\eta_2)$ topology. The $A_{LL}$ uncertainty contribution of $0.0007$ from uncertainty in the relative luminosity measurement and $6.1\%$ from the beam polarization uncertainty, which are common to all the data points, are separated from the listed values.

The correlation matrix for the point-to-point uncertainties coupling inclusive jet measurements with jets in the $|\eta|<0.5$ region with dijet measurements with the the $\textrm{sign}(\eta_1) = \textrm{sign}(\eta_2)$ topology. The $A_{LL}$ uncertainty contribution of $0.0007$ from uncertainty in the relative luminosity measurement and $6.1\%$ from the beam polarization uncertainty, which are common to all the data points, are separated from the listed values.

The correlation matrix for the point-to-point uncertainties coupling inclusive jet measurements with jets in the $|\eta|<0.5$ region with dijet measurements with the the $\textrm{sign}(\eta_1) \neq \textrm{sign}(\eta_2)$ topology. The $A_{LL}$ uncertainty contribution of $0.0007$ from uncertainty in the relative luminosity measurement and $6.1\%$ from the beam polarization uncertainty, which are common to all the data points, are separated from the listed values.

The correlation matrix for the point-to-point uncertainties coupling inclusive jet measurements with jets in the $|\eta|<1$ region with dijet measurements with the the $\textrm{sign}(\eta_1) = \textrm{sign}(\eta_2)$ topology. The $A_{LL}$ uncertainty contribution of $0.0007$ from uncertainty in the relative luminosity measurement and $6.1\%$ from the beam polarization uncertainty, which are common to all the data points, are separated from the listed values.

The correlation matrix for the point-to-point uncertainties coupling inclusive jet measurements with jets in the $|\eta|<1$ region with dijet measurements with the the $\textrm{sign}(\eta_1) \neq \textrm{sign}(\eta_2)$ topology. The $A_{LL}$ uncertainty contribution of $0.0007$ from uncertainty in the relative luminosity measurement and $6.1\%$ from the beam polarization uncertainty, which are common to all the data points, are separated from the listed values.

The correlation matrix for the point-to-point uncertainties for dijet measurements with the $\textrm{sign}(\eta_1) = \textrm{sign}(\eta_2)$ topology. The $A_{LL}$ uncertainty contribution of $0.0007$ from uncertainty in the relative luminosity measurement and $6.1\%$ from the beam polarization uncertainty, which are common to all the data points, are separated from the listed values.

The correlation matrix for the point-to-point uncertainties for dijet measurements with the $\textrm{sign}(\eta_1) \neq \textrm{sign}(\eta_2)$ topology. The $A_{LL}$ uncertainty contribution of $0.0007$ from uncertainty in the relative luminosity measurement and $6.1\%$ from the beam polarization uncertainty, which are common to all the data points, are separated from the listed values.

The correlation matrix for the point-to-point uncertainties coupling dijet measurements with the $\textrm{sign}(\eta_1) \neq \textrm{sign}(\eta_2)$ topology and the $\textrm{sign}(\eta_1) = \textrm{sign}(\eta_2)$ topology. The $A_{LL}$ uncertainty contribution of $0.0007$ from uncertainty in the relative luminosity measurement and $6.1\%$ from the beam polarization uncertainty, which are common to all the data points, are separated from the listed values.

$D_s^{\pm}$ invariant yield as a function of $p_{T}$ in 40-80% centrality bin of Au+Au collisions at $\sqrt{s_{_{\rm NN}}}$ = 200 GeV. The $p_T$ bins are 1.5 < $p_T$ < 2.5 GeV/c, 2.5 < $p_T$ < 3.5 GeV/c, 3.5 < $p_T$ < 5.0 GeV/c and 5.0 < $p_T$ < 8.0 GeV/c.

$D_{s}/D^{0}$ yield ratio as a function of $p_{T}$ in 0-10% centrality bin of Au+Au collisions at $\sqrt{s_{_{\rm NN}}}$ = 200 GeV. The $p_T$ bins are 1.5 < $p_T$ < 2.5 GeV/c, 2.5 < $p_T$ < 3.5 GeV/c, 3.5 < $p_T$ < 5.0 GeV/c and 5.0 < $p_T$ < 8.0 GeV/c.

$D_{s}/D^{0}$ yield ratio as a function of $p_{T}$ in 10-40% centrality bin of Au+Au collisions at $\sqrt{s_{_{\rm NN}}}$ = 200 GeV. The $p_T$ bins are 1.0 < $p_T$ < 2.0 GeV/c, 2.0 < $p_T$ < 2.5 GeV/c, 2.5 < $p_T$ < 3.5 GeV/c, 3.5 < $p_T$ < 5.0 GeV/c and 5.0 < $p_T$ < 8.0 GeV/c.

$D_{s}/D^{0}$ yield ratio as a function of $p_{T}$ in 40-80% centrality bin of Au+Au collisions at $\sqrt{s_{_{\rm NN}}}$ = 200 GeV. The $p_T$ bins are 1.5 < $p_T$ < 2.5 GeV/c, 2.5 < $p_T$ < 3.5 GeV/c, 3.5 < $p_T$ < 5.0 GeV/c and 5.0 < $p_T$ < 8.0 GeV/c.

$D_{s}/D^{0}$ yield ratio as a function of $p_{T}$ in 10-20% centrality bin of Au+Au collisions at $\sqrt{s_{_{\rm NN}}}$ = 200 GeV. The $p_T$ bins are 1.5 < $p_T$ < 2.5 GeV/c, 2.5 < $p_T$ < 3.5 GeV/c, 3.5 < $p_T$ < 5.0 GeV/c and 5.0 < $p_T$ < 8.0 GeV/c.

$D_{s}/D^{0}$ yield ratio as a function of $p_{T}$ in 20-40% centrality bin of Au+Au collisions at $\sqrt{s_{_{\rm NN}}}$ = 200 GeV. The $p_T$ bins are 1.5 < $p_T$ < 2.5 GeV/c, 2.5 < $p_T$ < 3.5 GeV/c, 3.5 < $p_T$ < 5.0 GeV/c and 5.0 < $p_T$ < 8.0 GeV/c.

The integrated $D_{s}/D^{0}$ yield ratio within 1.5 < $p_{T}$ < 5.0 GeV/c as a function of collision centrality. The centrality bins are 60-80%, 40-60%, 20-40%, 10-20% and 0-10%.

Centrality dependence of $\Lambda$ global polarization in Au+Au collisions at 200 GeV. Note that the results from previous measurement are rescaled using updated decay parameters (PDG2020), $\alpha_{\Lambda}$=0.732. The original data can be found in <a href="https://www.hepdata.net/record/ins1672785">this page</a>.

Rapidity dependence of directed flow, v1(y), for positive pions with transverse momentum pT > 0.2 GeV/c and total momentum magnitude |p| < 1.6 GeV/c from events within 10-30% centrality. Here, the BBC-based Event Plane method is used. Plotted error bars are statistical only, and systematic errors are of comparable size.

The rapidity dependence of the directed flow for the Lambda using the TPC event plane for 10-30% centrality. Plotted error bars are statistical only, while systematic errors are +/- 0.007.

The rapidity dependence of the directed flow for the K0s using the TPC event plane for 10-30% centrality. Plotted error bars are statistical only, while systematic errors are +/- 0.017.

Beam energy dependence of the directed flow slope dv1=dy at midrapidity for baryons and mesons measured by STAR.

v2 scaled by the number of constituent quarks (nq) for charged pions and protons for 0-30% central collisions. The values of v2 scaled with nq for pions and protons are consistent with each other within errors. For comparison, points from E895 are also shown.

The excitation function v2 for all charged particles or separately for protons and pions, measured by several experiments. The STAR FXT points for protons and for pions are near the region where a change in slope occurs.

Excitation function of Rout, Rside, and Rlong.

Transverse mass dependence of Rout, Rside, and Rlong for three experiments: E895, STAR, and E866. The corresponding center-of-mass collision energies are 4.3 GeV, 4.5 GeV, and 4.9 GeV, respectively. Pairs for the STAR points are created from negative pion tracks in the momentum range 0.15 < pT < 0.8 GeV/c from events in the 0-15% centrality range.

The centrality dependence of Rout, Rside, and Rlong. Errors are statistical only. Here positive and negative pion pairs in the momentum range 0.15 < pT < 0.8 GeV/c are used.

Directed flow slope at mid-rapidity, $dv_{1}/dy|_{y=0}$, as a function of beam energy in 10-40\% Au+Au collisions. Dots represent deuterons. Open squares are the published results in~\cite{Adamczyk_2018} for protons. The band denotes the results for deuterons from AMPT transport model. Statistical uncertainties (bars) and systematic uncertainties (horizontal lines) are shown separately. For visibility, the data points are staggered horizontally.

The $p_{T}$ dependence of $v_1/A$ in $|y|<0.6$ for protons (open squares) in 10-40\% Au+Au collisions at $\sqrt{s_{NN}}$ = 7.7, 11.5, 14.5, 19.6, 27, and 39 GeV. Statistical uncertainties (bars) and systematic uncertainties (horizontal lines) are shown separately.

The $d_{T}$ dependence of $v_1/A$ in $|y|<0.6$ for deuterons (solid circles) in 10-40\% Au+Au collisions at $\sqrt{s_{NN}}$ = 7.7, 11.5, 14.5, 19.6, 27, and 39 GeV. Statistical uncertainties (bars) and systematic uncertainties (horizontal lines) are shown separately.

The $p_{T}$ dependance of $v_{4}^{Inclusive}$, $v_{4}^{Linear}$, $v_{4}^{Mode-coupled}$, $\chi_{4,22}$ and $\rho_{4,22}$

The $p_{T}$ dependance of $v_{5}^{Inclusive}$, $v_{5}^{Linear}$, $v_{5}^{Mode-coupled}$, $\chi_{5,23}$ and $\rho_{5,23}$

$m_{\rm inv}$ dependences of $\Delta\gamma$ in 20-50$\%$ Au+Au collisions at 200 GeV.

The $\pi$ pair $\Delta\gamma$ at $m_{\rm inv} > m_{\rm inv}^{\rm low}$ for the data identified by TPC and more extended $p_T$ range identified by TPC and TOF.

The elliptic flow $v_2$ in bins of $q_2$ for the three centrality bins in the 20-50$\%$ range of Au+Au collisions at 200 GeV.

$m_{\rm inv}$ dependences of the $\Delta\gamma$ from ESE-selected event samples A (50$\%$ largest $q_{2}$) and B (50$\%$ smallest $q_{2}$) in 20-50$\%$ Au+Au collisions at 200 GeV.

$m_{\rm inv}$ dependences of the $\rm \Delta\gamma_{A} - \Delta\gamma_{B}$ in 20-50$\%$ Au+Au collisions at 200 GeV.

$\Delta\gamma_{\rm A}$ versus $\Delta\gamma_{\rm B}$ in 20-50$\%$ Au+Au collisions at 200 GeV. (The data are arranged as function of mass)