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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

ggH default MC + XH predictions

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

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

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

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

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

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

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

Fractions of selected events for several signal MC samples with a mass difference $\Delta m = 100$ GeV, illustrating how $\mathcal{A}\times\varepsilon$ varies with the model parameters.

Two-dimensional distribution of $m_{\mathrm{DV}}$ and track multiplicity for DVs in data events and events of a $R$-hadron signal sample with $m_{\tilde{g}} = 1.4$ TeV, $m_{\tilde{\chi}_{1}^{0}} = 100$ GeV and $\tau_{\tilde{g}} = 1$ ns that satisfy all signal region event selection criteria.

Two-dimensional distribution of $m_{\mathrm{DV}}$ and track multiplicity for DVs in data events and events of a $R$-hadron signal sample with $m_{\tilde{g}} = 1.4$ TeV, $m_{\tilde{\chi}_{1}^{0}} = 1.32$ TeV and $\tau_{\tilde{g}} = 1$ ns that satisfy all signal region event selection criteria.

Observed cross section upper 95% CL limits as a function of $m_{\tilde{g}}$ and $\tau$ for $m_{\tilde{\chi}_{1}^{0}}=100$ GeV. For the mass limits see the entry of Figure 8b.

Upper 95% CL limits on the signal cross section for $m_{\tilde{g}}=1.4$ TeV and fixed $m_{\tilde{\chi}_{1}^{0}}=100$ GeV as a function of lifetime $\tau$.

Vertex reconstruction efficiency as a function of radial position $R$ with and without the special LRT processing for one $R$-hadron signal sample with $m_{\tilde{g}} = 1.2$ TeV, $m_{\tilde{\chi}_{1}^{0}} = 100$ GeV and $\tau_{\tilde{g}} = 1$ ns. The efficiency is defined as the probability for a true LLP decay to be matched with a reconstructed DV fulfilling the vertex preselection criteria in events with a reconstructed primary vertex.

Upper 95% CL limits on the signal cross section for $m_{\tilde{g}}=2.0$ TeV and fixed $m_{\tilde{\chi}_{1}^{0}}=100$ GeV as a function of lifetime $\tau$.

Vertex reconstruction efficiency as a function of radial position $R$ for two $R$-hadron signal samples with $m_{\tilde{g}} = 1.2$ TeV, $\tau_{\tilde{g}} = 1$ ns and different neutralino masses. The efficiency is defined as the probability for a true LLP decay to be matched with a reconstructed DV fulfilling the vertex preselection criteria in events with a reconstructed primary vertex.

Lower 95% CL limits on $m_{\tilde{g}}$ for fixed $m_{\tilde{\chi}_{1}^{0}}=100$ GeV as a function of lifetime $\tau$.

Fractions of selected events for several signal MC samples with a gluino lifetime $\tau = 1$ ns, illustrating how $\mathcal{A}\times\varepsilon$ varies with the model parameters.

Upper 95% CL limits on the signal cross section for $m_{\tilde{g}}=1.4$ TeV and fixed $\Delta m=100$ GeV as a function of lifetime $\tau$.

Fractions of selected events for several signal MC samples with a mass difference $\Delta m = 100$ GeV, illustrating how $\mathcal{A}\times\varepsilon$ varies with the model parameters.

Upper 95% CL limits on the signal cross section for $m_{\tilde{g}}=2.0$ TeV and fixed $\Delta m=100$ GeV as a function of lifetime $\tau$.

Two-dimensional distribution of $m_{\mathrm{DV}}$ and track multiplicity for DVs in data events and events of a $R$-hadron signal sample with $m_{\tilde{g}} = 1.4$ TeV, $m_{\tilde{\chi}_{1}^{0}} = 100$ GeV and $\tau_{\tilde{g}} = 1$ ns that satisfy all signal region event selection criteria.

Lower 95% CL limit on $m_{\tilde{g}}$ for fixed $\Delta m=100$ GeV as a function of lifetime $\tau$.

Two-dimensional distribution of $m_{\mathrm{DV}}$ and track multiplicity for DVs in data events and events of a $R$-hadron signal sample with $m_{\tilde{g}} = 1.4$ TeV, $m_{\tilde{\chi}_{1}^{0}} = 1.32$ TeV and $\tau_{\tilde{g}} = 1$ ns that satisfy all signal region event selection criteria.

Upper 95% CL limits on the signal cross section for $m_{\tilde{g}}=1.4$ TeV and fixed $\tau=1$ ns as a function of $m_{\tilde{\chi}_{1}^{0}}$.

Upper 95% CL limits on the signal cross section for $m_{\tilde{g}}=1.4$ TeV and fixed $\Delta m=100$ GeV as a function of lifetime $\tau$.

Upper 95% CL limits on the signal cross section for $m_{\tilde{g}}=2.0$ TeV and fixed $\tau=1$ ns as a function of $m_{\tilde{\chi}_{1}^{0}}$.

Upper 95% CL limits on the signal cross section for $m_{\tilde{g}}=2.0$ TeV and fixed $\Delta m=100$ GeV as a function of lifetime $\tau$.

Observed 95% CL limit as a function of $m_{\tilde{g}}$ and $m_{\tilde{\chi}_{1}^{0}}$ for fixed $\tau=1$ ns.

Lower 95% CL limit on $m_{\tilde{g}}$ for fixed $\Delta m=100$ GeV as a function of lifetime $\tau$.

Two-dimensional distributions of $x$-$y$ positions of vertices observed in the data passing the vertex pre-selection and satisfying all signal region event-level requirements.

Two-dimensional distributions of $x$-$y$ positions of vertices observed in the data passing the vertex pre-selection and satisfying all signal region event-level requirements.

Distribution of the mass $m_{\mathrm{DV}}$ for vertices in data events and in events of five $R$-hadron signal samples with $m_{\tilde{g}} = 1.2$ TeV, $m_{\tilde{\chi}_{1}^{0}} = 100$ GeV and different $\tau_{\tilde{g}}$ that satisfy the signal region event requirements. All DV selections are applied except for the $m_{\mathrm{DV}}$ and track multiplicity requirements.

Distribution of the mass $m_{\mathrm{DV}}$ for vertices in data events and in events of five $R$-hadron signal samples with $m_{\tilde{g}} = 1.2$ TeV, $m_{\tilde{\chi}_{1}^{0}} = 100$ GeV and different $\tau_{\tilde{g}}$ that satisfy the signal region event requirements. All DV selections are applied except for the $m_{\mathrm{DV}}$ and track multiplicity requirements.

Distribution of the track multiplicity $n_{\mathrm{Tracks}}$ for vertices in data events and events of five $R$-hadron signal samples with $m_{\tilde{g}} = 1.2$ TeV, $m_{\tilde{\chi}_{1}^{0}} = 100$ GeV and and different $\tau_{\tilde{g}}$ that satisfy the signal region event requirements. All DV selections are applied except for the $m_{\mathrm{DV}}$ and track multiplicity requirements. The track multiplicity distribution requires vertices to have $m_{\mathrm{DV}}>3$ GeV.

Distribution of the track multiplicity $n_{\mathrm{Tracks}}$ for vertices in data events and events of five $R$-hadron signal samples with $m_{\tilde{g}} = 1.2$ TeV, $m_{\tilde{\chi}_{1}^{0}} = 100$ GeV and and different $\tau_{\tilde{g}}$ that satisfy the signal region event requirements. All DV selections are applied except for the $m_{\mathrm{DV}}$ and track multiplicity requirements. The track multiplicity distribution requires vertices to have $m_{\mathrm{DV}}>3$ GeV.

Observed cross section upper 95% CL limits as a function of $m_{\tilde{g}}$ and $\tau$ for $m_{\tilde{\chi}_{1}^{0}}=100$ GeV. For the mass limits see the entry of Figure 8b.

Observed cross section upper 95% CL limits as a function of $m_{\tilde{g}}$ and $\tau$ for $\Delta m=100$ GeV. For the mass limits see the entry of Figure 9b.

Observed cross section upper 95% CL limits as a function of $m_{\tilde{g}}$ and $\tau$ for $\Delta m=100$ GeV. For the mass limits see the entry of Figure 9b.

Observed cross section upper 95% CL limits as a function of $m_{\tilde{\chi}_{1}^{0}}$ and $m_{\tilde{g}}$ for $\tau = 1$ ns. For the mass limits see the entry of Figure 10b.

Observed cross section upper 95% CL limits as a function of $m_{\tilde{\chi}_{1}^{0}}$ and $m_{\tilde{g}}$ for $\tau = 1$ ns. For the mass limits see the entry of Figure 10b.

Parameterized event selection efficiencies as a function of truth MET for events which have all truth decay vertices occurring before the start of the ATLAS calorimeter. Event-level efficiencies are evaluated for events that have truth MET $> 200$ GeV, pass the trackless jet requirement, and have at least one displaced truth decay within the fiducial volume. To satisfy the event-level efficiency, events must then pass the full event selection.

Parameterized event selection efficiencies as a function of truth MET for events which have all truth decay vertices occurring before the start of the ATLAS calorimeter. Event-level efficiencies are evaluated for events that have truth MET $> 200$ GeV, pass the trackless jet requirement, and have at least one displaced truth decay within the fiducial volume. To satisfy the event-level efficiency, events must then pass the full event selection.

Parameterized event selection efficiencies as a function of truth MET for events which have the furthest truth decay occurring inside the calorimeter. Event-level efficiencies are evaluated for events that have truth MET $> 200$ GeV, pass the trackless jet requirement, and have at least one displaced truth decay within the fiducial volume. To satisfy the event-level efficiency, events must then pass the full event selection.

Parameterized event selection efficiencies as a function of truth MET for events which have the furthest truth decay occurring inside the calorimeter. Event-level efficiencies are evaluated for events that have truth MET $> 200$ GeV, pass the trackless jet requirement, and have at least one displaced truth decay within the fiducial volume. To satisfy the event-level efficiency, events must then pass the full event selection.

Parameterized event selection efficiencies as a function of truth MET for events which have the furthest truth decay occurring after the end of the ATLAS calorimeter. Event-level efficiencies are evaluated for events that have truth MET $> 200$ GeV, pass the trackless jet requirement, and have at least one displaced truth decay within the fiducial volume. To satisfy the event-level efficiency, events must then pass the full event selection.

Parameterized event selection efficiencies as a function of truth MET for events which have the furthest truth decay occurring after the end of the ATLAS calorimeter. Event-level efficiencies are evaluated for events that have truth MET $> 200$ GeV, pass the trackless jet requirement, and have at least one displaced truth decay within the fiducial volume. To satisfy the event-level efficiency, events must then pass the full event selection.

Parameterized vertex level efficiencies as a function of number of particles associated to a truth decay vertex, and the vertex invariant mass for truth decays with $4$ mm $< R_{\mathrm{decay}} < 22$ mm. Selected particles are required to have nonzero electric charge, $p_{T}(|Q|=1) > 1$ GeV, and $d_0 > 2$ mm. The per-vertex efficiency is evaluated only for truth vertices that have at least 5 associated tracks, an invariant mass $> 10$ GeV, and are in the region $4$ mm $< R_{\mathrm{decay}} < 300$ mm, and $|Z_{\mathrm{decay}}| < 300$ mm. A truth vertex satisfies the vertex level efficiency if it can be matched to a reconstructed DV which passes the final vertex selection.

Parameterized vertex level efficiencies as a function of number of particles associated to a truth decay vertex, and the vertex invariant mass for truth decays with $4$ mm $< R_{\mathrm{decay}} < 22$ mm. Selected particles are required to have nonzero electric charge, $p_{T}(|Q|=1) > 1$ GeV, and $d_0 > 2$ mm. The per-vertex efficiency is evaluated only for truth vertices that have at least 5 associated tracks, an invariant mass $> 10$ GeV, and are in the region $4$ mm $< R_{\mathrm{decay}} < 300$ mm, and $|Z_{\mathrm{decay}}| < 300$ mm. A truth vertex satisfies the vertex level efficiency if it can be matched to a reconstructed DV which passes the final vertex selection.

Parameterized vertex level efficiencies as a function of number of particles associated to a truth decay vertex, and the vertex invariant mass for truth decays with $22$ mm $< R_{\mathrm{decay}} < 25$ mm. Selected particles are required to have nonzero electric charge, $p_{T}(|Q|=1) > 1$ GeV, and $d_0 > 2$ mm. The per-vertex efficiency is evaluated only for truth vertices that have at least 5 associated tracks, an invariant mass $> 10$ GeV, and are in the region $4$ mm $< R_{\mathrm{decay}} < 300$ mm, and $|Z_{\mathrm{decay}}| < 300$ mm. A truth vertex satisfies the vertex level efficiency if it can be matched to a reconstructed DV which passes the final vertex selection.

Parameterized vertex level efficiencies as a function of number of particles associated to a truth decay vertex, and the vertex invariant mass for truth decays with $22$ mm $< R_{\mathrm{decay}} < 25$ mm. Selected particles are required to have nonzero electric charge, $p_{T}(|Q|=1) > 1$ GeV, and $d_0 > 2$ mm. The per-vertex efficiency is evaluated only for truth vertices that have at least 5 associated tracks, an invariant mass $> 10$ GeV, and are in the region $4$ mm $< R_{\mathrm{decay}} < 300$ mm, and $|Z_{\mathrm{decay}}| < 300$ mm. A truth vertex satisfies the vertex level efficiency if it can be matched to a reconstructed DV which passes the final vertex selection.

Parameterized vertex level efficiencies as a function of number of particles associated to a truth decay vertex, and the vertex invariant mass for truth decays with $25$ mm $< R_{\mathrm{decay}} < 29$ mm. Selected particles are required to have nonzero electric charge, $p_{T}(|Q|=1) > 1$ GeV, and $d_0 > 2$ mm. The per-vertex efficiency is evaluated only for truth vertices that have at least 5 associated tracks, an invariant mass $> 10$ GeV, and are in the region $4$ mm $< R_{\mathrm{decay}} < 300$ mm, and $|Z_{\mathrm{decay}}| < 300$ mm. A truth vertex satisfies the vertex level efficiency if it can be matched to a reconstructed DV which passes the final vertex selection.

Parameterized vertex level efficiencies as a function of number of particles associated to a truth decay vertex, and the vertex invariant mass for truth decays with $25$ mm $< R_{\mathrm{decay}} < 29$ mm. Selected particles are required to have nonzero electric charge, $p_{T}(|Q|=1) > 1$ GeV, and $d_0 > 2$ mm. The per-vertex efficiency is evaluated only for truth vertices that have at least 5 associated tracks, an invariant mass $> 10$ GeV, and are in the region $4$ mm $< R_{\mathrm{decay}} < 300$ mm, and $|Z_{\mathrm{decay}}| < 300$ mm. A truth vertex satisfies the vertex level efficiency if it can be matched to a reconstructed DV which passes the final vertex selection.

Upper 95% CL limits on the signal cross section for $m_{\tilde{g}}=1.4$ TeV and fixed $m_{\tilde{\chi}_{1}^{0}}=100$ GeV as a function of lifetime $\tau$.

Parameterized vertex level efficiencies as a function of number of particles associated to a truth decay vertex, and the vertex invariant mass for truth decays with $29$ mm $< R_{\mathrm{decay}} < 38$ mm. Selected particles are required to have nonzero electric charge, $p_{T}(|Q|=1) > 1$ GeV, and $d_0 > 2$ mm. The per-vertex efficiency is evaluated only for truth vertices that have at least 5 associated tracks, an invariant mass $> 10$ GeV, and are in the region $4$ mm $< R_{\mathrm{decay}} < 300$ mm, and $|Z_{\mathrm{decay}}| < 300$ mm. A truth vertex satisfies the vertex level efficiency if it can be matched to a reconstructed DV which passes the final vertex selection.

Parameterized vertex level efficiencies as a function of number of particles associated to a truth decay vertex, and the vertex invariant mass for truth decays with $29$ mm $< R_{\mathrm{decay}} < 38$ mm. Selected particles are required to have nonzero electric charge, $p_{T}(|Q|=1) > 1$ GeV, and $d_0 > 2$ mm. The per-vertex efficiency is evaluated only for truth vertices that have at least 5 associated tracks, an invariant mass $> 10$ GeV, and are in the region $4$ mm $< R_{\mathrm{decay}} < 300$ mm, and $|Z_{\mathrm{decay}}| < 300$ mm. A truth vertex satisfies the vertex level efficiency if it can be matched to a reconstructed DV which passes the final vertex selection.

Upper 95% CL limits on the signal cross section for $m_{\tilde{g}}=2.0$ TeV and fixed $m_{\tilde{\chi}_{1}^{0}}=100$ GeV as a function of lifetime $\tau$.

Parameterized vertex level efficiencies as a function of number of particles associated to a truth decay vertex, and the vertex invariant mass for truth decays with $38$ mm $< R_{\mathrm{decay}} < 46$ mm. Selected particles are required to have nonzero electric charge, $p_{T}(|Q|=1) > 1$ GeV, and $d_0 > 2$ mm. The per-vertex efficiency is evaluated only for truth vertices that have at least 5 associated tracks, an invariant mass $> 10$ GeV, and are in the region $4$ mm $< R_{\mathrm{decay}} < 300$ mm, and $|Z_{\mathrm{decay}}| < 300$ mm. A truth vertex satisfies the vertex level efficiency if it can be matched to a reconstructed DV which passes the final vertex selection.

Parameterized vertex level efficiencies as a function of number of particles associated to a truth decay vertex, and the vertex invariant mass for truth decays with $38$ mm $< R_{\mathrm{decay}} < 46$ mm. Selected particles are required to have nonzero electric charge, $p_{T}(|Q|=1) > 1$ GeV, and $d_0 > 2$ mm. The per-vertex efficiency is evaluated only for truth vertices that have at least 5 associated tracks, an invariant mass $> 10$ GeV, and are in the region $4$ mm $< R_{\mathrm{decay}} < 300$ mm, and $|Z_{\mathrm{decay}}| < 300$ mm. A truth vertex satisfies the vertex level efficiency if it can be matched to a reconstructed DV which passes the final vertex selection.

Lower 95% CL limits on $m_{\tilde{g}}$ for fixed $m_{\tilde{\chi}_{1}^{0}}=100$ GeV as a function of lifetime $\tau$.

Parameterized vertex level efficiencies as a function of number of particles associated to a truth decay vertex, and the vertex invariant mass for truth decays with $46$ mm $< R_{\mathrm{decay}} < 73$ mm. Selected particles are required to have nonzero electric charge, $p_{T}(|Q|=1) > 1$ GeV, and $d_0 > 2$ mm. The per-vertex efficiency is evaluated only for truth vertices that have at least 5 associated tracks, an invariant mass $> 10$ GeV, and are in the region $4$ mm $< R_{\mathrm{decay}} < 300$ mm, and $|Z_{\mathrm{decay}}| < 300$ mm. A truth vertex satisfies the vertex level efficiency if it can be matched to a reconstructed DV which passes the final vertex selection.

Parameterized vertex level efficiencies as a function of number of particles associated to a truth decay vertex, and the vertex invariant mass for truth decays with $46$ mm $< R_{\mathrm{decay}} < 73$ mm. Selected particles are required to have nonzero electric charge, $p_{T}(|Q|=1) > 1$ GeV, and $d_0 > 2$ mm. The per-vertex efficiency is evaluated only for truth vertices that have at least 5 associated tracks, an invariant mass $> 10$ GeV, and are in the region $4$ mm $< R_{\mathrm{decay}} < 300$ mm, and $|Z_{\mathrm{decay}}| < 300$ mm. A truth vertex satisfies the vertex level efficiency if it can be matched to a reconstructed DV which passes the final vertex selection.

Parameterized vertex level efficiencies as a function of number of particles associated to a truth decay vertex, and the vertex invariant mass for truth decays with $73$ mm $< R_{\mathrm{decay}} < 84$ mm. Selected particles are required to have nonzero electric charge, $p_{T}(|Q|=1) > 1$ GeV, and $d_0 > 2$ mm. The per-vertex efficiency is evaluated only for truth vertices that have at least 5 associated tracks, an invariant mass $> 10$ GeV, and are in the region $4$ mm $< R_{\mathrm{decay}} < 300$ mm, and $|Z_{\mathrm{decay}}| < 300$ mm. A truth vertex satisfies the vertex level efficiency if it can be matched to a reconstructed DV which passes the final vertex selection.

Parameterized vertex level efficiencies as a function of number of particles associated to a truth decay vertex, and the vertex invariant mass for truth decays with $73$ mm $< R_{\mathrm{decay}} < 84$ mm. Selected particles are required to have nonzero electric charge, $p_{T}(|Q|=1) > 1$ GeV, and $d_0 > 2$ mm. The per-vertex efficiency is evaluated only for truth vertices that have at least 5 associated tracks, an invariant mass $> 10$ GeV, and are in the region $4$ mm $< R_{\mathrm{decay}} < 300$ mm, and $|Z_{\mathrm{decay}}| < 300$ mm. A truth vertex satisfies the vertex level efficiency if it can be matched to a reconstructed DV which passes the final vertex selection.

Parameterized vertex level efficiencies as a function of number of particles associated to a truth decay vertex, and the vertex invariant mass for truth decays with $84$ mm $< R_{\mathrm{decay}} < 111$ mm. Selected particles are required to have nonzero electric charge, $p_{T}(|Q|=1) > 1$ GeV, and $d_0 > 2$ mm. The per-vertex efficiency is evaluated only for truth vertices that have at least 5 associated tracks, an invariant mass $> 10$ GeV, and are in the region $4$ mm $< R_{\mathrm{decay}} < 300$ mm, and $|Z_{\mathrm{decay}}| < 300$ mm. A truth vertex satisfies the vertex level efficiency if it can be matched to a reconstructed DV which passes the final vertex selection.

Parameterized vertex level efficiencies as a function of number of particles associated to a truth decay vertex, and the vertex invariant mass for truth decays with $84$ mm $< R_{\mathrm{decay}} < 111$ mm. Selected particles are required to have nonzero electric charge, $p_{T}(|Q|=1) > 1$ GeV, and $d_0 > 2$ mm. The per-vertex efficiency is evaluated only for truth vertices that have at least 5 associated tracks, an invariant mass $> 10$ GeV, and are in the region $4$ mm $< R_{\mathrm{decay}} < 300$ mm, and $|Z_{\mathrm{decay}}| < 300$ mm. A truth vertex satisfies the vertex level efficiency if it can be matched to a reconstructed DV which passes the final vertex selection.

Parameterized vertex level efficiencies as a function of number of particles associated to a truth decay vertex, and the vertex invariant mass for truth decays with $111$ mm $< R_{\mathrm{decay}} < 120$ mm. Selected particles are required to have nonzero electric charge, $p_{T}(|Q|=1) > 1$ GeV, and $d_0 > 2$ mm. The per-vertex efficiency is evaluated only for truth vertices that have at least 5 associated tracks, an invariant mass $> 10$ GeV, and are in the region $4$ mm $< R_{\mathrm{decay}} < 300$ mm, and $|Z_{\mathrm{decay}}| < 300$ mm. A truth vertex satisfies the vertex level efficiency if it can be matched to a reconstructed DV which passes the final vertex selection.

Parameterized vertex level efficiencies as a function of number of particles associated to a truth decay vertex, and the vertex invariant mass for truth decays with $111$ mm $< R_{\mathrm{decay}} < 120$ mm. Selected particles are required to have nonzero electric charge, $p_{T}(|Q|=1) > 1$ GeV, and $d_0 > 2$ mm. The per-vertex efficiency is evaluated only for truth vertices that have at least 5 associated tracks, an invariant mass $> 10$ GeV, and are in the region $4$ mm $< R_{\mathrm{decay}} < 300$ mm, and $|Z_{\mathrm{decay}}| < 300$ mm. A truth vertex satisfies the vertex level efficiency if it can be matched to a reconstructed DV which passes the final vertex selection.

Upper 95% CL limits on the signal cross section for $m_{\tilde{g}}=1.4$ TeV and fixed $\tau=1$ ns as a function of $m_{\tilde{\chi}_{1}^{0}}$.

Parameterized vertex level efficiencies as a function of number of particles associated to a truth decay vertex, and the vertex invariant mass for truth decays with $120$ mm $< R_{\mathrm{decay}} < 145$ mm. Selected particles are required to have nonzero electric charge, $p_{T}(|Q|=1) > 1$ GeV, and $d_0 > 2$ mm. The per-vertex efficiency is evaluated only for truth vertices that have at least 5 associated tracks, an invariant mass $> 10$ GeV, and are in the region $4$ mm $< R_{\mathrm{decay}} < 300$ mm, and $|Z_{\mathrm{decay}}| < 300$ mm. A truth vertex satisfies the vertex level efficiency if it can be matched to a reconstructed DV which passes the final vertex selection.

Parameterized vertex level efficiencies as a function of number of particles associated to a truth decay vertex, and the vertex invariant mass for truth decays with $120$ mm $< R_{\mathrm{decay}} < 145$ mm. Selected particles are required to have nonzero electric charge, $p_{T}(|Q|=1) > 1$ GeV, and $d_0 > 2$ mm. The per-vertex efficiency is evaluated only for truth vertices that have at least 5 associated tracks, an invariant mass $> 10$ GeV, and are in the region $4$ mm $< R_{\mathrm{decay}} < 300$ mm, and $|Z_{\mathrm{decay}}| < 300$ mm. A truth vertex satisfies the vertex level efficiency if it can be matched to a reconstructed DV which passes the final vertex selection.

Upper 95% CL limits on the signal cross section for $m_{\tilde{g}}=2.0$ TeV and fixed $\tau=1$ ns as a function of $m_{\tilde{\chi}_{1}^{0}}$.

Parameterized vertex level efficiencies as a function of number of particles associated to a truth decay vertex, and the vertex invariant mass for truth decays with $145$ mm $< R_{\mathrm{decay}} < 180$ mm. Selected particles are required to have nonzero electric charge, $p_{T}(|Q|=1) > 1$ GeV, and $d_0 > 2$ mm. The per-vertex efficiency is evaluated only for truth vertices that have at least 5 associated tracks, an invariant mass $> 10$ GeV, and are in the region $4$ mm $< R_{\mathrm{decay}} < 300$ mm, and $|Z_{\mathrm{decay}}| < 300$ mm. A truth vertex satisfies the vertex level efficiency if it can be matched to a reconstructed DV which passes the final vertex selection.

Parameterized vertex level efficiencies as a function of number of particles associated to a truth decay vertex, and the vertex invariant mass for truth decays with $145$ mm $< R_{\mathrm{decay}} < 180$ mm. Selected particles are required to have nonzero electric charge, $p_{T}(|Q|=1) > 1$ GeV, and $d_0 > 2$ mm. The per-vertex efficiency is evaluated only for truth vertices that have at least 5 associated tracks, an invariant mass $> 10$ GeV, and are in the region $4$ mm $< R_{\mathrm{decay}} < 300$ mm, and $|Z_{\mathrm{decay}}| < 300$ mm. A truth vertex satisfies the vertex level efficiency if it can be matched to a reconstructed DV which passes the final vertex selection.

Observed 95% CL limit as a function of $m_{\tilde{g}}$ and $m_{\tilde{\chi}_{1}^{0}}$ for fixed $\tau=1$ ns.

Parameterized vertex level efficiencies as a function of number of particles associated to a truth decay vertex, and the vertex invariant mass for truth decays with $180$ mm $< R_{\mathrm{decay}} < 300$ mm. Selected particles are required to have nonzero electric charge, $p_{T}(|Q|=1) > 1$ GeV, and $d_0 > 2$ mm. The per-vertex efficiency is evaluated only for truth vertices that have at least 5 associated tracks, an invariant mass $> 10$ GeV, and are in the region $4$ mm $< R_{\mathrm{decay}} < 300$ mm, and $|Z_{\mathrm{decay}}| < 300$ mm. A truth vertex satisfies the vertex level efficiency if it can be matched to a reconstructed DV which passes the final vertex selection.

Parameterized vertex level efficiencies as a function of number of particles associated to a truth decay vertex, and the vertex invariant mass for truth decays with $180$ mm $< R_{\mathrm{decay}} < 300$ mm. Selected particles are required to have nonzero electric charge, $p_{T}(|Q|=1) > 1$ GeV, and $d_0 > 2$ mm. The per-vertex efficiency is evaluated only for truth vertices that have at least 5 associated tracks, an invariant mass $> 10$ GeV, and are in the region $4$ mm $< R_{\mathrm{decay}} < 300$ mm, and $|Z_{\mathrm{decay}}| < 300$ mm. A truth vertex satisfies the vertex level efficiency if it can be matched to a reconstructed DV which passes the final vertex selection.

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

The observed exclusion contour at 95% CL as a function of the $m_{H^{\pm}}$ vs. tan$\beta$ and assuming $m_a$ = 150 $\mathrm{GeV}$, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.7$. Masses that are within the contours are excluded. Only signals simulating the tW+DM final states are considered in this contour.

The expected exclusion contour at 95% CL as a function of the $m_{H^{\pm}}$ vs. tan$\beta$ and assuming $m_a$ = 150 $\mathrm{GeV}$, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.7$. Masses that are within the contours are excluded. Only signals simulating the tW+DM final states are considered in this contour.

The observed exclusion contour at 95% CL as a function of the $m_{H^{\pm}}$ vs. tan$\beta$ and assuming $m_a$ = 250 $\mathrm{GeV}$, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.7$. Masses that are within the contours are excluded. Only signals simulating the tW+DM final states are considered in this contour.

The expected exclusion contour at 95% CL as a function of the $m_{H^{\pm}}$ vs. tan$\beta$ and assuming $m_a$ = 250 $\mathrm{GeV}$, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.7$. Masses that are within the contours are excluded. Only signals simulating the tW+DM final states are considered in this contour.

The observed exclusion contour at 95% CL as a function of the $m_a$ vs. $m_{H^{\pm}}$ and assuming tan$\beta$ = 1, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.7$. Masses that are within the contours are excluded. Signals simulating the tW+DM + tt+DM final states are considered in this contour.

The expected exclusion contour at 95% CL as a function of the $m_a$ vs. $m_{H^{\pm}}$ and assuming tan$\beta$ = 1, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.7$. Masses that are within the contours are excluded. Signals simulating the tW+DM + tt+DM final states are considered in this contour.

The observed exclusion contour at 95% CL as a function of the $m_{H^{\pm}}$ vs. tan$\beta$ and assuming $m_a$ = 150 $\mathrm{GeV}$, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.7$. Masses that are within the contours are excluded. Signals simulating the tW+DM + tt+DM final states are considered in this contour.

The expected exclusion contour at 95% CL as a function of the $m_{H^{\pm}}$ vs. tan$\beta$ and assuming $m_a$ = 150 $\mathrm{GeV}$, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.7$. Masses that are within the contours are excluded. Signals simulating the tW+DM + tt+DM final states are considered in this contour.

The observed exclusion contour at 95% CL as a function of the $m_{H^{\pm}}$ vs. tan$\beta$ and assuming $m_a$ = 250 $\mathrm{GeV}$, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.7$. Masses that are within the contours are excluded. Signals simulating the tW+DM + tt+DM final states are considered in this contour.

The expected exclusion contour at 95% CL as a function of the $m_{H^{\pm}}$ vs. tan$\beta$ and assuming $m_a$ = 250 $\mathrm{GeV}$, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.7$. Masses that are within the contours are excluded. Signals simulating the tW+DM + tt+DM final states are considered in this contour.

The observed exclusion contour at 95% CL as a function of the $m_a$ vs. $m_{H^{\pm}}$ and assuming tan$\beta$ = 1, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.35$. Masses that are within the contours are excluded. Only signals simulating the tW+DM final states are considered in this contour.

The expected exclusion contour at 95% CL as a function of the $m_a$ vs. $m_{H^{\pm}}$ and assuming tan$\beta$ = 1, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.35$. Masses that are within the contours are excluded. Only signals simulating the tW+DM final states are considered in this contour.

The observed exclusion contour at 95% CL as a function of the $m_{H^{\pm}}$ vs. tan$\beta$ and assuming $m_a$ = 150 $\mathrm{GeV}$, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.35$. Masses that are within the contours are excluded. Only signals simulating the tW+DM final states are considered in this contour.

The expected exclusion contour at 95% CL as a function of the $m_{H^{\pm}}$ vs. tan$\beta$ and assuming $m_a$ = 150 $\mathrm{GeV}$, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.35$. Masses that are within the contours are excluded. Only signals simulating the tW+DM final states are considered in this contour.

The observed exclusion contour at 95% CL as a function of the $m_{H^{\pm}}$ vs. tan$\beta$ and assuming $m_a$ = 250 $\mathrm{GeV}$, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.35$. Masses that are within the contours are excluded. Only signals simulating the tW+DM final states are considered in this contour.

The expected exclusion contour at 95% CL as a function of the $m_{H^{\pm}}$ vs. tan$\beta$ and assuming $m_a$ = 250 $\mathrm{GeV}$, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.35$. Masses that are within the contours are excluded. Only signals simulating the tW+DM final states are considered in this contour.

The observed exclusion contour at 95% CL as a function of the $m_a$ vs. $m_{H^{\pm}}$ and assuming tan$\beta$ = 1, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.35$. Masses that are within the contours are excluded. Signals simulating the tW+DM + tt+DM final states are considered in this contour.

The expected exclusion contour at 95% CL as a function of the $m_a$ vs. $m_{H^{\pm}}$ and assuming tan$\beta$ = 1, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.35$. Masses that are within the contours are excluded. Signals simulating the tW+DM + tt+DM final states are considered in this contour.

The observed exclusion contour at 95% CL as a function of the $m_{H^{\pm}}$ vs. tan$\beta$ and assuming $m_a$ = 150 $\mathrm{GeV}$, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.35$. Masses that are within the contours are excluded. Signals simulating the tW+DM + tt+DM final states are considered in this contour.

The expected exclusion contour at 95% CL as a function of the $m_{H^{\pm}}$ vs. tan$\beta$ and assuming $m_a$ = 150 $\mathrm{GeV}$, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.35$. Masses that are within the contours are excluded. Signals simulating the tW+DM + tt+DM final states are considered in this contour.

The observed exclusion contour at 95% CL as a function of the $m_{H^{\pm}}$ vs. tan$\beta$ and assuming $m_a$ = 250 $\mathrm{GeV}$, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.35$. Masses that are within the contours are excluded. Signals simulating the tW+DM + tt+DM final states are considered in this contour.

The expected exclusion contour at 95% CL as a function of the $m_{H^{\pm}}$ vs. tan$\beta$ and assuming $m_a$ = 250 $\mathrm{GeV}$, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.35$. Masses that are within the contours are excluded. Signals simulating the tW+DM + tt+DM final states are considered in this contour.

The observed exclusion contour at 95% CL as a function of the $m_a$ vs. $m_{H^{\pm}}$ and assuming tan$\beta$ = 1, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.7$. Masses that are within the contours are excluded. Only signals simulating the tW+DM final states are considered in this contour. These exclusion contours are derived using the 0L channel only.

The expected exclusion contour at 95% CL as a function of the $m_a$ vs. $m_{H^{\pm}}$ and assuming tan$\beta$ = 1, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.7$. Masses that are within the contours are excluded. Only signals simulating the tW+DM final states are considered in this contour. These exclusion contours are derived using the 0L channel only.

The observed exclusion contour at 95% CL as a function of the $m_{H^{\pm}}$ vs. tan$\beta$ and assuming $m_a$ = 150 $\mathrm{GeV}$, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.7$. Masses that are within the contours are excluded. Only signals simulating the tW+DM final states are considered in this contour. These exclusion contours are derived using the 0L channel only.

The expected exclusion contour at 95% CL as a function of the $m_{H^{\pm}}$ vs. tan$\beta$ and assuming $m_a$ = 150 $\mathrm{GeV}$, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.7$. Masses that are within the contours are excluded. Only signals simulating the tW+DM final states are considered in this contour. These exclusion contours are derived using the 0L channel only.

The observed exclusion contour at 95% CL as a function of the $m_{H^{\pm}}$ vs. tan$\beta$ and assuming $m_a$ = 250 $\mathrm{GeV}$, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.7$. Masses that are within the contours are excluded. Only signals simulating the tW+DM final states are considered in this contour. These exclusion contours are derived using the 0L channel only.

The expected exclusion contour at 95% CL as a function of the $m_{H^{\pm}}$ vs. tan$\beta$ and assuming $m_a$ = 250 $\mathrm{GeV}$, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.7$. Masses that are within the contours are excluded. Only signals simulating the tW+DM final states are considered in this contour. These exclusion contours are derived using the 0L channel only.

The observed exclusion contour at 95% CL as a function of the $m_a$ vs. $m_{H^{\pm}}$ and assuming tan$\beta$ = 1, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.7$. Masses that are within the contours are excluded. Signals simulating the tW+DM + tt+DM final states are considered in this contour. These exclusion contours are derived using the 0L channel only.

The expected exclusion contour at 95% CL as a function of the $m_a$ vs. $m_{H^{\pm}}$ and assuming tan$\beta$ = 1, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.7$. Masses that are within the contours are excluded. Signals simulating the tW+DM + tt+DM final states are considered in this contour. These exclusion contours are derived using the 0L channel only.

The observed exclusion contour at 95% CL as a function of the $m_{H^{\pm}}$ vs. tan$\beta$ and assuming $m_a$ = 150 $\mathrm{GeV}$, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.7$. Masses that are within the contours are excluded. Signals simulating the tW+DM + tt+DM final states are considered in this contour. These exclusion contours are derived using the 0L channel only.

The expected exclusion contour at 95% CL as a function of the $m_{H^{\pm}}$ vs. tan$\beta$ and assuming $m_a$ = 150 $\mathrm{GeV}$, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.7$. Masses that are within the contours are excluded. Signals simulating the tW+DM + tt+DM final states are considered in this contour. These exclusion contours are derived using the 0L channel only.

The observed exclusion contour at 95% CL as a function of the $m_{H^{\pm}}$ vs. tan$\beta$ and assuming $m_a$ = 250 $\mathrm{GeV}$, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.7$. Masses that are within the contours are excluded. Signals simulating the tW+DM + tt+DM final states are considered in this contour. These exclusion contours are derived using the 0L channel only.

The expected exclusion contour at 95% CL as a function of the $m_{H^{\pm}}$ vs. tan$\beta$ and assuming $m_a$ = 250 $\mathrm{GeV}$, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.7$. Masses that are within the contours are excluded. Signals simulating the tW+DM + tt+DM final states are considered in this contour. These exclusion contours are derived using the 0L channel only.

The observed exclusion contour at 95% CL as a function of the $m_a$ vs. $m_{H^{\pm}}$ and assuming tan$\beta$ = 1, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.35$. Masses that are within the contours are excluded. Only signals simulating the tW+DM final states are considered in this contour. These exclusion contours are derived using the 0L channel only.

The expected exclusion contour at 95% CL as a function of the $m_a$ vs. $m_{H^{\pm}}$ and assuming tan$\beta$ = 1, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.35$. Masses that are within the contours are excluded. Only signals simulating the tW+DM final states are considered in this contour. These exclusion contours are derived using the 0L channel only.

The observed exclusion contour at 95% CL as a function of the $m_{H^{\pm}}$ vs. tan$\beta$ and assuming $m_a$ = 150 $\mathrm{GeV}$, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.35$. Masses that are within the contours are excluded. Only signals simulating the tW+DM final states are considered in this contour. These exclusion contours are derived using the 0L channel only.

The expected exclusion contour at 95% CL as a function of the $m_{H^{\pm}}$ vs. tan$\beta$ and assuming $m_a$ = 150 $\mathrm{GeV}$, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.35$. Masses that are within the contours are excluded. Only signals simulating the tW+DM final states are considered in this contour. These exclusion contours are derived using the 0L channel only.

The observed exclusion contour at 95% CL as a function of the $m_{H^{\pm}}$ vs. tan$\beta$ and assuming $m_a$ = 250 $\mathrm{GeV}$, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.35$. Masses that are within the contours are excluded. Only signals simulating the tW+DM final states are considered in this contour. These exclusion contours are derived using the 0L channel only.

The expected exclusion contour at 95% CL as a function of the $m_{H^{\pm}}$ vs. tan$\beta$ and assuming $m_a$ = 250 $\mathrm{GeV}$, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.35$. Masses that are within the contours are excluded. Only signals simulating the tW+DM final states are considered in this contour. These exclusion contours are derived using the 0L channel only.

The observed exclusion contour at 95% CL as a function of the $m_a$ vs. $m_{H^{\pm}}$ and assuming tan$\beta$ = 1, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.35$. Masses that are within the contours are excluded. Signals simulating the tW+DM + tt+DM final states are considered in this contour. These exclusion contours are derived using the 0L channel only.

The expected exclusion contour at 95% CL as a function of the $m_a$ vs. $m_{H^{\pm}}$ and assuming tan$\beta$ = 1, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.35$. Masses that are within the contours are excluded. Signals simulating the tW+DM + tt+DM final states are considered in this contour. These exclusion contours are derived using the 0L channel only.

The observed exclusion contour at 95% CL as a function of the $m_{H^{\pm}}$ vs. tan$\beta$ and assuming $m_a$ = 150 $\mathrm{GeV}$, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.35$. Masses that are within the contours are excluded. Signals simulating the tW+DM + tt+DM final states are considered in this contour. These exclusion contours are derived using the 0L channel only.

The expected exclusion contour at 95% CL as a function of the $m_{H^{\pm}}$ vs. tan$\beta$ and assuming $m_a$ = 150 $\mathrm{GeV}$, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.35$. Masses that are within the contours are excluded. Signals simulating the tW+DM + tt+DM final states are considered in this contour. These exclusion contours are derived using the 0L channel only.

The observed exclusion contour at 95% CL as a function of the $m_{H^{\pm}}$ vs. tan$\beta$ and assuming $m_a$ = 250 $\mathrm{GeV}$, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.35$. Masses that are within the contours are excluded. Signals simulating the tW+DM + tt+DM final states are considered in this contour. These exclusion contours are derived using the 0L channel only.

The expected exclusion contour at 95% CL as a function of the $m_{H^{\pm}}$ vs. tan$\beta$ and assuming $m_a$ = 250 $\mathrm{GeV}$, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.35$. Masses that are within the contours are excluded. Signals simulating the tW+DM + tt+DM final states are considered in this contour. These exclusion contours are derived using the 0L channel only.

The observed exclusion contour at 95% CL as a function of the $m_a$ vs. $m_{H^{\pm}}$ and assuming tan$\beta$ = 1, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.7$. Masses that are within the contours are excluded. Only signals simulating the tW+DM final states are considered in this contour. These exclusion contours are derived using the 1L channel only.

The expected exclusion contour at 95% CL as a function of the $m_a$ vs. $m_{H^{\pm}}$ and assuming tan$\beta$ = 1, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.7$. Masses that are within the contours are excluded. Only signals simulating the tW+DM final states are considered in this contour. These exclusion contours are derived using the 1L channel only.

The observed exclusion contour at 95% CL as a function of the $m_{H^{\pm}}$ vs. tan$\beta$ and assuming $m_a$ = 150 $\mathrm{GeV}$, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.7$. Masses that are within the contours are excluded. Only signals simulating the tW+DM final states are considered in this contour. These exclusion contours are derived using the 1L channel only.

The expected exclusion contour at 95% CL as a function of the $m_{H^{\pm}}$ vs. tan$\beta$ and assuming $m_a$ = 150 $\mathrm{GeV}$, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.7$. Masses that are within the contours are excluded. Only signals simulating the tW+DM final states are considered in this contour. These exclusion contours are derived using the 1L channel only.

The observed exclusion contour at 95% CL as a function of the $m_{H^{\pm}}$ vs. tan$\beta$ and assuming $m_a$ = 250 $\mathrm{GeV}$, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.7$. Masses that are within the contours are excluded. Only signals simulating the tW+DM final states are considered in this contour. These exclusion contours are derived using the 1L channel only.

The expected exclusion contour at 95% CL as a function of the $m_{H^{\pm}}$ vs. tan$\beta$ and assuming $m_a$ = 250 $\mathrm{GeV}$, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.7$. Masses that are within the contours are excluded. Only signals simulating the tW+DM final states are considered in this contour. These exclusion contours are derived using the 1L channel only.

The observed exclusion contour at 95% CL as a function of the $m_a$ vs. $m_{H^{\pm}}$ and assuming tan$\beta$ = 1, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.7$. Masses that are within the contours are excluded. Signals simulating the tW+DM + tt+DM final states are considered in this contour. These exclusion contours are derived using the 1L channel only.

The expected exclusion contour at 95% CL as a function of the $m_a$ vs. $m_{H^{\pm}}$ and assuming tan$\beta$ = 1, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.7$. Masses that are within the contours are excluded. Signals simulating the tW+DM + tt+DM final states are considered in this contour. These exclusion contours are derived using the 1L channel only.

The observed exclusion contour at 95% CL as a function of the $m_{H^{\pm}}$ vs. tan$\beta$ and assuming $m_a$ = 150 $\mathrm{GeV}$, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.7$. Masses that are within the contours are excluded. Signals simulating the tW+DM + tt+DM final states are considered in this contour. These exclusion contours are derived using the 1L channel only.

The expected exclusion contour at 95% CL as a function of the $m_{H^{\pm}}$ vs. tan$\beta$ and assuming $m_a$ = 150 $\mathrm{GeV}$, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.7$. Masses that are within the contours are excluded. Signals simulating the tW+DM + tt+DM final states are considered in this contour. These exclusion contours are derived using the 1L channel only.

The observed exclusion contour at 95% CL as a function of the $m_{H^{\pm}}$ vs. tan$\beta$ and assuming $m_a$ = 250 $\mathrm{GeV}$, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.7$. Masses that are within the contours are excluded. Signals simulating the tW+DM + tt+DM final states are considered in this contour. These exclusion contours are derived using the 1L channel only.

The expected exclusion contour at 95% CL as a function of the $m_{H^{\pm}}$ vs. tan$\beta$ and assuming $m_a$ = 250 $\mathrm{GeV}$, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.7$. Masses that are within the contours are excluded. Signals simulating the tW+DM + tt+DM final states are considered in this contour. These exclusion contours are derived using the 1L channel only.

The observed exclusion contour at 95% CL as a function of the $m_a$ vs. $m_{H^{\pm}}$ and assuming tan$\beta$ = 1, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.35$. Masses that are within the contours are excluded. Only signals simulating the tW+DM final states are considered in this contour. These exclusion contours are derived using the 1L channel only.

The expected exclusion contour at 95% CL as a function of the $m_a$ vs. $m_{H^{\pm}}$ and assuming tan$\beta$ = 1, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.35$. Masses that are within the contours are excluded. Only signals simulating the tW+DM final states are considered in this contour. These exclusion contours are derived using the 1L channel only.

The observed exclusion contour at 95% CL as a function of the $m_{H^{\pm}}$ vs. tan$\beta$ and assuming $m_a$ = 150 $\mathrm{GeV}$, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.35$. Masses that are within the contours are excluded. Only signals simulating the tW+DM final states are considered in this contour. These exclusion contours are derived using the 1L channel only.

The expected exclusion contour at 95% CL as a function of the $m_{H^{\pm}}$ vs. tan$\beta$ and assuming $m_a$ = 150 $\mathrm{GeV}$, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.35$. Masses that are within the contours are excluded. Only signals simulating the tW+DM final states are considered in this contour. These exclusion contours are derived using the 1L channel only.

The expected exclusion contour at 95% CL as a function of the $m_{H^{\pm}}$ vs. tan$\beta$ and assuming $m_a$ = 250 $\mathrm{GeV}$, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.35$. Masses that are within the contours are excluded. Only signals simulating the tW+DM final states are considered in this contour. These exclusion contours are derived using the 1L channel only.

The observed exclusion contour at 95% CL as a function of the $m_a$ vs. $m_{H^{\pm}}$ and assuming tan$\beta$ = 1, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.35$. Masses that are within the contours are excluded. Signals simulating the tW+DM + tt+DM final states are considered in this contour. These exclusion contours are derived using the 1L channel only.

The expected exclusion contour at 95% CL as a function of the $m_a$ vs. $m_{H^{\pm}}$ and assuming tan$\beta$ = 1, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.35$. Masses that are within the contours are excluded. Signals simulating the tW+DM + tt+DM final states are considered in this contour. These exclusion contours are derived using the 1L channel only.

The observed exclusion contour at 95% CL as a function of the $m_{H^{\pm}}$ vs. tan$\beta$ and assuming $m_a$ = 150 $\mathrm{GeV}$, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.35$. Masses that are within the contours are excluded. Signals simulating the tW+DM + tt+DM final states are considered in this contour. These exclusion contours are derived using the 1L channel only.

The expected exclusion contour at 95% CL as a function of the $m_{H^{\pm}}$ vs. tan$\beta$ and assuming $m_a$ = 150 $\mathrm{GeV}$, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.35$. Masses that are within the contours are excluded. Signals simulating the tW+DM + tt+DM final states are considered in this contour. These exclusion contours are derived using the 1L channel only.

The observed exclusion contour at 95% CL as a function of the $m_{H^{\pm}}$ vs. tan$\beta$ and assuming $m_a$ = 250 $\mathrm{GeV}$, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.35$. Masses that are within the contours are excluded. Signals simulating the tW+DM + tt+DM final states are considered in this contour. These exclusion contours are derived using the 1L channel only.

The expected exclusion contour at 95% CL as a function of the $m_{H^{\pm}}$ vs. tan$\beta$ and assuming $m_a$ = 250 $\mathrm{GeV}$, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.35$. Masses that are within the contours are excluded. Signals simulating the tW+DM + tt+DM final states are considered in this contour. These exclusion contours are derived using the 1L channel only.

The observed exclusion contour at 95% CL as a function of the $m_a$ vs. $m_{H^{\pm}}$ and assuming tan$\beta$ = 1, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.7$. Masses that are within the contours are excluded. Only signals simulating the tW+DM final states are considered in this contour. These exclusion contours are derived using the 2L channel only.

The expected exclusion contour at 95% CL as a function of the $m_a$ vs. $m_{H^{\pm}}$ and assuming tan$\beta$ = 1, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.7$. Masses that are within the contours are excluded. Only signals simulating the tW+DM final states are considered in this contour. These exclusion contours are derived using the 2L channel only.

The observed exclusion contour at 95% CL as a function of the $m_{H^{\pm}}$ vs. tan$\beta$ and assuming $m_a$ = 150 $\mathrm{GeV}$, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.7$. Masses that are within the contours are excluded. Only signals simulating the tW+DM final states are considered in this contour. These exclusion contours are derived using the 2L channel only.

The expected exclusion contour at 95% CL as a function of the $m_{H^{\pm}}$ vs. tan$\beta$ and assuming $m_a$ = 150 $\mathrm{GeV}$, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.7$. Masses that are within the contours are excluded. Only signals simulating the tW+DM final states are considered in this contour. These exclusion contours are derived using the 2L channel only.

The observed exclusion contour at 95% CL as a function of the $m_{H^{\pm}}$ vs. tan$\beta$ and assuming $m_a$ = 250 $\mathrm{GeV}$, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.7$. Masses that are within the contours are excluded. Only signals simulating the tW+DM final states are considered in this contour. These exclusion contours are derived using the 2L channel only.

The expected exclusion contour at 95% CL as a function of the $m_{H^{\pm}}$ vs. tan$\beta$ and assuming $m_a$ = 250 $\mathrm{GeV}$, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.7$. Masses that are within the contours are excluded. Only signals simulating the tW+DM final states are considered in this contour. These exclusion contours are derived using the 2L channel only.

The observed exclusion contour at 95% CL as a function of the $m_a$ vs. $m_{H^{\pm}}$ and assuming tan$\beta$ = 1, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.7$. Masses that are within the contours are excluded. Signals simulating the tW+DM + tt+DM final states are considered in this contour. These exclusion contours are derived using the 2L channel only.

The expected exclusion contour at 95% CL as a function of the $m_a$ vs. $m_{H^{\pm}}$ and assuming tan$\beta$ = 1, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.7$. Masses that are within the contours are excluded. Signals simulating the tW+DM + tt+DM final states are considered in this contour. These exclusion contours are derived using the 2L channel only.

The observed exclusion contour at 95% CL as a function of the $m_{H^{\pm}}$ vs. tan$\beta$ and assuming $m_a$ = 150 $\mathrm{GeV}$, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.7$. Masses that are within the contours are excluded. Signals simulating the tW+DM + tt+DM final states are considered in this contour. These exclusion contours are derived using the 2L channel only.

The expected exclusion contour at 95% CL as a function of the $m_{H^{\pm}}$ vs. tan$\beta$ and assuming $m_a$ = 150 $\mathrm{GeV}$, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.7$. Masses that are within the contours are excluded. Signals simulating the tW+DM + tt+DM final states are considered in this contour. These exclusion contours are derived using the 2L channel only.

The observed exclusion contour at 95% CL as a function of the $m_{H^{\pm}}$ vs. tan$\beta$ and assuming $m_a$ = 250 $\mathrm{GeV}$, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.7$. Masses that are within the contours are excluded. Signals simulating the tW+DM + tt+DM final states are considered in this contour. These exclusion contours are derived using the 2L channel only.

The expected exclusion contour at 95% CL as a function of the $m_{H^{\pm}}$ vs. tan$\beta$ and assuming $m_a$ = 250 $\mathrm{GeV}$, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.7$. Masses that are within the contours are excluded. Signals simulating the tW+DM + tt+DM final states are considered in this contour. These exclusion contours are derived using the 2L channel only.

The expected exclusion contour at 95% CL as a function of the $m_a$ vs. $m_{H^{\pm}}$ and assuming tan$\beta$ = 1, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.35$. Masses that are within the contours are excluded. Only signals simulating the tW+DM final states are considered in this contour. These exclusion contours are derived using the 2L channel only.

The expected exclusion contour at 95% CL as a function of the $m_{H^{\pm}}$ vs. tan$\beta$ and assuming $m_a$ = 150 $\mathrm{GeV}$, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.35$. Masses that are within the contours are excluded. Only signals simulating the tW+DM final states are considered in this contour. These exclusion contours are derived using the 2L channel only.

The expected exclusion contour at 95% CL as a function of the $m_{H^{\pm}}$ vs. tan$\beta$ and assuming $m_a$ = 250 $\mathrm{GeV}$, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.35$. Masses that are within the contours are excluded. Only signals simulating the tW+DM final states are considered in this contour. These exclusion contours are derived using the 2L channel only.

The expected exclusion contour at 95% CL as a function of the $m_a$ vs. $m_{H^{\pm}}$ and assuming tan$\beta$ = 1, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.35$. Masses that are within the contours are excluded. Signals simulating the tW+DM + tt+DM final states are considered in this contour. These exclusion contours are derived using the 2L channel only.

The observed exclusion contour at 95% CL as a function of the $m_{H^{\pm}}$ vs. tan$\beta$ and assuming $m_a$ = 150 $\mathrm{GeV}$, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.35$. Masses that are within the contours are excluded. Signals simulating the tW+DM + tt+DM final states are considered in this contour. These exclusion contours are derived using the 2L channel only.

The expected exclusion contour at 95% CL as a function of the $m_{H^{\pm}}$ vs. tan$\beta$ and assuming $m_a$ = 150 $\mathrm{GeV}$, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.35$. Masses that are within the contours are excluded. Signals simulating the tW+DM + tt+DM final states are considered in this contour. These exclusion contours are derived using the 2L channel only.

The observed exclusion contour at 95% CL as a function of the $m_{H^{\pm}}$ vs. tan$\beta$ and assuming $m_a$ = 250 $\mathrm{GeV}$, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.35$. Masses that are within the contours are excluded. Signals simulating the tW+DM + tt+DM final states are considered in this contour. These exclusion contours are derived using the 2L channel only.

The expected exclusion contour at 95% CL as a function of the $m_{H^{\pm}}$ vs. tan$\beta$ and assuming $m_a$ = 250 $\mathrm{GeV}$, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.35$. Masses that are within the contours are excluded. Signals simulating the tW+DM + tt+DM final states are considered in this contour. These exclusion contours are derived using the 2L channel only.

Model dependent upper limit on the cross section for the $m_a$ vs. $ m_{H^{\pm}}$ signal grid assuming tan$\beta$ = 1, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.7$. Only signals simulating the tW+DM final states are considered. Upper limits with large $\mu_{\mathrm{sig}}$ for the observed limit are capped at 500.

Model dependent upper limit on the cross section for the $m_{H^{\pm}}$ vs. tan$\beta$ signal grid assuming $m_a$ = 150 GeV, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.7$. Only signals simulating the tW+DM final states are considered. Upper limits with large $\mu_{\mathrm{sig}}$ for the observed limit are capped at 500.

Model dependent upper limit on the cross section for the $m_{H^{\pm}}$ vs. tan$\beta$ signal grid assuming $m_a$ = 250 GeV, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.7$. Only signals simulating the tW+DM final states are considered. Upper limits with large $\mu_{\mathrm{sig}}$ for the observed limit are capped at 500.

Model dependent upper limit on the cross section for the $m_a$ vs. $ m_{H^{\pm}}$ signal grid assuming tan$\beta$ = 1, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.7$. Signals simulating the tW+DM + tt+DM final states are considered. Upper limits with large $\mu_{\mathrm{sig}}$ for the observed limit are capped at 500.

Model dependent upper limit on the cross section for the $m_{H^{\pm}}$ vs. tan$\beta$ signal grid assuming $m_a$ = 150 GeV, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.7$. Signals simulating the tW+DM + tt+DM final states are considered. Upper limits with large $\mu_{\mathrm{sig}}$ for the observed limit are capped at 500.

Model dependent upper limit on the cross section for the $m_{H^{\pm}}$ vs. tan$\beta$ signal grid assuming $m_a$ = 250 GeV, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.7$. Signals simulating the tW+DM + tt+DM final states are considered. Upper limits with large $\mu_{\mathrm{sig}}$ for the observed limit are capped at 500.

Model dependent upper limit on the cross section for the $m_a$ vs. $ m_{H^{\pm}}$ signal grid assuming tan$\beta$ = 1, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.35$. Only signals simulating the tW+DM final states are considered. Upper limits with large $\mu_{\mathrm{sig}}$ for the observed limit are capped at 500.

Model dependent upper limit on the cross section for the $m_{H^{\pm}}$ vs. tan$\beta$ signal grid assuming $m_a$ = 150 GeV, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.35$. Only signals simulating the tW+DM final states are considered. Upper limits with large $\mu_{\mathrm{sig}}$ for the observed limit are capped at 500.

Model dependent upper limit on the cross section for the $m_{H^{\pm}}$ vs. tan$\beta$ signal grid assuming $m_a$ = 250 GeV, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.35$. Only signals simulating the tW+DM final states are considered. Upper limits with large $\mu_{\mathrm{sig}}$ for the observed limit are capped at 500.

Model dependent upper limit on the cross section for the $m_a$ vs. $ m_{H^{\pm}}$ signal grid assuming tan$\beta$ = 1, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.35$. Signals simulating the tW+DM + tt+DM final states are considered. Upper limits with large $\mu_{\mathrm{sig}}$ for the observed limit are capped at 500.

Model dependent upper limit on the cross section for the $m_{H^{\pm}}$ vs. tan$\beta$ signal grid assuming $m_a$ = 150 GeV, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.35$. Signals simulating the tW+DM + tt+DM final states are considered. Upper limits with large $\mu_{\mathrm{sig}}$ for the observed limit are capped at 500.

Model dependent upper limit on the cross section for the $m_{H^{\pm}}$ vs. tan$\beta$ signal grid assuming $m_a$ = 250 GeV, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.35$. Signals simulating the tW+DM + tt+DM final states are considered. Upper limits with large $\mu_{\mathrm{sig}}$ for the observed limit are capped at 500.

Model dependent upper limit on the cross section for the $m_a$ vs. $ m_{H^{\pm}}$ signal grid assuming tan$\beta$ = 1, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.7$. Only signals simulating the tW+DM final states are considered. Upper limits with large $\mu_{\mathrm{sig}}$ for the observed limit are capped at 500.

Model dependent upper limit on the cross section for the $m_{H^{\pm}}$ vs. tan$\beta$ signal grid assuming $m_a$ = 150 GeV, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.7$. Only signals simulating the tW+DM final states are considered. Upper limits with large $\mu_{\mathrm{sig}}$ for the observed limit are capped at 500.

Model dependent upper limit on the cross section for the $m_{H^{\pm}}$ vs. tan$\beta$ signal grid assuming $m_a$ = 250 GeV, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.7$. Only signals simulating the tW+DM final states are considered. Upper limits with large $\mu_{\mathrm{sig}}$ for the observed limit are capped at 500.

Model dependent upper limit on the cross section for the $m_a$ vs. $ m_{H^{\pm}}$ signal grid assuming tan$\beta$ = 1, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.7$. Signals simulating the tW+DM + tt+DM final states are considered. Upper limits with large $\mu_{\mathrm{sig}}$ for the observed limit are capped at 500.

Model dependent upper limit on the cross section for the $m_{H^{\pm}}$ vs. tan$\beta$ signal grid assuming $m_a$ = 150 GeV, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.7$. Signals simulating the tW+DM + tt+DM final states are considered. Upper limits with large $\mu_{\mathrm{sig}}$ for the observed limit are capped at 500.

Model dependent upper limit on the cross section for the $m_{H^{\pm}}$ vs. tan$\beta$ signal grid assuming $m_a$ = 250 GeV, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.7$. Signals simulating the tW+DM + tt+DM final states are considered. Upper limits with large $\mu_{\mathrm{sig}}$ for the observed limit are capped at 500.

Model dependent upper limit on the cross section for the $m_a$ vs. $ m_{H^{\pm}}$ signal grid assuming tan$\beta$ = 1, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.35$. Only signals simulating the tW+DM final states are considered. Upper limits with large $\mu_{\mathrm{sig}}$ for the observed limit are capped at 500.

Model dependent upper limit on the cross section for the $m_{H^{\pm}}$ vs. tan$\beta$ signal grid assuming $m_a$ = 150 GeV, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.35$. Only signals simulating the tW+DM final states are considered. Upper limits with large $\mu_{\mathrm{sig}}$ for the observed limit are capped at 500.

Model dependent upper limit on the cross section for the $m_{H^{\pm}}$ vs. tan$\beta$ signal grid assuming $m_a$ = 250 GeV, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.35$. Only signals simulating the tW+DM final states are considered. Upper limits with large $\mu_{\mathrm{sig}}$ for the observed limit are capped at 500.

Model dependent upper limit on the cross section for the $m_a$ vs. $ m_{H^{\pm}}$ signal grid assuming tan$\beta$ = 1, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.35$. Signals simulating the tW+DM + tt+DM final states are considered. Upper limits with large $\mu_{\mathrm{sig}}$ for the observed limit are capped at 500.

Model dependent upper limit on the cross section for the $m_{H^{\pm}}$ vs. tan$\beta$ signal grid assuming $m_a$ = 150 GeV, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.35$. Signals simulating the tW+DM + tt+DM final states are considered. Upper limits with large $\mu_{\mathrm{sig}}$ for the observed limit are capped at 500.

Model dependent upper limit on the cross section for the $m_{H^{\pm}}$ vs. tan$\beta$ signal grid assuming $m_a$ = 250 GeV, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.35$. Signals simulating the tW+DM + tt+DM final states are considered. Upper limits with large $\mu_{\mathrm{sig}}$ for the observed limit are capped at 500.

Model dependent upper limit on the cross section for the $m_a$ vs. $ m_{H^{\pm}}$ signal grid assuming tan$\beta$ = 1, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.7$. Only signals simulating the tW+DM final states are considered. Upper limits with large $\mu_{\mathrm{sig}}$ for the observed limit are capped at 500.

Model dependent upper limit on the cross section for the $m_{H^{\pm}}$ vs. tan$\beta$ signal grid assuming $m_a$ = 150 GeV, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.7$. Only signals simulating the tW+DM final states are considered. Upper limits with large $\mu_{\mathrm{sig}}$ for the observed limit are capped at 500.

Model dependent upper limit on the cross section for the $m_{H^{\pm}}$ vs. tan$\beta$ signal grid assuming $m_a$ = 250 GeV, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.7$. Only signals simulating the tW+DM final states are considered. Upper limits with large $\mu_{\mathrm{sig}}$ for the observed limit are capped at 500.

Model dependent upper limit on the cross section for the $m_a$ vs. $ m_{H^{\pm}}$ signal grid assuming tan$\beta$ = 1, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.7$. Signals simulating the tW+DM + tt+DM final states are considered. Upper limits with large $\mu_{\mathrm{sig}}$ for the observed limit are capped at 500.

Model dependent upper limit on the cross section for the $m_{H^{\pm}}$ vs. tan$\beta$ signal grid assuming $m_a$ = 150 GeV, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.7$. Signals simulating the tW+DM + tt+DM final states are considered. Upper limits with large $\mu_{\mathrm{sig}}$ for the observed limit are capped at 500.

Model dependent upper limit on the cross section for the $m_{H^{\pm}}$ vs. tan$\beta$ signal grid assuming $m_a$ = 250 GeV, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.7$. Signals simulating the tW+DM + tt+DM final states are considered. Upper limits with large $\mu_{\mathrm{sig}}$ for the observed limit are capped at 500.

Model dependent upper limit on the cross section for the $m_a$ vs. $ m_{H^{\pm}}$ signal grid assuming tan$\beta$ = 1, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.35$. Only signals simulating the tW+DM final states are considered. Upper limits with large $\mu_{\mathrm{sig}}$ for the observed limit are capped at 500.

Model dependent upper limit on the cross section for the $m_{H^{\pm}}$ vs. tan$\beta$ signal grid assuming $m_a$ = 150 GeV, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.35$. Only signals simulating the tW+DM final states are considered. Upper limits with large $\mu_{\mathrm{sig}}$ for the observed limit are capped at 500.

Model dependent upper limit on the cross section for the $m_{H^{\pm}}$ vs. tan$\beta$ signal grid assuming $m_a$ = 250 GeV, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.35$. Only signals simulating the tW+DM final states are considered. Upper limits with large $\mu_{\mathrm{sig}}$ for the observed limit are capped at 500.

Model dependent upper limit on the cross section for the $m_a$ vs. $ m_{H^{\pm}}$ signal grid assuming tan$\beta$ = 1, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.35$. Signals simulating the tW+DM + tt+DM final states are considered. Upper limits with large $\mu_{\mathrm{sig}}$ for the observed limit are capped at 500.

Model dependent upper limit on the cross section for the $m_{H^{\pm}}$ vs. tan$\beta$ signal grid assuming $m_a$ = 150 GeV, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.35$. Signals simulating the tW+DM + tt+DM final states are considered. Upper limits with large $\mu_{\mathrm{sig}}$ for the observed limit are capped at 500.

Model dependent upper limit on the cross section for the $m_{H^{\pm}}$ vs. tan$\beta$ signal grid assuming $m_a$ = 250 GeV, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.35$. Signals simulating the tW+DM + tt+DM final states are considered. Upper limits with large $\mu_{\mathrm{sig}}$ for the observed limit are capped at 500.

The distributions of $m_{\mathrm{b1},\mathrm{W-tagged}}$ in the 0L inclusive signal region. For each bin yields for the data and total SM prediction are provided. The SM prediction is provided with the total uncertainty, including the MC statistical uncertainty, detector-related systematic uncertainties and theoretical uncertainties. The rightmost bin includes overflow events.

The distributions of $m_{\mathrm{T}}^{\mathrm{b,E_{\mathrm{T}^{\mathrm{miss}}}}}$ in the 0L inclusive signal region. For each bin yields for the data and total SM prediction are provided. The SM prediction is provided with the total uncertainty, including the MC statistical uncertainty, detector-related systematic uncertainties and theoretical uncertainties. The rightmost bin includes overflow events.

The distributions of $N_{\mathrm{W-tagged}}$ in the 0L inclusive signal region. For each bin yields for the data and total SM prediction are provided. The SM prediction is provided with the total uncertainty, including the MC statistical uncertainty, detector-related systematic uncertainties and theoretical uncertainties. The rightmost bin includes overflow events.

The distributions of $m_{\mathrm{b1},\mathrm{\cancel{b1}}}$ in the hadronic top inclusive signal region. For each bin yields for the data and total SM prediction are provided. The SM prediction is provided with the total uncertainty, including the MC statistical uncertainty, detector-related systematic uncertainties and theoretical uncertainties. The rightmost bin includes overflow events.

The distributions of $m_{\mathrm{b1},\mathrm{\cancel{b1}}}$ in the leptonic top inclusive signal region. For each bin yields for the data and total SM prediction are provided. The SM prediction is provided with the total uncertainty, including the MC statistical uncertainty, detector-related systematic uncertainties and theoretical uncertainties. The rightmost bin includes overflow events.

The distributions of $m_{\mathrm{b1},\mathrm{\cancel{b1}}}$ in the leptonic top inclusive signal region. For each bin yields for the data and total SM prediction are provided. The SM prediction is provided with the total uncertainty, including the MC statistical uncertainty, detector-related systematic uncertainties and theoretical uncertainties. The rightmost bin includes overflow events.

Cutflow for the reference point $(\it{m}_{\mathrm{H^{\pm}}}, \it{m}_{a}, tan\beta, sin\theta )=$ (500,100,1,0.7) , (800,150,20,0.7), (600,250,30,0.7), (1000,400,1,0.7) in 0L regions. Results are shown including all correction factors applied to simulation, and is normalised to 139 fb$^{-1}$.

Cutflow for the reference point $(\it{m}_{\mathrm{H^{\pm}}}, \it{m}_{a}, tan\beta, sin\theta )=$ (500,100,1,0.7) , (800,150,20,0.7), (600,250,30,0.7), (1000,400,1,0.7) in 1L leptonic top regions. Results are shown including all correction factors applied to simulation, and is normalised to 139 fb$^{-1}$.

Cutflow for the reference point $(\it{m}_{\mathrm{H^{\pm}}}, \it{m}_{a}, tan\beta, sin\theta )=$ (500,100,1,0.7) , (800,150,20,0.7), (600,250,30,0.7), (1000,400,1,0.7) in 1L hadronic top regions. Results are shown including all correction factors applied to simulation, and is normalised to 139 fb$^{-1}$.

Signal acceptance in the 0L region for 2HDM+a model DM signals on the plane defined by m$_a$--m$_{H^{\pm}}$ assuming tan$\beta$ = 1, m$_{\chi}$= 10 GeV and sin$\theta$ = 0.7. Please mind that the acceptance given in the table is multiplied by factor of $10^{3}$

Signal efficiency in the 0L region for 2HDM+a model DM signals on the plane defined by m$_a$--m$_{H^{\pm}}$ assuming tan$\beta$ = 1, m$_{\chi}$= 10 GeV and sin$\theta$ = 0.7. Please mind that the efficiency given in the table is multiplied by factor of $10^{2}$

Signal acceptance in the 0L region for 2HDM+a model DM signals on the plane defined by m$_{H^{\pm}}$--tan$\beta$ assuming m$_a$ = 150 GeV, m$_{\chi}$= 10 GeV and sin$\theta$ = 0.7. Please mind that the acceptance given in the table is multiplied by factor of $10^{3}$

Signal efficiency in the 0L region for 2HDM+a model DM signals on the plane defined by m$_{H^{\pm}}$--tan$\beta$ assuming m$_a$ = 150 GeV, m$_{\chi}$= 10 GeV and sin$\theta$ = 0.7. Please mind that the efficiency given in the table is multiplied by factor of $10^{2}$

Signal acceptance in the 0L region for 2HDM+a model DM signals on the plane defined by m$_{H^{\pm}}$--tan$\beta$ assuming m$_a$ = 250 GeV, m$_{\chi}$= 10 GeV and sin$\theta$ = 0.7. Please mind that the acceptance given in the table is multiplied by factor of $10^{3}$

Signal efficiency in the 0L region for 2HDM+a model DM signals on the plane defined by m$_{H^{\pm}}$--tan$\beta$ assuming m$_a$ = 250 GeV, m$_{\chi}$= 10 GeV and sin$\theta$ = 0.7. Please mind that the efficiency given in the table is multiplied by factor of $10^{2}$

Signal acceptance in the 1L region for 2HDM+a model DM signals on the plane defined by m$_a$--m$_{H^{\pm}}$ assuming tan$\beta$ = 1, m$_{\chi}$= 10 GeV and sin$\theta$ = 0.7. Please mind that the acceptance given in the table is multiplied by factor of $10^{3}$

Signal efficiency in the 1L region for 2HDM+a model DM signals on the plane defined by m$_a$--m$_{H^{\pm}}$ assuming tan$\beta$ = 1, m$_{\chi}$= 10 GeV and sin$\theta$ = 0.7. Please mind that the efficiency given in the table is multiplied by factor of $10^{2}$

Signal acceptance in the 1L region for 2HDM+a model DM signals on the plane defined by m$_{H^{\pm}}$--tan$\beta$ assuming m$_a$ = 150 GeV, m$_{\chi}$= 10 GeV and sin$\theta$ = 0.7. Please mind that the acceptance given in the table is multiplied by factor of $10^{3}$

Signal efficiency in the 1L region for 2HDM+a model DM signals on the plane defined by m$_{H^{\pm}}$--tan$\beta$ assuming m$_a$ = 150 GeV, m$_{\chi}$= 10 GeV and sin$\theta$ = 0.7. Please mind that the efficiency given in the table is multiplied by factor of $10^{2}$

Signal acceptance in the 1L region for 2HDM+a model DM signals on the plane defined by m$_{H^{\pm}}$--tan$\beta$ assuming m$_a$ = 250 GeV, m$_{\chi}$= 10 GeV and sin$\theta$ = 0.7. Please mind that the acceptance given in the table is multiplied by factor of $10^{3}$

Signal efficiency in the 1L region for 2HDM+a model DM signals on the plane defined by m$_{H^{\pm}}$--tan$\beta$ assuming m$_a$ = 250 GeV, m$_{\chi}$= 10 GeV and sin$\theta$ = 0.7. Please mind that the efficiency given in the table is multiplied by factor of $10^{2}$

Signal acceptance in the 0L region for 2HDM+a model DM signals on the plane defined by m$_a$--m$_{H^{\pm}}$ assuming tan$\beta$ = 1, m$_{\chi}$= 10 GeV and sin$\theta$ = 0.35. Please mind that the acceptance given in the table is multiplied by factor of $10^{3}$

Signal efficiency in the 0L region for 2HDM+a model DM signals on the plane defined by m$_a$--m$_{H^{\pm}}$ assuming tan$\beta$ = 1, m$_{\chi}$= 10 GeV and sin$\theta$ = 0.35. Please mind that the efficiency given in the table is multiplied by factor of $10^{2}$

Signal acceptance in the 0L region for 2HDM+a model DM signals on the plane defined by m$_{H^{\pm}}$--tan$\beta$ assuming m$_a$ = 150 GeV, m$_{\chi}$= 10 GeV and sin$\theta$ = 0.35. Please mind that the acceptance given in the table is multiplied by factor of $10^{3}$

Signal efficiency in the 0L region for 2HDM+a model DM signals on the plane defined by m$_{H^{\pm}}$--tan$\beta$ assuming m$_a$ = 150 GeV, m$_{\chi}$= 10 GeV and sin$\theta$ = 0.35. Please mind that the efficiency given in the table is multiplied by factor of $10^{2}$

Signal acceptance in the 0L region for 2HDM+a model DM signals on the plane defined by m$_{H^{\pm}}$--tan$\beta$ assuming m$_a$ = 250 GeV, m$_{\chi}$= 10 GeV and sin$\theta$ = 0.35. Please mind that the acceptance given in the table is multiplied by factor of $10^{3}$

Signal efficiency in the 0L region for 2HDM+a model DM signals on the plane defined by m$_{H^{\pm}}$--tan$\beta$ assuming m$_a$ = 250 GeV, m$_{\chi}$= 10 GeV and sin$\theta$ = 0.35. Please mind that the efficiency given in the table is multiplied by factor of $10^{2}$

Signal acceptance in the 1L region for 2HDM+a model DM signals on the plane defined by m$_a$--m$_{H^{\pm}}$ assuming tan$\beta$ = 1, m$_{\chi}$= 10 GeV and sin$\theta$ = 0.35. Please mind that the acceptance given in the table is multiplied by factor of $10^{3}$

Signal efficiency in the 1L region for 2HDM+a model DM signals on the plane defined by m$_a$--m$_{H^{\pm}}$ assuming tan$\beta$ = 1, m$_{\chi}$= 10 GeV and sin$\theta$ = 0.35. Please mind that the efficiency given in the table is multiplied by factor of $10^{2}$

Signal acceptance in the 1L region for 2HDM+a model DM signals on the plane defined by m$_{H^{\pm}}$--tan$\beta$ assuming m$_a$ = 150 GeV, m$_{\chi}$= 10 GeV and sin$\theta$ = 0.35. Please mind that the acceptance given in the table is multiplied by factor of $10^{3}$

Signal efficiency in the 1L region for 2HDM+a model DM signals on the plane defined by m$_{H^{\pm}}$--tan$\beta$ assuming m$_a$ = 150 GeV, m$_{\chi}$= 10 GeV and sin$\theta$ = 0.35. Please mind that the efficiency given in the table is multiplied by factor of $10^{2}$

Signal acceptance in the 1L region for 2HDM+a model DM signals on the plane defined by m$_{H^{\pm}}$--tan$\beta$ assuming m$_a$ = 250 GeV, m$_{\chi}$= 10 GeV and sin$\theta$ = 0.35. Please mind that the acceptance given in the table is multiplied by factor of $10^{3}$

Signal efficiency in the 1L region for 2HDM+a model DM signals on the plane defined by m$_{H^{\pm}}$--tan$\beta$ assuming m$_a$ = 250 GeV, m$_{\chi}$= 10 GeV and sin$\theta$ = 0.35. Please mind that the efficiency given in the table is multiplied by factor of $10^{2}$

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.

Three-body selection. Distributions of $M_{\Delta}^R$ in (a,b) $SR_{W}^{3-body}$ and (c,d) $SR_{T}^{3-body}$ for (left) same-flavour and (right) different-flavour events satisfying the selection criteria of the given SR, except the one for the presented variable, after the background fit. The contributions from all SM backgrounds are shown as a histogram stack. ''Others'' includes contributions from $VVV$, $t\bar{t} t$, $t\bar{t}t\bar{t}$, $t\bar{t} W$, $t\bar{t} WW$, $t\bar{t} WZ$, $t\bar{t} H$, and $tZ$ processes. The hatched bands represent the total statistical and systematic uncertainty. The rightmost bin of each plot includes overflow events. Reference top squark pair production signal models are overlayed for comparison. Red arrows in the upper panels indicate the signal region selection criteria. The bottom panels show the ratio of the observed data to the total SM background prediction, with hatched bands representing the total uncertainty in the background prediction; red arrows show data outside the vertical-axis range.

Three-body selection. Distributions of $M_{\Delta}^R$ in (a,b) $SR_{W}^{3-body}$ and (c,d) $SR_{T}^{3-body}$ for (left) same-flavour and (right) different-flavour events satisfying the selection criteria of the given SR, except the one for the presented variable, after the background fit. The contributions from all SM backgrounds are shown as a histogram stack. ''Others'' includes contributions from $VVV$, $t\bar{t} t$, $t\bar{t}t\bar{t}$, $t\bar{t} W$, $t\bar{t} WW$, $t\bar{t} WZ$, $t\bar{t} H$, and $tZ$ processes. The hatched bands represent the total statistical and systematic uncertainty. The rightmost bin of each plot includes overflow events. Reference top squark pair production signal models are overlayed for comparison. Red arrows in the upper panels indicate the signal region selection criteria. The bottom panels show the ratio of the observed data to the total SM background prediction, with hatched bands representing the total uncertainty in the background prediction; red arrows show data outside the vertical-axis range.

Three-body selection. Distributions of $M_{\Delta}^R$ in (a,b) $SR_{W}^{3-body}$ and (c,d) $SR_{T}^{3-body}$ for (left) same-flavour and (right) different-flavour events satisfying the selection criteria of the given SR, except the one for the presented variable, after the background fit. The contributions from all SM backgrounds are shown as a histogram stack. ''Others'' includes contributions from $VVV$, $t\bar{t} t$, $t\bar{t}t\bar{t}$, $t\bar{t} W$, $t\bar{t} WW$, $t\bar{t} WZ$, $t\bar{t} H$, and $tZ$ processes. The hatched bands represent the total statistical and systematic uncertainty. The rightmost bin of each plot includes overflow events. Reference top squark pair production signal models are overlayed for comparison. Red arrows in the upper panels indicate the signal region selection criteria. The bottom panels show the ratio of the observed data to the total SM background prediction, with hatched bands representing the total uncertainty in the background prediction; red arrows show data outside the vertical-axis range.

Four-body selection. (a) distributions of $E_{T}^{miss}$ in $SR^{4-body}_{Small\,\Delta m}$ and (b) distribution of $R_{2\ell 4j}$ in $SR^{4-body}_{Large\,\Delta m}$ for events satisfying the selection criteria of the given SR, except the one for the presented variable, after the background fit. The contributions from all SM backgrounds are shown as a histogram stack. ''Others'' includes contributions from $VVV$, $t\bar{t} t$, $t\bar{t}$, $t\bar{t} W$, $t\bar{t} WW$, $t\bar{t} WZ$, $t\bar{t} H$, and $tZ$ processes. The hatched bands represent the total statistical and systematic uncertainty. The rightmost bin of each plot includes overflow events. Reference top squark pair production signal models are overlayed for comparison. Red arrows in the upper panel indicate the signal region selection criteria. The bottom panels show the ratio of the observed data to the total SM background prediction, with hatched bands representing the total uncertainty in the background prediction; red arrows show data outside the vertical-axis range.

Four-body selection. (a) distributions of $E_{T}^{miss}$ in $SR^{4-body}_{Small\,\Delta m}$ and (b) distribution of $R_{2\ell 4j}$ in $SR^{4-body}_{Large\,\Delta m}$ for events satisfying the selection criteria of the given SR, except the one for the presented variable, after the background fit. The contributions from all SM backgrounds are shown as a histogram stack. ''Others'' includes contributions from $VVV$, $t\bar{t} t$, $t\bar{t}$, $t\bar{t} W$, $t\bar{t} WW$, $t\bar{t} WZ$, $t\bar{t} H$, and $tZ$ processes. The hatched bands represent the total statistical and systematic uncertainty. The rightmost bin of each plot includes overflow events. Reference top squark pair production signal models are overlayed for comparison. Red arrows in the upper panel indicate the signal region selection criteria. The bottom panels show the ratio of the observed data to the total SM background prediction, with hatched bands representing the total uncertainty in the background prediction; red arrows show data outside the vertical-axis range.

Exclusion limit contour (95% CL) for a simplified model assuming $\tilde{t}_1$ pair production, decaying via $\tilde{t}_1 \rightarrow t^{(*)}\tilde{\chi}_1^0$ with 100% branching ratio, in the (a) $m(\tilde{t}_1)$--$m(\tilde{\chi}_1^0)$ and (b) $m(\tilde{t}_1)$--$\Delta m(\tilde{t}_1,\tilde{\chi}_1^0)$ planes. The dashed lines and the shaded bands are the expected limits and their $\pm1\sigma$ uncertainties. The thick solid lines are the observed limits for the central value of the signal cross-section. The expected and observed limits do not include the effect of the theoretical uncertainties in the signal cross-section. The dotted lines show the effect on the observed limit when varying the signal cross-section by $\pm1\sigma$ of the theoretical uncertainty.

Exclusion limit contour (95% CL) for a simplified model assuming $\tilde{t}_1$ pair production, decaying via $\tilde{t}_1 \rightarrow t^{(*)}\tilde{\chi}_1^0$ with 100% branching ratio, in the (a) $m(\tilde{t}_1)$--$m(\tilde{\chi}_1^0)$ and (b) $m(\tilde{t}_1)$--$\Delta m(\tilde{t}_1,\tilde{\chi}_1^0)$ planes. The dashed lines and the shaded bands are the expected limits and their $\pm1\sigma$ uncertainties. The thick solid lines are the observed limits for the central value of the signal cross-section. The expected and observed limits do not include the effect of the theoretical uncertainties in the signal cross-section. The dotted lines show the effect on the observed limit when varying the signal cross-section by $\pm1\sigma$ of the theoretical uncertainty.

Exclusion limit contour (95% CL) for a simplified model assuming $\tilde{t}_1$ pair production, decaying via $\tilde{t}_1 \rightarrow t^{(*)}\tilde{\chi}_1^0$ with 100% branching ratio, in the (a) $m(\tilde{t}_1)$--$m(\tilde{\chi}_1^0)$ and (b) $m(\tilde{t}_1)$--$\Delta m(\tilde{t}_1,\tilde{\chi}_1^0)$ planes. The dashed lines and the shaded bands are the expected limits and their $\pm1\sigma$ uncertainties. The thick solid lines are the observed limits for the central value of the signal cross-section. The expected and observed limits do not include the effect of the theoretical uncertainties in the signal cross-section. The dotted lines show the effect on the observed limit when varying the signal cross-section by $\pm1\sigma$ of the theoretical uncertainty.

Exclusion limit contour (95% CL) for a simplified model assuming $\tilde{t}_1$ pair production, decaying via $\tilde{t}_1 \rightarrow t^{(*)}\tilde{\chi}_1^0$ with 100% branching ratio, in the (a) $m(\tilde{t}_1)$--$m(\tilde{\chi}_1^0)$ and (b) $m(\tilde{t}_1)$--$\Delta m(\tilde{t}_1,\tilde{\chi}_1^0)$ planes. The dashed lines and the shaded bands are the Observed limits and their $\pm1\sigma$ uncertainties. The thick solid lines are the observed limits for the central value of the signal cross-section. The expected and observed limits do not include the effect of the theoretical uncertainties in the signal cross-section. The dotted lines show the effect on the observed limit when varying the signal cross-section by $\pm1\sigma$ of the theoretical uncertainty.

Exclusion limit contour (95% CL) for a simplified model assuming $\tilde{t}_1$ pair production, decaying via $\tilde{t}_1 \rightarrow t^{(*)}\tilde{\chi}_1^0$ with 100% branching ratio, in the (a) $m(\tilde{t}_1)$--$m(\tilde{\chi}_1^0)$ and (b) $m(\tilde{t}_1)$--$\Delta m(\tilde{t}_1,\tilde{\chi}_1^0)$ planes. The dashed lines and the shaded bands are the expected limits and their $\pm1\sigma$ uncertainties. The thick solid lines are the observed limits for the central value of the signal cross-section. The expected and observed limits do not include the effect of the theoretical uncertainties in the signal cross-section. The dotted lines show the effect on the observed limit when varying the signal cross-section by $\pm1\sigma$ of the theoretical uncertainty.

Exclusion limit contour (95% CL) for a simplified model assuming $\tilde{t}_1$ pair production, decaying via $\tilde{t}_1 \rightarrow t^{(*)}\tilde{\chi}_1^0$ with 100% branching ratio, in the (a) $m(\tilde{t}_1)$--$m(\tilde{\chi}_1^0)$ and (b) $m(\tilde{t}_1)$--$\Delta m(\tilde{t}_1,\tilde{\chi}_1^0)$ planes. The dashed lines and the shaded bands are the expected limits and their $\pm1\sigma$ uncertainties. The thick solid lines are the observed limits for the central value of the signal cross-section. The expected and observed limits do not include the effect of the theoretical uncertainties in the signal cross-section. The dotted lines show the effect on the observed limit when varying the signal cross-section by $\pm1\sigma$ of the theoretical uncertainty.

Exclusion limit contour (95% CL) for a simplified model assuming $\tilde{t}_1$ pair production, decaying via $\tilde{t}_1 \rightarrow t^{(*)}\tilde{\chi}_1^0$ with 100\% branching ratio, in the (a) $m(\tilde{t}_1)$--$m(\tilde{\chi}_1^0)$ and (b) $m(\tilde{t}_1)$--$\Delta m(\tilde{t}_1,\tilde{\chi}_1^0)$ planes. The dashed lines and the shaded bands are the expected limits and their $\pm1\sigma$ uncertainties. The thick solid lines are the observed limits for the central value of the signal cross-section. The expected and observed limits do not include the effect of the theoretical uncertainties in the signal cross-section. The dotted lines show the effect on the observed limit when varying the signal cross-section by $\pm1\sigma$ of the theoretical uncertainty.

Exclusion limit contour (95% CL) for a simplified model assuming $\tilde{t}_1$ pair production, decaying via $\tilde{t}_1 \rightarrow t^{(*)}\tilde{\chi}_1^0$ with 100\% branching ratio, in the (a) $m(\tilde{t}_1)$--$m(\tilde{\chi}_1^0)$ and (b) $m(\tilde{t}_1)$--$\Delta m(\tilde{t}_1,\tilde{\chi}_1^0)$ planes. The dashed lines and the shaded bands are the expected limits and their $\pm1\sigma$ uncertainties. The thick solid lines are the observed limits for the central value of the signal cross-section. The expected and observed limits do not include the effect of the theoretical uncertainties in the signal cross-section. The dotted lines show the effect on the observed limit when varying the signal cross-section by $\pm1\sigma$ of the theoretical uncertainty.

Exclusion limit contour (95% CL) for a simplified model assuming $\tilde{t}_1$ pair production, decaying via $\tilde{t}_1 \rightarrow t^{(*)}\tilde{\chi}_1^0$ with 100\% branching ratio, in the (a) $m(\tilde{t}_1)$--$m(\tilde{\chi}_1^0)$ and (b) $m(\tilde{t}_1)$--$\Delta m(\tilde{t}_1,\tilde{\chi}_1^0)$ planes. The dashed lines and the shaded bands are the expected limits and their $\pm1\sigma$ uncertainties. The thick solid lines are the observed limits for the central value of the signal cross-section. The expected and observed limits do not include the effect of the theoretical uncertainties in the signal cross-section. The dotted lines show the effect on the observed limit when varying the signal cross-section by $\pm1\sigma$ of the theoretical uncertainty.

Exclusion limit contour (95% CL) for a simplified model assuming $\tilde{t}_1$ pair production, decaying via $\tilde{t}_1 \rightarrow t^{(*)}\tilde{\chi}_1^0$ with 100\% branching ratio, in the (a) $m(\tilde{t}_1)$--$m(\tilde{\chi}_1^0)$ and (b) $m(\tilde{t}_1)$--$\Delta m(\tilde{t}_1,\tilde{\chi}_1^0)$ planes. The dashed lines and the shaded bands are the expected limits and their $\pm1\sigma$ uncertainties. The thick solid lines are the observed limits for the central value of the signal cross-section. The expected and observed limits do not include the effect of the theoretical uncertainties in the signal cross-section. The dotted lines show the effect on the observed limit when varying the signal cross-section by $\pm1\sigma$ of the theoretical uncertainty.

Exclusion limit contour (95% CL) for a simplified model assuming $\tilde{t}_1$ pair production, decaying via $\tilde{t}_1 \rightarrow t^{(*)}\tilde{\chi}_1^0$ with 100\% branching ratio, in the (a) $m(\tilde{t}_1)$--$m(\tilde{\chi}_1^0)$ and (b) $m(\tilde{t}_1)$--$\Delta m(\tilde{t}_1,\tilde{\chi}_1^0)$ planes. The dashed lines and the shaded bands are the expected limits and their $\pm1\sigma$ uncertainties. The thick solid lines are the observed limits for the central value of the signal cross-section. The expected and observed limits do not include the effect of the theoretical uncertainties in the signal cross-section. The dotted lines show the effect on the observed limit when varying the signal cross-section by $\pm1\sigma$ of the theoretical uncertainty.

Exclusion limit contour (95% CL) for a simplified model assuming $\tilde{t}_1$ pair production, decaying via $\tilde{t}_1 \rightarrow t^{(*)}\tilde{\chi}_1^0$ with 100\% branching ratio, in the (a) $m(\tilde{t}_1)$--$m(\tilde{\chi}_1^0)$ and (b) $m(\tilde{t}_1)$--$\Delta m(\tilde{t}_1,\tilde{\chi}_1^0)$ planes. The dashed lines and the shaded bands are the expected limits and their $\pm1\sigma$ uncertainties. The thick solid lines are the observed limits for the central value of the signal cross-section. The expected and observed limits do not include the effect of the theoretical uncertainties in the signal cross-section. The dotted lines show the effect on the observed limit when varying the signal cross-section by $\pm1\sigma$ of the theoretical uncertainty.

Exclusion limits for (a) $t\bar{t} + \phi $ scalar and (b) $t\bar{t} + a $ pseudoscalar models as a function of the mediator mass for a DM particle mass of $m(\chi)=1$ GeV. 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 shows the observed (expected) exclusion limits.

Exclusion limits for (a) $t\bar{t} + \phi $ scalar and (b) $t\bar{t} + a $ pseudoscalar models as a function of the mediator mass for a DM particle mass of $m(\chi)=1$ GeV. 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 shows the observed (expected) exclusion limits.

Exclusion limits for (a) $t\bar{t} + \phi $ scalar and (b) $t\bar{t} + a $ pseudoscalar models as a function of the mediator mass for a DM particle mass of $m(\chi)=1$ GeV. 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 shows the observed (expected) exclusion limits.

Exclusion limits for (a) $t\bar{t} + \phi $ scalar and (b) $t\bar{t} + a $ pseudoscalar models as a function of the mediator mass for a DM particle mass of $m(\chi)=1$ GeV. 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 shows the observed (expected) exclusion limits.

Two-body selection. Background fit results for $\mathrm{CR}^{\mathrm{2-body}}_{t\bar{t}}$, $\mathrm{CR}^{\mathrm{2-body}}_{t\bar{t}Z}$, $\mathrm{VR}^{\mathrm{2-body}}_{t\bar{t}, DF}$, $\mathrm{VR}^{\mathrm{2-body}}_{t\bar{t}, SF}$ and $\mathrm{VR}^{\mathrm{2-body}}_{t\bar{t} Z}$. ''Others'' includes contributions from $VVV$, $t\bar{t} t$, $t\bar{t}t\bar{t}$, $t\bar{t} W$, $t\bar{t} WW$, $t\bar{t} WZ$, $t\bar{t} H$, and $tZ$ processes. Combined statistical and systematic uncertainties are given. Entries marked `--' indicate a negligible background contribution (less than 0.001 events). The individual uncertainties can be correlated, and do not necessarily add up in quadrature to the total background uncertainty.

Three-body selection. Background fit results for $\mathrm{CR}^{\mathrm{3-body}}_{t\bar{t}}$, $\mathrm{CR}^{\mathrm{3-body}}_{VV}$, $\mathrm{CR}^{\mathrm{2-body}}_{t\bar{t}Z}$, $\mathrm{VR}^{\mathrm{3-body}}_{VV}$, $\mathrm{VR(1)}^{\mathrm{3-body}}_{t\bar{t}}$ and $\mathrm{VR(2)}^{\mathrm{3-body}}_{t\bar{t}}$. ''Others'' includes contributions from $VVV$, $t\bar{t} t$, $t\bar{t}t\bar{t}$, $t\bar{t} W$, $t\bar{t} WW$, $t\bar{t} WZ$, $t\bar{t} H$, and $tZ$ processes. Combined statistical and systematic uncertainties are given. Entries marked `--' indicate a negligible background contribution (less than 0.001 events). The individual uncertainties can be correlated, and do not necessarily add up in quadrature to the total background uncertainty.

Four-body selection. Background fit results for $\mathrm{CR}^{\mathrm{4-body}}_{t\bar{t}}$,$\mathrm{CR}^{\mathrm{4-body}}_{VV}$, $\mathrm{VR}^{\mathrm{4-body}}_{t\bar{t}}$, $VR^{4-body}_{VV}$ and $\mathrm{VR}^{\mathrm{4-body}}_{VV,lll}$. The ''Others'' category contains the contributions from $VVV$, $t\bar{t} t$, $t\bar{t}t\bar{t}$, $t\bar{t} W$, $t\bar{t} WW$, $t\bar{t} WZ$, $t\bar{t} H$, and $tZ$. Combined statistical and systematic uncertainties are given. Entries marked `--' indicate a negligible background contribution (less than 0.001 events). The individual uncertainties can be correlated, and do not necessarily add up in quadrature to the total background uncertainty.

Two-body selection. Background fit results for the different-flavour leptons binned SRs. The ''Others'' category contains the contributions from $VVV$, $t\bar{t} t$, $t\bar{t}t\bar{t}$, $t\bar{t} W$, $t\bar{t} WW$, $t\bar{t} WZ$, $t\bar{t} H$, and $tZ$. Combined statistical and systematic uncertainties are given. Entries marked `--' indicate a negligible background contribution (less than 0.001 events). The individual uncertainties can be correlated, and do not necessarily add up in quadrature to the total background uncertainty.

Two-body selection. Background fit results for the same-flavour leptons binned SRs. The ''Others'' category contains the contributions from $VVV$, $t\bar{t} t$, $t\bar{t}t\bar{t}$, $t\bar{t} W$, $t\bar{t} WW$, $t\bar{t} WZ$, $t\bar{t} H$, and $tZ$. Combined statistical and systematic uncertainties are given. The individual uncertainties can be correlated, and do not necessarily add up in quadrature to the total background uncertainty.

Three-body selection. Observed event yields and background fit results for the three-body selection SRs. The ''Others'' category contains contributions from $VVV$, $t\bar{t} t$, $t\bar{t}t\bar{t}$, $t\bar{t} W$, $t\bar{t} WW$, $t\bar{t} WZ$, $t\bar{t} H$, and $tZ$. Combined statistical and systematic uncertainties are given. Entries marked `--' indicate a negligible background contribution (less than 0.001 events). The individual uncertainties can be correlated, and do not necessarily add up in quadrature to the total background uncertainty.

Four-body selection. Observed event yields and background fit results for SR$^{\mathrm{4-body}}_{\mathrm{Small}\,\Delta m}$ and SR$^{\mathrm{4-body}}_{\mathrm{Large}\,\Delta m}$. The ''Others'' category contains the contributions from $VVV$, $t\bar{t} t$, $t\bar{t}t\bar{t}$, $t\bar{t} W$, $t\bar{t} WW$, $t\bar{t} WZ$, $t\bar{t} H$, and $tZ$. Combined statistical and systematic uncertainties are given. The individual uncertainties can be correlated, and do not necessarily add up in quadrature to the total background uncertainty.

Exclusion limits contours (95% CL) for a simplified model assuming $\tilde{t}_1$ pair production, decaying via $\tilde{t}_1 \rightarrow t^{(*)}\tilde{\chi}^0_1$ with 100% branching ratio in $\tilde{t}_1--\tilde{\chi}^0_1$ masses planes. The dashed lines and the shaded bands are the expected limit and its $\pm 1\sigma$ uncertainty. The thick solid lines are the observed limits for the central value of the signal cross-section. The expected and observed limits do not include the effect of the theoretical uncertainties in the signal cross-section. The dotted lines show the effect on the observed limit when varying the signal cross-section by $\pm1\sigma$ of the theoretical uncertainty. The exclusion limits contours for the two-body, three-body and four-body selections are respectively shown in blue, green and red.

Exclusion limits contours (95% CL) for a simplified model assuming $\tilde{t}_1$ pair production, decaying via $\tilde{t}_1 \rightarrow t^{(*)}\tilde{\chi}^0_1$ with 100% branching ratio in $\tilde{t}_1--\tilde{\chi}^0_1$ masses planes. The dashed lines and the shaded bands are the expected limit and its $\pm 1\sigma$ uncertainty. The thick solid lines are the observed limits for the central value of the signal cross-section. The expected and observed limits do not include the effect of the theoretical uncertainties in the signal cross-section. The dotted lines show the effect on the observed limit when varying the signal cross-section by $\pm1\sigma$ of the theoretical uncertainty. The exclusion limits contours for the two-body, three-body and four-body selections are respectively shown in blue, green and red.

Exclusion limit contour (95% CL) for a simplified model assuming $\tilde{t}_1$ pair production, decaying via $\tilde{t}_1 \rightarrow t \tilde{\chi}_1^0$ with 100% branching ratio, in $\tilde{t}_1$--$\tilde{\chi}_1^0$ masses plane. The dashed lines and the shaded bands are the expected limit and its $\pm1\sigma$ uncertainty. The thick solid lines are the observed limits for the central value of the signal cross-section. The expected and observed limits do not include the effect of the theoretical uncertainties in the signal cross-section. The dotted lines show the effect on the observed limit when varying the signal cross-section by $\pm1\sigma$ of the theoretical uncertainty. The observed (a) and expected (b) CLs values are respectively shown.

Exclusion limit contour (95% CL) for a simplified model assuming $\tilde{t}_1$ pair production, decaying via $\tilde{t}_1 \rightarrow t \tilde{\chi}_1^0$ with 100% branching ratio, in $\tilde{t}_1$--$\tilde{\chi}_1^0$ masses plane. The dashed lines and the shaded bands are the expected limit and its $\pm1\sigma$ uncertainty. The thick solid lines are the observed limits for the central value of the signal cross-section. The expected and observed limits do not include the effect of the theoretical uncertainties in the signal cross-section. The dotted lines show the effect on the observed limit when varying the signal cross-section by $\pm1\sigma$ of the theoretical uncertainty. The observed (a) and expected (b) CLs values are respectively shown.

Exclusion limit contour (95% CL) for a simplified model assuming $\tilde{t}_1$ pair production, decaying via $\tilde{t}_1 \rightarrow t \tilde{\chi}_1^0$ with 100% branching ratio, in $\tilde{t}_1$--$\tilde{\chi}_1^0$ masses plane. The dashed lines and the shaded bands are the expected limit and its $\pm1\sigma$ uncertainty. The thick solid lines are the observed limits for the central value of the signal cross-section. The expected and observed limits do not include the effect of the theoretical uncertainties in the signal cross-section. The dotted lines show the effect on the observed limit when varying the signal cross-section by $\pm1\sigma$ of the theoretical uncertainty. The observed (a) and expected (b) CLs values are respectively shown.

Exclusion limit contour (95% CL) for a simplified model assuming $\tilde{t}_1$ pair production, decaying via $\tilde{t}_1 \rightarrow t \tilde{\chi}_1^0$ with 100% branching ratio, in $\tilde{t}_1$--$\tilde{\chi}_1^0$ masses plane. The dashed lines and the shaded bands are the expected limit and its $\pm1\sigma$ uncertainty. The thick solid lines are the observed limits for the central value of the signal cross-section. The expected and observed limits do not include the effect of the theoretical uncertainties in the signal cross-section. The dotted lines show the effect on the observed limit when varying the signal cross-section by $\pm1\sigma$ of the theoretical uncertainty. The observed (a) and expected (b) CLs values are respectively shown.

Exclusion limit contour (95% CL) for a simplified model assuming $\tilde{t}_1$ pair production, decaying via $\tilde{t}_1 \rightarrow b W \tilde{\chi}_1^0$ with 100% branching ratio, in $\tilde{t}_1$--$\tilde{\chi}_1^0$ masses plane. The dashed lines and the shaded bands are the expected limit and its $\pm1\sigma$ uncertainty. The thick solid lines are the observed limits for the central value of the signal cross-section. The expected and observed limits do not include the effect of the theoretical uncertainties in the signal cross-section. The dotted lines show the effect on the observed limit when varying the signal cross-section by $\pm1\sigma$ of the theoretical uncertainty. The observed (a) and expected (b) CLs values are respectively shown.

Exclusion limit contour (95% CL) for a simplified model assuming $\tilde{t}_1$ pair production, decaying via $\tilde{t}_1 \rightarrow b W \tilde{\chi}_1^0$ with 100% branching ratio, in $\tilde{t}_1$--$\tilde{\chi}_1^0$ masses plane. The dashed lines and the shaded bands are the expected limit and its $\pm1\sigma$ uncertainty. The thick solid lines are the observed limits for the central value of the signal cross-section. The expected and observed limits do not include the effect of the theoretical uncertainties in the signal cross-section. The dotted lines show the effect on the observed limit when varying the signal cross-section by $\pm1\sigma$ of the theoretical uncertainty. The observed (a) and expected (b) CLs values are respectively shown.

Exclusion limit contour (95% CL) for a simplified model assuming $\tilde{t}_1$ pair production, decaying via $\tilde{t}_1 \rightarrow b W \tilde{\chi}_1^0$ with 100% branching ratio, in $\tilde{t}_1$--$\tilde{\chi}_1^0$ masses plane. The dashed lines and the shaded bands are the expected limit and its $\pm1\sigma$ uncertainty. The thick solid lines are the observed limits for the central value of the signal cross-section. The expected and observed limits do not include the effect of the theoretical uncertainties in the signal cross-section. The dotted lines show the effect on the observed limit when varying the signal cross-section by $\pm 1\sigma$ of the theoretical uncertainty. The observed (a) and expected (b) CLs values are respectively shown.

Exclusion limit contour (95% CL) for a simplified model assuming $\tilde{t}_1$ pair production, decaying via $\tilde{t}_1 \rightarrow b W \tilde{\chi}_1^0$ with 100% branching ratio, in $\tilde{t}_1$--$\tilde{\chi}_1^0$ masses plane. The dashed lines and the shaded bands are the expected limit and its $\pm1\sigma$ uncertainty. The thick solid lines are the observed limits for the central value of the signal cross-section. The expected and observed limits do not include the effect of the theoretical uncertainties in the signal cross-section. The dotted lines show the effect on the observed limit when varying the signal cross-section by $\pm1\sigma$ of the theoretical uncertainty. The observed (a) and expected (b) CLs values are respectively shown.

Exclusion limit contour (95% CL) for a simplified model assuming $\tilde{t}_1$ pair production, decaying via $\tilde{t}_1 \rightarrow b l \nu \tilde{\chi}_1^0$ with 100% branching ratio, in $\tilde{t}_1$--$\tilde{\chi}_1^0$ masses plane. The dashed lines and the shaded bands are the expected limit and its $\pm1\sigma$ uncertainty. The thick solid lines are the observed limits for the central value of the signal cross-section. The expected and observed limits do not include the effect of the theoretical uncertainties in the signal cross-section. The dotted lines show the effect on the observed limit when varying the signal cross-section by $\pm1\sigma$ of the theoretical uncertainty. The observed (a) and expected (b) CLs values are respectively shown.

Exclusion limit contour (95% CL) for a simplified model assuming $\tilde{t}_1$ pair production, decaying via $\tilde{t}_1 \rightarrow b l \nu \tilde{\chi}_1^0$ with 100% branching ratio, in $\tilde{t}_1$--$\tilde{\chi}_1^0$ masses plane. The dashed lines and the shaded bands are the expected limit and its $\pm1\sigma$ uncertainty. The thick solid lines are the observed limits for the central value of the signal cross-section. The expected and observed limits do not include the effect of the theoretical uncertainties in the signal cross-section. The dotted lines show the effect on the observed limit when varying the signal cross-section by $\pm1\sigma$ of the theoretical uncertainty. The observed (a) and expected (b) CLs values are respectively shown.

Exclusion limit contour (95% CL) for a simplified model assuming $\tilde{t}_1$ pair production, decaying via $\tilde{t}_1 \rightarrow b l \nu \tilde{\chi}_1^0$ with 100% branching ratio, in $\tilde{t}_1$--$\tilde{\chi}_1^0$ masses plane. The dashed lines and the shaded bands are the expected limit and its $\pm1\sigma$ uncertainty.The thick solid lines are the observed limits for the central value of the signal cross-section. The expected and observed limits do not include the effect of the theoretical uncertainties in the signal cross-section. The dotted lines show the effect on the observed limit when varying the signal cross-section by $\pm1\sigma$ of the theoretical uncertainty. The observed (a) and expected (b) CLs values are respectively shown.

Exclusion limit contour (95% CL) for a simplified model assuming $\tilde{t}_1$ pair production, decaying via $\tilde{t}_1 \rightarrow b l \nu \tilde{\chi}_1^0$ with 100% branching ratio, in $\tilde{t}_1$--$\tilde{\chi}_1^0$ masses plane. The dashed lines and the shaded bands are the expected limit and its $\pm1\sigma$ uncertainty.The thick solid lines are the observed limits for the central value of the signal cross-section. The expected and observed limits do not include the effect of the theoretical uncertainties in the signal cross-section. The dotted lines show the effect on the observed limit when varying the signal cross-section by $\pm1\sigma$ of the theoretical uncertainty. The observed (a) and expected (b) CLs values are respectively shown.

Exclusion limits for (a) $t\bar{t} + \phi $ scalar and (b) $t\bar{t} + a $ pseudoscalar models as a function of the DM particle mass for a mediator mass of 10 GeV. 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 shows the observed (expected) exclusion limits.

Exclusion limits for (a) $t\bar{t} + \phi $ scalar and (b) $t\bar{t} + a $ pseudoscalar models as a function of the DM particle mass for a mediator mass of 10 GeV. 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 shows the observed (expected) exclusion limits.

Exclusion limits for (a) $t\bar{t} + \phi $ scalar and (b) $t\bar{t} + a $ pseudoscalar models as a function of the DM particle mass for a mediator mass of 10 GeV. 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 shows the observed (expected) exclusion limits.

Exclusion limits for (a) $t\bar{t} + \phi $ scalar and (b) $t\bar{t} + a $ pseudoscalar models as a function of the DM particle mass for a mediator mass of 10 GeV. 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 shows the observed (expected) exclusion limits.

Two-body selection efficiency (a) SR-DF$^{2-body}_{[110,120)}$, (b) SR-DF1$^{2-body}_{[120,140)}$, (c) SR-DF2$^{2-body}_{[140,160)}$, (d) SR-DF3$^{2-body}_{[160,180)}$, (e) SR-DF4$^{2-body}_{[180,220)}$, (f) SR-DF5$^{2-body}_{[220,\infty)}$ for a simplified model assuming $ \tilde{t}_1$ pair production.

Two-body selection efficiency (a) SR-DF$^{2-body}_{[110,120)}$, (b) SR-DF1$^{2-body}_{[120,140)}$, (c) SR-DF2$^{2-body}_{[140,160)}$, (d) SR-DF3$^{2-body}_{[160,180)}$, (e) SR-DF4$^{2-body}_{[180,220)}$, (f) SR-DF5$^{2-body}_{[220,\infty)}$ for a simplified model assuming $ \tilde{t}_1$ pair production.

Two-body selection efficiency (a) SR-DF$^{2-body}_{[110,120)}$, (b) SR-DF1$^{2-body}_{[120,140)}$, (c) SR-DF2$^{2-body}_{[140,160)}$, (d) SR-DF3$^{2-body}_{[160,180)}$, (e) SR-DF4$^{2-body}_{[180,220)}$, (f) SR-DF5$^{2-body}_{[220,\infty)}$ for a simplified model assuming $ \tilde{t}_1$ pair production.

Two-body selection efficiency (a) SR-DF$^{2-body}_{[110,120)}$, (b) SR-DF1$^{2-body}_{[120,140)}$, (c) SR-DF2$^{2-body}_{[140,160)}$, (d) SR-DF3$^{2-body}_{[160,180)}$, (e) SR-DF4$^{2-body}_{[180,220)}$, (f) SR-DF5$^{2-body}_{[220,\infty)}$ for a simplified model assuming $ \tilde{t}_1$ pair production.

Two-body selection efficiency (a) SR-DF$^{2-body}_{[110,120)}$, (b) SR-DF1$^{2-body}_{[120,140)}$, (c) SR-DF2$^{2-body}_{[140,160)}$, (d) SR-DF3$^{2-body}_{[160,180)}$, (e) SR-DF4$^{2-body}_{[180,220)}$, (f) SR-DF5$^{2-body}_{[220,\infty)}$ for a simplified model assuming $ \tilde{t}_1$ pair production.

Two-body selection efficiency (a) SR-DF$^{2-body}_{[110,120)}$, (b) SR-DF1$^{2-body}_{[120,140)}$, (c) SR-DF2$^{2-body}_{[140,160)}$, (d) SR-DF3$^{2-body}_{[160,180)}$, (e) SR-DF4$^{2-body}_{[180,220)}$, (f) SR-DF5$^{2-body}_{[220,\infty)}$ for a simplified model assuming $ \tilde{t}_1$ pair production.

Two-body selection efficiency (a) SR-SF$^{2-body}_{[110,120)}$, (b) SR-SF1$^{2-body}_{[120,140)}$, (c) SR-SF2$^{2-body}_{[140,160)}$, (d) SR-SF3$^{2-body}_{[160,180)}$, (e) SR-SF4$^{2-body}_{[180,220)}$, (f) SR-SF5$^{2-body}_{[220,\infty)}$ for a simplified model assuming $ \tilde{t}_1$ pair production.

Two-body selection efficiency (a) SR-SF$^{2-body}_{[110,120)}$, (b) SR-SF1$^{2-body}_{[120,140)}$, (c) SR-SF2$^{2-body}_{[140,160)}$, (d) SR-SF3$^{2-body}_{[160,180)}$, (e) SR-SF4$^{2-body}_{[180,220)}$, (f) SR-SF5$^{2-body}_{[220,\infty)}$ for a simplified model assuming $ \tilde{t}_1$ pair production.

Two-body selection efficiency (a) SR-SF$^{2-body}_{[110,120)}$, (b) SR-SF1$^{2-body}_{[120,140)}$, (c) SR-SF2$^{2-body}_{[140,160)}$, (d) SR-SF3$^{2-body}_{[160,180)}$, (e) SR-SF4$^{2-body}_{[180,220)}$, (f) SR-SF5$^{2-body}_{[220,\infty)}$ for a simplified model assuming $ \tilde{t}_1$ pair production.

Two-body selection efficiency (a) SR-SF$^{2-body}_{[110,120)}$, (b) SR-SF1$^{2-body}_{[120,140)}$, (c) SR-SF2$^{2-body}_{[140,160)}$, (d) SR-SF3$^{2-body}_{[160,180)}$, (e) SR-SF4$^{2-body}_{[180,220)}$, (f) SR-SF5$^{2-body}_{[220,\infty)}$ for a simplified model assuming $ \tilde{t}_1$ pair production.

Two-body selection efficiency (a) SR-SF$^{2-body}_{[110,120)}$, (b) SR-SF1$^{2-body}_{[120,140)}$, (c) SR-SF2$^{2-body}_{[140,160)}$, (d) SR-SF3$^{2-body}_{[160,180)}$, (e) SR-SF4$^{2-body}_{[180,220)}$, (f) SR-SF5$^{2-body}_{[220,\infty)}$ for a simplified model assuming $ \tilde{t}_1$ pair production.

Two-body selection efficiency (a) SR-SF$^{2-body}_{[110,120)}$, (b) SR-SF1$^{2-body}_{[120,140)}$, (c) SR-SF2$^{2-body}_{[140,160)}$, (d) SR-SF3$^{2-body}_{[160,180)}$, (e) SR-SF4$^{2-body}_{[180,220)}$, (f) SR-SF5$^{2-body}_{[220,\infty)}$ for a simplified model assuming $ \tilde{t}_1$ pair production.

Two-body selection efficiency (a) $SR^{2-body}_{[110,\infty)}$ , (b) $SR^{2-body}_{[120,\infty)}$ , (c) $SR^{2-body}_{[140,\infty)}$ , (d) $SR^{2-body}_{[160,\infty)}$ , (e) $SR^{2-body}_{[180,\infty)}$ , (f) $SR^{2-body}_{[200,\infty)}$ , (g) $SR^{2-body}_{[220,\infty)}$ for a simplified model assuming $ \tilde{t}_1$ pair production.

Two-body selection efficiency (a) $SR^{2-body}_{[110,\infty)}$ , (b) $SR^{2-body}_{[120,\infty)}$ , (c) $SR^{2-body}_{[140,\infty)}$ , (d) $SR^{2-body}_{[160,\infty)}$ , (e) $SR^{2-body}_{[180,\infty)}$ , (f) $SR^{2-body}_{[200,\infty)}$ , (g) $SR^{2-body}_{[220,\infty)}$ for a simplified model assuming $ \tilde{t}_1$ pair production.

Two-body selection efficiency (a) $SR^{2-body}_{[110,\infty)}$ , (b) $SR^{2-body}_{[120,\infty)}$ , (c) $SR^{2-body}_{[140,\infty)}$ , (d) $SR^{2-body}_{[160,\infty)}$ , (e) $SR^{2-body}_{[180,\infty)}$ , (f) $SR^{2-body}_{[200,\infty)}$ , (g) $SR^{2-body}_{[220,\infty)}$ for a simplified model assuming $ \tilde{t}_1$ pair production.

Two-body selection efficiency (a) $SR^{2-body}_{[110,\infty)}$ , (b) $SR^{2-body}_{[120,\infty)}$ , (c) $SR^{2-body}_{[140,\infty)}$ , (d) $SR^{2-body}_{[160,\infty)}$ , (e) $SR^{2-body}_{[180,\infty)}$ , (f) $SR^{2-body}_{[200,\infty)}$ , (g) $SR^{2-body}_{[220,\infty)}$ for a simplified model assuming $ \tilde{t}_1$ pair production.

Two-body selection efficiency (a) $SR^{2-body}_{[110,\infty)}$ , (b) $SR^{2-body}_{[120,\infty)}$ , (c) $SR^{2-body}_{[140,\infty)}$ , (d) $SR^{2-body}_{[160,\infty)}$ , (e) $SR^{2-body}_{[180,\infty)}$ , (f) $SR^{2-body}_{[200,\infty)}$ , (g) $SR^{2-body}_{[220,\infty)}$ for a simplified model assuming $ \tilde{t}_1$ pair production.

Two-body selection efficiency (a) $SR^{2-body}_{[110,\infty)}$ , (b) $SR^{2-body}_{[120,\infty)}$ , (c) $SR^{2-body}_{[140,\infty)}$ , (d) $SR^{2-body}_{[160,\infty)}$ , (e) $SR^{2-body}_{[180,\infty)}$ , (f) $SR^{2-body}_{[200,\infty)}$ , (g) $SR^{2-body}_{[220,\infty)}$ for a simplified model assuming $ \tilde{t}_1$ pair production.

Two-body selection efficiency (a) $SR^{2-body}_{[110,\infty)}$ , (b) $SR^{2-body}_{[120,\infty)}$ , (c) $SR^{2-body}_{[140,\infty)}$ , (d) $SR^{2-body}_{[160,\infty)}$ , (e) $SR^{2-body}_{[180,\infty)}$ , (f) $SR^{2-body}_{[200,\infty)}$ , (g) $SR^{2-body}_{[220,\infty)}$ for a simplified model assuming $ \tilde{t}_1$ pair production.

Two-body selection efficiency (a) SR-DF$^{2-body}_{[110,120)}$, (b) SR-DF1$^{2-body}_{[120,140)}$, (c) SR-DF2$^{2-body}_{[140,160)}$, (d) SR-DF3$^{2-body}_{[160,180)}$, (e) SR-DF4$^{2-body}_{[180,220)}$, (f) SR-DF5$^{2-body}_{[220,\infty)}$ for a simplified model assuming $t \tilde{t} + \phi$ pair production.

Two-body selection efficiency (a) SR-DF$^{2-body}_{[110,120)}$, (b) SR-DF1$^{2-body}_{[120,140)}$, (c) SR-DF2$^{2-body}_{[140,160)}$, (d) SR-DF3$^{2-body}_{[160,180)}$, (e) SR-DF4$^{2-body}_{[180,220)}$, (f) SR-DF5$^{2-body}_{[220,\infty)}$ for a simplified model assuming $t \tilde{t} + \phi$ pair production.

Two-body selection efficiency (a) SR-DF$^{2-body}_{[110,120)}$, (b) SR-DF1$^{2-body}_{[120,140)}$, (c) SR-DF2$^{2-body}_{[140,160)}$, (d) SR-DF3$^{2-body}_{[160,180)}$, (e) SR-DF4$^{2-body}_{[180,220)}$, (f) SR-DF5$^{2-body}_{[220,\infty)}$ for a simplified model assuming $t \tilde{t} + \phi$ pair production.

Two-body selection efficiency (a) SR-DF$^{2-body}_{[110,120)}$, (b) SR-DF1$^{2-body}_{[120,140)}$, (c) SR-DF2$^{2-body}_{[140,160)}$, (d) SR-DF3$^{2-body}_{[160,180)}$, (e) SR-DF4$^{2-body}_{[180,220)}$, (f) SR-DF5$^{2-body}_{[220,\infty)}$ for a simplified model assuming $t \tilde{t} + \phi$ pair production.

Two-body selection efficiency (a) SR-DF$^{2-body}_{[110,120)}$, (b) SR-DF1$^{2-body}_{[120,140)}$, (c) SR-DF2$^{2-body}_{[140,160)}$, (d) SR-DF3$^{2-body}_{[160,180)}$, (e) SR-DF4$^{2-body}_{[180,220)}$, (f) SR-DF5$^{2-body}_{[220,\infty)}$ for a simplified model assuming $t \tilde{t} + \phi$ pair production.

Two-body selection efficiency (a) SR-DF$^{2-body}_{[110,120)}$, (b) SR-DF1$^{2-body}_{[120,140)}$, (c) SR-DF2$^{2-body}_{[140,160)}$, (d) SR-DF3$^{2-body}_{[160,180)}$, (e) SR-DF4$^{2-body}_{[180,220)}$, (f) SR-DF5$^{2-body}_{[220,\infty)}$ for a simplified model assuming $t \tilde{t} + \phi$ pair production.

Two-body selection efficiency (a) SR-SF$^{2-body}_{[110,120)}$, (b) SR-SF1$^{2-body}_{[120,140)}$, (c) SR-SF2$^{2-body}_{[140,160)}$, (d) SR-SF3$^{2-body}_{[160,180)}$, (e) SR-SF4$^{2-body}_{[180,220)}$, (f) SR-SF5$^{2-body}_{[220,\infty)}$ for a simplified model assuming $t \tilde{t} + \phi$ pair production.

Two-body selection efficiency (a) SR-SF$^{2-body}_{[110,120)}$, (b) SR-SF1$^{2-body}_{[120,140)}$, (c) SR-SF2$^{2-body}_{[140,160)}$, (d) SR-SF3$^{2-body}_{[160,180)}$, (e) SR-SF4$^{2-body}_{[180,220)}$, (f) SR-SF5$^{2-body}_{[220,\infty)}$ for a simplified model assuming $t \tilde{t} + \phi$ pair production.

Two-body selection efficiency (a) SR-SF$^{2-body}_{[110,120)}$, (b) SR-SF1$^{2-body}_{[120,140)}$, (c) SR-SF2$^{2-body}_{[140,160)}$, (d) SR-SF3$^{2-body}_{[160,180)}$, (e) SR-SF4$^{2-body}_{[180,220)}$, (f) SR-SF5$^{2-body}_{[220,\infty)}$ for a simplified model assuming $t \tilde{t} + \phi$ pair production.

Two-body selection efficiency (a) SR-SF$^{2-body}_{[110,120)}$, (b) SR-SF1$^{2-body}_{[120,140)}$, (c) SR-SF2$^{2-body}_{[140,160)}$, (d) SR-SF3$^{2-body}_{[160,180)}$, (e) SR-SF4$^{2-body}_{[180,220)}$, (f) SR-SF5$^{2-body}_{[220,\infty)}$ for a simplified model assuming $t \tilde{t} + \phi$ pair production.

Two-body selection efficiency (a) SR-SF$^{2-body}_{[110,120)}$, (b) SR-SF1$^{2-body}_{[120,140)}$, (c) SR-SF2$^{2-body}_{[140,160)}$, (d) SR-SF3$^{2-body}_{[160,180)}$, (e) SR-SF4$^{2-body}_{[180,220)}$, (f) SR-SF5$^{2-body}_{[220,\infty)}$ for a simplified model assuming $t \tilde{t} + \phi$ pair production.

Two-body selection efficiency (a) SR-SF$^{2-body}_{[110,120)}$, (b) SR-SF1$^{2-body}_{[120,140)}$, (c) SR-SF2$^{2-body}_{[140,160)}$, (d) SR-SF3$^{2-body}_{[160,180)}$, (e) SR-SF4$^{2-body}_{[180,220)}$, (f) SR-SF5$^{2-body}_{[220,\infty)}$ for a simplified model assuming $t \tilde{t} + \phi$ pair production.

Two-body selection efficiency (a) $SR^{2-body}_{[110,\infty)}$ , (b) $SR^{2-body}_{[120,\infty)}$ , (c) $SR^{2-body}_{[140,\infty)}$ , (d) $SR^{2-body}_{[160,\infty)}$ , (e) $SR^{2-body}_{[180,\infty)}$ , (f) $SR^{2-body}_{[200,\infty)}$ , (g) $SR^{2-body}_{[220,\infty)}$ for a simplified model assuming $t \tilde{t} + \phi$ pair production.

Two-body selection efficiency (a) $SR^{2-body}_{[110,\infty)}$ , (b) $SR^{2-body}_{[120,\infty)}$ , (c) $SR^{2-body}_{[140,\infty)}$ , (d) $SR^{2-body}_{[160,\infty)}$ , (e) $SR^{2-body}_{[180,\infty)}$ , (f) $SR^{2-body}_{[200,\infty)}$ , (g) $SR^{2-body}_{[220,\infty)}$ for a simplified model assuming $t \tilde{t} + \phi$ pair production.

Two-body selection efficiency (a) $SR^{2-body}_{[110,\infty)}$ , (b) $SR^{2-body}_{[120,\infty)}$ , (c) $SR^{2-body}_{[140,\infty)}$ , (d) $SR^{2-body}_{[160,\infty)}$ , (e) $SR^{2-body}_{[180,\infty)}$ , (f) $SR^{2-body}_{[200,\infty)}$ , (g) $SR^{2-body}_{[220,\infty)}$ for a simplified model assuming $t \tilde{t} + \phi$ pair production.

Two-body selection efficiency (a) $SR^{2-body}_{[110,\infty)}$ , (b) $SR^{2-body}_{[120,\infty)}$ , (c) $SR^{2-body}_{[140,\infty)}$ , (d) $SR^{2-body}_{[160,\infty)}$ , (e) $SR^{2-body}_{[180,\infty)}$ , (f) $SR^{2-body}_{[200,\infty)}$ , (g) $SR^{2-body}_{[220,\infty)}$ for a simplified model assuming $t \tilde{t} + \phi$ pair production.

Two-body selection efficiency (a) $SR^{2-body}_{[110,\infty)}$ , (b) $SR^{2-body}_{[120,\infty)}$ , (c) $SR^{2-body}_{[140,\infty)}$ , (d) $SR^{2-body}_{[160,\infty)}$ , (e) $SR^{2-body}_{[180,\infty)}$ , (f) $SR^{2-body}_{[200,\infty)}$ , (g) $SR^{2-body}_{[220,\infty)}$ for a simplified model assuming $t \tilde{t} + \phi$ pair production.

Two-body selection efficiency (a) $SR^{2-body}_{[110,\infty)}$ , (b) $SR^{2-body}_{[120,\infty)}$ , (c) $SR^{2-body}_{[140,\infty)}$ , (d) $SR^{2-body}_{[160,\infty)}$ , (e) $SR^{2-body}_{[180,\infty)}$ , (f) $SR^{2-body}_{[200,\infty)}$ , (g) $SR^{2-body}_{[220,\infty)}$ for a simplified model assuming $t \tilde{t} + \phi$ pair production.

Two-body selection efficiency (a) $SR^{2-body}_{[110,\infty)}$ , (b) $SR^{2-body}_{[120,\infty)}$ , (c) $SR^{2-body}_{[140,\infty)}$ , (d) $SR^{2-body}_{[160,\infty)}$ , (e) $SR^{2-body}_{[180,\infty)}$ , (f) $SR^{2-body}_{[200,\infty)}$ , (g) $SR^{2-body}_{[220,\infty)}$ for a simplified model assuming $t \tilde{t} + \phi$ pair production.

Two-body selection efficiency (a) SR-DF$^{2-body}_{[110,120)}$, (b) SR-DF1$^{2-body}_{[120,140)}$, (c) SR-DF2$^{2-body}_{[140,160)}$, (d) SR-DF3$^{2-body}_{[160,180)}$, (e) SR-DF4$^{2-body}_{[180,220)}$, (f) SR-DF5$^{2-body}_{[220,\infty)}$ for a simplified model assuming $t \tilde{t} + a$ pair production.

Two-body selection efficiency (a) SR-DF$^{2-body}_{[110,120)}$, (b) SR-DF1$^{2-body}_{[120,140)}$, (c) SR-DF2$^{2-body}_{[140,160)}$, (d) SR-DF3$^{2-body}_{[160,180)}$, (e) SR-DF4$^{2-body}_{[180,220)}$, (f) SR-DF5$^{2-body}_{[220,\infty)}$ for a simplified model assuming $t \tilde{t} + a$ pair production.

Two-body selection efficiency (a) SR-DF$^{2-body}_{[110,120)}$, (b) SR-DF1$^{2-body}_{[120,140)}$, (c) SR-DF2$^{2-body}_{[140,160)}$, (d) SR-DF3$^{2-body}_{[160,180)}$, (e) SR-DF4$^{2-body}_{[180,220)}$, (f) SR-DF5$^{2-body}_{[220,\infty)}$ for a simplified model assuming $t \tilde{t} + a$ pair production.

Two-body selection efficiency (a) SR-DF$^{2-body}_{[110,120)}$, (b) SR-DF1$^{2-body}_{[120,140)}$, (c) SR-DF2$^{2-body}_{[140,160)}$, (d) SR-DF3$^{2-body}_{[160,180)}$, (e) SR-DF4$^{2-body}_{[180,220)}$, (f) SR-DF5$^{2-body}_{[220,\infty)}$ for a simplified model assuming $t \tilde{t} + a$ pair production.

Two-body selection efficiency (a) SR-DF$^{2-body}_{[110,120)}$, (b) SR-DF1$^{2-body}_{[120,140)}$, (c) SR-DF2$^{2-body}_{[140,160)}$, (d) SR-DF3$^{2-body}_{[160,180)}$, (e) SR-DF4$^{2-body}_{[180,220)}$, (f) SR-DF5$^{2-body}_{[220,\infty)}$ for a simplified model assuming $t \tilde{t} + a$ pair production.

Two-body selection efficiency (a) SR-DF$^{2-body}_{[110,120)}$, (b) SR-DF1$^{2-body}_{[120,140)}$, (c) SR-DF2$^{2-body}_{[140,160)}$, (d) SR-DF3$^{2-body}_{[160,180)}$, (e) SR-DF4$^{2-body}_{[180,220)}$, (f) SR-DF5$^{2-body}_{[220,\infty)}$ for a simplified model assuming $t \tilde{t} + a$ pair production.

Two-body selection efficiency (a) SR-SF$^{2-body}_{[110,120)}$, (b) SR-SF1$^{2-body}_{[120,140)}$, (c) SR-SF2$^{2-body}_{[140,160)}$, (d) SR-SF3$^{2-body}_{[160,180)}$, (e) SR-SF4$^{2-body}_{[180,220)}$, (f) SR-SF5$^{2-body}_{[220,\infty)}$ for a simplified model assuming $t \tilde{t} + a$ pair production.

Two-body selection efficiency (a) SR-SF$^{2-body}_{[110,120)}$, (b) SR-SF1$^{2-body}_{[120,140)}$, (c) SR-SF2$^{2-body}_{[140,160)}$, (d) SR-SF3$^{2-body}_{[160,180)}$, (e) SR-SF4$^{2-body}_{[180,220)}$, (f) SR-SF5$^{2-body}_{[220,\infty)}$ for a simplified model assuming $t \tilde{t} + a$ pair production.

Two-body selection efficiency (a) SR-SF$^{2-body}_{[110,120)}$, (b) SR-SF1$^{2-body}_{[120,140)}$, (c) SR-SF2$^{2-body}_{[140,160)}$, (d) SR-SF3$^{2-body}_{[160,180)}$, (e) SR-SF4$^{2-body}_{[180,220)}$, (f) SR-SF5$^{2-body}_{[220,\infty)}$ for a simplified model assuming $t \tilde{t} + a$ pair production.

Two-body selection efficiency (a) SR-SF$^{2-body}_{[110,120)}$, (b) SR-SF1$^{2-body}_{[120,140)}$, (c) SR-SF2$^{2-body}_{[140,160)}$, (d) SR-SF3$^{2-body}_{[160,180)}$, (e) SR-SF4$^{2-body}_{[180,220)}$, (f) SR-SF5$^{2-body}_{[220,\infty)}$ for a simplified model assuming $t \tilde{t} + a$ pair production.

Two-body selection efficiency (a) SR-SF$^{2-body}_{[110,120)}$, (b) SR-SF1$^{2-body}_{[120,140)}$, (c) SR-SF2$^{2-body}_{[140,160)}$, (d) SR-SF3$^{2-body}_{[160,180)}$, (e) SR-SF4$^{2-body}_{[180,220)}$, (f) SR-SF5$^{2-body}_{[220,\infty)}$ for a simplified model assuming $t \tilde{t} + a$ pair production.

Two-body selection efficiency (a) SR-SF$^{2-body}_{[110,120)}$, (b) SR-SF1$^{2-body}_{[120,140)}$, (c) SR-SF2$^{2-body}_{[140,160)}$, (d) SR-SF3$^{2-body}_{[160,180)}$, (e) SR-SF4$^{2-body}_{[180,220)}$, (f) SR-SF5$^{2-body}_{[220,\infty)}$ for a simplified model assuming $t \tilde{t} + a$ pair production.

Two-body selection efficiency (a) $SR^{2-body}_{[110,\infty)}$ , (b) $SR^{2-body}_{[120,\infty)}$ , (c) $SR^{2-body}_{[140,\infty)}$ , (d) $SR^{2-body}_{[160,\infty)}$ , (e) $SR^{2-body}_{[180,\infty)}$ , (f) $SR^{2-body}_{[200,\infty)}$ , (g) $SR^{2-body}_{[220,\infty)}$ for a simplified model assuming $t \tilde{t} + a$ pair production.

Two-body selection efficiency (a) $SR^{2-body}_{[110,\infty)}$ , (b) $SR^{2-body}_{[120,\infty)}$ , (c) $SR^{2-body}_{[140,\infty)}$ , (d) $SR^{2-body}_{[160,\infty)}$ , (e) $SR^{2-body}_{[180,\infty)}$ , (f) $SR^{2-body}_{[200,\infty)}$ , (g) $SR^{2-body}_{[220,\infty)}$ for a simplified model assuming $t \tilde{t} + a$ pair production.

Two-body selection efficiency (a) $SR^{2-body}_{[110,\infty)}$ , (b) $SR^{2-body}_{[120,\infty)}$ , (c) $SR^{2-body}_{[140,\infty)}$ , (d) $SR^{2-body}_{[160,\infty)}$ , (e) $SR^{2-body}_{[180,\infty)}$ , (f) $SR^{2-body}_{[200,\infty)}$ , (g) $SR^{2-body}_{[220,\infty)}$ for a simplified model assuming $t \tilde{t} + a$ pair production.

Two-body selection efficiency (a) $SR^{2-body}_{[110,\infty)}$ , (b) $SR^{2-body}_{[120,\infty)}$ , (c) $SR^{2-body}_{[140,\infty)}$ , (d) $SR^{2-body}_{[160,\infty)}$ , (e) $SR^{2-body}_{[180,\infty)}$ , (f) $SR^{2-body}_{[200,\infty)}$ , (g) $SR^{2-body}_{[220,\infty)}$ for a simplified model assuming $t \tilde{t} + a$ pair production.

Two-body selection efficiency (a) $SR^{2-body}_{[110,\infty)}$ , (b) $SR^{2-body}_{[120,\infty)}$ , (c) $SR^{2-body}_{[140,\infty)}$ , (d) $SR^{2-body}_{[160,\infty)}$ , (e) $SR^{2-body}_{[180,\infty)}$ , (f) $SR^{2-body}_{[200,\infty)}$ , (g) $SR^{2-body}_{[220,\infty)}$ for a simplified model assuming $t \tilde{t} + a$ pair production.

Two-body selection efficiency (a) $SR^{2-body}_{[110,\infty)}$ , (b) $SR^{2-body}_{[120,\infty)}$ , (c) $SR^{2-body}_{[140,\infty)}$ , (d) $SR^{2-body}_{[160,\infty)}$ , (e) $SR^{2-body}_{[180,\infty)}$ , (f) $SR^{2-body}_{[200,\infty)}$ , (g) $SR^{2-body}_{[220,\infty)}$ for a simplified model assuming $t \tilde{t} + a$ pair production.

Two-body selection efficiency (a) $SR^{2-body}_{[110,\infty)}$ , (b) $SR^{2-body}_{[120,\infty)}$ , (c) $SR^{2-body}_{[140,\infty)}$ , (d) $SR^{2-body}_{[160,\infty)}$ , (e) $SR^{2-body}_{[180,\infty)}$ , (f) $SR^{2-body}_{[200,\infty)}$ , (g) $SR^{2-body}_{[220,\infty)}$ for a simplified model assuming $t \tilde{t} + a$ pair production.

Three-body selection efficiency (a) SR-DF$^{3-body}_{t}$, (b) SR-SF$^{3-body}_{t}$, (c) SR-DF$^{3-body}_{W}$, (d) SR-SF$^{3-body}_{W}$ for a simplified model assuming $ \tilde{t}_1$ pair production.

Three-body selection efficiency (a) SR-DF$^{3-body}_{t}$, (b) SR-SF$^{3-body}_{t}$, (c) SR-DF$^{3-body}_{W}$, (d) SR-SF$^{3-body}_{W}$ for a simplified model assuming $ \tilde{t}_1$ pair production.

Three-body selection efficiency (a) SR-DF$^{3-body}_{t}$, (b) SR-SF$^{3-body}_{t}$, (c) SR-DF$^{3-body}_{W}$, (d) SR-SF$^{3-body}_{W}$ for a simplified model assuming $ \tilde{t}_1$ pair production.

Three-body selection efficiency (a) SR-DF$^{3-body}_{t}$, (b) SR-SF$^{3-body}_{t}$, (c) SR-DF$^{3-body}_{W}$, (d) SR-SF$^{3-body}_{W}$ for a simplified model assuming $ \tilde{t}_1$ pair production.

Four-body selection Efficiency (a) SR$^{4-body}_{Small \Delta m}$ , (b) $SR^{4-body}_{Large \Delta m}$ for a simplified model assuming $\tilde{t}_1$ pair production.

Four-body selection Efficiency (a) SR$^{4-body}_{Small \Delta m}$ , (b) $SR^{4-body}_{Large \Delta\ m}$ for a simplified model assuming $\tilde{t}_1$ pair production.

Two-body selection acceptance (a) SR-DF$^{2-body}_{[110,120)}$, (b) SR-DF1$^{2-body}_{[120,140)}$, (c) SR-DF2$^{2-body}_{[140,160)}$, (d) SR-DF3$^{2-body}_{[160,180)}$, (e) SR-DF4$^{2-body}_{[180,220)}$, (f) SR-DF5$^{2-body}_{[220,\infty)}$ for a simplified model assuming $ \tilde{t}_1$ pair production.

Two-body selection acceptance (a) SR-DF$^{2-body}_{[110,120)}$, (b) SR-DF1$^{2-body}_{[120,140)}$, (c) SR-DF2$^{2-body}_{[140,160)}$, (d) SR-DF3$^{2-body}_{[160,180)}$, (e) SR-DF4$^{2-body}_{[180,220)}$, (f) SR-DF5$^{2-body}_{[220,\infty)}$ for a simplified model assuming $ \tilde{t}_1$ pair production.

Two-body selection acceptance (a) SR-DF$^{2-body}_{[110,120)}$, (b) SR-DF1$^{2-body}_{[120,140)}$, (c) SR-DF2$^{2-body}_{[140,160)}$, (d) SR-DF3$^{2-body}_{[160,180)}$, (e) SR-DF4$^{2-body}_{[180,220)}$, (f) SR-DF5$^{2-body}_{[220,\infty)}$ for a simplified model assuming $ \tilde{t}_1$ pair production.

Two-body selection acceptance (a) SR-DF$^{2-body}_{[110,120)}$, (b) SR-DF1$^{2-body}_{[120,140)}$, (c) SR-DF2$^{2-body}_{[140,160)}$, (d) SR-DF3$^{2-body}_{[160,180)}$, (e) SR-DF4$^{2-body}_{[180,220)}$, (f) SR-DF5$^{2-body}_{[220,\infty)}$ for a simplified model assuming $ \tilde{t}_1$ pair production.

Two-body selection acceptance (a) SR-DF$^{2-body}_{[110,120)}$, (b) SR-DF1$^{2-body}_{[120,140)}$, (c) SR-DF2$^{2-body}_{[140,160)}$, (d) SR-DF3$^{2-body}_{[160,180)}$, (e) SR-DF4$^{2-body}_{[180,220)}$, (f) SR-DF5$^{2-body}_{[220,\infty)}$ for a simplified model assuming $ \tilde{t}_1$ pair production.

Two-body selection acceptance (a) SR-DF$^{2-body}_{[110,120)}$, (b) SR-DF1$^{2-body}_{[120,140)}$, (c) SR-DF2$^{2-body}_{[140,160)}$, (d) SR-DF3$^{2-body}_{[160,180)}$, (e) SR-DF4$^{2-body}_{[180,220)}$, (f) SR-DF5$^{2-body}_{[220,\infty)}$ for a simplified model assuming $ \tilde{t}_1$ pair production.

Two-body selection acceptance (a) SR-SF$^{2-body}_{[110,120)}$, (b) SR-SF1$^{2-body}_{[120,140)}$, (c) SR-SF2$^{2-body}_{[140,160)}$, (d) SR-SF3$^{2-body}_{[160,180)}$, (e) SR-SF4$^{2-body}_{[180,220)}$, (f) SR-SF5$^{2-body}_{[220,\infty)}$ for a simplified model assuming $ \tilde{t}_1$ pair production.

Two-body selection acceptance (a) SR-SF$^{2-body}_{[110,120)}$, (b) SR-SF1$^{2-body}_{[120,140)}$, (c) SR-SF2$^{2-body}_{[140,160)}$, (d) SR-SF3$^{2-body}_{[160,180)}$, (e) SR-SF4$^{2-body}_{[180,220)}$, (f) SR-SF5$^{2-body}_{[220,\infty)}$ for a simplified model assuming $ \tilde{t}_1$ pair production.

Two-body selection acceptance (a) SR-SF$^{2-body}_{[110,120)}$, (b) SR-SF1$^{2-body}_{[120,140)}$, (c) SR-SF2$^{2-body}_{[140,160)}$, (d) SR-SF3$^{2-body}_{[160,180)}$, (e) SR-SF4$^{2-body}_{[180,220)}$, (f) SR-SF5$^{2-body}_{[220,\infty)}$ for a simplified model assuming $ \tilde{t}_1$ pair production.

Two-body selection acceptance (a) SR-SF$^{2-body}_{[110,120)}$, (b) SR-SF1$^{2-body}_{[120,140)}$, (c) SR-SF2$^{2-body}_{[140,160)}$, (d) SR-SF3$^{2-body}_{[160,180)}$, (e) SR-SF4$^{2-body}_{[180,220)}$, (f) SR-SF5$^{2-body}_{[220,\infty)}$ for a simplified model assuming $ \tilde{t}_1$ pair production.

Two-body selection acceptance (a) SR-SF$^{2-body}_{[110,120)}$, (b) SR-SF1$^{2-body}_{[120,140)}$, (c) SR-SF2$^{2-body}_{[140,160)}$, (d) SR-SF3$^{2-body}_{[160,180)}$, (e) SR-SF4$^{2-body}_{[180,220)}$, (f) SR-SF5$^{2-body}_{[220,\infty)}$ for a simplified model assuming $ \tilde{t}_1$ pair production.

Two-body selection acceptance (a) SR-SF$^{2-body}_{[110,120)}$, (b) SR-SF1$^{2-body}_{[120,140)}$, (c) SR-SF2$^{2-body}_{[140,160)}$, (d) SR-SF3$^{2-body}_{[160,180)}$, (e) SR-SF4$^{2-body}_{[180,220)}$, (f) SR-SF5$^{2-body}_{[220,\infty)}$ for a simplified model assuming $ \tilde{t}_1$ pair production.

Two-body selection acceptance (a) $SR^{2-body}_{[110,\infty)}$ , (b) $SR^{2-body}_{[120,\infty)}$ , (c) $SR^{2-body}_{[140,\infty)}$ , (d) $SR^{2-body}_{[160,\infty)}$ , (e) $SR^{2-body}_{[180,\infty)}$ , (f) $SR^{2-body}_{[200,\infty)}$ , (g) $SR^{2-body}_{[220,\infty)}$ for a simplified model assuming $ \tilde{t}_1$ pair production.

Two-body selection acceptance (a) $SR^{2-body}_{[110,\infty)}$ , (b) $SR^{2-body}_{[120,\infty)}$ , (c) $SR^{2-body}_{[140,\infty)}$ , (d) $SR^{2-body}_{[160,\infty)}$ , (e) $SR^{2-body}_{[180,\infty)}$ , (f) $SR^{2-body}_{[200,\infty)}$ , (g) $SR^{2-body}_{[220,\infty)}$ for a simplified model assuming $ \tilde{t}_1$ pair production.

Two-body selection acceptance (a) $SR^{2-body}_{[110,\infty)}$ , (b) $SR^{2-body}_{[120,\infty)}$ , (c) $SR^{2-body}_{[140,\infty)}$ , (d) $SR^{2-body}_{[160,\infty)}$ , (e) $SR^{2-body}_{[180,\infty)}$ , (f) $SR^{2-body}_{[200,\infty)}$ , (g) $SR^{2-body}_{[220,\infty)}$ for a simplified model assuming $ \tilde{t}_1$ pair production.

Two-body selection acceptance (a) $SR^{2-body}_{[110,\infty)}$ , (b) $SR^{2-body}_{[120,\infty)}$ , (c) $SR^{2-body}_{[140,\infty)}$ , (d) $SR^{2-body}_{[160,\infty)}$ , (e) $SR^{2-body}_{[180,\infty)}$ , (f) $SR^{2-body}_{[200,\infty)}$ , (g) $SR^{2-body}_{[220,\infty)}$ for a simplified model assuming $ \tilde{t}_1$ pair production.

Two-body selection acceptance (a) $SR^{2-body}_{[110,\infty)}$ , (b) $SR^{2-body}_{[120,\infty)}$ , (c) $SR^{2-body}_{[140,\infty)}$ , (d) $SR^{2-body}_{[160,\infty)}$ , (e) $SR^{2-body}_{[180,\infty)}$ , (f) $SR^{2-body}_{[200,\infty)}$ , (g) $SR^{2-body}_{[220,\infty)}$ for a simplified model assuming $ \tilde{t}_1$ pair production.

Two-body selection acceptance (a) $SR^{2-body}_{[110,\infty)}$ , (b) $SR^{2-body}_{[120,\infty)}$ , (c) $SR^{2-body}_{[140,\infty)}$ , (d) $SR^{2-body}_{[160,\infty)}$ , (e) $SR^{2-body}_{[180,\infty)}$ , (f) $SR^{2-body}_{[200,\infty)}$ , (g) $SR^{2-body}_{[220,\infty)}$ for a simplified model assuming $ \tilde{t}_1$ pair production.

Two-body selection acceptance (a) $SR^{2-body}_{[110,\infty)}$ , (b) $SR^{2-body}_{[120,\infty)}$ , (c) $SR^{2-body}_{[140,\infty)}$ , (d) $SR^{2-body}_{[160,\infty)}$ , (e) $SR^{2-body}_{[180,\infty)}$ , (f) $SR^{2-body}_{[200,\infty)}$ , (g) $SR^{2-body}_{[220,\infty)}$ for a simplified model assuming $ \tilde{t}_1$ pair production.

Two-body selection acceptance (a) SR-DF$^{2-body}_{[110,120)}$, (b) SR-DF1$^{2-body}_{[120,140)}$, (c) SR-DF2$^{2-body}_{[140,160)}$, (d) SR-DF3$^{2-body}_{[160,180)}$, (e) SR-DF4$^{2-body}_{[180,220)}$, (f) SR-DF5$^{2-body}_{[220,\infty)}$ for a simplified model assuming $t \tilde{t} + \phi$ pair production.

Two-body selection acceptance (a) SR-DF$^{2-body}_{[110,120)}$, (b) SR-DF1$^{2-body}_{[120,140)}$, (c) SR-DF2$^{2-body}_{[140,160)}$, (d) SR-DF3$^{2-body}_{[160,180)}$, (e) SR-DF4$^{2-body}_{[180,220)}$, (f) SR-DF5$^{2-body}_{[220,\infty)}$ for a simplified model assuming $t \tilde{t} + \phi$ pair production.

Two-body selection acceptance (a) SR-DF$^{2-body}_{[110,120)}$, (b) SR-DF1$^{2-body}_{[120,140)}$, (c) SR-DF2$^{2-body}_{[140,160)}$, (d) SR-DF3$^{2-body}_{[160,180)}$, (e) SR-DF4$^{2-body}_{[180,220)}$, (f) SR-DF5$^{2-body}_{[220,\infty)}$ for a simplified model assuming $t \tilde{t} + \phi$ pair production.

Two-body selection acceptance (a) SR-DF$^{2-body}_{[110,120)}$, (b) SR-DF1$^{2-body}_{[120,140)}$, (c) SR-DF2$^{2-body}_{[140,160)}$, (d) SR-DF3$^{2-body}_{[160,180)}$, (e) SR-DF4$^{2-body}_{[180,220)}$, (f) SR-DF5$^{2-body}_{[220,\infty)}$ for a simplified model assuming $t \tilde{t} + \phi$ pair production.

Two-body selection acceptance (a) SR-DF$^{2-body}_{[110,120)}$, (b) SR-DF1$^{2-body}_{[120,140)}$, (c) SR-DF2$^{2-body}_{[140,160)}$, (d) SR-DF3$^{2-body}_{[160,180)}$, (e) SR-DF4$^{2-body}_{[180,220)}$, (f) SR-DF5$^{2-body}_{[220,\infty)}$ for a simplified model assuming $t \tilde{t} + \phi$ pair production.

Two-body selection acceptance (a) SR-DF$^{2-body}_{[110,120)}$, (b) SR-DF1$^{2-body}_{[120,140)}$, (c) SR-DF2$^{2-body}_{[140,160)}$, (d) SR-DF3$^{2-body}_{[160,180)}$, (e) SR-DF4$^{2-body}_{[180,220)}$, (f) SR-DF5$^{2-body}_{[220,\infty)}$ for a simplified model assuming $t \tilde{t} + \phi$ pair production.

Two-body selection acceptance (a) SR-SF$^{2-body}_{[110,120)}$, (b) SR-SF1$^{2-body}_{[120,140)}$, (c) SR-SF2$^{2-body}_{[140,160)}$, (d) SR-SF3$^{2-body}_{[160,180)}$, (e) SR-SF4$^{2-body}_{[180,220)}$, (f) SR-SF5$^{2-body}_{[220,\infty)}$ for a simplified model assuming $t \tilde{t} +\phi$ pair production.

Two-body selection acceptance (a) SR-SF$^{2-body}_{[110,120)}$, (b) SR-SF1$^{2-body}_{[120,140)}$, (c) SR-SF2$^{2-body}_{[140,160)}$, (d) SR-SF3$^{2-body}_{[160,180)}$, (e) SR-SF4$^{2-body}_{[180,220)}$, (f) SR-SF5$^{2-body}_{[220,\infty)}$ for a simplified model assuming $ t \tilde{t} +\phi$ pair production.

Two-body selection acceptance (a) SR-SF$^{2-body}_{[110,120)}$, (b) SR-SF1$^{2-body}_{[120,140)}$, (c) SR-SF2$^{2-body}_{[140,160)}$, (d) SR-SF3$^{2-body}_{[160,180)}$, (e) SR-SF4$^{2-body}_{[180,220)}$, (f) SR-SF5$^{2-body}_{[220,\infty)}$ for a simplified model assuming $ t \tilde{t} +\phi$ pair production.

Two-body selection acceptance (a) SR-SF$^{2-body}_{[110,120)}$, (b) SR-SF1$^{2-body}_{[120,140)}$, (c) SR-SF2$^{2-body}_{[140,160)}$, (d) SR-SF3$^{2-body}_{[160,180)}$, (e) SR-SF4$^{2-body}_{[180,220)}$, (f) SR-SF5$^{2-body}_{[220,\infty)}$ for a simplified model assuming $ t \tilde{t} +\phi$ pair production.

Two-body selection acceptance (a) SR-SF$^{2-body}_{[110,120)}$, (b) SR-SF1$^{2-body}_{[120,140)}$, (c) SR-SF2$^{2-body}_{[140,160)}$, (d) SR-SF3$^{2-body}_{[160,180)}$, (e) SR-SF4$^{2-body}_{[180,220)}$, (f) SR-SF5$^{2-body}_{[220,\infty)}$ for a simplified model assuming $ t \tilde{t} +\phi$ pair production.

Two-body selection acceptance (a) SR-SF$^{2-body}_{[110,120)}$, (b) SR-SF1$^{2-body}_{[120,140)}$, (c) SR-SF2$^{2-body}_{[140,160)}$, (d) SR-SF3$^{2-body}_{[160,180)}$, (e) SR-SF4$^{2-body}_{[180,220)}$, (f) SR-SF5$^{2-body}_{[220,\infty)}$ for a simplified model assuming $ t \tilde{t} +\phi$ pair production.

Two-body selection acceptance (a) $SR^{2-body}_{[110,\infty)}$ , (b) $SR^{2-body}_{[120,\infty)}$ , (c) $SR^{2-body}_{[140,\infty)}$ , (d) $SR^{2-body}_{[160,\infty)}$ , (e) $SR^{2-body}_{[180,\infty)}$ , (f) $SR^{2-body}_{[200,\infty)}$ , (g) $SR^{2-body}_{[220,\infty)}$ for a simplified model assuming $t \tilde{t} + \phi$ pair production.

Two-body selection acceptance (a) $SR^{2-body}_{[110,\infty)}$ , (b) $SR^{2-body}_{[120,\infty)}$ , (c) $SR^{2-body}_{[140,\infty)}$ , (d) $SR^{2-body}_{[160,\infty)}$ , (e) $SR^{2-body}_{[180,\infty)}$ , (f) $SR^{2-body}_{[200,\infty)}$ , (g) $SR^{2-body}_{[220,\infty)}$ for a simplified model assuming $t \tilde{t} + \phi$ pair production.

Two-body selection acceptance (a) $SR^{2-body}_{[110,\infty)}$ , (b) $SR^{2-body}_{[120,\infty)}$ , (c) $SR^{2-body}_{[140,\infty)}$ , (d) $SR^{2-body}_{[160,\infty)}$ , (e) $SR^{2-body}_{[180,\infty)}$ , (f) $SR^{2-body}_{[200,\infty)}$ , (g) $SR^{2-body}_{[220,\infty)}$ for a simplified model assuming $t \tilde{t} + \phi$ pair production.

Two-body selection acceptance (a) $SR^{2-body}_{[110,\infty)}$ , (b) $SR^{2-body}_{[120,\infty)}$ , (c) $SR^{2-body}_{[140,\infty)}$ , (d) $SR^{2-body}_{[160,\infty)}$ , (e) $SR^{2-body}_{[180,\infty)}$ , (f) $SR^{2-body}_{[200,\infty)}$ , (g) $SR^{2-body}_{[220,\infty)}$ for a simplified model assuming $t \tilde{t} + \phi$ pair production.

Two-body selection acceptance (a) $SR^{2-body}_{[110,\infty)}$ , (b) $SR^{2-body}_{[120,\infty)}$ , (c) $SR^{2-body}_{[140,\infty)}$ , (d) $SR^{2-body}_{[160,\infty)}$ , (e) $SR^{2-body}_{[180,\infty)}$ , (f) $SR^{2-body}_{[200,\infty)}$ , (g) $SR^{2-body}_{[220,\infty)}$ for a simplified model assuming $t \tilde{t} + \phi$ pair production.

Two-body selection acceptance (a) $SR^{2-body}_{[110,\infty)}$ , (b) $SR^{2-body}_{[120,\infty)}$ , (c) $SR^{2-body}_{[140,\infty)}$ , (d) $SR^{2-body}_{[160,\infty)}$ , (e) $SR^{2-body}_{[180,\infty)}$ , (f) $SR^{2-body}_{[200,\infty)}$ , (g) $SR^{2-body}_{[220,\infty)}$ for a simplified model assuming $t \tilde{t} + \phi$ pair production.

Two-body selection acceptance (a) $SR^{2-body}_{[110,\infty)}$ , (b) $SR^{2-body}_{[120,\infty)}$ , (c) $SR^{2-body}_{[140,\infty)}$ , (d) $SR^{2-body}_{[160,\infty)}$ , (e) $SR^{2-body}_{[180,\infty)}$ , (f) $SR^{2-body}_{[200,\infty)}$ , (g) $SR^{2-body}_{[220,\infty)}$ for a simplified model assuming $t \tilde{t} + \phi$ pair production.

Two-body selection acceptance (a) SR-DF$^{2-body}_{[110,120)}$, (b) SR-DF1$^{2-body}_{[120,140)}$, (c) SR-DF2$^{2-body}_{[140,160)}$, (d) SR-DF3$^{2-body}_{[160,180)}$, (e) SR-DF4$^{2-body}_{[180,220)}$, (f) SR-DF5$^{2-body}_{[220,\infty)}$ for a simplified model assuming $t \tilde{t} + a$ pair production.

Two-body selection acceptance (a) SR-DF$^{2-body}_{[110,120)}$, (b) SR-DF1$^{2-body}_{[120,140)}$, (c) SR-DF2$^{2-body}_{[140,160)}$, (d) SR-DF3$^{2-body}_{[160,180)}$, (e) SR-DF4$^{2-body}_{[180,220)}$, (f) SR-DF5$^{2-body}_{[220,\infty)}$ for a simplified model assuming $t \tilde{t} + a$ pair production.

Two-body selection acceptance (a) SR-DF$^{2-body}_{[110,120)}$, (b) SR-DF1$^{2-body}_{[120,140)}$, (c) SR-DF2$^{2-body}_{[140,160)}$, (d) SR-DF3$^{2-body}_{[160,180)}$, (e) SR-DF4$^{2-body}_{[180,220)}$, (f) SR-DF5$^{2-body}_{[220,\infty)}$ for a simplified model assuming $t \tilde{t} + a$ pair production.

Two-body selection acceptance (a) SR-DF$^{2-body}_{[110,120)}$, (b) SR-DF1$^{2-body}_{[120,140)}$, (c) SR-DF2$^{2-body}_{[140,160)}$, (d) SR-DF3$^{2-body}_{[160,180)}$, (e) SR-DF4$^{2-body}_{[180,220)}$, (f) SR-DF5$^{2-body}_{[220,\infty)}$ for a simplified model assuming $t \tilde{t} + a$ pair production.

Two-body selection acceptance (a) SR-DF$^{2-body}_{[110,120)}$, (b) SR-DF1$^{2-body}_{[120,140)}$, (c) SR-DF2$^{2-body}_{[140,160)}$, (d) SR-DF3$^{2-body}_{[160,180)}$, (e) SR-DF4$^{2-body}_{[180,220)}$, (f) SR-DF5$^{2-body}_{[220,\infty)}$ for a simplified model assuming $t \tilde{t} + a$ pair production.

Two-body selection acceptance (a) SR-DF$^{2-body}_{[110,120)}$, (b) SR-DF1$^{2-body}_{[120,140)}$, (c) SR-DF2$^{2-body}_{[140,160)}$, (d) SR-DF3$^{2-body}_{[160,180)}$, (e) SR-DF4$^{2-body}_{[180,220)}$, (f) SR-DF5$^{2-body}_{[220,\infty)}$ for a simplified model assuming $t \tilde{t} + a$ pair production.

Two-body selection acceptance (a) SR-SF$^{2-body}_{[110,120)}$, (b) SR-SF1$^{2-body}_{[120,140)}$, (c) SR-SF2$^{2-body}_{[140,160)}$, (d) SR-SF3$^{2-body}_{[160,180)}$, (e) SR-SF4$^{2-body}_{[180,220)}$, (f) SR-SF5$^{2-body}_{[220,\infty)}$ for a simplified model assuming t \tilde{t} + a$ pair production.

Two-body selection acceptance (a) SR-SF$^{2-body}_{[110,120)}$, (b) SR-SF1$^{2-body}_{[120,140)}$, (c) SR-SF2$^{2-body}_{[140,160)}$, (d) SR-SF3$^{2-body}_{[160,180)}$, (e) SR-SF4$^{2-body}_{[180,220)}$, (f) SR-SF5$^{2-body}_{[220,\infty)}$ for a simplified model assuming $t \tilde{t} + a$ pair production.

Two-body selection acceptance (a) SR-SF$^{2-body}_{[110,120)}$, (b) SR-SF1$^{2-body}_{[120,140)}$, (c) SR-SF2$^{2-body}_{[140,160)}$, (d) SR-SF3$^{2-body}_{[160,180)}$, (e) SR-SF4$^{2-body}_{[180,220)}$, (f) SR-SF5$^{2-body}_{[220,\infty)}$ for a simplified model assuming $t \tilde{t} + a$ pair production.

Two-body selection acceptance (a) SR-SF$^{2-body}_{[110,120)}$, (b) SR-SF1$^{2-body}_{[120,140)}$, (c) SR-SF2$^{2-body}_{[140,160)}$, (d) SR-SF3$^{2-body}_{[160,180)}$, (e) SR-SF4$^{2-body}_{[180,220)}$, (f) SR-SF5$^{2-body}_{[220,\infty)}$ for a simplified model assuming $t \tilde{t} + a$ pair production.

Two-body selection acceptance (a) SR-SF$^{2-body}_{[110,120)}$, (b) SR-SF1$^{2-body}_{[120,140)}$, (c) SR-SF2$^{2-body}_{[140,160)}$, (d) SR-SF3$^{2-body}_{[160,180)}$, (e) SR-SF4$^{2-body}_{[180,220)}$, (f) SR-SF5$^{2-body}_{[220,\infty)}$ for a simplified model assuming $t \tilde{t} + a$ pair production.

Two-body selection acceptance (a) SR-SF$^{2-body}_{[110,120)}$, (b) SR-SF1$^{2-body}_{[120,140)}$, (c) SR-SF2$^{2-body}_{[140,160)}$, (d) SR-SF3$^{2-body}_{[160,180)}$, (e) SR-SF4$^{2-body}_{[180,220)}$, (f) SR-SF5$^{2-body}_{[220,\infty)}$ for a simplified model assuming $t \tilde{t} + a$ pair production.

Two-body selection acceptance (a) $SR^{2-body}_{[110,\infty)}$ , (b) $SR^{2-body}_{[120,\infty)}$ , (c) $SR^{2-body}_{[140,\infty)}$ , (d) $SR^{2-body}_{[160,\infty)}$ , (e) $SR^{2-body}_{[180,\infty)}$ , (f) $SR^{2-body}_{[200,\infty)}$ , (g) $SR^{2-body}_{[220,\infty)}$ for a simplified model assuming $ \tilde{t}_1$ pair production.

Two-body selection acceptance (a) $SR^{2-body}_{[110,\infty)}$ , (b) $SR^{2-body}_{[120,\infty)}$ , (c) $SR^{2-body}_{[140,\infty)}$ , (d) $SR^{2-body}_{[160,\infty)}$ , (e) $SR^{2-body}_{[180,\infty)}$ , (f) $SR^{2-body}_{[200,\infty)}$ , (g) $SR^{2-body}_{[220,\infty)}$ for a simplified model assuming $ \tilde{t}_1$ pair production.

Two-body selection acceptance (a) $SR^{2-body}_{[110,\infty)}$ , (b) $SR^{2-body}_{[120,\infty)}$ , (c) $SR^{2-body}_{[140,\infty)}$ , (d) $SR^{2-body}_{[160,\infty)}$ , (e) $SR^{2-body}_{[180,\infty)}$ , (f) $SR^{2-body}_{[200,\infty)}$ , (g) $SR^{2-body}_{[220,\infty)}$ for a simplified model assuming $ \tilde{t}_1$ pair production.

Two-body selection acceptance (a) $SR^{2-body}_{[110,\infty)}$ , (b) $SR^{2-body}_{[120,\infty)}$ , (c) $SR^{2-body}_{[140,\infty)}$ , (d) $SR^{2-body}_{[160,\infty)}$ , (e) $SR^{2-body}_{[180,\infty)}$ , (f) $SR^{2-body}_{[200,\infty)}$ , (g) $SR^{2-body}_{[220,\infty)}$ for a simplified model assuming $ \tilde{t}_1$ pair production.

Two-body selection acceptance (a) $SR^{2-body}_{[110,\infty)}$ , (b) $SR^{2-body}_{[120,\infty)}$ , (c) $SR^{2-body}_{[140,\infty)}$ , (d) $SR^{2-body}_{[160,\infty)}$ , (e) $SR^{2-body}_{[180,\infty)}$ , (f) $SR^{2-body}_{[200,\infty)}$ , (g) $SR^{2-body}_{[220,\infty)}$ for a simplified model assuming $ \tilde{t}_1$ pair production.

Two-body selection acceptance (a) $SR^{2-body}_{[110,\infty)}$ , (b) $SR^{2-body}_{[120,\infty)}$ , (c) $SR^{2-body}_{[140,\infty)}$ , (d) $SR^{2-body}_{[160,\infty)}$ , (e) $SR^{2-body}_{[180,\infty)}$ , (f) $SR^{2-body}_{[200,\infty)}$ , (g) $SR^{2-body}_{[220,\infty)}$ for a simplified model assuming $ \tilde{t}_1$ pair production.

Two-body selection acceptance (a) $SR^{2-body}_{[110,\infty)}$ , (b) $SR^{2-body}_{[120,\infty)}$ , (c) $SR^{2-body}_{[140,\infty)}$ , (d) $SR^{2-body}_{[160,\infty)}$ , (e) $SR^{2-body}_{[180,\infty)}$ , (f) $SR^{2-body}_{[200,\infty)}$ , (g) $SR^{2-body}_{[220,\infty)}$ for a simplified model assuming $ \tilde{t}_1$ pair production.

Three-body selection acceptance (a) SR-DF$^{3-body}_{t}$, (b) SR-SF$^{3-body}_{t}$, (c) SR-DF$^{3-body}_{W}$, (d) SR-SF$^{3-body}_{W}$ for a simplified model assuming $ \tilde{t}_1$ pair production.

Three-body selection acceptance (a) SR-DF$^{3-body}_{t}$, (b) SR-SF$^{3-body}_{t}$, (c) SR-DF$^{3-body}_{W}$, (d) SR-SF$^{3-body}_{W}$ for a simplified model assuming $ \tilde{t}_1$ pair production.

Three-body selection acceptance (a) SR-DF$^{3-body}_{t}$, (b) SR-SF$^{3-body}_{t}$, (c) SR-DF$^{3-body}_{W}$, (d) SR-SF$^{3-body}_{W}$ for a simplified model assuming $ \tilde{t}_1$ pair production.

Three-body selection acceptance (a) SR-DF$^{3-body}_{t}$, (b) SR-SF$^{3-body}_{t}$, (c) SR-DF$^{3-body}_{W}$, (d) SR-SF$^{3-body}_{W}$ for a simplified model assuming $ \tilde{t}_1$ pair production.

Four-body selection acceptance (a) SR$^{4-body}_{Small \Delta m}$ , (b) $SR^{4-body}_{Large \Delta m}$ for a simplified model assuming $\tilde{t}_1$ pair production.

Four-body selection acceptance (a) SR$^{4-body}_{Small \Delta m}$ , (b) $SR^{4-body}_{Large \Delta m}$ for a simplified model assuming $\tilde{t}_1$ pair production.

Two-body selection The numbers indicate the observed upper limits on the signal strenght for (a) a simplified model assuming $\tilde{t}_1$ pair production, (b) for $t\tilde{t} + a $ pseudoscalar models, (c) for $t\tilde{t} + \phi $ scalar models. In Figure (a), the red line corresponds to the observed limit.

Two-body selection The numbers indicate the observed upper limits on the signal strenght for (a) a simplified model assuming $\tilde{t}_1$ pair production, (b) for $t\tilde{t} + a $ pseudoscalar models, (c) for $t\tilde{t} + \phi $ scalar models. In Figure (a), the red line corresponds to the observed limit.

Two-body selection The numbers indicate the observed upper limits on the signal strenght for (a) a simplified model assuming $\tilde{t}_1$ pair production, (b) for $t\tilde{t} + a $ pseudoscalar models, (c) for $t\tilde{t} + \phi $ scalar models. In Figure (a), the red line corresponds to the observed limit.

Three-body selection The numbers indicate the upper limits on the signal strenght for a simplified model assuming $\tilde{t}_1$ pair production. For comparison, the red line corresponds to the observed limit.

Four-body selection The numbers indicate the upper limits on the signal strenght for a simplified model assuming $\tilde{t}_1$ pair production. For comparison, the red line corresponds to the observed limit.

Two-body selection The numbers indicate the upper limits on the signal cross-section for (a) a simplified model assuming $\tilde{t}_1$ pair production, (b) for $t\tilde{t} + a $ pseudoscalar models, (c) for $t\tilde{t} + \phi $ scalar models. In Figure (a), the red line corresponds to the observed limit.

Two-body selection The numbers indicate the upper limits on the signal cross-section for (a) a simplified model assuming $\tilde{t}_1$ pair production, (b) for $t\tilde{t} + a $ pseudoscalar models, (c) for $t\tilde{t} + \phi $ scalar models. In Figure (a), the red line corresponds to the observed limit.

Two-body selection The numbers indicate the upper limits on the signal cross-section for (a) a simplified model assuming $\tilde{t}_1$ pair production, (b) for $t\tilde{t} + a $ pseudoscalar models, (c) for $t\tilde{t} + \phi $ scalar models. In Figure (a), the red line corresponds to the observed limit.

Three-body selection The numbers indicate the upper limits on the signal cross-section for a simplified model assuming $\tilde{t}_1$ pair production. For comparison, the red line corresponds to the observed limit.

Four-body selection The numbers indicate the upper limits on the signal cross-section for a simplified model assuming $\tilde{t}_1$ pair production. For comparison, the red line corresponds to the observed limit.

Two-body selection. Background fit results for the $inclusive$ SRs. The Others category contains the contributions from $VVV$, $t\bar{t} t$, $t\bar{t}t\bar{t}$, $t\bar{t} W$, $t\bar{t} WW$, $t\bar{t} WZ$, $t\bar{t} H$, and $tZ$. Combined statistical and systematic uncertainties are given. Note that the individual uncertainties can be correlated, and do not necessarily add up quadratically to the total background uncertainty.

Cut flow for the simplified signal model $\tilde{t}_1 \rightarrow t^{(*)}\tilde{\chi}^0_1$ with $m(\tilde{t}_1)=600~ GeV$ and $m(\tilde{\chi}^0_1)=400~ GeV$ in the SRs for the two-body selection. The number of events is normalized to the cross-section and to an integrated luminosity of $139~fb^{-1}$.

Cut flow for the scalar signal model $t\bar{t} + \phi $ with $m(\phi)=150~ GeV$ and $m(\chi)=1~ GeV$ in the SRs for the two-body selection. The number of events is normalized to the cross-section and to an integrated luminosity of $139~fb^{-1}$.

Cut flow for the pseudoscalar signal model $t\bar{t} + a $ with $m(a)=150~ GeV$ and $m(\chi)=1~ GeV$ in the SRs for the two-body selection. The number of events is normalized to the cross-section and to an integrated luminosity of $139~fb^{-1}$.

Cut flow for the simplified signal model $\tilde{t}_1 \rightarrow bW\tilde{\chi}^0_1$ with $m(\tilde{t}_1)=550~ GeV$ and $m(\tilde{\chi}^0_1)=385~ GeV$ in the SRs for the three-body selection. The number of events is normalized to the cross-section and to an integrated luminosity of $139~fb^{-1}$.

Cut flow for the simplified signal model $\tilde{t}_1 \rightarrow bW\tilde{\chi}^0_1$ with $m(\tilde{t}_1)=550~ GeV$ and $m(\tilde{\chi}^0_1)=400~ GeV$ in the SRs for the three-body selection. The number of events is normalized to the cross-section and to an integrated luminosity of $139~fb^{-1}$.

Cut flow for the simplified signal model $\tilde{t}_1 \rightarrow bW\tilde{\chi}^0_1$ with $m(\tilde{t}_1)=550~ GeV$ and $m(\tilde{\chi}^0_1)=430~ GeV$ in the SRs for the three-body selection. The number of events is normalized to the cross-section and to an integrated luminosity of $139~fb^{-1}$.

Cut flow for the simplified signal model $\tilde{t}_1 \rightarrow bW\tilde{\chi}^0_1$ with $m(\tilde{t}_1)=550~ GeV$ and $m(\tilde{\chi}^0_1)=460~ GeV$ in the SRs for the three-body selection. The number of events is normalized to the cross-section and to an integrated luminosity of $139~fb^{-1}$.

Cut flow for the simplified signal model $\tilde{t}_1 \rightarrow b l \nu \tilde{\chi}^0_1$ with $m(\tilde{t}_1)=400~ GeV$ and $m(\tilde{\chi}^0_1)=380~ GeV$ in the SRs for the four-body selection. The number of events is normalized to the cross-section and to an integrated luminosity of $139~fb^{-1}$.

Cut flow for the simplified signal model $\tilde{t}_1 \rightarrow b l \nu \tilde{\chi}^0_1$ with $m(\tilde{t}_1)=460~ GeV$ and $m(\tilde{\chi}^0_1)=415~ GeV$ in the SRs for the four-body selection. The number of events is normalized to the cross-section and to an integrated luminosity of $139~fb^{-1}$.

Cut flow for the simplified signal model $\tilde{t}_1 \rightarrow b l \nu \tilde{\chi}^0_1$ with $m(\tilde{t}_1)=400~ GeV$ and $m(\tilde{\chi}^0_1)=320~ GeV$ in the SRs for the four-body selection. The number of events is normalized to the cross-section and to an integrated luminosity of $139~fb^{-1}$.

The ratios (&mu;) of the 95&#37; C.L. upper limits on the combined s&rarr; W^&plusmn;W^&#8723; and s&rarr; ZZ cross section to simplified model expectations for the m_Z'=1 TeV scenario, for various m_s hypotheses. The observed limits (solid line) are consistent with the expectation under the SM-only hypothesis (dashed line) within uncertainties (filled band), except for a small excess for m_s=160 GeV, discussed in the text.

The ratios (&mu;) of the 95&#37; C.L. upper limits on the combined s&rarr; W^&plusmn;W^&#8723; and s&rarr; ZZ cross section to simplified model expectations for the m_Z'=1.7 TeV scenario, for various m_s hypotheses. The observed limits (solid line) are consistent with the expectation under the SM-only hypothesis (dashed line) within uncertainties (filled band), except for a small excess for m_s=160 GeV, discussed in the text.

Observed upper limits at 95&#37; C.L. on &sigma;(pp &rarr; s &chi;&chi;) × B(s&rarr; VV) for m_Z'=0.5 TeV signal points. The expected limits, varied up and down by one and two standard deviations, are shown as green and yellow bands, respectively. The observed and expected limits are compared to the theoretical LO cross section for the &sigma;(pp &rarr; s &chi;&chi;) × B(s&rarr; VV) process for m_Z'=0.5 TeV, shown in dashed blue.

Observed upper limits at 95&#37; C.L. on &sigma;(pp &rarr; s &chi;&chi;) × B(s&rarr; VV) for m_Z'=1 TeV signal points. The expected limits, varied up and down by one and two standard deviations, are shown as green and yellow bands, respectively. The observed and expected limits are compared to the theoretical LO cross section for the &sigma;(pp &rarr; s &chi;&chi;) × B(s&rarr; VV) process for m_Z'=1 TeV, shown in dashed blue.

Observed upper limits at 95&#37; C.L. on &sigma;(pp &rarr; s &chi;&chi;) × B(s&rarr; VV) for m_Z'=1.7 TeV signal points. The expected limits, varied up and down by one and two standard deviations, are shown as green and yellow bands, respectively. The observed and expected limits are compared to the theoretical LO cross section for the &sigma;(pp &rarr; s &chi;&chi;) × B(s&rarr; VV) process for m_Z'=1.7 TeV, shown in dashed blue.

SM background post-fit yields stacked in each SR and CR category and E_T^miss bin and data overlaid with the maximum likelihood estimators set to the conditional values of the combined signal and control region fit. The hatched uncertainty band shown includes simulation statistics uncertainties, experimental systematic uncertainties, and V+jets theory modelling systematic uncertainties. Pre-fit uncertainties cover differences between the data and pre-fit background prediction.

Cumulative efficiencies for the merged category for signal samples with m_s=160 GeV (a), m_s=235 GeV (b) and m_s=310 GeV (c), each with m_Z'=1 TeV. The dark Higgs candidate selection includes stringent jet substructure requirements and typically at most one candidate is present in signal events. Here, &Delta; &phi;_{jets_1,2,3 E_T^miss} is the smallest azimuthal angle between the E_T^miss and any of the three highest-p_T (leading) small-R jets.

Cumulative efficiencies for the intermediate category for signal samples with m_s=160 GeV (a), m_s=235 GeV (b) and m_s=310 GeV (c), each with m_Z'=1 TeV. The TAR+Comb algorithm reconstructs the dark Higgs candidate from a TAR jet with m^TAR&gt;60 GeV that is supplemented by up to two additional small-R jets within &Delta;R_cone=2.5 of the TAR jet. Here, &Delta; &phi;_{jets_1,2,3 E_T^miss} is the smallest azimuthal angle between the E_T^miss and any of the three highest-p_T (leading) small-R jets. For details see text.

The product of acceptance and efficiency (A &times; &varepsilon;), defined as the number of signal events satisfying the full set of selection criteria in the merged or intermediate signal regions, divided by the total number of generated signal events, for the s(W^&plusmn;W^&#8723;) dark Higgs signal points with dark Higgs boson mass m_s and Z' boson mass m_Z'.

The product of acceptance and efficiency (A &times; &varepsilon;), defined as the number of signal events satisfying the full set of selection criteria in the merged or intermediate signal regions, divided by the total number of generated signal events, for the s(ZZ) dark Higgs signal points with dark Higgs boson mass m_s and Z' boson mass m_Z'.

Observed excluded cross-section at 95% CL for the WZ-mediated simplified model of wino $\tilde{\chi}^{\pm}_{1}/\tilde{\chi}^{0}_{2}$ production from Fig 8(aux).

Expected exclusion limits at 95% CL for the WZ-mediated simplified model of wino $\tilde{\chi}^{\pm}_{1}/\tilde{\chi}^{0}_{2}$ production.

Observed exclusion limits at 95% CL for the Wh-mediated simplified model of wino $\tilde{\chi}^{\pm}_{1}/\tilde{\chi}^{0}_{2}$ production from from Fig 13(a) and from Fig 7 and Fig 10(aux).

Observed excluded cross-section at 95% CL for the Wh-mediated simplified model of wino $\tilde{\chi}^{\pm}_{1}/\tilde{\chi}^{0}_{2}$ production from Fig 7(aux) and Fig 10(aux).

positive one $\sigma$ observed exclusion limits at 95% CL for the Wh-mediated simplified model of wino $\tilde{\chi}^{\pm}_{1}/\tilde{\chi}^{0}_{2}$ production from from Fig 13(a) and from Fig 7 and Fig 10(aux).

negative one $\sigma$ observed exclusion limits at 95% CL for the Wh-mediated simplified model of wino $\tilde{\chi}^{\pm}_{1}/\tilde{\chi}^{0}_{2}$ production from from Fig 13(a) and from Fig 7 and Fig 10(aux).

Expected exclusion limits at 95% CL for the Wh-mediated simplified model of wino $\tilde{\chi}^{\pm}_{1}/\tilde{\chi}^{0}_{2}$ production.

Expected exclusion limits at 95% CL for the Wh-mediated simplified model of wino $\tilde{\chi}^{\pm}_{1}/\tilde{\chi}^{0}_{2}$ production.

Expected exclusion limits at 95% CL for the Wh-mediated simplified model of wino $\tilde{\chi}^{\pm}_{1}/\tilde{\chi}^{0}_{2}$ production.

Number of signal events expected for 139 fb$^{-1}$ at different stages of the event selection for the signal region $SR^{bRPV}_{2l-SS}$. in a susy scenario where $\tilde{\chi}^{0} _{1} and \tilde{\chi}^{0} _{2}$ are produced in pairs and decay to all possible allowed bRPV decays. The masses of the superpartners involved in the process are set to $m(\tilde{\chi}^{0} _{1}/\tilde{\chi}^{0} _{2})$ = 200 GeV, tan$\beta$=5. Only statistical uncertainties are shown.

Number of signal events expected for 139 fb$^{-1}$ at different stages of the event selection for the signal region $SR^{bRPV}_{3l}$. in a susy scenario where $\tilde{\chi}^{0} _{1} and \tilde{\chi}^{0} _{2}$ are produced in pairs and decay to all possible allowed bRPV decays. The masses of the superpartners involved in the process are set to $m(\tilde{\chi}^{0} _{1}/\tilde{\chi}^{0} _{2})$ = 200 GeV, tan$\beta$=5. Only statistical uncertainties are shown.

Number of signal events expected for 139 fb$^{-1}$ at different stages of the event selection for the signal region $SR^{WZ}_{high-m_{T2}}$. The wino-like doublet pair ($\tilde{\chi}^{\pm} _{1} and \tilde{\chi}^{0} _{2}$) were produced and then decays into $bino-like \tilde{\chi}^{0} _{1}$ which is the lightest SUSY particle (LSP) accompanied by mass on-shell or mass off-shell W and Z bosons. The masses of the superpartners involved in the process are set to $m(\tilde{\chi}^{\pm} _{1}/\tilde{\chi}^{0} _{2})$ = 150 GeV, $m(\tilde{\chi}^{0} _{1})$ = 50 GeV. Only statistical uncertainties are shown.

Number of signal events expected for 139 fb$^{-1}$ at different stages of the event selection for the signal region $SR^{WZ}_{low-m_{T2}}$. The wino-like doublet pair ($\tilde{\chi}^{\pm} _{1} and \tilde{\chi}^{0} _{2}$) were produced and then decays into $bino-like \tilde{\chi}^{0} _{1}$ which is the lightest SUSY particle (LSP) accompanied by mass on-shell or mass off-shell W and Z bosons. The masses of the superpartners involved in the process are set to $m(\tilde{\chi}^{\pm} _{1}/\tilde{\chi}^{0} _{2})$ = 150 GeV, $m(\tilde{\chi}^{0} _{1})$ = 50 GeV. Only statistical uncertainties are shown.

Number of signal events expected for 139 fb$^{-1}$ at different stages of the event selection for the low mass $SR^{RPV}_{2l1b}$, where the $\tilde{\chi}^{0} _{1} and \tilde{\chi}^{0} _{2}$ are directly produced and undergoes prompt RPV decays. The masses of the superpartners involved in the process are set to $m(\tilde{\chi}^{0} _{1}/\tilde{\chi}^{0} _{2})$ = 200 GeV. Only statistical uncertainties are shown.

Number of signal events expected for 139 fb$^{-1}$ at different stages of the event selection for the medium mass $SR^{RPV}_{2l1b}$, where the $\tilde{\chi}^{0} _{1} and \tilde{\chi}^{0} _{2}$ are directly produced and undergoes prompt RPV decays. The masses of the superpartners involved in the process are set to $m(\tilde{\chi}^{0} _{1}/\tilde{\chi}^{0} _{2})$ = 200 GeV. Only statistical uncertainties are shown.

Number of signal events expected for 139 fb$^{-1}$ at different stages of the event selection for the low mass $SR^{RPV}_{2l2b}$, where the $\tilde{\chi}^{0} _{1} and \tilde{\chi}^{0} _{2}$ are directly produced and undergoes prompt RPV decays. The masses of the superpartners involved in the process are set to $m(\tilde{\chi}^{0} _{1}/\tilde{\chi}^{0} _{2})$ = 200 GeV. Only statistical uncertainties are shown.

Number of signal events expected for 139 fb$^{-1}$ at different stages of the event selection for the medium mass $SR^{RPV}_{2l2b}$, where the $\tilde{\chi}^{0} _{1} and \tilde{\chi}^{0} _{2}$ are directly produced and undergoes prompt RPV decays. The masses of the superpartners involved in the process are set to $m(\tilde{\chi}^{0} _{1}/\tilde{\chi}^{0} _{2})$ = 200 GeV. Only statistical uncertainties are shown.

Number of signal events expected for 139 fb$^{-1}$ at different stages of the event selection for the high mass $SR^{RPV}_{2l2b}$, where the $\tilde{\chi}^{0} _{1} and \tilde{\chi}^{0} _{2}$ are directly produced and undergoes prompt RPV decays. The masses of the superpartners involved in the process are set to $m(\tilde{\chi}^{0} _{1}/\tilde{\chi}^{0} _{2})$ = 200 GeV. Only statistical uncertainties are shown.

Number of signal events expected for 139 fb$^{-1}$ at different stages of the event selection for the low mass $SR^{RPV}_{2l3b}$, where the $\tilde{\chi}^{0} _{1} and \tilde{\chi}^{0} _{2}$ are directly produced and undergoes prompt RPV decays. The masses of the superpartners involved in the process are set to $m(\tilde{\chi}^{0} _{1}/\tilde{\chi}^{0} _{2})$ = 200 GeV. Only statistical uncertainties are shown.

Number of signal events expected for 139 fb$^{-1}$ at different stages of the event selection for the medium mass $SR^{RPV}_{2l3b}$, where the $\tilde{\chi}^{0} _{1} and \tilde{\chi}^{0} _{2}$ are directly produced and undergoes prompt RPV decays. The masses of the superpartners involved in the process are set to $m(\tilde{\chi}^{0} _{1}/\tilde{\chi}^{0} _{2})$ = 200 GeV. Only statistical uncertainties are shown.

Number of signal events expected for 139 fb$^{-1}$ at different stages of the event selection for the high mass $SR^{RPV}_{2l3b}$, where the $\tilde{\chi}^{0} _{1} and \tilde{\chi}^{0} _{2}$ are directly produced and undergoes prompt RPV decays. The masses of the superpartners involved in the process are set to $m(\tilde{\chi}^{0} _{1}/\tilde{\chi}^{0} _{2})$ = 200 GeV. Only statistical uncertainties are shown.

Number of signal events expected for 139 fb$^{-1}$ at different stages of the event selection for the $SR^{Wh}_{low-m_{T2} }$. The wino-like doublet pair ($\tilde{\chi}^{\pm} _{1} and \tilde{\chi}^{0} _{2}$) were produced and then decays into $bino-like \tilde{\chi}^{0} _{1}$ which is the lightest SUSY particle (LSP) accompanied by mass on-shell or mass off-shell W and Higgs bosons. The masses of the superpartners involved in the process are set to $m(\tilde{\chi}^{\pm} _{1}/\tilde{\chi}^{0} _{2})$ = 300 GeV, $m(\tilde{\chi}^{0} _{1})$ = 100 GeV. Only statistical uncertainties are shown.

Number of signal events expected for 139 fb$^{-1}$ at different stages of the event selection for the $SR^{Wh}_{high-m_{T2} }$. The wino-like doublet pair ($\tilde{\chi}^{\pm} _{1} and \tilde{\chi}^{0} _{2}$) were produced and then decays into $bino-like \tilde{\chi}^{0} _{1}$ which is the lightest SUSY particle (LSP) accompanied by mass on-shell or mass off-shell W and Higgs bosons. The masses of the superpartners involved in the process are set to $m(\tilde{\chi}^{\pm} _{1}/\tilde{\chi}^{0} _{2})$ = 300 GeV, $m(\tilde{\chi}^{0} _{1})$ = 100 GeV. Only statistical uncertainties are shown.

Signal Hepdataeptance for $SR^{bRPV}_{2l-SS}$ signal region from Fig 13(a)(aux) in a SUSY scenario where $\tilde{\chi}^{0} _{1} and \tilde{\chi}^{0} _{2}$ are produced in pairs and decay to all possible allowed bRPV decays.

Signal Hepdataeptance for $SR^{bRPV}_{3l}$ signal region from Fig 13(b)(aux) in a SUSY scenario where $\tilde{\chi}^{0} _{1} and \tilde{\chi}^{0} _{2}$ are produced in pairs and decay to all possible allowed bRPV decays.

Signal acceptance for $SR^{WZ}_{high-m_{T2}}$ in a SUSY scenario where the wino-like doublet pair ($\tilde{\chi}^{\pm} _{1} and \tilde{\chi}^{0} _{2}$) were produced and then decays into $bino-like \tilde{\chi}^{0} _{1}$ which is the lightest SUSY particle (LSP) accompanied by mass on-shell or mass off-shell W and Z bosons.

Signal acceptance for $SR^{WZ}_{low-m_{T2}}$ in a SUSY scenario where the wino-like doublet pair ($\tilde{\chi}^{\pm} _{1} and \tilde{\chi}^{0} _{2}$) were produced and then decays into $bino-like \tilde{\chi}^{0} _{1}$ which is the lightest SUSY particle (LSP) accompanied by mass on-shell or mass off-shell W and Z bosons.

Signal acceptance for $SR^{RPV}_{2l1b}-L$ signal region in a SUSY scenario where the $\tilde{\chi}^{0} _{1} and \tilde{\chi}^{0} _{2}$ are directly produced and undergoes prompt RPV decays.

Signal acceptance for $SR^{RPV}_{2l1b}-M$ signal region in a SUSY scenario where the $\tilde{\chi}^{0} _{1} and \tilde{\chi}^{0} _{2}$ are directly produced and undergoes prompt RPV decays.

Signal acceptance for $SR^{RPV}_{2l2b}-L$ signal region in a SUSY scenario where the $\tilde{\chi}^{0} _{1} and \tilde{\chi}^{0} _{2}$ are directly produced and undergoes prompt RPV decays.

Signal acceptance for $SR^{RPV}_{2l2b}-M$ signal region in a SUSY scenario where the $\tilde{\chi}^{0} _{1} and \tilde{\chi}^{0} _{2}$ are directly produced and undergoes prompt RPV decays.

Signal acceptance for $SR^{RPV}_{2l2b}-H$ signal region in a SUSY scenario where the $\tilde{\chi}^{0} _{1} and \tilde{\chi}^{0} _{2}$ are directly produced and undergoes prompt RPV decays.

Signal acceptance for $SR^{RPV}_{2l3b}-L$ signal region in a SUSY scenario where the $\tilde{\chi}^{0} _{1} and \tilde{\chi}^{0} _{2}$ are directly produced and undergoes prompt RPV decays.

Signal acceptance for $SR^{RPV}_{2l3b}-M$ signal region in a SUSY scenario where the $\tilde{\chi}^{0} _{1} and \tilde{\chi}^{0} _{2}$ are directly produced and undergoes prompt RPV decays.

Signal acceptance for $SR^{RPV}_{2l3b}-H$ signal region in a SUSY scenario where the $\tilde{\chi}^{0} _{1} and \tilde{\chi}^{0} _{2}$ are directly produced and undergoes prompt RPV decays.

Signal efficiency for $SR^{bRPV}_{2l-SS}$ signal region in a SUSY scenario where $\tilde{\chi}^{0} _{1} and \tilde{\chi}^{0} _{2}$ are produced in pairs and decay to all possible allowed bRPV decays.

Signal efficiency for $SR^{bRPV}_{3l}$ signal region in a SUSY scenario where $\tilde{\chi}^{0} _{1} and \tilde{\chi}^{0} _{2}$ are produced in pairs and decay to all possible allowed bRPV decays.

Signal efficiency for $SR^{WZ}_{high-m_{T2}}$ in a SUSY scenario where the wino-like doublet pair ($\tilde{\chi}^{\pm} _{1} and \tilde{\chi}^{0} _{2}$) were produced and then decays into $bino-like \tilde{\chi}^{0} _{1}$ which is the lightest SUSY particle (LSP) accompanied by mass on-shell or mass off-shell W and Z bosons.

Signal efficiency for $SR^{WZ}_{low-m_{T2}}$ in a SUSY scenario where the wino-like doublet pair ($\tilde{\chi}^{\pm} _{1} and \tilde{\chi}^{0} _{2}$) were produced and then decays into $bino-like \tilde{\chi}^{0} _{1}$ which is the lightest SUSY particle (LSP) accompanied by mass on-shell or mass off-shell W and Z bosons.

Signal efficiency for $SR^{RPV}_{2l1b}-L$ signal region in a SUSY scenario where the $\tilde{\chi}^{0} _{1} and \tilde{\chi}^{0} _{2}$ are directly produced and undergoes prompt RPV decays.

Signal efficiency for $SR^{RPV}_{2l1b}-M$ signal region in a SUSY scenario where the $\tilde{\chi}^{0} _{1} and \tilde{\chi}^{0} _{2}$ are directly produced and undergoes prompt RPV decays.

Signal efficiency for $SR^{RPV}_{2l2b}-L$ signal region in a SUSY scenario where the $\tilde{\chi}^{0} _{1} and \tilde{\chi}^{0} _{2}$ are directly produced and undergoes prompt RPV decays.

Signal efficiency for $SR^{RPV}_{2l2b}-M$ signal region in a SUSY scenario where the $\tilde{\chi}^{0} _{1} and \tilde{\chi}^{0} _{2}$ are directly produced and undergoes prompt RPV decays.

Signal efficiency for $SR^{RPV}_{2l2b}-H$ signal region in a SUSY scenario where the $\tilde{\chi}^{0} _{1} and \tilde{\chi}^{0} _{2}$ are directly produced and undergoes prompt RPV decays.

Signal efficiency for $SR^{RPV}_{2l3b}-L$ signal region in a SUSY scenario where the $\tilde{\chi}^{0} _{1} and \tilde{\chi}^{0} _{2}$ are directly produced and undergoes prompt RPV decays.

Signal efficiency for $SR^{RPV}_{2l3b}-M$ signal region in a SUSY scenario where the $\tilde{\chi}^{0} _{1} and \tilde{\chi}^{0} _{2}$ are directly produced and undergoes prompt RPV decays.

Signal efficiency for $SR^{RPV}_{2l3b}-H$ signal region in a SUSY scenario where the $\tilde{\chi}^{0} _{1} and \tilde{\chi}^{0} _{2}$ are directly produced and undergoes prompt RPV decays.

Signal acceptance for $SR^{Wh}_{high-m_{T2} }$ signal region from Fig 11(a)(aux) in a SUSY scenario where direct production of a lightest $\tilde{\chi}^{\pm} _{1} and \tilde{\chi}^{0} _{2}$ , decay with 100% branching ratio to a final state with a same sign light lepton (e or $\mu$) pair and two lightest neutralino1, via the on-shell emission of SM W and Higgs bosons,

Signal acceptance for $SR^{Wh}_{low-m_{T2} }$ signal region from Fig 11(b)(aux) in a SUSY scenario where direct production of a lightest $\tilde{\chi}^{\pm} _{1} and \tilde{\chi}^{0} _{2}$ , decay with 100% branching ratio to a final state with a same sign light lepton (e or $\mu$) pair and two lightest neutralino1, via the on-shell emission of SM W and Higgs bosons,

Signal efficiency for $SR^{Wh}_{high-m_{T2} }$ signal region from Fig 15(a)(aux) in a SUSY scenario where direct production of a lightest $\tilde{\chi}^{\pm} _{1} and \tilde{\chi}^{0} _{2}$ , decay with 100% branching ratio to a final state with a same sign light lepton (e or $\mu$) pair and two lightest neutralino1, via the on-shell emission of SM W and Higgs bosons,

Signal efficiency for $SR^{Wh}_{low-m_{T2} }$ signal region from Fig 15(b)(aux) in a SUSY scenario where direct production of a lightest $\tilde{\chi}^{\pm} _{1} and \tilde{\chi}^{0} _{2}$ , decay with 100% branching ratio to a final state with a same sign light lepton (e or $\mu$) pair and two lightest neutralino1, via the on-shell emission of SM W and Higgs bosons,

Observed 95% X-section upper limits as a function of higgsino $\tilde{\chi}^{\pm}_{1}/\tilde{\chi}^{0}_{1}/\tilde{\chi}^{0}_{2}$ mass in the bilinear RPV model from Fig 14.

Observed 95% X-section upper limits as a function of higgsino $\tilde{\chi}^{0}_{1}/\tilde{\chi}^{0}_{2}$ mass in the UDD RPV model from Fig 18.

Observed 95% X-section upper limits as a function of wino $\tilde{\chi}^{\pm}_{1}/\tilde{\chi}^{0}_{2}$ mass in the WZ-mediated simplified model of wino $\tilde{\chi}^{\pm}_{1}/\tilde{\chi}^{0}_{2}$ production from Fig 9(aux).

N-1 distributions for $m_{T2}$ of observed data and expected background towards $SR^{WZ}_{high-m_{T2}}$ from publication's Figure 11(a) . The last bin is inclusive.

N-1 distributions for $m_{T2}$ of observed data and expected background towards $SR^{WZ}_{low-m_{T2}}$ from publication's Figure 11(b) . The last bin is inclusive.

N-1 distributions for $m_{T2}$ of observed data and expected background towards $SR^{bRPV}_{2l-SS}$ from publication's Figure 11(c) . The last bin is inclusive.

N-1 distributions for $m_{T2}$ of observed data and expected background towards $SR^{bRPV}_{3l}$ from publication's Figure 11(d) . The last bin is inclusive.

N-1 distributions for $\sum p^{b-jet}_{T}/\sum p^{jet}_{T}$ of observed data and expected background towards $SR^{RPV}_{2l1b}-L$ from publication's Figure 16(a) . The last bin is inclusive.

N-1 distributions for $\sum p^{b-jet}_{T}/\sum p^{jet}_{T}$ of observed data and expected background towards $SR^{RPV}_{2l2b}-M$ from publication's Figure 16(b) . The last bin is inclusive.

N-1 distributions for $\sum p^{b-jet}_{T}/\sum p^{jet}_{T}$ of observed data and expected background towards $SR^{RPV}_{2l3b}-H$ from publication's Figure 16(c) . The last bin is inclusive.

N-1 distribution for $E_{T}^{miss}$ in $SR^{Wh}_{high-m_{T2} }$ in ee channel

N-1 distribution for $E_{T}^{miss}$ in $SR^{Wh}_{high-m_{T2} }$ in e$\mu$ channel

N-1 distribution for $E_{T}^{miss}$ in $SR^{Wh}_{high-m_{T2} }$ in $\mu\mu$ channel

N-1 distribution for $\mathcal{S}(E_{T}^{miss})$ in $SR^{Wh}_{low-m_{T2} }$ in ee channel

N-1 distribution for $\mathcal{S}(E_{T}^{miss})$ in $SR^{Wh}_{low-m_{T2} }$ in e$\mu$ channel

N-1 distribution for $\mathcal{S}(E_{T}^{miss})$ in $SR^{Wh}_{low-m_{T2} }$ in $\mu\mu$ channel

Trigger efficiency for simulated signal events as a function of the LLP decay position in the x&ndash;y plane (L_xy) for LLPs decaying in the barrel (|&eta;|&lt;1.4) for one of the low-E_T signal samples for HLT CalRatio triggers seeded by the high-E_T L1 triggers with E_T thresholds of 60 GeV and 100 GeV and by the two versions of the low-E_T L1 triggers. Only statistical uncertainties are shown.

Trigger efficiency for simulated signal events as a function of the LLP decay position in the x&ndash;y plane (L_xy) for LLPs decaying in the barrel (|&eta;|&lt;1.4) for one of the high-E_T signal samples for HLT CalRatio triggers seeded by the high-E_T L1 triggers with E_T thresholds of 60 GeV and 100 GeV and by the two versions of the low-E_T L1 triggers. Only statistical uncertainties are shown.

The NN output scores in the dijet control region for the low-E_T training with no adversary network. Statistical uncertainties are shown in all plots. In cases where training with adversary networks is considered, the systematic uncertainty related to modelling discrepancies is included as well.

The NN output scores in the dijet control region for the low-E_T training with an adversary network included. Statistical uncertainties are shown in all plots. In cases where training with adversary networks is considered, the systematic uncertainty related to modelling discrepancies is included as well.

The NN output scores in the dijet control region for the high-E_T training with no adversary network. Statistical uncertainties are shown in all plots. In cases where training with adversary networks is considered, the systematic uncertainty related to modelling discrepancies is included as well.

The NN output scores in the dijet control region for the high-E_T training with an adversary network included. Statistical uncertainties are shown in all plots. In cases where training with adversary networks is considered, the systematic uncertainty related to modelling discrepancies is included as well.

Distribution of the low-E_T per-event BDT in main data, BIB data and some of the benchmark signal samples after preselection. Only statistical uncertainties are shown.

Distribution of the high-E_T per-event BDT outputs in main data, BIB data and some of the benchmark signal samples after preselection. Only statistical uncertainties are shown.

Sequential impact of each requirement on the number of events passing the selection for the high-E_T selections. The signal columns represent the cumulative fraction of events passing the selection than the number of events.

Sequential impact of each requirement on the number of events passing the selection for the low-E_T selections. The signal columns represent the cumulative fraction of events passing the selection than the number of events.

Application of the modified ABCD method to the final high-E_T selections. The a priori estimate refers to the "pre-unblinding" case, where the data in region A are ignored by removing the Poisson constraint in that region and the signal strength is fixed to zero. This matches the simple N^bkg_A=(N^bkg_B&middot; N^bkg_C)/N^bkg_D relation. The a posteriori estimate refers to the "post-unblinding" case, including the observed data in region A in the background-only global fit, obtained by fixing the signal strength to 0 (background-only fit) or allowing it to float (signal-plus-background fit). The table also shows one set of representative signal yields in each selection for the signal-plus-background fit. Only statistical uncertainties are included in the quoted error of the background, while the uncertainties in the signal include those from both statistical and experimental sources.

Application of the modified ABCD method to the final low-E_T selections. The a priori estimate refers to the "pre-unblinding" case, where the data in region A are ignored by removing the Poisson constraint in that region and the signal strength is fixed to zero. This matches the simple N^bkg_A=(N^bkg_B&middot; N^bkg_C)/N^bkg_D relation. The a posteriori estimate refers to the "post-unblinding" case, including the observed data in region A in the background-only global fit, obtained by fixing the signal strength to 0 (background-only fit) or allowing it to float (signal-plus-background fit). The table also shows one set of representative signal yields in each selection for the signal-plus-background fit. Only statistical uncertainties are included in the quoted error of the background, while the uncertainties in the signal include those from both statistical and experimental sources.

95% CL expected and observed limits on the BR of SM Higgs bosons to pairs of neutral LLPs (B_{H&rarr; ss}), showing the &plusmn; 1 &sigma; (green) and &plusmn; 2 &sigma; (yellow) expected limit bands, as well as a comparison with the results from previous ATLAS searches [36,78]. The cross-section for SM Higgs boson gluon--gluon fusion production is assumed to be 48.6 pb.

Trigger efficiency of simulated signal events as a function of the LLP decay position in the z direction for LLPs decaying in the calorimeter endcaps (1.4 &le; |&eta;| &lt; 2.5) for one of the low-E_T signal samples for HLT CalRatio triggers seeded by the 60 GeV-high-E_T trigger and by the two versions of the low-E_T triggers.

Trigger efficiency of simulated signal events as a function of the LLP decay position in the z direction for LLPs decaying in the calorimeter endcaps (1.4 &le; |&eta;| &lt; 2.5) for one of the high-E_T signal samples for HLT CalRatio triggers seeded by the 60 GeV-high-E_T trigger and by the two versions of the low-E_T triggers.

The simulated efficiencies as a function of c times the mean proper lifetime (c&tau;) of s for several different MC samples. A weight-based extrapolation procedure is used to determine the efficiency at a given mean proper lifetime.

The simulated efficiencies as a function of c times the mean proper lifetime (c&tau;) of s for several different MC samples. A weight-based extrapolation procedure is used to determine the efficiency at a given mean proper lifetime.

The simulated efficiencies as a function of c times the mean proper lifetime (c&tau;) of s for several different MC samples. A weight-based extrapolation procedure is used to determine the efficiency at a given mean proper lifetime.

The simulated efficiencies as a function of c times the mean proper lifetime (c&tau;) of s for several different MC samples. A weight-based extrapolation procedure is used to determine the efficiency at a given mean proper lifetime.

The simulated efficiencies as a function of c times the mean proper lifetime (c&tau;) of s for several different MC samples. A weight-based extrapolation procedure is used to determine the efficiency at a given mean proper lifetime.

The simulated efficiencies as a function of c times the mean proper lifetime (c&tau;) of s for several different MC samples. A weight-based extrapolation procedure is used to determine the efficiency at a given mean proper lifetime.

The simulated efficiencies as a function of c times the mean proper lifetime (c&tau;) of s for several different MC samples. A weight-based extrapolation procedure is used to determine the efficiency at a given mean proper lifetime.

The simulated efficiencies as a function of c times the mean proper lifetime (c&tau;) of s for several different MC samples. A weight-based extrapolation procedure is used to determine the efficiency at a given mean proper lifetime.

The simulated efficiencies as a function of c times the mean proper lifetime (c&tau;) of s for several different MC samples. A weight-based extrapolation procedure is used to determine the efficiency at a given mean proper lifetime.

The simulated efficiencies as a function of c times the mean proper lifetime (c&tau;) of s for several different MC samples. A weight-based extrapolation procedure is used to determine the efficiency at a given mean proper lifetime.

The simulated efficiencies as a function of c times the mean proper lifetime (c&tau;) of s for several different MC samples. A weight-based extrapolation procedure is used to determine the efficiency at a given mean proper lifetime.

The simulated efficiencies as a function of c times the mean proper lifetime (c&tau;) of s for several different MC samples. A weight-based extrapolation procedure is used to determine the efficiency at a given mean proper lifetime.

The simulated efficiencies as a function of c times the mean proper lifetime (c&tau;) of s for several different MC samples. A weight-based extrapolation procedure is used to determine the efficiency at a given mean proper lifetime.

The simulated efficiencies as a function of c times the mean proper lifetime (c&tau;) of s for several different MC samples. A weight-based extrapolation procedure is used to determine the efficiency at a given mean proper lifetime.

The event BDT in the dijet control region for the low-E_T training training. The uncertainties are combined statistical and ML modelling systematic uncertainty.

The event BDT in the dijet control region for the high-E_T training. The uncertainties are combined statistical and ML modelling systematic uncertainty.

The 95% CL observed limits, expected limits and &plusmn; 1 &sigma; and 2 &sigma; bands for the &Phi; mass of 60 GeV compared to the results from the 2016-data analysis and subsequent combinations, where available.

The 95% CL observed limits, expected limits and &plusmn; 1 &sigma; and 2 &sigma; bands for the &Phi; mass of 60 GeV compared to the results from the 2016-data analysis and subsequent combinations, where available.

The 95% CL observed limits, expected limits and &plusmn; 1 &sigma; and 2 &sigma; bands for the &Phi; mass 125 GeV, compared to the results from the 2016-data analysis and subsequent combinations, where available.

The 95% CL observed limits, expected limits and &plusmn; 1 &sigma; and 2 &sigma; bands for the &Phi; mass 125 GeV, compared to the results from the 2016-data analysis and subsequent combinations, where available. The 125 GeV mediator is assumed to be the SM Higgs boson.

The 95% CL observed limits, expected limits and &plusmn; 1 &sigma; and 2 &sigma; bands for the &Phi; mass 125 GeV, compared to the results from the 2016-data analysis and subsequent combinations, where available. The 125 GeV mediator is assumed to be the SM Higgs boson.

The 95% CL observed limits, expected limits and &plusmn; 1 &sigma; and 2 &sigma; bands for the &Phi; mass of 200 GeV, compared to the results from the 2016-data analysis and subsequent combinations, where available.

The 95% CL observed limits, expected limits and &plusmn; 1 &sigma; and 2 &sigma; bands for the &Phi; mass of 400 GeV, compared to the results from the 2016-data analysis and subsequent combinations, where available.

The 95% CL observed limits, expected limits and &plusmn; 1 &sigma; and 2 &sigma; bands for the &Phi; mass 600 GeV, compared to the results from the 2016-data analysis and subsequent combinations, where available.

The 95% CL observed limits, expected limits and &plusmn; 1 &sigma; and 2 &sigma; bands for the &Phi; mass of 600 GeV, compared to the results from the 2016-data analysis and subsequent combinations, where available.

The 95% CL observed limits, expected limits and &plusmn; 1 &sigma; and 2 &sigma; bands for the &Phi; masses of of 1000 GeV, compared to the results from the 2016-data analysis and subsequent combinations, where available.

The 95% CL observed limits, expected limits and &plusmn; 1 &sigma; and 2 &sigma; bands for the &Phi; masses of of 1000 GeV, compared to the results from the 2016-data analysis and subsequent combinations, where available.

The 95% CL observed limits, expected limits and &plusmn; 1 &sigma; and 2 &sigma; bands for the &Phi; masses of of 1000 GeV, compared to the results from the 2016-data analysis and subsequent combinations, where available.

The 95% CL observed limits, expected limits and &plusmn; 1 &sigma; and 2 &sigma; bands for the &Phi; masses of of 1000 GeV, compared to the results from the 2016-data analysis and subsequent combinations, where available.

Efficiency for an event to enter Region A of the high-E_T selection, as a function of the two LLP kinematics, decay type and decay position. LLP p_T is binned in the ranges of [0, 50 , 100, 200 , 400, 1600] GeV (5 bins), the LLP decay position is binned in decay position in L_xy in [0, 1.5, 2, 2.5, 3, 3.5, 3.9, &infin;] m for LLPs with |&eta;| &lt; 1.5 and L_z in [0, 3.6, 4.2, 4.8, 5.5, 6, &infin; ] m for LLPs with |&eta;|geq1.5, 13 bins in all. Finally, four decay types are considered: LLPs decaying to pairs of c, b, t, or &tau; in bins 0, 1, 2, 3 respectively. The efficiency is presented as a function of "Bin Index", which is calculated as follows: Bin Index = (decay position bin index ) &times; (number of p_T bins &times; number of decay type bins) + p_T bin index * (number of decay type bin) + decay type bin index. The efficiency for a given pair-produced LLP sample can be obtained by summing the efficiency values for each event as extracted from this map, and dividing it by the total number of events in the sample. The efficiency map is symmetric between the LLPs, so the choice of LLP1 and LLP2 is arbitrary. For the high-E_T selections, for overall efficiencies above 0.5%, the results are typically accurate to around 25%, but below this the efficiency can be overestimated and therefore should not be used for re-interpration.

Efficiency for an event to enter Region A of the low-E_T selection, as a function of the two LLP kinematics, decay type and decay position. LLP p_T is binned in the ranges of [0, 50 , 100, 200 , 400, 1600] GeV (5 bins), the LLP decay position is binned in L_xy in [0, 1.5, 2, 2.5, 3, 3.5, 3.9, &infin;] m for LLPs with |&eta;| &lt; 1.5 and L_z in [0, 3.6, 4.2, 4.8, 5.5, 6, &infin; ] m for LLPs with |&eta;|geq1.5, 13 bins in all. Finally, four decay types are considered: LLPs decaying to pairs of c, b, t, or &tau; in bins 0, 1, 2, 3 respectively. The efficiency is presented as a function of "Bin Index", which is calculated as follows: Bin Index = (decay position bin index ) &times; (number of p_T bins &times; number of decay type bins) + p_T bin index * (number of decay type bin) + decay type bin index. The efficiency for a given pair-produced LLP sample can be obtained by summing the efficiency values for each event as extracted from this map, and dividing it by the total number of events in the sample. The efficiency map is symmetric between the LLPs, so the choice of LLP1 and LLP2 is arbitrary. For the low-E_T selections, for overall efficiencies above 0.15%, the results are typically accurate to around 33%, and below this the efficiency is typically accurate up to a factor of 3.