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