$E_{\mathrm{T}}^{\mathrm{miss}} / \sqrt{H_{\mathrm{T}}}$ distribution in signal region 10j50 2b. Two benchmark signal models are overlaid on the plot for comparison. Labelled `pMSSM' and `2-step', they show signal distributions from the example SUSY models (as described in the paper): a pMSSM slice model with ($m \tilde{g}$, $m \tilde{\chi_{1}^{\pm}}$) = (1300, 200) GeV and a cascade decay model with ($m \tilde{g}$, $m \tilde{\chi_{1}^{0}}$) = (1300, 200) GeV.

$E_{\mathrm{T}}^{\mathrm{miss}} / \sqrt{H_{\mathrm{T}}}$ distribution in signal region 8j80 0b. Two benchmark signal models are overlaid on the plot for comparison. Labelled `pMSSM' and `2-step', they show signal distributions from the example SUSY models (as described in the paper): a pMSSM slice model with ($m \tilde{g}$, $m \tilde{\chi_{1}^{\pm}}$) = (1300, 200) GeV and a cascade decay model with ($m \tilde{g}$, $m \tilde{\chi_{1}^{0}}$) = (1300, 200) GeV.

$E_{\mathrm{T}}^{\mathrm{miss}} / \sqrt{H_{\mathrm{T}}}$ distribution in signal region 8j80 2b. Two benchmark signal models are overlaid on the plot for comparison. Labelled `pMSSM' and `2-step', they show signal distributions from the example SUSY models (as described in the paper): a pMSSM slice model with ($m \tilde{g}$, $m \tilde{\chi_{1}^{\pm}}$) = (1300, 200) GeV and a cascade decay model with ($m \tilde{g}$, $m \tilde{\chi_{1}^{0}}$) = (1300, 200) GeV.

Observed 95% CL limit for the pMSSM grid.

Observed 95% CL limit for the pMSSM grid when the signal cross section is increased by one standard deviation.

Observed 95% CL limit for the pMSSM grid when the signal cross section is decreased by one standard deviation.

Expected 95% CL limit for the pMSSM grid.

+1 sigma excursion of the expected 95% CL limit for the pMSSM grid.

-1 sigma excursion of the expected 95% CL limit for the pMSSM grid.

Observed 95% CL limit for the 2Step grid.

Observed 95% CL limit for the 2Step grid when the signal cross section is increased by one standard deviation.

Observed 95% CL limit for the 2Step grid when the signal cross section is decreased by one standard deviation.

Expected 95% CL limit for the 2Step grid.

+1 sigma excursion of the expected 95% CL limit for the 2Step grid.

-1 sigma excursion of the expected 95% CL limit for the 2Step grid.

$E_{\mathrm{T}}^{\mathrm{miss}} / \sqrt{H_{\mathrm{T}}}$ distribution in signal region 8j50 0b. Two benchmark signal models are overlaid on the plot for comparison. Labelled `pMSSM' and `2-step', they show signal distributions from the example SUSY models (as described in the paper): a pMSSM slice model with ($m \tilde{g}$, $m \tilde{\chi_{1}^{\pm}}$) = (1300, 200) GeV and a cascade decay model with ($m \tilde{g}$, $m \tilde{\chi_{1}^{0}}$) = (1300, 200) GeV.

$E_{\mathrm{T}}^{\mathrm{miss}} / \sqrt{H_{\mathrm{T}}}$ distribution in signal region 8j50 1b. Two benchmark signal models are overlaid on the plot for comparison. Labelled `pMSSM' and `2-step', they show signal distributions from the example SUSY models (as described in the paper): a pMSSM slice model with ($m \tilde{g}$, $m \tilde{\chi_{1}^{\pm}}$) = (1300, 200) GeV and a cascade decay model with ($m \tilde{g}$, $m \tilde{\chi_{1}^{0}}$) = (1300, 200) GeV.

$E_{\mathrm{T}}^{\mathrm{miss}} / \sqrt{H_{\mathrm{T}}}$ distribution in signal region 8j50 2b. Two benchmark signal models are overlaid on the plot for comparison. Labelled `pMSSM' and `2-step', they show signal distributions from the example SUSY models (as described in the paper): a pMSSM slice model with ($m \tilde{g}$, $m \tilde{\chi_{1}^{\pm}}$) = (1300, 200) GeV and a cascade decay model with ($m \tilde{g}$, $m \tilde{\chi_{1}^{0}}$) = (1300, 200) GeV.

$E_{\mathrm{T}}^{\mathrm{miss}} / \sqrt{H_{\mathrm{T}}}$ distribution in signal region 9j50 0b. Two benchmark signal models are overlaid on the plot for comparison. Labelled `pMSSM' and `2-step', they show signal distributions from the example SUSY models (as described in the paper): a pMSSM slice model with ($m \tilde{g}$, $m \tilde{\chi_{1}^{\pm}}$) = (1300, 200) GeV and a cascade decay model with ($m \tilde{g}$, $m \tilde{\chi_{1}^{0}}$) = (1300, 200) GeV.

$E_{\mathrm{T}}^{\mathrm{miss}} / \sqrt{H_{\mathrm{T}}}$ distribution in signal region 9j50 1b. Two benchmark signal models are overlaid on the plot for comparison. Labelled `pMSSM' and `2-step', they show signal distributions from the example SUSY models (as described in the paper): a pMSSM slice model with ($m \tilde{g}$, $m \tilde{\chi_{1}^{\pm}}$) = (1300, 200) GeV and a cascade decay model with ($m \tilde{g}$, $m \tilde{\chi_{1}^{0}}$) = (1300, 200) GeV.

$E_{\mathrm{T}}^{\mathrm{miss}} / \sqrt{H_{\mathrm{T}}}$ distribution in signal region 9j50 2b. Two benchmark signal models are overlaid on the plot for comparison. Labelled `pMSSM' and `2-step', they show signal distributions from the example SUSY models (as described in the paper): a pMSSM slice model with ($m \tilde{g}$, $m \tilde{\chi_{1}^{\pm}}$) = (1300, 200) GeV and a cascade decay model with ($m \tilde{g}$, $m \tilde{\chi_{1}^{0}}$) = (1300, 200) GeV.

$E_{\mathrm{T}}^{\mathrm{miss}} / \sqrt{H_{\mathrm{T}}}$ distribution in signal region 10j50 0b. Two benchmark signal models are overlaid on the plot for comparison. Labelled `pMSSM' and `2-step', they show signal distributions from the example SUSY models (as described in the paper): a pMSSM slice model with ($m \tilde{g}$, $m \tilde{\chi_{1}^{\pm}}$) = (1300, 200) GeV and a cascade decay model with ($m \tilde{g}$, $m \tilde{\chi_{1}^{0}}$) = (1300, 200) GeV.

$E_{\mathrm{T}}^{\mathrm{miss}} / \sqrt{H_{\mathrm{T}}}$ distribution in signal region 10j50 1b. Two benchmark signal models are overlaid on the plot for comparison. Labelled `pMSSM' and `2-step', they show signal distributions from the example SUSY models (as described in the paper): a pMSSM slice model with ($m \tilde{g}$, $m \tilde{\chi_{1}^{\pm}}$) = (1300, 200) GeV and a cascade decay model with ($m \tilde{g}$, $m \tilde{\chi_{1}^{0}}$) = (1300, 200) GeV.

$E_{\mathrm{T}}^{\mathrm{miss}} / \sqrt{H_{\mathrm{T}}}$ distribution in signal region 10j50 2b. Two benchmark signal models are overlaid on the plot for comparison. Labelled `pMSSM' and `2-step', they show signal distributions from the example SUSY models (as described in the paper): a pMSSM slice model with ($m \tilde{g}$, $m \tilde{\chi_{1}^{\pm}}$) = (1300, 200) GeV and a cascade decay model with ($m \tilde{g}$, $m \tilde{\chi_{1}^{0}}$) = (1300, 200) GeV.

$E_{\mathrm{T}}^{\mathrm{miss}} / \sqrt{H_{\mathrm{T}}}$ distribution in signal region 7j80 0b. Two benchmark signal models are overlaid on the plot for comparison. Labelled `pMSSM' and `2-step', they show signal distributions from the example SUSY models (as described in the paper): a pMSSM slice model with ($m \tilde{g}$, $m \tilde{\chi_{1}^{\pm}}$) = (1300, 200) GeV and a cascade decay model with ($m \tilde{g}$, $m \tilde{\chi_{1}^{0}}$) = (1300, 200) GeV.

$E_{\mathrm{T}}^{\mathrm{miss}} / \sqrt{H_{\mathrm{T}}}$ distribution in signal region 7j80 1b. Two benchmark signal models are overlaid on the plot for comparison. Labelled `pMSSM' and `2-step', they show signal distributions from the example SUSY models (as described in the paper): a pMSSM slice model with ($m \tilde{g}$, $m \tilde{\chi_{1}^{\pm}}$) = (1300, 200) GeV and a cascade decay model with ($m \tilde{g}$, $m \tilde{\chi_{1}^{0}}$) = (1300, 200) GeV.

$E_{\mathrm{T}}^{\mathrm{miss}} / \sqrt{H_{\mathrm{T}}}$ distribution in signal region 7j80 2b. Two benchmark signal models are overlaid on the plot for comparison. Labelled `pMSSM' and `2-step', they show signal distributions from the example SUSY models (as described in the paper): a pMSSM slice model with ($m \tilde{g}$, $m \tilde{\chi_{1}^{\pm}}$) = (1300, 200) GeV and a cascade decay model with ($m \tilde{g}$, $m \tilde{\chi_{1}^{0}}$) = (1300, 200) GeV.

$E_{\mathrm{T}}^{\mathrm{miss}} / \sqrt{H_{\mathrm{T}}}$ distribution in signal region 8j80 0b. Two benchmark signal models are overlaid on the plot for comparison. Labelled `pMSSM' and `2-step', they show signal distributions from the example SUSY models (as described in the paper): a pMSSM slice model with ($m \tilde{g}$, $m \tilde{\chi_{1}^{\pm}}$) = (1300, 200) GeV and a cascade decay model with ($m \tilde{g}$, $m \tilde{\chi_{1}^{0}}$) = (1300, 200) GeV.

$E_{\mathrm{T}}^{\mathrm{miss}} / \sqrt{H_{\mathrm{T}}}$ distribution in signal region 8j80 1b. Two benchmark signal models are overlaid on the plot for comparison. Labelled `pMSSM' and `2-step', they show signal distributions from the example SUSY models (as described in the paper): a pMSSM slice model with ($m \tilde{g}$, $m \tilde{\chi_{1}^{\pm}}$) = (1300, 200) GeV and a cascade decay model with ($m \tilde{g}$, $m \tilde{\chi_{1}^{0}}$) = (1300, 200) GeV.

$E_{\mathrm{T}}^{\mathrm{miss}} / \sqrt{H_{\mathrm{T}}}$ distribution in signal region 8j80 2b. Two benchmark signal models are overlaid on the plot for comparison. Labelled `pMSSM' and `2-step', they show signal distributions from the example SUSY models (as described in the paper): a pMSSM slice model with ($m \tilde{g}$, $m \tilde{\chi_{1}^{\pm}}$) = (1300, 200) GeV and a cascade decay model with ($m \tilde{g}$, $m \tilde{\chi_{1}^{0}}$) = (1300, 200) GeV.

Degree of multijet closure for signal and vaidation regions with at no b-jet requirement. The solid lines are the pre-fit predicted numbers of events and the points are the observed numbers. The blue hatched band shows only the statistical (MC and data) uncertainty on the background estimate. The bins labelled in bold are signal regions, while the others are validation regions. The template closure uncertainty for each SR bin is given by the maximal deviation of data from prediction in any non-SR bin to its left on this plot (although those for 80 GeV regions are independent of deviations in 50 GeV regions).

Degree of multijet closure for signal and vaidation regions with at least 1 b-jet. The solid lines are the pre-fit predicted numbers of events and the points are the observed numbers. The blue hatched band shows only the statistical (MC and data) uncertainty on the background estimate. The bins labelled in bold are signal regions, while the others are validation regions. The template closure uncertainty for each SR bin is given by the maximal deviation of data from prediction in any non-SR bin to its left on this plot (although those for 80 GeV regions are independent of deviations in 50 GeV regions).

Degree of multijet closure for signal and vaidation regions with at least 2 b-jets. The solid lines are the pre-fit predicted numbers of events and the points are the observed numbers. The blue hatched band shows only the statistical (MC and data) uncertainty on the background estimate. The bins labelled in bold are signal regions, while the others are validation regions. The template closure uncertainty for each SR bin is given by the maximal deviation of data from prediction in any non-SR bin to its left on this plot (although those for 80 GeV regions are independent of deviations in 50 GeV regions).

Summary of all 15 signal regions (post-fit).

Signal region yielding the best-expected CLs value, the best expected CLs value, and the corresponding observed CLs value for the 2Step grid.

Signal region yielding the best-expected CLs value, the best expected CLs value, and the corresponding observed CLs value for the pMSSM grid.

95% CLs observed upper limit on model cross-section for 2-step signal points for the best-expected signal region.

Performance of the 8j50-0b selection for the pMSSM grid: number of generated signal events; total signal cross-section; acceptance; efficiency (fractional); observed CL using this region alone; expected CL using this region alone.

Performance of the 8j50-1b selection for the pMSSM grid: number of generated signal events; total signal cross-section; acceptance; efficiency (fractional); observed CL using this region alone; expected CL using this region alone.

Performance of the 8j50-2b selection for the pMSSM grid: number of generated signal events; total signal cross-section; acceptance; efficiency (fractional); observed CL using this region alone; expected CL using this region alone.

Performance of the 9j50-0b selection for the pMSSM grid: number of generated signal events; total signal cross-section; acceptance; efficiency (fractional); observed CL using this region alone; expected CL using this region alone.

Performance of the 9j50-1b selection for the pMSSM grid: number of generated signal events; total signal cross-section; acceptance; efficiency (fractional); observed CL using this region alone; expected CL using this region alone.

Performance of the 9j50-2b selection for the pMSSM grid: number of generated signal events; total signal cross-section; acceptance; efficiency (fractional); observed CL using this region alone; expected CL using this region alone.

Performance of the 10j50-0b selection for the pMSSM grid: number of generated signal events; total signal cross-section; acceptance; efficiency (fractional); observed CL using this region alone; expected CL using this region alone.

Performance of the 10j50-1b selection for the pMSSM grid: number of generated signal events; total signal cross-section; acceptance; efficiency (fractional); observed CL using this region alone; expected CL using this region alone.

Performance of the 10j50-2b selection for the pMSSM grid: number of generated signal events; total signal cross-section; acceptance; efficiency (fractional); observed CL using this region alone; expected CL using this region alone.

Performance of the 7j80-0b selection for the pMSSM grid: number of generated signal events; total signal cross-section; acceptance; efficiency (fractional); observed CL using this region alone; expected CL using this region alone.

Performance of the 7j80-1b selection for the pMSSM grid: number of generated signal events; total signal cross-section; acceptance; efficiency (fractional); observed CL using this region alone; expected CL using this region alone.

Performance of the 7j80-2b selection for the pMSSM grid: number of generated signal events; total signal cross-section; acceptance; efficiency (fractional); observed CL using this region alone; expected CL using this region alone.

Performance of the 8j80-0b selection for the pMSSM grid: number of generated signal events; total signal cross-section; acceptance; efficiency (fractional); observed CL using this region alone; expected CL using this region alone.

Performance of the 8j80-1b selection for the pMSSM grid: number of generated signal events; total signal cross-section; acceptance; efficiency (fractional); observed CL using this region alone; expected CL using this region alone.

Performance of the 8j80-2b selection for the pMSSM grid: number of generated signal events; total signal cross-section; acceptance; efficiency (fractional); observed CL using this region alone; expected CL using this region alone.

Performance of the 8j50-0b selection for the 2Step grid: number of generated signal events; total signal cross-section; acceptance; efficiency (fractional); observed CL using this region alone; expected CL using this region alone.

Performance of the 8j50-1b selection for the 2Step grid: number of generated signal events; total signal cross-section; acceptance; efficiency (fractional); observed CL using this region alone; expected CL using this region alone.

Performance of the 8j50-2b selection for the 2Step grid: number of generated signal events; total signal cross-section; acceptance; efficiency (fractional); observed CL using this region alone; expected CL using this region alone.

Performance of the 9j50-0b selection for the 2Step grid: number of generated signal events; total signal cross-section; acceptance; efficiency (fractional); observed CL using this region alone; expected CL using this region alone.

Performance of the 9j50-1b selection for the 2Step grid: number of generated signal events; total signal cross-section; acceptance; efficiency (fractional); observed CL using this region alone; expected CL using this region alone.

Performance of the 9j50-2b selection for the 2Step grid: number of generated signal events; total signal cross-section; acceptance; efficiency (fractional); observed CL using this region alone; expected CL using this region alone.

Performance of the 10j50-0b selection for the 2Step grid: number of generated signal events; total signal cross-section; acceptance; efficiency (fractional); observed CL using this region alone; expected CL using this region alone.

Performance of the 10j50-1b selection for the 2Step grid: number of generated signal events; total signal cross-section; acceptance; efficiency (fractional); observed CL using this region alone; expected CL using this region alone.

Performance of the 10j50-2b selection for the 2Step grid: number of generated signal events; total signal cross-section; acceptance; efficiency (fractional); observed CL using this region alone; expected CL using this region alone.

Performance of the 7j80-0b selection for the 2Step grid: number of generated signal events; total signal cross-section; acceptance; efficiency (fractional); observed CL using this region alone; expected CL using this region alone.

Performance of the 7j80-1b selection for the 2Step grid: number of generated signal events; total signal cross-section; acceptance; efficiency (fractional); observed CL using this region alone; expected CL using this region alone.

Performance of the 7j80-2b selection for the 2Step grid: number of generated signal events; total signal cross-section; acceptance; efficiency (fractional); observed CL using this region alone; expected CL using this region alone.

Performance of the 8j80-0b selection for the 2Step grid: number of generated signal events; total signal cross-section; acceptance; efficiency (fractional); observed CL using this region alone; expected CL using this region alone.

Performance of the 8j80-1b selection for the 2Step grid: number of generated signal events; total signal cross-section; acceptance; efficiency (fractional); observed CL using this region alone; expected CL using this region alone.

Performance of the 8j80-2b selection for the 2Step grid: number of generated signal events; total signal cross-section; acceptance; efficiency (fractional); observed CL using this region alone; expected CL using this region alone.

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]].

Detector response matrix (PROB) for the acoplanarity variable (ACO) for mu+ mu- channel (empty bins are not reported).

Acoplanarity (ACO) distributions unfolded for detector resolution, and lepton pair trigger, reconstruction and identification efficiencies for e+ e- channel (empty bins are not reported).

Acoplanarity (ACO) distributions unfolded for detector resolution, and lepton pair trigger, reconstruction and identification efficiencies for mu+ mu- channel (empty bins are not reported).

The lepton pair exclusive selection efficiency in the fiducial region for signal MC events as a function of M and Y for e+ e- channel.

The lepton pair exclusive selection efficiency in the fiducial region for signal MC events as a function of M and Y for mu+ mu- channel.

The lepton pair exclusive selection efficiency in the fiducial region for single-dissociative MC events as a function of M and Y for e+ e- channel.

The lepton pair exclusive selection efficiency in the fiducial region for single-dissociative MC events as a function of M and Y for mu+ mu- channel.

Acceptance times efficiency for the tight channel of the Razor analysis for the direct squark-squark model.

Acceptance for the 1tau channel of the Taus+jets+MET analysis for the gluino simplified model.

Efficiency for the 1tau channel of the Taus+jets+MET analysis for the gluino simplified model.

Acceptance for the 2tau channel of the Taus+jets+MET analysis for the gluino simplified model.

Efficiency for the 2tau channel of the Taus+jets+MET analysis for the gluino simplified model.

Acceptance for the 1tau channel of the Taus+jets+MET analysis for the squark simplified model.

Efficiency for the 1tau channel of the Taus+jets+MET analysis for the squark simplified model.

Acceptance for the 2tau channel of the Taus+jets+MET analysis for the squark simplified model.

Efficiency for the 2tau channel of the Taus+jets+MET analysis for the squark simplified model.

Observed limit on the cross-section times branching ratio in fb at 95% CL for the direct squark-squark model.

Observed limit on the cross-section times branching ratio in fb at 95% CL for the gluino simplified model, obtained by the Taus+jets+MET analysis.

Observed limit on the cross-section times branching ratio in fb at 95% CL for the squark simplified model, obtained by the Taus+jets+MET analysis.

Observed 95% CL cross section times branching ratio upper limit on pMSSM qLqL production with M1=60 GeV and $m_{\tilde g}=$1.6 TeV interpreted by the 0L analysis including its two dedicated signal regions.

Observed 95% CL cross section times branching ratio upper limit on pMSSM qLqL production with M2=(M1+mqL)/2 and $m_{\tilde g}=$1.6 TeV interpreted by the 0L analysis including its two dedicated signal regions.

Observed 95% CL cross section times branching ratio upper limit on pMSSM qLqL production with M1=60 GeV with $m_{\tilde g}=$2.2 TeV interpreted by the 0L analysis including its two dedicated signal regions.

Observed 95% CL cross section times branching ratio upper limit on pMSSM qLqL production with M2=(M1+mqL)/2 and $m_{\tilde g}=$2.2 TeV interpreted by the 0L analysis including its two dedicated signal regions.

Observed 95% CL cross section times branching ratio upper limit on pMSSM qLqL production with M1=60 GeV and $m_{\tilde g}=$3.0 TeV interpreted by the 0L analysis including its two dedicated signal regions.

Observed 95% CL cross section times branching ratio upper limit on pMSSM qLqL production with M2=(M1+mqL)/2 and $m_{\tilde g}=$3.0 TeV interpreted by the 0L analysis including its two dedicated signal regions.

Observed and expected limits on the cross-section times branching ratio in fb at 95% CL for the gluino-mediated stop model with off-shell stops, for mass points with four- or five-body gluino decays, obtained by the SS/3L and 0/1L3B analyses.

: Cut flow for signal regions 0L_4jt+ and 0L_5jt of baseline signal points normalized to 20.3 fb$^\mathrm{-1}$. The point used for the 0L_4jt+ cutflow is mqL = 1150 GeV, M1 = 60 GeV, M2 = 1000 GeV, m(gluino) = 2200 GeV. The point used for the 0L_5jt cutflow is mqL = 950 GeV, M1 = 100 GeV, M2 = 525 GeV, m(gluino) = 2200 GeV.

Cut flow for the 1L(H)_7-jet SR (electron channel) for the main background from the $t\bar{t}$ production. The events are normalized to 20.3 fb$^\mathrm{-1}$.

Cut flow for the 1L(H)_7-jet SR (muon channel) for the main background from the $t\bar{t}$ production. The events are normalized to 20.3 fb$^\mathrm{-1}$.

Cut flow for 1L(H)_7-jet SR (electron channel). The signal point is taken from the gluino simplified model with 1-step decay via charginos and parameters $m_{\tilde{g}}$ = 1145 GeV, $m_{{\tilde \chi}_{1}^{\pm}}$ = 785 GeV, $m_{\tilde{\chi}_{1}^{0}}$ = 425 GeV; 20000 events were generated for this point. The events are normalized to 20.3 fb$^\mathrm{-1}$.

Cut flow for 1L(H)_7-jet SR (muon channel). The signal point is taken from the gluino simplified model with 1-step decay via charginos and parameters $m_{\tilde{g}}$ = 1145 GeV, $m_{\tilde{\chi}_{1}^{\pm}}$ = 785 GeV, $m_{\tilde{\chi}_{1}^{0}}$ = 425 GeV; 20000 events were generated for this point. The events are normalized to 20.3 fb$^\mathrm{-1}$.

Cut flow for 1L(H)_7-jet SR (electron channel). The signal point is taken from the gluino simplified model with 1-step decay via charginos and parameters $m_{\tilde{g}}$ = 1025 GeV, $m_{\tilde{\chi}_{1}^{\pm}}$ = 545 GeV, $m_{\tilde{\chi}_{1}^{0}}$ = 65 GeV; 20000 events were generated for this point. The events are normalized to 20.3 fb$^\mathrm{-1}$.

Cut flow for 1L(H)_7-jet SR (muon channel). The signal point is taken from the gluino simplified model with 1-step decay via charginos and parameters $m_{\tilde{g}}$ = 1025 GeV, $m_{\tilde{\chi}_{1}^{\pm}}$ = 545 GeV, $m_{\tilde{\chi}_{1}^{0}}$ = 65 GeV; 20000 events were generated for this point. The events are normalized to 20.3 fb$^\mathrm{-1}$.

Cut flow for 0LRaz_SR$_{\rm loose}$ and 0LRaz_SR$_{\rm tight}$ signal regions of baseline signal point ($m(\tilde{q}, \tilde{\chi}_1^0)=(850, 100)$ GeV) normalized to 20.3 fb$^{-1}$. 5000 events were generated for ths point.

The $\delta m=m_{3\ell}-m_{ell^+\ell^-}$ distributions for the $3\ell+jj$ category and $Z+\mu$ flavor channel.

The $\delta m=m_{3\ell}-m_{ell^+\ell^-}$ distributions for the $3\ell$-only category and $Z+e$ flavor channel.

The $\delta m=m_{3\ell}-m_{ell^+\ell^-}$ distributions for the $3\ell$-only category and $Z+\mu$ flavor channel.

95\% CL upper limits on the vector-like lepton cross section, for heavy vector-like leptons coupling exclusively to electrons.

95\% CL upper limits on the vector-like lepton cross section, for heavy vector-like leptons coupling exclusively to muons.

95% CL upper limits on the type-III seesaw production cross section, for the seesaw leptons coupling exclusively to electrons.

95% CL upper limits on the type-III seesaw production cross section, for the seesaw leptons coupling exclusively to muons.

95% CL upper limits on $\sigma_{\mathrm{vis}}$ for the $Z+e$ flavour channel, derived without dividing events into the three categories.

95% CL upper limits on $\sigma_{\mathrm{vis}}$ for the $Z+\mu$ flavour channel, derived without dividing events into the three categories.

95% CL upper limits on $\sigma_{\mathrm{vis}}$ for the $Z+e$ flavour channel, for events in the four lepton signal region.:.

95% CL upper limits on $\sigma_{\mathrm{vis}}$ for the $Z+\mu$ flavour channel, for events in the four lepton signal region.

95% CL upper limits on $\sigma_{\mathrm{vis}}$ for the $Z+e$ flavour channel, for events in the two-jet signal region.

95% CL upper limits on $\sigma_{\mathrm{vis}}$ for the $Z+\mu$ flavour channel, for events in the two-jet signal region.

Efficiencies for reconstructing and correctly identifying the $L^{\pm}\rightarrow Z(\ell\ell)e^{\pm}$ decay in the inclusive category and $Z+e$ flavor channel, for events within the fiducial volume for the type-III seesaw model. One heavy lepton is required to decay to $Z+e$, with the $Z$ boson decaying leptonically. The efficiencies are shown separately for each decay mode of the second heavy lepton.

Efficiencies for reconstructing and correctly identifying the $L^{\pm}\rightarrow Z(\ell\ell)\mu^{\pm}$ decay in the inclusive category and $Z+\mu$ flavor channel, for events within the fiducial volume for the type-III seesaw model. One heavy lepton is required to decay to $Z+\mu$, with the $Z$ boson decaying leptonically. The efficiencies are shown separately for each decay mode of the second heavy lepton.

Efficiencies for reconstructing and correctly identifying the $L^{\pm}\rightarrow Z(\ell\ell)e^{\pm}$ decay in the $4\ell$ category and $Z+e$ flavor channel, for events within the fiducial volume for the type-III seesaw model. One heavy lepton is required to decay to $Z+e$, with the $Z$ boson decaying leptonically. The efficiencies are shown separately for each decay mode of the second heavy lepton.

Efficiencies for reconstructing and correctly identifying the $L^{\pm}\rightarrow Z(\ell\ell)\mu^{\pm}$ decay in the $4\ell$ category and $Z+\mu$ flavor channel, for events within the fiducial volume for the type-III seesaw model. One heavy lepton is required to decay to $Z+\mu$, with the $Z$ boson decaying leptonically. The efficiencies are shown separately for each decay mode of the second heavy lepton.

Efficiencies for reconstructing and correctly identifying the $L^{\pm}\rightarrow Z(\ell\ell)e^{\pm}$ decay in the $3\ell+jj$ category and $Z+e$ flavor channel, for events within the fiducial volume for the type-III seesaw model. One heavy lepton is required to decay to $Z+e$, with the $Z$ boson decaying leptonically. The efficiencies are shown separately for each decay mode of the second heavy lepton.

Efficiencies for reconstructing and correctly identifying the $L^{\pm}\rightarrow Z(\ell\ell)\mu^{\pm}$ decay in the $3\ell+jj$ category and $Z+\mu$ flavor channel, for events within the fiducial volume for the type-III seesaw model. One heavy lepton is required to decay to $Z+\mu$, with the $Z$ boson decaying leptonically. The efficiencies are shown separately for each decay mode of the second heavy lepton.

Efficiencies for reconstructing and correctly identifying the $L^{\pm}\rightarrow Z(\ell\ell)e^{\pm}$ decay in the $3\ell$-only category and $Z+e$ flavor channel, for events within the fiducial volume for the type-III seesaw model. One heavy lepton is required to decay to $Z+e$, with the $Z$ boson decaying leptonically. The efficiencies are shown separately for each decay mode of the second heavy lepton.

Efficiencies for reconstructing and correctly identifying the $L^{\pm}\rightarrow Z(\ell\ell)\mu^{\pm}$ decay in the $3\ell$-only category and $Z+\mu$ flavor channel, for events within the fiducial volume for the type-III seesaw model. One heavy lepton is required to decay to $Z+\mu$, with the $Z$ boson decaying leptonically. The efficiencies are shown separately for each decay mode of the second heavy lepton.

Efficiency for the calorimetric MET>80 GeV trigger as a function of the metastable R-hadron mass. The R-hadron decays to g/qq plus neutralino of mass = m(gluino) - 100 GeV with a lifetime of 1 ns.

Efficiency for the calorimetric MET>80 GeV trigger as a function of the metastable R-hadron mass. The R-hadron decays to g/qq plus neutralino of mass 100 GeV with a lifetime of 1 ns.

Efficiency for the calorimetric MET>80 GeV trigger as a function of the stable chargino mass.

Total selection efficiency as a function of the stable R-hadron mass.

Total selection efficiency as a function of the metastable R-hadron mass. The R-hadron decays to g/qq plus neutralino of mass 100 GeV with a lifetime of 10 ns.

Total selection efficiency as a function of the metastable R-hadron mass. The R-hadron decays to g/qq plus neutralino of mass = m(gluino) - 100 GeV with a lifetime of 10 ns.

Total selection efficiency as a function of the metastable R-hadron mass. The R-hadron decays to g/qq plus neutralino of mass 100 GeV with a lifetime of 1 ns.

Total selection efficiency as a function of the stable chargino mass.

Ionisation distribution of all the CR2 tracks, and those not matched to a reconstructed muon. The two distributions are normalised to their total number of entries.

Distribution of the mass of selected candidates, derived from the specific ionisation loss, for an example of gluino R-hadron signal, for searches for stable particles. The signal distributions are stacked on the expected background, and a narrower binning is used for them to allow the signal shape to be seen more clearly. The number of signal events is that expected according to the theoretical cross sections.

Distribution of the mass of selected candidates, derived from the specific ionisation loss, for one example of chargino signal, for searches for stable particles. The signal distributions are stacked on the expected background, and a narrower binning is used for them to allow the signal shape to be seen more clearly. The number of signal events is that expected according to the theoretical cross sections.

Distribution of the mass of selected candidates, derived from the specific ionisation loss, for background and data, for searches for stable particles. The expected background is shown with its total uncertainty (sum in quadrature of statistical, normalisation and systematic errors).

Distribution of the mass of selected candidates, derived from the specific ionisation loss, for an example of gluino R-hadron signal, for searches for metastable particles. The signal distributions are stacked on the expected background, and a narrower binning is used for them to allow the signal shape to be seen more clearly. The number of signal events is that expected according to the theoretical cross sections.

Distribution of the mass of selected candidates, derived from the specific ionisation loss, for an example of chargino signal, for searches for metastable particles. The signal distributions are stacked on the expected background, and a narrower binning is used for them to allow the signal shape to be seen more clearly. The number of signal events is that expected according to the theoretical cross sections.

Distribution of the mass of selected candidates, derived from the specific ionisation loss, for background and data. The expected background is shown with its total uncertainty (sum in quadrature of statistical, normalisation and systematic errors).

Theoretical values for the cross section of gluino pairs production with their uncertainty.

Expected upper limits on the production cross section as a function of mass for metastable gluino R-hadrons, with lifetime tau =10 ns, decaying into g/qq plus a light neutralino of mass 100 GeV, in the background-only case, with its 1 sigma band.

Observed 95 PCT upper limits on the production cross section as a function of mass for metastable gluino R-hadrons, with lifetime tau =10 ns, decaying into g/qq plus a light neutralino of mass 100 GeV.

Expected upper limits on the production cross section as a function of mass for metastable gluino R-hadrons, with lifetime tau =10 ns, decaying into g/qq plus a heavy neutralino of mass(gluino) - 100 GeV, in the background-only case, with its 1 sigma band.

Observed 95 PCT upper limits on the production cross section as a function of mass for metastable gluino R-hadrons, with lifetime tau =10 ns, decaying into g/qq plus a heavy neutralino of mass(gluino) - 100 GeV.

The expected excluded range of lifetimes as a function of gluino mass for gluino R-hadrons decaying into g/qq plus a light neutralino of mass 100 GeV, with respect to the nominal theoretical cross section.

The expected excluded range of lifetimes as a function of gluino mass for gluino R-hadrons decaying into g/qq plus a light neutralino of mass 100 GeV, with respect to the nominal theoretical cross section, plus 1 experimental sigma.

The expected excluded range of lifetimes as a function of gluino mass for gluino R-hadrons decaying into g/qq plus a light neutralino of mass 100 GeV, with respect to the nominal theoretical cross section, minus 1 experimental sigma.

The observed excluded range of lifetimes as a function of gluino mass for gluino R-hadrons decaying into g/qq plus a light neutralino of mass 100 GeV, with respect to the nominal theoretical cross section.

The observed excluded range of lifetimes as a function of gluino mass for gluino R-hadrons decaying into g/qq plus a light neutralino of mass 100 GeV, with respect to the nominal theoretical cross section minus its uncertainty.

The observed excluded range of lifetimes as a function of gluino mass for gluino R-hadrons decaying into g/qq plus a light neutralino of mass 100 GeV, with respect to the nominal theoretical cross section plus its uncertainty.

The expected excluded range of lifetimes as a function of gluino mass for gluino R-hadrons decaying into g/qq plus a heavy neutralino of mass = m(gluino) - 100 GeV, with respect to the nominal theoretical cross section.

The expected excluded range of lifetimes as a function of gluino mass for gluino R-hadrons decaying into g/qq plus a heavy neutralino of mass = m(gluino) - 100 GeV, with respect to the nominal theoretical cross section, plus 1 experimental sigma.

The expected excluded range of lifetimes as a function of gluino mass for gluino R-hadrons decaying into g/qq plus a heavy neutralino of mass = m(gluino) - 100 GeV, with respect to the nominal theoretical cross section, minus 1 experimental sigma.

The observed excluded range of lifetimes as a function of gluino mass for gluino R-hadrons decaying into g/qq plus a heavy neutralino of mass = m(gluino) - 100 GeV, with respect to the nominal theoretical cross section.

The observed excluded range of lifetimes as a function of gluino mass for gluino R-hadrons decaying into g/qq plus a heavy neutralino of mass = m(gluino) - 100 GeV, with respect to the nominal theoretical cross section minus its uncertainty.

The observed excluded range of lifetimes as a function of gluino mass for gluino R-hadrons decaying into g/qq plus a heavy neutralino of mass = m(gluino) - 100 GeV, with respect to the nominal theoretical cross section plus its uncertainty.

Expected upper limits on the production cross section as a function of mass for metastable gluino R-hadrons, with lifetime tau =10 ns, decaying into tt plus a light neutralino of mass 100 GeV, in the background-only case, with its 1 sigma band.

Observed 95 PCT upper limits on the production cross section as a function of mass for metastable gluino R-hadrons, with lifetime tau =10 ns, decaying into tt plus a light neutralino of mass 100 GeV.

Expected upper limits on the production cross section as a function of mass for metastable gluino R-hadrons, with lifetime tau =10 ns, decaying into tt plus a heavy neutralino of mass(gluino) - 100 GeV, in the background-only case, with its 1 sigma band.

Observed 95 PCT upper limits on the production cross section as a function of mass for metastable gluino R-hadrons, with lifetime tau =10 ns, decaying into tt plus a heavy neutralino of mass(gluino) - 100 GeV.

The expected excluded range of lifetimes as a function of gluino mass for gluino R-hadrons decaying into tt plus a light neutralino of mass 100 GeV, with respect to the nominal theoretical cross section.

The expected excluded range of lifetimes as a function of gluino mass for gluino R-hadrons decaying into tt plus a light neutralino of mass 100 GeV, with respect to the nominal theoretical cross section, plus 1 experimental sigma.

The expected excluded range of lifetimes as a function of gluino mass for gluino R-hadrons decaying into tt plus a light neutralino of mass 100 GeV, with respect to the nominal theoretical cross section, minus 1 experimental sigma.

The observed excluded range of lifetimes as a function of gluino mass for gluino R-hadrons decaying into tt plus a light neutralino of mass 100 GeV, with respect to the nominal theoretical cross section.

The observed excluded range of lifetimes as a function of gluino mass for gluino R-hadrons decaying into tt plus a light neutralino of mass 100 GeV, with respect to the nominal theoretical cross section minus its uncertainty.

The observed excluded range of lifetimes as a function of gluino mass for gluino R-hadrons decaying into tt plus a light neutralino of mass 100 GeV, with respect to the nominal theoretical cross section plus its uncertainty.

The expected excluded range of lifetimes as a function of gluino mass for gluino R-hadrons decaying into tt plus a heavy neutralino of mass = m(gluino) - 100 GeV, with respect to the nominal theoretical cross section.

The expected excluded range of lifetimes as a function of gluino mass for gluino R-hadrons decaying into tt plus a heavy neutralino of mass = m(gluino) - 100 GeV, with respect to the nominal theoretical cross section, plus 1 experimental sigma.

The expected excluded range of lifetimes as a function of gluino mass for gluino R-hadrons decaying into tt plus a heavy neutralino of mass = m(gluino) - 100 GeV, with respect to the nominal theoretical cross section, minus 1 experimental sigma.

The observed excluded range of lifetimes as a function of gluino mass for gluino R-hadrons decaying into tt plus a heavy neutralino of mass = m(gluino) - 100 GeV, with respect to the nominal theoretical cross section.

The observed excluded range of lifetimes as a function of gluino mass for gluino R-hadrons decaying into tt plus a heavy neutralino of mass = m(gluino) - 100 GeV, with respect to the nominal theoretical cross section minus its uncertainty.

The observed excluded range of lifetimes as a function of gluino mass for gluino R-hadrons decaying into tt plus a heavy neutralino of mass = m(gluino) - 100 GeV, with respect to the nominal theoretical cross section plus its uncertainty.

Theoretical values for the production cross section of charginos or chargino/neutralino pairs, with their uncertainty.

Expected upper limits on the production cross section as a function of mass for metastable charginos, with lifetime tau =1.0 ns, decaying into neutralino + pion, in the background-only case, with its 1 sigma band.

Observed 95 PCT upper limits on the production cross section as a function of mass for metastable charginos, with lifetime tau =1.0 ns, decaying into neutralino + pion.

The expected excluded range of lifetimes as a function of chargino mass for charginos decaying into neutralino plus pion, with respect to the nominal theoretical cross section.

The expected excluded range of lifetimes as a function of chargino mass for charginos decaying into neutralino plus pion, with respect to the nominal theoretical cross section, plus 1 experimental sigma.

The expected excluded range of lifetimes as a function of chargino mass for charginos decaying into neutralino plus pion, with respect to the nominal theoretical cross section, minus 1 experimental sigma.

The observed excluded range of lifetimes as a function of chargino mass for charginos decaying into neutralino plus pion, with respect to the nominal theoretical cross section.

The observed excluded range of lifetimes as a function of chargino mass for charginos decaying into neutralino plus pion, with respect to the nominal theoretical cross section minus its uncertainty.

The observed excluded range of lifetimes as a function of gluino mass for chargino mass for charginos decaying into neutralino plus pion, with respect to the nominal theoretical cross section plus its uncertainty.

dEdx ionization for data, 1 TeV gluino R-hadrons stable and decaying in 100 GeV neutralinos with a 10 ns lifetime and for charginos of 350 GeV. Tracks that fulfil all the requirements up to including the High-m_T (see Tab.1 in the paper) are considered at this stage and normalised to an integrated luminosity of 18.4 fb^-1.

Expected upper limits on the production cross section as a function of mass for stable gluino R-hadrons, in case of background only, with its 1 sigma band.

Observed 95 PCT upper limits on the production cross section as a function of mass for stable gluino R-hadrons.

Theoretical values for the cross section of squark pairs production with their uncertainty.

Expected upper limits on the production cross section as a function of mass for stable sbottom R-hadrons, in case of background only, with its 1 sigma band.

Observed 95 PCT upper limits on the production cross section as a function of mass for stable sbottom $R$-hadrons. Cross section IN PB.

Expected upper limits on the production cross section as a function of mass for stop R-hadrons, in case of background only, with its 1 sigma band.

Observed 95 PCT upper limits on the production cross section as a function of mass for stop R-hadrons.

Expected upper limits on the production cross section as a function of mass for stable charginos, in case of background only, with its 1 sigma band.

Observed 95 PCT upper limits on the production cross section as a function of mass for stable charginos.

Expected upper limits on the production cross section as a function of mass for metastable gluino R-hadrons, with lifetime tau=1.0 ns, decaying to g/qq plus a light neutralino of mass 100 GeV.

Observed 95 PCT upper limits on the production cross section as a function of mass for metastable gluino R-hadrons, with lifetime tau=1.0 ns, decaying to g/qq plus a light neutralino of mass 100 GeV.

Expected upper limits on the production cross section as a function of mass for metastable gluino R-hadrons, with lifetime tau=1.0 ns, decaying to g/qq plus a heavy neutralino of mass = m(gluino) - 100 GeV.

Observed 95 PCT upper limits on the production cross section as a function of mass for metastable gluino R-hadrons, with lifetime tau=1.0 ns, decaying to g/qq plus a heavy neutralino of mass = m(gluino) - 100 GeV.

Expected upper limits on the production cross section as a function of mass for metastable gluino R-hadrons, with lifetime tau=1.0 ns, decaying to tt plus a light neutralino of mass 100 GeV.

Observed 95 PCT upper limits on the production cross section as a function of mass for metastable gluino R-hadrons, with lifetime tau=1.0 ns, decaying to tt plus a light neutralino of mass 100 GeV.

Expected upper limits on the production cross section as a function of mass for metastable gluino R-hadrons, with lifetime tau=1.0 ns, decaying to tt plus a heavy neutralino of mass = m(gluino) - 100 GeV.

Observed 95 PCT upper limits on the production cross section as a function of mass for metastable gluino R-hadrons, with lifetime tau=1.0 ns, decaying to tt plus a heavy neutralino of mass = m(gluino) - 100 GeV.

Expected upper limits on the production cross section as a function of mass for metastable charginos, with lifetime tau =15 ns, decaying to neutralino and pion, in case of background only, with its 1 sigma band.

Observed 95 PCT upper limits on the production cross section as a function of mass for metastable charginos, with lifetime tau =15 ns, decaying to neutralino and pion, in case of background only, with its 1 sigma band.

HT inclusive distribution in SRZ.

Jet multiplicity distribution in SRZ.

B-tagged jet multiplicity distribution in SRZ.

Dilepton invariant mass distribution in SRMedium.

Dilepton invariant mass distribution in SRMedium.

Dilepton invariant mass in SRHigh.

Expected 95% exclusion contour for the gluino simplified model with fixed N1 mass and decays to quarks and Z bosons in SRZ.

Observed 95% exclusion contour for the gluino simplified model with fixed N1 mass and decays to quarks and Z bosons in SRZ.

Expected 95% exclusion contour for the squark simplified model with fixed N1 mass and decays to quarks and Z bosons in SRZ.

Observed 95% exclusion contour for the squark simplified model with fixed N1 mass and decays to quarks and Z bosons in SRZ.

Expected 95% exclusion contour for the gluino simplified model with varying N1 mass and decays to quarks and Z bosons in SRZ.

Observed 95% exclusion contour for the gluino simplified model with varying N1 mass and decays to quarks and Z bosons in SRZ.

Expected 95% exclusion contour for the gluino simplified model with sleptons in the edge SRs.

Observed 95% exclusion contour for the gluino simplified model with sleptons in the edge SRs.

Expected 95% exclusion contour for the gluino simplified model with decays to on- or off-shell Z bosons in the edge SRs.

Observed 95% exclusion contour for the gluino simplified model with decays to on- or off-shell Z bosons in the edge SRs.

Expected 95% exclusion contour for the gluino simplified model with decays to on- or off-shell Z bosons in SRZ.

Observed 95% exclusion contour for the gluino simplified model with decays to on- or off-shell Z bosons in SRZ.

Cutflow table for three benchmark signal points from the $\tilde{g}$--$\chi_{2}^{0}$ on-shell grid in SRZ for the $ee$ and $\mu\mu$ channels separately. The number of generated events is N RAW, while N WEIGHTED corresponds to the number of events expected in 14.7/fb, after all MC simulation corrections have been applied.

Cutflow table for three benchmark signal points for the slepton model in SR-low, SR-medium and SR-high for the $ee$ and $\mu\mu$ channels separately. The number of generated events is N RAW, while N WEIGHTED corresponds to the number of events expected in 14.7/fb, after all MC simulation corrections have been applied. The event generation was filtered to require at least two leptons (electrons or muons) with $p_{T}>5$ GeV and $|\eta|<2.8$, with 20000 events required to pass the filter. The corresponding number of unfiltered events is provided in the ``Pre-filter events'' row of the table.

Signal acceptance for the gluino simplified model with decays to quarks and Z bosons in SRZ.

Signal efficiency for the gluino simplified model with decays to quarks and Z bosons in SRZ.

Signal acceptance for the gluino simplified model with decays via sleptons in SRLow.

Signal efficiency for the gluino simplified model with decays via sleptons in SRLow.

Signal acceptance for the gluino simplified model with decays via sleptons in SRMedium.

Signal efficiency for the gluino simplified model with decays via sleptons in SRMedium.

Signal acceptance for the gluino simplified model with decays via sleptons in SRHigh.

Signal efficiency for the gluino simplified model with decays via sleptons in SRHigh.

Upper limits on the signal cross-section at 95% CL for the gluino simplified model with fixed N1 mass and decays to quarks and Z bosons in SRZ.

Upper limits on the signal cross-section at 95% CL for the squark simplified model with fixed N1 mass and decays to quarks and Z bosons in SRZ.

Upper limits on the signal cross-section at 95% CL for the gluino simplified model with varying N1 mass and decays to quarks and Z bosons in SRZ.

Upper limits on the signal cross-section at 95% CL for the gluino simplified model with decays via sleptons.

Upper limits on the signal cross-section at 95% CL for the gluino simplified model with decays to on-shell or off-shell Z bosons.