The expected 95% CL exclusion contours for the combined fit of SR-lowMass and SR-highMass for simplified models with combined ${\tilde{\tau}}^{+}_{R,L} {\tilde{\tau}}^{-}_{R,L}$ production.

The observed 95% CL exclusion contours for the combined fit of SR-lowMass and SR-highMass for simplified models with ${\tilde{\tau}}_L {\tilde{\tau}}_L$ only production.

The expected 95% CL exclusion contours for the combined fit of SR-lowMass and SR-highMass for simplified models with ${\tilde{\tau}}_L {\tilde{\tau}}_L$ only production.

Observed 95% CL exclusion limits for simplified models with direct stau pair production in SR-lowMass.

Expected 95% CL exclusion limits for simplified models with direct stau pair production in SR-lowMass.

Observed 95% CL exclusion limits for simplified models with direct stau pair production in SR-highMass.

Expected 95% CL exclusion limits for simplified models with direct stau pair production in SR-highMass.

Signal acceptance in SR highMass for combined stau final states

Signal acceptance in SR lowMass for combined stau final states

Signal efficiency in SR highMass for combined stau final states

Signal efficiency in SR lowMass for combined stau final states

Signal acceptance*efficiency in SR highMass for combined stau final states

Signal acceptance*efficiency in SR lowMass for combined stau final states

Cutflow for two reference points (${\tilde{\tau}}^{+}_{R,L} {\tilde{\tau}}^{-}_{R,L}$ production) in SR. The column labelled $N_{weighted}$ shows the results including all correction factors applied to simulation, and is normalised to 139 fb$^{-1}$, while $N_{raw}$ in brackets shows the results for the generated number of events. The quoted uncertainties are statistical only. The "Generator filter" includes the requirements that two $\tau$ in the event have ${p}_{T} > 15$ GeV and $|\eta| <$ 2.6. The "Baseline Cut" includes the requirement of two baseline $\tau$ with a minimum value at 0.01 of the boosted decision tree discriminant (JetBDTSigTransMin $>$ 0.01) and ${p}_{T, \tau_{1}} > 50$ GeV and ${p}_{T, \tau_{2}} > 40$ GeV. At the step "Trigger & offline cuts", the following requirements are applied: the event is recorded using the asymmetric di-$\tau$ trigger (di-$\tau$ $E_{T}^{miss}$ trigger) in SR-lowMass (SR-highMass), and the lepton $p_{T}$ and $E_{T}^{miss}$ are required at plateau.

Observed and expected numbers of events in the control and signal regions where all control and signal region bins are included as constraints in the likelihood. The expected event yields of SM processes are given after the background-only fit. The entries marked as "--" are negligible. The uncertainties correspond to the sum in quadrature of statistical and systematic uncertainties. The correlation of systematic uncertainties among control regions and among background processes is fully taken into account.

The post-fit $m_{T2}$ distribution for SR-lowMass. The stacked histograms show the expected SM backgrounds. The multi-jet contribution is estimated from data using the ABCD method. The contributions of multi-jet and $W$+jets events are scaled with the corresponding normalization factors derived from the background-only fit. The hatched bands represent the sum in quadrature of systematic and statistical uncertainties of the total SM background. For illustration, the distributions from the SUSY reference points are also shown as dashed lines. The last bin includes the overflow events.

The post-fit $m_{T2}$ distribution for SR-highMass. The stacked histograms show the expected SM backgrounds. The multi-jet contribution is estimated from data using the ABCD method. The contributions of multi-jet and $W$+jets events are scaled with the corresponding normalization factors derived from the background-only fit. The hatched bands represent the sum in quadrature of systematic and statistical uncertainties of the total SM background. For illustration, the distributions from the SUSY reference points are also shown as dashed lines. The last bin includes the overflow events.

The $m_{T2}$ post-fit distributions in the multi-jet background validation region VR-F (lowMass). The stacked histograms show the contribution of each relevant SM process. The multi-jet shape is taken from VR-E in the ABCD method and the normalization is determined by the transfer factor $T$ and rescaled by a correction factor determined by the fit. The hatched bands represent the combined statistical and systematic uncertainties in the sum of the SM backgrounds shown. For illustration, the distributions of the SUSY reference points are also shown as dashed lines.

The $m_{T2}$ post-fit distributions in the multi-jet background validation VR-F (highMass). The stacked histograms show the contribution of each relevant SM process. The multi-jet shape is taken from VR-E in the ABCD method and the normalization is determined by the transfer factor $T$ and rescaled by a correction factor determined by the fit. The hatched bands represent the combined statistical and systematic uncertainties in the sum of the SM backgrounds shown. For illustration, the distributions of the SUSY reference points are also shown as dashed lines.

The $E_{T}^{miss}$ post-fit distributions in the multi-jet background validation region VR-F (lowMass). The stacked histograms show the contribution of each relevant SM process. The multi-jet shape is taken from VR-E in the ABCD method and the normalization is determined by the transfer factor $T$ and rescaled by a correction factor determined by the fit. The hatched bands represent the combined statistical and systematic uncertainties in the sum of the SM backgrounds shown. For illustration, the distributions of the SUSY reference points are also shown as dashed lines.

The $E_{T}^{miss}$ post-fit distributions in the multi-jet background validation region VR-F (highMass). The stacked histograms show the contribution of each relevant SM process. The multi-jet shape is taken from VR-E in the ABCD method and the normalization is determined by the transfer factor $T$ and rescaled by a correction factor determined by the fit. The hatched bands represent the combined statistical and systematic uncertainties in the sum of the SM backgrounds shown. For illustration, the distributions of the SUSY reference points are also shown as dashed lines.

The pre-fit $m_{T2}$ distribution in the $WCR$. The SM backgrounds other than multi-jet production are estimated from MC simulation. The multi-jet contribution is estimated from data using the OS--SS method. The hatched bands represent the combined statistical and systematic uncertainties of the total SM background. For illustration, the distributions of the SUSY reference points are also shown as dashed lines.

The post-fit yields in the $WVR$, $TVRs$, $ZVRs$ and $VVVRs$. The SM backgrounds other than multi-jet production are estimated from MC simulation. The multi-jet contribution is negligible and is estimated from data using the ABCD method, using CRs obtained with the same technique used for the SRs. The hatched bands represent the combined statistical and systematic uncertainties of the total SM background. For illustration, the distributions of the SUSY reference points are also shown as dashed lines. The lower panels show the ratio of data to the SM background estimate.

Values of the predicted signal yields in each of the 15 signal regions (for $ x=0.75 $).

Pre-fit background-only covariance among different search bins.

Expected $95\%$ CL upper limit on the production cross section of top squarks (for $ x=0.25 $).

Top squark - neutralino mass points lying on the expected exclusion contour (for $ x=0.25 $).

Top squark - neutralino mass points lying on the expected-$1\sigma_\text{expected}$ exclusion contour (for $ x=0.25 $).

Top squark - neutralino mass points lying on the expected+$1\sigma_\text{expected}$ exclusion contour (for $ x=0.25 $).

Observed $95\%$ CL upper limit on the production cross section of top squarks (for $ x=0.25 $).

Top squark - neutralino mass points lying on the observed exclusion contour (for $ x=0.25 $).

Top squark - neutralino mass points lying on the observed-$1\sigma_\text{observed}$ exclusion contour (for $ x=0.25 $).

Top squark - neutralino mass points lying on the observed+$1\sigma_\text{observed}$ exclusion contour (for $ x=0.25 $).

Expected $95\%$ CL upper limit on the production cross section of top squarks (for $ x=0.5 $).

Top squark - neutralino mass points lying on the expected exclusion contour (for $ x=0.5 $).

Top squark - neutralino mass points lying on the expected-$1\sigma_\text{expected}$ exclusion contour (for $ x=0.5 $).

Top squark - neutralino mass points lying on the expected+$1\sigma_\text{expected}$ exclusion contour (for $ x=0.5 $).

Observed $95\%$ CL upper limit on the production cross section of top squarks (for $ x=0.5 $).

Top squark - neutralino mass points lying on the observed exclusion contour (for $ x=0.5 $).

Top squark - neutralino mass points lying on the observed-$1\sigma_\text{observed}$ exclusion contour (for $ x=0.5 $).

Top squark - neutralino mass points lying on the observed+$1\sigma_\text{observed}$ exclusion contour (for $ x=0.5 $).

Expected $95\%$ CL upper limit on the production cross section of top squarks (for $ x=0.75 $).

Top squark - neutralino mass points lying on the expected exclusion contour (for $ x=0.75 $).

Top squark - neutralino mass points lying on the expected-$1\sigma_\text{expected}$ exclusion contour (for $ x=0.75 $).

Top squark - neutralino mass points lying on the expected+$1\sigma_\text{expected}$ exclusion contour (for $ x=0.75 $).

Observed $95\%$ CL upper limit on the production cross section of top squarks (for $ x=0.75 $).

Top squark - neutralino mass points lying on the observed exclusion contour (for $ x=0.75 $).

Top squark - neutralino mass points lying on the observed-$1\sigma_\text{observed}$ exclusion contour (for $ x=0.75 $).

Top squark - neutralino mass points lying on the observed+$1\sigma_\text{observed}$ exclusion contour (for $ x=0.75 $).

Exclusion limits at 95% CL for direct gluino pair production, where the gluinos decay to light-flavor quarks and a $\tilde{\chi}_1^\pm$ that decays to $W^\pm\tilde{\chi}_1^0$. Signal cross sections are calculated at approximately NNLO+NNLL order in $\alpha_S$, assuming unity branching fraction to $q_i\bar{q}_j W\pm\tilde{\chi}_1^0$.

Exclusion limits at 95% CL for direct gluino pair production, where the gluinos decay to bottom quarks ($\tilde{g}\to b\bar{b}\tilde{\chi}_1^0$). Signal cross sections are calculated at approximately NNLO+NNLL order in $\alpha_S$, assuming unity branching fraction for the given decay.

Exclusion limits at 95% CL for direct gluino pair production, where the gluinos decay to top quarks ($\tilde{g}\to t\bar{t}\tilde{\chi}_1^0$). Signal cross sections are calculated at approximately NNLO+NNLL order in $\alpha_S$, assuming unity branching fraction for the given decay.

Exclusion limits at 95% CL for light-flavor squark pair production, where the squarks decay to $q\tilde{\chi}_1^0$. Signal cross sections are calculated at approximately NNLO+NNLL order in $\alpha_S$, assuming unity branching fraction for the given decay and 1-fold degeneracy in the light-flavor squarks (corresponding to the inner set of curves in the limit plot). To get the theory cross section for other N-fold degeneracy assumptions (e.g. 8-fold for the outer curves in the limit plot), just multiply by N.

Exclusion limits at 95% CL for bottom squark pair production, where the squarks decay to $b\tilde{\chi}_1^0$. Signal cross sections are calculated at approximately NNLO+NNLL order in $\alpha_S$, assuming unity branching fraction for the given decay.

Exclusion limits at 95% CL for top squark pair production, where the squarks decay to $t\tilde{\chi}_1^0$. Signal cross sections are calculated at approximately NNLO+NNLL order in $\alpha_S$, assuming unity branching fraction for the given decay.

Exclusion limits at 95% CL for top squark pair production, where the squarks decay to $b\tilde{\chi}_1^\pm$ and the $\tilde{\chi}_1^0$ decay to $W^\pm\tilde{\chi}_1^0$. Signal cross sections are calculated at approximately NNLO+NNLL order in $\alpha_S$, assuming unity branching fraction for the given decay.

Exclusion limits at 95% CL for top squark pair production, where the squarks decay either to $b\tilde{\chi}_1^\pm\to bW^\pm\tilde{\chi}_1^0$ or to $t\tilde{\chi}_1^0$. Signal cross sections are calculated at approximately NNLO+NNLL order in $\alpha_S$, assuming unity branching fraction for the given decay.

Exclusion limits at 95% CL for top squark pair production, where the squarks decay to $c\tilde{\chi}_1^0$. Signal cross sections are calculated at approximately NNLO+NNLL order in $\alpha_S$, assuming unity branching fraction for the given decay.

Exclusion limits at 95% CL for the mono-$\phi$ model, in which a resonantly-produced colored scalar decays to a massive Dirac fermion and a quark. Signal cross sections are calculated at leading order in $\alpha_S$, assuming unity branching fraction for the given decay.

Cross section limits for $\mathrm{LQ}\to\mathrm{q}\nu$, where $q=u,\,d,\,s,\,\mathrm{or}\,c$. Limits are at the 95% confidence level. Theory cross sections are LO for vector LQ, and NLO for scalar LQ. Branching ratio is assumed to be 100% to $\mathrm{q}\nu$.

Cross section limits for $\mathrm{LQ}\to\mathrm{b}\nu$. Limits are at the 95% confidence level. Theory cross sections are LO for vector LQ, and NLO for scalar LQ. Branching ratio is assumed to be 100% to $\mathrm{b}\nu$.

Cross section limits for $\mathrm{LQ}\to\mathrm{t}\nu$. Limits are at the 95% confidence level. Theory cross sections are LO for vector LQ, and NLO for scalar LQ. Branching ratios are assumed to be $\mathcal{B}(\mathrm{LQ}\to\mathrm{t}\nu)=1-\beta$, and $\mathcal{B}(\mathrm{LQ}\to\mathrm{b}\tau)=\beta$.

Predictions and observations for monojet signal regions

Predictions and observations for signal regions with $250 \leq H_\mathrm{T} < 450$ GeV

Predictions and observations for signal regions with $450 \leq H_\mathrm{T} < 575$ GeV and $N_\mathrm{j}<7$

Predictions and observations for signal regions with $450 \leq H_\mathrm{T} < 575$ GeV and $N_\mathrm{j}\geq7$

Predictions and observations for signal regions with $575 \leq H_\mathrm{T} < 1200$ GeV and $N_\mathrm{j}^\mathrm{hi}<4$

Predictions and observations for signal regions with $575 \leq H_\mathrm{T} < 1200$ GeV and $4\leq N_\mathrm{j}^\mathrm{hi}<7$

Predictions and observations for signal regions with $575 \leq H_\mathrm{T} < 1200$ GeV and $N_\mathrm{j}\geq7$

Predictions and observations for signal regions with $1200 \leq H_\mathrm{T} < 1500$ GeV and $N_\mathrm{j}^\mathrm{hi}<4$

Predictions and observations for signal regions with $1200 \leq H_\mathrm{T} < 1500$ GeV and $4\leq N_\mathrm{j}^\mathrm{hi}<7$

Predictions and observations for signal regions with $1200 \leq H_\mathrm{T} < 1500$ GeV and $N_\mathrm{j}\geq7$

Predictions and observations for signal regions with $H_\mathrm{T} \geq 1500$ GeV and $N_\mathrm{j}<7$

Predictions and observations for signal regions with $H_\mathrm{T} \geq 1500$ GeV and $N_\mathrm{j}\geq7$

Covariance matrix for the 282 signal regions of the inclusive $M_\mathrm{T2}$ search

Correlation matrix for the 282 signal regions of the inclusive $M_\mathrm{T2}$ search

Bin number definitions for the $M_\mathrm{T2}$ covariance and correlation matrices

Exclusion limits at 95% CL for direct gluino pair production, where the gluinos decay to light-flavor quarks and either the lightest neutralino, or the lightest chargino, and the chargino is long-lived with $c\tau_0 = 10$ cm and mass O(100) MeV greater than the neutralino's mass. Signal cross sections are calculated at approximately NNLO+NNLL order in $\alpha_S$.

Exclusion limits at 95% CL for direct gluino pair production, where the gluinos decay to light-flavor quarks and either the lightest neutralino, or the lightest chargino, and the chargino is long-lived with $c\tau_0 = 50$ cm and mass O(100) MeV greater than the neutralino's mass. Signal cross sections are calculated at approximately NNLO+NNLL order in $\alpha_S$.

Exclusion limits at 95% CL for direct gluino pair production, where the gluinos decay to light-flavor quarks and either the lightest neutralino, or the lightest chargino, and the chargino is long-lived with $c\tau_0 = 200$ cm and mass O(100) MeV greater than the neutralino's mass. Signal cross sections are calculated at approximately NNLO+NNLL order in $\alpha_S$.

Exclusion limits at 95% CL for direct light squark pair production, where the squarks decay to light-flavor quarks and either the lightest neutralino, or the lightest chargino, and the chargino is long-lived with $c\tau_0 = 10$ cm and mass O(100) MeV greater than the neutralino's mass. Signal cross sections are calculated at approximately NNLO+NNLL order in $\alpha_S$, for a single light squark.

Exclusion limits at 95% CL for direct light squark pair production, where the squarks decay to light-flavor quarks and either the lightest neutralino, or the lightest chargino, and the chargino is long-lived with $c\tau_0 = 50$ cm and mass O(100) MeV greater than the neutralino's mass. Signal cross sections are calculated at approximately NNLO+NNLL order in $\alpha_S$, for a single light squark.

Exclusion limits at 95% CL for direct light squark pair production, where the squarks decay to light-flavor quarks and either the lightest neutralino, or the lightest chargino, and the chargino is long-lived with $c\tau_0 = 200$ cm and mass O(100) MeV greater than the neutralino's mass. Signal cross sections are calculated at approximately NNLO+NNLL order in $\alpha_S$, for a single light squark.

Exclusion limits at 95% CL for direct stop pair production, where the stops decay to either a top and the lightest neutralino, or a bottom and the lightest chargino, and the chargino is long-lived with $c\tau_0 = 10$ cm and mass O(100) MeV greater than the neutralino's mass. Signal cross sections are calculated at approximately NNLO+NNLL order in $\alpha_S$.

Exclusion limits at 95% CL for direct stop pair production, where the stops decay to either a top and the lightest neutralino, or a bottom and the lightest chargino, and the chargino is long-lived with $c\tau_0 = 50$ cm and mass O(100) MeV greater than the neutralino's mass. Signal cross sections are calculated at approximately NNLO+NNLL order in $\alpha_S$.

Exclusion limits at 95% CL for direct stop pair production, where the stops decay to either a top and the lightest neutralino, or a bottom and the lightest chargino, and the chargino is long-lived with $c\tau_0 = 200$ cm and mass O(100) MeV greater than the neutralino's mass. Signal cross sections are calculated at approximately NNLO+NNLL order in $\alpha_S$.

The maximum chargino mass excluded at 95% CL for direct gluino pair production, where the gluinos decay to light-flavor quarks and either the lightest neutralino, or the lightest chargino, and the chargino mass is O(100) MeV greater than the neutralino's mass. The chargino's lifetime is varied from $c\tau_{0} = 1$ to 2000 cm while the gluino mass is fixed to 1900 GeV. Signal cross sections are calculated at approximately NNLO+NNLL order in $\alpha_S$. If all kinematically allowed chargino masses are excluded, the curves, including 68 and 95% expected, tend to overlap. At short decay lengths, horizontal exclusion lines are obtained from the inclusive analysis, as this is not affected by track reconstruction inefficiencies, which may arise when the chargino decays before the CMS tracker, and therefore shows better sensitivity to scenarios with very small lifetime compared to the disappearing track search, based on median expected limits.

The maximum chargino mass excluded at 95% CL for direct squark pair production, where the squarks decay to light-flavor quarks and either the lightest neutralino, or the lightest chargino, and the chargino mass is O(100) MeV greater than the neutralino's mass. The chargino's lifetime is varied from $c\tau_{0} = 1$ to 2000 cm while the squark mass is fixed to 900 GeV. Signal cross sections are calculated at approximately NNLO+NNLL order in $\alpha_S$, for a single light squark. If all kinematically allowed chargino masses are excluded, the curves, including 68 and 95% expected, tend to overlap. At short decay lengths, horizontal exclusion lines are obtained from the inclusive analysis, as this is not affected by track reconstruction inefficiencies, which may arise when the chargino decays before the CMS tracker, and therefore shows better sensitivity to scenarios with very small lifetime compared to the disappearing track search, based on median expected limits.

The maximum chargino mass excluded at 95% CL for direct squark pair production, where the squarks decay to light-flavor quarks and either the lightest neutralino, or the lightest chargino, and the chargino mass is O(100) MeV greater than the neutralino's mass. The chargino's lifetime is varied from $c\tau_{0} = 1$ to 2000 cm while the squark mass is fixed to 1500 GeV. Signal cross sections are calculated at approximately NNLO+NNLL order in $\alpha_S$, and the eight light squarks' masses are assumed to be degenerate. If all kinematically allowed chargino masses are excluded, the curves, including 68 and 95% expected, tend to overlap. At short decay lengths, horizontal exclusion lines are obtained from the inclusive analysis, as this is not affected by track reconstruction inefficiencies, which may arise when the chargino decays before the CMS tracker, and therefore shows better sensitivity to scenarios with very small lifetime compared to the disappearing track search, based on median expected limits.

Exclusion limits at 95% CL for direct stop pair production, where the stops decay to either a top and the lightest neutralino, or a bottom and the lightest chargino, and the chargino mass is O(100) MeV greater than the neutralino's mass. The chargino's lifetime is varied from $c\tau_{0} = 1$ to 2000 cm while the stop mass is fixed to 1000 GeV. Signal cross sections are calculated at approximately NNLO+NNLL order in $\alpha_S$. If all kinematically allowed chargino masses are excluded, the curves, including 68 and 95% expected, tend to overlap. At short decay lengths, horizontal exclusion lines are obtained from the inclusive analysis, as this is not affected by track reconstruction inefficiencies, which may arise when the chargino decays before the CMS tracker, and therefore shows better sensitivity to scenarios with very small lifetime compared to the disappearing track search, based on median expected limits.

Exclusion limits at 95% CL for direct gluino pair production, where the gluinos decay to light-flavor quarks and either the lightest neutralino, or the lightest chargino, and the chargino mass is O(100) MeV greater than the neutralino's mass. The chargino's lifetime is varied from $c\tau_{0} = 5$ to 1000 cm while the gluino mass is fixed to 1600 GeV and the neutralino's mass is fixed to 1575 GeV. Signal cross sections are calculated at approximately NNLO+NNLL order in $\alpha_S$.

Exclusion limits at 95% CL for direct squark pair production, where the squarks decay to light-flavor quarks and either the lightest neutralino, or the lightest chargino, and the chargino mass is O(100) MeV greater than the neutralino's mass. The chargino's lifetime is varied from $c\tau_{0} = 5$ to 1000 cm while the squark mass is fixed to 2000 GeV and the neutralino's mass is fixed to 1000 GeV. Signal cross sections are calculated at approximately NNLO+NNLL order in $\alpha_S$, and the eight light squarks' masses are assumed to be degenerate.

Exclusion limits at 95% CL for direct stop pair production, where the stops decay to either a top and the lightest neutralino, or a bottom and the lightest chargino, and the chargino mass is O(100) MeV greater than the neutralino's mass. The chargino's lifetime is varied from $c\tau_{0} = 5$ to 1000 cm while the stop mass is fixed to 1100 GeV and the neutralino's mass is fixed to 1000 GeV. Signal cross sections are calculated at approximately NNLO+NNLL order in $\alpha_S$.

Predictions and observations for 2016 disappearing track signal regions

Predictions and observations for 2017-2018 pixel track signal regions

Predictions and observations for 2017-2018 medium (M) and long (L) length track signal regions

Covariance matrix for the 68 signal regions of the disappearing tracks $M_\mathrm{T2}$ search

Correlation matrix for the 68 signal regions of the disappearing tracks $M_\mathrm{T2}$ search

Distributions of $m_{T2}$ in VR-top-low for data and the estimated SM backgrounds. The normalisation factors extracted from the corresponding CRs are used to rescale the $t\bar t$, single top, WW, WZ and ZZ backgrounds. The fake and non-prompt leptons background (FNP) is calculated using the data-driven matrix method. The uncertainty band includes all sources of systematic and statistical errors and the last bin includes the overflow.

Distributions of $m_{T2}$ in VR-top-high for data and the estimated SM backgrounds. The normalisation factors extracted from the corresponding CRs are used to rescale the $t\bar t$, single top, WW, WZ and ZZ backgrounds. The fake and non-prompt leptons background (FNP) is calculated using the data-driven matrix method. The uncertainty band includes all sources of systematic and statistical errors and the last bin includes the overflow.

Distributions of $E_T^{miss}$ in VR-WW-0J for data and the estimated SM backgrounds. The normalisation factors extracted from the corresponding CRs are used to rescale the $t\bar t$, single top, WW, WZ and ZZ backgrounds. The fake and non-prompt leptons background (FNP) is calculated using the data-driven matrix method. The uncertainty band includes all sources of systematic and statistical errors and the last bin includes the overflow.

Distributions of $E_T^{miss}$ in VR-WW-1J for data and the estimated SM backgrounds. The normalisation factors extracted from the corresponding CRs are used to rescale the $t\bar t$, single top, WW, WZ and ZZ backgrounds. The fake and non-prompt leptons background (FNP) is calculated using the data-driven matrix method. The uncertainty band includes all sources of systematic and statistical errors and the last bin includes the overflow.

Distributions of $E_T^{miss}$ significance in VR-VZ for data and the estimated SM backgrounds. The normalisation factors extracted from the corresponding CRs are used to rescale the $t\bar t$, single top, WW, WZ and ZZ backgrounds. The fake and non-prompt leptons background (FNP) is calculated using the data-driven matrix method. The uncertainty band includes all sources of systematic and statistical errors and the last bin includes the overflow.

Distributions of $E_T^{miss}$ significance in VR-top-WW for data and the estimated SM backgrounds. The normalisation factors extracted from the corresponding CRs are used to rescale the $t\bar t$, single top, WW, WZ and ZZ backgrounds. The fake and non-prompt leptons background (FNP) is calculated using the data-driven matrix method. The uncertainty band includes all sources of systematic and statistical errors and the last bin includes the overflow.

Breakdown of the dominant systematic uncertainties on background estimates in the inclusive SRs requiring $m_{T2}$>100 GeV after performing the profile likelihood fit. Note that the individual uncertainties can be correlated, and do not necessarily add up quadratically to the total background uncertainty. The percentages show the size of the uncertainty relative to the total expected background. "Top theoretical uncertainties" refers to $t\bar t$ theoretical uncertainties and the uncertainty associated to $Wt-t\bar t$ interference added quadratically.

Observed events and predicted post-fit background yields for the DF inclusive SRs. The model independent upper limits at 95% confidence level (CL) on the observed and expected number of beyond the SM events $S^{0.95}_{obs/exp}$ and the effective beyond the SM cross-section $\sigma^{0.95}_{obs}$ are also reported. The last row reports the $p_0$-value of the SM-only hypothesis. For SRs where the data yield is smaller than expected, the $p$-value is truncated at 0.50. The background denoted as "Other" in the Table includes the non-dominant background sources for this analysis, i.e. Z+jets, $t\bar t$+V, Higgs and Drell-Yan events. A "–" symbol indicates that the background contribution is negligible.

Observed events and predicted post-fit background yields for the DF inclusive SRs. The model independent upper limits at 95% confidence level (CL) on the observed and expected number of beyond the SM events $S^{0.95}_{obs/exp}$ and the effective beyond the SM cross-section $\sigma^{0.95}_{obs}$ are also reported. The last row reports the $p_0$-value of the SM-only hypothesis. For SRs where the data yield is smaller than expected, the $p$-value is truncated at 0.50. The background denoted as "Other" in the Table includes the non-dominant background sources for this analysis, i.e. Z+jets, $t\bar t$+V, Higgs and Drell-Yan events. A "–" symbol indicates that the background contribution is negligible.

Observed events and predicted post-fit background yields for the SF inclusive SRs. The model independent upper limits at 95% confidence level (CL) on the observed and expected number of beyond the SM events $S^{0.95}_{obs/exp}$ and the effective beyond the SM cross-section $\sigma^{0.95}_{obs}$ are also reported. The last row reports the $p_0$-value of the SM-only hypothesis. For SRs where the data yield is smaller than expected, the $p$-value is truncated at 0.50. The background denoted as "Other" in the Table includes the non-dominant background sources for this analysis, i.e. Z+jets, $t\bar t$+V, Higgs and Drell-Yan events. A "–" symbol indicates that the background contribution is negligible.

Observed events and predicted post-fit background yields for the SF inclusive SRs. The model independent upper limits at 95% confidence level (CL) on the observed and expected number of beyond the SM events $S^{0.95}_{obs/exp}$ and the effective beyond the SM cross-section $\sigma^{0.95}_{obs}$ are also reported. The last row reports the $p_0$-value of the SM-only hypothesis. For SRs where the data yield is smaller than expected, the $p$-value is truncated at 0.50. The background denoted as "Other" in the Table includes the non-dominant background sources for this analysis, i.e. Z+jets, $t\bar t$+V, Higgs and Drell-Yan events. A "–" symbol indicates that the background contribution is negligible.

Distributions of $m_{T2}$ in SRSF-0J for data and the estimated SM backgrounds. The normalisation factors extracted from the corresponding CRs are used to rescale the $t\bar t$, single top, WW, WZ and ZZ backgrounds. The fake and non-prompt leptons background (FNP) is calculated using the data-driven matrix method. The uncertainty band includes all sources of systematic and statistical errors and the last bin includes the overflow. Distributions for three benchmark signal points are overlaid for comparison.

Distributions of $m_{T2}$ in SRSF-1J for data and the estimated SM backgrounds. The normalisation factors extracted from the corresponding CRs are used to rescale the $t\bar t$, single top, WW, WZ and ZZ backgrounds. The fake and non-prompt leptons background (FNP) is calculated using the data-driven matrix method. The uncertainty band includes all sources of systematic and statistical errors and the last bin includes the overflow. Distributions for three benchmark signal points are overlaid for comparison.

Distributions of $m_{T2}$ in SRDF-0J for data and the estimated SM backgrounds. The normalisation factors extracted from the corresponding CRs are used to rescale the $t\bar t$, single top, WW, WZ and ZZ backgrounds. The fake and non-prompt leptons background (FNP) is calculated using the data-driven matrix method. The uncertainty band includes all sources of systematic and statistical errors and the last bin includes the overflow. Distributions for three benchmark signal points are overlaid for comparison.

Distributions of $m_{T2}$ in SRDF-1J for data and the estimated SM backgrounds. The normalisation factors extracted from the corresponding CRs are used to rescale the $t\bar t$, single top, WW, WZ and ZZ backgrounds. The fake and non-prompt leptons background (FNP) is calculated using the data-driven matrix method. The uncertainty band includes all sources of systematic and statistical errors and the last bin includes the overflow. Distributions for three benchmark signal points are overlaid for comparison.

Observed exclusion limits on SUSY simplified models for chargino-pair production with $W$ boson mediated decays. All limits are computed at 95% CL.

Expected exclusion limits on SUSY simplified models for chargino-pair production with $W$ boson mediated decays. All limits are computed at 95% CL.

Observed exclusion limits on SUSY simplified models for chargino-pair production with slepton/sneutrino mediated mediated decays. The mass relation $m(\tilde{l}_L)=\frac{1}{2}[m(\tilde{\chi}^{\pm}_1 + m(\tilde{\chi}^{0}_1)]$ is assumed. All limits are computed at 95% CL.

Expected exclusion limits on SUSY simplified models for chargino-pair production with slepton/sneutrino mediated mediated decays. The mass relation $m(\tilde{l}_L)=\frac{1}{2}[m(\tilde{\chi}^{\pm}_1 + m(\tilde{\chi}^{0}_1)]$ is assumed. All limits are computed at 95% CL.

Observed exclusion limits on SUSY simplified models for slepton-pair production. All limits are computed at 95% CL.

Expected exclusion limits on SUSY simplified models for slepton-pair production. All limits are computed at 95% CL.

Observed exclusion limits on SUSY simplified models for selectron-pair production, with left and right handed selectron production combined. All limits are computed at 95% CL.

Expected exclusion limits on SUSY simplified models for selectron-pair production, with left and right handed selectron production combined. All limits are computed at 95% CL.

Observed exclusion limits on SUSY simplified models for left-handed selectron-pair production. All limits are computed at 95% CL.

Expected exclusion limits on SUSY simplified models for left-handed selectron-pair production. All limits are computed at 95% CL.

Observed exclusion limits on SUSY simplified models for right-handed selectron-pair production. All limits are computed at 95% CL.

Expected exclusion limits on SUSY simplified models for right-handed selectron-pair production. All limits are computed at 95% CL.

Observed exclusion limits on SUSY simplified models for smuon-pair production, with left and right handed smuon production combined. All limits are computed at 95% CL.

Expected exclusion limits on SUSY simplified models for smuon-pair production, with left and right handed smuon production combined. All limits are computed at 95% CL.

Observed exclusion limits on SUSY simplified models for left-handed smuon-pair production. All limits are computed at 95% CL.

Expected exclusion limits on SUSY simplified models for left-handed smuon-pair production. All limits are computed at 95% CL.

Observed exclusion limits on SUSY simplified models for right-handed smuon-pair production. All limits are computed at 95% CL.

Expected exclusion limits on SUSY simplified models for right-handed smuon-pair production. All limits are computed at 95% CL.

Observed events and predicted background yields from the fit in the binned DF SRs with $n_{non-b-tagged jets}=0$. For backgrounds whose normalisation is extracted from the fit in the CRs, the yield expected from the simulation before the fit is also reported. The background denoted as "Other" in the Table includes the non-dominant background sources for this analysis, i.e. Z+jets, $t\bar t$ +V, Higgs and Drell-Yan events. A "–" symbol indicates that the background contribution is negligible.

Observed events and predicted background yields from the fit in the binned DF SRs with $n_{non-b-tagged jets}=1$. For backgrounds whose normalisation is extracted from the fit in the CRs, the yield expected from the simulation before the fit is also reported. The background denoted as "Other" in the Table includes the non-dominant background sources for this analysis, i.e. Z+jets, $t\bar t$ +V, Higgs and Drell-Yan events. A "–" symbol indicates that the background contribution is negligible.

Observed events and predicted background yields from the fit in the binned SF SRs with $n_{non-b-tagged jets}=0$. For backgrounds whose normalisation is extracted from the fit in the CRs, the yield expected from the simulation before the fit is also reported. The background denoted as "Other" in the Table includes the non-dominant background sources for this analysis, i.e. Z+jets, $t\bar t$ +V, Higgs and Drell-Yan events. A "–" symbol indicates that the background contribution is negligible.

Observed events and predicted background yields from the fit in the binned SF SRs with $n_{non-b-tagged jets}=1$. For backgrounds whose normalisation is extracted from the fit in the CRs, the yield expected from the simulation before the fit is also reported. The background denoted as "Other" in the Table includes the non-dominant background sources for this analysis, i.e. Z+jets, $t\bar t$ +V, Higgs and Drell-Yan events. A "–" symbol indicates that the background contribution is negligible.

Observed exclusion limits on SUSY simplified models for left-handed slepton-pair production. All limits are computed at 95% CL.

Expected exclusion limits on SUSY simplified models for left-handed slepton-pair production. All limits are computed at 95% CL.

Observed exclusion limits on SUSY simplified models for right-handed slepton-pair production. All limits are computed at 95% CL.

Expected exclusion limits on SUSY simplified models for right-handed slepton-pair production. All limits are computed at 95% CL.

Upper limits on signal cross-section (fb) for chargino-pair production with W -boson-mediated decays.

Upper limits on signal cross-section (fb) for chargino-pair production with slepton/sneutrino-mediated decays. The mass relation $m(\tilde{l}_L)=\frac{1}{2}[m(\tilde{\chi}^{\pm}_1 + m(\tilde{\chi}^{0}_1)]$ is assumed.

Upper limits on signal cross-section (fb) for slepton-pair production.

Signal Acceptance for direct chargino-pair production with W-boson mediated decays in SR-SF-0J-[100,inf).

Signal Efficiency for direct chargino-pair production with W-boson mediated decays in SR-SF-0J-[100,inf).

Signal Acceptance for direct chargino-pair production with W-boson mediated decays in SR-SF-1J-[100,inf).

Signal Efficiency for direct chargino-pair production with W-boson mediated decays in SR-SF-1J-[100,inf).

Signal Acceptance for direct chargino-pair production with W-boson mediated decays in SR-DF-0J-[100,inf).

Signal Efficiency for direct chargino-pair production with W-boson mediated decays in SR-DF-0J-[100,inf).

Signal Acceptance for direct chargino-pair production with W-boson mediated decays in SR-DF-1J-[100,inf).

Signal Efficiency for direct chargino-pair production with W-boson mediated decays in SR-DF-1J-[100,inf).

Signal Acceptance for direct chargino-pair production with W-boson mediated decays in SR-SF-0J-[160,inf).

Signal Efficiency for direct chargino-pair production with W-boson mediated decays in SR-SF-0J-[160,inf).

Signal Acceptance for direct chargino-pair production with W-boson mediated decays in SR-SF-1J-[160,inf).

Signal Efficiency for direct chargino-pair production with W-boson mediated decays in SR-SF-1J-[160,inf).

Signal Acceptance for direct chargino-pair production with W-boson mediated decays in SR-DF-0J-[160,inf).

Signal Efficiency for direct chargino-pair production with W-boson mediated decays in SR-DF-0J-[160,inf).

Signal Acceptance for direct chargino-pair production with W-boson mediated decays in SR-DF-1J-[160,inf).

Signal Efficiency for direct chargino-pair production with W-boson mediated decays in SR-DF-1J-[160,inf).

Signal Acceptance for direct chargino-pair production with W-boson mediated decays in SR-SF-0J-[100,120).

Signal Efficiency for direct chargino-pair production with W-boson mediated decays in SR-SF-0J-[100,120).

Signal Acceptance for direct chargino-pair production with W-boson mediated decays in SR-SF-1J-[100,120).

Signal Efficiency for direct chargino-pair production with W-boson mediated decays in SR-SF-1J-[100,120).

Signal Acceptance for direct chargino-pair production with W-boson mediated decays in SR-DF-0J-[100,120).

Signal Efficiency for direct chargino-pair production with W-boson mediated decays in SR-DF-0J-[100,120).

Signal Acceptance for direct chargino-pair production with W-boson mediated decays in SR-DF-1J-[100,120).

Signal Efficiency for direct chargino-pair production with W-boson mediated decays in SR-DF-1J-[100,120).

Signal Acceptance for direct chargino-pair production with W-boson mediated decays in SR-SF-0J-[120,160).

Signal Efficiency for direct chargino-pair production with W-boson mediated decays in SR-SF-0J-[120,160).

Signal Acceptance for direct chargino-pair production with W-boson mediated decays in SR-SF-1J-[120,160).

Signal Efficiency for direct chargino-pair production with W-boson mediated decays in SR-SF-1J-[120,160).

Signal Acceptance for direct chargino-pair production with W-boson mediated decays in SR-DF-0J-[120,160).

Signal Efficiency for direct chargino-pair production with W-boson mediated decays in SR-DF-0J-[120,160).

Signal Acceptance for direct chargino-pair production with W-boson mediated decays in SR-DF-1J-[120,160).

Signal Efficiency for direct chargino-pair production with W-boson mediated decays in SR-DF-1J-[120,160).

Signal Acceptance for direct chargino-pair production with W-boson mediated decays in SR-SF-0J-[100,105).

Signal Efficiency for direct chargino-pair production with W-boson mediated decays in SR-SF-0J-[100,105).

Signal Acceptance for direct chargino-pair production with W-boson mediated decays in SR-SF-1J-[100,105).

Signal Efficiency for direct chargino-pair production with W-boson mediated decays in SR-SF-1J-[100,105).

Signal Acceptance for direct chargino-pair production with W-boson mediated decays in SR-DF-0J-[100,105).

Signal Efficiency for direct chargino-pair production with W-boson mediated decays in SR-DF-0J-[100,105).

Signal Acceptance for direct chargino-pair production with W-boson mediated decays in SR-DF-1J-[100,105).

Signal Efficiency for direct chargino-pair production with W-boson mediated decays in SR-DF-1J-[100,105).

Signal Acceptance for direct chargino-pair production with W-boson mediated decays in SR-SF-0J-[105,110).

Signal Efficiency for direct chargino-pair production with W-boson mediated decays in SR-SF-0J-[105,110).

Signal Acceptance for direct chargino-pair production with W-boson mediated decays in SR-SF-1J-[105,110).

Signal Efficiency for direct chargino-pair production with W-boson mediated decays in SR-SF-1J-[105,110).

Signal Acceptance for direct chargino-pair production with W-boson mediated decays in SR-DF-0J-[105,110).

Signal Efficiency for direct chargino-pair production with W-boson mediated decays in SR-DF-0J-[105,110).

Signal Acceptance for direct chargino-pair production with W-boson mediated decays in SR-DF-1J-[105,110).

Signal Efficiency for direct chargino-pair production with W-boson mediated decays in SR-DF-1J-[105,110).

Signal Acceptance for direct chargino-pair production with W-boson mediated decays in SR-SF-0J-[110,120).

Signal Efficiency for direct chargino-pair production with W-boson mediated decays in SR-SF-0J-[110,120).

Signal Acceptance for direct chargino-pair production with W-boson mediated decays in SR-SF-1J-[110,120).

Signal Efficiency for direct chargino-pair production with W-boson mediated decays in SR-SF-1J-[110,120).

Signal Acceptance for direct chargino-pair production with W-boson mediated decays in SR-DF-0J-[110,120).

Signal Efficiency for direct chargino-pair production with W-boson mediated decays in SR-DF-0J-[110,120).

Signal Acceptance for direct chargino-pair production with W-boson mediated decays in SR-DF-1J-[110,120).

Signal Efficiency for direct chargino-pair production with W-boson mediated decays in SR-DF-1J-[110,120).

Signal Acceptance for direct chargino-pair production with W-boson mediated decays in SR-SF-0J-[120,140).

Signal Efficiency for direct chargino-pair production with W-boson mediated decays in SR-SF-0J-[120,140).

Signal Acceptance for direct chargino-pair production with W-boson mediated decays in SR-SF-1J-[120,140).

Signal Efficiency for direct chargino-pair production with W-boson mediated decays in SR-SF-1J-[120,140).

Signal Acceptance for direct chargino-pair production with W-boson mediated decays in SR-DF-0J-[120,140).

Signal Efficiency for direct chargino-pair production with W-boson mediated decays in SR-DF-0J-[120,140).

Signal Acceptance for direct chargino-pair production with W-boson mediated decays in SR-DF-1J-[120,140).

Signal Efficiency for direct chargino-pair production with W-boson mediated decays in SR-DF-1J-[120,140).

Signal Acceptance for direct chargino-pair production with W-boson mediated decays in SR-SF-0J-[140,160).

Signal Efficiency for direct chargino-pair production with W-boson mediated decays in SR-SF-0J-[140,160).

Signal Acceptance for direct chargino-pair production with W-boson mediated decays in SR-SF-1J-[140,160).

Signal Efficiency for direct chargino-pair production with W-boson mediated decays in SR-SF-1J-[140,160).

Signal Acceptance for direct chargino-pair production with W-boson mediated decays in SR-DF-0J-[140,160).

Signal Efficiency for direct chargino-pair production with W-boson mediated decays in SR-DF-0J-[140,160).

Signal Acceptance for direct chargino-pair production with W-boson mediated decays in SR-DF-1J-[140,160).

Signal Efficiency for direct chargino-pair production with W-boson mediated decays in SR-DF-1J-[140,160).

Signal Acceptance for direct chargino-pair production with W-boson mediated decays in SR-SF-0J-[160,180).

Signal Efficiency for direct chargino-pair production with W-boson mediated decays in SR-SF-0J-[160,180).

Signal Acceptance for direct chargino-pair production with W-boson mediated decays in SR-SF-1J-[160,180).

Signal Efficiency for direct chargino-pair production with W-boson mediated decays in SR-SF-1J-[160,180).

Signal Acceptance for direct chargino-pair production with W-boson mediated decays in SR-DF-0J-[160,180).

Signal Efficiency for direct chargino-pair production with W-boson mediated decays in SR-DF-0J-[160,180).

Signal Acceptance for direct chargino-pair production with W-boson mediated decays in SR-DF-1J-[160,180).

Signal Efficiency for direct chargino-pair production with W-boson mediated decays in SR-DF-1J-[160,180).

Signal Acceptance for direct chargino-pair production with W-boson mediated decays in SR-SF-0J-[180,220).

Signal Efficiency for direct chargino-pair production with W-boson mediated decays in SR-SF-0J-[180,220).

Signal Acceptance for direct chargino-pair production with W-boson mediated decays in SR-SF-1J-[180,220).

Signal Efficiency for direct chargino-pair production with W-boson mediated decays in SR-SF-1J-[180,220).

Signal Acceptance for direct chargino-pair production with W-boson mediated decays in SR-DF-0J-[180,220).

Signal Efficiency for direct chargino-pair production with W-boson mediated decays in SR-DF-0J-[180,220).

Signal Acceptance for direct chargino-pair production with W-boson mediated decays in SR-DF-1J-[180,220).

Signal Efficiency for direct chargino-pair production with W-boson mediated decays in SR-DF-1J-[180,220).

Signal Acceptance for direct chargino-pair production with W-boson mediated decays in SR-SF-0J-[220,260).

Signal Efficiency for direct chargino-pair production with W-boson mediated decays in SR-SF-0J-[220,260).

Signal Acceptance for direct chargino-pair production with W-boson mediated decays in SR-SF-1J-[220,260).

Signal Efficiency for direct chargino-pair production with W-boson mediated decays in SR-SF-1J-[220,260).

Signal Acceptance for direct chargino-pair production with W-boson mediated decays in SR-DF-0J-[220,260).

Signal Efficiency for direct chargino-pair production with W-boson mediated decays in SR-DF-0J-[220,260).

Signal Acceptance for direct chargino-pair production with W-boson mediated decays in SR-DF-1J-[220,260).

Signal Efficiency for direct chargino-pair production with W-boson mediated decays in SR-DF-1J-[220,260).

Signal Acceptance for direct chargino-pair production with W-boson mediated decays in SR-SF-0J-[260,inf).

Signal Efficiency for direct chargino-pair production with W-boson mediated decays in SR-SF-0J-[260,inf).

Signal Acceptance for direct chargino-pair production with W-boson mediated decays in SR-SF-1J-[260,inf).

Signal Efficiency for direct chargino-pair production with W-boson mediated decays in SR-SF-1J-[260,inf).

Signal Acceptance for direct chargino-pair production with W-boson mediated decays in SR-DF-0J-[260,inf).

Signal Efficiency for direct chargino-pair production with W-boson mediated decays in SR-DF-0J-[260,inf).

Signal Acceptance for direct chargino-pair production with W-boson mediated decays in SR-DF-1J-[260,inf).

Signal Efficiency for direct chargino-pair production with W-boson mediated decays in SR-DF-1J-[260,inf).

Cutflow for supersymmetric model where $\tilde{\chi}_1^{\pm}\tilde{\chi}_1^{\mp}$ decay via $W^{\pm}W^{\mp}$. The masses of the two charginos are 300 GeV, while the mass of $\tilde{\chi}_1^{0}$ is 50 GeV. The numbers are normalised to the luminosity of 139~fb$^{-1}$.

Cutflow for supersymmetric model where $\tilde{\chi}_1^{\pm}\tilde{\chi}_1^{\mp}$ decay via slepton-neutrino/sneutrino-lepton pair. The masses of the two charginos are 600 GeV, while the mass of $\tilde{\chi}_1^{0}$ is 1 GeV. The slepton/sneutrino masses are 300 GeV. The numbers are normalised to the luminosity of 139~fb$^{-1}$.

Cutflow for supersymmetric model where $\tilde\ell\tilde\ell$ are produced. Only $\tilde{e}$ and $\tilde{\mu}$ are considered in this model. The masses of the two sleptons are 400 GeV, while the mass of $\tilde{\chi}_1^{0}$ is 200 GeV. The numbers are normalised to the luminosity of 139~fb$^{-1}$.

Distributions of MaxMin alternative algorithm $m(h_{\mathrm{cand1}},h_{\mathrm{cand2}})_{\mathrm{avg}}$ after the background-only fit. The backgrounds which contribute only a small amount (diboson, W+jets and ttbar+W/Z/h) are grouped and labelled as `Other'.

Distributions of Leading jet $p_T$ after the background-only fit. The backgrounds which contribute only a small amount (diboson, W+jets and ttbar+W/Z/h) are grouped and labelled as `Other'.

Distributions of MaxMin algorithm $m_{hcand}$ after the background-only fit. The backgrounds which contribute only a small amount (diboson, W+jets and ttbar+W/Z/h) are grouped and labelled as `Other'.

Signal efficiency in SRA_M_m60 for simplified models with '$\widetilde{b}\widetilde{b}$ $\rightarrow$ $b\bar{b} \widetilde{\chi}_2^0 \widetilde{\chi}_2^0$ $\rightarrow$ hh$\widetilde{\chi}_1^0 \widetilde{\chi}_1^0$, h $\rightarrow$ $b\bar{b}$' production

Signal acceptance in SRC_28 for simplified models with $\widetilde{b}\widetilde{b}$ $\rightarrow$ $b\bar{b} \widetilde{\chi}_2^0 \widetilde{\chi}_2^0$ $\rightarrow$ hh$\widetilde{\chi}_1^0 \widetilde{\chi}_1^0$, h $\rightarrow$ $b\bar{b}$ production

Signal acceptance in SRC_26 for simplified models with $\widetilde{b}\widetilde{b}$ $\rightarrow$ $b\bar{b} \widetilde{\chi}_2^0 \widetilde{\chi}_2^0$ $\rightarrow$ hh$\widetilde{\chi}_1^0 \widetilde{\chi}_1^0$, h $\rightarrow$ $b\bar{b}$ production

Signal acceptance in SRC_24 for simplified models with $\widetilde{b}\widetilde{b}$ $\rightarrow$ $b\bar{b} \widetilde{\chi}_2^0 \widetilde{\chi}_2^0$ $\rightarrow$ hh$\widetilde{\chi}_1^0 \widetilde{\chi}_1^0$, h $\rightarrow$ $b\bar{b}$ production

Signal acceptance in SRA_M_dm130 for simplified models with $\widetilde{b}\widetilde{b}$ $\rightarrow$ $b\bar{b} \widetilde{\chi}_2^0 \widetilde{\chi}_2^0$ $\rightarrow$ hh$\widetilde{\chi}_1^0 \widetilde{\chi}_1^0$, h $\rightarrow$ $b\bar{b}$ production

Signal acceptance in SRB for simplified models with $\widetilde{b}\widetilde{b}$ $\rightarrow$ $b\bar{b} \widetilde{\chi}_2^0 \widetilde{\chi}_2^0$ $\rightarrow$ hh$\widetilde{\chi}_1^0 \widetilde{\chi}_1^0$, h $\rightarrow$ $b\bar{b}$ production

Signal acceptance in SRA_L_dm130 for simplified models with $\widetilde{b}\widetilde{b}$ $\rightarrow$ $b\bar{b} \widetilde{\chi}_2^0 \widetilde{\chi}_2^0$ $\rightarrow$ hh$\widetilde{\chi}_1^0 \widetilde{\chi}_1^0$, h $\rightarrow$ $b\bar{b}$ production

Signal acceptance in SRC_incl for simplified models with $\widetilde{b}\widetilde{b}$ $\rightarrow$ $b\bar{b} \widetilde{\chi}_2^0 \widetilde{\chi}_2^0$ $\rightarrow$ hh$\widetilde{\chi}_1^0 \widetilde{\chi}_1^0$, h $\rightarrow$ $b\bar{b}$ production

Signal acceptance in SRA_L_m60 for simplified models with $\widetilde{b}\widetilde{b}$ $\rightarrow$ $b\bar{b} \widetilde{\chi}_2^0 \widetilde{\chi}_2^0$ $\rightarrow$ hh$\widetilde{\chi}_1^0 \widetilde{\chi}_1^0$, h $\rightarrow$ $b\bar{b}$ production

Signal acceptance in SRA_incl_dm130 for simplified models with $\widetilde{b}\widetilde{b}$ $\rightarrow$ $b\bar{b} \widetilde{\chi}_2^0 \widetilde{\chi}_2^0$ $\rightarrow$ hh$\widetilde{\chi}_1^0 \widetilde{\chi}_1^0$, h $\rightarrow$ $b\bar{b}$ production

Signal acceptance in SRA_incl_m60 for simplified models with $\widetilde{b}\widetilde{b}$ $\rightarrow$ $b\bar{b} \widetilde{\chi}_2^0 \widetilde{\chi}_2^0$ $\rightarrow$ hh$\widetilde{\chi}_1^0 \widetilde{\chi}_1^0$, h $\rightarrow$ $b\bar{b}$ production

Signal efficiency in SRA_H_m60 for simplified models with '$\widetilde{b}\widetilde{b}$ $\rightarrow$ $b\bar{b} \widetilde{\chi}_2^0 \widetilde{\chi}_2^0$ $\rightarrow$ hh$\widetilde{\chi}_1^0 \widetilde{\chi}_1^0$, h $\rightarrow$ $b\bar{b}$' production

Signal efficiency in SRA_L_dm130 for simplified models with '$\widetilde{b}\widetilde{b}$ $\rightarrow$ $b\bar{b} \widetilde{\chi}_2^0 \widetilde{\chi}_2^0$ $\rightarrow$ hh$\widetilde{\chi}_1^0 \widetilde{\chi}_1^0$, h $\rightarrow$ $b\bar{b}$' production

Signal efficiency in SRB for simplified models with '$\widetilde{b}\widetilde{b}$ $\rightarrow$ $b\bar{b} \widetilde{\chi}_2^0 \widetilde{\chi}_2^0$ $\rightarrow$ hh$\widetilde{\chi}_1^0 \widetilde{\chi}_1^0$, h $\rightarrow$ $b\bar{b}$' production

Signal acceptance in SRC_22 for simplified models with $\widetilde{b}\widetilde{b}$ $\rightarrow$ $b\bar{b} \widetilde{\chi}_2^0 \widetilde{\chi}_2^0$ $\rightarrow$ hh$\widetilde{\chi}_1^0 \widetilde{\chi}_1^0$, h $\rightarrow$ $b\bar{b}$ production

Signal efficiency in SRA_H_dm130 for simplified models with '$\widetilde{b}\widetilde{b}$ $\rightarrow$ $b\bar{b} \widetilde{\chi}_2^0 \widetilde{\chi}_2^0$ $\rightarrow$ hh$\widetilde{\chi}_1^0 \widetilde{\chi}_1^0$, h $\rightarrow$ $b\bar{b}$' production

Signal efficiency in SRC_24 for simplified models with '$\widetilde{b}\widetilde{b}$ $\rightarrow$ $b\bar{b} \widetilde{\chi}_2^0 \widetilde{\chi}_2^0$ $\rightarrow$ hh$\widetilde{\chi}_1^0 \widetilde{\chi}_1^0$, h $\rightarrow$ $b\bar{b}$' production

Signal efficiency in SRC_26 for simplified models with '$\widetilde{b}\widetilde{b}$ $\rightarrow$ $b\bar{b} \widetilde{\chi}_2^0 \widetilde{\chi}_2^0$ $\rightarrow$ hh$\widetilde{\chi}_1^0 \widetilde{\chi}_1^0$, h $\rightarrow$ $b\bar{b}$' production

Signal acceptance in SRA_H_m60 for simplified models with $\widetilde{b}\widetilde{b}$ $\rightarrow$ $b\bar{b} \widetilde{\chi}_2^0 \widetilde{\chi}_2^0$ $\rightarrow$ hh$\widetilde{\chi}_1^0 \widetilde{\chi}_1^0$, h $\rightarrow$ $b\bar{b}$ production

Signal efficiency in SRA_incl_m60 for simplified models with '$\widetilde{b}\widetilde{b}$ $\rightarrow$ $b\bar{b} \widetilde{\chi}_2^0 \widetilde{\chi}_2^0$ $\rightarrow$ hh$\widetilde{\chi}_1^0 \widetilde{\chi}_1^0$, h $\rightarrow$ $b\bar{b}$' production

Signal efficiency in SRC_22 for simplified models with '$\widetilde{b}\widetilde{b}$ $\rightarrow$ $b\bar{b} \widetilde{\chi}_2^0 \widetilde{\chi}_2^0$ $\rightarrow$ hh$\widetilde{\chi}_1^0 \widetilde{\chi}_1^0$, h $\rightarrow$ $b\bar{b}$' production

Signal acceptance in SRA_M_m60 for simplified models with $\widetilde{b}\widetilde{b}$ $\rightarrow$ $b\bar{b} \widetilde{\chi}_2^0 \widetilde{\chi}_2^0$ $\rightarrow$ hh$\widetilde{\chi}_1^0 \widetilde{\chi}_1^0$, h $\rightarrow$ $b\bar{b}$ production

Signal efficiency in SRC_28 for simplified models with '$\widetilde{b}\widetilde{b}$ $\rightarrow$ $b\bar{b} \widetilde{\chi}_2^0 \widetilde{\chi}_2^0$ $\rightarrow$ hh$\widetilde{\chi}_1^0 \widetilde{\chi}_1^0$, h $\rightarrow$ $b\bar{b}$' production

Signal acceptance in SRA_H_dm130 for simplified models with $\widetilde{b}\widetilde{b}$ $\rightarrow$ $b\bar{b} \widetilde{\chi}_2^0 \widetilde{\chi}_2^0$ $\rightarrow$ hh$\widetilde{\chi}_1^0 \widetilde{\chi}_1^0$, h $\rightarrow$ $b\bar{b}$ production

Signal efficiency in SRA_incl_dm130 for simplified models with '$\widetilde{b}\widetilde{b}$ $\rightarrow$ $b\bar{b} \widetilde{\chi}_2^0 \widetilde{\chi}_2^0$ $\rightarrow$ hh$\widetilde{\chi}_1^0 \widetilde{\chi}_1^0$, h $\rightarrow$ $b\bar{b}$' production

Signal efficiency in SRA_L_m60 for simplified models with '$\widetilde{b}\widetilde{b}$ $\rightarrow$ $b\bar{b} \widetilde{\chi}_2^0 \widetilde{\chi}_2^0$ $\rightarrow$ hh$\widetilde{\chi}_1^0 \widetilde{\chi}_1^0$, h $\rightarrow$ $b\bar{b}$' production

Signal efficiency in SRA_M_dm130 for simplified models with '$\widetilde{b}\widetilde{b}$ $\rightarrow$ $b\bar{b} \widetilde{\chi}_2^0 \widetilde{\chi}_2^0$ $\rightarrow$ hh$\widetilde{\chi}_1^0 \widetilde{\chi}_1^0$, h $\rightarrow$ $b\bar{b}$' production

Signal efficiency in SRC_incl for simplified models with '$\widetilde{b}\widetilde{b}$ $\rightarrow$ $b\bar{b} \widetilde{\chi}_2^0 \widetilde{\chi}_2^0$ $\rightarrow$ hh$\widetilde{\chi}_1^0 \widetilde{\chi}_1^0$, h $\rightarrow$ $b\bar{b}$' production

Observed 95% CLs exclusion limit for the $\Delta M(\widetilde{\chi}_{2}^{0},\widetilde{\chi}_{1}^{0})$=130GeV signal grid for the best combined signal regions.

Expected 95% CLs exclusion limit for the $\Delta M(\widetilde{\chi}_{2}^{0},\widetilde{\chi}_{1}^{0})$=130GeV signal grid for the best combined signal regions.

Observed 95% CLs exclusion limit for the $M(\widetilde{\chi}_{1}^{0})$=60GeV signal grid for the best combined signal regions.

Expected 95% CLs exclusion limit for the $M(\widetilde{\chi}_{1}^{0})$=60GeV signal grid for the best combined signal regions.

Model dependent upper limit on the best combined signal regions considered in the $\Delta M(\widetilde{\chi}_{2}^{0},\widetilde{\chi}_{1}^{0})$=130GeV signal grid

Model dependet upper limits on the best combined signal regions considered in the $\Delta M(\widetilde{\chi}_{2}^{0},\widetilde{\chi}_{1}^{0})$=130GeV signal grid

Result of background only fit applied to signal regions. Event yields from the signal regions compared with SM MC predictions for the 3 highest contributing backgrounds separately and combined minor backgrounds.

Expected background event yields and dominant systematic uncertainties on background estimates in the A-type (inclusive), B-type and C-type (inclusive) regions.

Background-only fit results for the A- and B-type regions performed using 139$fb^{-1}$ of data. The quoted uncertainties on the fitted SM background include both the statistical and systematic uncertainties.

Background-only fit results for the C-type region performed using 139$fb^{-1}$ of data. The quoted uncertainties on the fitted SM background include both the statistical and systematic uncertainties.

Observed 95% CL upper limits on the visible cross sections σvis, the observed (S95obs) and expected (S95exp) 95% CL upper limits on the number of signal events with ± 1 σ excursions of the expectation, the CL of the background-only hypothesis, CLB, the discovery p-value (p0), truncated at 0.5, and the associated significance.

Cutflow of the MC events scaled to 139 $fb^{-1}$ for the SRA selections, with a scalar bottom signal of m$(\widetilde{b}_{1},\widetilde{\chi}_2^0,\widetilde{\chi}_1^0) = (1100, 330, 200)$ GeV, considered.

Cutflow of the MC events scaled to 139 $fb^{-1}$ for the SRB selections, with a scalar bottom signal of m$(\widetilde{b}_{1},\widetilde{\chi}_2^0,\widetilde{\chi}_1^0) = (700, 680, 550)$ GeV, considered.

Cutflow of the MC events scaled to 139 $fb^{-1}$ for the SRC selections, with a scalar bottom signal of m$(\widetilde{b}_{1},\widetilde{\chi}_2^0,\widetilde{\chi}_1^0) = (1200, 1150, 60)$ GeV, considered.

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

Expected limit contour for GGM scenario

Observed limit contour for GGM scenario

Expected limit contour for GGM scenario

Observed limit contour for GGM scenario

Expected limit contour for Neutralino BF scenario

Observed limit contour for Neutralino BF scenario

Expected limit contour for Chargino BF scenario

Observed limit contour for Chargino BF scenario

Expected limit contour for T5Wg

Observed limit contour for T5Wg

Expected limit contour for Gluino BF scenario

Observed limit contour for Gluino BF scenario

Covariance matrix for all bins (additional material)

Correlation matrix for all bins (additional material)

The observed and expected 95% CL upper limits on model-independent visible cross-sections, along with the observed $p0$ values, for the stable signal region, as a function of different mass windows, for which the lower bound is shown. The upper boundary on the mass window is 5 TeV for all windows.

The observed and expected 95% CL upper limits on model-independent visible cross-sections, along with the observed $p0$ values, for the metastable signal region, as a function of different mass windows, for which the lower bound is shown. The upper boundary on the mass window is 5 TeV for all windows.

For each gluino lifetime and mass in the signal samples, the lower boundary of the mass window in which at least $70\%$ of the reconstructed signal appears. The upper boundary for all mass windows is 5 TeV.

Acceptance and efficiency for a representative set of pair-produced gluino signal samples. The mass of the gluino ($m(\tilde{g})$), its lifetime ($\tau(\tilde{g})$) and the mass of the neutralino ($m(\tilde{\chi}^{0})$) are given in the first three columns. The Pythia 6.4.27 signal samples shown in this table are not reweighted to match the transverse momentum of the gluino-gluino system as simulated by MadGraph5_aMC@NLO. The acceptance is defined as the fraction of events passing a loose set of fiducial requirements. The full simulation efficiency (Full sim. $\epsilon$) is defined as the ratio of the number of reconstructed events, as expected by the full ATLAS simulation, and the number of events passing the fiducial requirements. The parameterised simulation efficiency (Param. sim. $\epsilon$) is defined as the ratio of the number of events estimated using a set of parametrised efficiencies (see auxiliary Figures 9,10,11,12) and the number of events passing the fiducial requirements alone.

The reconstructed candidate track mass distributions for observed data, predicted background, and the expected contribution from two signal models in the metastable R-hadron signal region. The yellow band around the background estimation includes both the statistical and systematic uncertainties.

The reconstructed candidate track mass distributions for observed data, predicted background, and the expected contribution from two signal models in the stable R-hadron signal region. The yellow band around the background estimation includes both the statistical and systematic uncertainties.

The 95% CL upper limit on the cross-section as a function of mass for gluinos with lifetime $\tau = 10$ ns decaying into $q\bar{q}$ and a 100 GeV neutralino, with the observed limit shown as a solid black line. The predicted production cross-section values are shown in purple along with their uncertainty. The expected upper limit in the case of only background is shown by the dashed black line, with a green $\pm 1\sigma$ and a yellow $\pm 2\sigma$ band.

The 95% CL upper limit on the cross-section as a function of mass for stable gluino $R$-hadrons, with the observed limit shown as a solid black line. The predicted production cross-section values are shown in purple along with their uncertainty. The expected upper limit in the case of only background is shown by the dashed black line, with a green $\pm 1\sigma$ and a yellow $\pm 2\sigma$ band.

Observed 95% lower limits on the gluino mass in the gluino lifetime--mass plane. The excluded area is to the left of the curves.

Expected 95% lower limits on the gluino mass in the gluino lifetime--mass plane. The excluded area is to the left of the curves.

The 95% CL upper limit on the cross-section as a function of mass for gluinos with lifetime $\tau = 1$ ns decaying into $q\bar{q}$ and a 100 GeV neutralino, with the observed limit shown as a solid black line. The predicted production cross-section values are shown in purple along with their uncertainty. The expected upper limit in the case of only background is shown by the dashed black line, with a green $\pm 1\sigma$ and a yellow $\pm 2\sigma$ band.

The 95% CL upper limit on the cross-section as a function of mass for gluinos with lifetime $\tau = 3$ ns decaying into $q\bar{q}$ and a 100 GeV neutralino, with the observed limit shown as a solid black line. The predicted production cross-section values are shown in purple along with their uncertainty. The expected upper limit in the case of only background is shown by the dashed black line, with a green $\pm 1\sigma$ and a yellow $\pm 2\sigma$ band.

The 95% CL upper limit on the cross-section as a function of mass for gluinos with lifetime $\tau = 30$ ns decaying into $q\bar{q}$ and a 100 GeV neutralino, with the observed limit shown as a solid black line. The predicted production cross-section values are shown in purple along with their uncertainty. The expected upper limit in the case of only background is shown by the dashed black line, with a green $\pm 1\sigma$ and a yellow $\pm 2\sigma$ band.

The 95% CL upper limit on the cross-section as a function of mass for gluinos with lifetime $\tau = 50$ ns decaying into $q\bar{q}$ and a 100 GeV neutralino, with the observed limit shown as a solid black line. The predicted production cross-section values are shown in purple along with their uncertainty. The expected upper limit in the case of only background is shown by the dashed black line, with a green $\pm 1\sigma$ and a yellow $\pm 2\sigma$ band.

The relationship between generated and reconstructed mass for gluino $R$-hadrons. Above 1500 GeV, the reconstructed mass falls below the generated mass due to bias in the reconstructed momentum. The uncertainty on the reconstructed mass is dominated by momentum uncertainty. The black dots represent the reconstructed mass computed as the most probable value of a Gaussian fit function, with the error bars showing its statistical uncertainty, while the orange band is the full-width at half maximum of the reconstructed mass distribution.

The parameterised efficiency for events to pass metastable event selections (including trigger, E$_{T}^{miss}$, and event cleaning requirements) as a function of the true E$_{T}^{miss}$ in the system, which is calculated at generator level. Event-level efficiencies are evaluated for events which have at least true E$_{T}^{miss} > 50$ GeV. The metastable event efficiencies are evaluated for different radial regions depending on the smallest radial distance, R, at which an R-hadron decays in the detector.

The parameterised efficiency for events to pass metastable event selections (including trigger, E$_{T}^{miss}$, and event cleaning requirements) as a function of the true E$_{T}^{miss}$ in the system, which is calculated at generator level. Event-level efficiencies are evaluated for events which have at least true E$_{T}^{miss} > 50$ GeV. The stable event efficiencies are evaluated for samples in which no R-hadron decays within the detector.

The parameterised efficiency for particles to pass full track selections in the metastable signal region, as function of the particle’s $\beta$, in different bins of transverse momentum, $p_{T}$, and for different radial decay positions of the particle. The efficiency is evaluated for particles which pass a loose set of fiducial requirements at generator level.

The parameterised efficiency for particles to pass full track selections in the metastable signal region, as function of the particle’s $\beta$, in different bins of transverse momentum, $p_{T}$, and for different radial decay positions of the particle. The efficiency is evaluated for particles which pass a loose set of fiducial requirements at generator level.

The parameterised efficiency for particles to pass full track selections in the metastable signal region, as function of the particle’s $\beta$, in different bins of transverse momentum, $p_{T}$, and for different radial decay positions of the particle. The efficiency is evaluated for particles which pass a loose set of fiducial requirements at generator level.

The parameterised efficiency for particles to pass full track selections in the metastable signal region, as function of the particle’s $\beta$, in different bins of transverse momentum, $p_{T}$, and for different radial decay positions of the particle. The efficiency is evaluated for particles which pass a loose set of fiducial requirements at generator level.

The parameterised efficiency for particles to pass full track selections in the metastable signal region, as function of the particle’s $\beta$, in different bins of transverse momentum, $p_{T}$, and for different radial decay positions of the particle. The efficiency is evaluated for particles which pass a loose set of fiducial requirements at generator level.

The parameterised efficiency for particles to pass full track selections in the metastable signal region, as function of the particle’s $\beta$, in different bins of transverse momentum, $p_{T}$. The stable efficiency is evaluated for samples which do not decay within the detector. The efficiency is evaluated for particles which pass a loose set of fiducial requirements at generator level.

The $m_{\mathrm{eff}}$ distribution for events passing the signal region requirements except the $m_{\mathrm{eff}}$ requirement in SR2. Distributions for data, the estimated SM backgrounds, and an example SUSY scenario are shown. "Other" is the sum of the $tWZ$, $t\bar{t}WW$, and $t\bar{t}t\bar{t}$ backgrounds. The last bin captures the overflow events. Both the statistical and systematic uncertainties in the SM background are included in the shaded band. The red arrows indicate the $m_{\mathrm{eff}}$ selections in the signal region.

Expected 95% CL exclusion limits on wino $W/Z$ NLSP pair production with RPV $\tilde{\chi}_1^0$ decays via $\lambda_{12k} \neq 0$ where $i,k \in{1,2}$. The limits are set using the statistical combination of disjoint signal regions. Where the signal regions are not mutually exclusive, the observed $\mathrm{CL}_s$ value is taken from the signal region with the better expected $\mathrm{CL}_s$ value.

Observed 95% CL exclusion limits on wino $W/Z$ NLSP pair production with RPV $\tilde{\chi}_1^0$ decays via $\lambda_{12k} \neq 0$ where $i,k \in{1,2}$. The limits are set using the statistical combination of disjoint signal regions. Where the signal regions are not mutually exclusive, the observed $\mathrm{CL}_s$ value is taken from the signal region with the better expected $\mathrm{CL}_s$ value.

Expected 95% CL exclusion limits on wino $W/Z$ NLSP pair production with RPV $\tilde{\chi}_1^0$ decays via $\lambda_{i33} \neq 0$ where $i,k \in{1,2}$. The limits are set using the statistical combination of disjoint signal regions. Where the signal regions are not mutually exclusive, the observed $\mathrm{CL}_s$ value is taken from the signal region with the better expected $\mathrm{CL}_s$ value.

Observed 95% CL exclusion limits on wino $W/Z$ NLSP pair production with RPV $\tilde{\chi}_1^0$ decays via $\lambda_{i33} \neq 0$ where $i,k \in{1,2}$. The limits are set using the statistical combination of disjoint signal regions. Where the signal regions are not mutually exclusive, the observed $\mathrm{CL}_s$ value is taken from the signal region with the better expected $\mathrm{CL}_s$ value.

Expected 95% CL exclusion limits on wino $W/h$ NLSP pair production with RPV $\tilde{\chi}_1^0$ decays via $\lambda_{12k} \neq 0$ where $i,k \in{1,2}$. The limits are set using the statistical combination of disjoint signal regions. Where the signal regions are not mutually exclusive, the observed $\mathrm{CL}_s$ value is taken from the signal region with the better expected $\mathrm{CL}_s$ value.

Observed 95% CL exclusion limits on wino $W/h$ NLSP pair production with RPV $\tilde{\chi}_1^0$ decays via $\lambda_{12k} \neq 0$ where $i,k \in{1,2}$. The limits are set using the statistical combination of disjoint signal regions. Where the signal regions are not mutually exclusive, the observed $\mathrm{CL}_s$ value is taken from the signal region with the better expected $\mathrm{CL}_s$ value.

Expected 95% CL exclusion limits on wino $W/h$ NLSP pair production with RPV $\tilde{\chi}_1^0$ decays via $\lambda_{i33} \neq 0$ where $i,k \in{1,2}$. The limits are set using the statistical combination of disjoint signal regions. Where the signal regions are not mutually exclusive, the observed $\mathrm{CL}_s$ value is taken from the signal region with the better expected $\mathrm{CL}_s$ value.

Observed 95% CL exclusion limits on wino $W/h$ NLSP pair production with RPV $\tilde{\chi}_1^0$ decays via $\lambda_{i33} \neq 0$ where $i,k \in{1,2}$. The limits are set using the statistical combination of disjoint signal regions. Where the signal regions are not mutually exclusive, the observed $\mathrm{CL}_s$ value is taken from the signal region with the better expected $\mathrm{CL}_s$ value.

Expected 95% CL exclusion limits on $\tilde{\ell}/\tilde{\nu}$ NLSP pair production with RPV $\tilde{\chi}_1^0$ decays via $\lambda_{12k} \neq 0$ where $i,k \in{1,2}$. The limits are set using the statistical combination of disjoint signal regions. Where the signal regions are not mutually exclusive, the observed $\mathrm{CL}_s$ value is taken from the signal region with the better expected $\mathrm{CL}_s$ value.

Observed 95% CL exclusion limits on $\tilde{\ell}/\tilde{\nu}$ NLSP pair production with RPV $\tilde{\chi}_1^0$ decays via $\lambda_{12k} \neq 0$ where $i,k \in{1,2}$. The limits are set using the statistical combination of disjoint signal regions. Where the signal regions are not mutually exclusive, the observed $\mathrm{CL}_s$ value is taken from the signal region with the better expected $\mathrm{CL}_s$ value.

Expected 95% CL exclusion limits on $\tilde{\ell}/\tilde{\nu}$ NLSP pair production with RPV $\tilde{\chi}_1^0$ decays via $\lambda_{i33} \neq 0$ where $i,k \in{1,2}$. The limits are set using the statistical combination of disjoint signal regions. Where the signal regions are not mutually exclusive, the observed $\mathrm{CL}_s$ value is taken from the signal region with the better expected $\mathrm{CL}_s$ value.

Observed 95% CL exclusion limits on $\tilde{\ell}/\tilde{\nu}$ NLSP pair production with RPV $\tilde{\chi}_1^0$ decays via $\lambda_{i33} \neq 0$ where $i,k \in{1,2}$. The limits are set using the statistical combination of disjoint signal regions. Where the signal regions are not mutually exclusive, the observed $\mathrm{CL}_s$ value is taken from the signal region with the better expected $\mathrm{CL}_s$ value.

Expected 95% CL exclusion limits on gluino NSLP pair production with RPV $\tilde{\chi}_1^0$ decays via $\lambda_{12k} \neq 0$ where $i,k \in{1,2}$. The limits are set using the statistical combination of disjoint signal regions. Where the signal regions are not mutually exclusive, the observed $\mathrm{CL}_s$ value is taken from the signal region with the better expected $\mathrm{CL}_s$ value.

Observed 95% CL exclusion limits on gluino NSLP pair production with RPV $\tilde{\chi}_1^0$ decays via $\lambda_{12k} \neq 0$ where $i,k \in{1,2}$. The limits are set using the statistical combination of disjoint signal regions. Where the signal regions are not mutually exclusive, the observed $\mathrm{CL}_s$ value is taken from the signal region with the better expected $\mathrm{CL}_s$ value.

Expected 95% CL exclusion limits on gluino NSLP pair production with RPV $\tilde{\chi}_1^0$ decays via $\lambda_{i33} \neq 0$ where $i,k \in{1,2}$. The limits are set using the statistical combination of disjoint signal regions. Where the signal regions are not mutually exclusive, the observed $\mathrm{CL}_s$ value is taken from the signal region with the better expected $\mathrm{CL}_s$ value.

Observed 95% CL exclusion limits on gluino NSLP pair production with RPV $\tilde{\chi}_1^0$ decays via $\lambda_{i33} \neq 0$ where $i,k \in{1,2}$. The limits are set using the statistical combination of disjoint signal regions. Where the signal regions are not mutually exclusive, the observed $\mathrm{CL}_s$ value is taken from the signal region with the better expected $\mathrm{CL}_s$ value.

Expected 95% CL exclusion limits on the higgsino GGM models. The limits are set using the statistical combination of disjoint signal regions. Where the signal regions are not mutually exclusive, the observed $\mathrm{CL}_s$ value is taken from the signal region with the better expected $\mathrm{CL}_s$ value.

Observed 95% CL exclusion limits on the higgsino GGM models. The limits are set using the statistical combination of disjoint signal regions. Where the signal regions are not mutually exclusive, the observed $\mathrm{CL}_s$ value is taken from the signal region with the better expected $\mathrm{CL}_s$ value.

Observed exclusion limits on the $\tilde{\chi}_1^{\pm}/\tilde{\chi}_2^0$ masses in the context of the wino $W/Z$ RPV scenario with $\lambda_{12k} \neq 0$ couplings. The signal region used to obtain the limits is SR0A. All limits are computed at 95% CL. The grey numbers represent the 95% CL upper limits on the production cross-section (in fb) obtained using the signal efficiency and acceptance specific to this model.

Observed exclusion limits on the $\tilde{\chi}_1^{\pm}/\tilde{\chi}_2^0$ masses in the context of the wino $W/Z$ RPV scenario with $\lambda_{i33} \neq 0$ couplings. The signal region used to obtain the limits is the combination of SR0A, SR1 and SR2. All limits are computed at 95% CL. The grey numbers represent the 95% CL upper limits on the production cross-section (in fb) obtained using the signal efficiency and acceptance specific to this model.

Observed exclusion limits on the $\tilde{\chi}_1^{\pm}/\tilde{\chi}_2^0$ masses in the context of the wino $W/Z$ RPV scenario with $\lambda_{12k} \neq 0$ couplings. The signal region used to obtain the limits is SR0B. All limits are computed at 95% CL. The grey numbers represent the 95% CL upper limits on the production cross-section (in fb) obtained using the signal efficiency and acceptance specific to this model.

Observed exclusion limits on the $\tilde{\chi}_1^{\pm}/\tilde{\chi}_2^0$ masses in the context of the wino $W/Z$ RPV scenario with $\lambda_{i33} \neq 0$ couplings. The signal region used to obtain the limits is the combination of SR0B, SR1 and SR2. All limits are computed at 95% CL. The grey numbers represent the 95% CL upper limits on the production cross-section (in fb) obtained using the signal efficiency and acceptance specific to this model.

Observed exclusion limits on the $\tilde{\chi}_1^{\pm}/\tilde{\chi}_2^0$ masses in the context of the wino $W/h$ RPV scenario with $\lambda_{12k} \neq 0$ couplings. The signal region used to obtain the limits is SR0A. All limits are computed at 95% CL. The grey numbers represent the 95% CL upper limits on the production cross-section (in fb) obtained using the signal efficiency and acceptance specific to this model.

Observed exclusion limits on the $\tilde{\chi}_1^{\pm}/\tilde{\chi}_2^0$ masses in the context of the wino $W/h$ RPV scenario with $\lambda_{i33} \neq 0$ couplings. The signal region used to obtain the limits is the combination of SR0A, SR1 and SR2. All limits are computed at 95% CL. The grey numbers represent the 95% CL upper limits on the production cross-section (in fb) obtained using the signal efficiency and acceptance specific to this model.

Observed exclusion limits on the $\tilde{\chi}_1^{\pm}/\tilde{\chi}_2^0$ masses in the context of the wino $W/h$ RPV scenario with $\lambda_{12k} \neq 0$ couplings. The signal region used to obtain the limits is SR0B. All limits are computed at 95% CL. The grey numbers represent the 95% CL upper limits on the production cross-section (in fb) obtained using the signal efficiency and acceptance specific to this model.

Observed exclusion limits on the $\tilde{\chi}_1^{\pm}/\tilde{\chi}_2^0$ masses in the context of the wino $W/h$ RPV scenario with $\lambda_{i33} \neq 0$ couplings. The signal region used to obtain the limits is the combination of SR0B, SR1 and SR2. All limits are computed at 95% CL. The grey numbers represent the 95% CL upper limits on the production cross-section (in fb) obtained using the signal efficiency and acceptance specific to this model.

Observed exclusion limits on the slepton/sneutrino masses in the context of the slepton/sneutrino RPV scenario with $\lambda_{12k} \neq 0$ couplings. The signal region used to obtain the limits is SR0A. All limits are computed at 95% CL. The grey numbers represent the 95% CL upper limits on the production cross-section (in fb) obtained using the signal efficiency and acceptance specific to this model.

Observed exclusion limits on the slepton/sneutrino masses in the context of the slepton/sneutrino RPV scenario with $\lambda_{i33} \neq 0$ couplings. The signal region used to obtain the limits is the combination of SR0A, SR1 and SR2. All limits are computed at 95% CL. The grey numbers represent the 95% CL upper limits on the production cross-section (in fb) obtained using the signal efficiency and acceptance specific to this model.

Observed exclusion limits on the slepton/sneutrino masses in the context of the slepton/sneutrino RPV scenario with $\lambda_{12k} \neq 0$ couplings. The signal region used to obtain the limits is SR0B. All limits are computed at 95% CL. The grey numbers represent the 95% CL upper limits on the production cross-section (in fb) obtained using the signal efficiency and acceptance specific to this model.

Observed exclusion limits on the slepton/sneutrino masses in the context of the slepton/sneutrino RPV scenario with $\lambda_{i33} \neq 0$ couplings. The signal region used to obtain the limits is the combination of SR0B, SR1 and SR2. All limits are computed at 95% CL. The grey numbers represent the 95% CL upper limits on the production cross-section (in fb) obtained using the signal efficiency and acceptance specific to this model.

Observed exclusion limits on the gluino masses in the context of the gluino RPV scenario with $\lambda_{12k} \neq 0$ couplings. The signal region used to obtain the limits is SR0A. All limits are computed at 95% CL. The grey numbers represent the 95% CL upper limits on the production cross-section (in fb) obtained using the signal efficiency and acceptance specific to this model.

Observed exclusion limits on the gluino masses in the context of the gluino RPV scenario with $\lambda_{i33} \neq 0$ couplings. The signal region used to obtain the limits is the combination of SR0A, SR1 and SR2. All limits are computed at 95% CL. The grey numbers represent the 95% CL upper limits on the production cross-section (in fb) obtained using the signal efficiency and acceptance specific to this model.

Observed exclusion limits on the gluino masses in the context of the gluino RPV scenario with $\lambda_{12k} \neq 0$ couplings. The signal region used to obtain the limits is SR0B. All limits are computed at 95% CL. The grey numbers represent the 95% CL upper limits on the production cross-section (in fb) obtained using the signal efficiency and acceptance specific to this model.

Observed exclusion limits on the gluino masses in the context of the gluino RPV scenario with $\lambda_{i33} \neq 0$ couplings. The signal region used to obtain the limits is the combination of SR0B, SR1 and SR2. All limits are computed at 95% CL. The grey numbers represent the 95% CL upper limits on the production cross-section (in fb) obtained using the signal efficiency and acceptance specific to this model.

Observed exclusion limits on the $\tilde{\chi}_1^{\pm}\tilde{\chi}_2^{0}\tilde{\chi}_1^{0}$ masses in the context of the higgsino GGM scenario in SR0C. All limits are computed at 95% CL. The grey numbers represent the 95% CL upper limits on the production cross-section (in fb) obtained using the signal efficiency and acceptance specific to this model.

Observed exclusion limits on the $\tilde{\chi}_1^{\pm}\tilde{\chi}_2^{0}\tilde{\chi}_1^{0}$ masses in the context of the higgsino GGM scenario in SR0D. All limits are computed at 95% CL. The grey numbers represent the 95% CL upper limits on the production cross-section (in fb) obtained using the signal efficiency and acceptance specific to this model.

The best expected exclusion power between signal regions at each signal point as is adopted in the limit calculation of the wino $W/Z$ model for $\lambda_{12k} \neq 0$ RPV couplings. For the $\lambda_{12k} \neq 0$ case, the results from SR0B are adopted everywhere in the final exclusion limit contours, as is found to be the most powerful signal region among SR0A and SR0B in the majority of the signal grid points of this model.

The best expected exclusion power between signal regions at each signal point as is adopted in the limit calculation of the wino $W/Z$ model for $\lambda_{i33} \neq 0$ RPV couplings. A combination of SR0A + SR1 + SR2 or SR0B + SR1 + SR2 is tested to detect the best expected limit; shown on the plot is the region (SR0A or SR0B) which is chosen to provide the best expected limit in combination with SR1 and SR2. The results from the combination SR0B + SR1 + SR2 are finally adopted everywhere in the final exclusion limit contour since they provide the best expected limit.

The best expected exclusion power between signal regions at each signal point as is adopted in the limit calculation of the wino $W/h$ model for $\lambda_{12k} \neq 0$ RPV couplings. For the $\lambda_{12k} \neq 0$ case, the results from SR0B are adopted everywhere in the final exclusion limit contours, as is found to be the most powerful signal region among SR0A and SR0B in the majority of the signal grid points of this model.

The best expected exclusion power between signal regions at each signal point as is adopted in the limit calculation of the wino $W/h$ model for $\lambda_{i33} \neq 0$ RPV couplings. A combination of SR0A + SR1 + SR2 or SR0B + SR1 + SR2 is tested to detect the best expected limit; shown on the plot is the region (SR0A or SR0B) which is chosen to provide the best expected limit in combination with SR1 and SR2. The results from the combination SR0B + SR1 + SR2 are finally adopted everywhere in the final exclusion limit contour since they provide the best expected limit.

The best expected exclusion power between signal regions at each signal point as is adopted in the limit calculation of the slepton/sneutrino model for $\lambda_{12k} \neq 0$ RPV couplings.

The best expected exclusion power between signal regions at each signal point as is adopted in the limit calculation of the slepton/sneutrino model for $\lambda_{i33} \neq 0$ RPV couplings. A combination of SR0A + SR1 + SR2 or SR0B + SR1 + SR2 is tested to detect the best expected limit; shown on the plot is the region (SR0A or SR0B) which is chosen to provide the best expected limit in combination with SR1 and SR2.

The best expected exclusion power between signal regions at each signal point as is adopted in the limit calculation of the gluino model for $\lambda_{12k} \neq 0$ RPV couplings.

The best expected exclusion power between signal regions at each signal point as is adopted in the limit calculation of the gluino model for $\lambda_{i33} \neq 0$ RPV couplings. A combination of SR0A + SR1 + SR2 or SR0B + SR1 + SR2 is tested to detect the best expected limit; shown on the plot is the region (SR0A or SR0B) which is chosen to provide the best expected limit in combination with SR1 and SR2.

The best expected exclusion power between the overlapping SR0C and SR0D at each signal point as is adopted in the limit combination of the GGM higgsino model.

Acceptance ($A$) at particle level and selection efficiency ($\epsilon$) at reconstruction level for the Wino $\tilde{\chi}_1^{+}\tilde{\chi}_1^{-}$ $\lambda_{12k} \neq 0$ model fulfilling the selection criteria of SR0A. The implementation of "simple" cut selection at truth level uses the SimpleAnalysis framework and the code can be found in the Resources of the HepData entry for this paper.

Acceptance ($A$) at particle level and selection efficiency ($\epsilon$) at reconstruction level for the Wino $\tilde{\chi}_1^{\pm}\tilde{\chi}_2^{0}$ $W/Z$ $\lambda_{12k} \neq 0$ model fulfilling the selection criteria of SR0A. The implementation of "simple" cut selection at truth level uses the SimpleAnalysis framework and the code can be found in the Resources of the HepData entry for this paper.

Acceptance ($A$) at particle level and selection efficiency ($\epsilon$) at reconstruction level for the Wino $\tilde{\chi}_1^{\pm}\tilde{\chi}_2^{0}$ $W/h$ $\lambda_{12k} \neq 0$ model fulfilling the selection criteria of SR0A. The implementation of "simple" cut selection at truth level uses the SimpleAnalysis framework and the code can be found in the Resources of the HepData entry for this paper.

Acceptance ($A$) at particle level and selection efficiency ($\epsilon$) at reconstruction level for the slepton/sneutrino $\lambda_{12k} \neq 0$ model fulfilling the selection criteria of SR0A. The implementation of "simple" cut selection at truth level uses the SimpleAnalysis framework and the code can be found in the Resources of the HepData entry for this paper.

Acceptance ($A$) at particle level and selection efficiency ($\epsilon$) at reconstruction level for the gluino $\lambda_{12k} \neq 0$ model fulfilling the selection criteria of SR0A. The implementation of "simple" cut selection at truth level uses the SimpleAnalysis framework and the code can be found in the Resources of the HepData entry for this paper.

Acceptance ($A$) at particle level and selection efficiency ($\epsilon$) at reconstruction level for the Wino $\tilde{\chi}_1^{+}\tilde{\chi}_1^{-}$ $\lambda_{12k} \neq 0$ model fulfilling the selection criteria of SR0B. The implementation of "simple" cut selection at truth level uses the SimpleAnalysis framework and the code can be found in the Resources of the HepData entry for this paper.

Acceptance ($A$) at particle level and selection efficiency ($\epsilon$) at reconstruction level for the Wino $\tilde{\chi}_1^{\pm}\tilde{\chi}_2^{0}$ $W/Z$ $\lambda_{12k} \neq 0$ model fulfilling the selection criteria of SR0B. The implementation of "simple" cut selection at truth level uses the SimpleAnalysis framework and the code can be found in the Resources of the HepData entry for this paper.

Acceptance ($A$) at particle level and selection efficiency ($\epsilon$) at reconstruction level for the Wino $\tilde{\chi}_1^{\pm}\tilde{\chi}_2^{0}$ $W/h$ $\lambda_{12k} \neq 0$ model fulfilling the selection criteria of SR0B. The implementation of "simple" cut selection at truth level uses the SimpleAnalysis framework and the code can be found in the Resources of the HepData entry for this paper.

Acceptance ($A$) at particle level and selection efficiency ($\epsilon$) at reconstruction level for the slepton/sneutrino $\lambda_{12k} \neq 0$ model fulfilling the selection criteria of SR0B. The implementation of "simple" cut selection at truth level uses the SimpleAnalysis framework and the code can be found in the Resources of the HepData entry for this paper.

Acceptance ($A$) at particle level and selection efficiency ($\epsilon$) at reconstruction level for the gluino $\lambda_{12k} \neq 0$ model fulfilling the selection criteria of SR0B. The implementation of "simple" cut selection at truth level uses the SimpleAnalysis framework and the code can be found in the Resources of the HepData entry for this paper.

Acceptance ($A$) at particle level and selection efficiency ($\epsilon$) at reconstruction level for the GGM Higgsino models with BR($\tilde{\chi}_1^0 \rightarrow Z \tilde{G}$)=100% and BR($\tilde{\chi}_1^0 \rightarrow Z \tilde{G}$)=50% fulfilling the selection criteria of SR0C. The implementation of "simple" cut selection at truth level uses the SimpleAnalysis framework and the code can be found in the Resources of the HepData entry for this paper.

Acceptance ($A$) at particle level and selection efficiency ($\epsilon$) at reconstruction level for the GGM Higgsino models with BR($\tilde{\chi}_1^0 \rightarrow Z \tilde{G}$)=100% and BR($\tilde{\chi}_1^0 \rightarrow Z \tilde{G}$)=50% fulfilling the selection criteria of SR0D. The implementation of "simple" cut selection at truth level uses the SimpleAnalysis framework and the code can be found in the Resources of the HepData entry for this paper.

Acceptance ($A$) at particle level and selection efficiency ($\epsilon$) at reconstruction level for the Wino $\tilde{\chi}_1^{+}\tilde{\chi}_1^{-}$ $\lambda_{i33} \neq 0$ model fulfilling the selection criteria of SR1. The implementation of "simple" cut selection at truth level uses the SimpleAnalysis framework and the code can be found in the Resources of the HepData entry for this paper.

Acceptance ($A$) at particle level and selection efficiency ($\epsilon$) at reconstruction level for the Wino $\tilde{\chi}_1^{\pm}\tilde{\chi}_2^{0}$ $W/Z$ $\lambda_{i33} \neq 0$ model fulfilling the selection criteria of SR1. The implementation of "simple" cut selection at truth level uses the SimpleAnalysis framework and the code can be found in the Resources of the HepData entry for this paper.

Acceptance ($A$) at particle level and selection efficiency ($\epsilon$) at reconstruction level for the Wino $\tilde{\chi}_1^{\pm}\tilde{\chi}_2^{0}$ $W/h$ $\lambda_{i33} \neq 0$ model fulfilling the selection criteria of SR1. The implementation of "simple" cut selection at truth level uses the SimpleAnalysis framework and the code can be found in the Resources of the HepData entry for this paper.

Acceptance ($A$) at particle level and selection efficiency ($\epsilon$) at reconstruction level for the slepton/sneutrino $\lambda_{i33} \neq 0$ model fulfilling the selection criteria of SR1. The implementation of "simple" cut selection at truth level uses the SimpleAnalysis framework and the code can be found in the Resources of the HepData entry for this paper.

Acceptance ($A$) at particle level and selection efficiency ($\epsilon$) at reconstruction level for the gluino $\lambda_{i33} \neq 0$ model fulfilling the selection criteria of SR1. The implementation of "simple" cut selection at truth level uses the SimpleAnalysis framework and the code can be found in the Resources of the HepData entry for this paper.

Acceptance ($A$) at particle level and selection efficiency ($\epsilon$) at reconstruction level for the Wino $\tilde{\chi}_1^{+}\tilde{\chi}_1^{-}$ $\lambda_{i33} \neq 0$ model fulfilling the selection criteria of SR2. The implementation of "simple" cut selection at truth level uses the SimpleAnalysis framework and the code can be found in the Resources of the HepData entry for this paper.

Acceptance ($A$) at particle level and selection efficiency ($\epsilon$) at reconstruction level for the Wino $\tilde{\chi}_1^{\pm}\tilde{\chi}_2^{0}$ $W/Z$ $\lambda_{i33} \neq 0$ model fulfilling the selection criteria of SR2. The implementation of "simple" cut selection at truth level uses the SimpleAnalysis framework and the code can be found in the Resources of the HepData entry for this paper.

Acceptance ($A$) at particle level and selection efficiency ($\epsilon$) at reconstruction level for the Wino $\tilde{\chi}_1^{\pm}\tilde{\chi}_2^{0}$ $W/h$ $\lambda_{i33} \neq 0$ model fulfilling the selection criteria of SR2. The implementation of "simple" cut selection at truth level uses the SimpleAnalysis framework and the code can be found in the Resources of the HepData entry for this paper.

Acceptance ($A$) at particle level and selection efficiency ($\epsilon$) at reconstruction level for the slepton/sneutrino $\lambda_{i33} \neq 0$ model fulfilling the selection criteria of SR2. The implementation of "simple" cut selection at truth level uses the SimpleAnalysis framework and the code can be found in the Resources of the HepData entry for this paper.

Acceptance ($A$) at particle level and selection efficiency ($\epsilon$) at reconstruction level for the gluino $\lambda_{i33} \neq 0$ model fulfilling the selection criteria of SR2. The implementation of "simple" cut selection at truth level uses the SimpleAnalysis framework and the code can be found in the Resources of the HepData entry for this paper.

Cutflow event yields in regions SR0A and SR0B for RPV models with the $\lambda_{12k} \neq 0$ coupling. All yields correspond to weighted events, so that effects from lepton reconstruction efficiencies, trigger corrections, pileup reweighting, etc., are included. They are normalized to the integrated luminosity of the data sample, $\int L dt = 36.1$ fb$^{-1}$. The "Initial" number indicates the weighted number of events before any selection cut is applied. The last entries show the normalized number of events surviving the selection requirements of SR0A and SR0B.

Cutflow event yields in regions SR0C and SR0D for GGM Higgsino models with BR($\tilde{\chi}_1^0 \rightarrow Z \tilde{G}$)=100% or BR($\tilde{\chi}_1^0 \rightarrow Z \tilde{G}$)=50%, and $\tilde{\chi}_1^{\pm}\tilde{\chi}_2^0\tilde{\chi}_1^0$ mass of 400 GeV. All yields correspond to weighted events, so that effects from lepton reconstruction efficiencies, trigger corrections, pileup reweighting, etc., are included. They are normalized to the integrated luminosity of the data sample, $\int L dt = 36.1$ fb$^{-1}$. The "Initial" number indicates the weighted number of events before any selection cut is applied. The "Generator Filter" step is applied during the MC generation of the simulated events; the BR=100% sample has a generator filter of $\geq 4e/\mu$ leptons with $p_{\mathrm{T}}>4$ GeV and $|\eta|<2.8$, and the BR=50% sample has a generator filter of $\geq 4 e/\mu/\tau_{\mathrm{had-vis}}$ leptons with $p_{\mathrm{T}}(e,\mu)>4$ GeV, $p_{\mathrm{T}}(\tau_{\mathrm{had-vis}}^{\mathrm{visible}})>15$ GeV and $|\eta|<2.8$. The $ZZ$ selection cutflow step refers to the mass window cut for the leading and subleading $Z$ boson candidate between $81.2-101.2$ GeV and $61.2-101.2$ GeV, respectively. The last entries show the efficiency of events surviving the selection requirements defined in SR0C and SR0D.

Cutflow event yields in region SR1 for RPV models with the $\lambda_{i33} \neq 0$ coupling. All yields correspond to weighted events, so that effects from lepton reconstruction efficiencies, trigger corrections, pileup reweighting, etc., are included. They are normalized to the integrated luminosity of the data sample, $\int L dt = 36.1$ fb$^{-1}$. The "Initial" number indicates the weighted number of events before any selection cut is applied. The last entries show the normalized number of events surviving the selection requirements of SR1.

Cutflow event yields in region SR2 for RPV models with the $\lambda_{i33} \neq 0$ coupling. All yields correspond to weighted events, so that effects from lepton reconstruction efficiencies, trigger corrections, pileup reweighting, etc., are included. They are normalized to the integrated luminosity of the data sample, $\int L dt = 36.1$ fb$^{-1}$. The "Initial" number indicates the weighted number of events before any selection cut is applied. The last entries show the normalized number of events surviving the selection requirements of SR2.

Cross sections for the slepton/snuetrino model for different NLSP masses.

Scan of the likelihood function for the search in the $\tau\tau$ final state for a single narrow resonance, $\phi$, produced via gluon fusion ($gg\phi$) or in association with b quarks ($bb\phi$). The scan is performed in 40000 points of the ($\sigma(gg\phi)\cdot B(\phi\rightarrow\tau\tau)$, $\sigma(bb\phi)\cdot B(\phi\rightarrow\tau\tau)$) plane. An asimov dataset constructed from the expectation of all backgrounds is tested against a background hypothesis without the SM Higgs boson. For further details and instructions, please have a look into the following README file http://cms-results.web.cern.ch/cms-results/public-results/publications/HIG-17-020/2D-likelihood-scans/README.txt. Selected examples of such a likelihood scan are given in Figure 8 of the paper.

Scan of the likelihood function for the search in the $\tau\tau$ final state for a single narrow resonance, $\phi$, produced via gluon fusion ($gg\phi$) or in association with b quarks ($bb\phi$). The scan is performed in 40000 points of the ($\sigma(gg\phi)\cdot B(\phi\rightarrow\tau\tau)$, $\sigma(bb\phi)\cdot B(\phi\rightarrow\tau\tau)$) plane. The observation in data is tested against a background hypothesis including the SM Higgs boson. For further details and instructions, please have a look into the following README file http://cms-results.web.cern.ch/cms-results/public-results/publications/HIG-17-020/2D-likelihood-scans/README.txt. Selected examples of such a likelihood scan are given in Figure 8 of the paper.

Scan of the likelihood function for the search in the $\tau\tau$ final state for a single narrow resonance, $\phi$, produced via gluon fusion ($gg\phi$) or in association with b quarks ($bb\phi$). The scan is performed in 40000 points of the ($\sigma(gg\phi)\cdot B(\phi\rightarrow\tau\tau)$, $\sigma(bb\phi)\cdot B(\phi\rightarrow\tau\tau)$) plane. The observation in data is tested against a background hypothesis without the SM Higgs boson. For further details and instructions, please have a look into the following README file http://cms-results.web.cern.ch/cms-results/public-results/publications/HIG-17-020/2D-likelihood-scans/README.txt. Selected examples of such a likelihood scan are given in Figure 8 of the paper.