The expected limit contour in the GLUINO-NEUTRALINO plane.

The observed limit contour in the SQUARK-NEUTRALINO plane.

The expected limit contour in the SQUARK-NEUTRALINO plane.

The stop pair production cross section and its uncertainty due to PDF, factorization and the renormalization scales.

(1) Number of generated MC events for the spin 1/2 top partner signal grid (2) Relative Cross section uncertainties for the spin 1/2 top partner signal grid.

(1) Acceptance of the same flavour selection for the spin 1/2 top partner signal grid (2) Selection efficiency of the same flavour selection for the spin 1/2 top partner signal grid (3) Product of the acceptance and efficiency of the same flavour selection for the spin 1/2 top partner signal grid (4) Relative experimental uncertainties on the acceptance times efficiency of the same flavour selection for the spin 1/2 top partner signal grid.

(1) Acceptance of the different flavour selection for the spin 1/2 top partner signal grid (2) Selection efficiency of the different flavour selection for the spin 1/2 top partner signal grid (3) Product of the acceptance and efficiency of the different flavour selection for the spin 1/2 top partner signal grid (4) Relative experimental uncertainties on the acceptance times efficiency of the different flavour selection for the spin 1/2 top partner signal grid.

(1) Observed CLs values for the scalar top signal grid (2) Expected CLs values for the scalar top signal grid.

(1) Observed CLs values for the spin 1/2 top partner signal grid (2) Expected CLs values for the spin 1/2 top partner signal grid.

Cross section limits [pb] for the scalar top signal grid.

Cross section limits [pb] for the spin 1/2 top partner signal grid.

Observed 95% CL limit for stop grid as a function of the scalar top and neutralino masses.

Observed 95% CL limit for stop grid as a function of the scalar top and neutralino masses, varying signal cross section of +1sigma.

Observed 95% CL limit for stop grid as a function of the scalar top and neutralino masses, varying signal cross section of -1sigma.

Expected 95% CL limit for stop grid as a function of the scalar top and neutralino masses.

Expected 95% CL limit for stop grid as a function of the scalar top and neutralino masses, varying the uncertainty of +1sigma.

Expected 95% CL limit for stop grid as a function of the scalar top and neutralino masses, varying the uncertainty of -1sigma.

Observed 95% CL limit for top partner grid as a function of the top partner and neutralino masses.

Observed 95% CL limit for top partner grid as a function of the top partner and neutralino masses, varying signal cross section of +1sigma.

Observed 95% CL limit for top partner grid as a function of the top partner and neutralino masses, varying signal cross section of -1sigma.

Expected 95% CL limit for top partner grid as a function of the top partner and neutralino masses.

Expected 95% CL limit for top partner grid as a function of the top partner and neutralino masses, varying the uncertainty of +1sigma.

Expected 95% CL limit for top partner grid as a function of the top partner and neutralino masses, varying the uncertainty of -1sigma.

Transverse momentum distribution for the second leading lepton for events in the SR2 signal region for DATA and SM predictions.

Transverse momentum distribution for the third leading lepton for events in the SR1 signal region for DATA and SM predictions.

Transverse momentum distribution for the third leading lepton for events in the SR2 signal region for DATA and SM predictions.

Missing transverse energy for events in the SR1 signal region for DATA and SM predictions.

Missing transverse energy for events in the SR2 signal region for DATA and SM predictions.

Invariant mass of the same-flavour-opposite-sign (SFOS) lepton pair for events in the SR1 signal region for DATA and SM predictions.

The Cross Section in signal region SR1 for the SUSY pMSSM model with M1=100 GeV grid.

The Cross Section in signal region SR1 for the SUSY simplified model grid.

The Number of generated Events in signal region SR1 for the SUSY pMSSM model with M1=100 GeV grid.

The Number of generated Events in signal region SR1 for the SUSY simplified model grid.

The Efficiency in signal region SR1 for the SUSY pMSSM model with M1=100 GeV grid.

The Efficiency in signal region SR1 for the SUSY simplified model grid.

The Acceptance in signal region SR1 for the SUSY pMSSM model with M1=100 GeV grid.

The Acceptance in signal region SR1 for the SUSY simplified model grid.

The Acceptance*Efficiency in signal region SR1 for the SUSY pMSSM model with M1=100 GeV grid.

The Acceptance*Efficiency in signal region SR1 for the SUSY simplified model grid.

The Systematic Uncertainty of the data (excluding the Monte Carlo) in signal region SR1 for the SUSY pMSSM model with M1=100 GeV grid.

The Systematic Uncertainty of the data (excluding the Monte Carlo) in signal region SR1 for the SUSY simplified model grid.

CL values for the pMSSM with M1=100 GeV model grid for the SR1 signal region.

CL values for the simplified model model grid for the SR1 signal region.

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.

Parametrization of the $A_{X}$ factor, defined as the fraction of diphoton resonances satisfying the fiducial acceptance at the particle level, as function of $m_{X}$ obtained from the narrow width simulated signal samples produced in gluon fusion.

The correction factor, $C_{X}$, defined as the ratio between the number of reconstructed signal events passing the analysis cuts and the number of signal events at the particle level generated within the fiducial volume, and acceptance correction factor, $A_{X}$, defined as the fraction of diphoton resonances satisfying the fiducial acceptance at the particle level. Both are computed for NWA spin-0 models as a function of $m_{X}$.

Effect of event selections on a scalar MC signal sample generated for $m_{X}$ = 15 GeV and on the data. For the MC sample, the efficiencies are shown after applying event weights and a truth level filter that requires two photons with $p^{\gamma\gamma}_{T}>40$ GeV; for the data, the absolute yields are shown. The initial yields for data include a trigger preselection that is the OR of a list of single photon and diphoton triggers. The "2 $loose$ photons" step includes the kinematic acceptance cuts.

Parameterization of the Double Sided Crystal Ball function parameters describing the scalar mass resolution model as a function of $m_{X}$ [GeV].

Numbers of observed and background events for SR3L low for each bin of the distribution in Meff. The table corresponds to Fig. 4(d). The statistical and systematic uncertainties are combined for the predicted numbers.

Numbers of observed and background events for SR3L high for each bin of the distribution in Meff. The table corresponds to Fig. 4(e). The statistical and systematic uncertainties are combined for the predicted numbers.

The efficiencies are calculated for all simplified models. For each model, the values are given for the five signal regions and their combination. The simplified model is for direct production of squarks that decay into two steps into q q W Z W Z chi1^0 chi1^0 (see Fig. 6c in the paper).

The efficiencies are calculated for all simplified extra dimension model (see Fig. 8d in the paper). For each model, the values are given for the five signal regions and their combination.

The efficiencies are calculated for all simplified models. For each model, the values are given for the five signal regions and their combination. The simplified model is for direct production of gluinos that decay via sleptons into q q q q l l (l l) chi1^0 chi1^0 + neutrinos (see Fig. 6d in the paper).

The efficiencies are calculated for all simplified models. For each model, the values are given for the five signal regions and their combination. The simplified model is for direct production of gluinos that decay into q q q q W W chi1^0 chi1^0 (see Fig. 6a in the paper).

The efficiencies are calculated for all simplified models. For each model, the values are given for the five signal regions and their combination. The simplified model is for direct production of gluinos that decay into t tbar t tbar chi1^0 chi1^0 (see Fig. 5a in the paper). This particular model assumes that top quark is much heavier than gluino.

The efficiencies are calculated for all mSUGRA models (see Fig. 8a in the paper). For each model, the values are given for the five signal regions and their combination. The model assumes tan(beta)=30, A0=2m0, and mu>0.

The efficiencies are calculated for all simplified models. For each model, the values are given for the five signal regions and their combination. The simplified model is for direct production of gluinos. A gluino decays into t c chi1^0 (see Fig. 5c in the paper). This particular model assumes that m(chi1^0) = m(stop) - 20 GeV.

The efficiencies are calculated for all GMSB models (see Fig. 8c in the paper). For each model, the values are given for the five signal regions and their combination. The model assumes mmess=250 TeV, m5=3, mu>0, and Cgrav=1.

The efficiencies are calculated for all simplified models. For each model, the values are given for the five signal regions and their combination. The simplified model is for direct production of bottom squarks. A bottom squark decays into t chi1^(+-) and chi1^(+-) --> W^(+-) chi1^0 (see Fig. 7a in the paper). This particular model assumes that m(chi1^0)=60 GeV.

The efficiencies are calculated for all simplified models. For each model, the values are given for the five signal regions and their combination. The simplified model is for direct production of gluinos and top squarks. Top squarks undergo R-parity violating decays into b s and gluinos decay into t stop (see Fig. 5d in the paper).

The efficiencies are calculated for all simplified models. For each model, the values are given for the five signal regions and their combination. The simplified model is for direct production of bottom squarks. A bottom squark decays into t chi1^(+-) and chi1^(+-) --> W^(+-) chi1^0 (see Fig. 7b in the paper). This particular model assumes that m(chi1^0)=2(chi1^0).

The efficiencies are calculated for all mSUGRA/CMSSM models with bRPV (see Fig. 8b in the paper). For each model, the values are given for the five signal regions and their combination. The model assumes tan(beta)=30, A0=2m0, mu>0, and bRPV.

The efficiencies are calculated for all simplified models. For each model, the values are given for the five signal regions and their combination. The simplified model is for direct production of squarks. Squarks decay into q q l l (l l) chi1^0 chi1^0 + neutrinos (see Fig. 6e in the paper).

The efficiencies are calculated for all simplified models. For each model, the values are given for the five signal regions and their combination. The simplified model is for direct pair-production of gluinos that decay via a two-step process into q q q q W Z W Z chi1^0 chi1^0 (see Fig. 6b in the paper).

The efficiencies are calculated for all simplified models. For each model, the values are given for the five signal regions and their combination. The simplified model is for direct pair production of gluinos. A gluino decays into t stop. Consequently, a top squark squark decays into b chi1^(+-) and chi1^(+-) --> W^(+-) chi1^0 (see Fig. 5b in the paper). This particular model assumes that m(stop) < m(gluino), m(chi1^0)=6 GeV, and m(chi1^(+-))=118 GeV.

The acceptances (in percent, %) are calculated for all simplified models. For each model, the values are given for the five signal regions and their combination. The simplified model is for direct production of squarks that decay into two steps into q q W Z W Z chi1^0 chi1^0 (see Fig. 6c in the paper).

The acceptances (in percent, %) are calculated for all simplified extra dimension model (see Fig. 8d in the paper). For each model, the values are given for the five signal regions and their combination.

The acceptances (in percent, %) are calculated for all simplified models. For each model, the values are given for the five signal regions and their combination. The simplified model is for direct production of gluinos that decay via sleptons into q q q q l l (l l) chi1^0 chi1^0 + neutrinos (see Fig. 6d in the paper).

The acceptances (in percent, %) are calculated for all simplified models. For each model, the values are given for the five signal regions and their combination. The simplified model is for direct production of gluinos that decay into q q q q W W chi1^0 chi1^0 (see Fig. 6a in the paper).

The acceptances (in percent, %) are calculated for all simplified models. For each model, the values are given for the five signal regions and their combination. The simplified model is for direct production of gluinos that decay into t tbar t tbar chi1^0 chi1^0 (see Fig. 5a in the paper). This particular model assumes that top quark is much heavier than gluino.

The acceptances (in percent, %) are calculated for all mSUGRA models (see Fig. 8a in the paper). For each model, the values are given for the five signal regions and their combination. The model assumes tan(beta)=30, A0=2m0, and mu>0.

The acceptances (in percent, %) are calculated for all simplified models. For each model, the values are given for the five signal regions and their combination. The simplified model is for direct production of gluinos. A gluino decays into t c chi1^0 (see Fig. 5c in the paper). This particular model assumes that m(chi1^0) = m(stop) - 20 GeV.

The acceptances (in percent, %) are calculated for all GMSB models (see Fig. 8c in the paper). For each model, the values are given for the five signal regions and their combination. The model assumes mmess=250 TeV, m5=3, mu>0, and Cgrav=1.

The acceptances (in percent, %) are calculated for all simplified models. For each model, the values are given for the five signal regions and their combination. The simplified model is for direct production of bottom squarks. A bottom squark decays into t chi1^(+-) and chi1^(+-) --> W^(+-) chi1^0 (see Fig. 7a in the paper). This particular model assumes that m(chi1^0)=60 GeV.

The acceptances (in percent, %) are calculated for all simplified models. For each model, the values are given for the five signal regions and their combination. The simplified model is for direct production of gluinos and top squarks. Top squarks undergo R-parity violating decays into bs and gluinos decay into t stop (see Fig. 5d in the paper).

The acceptances (in percent, %) are calculated for all simplified models. For each model, the values are given for the five signal regions and their combination. The simplified model is for direct production of bottom squarks. A bottom squark decays into t chi1^(+-) and chi1^(+-) --> W chi1^0 (see Fig. 7b in the paper). This particular model assumes that m(chi1^0)=2(chi1^0).

The acceptances (in percent, %) are calculated for all mSUGRA/CMSSM models with bRPV (see Fig. 8b in the paper). For each model, the values are given for the five signal regions and their combination. The model assumes tan(beta)=30, A0=2m0, mu>0, and bRPV.

The acceptances (in percent, %) are calculated for all simplified models. For each model, the values are given for the five signal regions and their combination. The simplified model is for direct production of squarks. Squarks decay into q q l l (l l) chi1^0 chi1^0 + neutrinos (see Fig. 6e in the paper).

The acceptances (in percent, %) are calculated for all simplified models. For each model, the values are given for the five signal regions and their combination. The simplified model is for direct pair-production of gluinos that decay via a two-step process into q q q q W Z W Z chi1^0 chi1^0 (see Fig. 6b in the paper).

The acceptances (in percent, %) are calculated for all simplified models. For each model, the values are given for the five signal regions and their combination. The simplified model is for direct pair production of gluinos. A gluino decays into t stop. Consequently, a top squark squark decays into b chi1^(+-) and chi1^(+-) --> W^(+-) chi1^0 (see Fig. 5b in the paper). This particular model assumes that m(stop) < m(gluino), m(chi1^0)=6 GeV, and m(chi1^(+-))=118 GeV.

The limits on observed cross section are calculated for all simplified models. The simplified models are for direct pair production of squarks that decay into two steps into q q W Z W Z chi1^0 chi1^0 (see Fig. 6c in the paper).

The limits on observed cross sections are calculated for all simplified models. The simplified models are for direct pair-production of gluinos that decay via sleptons into q q q q l l (l l) chi1^0 chi1^0 + neutrinos (see Fig. 6d in the paper).

The limits on observed cross sections are calculated for all simplified models. The simplified models are for direct production of gluinos that decay into q q q q W W chi1^0 chi1^0 (see Fig. 6a in the paper).

The limits on observed cross sections are calculated for all simplified models. The simplified models are for direct production of gluinos that decay into t tbar t tbar chi1^0 chi1^0 (see Fig. 5a in the paper). This particular model assumes that top quark is much heavier than gluino.

The limits on observed cross sections are calculated for all simplified models. The simplified models are for direct pair production of gluinos. A gluino decays into t c chi1^0 (see Fig. 5c in the paper). This particular model assumes that m(chi1^0) = m(stop) - 20 GeV.

The limits on observed cross sections are calculated for all simplified models. The simplified models are for direct production of bottom squarks. A bottom squark decays into t chi1^(+-) and chi1^(+-) --> W^(+-) chi1^0 (see Fig. 7a in the paper). This particular model assumes that m(chi1^0)=60 GeV.

The limits on observed cross sections are calculated for all simplified models. The simplified models are for direct production of gluinos and top squarks. Top squarks undergo R-parity violating decays into bs and gluinos decay into t stop (see Fig. 5d in the paper).

The limits on observed cross sections are calculated for all simplified models. The simplified models are for direct production of bottom squarks. A bottom squark decays into t chi1^(+-) and chi1^(+-) --> W^(+-) chi1^0 (see Fig. 7b in the paper). This particular model assumes that m(chi1^0)=2(chi1^0).

The limits on observed cross sections are calculated for all simplified models. The simplified models are for direct production of squarks. Squarks decay into q q l l (l l) chi1^0 chi1^0 + neutrinos (see Fig. 6e in the paper).

The limits on observed cross sections are calculated for all simplified models. The simplified models are for direct pair-production of gluinos that decay via a two-step process into q q q q W Z W Z chi1^0 chi1^0 (see Fig. 6b in the paper).

The limits on observed cross sections are calculated for all simplified models. The simplified models are for direct pair production of gluinos. A gluino decays into t stop. Consequently, a top squark squark decays into b chi1^(+-) and chi1^(+-) --> W^(+-) chi1^0 (see Fig. 5b in the paper). This particular model assumes that m(stop) < m(gluino), m(chi1^0)=6 GeV, and m(chi1^(+-))=118 GeV.

The signal event yields are calculated for all simplified models. For each model, the values are given for the five signal regions and their combination. The simplified model is for direct production of squarks that decay into two steps into q q W Z W Z chi1^0 chi1^0 (see Fig. 6c in the paper).

The signal event yields are calculated for all simplified extra dimension model (see Fig. 8d in the paper). For each model, the values are given for the five signal regions and their combination.

The signal event yields are calculated for all simplified models. For each model, the values are given for the five signal regions and their combination. The simplified model is for direct production of gluinos that decay via sleptons into q q q q l l (l l) chi1^0 chi1^0 + neutrinos (see Fig. 6d in the paper).

The signal event yields are calculated for all simplified models. For each model, the values are given for the five signal regions and their combination. The simplified model is for direct production of gluinos that decay into q q q q W W chi1^0 chi1^0 (see Fig. 6a in the paper).

The signal event yields are calculated for all simplified models. For each model, the values are given for the five signal regions and their combination. The simplified model is for direct production of gluinos that decay into t tbar t tbar chi1^0 chi1^0 (see Fig. 5a in the paper). This particular model assumes that top quark is much heavier than gluino.

The signal event yields are calculated for all mSUGRA models (see Fig. 8a in the paper). For each model, the values are given for the five signal regions and their combination. The model assumes tan(beta)=30, A0=2m0, and mu>0.

The signal event yields are calculated for all simplified models. For each model, the values are given for the five signal regions and their combination. The simplified model is for direct production of gluinos. A gluino decays into t c chi1^0 (see Fig. 5c in the paper). This particular model assumes that m(chi1^0) = m(stop)-20 GeV.

The signal event yields are calculated for all GMSB models (see Fig. 8c in the paper). For each model, the values are given for the five signal regions and their combination. The model assumes mmess=250 TeV, m5=3, mu>0, and Cgrav=1.

The signal event yields are calculated for all simplified models. For each model, the values are given for the five signal regions and their combination. The simplified model is for direct production of bottom squarks. A bottom squark decays into t chi1^(+-) and chi1^(+-) --> W^(+-) chi1^0 (see Fig. 7a in the paper). This particular model assumes that m(chi1^0)=60 GeV.

The signal event yields are calculated for all simplified models. For each model, the values are given for the five signal regions and their combination. The simplified model is for direct production of gluinos and top squarks. Top squarks undergo R-parity violating decays into bs and gluinos decay into t stop (see Fig. 5d in the paper).

The signal event yields are calculated for all simplified models. For each model, the values are given for the five signal regions and their combination. The simplified model is for direct production of bottom squarks. A bottom squark decays into t chi1^(+-) and chi1^(+-) --> W^(+-) chi1^0 (see Fig. 7b in the paper). This particular model assumes that m(chi1^0)=2(chi1^0).

The signal event yields are calculated for all mSUGRA/CMSSM models with bRPV (see Fig. 8b in the paper). For each model, the values are given for the five signal regions and their combination. The model assumes tan(beta)=30, A0=2m0, mu>0, and bRPV.

The signal event yields are calculated for all simplified models. For each model, the values are given for the five signal regions and their combination. The simplified model is for direct production of squarks. Squarks decay into q q l l (l l) chi1^0 chi1^0 + neutrinos (see Fig. 6e in the paper).

The signal event yields are calculated for all simplified models. For each model, the values are given for the five signal regions and their combination. The simplified model is for direct pair-production of gluinos that decay via a two-step process into q q q q W Z W Z chi1^0 chi1^0 (see Fig. 6b in the paper).

The signal event yields are calculated for all simplified models. For each model, the values are given for the five signal regions and their combination. The simplified model is for direct pair-production of gluinos. A gluino decays into t stop. Consequently, a top squark squark decays into b chi1^(+-) and chi1^(+-) --> W^(+-) chi1^0 (see Fig. 5b in the paper). This particular model assumes that m(stop) < m(gluino), m(chi1^0)=6 GeV, and m(chi1^(+-))=118 GeV.

Experimental uncertainties on the signal event yields are calculated for all simplified models. For each model, the values are given for the five signal regions and their combination. The simplified model is for direct production of squarks that decay into two steps into q q W Z W Z chi1^0 chi1^0 (see Fig. 6c in the paper).

Experimental uncertainties on the signal event yields are calculated for all simplified extra dimension model (see Fig. 8d in the paper). For each model, the values are given for the five signal regions and their combination.

Experimental uncertainties on the signal event yields are calculated for all simplified models. For each model, the values are given for the five signal regions and their combination. The simplified model is for direct production of gluinos that decay via sleptons into q q q q l l (l l) chi1^0 chi1^0 + neutrinos (see Fig. 6d in the paper).

Experimental uncertainties on the signal event yields are calculated for all simplified models. For each model, the values are given for the five signal regions and their combination. The simplified model is for direct production of gluinos that decay into q q q q W W chi1^0 chi1^0 (see Fig. 6a in the paper).

Experimental uncertainties on the signal event yields are calculated for all simplified models. For each model, the values are given for the five signal regions and their combination. The simplified model is for direct production of gluinos that decay into t tbar t tbar chi1^0 chi1^0 (see Fig. 5a in the paper). This particular model assumes that top quark is much heavier than gluino.

Experimental uncertainties on the signal event yields are calculated for all mSUGRA models (see Fig. 8a in the paper). For each model, the values are given for the five signal regions and their combination. The model assumes tan(beta)=30, A0=2m0, and mu>0.

Experimental uncertainties on the signal event yields are calculated for all simplified models. For each model, the values are given for the five signal regions and their combination. The simplified model is for direct production of gluinos. A gluino decays into t c chi1^0 (see Fig. 5c in the paper). This particular model assumes that m(chi1^0) = m(stop) - 20 GeV.

Experimental uncertainties on the signal event yields are calculated for all GMSB models (see Fig. 8c in the paper). For each model, the values are given for the five signal regions and their combination. The model assumes mmess=250 TeV, m5=3, mu>0, and Cgrav=1.

Experimental uncertainties on the signal event yields are calculated for all simplified models. For each model, the values are given for the five signal regions and their combination. The simplified model is for direct production of bottom squarks. A bottom squark decays into t chi1^(+-) and chi1^(+-) --> W^(+-) chi1^0 (see Fig. 7a in the paper). This particular model assumes that m(chi1^0)=60 GeV.

Experimental uncertainties on the signal event yields are calculated for all simplified models. For each model, the values are given for the five signal regions and their combination. The simplified model is for direct production of gluinos and top squarks. Top squarks undergo R-parity violating decays into bs and gluinos decay into t stop (see Fig. 5d in the paper).

Experimental uncertainties on the signal event yields are calculated for all simplified models. For each model, the values are given for the five signal regions and their combination. The simplified model is for direct production of bottom squarks. A bottom squark decays into t chi1^(+-) and chi1^(+-) --> W^(+-) chi1^0 (see Fig. 7b in the paper). This particular model assumes that m(chi1^0)=2(chi1^0).

Experimental uncertainties on the signal event yields are calculated for all mSUGRA/CMSSM models with bRPV (see Fig. 8b in the paper). For each model, the values are given for the five signal regions and their combination. The model assumes tan(beta)=30, A0=2m0, mu>0, and bRPV.

Experimental uncertainties on the signal event yields are calculated for all simplified models. For each model, the values are given for the five signal regions and their combination. The simplified model is for direct production of squarks. Squarks decay into q q l l (l l) chi1^0 chi1^0 + neutrinos (see Fig. 6e in the paper).

Experimental uncertainties on the signal event yields are calculated for all simplified models. For each model, the values are given for the five signal regions and their combination. The simplified model is for direct pair-production of gluinos that decay via a two-step process into q q q q W Z W Z chi1^0 chi1^0 (see Fig. 6b in the paper).

Experimental uncertainties on the signal event yields are calculated for all simplified models. For each model, the values are given for the five signal regions and their combination. The simplified model is for direct pair-production of gluinos. A gluino decays into t stop. Consequently, a top squark squark decays into b chi1^(+-) and chi1^(+-) --> W^(+-) chi1^0 (see Fig. 5b in the paper). This particular model assumes that m(stop) < m(gluino), m(chi1^0)=6 GeV, and m(chi1^(+-))=118 GeV.

Statistical uncertainties on the signal event yields are calculated for all simplified models. For each model, the values are given for the five signal regions and their combination. The simplified model is for direct production of squarks that decay into two steps into q q W Z W Z chi1^0 chi1^0 (see Fig. 6c in the paper).

Statistical uncertainties on the signal event yields are calculated for all simplified extra dimension model (see Fig. 8d in the paper). For each model, the values are given for the five signal regions and their combination.

Statistical uncertainties on the signal event yields are calculated for all simplified models. For each model, the values are given for the five signal regions and their combination. The simplified model is for direct production of gluinos that decay via sleptons into q q q q l l (l l) chi1^0 chi1^0 + neutrinos (see Fig. 6d in the paper).

Statistical uncertainties on the signal event yields are calculated for all simplified models. For each model, the values are given for the five signal regions and their combination. The simplified model is for direct production of gluinos that decay into q q q q W W chi1^0 chi1^0 (see Fig. 6a in the paper).

Statistical uncertainties on the signal event yields are calculated for all simplified models. For each model, the values are given for the five signal regions and their combination. The simplified model is for direct production of gluinos that decay into t tbar t tbar chi1^0 chi1^0 (see Fig. 5a in the paper). This particular model assumes that top quark is much heavier than gluino.

Statistical uncertainties on the signal event yields are calculated for all mSUGRA models (see Fig. 8a in the paper). For each model, the values are given for the five signal regions and their combination. The model assumes tan(beta)=30, A0=2m0, and mu>0.

Statistical uncertainties on the signal event yields are calculated for all simplified models. For each model, the values are given for the five signal regions and their combination. The simplified model is for direct production of gluinos. A gluino decays into t c chi1^0 (see Fig. 5c in the paper). This particular model assumes that m(chi1^0) = m(stop) - 20 GeV.

Statistical uncertainties on the signal event yields are calculated for all GMSB models (see Fig. 8c in the paper). For each model, the values are given for the five signal regions and their combination. The model assumes mmess=250 TeV, m5=3, mu>0, and Cgrav=1.

Statistical uncertainties on the signal event yields are calculated for all simplified models. For each model, the values are given for the five signal regions and their combination. The simplified model is for direct production of bottom squarks. A bottom squark decays into t chi1^(+-) and chi1^(+-) --> W^(+-) chi1^0 (see Fig. 7a in the paper). This particular model assumes that m(chi1^0)=60 GeV.

Statistical uncertainties on the signal event yields are calculated for all simplified models. For each model, the values are given for the five signal regions and their combination. The simplified model is for direct production of gluinos and top squarks. Top squarks undergo R-parity violating decays into bs and gluinos decay into t stop (see Fig. 5d in the paper).

Statistical uncertainties on the signal event yields are calculated for all simplified models. For each model, the values are given for the five signal regions and their combination. The simplified model is for direct production of bottom squarks. A bottom squark decays into t chi1^(+-) and chi1^(+-) --> W^(+-) chi1^0 (see Fig. 7b in the paper). This particular model assumes that m(chi1^0)=2(chi1^0).

Statistical uncertainties on the signal event yields are calculated for all mSUGRA/CMSSM models with bRPV (see Fig. 8b in the paper). For each model, the values are given for the five signal regions and their combination. The model assumes tan(beta)=30, A0=2m0, mu>0, and bRPV.

Statistical uncertainties on the signal event yields are calculated for all simplified models. For each model, the values are given for the five signal regions and their combination. The simplified model is for direct production of squarks. Squarks decay into q q l l (l l) chi1^0 chi1^0 + neutrinos (see Fig. 6e in the paper).

Statistical uncertainties on the signal event yields are calculated for all simplified models. For each model, the values are given for the five signal regions and their combination. The simplified model is for direct pair-production of gluinos that decay via a two-step process into q q q q W Z W Z chi1^0 chi1^0 (see Fig. 6b in the paper).

Statistical uncertainties on the signal event yields are calculated for all simplified models. For each model, the values are given for the five signal regions and their combination. The simplified model is for direct pair-production of gluinos. A gluino decays into t stop. Consequently, a top squark squark decays into b chi1^(+-) and chi1^(+-) --> W ^(+-) chi1^0 (see Fig. 5b in the paper). This particular model assumes that m(stop) < m(gluino), m(chi1^0)=6 GeV, and m(chi1^(+-))=118 GeV.

The confidence levels are calculated for all simplified models. For each model, the observed and expected values are given. The simplified model is for direct production of gluinos that decay into t tbar t tbar chi1^0 chi1^0 (see Fig. 5a in the paper). This particular model assumes that top quark is much heavier than gluino.

The confidence levels are calculated for all simplified models. For each model, the observed and expected values are given. The simplified model is for direct production of squarks that decay into two steps into q q W Z W Z chi1^0 chi1^0 (see Fig. 6c in the paper).

The confidence levels are calculated for all simplified models. For each model, the values are given for the five signal regions and their combination. The simplified model is for direct pair-production of gluinos that decay via a two-step process into q q q q W Z W Z chi1^0 chi1^0 (see Fig. 6b in the paper).

The confidence levels are calculated for all simplified models. For each model, the expected and observed values are given. The simplified model is for direct production of gluinos that decay via sleptons into q q q q l l (l l) chi1^0 chi1^0 + neutrinos (see Fig. 6d in the paper).

The confidence levels are calculated for all simplified models. For each model, the expected and observed values are given. The simplified model is for direct pair-production of gluinos. A gluino decays into t stop. Consequently, a top squark squark decays into b chi1^(+-) and chi1^(+-) --> W^(+-) chi1^0 (see Fig. 5b in the paper). This particular model assumes that m(stop) < m(gluion), m(chi1^0)=6 GeV, and m(chi1^(+-))=118 GeV.

The confidence levels are calculated for all simplified models. For each model, the expected and observed values are given. The simplified model is for direct production of gluinos. A gluino decays into t c chi1^0 (see Fig. 5c in the paper). This particular model assumes that m(chi1^0) = m(stop) - 20 GeV.

The confidence levels are calculated for all simplified models. For each model, the expected and observed values are given. The simplified model is for direct production of bottom squarks. A bottom squark decays into t chi1^(+-) and chi1^(+-) --> W^(+-) chi1^0 (see Fig. 7b in the paper). This particular model assumes that m(chi1^0)=2(chi1^0).

The confidence levels are calculated for all simplified models. For each model, the expected and observed values are given. The simplified model is for direct production of bottom squarks. A bottom squark decays into t chi1^(+-) and chi1^(+-) --> W^(+-) chi1^0 (see Fig. 7a in the paper). This particular model assumes that m(chi1^0)=60 GeV.

The confidence levels are calculated for all simplified models. For each model, the expected and observed values are given. The simplified model is for direct production of squarks. Squarks decay into q q l l (l l) chi1^0 chi1^0 + neutrinos (see Fig. 6e in the paper).

The confidence levels are calculated for all GMSB models (see Fig. 8c in the paper). For each model, the expected and observed values are given. The model assumes mmess=250 TeV, m5=3, mu>0, and Cgrav=1.

The confidence levels are calculated for all simplified models. For each model, the expected and observed values are given. The simplified model is for direct production of gluinos and top squarks. Top squarks undergo R-parity violating decays into bs and gluinos decay into t stop (see Fig. 5d in the paper).

The confidence levels are calculated for all mSUGRA/CMSSM models with bRPV (see Fig. 8b in the paper). For each model, the expected and observed values are given. The model assumes tan(beta)=30, A0=2m0, mu>0, and bRPV.

The confidence levels are calculated for all simplified extra dimension model (see Fig. 8d in the paper). For each model, the expected and observed values are given.

The confidence levels are calculated for all simplified models. For each model, the expected and observed values are given. The simplified model is for direct production of gluinos that decay into q q q q W W chi1^0 chi1^0 (see Fig. 6a in the paper).

The confidence levels are calculated for all mSUGRA models (see Fig. 8a in the paper). For each model, the expected and observed values are given. The model assumes tan(beta)=30, A0=2m0, and mu>0.

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.

1$\tau$ MediumMass SR eff.

2$\tau$ Compressed SR eff.

2$\tau$ Compressed SR eff.

2$\tau$ HighMass SR eff.

2$\tau$ HighMass SR eff.

2$\tau$ multibin SR eff.

2$\tau$ multibin SR eff.

2$\tau$ GMSB SR eff.

2$\tau$ GMSB SR eff.

1$\tau$ Compressed SR eff.

1$\tau$ Compressed SR eff.

1$\tau$ MediumMass SR eff.

1$\tau$ MediumMass SR eff.

2$\tau$ Compressed SR eff.

2$\tau$ Compressed SR eff.

2$\tau$ HighMass SR eff.

2$\tau$ HighMass SR eff.

2$\tau$ multibin SR eff.

2$\tau$ multibin SR eff.

2$\tau$ GMSB SR eff.

2$\tau$ GMSB SR eff.

1$\tau$ Compressed SR acceptance.

1$\tau$ Compressed SR acceptance.

1$\tau$ MediumMass SR acceptance.

1$\tau$ MediumMass SR acceptance.

2$\tau$ Compressed SR acceptance.

2$\tau$ Compressed SR acceptance.

2$\tau$ HighMass SR acceptance.

2$\tau$ HighMass SR acceptance.

2$\tau$ multibin SR acceptance.

2$\tau$ multibin SR acceptance.

2$\tau$ GMSB SR acceptance.

2$\tau$ GMSB SR acceptance.

1$\tau$ Compressed SR acceptance.

1$\tau$ Compressed SR acceptance.

1$\tau$ MediumMass SR acceptance.

1$\tau$ MediumMass SR acceptance.

2$\tau$ Compressed SR acceptance.

2$\tau$ Compressed SR acceptance.

2$\tau$ HighMass SR acceptance.

2$\tau$ HighMass SR acceptance.

2$\tau$ multibin SR acceptance.

2$\tau$ multibin SR acceptance.

2$\tau$ GMSB SR acceptance.

2$\tau$ GMSB SR acceptance.

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

Summary of the total uncertainty in the background prediction for each SR of the tt0L-low, tt0L-high, tt1L and tt2L analysis channels in the statistical combination. Their dominant contributions are indicated by individual lines. Individual uncertainties can be correlated, and do not necessarily add up in quadrature to the total background uncertainty.

Exclusion limits for colour-neutral scalar mediator dark matter models as a function of the mediator mass $m(\phi)$ for a DM mass $m_{\chi} = 1$ GeV. Associated production of DM with both single top quarks ($tW$ and $tj$ channels) and top quark pairs is considered. The limits are calculated at 95% CL and are expressed in terms of the ratio of the excluded cross section to the cross section for a coupling assumption of $g = g_q = g_{\chi} = 1$. The solid (dashed) lines show the observed (expected) exclusion limits for each individual channel and their statistical combination.

Exclusion limits for colour-neutral pseudoscalar mediator dark matter models as a function of the mediator mass $m(a)$ for a DM mass $m_{\chi} = 1$ GeV. Associated production of DM with both single top quarks ($tW$ and $tj$ channels) and top quark pairs is considered. The limits are calculated at 95% CL and are expressed in terms of the ratio of the excluded cross section to the cross section for a coupling assumption of $g = g_q = g_{\chi} = 1$. The solid (dashed) lines show the observed (expected) exclusion limits for each individual channel and their statistical combination.

$E_{\text{T}}^{\text{miss}}$ distribution in SR0X for the tt0L-low analysis. The contributions from all SM backgrounds are shown after the profile likelihood simultaneous fit to all tt0L-low CRs, with the hatched bands representing the total uncertainty. The category '$t\bar{t}$ (other)' represents $t\bar{t}$ events without extra jets or events with extra light-flavour jets. 'Other' includes contributions from $t\bar{t}W$, $tZ$ and $tWZ$ processes. The expected distributions for selected signal models are shown as dashed lines. The overflow events are included in the last bin. The bottom panels show the ratio of the observed data to the total SM background prediction, with the hatched area representing the total uncertainty in the background prediction and the red arrows marking data outside the vertical-axis range.

$E_{\text{T}}^{\text{miss}}$ distribution in SRWX for the tt0L-low analysis. The contributions from all SM backgrounds are shown after the profile likelihood simultaneous fit to all tt0L-low CRs, with the hatched bands representing the total uncertainty. The category '$t\bar{t}$ (other)' represents $t\bar{t}$ events without extra jets or events with extra light-flavour jets. 'Other' includes contributions from $t\bar{t}W$, $tZ$ and $tWZ$ processes. The expected distributions for selected signal models are shown as dashed lines. The overflow events are included in the last bin. The bottom panels show the ratio of the observed data to the total SM background prediction, with the hatched area representing the total uncertainty in the background prediction and the red arrows marking data outside the vertical-axis range.

$E_{\text{T}}^{\text{miss}}$ distribution in SRTX for the tt0L-low analysis. The contributions from all SM backgrounds are shown after the profile likelihood simultaneous fit to all tt0L-low CRs, with the hatched bands representing the total uncertainty. The category '$t\bar{t}$ (other)' represents $t\bar{t}$ events without extra jets or events with extra light-flavour jets. 'Other' includes contributions from $t\bar{t}W$, $tZ$ and $tWZ$ processes. The expected distributions for selected signal models are shown as dashed lines. The overflow events are included in the last bin. The bottom panels show the ratio of the observed data to the total SM background prediction, with the hatched area representing the total uncertainty in the background prediction and the red arrows marking data outside the vertical-axis range.

Exclusion limits for colour-neutral scalar mediator dark matter models as a function of the mediator mass $m(\phi)$ for a DM mass $m_{\chi} = 1$ GeV. Associated production of DM with both single top quarks ($tW$ and $tj$ channels) and top quark pairs is considered. The limits are calculated at 95% CL and are expressed in terms of the ratio of the excluded cross section to the nominal cross section for a coupling assumption of $g = g_q = g_{\chi} = 1$. The solid (dashed) lines show the observed (expected) exclusion limits for the tt0L-high and tt0L-low analyses and their statistical combination.

Exclusion limits for colour-neutral pseudoscalar mediator dark matter models as a function of the mediator mass $m(a)$ for a DM mass $m_{\chi} = 1$ GeV. Associated production of DM with both single top quarks ($tW$ and $tj$ channels) and top quark pairs is considered. The limits are calculated at 95% CL and are expressed in terms of the ratio of the excluded cross section to the nominal cross section for a coupling assumption of $g = g_q = g_{\chi} = 1$. The solid (dashed) lines show the observed (expected) exclusion limits for the tt0L-high and tt0L-low analyses and their statistical combination.

Exclusion limits for colour-neutral scalar mediator dark matter models as a function of the mediator mass $m(\phi)$ for a DM mass $m_{\chi} = 1$ GeV. Only associated production of DM with top quark pairs is considered for this interpretation. The limits are calculated at 95% CL and are expressed in terms of the ratio of the excluded cross section to the cross section for a coupling assumption of $g = g_q = g_{\chi} = 1$. The solid (dashed) lines show the observed (expected) exclusion limits for each individual channel and their statistical combination.

Exclusion limits for colour-neutral pseudoscalar mediator dark matter models as a function of the mediator mass $m(a)$ for a DM mass $m_{\chi} = 1$ GeV. Only associated production of DM with top quark pairs is considered for this interpretation. The limits are calculated at 95% CL and are expressed in terms of the ratio of the excluded cross section to the cross section for a coupling assumption of $g = g_q = g_{\chi} = 1$. The solid (dashed) lines show the observed (expected) exclusion limits for each individual channel and their statistical combination.

Exclusion limits for colour-neutral scalar mediator dark matter models as a function of the mediator mass $m(\phi)$ for a DM mass $m_{\chi} = 1$ GeV. Only associated production of DM with top quark pairs is considered for this interpretation. The limits are calculated at 95% CL and are expressed in terms of the ratio of the excluded cross section to the nominal cross section for a coupling assumption of $g = g_q = g_{\chi} = 1$. The solid (dashed) lines show the observed (expected) exclusion limits for the tt0L-high and tt0L-low analyses and their statistical combination.

Exclusion limits for colour-neutral pseudoscalar mediator dark matter models as a function of the mediator mass $m(a)$ for a DM mass $m_{\chi} = 1$ GeV. Only associated production of DM with top quark pairs is considered for this interpretation. The limits are calculated at 95% CL and are expressed in terms of the ratio of the excluded cross section to the nominal cross section for a coupling assumption of $g = g_q = g_{\chi} = 1$. The solid (dashed) lines show the observed (expected) exclusion limits for the tt0L-high and tt0L-low analyses and their statistical combination.

Representative fit distribution in the different flavour leptons signal region for the tt2L analysis: each bin of such distribution, starting from the red arrow, corresponds to a single SR included in the fit. 'FNP' includes the contribution from fake/non-prompt lepton background arising from jets (mainly $\pi/K$, heavy-flavour hadron decays and photon conversion) misidentified as leptons, estimated in a purely data-driven way. 'Other' includes contributions from $t\bar{t}W$, $tZ$ and $tWZ$ processes. The total uncertainty in the SM expectation is represented with hatched bands and the expected distributions for selected signal models are shown as dashed lines.

Signal acceptance in SR0X, SRWX and SRTX for simplified DM+$t\bar{t}$ model, defined as the number of accepted events at generator level in signal Monte Carlo simulation divided by the total number of events in the sample.

Signal acceptance in SR0X, SRWX and SRTX for simplified DM+$tW$ model, defined as the number of accepted events at generator level in signal Monte Carlo simulation divided by the total number of events in the sample.

Signal acceptance in SR0X, SRWX and SRTX for simplified DM+$tj$ model, defined as the number of accepted events at generator level in signal Monte Carlo simulation divided by the total number of events in the sample.

Signal efficiency in SR0X, SRWX and SRTX for simplified DM+$t\bar{t}$ model, defined as the number of selected reconstructed events divided by the acceptance.

Signal efficiency in SR0X, SRWX and SRTX for simplified DM+$tW$ model, defined as the number of selected reconstructed events divided by the acceptance.

Signal efficiency in SR0X, SRWX and SRTX for simplified DM+$tj$ model, defined as the number of selected reconstructed events divided by the acceptance.

Cutflow for the reference point DM+$t\bar{t}$ $m(\phi, \chi) = (10, 1)$ GeV in signal region SR0X. The column labelled 'weighted' shows the event yield including all correction factors applied to simulation, and is normalised to 139 fb$^{-1}$. A notable exception concerns the 'weighted' numbers in the first and the second row, labelled 'Total' and 'Filtered', which correspond to $\mathcal{L}\cdot\sigma$ and $\mathcal{L}\cdot\sigma\cdot\epsilon$ expected, respectively. The 'Skim' selection requires the $p_{\text{T}}$ of the leading four jets to be above (80, 60, 40, 40) GeV, the missing transverse momentum $E_{\text{T}}^{\text{miss}} > 140$ GeV, the missing momentum significance $\mathcal{S} > 8$, $\Delta\phi_{\min}(\vec{p}_{\text{T,1-4}},\vec{p}_{\text{T}}^{\text{miss}}) > 0.4$ and a lepton veto. The 'Orthogonalisation' selection is defined in the main body. In total 2045000 raw MC events were generated prior to the specified cuts, with the column 'Unweighted yield' collecting the numbers after each cut.

Cutflow for the reference point DM+$t\bar{t}$ $m(\phi, \chi) = (10, 1)$ GeV in signal region SRWX. The column labelled 'weighted' shows the event yield including all correction factors applied to simulation, and is normalised to 139 fb$^{-1}$. A notable exception concerns the 'weighted' numbers in the first and the second row, labelled 'Total' and 'Filtered', which correspond to $\mathcal{L}\cdot\sigma$ and $\mathcal{L}\cdot\sigma\cdot\epsilon$ expected, respectively. The 'Skim' selection requires the $p_{\text{T}}$ of the leading four jets to be above (80, 60, 40, 40) GeV, the missing transverse momentum $E_{\text{T}}^{\text{miss}} > 140$ GeV, the missing momentum significance $\mathcal{S} > 8$, $\Delta\phi_{\min}(\vec{p}_{\text{T,1-4}},\vec{p}_{\text{T}}^{\text{miss}}) > 0.4$ and a lepton veto. The 'Orthogonalisation' selection is defined in the main body. In total 2045000 raw MC events were generated prior to the specified cuts, with the column 'Unweighted yield' collecting the numbers after each cut.

Cutflow for the reference point DM+$t\bar{t}$ $m(\phi, \chi) = (10, 1)$ GeV in signal region SRTX. The column labelled 'weighted' shows the event yield including all correction factors applied to simulation, and is normalised to 139 fb$^{-1}$. A notable exception concerns the 'weighted' numbers in the first and the second row, labelled 'Total' and 'Filtered', which correspond to $\mathcal{L}\cdot\sigma$ and $\mathcal{L}\cdot\sigma\cdot\epsilon$ expected, respectively. The 'Skim' selection requires the $p_{\text{T}}$ of the leading four jets to be above (80, 60, 40, 40) GeV, the missing transverse momentum $E_{\text{T}}^{\text{miss}} > 140$ GeV, the missing momentum significance $\mathcal{S} > 8$, $\Delta\phi_{\min}(\vec{p}_{\text{T,1-4}},\vec{p}_{\text{T}}^{\text{miss}}) > 0.4$ and a lepton veto. The 'Orthogonalisation' selection is defined in the main body. In total 2045000 raw MC events were generated prior to the specified cuts, with the column 'Unweighted yield' collecting the numbers after each cut.

Cutflow for the reference point DM+$t\bar{t}$ $m(a, \chi) = (10, 1)$ GeV in signal region SR0X. The column labelled 'weighted' shows the event yield including all correction factors applied to simulation, and is normalised to 139 fb$^{-1}$. A notable exception concerns the 'weighted' numbers in the first and the second row, labelled 'Total' and 'Filtered', which correspond to $\mathcal{L}\cdot\sigma$ and $\mathcal{L}\cdot\sigma\cdot\epsilon$ expected, respectively. The 'Skim' selection requires the $p_{\text{T}}$ of the leading four jets to be above (80, 60, 40, 40) GeV, the missing transverse momentum $E_{\text{T}}^{\text{miss}} > 140$ GeV, the missing momentum significance $\mathcal{S} > 8$, $\Delta\phi_{\min}(\vec{p}_{\text{T,1-4}},\vec{p}_{\text{T}}^{\text{miss}}) > 0.4$ and a lepton veto. The 'Orthogonalisation' selection is defined in the main body. In total 400000 raw MC events were generated prior to the specified cuts, with the column 'Unweighted yield' collecting the numbers after each cut.

Cutflow for the reference point DM+$t\bar{t}$ $m(a, \chi) = (10, 1)$ GeV in signal region SRWX. The column labelled 'weighted' shows the event yield including all correction factors applied to simulation, and is normalised to 139 fb$^{-1}$. A notable exception concerns the 'weighted' numbers in the first and the second row, labelled 'Total' and 'Filtered', which correspond to $\mathcal{L}\cdot\sigma$ and $\mathcal{L}\cdot\sigma\cdot\epsilon$ expected, respectively. The 'Skim' selection requires the $p_{\text{T}}$ of the leading four jets to be above (80, 60, 40, 40) GeV, the missing transverse momentum $E_{\text{T}}^{\text{miss}} > 140$ GeV, the missing momentum significance $\mathcal{S} > 8$, $\Delta\phi_{\min}(\vec{p}_{\text{T,1-4}},\vec{p}_{\text{T}}^{\text{miss}}) > 0.4$ and a lepton veto. The 'Orthogonalisation' selection is defined in the main body. In total 400000 raw MC events were generated prior to the specified cuts, with the column 'Unweighted yield' collecting the numbers after each cut.

Cutflow for the reference point DM+$t\bar{t}$ $m(a, \chi) = (10, 1)$ GeV in signal region SRTX. The column labelled 'weighted' shows the event yield including all correction factors applied to simulation, and is normalised to 139 fb$^{-1}$. A notable exception concerns the 'weighted' numbers in the first and the second row, labelled 'Total' and 'Filtered', which correspond to $\mathcal{L}\cdot\sigma$ and $\mathcal{L}\cdot\sigma\cdot\epsilon$ expected, respectively. The 'Skim' selection requires the $p_{\text{T}}$ of the leading four jets to be above (80, 60, 40, 40) GeV, the missing transverse momentum $E_{\text{T}}^{\text{miss}} > 140$ GeV, the missing momentum significance $\mathcal{S} > 8$, $\Delta\phi_{\min}(\vec{p}_{\text{T,1-4}},\vec{p}_{\text{T}}^{\text{miss}}) > 0.4$ and a lepton veto. The 'Orthogonalisation' selection is defined in the main body. In total 400000 raw MC events were generated prior to the specified cuts, with the column 'Unweighted yield' collecting the numbers after each cut.

Cutflow for the reference point DM+$tW$ $m(\phi, \chi) = (10, 1)$ GeV in signal region SR0X. The column labelled 'weighted' shows the event yield including all correction factors applied to simulation, and is normalised to 139 fb$^{-1}$. A notable exception concerns the 'weighted' numbers in the first and the second row, labelled 'Total' and 'Filtered', which correspond to $\mathcal{L}\cdot\sigma$ and $\mathcal{L}\cdot\sigma\cdot\epsilon$ expected, respectively. The 'Skim' selection requires the $p_{\text{T}}$ of the leading four jets to be above (80, 60, 40, 40) GeV, the missing transverse momentum $E_{\text{T}}^{\text{miss}} > 140$ GeV, the missing momentum significance $\mathcal{S} > 8$, $\Delta\phi_{\min}(\vec{p}_{\text{T,1-4}},\vec{p}_{\text{T}}^{\text{miss}}) > 0.4$ and a lepton veto. The 'Orthogonalisation' selection is defined in the main body. In total 120000 raw MC events were generated prior to the specified cuts, with the column 'Unweighted yield' collecting the numbers after each cut.

Cutflow for the reference point DM+$tW$ $m(\phi, \chi) = (10, 1)$ GeV in signal region SRWX. The column labelled 'weighted' shows the event yield including all correction factors applied to simulation, and is normalised to 139 fb$^{-1}$. A notable exception concerns the 'weighted' numbers in the first and the second row, labelled 'Total' and 'Filtered', which correspond to $\mathcal{L}\cdot\sigma$ and $\mathcal{L}\cdot\sigma\cdot\epsilon$ expected, respectively. The 'Skim' selection requires the $p_{\text{T}}$ of the leading four jets to be above (80, 60, 40, 40) GeV, the missing transverse momentum $E_{\text{T}}^{\text{miss}} > 140$ GeV, the missing momentum significance $\mathcal{S} > 8$, $\Delta\phi_{\min}(\vec{p}_{\text{T,1-4}},\vec{p}_{\text{T}}^{\text{miss}}) > 0.4$ and a lepton veto. The 'Orthogonalisation' selection is defined in the main body. In total 120000 raw MC events were generated prior to the specified cuts, with the column 'Unweighted yield' collecting the numbers after each cut.

Cutflow for the reference point DM+$tW$ $m(\phi, \chi) = (10, 1)$ GeV in signal region SRTX. The column labelled 'weighted' shows the event yield including all correction factors applied to simulation, and is normalised to 139 fb$^{-1}$. A notable exception concerns the 'weighted' numbers in the first and the second row, labelled 'Total' and 'Filtered', which correspond to $\mathcal{L}\cdot\sigma$ and $\mathcal{L}\cdot\sigma\cdot\epsilon$ expected, respectively. The 'Skim' selection requires the $p_{\text{T}}$ of the leading four jets to be above (80, 60, 40, 40) GeV, the missing transverse momentum $E_{\text{T}}^{\text{miss}} > 140$ GeV, the missing momentum significance $\mathcal{S} > 8$, $\Delta\phi_{\min}(\vec{p}_{\text{T,1-4}},\vec{p}_{\text{T}}^{\text{miss}}) > 0.4$ and a lepton veto. The 'Orthogonalisation' selection is defined in the main body. In total 120000 raw MC events were generated prior to the specified cuts, with the column 'Unweighted yield' collecting the numbers after each cut.

Cutflow for the reference point DM+$tW$ $m(a, \chi) = (10, 1)$ GeV in signal region SR0X. The column labelled 'weighted' shows the event yield including all correction factors applied to simulation, and is normalised to 139 fb$^{-1}$. A notable exception concerns the 'weighted' numbers in the first and the second row, labelled 'Total' and 'Filtered', which correspond to $\mathcal{L}\cdot\sigma$ and $\mathcal{L}\cdot\sigma\cdot\epsilon$ expected, respectively. The 'Skim' selection requires the $p_{\text{T}}$ of the leading four jets to be above (80, 60, 40, 40) GeV, the missing transverse momentum $E_{\text{T}}^{\text{miss}} > 140$ GeV, the missing momentum significance $\mathcal{S} > 8$, $\Delta\phi_{\min}(\vec{p}_{\text{T,1-4}},\vec{p}_{\text{T}}^{\text{miss}}) > 0.4$ and a lepton veto. The 'Orthogonalisation' selection is defined in the main body. In total 100000 raw MC events were generated prior to the specified cuts, with the column 'Unweighted yield' collecting the numbers after each cut.

Cutflow for the reference point DM+$tW$ $m(a, \chi) = (10, 1)$ GeV in signal region SRWX. The column labelled 'weighted' shows the event yield including all correction factors applied to simulation, and is normalised to 139 fb$^{-1}$. A notable exception concerns the 'weighted' numbers in the first and the second row, labelled 'Total' and 'Filtered', which correspond to $\mathcal{L}\cdot\sigma$ and $\mathcal{L}\cdot\sigma\cdot\epsilon$ expected, respectively. The 'Skim' selection requires the $p_{\text{T}}$ of the leading four jets to be above (80, 60, 40, 40) GeV, the missing transverse momentum $E_{\text{T}}^{\text{miss}} > 140$ GeV, the missing momentum significance $\mathcal{S} > 8$, $\Delta\phi_{\min}(\vec{p}_{\text{T,1-4}},\vec{p}_{\text{T}}^{\text{miss}}) > 0.4$ and a lepton veto. The 'Orthogonalisation' selection is defined in the main body. In total 100000 raw MC events were generated prior to the specified cuts, with the column 'Unweighted yield' collecting the numbers after each cut.

Cutflow for the reference point DM+$tW$ $m(a, \chi) = (10, 1)$ GeV in signal region SRTX. The column labelled 'weighted' shows the event yield including all correction factors applied to simulation, and is normalised to 139 fb$^{-1}$. A notable exception concerns the 'weighted' numbers in the first and the second row, labelled 'Total' and 'Filtered', which correspond to $\mathcal{L}\cdot\sigma$ and $\mathcal{L}\cdot\sigma\cdot\epsilon$ expected, respectively. The 'Skim' selection requires the $p_{\text{T}}$ of the leading four jets to be above (80, 60, 40, 40) GeV, the missing transverse momentum $E_{\text{T}}^{\text{miss}} > 140$ GeV, the missing momentum significance $\mathcal{S} > 8$, $\Delta\phi_{\min}(\vec{p}_{\text{T,1-4}},\vec{p}_{\text{T}}^{\text{miss}}) > 0.4$ and a lepton veto. The 'Orthogonalisation' selection is defined in the main body. In total 100000 raw MC events were generated prior to the specified cuts, with the column 'Unweighted yield' collecting the numbers after each cut.

Cutflow for the reference point DM+$tj$ $m(\phi, \chi) = (10, 1)$ GeV in signal region SR0X. The column labelled 'weighted' shows the event yield including all correction factors applied to simulation, and is normalised to 139 fb$^{-1}$. A notable exception concerns the 'weighted' numbers in the first and the second row, labelled 'Total' and 'Filtered', which correspond to $\mathcal{L}\cdot\sigma$ and $\mathcal{L}\cdot\sigma\cdot\epsilon$ expected, respectively. The 'Skim' selection requires the $p_{\text{T}}$ of the leading four jets to be above (80, 60, 40, 40) GeV, the missing transverse momentum $E_{\text{T}}^{\text{miss}} > 140$ GeV, the missing momentum significance $\mathcal{S} > 8$, $\Delta\phi_{\min}(\vec{p}_{\text{T,1-4}},\vec{p}_{\text{T}}^{\text{miss}}) > 0.4$ and a lepton veto. The 'Orthogonalisation' selection is defined in the main body. In total 169000 raw MC events were generated prior to the specified cuts, with the column 'Unweighted yield' collecting the numbers after each cut.

Cutflow for the reference point DM+$tj$ $m(\phi, \chi) = (10, 1)$ GeV in signal region SRWX. The column labelled 'weighted' shows the event yield including all correction factors applied to simulation, and is normalised to 139 fb$^{-1}$. A notable exception concerns the 'weighted' numbers in the first and the second row, labelled 'Total' and 'Filtered', which correspond to $\mathcal{L}\cdot\sigma$ and $\mathcal{L}\cdot\sigma\cdot\epsilon$ expected, respectively. The 'Skim' selection requires the $p_{\text{T}}$ of the leading four jets to be above (80, 60, 40, 40) GeV, the missing transverse momentum $E_{\text{T}}^{\text{miss}} > 140$ GeV, the missing momentum significance $\mathcal{S} > 8$, $\Delta\phi_{\min}(\vec{p}_{\text{T,1-4}},\vec{p}_{\text{T}}^{\text{miss}}) > 0.4$ and a lepton veto. The 'Orthogonalisation' selection is defined in the main body. In total 169000 raw MC events were generated prior to the specified cuts, with the column 'Unweighted yield' collecting the numbers after each cut.

Cutflow for the reference point DM+$tj$ $m(\phi, \chi) = (10, 1)$ GeV in signal region SRTX. The column labelled 'weighted' shows the event yield including all correction factors applied to simulation, and is normalised to 139 fb$^{-1}$. A notable exception concerns the 'weighted' numbers in the first and the second row, labelled 'Total' and 'Filtered', which correspond to $\mathcal{L}\cdot\sigma$ and $\mathcal{L}\cdot\sigma\cdot\epsilon$ expected, respectively. The 'Skim' selection requires the $p_{\text{T}}$ of the leading four jets to be above (80, 60, 40, 40) GeV, the missing transverse momentum $E_{\text{T}}^{\text{miss}} > 140$ GeV, the missing momentum significance $\mathcal{S} > 8$, $\Delta\phi_{\min}(\vec{p}_{\text{T,1-4}},\vec{p}_{\text{T}}^{\text{miss}}) > 0.4$ and a lepton veto. The 'Orthogonalisation' selection is defined in the main body. In total 169000 raw MC events were generated prior to the specified cuts, with the column 'Unweighted yield' collecting the numbers after each cut.

Cutflow for the reference point DM+$tj$ $m(a, \chi) = (10, 1)$ GeV in signal region SR0X. The column labelled 'weighted' shows the event yield including all correction factors applied to simulation, and is normalised to 139 fb$^{-1}$. A notable exception concerns the 'weighted' numbers in the first and the second row, labelled 'Total' and 'Filtered', which correspond to $\mathcal{L}\cdot\sigma$ and $\mathcal{L}\cdot\sigma\cdot\epsilon$ expected, respectively. The 'Skim' selection requires the $p_{\text{T}}$ of the leading four jets to be above (80, 60, 40, 40) GeV, the missing transverse momentum $E_{\text{T}}^{\text{miss}} > 140$ GeV, the missing momentum significance $\mathcal{S} > 8$, $\Delta\phi_{\min}(\vec{p}_{\text{T,1-4}},\vec{p}_{\text{T}}^{\text{miss}}) > 0.4$ and a lepton veto. The 'Orthogonalisation' selection is defined in the main body. In total 140000 raw MC events were generated prior to the specified cuts, with the column 'Unweighted yield' collecting the numbers after each cut.

Cutflow for the reference point DM+$tj$ $m(a, \chi) = (10, 1)$ GeV in signal region SRWX. The column labelled 'weighted' shows the event yield including all correction factors applied to simulation, and is normalised to 139 fb$^{-1}$. A notable exception concerns the 'weighted' numbers in the first and the second row, labelled 'Total' and 'Filtered', which correspond to $\mathcal{L}\cdot\sigma$ and $\mathcal{L}\cdot\sigma\cdot\epsilon$ expected, respectively. The 'Skim' selection requires the $p_{\text{T}}$ of the leading four jets to be above (80, 60, 40, 40) GeV, the missing transverse momentum $E_{\text{T}}^{\text{miss}} > 140$ GeV, the missing momentum significance $\mathcal{S} > 8$, $\Delta\phi_{\min}(\vec{p}_{\text{T,1-4}},\vec{p}_{\text{T}}^{\text{miss}}) > 0.4$ and a lepton veto. The 'Orthogonalisation' selection is defined in the main body. In total 140000 raw MC events were generated prior to the specified cuts, with the column 'Unweighted yield' collecting the numbers after each cut.

Cutflow for the reference point DM+$tj$ $m(a, \chi) = (10, 1)$ GeV in signal region SRTX. The column labelled 'weighted' shows the event yield including all correction factors applied to simulation, and is normalised to 139 fb$^{-1}$. A notable exception concerns the 'weighted' numbers in the first and the second row, labelled 'Total' and 'Filtered', which correspond to $\mathcal{L}\cdot\sigma$ and $\mathcal{L}\cdot\sigma\cdot\epsilon$ expected, respectively. The 'Skim' selection requires the $p_{\text{T}}$ of the leading four jets to be above (80, 60, 40, 40) GeV, the missing transverse momentum $E_{\text{T}}^{\text{miss}} > 140$ GeV, the missing momentum significance $\mathcal{S} > 8$, $\Delta\phi_{\min}(\vec{p}_{\text{T,1-4}},\vec{p}_{\text{T}}^{\text{miss}}) > 0.4$ and a lepton veto. The 'Orthogonalisation' selection is defined in the main body. In total 140000 raw MC events were generated prior to the specified cuts, with the column 'Unweighted yield' collecting the numbers after each cut.