The observed and expected $E_T^{miss}$ distribution in the dimuon SR-Z. The negigible estimated contribution from Z+jets is omitted in these distributions. The last bin contains the overflow.

The observed and expected dielectron invariant mass distribution in SR-loose. The last bin contains the overflow.

The observed and expected dimuon invariant mass distribution in SR-loose. The last bin contains the overflow.

The observed and expected dielectron invariant mass distribution in the two-jet $b$-veto SR. The last bin contains the overflow.

The observed and expected dimuon invariant mass distribution in the two-jet $b$-veto SR. The last bin contains the overflow.

The observed and expected dielectron invariant mass distribution in the four jet b-veto SR. The last bin contains the overflow.

The observed and expected dimuon invariant mass distribution in the four-jet $b$-veto SR. The last bin contains the overflow.

The observed and expected dielectron invariant mass distribution in the two-jet $b$-tag SR. The last bin contains the overflow.

The observed and expected dimuon invariant mass distribution in two-jet $b$-tag SR. The last bin contains the overflow.

The observed and expected dielectron invariant mass distribution in the four-jet $b$-tag SR. The last bin contains the overflow.

The observed and expected dimuon invariant mass distribution in the four-jet $b$-tag SR. The last bin contains the overflow.

Expected 95% exclusion contour for the GGM model with $\tan(\beta)=1.5$ in SR-Z.

Observed 95% exclusion contour for the GGM model with $\tan(\beta)=1.5$ in SR-Z.

Expected 95% exclusion contour for the GGM model with $\tan(\beta)=30$ in SR-Z.

Observed 95% exclusion contour for the GGM model with $\tan(\beta)=30$ in SR-Z.

Expected 95% exclusion contour for the two-step first- and second-generation squark simplified model with sleptons in the two-jet $b$-veto SR.

Observed 95% exclusion contour for the two-step first- and second-generation squark simplified model with sleptons in the two-jet $b$-veto SR.

Expected 95% exclusion contour for the two-step gluino simplified model with sleptons in the four-jet $b$-veto SR.

Observed 95% exclusion contour for the two-step gluino simplified model with sleptons in the four-jet $b$-veto SR.

Number of generated events in the two-step gluino simplified model with sleptons.

Production cross-section in the two-step gluino simplified model with sleptons.

Number of generated events in the two-step first- and second-generation squark simplified model with sleptons.

Production cross-section in the two-step first- and second-generation squark simplified model with sleptons.

Number of generated events in the GGM model with $\tan(\beta)=1.5$.

Production cross-section in the GGM model with $\tan(\beta)=1.5$.

Number of generated events in the GGM model with $\tan(\beta)=30$.

Production cross-section in the GGM model with $\tan(\beta)=30$.

Total experimental uncertainty [%] for the two-step gluino simplified model with sleptons.

Total experimental uncertainty [%] for the GGM model with $\tan(\beta)=1.5$.

Total experimental uncertainty for the GGM model with $\tan(\beta)=30$.

Signal acceptance for the GGM model with $\tan(\beta)=1.5$ in the combined electron and muon SR-Z.

Signal acceptance for the GGM model with $\tan(\beta)=30$ in the combined electron and muon SR-Z.

Signal efficiency for the GGM model with $\tan(\beta)=1.5$ in the dielectron SR-Z.

Signal efficiency for the GGM model with $\tan(\beta)=30$ in the dielectron SR-Z.

Signal efficiency for the GGM model with $\tan(\beta)=1.5$ in the dimuon SR-Z.

Signal efficiency for the GGM model with $\tan(\beta)=30$ in the dimuon SR-Z.

Signal efficiency for the GGM model with $\tan(\beta)=1.5$ in the electron and muon combined SR-Z.

Signal efficiency for the GGM model with $\tan(\beta)=30$ in the the electron and muon combined SR-Z.

Signal acceptance for the two-step first- and second-generation squarks simplified model with sleptons in the two-jet $b$-veto SR.

Signal acceptance for the two-step gluino simplified model with sleptons in the four-jet $b$-veto SR.

Signal efficiency for the two-step first- and second-generation squarks simplified model with sleptons in the two-jet $b$-veto SR.

Signal efficiency for the two-step gluino simplified model with sleptons in the four-jet $b$-veto SR.

Upper limits on the signal cross-section at 95% CL for the GGM model with $\tan(\beta)=1.5$.

Observed CL$_{\text{S}}$ for the GGM model with $\tan(\beta)=1.5$.

Expected CL$_{\text{S}}$ for the GGM model with $\tan(\beta)=1.5$.

Upper limits on the signal cross-section at 95% CL for the GGM model with $\tan(\beta)=30$.

Observed CL$_{\text{S}}$ for the GGM model with $\tan(\beta)=30$.

Expected CL$_{\text{S}}$ for the GGM model with $\tan(\beta)=30$.

Upper limits on the signal stength at 95% CL for the two-step first- and second-generation squark simplified model with sleptons. The excluded signal strength is defined as the ratio of the observed excluded production cross section to the expected production cross section calculated at NLO+NLL.

Upper limits on the signal stength at 95% CL for the two-step gluino simplified model with sleptons. The excluded signal strength is defined as the ratio of the observed excluded production cross section to the expected production cross section calculated at NLO+NLL.

Observed CL$_{\text{S}}$ for the two-step first- and second-generation squark simplified model with sleptons.

Expected CL$_{\text{S}}$ for the two-step first- and second-generation squark simplified model with sleptons.

Observed CL$_{\text{S}}$ for the two-step gluino simplified model with sleptons.

Expected CL$_{\text{S}}$ for the two-step gluino simplified model with sleptons.

Cutflow table for three benchmark signal points in SR-Z for the $ee$ and $\mu\mu$ channels separately. The three signal points are taken from the $\tan\beta = 1.5$ grid. 100000 events were generated for each of these points. Shown here are both the unweighted number of events and the number of events normalised to 20.3 fb$^{-1}$. The total experimental systematic uncertainty on the signal yields is indicated at the last cut, along with the corresponding observed and expected $CL_S$ values.

Cutflow table for three benchmark signal points in the two jet b-veto SR of the off-$Z$ search for the $ee$ and $\mu\mu$ channels separately. Shown here are both the unweighted number of events and the number of events normalised to 20.3$^{-1}$. Except for the last two rows indicating the dilepton mass requirements, quoted event yields include all requirements from the top of the table down to the given row.

Cutflow table for three benchmark signal points in the four jet b-veto SR of the off-$Z$ search for the $ee$ and $\mu\mu$ channels separately. Shown here are both the unweighted number of events and the number of events normalised to 20.3 fb$^{-1}$. Except for the last two rows indicating the dilepton mass requirements, quoted event yields include all requirements from the top of the table down to the given row.

Comparison of data with estimated backgrounds in the $\tilde{E}_\text{T}^\text{miss}$ distribution with the TZCR1 event selection except for the requirement on $\tilde{E}_\text{T}^\text{miss}$. The variables $\tilde{E}_\text{T}^\text{miss}$ and $\tilde{m}_\text{T}$ are constructed in the same way as $E_\text{T}^\text{miss}$ and $m_\text{T}$ but treating the leading photon transverse momentum as invisible. The predicted backgrounds are scaled with normalization factors. The uncertainty band includes statistical and all experimental systematic uncertainties. The last bin includes overflow.

Comparison of data with estimated backgrounds in the $\tilde{m}_\text{T}$ distribution with the TZCR1 event selection except for the requirement on $\tilde{m}_\text{T}$. The variables $\tilde{E}_\text{T}^\text{miss}$ and $\tilde{m}_\text{T}$ are constructed in the same way as $E_\text{T}^\text{miss}$ and $m_\text{T}$ but treating the leading photon transverse momentum as invisible. The predicted backgrounds are scaled with normalization factors. The uncertainty band includes statistical and all experimental systematic uncertainties. The last bin includes overflow.

Comparison of the observed data ($n_\text{obs}$) with the predicted background ($n_\text{exp}$) in the validation and signal regions. The background predictions are obtained using the background-only fit configuration. The bottom panel shows the significance of the difference between data and predicted background, where the significance is based on the total uncertainty ($\sigma_\text{tot}$).

Jet multiplicity distributions for events where exactly two signal leptons are selected. No correction factors are included in the background normalizations. The uncertainty band includes statistical and all experimental systematic uncertainties. The last bin includes overflow.

Jet multiplicity distributions for events where exactly one lepton plus one $\tau$ candidate are selected. No correction factors are included in the background normalizations. The uncertainty band includes statistical and all experimental systematic uncertainties. The last bin includes overflow.

The $E_\text{T}^\text{miss}$ distribution in SR1. In the plot, the full event selection in the corresponding signal region is applied, except for the requirement on $E_\text{T}^\text{miss}$. The predicted backgrounds are scaled with normalization factors. The uncertainty band includes statistical and all experimental systematic uncertainties. The last bin contains the overflow. Benchmark signal models are overlaid for comparison. The benchmark models are specified by the gluino and stop masses, given in TeV in the table.

The $m_\text{T}$ distribution in SR1. In the plot, the full event selection in the corresponding signal region is applied, except for the requirement on $m_\text{T}$. The predicted backgrounds are scaled with normalization factors. The uncertainty band includes statistical and all experimental systematic uncertainties. The last bin contains the overflow. Benchmark signal models are overlaid for comparison. The benchmark models are specified by the gluino and stop masses, given in TeV in the table.

Expected (black dashed) 95% excluded regions in the plane of $m_{\tilde{g}}$ versus $m_{\tilde{t}_1}$ for gluino-mediated stop production.

Observed (red solid) 95% excluded regions in the plane of $m_{\tilde{g}}$ versus $m_{\tilde{t}_1}$ for gluino-mediated stop production.

Expected (black dashed) 95% excluded regions in the plane of $m_{\tilde{t}_1}$ versus $m_{\tilde{\chi}_1^0}$ for direct stop production.

Observed (red solid) 95% excluded regions in the plane of $m_{\tilde{t}_1}$ versus $m_{\tilde{\chi}_1^0}$ for direct stop production.

The expected upper limits on $T$ quark pair production times the squared branching ratio for $T \rightarrow tZ$ as a function of the $T$ quark mass.

The observed upper limits on $T$ quark pair production times the squared branching ratio for $T \rightarrow tZ$ as a function of the $T$ quark mass.

The expected limits on $T$ quarks as a function of the branching ratios $B\left(T \rightarrow bW\right)$ and $B\left(T \rightarrow tH\right)$ for a $T$ quark with a mass of 800 GeV. The $T$ is assumed to decay in three possible ways: $T \to tZ$, $T \to tH$, and $T \to bW$.

The observed limits on $T$ quarks as a function of the branching ratios $B\left(T \rightarrow bW\right)$ and $B\left(T \rightarrow tH\right)$ for a $T$ quark with a mass of 800 GeV. The $T$ is assumed to decay in three possible ways: $T \to tZ$, $T \to tH$, and $T \to bW$.

The $m_\text{T}$ distribution in the WVR2-tail validation region which has the same preselection and jet $p_\text{T}$ requirements as SR2.

The $am_\text{T2}$ distribution in the WVR2-tail validation region which has the same preselection and jet $p_\text{T}$ requirements as SR2.

Large-radius jet mass ($R=1.2$), decomposed into the number of small-radius jet constituents. The lower panel shows the ratio of the total data to the total prediction (summed over all jet multiplicities). Events are required to have one lepton, four jets with $p_\text{T}>80,50,40,40$ GeV, at least one $b$-tagged jet, $E_\text{T}^\text{miss}>200$ GeV, and $m_\text{T}>30$ GeV.

Distribution of $m_\text{T2}^\tau$ in data for a selection enriched in $t\bar{t}$ events with one hadronically decaying $\tau$. Events that have no hadronic $\tau$ candidate (that passes the Loose identification criteria, as well as other requirements) are not shown in the plot.

Upper limits on the model cross-section in units of pb for the gluino-mediated stop models.

Upper limits on the model cross-section in units of pb for the models with direct stop pair production.

Illustration of the best expected signal region per signal grid point for the gluino-mediated stop models. This mapping is used for the final combined exclusion limits.

Illustration of the best expected signal region per signal grid point for models with direct stop pair production. This mapping is used for the final combined exclusion limits.

Expected $CL_s$ values for the gluino-mediated stop models.

Observed $CL_s$ values for the gluino-mediated stop models.

Expected $CL_s$ values for the direct stop pair production models.

Observed $CL_s$ values for the direct stop pair production models.

Expected limit using SR1 for models with direct stop pair production and an unpolarized stop (and bino LSP).

Expected limit using SR1 for models with direct stop pair production with $\tilde{t}_1=\tilde{t}_L$ (and bino LSP).

Expected limit using SR1 for models with direct stop pair production with $\tilde{t}_1\sim\tilde{t}_R$ (and bino LSP).

Observed limit using SR1 for models with direct stop pair production and an unpolarized stop (and bino LSP).

Observed limit using SR1 for models with direct stop pair production with $\tilde{t}_1=\tilde{t}_L$ (and bino LSP).

Observed limit using SR1 for models with direct stop pair production with $\tilde{t}_1\sim\tilde{t}_R$ (and bino LSP).

Expected limit using SR2 for models with direct stop pair production and an unpolarized stop (and bino LSP).

Expected limit using SR2 for models with direct stop pair production with $\tilde{t}_1=\tilde{t}_L$ (and bino LSP).

Expected limit using SR2 for models with direct stop pair production with $\tilde{t}_1\sim\tilde{t}_R$ (and bino LSP).

Observed limit using SR2 for models with direct stop pair production and an unpolarized stop (and bino LSP).

Observed limit using SR2 for models with direct stop pair production with $\tilde{t}_1=\tilde{t}_L$ (and bino LSP).

Observed limit using SR2 for models with direct stop pair production with $\tilde{t}_1\sim\tilde{t}_R$ (and bino LSP).

Expected limit using SR1+SR2 (best expected) for models with direct stop pair production and an unpolarized stop (and bino LSP).

Expected limit using SR1+SR2 (best expected) for models with direct stop pair production with $\tilde{t}_1=\tilde{t}_L$ (and bino LSP).

Expected limit using SR1+SR2 (best expected) for models with direct stop pair production with $\tilde{t}_1\sim\tilde{t}_R$ (and bino LSP).

Observed limit using SR1+SR2 (best expected) for models with direct stop pair production and an unpolarized stop (and bino LSP).

Observed limit using SR1+SR2 (best expected) for models with direct stop pair production with $\tilde{t}_1=\tilde{t}_L$ (and bino LSP).

Observed limit using SR1+SR2 (best expected) for models with direct stop pair production with $\tilde{t}_1\sim\tilde{t}_R$ (and bino LSP).

Acceptance for SR1 in the gluino-mediated stop models. The acceptance is defined as the fraction of signal events that pass the analysis selection performed on generator-level objects, therefore emulating an ideal detector with perfect particle identification and no measurement resolution effects.

Acceptance for SR1 in the direct stop pair production. The acceptance is defined as the fraction of signal events that pass the analysis selection performed on generator-level objects, therefore emulating an ideal detector with perfect particle identification and no measurement resolution effects.

Acceptance for SR2 in the gluino-mediated stop models. The acceptance is defined as the fraction of signal events that pass the analysis selection performed on generator-level objects, therefore emulating an ideal detector with perfect particle identification and no measurement resolution effects.

Acceptance for SR2 in the direct stop pair production. The acceptance is defined as the fraction of signal events that pass the analysis selection performed on generator-level objects, therefore emulating an ideal detector with perfect particle identification and no measurement resolution effects.

Acceptance for SR3 in the gluino-mediated stop models. The acceptance is defined as the fraction of signal events that pass the analysis selection performed on generator-level objects, therefore emulating an ideal detector with perfect particle identification and no measurement resolution effects.

Acceptance for SR3 in the direct stop pair production. The acceptance is defined as the fraction of signal events that pass the analysis selection performed on generator-level objects, therefore emulating an ideal detector with perfect particle identification and no measurement resolution effects.

Efficiency for SR1 in the gluino-mediated stop models. The efficiency is the ratio between the expected signal rate calculated with simulated data passing all the reconstruction level cuts applied to reconstructed objects, and the signal rate for an ideal detector (with perfect particle identification and no measurement resolution effects).

Efficiency for SR1 in the direct stop pair production. The efficiency is the ratio between the expected signal rate calculated with simulated data passing all the reconstruction level cuts applied to reconstructed objects, and the signal rate for an ideal detector (with perfect particle identification and no measurement resolution effects).

Efficiency for SR2 in the gluino-mediated stop models. The efficiency is the ratio between the expected signal rate calculated with simulated data passing all the reconstruction level cuts applied to reconstructed objects, and the signal rate for an ideal detector (with perfect particle identification and no measurement resolution effects).

Efficiency for SR2 in the direct stop pair production. The efficiency is the ratio between the expected signal rate calculated with simulated data passing all the reconstruction level cuts applied to reconstructed objects, and the signal rate for an ideal detector (with perfect particle identification and no measurement resolution effects).

Efficiency for SR3 in the gluino-mediated stop models. The efficiency is the ratio between the expected signal rate calculated with simulated data passing all the reconstruction level cuts applied to reconstructed objects, and the signal rate for an ideal detector (with perfect particle identification and no measurement resolution effects).

Efficiency for SR3 in the direct stop pair production. The efficiency is the ratio between the expected signal rate calculated with simulated data passing all the reconstruction level cuts applied to reconstructed objects, and the signal rate for an ideal detector (with perfect particle identification and no measurement resolution effects).

Cross-section upper limit in the $m(\tilde b_1)$-$m(\tilde\chi^0_1)$ plane for the best expected signal region.

Cross-section upper limit in the $m(\tilde b_1)$-$m(\tilde\chi^0_1)$ plane for the SRA250 signal region.

Cross-section upper limit in the $m(\tilde b_1)$-$m(\tilde\chi^0_1)$ plane for the SRA350 signal region.

Cross-section upper limit in the $m(\tilde b_1)$-$m(\tilde\chi^0_1)$ plane for the SRA450 signal region.

Cross-section upper limit in the $m(\tilde b_1)$-$m(\tilde\chi^0_1)$ plane for the SRB signal region.

Expected CLs values in the $m(\tilde b_1)$-$m(\tilde\chi^0_1)$ plane for the best expected signal region.

Expected CLs values in the $m(\tilde b_1)$-$m(\tilde\chi^0_1)$ plane for the SRA250 signal region.

Expected exclusion limit at 95% CL in the $m(\tilde b_1)$-$m(\tilde\chi^0_1)$ plane for the sbottom pair production scenario, for signal region SRA250.

Observed exclusion limit at 95% CL in the $m(\tilde b_1)$-$m(\tilde\chi^0_1)$ plane for the sbottom pair production scenario, for signal region SRA250.

Expected CLs values in the $m(\tilde b_1)$-$m(\tilde\chi^0_1)$ plane for the SRA350 signal region.

Expected exclusion limit at 95% CL in the $m(\tilde b_1)$-$m(\tilde\chi^0_1)$ plane for the sbottom pair production scenario, for signal region SRA350.

Observed exclusion limit at 95% CL in the $m(\tilde b_1)$-$m(\tilde\chi^0_1)$ plane for the sbottom pair production scenario, for signal region SRA350.

Expected CLs values in the $m(\tilde b_1)$-$m(\tilde\chi^0_1)$ plane for the SRA450 signal region.

Expected exclusion limit at 95% CL in the $m(\tilde b_1)$-$m(\tilde\chi^0_1)$ plane for the sbottom pair production scenario, for signal region SRA450.

Observed exclusion limit at 95% CL in the $m(\tilde b_1)$-$m(\tilde\chi^0_1)$ plane for the sbottom pair production scenario, for signal region SRA450.

Expected CLs values in the $m(\tilde b_1)$-$m(\tilde\chi^0_1)$ plane for the SRB signal region.

Expected exclusion limit at 95% CL in the $m(\tilde b_1)$-$m(\tilde\chi^0_1)$ plane for the sbottom pair production scenario, for signal region SRB.

Observed exclusion limit at 95% CL in the $m(\tilde b_1)$-$m(\tilde\chi^0_1)$ plane for the sbottom pair production scenario, for signal region SRB.

Observed CLs values in the $m(\tilde b_1)$-$m(\tilde\chi^0_1)$ plane for the best expected signal region.

Observed CLs values in the $m(\tilde b_1)$-$m(\tilde\chi^0_1)$ plane for the SRA250 signal region.

Observed CLs values in the $m(\tilde b_1)$-$m(\tilde\chi^0_1)$ plane for the SRA350 signal region.

Observed CLs values in the $m(\tilde b_1)$-$m(\tilde\chi^0_1)$ plane for the SRA450 signal region.

Observed CLs values in the $m(\tilde b_1)$-$m(\tilde\chi^0_1)$ plane for the SRB signal region.

Signal efficiency (in %) in the $m(\tilde b_1)$-$m(\tilde\chi^0_1)$ plane for the sbottom pair production scenario, for the best expected signal region.

Signal efficiency (in %) in the $m(\tilde b_1)$-$m(\tilde\chi^0_1)$ plane for the sbottom pair production scenario, for the SRA250 signal region.

Signal efficiency (in %) in the $m(\tilde b_1)$-$m(\tilde\chi^0_1)$ plane for the sbottom pair production scenario, for the SRA350 signal region.

Signal efficiency (in %) in the $m(\tilde b_1)$-$m(\tilde\chi^0_1)$ plane for the sbottom pair production scenario, for the SRA450 signal region.

Signal efficiency (in %) in the $m(\tilde b_1)$-$m(\tilde\chi^0_1)$ plane for the sbottom pair production scenario, for the SRB signal region.

Signal acceptance (in %) in the $m(\tilde b_1)$-$m(\tilde\chi^0_1)$ plane for the sbottom pair production scenario, for the best expected signal region.

Signal acceptance (in %) in the $m(\tilde b_1)$-$m(\tilde\chi^0_1)$ plane for the sbottom pair production scenario, for the SRA250 signal region.

Signal acceptance (in %) in the $m(\tilde b_1)$-$m(\tilde\chi^0_1)$ plane for the sbottom pair production scenario, for the SRA350 signal region.

Signal acceptance (in %) in the $m(\tilde b_1)$-$m(\tilde\chi^0_1)$ plane for the sbottom pair production scenario, for the SRA450 signal region.

Signal acceptance (in %) in the $m(\tilde b_1)$-$m(\tilde\chi^0_1)$ plane for the sbottom pair production scenario, for the SRB signal region.

Total experimental systematic uncertainty in percent on the signal efficiency times acceptance in the $m(\tilde b_1)$-$m(\tilde\chi^0_1)$ plane. The best expected signal region selection is used per point.

Missing transverse momentum distribution after SR3b selection, beside the $E_\mathrm{T}^\mathrm{miss}$ requirement. The results in the signal region correspond to the last inclusive bin. The systematic uncertainties include theory uncertainties for the backgrounds with prompt SS/3L and the full systematic uncertainties for data-driven backgrounds. For illustration the distribution for a benchmark SUSY scenario ($pp\to \tilde g\tilde g$, $\tilde g\to t\bar t\tilde\chi_1^0$, $m_{\tilde g}=1.2$ TeV, $m_{\tilde\chi_1^0}=0.7$ TeV) is also shown.

Observed exclusion limits on the $\tilde g$ and $\tilde\chi_1^0$ masses in the context of SUSY scenarios with simplified mass spectra featuring $\tilde g\tilde g$ pair production with exclusive $\tilde g\to qq(\tilde\ell\ell/\tilde\nu\nu)$ decays. All limits are computed at 95% CL.

Expected exclusion limits on the $\tilde g$ and $\tilde\chi_1^0$ masses in the context of SUSY scenarios with simplified mass spectra featuring $\tilde g\tilde g$ pair production with exclusive $\tilde g\to qq(\tilde\ell\ell/\tilde\nu\nu)$ decays. All limits are computed at 95% CL.

Upper limits on signal cross-sections as function of the $\tilde g$ and $\tilde\chi_1^0$ masses in the context of SUSY scenarios with simplified mass spectra featuring $\tilde g\tilde g$ pair production with exclusive $\tilde g\to qq(\tilde\ell\ell/\tilde\nu\nu)$ decays, obtained using the signal efficiency and acceptance specific to each model. All limits are computed at 95% CL.

Observed exclusion limits on the $\tilde g$ and $\tilde\chi_1^0$ masses in the context of SUSY scenarios with simplified mass spectra featuring $\tilde g\tilde g$ pair production with exclusive $\tilde g\to qqWZ\tilde\chi_1^0$ decays. All limits are computed at 95% CL.

Expected exclusion limits on the $\tilde g$ and $\tilde\chi_1^0$ masses in the context of SUSY scenarios with simplified mass spectra featuring $\tilde g\tilde g$ pair production with exclusive $\tilde g\to qqWZ\tilde\chi_1^0$ decays. All limits are computed at 95% CL.

Upper limits on signal cross-sections as function of the $\tilde g$ and $\tilde\chi_1^0$ masses in the context of SUSY scenarios with simplified mass spectra featuring $\tilde g\tilde g$ pair production with exclusive $\tilde g\to qqWZ\tilde\chi_1^0$ decays, obtained using the signal efficiency and acceptance specific to each model. All limits are computed at 95% CL.

Observed exclusion limits on the $\tilde b_1$ and $\tilde\chi_1^0$ masses in the context of SUSY scenarios with simplified mass spectra featuring $\tilde b_1\tilde b_1^*$ pair production with exclusive $\tilde b_1\to t\tilde\chi_1^-$ decays. All limits are computed at 95% CL.

Expected exclusion limits on the $\tilde b_1$ and $\tilde\chi_1^0$ masses in the context of SUSY scenarios with simplified mass spectra featuring $\tilde b_1\tilde b_1^*$ pair production with exclusive $\tilde b_1\to t\tilde\chi_1^-$ decays. All limits are computed at 95% CL.

Upper limits on signal cross-sections as function of the $\tilde b_1$ and $\tilde\chi_1^0$ masses in the context of SUSY scenarios with simplified mass spectra featuring $\tilde b_1\tilde b_1^*$ pair production with exclusive $\tilde b_1\to t\tilde\chi_1^-$ decays, obtained using the signal efficiency and acceptance specific to each model. All limits are computed at 95% CL.

Observed exclusion limits on the $\tilde g$ and $\tilde\chi_1^0$ masses in the context of SUSY scenarios with simplified mass spectra featuring $\tilde g\tilde g$ pair production with exclusive $\tilde g\to t\bar t\tilde\chi_1^0$ decays. All limits are computed at 95% CL.

Expected exclusion limits on the $\tilde g$ and $\tilde\chi_1^0$ masses in the context of SUSY scenarios with simplified mass spectra featuring $\tilde g\tilde g$ pair production with exclusive $\tilde g\to t\bar t\tilde\chi_1^0$ decays. All limits are computed at 95% CL.

Upper limits on signal cross-sections as function of the $\tilde g$ and $\tilde\chi_1^0$ masses in the context of SUSY scenarios with simplified mass spectra featuring $\tilde g\tilde g$ pair production with exclusive $\tilde g\to t\bar t\tilde\chi_1^0$ decays, obtained using the signal efficiency and acceptance specific to each model. All limits are computed at 95% CL.

SUSY scenario with $\tilde g\tilde g$ production and $\tilde g\to q\bar q(\tilde\ell\ell/\tilde\nu\nu)$ decay: signal acceptance (in %) in the signal region SR0b3j. The benchmark scenarios used to set exclusion limits are materialized by black dot markers. Acceptance and efficiency are defined as in appendix A of [JHEP 06 (2014) 124, arXiv: 1403.4853v1 [hep-ex]].

SUSY scenario with $\tilde g\tilde g$ production and $\tilde g\to q\bar q(\tilde\ell\ell/\tilde\nu\nu)$ decay: reconstruction efficiency (in %) in the signal region SR0b3j. The benchmark scenarios used to set exclusion limits are materialized by black dot markers. Acceptance and efficiency are defined as in appendix A of [JHEP 06 (2014) 124, arXiv: 1403.4853v1 [hep-ex]].

SUSY scenario with $\tilde g\tilde g$ production and $\tilde g\to q\bar qWZ\tilde\chi_1^0$ decay: signal acceptance (in %) in the signal region SR0b5j. The benchmark scenarios used to set exclusion limits are materialized by black dot markers. Acceptance and efficiency are defined as in appendix A of [JHEP 06 (2014) 124, arXiv: 1403.4853v1 [hep-ex]].

SUSY scenario with $\tilde g\tilde g$ production and $\tilde g\to q\bar qWZ\tilde\chi_1^0$ decay: reconstruction efficiency (in %) in the signal region SR0b5j. The benchmark scenarios used to set exclusion limits are materialized by black dot markers. Acceptance and efficiency are defined as in appendix A of [JHEP 06 (2014) 124, arXiv: 1403.4853v1 [hep-ex]].

SUSY scenario with $\tilde b_1\tilde b_1^*$ production and $\tilde b_1\to tW\tilde\chi_1^0$ decay: signal acceptance (in %) in the signal region SR1b. The benchmark scenarios used to set exclusion limits are materialized by black dot markers. Acceptance and efficiency are defined as in appendix A of [JHEP 06 (2014) 124, arXiv: 1403.4853v1 [hep-ex]].

SUSY scenario with $\tilde b_1\tilde b_1^*$ production and $\tilde b_1\to tW\tilde\chi_1^0$ decay: reconstruction efficiency (in %) in the signal region SR1b. The benchmark scenarios used to set exclusion limits are materialized by black dot markers. Acceptance and efficiency are defined as in appendix A of [JHEP 06 (2014) 124, arXiv: 1403.4853v1 [hep-ex]].

SUSY scenario with $\tilde g\tilde g$ production and $\tilde g\to t\bar t\tilde\chi_1^0$ decay: signal acceptance (in %) in the signal region SR3b. The benchmark scenarios used to set exclusion limits are materialized by black dot markers. Acceptance and efficiency are defined as in appendix A of [JHEP 06 (2014) 124, arXiv: 1403.4853v1 [hep-ex]].

SUSY scenario with $\tilde g\tilde g$ production and $\tilde g\to t\bar t\tilde\chi_1^0$ decay: reconstruction efficiency (in %) in the signal region SR3b. The benchmark scenarios used to set exclusion limits are materialized by black dot markers. Acceptance and efficiency are defined as in appendix A of [JHEP 06 (2014) 124, arXiv: 1403.4853v1 [hep-ex]].

Data (bold boxes) and background estimates (colour fill) for m_beta vs. m_betagamma for the gluino R-hadron search (1000 GeV). The blue thin-line boxes illustrate the expected signal (on top of background) for the given R-hadron mass hypothesis. The black dashed vertical/horizontal lines at 500 GeV show the mass selection (signal region in the top-right). Two events pass this selection.

Expected (dashed black line) and observed (solid red line) 95% CL upper limits on the cross section as a function of mass for the production of long-lived gluino R-hadrons. The theory prediction along with its +-1sigma uncertainty is show as a black line and a blue band, respectively. The observed 8 TeV Run-1 limit and theory prediction [arXiv:1411.6795] are shown in dash-dotted and dotted lines, respectively.

Expected (dashed black line) and observed (solid red line) 95% CL upper limits on the cross section as a function of mass for the production of bottom-squark R-hadrons. The theory prediction along with its +-1sigma uncertainty is show as a black line and a blue band, respectively. The observed 8 TeV Run-1 limit and theory prediction [arXiv:1411.6795] are shown in dash-dotted and dotted lines, respectively.

Expected (dashed black line) and observed (solid red line) 95% CL upper limits on the cross section as a function of mass for the production of top-squark R-hadrons. The theory prediction along with its +-1sigma uncertainty is show as a black line and a blue band, respectively. The observed 8 TeV Run-1 limit and theory prediction [arXiv:1411.6795] are shown in dash-dotted and dotted lines, respectively.

Final selection requirements as a function of the simulated R-hadron mass.

Summary of all studied systematic uncertainties. Ranges indicate a dependency on the R-hadron mass hypothesis (from low to high masses).

Expected signal yield (Nsig) and efficiency (eff.), estimated background (Nbkg) and observed number of events in data (Nobs) for the full mass range after the final selection using 3.2/fb of data. The stated uncertainties include both the statistical and systematic contribution.

Distribution of the truth-level beta for gluino R-hadrons in exemplary signal MC samples and muons in a Zmumu MC sample. All distributions have been normalised to one. The last bin contains the overflow of the histograms. The distributions illustrate the good discriminating power of the variables.

Distribution of the truth-level betagamma for gluino R-hadrons in exemplary signal MC samples and muons in a Zmumu MC sample. All distributions have been normalised to one. The last bin contains the overflow of the histograms. The distributions illustrate the good discriminating power of the variables.

Expected (dashed black line) and observed (solid red line) 95% confidence level upper limits on the cross section as a function of mass for the production of long-lived gluino R-hadrons. The theory prediction along with its +-1sigma uncertainty is show as a black line and a blue band, respectively. For meta-stable gluinos with a lifetime of 50 ns. (mass exclusion: about 1660 GeV expected, 1520 GeV observed).

Expected (dashed black line) and observed (solid red line) 95% confidence level upper limits on the cross section as a function of mass for the production of long-lived gluino R-hadrons. The theory prediction along with its +-1sigma uncertainty is show as a black line and a blue band, respectively. For meta-stable gluinos with a lifetime of 30 ns. (mass exclusion: about 1660 GeV expected, 1520 GeV observed).

Expected (dashed black line) and observed (solid red line) 95% confidence level upper limits on the cross section as a function of mass for the production of long-lived gluino R-hadrons. The theory prediction along with its +-1sigma uncertainty is show as a black line and a blue band, respectively. For meta-stable gluinos with a lifetime of 10 ns. (mass exclusion: about 1660 GeV expected, 1520 GeV observed).

Object-quality selection cut-flow with observed data and exemplary expected events (scaled to 3.2/fb for MC) in the gluino R-hadron search.

Object-quality selection cut-flow with observed data and exemplary expected events (scaled to 3.2/fb for MC) in the squark R-hadron search.

Expected signal yield (Nsig) and efficiency (eff.), estimated background (Nbkg) and observed number of events in data (Nobs) for the full mass range in the meta-stable gluino R-hadron search using 3.2/fb of data. The stated uncertainties include both the statistical and systematic contribution.

The distribution of the missing transverse momentum is shown in soft-lepton 2-jet W+jets control regions after normalising the ttbar and W+jets background processes in the simultaneous fit.

Expected background yields as obtained in the background-only fits in all hard-lepton and soft-lepton validation together with observed data are given. Uncertainties in the fitted background estimates combine statistical (in the simulated event yields) and systematic uncertainties.

Expected background yields as obtained in the background-only fits in all hard-lepton and soft-lepton signal together with observed data are given. Uncertainties in the fitted background estimates combine statistical (in the simulated event yields) and systematic uncertainties.

Distributions of mt for the hard-lepton 4-jet low-x signal region. The requirement on the variable plotted is removed from the definitions of the signal regions, where the arrow indicates the position of the cut in the signal region. The lower panels of the plots show the ratio of the observed data to the total background prediction as derived in the background-only fit. The uncertainty bands plotted include all statistical and systematic uncertainties as discussed in Section 7. The component `Others' is the sum of Z+jets and ttbar+V. The last bin includes the overflow.

Distributions of met/meff for the 4-jet high-x signal region. The requirement on the variable plotted is removed from the definitions of the signal regions, where the arrow indicates the position of the cut in the signal region. The lower panels of the plots show the ratio of the observed data to the total background prediction as derived in the background-only fit. The uncertainty bands plotted include all statistical and systematic uncertainties as discussed in Section 7. The component `Others' is the sum of Z+jets and ttbar+V. The last bin includes the overflow.

Distributions of mt for the hard-lepton 5-jet signal region. The requirement on the variable plotted is removed from the definitions of the signal regions, where the arrow indicates the position of the cut in the signal region. The lower panels of the plots show the ratio of the observed data to the total background prediction as derived in the background-only fit. The uncertainty bands plotted include all statistical and systematic uncertainties as discussed in Section 7. The component `Others' is the sum of Z+jets and ttbar+V. The last bin includes the overflow.

Distributions of mt for the hard-lepton 6-jet signal region. The requirement on the variable plotted is removed from the definitions of the signal regions, where the arrow indicates the position of the cut in the signal region. The lower panels of the plots show the ratio of the observed data to the total background prediction as derived in the background-only fit. The uncertainty bands plotted include all statistical and systematic uncertainties as discussed in Section 7. The component `Others' is the sum of Z+jets and ttbar+V. The last bin includes the overflow.

Distributions of met for the soft-lepton 2-jet signal region. The requirement on the variable plotted is removed from the definitions of the signal regions, where the arrow indicates the position of the cut in the signal region. The lower panels of the plots show the ratio of the observed data to the total background prediction as derived in the background-only fit. The uncertainty bands plotted include all statistical and systematic uncertainties as discussed in Section 7. The component `Others' is the sum of Z+jets and ttbar+V. The last bin includes the overflow.

Distributions of met for the soft-lepton 5-jet signal region. The requirement on the variable plotted is removed from the definitions of the signal regions, where the arrow indicates the position of the cut in the signal region. The lower panels of the plots show the ratio of the observed data to the total background prediction as derived in the background-only fit. The uncertainty bands plotted include all statistical and systematic uncertainties as discussed in Section 7. The component `Others' is the sum of Z+jets and ttbar+V. The last bin includes the overflow.

The observed combined 95% CL exclusion limits in the the gluino simplified models using for each model point the signal region with the best expected sensitivity. The limits are presented in the (gluino, chargino) mass plane for the scenario where the mass of the chargino is fixed to $x=(m(\tilde\chi^\pm_1)-m(\tilde\chi^0_1))/(m(\tilde g) - m(\tilde\chi^0_1)) = 1/2$ models.

The expected combined 95% CL exclusion limits in the the gluino simplified models using for each model point the signal region with the best expected sensitivity. The limits are presented in the (gluino, chargino) mass plane for the scenario where the mass of the chargino is fixed to $x=(m(\tilde\chi^\pm_1)-m(\tilde\chi^0_1))/(m(\tilde g) - m(\tilde\chi^0_1)) = 1/2$ models.

The yellow band ($+ 1 \sigma$) of the combined 95% CL exclusion limits in the the gluino simplified models using for each model point the signal region with the best expected sensitivity. The limits are presented in the (gluino, chargino) mass plane for the scenario where the mass of the chargino is fixed to $x=(m(\tilde\chi^\pm_1)-m(\tilde\chi^0_1))/(m(\tilde g) - m(\tilde\chi^0_1)) = 1/2$ models. The yellow band represents the $\pm 1 \sigma$ variation of the median expected limit due to the experimental and theoretical uncertainties.

The yellow band ($- 1 \sigma$) of the combined 95% CL exclusion limits in the the gluino simplified models using for each model point the signal region with the best expected sensitivity. The limits are presented in the (gluino, chargino) mass plane for the scenario where the mass of the chargino is fixed to $x=(m(\tilde\chi^\pm_1)-m(\tilde\chi^0_1))/(m(\tilde g) - m(\tilde\chi^0_1)) = 1/2$ models. The yellow band represents the $\pm 1 \sigma$ variation of the median expected limit due to the experimental and theoretical uncertainties.

The observed combined 95% CL exclusion limits in the the gluino simplified models using for each model point the signal region with the best expected sensitivity. The limits are presented in the (gluino, x) plane for the chargino = 60 GeV models where $x=(m(\tilde\chi^\pm_1)-m(\tilde\chi^0_1))/(m(\tilde g) - m(\tilde\chi^0_1))$.

The expected combined 95% CL exclusion limits in the the gluino simplified models using for each model point the signal region with the best expected sensitivity. The limits are presented in the (gluino, x) plane for the chargino = 60 GeV models where $x=(m(\tilde\chi^\pm_1)-m(\tilde\chi^0_1))/(m(\tilde g) - m(\tilde\chi^0_1))$.

The yellow band ($+ 1 \sigma$) of the combined 95% CL exclusion limits in the the gluino simplified models using for each model point the signal region with the best expected sensitivity. The limits are presented in the (gluino, x) plane for the chargino = 60 GeV models where $x=(m(\tilde\chi^\pm_1)-m(\tilde\chi^0_1))/(m(\tilde g) - m(\tilde\chi^0_1))$. The yellow band represents the $\pm 1 \sigma$ variation of the median expected limit due to the experimental and theoretical uncertainties.

The yellow band ($- 1 \sigma$) of the combined 95% CL exclusion limits in the the gluino simplified models using for each model point the signal region with the best expected sensitivity. The limits are presented in the (gluino, x) plane for the chargino = 60 GeV models where $x=(m(\tilde\chi^\pm_1)-m(\tilde\chi^0_1))/(m(\tilde g) - m(\tilde\chi^0_1))$. The yellow band represents the $\pm 1 \sigma$ variation of the median expected limit due to the experimental and theoretical uncertainties.

The observed limits for the soft-lepton 2-jet signal region. The limits are presented in the (gluino, chargino) mass plane for the scenario where the mass of the chargino is fixed to $x=(m(\tilde\chi^\pm_1)-m(\tilde\chi^0_1))/(m(\tilde g) - m(\tilde\chi^0_1)) = 1/2$ models.

The expected limits for the soft-lepton 2-jet signal region. The limits are presented in the (gluino, chargino) mass plane for the scenario where the mass of the chargino is fixed to $x=(m(\tilde\chi^\pm_1)-m(\tilde\chi^0_1))/(m(\tilde g) - m(\tilde\chi^0_1)) = 1/2$ models.

The observed limits for the hard-lepton 4-jet low-x signal region. The limits are presented in the (gluino, chargino) mass plane for the scenario where the mass of the chargino is fixed to $x=(m(\tilde\chi^\pm_1)-m(\tilde\chi^0_1))/(m(\tilde g) - m(\tilde\chi^0_1)) = 1/2$ models.

The expected limits for the hard-lepton 4-jet low-x signal region. The limits are presented in the (gluino, chargino) mass plane for the scenario where the mass of the chargino is fixed to $x=(m(\tilde\chi^\pm_1)-m(\tilde\chi^0_1))/(m(\tilde g) - m(\tilde\chi^0_1)) = 1/2$ models.

The observed limits for the hard-lepton 5-jet signal region. The limits are presented in the (gluino, chargino) mass plane for the scenario where the mass of the chargino is fixed to $x=(m(\tilde\chi^\pm_1)-m(\tilde\chi^0_1))/(m(\tilde g) - m(\tilde\chi^0_1)) = 1/2$ models.

The expected limits for the hard-lepton 5-jet signal region. The limits are presented in the (gluino, chargino) mass plane for the scenario where the mass of the chargino is fixed to $x=(m(\tilde\chi^\pm_1)-m(\tilde\chi^0_1))/(m(\tilde g) - m(\tilde\chi^0_1)) = 1/2$ models.

The observed limits for the hard-lepton 6-jet signal region. The limits are presented in the (gluino, chargino) mass plane for the scenario where the mass of the chargino is fixed to $x=(m(\tilde\chi^\pm_1)-m(\tilde\chi^0_1))/(m(\tilde g) - m(\tilde\chi^0_1)) = 1/2$ models.

The expected limits for the hard-lepton 6-jet signal region. The limits are presented in the (gluino, chargino) mass plane for the scenario where the mass of the chargino is fixed to $x=(m(\tilde\chi^\pm_1)-m(\tilde\chi^0_1))/(m(\tilde g) - m(\tilde\chi^0_1)) = 1/2$ models.

The observed limits for the soft-lepton 5-jet signal region. The limits are presented in the (gluino, x) plane for the chargino = 60 GeV models where $x=(m(\tilde\chi^\pm_1)-m(\tilde\chi^0_1))/(m(\tilde g) - m(\tilde\chi^0_1))$.

The expected limits for the soft-lepton 5-jet signal region. The limits are presented in the (gluino, x) plane for the chargino = 60 GeV models where $x=(m(\tilde\chi^\pm_1)-m(\tilde\chi^0_1))/(m(\tilde g) - m(\tilde\chi^0_1))$.

The observed limits for the hard-lepton 4-jet low-x signal region. The limits are presented in the (gluino, x) plane for the chargino = 60 GeV models where $x=(m(\tilde\chi^\pm_1)-m(\tilde\chi^0_1))/(m(\tilde g) - m(\tilde\chi^0_1))$.

The expected limits for the hard-lepton 4-jet low-x signal region. The limits are presented in the (gluino, x) plane for the chargino = 60 GeV models where $x=(m(\tilde\chi^\pm_1)-m(\tilde\chi^0_1))/(m(\tilde g) - m(\tilde\chi^0_1))$.

The observed limits for the hard-lepton 4-jet high-x signal region. The limits are presented in the (gluino, x) plane for the chargino = 60 GeV models where $x=(m(\tilde\chi^\pm_1)-m(\tilde\chi^0_1))/(m(\tilde g) - m(\tilde\chi^0_1))$.

The expected limits for the hard-lepton 4-jet high-x signal region. The limits are presented in the (gluino, x) plane for the chargino = 60 GeV models where $x=(m(\tilde\chi^\pm_1)-m(\tilde\chi^0_1))/(m(\tilde g) - m(\tilde\chi^0_1))$.

The observed limits for the hard-lepton 5-jet signal region. The limits are presented in the (gluino, x) plane for the chargino = 60 GeV models where $x=(m(\tilde\chi^\pm_1)-m(\tilde\chi^0_1))/(m(\tilde g) - m(\tilde\chi^0_1))$.

The expected limits for the hard-lepton 5-jet signal region. The limits are presented in the (gluino, x) plane for the chargino = 60 GeV models where $x=(m(\tilde\chi^\pm_1)-m(\tilde\chi^0_1))/(m(\tilde g) - m(\tilde\chi^0_1))$.

The observed limits for the hard-lepton 6-jet signal region. The limits are presented in the (gluino, x) plane for the chargino = 60 GeV models where $x=(m(\tilde\chi^\pm_1)-m(\tilde\chi^0_1))/(m(\tilde g) - m(\tilde\chi^0_1))$.

The expected limits for the hard-lepton 6-jet signal region. The limits are presented in the (gluino, x) plane for the chargino = 60 GeV models where $x=(m(\tilde\chi^\pm_1)-m(\tilde\chi^0_1))/(m(\tilde g) - m(\tilde\chi^0_1))$.

Number of generated events in the (gluino, chargino) mass plane for the scenario where the mass of the chargino is fixed to $x=(m(\tilde\chi^\pm_1)-m(\tilde\chi^0_1))/(m(\tilde g) - m(\tilde\chi^0_1)) = 1/2$.

Number of generated events in the (gluino, x) plane for the chargino = 60 GeV models.

Production cross-section in the (gluino, chargino) mass plane for the scenario where the mass of the chargino is fixed to $x=(m(\tilde\chi^\pm_1)-m(\tilde\chi^0_1))/(m(\tilde g) - m(\tilde\chi^0_1)) = 1/2$.

Production cross-section in the (gluino, x) plane for the chargino = 60 GeV models.

Acceptance times efficiency obtained in the different signal regions in the (gluino, chargino) mass plane for the scenario where the mass of the chargino is fixed to $x=(m(\tilde\chi^\pm_1)-m(\tilde\chi^0_1))/(m(\tilde g) - m(\tilde\chi^0_1)) = 1/2$ (hard-lepton 4-jet low-x).

Acceptance times efficiency in the different signal regions in the (gluino, chargino) mass plane for the scenario where the mass of the chargino is fixed to $x=(m(\tilde\chi^\pm_1)-m(\tilde\chi^0_1))/(m(\tilde g) - m(\tilde\chi^0_1)) = 1/2$ (hard-lepton 4-jet high-x).

Acceptance times efficiency in the different signal regions in the (gluino, chargino) mass plane for the scenario where the mass of the chargino is fixed to $x=(m(\tilde\chi^\pm_1)-m(\tilde\chi^0_1))/(m(\tilde g) - m(\tilde\chi^0_1)) = 1/2$ (hard-lepton 5-jet).

Acceptance times efficiency in the different signal regions in the (gluino, chargino) mass plane for the scenario where the mass of the chargino is fixed to $x=(m(\tilde\chi^\pm_1)-m(\tilde\chi^0_1))/(m(\tilde g) - m(\tilde\chi^0_1)) = 1/2$ (hard-lepton 6-jet).

Acceptance times efficiency in the different signal regions in the (gluino, chargino) mass plane for the scenario where the mass of the chargino is fixed to $x=(m(\tilde\chi^\pm_1)-m(\tilde\chi^0_1))/(m(\tilde g) - m(\tilde\chi^0_1)) = 1/2$ (soft-lepton 2-jet).

Acceptance times efficiency obtained in the different signal regions in the (gluino, x) plane for the chargino = 60 GeV models (hard-lepton 4-jet low-x).

Acceptance times efficiency in the different signal regions in the (gluino, x) plane for the chargino = 60 GeV models (hard-lepton 4-jet high-x).

Acceptance times efficiency in the different signal regions in the (gluino, x) plane for the chargino = 60 GeV models (hard-lepton 5-jet).

Acceptance times efficiency in the different signal regions in the (gluino, x) plane for the chargino = 60 GeV models (hard-lepton 6-jet).

Acceptance times efficiency in the different signal regions in the (gluino, x) plane for the chargino = 60 GeV models (soft-lepton 5-jet).

The observed CLs values as obtained in the different signal regions in the (gluino, chargino) mass plane for the scenario where the mass of the chargino is fixed to $x=(m(\tilde\chi^\pm_1)-m(\tilde\chi^0_1))/(m(\tilde g) - m(\tilde\chi^0_1)) = 1/2$ (hard-lepton 4-jet low-x).

The observed CLs values as obtained in the different signal regions in the (gluino, chargino) mass plane for the scenario where the mass of the chargino is fixed to $x=(m(\tilde\chi^\pm_1)-m(\tilde\chi^0_1))/(m(\tilde g) - m(\tilde\chi^0_1)) = 1/2$ (hard-lepton 4-jet high-x).

The observed CLs values as obtained in the different signal regions in the (gluino, chargino) mass plane for the scenario where the mass of the chargino is fixed to $x=(m(\tilde\chi^\pm_1)-m(\tilde\chi^0_1))/(m(\tilde g) - m(\tilde\chi^0_1)) = 1/2$ (hard-lepton 5-jet).

The observed CLs values as obtained in the different signal regions in the (gluino, chargino) mass plane for the scenario where the mass of the chargino is fixed to $x=(m(\tilde\chi^\pm_1)-m(\tilde\chi^0_1))/(m(\tilde g) - m(\tilde\chi^0_1)) = 1/2$ (hard-lepton 6-jet).

The observed CLs values as obtained in the different signal regions in the (gluino, chargino) mass plane for the scenario where the mass of the chargino is fixed to $x=(m(\tilde\chi^\pm_1)-m(\tilde\chi^0_1))/(m(\tilde g) - m(\tilde\chi^0_1)) = 1/2$ (soft-lepton 2-jet).

The expected CLs values as obtained in the different signal regions in the (gluino, chargino) mass plane for the scenario where the mass of the chargino is fixed to $x=(m(\tilde\chi^\pm_1)-m(\tilde\chi^0_1))/(m(\tilde g) - m(\tilde\chi^0_1)) = 1/2$ (hard-lepton 4-jet low-x).

The expected CLs values as obtained in the different signal regions in the (gluino, chargino) mass plane for the scenario where the mass of the chargino is fixed to $x=(m(\tilde\chi^\pm_1)-m(\tilde\chi^0_1))/(m(\tilde g) - m(\tilde\chi^0_1)) = 1/2$ (hard-lepton 4-jet high-x).

The expected CLs values as obtained in the different signal regions in the (gluino, chargino) mass plane for the scenario where the mass of the chargino is fixed to $x=(m(\tilde\chi^\pm_1)-m(\tilde\chi^0_1))/(m(\tilde g) - m(\tilde\chi^0_1)) = 1/2$ (hard-lepton 5-jet).

The expected CLs values as obtained in the different signal regions in the (gluino, chargino) mass plane for the scenario where the mass of the chargino is fixed to $x=(m(\tilde\chi^\pm_1)-m(\tilde\chi^0_1))/(m(\tilde g) - m(\tilde\chi^0_1)) = 1/2$ (hard-lepton 6-jet).

The expected CLs values as obtained in the different signal regions in the (gluino, chargino) mass plane for the scenario where the mass of the chargino is fixed to $x=(m(\tilde\chi^\pm_1)-m(\tilde\chi^0_1))/(m(\tilde g) - m(\tilde\chi^0_1)) = 1/2$ (soft-lepton 2-jet).

The observed CLs values as obtained in the different signal regions in the (gluino, x) plane for the chargino = 60 GeV models (hard-lepton 4-jet low-x).

The observed CLs values as obtained in the different signal regions in the (gluino, x) plane for the chargino = 60 GeV models (hard-lepton 4-jet high-x).

The observed CLs values as obtained in the different signal regions in the (gluino, x) plane for the chargino = 60 GeV models (hard-lepton 5-jet).

The observed CLs values as obtained in the different signal regions in the (gluino, x) plane for the chargino = 60 GeV models (hard-lepton 6-jet).

The observed CLs values as obtained in the different signal regions in the (gluino, x) plane for the chargino = 60 GeV models (soft-lepton 5-jet).

The expected CLs values as obtained in the different signal regions in the (gluino, x) plane for the chargino = 60 GeV models (hard-lepton 4-jet low-x).

The expected CLs values as obtained in the different signal regions in the (gluino, x) plane for the chargino = 60 GeV models (hard-lepton 4-jet high-x).

The expected CLs values as obtained in the different signal regions in the (gluino, x) plane for the chargino = 60 GeV models (hard-lepton 5-jet).

The expected CLs values as obtained in the different signal regions in the (gluino, x) plane for the chargino = 60 GeV models (hard-lepton 6-jet).

The expected CLs values as obtained in the different signal regions in the (gluino, x) plane for the chargino = 60 GeV models (soft-lepton 5-jet).

The signal region yielding in the best expected limit is indicated for every signal point used in the the gluino simplified models for the (gluino, chargino) mass plane for the scenario where the mass of the chargino is fixed to $x=(m(\tilde\chi^\pm_1)-m(\tilde\chi^0_1))/(m(\tilde g) - m(\tilde\chi^0_1)) = 1/2$.

The signal region yielding in the best expected limit is indicated for every signal point used in the the gluino simplified models for the (gluino, x) mass plane where for the chargino = 60 GeV and $x=(m(\tilde\chi^\pm_1)-m(\tilde\chi^0_1))/(m(\tilde g) - m(\tilde\chi^0_1))$.

Model-dependent 95% CL upper limits on the visible cross-section in addition to the observed and expected exclusion limits for the (gluino, chargino) mass plane for the scenario where the mass of the chargino is fixed to $x=(m(\tilde\chi^\pm_1)-m(\tilde\chi^0_1))/(m(\tilde g) - m(\tilde\chi^0_1)) = 1/2$.

Model-dependent 95% CL upper limits on the visible cross-section in addition to the observed and expected exclusion limits for the (gluino, x) mass plane where for the chargino = 60 GeV and $x=(m(\tilde\chi^\pm_1)-m(\tilde\chi^0_1))/(m(\tilde g) - m(\tilde\chi^0_1))$.

Simulated background event samples: the corresponding generator, parton shower, cross-section normalisation, PDF set and underlying-event tune are shown.

Overview of the selection criteria for the soft-lepton signal regions. The symbol $p_{T}^{l}$ refers to signal leptons.

Overview of the selection criteria for the hard-lepton signal regions. The symbol $p_{T}^{l}$ refers to signal leptons.

Background fit results for the hard-lepton and soft-lepton signal regions, for an integrated luminosity of 3.2 fb-1. Uncertainties in the fitted background estimates combine statistical (in the simulated event yields) and systematic uncertainties. The uncertainties in this table are symmetrised for propagation purposes but truncated at zero to remain within the physical boundaries.

Breakdown of upper limits. The columns show from left to right: the name of the respective signal region; the 95% confidence level (CL) upper limits on the visible cross-section and on the number of signal events the 95% CL upper limit on the number of signal events, given the expected number (and $\pm 1 \sigma$ variations on the expectation) of background events; the two-sided CLb value, i.e. the confidence level observed for the background-only hypothesis and the one-sided discovery p-value (p(s = 0)). The discovery p-values are capped to 0.5 in the case of observing less events than the fitted background estimates.

Table shows the data, fitted background and expected signal event counts for a benchmark signal point in each bin of the mt distribution shown in figure 5 (top left). The fit results are shown for an integrated luminosity of 3.2 fb-1.

Table shows the data, fitted background and expected signal event counts for a benchmark signal point in each bin of the met/meff distribution shown in figure 5 (top right). The fit results are shown for an integrated luminosity of 3.2 fb-1.

Table shows the data, fitted background and expected signal event counts for a benchmark signal point in each bin of the mt distribution shown in figure 5 (middle left). The fit results are shown for an integrated luminosity of 3.2 fb-1.

Table shows the data, fitted background and expected signal event counts for a benchmark signal point in each bin of the mt distribution shown in figure 5 (middle right). The fit results are shown for an integrated luminosity of 3.2 fb-1.

Table shows the data, fitted background and expected signal event counts for a benchmark signal point in each bin of the met distribution shown in figure 5 (bottom left). The fit results are shown for an integrated luminosity of 3.2 fb-1.

Table shows the data, fitted background and expected signal event counts for a benchmark signal point in each bin of the met distribution shown in figure 5 (bottom right). The fit results are shown for an integrated luminosity of 3.2 fb-1.

Cutflow table for the hard-lepton signal regions with representative target signal models. The weighted numbers are normalized to 3.2 fb-1 and rounded to the statistical error.

Cutflow table for the hard-lepton signal regions with representative target signal models. The weighted numbers are normalized to 3.2 fb-1 and rounded to the statistical error.

Kinematic distributions for "extended SR selections" of the 2-tau channel, for mTsum in the Compressed SR selection without the mTsum>1400 GeV requirement. The last bin includes overflow events. Cited uncertainties are statistical uncertainties only. Signal predictions are overlaid for several benchmark models, normalised to their predicted cross sections. For the simplified model, "MM" refers to a medium mass splitting, with m(gluino)=1145 GeV and m(neutralino)=265 GeV; "HM" denotes a high mass splitting scenario, with m(gluino)=1305 GeV and m(neutralino)=105 GeV. The GMSB benchmark model corresponds to Lambda = 90 TeV and tanbeta = 40.

Kinematic distributions for "extended SR selections" of the 2-tau channel, for mTtau1+mTtau2 in the High-Mass SR selection without the mTtau1+mTtau2>350GeV requirement. The last bin includes overflow events. Cited uncertainties are statistical uncertainties only. Signal predictions are overlaid for several benchmark models, normalised to their predicted cross sections. For the simplified model, "MM" refers to a medium mass splitting, with m(gluino)=1145 GeV and m(neutralino)=265 GeV; "HM" denotes a high mass splitting scenario, with m(gluino)=1305 GeV and m(neutralino)=105 GeV. The GMSB benchmark model corresponds to Lambda = 90 TeV and tanbeta = 40.

Kinematic distributions for "extended SR selections" of the 2-tau channel, for HT in the GMSB SR selection without the HT > 1700 GeV requirement. The last bin includes overflow events. Cited uncertainties are statistical uncertainties only. Signal predictions are overlaid for several benchmark models, normalised to their predicted cross sections. For the simplified model, "MM" refers to a medium mass splitting, with m(gluino)=1145 GeV and m(neutralino)=265 GeV; "HM" denotes a high mass splitting scenario, with m(gluino)=1305 GeV and m(neutralino)=105 GeV. The GMSB benchmark model corresponds to Lambda = 90 TeV and tanbeta = 40.

Expected exclusion contour at the 95% confidence level for the simplified model of gluino pair production, based on the combined results from the 1tau and 2tau channel. The result is obtained using 3.2 fb-1 of sqrt(s) = 13 TeV ATLAS data.

Observed exclusion contour at the 95% confidence level for the simplified model of gluino pair production, based on the combined results from the 1tau and 2tau channel. The result is obtained using 3.2 fb-1 of sqrt(s) = 13 TeV ATLAS data.

Expected exclusion contour at the 95% confidence level for the simplified model of gluino pair production, based on results from the 2tau channel. The result is obtained using 3.2 fb-1 of sqrt(s) = 13 TeV ATLAS data.

Expected exclusion contour at the 95% confidence level for the simplified model of gluino pair production, based on results from the 1tau channel. The result is obtained using 3.2 fb-1 of sqrt(s) = 13 TeV ATLAS data.

Observed exclusion contours at the 95% confidence level for the gauge-mediated supersymmetry-breaking model, based on results from the 2 tau channel. The result is obtained using 3.2 fb-1 of sqrt(s) = 13 TeV ATLAS data. Additional model parameters are M(mess) = 250 TeV, N5 = 3, mu>0 and Cgrav =1.

Expected exclusion contours at the 95% confidence level for the gauge-mediated supersymmetry-breaking model, based on results from the 2 tau channel. The result is obtained using 3.2 fb-1 of sqrt(s) = 13 TeV ATLAS data. Additional model parameters are M(mess) = 250 TeV, N5 = 3, mu>0 and Cgrav =1.

Observed upper cross section limits in pb for the simplified model of gluino pair production for the combination of all SRs.

Best expected signal region for the simplified model of gluino pair production. The respective SR has been used in the combination of the results.

Acceptance for the gluino production simplified model grid in the Compressed 1tau signal region.

Efficiency for the gluino production simplified model grid in the Compressed 1tau signal region.

Acceptance times Efficiency for the gluino production simplified model grid in the Compressed 1tau signal region.

Acceptance for the gluino production simplified model grid in the medium mass 1tau signal region.

Efficiency for the gluino production simplified model grid in the medium mass 1tau signal region.

Acceptance times Efficiency for the gluino production simplified model grid in the medium mass 1tau signal region.

Acceptance for the gluino production simplified model grid in the high mass 1tau signal region.

Efficiency for the gluino production simplified model grid in the high mass 1tau signal region.

Acceptance times Efficiency for the gluino production simplified model grid in the high mass 1tau signal region.

Acceptance for the gluino production simplified model grid in the compressed 2tau signal region.

Efficiency for the gluino production simplified model grid in the compressed 2tau signal region.

Acceptance times Efficiency for the gluino production simplified model grid in the compressed 2tau signal region.

Acceptance for the gluino production simplified model grid in the high mass 2tau signal region.

Efficiency for the gluino production simplified model grid in the high mass 2tau signal region.

Acceptance times Efficiency for the gluino production simplified model grid in the high mass 2tau signal region.

Acceptance for the GMSB model grid in the 2tau signal region.

Efficiency for the GMSB model grid in the 2tau signal region.

Acceptance times Efficiency for the GMSB model grid in the 2tau signal region.

Fractions of selected events for several signal MC samples with a mass difference $\Delta m = 100$ GeV, illustrating how $\mathcal{A}\times\varepsilon$ varies with the model parameters.

Two-dimensional distribution of $m_{\mathrm{DV}}$ and track multiplicity for DVs in data events and events of a $R$-hadron signal sample with $m_{\tilde{g}} = 1.4$ TeV, $m_{\tilde{\chi}_{1}^{0}} = 100$ GeV and $\tau_{\tilde{g}} = 1$ ns that satisfy all signal region event selection criteria.

Two-dimensional distribution of $m_{\mathrm{DV}}$ and track multiplicity for DVs in data events and events of a $R$-hadron signal sample with $m_{\tilde{g}} = 1.4$ TeV, $m_{\tilde{\chi}_{1}^{0}} = 1.32$ TeV and $\tau_{\tilde{g}} = 1$ ns that satisfy all signal region event selection criteria.

Observed cross section upper 95% CL limits as a function of $m_{\tilde{g}}$ and $\tau$ for $m_{\tilde{\chi}_{1}^{0}}=100$ GeV. For the mass limits see the entry of Figure 8b.

Upper 95% CL limits on the signal cross section for $m_{\tilde{g}}=1.4$ TeV and fixed $m_{\tilde{\chi}_{1}^{0}}=100$ GeV as a function of lifetime $\tau$.

Vertex reconstruction efficiency as a function of radial position $R$ with and without the special LRT processing for one $R$-hadron signal sample with $m_{\tilde{g}} = 1.2$ TeV, $m_{\tilde{\chi}_{1}^{0}} = 100$ GeV and $\tau_{\tilde{g}} = 1$ ns. The efficiency is defined as the probability for a true LLP decay to be matched with a reconstructed DV fulfilling the vertex preselection criteria in events with a reconstructed primary vertex.

Upper 95% CL limits on the signal cross section for $m_{\tilde{g}}=2.0$ TeV and fixed $m_{\tilde{\chi}_{1}^{0}}=100$ GeV as a function of lifetime $\tau$.

Vertex reconstruction efficiency as a function of radial position $R$ for two $R$-hadron signal samples with $m_{\tilde{g}} = 1.2$ TeV, $\tau_{\tilde{g}} = 1$ ns and different neutralino masses. The efficiency is defined as the probability for a true LLP decay to be matched with a reconstructed DV fulfilling the vertex preselection criteria in events with a reconstructed primary vertex.

Lower 95% CL limits on $m_{\tilde{g}}$ for fixed $m_{\tilde{\chi}_{1}^{0}}=100$ GeV as a function of lifetime $\tau$.

Fractions of selected events for several signal MC samples with a gluino lifetime $\tau = 1$ ns, illustrating how $\mathcal{A}\times\varepsilon$ varies with the model parameters.

Upper 95% CL limits on the signal cross section for $m_{\tilde{g}}=1.4$ TeV and fixed $\Delta m=100$ GeV as a function of lifetime $\tau$.

Fractions of selected events for several signal MC samples with a mass difference $\Delta m = 100$ GeV, illustrating how $\mathcal{A}\times\varepsilon$ varies with the model parameters.

Upper 95% CL limits on the signal cross section for $m_{\tilde{g}}=2.0$ TeV and fixed $\Delta m=100$ GeV as a function of lifetime $\tau$.

Two-dimensional distribution of $m_{\mathrm{DV}}$ and track multiplicity for DVs in data events and events of a $R$-hadron signal sample with $m_{\tilde{g}} = 1.4$ TeV, $m_{\tilde{\chi}_{1}^{0}} = 100$ GeV and $\tau_{\tilde{g}} = 1$ ns that satisfy all signal region event selection criteria.

Lower 95% CL limit on $m_{\tilde{g}}$ for fixed $\Delta m=100$ GeV as a function of lifetime $\tau$.

Two-dimensional distribution of $m_{\mathrm{DV}}$ and track multiplicity for DVs in data events and events of a $R$-hadron signal sample with $m_{\tilde{g}} = 1.4$ TeV, $m_{\tilde{\chi}_{1}^{0}} = 1.32$ TeV and $\tau_{\tilde{g}} = 1$ ns that satisfy all signal region event selection criteria.

Upper 95% CL limits on the signal cross section for $m_{\tilde{g}}=1.4$ TeV and fixed $\tau=1$ ns as a function of $m_{\tilde{\chi}_{1}^{0}}$.

Upper 95% CL limits on the signal cross section for $m_{\tilde{g}}=1.4$ TeV and fixed $\Delta m=100$ GeV as a function of lifetime $\tau$.

Upper 95% CL limits on the signal cross section for $m_{\tilde{g}}=2.0$ TeV and fixed $\tau=1$ ns as a function of $m_{\tilde{\chi}_{1}^{0}}$.

Upper 95% CL limits on the signal cross section for $m_{\tilde{g}}=2.0$ TeV and fixed $\Delta m=100$ GeV as a function of lifetime $\tau$.

Observed 95% CL limit as a function of $m_{\tilde{g}}$ and $m_{\tilde{\chi}_{1}^{0}}$ for fixed $\tau=1$ ns.

Lower 95% CL limit on $m_{\tilde{g}}$ for fixed $\Delta m=100$ GeV as a function of lifetime $\tau$.

Two-dimensional distributions of $x$-$y$ positions of vertices observed in the data passing the vertex pre-selection and satisfying all signal region event-level requirements.

Two-dimensional distributions of $x$-$y$ positions of vertices observed in the data passing the vertex pre-selection and satisfying all signal region event-level requirements.

Distribution of the mass $m_{\mathrm{DV}}$ for vertices in data events and in events of five $R$-hadron signal samples with $m_{\tilde{g}} = 1.2$ TeV, $m_{\tilde{\chi}_{1}^{0}} = 100$ GeV and different $\tau_{\tilde{g}}$ that satisfy the signal region event requirements. All DV selections are applied except for the $m_{\mathrm{DV}}$ and track multiplicity requirements.

Distribution of the mass $m_{\mathrm{DV}}$ for vertices in data events and in events of five $R$-hadron signal samples with $m_{\tilde{g}} = 1.2$ TeV, $m_{\tilde{\chi}_{1}^{0}} = 100$ GeV and different $\tau_{\tilde{g}}$ that satisfy the signal region event requirements. All DV selections are applied except for the $m_{\mathrm{DV}}$ and track multiplicity requirements.

Distribution of the track multiplicity $n_{\mathrm{Tracks}}$ for vertices in data events and events of five $R$-hadron signal samples with $m_{\tilde{g}} = 1.2$ TeV, $m_{\tilde{\chi}_{1}^{0}} = 100$ GeV and and different $\tau_{\tilde{g}}$ that satisfy the signal region event requirements. All DV selections are applied except for the $m_{\mathrm{DV}}$ and track multiplicity requirements. The track multiplicity distribution requires vertices to have $m_{\mathrm{DV}}>3$ GeV.

Distribution of the track multiplicity $n_{\mathrm{Tracks}}$ for vertices in data events and events of five $R$-hadron signal samples with $m_{\tilde{g}} = 1.2$ TeV, $m_{\tilde{\chi}_{1}^{0}} = 100$ GeV and and different $\tau_{\tilde{g}}$ that satisfy the signal region event requirements. All DV selections are applied except for the $m_{\mathrm{DV}}$ and track multiplicity requirements. The track multiplicity distribution requires vertices to have $m_{\mathrm{DV}}>3$ GeV.

Observed cross section upper 95% CL limits as a function of $m_{\tilde{g}}$ and $\tau$ for $m_{\tilde{\chi}_{1}^{0}}=100$ GeV. For the mass limits see the entry of Figure 8b.

Observed cross section upper 95% CL limits as a function of $m_{\tilde{g}}$ and $\tau$ for $\Delta m=100$ GeV. For the mass limits see the entry of Figure 9b.

Observed cross section upper 95% CL limits as a function of $m_{\tilde{g}}$ and $\tau$ for $\Delta m=100$ GeV. For the mass limits see the entry of Figure 9b.

Observed cross section upper 95% CL limits as a function of $m_{\tilde{\chi}_{1}^{0}}$ and $m_{\tilde{g}}$ for $\tau = 1$ ns. For the mass limits see the entry of Figure 10b.

Observed cross section upper 95% CL limits as a function of $m_{\tilde{\chi}_{1}^{0}}$ and $m_{\tilde{g}}$ for $\tau = 1$ ns. For the mass limits see the entry of Figure 10b.

Parameterized event selection efficiencies as a function of truth MET for events which have all truth decay vertices occurring before the start of the ATLAS calorimeter. Event-level efficiencies are evaluated for events that have truth MET $> 200$ GeV, pass the trackless jet requirement, and have at least one displaced truth decay within the fiducial volume. To satisfy the event-level efficiency, events must then pass the full event selection.

Parameterized event selection efficiencies as a function of truth MET for events which have all truth decay vertices occurring before the start of the ATLAS calorimeter. Event-level efficiencies are evaluated for events that have truth MET $> 200$ GeV, pass the trackless jet requirement, and have at least one displaced truth decay within the fiducial volume. To satisfy the event-level efficiency, events must then pass the full event selection.

Parameterized event selection efficiencies as a function of truth MET for events which have the furthest truth decay occurring inside the calorimeter. Event-level efficiencies are evaluated for events that have truth MET $> 200$ GeV, pass the trackless jet requirement, and have at least one displaced truth decay within the fiducial volume. To satisfy the event-level efficiency, events must then pass the full event selection.

Parameterized event selection efficiencies as a function of truth MET for events which have the furthest truth decay occurring inside the calorimeter. Event-level efficiencies are evaluated for events that have truth MET $> 200$ GeV, pass the trackless jet requirement, and have at least one displaced truth decay within the fiducial volume. To satisfy the event-level efficiency, events must then pass the full event selection.

Parameterized event selection efficiencies as a function of truth MET for events which have the furthest truth decay occurring after the end of the ATLAS calorimeter. Event-level efficiencies are evaluated for events that have truth MET $> 200$ GeV, pass the trackless jet requirement, and have at least one displaced truth decay within the fiducial volume. To satisfy the event-level efficiency, events must then pass the full event selection.

Parameterized event selection efficiencies as a function of truth MET for events which have the furthest truth decay occurring after the end of the ATLAS calorimeter. Event-level efficiencies are evaluated for events that have truth MET $> 200$ GeV, pass the trackless jet requirement, and have at least one displaced truth decay within the fiducial volume. To satisfy the event-level efficiency, events must then pass the full event selection.

Parameterized vertex level efficiencies as a function of number of particles associated to a truth decay vertex, and the vertex invariant mass for truth decays with $4$ mm $< R_{\mathrm{decay}} < 22$ mm. Selected particles are required to have nonzero electric charge, $p_{T}(|Q|=1) > 1$ GeV, and $d_0 > 2$ mm. The per-vertex efficiency is evaluated only for truth vertices that have at least 5 associated tracks, an invariant mass $> 10$ GeV, and are in the region $4$ mm $< R_{\mathrm{decay}} < 300$ mm, and $|Z_{\mathrm{decay}}| < 300$ mm. A truth vertex satisfies the vertex level efficiency if it can be matched to a reconstructed DV which passes the final vertex selection.

Parameterized vertex level efficiencies as a function of number of particles associated to a truth decay vertex, and the vertex invariant mass for truth decays with $4$ mm $< R_{\mathrm{decay}} < 22$ mm. Selected particles are required to have nonzero electric charge, $p_{T}(|Q|=1) > 1$ GeV, and $d_0 > 2$ mm. The per-vertex efficiency is evaluated only for truth vertices that have at least 5 associated tracks, an invariant mass $> 10$ GeV, and are in the region $4$ mm $< R_{\mathrm{decay}} < 300$ mm, and $|Z_{\mathrm{decay}}| < 300$ mm. A truth vertex satisfies the vertex level efficiency if it can be matched to a reconstructed DV which passes the final vertex selection.

Parameterized vertex level efficiencies as a function of number of particles associated to a truth decay vertex, and the vertex invariant mass for truth decays with $22$ mm $< R_{\mathrm{decay}} < 25$ mm. Selected particles are required to have nonzero electric charge, $p_{T}(|Q|=1) > 1$ GeV, and $d_0 > 2$ mm. The per-vertex efficiency is evaluated only for truth vertices that have at least 5 associated tracks, an invariant mass $> 10$ GeV, and are in the region $4$ mm $< R_{\mathrm{decay}} < 300$ mm, and $|Z_{\mathrm{decay}}| < 300$ mm. A truth vertex satisfies the vertex level efficiency if it can be matched to a reconstructed DV which passes the final vertex selection.

Parameterized vertex level efficiencies as a function of number of particles associated to a truth decay vertex, and the vertex invariant mass for truth decays with $22$ mm $< R_{\mathrm{decay}} < 25$ mm. Selected particles are required to have nonzero electric charge, $p_{T}(|Q|=1) > 1$ GeV, and $d_0 > 2$ mm. The per-vertex efficiency is evaluated only for truth vertices that have at least 5 associated tracks, an invariant mass $> 10$ GeV, and are in the region $4$ mm $< R_{\mathrm{decay}} < 300$ mm, and $|Z_{\mathrm{decay}}| < 300$ mm. A truth vertex satisfies the vertex level efficiency if it can be matched to a reconstructed DV which passes the final vertex selection.

Parameterized vertex level efficiencies as a function of number of particles associated to a truth decay vertex, and the vertex invariant mass for truth decays with $25$ mm $< R_{\mathrm{decay}} < 29$ mm. Selected particles are required to have nonzero electric charge, $p_{T}(|Q|=1) > 1$ GeV, and $d_0 > 2$ mm. The per-vertex efficiency is evaluated only for truth vertices that have at least 5 associated tracks, an invariant mass $> 10$ GeV, and are in the region $4$ mm $< R_{\mathrm{decay}} < 300$ mm, and $|Z_{\mathrm{decay}}| < 300$ mm. A truth vertex satisfies the vertex level efficiency if it can be matched to a reconstructed DV which passes the final vertex selection.

Parameterized vertex level efficiencies as a function of number of particles associated to a truth decay vertex, and the vertex invariant mass for truth decays with $25$ mm $< R_{\mathrm{decay}} < 29$ mm. Selected particles are required to have nonzero electric charge, $p_{T}(|Q|=1) > 1$ GeV, and $d_0 > 2$ mm. The per-vertex efficiency is evaluated only for truth vertices that have at least 5 associated tracks, an invariant mass $> 10$ GeV, and are in the region $4$ mm $< R_{\mathrm{decay}} < 300$ mm, and $|Z_{\mathrm{decay}}| < 300$ mm. A truth vertex satisfies the vertex level efficiency if it can be matched to a reconstructed DV which passes the final vertex selection.

Upper 95% CL limits on the signal cross section for $m_{\tilde{g}}=1.4$ TeV and fixed $m_{\tilde{\chi}_{1}^{0}}=100$ GeV as a function of lifetime $\tau$.

Parameterized vertex level efficiencies as a function of number of particles associated to a truth decay vertex, and the vertex invariant mass for truth decays with $29$ mm $< R_{\mathrm{decay}} < 38$ mm. Selected particles are required to have nonzero electric charge, $p_{T}(|Q|=1) > 1$ GeV, and $d_0 > 2$ mm. The per-vertex efficiency is evaluated only for truth vertices that have at least 5 associated tracks, an invariant mass $> 10$ GeV, and are in the region $4$ mm $< R_{\mathrm{decay}} < 300$ mm, and $|Z_{\mathrm{decay}}| < 300$ mm. A truth vertex satisfies the vertex level efficiency if it can be matched to a reconstructed DV which passes the final vertex selection.

Parameterized vertex level efficiencies as a function of number of particles associated to a truth decay vertex, and the vertex invariant mass for truth decays with $29$ mm $< R_{\mathrm{decay}} < 38$ mm. Selected particles are required to have nonzero electric charge, $p_{T}(|Q|=1) > 1$ GeV, and $d_0 > 2$ mm. The per-vertex efficiency is evaluated only for truth vertices that have at least 5 associated tracks, an invariant mass $> 10$ GeV, and are in the region $4$ mm $< R_{\mathrm{decay}} < 300$ mm, and $|Z_{\mathrm{decay}}| < 300$ mm. A truth vertex satisfies the vertex level efficiency if it can be matched to a reconstructed DV which passes the final vertex selection.

Upper 95% CL limits on the signal cross section for $m_{\tilde{g}}=2.0$ TeV and fixed $m_{\tilde{\chi}_{1}^{0}}=100$ GeV as a function of lifetime $\tau$.

Parameterized vertex level efficiencies as a function of number of particles associated to a truth decay vertex, and the vertex invariant mass for truth decays with $38$ mm $< R_{\mathrm{decay}} < 46$ mm. Selected particles are required to have nonzero electric charge, $p_{T}(|Q|=1) > 1$ GeV, and $d_0 > 2$ mm. The per-vertex efficiency is evaluated only for truth vertices that have at least 5 associated tracks, an invariant mass $> 10$ GeV, and are in the region $4$ mm $< R_{\mathrm{decay}} < 300$ mm, and $|Z_{\mathrm{decay}}| < 300$ mm. A truth vertex satisfies the vertex level efficiency if it can be matched to a reconstructed DV which passes the final vertex selection.

Parameterized vertex level efficiencies as a function of number of particles associated to a truth decay vertex, and the vertex invariant mass for truth decays with $38$ mm $< R_{\mathrm{decay}} < 46$ mm. Selected particles are required to have nonzero electric charge, $p_{T}(|Q|=1) > 1$ GeV, and $d_0 > 2$ mm. The per-vertex efficiency is evaluated only for truth vertices that have at least 5 associated tracks, an invariant mass $> 10$ GeV, and are in the region $4$ mm $< R_{\mathrm{decay}} < 300$ mm, and $|Z_{\mathrm{decay}}| < 300$ mm. A truth vertex satisfies the vertex level efficiency if it can be matched to a reconstructed DV which passes the final vertex selection.

Lower 95% CL limits on $m_{\tilde{g}}$ for fixed $m_{\tilde{\chi}_{1}^{0}}=100$ GeV as a function of lifetime $\tau$.

Parameterized vertex level efficiencies as a function of number of particles associated to a truth decay vertex, and the vertex invariant mass for truth decays with $46$ mm $< R_{\mathrm{decay}} < 73$ mm. Selected particles are required to have nonzero electric charge, $p_{T}(|Q|=1) > 1$ GeV, and $d_0 > 2$ mm. The per-vertex efficiency is evaluated only for truth vertices that have at least 5 associated tracks, an invariant mass $> 10$ GeV, and are in the region $4$ mm $< R_{\mathrm{decay}} < 300$ mm, and $|Z_{\mathrm{decay}}| < 300$ mm. A truth vertex satisfies the vertex level efficiency if it can be matched to a reconstructed DV which passes the final vertex selection.

Parameterized vertex level efficiencies as a function of number of particles associated to a truth decay vertex, and the vertex invariant mass for truth decays with $46$ mm $< R_{\mathrm{decay}} < 73$ mm. Selected particles are required to have nonzero electric charge, $p_{T}(|Q|=1) > 1$ GeV, and $d_0 > 2$ mm. The per-vertex efficiency is evaluated only for truth vertices that have at least 5 associated tracks, an invariant mass $> 10$ GeV, and are in the region $4$ mm $< R_{\mathrm{decay}} < 300$ mm, and $|Z_{\mathrm{decay}}| < 300$ mm. A truth vertex satisfies the vertex level efficiency if it can be matched to a reconstructed DV which passes the final vertex selection.

Parameterized vertex level efficiencies as a function of number of particles associated to a truth decay vertex, and the vertex invariant mass for truth decays with $73$ mm $< R_{\mathrm{decay}} < 84$ mm. Selected particles are required to have nonzero electric charge, $p_{T}(|Q|=1) > 1$ GeV, and $d_0 > 2$ mm. The per-vertex efficiency is evaluated only for truth vertices that have at least 5 associated tracks, an invariant mass $> 10$ GeV, and are in the region $4$ mm $< R_{\mathrm{decay}} < 300$ mm, and $|Z_{\mathrm{decay}}| < 300$ mm. A truth vertex satisfies the vertex level efficiency if it can be matched to a reconstructed DV which passes the final vertex selection.

Parameterized vertex level efficiencies as a function of number of particles associated to a truth decay vertex, and the vertex invariant mass for truth decays with $73$ mm $< R_{\mathrm{decay}} < 84$ mm. Selected particles are required to have nonzero electric charge, $p_{T}(|Q|=1) > 1$ GeV, and $d_0 > 2$ mm. The per-vertex efficiency is evaluated only for truth vertices that have at least 5 associated tracks, an invariant mass $> 10$ GeV, and are in the region $4$ mm $< R_{\mathrm{decay}} < 300$ mm, and $|Z_{\mathrm{decay}}| < 300$ mm. A truth vertex satisfies the vertex level efficiency if it can be matched to a reconstructed DV which passes the final vertex selection.

Parameterized vertex level efficiencies as a function of number of particles associated to a truth decay vertex, and the vertex invariant mass for truth decays with $84$ mm $< R_{\mathrm{decay}} < 111$ mm. Selected particles are required to have nonzero electric charge, $p_{T}(|Q|=1) > 1$ GeV, and $d_0 > 2$ mm. The per-vertex efficiency is evaluated only for truth vertices that have at least 5 associated tracks, an invariant mass $> 10$ GeV, and are in the region $4$ mm $< R_{\mathrm{decay}} < 300$ mm, and $|Z_{\mathrm{decay}}| < 300$ mm. A truth vertex satisfies the vertex level efficiency if it can be matched to a reconstructed DV which passes the final vertex selection.

Parameterized vertex level efficiencies as a function of number of particles associated to a truth decay vertex, and the vertex invariant mass for truth decays with $84$ mm $< R_{\mathrm{decay}} < 111$ mm. Selected particles are required to have nonzero electric charge, $p_{T}(|Q|=1) > 1$ GeV, and $d_0 > 2$ mm. The per-vertex efficiency is evaluated only for truth vertices that have at least 5 associated tracks, an invariant mass $> 10$ GeV, and are in the region $4$ mm $< R_{\mathrm{decay}} < 300$ mm, and $|Z_{\mathrm{decay}}| < 300$ mm. A truth vertex satisfies the vertex level efficiency if it can be matched to a reconstructed DV which passes the final vertex selection.

Parameterized vertex level efficiencies as a function of number of particles associated to a truth decay vertex, and the vertex invariant mass for truth decays with $111$ mm $< R_{\mathrm{decay}} < 120$ mm. Selected particles are required to have nonzero electric charge, $p_{T}(|Q|=1) > 1$ GeV, and $d_0 > 2$ mm. The per-vertex efficiency is evaluated only for truth vertices that have at least 5 associated tracks, an invariant mass $> 10$ GeV, and are in the region $4$ mm $< R_{\mathrm{decay}} < 300$ mm, and $|Z_{\mathrm{decay}}| < 300$ mm. A truth vertex satisfies the vertex level efficiency if it can be matched to a reconstructed DV which passes the final vertex selection.

Parameterized vertex level efficiencies as a function of number of particles associated to a truth decay vertex, and the vertex invariant mass for truth decays with $111$ mm $< R_{\mathrm{decay}} < 120$ mm. Selected particles are required to have nonzero electric charge, $p_{T}(|Q|=1) > 1$ GeV, and $d_0 > 2$ mm. The per-vertex efficiency is evaluated only for truth vertices that have at least 5 associated tracks, an invariant mass $> 10$ GeV, and are in the region $4$ mm $< R_{\mathrm{decay}} < 300$ mm, and $|Z_{\mathrm{decay}}| < 300$ mm. A truth vertex satisfies the vertex level efficiency if it can be matched to a reconstructed DV which passes the final vertex selection.

Upper 95% CL limits on the signal cross section for $m_{\tilde{g}}=1.4$ TeV and fixed $\tau=1$ ns as a function of $m_{\tilde{\chi}_{1}^{0}}$.

Parameterized vertex level efficiencies as a function of number of particles associated to a truth decay vertex, and the vertex invariant mass for truth decays with $120$ mm $< R_{\mathrm{decay}} < 145$ mm. Selected particles are required to have nonzero electric charge, $p_{T}(|Q|=1) > 1$ GeV, and $d_0 > 2$ mm. The per-vertex efficiency is evaluated only for truth vertices that have at least 5 associated tracks, an invariant mass $> 10$ GeV, and are in the region $4$ mm $< R_{\mathrm{decay}} < 300$ mm, and $|Z_{\mathrm{decay}}| < 300$ mm. A truth vertex satisfies the vertex level efficiency if it can be matched to a reconstructed DV which passes the final vertex selection.

Parameterized vertex level efficiencies as a function of number of particles associated to a truth decay vertex, and the vertex invariant mass for truth decays with $120$ mm $< R_{\mathrm{decay}} < 145$ mm. Selected particles are required to have nonzero electric charge, $p_{T}(|Q|=1) > 1$ GeV, and $d_0 > 2$ mm. The per-vertex efficiency is evaluated only for truth vertices that have at least 5 associated tracks, an invariant mass $> 10$ GeV, and are in the region $4$ mm $< R_{\mathrm{decay}} < 300$ mm, and $|Z_{\mathrm{decay}}| < 300$ mm. A truth vertex satisfies the vertex level efficiency if it can be matched to a reconstructed DV which passes the final vertex selection.

Upper 95% CL limits on the signal cross section for $m_{\tilde{g}}=2.0$ TeV and fixed $\tau=1$ ns as a function of $m_{\tilde{\chi}_{1}^{0}}$.

Parameterized vertex level efficiencies as a function of number of particles associated to a truth decay vertex, and the vertex invariant mass for truth decays with $145$ mm $< R_{\mathrm{decay}} < 180$ mm. Selected particles are required to have nonzero electric charge, $p_{T}(|Q|=1) > 1$ GeV, and $d_0 > 2$ mm. The per-vertex efficiency is evaluated only for truth vertices that have at least 5 associated tracks, an invariant mass $> 10$ GeV, and are in the region $4$ mm $< R_{\mathrm{decay}} < 300$ mm, and $|Z_{\mathrm{decay}}| < 300$ mm. A truth vertex satisfies the vertex level efficiency if it can be matched to a reconstructed DV which passes the final vertex selection.

Parameterized vertex level efficiencies as a function of number of particles associated to a truth decay vertex, and the vertex invariant mass for truth decays with $145$ mm $< R_{\mathrm{decay}} < 180$ mm. Selected particles are required to have nonzero electric charge, $p_{T}(|Q|=1) > 1$ GeV, and $d_0 > 2$ mm. The per-vertex efficiency is evaluated only for truth vertices that have at least 5 associated tracks, an invariant mass $> 10$ GeV, and are in the region $4$ mm $< R_{\mathrm{decay}} < 300$ mm, and $|Z_{\mathrm{decay}}| < 300$ mm. A truth vertex satisfies the vertex level efficiency if it can be matched to a reconstructed DV which passes the final vertex selection.

Observed 95% CL limit as a function of $m_{\tilde{g}}$ and $m_{\tilde{\chi}_{1}^{0}}$ for fixed $\tau=1$ ns.

Parameterized vertex level efficiencies as a function of number of particles associated to a truth decay vertex, and the vertex invariant mass for truth decays with $180$ mm $< R_{\mathrm{decay}} < 300$ mm. Selected particles are required to have nonzero electric charge, $p_{T}(|Q|=1) > 1$ GeV, and $d_0 > 2$ mm. The per-vertex efficiency is evaluated only for truth vertices that have at least 5 associated tracks, an invariant mass $> 10$ GeV, and are in the region $4$ mm $< R_{\mathrm{decay}} < 300$ mm, and $|Z_{\mathrm{decay}}| < 300$ mm. A truth vertex satisfies the vertex level efficiency if it can be matched to a reconstructed DV which passes the final vertex selection.

Parameterized vertex level efficiencies as a function of number of particles associated to a truth decay vertex, and the vertex invariant mass for truth decays with $180$ mm $< R_{\mathrm{decay}} < 300$ mm. Selected particles are required to have nonzero electric charge, $p_{T}(|Q|=1) > 1$ GeV, and $d_0 > 2$ mm. The per-vertex efficiency is evaluated only for truth vertices that have at least 5 associated tracks, an invariant mass $> 10$ GeV, and are in the region $4$ mm $< R_{\mathrm{decay}} < 300$ mm, and $|Z_{\mathrm{decay}}| < 300$ mm. A truth vertex satisfies the vertex level efficiency if it can be matched to a reconstructed DV which passes the final vertex selection.

Lower mass requirement for signal regions.

Lower mass requirement for signal regions.

Lower mass requirement for signal regions.

Expected and observed events in the 16 discovery regions along with the according control regions.

Expected and observed events in the 16 discovery regions along with the according control regions.

Expected signal yield and acceptance x efficiency, estimated background and observed number of events in data for the full range of simulated masses in the MS-agnostic R-hadron search.

Expected signal yield and acceptance x efficiency, estimated background and observed number of events in data for the full range of simulated masses in the MS-agnostic R-hadron search.

Expected signal yield and acceptance x efficiency, estimated background and observed number of events in data for the full range of simulated masses in the full-detector R-hadron search.

Expected signal yield and acceptance x efficiency, estimated background and observed number of events in data for the full range of simulated masses in the full-detector R-hadron search.

p0-values and model-independent upper limits on cross-section x acceptance x efficiency for the 16 discovery regions.

p0-values and model-independent upper limits on cross-section x acceptance x efficiency for the 16 discovery regions.

Expected signal yield and acceptance x efficiency, estimated background and observed number of events in data for the full range of simulated masses in the MS-agnostic search for metastable gluino R-hadrons.

Expected signal yield and acceptance x efficiency, estimated background and observed number of events in data for the full range of simulated masses in the MS-agnostic search for metastable gluino R-hadrons.

Expected signal yield and acceptance x efficiency, estimated background and observed number of events in data for the full range of simulated masses in the full-detector direct-stau search.

Expected signal yield and acceptance x efficiency, estimated background and observed number of events in data for the full range of simulated masses in the full-detector direct-stau search.

Expected signal yield and acceptance x efficiency, estimated background and observed number of events in data for the full range of simulated masses in the full-detector chargino search.

Expected signal yield and acceptance x efficiency, estimated background and observed number of events in data for the full range of simulated masses in the full-detector chargino search.

Upper cross-section limit in gluino R-hadron search.

Upper cross-section limit in gluino R-hadron search.

Upper cross-section limit in sbottom R-hadron search.

Upper cross-section limit in sbottom R-hadron search.

Upper cross-section limit in stop R-hadron search.

Upper cross-section limit in stop R-hadron search.

Upper cross-section limit in stau search.

Upper cross-section limit in stau search.

Upper cross-section limit in chargino search.

Upper cross-section limit in chargino search.

Lower mass limit as function of gluino lifetime.

Lower mass limit as function of gluino lifetime.

Acceptance x efficiency, acceptance and efficiency for the full range of simulated masses in the MS-agnostic R-hadron search.

Acceptance x efficiency, acceptance and efficiency for the full range of simulated masses in the MS-agnostic R-hadron search.

Upper cross-section limit in meta-stable gluino R-hadron search.

Upper cross-section limit in meta-stable gluino R-hadron search.

Flavor composition of 800 GeV stop R-hadrons simulated using the generic model as a function of radial distance from the interaction point.

Flavor composition of 800 GeV stop R-hadrons simulated using the generic model as a function of radial distance from the interaction point.

Flavor composition of 800 GeV anti-stop R-hadrons simulated using the generic model as a function of radial distance from the interaction point.

Flavor composition of 800 GeV anti-stop R-hadrons simulated using the generic model as a function of radial distance from the interaction point.

Flavor composition of 800 GeV stop R-hadrons simulated using the Regge model as a function of radial distance from the interaction point.

Flavor composition of 800 GeV stop R-hadrons simulated using the Regge model as a function of radial distance from the interaction point.

Flavor composition of 800 GeV anti-stop R-hadrons simulated using the Regge model as a function of radial distance from the interaction point.

Flavor composition of 800 GeV anti-stop R-hadrons simulated using the Regge model as a function of radial distance from the interaction point.

ETmiss trigger efficiency as function of true ETmiss (EtmissTurnOn).

ETmiss trigger efficiency as function of true ETmiss (EtmissTurnOn).

Single-muon trigger efficiency as function of $|\eta|$ and $\beta$ (SingleMuTurnOn).

Single-muon trigger efficiency as function of $|\eta|$ and $\beta$ (SingleMuTurnOn).

Candidate reconstruction efficiency for ID+Calo selection (IDCaloEff).

Candidate reconstruction efficiency for ID+Calo selection (IDCaloEff).

Candidate reconstruction efficiency for loose selection (LooseEff).

Candidate reconstruction efficiency for loose selection (LooseEff).

Efficiency for a loose candidate to be promoted to a tight candidate (TightPromotionEff).

Efficiency for a loose candidate to be promoted to a tight candidate (TightPromotionEff).

Resolution and average of reconstructed dE/dx mass for a given simulated mass for ID+calo candidates.

Resolution and average of reconstructed dE/dx mass for a given simulated mass for ID+calo candidates.

Resolution and average of reconstructed ToF mass for a given simulated mass for ID+calo candidates.

Resolution and average of reconstructed ToF mass for a given simulated mass for ID+calo candidates.

Resolution and average of reconstructed ToF mass for a given simulated mass for FullDet candidates.

Resolution and average of reconstructed ToF mass for a given simulated mass for FullDet candidates.

Summary of the total uncertainty in the background prediction for each SR of the tt0L-low, tt0L-high, tt1L and tt2L analysis channels in the statistical combination. Their dominant contributions are indicated by individual lines. Individual uncertainties can be correlated, and do not necessarily add up in quadrature to the total background uncertainty.

Exclusion limits for colour-neutral scalar mediator dark matter models as a function of the mediator mass $m(\phi)$ for a DM mass $m_{\chi} = 1$ GeV. Associated production of DM with both single top quarks ($tW$ and $tj$ channels) and top quark pairs is considered. The limits are calculated at 95% CL and are expressed in terms of the ratio of the excluded cross section to the cross section for a coupling assumption of $g = g_q = g_{\chi} = 1$. The solid (dashed) lines show the observed (expected) exclusion limits for each individual channel and their statistical combination.

Exclusion limits for colour-neutral pseudoscalar mediator dark matter models as a function of the mediator mass $m(a)$ for a DM mass $m_{\chi} = 1$ GeV. Associated production of DM with both single top quarks ($tW$ and $tj$ channels) and top quark pairs is considered. The limits are calculated at 95% CL and are expressed in terms of the ratio of the excluded cross section to the cross section for a coupling assumption of $g = g_q = g_{\chi} = 1$. The solid (dashed) lines show the observed (expected) exclusion limits for each individual channel and their statistical combination.

$E_{\text{T}}^{\text{miss}}$ distribution in SR0X for the tt0L-low analysis. The contributions from all SM backgrounds are shown after the profile likelihood simultaneous fit to all tt0L-low CRs, with the hatched bands representing the total uncertainty. The category '$t\bar{t}$ (other)' represents $t\bar{t}$ events without extra jets or events with extra light-flavour jets. 'Other' includes contributions from $t\bar{t}W$, $tZ$ and $tWZ$ processes. The expected distributions for selected signal models are shown as dashed lines. The overflow events are included in the last bin. The bottom panels show the ratio of the observed data to the total SM background prediction, with the hatched area representing the total uncertainty in the background prediction and the red arrows marking data outside the vertical-axis range.

$E_{\text{T}}^{\text{miss}}$ distribution in SRWX for the tt0L-low analysis. The contributions from all SM backgrounds are shown after the profile likelihood simultaneous fit to all tt0L-low CRs, with the hatched bands representing the total uncertainty. The category '$t\bar{t}$ (other)' represents $t\bar{t}$ events without extra jets or events with extra light-flavour jets. 'Other' includes contributions from $t\bar{t}W$, $tZ$ and $tWZ$ processes. The expected distributions for selected signal models are shown as dashed lines. The overflow events are included in the last bin. The bottom panels show the ratio of the observed data to the total SM background prediction, with the hatched area representing the total uncertainty in the background prediction and the red arrows marking data outside the vertical-axis range.

$E_{\text{T}}^{\text{miss}}$ distribution in SRTX for the tt0L-low analysis. The contributions from all SM backgrounds are shown after the profile likelihood simultaneous fit to all tt0L-low CRs, with the hatched bands representing the total uncertainty. The category '$t\bar{t}$ (other)' represents $t\bar{t}$ events without extra jets or events with extra light-flavour jets. 'Other' includes contributions from $t\bar{t}W$, $tZ$ and $tWZ$ processes. The expected distributions for selected signal models are shown as dashed lines. The overflow events are included in the last bin. The bottom panels show the ratio of the observed data to the total SM background prediction, with the hatched area representing the total uncertainty in the background prediction and the red arrows marking data outside the vertical-axis range.

Exclusion limits for colour-neutral scalar mediator dark matter models as a function of the mediator mass $m(\phi)$ for a DM mass $m_{\chi} = 1$ GeV. Associated production of DM with both single top quarks ($tW$ and $tj$ channels) and top quark pairs is considered. The limits are calculated at 95% CL and are expressed in terms of the ratio of the excluded cross section to the nominal cross section for a coupling assumption of $g = g_q = g_{\chi} = 1$. The solid (dashed) lines show the observed (expected) exclusion limits for the tt0L-high and tt0L-low analyses and their statistical combination.

Exclusion limits for colour-neutral pseudoscalar mediator dark matter models as a function of the mediator mass $m(a)$ for a DM mass $m_{\chi} = 1$ GeV. Associated production of DM with both single top quarks ($tW$ and $tj$ channels) and top quark pairs is considered. The limits are calculated at 95% CL and are expressed in terms of the ratio of the excluded cross section to the nominal cross section for a coupling assumption of $g = g_q = g_{\chi} = 1$. The solid (dashed) lines show the observed (expected) exclusion limits for the tt0L-high and tt0L-low analyses and their statistical combination.

Exclusion limits for colour-neutral scalar mediator dark matter models as a function of the mediator mass $m(\phi)$ for a DM mass $m_{\chi} = 1$ GeV. Only associated production of DM with top quark pairs is considered for this interpretation. The limits are calculated at 95% CL and are expressed in terms of the ratio of the excluded cross section to the cross section for a coupling assumption of $g = g_q = g_{\chi} = 1$. The solid (dashed) lines show the observed (expected) exclusion limits for each individual channel and their statistical combination.

Exclusion limits for colour-neutral pseudoscalar mediator dark matter models as a function of the mediator mass $m(a)$ for a DM mass $m_{\chi} = 1$ GeV. Only associated production of DM with top quark pairs is considered for this interpretation. The limits are calculated at 95% CL and are expressed in terms of the ratio of the excluded cross section to the cross section for a coupling assumption of $g = g_q = g_{\chi} = 1$. The solid (dashed) lines show the observed (expected) exclusion limits for each individual channel and their statistical combination.

Exclusion limits for colour-neutral scalar mediator dark matter models as a function of the mediator mass $m(\phi)$ for a DM mass $m_{\chi} = 1$ GeV. Only associated production of DM with top quark pairs is considered for this interpretation. The limits are calculated at 95% CL and are expressed in terms of the ratio of the excluded cross section to the nominal cross section for a coupling assumption of $g = g_q = g_{\chi} = 1$. The solid (dashed) lines show the observed (expected) exclusion limits for the tt0L-high and tt0L-low analyses and their statistical combination.

Exclusion limits for colour-neutral pseudoscalar mediator dark matter models as a function of the mediator mass $m(a)$ for a DM mass $m_{\chi} = 1$ GeV. Only associated production of DM with top quark pairs is considered for this interpretation. The limits are calculated at 95% CL and are expressed in terms of the ratio of the excluded cross section to the nominal cross section for a coupling assumption of $g = g_q = g_{\chi} = 1$. The solid (dashed) lines show the observed (expected) exclusion limits for the tt0L-high and tt0L-low analyses and their statistical combination.

Representative fit distribution in the different flavour leptons signal region for the tt2L analysis: each bin of such distribution, starting from the red arrow, corresponds to a single SR included in the fit. 'FNP' includes the contribution from fake/non-prompt lepton background arising from jets (mainly $\pi/K$, heavy-flavour hadron decays and photon conversion) misidentified as leptons, estimated in a purely data-driven way. 'Other' includes contributions from $t\bar{t}W$, $tZ$ and $tWZ$ processes. The total uncertainty in the SM expectation is represented with hatched bands and the expected distributions for selected signal models are shown as dashed lines.

Signal acceptance in SR0X, SRWX and SRTX for simplified DM+$t\bar{t}$ model, defined as the number of accepted events at generator level in signal Monte Carlo simulation divided by the total number of events in the sample.

Signal acceptance in SR0X, SRWX and SRTX for simplified DM+$tW$ model, defined as the number of accepted events at generator level in signal Monte Carlo simulation divided by the total number of events in the sample.

Signal acceptance in SR0X, SRWX and SRTX for simplified DM+$tj$ model, defined as the number of accepted events at generator level in signal Monte Carlo simulation divided by the total number of events in the sample.

Signal efficiency in SR0X, SRWX and SRTX for simplified DM+$t\bar{t}$ model, defined as the number of selected reconstructed events divided by the acceptance.

Signal efficiency in SR0X, SRWX and SRTX for simplified DM+$tW$ model, defined as the number of selected reconstructed events divided by the acceptance.

Signal efficiency in SR0X, SRWX and SRTX for simplified DM+$tj$ model, defined as the number of selected reconstructed events divided by the acceptance.

Cutflow for the reference point DM+$t\bar{t}$ $m(\phi, \chi) = (10, 1)$ GeV in signal region SR0X. The column labelled 'weighted' shows the event yield including all correction factors applied to simulation, and is normalised to 139 fb$^{-1}$. A notable exception concerns the 'weighted' numbers in the first and the second row, labelled 'Total' and 'Filtered', which correspond to $\mathcal{L}\cdot\sigma$ and $\mathcal{L}\cdot\sigma\cdot\epsilon$ expected, respectively. The 'Skim' selection requires the $p_{\text{T}}$ of the leading four jets to be above (80, 60, 40, 40) GeV, the missing transverse momentum $E_{\text{T}}^{\text{miss}} > 140$ GeV, the missing momentum significance $\mathcal{S} > 8$, $\Delta\phi_{\min}(\vec{p}_{\text{T,1-4}},\vec{p}_{\text{T}}^{\text{miss}}) > 0.4$ and a lepton veto. The 'Orthogonalisation' selection is defined in the main body. In total 2045000 raw MC events were generated prior to the specified cuts, with the column 'Unweighted yield' collecting the numbers after each cut.

Cutflow for the reference point DM+$t\bar{t}$ $m(\phi, \chi) = (10, 1)$ GeV in signal region SRWX. The column labelled 'weighted' shows the event yield including all correction factors applied to simulation, and is normalised to 139 fb$^{-1}$. A notable exception concerns the 'weighted' numbers in the first and the second row, labelled 'Total' and 'Filtered', which correspond to $\mathcal{L}\cdot\sigma$ and $\mathcal{L}\cdot\sigma\cdot\epsilon$ expected, respectively. The 'Skim' selection requires the $p_{\text{T}}$ of the leading four jets to be above (80, 60, 40, 40) GeV, the missing transverse momentum $E_{\text{T}}^{\text{miss}} > 140$ GeV, the missing momentum significance $\mathcal{S} > 8$, $\Delta\phi_{\min}(\vec{p}_{\text{T,1-4}},\vec{p}_{\text{T}}^{\text{miss}}) > 0.4$ and a lepton veto. The 'Orthogonalisation' selection is defined in the main body. In total 2045000 raw MC events were generated prior to the specified cuts, with the column 'Unweighted yield' collecting the numbers after each cut.

Cutflow for the reference point DM+$t\bar{t}$ $m(\phi, \chi) = (10, 1)$ GeV in signal region SRTX. The column labelled 'weighted' shows the event yield including all correction factors applied to simulation, and is normalised to 139 fb$^{-1}$. A notable exception concerns the 'weighted' numbers in the first and the second row, labelled 'Total' and 'Filtered', which correspond to $\mathcal{L}\cdot\sigma$ and $\mathcal{L}\cdot\sigma\cdot\epsilon$ expected, respectively. The 'Skim' selection requires the $p_{\text{T}}$ of the leading four jets to be above (80, 60, 40, 40) GeV, the missing transverse momentum $E_{\text{T}}^{\text{miss}} > 140$ GeV, the missing momentum significance $\mathcal{S} > 8$, $\Delta\phi_{\min}(\vec{p}_{\text{T,1-4}},\vec{p}_{\text{T}}^{\text{miss}}) > 0.4$ and a lepton veto. The 'Orthogonalisation' selection is defined in the main body. In total 2045000 raw MC events were generated prior to the specified cuts, with the column 'Unweighted yield' collecting the numbers after each cut.

Cutflow for the reference point DM+$t\bar{t}$ $m(a, \chi) = (10, 1)$ GeV in signal region SR0X. The column labelled 'weighted' shows the event yield including all correction factors applied to simulation, and is normalised to 139 fb$^{-1}$. A notable exception concerns the 'weighted' numbers in the first and the second row, labelled 'Total' and 'Filtered', which correspond to $\mathcal{L}\cdot\sigma$ and $\mathcal{L}\cdot\sigma\cdot\epsilon$ expected, respectively. The 'Skim' selection requires the $p_{\text{T}}$ of the leading four jets to be above (80, 60, 40, 40) GeV, the missing transverse momentum $E_{\text{T}}^{\text{miss}} > 140$ GeV, the missing momentum significance $\mathcal{S} > 8$, $\Delta\phi_{\min}(\vec{p}_{\text{T,1-4}},\vec{p}_{\text{T}}^{\text{miss}}) > 0.4$ and a lepton veto. The 'Orthogonalisation' selection is defined in the main body. In total 400000 raw MC events were generated prior to the specified cuts, with the column 'Unweighted yield' collecting the numbers after each cut.

Cutflow for the reference point DM+$t\bar{t}$ $m(a, \chi) = (10, 1)$ GeV in signal region SRWX. The column labelled 'weighted' shows the event yield including all correction factors applied to simulation, and is normalised to 139 fb$^{-1}$. A notable exception concerns the 'weighted' numbers in the first and the second row, labelled 'Total' and 'Filtered', which correspond to $\mathcal{L}\cdot\sigma$ and $\mathcal{L}\cdot\sigma\cdot\epsilon$ expected, respectively. The 'Skim' selection requires the $p_{\text{T}}$ of the leading four jets to be above (80, 60, 40, 40) GeV, the missing transverse momentum $E_{\text{T}}^{\text{miss}} > 140$ GeV, the missing momentum significance $\mathcal{S} > 8$, $\Delta\phi_{\min}(\vec{p}_{\text{T,1-4}},\vec{p}_{\text{T}}^{\text{miss}}) > 0.4$ and a lepton veto. The 'Orthogonalisation' selection is defined in the main body. In total 400000 raw MC events were generated prior to the specified cuts, with the column 'Unweighted yield' collecting the numbers after each cut.

Cutflow for the reference point DM+$t\bar{t}$ $m(a, \chi) = (10, 1)$ GeV in signal region SRTX. The column labelled 'weighted' shows the event yield including all correction factors applied to simulation, and is normalised to 139 fb$^{-1}$. A notable exception concerns the 'weighted' numbers in the first and the second row, labelled 'Total' and 'Filtered', which correspond to $\mathcal{L}\cdot\sigma$ and $\mathcal{L}\cdot\sigma\cdot\epsilon$ expected, respectively. The 'Skim' selection requires the $p_{\text{T}}$ of the leading four jets to be above (80, 60, 40, 40) GeV, the missing transverse momentum $E_{\text{T}}^{\text{miss}} > 140$ GeV, the missing momentum significance $\mathcal{S} > 8$, $\Delta\phi_{\min}(\vec{p}_{\text{T,1-4}},\vec{p}_{\text{T}}^{\text{miss}}) > 0.4$ and a lepton veto. The 'Orthogonalisation' selection is defined in the main body. In total 400000 raw MC events were generated prior to the specified cuts, with the column 'Unweighted yield' collecting the numbers after each cut.

Cutflow for the reference point DM+$tW$ $m(\phi, \chi) = (10, 1)$ GeV in signal region SR0X. The column labelled 'weighted' shows the event yield including all correction factors applied to simulation, and is normalised to 139 fb$^{-1}$. A notable exception concerns the 'weighted' numbers in the first and the second row, labelled 'Total' and 'Filtered', which correspond to $\mathcal{L}\cdot\sigma$ and $\mathcal{L}\cdot\sigma\cdot\epsilon$ expected, respectively. The 'Skim' selection requires the $p_{\text{T}}$ of the leading four jets to be above (80, 60, 40, 40) GeV, the missing transverse momentum $E_{\text{T}}^{\text{miss}} > 140$ GeV, the missing momentum significance $\mathcal{S} > 8$, $\Delta\phi_{\min}(\vec{p}_{\text{T,1-4}},\vec{p}_{\text{T}}^{\text{miss}}) > 0.4$ and a lepton veto. The 'Orthogonalisation' selection is defined in the main body. In total 120000 raw MC events were generated prior to the specified cuts, with the column 'Unweighted yield' collecting the numbers after each cut.

Cutflow for the reference point DM+$tW$ $m(\phi, \chi) = (10, 1)$ GeV in signal region SRWX. The column labelled 'weighted' shows the event yield including all correction factors applied to simulation, and is normalised to 139 fb$^{-1}$. A notable exception concerns the 'weighted' numbers in the first and the second row, labelled 'Total' and 'Filtered', which correspond to $\mathcal{L}\cdot\sigma$ and $\mathcal{L}\cdot\sigma\cdot\epsilon$ expected, respectively. The 'Skim' selection requires the $p_{\text{T}}$ of the leading four jets to be above (80, 60, 40, 40) GeV, the missing transverse momentum $E_{\text{T}}^{\text{miss}} > 140$ GeV, the missing momentum significance $\mathcal{S} > 8$, $\Delta\phi_{\min}(\vec{p}_{\text{T,1-4}},\vec{p}_{\text{T}}^{\text{miss}}) > 0.4$ and a lepton veto. The 'Orthogonalisation' selection is defined in the main body. In total 120000 raw MC events were generated prior to the specified cuts, with the column 'Unweighted yield' collecting the numbers after each cut.

Cutflow for the reference point DM+$tW$ $m(\phi, \chi) = (10, 1)$ GeV in signal region SRTX. The column labelled 'weighted' shows the event yield including all correction factors applied to simulation, and is normalised to 139 fb$^{-1}$. A notable exception concerns the 'weighted' numbers in the first and the second row, labelled 'Total' and 'Filtered', which correspond to $\mathcal{L}\cdot\sigma$ and $\mathcal{L}\cdot\sigma\cdot\epsilon$ expected, respectively. The 'Skim' selection requires the $p_{\text{T}}$ of the leading four jets to be above (80, 60, 40, 40) GeV, the missing transverse momentum $E_{\text{T}}^{\text{miss}} > 140$ GeV, the missing momentum significance $\mathcal{S} > 8$, $\Delta\phi_{\min}(\vec{p}_{\text{T,1-4}},\vec{p}_{\text{T}}^{\text{miss}}) > 0.4$ and a lepton veto. The 'Orthogonalisation' selection is defined in the main body. In total 120000 raw MC events were generated prior to the specified cuts, with the column 'Unweighted yield' collecting the numbers after each cut.

Cutflow for the reference point DM+$tW$ $m(a, \chi) = (10, 1)$ GeV in signal region SR0X. The column labelled 'weighted' shows the event yield including all correction factors applied to simulation, and is normalised to 139 fb$^{-1}$. A notable exception concerns the 'weighted' numbers in the first and the second row, labelled 'Total' and 'Filtered', which correspond to $\mathcal{L}\cdot\sigma$ and $\mathcal{L}\cdot\sigma\cdot\epsilon$ expected, respectively. The 'Skim' selection requires the $p_{\text{T}}$ of the leading four jets to be above (80, 60, 40, 40) GeV, the missing transverse momentum $E_{\text{T}}^{\text{miss}} > 140$ GeV, the missing momentum significance $\mathcal{S} > 8$, $\Delta\phi_{\min}(\vec{p}_{\text{T,1-4}},\vec{p}_{\text{T}}^{\text{miss}}) > 0.4$ and a lepton veto. The 'Orthogonalisation' selection is defined in the main body. In total 100000 raw MC events were generated prior to the specified cuts, with the column 'Unweighted yield' collecting the numbers after each cut.

Cutflow for the reference point DM+$tW$ $m(a, \chi) = (10, 1)$ GeV in signal region SRWX. The column labelled 'weighted' shows the event yield including all correction factors applied to simulation, and is normalised to 139 fb$^{-1}$. A notable exception concerns the 'weighted' numbers in the first and the second row, labelled 'Total' and 'Filtered', which correspond to $\mathcal{L}\cdot\sigma$ and $\mathcal{L}\cdot\sigma\cdot\epsilon$ expected, respectively. The 'Skim' selection requires the $p_{\text{T}}$ of the leading four jets to be above (80, 60, 40, 40) GeV, the missing transverse momentum $E_{\text{T}}^{\text{miss}} > 140$ GeV, the missing momentum significance $\mathcal{S} > 8$, $\Delta\phi_{\min}(\vec{p}_{\text{T,1-4}},\vec{p}_{\text{T}}^{\text{miss}}) > 0.4$ and a lepton veto. The 'Orthogonalisation' selection is defined in the main body. In total 100000 raw MC events were generated prior to the specified cuts, with the column 'Unweighted yield' collecting the numbers after each cut.

Cutflow for the reference point DM+$tW$ $m(a, \chi) = (10, 1)$ GeV in signal region SRTX. The column labelled 'weighted' shows the event yield including all correction factors applied to simulation, and is normalised to 139 fb$^{-1}$. A notable exception concerns the 'weighted' numbers in the first and the second row, labelled 'Total' and 'Filtered', which correspond to $\mathcal{L}\cdot\sigma$ and $\mathcal{L}\cdot\sigma\cdot\epsilon$ expected, respectively. The 'Skim' selection requires the $p_{\text{T}}$ of the leading four jets to be above (80, 60, 40, 40) GeV, the missing transverse momentum $E_{\text{T}}^{\text{miss}} > 140$ GeV, the missing momentum significance $\mathcal{S} > 8$, $\Delta\phi_{\min}(\vec{p}_{\text{T,1-4}},\vec{p}_{\text{T}}^{\text{miss}}) > 0.4$ and a lepton veto. The 'Orthogonalisation' selection is defined in the main body. In total 100000 raw MC events were generated prior to the specified cuts, with the column 'Unweighted yield' collecting the numbers after each cut.

Cutflow for the reference point DM+$tj$ $m(\phi, \chi) = (10, 1)$ GeV in signal region SR0X. The column labelled 'weighted' shows the event yield including all correction factors applied to simulation, and is normalised to 139 fb$^{-1}$. A notable exception concerns the 'weighted' numbers in the first and the second row, labelled 'Total' and 'Filtered', which correspond to $\mathcal{L}\cdot\sigma$ and $\mathcal{L}\cdot\sigma\cdot\epsilon$ expected, respectively. The 'Skim' selection requires the $p_{\text{T}}$ of the leading four jets to be above (80, 60, 40, 40) GeV, the missing transverse momentum $E_{\text{T}}^{\text{miss}} > 140$ GeV, the missing momentum significance $\mathcal{S} > 8$, $\Delta\phi_{\min}(\vec{p}_{\text{T,1-4}},\vec{p}_{\text{T}}^{\text{miss}}) > 0.4$ and a lepton veto. The 'Orthogonalisation' selection is defined in the main body. In total 169000 raw MC events were generated prior to the specified cuts, with the column 'Unweighted yield' collecting the numbers after each cut.

Cutflow for the reference point DM+$tj$ $m(\phi, \chi) = (10, 1)$ GeV in signal region SRWX. The column labelled 'weighted' shows the event yield including all correction factors applied to simulation, and is normalised to 139 fb$^{-1}$. A notable exception concerns the 'weighted' numbers in the first and the second row, labelled 'Total' and 'Filtered', which correspond to $\mathcal{L}\cdot\sigma$ and $\mathcal{L}\cdot\sigma\cdot\epsilon$ expected, respectively. The 'Skim' selection requires the $p_{\text{T}}$ of the leading four jets to be above (80, 60, 40, 40) GeV, the missing transverse momentum $E_{\text{T}}^{\text{miss}} > 140$ GeV, the missing momentum significance $\mathcal{S} > 8$, $\Delta\phi_{\min}(\vec{p}_{\text{T,1-4}},\vec{p}_{\text{T}}^{\text{miss}}) > 0.4$ and a lepton veto. The 'Orthogonalisation' selection is defined in the main body. In total 169000 raw MC events were generated prior to the specified cuts, with the column 'Unweighted yield' collecting the numbers after each cut.

Cutflow for the reference point DM+$tj$ $m(\phi, \chi) = (10, 1)$ GeV in signal region SRTX. The column labelled 'weighted' shows the event yield including all correction factors applied to simulation, and is normalised to 139 fb$^{-1}$. A notable exception concerns the 'weighted' numbers in the first and the second row, labelled 'Total' and 'Filtered', which correspond to $\mathcal{L}\cdot\sigma$ and $\mathcal{L}\cdot\sigma\cdot\epsilon$ expected, respectively. The 'Skim' selection requires the $p_{\text{T}}$ of the leading four jets to be above (80, 60, 40, 40) GeV, the missing transverse momentum $E_{\text{T}}^{\text{miss}} > 140$ GeV, the missing momentum significance $\mathcal{S} > 8$, $\Delta\phi_{\min}(\vec{p}_{\text{T,1-4}},\vec{p}_{\text{T}}^{\text{miss}}) > 0.4$ and a lepton veto. The 'Orthogonalisation' selection is defined in the main body. In total 169000 raw MC events were generated prior to the specified cuts, with the column 'Unweighted yield' collecting the numbers after each cut.

Cutflow for the reference point DM+$tj$ $m(a, \chi) = (10, 1)$ GeV in signal region SR0X. The column labelled 'weighted' shows the event yield including all correction factors applied to simulation, and is normalised to 139 fb$^{-1}$. A notable exception concerns the 'weighted' numbers in the first and the second row, labelled 'Total' and 'Filtered', which correspond to $\mathcal{L}\cdot\sigma$ and $\mathcal{L}\cdot\sigma\cdot\epsilon$ expected, respectively. The 'Skim' selection requires the $p_{\text{T}}$ of the leading four jets to be above (80, 60, 40, 40) GeV, the missing transverse momentum $E_{\text{T}}^{\text{miss}} > 140$ GeV, the missing momentum significance $\mathcal{S} > 8$, $\Delta\phi_{\min}(\vec{p}_{\text{T,1-4}},\vec{p}_{\text{T}}^{\text{miss}}) > 0.4$ and a lepton veto. The 'Orthogonalisation' selection is defined in the main body. In total 140000 raw MC events were generated prior to the specified cuts, with the column 'Unweighted yield' collecting the numbers after each cut.

Cutflow for the reference point DM+$tj$ $m(a, \chi) = (10, 1)$ GeV in signal region SRWX. The column labelled 'weighted' shows the event yield including all correction factors applied to simulation, and is normalised to 139 fb$^{-1}$. A notable exception concerns the 'weighted' numbers in the first and the second row, labelled 'Total' and 'Filtered', which correspond to $\mathcal{L}\cdot\sigma$ and $\mathcal{L}\cdot\sigma\cdot\epsilon$ expected, respectively. The 'Skim' selection requires the $p_{\text{T}}$ of the leading four jets to be above (80, 60, 40, 40) GeV, the missing transverse momentum $E_{\text{T}}^{\text{miss}} > 140$ GeV, the missing momentum significance $\mathcal{S} > 8$, $\Delta\phi_{\min}(\vec{p}_{\text{T,1-4}},\vec{p}_{\text{T}}^{\text{miss}}) > 0.4$ and a lepton veto. The 'Orthogonalisation' selection is defined in the main body. In total 140000 raw MC events were generated prior to the specified cuts, with the column 'Unweighted yield' collecting the numbers after each cut.

Cutflow for the reference point DM+$tj$ $m(a, \chi) = (10, 1)$ GeV in signal region SRTX. The column labelled 'weighted' shows the event yield including all correction factors applied to simulation, and is normalised to 139 fb$^{-1}$. A notable exception concerns the 'weighted' numbers in the first and the second row, labelled 'Total' and 'Filtered', which correspond to $\mathcal{L}\cdot\sigma$ and $\mathcal{L}\cdot\sigma\cdot\epsilon$ expected, respectively. The 'Skim' selection requires the $p_{\text{T}}$ of the leading four jets to be above (80, 60, 40, 40) GeV, the missing transverse momentum $E_{\text{T}}^{\text{miss}} > 140$ GeV, the missing momentum significance $\mathcal{S} > 8$, $\Delta\phi_{\min}(\vec{p}_{\text{T,1-4}},\vec{p}_{\text{T}}^{\text{miss}}) > 0.4$ and a lepton veto. The 'Orthogonalisation' selection is defined in the main body. In total 140000 raw MC events were generated prior to the specified cuts, with the column 'Unweighted yield' collecting the numbers after each cut.