The ratios (μ) of the 95% C.L. upper limits on the combined s→ W^±W^∓ and s→ 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 (μ) of the 95% C.L. upper limits on the combined s→ W^±W^∓ and s→ 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% C.L. on σ(pp → s χχ) × B(s→ 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 σ(pp → s χχ) × B(s→ VV) process for m_Z'=0.5 TeV, shown in dashed blue.

Observed upper limits at 95% C.L. on σ(pp → s χχ) × B(s→ 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 σ(pp → s χχ) × B(s→ VV) process for m_Z'=1 TeV, shown in dashed blue.

Observed upper limits at 95% C.L. on σ(pp → s χχ) × B(s→ 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 σ(pp → s χχ) × B(s→ 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, Δ φ_{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>60 GeV that is supplemented by up to two additional small-R jets within ΔR_cone=2.5 of the TAR jet. Here, Δ φ_{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 × ϵ), 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^±W^∓) dark Higgs signal points with dark Higgs boson mass m_s and Z' boson mass m_Z'.

The product of acceptance and efficiency (A × ϵ), 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'.

The observed exclusion contour at 95% CL as a function of the $\it{m}_{\tilde{\chi}^{0}_{1}}$ vs. $\it{m}_{\tilde{t}}$. Masses that are within the contours are excluded. Limits are shown in the case of a higgsino LSP. The results are constrained by the kinematic limits of the top-squark decay into a chargino and a bottom quark (upper diagonal line) and into a neutralino and a top quark (lower diagonal line), respectively.

The expected exclusion contour at 95% CL as a function of the $\it{m}_{\tilde{\chi}^{0}_{1}}$ vs. $\it{m}_{\tilde{t}}$. Masses that are within the contours are excluded. Limits are shown in the case of a higgsino LSP. The results are constrained by the kinematic limits of the top-squark decay into a chargino and a bottom quark (upper diagonal line) and into a neutralino and a top quark (lower diagonal line), respectively.

The observed exclusion contour at 95% CL as a function of the $\it{m}_{\tilde{\chi}^{0}_{1}}$ vs. $\it{m}_{\tilde{t}}$. Masses that are within the contours are excluded. Limits are shown for the region $m_{\tilde{t}} - m_{\tilde{\chi}^0_{1,2}, \tilde{\chi}^\pm_{1}} \geq m_\text{top}$ where $B(\tilde{t} \rightarrow b \chi^{+}_{1}) = B(\tilde{t} \rightarrow t \chi^{0}_{1,2}) = 0.5$.

The expected exclusion contour at 95% CL as a function of the $\it{m}_{\tilde{\chi}^{0}_{1}}$ vs. $\it{m}_{\tilde{t}}$. Masses that are within the contours are excluded. Limits are shown for the region $m_{\tilde{t}} - m_{\tilde{\chi}^0_{1,2}, \tilde{\chi}^\pm_{1}} \geq m_\text{top}$ where $B(\tilde{t} \rightarrow b \chi^{+}_{1}) = B(\tilde{t} \rightarrow t \chi^{0}_{1,2}) = 0.5$.

Observed model-dependent upper limit on the cross section for the $(\tilde{t},\tilde{\chi}^{\pm}_{1})$ signal grid. Limits are shown for $B(\tilde{t} \rightarrow b \chi^{+}_{1})$ equal to unity.

Observed model-dependent upper limit on the cross section for the $(\tilde{t},\tilde{\chi}^{\pm}_{1} / \tilde{\chi}^{0}_{1,2})$ signal grid. Limits are shown in the case of a higgsino LSP. The results are constrained by the kinematic limits of the top-squark decay into a chargino and a bottom quark (upper diagonal line) and into a neutralino and a top quark (lower diagonal line), respectively.

Expected model-dependent upper limit on the cross section for the $(\tilde{t},\tilde{\chi}^{\pm}_{1})$ signal grid. Limits are shown for $B(\tilde{t} \rightarrow b \chi^{+}_{1})$ equal to unity.

Expected model-dependent upper limit on the cross section for the $(\tilde{t},\tilde{\chi}^{\pm}_{1} / \tilde{\chi}^{0}_{1,2})$ signal grid. Limits are shown in the case of a higgsino LSP. The results are constrained by the kinematic limits of the top-squark decay into a chargino and a bottom quark (upper diagonal line) and into a neutralino and a top quark (lower diagonal line), respectively.

Expected background and observed number of events in different jet and $b$-tag multiplicity bins.

Cut flow for a model of top-squark pair production with the top squark decaying to a $b$-quark and a chargino. The chargino decays through the non-zero RPV coupling $\lambda^{''}_{323}$ via a virtual top squark to $bbs$ quark triplets ($m_{\tilde{t}}$ = 800 GeV, $m_{\tilde{\chi}^{\pm}_{1}}$ = 750 GeV). The multijet trigger consists of four jets satisfying $p_{\text{T}}\geq(100)120$ GeV for the 2015-2016 (2017-2018) data period. Selections with negligible inefficiencies on the given sample, such as data quality requirements, are not displayed. The numbers in $N_{\text{weighted}}$ are normalized by the integrated luminosity of 139 fb$^{-1}$.

Signal acceptance for $\tilde{t} \rightarrow b\tilde{\chi}^{+}_{1}(\tilde{\chi}^{+}_{1} \rightarrow \bar{b}\bar{b}\bar{s}) $ and c.c. model. Please mind that the acceptance given in the table is reported in %.

Signal acceptance for $\tilde{t} \rightarrow b\tilde{\chi}^{+}_{1}(\tilde{\chi}^{+}_{1} \rightarrow \bar{b}\bar{b}\bar{s}) $ and c.c. model. Please mind that the acceptance given in the table is reported in %.

Signal acceptance for $\tilde{t} \rightarrow b\tilde{\chi}^{+}_{1}(\tilde{\chi}^{+}_{1} \rightarrow \bar{b}\bar{b}\bar{s}) $ and c.c. model. Please mind that the acceptance given in the table is reported in %.

Signal acceptance for $\tilde{t} \rightarrow b\tilde{\chi}^{+}_{1}(\tilde{\chi}^{+}_{1} \rightarrow \bar{b}\bar{b}\bar{s}) $ and c.c. model. Please mind that the acceptance given in the table is reported in %.

Signal acceptance for $\tilde{t} \rightarrow b\tilde{\chi}^{+}_{1}(\tilde{\chi}^{+}_{1} \rightarrow \bar{b}\bar{b}\bar{s}) $ and c.c. model. Please mind that the acceptance given in the table is reported in %.

Signal acceptance for $\tilde{t} \rightarrow b\tilde{\chi}^{+}_{1}(\tilde{\chi}^{+}_{1} \rightarrow \bar{b}\bar{b}\bar{s}) $ and c.c. model. Please mind that the acceptance given in the table is reported in %.

Signal acceptance for $\tilde{t} \rightarrow b\tilde{\chi}^{+}_{1}(\tilde{\chi}^{+}_{1} \rightarrow \bar{b}\bar{b}\bar{s}) $ and c.c. model. Please mind that the acceptance given in the table is reported in %.

Signal acceptance for $\tilde{t} \rightarrow b\tilde{\chi}^{+}_{1}(\tilde{\chi}^{+}_{1} \rightarrow \bar{b}\bar{b}\bar{s}) $ and c.c. model. Please mind that the acceptance given in the table is reported in %.

Signal acceptance for $\tilde{t} \rightarrow b\tilde{\chi}^{+}_{1}(\tilde{\chi}^{+}_{1} \rightarrow \bar{b}\bar{b}\bar{s}) $ and c.c. model. Please mind that the acceptance given in the table is reported in %.

Signal efficiency for $\tilde{t} \rightarrow b\tilde{\chi}^{+}_{1}(\tilde{\chi}^{+}_{1} \rightarrow \bar{b}\bar{b}\bar{s}) $ and c.c. model. Please mind that the efficiency given in the table is reported in %.

Signal efficiency for $\tilde{t} \rightarrow b\tilde{\chi}^{+}_{1}(\tilde{\chi}^{+}_{1} \rightarrow \bar{b}\bar{b}\bar{s}) $ and c.c. model. Please mind that the efficiency given in the table is reported in %.

Signal efficiency for $\tilde{t} \rightarrow b\tilde{\chi}^{+}_{1}(\tilde{\chi}^{+}_{1} \rightarrow \bar{b}\bar{b}\bar{s}) $ and c.c. model. Please mind that the efficiency given in the table is reported in %.

Signal efficiency for $\tilde{t} \rightarrow b\tilde{\chi}^{+}_{1}(\tilde{\chi}^{+}_{1} \rightarrow \bar{b}\bar{b}\bar{s}) $ and c.c. model. Please mind that the efficiency given in the table is reported in %.

Signal efficiency for $\tilde{t} \rightarrow b\tilde{\chi}^{+}_{1}(\tilde{\chi}^{+}_{1} \rightarrow \bar{b}\bar{b}\bar{s}) $ and c.c. model. Please mind that the efficiency given in the table is reported in %.

Signal efficiency for $\tilde{t} \rightarrow b\tilde{\chi}^{+}_{1}(\tilde{\chi}^{+}_{1} \rightarrow \bar{b}\bar{b}\bar{s}) $ and c.c. model. Please mind that the efficiency given in the table is reported in %.

Signal efficiency for $\tilde{t} \rightarrow b\tilde{\chi}^{+}_{1}(\tilde{\chi}^{+}_{1} \rightarrow \bar{b}\bar{b}\bar{s}) $ and c.c. model. Please mind that the efficiency given in the table is reported in %.

Signal efficiency for $\tilde{t} \rightarrow b\tilde{\chi}^{+}_{1}(\tilde{\chi}^{+}_{1} \rightarrow \bar{b}\bar{b}\bar{s}) $ and c.c. model. Please mind that the efficiency given in the table is reported in %.

Signal efficiency for $\tilde{t} \rightarrow b\tilde{\chi}^{+}_{1}(\tilde{\chi}^{+}_{1} \rightarrow \bar{b}\bar{b}\bar{s}) $ and c.c. model. Please mind that the efficiency given in the table is reported in %.

Signal acceptance for $\tilde{t} \rightarrow t\tilde{\chi}^{0}_{1,2}(\tilde{\chi}^{0}_{1,2} \rightarrow tbs) / b\tilde{\chi}^{+}_{1}(\tilde{\chi}^{+}_{1} \rightarrow \bar{b}\bar{b}\bar{s}) $ and c.c. model. Please mind that the acceptance given in the table is reported in %.

Signal acceptance for $\tilde{t} \rightarrow t\tilde{\chi}^{0}_{1,2}(\tilde{\chi}^{0}_{1,2} \rightarrow tbs) / b\tilde{\chi}^{+}_{1}(\tilde{\chi}^{+}_{1} \rightarrow \bar{b}\bar{b}\bar{s}) $ and c.c. model. Please mind that the acceptance given in the table is reported in %.

Signal acceptance for $\tilde{t} \rightarrow t\tilde{\chi}^{0}_{1,2}(\tilde{\chi}^{0}_{1,2} \rightarrow tbs) / b\tilde{\chi}^{+}_{1}(\tilde{\chi}^{+}_{1} \rightarrow \bar{b}\bar{b}\bar{s}) $ and c.c. model. Please mind that the acceptance given in the table is reported in %.

Signal acceptance for $\tilde{t} \rightarrow t\tilde{\chi}^{0}_{1,2}(\tilde{\chi}^{0}_{1,2} \rightarrow tbs) / b\tilde{\chi}^{+}_{1}(\tilde{\chi}^{+}_{1} \rightarrow \bar{b}\bar{b}\bar{s}) $ and c.c. model. Please mind that the acceptance given in the table is reported in %.

Signal acceptance for $\tilde{t} \rightarrow t\tilde{\chi}^{0}_{1,2}(\tilde{\chi}^{0}_{1,2} \rightarrow tbs) / b\tilde{\chi}^{+}_{1}(\tilde{\chi}^{+}_{1} \rightarrow \bar{b}\bar{b}\bar{s}) $ and c.c. model. Please mind that the acceptance given in the table is reported in %.

Signal acceptance for $\tilde{t} \rightarrow t\tilde{\chi}^{0}_{1,2}(\tilde{\chi}^{0}_{1,2} \rightarrow tbs) / b\tilde{\chi}^{+}_{1}(\tilde{\chi}^{+}_{1} \rightarrow \bar{b}\bar{b}\bar{s}) $ and c.c. model. Please mind that the acceptance given in the table is reported in %.

Signal acceptance for $\tilde{t} \rightarrow t\tilde{\chi}^{0}_{1,2}(\tilde{\chi}^{0}_{1,2} \rightarrow tbs) / b\tilde{\chi}^{+}_{1}(\tilde{\chi}^{+}_{1} \rightarrow \bar{b}\bar{b}\bar{s}) $ and c.c. model. Please mind that the acceptance given in the table is reported in %.

Signal acceptance for $\tilde{t} \rightarrow t\tilde{\chi}^{0}_{1,2}(\tilde{\chi}^{0}_{1,2} \rightarrow tbs) / b\tilde{\chi}^{+}_{1}(\tilde{\chi}^{+}_{1} \rightarrow \bar{b}\bar{b}\bar{s}) $ and c.c. model. Please mind that the acceptance given in the table is reported in %.

Signal acceptance for $\tilde{t} \rightarrow t\tilde{\chi}^{0}_{1,2}(\tilde{\chi}^{0}_{1,2} \rightarrow tbs) / b\tilde{\chi}^{+}_{1}(\tilde{\chi}^{+}_{1} \rightarrow \bar{b}\bar{b}\bar{s}) $ and c.c. model. Please mind that the acceptance given in the table is reported in %.

Signal acceptance for $\tilde{t} \rightarrow t\tilde{\chi}^{0}_{1,2}(\tilde{\chi}^{0}_{1,2} \rightarrow tbs) / b\tilde{\chi}^{+}_{1}(\tilde{\chi}^{+}_{1} \rightarrow \bar{b}\bar{b}\bar{s}) $ and c.c. model. Please mind that the acceptance given in the table is reported in %.

Signal acceptance for $\tilde{t} \rightarrow t\tilde{\chi}^{0}_{1,2}(\tilde{\chi}^{0}_{1,2} \rightarrow tbs) / b\tilde{\chi}^{+}_{1}(\tilde{\chi}^{+}_{1} \rightarrow \bar{b}\bar{b}\bar{s}) $ and c.c. model. Please mind that the acceptance given in the table is reported in %.

Signal acceptance for $\tilde{t} \rightarrow t\tilde{\chi}^{0}_{1,2}(\tilde{\chi}^{0}_{1,2} \rightarrow tbs) / b\tilde{\chi}^{+}_{1}(\tilde{\chi}^{+}_{1} \rightarrow \bar{b}\bar{b}\bar{s}) $ and c.c. model. Please mind that the acceptance given in the table is reported in %.

Signal acceptance for $\tilde{t} \rightarrow t\tilde{\chi}^{0}_{1,2}(\tilde{\chi}^{0}_{1,2} \rightarrow tbs) / b\tilde{\chi}^{+}_{1}(\tilde{\chi}^{+}_{1} \rightarrow \bar{b}\bar{b}\bar{s}) $ and c.c. model. Please mind that the acceptance given in the table is reported in %.

Signal acceptance for $\tilde{t} \rightarrow t\tilde{\chi}^{0}_{1,2}(\tilde{\chi}^{0}_{1,2} \rightarrow tbs) / b\tilde{\chi}^{+}_{1}(\tilde{\chi}^{+}_{1} \rightarrow \bar{b}\bar{b}\bar{s}) $ and c.c. model. Please mind that the acceptance given in the table is reported in %.

Signal acceptance for $\tilde{t} \rightarrow t\tilde{\chi}^{0}_{1,2}(\tilde{\chi}^{0}_{1,2} \rightarrow tbs) / b\tilde{\chi}^{+}_{1}(\tilde{\chi}^{+}_{1} \rightarrow \bar{b}\bar{b}\bar{s}) $ and c.c. model. Please mind that the acceptance given in the table is reported in %.

Signal acceptance for $\tilde{t} \rightarrow t\tilde{\chi}^{0}_{1,2}(\tilde{\chi}^{0}_{1,2} \rightarrow tbs) / b\tilde{\chi}^{+}_{1}(\tilde{\chi}^{+}_{1} \rightarrow \bar{b}\bar{b}\bar{s}) $ and c.c. model. Please mind that the acceptance given in the table is reported in %.

Signal acceptance for $\tilde{t} \rightarrow t\tilde{\chi}^{0}_{1,2}(\tilde{\chi}^{0}_{1,2} \rightarrow tbs) / b\tilde{\chi}^{+}_{1}(\tilde{\chi}^{+}_{1} \rightarrow \bar{b}\bar{b}\bar{s}) $ and c.c. model. Please mind that the acceptance given in the table is reported in %.

Signal acceptance for $\tilde{t} \rightarrow t\tilde{\chi}^{0}_{1,2}(\tilde{\chi}^{0}_{1,2} \rightarrow tbs) / b\tilde{\chi}^{+}_{1}(\tilde{\chi}^{+}_{1} \rightarrow \bar{b}\bar{b}\bar{s}) $ and c.c. model. Please mind that the acceptance given in the table is reported in %.

Post-fit yields for data and prediction in each of the single-bin signal regions of the analysis.

Observed 95% confidence level limit for the two-step signal grid.

Observed 95% confidence level limit for the two-step signal grid with the signal cross section increased by one sigma.

Observed 95% confidence level limit for the two-step signal grid with the signal cross section decreased by one sigma.

Expected 95% confidence level limit for the two-step signal grid.

Expected 95% confidence level limit for the two-step signal grid plus one sigma from experimental systematics.

Expected 95% confidence level limit for the two-step signal grid minus one sigma from experimental systematics.

Observed 95% confidence level limit for the Gtt signal grid.

Observed 95% confidence level limit for the Gtt signal grid with the signal cross section increased by one sigma.

Observed 95% confidence level limit for the Gtt signal grid with the signal cross section decreased by one sigma.

Expected 95% confidence level limit for the Gtt signal grid.

Expected 95% confidence level limit for the Gtt signal grid plus one sigma from experimental systematics.

Expected 95% confidence level limit for the Gtt signal grid minus one sigma from experimental systematics.

Observed 95% confidence level limit for the RPV signal grid.

Observed 95% confidence level limit for the RPV signal grid with the signal cross section increased by one sigma.

Observed 95% confidence level limit for the RPV signal grid with the signal cross section decreased by one sigma.

Expected 95% confidence level limit for the RPV signal grid.

Expected 95% confidence level limit for the RPV signal grid plus one sigma from experimental systematics.

Expected 95% confidence level limit for the RPV signal grid minus one sigma from experimental systematics.

Observed 95% confidence level limit for the two-step signal grid.

Expected 95% confidence level limit for the two-step signal grid.

Observed 95% confidence level limit for the Gtt signal grid.

Expected 95% confidence level limit for the Gtt signal grid.

Observed 95% confidence level limit for the RPV signal grid.

Expected 95% confidence level limit for the RPV signal grid.

$\mathcal{S}(E_{\mathrm{T}}^{\mathrm{miss}})$ distribution in the signal region SR-10ij50-0ib-MJ340. Two benchmark signal models are shown along with the background yields. These models, each representing a single mass point, are labelled 'RPV' with $(m_{\tilde{g}}, m_{\tilde{t}}) = (1600, 600) \, \mathrm{GeV}$ and 'two-step' with $(m_{\tilde{g}}, m_{\tilde{\chi^{0}_{1}}}) = (1600, 100) \, \mathrm{GeV}$.

$\mathcal{S}(E_{\mathrm{T}}^{\mathrm{miss}})$ distribution in the signal region SR-12ij50-2ib. Two benchmark signal models are shown along with the background yields. These models, each representing a single mass point, are labelled 'RPV' with $(m_{\tilde{g}}, m_{\tilde{t}}) = (1600, 600) \, \mathrm{GeV}$ and 'two-step' with $(m_{\tilde{g}}, m_{\tilde{\chi^{0}_{1}}}) = (1600, 100) \, \mathrm{GeV}$.

$\mathcal{S}(E_{\mathrm{T}}^{\mathrm{miss}})$ distribution in the signal region SR-9ij80-0ib. Two benchmark signal models are shown along with the background yields. These models, each representing a single mass point, are labelled 'RPV' with $(m_{\tilde{g}}, m_{\tilde{t}}) = (1600, 600) \, \mathrm{GeV}$ and 'two-step' with $(m_{\tilde{g}}, m_{\tilde{\chi^{0}_{1}}}) = (1600, 100) \, \mathrm{GeV}$.

Number of signal events expected for $139 \, \mathrm{fb}^{-1} $ after different analysis selections in the signal region SR-8ij50-0ib-MJ500. This 'two-step' model requires that a strongly produced gluino decays into quarks, the W and Z bosons, and the lightest stable neutralino where $(m_{\tilde{g}}, m_{\tilde{\chi^{0}_{1}}}) = (1600, 100) \, \mathrm{GeV}$.

Number of signal events expected for $139 \, \mathrm{fb}^{-1} $ after different analysis selections in the signal region SR-9ij50-0ib-MJ340. This 'two-step' model requires that a strongly produced gluino decays into quarks, the W and Z bosons, and the lightest stable neutralino where $(m_{\tilde{g}}, m_{\tilde{\chi^{0}_{1}}}) = (1600, 100) \, \mathrm{GeV}$.

Number of signal events expected for $139 \, \mathrm{fb}^{-1} $ after different analysis selections in the signal region SR-10ij50-0ib-MJ340. This 'two-step' model requires that a strongly produced gluino decays into quarks, the W and Z bosons, and the lightest stable neutralino where $(m_{\tilde{g}}, m_{\tilde{\chi^{0}_{1}}}) = (1600, 100) \, \mathrm{GeV}$.

Number of signal events expected for $139 \, \mathrm{fb}^{-1} $ after different analysis selections in the signal region SR-10ij50-0ib-MJ500. This 'two-step' model requires that a strongly produced gluino decays into quarks, the W and Z bosons, and the lightest stable neutralino where $(m_{\tilde{g}}, m_{\tilde{\chi^{0}_{1}}}) = (1600, 100) \, \mathrm{GeV}$.

Number of signal events expected for $139 \, \mathrm{fb}^{-1} $ after different analysis selections in the signal region SR-10ij50-1ib-MJ500. This 'two-step' model requires that a strongly produced gluino decays into quarks, the W and Z bosons, and the lightest stable neutralino where $(m_{\tilde{g}}, m_{\tilde{\chi^{0}_{1}}}) = (1600, 100) \, \mathrm{GeV}$.

Number of signal events expected for $139 \, \mathrm{fb}^{-1} $ after different analysis selections in the signal region SR-11ij50-0ib. This 'two-step' model requires that a strongly produced gluino decays into quarks, the W and Z bosons, and the lightest stable neutralino where $(m_{\tilde{g}}, m_{\tilde{\chi^{0}_{1}}}) = (1600, 100) \, \mathrm{GeV}$.

Number of signal events expected for $139 \, \mathrm{fb}^{-1} $ after different analysis selections in the signal region SR-12ij50-2ib. This 'two-step' model requires that a strongly produced gluino decays into quarks, the W and Z bosons, and the lightest stable neutralino where $(m_{\tilde{g}}, m_{\tilde{\chi^{0}_{1}}}) = (1600, 100) \, \mathrm{GeV}$.

Number of signal events expected for $139 \, \mathrm{fb}^{-1} $ after different analysis selections in the signal region SR-9ij80-0ib. This 'two-step' model requires that a strongly produced gluino decays into quarks, the W and Z bosons, and the lightest stable neutralino where $(m_{\tilde{g}}, m_{\tilde{\chi^{0}_{1}}}) = (1600, 100) \, \mathrm{GeV}$.

Acceptance for the signal region SR-8ij50-0ib-MJ500 showing the acceptance for the complete two-step signal grid.

Efficiency for the signal region SR-8ij50-0ib-MJ500 showing the efficiency for the complete two-step signal grid.

Acceptance for the signal region SR-9ij50-0ib-MJ340 showing the acceptance for the complete two-step signal grid.

Efficiency for the signal region SR-9ij50-0ib-MJ340 showing the efficiency for the complete two-step signal grid.

Acceptance for the signal region SR-10ij50-0ib-MJ340 showing the acceptance for the complete two-step signal grid.

Efficiency for the signal region SR-10ij50-0ib-MJ340 showing the efficiency for the complete two-step signal grid.

Acceptance for the signal region SR-10ij50-0ib-MJ500 showing the acceptance for the complete two-step signal grid.

Efficiency for the signal region SR-10ij50-0ib-MJ500 showing the efficiency for the complete two-step signal grid.

Acceptance for the signal region SR-10ij50-1ib-MJ500 showing the acceptance for the complete two-step signal grid.

Efficiency for the signal region SR-10ij50-1ib-MJ500 showing the efficiency for the complete two-step signal grid.

Acceptance for the signal region SR-11ij50-0ib showing the acceptance for the complete two-step signal grid.

Efficiency for the signal region SR-11ij50-0ib showing the efficiency for the complete two-step signal grid.

Acceptance for the signal region SR-12ij50-2ib showing the acceptance for the complete two-step signal grid.

Efficiency for the signal region SR-12ij50-2ib showing the efficiency for the complete two-step signal grid.

Acceptance for the signal region SR-9ij80-0ib showing the acceptance for the complete two-step signal grid.

Efficiency for the signal region SR-9ij80-0ib showing the efficiency for the complete two-step signal grid.

The normalisation factors for the dominant backgrounds of the analysis in each of the multi-bin and single-bin regions.

Post-fit yields for data and prediction in each of the single-bin validation regions to test the $N_{\mathrm{jet}}$ extraction.

Post-fit yields for data and prediction in each of the single-bin validation regions to test the $\mathcal{S}(E_{\mathrm{T}}^{\mathrm{miss}})$ extrapolation.

Post-fit yields for data and prediction in each of the multi-bin validation regions to test the $N_{\mathrm{jet}}$ extraction.

Post-fit yields for data and prediction in each of the multi-bin validation regions to test the $\mathcal{S}(E_{\mathrm{T}}^{\mathrm{miss}})$ extrapolation.

The observed Cls from the best expected signal regions for the two-step decay.

The observed Cls from the best expected signal regions for the Gtt decay.

The observed Cls from the best expected signal regions for the RPV decay.

Number of events in each signal region broken down by background type and the number of observed data events.

From left to right; the $95\%$ CL upper limits on the visible cross section (${\langle \epsilon\sigma \rangle}^{95}_{obs}$) and on the number of signal events. Next is the $95\%$ CL upper limit on the number of signal events, given the expected number of background events. The last two columns show the confidence level for the background only hypothesis ($CL_{b}$) and the dicovery $p$-value along with the Gaussian significance (Z).

Visualisation of the highest jet multiplicity event selected in signal regions targeting long cascade decays of pair-produced gluinos. This event was recorded by ATLAS on 23 October 2016, and contains 16 jets, illustrated by cones. Yellow blocks represent the calorimeter energy measured in noise-suppressed clusters. Of the reconstructed jets, 13 (11) have transverse momenta above 50 GeV (80 GeV), with 3 (2) being b-tagged. The leading jet has a transverse momentum of 507 GeV, and the sum of jet transverse momenta $H_T=2.9$ TeV. A value of 343 GeV is observed for the $E_{T}^{miss}$, whose direction is shown by the dashed red line, producing a significance $S(E_{T}^{miss})=6.4$. The sum of the masses of large-radius jets is evaluated as $M_{J}^{\Sigma}=1070$ GeV.

Visualisation of the highest jet multiplicity event selected in a control region used to make predictions of the background from multijet production. This event was recorded by ATLAS on 18 July 2018, and contains 19 jets, illustrated by cones. Yellow blocks represent the calorimeter energy measured in in noise-suppressed clusters. Of the reconstructed jets, 16 (10) have transverse momenta above 50 GeV (80 GeV). No jets were b-tagged. The leading et has a transverse momentum of 371 GeV, and the sum of jet transverse momenta $H_T=2.2$ TeV. A value of 8 GeV is observed for the $E_{T}^{miss}$, whose direction is shown by the dashed red line, producing a significance $S(E_{T}^{miss})=0.2$. The sum of the masses of large-radius jets is evaluated as $M_{J}^{\Sigma}=767$ GeV.

Expected and observed 95% CL upper limits on the production of a heavy, narrow-width, scalar resonance decaying to a pair of Higgs bosons ($X\rightarrow HH$). The final state used in the search consists of a boosted $b\bar{b}$ pair and a boosted hadronically decaying $\tau^{+}\tau^{-}$ pair, and the SM braching ratio of the Higgs boson are assumed. The $\pm 1\sigma$ and $\pm 2\sigma$ variations about the expected limit are indicated by the error bands. Two different requirements are applied on the visible mass of the two boosted Higgs boson candidates for the resonance mass hypotheses of 1.6 TeV and 2.5 TeV, leading to discontinuities in the limits (at 1.6 TeV, the difference between imposing no requirement and $m^{vis}_{HH}>900$ GeV is less than 1% though).

Distribution of the resonance mass in the c-tag Signal Region of the $ ce$ channel for the post-fit background, the observed data, and the expected signal with $m_{LQ} = 1$ TeV.

Distribution of the resonance mass in the b-tag Signal Region of the $ ce$ channel for the post-fit background, the observed data, and the expected signal with $m_{LQ} = 1$ TeV.

Distribution of the resonance mass in the untagged Signal Region of the $ c\mu$ channel for the post-fit background, the observed data, and the expected signal with $m_{LQ} = 1$ TeV.

Distribution of the resonance mass in the c-tag Signal Region of the $ c\mu$ channel for the post-fit background, the observed data, and the expected signal with $m_{LQ} = 1$ TeV.

Distribution of the resonance mass in the b-tag Signal Region of the $ c\mu$ channel for the post-fit background, the observed data, and the expected signal with $m_{LQ} = 1$ TeV.

Distribution of the resonance mass in the 0-tag Signal Region of the $ be$ channel for the post-fit background, the observed data, and the expected signal with $m_{LQ} = 1$ TeV.

Distribution of the resonance mass in the 1-tag Signal Region of the $ be$ channel for the post-fit background, the observed data, and the expected signal with $m_{LQ} = 1$ TeV.

Distribution of the resonance mass in the 2-tag Signal Region of the $ be$ channel for the post-fit background, the observed data, and the expected signal with $m_{LQ} = 1$ TeV.

Distribution of the resonance mass in the 0-tag Signal Region of the $ b\mu$ channel for the post-fit background, the observed data, and the expected signal with $m_{LQ} = 1$ TeV.

Distribution of the resonance mass in the 1-tag Signal Region of the $ b\mu$ channel for the post-fit background, the observed data, and the expected signal with $m_{LQ} = 1$ TeV.

Distribution of the resonance mass in the 2-tag Signal Region of the $ b\mu$ channel for the post-fit background, the observed data, and the expected signal with $m_{LQ} = 1$ TeV.

The observed and expected limits on the leptoquark pair production cross-section at 95% CL for $\mathcal{B}=1$ into electrons, shown as a function of $m_{LQ}$ for the $qe$ channel.

The observed and expected limits on the leptoquark pair production cross-section at 95% CL for $\mathcal{B}=1$ into muons, shown as a function of $m_{LQ}$ for the $q\mu$ channel.

The observed and expected limits on the leptoquark pair production cross-section at 95% CL for $\mathcal{B}=1$ into electrons, shown as a function of $m_{LQ}$ for the $ce$ channel.

The observed and expected limits on the leptoquark pair production cross-section at 95% CL for $\mathcal{B}=1$ into muons, shown as a function of $m_{LQ}$ for the $c\mu$ channel.

The observed and expected limits on the leptoquark pair production cross-section at 95% CL for $\mathcal{B}=1$ into electrons, shown as a function of $m_{LQ}$ for the $be$ channel.

The observed and expected limits on the leptoquark pair production cross-section at 95% CL for $\mathcal{B}=1$ into muons, shown as a function of $m_{LQ}$ for the $b\mu$ channel.

The observed and expected limits on the leptoquark branching ratio at 95% CL, shown as a function of $m_{LQ}$ for the $qe$ channel.

The observed and expected limits on the leptoquark branching ratio at 95% CL, shown as a function of $m_{LQ}$ for the $q\mu$ channel.

The observed and expected limits on the leptoquark branching ratio at 95% CL, shown as a function of $m_{LQ}$ for the $ce$ channel.

The observed and expected limits on the leptoquark branching ratio at 95% CL, shown as a function of $m_{LQ}$ for the $c\mu$ channel.

The observed and expected limits on the leptoquark branching ratio at 95% CL, shown as a function of $m_{LQ}$ for the $be$ channel.

The observed and expected limits on the leptoquark branching ratio at 95% CL, shown as a function of $m_{LQ}$ for the $b\mu$ channel.

The signal selection efficiency x acceptance summed over all signal regions, for all masses and LQ decay channels considered.

The observed and expected limits for all masses and LQ decay channels considered.

Cutflow Table in the electron channel, considering signal samples with LQ mass of 1 TeV.

Cutflow Table in the muon channel, considering signal samples with LQ mass of 1 TeV.

Event yields for a simulated $ZH(\rightarrow aa\rightarrow (b\bar{b})(b\bar{b}))$ sample with $m_a = 17.5\,\text{GeV}$. The signal sample is produced with cross section equals to the standard model $pp\to ZH$, i.e. $0.88\,\text{pb}$. Cut 0 corresponds to the initial number of events. Cut 1 requires the single lepton trigger. Cut 2 requires 2 identified leptons. Cut 3 requires the Z-boson mass window. Cut 4 requires 2 reconstructed $a\to b\bar{b}$ candidates. Cut 5a requires 2 identified $a\to b\bar{b}$ candidates in the 1HPC1LPC region. Cut 6a requires the 2 $a\to b\bar{b}$ candidates in the 1HPC1LPC region to be inside the Higgs mass window. Cut 5b requires 2 identified $a\to b\bar{b}$ candidates in the 2HPC region. Cut 6b requires the 2 $a\to b\bar{b}$ candidates in the 2HPC region to be inside the Higgs mass window.

Background yield table for Z+jets, $t\bar{t}$, and rare sources. Observed data yield. Signal $ZH(\rightarrow aa\rightarrow (b\bar{b})(b\bar{b}))$ yield with $m_a = 20\,\text{GeV}$. The signal sample is produced with cross section equals to the standard model $pp\to ZH$, i.e. $0.88\,\text{pb}$. The table includes the yields in two signal regions with leptons consistent with an on-shell Z-boson decay, one with 2 $a\to b\bar{b}$ candidates in the 2HPC region and one with 2 $a\to b\bar{b}$ candidates in the 1HPC1LPC region. The table also includes the yields in four control regions, one with leptons consistent with an on-shell Z-boson decay and 2 $a\to b\bar{b}$ candidates in the Low Purity Category (LPC), and three others where the leptons are not consistent an on-shell Z-boson decay.

The minus $1\sigma$ variation of observed exclusion contour obtained by varying the signal cross section within its uncertainty. The contour is given as a function of the $\it{m}_{\tilde{\chi}^{0}_{1}}$ vs. $\it{m}_{\tilde{t}}$.

The plus $1\sigma$ variation of observed exclusion contour obtained by varying the signal cross section within its uncertainty. The contour is given as a function of the $\it{m}_{\tilde{\chi}^{0}_{1}}$ vs. $\it{m}_{\tilde{t}}$.

The minus $1\sigma$ variation of expected exclusion contour obtained by varying MC statistical uncertainties, detector-related systematic uncertainties, and theoretical uncertainties (excluding signal cross section uncertainties). The contour is given as a function of the $\it{m}_{\tilde{\chi}^{0}_{1}}$ vs. $\it{m}_{\tilde{t}}$.

The plus $1\sigma$ variation of expected exclusion contour obtained by varying MC statistical uncertainties, detector-related systematic uncertainties, and theoretical uncertainties (excluding signal cross section uncertainties). The contour is given as a function of the $\it{m}_{\tilde{\chi}^{0}_{1}}$ vs. $\it{m}_{\tilde{t}}$.

The observed exclusion contour at 95% CL as a function of the $\it{m}_{LQ_{3}^{u}}$ vs. $\mathrm{BR}(\it{m}_{LQ_{3}^{u}}\rightarrow b \tau)$. Points that are within the contours are excluded.

The expected exclusion contour at 95% CL as a function of the $\it{m}_{LQ_{3}^{u}}$ vs. $\mathrm{BR}(\it{m}_{LQ_{3}^{u}}\rightarrow b \tau)$. Points that are within the contours are excluded.

The minus $1\sigma$ variation of observed exclusion contour obtained by varying the signal cross section within its uncertainty. The contour is given as a function of the $\it{m}_{LQ_{3}^{u}}$ vs. $\mathrm{BR}(\it{m}_{LQ_{3}^{u}}\rightarrow b \tau)$

The plus $1\sigma$ variation of observed exclusion contour obtained by varying the signal cross section within its uncertainty. The contour is given as a function of the $\it{m}_{LQ_{3}^{u}}$ vs. $\mathrm{BR}(\it{m}_{LQ_{3}^{u}}\rightarrow b \tau)$

The plus $1\sigma$ variation of expected exclusion contour obtained by varying MC statistical uncertainties, detector-related systematic uncertainties, and theoretical uncertainties (excluding signal cross section uncertainties). The contour is given as a function of the $\it{m}_{LQ_{3}^{u}}$ vs. $\mathrm{BR}(\it{m}_{LQ_{3}^{u}}\rightarrow b \tau)$

The minus $1\sigma$ variation of expected exclusion contour obtained by varying MC statistical uncertainties, detector-related systematic uncertainties, and theoretical uncertainties (excluding signal cross section uncertainties). The contour is given as a function of the $\it{m}_{LQ_{3}^{u}}$ vs. $\mathrm{BR}(\it{m}_{LQ_{3}^{u}}\rightarrow b \tau)$

Model dependent upper limit on the cross section for the $(\tilde{t},\tilde{\chi}^{0}_{1})$ signal grid. The column titled 'Leading Region' stores information on which of the fit regions (SRA-B, SRC or SRD) is the dominant based on the expected CLs values.

Expected model dependent upper limit on the cross section for the $(\tilde{t},\tilde{\chi}^{0}_{1})$ signal grid. The column titled 'Leading Region' stores information on which of the fit regions (SRA-B, SRC or SRD) is the dominant based on the expected CLs values.

Model dependent upper limit on the cross section for the $LQ_{3}^{u}$ signal grid with $\mathrm{BR}(\it{m}_{LQ_{3}^{u}}\rightarrow b \tau))=0$ %. Only the SRA-B fit region is considered in this interpretation.

Expected model dependent upper limit on the cross section for the $LQ_{3}^{u}$ signal grid with $\mathrm{BR}(\it{m}_{LQ_{3}^{u}}\rightarrow b \tau))=0$ %. Only the SRA-B fit region is considered in this interpretation.

The distributions of $S$ in SRA-TW. For each bin yields for the data, total SM prediction and a representative signal point are provided. The SM prediction is provided with the MC statistical uncertainties, labeled 'stat', and the remaining uncertainties, labeled 'syst' that include detector-related systematic uncertainties and theoretical uncertainties. The signal predictions is provided with the MC statistical uncertainties only. The rightmost bin includes overflow events.

The distributions of $\it{m}^{\mathrm{R=1.2}}_{1}$ in SRB-TT. For each bin yields for the data, total SM prediction and a representative signal point are provided. The SM prediction is provided with the MC statistical uncertainties, labeled 'stat', and the remaining uncertainties, labeled 'syst' that include detector-related systematic uncertainties and theoretical uncertainties. The signal predictions is provided with the MC statistical uncertainties only. The rightmost bin includes overflow events.

The distributions of R$_{ISR}$ in SRC signal regions before R$_{ISR}$ cuts are applied. For each bin yields for the data, total SM prediction and a representative signal point are provided. The SM prediction is provided with the MC statistical uncertainties, labeled 'stat', and the remaining uncertainties, labeled 'syst' that include detector-related systematic uncertainties and theoretical uncertainties. The signal predictions is provided with the MC statistical uncertainties only. The rightmost bin includes overflow events.

The distributions of $E^{miss}_{T}/\sqrt{H_{T}}$ in SRD0. For each bin yields for the data, total SM prediction and a representative signal point are provided. The SM prediction is provided with the MC statistical uncertainties, labeled 'stat', and the remaining uncertainties, labeled 'syst' that include detector-related systematic uncertainties and theoretical uncertainties. The signal predictions is provided with the MC statistical uncertainties only. The rightmost bin includes overflow events.

The distributions of $E^{miss}_{T}/\sqrt{H_{T}}$ in SRD1. For each bin yields for the data, total SM prediction and a representative signal point are provided. The SM prediction is provided with the MC statistical uncertainties, labeled 'stat', and the remaining uncertainties, labeled 'syst' that include detector-related systematic uncertainties and theoretical uncertainties. The signal predictions is provided with the MC statistical uncertainties only. The rightmost bin includes overflow events.

The distributions of $E^{miss}_{T}/\sqrt{H_{T}}$ in SRD2. For each bin yields for the data, total SM prediction and a representative signal point are provided. The SM prediction is provided with the MC statistical uncertainties, labeled 'stat', and the remaining uncertainties, labeled 'syst' that include detector-related systematic uncertainties and theoretical uncertainties. The signal predictions is provided with the MC statistical uncertainties only. The rightmost bin includes overflow events.

Cutflow for the reference point $(\it{m}_{\tilde{t}}, \it{m}_{\tilde{\chi}^{0}_{1}})= (1300,1)\ \mathrm{GeV} $ in SRA-TT. The column labelled ''Weighted yield'' shows the results including all correction factors applied to simulation, and is normalised to 139 fb$^{-1}$. A notable exception concerns results in the first row, labelled ''Total'', that corresponds to plain $\sigma \cdot \mathcal{L}$ expected. The ''Derivation skim'' includes the requirements that $H_{T}$, the scalar sum of $p_{T}$ of jets and leptons, $H_{T}>150\ \mathrm{GeV}$ or that a ''baseline'' electron or muon has $p_{T}>20\ \mathrm{GeV}$. The definition of ''baseline'' electron/muons, lepton and $\tau$ vetos are described in the main body of the paper. In total 30000 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 $(\it{m}_{\tilde{t}}, \it{m}_{\tilde{\chi}^{0}_{1}})= (1300,1)\ \mathrm{GeV} $ in SRA-TW. The column labelled ''Weighted yield'' shows the results including all correction factors applied to simulation, and is normalised to 139 fb$^{-1}$. A notable exception concerns results in the first row, labelled ''Total'', that corresponds to plain $\sigma \cdot \mathcal{L}$ expected. The ''Derivation skim'' includes the requirements that $H_{T}$, the scalar sum of $p_{T}$ of jets and leptons, $H_{T}>150\ \mathrm{GeV}$ or that a ''baseline'' electron or muon has $p_{T}>20\ \mathrm{GeV}$. The definition of ''baseline'' electron/muons, lepton and $\tau$ vetos are described in the main body of the paper. In total 30000 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 $(\it{m}_{\tilde{t}}, \it{m}_{\tilde{\chi}^{0}_{1}})= (1300,1)\ \mathrm{GeV} $ in SRA-T0. The column labelled ''Weighted yield'' shows the results including all correction factors applied to simulation, and is normalised to 139 fb$^{-1}$. A notable exception concerns results in the first row, labelled ''Total'', that corresponds to plain $\sigma \cdot \mathcal{L}$ expected. The ''Derivation skim'' includes the requirements that $H_{T}$, the scalar sum of $p_{T}$ of jets and leptons, $H_{T}>150\ \mathrm{GeV}$ or that a ''baseline'' electron or muon has $p_{T}>20\ \mathrm{GeV}$. The definition of ''baseline'' electron/muons, lepton and $\tau$ vetos are described in the main body of the paper. In total 30000 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 $(\it{m}_{\tilde{t}}, \it{m}_{\tilde{\chi}^{0}_{1}})= (700,400)\ \mathrm{GeV} $ in signal regions SRB-TT, SRB-TW and SRB-T0. The regions differ by the last cut applied. The column labelled ''Weighted yield'' shows the results including all correction factors applied to simulation, and is normalised to 139 fb$^{-1}$. A notable exception concerns results in the first row, labelled ''Total'', that corresponds to plain $\sigma \cdot \mathcal{L}$ expected. The ''Derivation skim'' includes the requirements that $H_{T}$, the scalar sum of $p_{T}$ of jets and leptons, $H_{T}>150\ \mathrm{GeV}$ or that a ''baseline'' electron or muon has $p_{T}>20\ \mathrm{GeV}$. The definition of ''baseline'' electron/muons, lepton and $\tau$ vetos are described in the main body of the paper. In total 60000 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 $(\it{m}_{\tilde{t}}, \it{m}_{\tilde{\chi}^{0}_{1}})= (500,327)\ \mathrm{GeV} $ in regions SRC-1, SRC-2, SRC-3, SRC-4 and SRC-5. The regions differ by the last cut applied. The column labelled ''Weighted yield'' shows the results including all correction factors applied to simulation, and is normalised to 139 fb$^{-1}$. A notable exception concerns results in the first row, labelled ''Total'', that corresponds to plain $\sigma \cdot \mathcal{L}$ expected. The ''Derivation skim'' includes the requirements that $H_{T}$, the scalar sum of $p_{T}$ of jets and leptons, $H_{T}>150\ \mathrm{GeV}$ or that a ''baseline'' electron or muon has $p_{T}>20\ \mathrm{GeV}$. The definition of ''baseline'' electron/muons, lepton and $\tau$ vetos are described in the main body of the paper. In total 150000 raw MC events with filter efficiency of 0.384 were generated prior to the specified cuts, with the column ''Unweighted yield'' collecting the numbers after each cut.

Cutflow for the reference point $(\it{m}_{\tilde{t}}, \it{m}_{\tilde{\chi}^{0}_{1}})= (550,500)\ \mathrm{GeV} $ in SRD0. The column labelled ''Weighted yield'' shows the results including all correction factors applied to simulation, and is normalised to 139 fb$^{-1}$. A notable exception concerns results in the first row, labelled ''Total'', that corresponds to plain $\sigma \cdot \mathcal{L}$ expected. The ''Derivation skim'' includes the requirements that $H_{T}$, the scalar sum of $p_{T}$ of jets and leptons, $H_{T}>150\ \mathrm{GeV}$ or that a ''baseline'' electron or muon has $p_{T}>20\ \mathrm{GeV}$. The definition of ''baseline'' electron/muons, lepton and $\tau$ vetos are described in the main body of the paper. In total 90000 raw MC events with filter efficiency of 0.428 were generated prior to the specified cuts, with the column ''Unweighted yield'' collecting the numbers after each cut.

Cutflow for the reference point $(\it{m}_{\tilde{t}}, \it{m}_{\tilde{\chi}^{0}_{1}})= (550,500)\ \mathrm{GeV} $ in SRD1. The column labelled ''Weighted yield'' shows the results including all correction factors applied to simulation, and is normalised to 139 fb$^{-1}$. A notable exception concerns results in the first row, labelled ''Total'', that corresponds to plain $\sigma \cdot \mathcal{L}$ expected. The ''Derivation skim'' includes the requirements that $H_{T}$, the scalar sum of $p_{T}$ of jets and leptons, $H_{T}>150\ \mathrm{GeV}$ or that a ''baseline'' electron or muon has $p_{T}>20\ \mathrm{GeV}$. The definition of ''baseline'' electron/muons, lepton and $\tau$ vetos are described in the main body of the paper. In total 90000 raw MC events with filter efficiency of 0.428 were generated prior to the specified cuts, with the column ''Unweighted yield'' collecting the numbers after each cut.

Cutflow for the reference point $(\it{m}_{\tilde{t}}, \it{m}_{\tilde{\chi}^{0}_{1}})= (550,500)\ \mathrm{GeV} $ in SRD2. The column labelled ''Weighted yield'' shows the results including all correction factors applied to simulation, and is normalised to 139 fb$^{-1}$. A notable exception concerns results in the first row, labelled ''Total'', that corresponds to plain $\sigma \cdot \mathcal{L}$ expected. The ''Derivation skim'' includes the requirements that $H_{T}$, the scalar sum of $p_{T}$ of jets and leptons, $H_{T}>150\ \mathrm{GeV}$ or that a ''baseline'' electron or muon has $p_{T}>20\ \mathrm{GeV}$. The definition of ''baseline'' electron/muons, lepton and $\tau$ vetos are described in the main body of the paper. In total 90000 raw MC events with filter efficiency of 0.428 were generated prior to the specified cuts, with the column ''Unweighted yield'' collecting the numbers after each cut.

Signal acceptance in SRA-TT for simplified $(\tilde{t},\tilde{\chi^{0}_1})$ model. Please mind that the acceptance given in the table is multiplied by factor of $10^{3}$

Signal efficiency in SRA-TT for simplified $(\tilde{t},\tilde{\chi^{0}_1})$ model. Please mind that the efficiency in the table is reported in %.

Signal acceptance in SRA-TW for simplified $(\tilde{t},\tilde{\chi^{0}_1})$ model. Please mind that the acceptance given in the table is multiplied by factor of $10^{3}$

Signal efficiency in SRA-TW for simplified $(\tilde{t},\tilde{\chi^{0}_1})$ model. Please mind that the efficiency in the table is reported in %.

Signal acceptance in SRA-T0 for simplified $(\tilde{t},\tilde{\chi^{0}_1})$ model. Please mind that the acceptance given in the table is multiplied by factor of $10^{3}$

Signal efficiency in SRA-T0 for simplified $(\tilde{t},\tilde{\chi^{0}_1})$ model. Please mind that the efficiency in the table is reported in %.

Signal acceptance in SRB-TT for simplified $(\tilde{t},\tilde{\chi^{0}_1})$ model. Please mind that the acceptance given in the table is multiplied by factor of $10^{3}$

Signal efficiency in SRB-TT for simplified $(\tilde{t},\tilde{\chi^{0}_1})$ model. Please mind that the efficiency in the table is reported in %.

Signal acceptance in SRB-TW for simplified $(\tilde{t},\tilde{\chi^{0}_1})$ model. Please mind that the acceptance given in the table is multiplied by factor of $10^{3}$

Signal efficiency in SRB-TW for simplified $(\tilde{t},\tilde{\chi^{0}_1})$ model. Please mind that the efficiency in the table is reported in %.

Signal acceptance in SRB-T0 for simplified $(\tilde{t},\tilde{\chi^{0}_1})$ model. Please mind that the acceptance given in the table is multiplied by factor of $10^{3}$

Signal efficiency in SRB-T0 for simplified $(\tilde{t},\tilde{\chi^{0}_1})$ model. Please mind that the efficiency in the table is reported in %.

Signal acceptance in SRC1 for simplified $(\tilde{t},\tilde{\chi^{0}_1})$ model. Please mind that the acceptance given in the table is multiplied by factor of $10^{5}$ and the results are given here in the $\it{m}_{\tilde{t}}-\it{m}_{\tilde{\chi}^{0}_{1}}$ plane as opposed to the $\it{m}_{\tilde{t}}-\Delta(\it{m}_{\tilde{\chi}^{0}_{1}},\it{m}_{\tilde{t}})$ one showed in the paper plot.

Signal efficiency in SRC1 for simplified $(\tilde{t},\tilde{\chi^{0}_1})$ model. Please mind that the efficiency in the table is reported in % and the results are given here in the $\it{m}_{\tilde{t}}-\it{m}_{\tilde{\chi}^{0}_{1}}$ plane as opposed to the $\it{m}_{\tilde{t}}-\Delta(\it{m}_{\tilde{\chi}^{0}_{1}},\it{m}_{\tilde{t}})$ one showed in the paper plot.

Signal acceptance in SRC2 for simplified $(\tilde{t},\tilde{\chi^{0}_1})$ model. Please mind that the acceptance given in the table is multiplied by factor of $10^{5}$ and the results are given here in the $\it{m}_{\tilde{t}}-\it{m}_{\tilde{\chi}^{0}_{1}}$ plane as opposed to the $\it{m}_{\tilde{t}}-\Delta(\it{m}_{\tilde{\chi}^{0}_{1}},\it{m}_{\tilde{t}})$ one showed in the paper plot.

Signal efficiency in SRC2 for simplified $(\tilde{t},\tilde{\chi^{0}_1})$ model. Please mind that the efficiency in the table is reported in % and the results are given here in the $\it{m}_{\tilde{t}}-\it{m}_{\tilde{\chi}^{0}_{1}}$ plane as opposed to the $\it{m}_{\tilde{t}}-\Delta(\it{m}_{\tilde{\chi}^{0}_{1}},\it{m}_{\tilde{t}})$ one showed in the paper plot.

Signal acceptance in SRC3 for simplified $(\tilde{t},\tilde{\chi^{0}_1})$ model. Please mind that the acceptance given in the table is multiplied by factor of $10^{5}$ and the results are given here in the $\it{m}_{\tilde{t}}-\it{m}_{\tilde{\chi}^{0}_{1}}$ plane as opposed to the $\it{m}_{\tilde{t}}-\Delta(\it{m}_{\tilde{\chi}^{0}_{1}},\it{m}_{\tilde{t}})$ one showed in the paper plot.

Signal efficiency in SRC3 for simplified $(\tilde{t},\tilde{\chi^{0}_1})$ model. Please mind that the efficiency in the table is reported in % and the results are given here in the $\it{m}_{\tilde{t}}-\it{m}_{\tilde{\chi}^{0}_{1}}$ plane as opposed to the $\it{m}_{\tilde{t}}-\Delta(\it{m}_{\tilde{\chi}^{0}_{1}},\it{m}_{\tilde{t}})$ one showed in the paper plot.

Signal acceptance in SRC4 for simplified $(\tilde{t},\tilde{\chi^{0}_1})$ model. Please mind that the acceptance given in the table is multiplied by factor of $10^{5}$ and the results are given here in the $\it{m}_{\tilde{t}}-\it{m}_{\tilde{\chi}^{0}_{1}}$ plane as opposed to the $\it{m}_{\tilde{t}}-\Delta(\it{m}_{\tilde{\chi}^{0}_{1}},\it{m}_{\tilde{t}})$ one showed in the paper plot.

Signal efficiency in SRC4 for simplified $(\tilde{t},\tilde{\chi^{0}_1})$ model. Please mind that the efficiency in the table is reported in % and the results are given here in the $\it{m}_{\tilde{t}}-\it{m}_{\tilde{\chi}^{0}_{1}}$ plane as opposed to the $\it{m}_{\tilde{t}}-\Delta(\it{m}_{\tilde{\chi}^{0}_{1}},\it{m}_{\tilde{t}})$ plane showed in the paper plot.

Signal acceptance in SRC5 for simplified $(\tilde{t},\tilde{\chi^{0}_1})$ model. Please mind that the acceptance given in the table is multiplied by factor of $10^{5}$ and the results are given here in the $\it{m}_{\tilde{t}}-\it{m}_{\tilde{\chi}^{0}_{1}}$ plane as opposed to the $\it{m}_{\tilde{t}}-\Delta(\it{m}_{\tilde{\chi}^{0}_{1}},\it{m}_{\tilde{t}})$ plane showed in the paper plot.

Signal efficiency in SRC5 for simplified $(\tilde{t},\tilde{\chi^{0}_1})$ model. Please mind that the efficiency in the table is reported in % and the results are given here in the $\it{m}_{\tilde{t}}-\it{m}_{\tilde{\chi}^{0}_{1}}$ plane as opposed to the $\it{m}_{\tilde{t}}-\Delta(\it{m}_{\tilde{\chi}^{0}_{1}},\it{m}_{\tilde{t}})$ one showed in the paper plot.

Signal acceptance in SRD0 for simplified $(\tilde{t},\tilde{\chi^{0}_1})$ model. Please mind that the acceptance given in the table is multiplied by factor of $10^{5}$ and the results are given here in the $\it{m}_{\tilde{t}}-\it{m}_{\tilde{\chi}^{0}_{1}}$ plane as opposed to the $\it{m}_{\tilde{t}}-\Delta(\it{m}_{\tilde{\chi}^{0}_{1}},\it{m}_{\tilde{t}})$ one showed in the paper plot.

Signal efficiency in SRD0 for simplified $(\tilde{t},\tilde{\chi^{0}_1})$ model. Please mind that the efficiency in the table is reported in % and the results are given here in the $\it{m}_{\tilde{t}}-\it{m}_{\tilde{\chi}^{0}_{1}}$ plane as opposed to the $\it{m}_{\tilde{t}}-\Delta(\it{m}_{\tilde{\chi}^{0}_{1}},\it{m}_{\tilde{t}})$ one showed in the paper plot.

Signal acceptance in SRD1 for simplified $(\tilde{t},\tilde{\chi^{0}_1})$ model. Please mind that the acceptance given in the table is multiplied by factor of $10^{5}$ and the results are given here in the $\it{m}_{\tilde{t}}-\it{m}_{\tilde{\chi}^{0}_{1}}$ plane as opposed to the $\it{m}_{\tilde{t}}-\Delta(\it{m}_{\tilde{\chi}^{0}_{1}},\it{m}_{\tilde{t}})$ one showed in the paper plot.

Signal efficiency in SRD1 for simplified $(\tilde{t},\tilde{\chi^{0}_1})$ model. Please mind that the efficiency in the table is reported in % and the results are given here in the $\it{m}_{\tilde{t}}-\it{m}_{\tilde{\chi}^{0}_{1}}$ plane as opposed to the $\it{m}_{\tilde{t}}-\Delta(\it{m}_{\tilde{\chi}^{0}_{1}},\it{m}_{\tilde{t}})$ one showed in the paper plot.

Signal acceptance in SRD2 for simplified $(\tilde{t},\tilde{\chi^{0}_1})$ model. Please mind that the acceptance given in the table is multiplied by factor of $10^{5}$ and the results are given here in the $\it{m}_{\tilde{t}}-\it{m}_{\tilde{\chi}^{0}_{1}}$ plane as opposed to the $\it{m}_{\tilde{t}}-\Delta(\it{m}_{\tilde{\chi}^{0}_{1}},\it{m}_{\tilde{t}})$ one showed in the paper plot.

Signal efficiency in SRD2 for simplified $(\tilde{t},\tilde{\chi^{0}_1})$ model. Please mind that the efficiency in the table is reported in % and the results are given here in the $\it{m}_{\tilde{t}}-\it{m}_{\tilde{\chi}^{0}_{1}}$ plane as opposed to the $\it{m}_{\tilde{t}}-\Delta(\it{m}_{\tilde{\chi}^{0}_{1}},\it{m}_{\tilde{t}})$ one showed in the paper plot.

Observed and predicted mTtot distribution in the b-tag category of 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 combined prediction for A and H bosons with masses of 400, 1000 and 1500 GeV and $\tan\beta$ = 6, 12 and 25 respectively in the mh125 scenario are also provided. The combined prediction for A and H bosons with masses of 1000 and 1500 GeV is scaled by 100 in the paper figure, but not in the HepData table.

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 Higgs boson mass.

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.

The observed 95% CL upper limits on $\tan\beta$ as a function of $m_{A}$ in the $M_{h}^{125}$ scenario. The lowest value of $\tan\beta$ considered for the $M_{h}^{125}$ scenario is 0.5. The points in the region which is called "Not applicable" in the paper figure are kept in the HEPData table. Linear connection is applied in the range of signal mass points from 400 to 1000 GeV in the paper figure. The range of $\tan\beta$ shown in the paper figure and the HEPData is from 1 to 60. No theoretical uncertainty is considered when computing these limits.

The expected 95% CL upper limits on $\tan\beta$ as a function of $m_{A}$ in the $M_{h}^{125}$ scenario. The lowest value of $\tan\beta$ considered for the $M_{h}^{125}$ scenario is 0.5. The points in the region which is called "Not applicable" in the paper figure are kept in the HEPData table. Linear connection is applied in the range of signal mass points from 400 to 1000 GeV in the paper figure. The range of $\tan\beta$ shown in the paper figure and the HEPData is from 1 to 60. No theoretical uncertainty is considered when computing these limits.

The expected 95% CL upper limits with plus one sigma on $\tan\beta$ as a function of $m_{A}$ in the $M_{h}^{125}$ scenario. The lowest value of $\tan\beta$ considered for the $M_{h}^{125}$ scenario is 0.5. The points in the region which is called "Not applicable" in the paper figure are kept in the HEPData table. Linear connection is applied in the range of signal mass points from 400 to 1000 GeV in the paper figure. The range of $\tan\beta$ shown in the paper figure and the HEPData is from 1 to 60. No theoretical uncertainty is considered when computing these limits.

The expected 95% CL upper limits with minus one sigma on $\tan\beta$ as a function of $m_{A}$ in the $M_{h}^{125}$ scenario. The lowest value of $\tan\beta$ considered for the $M_{h}^{125}$ scenario is 0.5. The points in the region which is called "Not applicable" in the paper figure are kept in the HEPData table. Linear connection is applied in the range of signal mass points from 400 to 1000 GeV in the paper figure. The range of $\tan\beta$ shown in the paper figure and the HEPData is from 1 to 60. No theoretical uncertainty is considered when computing these limits.

The expected 95% CL upper limits with plus two sigma on $\tan\beta$ as a function of $m_{A}$ in the $M_{h}^{125}$ scenario. The lowest value of $\tan\beta$ considered for the $M_{h}^{125}$ scenario is 0.5. The points in the region which is called "Not applicable" in the paper figure are kept in the HEPData table. Linear connection is applied in the range of signal mass points from 400 to 1000 GeV in the paper figure. The range of $\tan\beta$ shown in the paper figure and the HEPData is from 1 to 60. No theoretical uncertainty is considered when computing these limits.

The expected 95% CL upper limits with minus two sigma on $\tan\beta$ as a function of $m_{A}$ in the $M_{h}^{125}$ scenario. The lowest value of $\tan\beta$ considered for the $M_{h}^{125}$ scenario is 0.5. The points in the region which is called "Not applicable" in the paper figure are kept in the HEPData table. Linear connection is applied in the range of signal mass points from 400 to 1000 GeV in the paper figure. The range of $\tan\beta$ shown in the paper figure and the HEPData is from 1 to 60. No theoretical uncertainty is considered when computing these limits.

The observed 95% CL upper limits on $\tan\beta$ as a function of $m_{A}$ in the hMSSM scenario. The lowest value of $\tan\beta$ considered by the hMSSM scenario is 0.8 and the highest value of mass is 2 TeV. The points in the region which is called "Not applicable" in the paper figure are kept in the HEPData table. Linear connection is applied in the range of signal mass points from 400 to 1000 GeV in the paper figure. The range of $\tan\beta$ shown in the paper figure and the HEPData is from 1 to 60. No theoretical uncertainty is considered when computing these limits.

The expected 95% CL upper limits on $\tan\beta$ as a function of $m_{A}$ in the hMSSM scenario. The lowest value of $\tan\beta$ considered by the hMSSM scenario is 0.8 and the highest value of mass is 2 TeV. The points in the region which is called "Not applicable" in the paper figure are kept in the HEPData table. Linear connection is applied in the range of signal mass points from 400 to 1000 GeV in the paper figure. The range of $\tan\beta$ shown in the paper figure and the HEPData is from 1 to 60. No theoretical uncertainty is considered when computing these limits.

The expected 95% CL upper limits with plus one sigma on $\tan\beta$ as a function of $m_{A}$ in the hMSSM scenario. The lowest value of $\tan\beta$ considered by the hMSSM scenario is 0.8 and the highest value of mass is 2 TeV. The points in the region which is called "Not applicable" in the paper figure are kept in the HEPData table. Linear connection is applied in the range of signal mass points from 400 to 1000 GeV in the paper figure. The range of $\tan\beta$ shown in the paper figure and the HEPData is from 1 to 60. No theoretical uncertainty is considered when computing these limits.

The expected 95% CL upper limits with minus one sigma on $\tan\beta$ as a function of $m_{A}$ in the hMSSM scenario. The lowest value of $\tan\beta$ considered by the hMSSM scenario is 0.8 and the highest value of mass is 2 TeV. The points in the region which is called "Not applicable" in the paper figure are kept in the HEPData table. Linear connection is applied in the range of signal mass points from 400 to 1000 GeV in the paper figure. The range of $\tan\beta$ shown in the paper figure and the HEPData is from 1 to 60. No theoretical uncertainty is considered when computing these limits.

The expected 95% CL upper limits with plus two sigma on $\tan\beta$ as a function of $m_{A}$ in the hMSSM scenario. The lowest value of $\tan\beta$ considered by the hMSSM scenario is 0.8 and the highest value of mass is 2 TeV. The points in the region which is called "Not applicable" in the paper figure are kept in the HEPData table. Linear connection is applied in the range of signal mass points from 400 to 1000 GeV in the paper figure. The range of $\tan\beta$ shown in the paper figure and the HEPData is from 1 to 60. No theoretical uncertainty is considered when computing these limits.

The expected 95% CL upper limits with minus two sigma on $\tan\beta$ as a function of $m_{A}$ in the hMSSM scenario. The lowest value of $\tan\beta$ considered by the hMSSM scenario is 0.8 and the highest value of mass is 2 TeV. The points in the region which is called "Not applicable" in the paper figure are kept in the HEPData table. Linear connection is applied in the range of signal mass points from 400 to 1000 GeV in the paper figure. The range of $\tan\beta$ shown in the paper figure and the HEPData is from 1 to 60. No theoretical uncertainty is considered when computing these limits.

Acceptance times efficiency for a scalar boson produced by gluon-gluon fusion as a function of the scalar boson mass.

Acceptance times efficiency for a scalar boson produced by b-associated production as a function of the scalar boson mass.

Observed 95% CL upper limits on the scalar boson production cross section times ditau branching fraction as a function of the scalar boson mass and the fraction of the b-associated production. The limits are calculated from a statistical combination of the 1l1tau_h and 2tau_h channels.

Expected 95% CL upper limits on the scalar boson production cross section times ditau branching fraction as a function of the scalar boson mass and the fraction of the b-associated production. The limits are calculated from a statistical combination of the 1l1tau_h and 2tau_h channels.

Observed two dimensional likelihood scan of the gluon-gluon fusion cross section times branching 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 scalar 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. 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. The value of $2\Delta(\mathrm{NLL})$ for 200 GeV signal mass point is shown in the HEPData table.

Observed two dimensional likelihood scan of the gluon-gluon fusion cross section times branching 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 scalar 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. 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. The value of $2\Delta(\mathrm{NLL})$ for 250 GeV signal mass point is shown in the HEPData table.

Observed two dimensional likelihood scan of the gluon-gluon fusion cross section times branching 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 scalar 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. 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. The value of $2\Delta(\mathrm{NLL})$ for 300 GeV signal mass point is shown in the HEPData table.

Observed two dimensional likelihood scan of the gluon-gluon fusion cross section times branching 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 scalar 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. 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. The value of $2\Delta(\mathrm{NLL})$ for 350 GeV signal mass point is shown in the HEPData table.

Observed two dimensional likelihood scan of the gluon-gluon fusion cross section times branching 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 scalar 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. 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. The value of $2\Delta(\mathrm{NLL})$ for 400 GeV signal mass point is shown in the HEPData table.

Observed two dimensional likelihood scan of the gluon-gluon fusion cross section times branching 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 scalar 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. 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. The value of $2\Delta(\mathrm{NLL})$ for 500 GeV signal mass point is shown in the HEPData table.

Observed two dimensional likelihood scan of the gluon-gluon fusion cross section times branching 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 scalar 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. 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. The value of $2\Delta(\mathrm{NLL})$ for 600 GeV signal mass point is shown in the HEPData table.

Observed two dimensional likelihood scan of the gluon-gluon fusion cross section times branching 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 scalar 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. 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. The value of $2\Delta(\mathrm{NLL})$ for 700 GeV signal mass point is shown in the HEPData table.

Observed two dimensional likelihood scan of the gluon-gluon fusion cross section times branching 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 scalar 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. 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. The value of $2\Delta(\mathrm{NLL})$ for 800 GeV signal mass point is shown in the HEPData table.

Observed two dimensional likelihood scan of the gluon-gluon fusion cross section times branching 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 scalar 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. 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. The value of $2\Delta(\mathrm{NLL})$ for 1000 GeV signal mass point is shown in the HEPData table.

Observed two dimensional likelihood scan of the gluon-gluon fusion cross section times branching 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 scalar 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. 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. The value of $2\Delta(\mathrm{NLL})$ for 1200 GeV signal mass point is shown in the HEPData table.

Observed two dimensional likelihood scan of the gluon-gluon fusion cross section times branching 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 scalar 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. 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. The value of $2\Delta(\mathrm{NLL})$ for 1500 GeV signal mass point is shown in the HEPData table.

Observed two dimensional likelihood scan of the gluon-gluon fusion cross section times branching 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 scalar 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. 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. The value of $2\Delta(\mathrm{NLL})$ for 2000 GeV signal mass point is shown in the HEPData table.

Observed two dimensional likelihood scan of the gluon-gluon fusion cross section times branching 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 scalar 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. 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. The value of $2\Delta(\mathrm{NLL})$ for 2500 GeV signal mass point is shown in the HEPData table.

Expected two dimensional likelihood scan of the gluon-gluon fusion cross section times branching 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 scalar 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. 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. The value of $2\Delta(\mathrm{NLL})$ for 200 GeV signal mass point is shown in the HEPData table.

Expected two dimensional likelihood scan of the gluon-gluon fusion cross section times branching 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 scalar 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. 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. The value of $2\Delta(\mathrm{NLL})$ for 250 GeV signal mass point is shown in the HEPData table.

Expected two dimensional likelihood scan of the gluon-gluon fusion cross section times branching 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 scalar 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. 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. The value of $2\Delta(\mathrm{NLL})$ for 300 GeV signal mass point is shown in the HEPData table.

Expected two dimensional likelihood scan of the gluon-gluon fusion cross section times branching 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 scalar 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. 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. The value of $2\Delta(\mathrm{NLL})$ for 350 GeV signal mass point is shown in the HEPData table.

Expected two dimensional likelihood scan of the gluon-gluon fusion cross section times branching 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 scalar 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. 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. The value of $2\Delta(\mathrm{NLL})$ for 400 GeV signal mass point is shown in the HEPData table.

Expected two dimensional likelihood scan of the gluon-gluon fusion cross section times branching 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 scalar 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. 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. The value of $2\Delta(\mathrm{NLL})$ for 500 GeV signal mass point is shown in the HEPData table.

Expected two dimensional likelihood scan of the gluon-gluon fusion cross section times branching 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 scalar 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. 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. The value of $2\Delta(\mathrm{NLL})$ for 600 GeV signal mass point is shown in the HEPData table.

Expected two dimensional likelihood scan of the gluon-gluon fusion cross section times branching 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 scalar 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. 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. The value of $2\Delta(\mathrm{NLL})$ for 700 GeV signal mass point is shown in the HEPData table.

Expected two dimensional likelihood scan of the gluon-gluon fusion cross section times branching 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 scalar 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. 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. The value of $2\Delta(\mathrm{NLL})$ for 800 GeV signal mass point is shown in the HEPData table.

Expected two dimensional likelihood scan of the gluon-gluon fusion cross section times branching 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 scalar 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. 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. The value of $2\Delta(\mathrm{NLL})$ for 1000 GeV signal mass point is shown in the HEPData table.

Expected two dimensional likelihood scan of the gluon-gluon fusion cross section times branching 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 scalar 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. 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. The value of $2\Delta(\mathrm{NLL})$ for 1200 GeV signal mass point is shown in the HEPData table.

Expected two dimensional likelihood scan of the gluon-gluon fusion cross section times branching 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 scalar 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. 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. The value of $2\Delta(\mathrm{NLL})$ for 1500 GeV signal mass point is shown in the HEPData table.

Expected two dimensional likelihood scan of the gluon-gluon fusion cross section times branching 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 scalar 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. 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. The value of $2\Delta(\mathrm{NLL})$ for 2000 GeV signal mass point is shown in the HEPData table.

Expected two dimensional likelihood scan of the gluon-gluon fusion cross section times branching 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 scalar 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. 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. The value of $2\Delta(\mathrm{NLL})$ for 2500 GeV signal mass point is shown in the HEPData table.

The observed 95% CL upper limits on $\tan\beta$ as a function of $m_{A}$ in the hMSSM scenario. The lowest value of $\tan\beta$ considered by the hMSSM scenario is 0.8 and the highest value of mass is 2 TeV. The points in the region which is called "Not applicable" in the paper figure are kept in the HEPData table. Linear connection is applied in the range of signal mass points from 400 to 1000 GeV in the paper figure. The range of $\tan\beta$ shown in the paper figure and the HEPData is from 1 to 60. The theoretical uncertainty of signal cross section is considered.

The expected 95% CL upper limits on $\tan\beta$ as a function of $m_{A}$ in the hMSSM scenario. The lowest value of $\tan\beta$ considered by the hMSSM scenario is 0.8 and the highest value of mass is 2 TeV. The points in the region which is called "Not applicable" in the paper figure are kept in the HEPData table. Linear connection is applied in the range of signal mass points from 400 to 1000 GeV in the paper figure. The range of $\tan\beta$ shown in the paper figure and the HEPData is from 1 to 60. The theoretical uncertainty of signal cross section is considered.

The expected 95% CL upper limits with plus one sigma on $\tan\beta$ as a function of $m_{A}$ in the hMSSM scenario. The lowest value of $\tan\beta$ considered by the hMSSM scenario is 0.8 and the highest value of mass is 2 TeV. The points in the region which is called "Not applicable" in the paper figure are kept in the HEPData table. Linear connection is applied in the range of signal mass points from 400 to 1000 GeV in the paper figure. The range of $\tan\beta$ shown in the paper figure and the HEPData is from 1 to 60. The theoretical uncertainty of signal cross section is considered.

The expected 95% CL upper limits with minus one sigma on $\tan\beta$ as a function of $m_{A}$ in the hMSSM scenario. The lowest value of $\tan\beta$ considered by the hMSSM scenario is 0.8 and the highest value of mass is 2 TeV. The points in the region which is called "Not applicable" in the paper figure are kept in the HEPData table. Linear connection is applied in the range of signal mass points from 400 to 1000 GeV in the paper figure. The range of $\tan\beta$ shown in the paper figure and the HEPData is from 1 to 60. The theoretical uncertainty of signal cross section is considered.

The expected 95% CL upper limits with plus two sigma on $\tan\beta$ as a function of $m_{A}$ in the hMSSM scenario. The lowest value of $\tan\beta$ considered by the hMSSM scenario is 0.8 and the highest value of mass is 2 TeV. The points in the region which is called "Not applicable" in the paper figure are kept in the HEPData table. Linear connection is applied in the range of signal mass points from 400 to 1000 GeV in the paper figure. The range of $\tan\beta$ shown in the paper figure and the HEPData is from 1 to 60. The theoretical uncertainty of signal cross section is considered.

The expected 95% CL upper limits with minus two sigma on $\tan\beta$ as a function of $m_{A}$ in the hMSSM scenario. The lowest value of $\tan\beta$ considered by the hMSSM scenario is 0.8 and the highest value of mass is 2 TeV. The points in the region which is called "Not applicable" in the paper figure are kept in the HEPData table. Linear connection is applied in the range of signal mass points from 400 to 1000 GeV in the paper figure. The range of $\tan\beta$ shown in the paper figure and the HEPData is from 1 to 60. The theoretical uncertainty of signal cross section is considered.

The observed 95% CL upper limits on $\tan\beta$ as a function of $m_{A}$ in the $M_{h}^{125}$ scenario. The lowest value of $\tan\beta$ considered by the $M_{h}^{125}$ scenario is 0.5. The points in the region which is called "Not applicable" in the paper figure are kept in the HEPData table. Linear connection is applied in the range of signal mass points from 400 to 1000 GeV in the paper figure. The range of $\tan\beta$ shown in the paper figure and the HEPData is from 1 to 60. The theoretical uncertainty of signal cross section is considered.

The expected 95% CL upper limits on $\tan\beta$ as a function of $m_{A}$ in the $M_{h}^{125}$ scenario. The lowest value of $\tan\beta$ considered by the $M_{h}^{125}$ scenario is 0.5. The points in the region which is called "Not applicable" in the paper figure are kept in the HEPData table. Linear connection is applied in the range of signal mass points from 400 to 1000 GeV in the paper figure. The range of $\tan\beta$ shown in the paper figure and the HEPData is from 1 to 60. The theoretical uncertainty of signal cross section is considered.

The expected 95% CL upper limits with plus one sigma on $\tan\beta$ as a function of $m_{A}$ in the $M_{h}^{125}$ scenario. The lowest value of $\tan\beta$ considered by the $M_{h}^{125}$ scenario is 0.5. The points in the region which is called "Not applicable" in the paper figure are kept in the HEPData table. Linear connection is applied in the range of signal mass points from 400 to 1000 GeV in the paper figure. The range of $\tan\beta$ shown in the paper figure and the HEPData is from 1 to 60. The theoretical uncertainty of signal cross section is considered.

The expected 95% CL upper limits with minus one sigma on $\tan\beta$ as a function of $m_{A}$ in the $M_{h}^{125}$ scenario. The lowest value of $\tan\beta$ considered by the $M_{h}^{125}$ scenario is 0.5. The points in the region which is called "Not applicable" in the paper figure are kept in the HEPData table. Linear connection is applied in the range of signal mass points from 400 to 1000 GeV in the paper figure. The range of $\tan\beta$ shown in the paper figure and the HEPData is from 1 to 60. The theoretical uncertainty of signal cross section is considered.

The expected 95% CL upper limits with plus two sigma on $\tan\beta$ as a function of $m_{A}$ in the $M_{h}^{125}$ scenario. The lowest value of $\tan\beta$ considered by the $M_{h}^{125}$ scenario is 0.5. The points in the region which is called "Not applicable" in the paper figure are kept in the HEPData table. Linear connection is applied in the range of signal mass points from 400 to 1000 GeV in the paper figure. The range of $\tan\beta$ shown in the paper figure and the HEPData is from 1 to 60. The theoretical uncertainty of signal cross section is considered.

The expected 95% CL upper limits with minus two sigma on $\tan\beta$ as a function of $m_{A}$ in the $M_{h}^{125}$ scenario. The lowest value of $\tan\beta$ considered by the $M_{h}^{125}$ scenario is 0.5. The points in the region which is called "Not applicable" in the paper figure are kept in the HEPData table. Linear connection is applied in the range of signal mass points from 400 to 1000 GeV in the paper figure. The range of $\tan\beta$ shown in the paper figure and the HEPData is from 1 to 60. The theoretical uncertainty of signal cross section is considered.

The observed 95% CL upper limits on $\tan\beta$ as a function of $m_{A}$ in the $M_{h}^{125}(\widetilde{\chi})$ scenario. The lowest value of $\tan\beta$ considered by the $M_{h}^{125}(\widetilde{\chi})$ scenario is 0.5. The points in the region which is called "Not applicable" in the paper figure are kept in the HEPData table. Linear connection is applied in the range of signal mass points from 400 to 1000 GeV in the paper figure. The range of $\tan\beta$ shown in the paper figure and the HEPData is from 1 to 60. The theoretical uncertainty of signal cross section is considered.

The expected 95% CL upper limits on $\tan\beta$ as a function of $m_{A}$ in the $M_{h}^{125}(\widetilde{\chi})$ scenario. The lowest value of $\tan\beta$ considered by the $M_{h}^{125}(\widetilde{\chi})$ scenario is 0.5. The points in the region which is called "Not applicable" in the paper figure are kept in the HEPData table. Linear connection is applied in the range of signal mass points from 400 to 1000 GeV in the paper figure. The range of $\tan\beta$ shown in the paper figure and the HEPData is from 1 to 60. The theoretical uncertainty of signal cross section is considered.

The expected 95% CL upper limits with plus one sigma on $\tan\beta$ as a function of $m_{A}$ in the $M_{h}^{125}(\widetilde{\chi})$ scenario. The lowest value of $\tan\beta$ considered by the $M_{h}^{125}(\widetilde{\chi})$ scenario is 0.5. The points in the region which is called "Not applicable" in the paper figure are kept in the HEPData table. Linear connection is applied in the range of signal mass points from 400 to 1000 GeV in the paper figure. The range of $\tan\beta$ shown in the paper figure and the HEPData is from 1 to 60. The theoretical uncertainty of signal cross section is considered.

The expected 95% CL upper limits with minus one sigma on $\tan\beta$ as a function of $m_{A}$ in the $M_{h}^{125}(\widetilde{\chi})$ scenario. The lowest value of $\tan\beta$ considered by the $M_{h}^{125}(\widetilde{\chi})$ scenario is 0.5. The points in the region which is called "Not applicable" in the paper figure are kept in the HEPData table. Linear connection is applied in the range of signal mass points from 400 to 1000 GeV in the paper figure. The range of $\tan\beta$ shown in the paper figure and the HEPData is from 1 to 60. The theoretical uncertainty of signal cross section is considered.

The expected 95% CL upper limits with plus two sigma on $\tan\beta$ as a function of $m_{A}$ in the $M_{h}^{125}(\widetilde{\chi})$ scenario. The lowest value of $\tan\beta$ considered by the $M_{h}^{125}(\widetilde{\chi})$ scenario is 0.5. The points in the region which is called "Not applicable" in the paper figure are kept in the HEPData table. Linear connection is applied in the range of signal mass points from 400 to 1000 GeV in the paper figure. The range of $\tan\beta$ shown in the paper figure and the HEPData is from 1 to 60. The theoretical uncertainty of signal cross section is considered.

The expected 95% CL upper limits with minus two sigma on $\tan\beta$ as a function of $m_{A}$ in the $M_{h}^{125}(\widetilde{\chi})$ scenario. The lowest value of $\tan\beta$ considered by the $M_{h}^{125}(\widetilde{\chi})$ scenario is 0.5. The points in the region which is called "Not applicable" in the paper figure are kept in the HEPData table. Linear connection is applied in the range of signal mass points from 400 to 1000 GeV in the paper figure. The range of $\tan\beta$ shown in the paper figure and the HEPData is from 1 to 60. The theoretical uncertainty of signal cross section is considered.

The observed 95% CL upper limits on $\tan\beta$ as a function of $m_{A}$ in the $M_{h}^{125}(\widetilde{\tau})$ scenario. The lowest value of $\tan\beta$ considered by the $M_{h}^{125}(\widetilde{\tau})$ scenario is 0.5. The points in the region which is called "Not applicable" in the paper figure are kept in the HEPData table. Linear connection is applied in the range of signal mass points from 400 to 1000 GeV in the paper figure. The range of $\tan\beta$ shown in the paper figure and the HEPData is from 1 to 60. The theoretical uncertainty of signal cross section is considered.

The expected 95% CL upper limits on $\tan\beta$ as a function of $m_{A}$ in the $M_{h}^{125}(\widetilde{\tau})$ scenario. The lowest value of $\tan\beta$ considered by the $M_{h}^{125}(\widetilde{\tau})$ scenario is 0.5. The points in the region which is called "Not applicable" in the paper figure are kept in the HEPData table. Linear connection is applied in the range of signal mass points from 400 to 1000 GeV in the paper figure. The range of $\tan\beta$ shown in the paper figure and the HEPData is from 1 to 60. The theoretical uncertainty of signal cross section is considered.

The expected 95% CL upper limits with plus one sigma on $\tan\beta$ as a function of $m_{A}$ in the $M_{h}^{125}(\widetilde{\tau})$ scenario. The lowest value of $\tan\beta$ considered by the $M_{h}^{125}(\widetilde{\tau})$ scenario is 0.5. The points in the region which is called "Not applicable" in the paper figure are kept in the HEPData table. Linear connection is applied in the range of signal mass points from 400 to 1000 GeV in the paper figure. The range of $\tan\beta$ shown in the paper figure and the HEPData is from 1 to 60. The theoretical uncertainty of signal cross section is considered.

The expected 95% CL upper limits with minus one sigma on $\tan\beta$ as a function of $m_{A}$ in the $M_{h}^{125}(\widetilde{\tau})$ scenario. The lowest value of $\tan\beta$ considered by the $M_{h}^{125}(\widetilde{\tau})$ scenario is 0.5. The points in the region which is called "Not applicable" in the paper figure are kept in the HEPData table. Linear connection is applied in the range of signal mass points from 400 to 1000 GeV in the paper figure. The range of $\tan\beta$ shown in the paper figure and the HEPData is from 1 to 60. The theoretical uncertainty of signal cross section is considered.

The expected 95% CL upper limits with plus two sigma on $\tan\beta$ as a function of $m_{A}$ in the $M_{h}^{125}(\widetilde{\tau})$ scenario. The lowest value of $\tan\beta$ considered by the $M_{h}^{125}(\widetilde{\tau})$ scenario is 0.5. The points in the region which is called "Not applicable" in the paper figure are kept in the HEPData table. Linear connection is applied in the range of signal mass points from 400 to 1000 GeV in the paper figure. The range of $\tan\beta$ shown in the paper figure and the HEPData is from 1 to 60. The theoretical uncertainty of signal cross section is considered.

The expected 95% CL upper limits with minus two sigma on $\tan\beta$ as a function of $m_{A}$ in the $M_{h}^{125}(\widetilde{\tau})$ scenario. The lowest value of $\tan\beta$ considered by the $M_{h}^{125}(\widetilde{\tau})$ scenario is 0.5. The points in the region which is called "Not applicable" in the paper figure are kept in the HEPData table. Linear connection is applied in the range of signal mass points from 400 to 1000 GeV in the paper figure. The range of $\tan\beta$ shown in the paper figure and the HEPData is from 1 to 60. The theoretical uncertainty of signal cross section is considered.

The observed 95% CL upper limits with one sigma on $\tan\beta$ as a function of $m_{A}$ in the $M_{h}^{125}(alignment)$ scenario. The lowest value of $\tan\beta$ considered by the $M_{h}^{125}(alignment)$ scenario is 1.0. The highest value of $\tan\beta$ considered by the $M_{h}^{125}(alignment)$ scenario is 20.0. The points in the region which is called "Not applicable" in the paper figure are kept in the HEPData table. Linear connection is applied in the range of signal mass points from 400 to 1000 GeV in the paper figure. The theoretical uncertainty of signal cross section is considered.

The expected 95% CL upper limits on $\tan\beta$ as a function of $m_{A}$ in the $M_{h}^{125}(alignment)$ scenario. The lowest value of $\tan\beta$ considered by the $M_{h}^{125}(alignment)$ scenario is 1.0. The highest value of $\tan\beta$ considered by the $M_{h}^{125}(alignment)$ scenario is 20.0. The points in the region which is called "Not applicable" in the paper figure are kept in the HEPData table. Linear connection is applied in the range of signal mass points from 400 to 1000 GeV in the paper figure. The theoretical uncertainty of signal cross section is considered.

The expected 95% CL upper limits with plus one sigma on $\tan\beta$ as a function of $m_{A}$ in the $M_{h}^{125}(alignment)$ scenario. The lowest value of $\tan\beta$ considered by the $M_{h}^{125}(alignment)$ scenario is 1.0. The highest value of $\tan\beta$ considered by the $M_{h}^{125}(alignment)$ scenario is 20.0. The points in the region which is called "Not applicable" in the paper figure are kept in the HEPData table. Linear connection is applied in the range of signal mass points from 400 to 1000 GeV in the paper figure. The theoretical uncertainty of signal cross section is considered.

The expected 95% CL upper limits with minus one sigma on $\tan\beta$ as a function of $m_{A}$ in the $M_{h}^{125}(alignment)$ scenario. The lowest value of $\tan\beta$ considered by the $M_{h}^{125}(alignment)$ scenario is 1.0. The highest value of $\tan\beta$ considered by the $M_{h}^{125}(alignment)$ scenario is 20.0. The points in the region which is called "Not applicable" in the paper figure are kept in the HEPData table. Linear connection is applied in the range of signal mass points from 400 to 1000 GeV in the paper figure. The theoretical uncertainty of signal cross section is considered.

The expected 95% CL upper limits with plus two sigma on $\tan\beta$ as a function of $m_{A}$ in the $M_{h}^{125}(alignment)$ scenario. The lowest value of $\tan\beta$ considered by the $M_{h}^{125}(alignment)$ scenario is 1.0. The highest value of $\tan\beta$ considered by the $M_{h}^{125}(alignment)$ scenario is 20.0. The points in the region which is called "Not applicable" in the paper figure are kept in the HEPData table. Linear connection is applied in the range of signal mass points from 400 to 1000 GeV in the paper figure. The theoretical uncertainty of signal cross section is considered.

The expected 95% CL upper limits with minus two sigma on $\tan\beta$ as a function of $m_{A}$ in the $M_{h}^{125}(alignment)$ scenario. The lowest value of $\tan\beta$ considered by the $M_{h}^{125}(alignment)$ scenario is 1.0. The highest value of $\tan\beta$ considered by the $M_{h}^{125}(alignment)$ scenario is 20.0. The points in the region which is called "Not applicable" in the paper figure are kept in the HEPData table. Linear connection is applied in the range of signal mass points from 400 to 1000 GeV in the paper figure. The theoretical uncertainty of signal cross section is considered.

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

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

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

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

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

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

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

Signal acceptance in SR highMass for combined stau final states

Signal acceptance in SR lowMass for combined stau final states

Signal efficiency in SR highMass for combined stau final states

Signal efficiency in SR lowMass for combined stau final states

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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