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

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.

The cross section times BR(A->ZH) times BR(H->bb) limits for an A boson in the type-I, tanbeta 1, 2HDM. For each signal point, characterised by the mass pair (mA, mH), two limits are provided, the observed and the expected. The correct width as predicted by this particular parameter choice of the 2HDM is used. For type-I only gluon-gluon fusion is used and the nb=2 category only.

The cross section times BR(A->ZH) times BR(H->bb) limits for an A boson in the type-I, tanbeta 5, 2HDM. For each signal point, characterised by the mass pair (mA, mH), two limits are provided, the observed and the expected. The correct width as predicted by this particular parameter choice of the 2HDM is used. For type-I only gluon-gluon fusion is used and the nb=2 category only.

The cross section times BR(A->ZH) times BR(H->bb) limits for an A boson in the type-I, tanbeta 10, 2HDM. For each signal point, characterised by the mass pair (mA, mH), two limits are provided, the observed and the expected. The correct width as predicted by this particular parameter choice of the 2HDM is used. For type-I only gluon-gluon fusion is used and the nb=2 category only.

The cross section times BR(A->ZH) times BR(H->bb) limits for an A boson in the type-I, tanbeta 20, 2HDM. For each signal point, characterised by the mass pair (mA, mH), two limits are provided, the observed and the expected. The correct width as predicted by this particular parameter choice of the 2HDM is used. For type-I only gluon-gluon fusion is used and the nb=2 category only.

The cross section times BR(A->ZH) times BR(H->bb) limits for an A boson in the type-II, tanbeta 1, 2HDM. For each signal point, characterised by the mass pair (mA, mH), two limits are provided, the observed and the expected. The correct width as predicted by this particular parameter choice of the 2HDM is used. For type-II both gluon-gluon fusion and b-associated production are used and the nb=2 and nb>=3 categories are combined.

The cross section times BR(A->ZH) times BR(H->bb) limits for an A boson in the type-II, tanbeta 5, 2HDM. For each signal point, characterised by the mass pair (mA, mH), two limits are provided, the observed and the expected. The correct width as predicted by this particular parameter choice of the 2HDM is used. For type-II both gluon-gluon fusion and b-associated production are used and the nb=2 and nb>=3 categories are combined.

The cross section times BR(A->ZH) times BR(H->bb) limits for an A boson in the type-II, tanbeta 10, 2HDM. For each signal point, characterised by the mass pair (mA, mH), two limits are provided, the observed and the expected. The correct width as predicted by this particular parameter choice of the 2HDM is used. For type-II both gluon-gluon fusion and b-associated production are used and the nb=2 and nb>=3 categories are combined.

The cross section times BR(A->ZH) times BR(H->bb) limits for an A boson in the type-II, tanbeta 20, 2HDM. For each signal point, characterised by the mass pair (mA, mH), two limits are provided, the observed and the expected. The correct width as predicted by this particular parameter choice of the 2HDM is used. For type-II both gluon-gluon fusion and b-associated production are used and the nb=2 and nb>=3 categories are combined.

The cross section times BR(A->ZH) times BR(H->bb) limits for an A boson in the lepton specific, tanbeta 1, 2HDM. For each signal point, characterised by the mass pair (mA, mH), two limits are provided, the observed and the expected. The correct width as predicted by this particular parameter choice of the 2HDM is used. For lepton specific only gluon-gluon fusion is used and the nb=2 category only.

The cross section times BR(A->ZH) times BR(H->bb) limits for an A boson in the lepton specific, tanbeta 2, 2HDM. For each signal point, characterised by the mass pair (mA, mH), two limits are provided, the observed and the expected. The correct width as predicted by this particular parameter choice of the 2HDM is used. For lepton specific only gluon-gluon fusion is used and the nb=2 category only.

The cross section times BR(A->ZH) times BR(H->bb) limits for an A boson in the lepton specific, tanbeta 3, 2HDM. For each signal point, characterised by the mass pair (mA, mH), two limits are provided, the observed and the expected. The correct width as predicted by this particular parameter choice of the 2HDM is used. For lepton specific only gluon-gluon fusion is used and the nb=2 category only.

The cross section times BR(A->ZH) times BR(H->bb) limits for an A boson in the flipped, tanbeta 1, 2HDM. For each signal point, characterised by the mass pair (mA, mH), two limits are provided, the observed and the expected. The correct width as predicted by this particular parameter choice of the 2HDM is used. For flipped both gluon-gluon fusion and b-associated production are used and the nb=2 and nb>=3 categories are combined.

The cross section times BR(A->ZH) times BR(H->bb) limits for an A boson in the flipped, tanbeta 5, 2HDM. For each signal point, characterised by the mass pair (mA, mH), two limits are provided, the observed and the expected. The correct width as predicted by this particular parameter choice of the 2HDM is used. For flipped both gluon-gluon fusion and b-associated production are used and the nb=2 and nb>=3 categories are combined.

The cross section times BR(A->ZH) times BR(H->bb) limits for an A boson in the flipped, tanbeta 10, 2HDM. For each signal point, characterised by the mass pair (mA, mH), two limits are provided, the observed and the expected. The correct width as predicted by this particular parameter choice of the 2HDM is used. For flipped both gluon-gluon fusion and b-associated production are used and the nb=2 and nb>=3 categories are combined.

The cross section times BR(A->ZH) times BR(H->bb) limits for an A boson in the flipped, tanbeta 20, 2HDM. For each signal point, characterised by the mass pair (mA, mH), two limits are provided, the observed and the expected. The correct width as predicted by this particular parameter choice of the 2HDM is used. For flipped both gluon-gluon fusion and b-associated production are used and the nb=2 and nb>=3 categories are combined.

The mass distribution of the llbb system used in the fit to derive the limits for the mbb window centered at 130 GeV for the nb = 2 (2 tag) category.

The mass distribution of the llbb system used in the fit to derive the limits for the mbb window centered at 130 GeV for the nb >= 3 (3 tag) category.

The mass distribution of the llbb system used in the fit to derive the limits for the mbb window centered at 140 GeV for the nb = 2 (2 tag) category.

The mass distribution of the llbb system used in the fit to derive the limits for the mbb window centered at 140 GeV for the nb >= 3 (3 tag) category.

The mass distribution of the llbb system used in the fit to derive the limits for the mbb window centered at 150 GeV for the nb = 2 (2 tag) category.

The mass distribution of the llbb system used in the fit to derive the limits for the mbb window centered at 150 GeV for the nb >= 3 (3 tag) category.

The mass distribution of the llbb system used in the fit to derive the limits for the mbb window centered at 160 GeV for the nb = 2 (2 tag) category.

The mass distribution of the llbb system used in the fit to derive the limits for the mbb window centered at 160 GeV for the nb >= 3 (3 tag) category.

The mass distribution of the llbb system used in the fit to derive the limits for the mbb window centered at 170 GeV for the nb = 2 (2 tag) category.

The mass distribution of the llbb system used in the fit to derive the limits for the mbb window centered at 170 GeV for the nb >= 3 (3 tag) category.

The mass distribution of the llbb system used in the fit to derive the limits for the mbb window centered at 180 GeV for the nb = 2 (2 tag) category.

The mass distribution of the llbb system used in the fit to derive the limits for the mbb window centered at 180 GeV for the nb >= 3 (3 tag) category.

The mass distribution of the llbb system used in the fit to derive the limits for the mbb window centered at 190 GeV for the nb = 2 (2 tag) category.

The mass distribution of the llbb system used in the fit to derive the limits for the mbb window centered at 190 GeV for the nb >= 3 (3 tag) category.

The mass distribution of the llbb system used in the fit to derive the limits for the mbb window centered at 210 GeV for the nb = 2 (2 tag) category.

The mass distribution of the llbb system used in the fit to derive the limits for the mbb window centered at 210 GeV for the nb >= 3 (3 tag) category.

The mass distribution of the llbb system used in the fit to derive the limits for the mbb window centered at 220 GeV for the nb = 2 (2 tag) category.

The mass distribution of the llbb system used in the fit to derive the limits for the mbb window centered at 220 GeV for the nb >= 3 (3 tag) category.

The mass distribution of the llbb system used in the fit to derive the limits for the mbb window centered at 230 GeV for the nb = 2 (2 tag) category.

The mass distribution of the llbb system used in the fit to derive the limits for the mbb window centered at 230 GeV for the nb >= 3 (3 tag) category.

The mass distribution of the llbb system used in the fit to derive the limits for the mbb window centered at 240 GeV for the nb = 2 (2 tag) category.

The mass distribution of the llbb system used in the fit to derive the limits for the mbb window centered at 240 GeV for the nb >= 3 (3 tag) category.

The mass distribution of the llbb system used in the fit to derive the limits for the mbb window centered at 250 GeV for the nb = 2 (2 tag) category.

The mass distribution of the llbb system used in the fit to derive the limits for the mbb window centered at 250 GeV for the nb >= 3 (3 tag) category.

The mass distribution of the llbb system used in the fit to derive the limits for the mbb window centered at 260 GeV for the nb = 2 (2 tag) category.

The mass distribution of the llbb system used in the fit to derive the limits for the mbb window centered at 260 GeV for the nb >= 3 (3 tag) category.

The mass distribution of the llbb system used in the fit to derive the limits for the mbb window centered at 270 GeV for the nb = 2 (2 tag) category.

The mass distribution of the llbb system used in the fit to derive the limits for the mbb window centered at 270 GeV for the nb >= 3 (3 tag) category.

The mass distribution of the llbb system used in the fit to derive the limits for the mbb window centered at 280 GeV for the nb = 2 (2 tag) category.

The mass distribution of the llbb system used in the fit to derive the limits for the mbb window centered at 280 GeV for the nb >= 3 (3 tag) category.

The mass distribution of the llbb system used in the fit to derive the limits for the mbb window centered at 290 GeV for the nb = 2 (2 tag) category.

The mass distribution of the llbb system used in the fit to derive the limits for the mbb window centered at 290 GeV for the nb >= 3 (3 tag) category.

The mass distribution of the llbb system used in the fit to derive the limits for the mbb window centered at 310 GeV for the nb = 2 (2 tag) category.

The mass distribution of the llbb system used in the fit to derive the limits for the mbb window centered at 310 GeV for the nb >= 3 (3 tag) category.

The mass distribution of the llbb system used in the fit to derive the limits for the mbb window centered at 320 GeV for the nb = 2 (2 tag) category.

The mass distribution of the llbb system used in the fit to derive the limits for the mbb window centered at 320 GeV for the nb >= 3 (3 tag) category.

The mass distribution of the llbb system used in the fit to derive the limits for the mbb window centered at 330 GeV for the nb = 2 (2 tag) category.

The mass distribution of the llbb system used in the fit to derive the limits for the mbb window centered at 330 GeV for the nb >= 3 (3 tag) category.

The mass distribution of the llbb system used in the fit to derive the limits for the mbb window centered at 340 GeV for the nb = 2 (2 tag) category.

The mass distribution of the llbb system used in the fit to derive the limits for the mbb window centered at 340 GeV for the nb >= 3 (3 tag) category.

The mass distribution of the llbb system used in the fit to derive the limits for the mbb window centered at 350 GeV for the nb = 2 (2 tag) category.

The mass distribution of the llbb system used in the fit to derive the limits for the mbb window centered at 350 GeV for the nb >= 3 (3 tag) category.

The mass distribution of the llbb system used in the fit to derive the limits for the mbb window centered at 360 GeV for the nb = 2 (2 tag) category.

The mass distribution of the llbb system used in the fit to derive the limits for the mbb window centered at 360 GeV for the nb >= 3 (3 tag) category.

The mass distribution of the llbb system used in the fit to derive the limits for the mbb window centered at 370 GeV for the nb = 2 (2 tag) category.

The mass distribution of the llbb system used in the fit to derive the limits for the mbb window centered at 370 GeV for the nb >= 3 (3 tag) category.

The mass distribution of the llbb system used in the fit to derive the limits for the mbb window centered at 380 GeV for the nb = 2 (2 tag) category.

The mass distribution of the llbb system used in the fit to derive the limits for the mbb window centered at 380 GeV for the nb >= 3 (3 tag) category.

The mass distribution of the llbb system used in the fit to derive the limits for the mbb window centered at 390 GeV for the nb = 2 (2 tag) category.

The mass distribution of the llbb system used in the fit to derive the limits for the mbb window centered at 390 GeV for the nb >= 3 (3 tag) category.

The mass distribution of the llbb system used in the fit to derive the limits for the mbb window centered at 400 GeV for the nb = 2 (2 tag) category.

The mass distribution of the llbb system used in the fit to derive the limits for the mbb window centered at 400 GeV for the nb >= 3 (3 tag) category.

The mass distribution of the llbb system used in the fit to derive the limits for the mbb window centered at 410 GeV for the nb = 2 (2 tag) category.

The mass distribution of the llbb system used in the fit to derive the limits for the mbb window centered at 410 GeV for the nb >= 3 (3 tag) category.

The mass distribution of the llbb system used in the fit to derive the limits for the mbb window centered at 420 GeV for the nb = 2 (2 tag) category.

The mass distribution of the llbb system used in the fit to derive the limits for the mbb window centered at 420 GeV for the nb >= 3 (3 tag) category.

The mass distribution of the llbb system used in the fit to derive the limits for the mbb window centered at 430 GeV for the nb = 2 (2 tag) category.

The mass distribution of the llbb system used in the fit to derive the limits for the mbb window centered at 430 GeV for the nb >= 3 (3 tag) category.

The mass distribution of the llbb system used in the fit to derive the limits for the mbb window centered at 440 GeV for the nb = 2 (2 tag) category.

The mass distribution of the llbb system used in the fit to derive the limits for the mbb window centered at 440 GeV for the nb >= 3 (3 tag) category.

The mass distribution of the llbb system used in the fit to derive the limits for the mbb window centered at 450 GeV for the nb = 2 (2 tag) category.

The mass distribution of the llbb system used in the fit to derive the limits for the mbb window centered at 450 GeV for the nb >= 3 (3 tag) category.

The mass distribution of the llbb system used in the fit to derive the limits for the mbb window centered at 460 GeV for the nb = 2 (2 tag) category.

The mass distribution of the llbb system used in the fit to derive the limits for the mbb window centered at 460 GeV for the nb >= 3 (3 tag) category.

The mass distribution of the llbb system used in the fit to derive the limits for the mbb window centered at 470 GeV for the nb = 2 (2 tag) category.

The mass distribution of the llbb system used in the fit to derive the limits for the mbb window centered at 470 GeV for the nb >= 3 (3 tag) category.

The mass distribution of the llbb system used in the fit to derive the limits for the mbb window centered at 480 GeV for the nb = 2 (2 tag) category.

The mass distribution of the llbb system used in the fit to derive the limits for the mbb window centered at 480 GeV for the nb >= 3 (3 tag) category.

The mass distribution of the llbb system used in the fit to derive the limits for the mbb window centered at 490 GeV for the nb = 2 (2 tag) category.

The mass distribution of the llbb system used in the fit to derive the limits for the mbb window centered at 490 GeV for the nb >= 3 (3 tag) category.

The mass distribution of the llbb system used in the fit to derive the limits for the mbb window centered at 510 GeV for the nb = 2 (2 tag) category.

The mass distribution of the llbb system used in the fit to derive the limits for the mbb window centered at 510 GeV for the nb >= 3 (3 tag) category.

The mass distribution of the llbb system used in the fit to derive the limits for the mbb window centered at 520 GeV for the nb = 2 (2 tag) category.

The mass distribution of the llbb system used in the fit to derive the limits for the mbb window centered at 520 GeV for the nb >= 3 (3 tag) category.

The mass distribution of the llbb system used in the fit to derive the limits for the mbb window centered at 530 GeV for the nb = 2 (2 tag) category.

The mass distribution of the llbb system used in the fit to derive the limits for the mbb window centered at 530 GeV for the nb >= 3 (3 tag) category.

The mass distribution of the llbb system used in the fit to derive the limits for the mbb window centered at 540 GeV for the nb = 2 (2 tag) category.

The mass distribution of the llbb system used in the fit to derive the limits for the mbb window centered at 540 GeV for the nb >= 3 (3 tag) category.

The mass distribution of the llbb system used in the fit to derive the limits for the mbb window centered at 550 GeV for the nb = 2 (2 tag) category.

The mass distribution of the llbb system used in the fit to derive the limits for the mbb window centered at 550 GeV for the nb >= 3 (3 tag) category.

The mass distribution of the llbb system used in the fit to derive the limits for the mbb window centered at 560 GeV for the nb = 2 (2 tag) category.

The mass distribution of the llbb system used in the fit to derive the limits for the mbb window centered at 560 GeV for the nb >= 3 (3 tag) category.

The mass distribution of the llbb system used in the fit to derive the limits for the mbb window centered at 570 GeV for the nb = 2 (2 tag) category.

The mass distribution of the llbb system used in the fit to derive the limits for the mbb window centered at 570 GeV for the nb >= 3 (3 tag) category.

The mass distribution of the llbb system used in the fit to derive the limits for the mbb window centered at 580 GeV for the nb = 2 (2 tag) category.

The mass distribution of the llbb system used in the fit to derive the limits for the mbb window centered at 580 GeV for the nb >= 3 (3 tag) category.

The mass distribution of the llbb system used in the fit to derive the limits for the mbb window centered at 590 GeV for the nb = 2 (2 tag) category.

The mass distribution of the llbb system used in the fit to derive the limits for the mbb window centered at 590 GeV for the nb >= 3 (3 tag) category.

The mass distribution of the llbb system used in the fit to derive the limits for the mbb window centered at 600 GeV for the nb = 2 (2 tag) category.

The mass distribution of the llbb system used in the fit to derive the limits for the mbb window centered at 600 GeV for the nb >= 3 (3 tag) category.

The mass distribution of the llbb system used in the fit to derive the limits for the mbb window centered at 610 GeV for the nb = 2 (2 tag) category.

The mass distribution of the llbb system used in the fit to derive the limits for the mbb window centered at 610 GeV for the nb >= 3 (3 tag) category.

The mass distribution of the llbb system used in the fit to derive the limits for the mbb window centered at 620 GeV for the nb = 2 (2 tag) category.

The mass distribution of the llbb system used in the fit to derive the limits for the mbb window centered at 620 GeV for the nb >= 3 (3 tag) category.

The mass distribution of the llbb system used in the fit to derive the limits for the mbb window centered at 630 GeV for the nb = 2 (2 tag) category.

The mass distribution of the llbb system used in the fit to derive the limits for the mbb window centered at 630 GeV for the nb >= 3 (3 tag) category.

The mass distribution of the llbb system used in the fit to derive the limits for the mbb window centered at 640 GeV for the nb = 2 (2 tag) category.

The mass distribution of the llbb system used in the fit to derive the limits for the mbb window centered at 640 GeV for the nb >= 3 (3 tag) category.

The mass distribution of the llbb system used in the fit to derive the limits for the mbb window centered at 650 GeV for the nb = 2 (2 tag) category.

The mass distribution of the llbb system used in the fit to derive the limits for the mbb window centered at 650 GeV for the nb >= 3 (3 tag) category.

The mass distribution of the llbb system used in the fit to derive the limits for the mbb window centered at 660 GeV for the nb = 2 (2 tag) category.

The mass distribution of the llbb system used in the fit to derive the limits for the mbb window centered at 660 GeV for the nb >= 3 (3 tag) category.

The mass distribution of the llbb system used in the fit to derive the limits for the mbb window centered at 670 GeV for the nb = 2 (2 tag) category.

The mass distribution of the llbb system used in the fit to derive the limits for the mbb window centered at 670 GeV for the nb >= 3 (3 tag) category.

The mass distribution of the llbb system used in the fit to derive the limits for the mbb window centered at 680 GeV for the nb = 2 (2 tag) category.

The mass distribution of the llbb system used in the fit to derive the limits for the mbb window centered at 680 GeV for the nb >= 3 (3 tag) category.

The mass distribution of the llbb system used in the fit to derive the limits for the mbb window centered at 690 GeV for the nb = 2 (2 tag) category.

The mass distribution of the llbb system used in the fit to derive the limits for the mbb window centered at 690 GeV for the nb >= 3 (3 tag) category.

The mass distribution of the llbb system used in the fit to derive the limits for the mbb window centered at 700 GeV for the nb = 2 (2 tag) category.

The mass distribution of the llbb system used in the fit to derive the limits for the mbb window centered at 700 GeV for the nb >= 3 (3 tag) category.

The mass distribution of the llbb system used in the fit to derive the limits for the mbb window centered at 200 GeV for the nb = 2 (2 tag) category.

The mass distribution of the llbb system used in the fit to derive the limits for the mbb window centered at 200 GeV for the nb >= 3 (3 tag) category.

The mass distribution of the llbb system used in the fit to derive the limits for the mbb window centered at 300 GeV for the nb = 2 (2 tag) category.

The mass distribution of the llbb system used in the fit to derive the limits for the mbb window centered at 300 GeV for the nb >= 3 (3 tag) category.

The mass distribution of the llbb system used in the fit to derive the limits for the mbb window centered at 500 GeV for the nb = 2 (2 tag) category.

The mass distribution of the llbb system used in the fit to derive the limits for the mbb window centered at 500 GeV for the nb >= 3 (3 tag) category.

Measured differential cross section with associated uncertainties as a function of COS(THETA*). Each systematic uncertainty sources is fully uncorrelated with the other sources and fully correlated across bins, except for the background modelling systematics for which an uncorrelated treatment across bins is more appropriate.

Measured differential cross section with associated uncertainties as a function of DELTAYRAP(2GAMMA). Each systematic uncertainty sources is fully uncorrelated with the other sources and fully correlated across bins, except for the background modelling systematics for which an uncorrelated treatment across bins is more appropriate.

Measured differential cross section with associated uncertainties as a function of MULT(JET,PT>30 GEV). Each systematic uncertainty sources is fully uncorrelated with the other sources and fully correlated across bins, except for the background modelling systematics for which an uncorrelated treatment across bins is more appropriate.

Measured differential cross section with associated uncertainties as a function of MULT(JET,PT>50 GEV). Each systematic uncertainty sources is fully uncorrelated with the other sources and fully correlated across bins, except for the background modelling systematics for which an uncorrelated treatment across bins is more appropriate.

Measured differential cross section with associated uncertainties as a function of PT(JET1). Each systematic uncertainty sources is fully uncorrelated with the other sources and fully correlated across bins, except for the background modelling systematics for which an uncorrelated treatment across bins is more appropriate.

Measured differential cross section with associated uncertainties as a function of PT(JET2). Each systematic uncertainty sources is fully uncorrelated with the other sources and fully correlated across bins, except for the background modelling systematics for which an uncorrelated treatment across bins is more appropriate.

Measured differential cross section with associated uncertainties as a function of HT. Each systematic uncertainty sources is fully uncorrelated with the other sources and fully correlated across bins, except for the background modelling systematics for which an uncorrelated treatment across bins is more appropriate.

Measured differential cross section with associated uncertainties as a function of YRAP(JET1). Each systematic uncertainty sources is fully uncorrelated with the other sources and fully correlated across bins, except for the background modelling systematics for which an uncorrelated treatment across bins is more appropriate.

Measured differential cross section with associated uncertainties as a function of YRAP(JET2). Each systematic uncertainty sources is fully uncorrelated with the other sources and fully correlated across bins, except for the background modelling systematics for which an uncorrelated treatment across bins is more appropriate.

Measured differential cross section with associated uncertainties as a function of M(2JET). Each systematic uncertainty sources is fully uncorrelated with the other sources and fully correlated across bins, except for the background modelling systematics for which an uncorrelated treatment across bins is more appropriate.

Measured differential cross section with associated uncertainties as a function of DELTAYRAP(2JET). Each systematic uncertainty sources is fully uncorrelated with the other sources and fully correlated across bins, except for the background modelling systematics for which an uncorrelated treatment across bins is more appropriate.

Measured differential cross section with associated uncertainties as a function of ABSDPHI(2JET). Each systematic uncertainty sources is fully uncorrelated with the other sources and fully correlated across bins, except for the background modelling systematics for which an uncorrelated treatment across bins is more appropriate.

Measured differential cross section with associated uncertainties as a function of DPHI(2JET). Each systematic uncertainty sources is fully uncorrelated with the other sources and fully correlated across bins, except for the background modelling systematics for which an uncorrelated treatment across bins is more appropriate.

Measured differential cross section with associated uncertainties as a function of PT(2GAMMA2JET). Each systematic uncertainty sources is fully uncorrelated with the other sources and fully correlated across bins, except for the background modelling systematics for which an uncorrelated treatment across bins is more appropriate.

Measured differential cross section with associated uncertainties as a function of DPHI(2GAMMA,2JET). Each systematic uncertainty sources is fully uncorrelated with the other sources and fully correlated across bins, except for the background modelling systematics for which an uncorrelated treatment across bins is more appropriate.

Measured differential cross section with associated uncertainties as a function of TAUJET. Each systematic uncertainty sources is fully uncorrelated with the other sources and fully correlated across bins, except for the background modelling systematics for which an uncorrelated treatment across bins is more appropriate.

Measured differential cross section with associated uncertainties as a function of SUM(TAUJET). Each systematic uncertainty sources is fully uncorrelated with the other sources and fully correlated across bins, except for the background modelling systematics for which an uncorrelated treatment across bins is more appropriate.

Measured differential cross section with associated uncertainties as a function of PT(2GAMMA) [NJET=0,PT>30 GEV]. Each systematic uncertainty sources is fully uncorrelated with the other sources and fully correlated across bins, except for the background modelling systematics for which an uncorrelated treatment across bins is more appropriate.

Measured differential cross section with associated uncertainties as a function of PT(2GAMMA) [NJET=1,PT>30 GEV]. Each systematic uncertainty sources is fully uncorrelated with the other sources and fully correlated across bins, except for the background modelling systematics for which an uncorrelated treatment across bins is more appropriate.

Measured differential cross section with associated uncertainties as a function of PT(2GAMMA) [NJET=2,PT>30 GEV]. Each systematic uncertainty sources is fully uncorrelated with the other sources and fully correlated across bins, except for the background modelling systematics for which an uncorrelated treatment across bins is more appropriate.

Measured differential cross section with associated uncertainties as a function of PT(2GAMMA) [NJET>=3,PT>30 GEV]. Each systematic uncertainty sources is fully uncorrelated with the other sources and fully correlated across bins, except for the background modelling systematics for which an uncorrelated treatment across bins is more appropriate.

The measured cross sections or cross section limits of the diphoton, VBF-enhanced, Nlepton $\geq$ 1, high $E_{T}^{miss}$, and ttH-enhanced fiducial regions are shown.

Measured differential cross section with associated uncertainties as a function of diphoton transverse momentum in bins of ABS(COS(THETA*)). Each systematic uncertainty sources is fully uncorrelated with the other sources and fully correlated across bins, except for the background modelling systematics for which an uncorrelated treatment across bins is more appropriate.Each systematic uncertainty sources is fully uncorrelated with the other sources.

Non-perturbative correction factors in percent accounting for the impact of hadronisation and the underlying event activity for all measured variables and fiducial regions. Regions of phase space where no reliable estimate could be obtained are listed as 100 without uncertainties. Uncertainties are evaluated by deriving these factors using different generators and tunes as described in the text. No factor are given for the Nlepton $\geq$ 1 and High-$E_{T}^{miss}$ fiducial regions as the gluon fusion contamination in both is negligible.

Isolation efficiencies in percent for gluon fusion $H\rightarrow\gamma\gamma$ for each fiducial region/variable bin measured in this analysis. The isolation efficiency is defined as the probability for both photons to fulfil the isolation criteria (as described in Section 9.1) for events that pass the diphoton kinematic criteria. Regions of phase space where no reliable estimate could be obtained are listed as 100 without uncertainties. Uncertainties are assigned in the same way as for the non-perturbative correction factors: by varying the fragmentation and underlying event modelling. These factors can be multiplied by the kinematic acceptance factors (see Table 29) to extrapolate an inclusive gluon fusion Higgs prediction to the fiducial volume used in this analysis. No factors for the Nlepton $\geq$ 1 and High $E_{T}^{miss}$ fiducial regions are provided as the gluon fusion contamination is negligible.

Combined non-perturbative (Table 27) and particle-level isolation correction factors (Table 28) in percent accounting for the impact of hadronisation and the underlying event activity for all measured variables and fiducial regions. Regions of phase space where no reliable estimate could be obtained are listed as 100 without uncertainties. The uncertainties on the combined values properly take into account the correlations between both multiplicative factors.

Diphoton kinematic acceptances in percent for gluon-gluon fusion for the diphoton fiducial region and all differential variable bins studied in this paper, defined as the probability to fulfill the diphoton kinematic criteria: $p_{T}$/$m_{\gamma\gamma}$ < 0.35 (0.25) for the leading (subleading) photon and $|\eta_{\gamma\gamma}|$ < 2.37. The factors are evaluated using the Powheg NNLOPSevent generator. Uncertainties are taken from PDF variations. QCD scale variations have a negligible impact on these factors. The range of each bin is given in Table 26.

observed statistical correlations between pTyy, Njets, mjj, |DeltaPhijj| and pTj1

ggH default MC + XH predictions

XH ( = VBF + VH + ttH + bbH ) MC predictions

Best-fit values and uncertainties of the production-mode cross sections times branching ratio.

Best-fit values and uncertainties of the simplified template cross sections times branching ratio.

Observed correlations between the measured simplified template cross sections, including both the statistical and systematic uncertainties.

Best-fit values and uncertainties of the simplified template cross sections times branching ratio.

Observed correlations between the measured simplified template cross sections, including both the statistical and systematic uncertainties.

Observed correlations between the measured simplified template cross sections, including both the statistical and systematic uncertainties.

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.

Observed and predicted mTtot distribution in the b-tag category of the 2tau_h channel. Despite listing this as an exclusive final state (as there must be at least one b-jets), there is no explicit selection on the presence of additional light-flavour jets. Please note that the bin content is divided by the bin width in the paper figure, but not in the HepData table. The last bin includes overflows. The combined prediction for A and H bosons with masses of 300, 500 and 800 GeV and $\tan\beta$ = 10 in the hMSSM scenario are also provided.

Observed and predicted mTtot distribution for the b-inclusive selection in the 1l1tau_h channel. Please note that the bin content is divided by the bin width in the paper figure, but not in the HepData table. In the paper, the first bin is cut off at 60 GeV for aesthetics but contains underflows down to 50 GeV as in the HepData table. The last bin includes overflows. The prediction for a SSM Zprime with masses of 1500, 2000 and 2500 GeV are also provided.

Observed and predicted mTtot distribution for the b-inclusive selection in the 2tau_h channel. Please note that the bin content is divided by the bin width in the paper figure, but not in the HepData table. The last bin includes overflows. The prediction for a SSM Zprime with masses of 1500, 2000 and 2500 GeV are also provided.

Observed and expected 95% CL upper limits on the b-associated Higgs boson production cross section times ditau branching fraction as a function of the boson mass.

Observed and expected 95% CL upper limits on the Drell Yan production cross section times ditau branching fraction as a function of the Zprime boson mass.

Observed and expected 95% CL upper limits on the Higgs boson production cross section times ditau branching fraction as a function of the boson mass and the relative strength of the b-associated production.

Ratio of the 95% CL upper limits on the production cross section times branching fraction for alternate Zprime models with respect to the SSM, both observed and expected are shown.

Acceptance, acceptance times efficiency and b-tag category fraction for a scalar boson produced by gluon-gluon fusion as a function of the scalar boson mass.

Acceptance, acceptance times efficiency and b-tag category fraction for a scalar boson produced by b-associated production as a function of the scalar boson mass.

Acceptance and acceptance times efficiency for a heavy gauge boson produced by Drell Yan as a function of the gauge boson mass.

Two dimensional likelihood scan of the gluon-gluon fusion cross section times braching fraction, $\sigma(gg\phi)\times B(\phi\to\tau\tau)$, vs the b-associated production times branching fraction, $\sigma(bb\phi)\times B(\phi\to\tau\tau)$ for the Higgs boson mass ($m_\phi$) indicated in the table. For each mass, 10000 points are scanned. At each point $\Delta(\mathrm{NLL})$ is calculated, defined as the negative-log-likelihood (NLL) of the conditional fit with $\sigma(gg\phi)$ and $\sigma(bb\phi)$ fixed to their values at the point and with the minimum NLL value at any point subtracted. Vaules are provided for the fit to the observed data and to the expected data, which is the sum of Standard Model contributions not including the SM Higgs boson. The best-fit point and the preferred 68% and 95% boundaries are found at $2\Delta(\mathrm{NLL})$ values of 0.0, 2.30 and 5.90, respectively.

Two dimensional likelihood scan of the gluon-gluon fusion cross section times braching fraction, $\sigma(gg\phi)\times B(\phi\to\tau\tau)$, vs the b-associated production times branching fraction, $\sigma(bb\phi)\times B(\phi\to\tau\tau)$ for the Higgs boson mass ($m_\phi$) indicated in the table. For each mass, 10000 points are scanned. At each point $\Delta(\mathrm{NLL})$ is calculated, defined as the negative-log-likelihood (NLL) of the conditional fit with $\sigma(gg\phi)$ and $\sigma(bb\phi)$ fixed to their values at the point and with the minimum NLL value at any point subtracted. Vaules are provided for the fit to the observed data and to the expected data, which is the sum of Standard Model contributions not including the SM Higgs boson. The best-fit point and the preferred 68% and 95% boundaries are found at $2\Delta(\mathrm{NLL})$ values of 0.0, 2.30 and 5.90, respectively.

Two dimensional likelihood scan of the gluon-gluon fusion cross section times braching fraction, $\sigma(gg\phi)\times B(\phi\to\tau\tau)$, vs the b-associated production times branching fraction, $\sigma(bb\phi)\times B(\phi\to\tau\tau)$ for the Higgs boson mass ($m_\phi$) indicated in the table. For each mass, 10000 points are scanned. At each point $\Delta(\mathrm{NLL})$ is calculated, defined as the negative-log-likelihood (NLL) of the conditional fit with $\sigma(gg\phi)$ and $\sigma(bb\phi)$ fixed to their values at the point and with the minimum NLL value at any point subtracted. Vaules are provided for the fit to the observed data and to the expected data, which is the sum of Standard Model contributions not including the SM Higgs boson. The best-fit point and the preferred 68% and 95% boundaries are found at $2\Delta(\mathrm{NLL})$ values of 0.0, 2.30 and 5.90, respectively.

Two dimensional likelihood scan of the gluon-gluon fusion cross section times braching fraction, $\sigma(gg\phi)\times B(\phi\to\tau\tau)$, vs the b-associated production times branching fraction, $\sigma(bb\phi)\times B(\phi\to\tau\tau)$ for the Higgs boson mass ($m_\phi$) indicated in the table. For each mass, 10000 points are scanned. At each point $\Delta(\mathrm{NLL})$ is calculated, defined as the negative-log-likelihood (NLL) of the conditional fit with $\sigma(gg\phi)$ and $\sigma(bb\phi)$ fixed to their values at the point and with the minimum NLL value at any point subtracted. Vaules are provided for the fit to the observed data and to the expected data, which is the sum of Standard Model contributions not including the SM Higgs boson. The best-fit point and the preferred 68% and 95% boundaries are found at $2\Delta(\mathrm{NLL})$ values of 0.0, 2.30 and 5.90, respectively.

Two dimensional likelihood scan of the gluon-gluon fusion cross section times braching fraction, $\sigma(gg\phi)\times B(\phi\to\tau\tau)$, vs the b-associated production times branching fraction, $\sigma(bb\phi)\times B(\phi\to\tau\tau)$ for the Higgs boson mass ($m_\phi$) indicated in the table. For each mass, 10000 points are scanned. At each point $\Delta(\mathrm{NLL})$ is calculated, defined as the negative-log-likelihood (NLL) of the conditional fit with $\sigma(gg\phi)$ and $\sigma(bb\phi)$ fixed to their values at the point and with the minimum NLL value at any point subtracted. Vaules are provided for the fit to the observed data and to the expected data, which is the sum of Standard Model contributions not including the SM Higgs boson. The best-fit point and the preferred 68% and 95% boundaries are found at $2\Delta(\mathrm{NLL})$ values of 0.0, 2.30 and 5.90, respectively.

Two dimensional likelihood scan of the gluon-gluon fusion cross section times braching fraction, $\sigma(gg\phi)\times B(\phi\to\tau\tau)$, vs the b-associated production times branching fraction, $\sigma(bb\phi)\times B(\phi\to\tau\tau)$ for the Higgs boson mass ($m_\phi$) indicated in the table. For each mass, 10000 points are scanned. At each point $\Delta(\mathrm{NLL})$ is calculated, defined as the negative-log-likelihood (NLL) of the conditional fit with $\sigma(gg\phi)$ and $\sigma(bb\phi)$ fixed to their values at the point and with the minimum NLL value at any point subtracted. Vaules are provided for the fit to the observed data and to the expected data, which is the sum of Standard Model contributions not including the SM Higgs boson. The best-fit point and the preferred 68% and 95% boundaries are found at $2\Delta(\mathrm{NLL})$ values of 0.0, 2.30 and 5.90, respectively.

Two dimensional likelihood scan of the gluon-gluon fusion cross section times braching fraction, $\sigma(gg\phi)\times B(\phi\to\tau\tau)$, vs the b-associated production times branching fraction, $\sigma(bb\phi)\times B(\phi\to\tau\tau)$ for the Higgs boson mass ($m_\phi$) indicated in the table. For each mass, 10000 points are scanned. At each point $\Delta(\mathrm{NLL})$ is calculated, defined as the negative-log-likelihood (NLL) of the conditional fit with $\sigma(gg\phi)$ and $\sigma(bb\phi)$ fixed to their values at the point and with the minimum NLL value at any point subtracted. Vaules are provided for the fit to the observed data and to the expected data, which is the sum of Standard Model contributions not including the SM Higgs boson. The best-fit point and the preferred 68% and 95% boundaries are found at $2\Delta(\mathrm{NLL})$ values of 0.0, 2.30 and 5.90, respectively.

Two dimensional likelihood scan of the gluon-gluon fusion cross section times braching fraction, $\sigma(gg\phi)\times B(\phi\to\tau\tau)$, vs the b-associated production times branching fraction, $\sigma(bb\phi)\times B(\phi\to\tau\tau)$ for the Higgs boson mass ($m_\phi$) indicated in the table. For each mass, 10000 points are scanned. At each point $\Delta(\mathrm{NLL})$ is calculated, defined as the negative-log-likelihood (NLL) of the conditional fit with $\sigma(gg\phi)$ and $\sigma(bb\phi)$ fixed to their values at the point and with the minimum NLL value at any point subtracted. Vaules are provided for the fit to the observed data and to the expected data, which is the sum of Standard Model contributions not including the SM Higgs boson. The best-fit point and the preferred 68% and 95% boundaries are found at $2\Delta(\mathrm{NLL})$ values of 0.0, 2.30 and 5.90, respectively.

Two dimensional likelihood scan of the gluon-gluon fusion cross section times braching fraction, $\sigma(gg\phi)\times B(\phi\to\tau\tau)$, vs the b-associated production times branching fraction, $\sigma(bb\phi)\times B(\phi\to\tau\tau)$ for the Higgs boson mass ($m_\phi$) indicated in the table. For each mass, 10000 points are scanned. At each point $\Delta(\mathrm{NLL})$ is calculated, defined as the negative-log-likelihood (NLL) of the conditional fit with $\sigma(gg\phi)$ and $\sigma(bb\phi)$ fixed to their values at the point and with the minimum NLL value at any point subtracted. Vaules are provided for the fit to the observed data and to the expected data, which is the sum of Standard Model contributions not including the SM Higgs boson. The best-fit point and the preferred 68% and 95% boundaries are found at $2\Delta(\mathrm{NLL})$ values of 0.0, 2.30 and 5.90, respectively.

Two dimensional likelihood scan of the gluon-gluon fusion cross section times braching fraction, $\sigma(gg\phi)\times B(\phi\to\tau\tau)$, vs the b-associated production times branching fraction, $\sigma(bb\phi)\times B(\phi\to\tau\tau)$ for the Higgs boson mass ($m_\phi$) indicated in the table. For each mass, 10000 points are scanned. At each point $\Delta(\mathrm{NLL})$ is calculated, defined as the negative-log-likelihood (NLL) of the conditional fit with $\sigma(gg\phi)$ and $\sigma(bb\phi)$ fixed to their values at the point and with the minimum NLL value at any point subtracted. Vaules are provided for the fit to the observed data and to the expected data, which is the sum of Standard Model contributions not including the SM Higgs boson. The best-fit point and the preferred 68% and 95% boundaries are found at $2\Delta(\mathrm{NLL})$ values of 0.0, 2.30 and 5.90, respectively.

Two dimensional likelihood scan of the gluon-gluon fusion cross section times braching fraction, $\sigma(gg\phi)\times B(\phi\to\tau\tau)$, vs the b-associated production times branching fraction, $\sigma(bb\phi)\times B(\phi\to\tau\tau)$ for the Higgs boson mass ($m_\phi$) indicated in the table. For each mass, 10000 points are scanned. At each point $\Delta(\mathrm{NLL})$ is calculated, defined as the negative-log-likelihood (NLL) of the conditional fit with $\sigma(gg\phi)$ and $\sigma(bb\phi)$ fixed to their values at the point and with the minimum NLL value at any point subtracted. Vaules are provided for the fit to the observed data and to the expected data, which is the sum of Standard Model contributions not including the SM Higgs boson. The best-fit point and the preferred 68% and 95% boundaries are found at $2\Delta(\mathrm{NLL})$ values of 0.0, 2.30 and 5.90, respectively.

Two dimensional likelihood scan of the gluon-gluon fusion cross section times braching fraction, $\sigma(gg\phi)\times B(\phi\to\tau\tau)$, vs the b-associated production times branching fraction, $\sigma(bb\phi)\times B(\phi\to\tau\tau)$ for the Higgs boson mass ($m_\phi$) indicated in the table. For each mass, 10000 points are scanned. At each point $\Delta(\mathrm{NLL})$ is calculated, defined as the negative-log-likelihood (NLL) of the conditional fit with $\sigma(gg\phi)$ and $\sigma(bb\phi)$ fixed to their values at the point and with the minimum NLL value at any point subtracted. Vaules are provided for the fit to the observed data and to the expected data, which is the sum of Standard Model contributions not including the SM Higgs boson. The best-fit point and the preferred 68% and 95% boundaries are found at $2\Delta(\mathrm{NLL})$ values of 0.0, 2.30 and 5.90, respectively.

Two dimensional likelihood scan of the gluon-gluon fusion cross section times braching fraction, $\sigma(gg\phi)\times B(\phi\to\tau\tau)$, vs the b-associated production times branching fraction, $\sigma(bb\phi)\times B(\phi\to\tau\tau)$ for the Higgs boson mass ($m_\phi$) indicated in the table. For each mass, 10000 points are scanned. At each point $\Delta(\mathrm{NLL})$ is calculated, defined as the negative-log-likelihood (NLL) of the conditional fit with $\sigma(gg\phi)$ and $\sigma(bb\phi)$ fixed to their values at the point and with the minimum NLL value at any point subtracted. Vaules are provided for the fit to the observed data and to the expected data, which is the sum of Standard Model contributions not including the SM Higgs boson. The best-fit point and the preferred 68% and 95% boundaries are found at $2\Delta(\mathrm{NLL})$ values of 0.0, 2.30 and 5.90, respectively.

Two dimensional likelihood scan of the gluon-gluon fusion cross section times braching fraction, $\sigma(gg\phi)\times B(\phi\to\tau\tau)$, vs the b-associated production times branching fraction, $\sigma(bb\phi)\times B(\phi\to\tau\tau)$ for the Higgs boson mass ($m_\phi$) indicated in the table. For each mass, 10000 points are scanned. At each point $\Delta(\mathrm{NLL})$ is calculated, defined as the negative-log-likelihood (NLL) of the conditional fit with $\sigma(gg\phi)$ and $\sigma(bb\phi)$ fixed to their values at the point and with the minimum NLL value at any point subtracted. Vaules are provided for the fit to the observed data and to the expected data, which is the sum of Standard Model contributions not including the SM Higgs boson. The best-fit point and the preferred 68% and 95% boundaries are found at $2\Delta(\mathrm{NLL})$ values of 0.0, 2.30 and 5.90, respectively.

Two dimensional likelihood scan of the gluon-gluon fusion cross section times braching fraction, $\sigma(gg\phi)\times B(\phi\to\tau\tau)$, vs the b-associated production times branching fraction, $\sigma(bb\phi)\times B(\phi\to\tau\tau)$ for the Higgs boson mass ($m_\phi$) indicated in the table. For each mass, 10000 points are scanned. At each point $\Delta(\mathrm{NLL})$ is calculated, defined as the negative-log-likelihood (NLL) of the conditional fit with $\sigma(gg\phi)$ and $\sigma(bb\phi)$ fixed to their values at the point and with the minimum NLL value at any point subtracted. Vaules are provided for the fit to the observed data and to the expected data, which is the sum of Standard Model contributions not including the SM Higgs boson. The best-fit point and the preferred 68% and 95% boundaries are found at $2\Delta(\mathrm{NLL})$ values of 0.0, 2.30 and 5.90, respectively.

Observed and expected 95% CL upper limits on the gluon-gluon fusion Higgs boson production cross section times ditau branching fraction as a function of the boson mass.

Transverse mass distribution for events satisfying all selection criteria in the muon channel.

Upper limits at the 95% CL on the cross section for SSM W' production and decay to the electron+neutrino channel as a function of the W' pole mass.

Upper limits at the 95% CL on the cross section for SSM W' production and decay to the electron+neutrino channel as a function of the W' pole mass.

Upper limits at the 95% CL on the cross section for SSM W' production and decay to the muon+neutrino channel as a function of the W' pole mass.

Upper limits at the 95% CL on the cross section for SSM W' production and decay to the muon+neutrino channel as a function of the W' pole mass.

Combined (electron and muon channels) upper limits at the 95% CL on the cross section for SSM W' production and decay to a single lepton generation as a function of the W' pole mass.

Combined (electron and muon channels) upper limits at the 95% CL on the cross section for SSM W' production and decay to a single lepton generation as a function of the W' pole mass.

Product of acceptance and efficiency for the electron and muon selections as a function of the W' pole mass.

Product of acceptance and efficiency for the electron and muon selections as a function of the W' pole mass.

Dijet mass spectra after the background only fit with the background prediction in the high-mass region with two b-tags.

Dijet mass spectra after the background only fit with the background prediction in the low-mass region with two b-tags.

The online b-tagging efficiency with respect to the offline b-tagging efficiency as a function of pT. The b-tagging online and offline working points correspond to an efficiency of 60% and 70%, respectively.

Observed and expected 95% credibility-level upper limits on the cross-section for the b* model in the high-mass region with inclusive b-jet selection.

Observed and expected 95% credibility-level upper limits on the cross-section times branching ratio for the SSM and leptophobic Z' models in the low- and high-mass region with two b-tags selection.

Observed and expected 95% credibility-level upper limits on the cross-section for DM Z' models in the low-mass region with two b-tags selection. The Z' is expected to decay to all five quark flavors other than the top quark and the mediator to SM quark coupling (gSM) equal to 0.1 is assumed.

Observed and expected 95% credibility-level upper limits on the cross-section times branching ratio for DM Z'->bb models in the high-mass region with two b-tags selection. The Z' is expected to decay to bb only and the mediator to SM quark coupling (gSM) equal to 0.25 is assumed.

Observed and expected 95% credibility-level upper limits on cross section times acceptance times branching ratio of X --> bb, including kinematic acceptance and b-tagging efficiencies, for resonances with intrinsic width smaller than the detector resolution. The width of the Gaussian reconstructed shape is dominated by the dijet mass resolution. The table shows the limits obtained from the high-mass inclusive one b-tag selection.

Observed and expected 95% credibility-level upper limits on cross section times acceptance times branching ratio of X --> bb, including kinematic acceptance and b-tagging efficiencies, for resonances with intrinsic width smaller than the detector resolution. The width of the Gaussian reconstructed shape is dominated by the dijet mass resolution. The table shows the limits obtained from the combined low- and high-mass two b-tags selection.

The mass distributions for the inclusive one b-tagged selection and two b-tagged selection using an integrated luminosity of 36.1 $fb^{-1}$. The inclusive one b-tagged Pythia8 MC distribution is normalized to the inclusive one b-tagged data. The two b-tagged Pythia8 MC distribution is normalized to the two b-tagged data. The systematic uncertainty band is dominated by the b-tagging scale factor and the b-jet energy scale.

Signal acceptance times efficiency for inclusive 1 b-tag and 2 b-tag categories as a function of the simulated b* and the Z' masses.

Signal acceptance times efficiency for 2 b-tag categories as a function of the simulated Z' masses.

The flavor composition of the simulated dijet background as a function of dijet mass before tagging. The fraction of the six combinations of the b-jet , c-jet and light-flavor jet are shown. All offline selections are applied.

The flavor composition of the simulated dijet background as a function of dijet mass with inclusive one b-tag. The fraction of the six combinations of the b-jet , c-jet and light-flavor jet are shown. All offline selections are applied.

The flavor composition of the simulated dijet background as a function of dijet mass with two b-tags. The fraction of the six combinations of the b-jet , c-jet and light-flavor jet are shown. All offline selections are applied.

Observed and expected 95% credibility-level upper limits on cross section times acceptance times branching ratio of X --> bb, including kinematic acceptance and b-tagging efficiencies, for resonances exhibiting a generic Gaussian shape at particle level. The table shows the limits obtained from the inclusive b-jet selection. The limits corresponding to Gaussian-shaped resonances with width of Γ(X)/m(X) = 3%.

Observed and expected 95% credibility-level upper limits on cross section times acceptance times branching ratio of X --> bb, including kinematic acceptance and b-tagging efficiencies, for resonances exhibiting a generic Gaussian shape at particle level. The table shows the limits obtained from the inclusive b-jet selection. The limits corresponding to Gaussian-shaped resonances with width of Γ(X)/m(X) = 7%.

Observed and expected 95% credibility-level upper limits on cross section times acceptance times branching ratio of X --> bb, including kinematic acceptance and b-tagging efficiencies, for resonances exhibiting a generic Gaussian shape at particle level. The table shows the limits obtained from the inclusive b-jet selection. The limits corresponding to Gaussian-shaped resonances with width of Γ(X)/m(X) = 10%.

Observed and expected 95% credibility-level upper limits on cross section times acceptance times branching ratio of X --> bb, including kinematic acceptance and b-tagging efficiencies, for resonances exhibiting a generic Gaussian shape at particle level. The table shows the limits obtained from the inclusive b-jet selection. The limits corresponding to Gaussian-shaped resonances with width of Γ(X)/m(X) = 15%.

Observed and expected 95% credibility-level upper limits on cross section times acceptance times branching ratio of X --> bb, including kinematic acceptance and b-tagging efficiencies, for resonances exhibiting a generic Gaussian shape at particle level. The table shows the limits obtained from the combined low- and high-mass two b-tags selection. The limits corresponding to Gaussian-shaped resonances with width of Γ(X)/m(X) = 3%.

Observed and expected 95% credibility-level upper limits on cross section times acceptance times branching ratio of X --> bb, including kinematic acceptance and b-tagging efficiencies, for resonances exhibiting a generic Gaussian shape at particle level. The table shows the limits obtained from the combined low- and high-mass two b-tags selection. The limits corresponding to Gaussian-shaped resonances with width of Γ(X)/m(X) = 7%.

Observed and expected 95% credibility-level upper limits on cross section times acceptance times branching ratio of X --> bb, including kinematic acceptance and b-tagging efficiencies, for resonances exhibiting a generic Gaussian shape at particle level. The table shows the limits obtained from the combined low- and high-mass two b-tags selection. The limits corresponding to Gaussian-shaped resonances with width of Γ(X)/m(X) = 10%.

Observed and expected 95% credibility-level upper limits on cross section times acceptance times branching ratio of X --> bb, including kinematic acceptance and b-tagging efficiencies, for resonances exhibiting a generic Gaussian shape at particle level. The table shows the limits obtained from the combined low- and high-mass two b-tags selection. The limits corresponding to Gaussian-shaped resonances with width of Γ(X)/m(X) = 15%.

1$\tau$ MediumMass SR eff.

2$\tau$ Compressed SR eff.

2$\tau$ Compressed SR eff.

2$\tau$ HighMass SR eff.

2$\tau$ HighMass SR eff.

2$\tau$ multibin SR eff.

2$\tau$ multibin SR eff.

2$\tau$ GMSB SR eff.

2$\tau$ GMSB SR eff.

1$\tau$ Compressed SR eff.

1$\tau$ Compressed SR eff.

1$\tau$ MediumMass SR eff.

1$\tau$ MediumMass SR eff.

2$\tau$ Compressed SR eff.

2$\tau$ Compressed SR eff.

2$\tau$ HighMass SR eff.

2$\tau$ HighMass SR eff.

2$\tau$ multibin SR eff.

2$\tau$ multibin SR eff.

2$\tau$ GMSB SR eff.

2$\tau$ GMSB SR eff.

1$\tau$ Compressed SR acceptance.

1$\tau$ Compressed SR acceptance.

1$\tau$ MediumMass SR acceptance.

1$\tau$ MediumMass SR acceptance.

2$\tau$ Compressed SR acceptance.

2$\tau$ Compressed SR acceptance.

2$\tau$ HighMass SR acceptance.

2$\tau$ HighMass SR acceptance.

2$\tau$ multibin SR acceptance.

2$\tau$ multibin SR acceptance.

2$\tau$ GMSB SR acceptance.

2$\tau$ GMSB SR acceptance.

1$\tau$ Compressed SR acceptance.

1$\tau$ Compressed SR acceptance.

1$\tau$ MediumMass SR acceptance.

1$\tau$ MediumMass SR acceptance.

2$\tau$ Compressed SR acceptance.

2$\tau$ Compressed SR acceptance.

2$\tau$ HighMass SR acceptance.

2$\tau$ HighMass SR acceptance.

2$\tau$ multibin SR acceptance.

2$\tau$ multibin SR acceptance.

2$\tau$ GMSB SR acceptance.

2$\tau$ GMSB SR acceptance.

Cutflow table of the $1\tau$ compressed SR for the four signal benchmark scenarios of low, medium, and high mass-splitting in the simplified model as well as the GMSB model.

Cutflow table of the $1\tau$ compressed SR for the four signal benchmark scenarios of low, medium, and high mass-splitting in the simplified model as well as the GMSB model.

Cutflow table of the $1\tau$ medium-mass SR for the four signal benchmark scenarios of low, medium, and high mass-splitting in the simplified model as well as the GMSB model.

Cutflow table of the $1\tau$ medium-mass SR for the four signal benchmark scenarios of low, medium, and high mass-splitting in the simplified model as well as the GMSB model.

Cutflow table of the $2\tau$ compressed SR for the four signal benchmark scenarios of low, medium, and high mass-splitting in the simplified model as well as the GMSB model.

Cutflow table of the $2\tau$ compressed SR for the four signal benchmark scenarios of low, medium, and high mass-splitting in the simplified model as well as the GMSB model.

Cutflow table of the $2\tau$ high-mass SR for the four signal benchmark scenarios of low, medium, and high mass-splitting in the simplified model as well as the GMSB model.

Cutflow table of the $2\tau$ high-mass SR for the four signal benchmark scenarios of low, medium, and high mass-splitting in the simplified model as well as the GMSB model.

Cutflow table of the $2\tau$ multibin SR for the four signal benchmark scenarios of low, medium, and high mass-splitting in the simplified model as well as the GMSB model.

Cutflow table of the $2\tau$ multibin SR for the four signal benchmark scenarios of low, medium, and high mass-splitting in the simplified model as well as the GMSB model.

Cutflow table of the $2\tau$ GMSB SR for the four signal benchmark scenarios of low, medium, and high mass-splitting in the simplified model as well as the GMSB model.

Cutflow table of the $2\tau$ GMSB SR for the four signal benchmark scenarios of low, medium, and high mass-splitting in the simplified model as well as the GMSB model.

Best performing fit setups entering the final combination as a function of the LSP mass and the gluino mass. 'S' marks the simultaneous fit of the four simplified model single-bin SRs, 'M' denotes the simultaneous fit of the two $1\tau$ SRs and the $2\tau$ multibin SR.

Best performing fit setups entering the final combination as a function of the LSP mass and the gluino mass. 'S' marks the simultaneous fit of the four simplified model single-bin SRs, 'M' denotes the simultaneous fit of the two $1\tau$ SRs and the $2\tau$ multibin SR.

Observed exclusion contour at 95% CL as a function of tanBeta and the SUSY-breaking mass scale Lambda.

Observed exclusion contour at 95% CL as a function of tanBeta and the SUSY-breaking mass scale Lambda.

Expected exclusion contour at 95% CL as a function of tanBeta and the SUSY-breaking mass scale Lambda.

Expected exclusion contour at 95% CL as a function of tanBeta and the SUSY-breaking mass scale Lambda.

Observed exclusion contour at 95% CL as a function of the LSP mass and the gluino mass.

Observed exclusion contour at 95% CL as a function of the LSP mass and the gluino mass.

Expected exclusion contour at 95% CL as a function of the LSP mass and the gluino mass.

Expected exclusion contour at 95% CL as a function of the LSP mass and the gluino mass.

Observed upper limits on the production cross section at 95% CL in pb as a function of tanBeta and SUSY breaking mass scale Lambda.

Observed upper limits on the production cross section at 95% CL in pb as a function of tanBeta and SUSY breaking mass scale Lambda.

Observed upper limits on the production cross section at 95% CL in pb as a function of the LSP mass and the gluino mass.

Observed upper limits on the production cross section at 95% CL in pb as a function of the LSP mass and the gluino mass.

Yields of the expected background from the SM in the bins of the multibin SR of the $2\tau$ channel with all bins being simultaneously used to constrain the background prediction. Expectation is given with the scalings computed in the combined fit applied. Uncertainties are statistial plus systematrics. Only the subsamples contributing the respective region are considered.

Yields of the expected background from the SM in the bins of the multibin SR of the $2\tau$ channel with all bins being simultaneously used to constrain the background prediction. Expectation is given with the scalings computed in the combined fit applied. Uncertainties are statistial plus systematrics. Only the subsamples contributing the respective region are considered.

$m_{\mathrm{T}}^{\tau}$ in the compressed $m_{\mathrm{T}}^{\tau}$ VR of the $1\tau$ channel, illustrating the background modeling after the fit. The last bin includes overflow events.

$m_{\mathrm{T}}^{\tau}$ in the compressed $m_{\mathrm{T}}^{\tau}$ VR of the $1\tau$ channel, illustrating the background modeling after the fit. The last bin includes overflow events.

$E_{\mathrm{T}}^{\mathrm{miss}}$ in the compressed $E_{\mathrm{T}}^{\mathrm{miss}}$ VR of the $1\tau$ channel, illustrating the background modeling after the fit. The last bin includes overflow events.

$E_{\mathrm{T}}^{\mathrm{miss}}$ in the compressed $E_{\mathrm{T}}^{\mathrm{miss}}$ VR of the $1\tau$ channel, illustrating the background modeling after the fit. The last bin includes overflow events.

$m_{\mathrm{T}}^{\tau}$ in the medium-mass $m_{\mathrm{T}}^{\tau}$ VR of the $1\tau$ channel, illustrating the background modeling after the fit. The last bin includes overflow events.

$m_{\mathrm{T}}^{\tau}$ in the medium-mass $m_{\mathrm{T}}^{\tau}$ VR of the $1\tau$ channel, illustrating the background modeling after the fit. The last bin includes overflow events.

$E_{\mathrm{T}}^{\mathrm{miss}}$ in the medium-mass $E_{\mathrm{T}}^{\mathrm{miss}}$ VR of the $1\tau$ channel, illustrating the background modeling after the fit. The last bin includes overflow events.

$E_{\mathrm{T}}^{\mathrm{miss}}$ in the medium-mass $E_{\mathrm{T}}^{\mathrm{miss}}$ VR of the $1\tau$ channel, illustrating the background modeling after the fit. The last bin includes overflow events.

$H_{\mathrm{T}}$ in the medium-mass $H_{\mathrm{T}}$ VR of the $1\tau$ channel, illustrating the background modeling after the fit. The last bin includes overflow events.

$H_{\mathrm{T}}$ in the medium-mass $H_{\mathrm{T}}$ VR of the $1\tau$ channel, illustrating the background modeling after the fit. The last bin includes overflow events.

$m_{\mathrm{T}}^{\tau_1}$ + $m_{\mathrm{T}}^{\tau_2}$ in the top VR of the $2\tau$ channel, illustrating the background modeling after the fit. The last bin includes overflow events.

$m_{\mathrm{T}}^{\tau_1}$ + $m_{\mathrm{T}}^{\tau_2}$ in the top VR of the $2\tau$ channel, illustrating the background modeling after the fit. The last bin includes overflow events.

$H_{\mathrm{T}}$ in the $W$ VR of the $2\tau$ channel, illustrating the background modeling after the fit. The last bin includes overflow events.

$H_{\mathrm{T}}$ in the $W$ VR of the $2\tau$ channel, illustrating the background modeling after the fit. The last bin includes overflow events.

$m_{\mathrm{T}}^{\tau_1}$ + $m_{\mathrm{T}}^{\tau_2}$ in the $Z$ VR of the $2\tau$ channel, illustrating the background modeling after the fit. The last bin includes overflow events.

$m_{\mathrm{T}}^{\tau_1}$ + $m_{\mathrm{T}}^{\tau_2}$ in the $Z$ VR of the $2\tau$ channel, illustrating the background modeling after the fit. The last bin includes overflow events.

$m_{\mathrm{T}}^{\tau}$ in the compressed SR of the $1\tau$ channel before application of the $m_{\mathrm{T}}^{\tau}$ > 80 GeV requirement. The last bin includes overflow events. Signal predictions corresponding to the simplified model scenarios of low (LM), medium (MM), and high mass-splitting (HM) as well as for the GMSB benchmark are given.

$m_{\mathrm{T}}^{\tau}$ in the compressed SR of the $1\tau$ channel before application of the $m_{\mathrm{T}}^{\tau}$ > 80 GeV requirement. The last bin includes overflow events. Signal predictions corresponding to the simplified model scenarios of low (LM), medium (MM), and high mass-splitting (HM) as well as for the GMSB benchmark are given.

$H_{\mathrm{T}}$ in the medium-mass SR of the $1\tau$ channel before application of the $H_{\mathrm{T}}$ > 1000 GeV requirement. The last bin includes overflow events. Signal predictions corresponding to the simplified model scenarios of low (LM), medium (MM), and high mass-splitting (HM) as well as for the GMSB benchmark are given.

$H_{\mathrm{T}}$ in the medium-mass SR of the $1\tau$ channel before application of the $H_{\mathrm{T}}$ > 1000 GeV requirement. The last bin includes overflow events. Signal predictions corresponding to the simplified model scenarios of low (LM), medium (MM), and high mass-splitting (HM) as well as for the GMSB benchmark are given.

$m_{\mathrm{T}}^{\mathrm{sum}}$ in the compressed SR of the $2\tau$ channel before application of the $m_{\mathrm{T}}^{\mathrm{sum}}$ > 1600 GeV requirement. The last bin includes overflow events. Signal predictions corresponding to the simplified model scenarios of low (LM), medium (MM), and high mass-splitting (HM) as well as for the GMSB benchmark are given.

$m_{\mathrm{T}}^{\mathrm{sum}}$ in the compressed SR of the $2\tau$ channel before application of the $m_{\mathrm{T}}^{\mathrm{sum}}$ > 1600 GeV requirement. The last bin includes overflow events. Signal predictions corresponding to the simplified model scenarios of low (LM), medium (MM), and high mass-splitting (HM) as well as for the GMSB benchmark are given.

$H_{\mathrm{T}}$ in the high-mass SR of the $2\tau$ channel before application of the $H_{\mathrm{T}}$ > 1100 GeV requirement. The last bin includes overflow events. Signal predictions corresponding to the simplified model scenarios of low (LM), medium (MM), and high mass-splitting (HM) as well as for the GMSB benchmark are given.

$H_{\mathrm{T}}$ in the high-mass SR of the $2\tau$ channel before application of the $H_{\mathrm{T}}$ > 1100 GeV requirement. The last bin includes overflow events. Signal predictions corresponding to the simplified model scenarios of low (LM), medium (MM), and high mass-splitting (HM) as well as for the GMSB benchmark are given.

mT(tau_1) + mT(tau_2) in the multibin SR of the 2T channel. The last bin includes overflow events. Signal predictions corresponding to the simplified model scenarios of low (LM), medium (MM), and high mass-splitting (HM) as well as for the GMSB benchmark are given.

mT(tau_1) + mT(tau_2) in the multibin SR of the 2T channel. The last bin includes overflow events. Signal predictions corresponding to the simplified model scenarios of low (LM), medium (MM), and high mass-splitting (HM) as well as for the GMSB benchmark are given.

$H_{\mathrm{T}}$ in the GMSB SR of the $2\tau$ channel before application of the $H_{\mathrm{T}}$ > 1900 GeV requirement. The last bin includes overflow events. Signal predictions corresponding to the simplified model scenarios of low (LM), medium (MM), and high mass-splitting (HM) as well as for the GMSB benchmark are given.

$H_{\mathrm{T}}$ in the GMSB SR of the $2\tau$ channel before application of the $H_{\mathrm{T}}$ > 1900 GeV requirement. The last bin includes overflow events. Signal predictions corresponding to the simplified model scenarios of low (LM), medium (MM), and high mass-splitting (HM) as well as for the GMSB benchmark are given.

Observed 95% confidence-level upper limits on the cross section for Drell-Yan spin-1/2 HECO production as a function of mass for various values of electric charge in the range $20\le|z|\le100$.

Selection efficiency as a function of transverse kinetic energy $E^\text{kin}_\text{T}=E_\text{kin}\sin\theta$ and pseudorapidity $|\eta|$ for $g=1g_\textrm{D}$ monopoles of mass 200 GeV.

Selection efficiency as a function of transverse kinetic energy $E^\text{kin}_\text{T}=E_\text{kin}\sin\theta$ and pseudorapidity $|\eta|$ for $g=1g_\textrm{D}$ monopoles of mass 500 GeV.

Selection efficiency as a function of transverse kinetic energy $E^\text{kin}_\text{T}=E_\text{kin}\sin\theta$ and pseudorapidity $|\eta|$ for $g=1g_\textrm{D}$ monopoles of mass 1000 GeV.

Selection efficiency as a function of transverse kinetic energy $E^\text{kin}_\text{T}=E_\text{kin}\sin\theta$ and pseudorapidity $|\eta|$ for $g=1g_\textrm{D}$ monopoles of mass 1500 GeV.

Selection efficiency as a function of transverse kinetic energy $E^\text{kin}_\text{T}=E_\text{kin}\sin\theta$ and pseudorapidity $|\eta|$ for $g=1g_\textrm{D}$ monopoles of mass 2000 GeV.

Selection efficiency as a function of transverse kinetic energy $E^\text{kin}_\text{T}=E_\text{kin}\sin\theta$ and pseudorapidity $|\eta|$ for $g=1g_\textrm{D}$ monopoles of mass 2500 GeV.

Selection efficiency as a function of transverse kinetic energy $E^\text{kin}_\text{T}=E_\text{kin}\sin\theta$ and pseudorapidity $|\eta|$ for $g=1g_\textrm{D}$ monopoles of mass 3000 GeV.

Selection efficiency as a function of transverse kinetic energy $E^\text{kin}_\text{T}=E_\text{kin}\sin\theta$ and pseudorapidity $|\eta|$ for $g=1g_\textrm{D}$ monopoles of mass 4000 GeV.

Selection efficiency as a function of transverse kinetic energy $E^\text{kin}_\text{T}=E_\text{kin}\sin\theta$ and pseudorapidity $|\eta|$ for $g=2g_\textrm{D}$ monopoles of mass 200 GeV.

Selection efficiency as a function of transverse kinetic energy $E^\text{kin}_\text{T}=E_\text{kin}\sin\theta$ and pseudorapidity $|\eta|$ for $g=2g_\textrm{D}$ monopoles of mass 500 GeV.

Selection efficiency as a function of transverse kinetic energy $E^\text{kin}_\text{T}=E_\text{kin}\sin\theta$ and pseudorapidity $|\eta|$ for $g=2g_\textrm{D}$ monopoles of mass 1000 GeV.

Selection efficiency as a function of transverse kinetic energy $E^\text{kin}_\text{T}=E_\text{kin}\sin\theta$ and pseudorapidity $|\eta|$ for $g=2g_\textrm{D}$ monopoles of mass 1500 GeV.

Selection efficiency as a function of transverse kinetic energy $E^\text{kin}_\text{T}=E_\text{kin}\sin\theta$ and pseudorapidity $|\eta|$ for $g=2g_\textrm{D}$ monopoles of mass 2000 GeV.

Selection efficiency as a function of transverse kinetic energy $E^\text{kin}_\text{T}=E_\text{kin}\sin\theta$ and pseudorapidity $|\eta|$ for $g=2g_\textrm{D}$ monopoles of mass 2500 GeV.

Selection efficiency as a function of transverse kinetic energy $E^\text{kin}_\text{T}=E_\text{kin}\sin\theta$ and pseudorapidity $|\eta|$ for $g=2g_\textrm{D}$ monopoles of mass 3000 GeV.

Selection efficiency as a function of transverse kinetic energy $E^\text{kin}_\text{T}=E_\text{kin}\sin\theta$ and pseudorapidity $|\eta|$ for $g=2g_\textrm{D}$ monopoles of mass 4000 GeV.

Total selection efficiency (i.e., the fraction of MC HECOs surviving the trigger and offline selection criteria) as a function of transverse kinetic energy $E^\text{kin}_\text{T}=E_\text{kin}\sin\theta$ and pseudorapidity $|\eta|$ for HECOs of charge $|z|=20$ of mass 200 GeV.

Total selection efficiency (i.e., the fraction of MC HECOs surviving the trigger and offline selection criteria) as a function of transverse kinetic energy $E^\text{kin}_\text{T}=E_\text{kin}\sin\theta$ and pseudorapidity $|\eta|$ for HECOs of charge $|z|=20$ of mass 500 GeV.

Total selection efficiency (i.e., the fraction of MC HECOs surviving the trigger and offline selection criteria) as a function of transverse kinetic energy $E^\text{kin}_\text{T}=E_\text{kin}\sin\theta$ and pseudorapidity $|\eta|$ for HECOs of charge $|z|=20$ of mass 1000 GeV.

Total selection efficiency (i.e., the fraction of MC HECOs surviving the trigger and offline selection criteria) as a function of transverse kinetic energy $E^\text{kin}_\text{T}=E_\text{kin}\sin\theta$ and pseudorapidity $|\eta|$ for HECOs of charge $|z|=20$ of mass 1500 GeV.

Total selection efficiency (i.e., the fraction of MC HECOs surviving the trigger and offline selection criteria) as a function of transverse kinetic energy $E^\text{kin}_\text{T}=E_\text{kin}\sin\theta$ and pseudorapidity $|\eta|$ for HECOs of charge $|z|=20$ of mass 2000 GeV.

Total selection efficiency (i.e., the fraction of MC HECOs surviving the trigger and offline selection criteria) as a function of transverse kinetic energy $E^\text{kin}_\text{T}=E_\text{kin}\sin\theta$ and pseudorapidity $|\eta|$ for HECOs of charge $|z|=20$ of mass 2500 GeV.

Total selection efficiency (i.e., the fraction of MC HECOs surviving the trigger and offline selection criteria) as a function of transverse kinetic energy $E^\text{kin}_\text{T}=E_\text{kin}\sin\theta$ and pseudorapidity $|\eta|$ for HECOs of charge $|z|=20$ of mass 3000 GeV.

Total selection efficiency (i.e., the fraction of MC HECOs surviving the trigger and offline selection criteria) as a function of transverse kinetic energy $E^\text{kin}_\text{T}=E_\text{kin}\sin\theta$ and pseudorapidity $|\eta|$ for HECOs of charge $|z|=20$ of mass 4000 GeV.

Total selection efficiency (i.e., the fraction of MC HECOs surviving the trigger and offline selection criteria) as a function of transverse kinetic energy $E^\text{kin}_\text{T}=E_\text{kin}\sin\theta$ and pseudorapidity $|\eta|$ for HECOs of charge $|z|=40$ of mass 200 GeV.

Total selection efficiency (i.e., the fraction of MC HECOs surviving the trigger and offline selection criteria) as a function of transverse kinetic energy $E^\text{kin}_\text{T}=E_\text{kin}\sin\theta$ and pseudorapidity $|\eta|$ for HECOs of charge $|z|=40$ of mass 500 GeV.

Total selection efficiency (i.e., the fraction of MC HECOs surviving the trigger and offline selection criteria) as a function of transverse kinetic energy $E^\text{kin}_\text{T}=E_\text{kin}\sin\theta$ and pseudorapidity $|\eta|$ for HECOs of charge $|z|=40$ of mass 1000 GeV.

Total selection efficiency (i.e., the fraction of MC HECOs surviving the trigger and offline selection criteria) as a function of transverse kinetic energy $E^\text{kin}_\text{T}=E_\text{kin}\sin\theta$ and pseudorapidity $|\eta|$ for HECOs of charge $|z|=40$ of mass 1500 GeV.

Total selection efficiency (i.e., the fraction of MC HECOs surviving the trigger and offline selection criteria) as a function of transverse kinetic energy $E^\text{kin}_\text{T}=E_\text{kin}\sin\theta$ and pseudorapidity $|\eta|$ for HECOs of charge $|z|=40$ of mass 2000 GeV.

Total selection efficiency (i.e., the fraction of MC HECOs surviving the trigger and offline selection criteria) as a function of transverse kinetic energy $E^\text{kin}_\text{T}=E_\text{kin}\sin\theta$ and pseudorapidity $|\eta|$ for HECOs of charge $|z|=40$ of mass 2500 GeV.

Total selection efficiency (i.e., the fraction of MC HECOs surviving the trigger and offline selection criteria) as a function of transverse kinetic energy $E^\text{kin}_\text{T}=E_\text{kin}\sin\theta$ and pseudorapidity $|\eta|$ for HECOs of charge $|z|=40$ of mass 3000 GeV.

Total selection efficiency (i.e., the fraction of MC HECOs surviving the trigger and offline selection criteria) as a function of transverse kinetic energy $E^\text{kin}_\text{T}=E_\text{kin}\sin\theta$ and pseudorapidity $|\eta|$ for HECOs of charge $|z|=40$ of mass 4000 GeV.

Total selection efficiency (i.e., the fraction of MC HECOs surviving the trigger and offline selection criteria) as a function of transverse kinetic energy $E^\text{kin}_\text{T}=E_\text{kin}\sin\theta$ and pseudorapidity $|\eta|$ for HECOs of charge $|z|=60$ of mass 200 GeV.

Total selection efficiency (i.e., the fraction of MC HECOs surviving the trigger and offline selection criteria) as a function of transverse kinetic energy $E^\text{kin}_\text{T}=E_\text{kin}\sin\theta$ and pseudorapidity $|\eta|$ for HECOs of charge $|z|=60$ of mass 500 GeV.

Total selection efficiency (i.e., the fraction of MC HECOs surviving the trigger and offline selection criteria) as a function of transverse kinetic energy $E^\text{kin}_\text{T}=E_\text{kin}\sin\theta$ and pseudorapidity $|\eta|$ for HECOs of charge $|z|=60$ of mass 1000 GeV.

Total selection efficiency (i.e., the fraction of MC HECOs surviving the trigger and offline selection criteria) as a function of transverse kinetic energy $E^\text{kin}_\text{T}=E_\text{kin}\sin\theta$ and pseudorapidity $|\eta|$ for HECOs of charge $|z|=60$ of mass 1500 GeV.

Total selection efficiency (i.e., the fraction of MC HECOs surviving the trigger and offline selection criteria) as a function of transverse kinetic energy $E^\text{kin}_\text{T}=E_\text{kin}\sin\theta$ and pseudorapidity $|\eta|$ for HECOs of charge $|z|=60$ of mass 2000 GeV.

Total selection efficiency (i.e., the fraction of MC HECOs surviving the trigger and offline selection criteria) as a function of transverse kinetic energy $E^\text{kin}_\text{T}=E_\text{kin}\sin\theta$ and pseudorapidity $|\eta|$ for HECOs of charge $|z|=60$ of mass 2500 GeV.

Total selection efficiency (i.e., the fraction of MC HECOs surviving the trigger and offline selection criteria) as a function of transverse kinetic energy $E^\text{kin}_\text{T}=E_\text{kin}\sin\theta$ and pseudorapidity $|\eta|$ for HECOs of charge $|z|=60$ of mass 3000 GeV.

Total selection efficiency (i.e., the fraction of MC HECOs surviving the trigger and offline selection criteria) as a function of transverse kinetic energy $E^\text{kin}_\text{T}=E_\text{kin}\sin\theta$ and pseudorapidity $|\eta|$ for HECOs of charge $|z|=60$ of mass 4000 GeV.

Total selection efficiency (i.e., the fraction of MC HECOs surviving the trigger and offline selection criteria) as a function of transverse kinetic energy $E^\text{kin}_\text{T}=E_\text{kin}\sin\theta$ and pseudorapidity $|\eta|$ for HECOs of charge $|z|=80$ of mass 200 GeV.

Total selection efficiency (i.e., the fraction of MC HECOs surviving the trigger and offline selection criteria) as a function of transverse kinetic energy $E^\text{kin}_\text{T}=E_\text{kin}\sin\theta$ and pseudorapidity $|\eta|$ for HECOs of charge $|z|=80$ of mass 500 GeV.

Total selection efficiency (i.e., the fraction of MC HECOs surviving the trigger and offline selection criteria) as a function of transverse kinetic energy $E^\text{kin}_\text{T}=E_\text{kin}\sin\theta$ and pseudorapidity $|\eta|$ for HECOs of charge $|z|=80$ of mass 1000 GeV.

Total selection efficiency (i.e., the fraction of MC HECOs surviving the trigger and offline selection criteria) as a function of transverse kinetic energy $E^\text{kin}_\text{T}=E_\text{kin}\sin\theta$ and pseudorapidity $|\eta|$ for HECOs of charge $|z|=80$ of mass 1500 GeV.

Total selection efficiency (i.e., the fraction of MC HECOs surviving the trigger and offline selection criteria) as a function of transverse kinetic energy $E^\text{kin}_\text{T}=E_\text{kin}\sin\theta$ and pseudorapidity $|\eta|$ for HECOs of charge $|z|=80$ of mass 2000 GeV.

Total selection efficiency (i.e., the fraction of MC HECOs surviving the trigger and offline selection criteria) as a function of transverse kinetic energy $E^\text{kin}_\text{T}=E_\text{kin}\sin\theta$ and pseudorapidity $|\eta|$ for HECOs of charge $|z|=80$ of mass 2500 GeV.

Total selection efficiency (i.e., the fraction of MC HECOs surviving the trigger and offline selection criteria) as a function of transverse kinetic energy $E^\text{kin}_\text{T}=E_\text{kin}\sin\theta$ and pseudorapidity $|\eta|$ for HECOs of charge $|z|=80$ of mass 3000 GeV.

Total selection efficiency (i.e., the fraction of MC HECOs surviving the trigger and offline selection criteria) as a function of transverse kinetic energy $E^\text{kin}_\text{T}=E_\text{kin}\sin\theta$ and pseudorapidity $|\eta|$ for HECOs of charge $|z|=80$ of mass 4000 GeV.

Total selection efficiency (i.e., the fraction of MC HECOs surviving the trigger and offline selection criteria) as a function of transverse kinetic energy $E^\text{kin}_\text{T}=E_\text{kin}\sin\theta$ and pseudorapidity $|\eta|$ for HECOs of charge $|z|=100$ of mass 200 GeV.

Total selection efficiency (i.e., the fraction of MC HECOs surviving the trigger and offline selection criteria) as a function of transverse kinetic energy $E^\text{kin}_\text{T}=E_\text{kin}\sin\theta$ and pseudorapidity $|\eta|$ for HECOs of charge $|z|=100$ of mass 500 GeV.

Total selection efficiency (i.e., the fraction of MC HECOs surviving the trigger and offline selection criteria) as a function of transverse kinetic energy $E^\text{kin}_\text{T}=E_\text{kin}\sin\theta$ and pseudorapidity $|\eta|$ for HECOs of charge $|z|=100$ of mass 1000 GeV.

Total selection efficiency (i.e., the fraction of MC HECOs surviving the trigger and offline selection criteria) as a function of transverse kinetic energy $E^\text{kin}_\text{T}=E_\text{kin}\sin\theta$ and pseudorapidity $|\eta|$ for HECOs of charge $|z|=100$ of mass 1500 GeV.

Total selection efficiency (i.e., the fraction of MC HECOs surviving the trigger and offline selection criteria) as a function of transverse kinetic energy $E^\text{kin}_\text{T}=E_\text{kin}\sin\theta$ and pseudorapidity $|\eta|$ for HECOs of charge $|z|=100$ of mass 2000 GeV.

Total selection efficiency (i.e., the fraction of MC HECOs surviving the trigger and offline selection criteria) as a function of transverse kinetic energy $E^\text{kin}_\text{T}=E_\text{kin}\sin\theta$ and pseudorapidity $|\eta|$ for HECOs of charge $|z|=100$ of mass 2500 GeV.

Total selection efficiency (i.e., the fraction of MC HECOs surviving the trigger and offline selection criteria) as a function of transverse kinetic energy $E^\text{kin}_\text{T}=E_\text{kin}\sin\theta$ and pseudorapidity $|\eta|$ for HECOs of charge $|z|=100$ of mass 3000 GeV.

Total selection efficiency (i.e., the fraction of MC HECOs surviving the trigger and offline selection criteria) as a function of transverse kinetic energy $E^\text{kin}_\text{T}=E_\text{kin}\sin\theta$ and pseudorapidity $|\eta|$ for HECOs of charge $|z|=100$ of mass 4000 GeV.

Upper limits at the 95% CL on the cross section for SSM $W^\prime$ production and decay to the muon+neutrino channel as a function of the $W^\prime$ pole mass.

Combined (electron and muon channels) upper limits at the 95% CL on the cross section for SSM $W^\prime$ production and decay to a single lepton generation as a function of the $W^\prime$ pole mass.

Observed upper limits at the 95% CL on the cross section for generic $W^\prime$ production and decay to the electron+neutrino channel as a function of the $W^\prime$ pole mass.

Observed upper limits at the 95% CL on the cross section for generic $W^\prime$ production and decay to the muon+neutrino channel as a function of the $W^\prime$ pole mass.

Combined (electron and muon channels) observed upper limits at the 95% CL on the cross section for generic $W^\prime$ production and decay to a single lepton generation as a function of the $W^\prime$ pole mass.

Observed upper limits at the 95% CL on the visible cross section in the electron+neutrino channel as a function of the transverse mass threshold.

Observed upper limits at the 95% CL on the visible cross section in the muon+neutrino channel as a function of the transverse mass threshold.

Product of acceptance and efficiency for the electron and muon selections as a function of the $W^\prime$ pole mass.

Expected upper limits at the 95% CL on the cross section for generic $W^\prime$ production and decay to the electron+neutrino channel as a function of the $W^\prime$ pole mass.

Expected upper limits at the 95% CL on the cross section for generic $W^\prime$ production and decay to the muon+neutrino channel as a function of the $W^\prime$ pole mass.

Combined (electron and muon channels) expected upper limits at the 95% CL on the cross section for generic $W^\prime$ production and decay to a single lepton generation as a function of the $W^\prime$ pole mass.

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.

Dijet invariant mass distribution for the category with at least one b-tagged jet.

Dijet invariant mass distribution for the category with both jets b-tagged.

The 95% CL upper limits on the cross-section times acceptance times branching ratio into two jets as a function of the mass of q* signal.

The 95% CL upper limits on the cross-section times acceptance times branching ratio into two jets as a function of the mass of QBH signal.

The 95% CL upper limits on the cross-section times acceptance times branching ratio into two jets as a function of the mass of W' signal.

The 95% CL upper limits on the cross-section times acceptance times branching ratio into two jets as a function of the mass of W* signal.

The upper limits on the DM mediator Z' signal at 95% CL from the inclusive category. The 95% CL upper limits are set on the universal quark coupling $g_q$ as a function of the Z' mass.

The 95% CL upper limit on the cross-section times acceptance times b-tagging efficiency times branching ratio as a function of the mass of the b* signal.

The 95% CL upper limit on the cross-section times acceptance times b-tagging efficiency times branching ratio as a function of the signal mass in the DM mediator Z' with $g_q$ = 0.25 model.

The 95% CL upper limit on the cross-section times acceptance times b-tagging efficiency times branching ratio as a function of the signal mass in the SSM Z' model.

The 95% CL upper limit on the cross-section times acceptance times b-tagging efficiency times branching ratio as a function of the signal mass in the graviton with $k/\overline{M}_{PL}=0.2$ model.

The 95% CL upper limit on the cross-section times kinematic acceptance times branching ratio for resonances with a generic Gaussian shape, as a function of the Gaussian mean mass in the inclusive category. Different widths, from 0% up to 15% of the signal mass, are considered. Gaussian-shape signals with 0% widths correspond to signal widths smaller than the experimental resolution. For a Gaussian-shaped signal with a relative width of 15%, the limits are truncated at high mass when the broad signal starts to overlap the upper end of the mass spectrum.

The 95% CL upper limit on the cross-section times kinematic acceptance times b-tagging efficiency times branching ratio for resonances with a generic Gaussian shape, as a function of the Gaussian mean mass in the 1 b category. Different widths, from 0% up to 15% of the signal mass, are considered. Gaussian-shape signals with 0% widths correspond to signal widths smaller than the experimental resolution. For a Gaussian-shaped signal with a relative width of 15%, the limits are truncated at high mass when the broad signal starts to overlap the upper end of the mass spectrum.

The 95% CL upper limit on the cross-section times kinematic acceptance times b-tagging efficiency times branching ratio for resonances with a generic Gaussian shape, as a function of the Gaussian mean mass in the 2 b category. Different widths, from 0% up to 15% of the signal mass, are considered. Gaussian-shape signals with 0% widths correspond to signal widths smaller than the experimental resolution. For a Gaussian-shaped signal with a relative width of 15%, the limits are truncated at high mass when the broad signal starts to overlap the upper end of the mass spectrum.

The expected 95% CL upper limits on the cross-section times acceptance times b-tagging efficiency times branching ratio as a function of the DM mediator Z' mass for the current and previous iterations of the analysis. The upper limit of the previous result was obtained with the Bayesian method and is also shown scaled to the 139 fb$^{-1}$ integrated luminosity of the current result to illustrate the effect of the analysis improvements.

Acceptance for the QBH, Q* and W' benchmark signal models in the inclusive category, as a function of the simulated mass.

Acceptance for the DM Z' benchmark signal model for various $g_q$ coupling parameters in the inclusive category, as a function of the simulated mass.

Acceptance for the W* benchmark signal model in the inclusive category, as a function of the simulated mass.

Acceptance times b-tagging efficiency for the b* benchmark signal model in the 1b category, as a function of the simulated mass.

Acceptance times b-tagging efficiency for the DM Z' with $g_q$=0.25, SSM Z' and graviton with $k/\overline{M}_{PL}=0.2$ benchmark signal models in the 2b category, as a function of the simulated mass.

The expected 95% CL upper limits on the cross-section times branching ratio as a function of the DM mediator Z' mass for the current and previous iterations of the analysis. The upper limit of the previous result was obtained with the Bayesian method and is also shown scaled to the 139 fb$^{-1}$ integrated luminosity of the current result to illustrate the effect of the analysis improvements. The current b-tagging requirement is tighter than the previous one for high-$p_T$ jets, resulting in a data sample with limited size for mass above 4 TeV. The background rejection, instead, has improved significantly across the entire mass spectrum inspected by the analysis.

Dependence of the overall signal efficiency on the pT of the long-lived Z' for Z' -> ee. The error bars indicate the total uncertainties.

Dependence of the overall signal efficiency on the transverse decay radius Rxy of the long-lived Z' for Z' -> emu. The error bars indicate the total uncertainties.

Dependence of the overall signal efficiency on the pT of the long-lived Z' for Z' -> emu. The error bars indicate the total uncertainties.

Dependence of the overall signal efficiency on the transverse decay radius Rxy of the long-lived Z' for Z' -> mumu. The error bars indicate the total uncertainties.

Dependence of the overall signal efficiency on the pT of the long-lived Z' for Z' -> mumu. The error bars indicate the total uncertainties.

Overall signal efficiency as a function of the mean proper lifetime (ctau) of the neutralino for the lambda121 scenario of the RPV SUSY model. The error bars indicate the total uncertainties.

Overall signal efficiency as a function of the mean proper lifetime (ctau) of the neutralino for the lambda122 scenario of the RPV SUSY model. The error bars indicate the total uncertainties.

95% CL upper limits on the squark-antisquark production cross-section as a function of the mean proper lifetime (ctau) of the neutralino for the lambda121 scenario of the RPV SUSY model and a 700 GeV squark. The uncertainties of the expected limit indicate the +-1sigma variations.

95% CL upper limits on the squark-antisquark production cross-section as a function of the mean proper lifetime (ctau) of the neutralino for the lambda122 scenario of the RPV SUSY model and a 700 GeV squark. The uncertainties of the expected limit indicate the +-1sigma variations.

95% CL upper limits on the squark-antisquark production cross-section as a function of the mean proper lifetime (ctau) of the neutralino for the lambda121 scenario of the RPV SUSY model and a 1600 GeV squark. The uncertainties of the expected limit indicate the +-1sigma variations.

95% CL upper limits on the squark-antisquark production cross-section as a function of the mean proper lifetime (ctau) of the neutralino for the lambda122 scenario of the RPV SUSY model and a 1600 GeV squark. The uncertainties of the expected limit indicate the +-1sigma variations.

Fraction of detector volume covered by the material veto as a function of z and Rxy of the displaced dilepton vertex.

Fraction of detector volume covered by the disabled pixel modules veto veto as a function of z and Rxy of the displaced dilepton vertex.

Observed Rxy distribution of vertices composed of two non-leptonic tracks in a control sample in the data and the predicted distribution obtained from the event mixing. The error bars indicate the statistical uncertainties.

Rxy distributions of Kshort vertices from the large radius tracking in the data and background MC samples. Data has been normalised such that the total number of Kshort from the standard tracking in the data agrees with the total number of Kshort from the standard tracking in the MC. The error bars indicate the statistical uncertainties.

Acceptance per decay as a function of the mean proper lifetime (ctau) of the neutralino for neutralino -> eenu. The error bars indicate the statistical uncertainties.

Acceptance per decay as a function of the mean proper lifetime (ctau) of the neutralino for neutralino -> emunu. The error bars indicate the statistical uncertainties.

Acceptance per decay as a function of the mean proper lifetime (ctau) of the neutralino for neutralino -> mumunu. The error bars indicate the statistical uncertainties.

Detection efficiency per decay as a function of the mean proper lifetime (ctau) of the neutralino for neutralino -> eenu. The error bars indicate the total uncertainties.

Detection efficiency per decay as a function of the mean proper lifetime (ctau) of the neutralino for neutralino -> emunu. The error bars indicate the total uncertainties.

Detection efficiency per decay as a function of the mean proper lifetime (ctau) of the neutralino for neutralino -> mumunu. The error bars indicate the total uncertainties.

Overall signal efficiency (acceptance times efficiency) per decay as a function of the mean proper lifetime (ctau) of the neutralino for neutralino -> eenu. The error bars indicate the total uncertainties.

Overall signal efficiency (acceptance times efficiency) per decay as a function of the mean proper lifetime (ctau) of the neutralino for neutralino -> emunu. The error bars indicate the total uncertainties.

Overall signal efficiency (acceptance times efficiency) per decay as a function of the mean proper lifetime (ctau) of the neutralino for neutralino -> mumunu. The error bars indicate the total uncertainties.

Detection efficiency per decay for Rxy < 22 mm as a function of the invariant mass and pT of the electron pair in LLP -> eeX.

Detection efficiency per decay for 22 <= Rxy < 38 mm as a function of the invariant mass and pT of the electron pair in LLP -> eeX.

Detection efficiency per decay for 38 <= Rxy < 73 mm as a function of the invariant mass and pT of the electron pair in LLP -> eeX.

Detection efficiency per decay for 73 <= Rxy < 111 mm as a function of the invariant mass and pT of the electron pair in LLP -> eeX.

Detection efficiency per decay for 111 <= Rxy < 145 mm as a function of the invariant mass and pT of the electron pair in LLP -> eeX.

Detection efficiency per decay for 145 <= Rxy < 300 mm as a function of the invariant mass and pT of the electron pair in LLP -> eeX.

Detection efficiency per decay for Rxy < 22 mm as a function of the invariant mass and pT of the electron and muon pair in LLP -> emuX.

Detection efficiency per decay for 22 <= Rxy < 38 mm as a function of the invariant mass and pT of the electron and muon pair in LLP -> emuX.

Detection efficiency per decay for 38 <= Rxy < 73 mm as a function of the invariant mass and pT of the electron and muon pair in LLP -> emuX.

Detection efficiency per decay for 73 <= Rxy < 111 mm as a function of the invariant mass and pT of the electron and muon pair in LLP -> emuX.

Detection efficiency per decay for 111 <= Rxy < 145 mm as a function of the invariant mass and pT of the electron and muon pair in LLP -> emuX.

Detection efficiency per decay for 145 <= Rxy < 300 mm as a function of the invariant mass and pT of the electron and muon pair in LLP -> emuX.

Detection efficiency per decay for Rxy < 22 mm as a function of the invariant mass and pT of the muon pair in LLP -> mumuX.

Detection efficiency per decay for 22 <= Rxy < 38 mm as a function of the invariant mass and pT of the muon pair in LLP -> mumuX.

Detection efficiency per decay for 38 <= Rxy < 73 mm as a function of the invariant mass and pT of the muon pair in LLP -> mumuX.

Detection efficiency per decay for 73 <= Rxy < 111 mm as a function of the invariant mass and pT of the muon pair in LLP -> mumuX.

Detection efficiency per decay for 111 <= Rxy < 145 mm as a function of the invariant mass and pT of the muon pair in LLP -> mumuX.

Detection efficiency per decay for 145 <= Rxy < 300 mm as a function of the invariant mass and pT of the muon pair in LLP -> mumuX.

Detection efficiency per decay as a function of the transverse decay radius Rxy and the dilepton pT for a LLP mass of 100 GeV and LLP -> ee.

Detection efficiency per decay as a function of the transverse decay radius Rxy and the dilepton pT for a LLP mass of 100 GeV and LLP -> emu.

Detection efficiency per decay as a function of the transverse decay radius Rxy and the dilepton pT for a LLP mass of 100 GeV and LLP -> mumu.

Detection efficiency per decay as a function of the transverse decay radius Rxy and the dilepton pT for a LLP mass of 250 GeV and LLP -> ee.

Detection efficiency per decay as a function of the transverse decay radius Rxy and the dilepton pT for a LLP mass of 250 GeV and LLP -> emu.

Detection efficiency per decay as a function of the transverse decay radius Rxy and the dilepton pT for a LLP mass of 250 GeV and LLP -> mumu.

Detection efficiency per decay as a function of the transverse decay radius Rxy and the dilepton pT for a LLP mass of 500 GeV and LLP -> ee.

Detection efficiency per decay as a function of the transverse decay radius Rxy and the dilepton pT for a LLP mass of 500 GeV and LLP -> emu.

Detection efficiency per decay as a function of the transverse decay radius Rxy and the dilepton pT for a LLP mass of 500 GeV and LLP -> mumu.

Detection efficiency per decay as a function of the transverse decay radius Rxy and the dilepton pT for a LLP mass of 750 GeV and LLP -> ee.

Detection efficiency per decay as a function of the transverse decay radius Rxy and the dilepton pT for a LLP mass of 750 GeV and LLP -> emu.

Detection efficiency per decay as a function of the transverse decay radius Rxy and the dilepton pT for a LLP mass of 750 GeV and LLP -> mumu.

Detection efficiency per decay as a function of the transverse decay radius Rxy and the dilepton pT for a LLP mass of 1000 GeV and LLP -> ee.

Detection efficiency per decay as a function of the transverse decay radius Rxy and the dilepton pT for a LLP mass of 1000 GeV and LLP -> emu.

Detection efficiency per decay as a function of the transverse decay radius Rxy and the dilepton pT for a LLP mass of 1000 GeV and LLP -> mumu.

Acceptance, detection efficiency, and overall signal efficiency in the Z' toy model for mZ' = 100 GeV, three mean proper lifetimes (ctau) and the three decay modes of the Z'.

Acceptance, detection efficiency, and overall signal efficiency in the Z' toy model for mZ' = 250 GeV, three mean proper lifetimes (ctau) and the three decay modes of the Z'.

Acceptance, detection efficiency, and overall signal efficiency in the Z' toy model for mZ' = 500 GeV, three mean proper lifetimes (ctau) and the three decay modes of the Z'.

Acceptance, detection efficiency, and overall signal efficiency in the Z' toy model for mZ' = 750 GeV, three mean proper lifetimes (ctau) and the three decay modes of the Z'.

Acceptance, detection efficiency, and overall signal efficiency in the Z' toy model for mZ' = 1000 GeV, three mean proper lifetimes (ctau) and the three decay modes of the Z'.

Post-fit yields for data and prediction in each of the single-bin signal regions of the analysis.

Observed 95% confidence level limit for the two-step signal grid.

Observed 95% confidence level limit for the two-step signal grid with the signal cross section increased by one sigma.

Observed 95% confidence level limit for the two-step signal grid with the signal cross section decreased by one sigma.

Expected 95% confidence level limit for the two-step signal grid.

Expected 95% confidence level limit for the two-step signal grid plus one sigma from experimental systematics.

Expected 95% confidence level limit for the two-step signal grid minus one sigma from experimental systematics.

Observed 95% confidence level limit for the Gtt signal grid.

Observed 95% confidence level limit for the Gtt signal grid with the signal cross section increased by one sigma.

Observed 95% confidence level limit for the Gtt signal grid with the signal cross section decreased by one sigma.

Expected 95% confidence level limit for the Gtt signal grid.

Expected 95% confidence level limit for the Gtt signal grid plus one sigma from experimental systematics.

Expected 95% confidence level limit for the Gtt signal grid minus one sigma from experimental systematics.

Observed 95% confidence level limit for the RPV signal grid.

Observed 95% confidence level limit for the RPV signal grid with the signal cross section increased by one sigma.

Observed 95% confidence level limit for the RPV signal grid with the signal cross section decreased by one sigma.

Expected 95% confidence level limit for the RPV signal grid.

Expected 95% confidence level limit for the RPV signal grid plus one sigma from experimental systematics.

Expected 95% confidence level limit for the RPV signal grid minus one sigma from experimental systematics.

Observed 95% confidence level limit for the two-step signal grid.

Expected 95% confidence level limit for the two-step signal grid.

Observed 95% confidence level limit for the Gtt signal grid.

Expected 95% confidence level limit for the Gtt signal grid.

Observed 95% confidence level limit for the RPV signal grid.

Expected 95% confidence level limit for the RPV signal grid.

$\mathcal{S}(E_{\mathrm{T}}^{\mathrm{miss}})$ distribution in the signal region SR-10ij50-0ib-MJ340. Two benchmark signal models are shown along with the background yields. These models, each representing a single mass point, are labelled 'RPV' with $(m_{\tilde{g}}, m_{\tilde{t}}) = (1600, 600) \, \mathrm{GeV}$ and 'two-step' with $(m_{\tilde{g}}, m_{\tilde{\chi^{0}_{1}}}) = (1600, 100) \, \mathrm{GeV}$.

$\mathcal{S}(E_{\mathrm{T}}^{\mathrm{miss}})$ distribution in the signal region SR-12ij50-2ib. Two benchmark signal models are shown along with the background yields. These models, each representing a single mass point, are labelled 'RPV' with $(m_{\tilde{g}}, m_{\tilde{t}}) = (1600, 600) \, \mathrm{GeV}$ and 'two-step' with $(m_{\tilde{g}}, m_{\tilde{\chi^{0}_{1}}}) = (1600, 100) \, \mathrm{GeV}$.

$\mathcal{S}(E_{\mathrm{T}}^{\mathrm{miss}})$ distribution in the signal region SR-9ij80-0ib. Two benchmark signal models are shown along with the background yields. These models, each representing a single mass point, are labelled 'RPV' with $(m_{\tilde{g}}, m_{\tilde{t}}) = (1600, 600) \, \mathrm{GeV}$ and 'two-step' with $(m_{\tilde{g}}, m_{\tilde{\chi^{0}_{1}}}) = (1600, 100) \, \mathrm{GeV}$.

Number of signal events expected for $139 \, \mathrm{fb}^{-1} $ after different analysis selections in the signal region SR-8ij50-0ib-MJ500. This 'two-step' model requires that a strongly produced gluino decays into quarks, the W and Z bosons, and the lightest stable neutralino where $(m_{\tilde{g}}, m_{\tilde{\chi^{0}_{1}}}) = (1600, 100) \, \mathrm{GeV}$.

Number of signal events expected for $139 \, \mathrm{fb}^{-1} $ after different analysis selections in the signal region SR-9ij50-0ib-MJ340. This 'two-step' model requires that a strongly produced gluino decays into quarks, the W and Z bosons, and the lightest stable neutralino where $(m_{\tilde{g}}, m_{\tilde{\chi^{0}_{1}}}) = (1600, 100) \, \mathrm{GeV}$.

Number of signal events expected for $139 \, \mathrm{fb}^{-1} $ after different analysis selections in the signal region SR-10ij50-0ib-MJ340. This 'two-step' model requires that a strongly produced gluino decays into quarks, the W and Z bosons, and the lightest stable neutralino where $(m_{\tilde{g}}, m_{\tilde{\chi^{0}_{1}}}) = (1600, 100) \, \mathrm{GeV}$.

Number of signal events expected for $139 \, \mathrm{fb}^{-1} $ after different analysis selections in the signal region SR-10ij50-0ib-MJ500. This 'two-step' model requires that a strongly produced gluino decays into quarks, the W and Z bosons, and the lightest stable neutralino where $(m_{\tilde{g}}, m_{\tilde{\chi^{0}_{1}}}) = (1600, 100) \, \mathrm{GeV}$.

Number of signal events expected for $139 \, \mathrm{fb}^{-1} $ after different analysis selections in the signal region SR-10ij50-1ib-MJ500. This 'two-step' model requires that a strongly produced gluino decays into quarks, the W and Z bosons, and the lightest stable neutralino where $(m_{\tilde{g}}, m_{\tilde{\chi^{0}_{1}}}) = (1600, 100) \, \mathrm{GeV}$.

Number of signal events expected for $139 \, \mathrm{fb}^{-1} $ after different analysis selections in the signal region SR-11ij50-0ib. This 'two-step' model requires that a strongly produced gluino decays into quarks, the W and Z bosons, and the lightest stable neutralino where $(m_{\tilde{g}}, m_{\tilde{\chi^{0}_{1}}}) = (1600, 100) \, \mathrm{GeV}$.

Number of signal events expected for $139 \, \mathrm{fb}^{-1} $ after different analysis selections in the signal region SR-12ij50-2ib. This 'two-step' model requires that a strongly produced gluino decays into quarks, the W and Z bosons, and the lightest stable neutralino where $(m_{\tilde{g}}, m_{\tilde{\chi^{0}_{1}}}) = (1600, 100) \, \mathrm{GeV}$.

Number of signal events expected for $139 \, \mathrm{fb}^{-1} $ after different analysis selections in the signal region SR-9ij80-0ib. This 'two-step' model requires that a strongly produced gluino decays into quarks, the W and Z bosons, and the lightest stable neutralino where $(m_{\tilde{g}}, m_{\tilde{\chi^{0}_{1}}}) = (1600, 100) \, \mathrm{GeV}$.

Acceptance for the signal region SR-8ij50-0ib-MJ500 showing the acceptance for the complete two-step signal grid.

Efficiency for the signal region SR-8ij50-0ib-MJ500 showing the efficiency for the complete two-step signal grid.

Acceptance for the signal region SR-9ij50-0ib-MJ340 showing the acceptance for the complete two-step signal grid.

Efficiency for the signal region SR-9ij50-0ib-MJ340 showing the efficiency for the complete two-step signal grid.

Acceptance for the signal region SR-10ij50-0ib-MJ340 showing the acceptance for the complete two-step signal grid.

Efficiency for the signal region SR-10ij50-0ib-MJ340 showing the efficiency for the complete two-step signal grid.

Acceptance for the signal region SR-10ij50-0ib-MJ500 showing the acceptance for the complete two-step signal grid.

Efficiency for the signal region SR-10ij50-0ib-MJ500 showing the efficiency for the complete two-step signal grid.

Acceptance for the signal region SR-10ij50-1ib-MJ500 showing the acceptance for the complete two-step signal grid.

Efficiency for the signal region SR-10ij50-1ib-MJ500 showing the efficiency for the complete two-step signal grid.

Acceptance for the signal region SR-11ij50-0ib showing the acceptance for the complete two-step signal grid.

Efficiency for the signal region SR-11ij50-0ib showing the efficiency for the complete two-step signal grid.

Acceptance for the signal region SR-12ij50-2ib showing the acceptance for the complete two-step signal grid.

Efficiency for the signal region SR-12ij50-2ib showing the efficiency for the complete two-step signal grid.

Acceptance for the signal region SR-9ij80-0ib showing the acceptance for the complete two-step signal grid.

Efficiency for the signal region SR-9ij80-0ib showing the efficiency for the complete two-step signal grid.

The normalisation factors for the dominant backgrounds of the analysis in each of the multi-bin and single-bin regions.

Post-fit yields for data and prediction in each of the single-bin validation regions to test the $N_{\mathrm{jet}}$ extraction.

Post-fit yields for data and prediction in each of the single-bin validation regions to test the $\mathcal{S}(E_{\mathrm{T}}^{\mathrm{miss}})$ extrapolation.