Distributions of $s_{\text{T}}$ in the single-tau one-bin SR. The stacked histograms show the various SM background contributions. The hatched band indicates the total statistical and systematic uncertainty of the SM background. The $t\bar{t}$ (2 real $\tau$) and $t\bar{t}$ (1 real $\tau$) as well as the single-top background contributions are scaled with the normalization factors obtained from the background-only fit. Minor backgrounds are grouped together and denoted as 'Other'. This includes $t\bar{t}$-fake, single top, and other top (di-tau channel) or $t\bar{t}$-fake, $t\bar{t}+H$, multiboson, and other top (single-tau channel). The overlaid dotted lines show the additional contributions for signal scenarios close to the expected exclusion contour with the particle type and the mass and $\beta$ parameters for the simplified models indicated in the legend. For the leptoquark signal model the shapes of the distributions for $\text{LQ}_{3}^{\text{d}}$ and $\text{LQ}_{3}^{\text{v}}$ (not shown) are similar to that of $\text{LQ}_{3}^{\text{u}}$. The rightmost bin includes the overflow.

Distributions of $m_{\text{T}}(\tau)$ in the single-tau one-bin SR. The stacked histograms show the various SM background contributions. The hatched band indicates the total statistical and systematic uncertainty of the SM background. The $t\bar{t}$ (2 real $\tau$) and $t\bar{t}$ (1 real $\tau$) as well as the single-top background contributions are scaled with the normalization factors obtained from the background-only fit. Minor backgrounds are grouped together and denoted as 'Other'. This includes $t\bar{t}$-fake, single top, and other top (di-tau channel) or $t\bar{t}$-fake, $t\bar{t}+H$, multiboson, and other top (single-tau channel). The overlaid dotted lines show the additional contributions for signal scenarios close to the expected exclusion contour with the particle type and the mass and $\beta$ parameters for the simplified models indicated in the legend. For the leptoquark signal model the shapes of the distributions for $\text{LQ}_{3}^{\text{d}}$ and $\text{LQ}_{3}^{\text{v}}$ (not shown) are similar to that of $\text{LQ}_{3}^{\text{u}}$. The rightmost bin includes the overflow.

Distributions of $\Sigma m_{\text{T}}(b_{1,2})$ in the single-tau $p_{\text{T}}(\tau)$-binned SR. The stacked histograms show the various SM background contributions. The hatched band indicates the total statistical and systematic uncertainty of the SM background. The $t\bar{t}$ (2 real $\tau$) and $t\bar{t}$ (1 real $\tau$) as well as the single-top background contributions are scaled with the normalization factors obtained from the background-only fit. Minor backgrounds are grouped together and denoted as 'Other'. This includes $t\bar{t}$-fake, single top, and other top (di-tau channel) or $t\bar{t}$-fake, $t\bar{t}+H$, multiboson, and other top (single-tau channel). The overlaid dotted lines show the additional contributions for signal scenarios close to the expected exclusion contour with the particle type and the mass and $\beta$ parameters for the simplified models indicated in the legend. For the leptoquark signal model the shapes of the distributions for $\text{LQ}_{3}^{\text{d}}$ and $\text{LQ}_{3}^{\text{v}}$ (not shown) are similar to that of $\text{LQ}_{3}^{\text{u}}$. The rightmost bin includes the overflow.

Distributions of $p_{\text{T}}(\tau)$ in the single-tau $p_{\text{T}}(\tau)$-binned SR. The stacked histograms show the various SM background contributions. The hatched band indicates the total statistical and systematic uncertainty of the SM background. The $t\bar{t}$ (2 real $\tau$) and $t\bar{t}$ (1 real $\tau$) as well as the single-top background contributions are scaled with the normalization factors obtained from the background-only fit. Minor backgrounds are grouped together and denoted as 'Other'. This includes $t\bar{t}$-fake, single top, and other top (di-tau channel) or $t\bar{t}$-fake, $t\bar{t}+H$, multiboson, and other top (single-tau channel). The overlaid dotted lines show the additional contributions for signal scenarios close to the expected exclusion contour with the particle type and the mass and $\beta$ parameters for the simplified models indicated in the legend. For the leptoquark signal model the shapes of the distributions for $\text{LQ}_{3}^{\text{d}}$ and $\text{LQ}_{3}^{\text{v}}$ (not shown) are similar to that of $\text{LQ}_{3}^{\text{u}}$. The rightmost bin includes the overflow.

Observed event yields in data ('Observed') and expected event yields for SM background processes obtained from the background-only fit ('Total bkg.' and rows below) in the signal regions of the di-tau and single-tau channels. The quoted uncertainties include both the statistical and systematic uncertainties and are truncated at zero yield. By construction, no $t\bar{t}$ (2 real $\tau$) events can pass the selections in the single-tau channel. As the individual uncertainties are correlated, they do not add in quadrature to equal the total background uncertainty.

From left to right: upper limits at the 95% confidence level (CL) on the visible cross section ($\sigma_\text{vis}$) and on the number of signal events ($S_{\text{obs}}^{95}$). The third column ($S_{\text{exp}}^{95}$) shows the upper limit at the 95% CL on the number of signal events, given the expected number (and $\pm 1\,\sigma$ excursions on the expectation) of background events. The last two columns indicate the confidence level observed for the background-only hypothesis ($\text{CL}_{b}$), the discovery $p$-value ($p(s=0)$) and the significance ($Z$). In the di-tau SR, where fewer events are observed than predicted by the fitted background estimate, the $p$-value is capped at 0.5.

Expected and observed exclusion contours at the 95% confidence level for the vector third-generation leptoquark signal model, as a function of the mass $m(\text{LQ}_{3}^{\text{v}})$ and the branching fraction $B(\text{LQ}_{3}^{\text{v}} \rightarrow b\tau)$ into a quark and a charged lepton. The plot shows the exclusion contour for the minimal-coupling scenario. The limits are derived from the binned single-tau signal region.

Expected and observed exclusion contours at the 95% confidence level for the vector third-generation leptoquark signal model, as a function of the mass $m(\text{LQ}_{3}^{\text{v}})$ and the branching fraction $B(\text{LQ}_{3}^{\text{v}} \rightarrow b\tau)$ into a quark and a charged lepton. The plot shows the exclusion contour for the minimal-coupling scenario. The limits are derived from the binned single-tau signal region.

Expected and observed exclusion contours at the 95% confidence level for the vector third-generation leptoquark signal model, as a function of the mass $m(\text{LQ}_{3}^{\text{v}})$ and the branching fraction $B(\text{LQ}_{3}^{\text{v}} \rightarrow b\tau)$ into a quark and a charged lepton. The plot shows the exclusion contour for vector leptoquarks with additional gauge couplings. The limits are derived from the binned single-tau signal region.

Expected and observed exclusion contours at the 95% confidence level for the vector third-generation leptoquark signal model, as a function of the mass $m(\text{LQ}_{3}^{\text{v}})$ and the branching fraction $B(\text{LQ}_{3}^{\text{v}} \rightarrow b\tau)$ into a quark and a charged lepton. The plot shows the exclusion contour for vector leptoquarks with additional gauge couplings. The limits are derived from the binned single-tau signal region.

Exclusion contours at the 95% confidence level for the stop-stau signal model as a function of the masses of the top squark $m(\tilde{t}_{1})$ and of the tau slepton $m(\tilde{\tau}_{1})$. Expected and observed limits are shown for the present search in comparison to observed limits from previous ATLAS analyses based on data from Run-1 of the LHC at $\sqrt{s} = 8$ TeV [Eur. Phys. J. C 76 (2016)] and on a partial dataset from Run 2 at $\sqrt{s} = 13$ TeV [Phys. Rev. D 98 (2018) 032008]. The green band indicates the limit on the mass of the tau slepton (for a massless LSP) from the LEP experiments.

Exclusion contours at the 95% confidence level for the stop-stau signal model as a function of the masses of the top squark $m(\tilde{t}_{1})$ and of the tau slepton $m(\tilde{\tau}_{1})$. Expected and observed limits are shown for the present search in comparison to observed limits from previous ATLAS analyses based on data from Run-1 of the LHC at $\sqrt{s} = 8$ TeV [Eur. Phys. J. C 76 (2016)] and on a partial dataset from Run 2 at $\sqrt{s} = 13$ TeV [Phys. Rev. D 98 (2018) 032008]. The green band indicates the limit on the mass of the tau slepton (for a massless LSP) from the LEP experiments.

Expected and observed exclusion contours at the 95% confidence level for the scalar third-generation leptoquark signal model, as a function of the mass $m(\text{LQ}_{3}^{\text{u}})$ and the branching fraction $B(\text{LQ}_{3}^{\text{u}} \rightarrow q\ell)$ into a quark and a charged lepton. The plot shows the exclusion contour for up-type leptoquarks $\text{LQ}_{3}^{\text{u}})$ with charge $+2/3e$. The limits are derived from the binned single-tau signal region. Shown in gray for comparison are the observed exclusion-limit contours from the previous ATLAS publication that targets the same leptoquark models but is based on a subset of the Run-2 data [JHEP 06 (2019) 144]. In this previous publication five different analyses are considered that target not only the final state studied here but also the final states that correspond to a branching fraction $B(\text{LQ}_{3}^{\text{u}} \rightarrow q\ell)$ of 0 or 1, leading to the concave shapes of the gray exclusion contours.

Expected and observed exclusion contours at the 95% confidence level for the scalar third-generation leptoquark signal model, as a function of the mass $m(\text{LQ}_{3}^{\text{u}})$ and the branching fraction $B(\text{LQ}_{3}^{\text{u}} \rightarrow q\ell)$ into a quark and a charged lepton. The plot shows the exclusion contour for up-type leptoquarks $\text{LQ}_{3}^{\text{u}})$ with charge $+2/3e$. The limits are derived from the binned single-tau signal region. Shown in gray for comparison are the observed exclusion-limit contours from the previous ATLAS publication that targets the same leptoquark models but is based on a subset of the Run-2 data [JHEP 06 (2019) 144]. In this previous publication five different analyses are considered that target not only the final state studied here but also the final states that correspond to a branching fraction $B(\text{LQ}_{3}^{\text{u}} \rightarrow q\ell)$ of 0 or 1, leading to the concave shapes of the gray exclusion contours.

Expected and observed exclusion contours at the 95% confidence level for the scalar third-generation leptoquark signal model, as a function of the mass $m(\text{LQ}_{3}^{\text{d}})$ and the branching fraction $B(\text{LQ}_{3}^{\text{d}} \rightarrow q\ell)$ into a quark and a charged lepton. The plot shows the exclusion contour for down-type leptoquarks $\text{LQ}_{3}^{\text{d}})$ with charge $-1/3e$. The limits are derived from the binned single-tau signal region. Shown in gray for comparison are the observed exclusion-limit contours from the previous ATLAS publication that targets the same leptoquark models but is based on a subset of the Run-2 data [JHEP 06 (2019) 144]. In this previous publication five different analyses are considered that target not only the final state studied here but also the final states that correspond to a branching fraction $B(\text{LQ}_{3}^{\text{d}} \rightarrow q\ell)$ of 0 or 1, leading to the concave shapes of the gray exclusion contours.

Expected and observed exclusion contours at the 95% confidence level for the scalar third-generation leptoquark signal model, as a function of the mass $m(\text{LQ}_{3}^{\text{d}})$ and the branching fraction $B(\text{LQ}_{3}^{\text{d}} \rightarrow q\ell)$ into a quark and a charged lepton. The plot shows the exclusion contour for down-type leptoquarks $\text{LQ}_{3}^{\text{d}})$ with charge $-1/3e$. The limits are derived from the binned single-tau signal region. Shown in gray for comparison are the observed exclusion-limit contours from the previous ATLAS publication that targets the same leptoquark models but is based on a subset of the Run-2 data [JHEP 06 (2019) 144]. In this previous publication five different analyses are considered that target not only the final state studied here but also the final states that correspond to a branching fraction $B(\text{LQ}_{3}^{\text{d}} \rightarrow q\ell)$ of 0 or 1, leading to the concave shapes of the gray exclusion contours.

Upper limits on the signal cross section at the 95 % confidence level for the stop-stau signal model.

Upper limits on the signal cross section at the 95 % confidence level for the scalar third-generation leptoquark signal model with up-type leptoquarks.

Upper limits on the signal cross section at the 95 % confidence level for the scalar third-generation leptoquark signal model with down-type leptoquarks.

Upper limits on the signal cross section at the 95 % confidence level for the vector third-generation leptoquark signal model with minimal coupling (MC).

Upper limits on the signal cross section at the 95 % confidence level for the vector third-generation leptoquark signal model with additional gauge couplings (YM).

Acceptance of the one-bin signal region of the single-tau channel for pair production of up-type leptoquarks $\text{LQ}_{3}^{\text{u}}$.

Efficiency of the one-bin signal region of the single-tau channel for pair production of up-type leptoquarks $\text{LQ}_{3}^{\text{u}}$. The plot does not show efficiencies for a branching fraction $B(\text{LQ}_{3}^{\text{u}} \rightarrow b\tau)$ of 0 or 1 because here the acceptance at generator level becomes zero and the efficiency is thus undefined.

Acceptance of the first bin of the multi-bin signal region (50 GeV $< p_{\text{T}}(\tau) <$ 100 GeV) of the single-tau channel for pair production of up-type leptoquarks $\text{LQ}_{3}^{\text{u}}$.

Efficiency of the first bin of the multi-bin signal region (50 GeV $< p_{\text{T}}(\tau) <$ 100 GeV) of the single-tau channel for pair production of up-type leptoquarks $\text{LQ}_{3}^{\text{u}}$. The plot does not show efficiencies for a branching fraction $B(\text{LQ}_{3}^{\text{u}} \rightarrow b\tau)$ of 0 or 1 because here the acceptance at generator level becomes zero and the efficiency is thus undefined.

Acceptance of the middle bin of the multi-bin signal region (100 GeV $< p_{\text{T}}(\tau) <$ 200 GeV) of the single-tau channel for pair production of up-type leptoquarks $\text{LQ}_{3}^{\text{u}}$.

Efficiency of the middle bin of the multi-bin signal region (100 GeV $< p_{\text{T}}(\tau) <$ 200 GeV) of the single-tau channel for pair production of up-type leptoquarks $\text{LQ}_{3}^{\text{u}}$. The plot does not show efficiencies for a branching fraction $B(\text{LQ}_{3}^{\text{u}} \rightarrow b\tau)$ of 0 or 1 because here the acceptance at generator level becomes zero and the efficiency is thus undefined.

Acceptance of the last bin of the multi-bin signal region (200 GeV $< p_{\text{T}}(\tau)$) of the single-tau channel for pair production of up-type leptoquarks $\text{LQ}_{3}^{\text{u}}$.

Efficiency of the last bin of the multi-bin signal region (200 GeV $< p_{\text{T}}(\tau)$) of the single-tau channel for pair production of up-type leptoquarks $\text{LQ}_{3}^{\text{u}}$. The plot does not show efficiencies for a branching fraction $B(\text{LQ}_{3}^{\text{u}} \rightarrow b\tau)$ of 0 or 1 because here the acceptance at generator level becomes zero and the efficiency is thus undefined.

Acceptance of the signal region of the di-tau channel for pair production of up-type leptoquarks $\text{LQ}_{3}^{\text{u}}$.

Efficiency of the signal region of the di-tau channel for pair production of up-type leptoquarks $\text{LQ}_{3}^{\text{u}}$. The plot does not show efficiencies for a branching fraction $B(\text{LQ}_{3}^{\text{u}} \rightarrow b\tau)$ of 0 because here the acceptance at generator level becomes zero and the efficiency is thus undefined.

Acceptance of the one-bin signal region of the single-tau channel for pair production of down-type leptoquarks $\text{LQ}_{3}^{\text{d}}$.

Efficiency of the one-bin signal region of the single-tau channel for pair production of down-type leptoquarks $\text{LQ}_{3}^{\text{d}}$. The plot does not show efficiencies for a branching fraction $B(\text{LQ}_{3}^{\text{d}} \rightarrow t\tau)$ of 0 or 1 because here the acceptance at generator level becomes zero and the efficiency is thus undefined.

Acceptance of the first bin of the multi-bin signal region (50 GeV $< p_{\text{T}}(\tau) <$ 100 GeV) of the single-tau channel for pair production of down-type leptoquarks $\text{LQ}_{3}^{\text{d}}$.

Efficiency of the first bin of the multi-bin signal region (50 GeV $< p_{\text{T}}(\tau) <$ 100 GeV) of the single-tau channel for pair production of down-type leptoquarks $\text{LQ}_{3}^{\text{d}}$. The plot does not show efficiencies for a branching fraction $B(\text{LQ}_{3}^{\text{d}} \rightarrow t\tau)$ of 0 or 1 because here the acceptance at generator level becomes zero and the efficiency is thus undefined.

Acceptance of the middle bin of the multi-bin signal region (100 GeV $< p_{\text{T}}(\tau) <$ 200 GeV) of the single-tau channel for pair production of down-type leptoquarks $\text{LQ}_{3}^{\text{d}}$.

Efficiency of the middle bin of the multi-bin signal region (100 GeV $< p_{\text{T}}(\tau) <$ 200 GeV) of the single-tau channel for pair production of down-type leptoquarks $\text{LQ}_{3}^{\text{d}}$. The plot does not show efficiencies for a branching fraction $B(\text{LQ}_{3}^{\text{d}} \rightarrow t\tau)$ of 0 or 1 because here the acceptance at generator level becomes zero and the efficiency is thus undefined.

Acceptance of the last bin of the multi-bin signal region (200 GeV $< p_{\text{T}}(\tau)$) of the single-tau channel for pair production of down-type leptoquarks $\text{LQ}_{3}^{\text{d}}$.

Efficiency of the last bin of the multi-bin signal region (200 GeV $< p_{\text{T}}(\tau)$) of the single-tau channel for pair production of down-type leptoquarks $\text{LQ}_{3}^{\text{d}}$. The plot does not show efficiencies for a branching fraction $B(\text{LQ}_{3}^{\text{d}} \rightarrow t\tau)$ of 0 or 1 because here the acceptance at generator level becomes zero and the efficiency is thus undefined.

Acceptance of the signal region of the di-tau channel for pair production of down-type leptoquarks $\text{LQ}_{3}^{\text{d}}$.

Efficiency of the signal region of the di-tau channel for pair production of down-type leptoquarks $\text{LQ}_{3}^{\text{d}}$. The plot does not show efficiencies for a branching fraction $B(\text{LQ}_{3}^{\text{d}} \rightarrow t\tau)$ of 0 because here the acceptance at generator level becomes zero and the efficiency is thus undefined.

Acceptance of the one-bin signal region of the single-tau channel for pair production of vector leptoquarks $\text{LQ}_{3}^{\text{v}}$ in the minimal-coupling scenario.

Efficiency of the one-bin signal region of the single-tau channel for pair production of vector leptoquarks $\text{LQ}_{3}^{\text{v}}$ in the minimal-coupling scenario. The plot does not show efficiencies for a branching fraction $B(\text{LQ}_{3}^{\text{v}} \rightarrow b\tau)$ of 0 or 1 because here the acceptance at generator level becomes zero and the efficiency is thus undefined.

Acceptance of the first bin of the multi-bin signal region (50 GeV $< p_{\text{T}}(\tau) <$ 100 GeV) of the single-tau channel for pair production of vector leptoquarks $\text{LQ}_{3}^{\text{v}}$ in the minimal-coupling scenario.

Efficiency of the first bin of the multi-bin signal region (50 GeV $< p_{\text{T}}(\tau) <$ 100 GeV) of the single-tau channel for pair production of vector leptoquarks $\text{LQ}_{3}^{\text{v}}$ in the minimal-coupling scenario. The plot does not show efficiencies for a branching fraction $B(\text{LQ}_{3}^{\text{v}} \rightarrow b\tau)$ of 0 or 1 because here the acceptance at generator level becomes zero and the efficiency is thus undefined.

Acceptance of the middle bin of the multi-bin signal region (100 GeV $< p_{\text{T}}(\tau) <$ 200 GeV) of the single-tau channel for pair production of vector leptoquarks $\text{LQ}_{3}^{\text{v}}$ in the minimal-coupling scenario.

Efficiency of the middle bin of the multi-bin signal region (100 GeV $< p_{\text{T}}(\tau) <$ 200 GeV) of the single-tau channel for pair production of vector leptoquarks $\text{LQ}_{3}^{\text{v}}$ in the minimal-coupling scenario. The plot does not show efficiencies for a branching fraction $B(\text{LQ}_{3}^{\text{v}} \rightarrow b\tau)$ of 0 or 1 because here the acceptance at generator level becomes zero and the efficiency is thus undefined.

Acceptance of the last bin of the multi-bin signal region (200 GeV $< p_{\text{T}}(\tau)$) of the single-tau channel for pair production of vector leptoquarks $\text{LQ}_{3}^{\text{v}}$ in the minimal-coupling scenario.

Efficiency of the last bin of the multi-bin signal region (200 GeV $< p_{\text{T}}(\tau)$) of the single-tau channel for pair production of vector leptoquarks $\text{LQ}_{3}^{\text{v}}$ in the minimal-coupling scenario. The plot does not show efficiencies for a branching fraction $B(\text{LQ}_{3}^{\text{v}} \rightarrow b\tau)$ of 0 or 1 because here the acceptance at generator level becomes zero and the efficiency is thus undefined.

Acceptance of the signal region of the di-tau channel for pair production of vector leptoquarks $\text{LQ}_{3}^{\text{v}}$ in the minimal-coupling scenario.

Efficiency of the signal region of the di-tau channel for pair production of vector leptoquarks $\text{LQ}_{3}^{\text{v}}$ in the minimal-coupling scenario. The plot does not show efficiencies for a branching fraction $B(\text{LQ}_{3}^{\text{v}} \rightarrow b\tau)$ of 0 because here the acceptance at generator level becomes zero and the efficiency is thus undefined.

Acceptance of the one-bin signal region of the single-tau channel for pair production of vector leptoquarks $\text{LQ}_{3}^{\text{v}}$ with additional gauge couplings.

Efficiency of the one-bin signal region of the single-tau channel for pair production of vector leptoquarks $\text{LQ}_{3}^{\text{v}}$ with additional gauge couplings. The plot does not show efficiencies for a branching fraction $B(\text{LQ}_{3}^{\text{v}} \rightarrow b\tau)$ of 0 or 1 because here the acceptance at generator level becomes zero and the efficiency is thus undefined.

Acceptance of the first bin of the multi-bin signal region (50 GeV $< p_{\text{T}}(\tau) <$ 100 GeV) of the single-tau channel for pair production of vector leptoquarks $\text{LQ}_{3}^{\text{v}}$ with additional gauge couplings.

Efficiency of the first bin of the multi-bin signal region (50 GeV $< p_{\text{T}}(\tau) <$ 100 GeV) of the single-tau channel for pair production of vector leptoquarks $\text{LQ}_{3}^{\text{v}}$ with additional gauge couplings. The plot does not show efficiencies for a branching fraction $B(\text{LQ}_{3}^{\text{v}} \rightarrow b\tau)$ of 0 or 1 because here the acceptance at generator level becomes zero and the efficiency is thus undefined.

Acceptance of the middle bin of the multi-bin signal region (100 GeV $< p_{\text{T}}(\tau) <$ 200 GeV) of the single-tau channel for pair production of vector leptoquarks $\text{LQ}_{3}^{\text{v}}$ with additional gauge couplings.

Efficiency of the middle bin of the multi-bin signal region (100 GeV $< p_{\text{T}}(\tau) <$ 200 GeV) of the single-tau channel for pair production of vector leptoquarks $\text{LQ}_{3}^{\text{v}}$ with additional gauge couplings. The plot does not show efficiencies for a branching fraction $B(\text{LQ}_{3}^{\text{v}} \rightarrow b\tau)$ of 0 or 1 because here the acceptance at generator level becomes zero and the efficiency is thus undefined.

Acceptance of the last bin of the multi-bin signal region (200 GeV $< p_{\text{T}}(\tau)$) of the single-tau channel for pair production of vector leptoquarks $\text{LQ}_{3}^{\text{v}}$ with additional gauge couplings.

Efficiency of the last bin of the multi-bin signal region (200 GeV $< p_{\text{T}}(\tau)$) of the single-tau channel for pair production of vector leptoquarks $\text{LQ}_{3}^{\text{v}}$ with additional gauge couplings. The plot does not show efficiencies for a branching fraction $B(\text{LQ}_{3}^{\text{v}} \rightarrow b\tau)$ of 0 or 1 because here the acceptance at generator level becomes zero and the efficiency is thus undefined.

Acceptance of the signal region of the di-tau channel for pair production of vector leptoquarks $\text{LQ}_{3}^{\text{v}}$ with additional gauge couplings.

Efficiency of the signal region of the di-tau channel for pair production of vector leptoquarks $\text{LQ}_{3}^{\text{v}}$ with additional gauge couplings. The plot does not show efficiencies for a branching fraction $B(\text{LQ}_{3}^{\text{v}} \rightarrow b\tau)$ of 0 because here the acceptance at generator level becomes zero and the efficiency is thus undefined.

Acceptance of the one-bin signal region of the single-tau channel for pair production of top squarks with decays via tau sleptons.

Efficiency of the one-bin signal region of the single-tau channel for pair production of top squarks with decays via tau sleptons.

Acceptance of the first bin of the multi-bin signal region (50 GeV $< p_{\text{T}}(\tau) <$ 100 GeV) of the single-tau channel for pair production of top squarks with decays via tau sleptons.

Efficiency of the first bin of the multi-bin signal region (50 GeV $< p_{\text{T}}(\tau) <$ 100 GeV) of the single-tau channel for pair production of top squarks with decays via tau sleptons.

Acceptance of the middle bin of the multi-bin signal region (100 GeV $< p_{\text{T}}(\tau) <$ 200 GeV) of the single-tau channel for pair production of top squarks with decays via tau sleptons.

Efficiency of the middle bin of the multi-bin signal region (100 GeV $< p_{\text{T}}(\tau) <$ 200 GeV) of the single-tau channel for pair production of top squarks with decays via tau sleptons.

Acceptance of the last bin of the multi-bin signal region (200 GeV $< p_{\text{T}}(\tau)$) of the single-tau channel for pair production of top squarks with decays via tau sleptons.

Efficiency of the last bin of the multi-bin signal region (200 GeV $< p_{\text{T}}(\tau)$) of the single-tau channel for pair production of top squarks with decays via tau sleptons.

Acceptance of the signal region of the di-tau channel for pair production of top squarks with decays via tau sleptons.

Efficiency of the signal region of the di-tau channel for pair production of top squarks with decays via tau sleptons.

Cutflow for the benchmark signal model $m(\tilde{t}_{1}) = 1350$ GeV, $m(\tilde{\tau}_{1}) = 1090$ GeV for the di-tau SR. The simulated sample contains 30,000 raw MC events. Weighted event yields are reported, normalized to an integrated luminosity of 139 fb$^{-1}$. 'Preselection' refers to the preselection for the di-tau channel.

Cutflow for the benchmark signal model $m(\tilde{t}_{1}) = 1350$ GeV, $m(\tilde{\tau}_{1}) = 1090$ GeV for the single-tau one-bin SR. The simulated sample contains 30,000 raw MC events. Weighted event yields are reported, normalized to an integrated luminosity of 139 fb$^{-1}$. 'Preselection' refers to the preselection for the single-tau channel.

Cutflow for the benchmark signal model $m(\tilde{t}_{1}) = 1350$ GeV, $m(\tilde{\tau}_{1}) = 1090$ GeV for the single-tau multi-bin SR. The simulated sample contains 30,000 raw MC events. Weighted event yields are reported, normalized to an integrated luminosity of 139 fb$^{-1}$. 'Preselection' refers to the preselection for the single-tau channel.

Cutflow for the benchmark signal model $m(\text{LQ}_{3}^{\text{u}}) = 1.2$ TeV, $\beta = 0.5$ for the di-tau SR. The simulated sample contains 210,000 raw MC events. Weighted event yields are reported, normalized to an integrated luminosity of 139 fb$^{-1}$. 'Preselection' refers to the preselection for the di-tau channel.

Cutflow for the benchmark signal model $m(\text{LQ}_{3}^{\text{u}}) = 1.2$ TeV, $\beta = 0.5$ for the single-tau one-bin SR. The simulated sample contains 210,000 raw MC events. Weighted event yields are reported, normalized to an integrated luminosity of 139 fb$^{-1}$. 'Preselection' refers to the preselection for the single-tau channel.

Cutflow for the benchmark signal model $m(\text{LQ}_{3}^{\text{u}}) = 1.2$ TeV, $\beta = 0.5$ for the single-tau multi-bin SR. The simulated sample contains 210,000 raw MC events. Weighted event yields are reported, normalized to an integrated luminosity of 139 fb$^{-1}$. 'Preselection' refers to the preselection for the single-tau channel.

Cutflow for the benchmark signal model $m(\text{LQ}_{3}^{\text{d}}) = 1.2$ TeV, $\beta = 0.5$ for the di-tau SR. The simulated sample contains 210,000 raw MC events. Weighted event yields are reported, normalized to an integrated luminosity of 139 fb$^{-1}$. 'Preselection' refers to the preselection for the di-tau channel.

Cutflow for the benchmark signal model $m(\text{LQ}_{3}^{\text{d}}) = 1.2$ TeV, $\beta = 0.5$ for the single-tau one-bin SR. The simulated sample contains 210,000 raw MC events. Weighted event yields are reported, normalized to an integrated luminosity of 139 fb$^{-1}$. 'Preselection' refers to the preselection for the single-tau channel.

Cutflow for the benchmark signal model $m(\text{LQ}_{3}^{\text{d}}) = 1.2$ TeV, $\beta = 0.5$ for the single-tau multi-bin SR. The simulated sample contains 210,000 raw MC events. Weighted event yields are reported, normalized to an integrated luminosity of 139 fb$^{-1}$. 'Preselection' refers to the preselection for the single-tau channel.

Cutflow for the benchmark signal model $m(\text{LQ}_{3}^{\text{v}}) = 1.4$ TeV, $\beta = 0.5$ in the minimal-coupling scenario for the di-tau SR. The simulated sample contains 50,000 raw MC events. Weighted event yields are reported, normalized to an integrated luminosity of 139 fb$^{-1}$. 'Preselection' refers to the preselection for the di-tau channel.

Cutflow for the benchmark signal model $m(\text{LQ}_{3}^{\text{v}}) = 1.4$ TeV, $\beta = 0.5$ in the minimal-coupling scenario for the single-tau one-bin SR. The simulated sample contains 50,000 raw MC events. Weighted event yields are reported, normalized to an integrated luminosity of 139 fb$^{-1}$. 'Preselection' refers to the preselection for the single-tau channel.

Cutflow for the benchmark signal model $m(\text{LQ}_{3}^{\text{v}}) = 1.4$ TeV, $\beta = 0.5$ in the minimal-coupling scenario for the single-tau multi-bin SR. The simulated sample contains 50,000 raw MC events. Weighted event yields are reported, normalized to an integrated luminosity of 139 fb$^{-1}$. 'Preselection' refers to the preselection for the single-tau channel.

Cutflow for the benchmark signal model $m(\text{LQ}_{3}^{\text{v}}) = 1.4$ TeV, $\beta = 0.5$ in the Yang--Mills scenario for the di-tau SR. The simulated sample contains 50,000 raw MC events. Weighted event yields are reported, normalized to an integrated luminosity of 139 fb$^{-1}$. 'Preselection' refers to the preselection for the di-tau channel.

Cutflow for the benchmark signal model $m(\text{LQ}_{3}^{\text{v}}) = 1.4$ TeV, $\beta = 0.5$ in the Yang--Mills scenario for the single-tau one-bin SR. The simulated sample contains 50,000 raw MC events. Weighted event yields are reported, normalized to an integrated luminosity of 139 fb$^{-1}$. 'Preselection' refers to the preselection for the single-tau channel.

Cutflow for the benchmark signal model $m(\text{LQ}_{3}^{\text{v}}) = 1.4$ TeV, $\beta = 0.5$ in the Yang--Mills scenario for the single-tau multi-bin SR. The simulated sample contains 50,000 raw MC events. Weighted event yields are reported, normalized to an integrated luminosity of 139 fb$^{-1}$. 'Preselection' refers to the preselection for the single-tau channel.

Distribution of the transverse mass of the $W$-candidate $m_{\rm T}^{W}$ for the one lepton and two $b$-jets channel in SRlbb-1 and SRlbb-2, $m_{bb}$ sidebands. The background histograms are obtained from the background-only fit, and their uncertainty represents the total background uncertainty after the fit. The last bin includes overflow.

Distribution of the number of $b$-jets for the one lepton and two $b$-jets channel in SRlbb-1 and SRlbb-2 without the b-jet multiplicity requirement, central $m_{bb}$ bin. The background histograms are obtained from the background-only fit, and their uncertainty represents the total background uncertainty after the fit.

Distribution of the invariant mass of the $b$-jets $m_{bb}$ for the one lepton and two $b$-jets channel in the SRlbb-1 and SRlbb-2. The background histograms are obtained from the background-only fit, and their uncertainty represents the total background uncertainty after the fit.

Distribution of missing transverse momentum $E_{\rm T}^{\rm miss}$ in the one lepton and two photons signal regions SR$\ell\gamma\gamma$-1 and SR$\ell\gamma\gamma$-2 for the Higgs-mass window ($120\lt m_{\gamma\gamma} \lt 130$ GeV). The final column (background) is a simulation-based cross check. The contributions from non-Higgs backgrounds are scaled by 10 GeV / 50 GeV = 0.2 from the $m_{\gamma\gamma}$ sideband ($100 \lt m_{\gamma\gamma} \lt 120$ GeV and $130 \lt m_{\gamma\gamma} \lt 160$ GeV) into the Higgs-mass window. The last bin includes overflow. The distributions of a signal hypothesis are also shown for $m(\tilde{\chi}^{\pm}_{1}\tilde{\chi}^{0}_{2},\tilde{\chi}^{0}_{1})=(165,35)$ GeV.

Distribution of the azimuth difference between the $W$ and Higgs boson candidates $\Delta\phi(W,h)$ in the one lepton and two photons signal regions SR$\ell\gamma\gamma$-1 and SR$\ell\gamma\gamma$-2 for the Higgs-mass window ($120\lt m_{\gamma\gamma} \lt 130$ GeV). The final column (background) is a simulation-based cross check. The contributions from non-Higgs backgrounds are scaled by 10 GeV / 50 GeV = 0.2 from the $m_{\gamma\gamma}$ sideband ($100 \lt m_{\gamma\gamma} \lt 120$ GeV and $130 \lt m_{\gamma\gamma} \lt 160$ GeV) into the Higgs-mass window. The distributions of a signal hypothesis are also shown for $m(\tilde{\chi}^{\pm}_{1}\tilde{\chi}^{0}_{2},\tilde{\chi}^{0}_{1})=(165,35)$ GeV.

Distribution of the transverse mass of the $W$ and photon system $m_{\rm{T}}^{W\gamma_1}$ in the one lepton and two photons signal regions SR$\ell\gamma\gamma$-1 and SR$\ell\gamma\gamma$-2 for the Higgs-mass window ($120\lt m_{\gamma\gamma} \lt 130$ GeV). The final column (background) is a simulation-based cross check. The contributions from non-Higgs backgrounds are scaled by 10 GeV / 50 GeV = 0.2 from the $m_{\gamma\gamma}$ sideband ($100 \lt m_{\gamma\gamma} \lt 120$ GeV and $130 \lt m_{\gamma\gamma} \lt 160$ GeV) into the Higgs-mass window. The last bin includes overflow. The distributions of a signal hypothesis are also shown for $m(\tilde{\chi}^{\pm}_{1}\tilde{\chi}^{0}_{2},\tilde{\chi}^{0}_{1})=(165,35)$ GeV.

Distribution of the transverse mass of the $W$ and photon system $m_{\rm{T}}^{W\gamma_2}$ in the one lepton and two photons signal regions SR$\ell\gamma\gamma$-1 and SR$\ell\gamma\gamma$-2 for the Higgs-mass window ($120\lt m_{\gamma\gamma} \lt 130$ GeV). The final column (background) is a simulation-based cross check. The contributions from non-Higgs backgrounds are scaled by 10 GeV / 50 GeV = 0.2 from the $m_{\gamma\gamma}$ sideband ($100 \lt m_{\gamma\gamma} \lt 120$ GeV and $130 \lt m_{\gamma\gamma} \lt 160$ GeV) into the Higgs-mass window. The last bin includes overflow. The distributions of a signal hypothesis are also shown for $m(\tilde{\chi}^{\pm}_{1}\tilde{\chi}^{0}_{2},\tilde{\chi}^{0}_{1})=(165,35)$ GeV.

Results of the background-only fit to the diphoton invariant mass, $m_{\gamma\gamma}$, distribution in the one lepton and two photons signal region SR$l\gamma\gamma$-1. The contributions from SM Higgs boson production are constrained to the MC prediction and associated systematic uncertainties. The fit is performed on events with 100 GeV $ \lt m_{\gamma\gamma} \lt $ 160 GeV, with events in SR$l\gamma\gamma$-1 or SR$l\gamma\gamma$-2 in the Higgs-mass window (120 GeV $\le m_{\gamma\gamma} \le$ 130 GeV) excluded from the fit.

Results of the background-only fit to the diphoton invariant mass, $m_{\gamma\gamma}$, distribution in the one lepton and two photons signal region SR$l\gamma\gamma$-2. The contributions from SM Higgs boson production are constrained to the MC prediction and associated systematic uncertainties. The fit is performed on events with 100 GeV $ \lt m_{\gamma\gamma} \lt $ 160 GeV, with events in SR$l\gamma\gamma$-1 or SR$l\gamma\gamma$-2 in the Higgs-mass window (120 GeV $\le m_{\gamma\gamma} \le$ 130 GeV) excluded from the fit.

Results of the background-only fit to the diphoton invariant mass, $m_{\gamma\gamma}$, distribution in the one lepton and two photons validation region VR$l\gamma\gamma$-1. The contributions from SM Higgs boson production are constrained to the MC prediction and associated systematic uncertainties. The fit is performed on events with 100 GeV $ \lt m_{\gamma\gamma} \lt $ 160 GeV, with events in SR$l\gamma\gamma$-1 or SR$l\gamma\gamma$-2 in the Higgs-mass window (120 GeV $\le m_{\gamma\gamma} \le$ 130 GeV) excluded from the fit.

Results of the background-only fit to the diphoton invariant mass, $m_{\gamma\gamma}$, distribution in the one lepton and two photons validation region VR$l\gamma\gamma$-2. The contributions from SM Higgs boson production are constrained to the MC prediction and associated systematic uncertainties. The fit is performed on events with 100 GeV $ \lt m_{\gamma\gamma} \lt $ 160 GeV, with events in SR$l\gamma\gamma$-1 or SR$l\gamma\gamma$-2 in the Higgs-mass window (120 GeV $\le m_{\gamma\gamma} \le$ 130 GeV) excluded from the fit.

Distribution of effective mass $m_{\rm eff}$ in the validation region of the same-sign $e\mu$ channel. This validation region is defined by requiring one, two, or three jets, and reversing the $m_{lj}$, $m_{ljj}$ criteria. The distribution of a signal hypothesis is also shown.

Distribution of effective mass $m_{\rm eff}$ for the same-sign dilepton channel in the signal regions with one jet SR$ll$-1. SR$ll$-1 is the sum of SR$ee$-1, SR$e\mu$-1, and SR$\mu\mu$-1. All selection criteria are applied, except for the one on $m_{\rm eff}$. The distributions of a signal hypothesis are also shown for $m(\tilde{\chi}^{\pm}_{1}\tilde{\chi}^{0}_{2},\tilde{\chi}^{0}_{1})=(130,0)$ GeV. The last bin includes overflow.

Distribution of effective mass $m_{\rm eff}$ for the same-sign dilepton channel in the signal regions with two or three jets SR$ll$-2. SR$ll$-2 is the sum of SR$ee$-2, SR$e\mu$-2, and SR$\mu\mu$-2. All selection criteria are applied, except for the one on $m_{\rm eff}$. The distributions of a signal hypothesis are also shown for $m(\tilde{\chi}^{\pm}_{1}\tilde{\chi}^{0}_{2},\tilde{\chi}^{0}_{1})=(130,0)$ GeV. The last bin includes overflow.

Distribution of largest transverse mass $m_{\rm T}^{\rm max}$ for the same-sign dilepton channel in the signal regions with one jet SR$ll$-1. SR$ll$-1 is the sum of SR$ee$-1, SR$e\mu$-1, and SR$\mu\mu$-1. All selection criteria are applied, except for the one on $m_{\rm T}^{\rm max}$. The distributions of a signal hypothesis are also shown for $m(\tilde{\chi}^{\pm}_{1}\tilde{\chi}^{0}_{2},\tilde{\chi}^{0}_{1})=(130,0)$ GeV. The last bin includes overflow.

Distribution of largest transverse mass $m_{\rm T}^{\rm max}$ for the same-sign dilepton channel in the signal regions with two or three jets SR$ll$-2. SR$ll$-2 is the sum of SR$ee$-2, SR$e\mu$-2, and SR$\mu\mu$-2. All selection criteria are applied, except for the one on $m_{\rm T}^{\rm max}$. The distributions of a signal hypothesis are also shown for $m(\tilde{\chi}^{\pm}_{1}\tilde{\chi}^{0}_{2},\tilde{\chi}^{0}_{1})=(130,0)$ GeV. The last bin includes overflow.

Distribution of invariant mass of lepton and jet $m_{lj}$ for the same-sign dilepton channel in the signal regions with one jet SR$ll$-1. SR$ll$-1 is the sum of SR$ee$-1, SR$e\mu$-1, and SR$\mu\mu$-1. All selection criteria are applied, except for the one on $m_{lj}$. The distributions of a signal hypothesis are also shown for $m(\tilde{\chi}^{\pm}_{1}\tilde{\chi}^{0}_{2},\tilde{\chi}^{0}_{1})=(130,0)$ GeV. The last bin includes overflow.

Distribution of invariant mass of lepton and di-jet $m_{ljj}$ for the same-sign dilepton channel in the signal regions with two or three jets SR$ll$-2. SR$ll$-2 is the sum of SR$ee$-2, SR$e\mu$-2, and SR$\mu\mu$-2. All selection criteria are applied, except for the one on $m_{ljj}$. The distributions of a signal hypothesis are also shown for $m(\tilde{\chi}^{\pm}_{1}\tilde{\chi}^{0}_{2},\tilde{\chi}^{0}_{1})=(130,0)$ GeV. The last bin includes overflow.

One lepton and two $b$-jets channel: observed and expected 95% CL upper limits on the cross section normalised by the simplified model prediction as a function of the common mass $m_{\tilde{\chi}_1^\pm \tilde{\chi}^0_2}$ for $m_{\tilde{\chi}^0_1}=0$.

One lepton and two photons channel: observed and expected 95% CL upper limits on the cross section normalised by the simplified model prediction as a function of the common mass $m_{\tilde{\chi}_1^\pm \tilde{\chi}^0_2}$ for $m_{\tilde{\chi}^0_1}=0$.

Same-sign dilepton channel: observed and expected 95% CL upper limits on the cross section normalised by the simplified model prediction as a function of the common mass $m_{\tilde{\chi}_1^\pm \tilde{\chi}^0_2}$ for $m_{\tilde{\chi}^0_1}=0$.

Combination: observed and expected 95% CL upper limits on the cross section normalised by the simplified model prediction as a function of the common mass $m_{\tilde{\chi}_1^\pm \tilde{\chi}^0_2}$ for $m_{\tilde{\chi}^0_1}=0$. This combination is obtained using the result from the ATLAS three-lepton search, J. High Energy Phys. 04 (2014) 169, in addition to the three channels reported in this paper.

One lepton and two b-jets channel: Expected 95% CL exclusion contour for chargino neutralino production via Wh.

One lepton and two b-jets channel: Observed 95% CL exclusion contour for chargino neutralino production via Wh.

One lepton and two photons channel: Expected 95% CL exclusion contour for chargino neutralino production via Wh.

One lepton and two photons channel: Observed 95% CL exclusion contour for chargino neutralino production via Wh.

Same-sign dilepton channel: Expected 95% CL exclusion contour for chargino neutralino production via Wh.

Same-sign dilepton channel: Observed 95% CL exclusion contour for chargino neutralino production via Wh.

Combination: Expected 95% CL exclusion contour for chargino neutralino production via Wh.

Combination: Observed 95% CL exclusion contour for chargino neutralino production via Wh.

Combination: excluded model cross-section at 95% CL IN PB.

Acceptance SR$\ell\gamma\gamma$-1.

Acceptance SR$\ell\gamma\gamma$-2.

Efficiency SR$\ell\gamma\gamma$-1.

Efficiency SR$\ell\gamma\gamma$-2.

$m_{\gamma\gamma}$ distribution in the SR$l\gamma\gamma$-1 region for the full $m_{\gamma\gamma}$ window. The last column (background) is from a simulation-based cross check. $Z\gamma$ events, with $Z\rightarrow{ee}$ and one electron mis-identified as a photon, are included in the $Z\gamma(\gamma)$ background. The distributions of a signal hypothesis are also shown for $m(\tilde{\chi}^{\pm}_{1}\tilde{\chi}^{0}_{2},\tilde{\chi}^{0}_{1})=(165,35)$ GeV.

$m_{\gamma\gamma}$ distribution in the SR$l\gamma\gamma$-2 region for the full $m_{\gamma\gamma}$ window. The last column (background) is from a simulation-based cross check. $Z\gamma$ events, with $Z\rightarrow{ee}$ and one electron mis-identified as a photon, are included in the $Z\gamma(\gamma)$ background. The distributions of a signal hypothesis are also shown for $m(\tilde{\chi}^{\pm}_{1}\tilde{\chi}^{0}_{2},\tilde{\chi}^{0}_{1})=(165,35)$ GeV.

$m_{\gamma\gamma}$ distribution in the VR$l\gamma\gamma$-1 region for the full $m_{\gamma\gamma}$ window. The last column (background) is from a simulation-based cross check. $Z\gamma$ events, with $Z\rightarrow{ee}$ and one electron mis-identified as a photon, are included in the $Z\gamma(\gamma)$ background. The distributions of a signal hypothesis are also shown for $m(\tilde{\chi}^{\pm}_{1}\tilde{\chi}^{0}_{2},\tilde{\chi}^{0}_{1})=(165,35)$ GeV.

$m_{\gamma\gamma}$ distribution in the VR$l\gamma\gamma$-2 region for the full $m_{\gamma\gamma}$ window. The last column (background) is from a simulation-based cross check. $Z\gamma$ events, with $Z\rightarrow{ee}$ and one electron mis-identified as a photon, are included in the $Z\gamma(\gamma)$ background. The distributions of a signal hypothesis are also shown for $m(\tilde{\chi}^{\pm}_{1}\tilde{\chi}^{0}_{2},\tilde{\chi}^{0}_{1})=(165,35)$ GeV.

Acceptance for SRlbb-1, central $m_{bb}$ bin.

Acceptance for SRlbb-1, $m_{bb}$ sideband.

Acceptance for SRlbb-2, central $m_{bb}$ bin.

Acceptance for SRlbb-2, $m_{bb}$ sideband.

Efficiency for SRlbb-1, central $m_{bb}$ bin.

Efficiency for SRlbb-1, $m_{bb}$ sideband.

Efficiency for SRlbb-2, central $m_{bb}$ bin.

Efficiency for SRlbb-2, $m_{bb}$ sideband.

Acceptance for the same-sign $ee$ channel with one jet.

Acceptance for the same-sign $ee$ channel with two or three jets.

Acceptance for the same-sign $e\mu$ channel with one jet.

Acceptance for the same-sign $e\mu$ channel with two or three jets.

Acceptance for the same-sign $\mu\mu$ channel with one jet.

Acceptance for the same-sign $\mu\mu$ channel with two or three jets.

Efficiency for the same-sign $ee$ channel with one jet.

Efficiency for the same-sign $ee$ channel with two or three jets.

Efficiency for the same-sign $e\mu$ channel with one jet.

Efficiency for the same-sign $e\mu$ channel with two or three jets.

Efficiency for the same-sign $\mu\mu$ channel with one jet.

Efficiency for the same-sign $\mu\mu$ channel with two or three jets.

Observed $m_{t\bar{t}}^{reco}$ distributions in data for SR$2b$, shown together with the result of the fit with the three-shape-parameter function. The error bars indicate the effect of the fit parameter uncertainty on the background prediction. The bin width of the distributions is chosen to be the same as that used in the background parameterization.

Observed and expected upper limits on the cross-section times branching fraction of the $Z'_{TC2} \rightarrow t\bar{t}$ as a function of the $Z'_{TC2}$ mass ($\Gamma/m = 1.2$%).

Fit parameter values obtained from fits to data in the two signal regions. The fit function is $f(x) = p_{0} (1-x)^{p_{1}} x^{p_{2} + p_{3} \log(x)}$, where $x = m^{reco}_{t\bar{t}}/\sqrt{s}$.

Selection acceptance times efficiency as a function of WPRIME mass at truth level for left- and right-handed WPRIME MC. The total efficiency curves correspond to the sum of the efficiencies of the one b-tag and two b-tag categories.

Cutflow (efficiency with respect to total number of events in %) for several WPRIME masses in the left-handed and in the right-handed model for hadronic top-quark decays.

Overview of the signal acceptance times efficiency (A x Eff) for the hadronic top-quark decay for both categories, the predicted cross section times branching ratio to TOP BOTTOM, as well as the observed 95% CL limit on the cross section times branching ratio for several WPRIME masses in the left-handed and in the right-handed model. The expected signal event yields in the two categories for 20.3 fb-1 of 8 TeV data are also shown.

Expected and observed $am_{T2}$ distribution for bCd_high1 SR, before applying the $am_{T2}>200$ GeV requirement. The uncertainty includes statistical and all experimental systematic uncertainties. The last bin includes overflows.

Expected and observed leading b-jet $p_T$ distribution for bCd_high2 SR, before applying the leading b-jet $p_T>170$ GeV requirement. The uncertainty includes statistical and all experimental systematic uncertainties. The last bin includes overflows.

Expected and observed $E_T^{miss}$ distribution for tNbC_mix SR, before applying the $E_T^{miss}>270$ GeV requirement. The uncertainty includes statistical and all experimental systematic uncertainties. The last bin includes overflows.

Expected and observed lepton $p_T$ distribution for bCa_low SR. The uncertainty includes statistical and all experimental systematic uncertainties. The last bin includes overflows.

Expected and observed lepton $p_T$ distribution for bCa_med SR. The uncertainty includes statistical and all experimental systematic uncertainties. The last bin includes overflows.

Expected and observed $am_T2$ distribution for bCb_med1 SR. The uncertainty includes statistical and all experimental systematic uncertainties. The last bin includes overflows.

Expected and observed $am_T2$ distribution for bCb_high SR. The uncertainty includes statistical and all experimental systematic uncertainties. The last bin includes overflows.

Best expected signal region for the $\tilde t_1\to t\chi^0_1$ scenario with $m_{\tilde t_1}>m_t+m_{\chi^0_1}$. This mapping is used for the final combined exclusion limits.

Best expected signal region for the $\tilde t_1$ three-body scenario ($\tilde t_1\to bW\chi^0_1$). This mapping is used for the final combined exclusion limits.

Best expected signal region for the $\tilde t_1$ four-body scenario ($\tilde t_1\to bff'\chi^0_1$). This mapping is used for the final combined exclusion limits.

Best expected signal region for the $\tilde t_1\to b\chi^\pm_1$ scenario with $m_{\chi^\pm_1}=2\times m_{\chi^0_1}$. This mapping is used for the final combined exclusion limits.

Best expected signal region for the $\tilde t_1\to b\chi^\pm_1$ scenario with $m_{\chi^\pm_1}=150$ GeV. This mapping is used for the final combined exclusion limits.

Best expected signal region for the $\tilde t_1\to b\chi^\pm_1$ scenario with $m_{\chi^\pm_1}=106$ GeV. This mapping is used for the final combined exclusion limits.

Best expected signal region for the $\tilde t_1\to b\chi^\pm_1$ scenario with $m_{\chi^\pm_1}=m_{\chi^0_1}+5$ GeV. This mapping is used for the final combined exclusion limits.

Best expected signal region for the $\tilde t_1\to b\chi^\pm_1$ scenario with $m_{\chi^\pm_1}=m_{\chi^0_1}+20$ GeV. This mapping is used for the final combined exclusion limits.

Best expected signal region for the $\tilde t_1\to b\chi^\pm_1$ scenario with $m_{\chi^\pm_1}=m_{\tilde t_1}-10$ GeV. This mapping is used for the final combined exclusion limits.

Best expected signal region for the $\tilde t_1\to b\chi^\pm_1$ scenario with $m_{\tilde t_1}=300$ GeV. This mapping is used for the final combined exclusion limits.

Upper limits on the model cross-section for the $\tilde t_1\to t\chi^0_1$ scenario with $m_{\tilde t_1}>m_t+m_{\chi^0_1}$.

Observed exclusion contour for the $\tilde t_1\to t\chi^0_1$ scenario with $m_{\tilde t_1}>m_t+m_{\chi^0_1}$.

Expected exclusion contour for the $\tilde t_1\to t\chi^0_1$ scenario with $m_{\tilde t_1}>m_t+m_{\chi^0_1}$.

Upper limit on signal events for the $\tilde t_1$ three-body scenario ($\tilde t_1\to bW\chi^0_1$).

Observed exclusion contour for the $\tilde t_1$ three-body scenario ($\tilde t_1\to bW\chi^0_1$).

Expected exclusion contour for the $\tilde t_1$ three-body scenario ($\tilde t_1\to bW\chi^0_1$).

Upper limit on signal events for the $\tilde t_1$ four-body scenario ($\tilde t_1\to bff'\chi^0_1$).

Observed exclusion contour for the $\tilde t_1$ four-body scenario ($\tilde t_1\to bff'\chi^0_1$).

Expected exclusion contour for the $\tilde t_1$ four-body scenario ($\tilde t_1\to bff'\chi^0_1$).

Upper limit on signal events for the $\tilde t_1\to b\chi^\pm_1$ scenario with $m_{\chi^\pm_1}=2\times m_{\chi^0_1}$.

Observed exclusion contour for the $\tilde t_1\to b\chi^\pm_1$ scenario with $m_{\chi^\pm_1}=2\times m_{\chi^0_1}$.

Expected exclusion contour for the $\tilde t_1\to b\chi^\pm_1$ scenario with $m_{\chi^\pm_1}=2\times m_{\chi^0_1}$.

Upper limit on signal events for the $\tilde t_1\to b\chi^\pm_1$ scenario with $m_{\chi^\pm_1}=150$ GeV.

Observed exclusion contour for the $\tilde t_1\to b\chi^\pm_1$ scenario with $m_{\chi^\pm_1}=150$ GeV.

Expected exclusion contour for the $\tilde t_1\to b\chi^\pm_1$ scenario with $m_{\chi^\pm_1}=150$ GeV.

Upper limit on signal events for the $\tilde t_1\to b\chi^\pm_1$ scenario with $m_{\chi^\pm_1}=106$ GeV.

Observed exclusion contour for the $\tilde t_1\to b\chi^\pm_1$ scenario with $m_{\chi^\pm_1}=106$ GeV.

Expected exclusion contour for the $\tilde t_1\to b\chi^\pm_1$ scenario with $m_{\chi^\pm_1}=106$ GeV.

Upper limit on signal events for the $\tilde t_1\to b\chi^\pm_1$ scenario with $m_{\chi^\pm_1}=m_{\chi^0_1}+5$ GeV.

Observed exclusion contour for the $\tilde t_1\to b\chi^\pm_1$ scenario with $m_{\chi^\pm_1}=m_{\chi^0_1}+5$ GeV.

Expected exclusion contour for the $\tilde t_1\to b\chi^\pm_1$ scenario with $m_{\chi^\pm_1}=m_{\chi^0_1}+5$ GeV.

Upper limit on signal events for the $\tilde t_1\to b\chi^\pm_1$ scenario with $m_{\chi^\pm_1}=m_{\chi^0_1}+20$ GeV.

Observed exclusion contour for the $\tilde t_1\to b\chi^\pm_1$ scenario with $m_{\chi^\pm_1}=m_{\chi^0_1}+20$ GeV.

Expected exclusion contour for the $\tilde t_1\to b\chi^\pm_1$ scenario with $m_{\chi^\pm_1}=m_{\chi^0_1}+20$ GeV.

Upper limit on signal events for the $\tilde t_1\to b\chi^\pm_1$ scenario with $m_{\chi^\pm_1}=m_{\tilde t_1}-10$ GeV.

Observed exclusion contour for the $\tilde t_1\to b\chi^\pm_1$ scenario with $m_{\chi^\pm_1}=m_{\tilde t_1}-10$ GeV.

Expected exclusion contour for the $\tilde t_1\to b\chi^\pm_1$ scenario with $m_{\chi^\pm_1}=m_{\tilde t_1}-10$ GeV.

Upper limit on signal events for the $\tilde t_1\to b\chi^\pm_1$ scenario with $m_{\tilde t_1}=300$ GeV.

Observed exclusion contour for the $\tilde t_1\to b\chi^\pm_1$ scenario with $m_{\tilde t_1}=300$ GeV.

Expected exclusion contour for the $\tilde t_1\to b\chi^\pm_1$ scenario with $m_{\tilde t_1}=300$ GeV.

Acceptance of tN_diag SR ($E_T^{miss}>150$ GeV, $m_T>140$ GeV) for the $\tilde t_1\to t\chi^0_1$ scenario with $m_{\tilde t_1}>m_t+m_{\chi^0_1}$. The acceptance is defined as the fraction of signal events that pass the analysis selection performed on generator-level objects, therefore emulating an ideal detector with perfect particle identification and no measurement resolution effects.

Acceptance of tN_med SR for the $\tilde t_1\to t\chi^0_1$ scenario with $m_{\tilde t_1}>m_t+m_{\chi^0_1}$. The acceptance is defined as the fraction of signal events that pass the analysis selection performed on generator-level objects, therefore emulating an ideal detector with perfect particle identification and no measurement resolution effects.

Acceptance of tN_boost SR for the $\tilde t_1\to t\chi^0_1$ scenario with $m_{\tilde t_1}>m_t+m_{\chi^0_1}$. The acceptance is defined as the fraction of signal events that pass the analysis selection performed on generator-level objects, therefore emulating an ideal detector with perfect particle identification and no measurement resolution effects.

Acceptance of bCb_med2 SR ($am_{T2}>250$ GeV, $m_T>60$ GeV) for the $\tilde t_1\to b\chi^\pm_1$ scenario with $m_{\chi^\pm_1}=2\times m_{\chi^0_1}$. The acceptance is defined as the fraction of signal events that pass the analysis selection performed on generator-level objects, therefore emulating an ideal detector with perfect particle identification and no measurement resolution effects.

Acceptance of bCc_diag SR for the $\tilde t_1\to b\chi^\pm_1$ scenario with $m_{\chi^\pm_1}=2\times m_{\chi^0_1}$. The acceptance is defined as the fraction of signal events that pass the analysis selection performed on generator-level objects, therefore emulating an ideal detector with perfect particle identification and no measurement resolution effects.

Acceptance of bCd_bulk SR ($am_{T2}>175$ GeV, $m_T>120$ GeV) for the $\tilde t_1\to b\chi^\pm_1$ scenario with $m_{\chi^\pm_1}=2\times m_{\chi^0_1}$. The acceptance is defined as the fraction of signal events that pass the analysis selection performed on generator-level objects, therefore emulating an ideal detector with perfect particle identification and no measurement resolution effects.

Acceptance of bCd_high1 SR for the $\tilde t_1\to b\chi^\pm_1$ scenario with $m_{\chi^\pm_1}=2\times m_{\chi^0_1}$. The acceptance is defined as the fraction of signal events that pass the analysis selection performed on generator-level objects, therefore emulating an ideal detector with perfect particle identification and no measurement resolution effects.

Acceptance of bCd_high2 SR for the $\tilde t_1\to b\chi^\pm_1$ scenario with $m_{\chi^\pm_1}=2\times m_{\chi^0_1}$. The acceptance is defined as the fraction of signal events that pass the analysis selection performed on generator-level objects, therefore emulating an ideal detector with perfect particle identification and no measurement resolution effects.

Acceptance of bCa_med for the $\tilde t_1\to b\chi^\pm_1$ scenario with $m_{\chi^\pm_1}=m_{\chi^0_1}+20$ GeV. The acceptance is defined as the fraction of signal events that pass the analysis selection performed on generator-level objects, therefore emulating an ideal detector with perfect particle identification and no measurement resolution effects.

Acceptance of bCa_low for the $\tilde t_1\to b\chi^\pm_1$ scenario with $m_{\chi^\pm_1}=m_{\chi^0_1}+20$ GeV. The acceptance is defined as the fraction of signal events that pass the analysis selection performed on generator-level objects, therefore emulating an ideal detector with perfect particle identification and no measurement resolution effects.

Acceptance of bCb_med1 for the $\tilde t_1\to b\chi^\pm_1$ scenario with $m_{\chi^\pm_1}=m_{\chi^0_1}+20$ GeV. The acceptance is defined as the fraction of signal events that pass the analysis selection performed on generator-level objects, therefore emulating an ideal detector with perfect particle identification and no measurement resolution effects.

Acceptance of bCb_high for the $\tilde t_1\to b\chi^\pm_1$ scenario with $m_{\chi^\pm_1}=m_{\chi^0_1}+20$ GeV. The acceptance is defined as the fraction of signal events that pass the analysis selection performed on generator-level objects, therefore emulating an ideal detector with perfect particle identification and no measurement resolution effects.

Acceptance of 3-body SR ($80<am_{T2}<90$ GeV, $m_T>120$ GeV) for the 3-body scenario ($\tilde t_1\to b W\chi^0_1$). The acceptance is defined as the fraction of signal events that pass the analysis selection performed on generator-level objects, therefore emulating an ideal detector with perfect particle identification and no measurement resolution effects.

Acceptance of tNbC_mix SR for the asymmetric scenario ($\tilde t_1$, $\tilde t_1\to t\chi^0_1$, b $\chi^\pm_1$) with $m_{\chi^\pm_1}=2\times m_{\chi^0_1}$. The acceptance is defined as the fraction of signal events that pass the analysis selection performed on generator-level objects, therefore emulating an ideal detector with perfect particle identification and no measurement resolution effects.

Efficiency of tN_diag SR ($E_T^{miss}>150$ GeV, $m_T>140$ GeV) for the $\tilde t_1\to t\chi^0_1$ scenario with $m_{\tilde t_1}>m_t+m_{\chi^0_1}$. The efficiency is the ratio between the expected signal rate calculated with simulated data passing all the reconstruction level cuts applied to reconstructed objects, and the signal rate for an ideal detector (with perfect particle identification and no measurement resolution effects).

Efficiency of tN_med SR for the $\tilde t_1\to t\chi^0_1$ scenario with $m_{\tilde t_1}>m_t+m_{\chi^0_1}$. The efficiency is the ratio between the expected signal rate calculated with simulated data passing all the reconstruction level cuts applied to reconstructed objects, and the signal rate for an ideal detector (with perfect particle identification and no measurement resolution effects).

Efficiency of tN_boost SR for the $\tilde t_1\to t\chi^0_1$ scenario with $m_{\tilde t_1}>m_t+m_{\chi^0_1}$. The efficiency is the ratio between the expected signal rate calculated with simulated data passing all the reconstruction level cuts applied to reconstructed objects, and the signal rate for an ideal detector (with perfect particle identification and no measurement resolution effects).

Efficiency of bCb_med2 SR ($am_{T2}>250$ GeV, $m_T>60$ GeV) for the $\tilde t_1\to b\chi^\pm_1$ scenario with $m_{\chi^\pm_1}=2\times m_{\chi^0_1}$. The efficiency is the ratio between the expected signal rate calculated with simulated data passing all the reconstruction level cuts applied to reconstructed objects, and the signal rate for an ideal detector (with perfect particle identification and no measurement resolution effects).

Efficiency of bCc_diag SR for the $\tilde t_1\to b\chi^\pm_1$ scenario with $m_{\chi^\pm_1}=2\times m_{\chi^0_1}$. The efficiency is the ratio between the expected signal rate calculated with simulated data passing all the reconstruction level cuts applied to reconstructed objects, and the signal rate for an ideal detector (with perfect particle identification and no measurement resolution effects).

Efficiency of bCd_bulk SR ($am_{T2}>175$ GeV, $m_T>120$ GeV) for the $\tilde t_1\to b\chi^\pm_1$ scenario with $m_{\chi^\pm_1}=2\times m_{\chi^0_1}$. The efficiency is the ratio between the expected signal rate calculated with simulated data passing all the reconstruction level cuts applied to reconstructed objects, and the signal rate for an ideal detector (with perfect particle identification and no measurement resolution effects).

Efficiency of bCd_high1 SR for the $\tilde t_1\to b\chi^\pm_1$ scenario with $m_{\chi^\pm_1}=2\times m_{\chi^0_1}$. The efficiency is the ratio between the expected signal rate calculated with simulated data passing all the reconstruction level cuts applied to reconstructed objects, and the signal rate for an ideal detector (with perfect particle identification and no measurement resolution effects).

Efficiency of bCd_high2 SR for the $\tilde t_1\to b\chi^\pm_1$ scenario with $m_{\chi^\pm_1}=2\times m_{\chi^0_1}$. The efficiency is the ratio between the expected signal rate calculated with simulated data passing all the reconstruction level cuts applied to reconstructed objects, and the signal rate for an ideal detector (with perfect particle identification and no measurement resolution effects).

Efficiency of bCa_med for the $\tilde t_1\to b\chi^\pm_1$ scenario with $m_{\chi^\pm_1}=m_{\chi^0_1}+20$ GeV. The efficiency is the ratio between the expected signal rate calculated with simulated data passing all the reconstruction level cuts applied to reconstructed objects, and the signal rate for an ideal detector (with perfect particle identification and no measurement resolution effects).

Efficiency of bCa_low for the $\tilde t_1\to b\chi^\pm_1$ scenario with $m_{\chi^\pm_1}=m_{\chi^0_1}+20$ GeV. The efficiency is the ratio between the expected signal rate calculated with simulated data passing all the reconstruction level cuts applied to reconstructed objects, and the signal rate for an ideal detector (with perfect particle identification and no measurement resolution effects).

Efficiency of bCb_med1 for the $\tilde t_1\to b\chi^\pm_1$ scenario with $m_{\chi^\pm_1}=m_{\chi^0_1}+20$ GeV. The efficiency is the ratio between the expected signal rate calculated with simulated data passing all the reconstruction level cuts applied to reconstructed objects, and the signal rate for an ideal detector (with perfect particle identification and no measurement resolution effects).

Efficiency of bCb_high for the $\tilde t_1\to b\chi^\pm_1$ scenario with $m_{\chi^\pm_1}=m_{\chi^0_1}+20$ GeV. The efficiency is the ratio between the expected signal rate calculated with simulated data passing all the reconstruction level cuts applied to reconstructed objects, and the signal rate for an ideal detector (with perfect particle identification and no measurement resolution effects).

Efficiency of 3-body SR ($80<am_{T2}<90$ GeV, $m_T>120$ GeV) for the 3-body scenario ($\tilde t_1\to b W\chi^0_1$). The efficiency is the ratio between the expected signal rate calculated with simulated data passing all the reconstruction level cuts applied to reconstructed objects, and the signal rate for an ideal detector (with perfect particle identification and no measurement resolution effects).

Efficiency of tNbC_mix SR for the asymmetric scenario ($\tilde t_1$, $\tilde t_1\to t\chi^0_1$, b $\chi^\pm_1$) with $m_{\chi^\pm_1}=2\times m_{\chi^0_1}$. The efficiency is the ratio between the expected signal rate calculated with simulated data passing all the reconstruction level cuts applied to reconstructed objects, and the signal rate for an ideal detector (with perfect particle identification and no measurement resolution effects).

Number of generated events for the $\tilde t_1\to t\chi^0_1$ scenario with $m_{\tilde t_1}>m_t+m_{\chi^0_1}$.

Number of generated events for the $\tilde t_1\to b\chi^\pm_1$ scenario with $m_{\chi^\pm_1}=2\times m_{\chi^0_1}$.

Number of generated events for the $\tilde t_1\to b\chi^\pm_1$ scenario with $m_{\chi^\pm_1}=m_{\chi^0_1}+20$ GeV; $E_T^{miss}$(gen)$>60$ GeV.

Number of generated events for the $\tilde t_1\to b\chi^\pm_1$ scenario with $m_{\chi^\pm_1}=m_{\chi^0_1}+20$ GeV; $E_T^{miss}$(gen)$>250$ GeV.

Number of generated events for the 3-body scenario ($\tilde t_1\to b W\chi^0_1$).

Number of generated events for the asymmetric scenario ($\tilde t_1$, $\tilde t_1\to t\chi^0_1$, b $\chi^\pm_1$) with $m_{\chi^\pm_1}=2\times m_{\chi^0_1}$.

Cross-section for the $\tilde t_1\to t\chi^0_1$ scenario with $m_{\tilde t_1}>m_t+m_{\chi^0_1}$.

Cross-section for the $\tilde t_1\to b\chi^\pm_1$ scenario with $m_{\chi^\pm_1}=2\times m_{\chi^0_1}$.

Cross-section for the $\tilde t_1\to b\chi^\pm_1$ scenario with $m_{\chi^\pm_1}=m_{\chi^0_1}+20$ GeV.

Cross-section for the 3-body scenario ($\tilde t_1\to b W\chi^0_1$).

Cross-section for the asymmetric scenario ($\tilde t_1$, $\tilde t_1\to t\chi^0_1$, b $\chi^\pm_1$) with $m_{\chi^\pm_1}=2\times m_{\chi^0_1}$.

Combined experimental systematic uncertainty of expected tN_diag SR yields for the $\tilde t_1\to t\chi^0_1$ scenario with $m_{\tilde t_1}>m_t+m_{\chi^0_1}$, using the 2 highest $E_T^{miss}$ and 2 highest $m_T$ bins.

Combined experimental systematic uncertainty of expected tN_med SR yields for the $\tilde t_1\to t\chi^0_1$ scenario with $m_{\tilde t_1}>m_t+m_{\chi^0_1}$.

Combined experimental systematic uncertainty of expected tN_boost SR yields for the $\tilde t_1\to t\chi^0_1$ scenario with $m_{\tilde t_1}>m_t+m_{\chi^0_1}$.

Combined experimental systematic uncertainty of expected bCb_med2 SR yields for the $\tilde t_1\to b\chi^\pm_1$ scenario with $m_{\chi^\pm_1}=2\times m_{\chi^0_1}$, using the 2 highest $am_{T2}$ and 2 highest $m_T$ bins.

Combined experimental systematic uncertainty of expected bCc_diag SR yields for the $\tilde t_1\to b\chi^\pm_1$ scenario with $m_{\chi^\pm_1}=2\times m_{\chi^0_1}$.

Combined experimental systematic uncertainty of expected bCd_bulk SR yields for the $\tilde t_1\to b\chi^\pm_1$ scenario with $m_{\chi^\pm_1}=2\times m_{\chi^0_1}$, using the 2 highest $am_{T2}$ and 2 highest $m_T$ bins.

Combined experimental systematic uncertainty of expected bCd_high1 SR yields for the $\tilde t_1\to b\chi^\pm_1$ scenario with $m_{\chi^\pm_1}=2\times m_{\chi^0_1}$.

Combined experimental systematic uncertainty of expected bCd_high2 SR yields for the $\tilde t_1\to b\chi^\pm_1$ scenario with $m_{\chi^\pm_1}=2\times m_{\chi^0_1}$.

Combined experimental systematic uncertainty of expected bCa_med SR yields for the $\tilde t_1\to b\chi^\pm_1$ scenario with $m_{\chi^\pm_1}=m_{\chi^0_1}+20$ GeV.

Combined experimental systematic uncertainty of expected bCa_low SR yields for the $\tilde t_1\to b\chi^\pm_1$ scenario with $m_{\chi^\pm_1}=m_{\chi^0_1}+20$ GeV.

Combined experimental systematic uncertainty of expected bCb_med1 SR yields for the $\tilde t_1\to b\chi^\pm_1$ scenario with $m_{\chi^\pm_1}=m_{\chi^0_1}+20$ GeV.

Combined experimental systematic uncertainty of expected bCb_high SR yields for the $\tilde t_1\to b\chi^\pm_1$ scenario with $m_{\chi^\pm_1}=m_{\chi^0_1}+20$ GeV.

Combined experimental systematic uncertainty of expected 3-body SR yields for the 3-body scenario ($\tilde t_1\to b W\chi^0_1$), using the 2 lowest $am_{T2}$ and 2 highest $m_T$ bins.

Combined experimental systematic uncertainty of expected tNbC_mix SR yields for the asymmetric scenario ($\tilde t_1$, $\tilde t_1\to t\chi^0_1$, b $\chi^\pm_1$) with $m_{\chi^\pm_1}=2\times m_{\chi^0_1}$.

Observed CLs in tN_diag SR for the $\tilde t_1\to t\chi^0_1$ scenario with $m_{\tilde t_1}>m_t+m_{\chi^0_1}$.

Observed CLs in tN_med SR for the $\tilde t_1\to t\chi^0_1$ scenario with $m_{\tilde t_1}>m_t+m_{\chi^0_1}$.

Observed CLs in tN_boost SR for the $\tilde t_1\to t\chi^0_1$ scenario with $m_{\tilde t_1}>m_t+m_{\chi^0_1}$.

Observed CLs in bCb_med2 SR for the $\tilde t_1\to b\chi^\pm_1$ scenario with $m_{\chi^\pm_1}=2\times m_{\chi^0_1}$.

Observed CLs in bCc_diag SR for the $\tilde t_1\to b\chi^\pm_1$ scenario with $m_{\chi^\pm_1}=2\times m_{\chi^0_1}$.

Observed CLs in bCd_bulk SR for the $\tilde t_1\to b\chi^\pm_1$ scenario with $m_{\chi^\pm_1}=2\times m_{\chi^0_1}$.

Observed CLs in bCd_high1 SR for the $\tilde t_1\to b\chi^\pm_1$ scenario with $m_{\chi^\pm_1}=2\times m_{\chi^0_1}$.

Observed CLs in bCd_high2 SR for the $\tilde t_1\to b\chi^\pm_1$ scenario with $m_{\chi^\pm_1}=2\times m_{\chi^0_1}$.

Observed CLs in bCa_med SR for the $\tilde t_1\to b\chi^\pm_1$ scenario with $m_{\chi^\pm_1}=m_{\chi^0_1}+20$ GeV.

Observed CLs in bCa_low SR for the $\tilde t_1\to b\chi^\pm_1$ scenario with $m_{\chi^\pm_1}=m_{\chi^0_1}+20$ GeV.

Observed CLs in bCb_med1 SR for the $\tilde t_1\to b\chi^\pm_1$ scenario with $m_{\chi^\pm_1}=m_{\chi^0_1}+20$ GeV.

Observed CLs in bCb_high SR for the $\tilde t_1\to b\chi^\pm_1$ scenario with $m_{\chi^\pm_1}=m_{\chi^0_1}+20$ GeV.

Observed CLs in 3-body SR for the 3-body scenario ($\tilde t_1\to b W\chi^0_1$).

Observed CLs in tNbC_mix SR for the mixed scenario (50% $\tilde t_1\to t\chi^0_1$, 50% $\tilde t_1\to b\chi^0_1$).

Expected CLs in tN_diag SR for the $\tilde t_1\to t\chi^0_1$ scenario with $m_{\tilde t_1}>m_t+m_{\chi^0_1}$.

Expected CLs in tN_med SR for the $\tilde t_1\to t\chi^0_1$ scenario with $m_{\tilde t_1}>m_t+m_{\chi^0_1}$.

Expected CLs in tN_boost SR for the $\tilde t_1\to t\chi^0_1$ scenario with $m_{\tilde t_1}>m_t+m_{\chi^0_1}$.

Expected CLs in bCb_med2 SR for the $\tilde t_1\to b\chi^\pm_1$ scenario with $m_{\chi^\pm_1}=2\times m_{\chi^0_1}$.

Expected CLs in bCc_diag SR for the $\tilde t_1\to b\chi^\pm_1$ scenario with $m_{\chi^\pm_1}=2\times m_{\chi^0_1}$.

Expected CLs in bCd_bulk SR for the $\tilde t_1\to b\chi^\pm_1$ scenario with $m_{\chi^\pm_1}=2\times m_{\chi^0_1}$.

Expected CLs in bCd_high1 SR for the $\tilde t_1\to b\chi^\pm_1$ scenario with $m_{\chi^\pm_1}=2\times m_{\chi^0_1}$.

Expected CLs in bCd_high2 SR for the $\tilde t_1\to b\chi^\pm_1$ scenario with $m_{\chi^\pm_1}=2\times m_{\chi^0_1}$.

Expected CLs in bCa_med SR for the $\tilde t_1\to b\chi^\pm_1$ scenario with $m_{\chi^\pm_1}=m_{\chi^0_1}+20$ GeV.

Expected CLs in bCa_low SR for the $\tilde t_1\to b\chi^\pm_1$ scenario with $m_{\chi^\pm_1}=m_{\chi^0_1}+20$ GeV.

Expected CLs in bCb_med1 SR for the $\tilde t_1\to b\chi^\pm_1$ scenario with $m_{\chi^\pm_1}=m_{\chi^0_1}+20$ GeV.

Expected CLs in bCb_high SR for the $\tilde t_1\to b\chi^\pm_1$ scenario with $m_{\chi^\pm_1}=m_{\chi^0_1}+20$ GeV.

Expected CLs in 3-body SR for the 3-body scenario ($\tilde t_1\to b W\chi^0_1$).

Expected CLs in tNbC_mix SR for the mixed scenario (50% $\tilde t_1\to t\chi^0_1$, 50% $\tilde t_1\to b\chi^\pm_1$).

The effective mass distribution in VR2Z.

The ETmiss distribution in SR0noZa.

The effective mass distribution in SR0noZa.

The ETmiss distribution in SR1noZa.

The effective mass distribution in SR1noZa.

The ETmiss distribution in SR2noZa.

The effective mass distribution in SR2noZa.

The ETmiss distribution in SR0noZb.

The effective mass distribution in SR0noZb.

The ETmiss distribution in SR1noZb.

The effective mass distribution in SR1noZb.

The ETmiss distribution in SR2noZb.

The effective mass distribution in SR2noZb.

The ETmiss distribution in SR0Z.

The effective mass distribution in SR0Z.

The ETmiss distribution in SR1Z.

The effective mass distribution in SR1Z.

The ETmiss distribution in SR2Z.

The effective mass distribution in SR2Z.

Observed 95% CL exclusion contour for the RPV chargino NLSP model with lambda_121 != 0.

Expected 95% CL exclusion contour for the RPV chargino NLSP model with lambda_121 != 0.

Observed 95% CL exclusion contour for the RPV chargino NLSP model with lambda_122 != 0.

Expected 95% CL exclusion contour for the RPV chargino NLSP model with lambda_122 != 0.

Observed 95% CL exclusion contour for the RPV chargino NLSP model with lambda_133 != 0.

Expected 95% CL exclusion contour for the RPV chargino NLSP model with lambda_133 != 0.

Observed 95% CL exclusion contour for the RPV chargino NLSP model with lambda_233 != 0.

Expected 95% CL exclusion contour for the RPV chargino NLSP model with lambda_233 != 0.

Observed 95% CL exclusion contour for the RPV gluino NLSP model with lambda_121 != 0.

Expected 95% CL exclusion contour for the RPV gluino NLSP model with lambda_121 != 0.

Observed 95% CL exclusion contour for the RPV gluino NLSP model with lambda_122 != 0.

Expected 95% CL exclusion contour for the RPV gluino NLSP model with lambda_122 != 0.

Observed 95% CL exclusion contour for the RPV gluino NLSP model with lambda_133 != 0.

Expected 95% CL exclusion contour for the RPV gluino NLSP model with lambda_133 != 0.

Observed 95% CL exclusion contour for the RPV gluino NLSP model with lambda_233 != 0.

Expected 95% CL exclusion contour for the RPV gluino NLSP model with lambda_233 != 0.

Observed 95% CL exclusion contour for the RPV Lslepton NLSP model with lambda_121 != 0.

Expected 95% CL exclusion contour for the RPV Lslepton NLSP model with lambda_121 != 0.

Observed 95% CL exclusion contour for the RPV Lslepton NLSP model with lambda_122 != 0.

Expected 95% CL exclusion contour for the RPV Lslepton NLSP model with lambda_122 != 0.

Observed 95% CL exclusion contour for the RPV Lslepton NLSP model with lambda_133 != 0.

Expected 95% CL exclusion contour for the RPV Lslepton NLSP model with lambda_133 != 0.

Observed 95% CL exclusion contour for the RPV Lslepton NLSP model with lambda_233 != 0.

Expected 95% CL exclusion contour for the RPV Lslepton NLSP model with lambda_233 != 0.

Observed 95% CL exclusion contour for the RPV Rslepton NLSP model with lambda_121 != 0.

Expected 95% CL exclusion contour for the RPV Rslepton NLSP model with lambda_121 != 0.

Observed 95% CL exclusion contour for the RPV Rslepton NLSP model with lambda_122 != 0.

Expected 95% CL exclusion contour for the RPV Rslepton NLSP model with lambda_122 != 0.

Observed 95% CL exclusion contour for the RPV Rslepton NLSP model with lambda_133 != 0.

Expected 95% CL exclusion contour for the RPV Rslepton NLSP model with lambda_133 != 0.

Observed 95% CL exclusion contour for the RPV Rslepton NLSP model with lambda_233 != 0.

Expected 95% CL exclusion contour for the RPV Rslepton NLSP model with lambda_233 != 0.

Observed 95% CL exclusion contour for the RPV sneutrino NLSP model with lambda_121 != 0.

Expected 95% CL exclusion contour for the RPV sneutrino NLSP model with lambda_121 != 0.

Observed 95% CL exclusion contour for the RPV sneutrino NLSP model with lambda_122 != 0.

Expected 95% CL exclusion contour for the RPV sneutrino NLSP model with lambda_122 != 0.

Observed 95% CL exclusion contour for the RPV sneutrino NLSP model with lambda_133 != 0.

Expected 95% CL exclusion contour for the RPV sneutrino NLSP model with lambda_133 != 0.

Observed 95% CL exclusion contour for the RPV sneutrino NLSP model with lambda_233 != 0.

Expected 95% CL exclusion contour for the RPV sneutrino NLSP model with lambda_233 != 0.

Observed 95% CL exclusion contour for the R-slepton RPC model.

Expected 95% CL exclusion contour for the R-slepton RPC model.

Observed and expected 95% CL cross-section upper limits for the Stau RPC model, together with the theoretically predicted cross-section.

Observed and expected 95% CL cross-section upper limits for the Z RPC model, together with the theoretically predicted cross-section.

Observed 95% CL exclusion contour for the GGM tan beta = 1.5 model.

Expected 95% CL exclusion contour for the GGM tan beta = 1.5 model.

Observed 95% CL exclusion contour for the GGM tan beta = 30 model.

Expected 95% CL exclusion contour for the GGM tan beta = 30 model.

Observed 95% CL cross-section upper limit for the RPV chargino NLSP models with lambda_121 != 0 and lambda_122 != 0, and the selection of Z-veto signal regions used to set limits in these models. The combination of regions used is ordered by the minimum number of hadronic taus required. For example, ``bba' means that the regions SR0noZb, SR1noZb and SR2noZa were used, in addition to the three Z-rich regions (SR0-2Z).

Observed 95% CL cross-section upper limit for the RPV chargino NLSP models with lambda_133 != 0 and lambda_233 != 0, and the selection of Z-veto signal regions used to set limits in these models. The combination of regions used is ordered by the minimum number of hadronic taus required. For example, ``bba' means that the regions SR0noZb, SR1noZb and SR2noZa were used, in addition to the three Z-rich regions (SR0-2Z).

Observed 95% CL cross-section upper limit for the RPV gluino NLSP models with lambda_121 != 0 and lambda_122 != 0, and the selection of Z-veto signal regions used to set limits in these models. The combination of regions used is ordered by the minimum number of hadronic taus required. For example, ``bba' means that the regions SR0noZb, SR1noZb and SR2noZa were used, in addition to the three Z-rich regions (SR0-2Z).

Observed 95% CL cross-section upper limit for the RPV gluino NLSP models with lambda_133 != 0 and lambda_233 != 0, and the selection of Z-veto signal regions used to set limits in these models. The combination of regions used is ordered by the minimum number of hadronic taus required. For example, ``bba' means that the regions SR0noZb, SR1noZb and SR2noZa were used, in addition to the three Z-rich regions (SR0-2Z).

Observed 95% CL cross-section upper limit for the RPV Lslepton NLSP models with lambda_121 != 0 and lambda_122 != 0, and the selection of Z-veto signal regions used to set limits in these models. The combination of regions used is ordered by the minimum number of hadronic taus required. For example, ``bba' means that the regions SR0noZb, SR1noZb and SR2noZa were used, in addition to the three Z-rich regions (SR0-2Z).

Observed 95% CL cross-section upper limit for the RPV Lslepton NLSP models with lambda_133 != 0 and lambda_233 != 0, and the selection of Z-veto signal regions used to set limits in these models. The combination of regions used is ordered by the minimum number of hadronic taus required. For example, ``bba' means that the regions SR0noZb, SR1noZb and SR2noZa were used, in addition to the three Z-rich regions (SR0-2Z).

Observed 95% CL cross-section upper limit for the RPV Rslepton NLSP models with lambda_121 != 0 and lambda_122 != 0, and the selection of Z-veto signal regions used to set limits in these models. The combination of regions used is ordered by the minimum number of hadronic taus required. For example, ``bba' means that the regions SR0noZb, SR1noZb and SR2noZa were used, in addition to the three Z-rich regions (SR0-2Z).

Observed 95% CL cross-section upper limit for the RPV Rslepton NLSP models with lambda_133 != 0 and lambda_233 != 0, and the selection of Z-veto signal regions used to set limits in these models. The combination of regions used is ordered by the minimum number of hadronic taus required. For example, ``bba' means that the regions SR0noZb, SR1noZb and SR2noZa were used, in addition to the three Z-rich regions (SR0-2Z).

Observed 95% CL cross-section upper limit for the RPV sneutrino NLSP models with lambda_121 != 0 and lambda_122 != 0, and the selection of Z-veto signal regions used to set limits in these models. The combination of regions used is ordered by the minimum number of hadronic taus required. For example, ``bba' means that the regions SR0noZb, SR1noZb and SR2noZa were used, in addition to the three Z-rich regions (SR0-2Z).

Observed 95% CL cross-section upper limit for the RPV sneutrino NLSP models with lambda_133 != 0 and lambda_233 != 0, and the selection of Z-veto signal regions used to set limits in these models. The combination of regions used is ordered by the minimum number of hadronic taus required. For example, ``bba' means that the regions SR0noZb, SR1noZb and SR2noZa were used, in addition to the three Z-rich regions (SR0-2Z).

Observed 95% CL cross-section upper limit for the R-slepton RPC model, and the selection of Z-veto signal regions used to set limits in this model. The combination of regions used is ordered by the minimum number of hadronic taus required. For example, ``bbb' means that the regions SR0noZb, SR1noZb and SR2noZb were used, in addition to the three Z-rich regions (SR0-2Z). For the RPC stau and Z models, the ``aaa' combination of regions was used throughout.

Performance of the SR0noZa selection in the R-slepton RPC model: number of generated signal events; total signal cross-section; acceptance; efficiency; total experimental systematic uncertainty, not including Monte Carlo statistics; observed CL using this region alone; expected CL using this region alone.

Performance of the SR0noZb selection in the RPV chargino NLSP model with lambda_121 != 0: number of generated signal events; total signal cross-section; acceptance; efficiency; total experimental systematic uncertainty, not including Monte Carlo statistics; observed CL using this region alone; expected CL using this region alone.

Performance of the SR1noZa selection in the RPV sneutrino NLSP model with lambda_233 != 0: number of generated signal events; total signal cross-section; acceptance; efficiency; total experimental systematic uncertainty, not including Monte Carlo statistics; observed CL using this region alone; expected CL using this region alone.

Performance of the SR1noZb selection in the RPV gluino NLSP model with lambda_133 != 0: number of generated signal events; total signal cross-section; acceptance; efficiency; total experimental systematic uncertainty, not including Monte Carlo statistics; observed CL using this region alone; expected CL using this region alone.

Performance of the SR2noZa selection in the RPV sneutrino NLSP model with lambda_233 != 0: number of generated signal events; total signal cross-section; acceptance; efficiency; total experimental systematic uncertainty, not including Monte Carlo statistics; observed CL using this region alone; expected CL using this region alone.

Performance of the SR2noZb selection in the RPV gluino NLSP model with lambda_133 != 0: number of generated signal events; total signal cross-section; acceptance; efficiency; total experimental systematic uncertainty, not including Monte Carlo statistics; observed CL using this region alone; expected CL using this region alone.

Performance of the SR0Z selection in the GGM tan beta = 30 model: number of generated signal events; total signal cross-section; acceptance; efficiency; total experimental systematic uncertainty, not including Monte Carlo statistics; observed CL using this region alone; expected CL using this region alone.

Cut flows for a representative selection of SUSY signal points in the Z-veto signal regions. In each case, m2 and m1 refer to the axes of the plots in Sec. XI, where m2 is the larger of the two masses. The number of events expected for a luminosity of 20.3 fb-1 is quoted at each step of the selection. The preselection requires four baseline leptons, at least two of which are light leptons; the signal lepton selection is made at the ``Lepton Multiplicity' stage. ``Event Cleaning' refers to the selection criteria applied to remove non-collision backgrounds and detector noise.

Cut flows for a representative selection of SUSY signal points in the Z-rich signal regions. In each case, m2 and m1 refer to the axes of the plots in Sec. XI, where m2 is the larger of the two masses (or the value of mu in the case of GGM models). The number of events expected for a luminosity of 20.3 fb-1 is quoted at each step of the selection. The preselection requires four baseline leptons, at least two of which are light leptons; the signal lepton selection is made at the ``Lepton Multiplicity' stage. ``Event Cleaning' refers to the selection criteria applied to remove non-collision backgrounds and detector noise.

Cut flows by lepton channel for a representative selection of SUSY signal points in the SR0noZa signal region. In each case, m2 and m1 refer to the axes of the plots in Sec. XI, where m2 is the larger of the two masses. The number of events expected for a luminosity of 20.3 fb-1 is quoted at each step of the selection. The preselection requires four baseline leptons, at least two of which are light leptons; the signal lepton selection is made at the ``Lepton Multiplicity' stage. ``Event Cleaning' refers to the selection criteria applied to remove non-collision backgrounds and detector noise. The RPC R-slepton model is used, with (m2,m1) = (450,300) GeV.

Cut flows by lepton channel for a representative selection of SUSY signal points in the SR1noZb signal region. In each case, m2 and m1 refer to the axes of the plots in Sec. XI, where m2 is the larger of the two masses. The number of events expected for a luminosity of 20.3 fb-1 is quoted at each step of the selection. The preselection requires four baseline leptons, at least two of which are light leptons; the signal lepton selection is made at the ``Lepton Multiplicity' stage. ``Event Cleaning' refers to the selection criteria applied to remove non-collision backgrounds and detector noise. The RPV gluino NLSP model is used, with lambda_133 != 0 and (m2,m1) = (800,400) GeV.

Cut flows by lepton channel for a representative selection of SUSY signal points in the SR0Z signal region. In each case, m2 and m1 refer to the axes of the plots in Sec. XI, where m2 is the value of mu. The number of events expected for a luminosity of 20.3 fb-1 is quoted at each step of the selection. The preselection requires four baseline leptons, at least two of which are light leptons; the signal lepton selection is made at the ``Lepton Multiplicity' stage. ``Event Cleaning' refers to the selection criteria applied to remove non-collision backgrounds and detector noise. The GGM tan beta = 30 model is used, with (m2,m1) = (200,1000) GeV.

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

Observed 95% CL limit for the pMSSM grid.

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

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

Expected 95% CL limit for the pMSSM grid.

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

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

Observed 95% CL limit for the 2Step grid.

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

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

Expected 95% CL limit for the 2Step grid.

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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