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.

Cross section upper limit vs m(GLUINO) for SMS model T5ZZ.

Cross section upper limit vs m(GLUINO) for SMS model T5ZZ.

Cross section upper limit vs m(GLUINO) for SMS model T5ZZ.

Cross section vs m(GLUINO) for SMS model T5ZZ.

Cross section vs m(GLUINO) for SMS model T5ZZ.

Cross section vs m(GLUINO) for SMS model T5ZZ.

Number of events in the pTmiss CR, transfer factor, background prediction, and observed yield in each of the six pTmiss bins. The systematic uncertainties in the background prediction include the shape uncertainties in addition to the uncertainty in T. Also listed in the last column is the number of expected signal events and corresponding statistical uncertainties for one example mass point, m(GLUINO) = 1700 GeV. The last pTmiss bin includes overflows.

Pre-fit background covariance matrix $\sigma_{xy}$ for the 6 pTmiss bins.

Pre-fit background correlation matrix $\rho_{xy}$ for the 6 pTmiss bins.

Observed significance (in standard deviations) of the T5ZZ signal model as a function of the gluino mass.

Cut flow for the signal for two example mass points of the T5ZZ signal model. Here the hadronic baseline is defined by $N_{\mathrm{jets}}\ge 2$, $p_T^{\mathrm{miss}}$ > 300 GeV, $H_T$ > 300 GeV, $\Delta\phi$ cuts, and the isolated photon, lepton and track veto.

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.

Track multiplicity $n_{Tracks}$ for preselected DVs in muon-triggered events with at least one muon passing the full selection. Along with the data shown with black markers, the stacked filled histograms represent the background estimates, and predictions for signal scenarios are overlaid with dashed lines. The errors include statistical and systematic uncertainties and are indicated by hatched bands. The DV full selection requirements, $n_{Tracks} \geq 3$ and $m_{DV} > 20$ GeV are visualized with a black arrow.

Invariant mass $m_{DV}$ for the highest-mass preselected DV with at least three associated tracks in MET-triggered events with at least one muon passing the full selection. Along with the data shown with black markers, the stacked filled histograms represent the background estimates, and predictions for signal scenarios are overlaid with dashed lines. The errors include statistical and systematic uncertainties and are indicated by hatched bands. The DV full selection requirements, $n_{Tracks} \geq 3$ and $m_{DV} > 20$ GeV are visualized with a black arrow.

Invariant mass $m_{DV}$ for the highest-mass preselected DV with at least three associated tracks in muon-triggered events with at least one muon passing the full selection. Along with the data shown with black markers, the stacked filled histograms represent the background estimates, and predictions for signal scenarios are overlaid with dashed lines. The errors include statistical and systematic uncertainties and are indicated by hatched bands. The DV full selection requirements, $n_{Tracks} \geq 3$ and $m_{DV} > 20$ GeV are visualized with a black arrow.

The observed event yields in the control, validation and signal regions are shown for the MET Trigger selections, along with the predicted background yields. The bottom panel shows the ratio of observed events to the total background yields. The errors represent the total uncertainty of the backgrounds prediction, including the statistical and systematic uncertainties added in quadrature.

The observed event yields in the control, validation and signal regions are shown for the Muon Trigger selections, along with the predicted background yields. The bottom panel shows the ratio of observed events to the total background yields. The errors represent the total uncertainty of the backgrounds prediction, including the statistical and systematic uncertainties added in quadrature.

Expected exclusion limits at 95% CL on m($\tilde{t}$) as a function of $\tau(\tilde{t})$.

Expected (1 sigma band) exclusion limits at 95% CL on m($\tilde{t}$) as a function of $\tau(\tilde{t})$.

Expected (2 sigma band) exclusion limits at 95% CL on m($\tilde{t}$) as a function of $\tau(\tilde{t})$.

Observed exclusion limits at 95% CL on m($\tilde{t}$) as a function of $\tau(\tilde{t})$.

Observed (+1 sigma) exclusion limits at 95% CL on m($\tilde{t}$) as a function of $\tau(\tilde{t})$.

Observed (-1 sigma) exclusion limits at 95% CL on m($\tilde{t}$) as a function of $\tau(\tilde{t})$.

Exclusion limits on the production cross section as a function of m($\tilde{t}$) are shown for several values of $\tau(\tilde{t})$ along with the nominal signal production cross section and its theoretical uncertainty.

Parameterized event selection efficiencies for the $E_{T}^{miss}$ Trigger SR. The event-level efficiencies for each SR are extracted for all events passing the acceptance of the corresponding SR.

Parameterized event selection efficiencies for the Muon Trigger SR. The event-level efficiencies for each SR are extracted for all events passing the acceptance of the corresponding SR.

Parameterized muon-level reconstruction efficiencies as a function of the muon $p_{T}$ and $d_{0}$. The muon-level efficiencies are extracted using muons passing the muon acceptance criteria.

Parameterized vertex-level reconstruction efficiencies as a function of the radial position of the truth vertex. The efficiency is calculated independent of the muons originating from this truth vertex.

Parameterized vertex-level reconstruction efficiencies as a function of the radial position of the truth vertex. The efficiency is calculated only for truth vertices which have a muon originating from them which is matched to a reconstructed muon.

The $p_{T}$ distribution of all muons originating from LLP decays in the samples used to calculate and validate the efficiencies.

The invariant mass and multiplicity of selected decay products of all truth vertices used in the calculation and validation of the reconstructed efficiencies.

Observed and expected limits at 95% CL on the cross-section of narrow resonance $X$ to $HH$ as a function of the mass $m_{X}$.

Observed and expected limits at 95% CL on the cross-section of broad resonance $X$ to $HH$ as a function of the mass $m_{X}$.

Observed and expected limits at 95% CL on the cross-sections of non-resonant $HH$ production via VBF as a function of the di-vector-boson-di-Higgs-boson coupling modifier $\kappa_{2V}$. The theory prediction of the cross-section as a function of $\kappa_{2V}$ is also shown. The minimum theory cross-section is 1.37 fb at $\kappa_{2V}=1.17$.

Distributions in SR-low of the data and post-fit background prediction for p_T^soft/(p_T^soft + m_eff^3ℓ). The SR-low event selections are applied for each distribution except for the variable shown, where the selection is indicated by a red arrow. The normalization factor for the WZ background is derived from the background-only estimation described in Section 7. The expected distribution for a benchmark signal model is included for comparison. The first (last) bin includes underflow (overflow). The "Top-quark like" category contains the tt̄, Wt, and WW processes while the "Others" category contains backgrounds from triboson production and processes that include a Higgs boson, 3 or more tops, and tops produced in association with W or Z bosons. The bottom panel shows the ratio of the data to the post-fit background prediction. The hatched bands indicate the combined theoretical, experimental, and MC statistical uncertainties.

Distributions in SR-ISR of the data and post-fit background prediction for m_T. The SR-ISR event selections are applied for each distribution except for the variable shown, where the selection is indicated by a red arrow. The normalization factor for the WZ background is derived from the background-only estimation described in Section 7. The expected distribution for a benchmark signal model is included for comparison. The first (last) bin includes underflow (overflow). The "Top-quark like" category contains the tt̄, Wt, and WW processes while the "Others" category contains backgrounds from triboson production and processes that include a Higgs boson, 3 or more tops, and tops produced in association with W or Z bosons. The bottom panel shows the ratio of the data to the post-fit background prediction. The hatched bands indicate the combined theoretical, experimental, and MC statistical uncertainties.

Distributions in SR-ISR of the data and post-fit background prediction for R(E_T^miss,jets). The SR-ISR event selections are applied for each distribution except for the variable shown, where the selection is indicated by a red arrow. The normalization factor for the WZ background is derived from the background-only estimation described in Section 7. The expected distribution for a benchmark signal model is included for comparison. The first (last) bin includes underflow (overflow). The "Top-quark like" category contains the tt̄, Wt, and WW processes while the "Others" category contains backgrounds from triboson production and processes that include a Higgs boson, 3 or more tops, and tops produced in association with W or Z bosons. The bottom panel shows the ratio of the data to the post-fit background prediction. The hatched bands indicate the combined theoretical, experimental, and MC statistical uncertainties.

Distributions in SR-ISR of the data and post-fit background prediction for p_T^soft. The SR-ISR event selections are applied for each distribution except for the variable shown, where the selection is indicated by a red arrow. The normalization factor for the WZ background is derived from the background-only estimation described in Section 7. The expected distribution for a benchmark signal model is included for comparison. The first (last) bin includes underflow (overflow). The "Top-quark like" category contains the tt̄, Wt, and WW processes while the "Others" category contains backgrounds from triboson production and processes that include a Higgs boson, 3 or more tops, and tops produced in association with W or Z bosons. The bottom panel shows the ratio of the data to the post-fit background prediction. The hatched bands indicate the combined theoretical, experimental, and MC statistical uncertainties.

Distributions in SR-ISR of the data and post-fit background prediction for p_T^jets. The SR-ISR event selections are applied for each distribution except for the variable shown, where the selection is indicated by a red arrow. The normalization factor for the WZ background is derived from the background-only estimation described in Section 7. The expected distribution for a benchmark signal model is included for comparison. The first (last) bin includes underflow (overflow). The "Top-quark like" category contains the tt̄, Wt, and WW processes while the "Others" category contains backgrounds from triboson production and processes that include a Higgs boson, 3 or more tops, and tops produced in association with W or Z bosons. The bottom panel shows the ratio of the data to the post-fit background prediction. The hatched bands indicate the combined theoretical, experimental, and MC statistical uncertainties.

Observed exclusion contour on C1N2 production assuming on-shell $W/Z$ decays as a function of the C1/N2 and N1 masses, and derived from the combined fit of low-mass and ISR regions.

Expected exclusion contour on C1N2 production assuming on-shell $W/Z$ decays as a function of the C1/N2 and N1 masses, and derived from the combined fit of low-mass and ISR regions.

Plus 1$\sigma$ uncertainty, varying the signal cross section within its uncertainty, on the observed exclusion contour on C1N2 production assuming on-shell $W/Z$ decays as a function of the C1/N2 and N1 masses, and derived from the combined fit of low-mass and ISR regions.

Minus 1$\sigma$ uncertainty, varying the signal cross section within its uncertainty, on the observed exclusion contour on C1N2 production assuming on-shell $W/Z$ decays as a function of the C1/N2 and N1 masses, and derived from the combined fit of low-mass and ISR regions.

Plus 1$\sigma$ uncertainty, due to uncertainties in the background prediction and experimental uncertainties affecting the signal, on the expected exclusion contour on C1N2 production assuming on-shell $W/Z$ decays as a function of the C1/N2 and N1 masses, and derived from the combined fit of low-mass and ISR regions.

Minus 1$\sigma$ uncertainty, due to uncertainties in the background prediction and experimental uncertainties affecting the signal, on the expected exclusion contour on C1N2 production assuming on-shell $W/Z$ decays as a function of the C1/N2 and N1 masses, and derived from the combined fit of low-mass and ISR regions.

Upper limits on observed wino-bino simplified model signal cross section $\sigma_\text{obs}^\text{95}$.

Upper limits on expected wino-bino simplified model signal cross section $\sigma_\text{exp}^\text{95}$.

Signal acceptance in SR-low, for signals with $m(\widetilde{\chi}^{\pm}_{1}/\widetilde{\chi}^{0}_{2}) - m\widetilde{\chi}^{0}_{1} \geq 100$ GeV.

Signal efficiency in SR-low, for signals with $m(\widetilde{\chi}^{\pm}_{1}/\widetilde{\chi}^{0}_{2}) - m\widetilde{\chi}^{0}_{1} \geq 100$ GeV.

Signal acceptance in SR-ISR, for signals with $m(\widetilde{\chi}^{\pm}_{1}/\widetilde{\chi}^{0}_{2}) - m\widetilde{\chi}^{0}_{1} \geq 100$ GeV.

Signal efficiency in SR-ISR, for signals with $m(\widetilde{\chi}^{\pm}_{1}/\widetilde{\chi}^{0}_{2}) - m\widetilde{\chi}^{0}_{1} \geq 100$ GeV.

Signal acceptance in SR-low, for signals with $m(\widetilde{\chi}^{\pm}_{1}/\widetilde{\chi}^{0}_{2}) - m\widetilde{\chi}^{0}_{1} < 100$ GeV.

Signal efficiency in SR-low, for signals with $m(\widetilde{\chi}^{\pm}_{1}/\widetilde{\chi}^{0}_{2}) - m\widetilde{\chi}^{0}_{1} < 100$ GeV.

Signal acceptance in SR-ISR, for signals with $m(\widetilde{\chi}^{\pm}_{1}/\widetilde{\chi}^{0}_{2}) - m\widetilde{\chi}^{0}_{1} < 100$ GeV.

Signal efficiency in SR-ISR, for signals with $m(\widetilde{\chi}^{\pm}_{1}/\widetilde{\chi}^{0}_{2}) - m\widetilde{\chi}^{0}_{1} < 100$ GeV.

The observed and expected yields after the background-only fit in the SRs. The normalization factors of the $WZ$ sample for the low-mass and ISR regions are different and are treated separately in the combined fit. \The "Top-quark like" category contains the tt&#772;, Wt, and WW processes while the "Others" category contains backgrounds from triboson production and processes that include a Higgs boson, 3 or more tops, and tops produced in association with W or Z bosons. Combined statistical and systematic uncertainties are presented. The individual uncertainties can be correlated and do not necessarily add in quadrature to equal the total background uncertainty.

Summary of the expected background and data yields in $\text{SR-low}$ and $\text{SR-ISR}$. The second and third columns show the data and total expected background with systematic uncertainties. The fourth column gives the model-independent upper limits at 95\% CL on the visible cross section ($\sigma_\text{vis}$). The fifth and sixth columns give the visible number of observed ($S^{95}_\text{obs}$) and expected ($S^{95}_\text{exp}$) events of a generic beyond-the-SM process, where uncertainties on $S^{95}_\text{exp}$ reflect the $\pm 1 \sigma$ uncertainties on the background estimation. The last column shows the discovery $p$-value and Gaussian significance $Z$ assuming no signal.

Upper limits on observed (expected) wino-bino simplified model signal cross section $\sigma_\text{obs(exp)}^\text{95}$.

Full list of event selections and MC generator-weighted yields and in $\text{SR-ISR}$ for the main $WZ$ background and a representative $\tilde{\chi}^{\pm}_{1}\tilde{\chi}^{0}_{2}$ signal point of mass 200 GeV and mass splitting $\Delta m = 100$ GeV normalized to 139 fb$^{-1}$. 40000 events were generated.

Full list of event selections and MC generator-weighted yields and in $\text{SR-low}$ for the main $WZ$ background and a representative $\tilde{\chi}^{\pm}_{1}\tilde{\chi}^{0}_{2}$ signal point of mass 200 GeV and mass splitting $\Delta m = 100$ GeV normalized to 139 fb$^{-1}$. 40000 events were generated.

Observed $CL_S$ limits on $\sigma \times BR$ for $m_{\Phi} = 200-400$ GeV.

Observed $CL_S$ limits on $\sigma \times BR$ for $m_{\Phi} = 600-1000$ GeV.

Combined limits from this analysis (ID) and the CR and MS analyses for $m_{H} = 125$ GeV, $m_s = 15$ GeV.

Combined limits from this analysis (ID) and the CR and MS analyses for $m_{H} = 125$ GeV, $m_s = 25$ GeV.

Combined limits from this analysis (ID) and the CR and MS analyses for $m_{H} = 125$ GeV, $m_s = 40$ GeV.

Combined limits from this analysis (ID) and the CR and MS analyses for $m_{H} = 125$ GeV, $m_s = 55$ GeV.

Combined limits from this analysis (ID) and the CR and MS analyses for $m_{\Phi} = 200$ GeV, $m_s = 25$ GeV.

Combined limits from this analysis (ID) and the CR and MS analyses for $m_{\Phi} = 200$ GeV, $m_s = 50$ GeV.

Comparison of the IDVx reconstruction and selection efficiency for a Higgs with a mass of 125 GeV decaying to an LLP with a mass of 8 GeV, using only standard tracking versus using standard and large radius tracking for all ID vertices in the signal MC sample.

Comparison of the IDVx reconstruction and selection efficiency for a Higgs with a mass of 125 GeV decaying to an LLP with a mass of 8 GeV, using only standard tracking versus using standard and large radius tracking for ID vertices passing the full IDVx selection criteria.

Comparison of the IDVx reconstruction and selection efficiency for a $\Phi$ with a mass of 1000 GeV decaying to an LLP with a mass of 150 GeV, using only standard tracking versus using standard and large radius tracking for all ID vertices in the signal MC sample.

Comparison of the IDVx reconstruction and selection efficiency for a $\Phi$ with a mass of 1000 GeV decaying to an LLP with a mass of 150 GeV, using only standard tracking versus using standard and large radius tracking for ID vertices passing the full IDVx selection criteria.

The IDVx selection efficiency as a function of long-lived particle decay $z$ position for MC signal samples with a 125 GeV Higgs boson decaying to LLPs with masses of 8, 25, and 55 GeV. The efficiency for the 55 GeV LLP sample is shown both with and without the material veto applied.

The IDVx selection efficiency as a function of long-lived particle decay $z$ position for MC signal samples with mediators of masses 200, 400, and 600 GeV decaying to LLPs with a mass of 50 GeV. The efficiency for the 200 GeV mediator sample is shown both with and without the material veto applied.

The impact of the IDVx selections on the selection efficiency for the MC signal sample with a Higgs with a mass of 125 GeV decaying to an LLP with a mass of 8 GeV.

The impact of the IDVx selections on the selection efficiency for the MC signal sample with a Higgs with a mass of 125 GeV decaying to an LLP with a mass of 8 GeV.

The impact of the IDVx selections on the selection efficiency for the MC signal sample with a Higgs with a mass of 125 GeV decaying to an LLP with a mass of 8 GeV.

The impact of the IDVx selections on the selection efficiency for the MC signal sample with a $\Phi$ with a mass of 400 GeV decaying to an LLP with a mass of 50 GeV.

The impact of the IDVx selections on the selection efficiency for the MC signal sample with a $\Phi$ with a mass of 400 GeV decaying to an LLP with a mass of 50 GeV.

The impact of the IDVx selections on the selection efficiency for the MC signal sample with a $\Phi$ with a mass of 400 GeV decaying to an LLP with a mass of 50 GeV.

The impact of the IDVx selections on the selection efficiency for the MC signal sample with a $\Phi$ with a mass of 1000 GeV decaying to an LLP with a mass of 150 GeV.

The impact of the IDVx selections on the selection efficiency for the MC signal sample with a $\Phi$ with a mass of 1000 GeV decaying to an LLP with a mass of 150 GeV.

The impact of the IDVx selections on the selection efficiency for the MC signal sample with a $\Phi$ with a mass of 1000 GeV decaying to an LLP with a mass of 150 GeV.

$CL_S$ limits on $B_{H\rightarrow ss}$ for $m_{H} = 125$ GeV, $m_s = 8$ GeV.

$CL_S$ limits on $B_{H\rightarrow ss}$ for $m_{H} = 125$ GeV, $m_s = 15$ GeV.

$CL_S$ limits on $B_{H\rightarrow ss}$ for $m_{H} = 125$ GeV, $m_s = 25$ GeV.

$CL_S$ limits on $B_{H\rightarrow ss}$ for $m_{H} = 125$ GeV, $m_s = 40$ GeV.

$CL_S$ limits on $B_{H\rightarrow ss}$ for $m_{H} = 125$ GeV, $m_s = 55$ GeV.

$CL_S$ limits on $\sigma \times B_{\Phi\rightarrow ss}$ for $m_{\Phi} = 200$ GeV, $m_s = 8$ GeV.

$CL_S$ limits on $\sigma \times B_{\Phi\rightarrow ss}$ for $m_{\Phi} = 200$ GeV, $m_s = 25$ GeV.

$CL_S$ limits on $\sigma \times B_{\Phi\rightarrow ss}$ for $m_{\Phi} = 200$ GeV, $m_s = 50$ GeV.

$CL_S$ limits on $\sigma \times B_{\Phi\rightarrow ss}$ for $m_{\Phi} = 400$ GeV, $m_s = 50$ GeV.

$CL_S$ limits on $\sigma \times B_{\Phi\rightarrow ss}$ for $m_{\Phi} = 400$ GeV, $m_s = 100$ GeV.

$CL_S$ limits on $\sigma \times B_{\Phi\rightarrow ss}$ for $m_{\Phi} = 600$ GeV, $m_s = 50$ GeV.

$CL_S$ limits on $\sigma \times B_{\Phi\rightarrow ss}$ for $m_{\Phi} = 600$ GeV, $m_s = 150$ GeV.

$CL_S$ limits on $\sigma \times B_{\Phi\rightarrow ss}$ for $m_{\Phi} = 1000$ GeV, $m_s = 50$ GeV.

$CL_S$ limits on $\sigma \times B_{\Phi\rightarrow ss}$ for $m_{\Phi} = 1000$ GeV, $m_s = 150$ GeV.

$CL_S$ limits on $\sigma \times B_{\Phi\rightarrow ss}$ for $m_{\Phi} = 1000$ GeV, $m_s = 400$ GeV.

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.

The rho distributions for the mixed-flavour category events corresponding to a signal hypothesis with mH = 261 GeV and mA = 150 GeV.

The rho distributions for the same-flavour category events corresponding to a signal hypothesis with mH = 442 GeV and mA = 193 GeV. The signal is normalised to its theoretical cross section.

The rho distributions for the mixed-flavour category events corresponding to a signal hypothesis with mH = 442 GeV and mA = 193 GeV.

Expected and observed 95% CL upper limits on the product of the production cross section and branching fraction for H(A) -> ZA(H) -> bbbar as a function of mA and mH. The limits are computed using the asymptotic CLs method, combining the ee and μμ channels.

Expected and observed 95% CL upper limits on the signal strength for the Type-II 2HDM benchmark (tan(beta)=1.5, cos(beta-alpha)=0.01) as a function of mA and mH . The limits are computed using the asymptotic CLs method, combining the ee and μμ channels.

Expected and observed 95% CL upper limits on the signal strength for the Type-II 2HDM benchmark (mH = 379 GeV and m A = 172 GeV) as a function of tan(beta) and cos(beta-alpha). The limits are computed using the asymptotic CLs method, combining the ee and μμ channels.

P-value for H(A) -> ZA(H) -> bbbar production as a function of mA and mH. The limits are computed using the asymptotic CLs method, combining the ee and μμ channels.

The post-fit $m_{CT}$ distribution is shown in the validation region VR-offLM after all the selection requirements are applied other than the $m_{CT}$ cut. The stacked histograms show the expected SM backgrounds. The hatched bands represent the sum in quadrature of systematic and statistical uncertainties of the total SM background. The red line with arrow indicates the $m_{CT}$ cut used in SR selection. The first and the last bin include the underflow and overflow events (where present), respectively.

The post-fit $m_{CT}$ distribution is shown in the validation region VR-offMM after all the selection requirements are applied other than the $m_{CT}$ cut. The stacked histograms show the expected SM backgrounds. The hatched bands represent the sum in quadrature of systematic and statistical uncertainties of the total SM background. The red line with arrow indicates the $m_{CT}$ cut used in SR selection. The first and the last bin include the underflow and overflow events (where present), respectively.

The post-fit $m_{CT}$ distribution is shown in the validation region VR-offHM after all the selection requirements are applied other than the $m_{CT}$ cut. The stacked histograms show the expected SM backgrounds. The hatched bands represent the sum in quadrature of systematic and statistical uncertainties of the total SM background. The red line with arrow indicates the $m_{CT}$ cut used in SR selection. The first and the last bin include the underflow and overflow events (where present), respectively.

The post-fit $m_{CT}$ distribution for SR-HM. The stacked histograms show the expected SM backgrounds. The hatched bands represent the sum in quadrature of systematic and statistical uncertainties of the total SM background. For illustration, the distribution of the SUSY reference points are also shown as dashed lines. The overflow events, where present, are included in the last bin.

The post-fit $m_{CT}$ distribution for SR-MM. The stacked histograms show the expected SM backgrounds. The hatched bands represent the sum in quadrature of systematic and statistical uncertainties of the total SM background. For illustration, the distribution of the SUSY reference points are also shown as dashed lines. The overflow events, where present, are included in the last bin.

The post-fit $m_{CT}$ distribution for SR-LM. The stacked histograms show the expected SM backgrounds. The hatched bands represent the sum in quadrature of systematic and statistical uncertainties of the total SM background. For illustration, the distribution of the SUSY reference points are also shown as dashed lines. The overflow events, where present, are included in the last bin.

The post-fit $m_{bb}$ distribution is shown in the signal region SR-HM after all the selection requirements are applied other than the $m_{bb}$ cut. The stacked histograms show the expected SM backgrounds. The hatched bands represent the sum in quadrature of systematic and statistical uncertainties of the total SM background. For illustration, the distribution of the SUSY reference points are also shown as dashed lines. The red line with arrow indicates the $m_{bb}$ cut used in SR selection.The overflow events, where present, are included in the last bin.

The post-fit $m_{bb}$ distribution is shown in the signal region SR-MM after all the selection requirements are applied other than the $m_{bb}$ cut. The hatched bands represent the sum in quadrature of systematic and statistical uncertainties of the total SM background. For illustration, the distribution of the SUSY reference points are also shown as dashed lines. The red line with arrow indicates the $m_{bb}$ cut used in SR selection. The overflow events, where present, are included in the last bin.

The post-fit $m_{bb}$ distribution is shown in the signal region SR-LM after all the selection requirements are applied other than the $m_{bb}$ cut. The hatched bands represent the sum in quadrature of systematic and statistical uncertainties of the total SM background. For illustration, the distribution of the SUSY reference points are also shown as dashed lines. The red line with arrow indicates the $m_{bb}$ cut used in SR selection. The overflow events, where present, are included in the last bin.

The observed exclusion for simplified models with $\tilde\chi^\pm_1 \tilde\chi^0_2 \rightarrow Wh\tilde\chi^0_1\tilde\chi^0_1, W \rightarrow l\nu, h \rightarrow b\bar{b}$ production. Experimental and theoretical systematic uncertainties are applied to background and signal samples and illustrated by the yellow band and the red dotted contour lines, respectively. The red dotted lines indicate the $\pm$ 1 standard-deviation variation on the observed exclusion limit due to theoretical uncertainties in the signal cross-section.

The observed exclusion up limit for simplified models with $\tilde\chi^\pm_1 \tilde\chi^0_2 \rightarrow Wh\tilde\chi^0_1\tilde\chi^0_1, W \rightarrow l\nu, h \rightarrow b\bar{b}$ production. The red dotted lines indicate the $\pm 1 \sigma$ on the observed exclusion limit due to the theoretical uncertainties in the signal cross-section.

The observed exclusion down limit for simplified models with $\tilde\chi^\pm_1 \tilde\chi^0_2 \rightarrow Wh\tilde\chi^0_1\tilde\chi^0_1, W \rightarrow l\nu, h \rightarrow b\bar{b}$ production. The red dotted lines indicate the $\pm 1 \sigma$ on the observed exclusion limit due to the theoretical uncertainties in the signal cross-section.

The expected exclusion for simplified models with $\tilde\chi^\pm_1 \tilde\chi^0_2 \rightarrow Wh\tilde\chi^0_1\tilde\chi^0_1, W \rightarrow l\nu, h \rightarrow b\bar{b}$ production. Experimental and theoretical systematic uncertainties are applied to background and signal samples and illustrated by the yellow band and the red dotted contour lines, respectively. The red dotted lines indicate the $\pm$ 1 standard-deviation variation on the observed exclusion limit due to theoretical uncertainties in the signal cross-section.

Upper limits on the cross sections for simplified models with $\tilde\chi^\pm_1 \tilde\chi^0_2 \rightarrow Wh\tilde\chi^0_1\tilde\chi^0_1, W \rightarrow l\nu, h \rightarrow b\bar{b}$ production.

Signal acceptance in SR-LM for simplified models with $\tilde\chi^\pm_1 \tilde\chi^0_2 \rightarrow Wh\tilde\chi^0_1\tilde\chi^0_1, W \rightarrow l\nu, h \rightarrow b\bar{b}$ production.

Signal acceptance in SR-LM low $m_{CT}$ for simplified models with $\tilde\chi^\pm_1 \tilde\chi^0_2 \rightarrow Wh\tilde\chi^0_1\tilde\chi^0_1, W \rightarrow l\nu, h \rightarrow b\bar{b}$ production.

Signal acceptance in SR-LM med. $m_{CT}$ for simplified models with $\tilde\chi^\pm_1 \tilde\chi^0_2 \rightarrow Wh\tilde\chi^0_1\tilde\chi^0_1, W \rightarrow l\nu, h \rightarrow b\bar{b}$ production.

Signal acceptance in SR-LM high $m_{CT}$ for simplified models with $\tilde\chi^\pm_1 \tilde\chi^0_2 \rightarrow Wh\tilde\chi^0_1\tilde\chi^0_1, W \rightarrow l\nu, h \rightarrow b\bar{b}$ production.

Signal acceptance in SR-MM for simplified models with $\tilde\chi^\pm_1 \tilde\chi^0_2 \rightarrow Wh\tilde\chi^0_1\tilde\chi^0_1, W \rightarrow l\nu, h \rightarrow b\bar{b}$ production.

Signal acceptance in SR-MM low $m_{CT}$ for simplified models with $\tilde\chi^\pm_1 \tilde\chi^0_2 \rightarrow Wh\tilde\chi^0_1\tilde\chi^0_1, W \rightarrow l\nu, h \rightarrow b\bar{b}$ production.

Signal acceptance in SR-MM med. $m_{CT}$ for simplified models with $\tilde\chi^\pm_1 \tilde\chi^0_2 \rightarrow Wh\tilde\chi^0_1\tilde\chi^0_1, W \rightarrow l\nu, h \rightarrow b\bar{b}$ production.

Signal acceptance in SR-MM high $m_{CT}$ for simplified models with $\tilde\chi^\pm_1 \tilde\chi^0_2 \rightarrow Wh\tilde\chi^0_1\tilde\chi^0_1, W \rightarrow l\nu, h \rightarrow b\bar{b}$ production.

Signal acceptance in SR-HM for simplified models with $\tilde\chi^\pm_1 \tilde\chi^0_2 \rightarrow Wh\tilde\chi^0_1\tilde\chi^0_1, W \rightarrow l\nu, h \rightarrow b\bar{b}$ production. 1lb\bar{b}$ production

Signal acceptance in SR-HM low $m_{CT}$ for simplified models with $\tilde\chi^\pm_1 \tilde\chi^0_2 \rightarrow Wh\tilde\chi^0_1\tilde\chi^0_1, W \rightarrow l\nu, h \rightarrow b\bar{b}$ production.

Signal acceptance in SR-HM med. $m_{CT}$ for simplified models with $\tilde\chi^\pm_1 \tilde\chi^0_2 \rightarrow Wh\tilde\chi^0_1\tilde\chi^0_1, W \rightarrow l\nu, h \rightarrow b\bar{b}$ production.

Signal acceptance in SR-HM high $m_{CT}$ for simplified models with $\tilde\chi^\pm_1 \tilde\chi^0_2 \rightarrow Wh\tilde\chi^0_1\tilde\chi^0_1, W \rightarrow l\nu, h \rightarrow b\bar{b}$ production.

Signal efficiency in SR-LM for simplified models with $\tilde\chi^\pm_1 \tilde\chi^0_2 \rightarrow Wh\tilde\chi^0_1\tilde\chi^0_1, W \rightarrow l\nu, h \rightarrow b\bar{b}$ production.

Signal efficiency in SR-LM low $m_{CT}$ for simplified models with $\tilde\chi^\pm_1 \tilde\chi^0_2 \rightarrow Wh\tilde\chi^0_1\tilde\chi^0_1, W \rightarrow l\nu, h \rightarrow b\bar{b}$ production.

Signal efficiency in SR-LM med. $m_{CT}$ for simplified models with $\tilde\chi^\pm_1 \tilde\chi^0_2 \rightarrow Wh\tilde\chi^0_1\tilde\chi^0_1, W \rightarrow l\nu, h \rightarrow b\bar{b}$ production.

Signal efficiency in SR-LM high $m_{CT}$ for simplified models with $\tilde\chi^\pm_1 \tilde\chi^0_2 \rightarrow Wh\tilde\chi^0_1\tilde\chi^0_1, W \rightarrow l\nu, h \rightarrow b\bar{b}$ production.

Signal efficiency in SR-MM for simplified models with $\tilde\chi^\pm_1 \tilde\chi^0_2 \rightarrow Wh\tilde\chi^0_1\tilde\chi^0_1, W \rightarrow l\nu, h \rightarrow b\bar{b}$ production.

Signal efficiency in SR-MM low $m_{CT}$ for simplified models with $\tilde\chi^\pm_1 \tilde\chi^0_2 \rightarrow Wh\tilde\chi^0_1\tilde\chi^0_1, W \rightarrow l\nu, h \rightarrow b\bar{b}$ production.

Signal efficiency in SR-MM med. $m_{CT}$ for simplified models with $\tilde\chi^\pm_1 \tilde\chi^0_2 \rightarrow Wh\tilde\chi^0_1\tilde\chi^0_1, W \rightarrow l\nu, h \rightarrow b\bar{b}$ production.

Signal efficiency in SR-MM high $m_{CT}$ for simplified models with $\tilde\chi^\pm_1 \tilde\chi^0_2 \rightarrow Wh\tilde\chi^0_1\tilde\chi^0_1, W \rightarrow l\nu, h \rightarrow b\bar{b}$ production.

Signal efficiency in SR-HM for simplified models with $\tilde\chi^\pm_1 \tilde\chi^0_2 \rightarrow Wh\tilde\chi^0_1\tilde\chi^0_1, W \rightarrow l\nu, h \rightarrow b\bar{b}$ production.

Signal efficiency in SR-HM low $m_{CT}$ for simplified models with $\tilde\chi^\pm_1 \tilde\chi^0_2 \rightarrow Wh\tilde\chi^0_1\tilde\chi^0_1, W \rightarrow l\nu, h \rightarrow b\bar{b}$ production.

Signal efficiency in SR-HM med. $m_{CT}$ for simplified models with $\tilde\chi^\pm_1 \tilde\chi^0_2 \rightarrow Wh\tilde\chi^0_1\tilde\chi^0_1, W \rightarrow l\nu, h \rightarrow b\bar{b}$ production.

Signal efficiency in SR-HM high $m_{CT}$ for simplified models with $\tilde\chi^\pm_1 \tilde\chi^0_2 \rightarrow Wh\tilde\chi^0_1\tilde\chi^0_1, W \rightarrow l\nu, h \rightarrow b\bar{b}$ production.

Event selection cutflow for a representative signal sample for the SR-LM low $m_{CT}$. The masses of next-lightest-neutralinos and LSPs are reported. While the first row of the table reports the total raw MC events produced, all subsequent rows show weighted events. Only statistical uncertainties are shown. Samples are produced with generator filters which selects $h\rightarrow b\bar{b}$ and $W\rightarrow\ell\nu$ decays.

Event selection cutflow for a representative signal sample for the SR-LM med. $m_{CT}$. The masses of next-lightest-neutralinos and LSPs are reported. While the first row of the table reports the total raw MC events produced, all subsequent rows show weighted events. Only statistical uncertainties are shown. Samples are produced with generator filters which selects $h\rightarrow b\bar{b}$ and $W\rightarrow\ell\nu$ decays.

Event selection cutflow for a representative signal sample for the SR-LM high $m_{CT}$. The masses of next-lightest-neutralinos and LSPs are reported. While the first row of the table reports the total raw MC events produced, all subsequent rows show weighted events. Only statistical uncertainties are shown. Samples are produced with generator filters which selects $h\rightarrow b\bar{b}$ and $W\rightarrow\ell\nu$ decays.

Event selection cutflow for a representative signal sample for the SR-MM low $m_{CT}$. The masses of next-lightest-neutralinos and LSPs are reported. While the first row of the table reports the total raw MC events produced, all subsequent rows show weighted events. Only statistical uncertainties are shown. Samples are produced with generator filters which selects $h\rightarrow b\bar{b}$ and $W\rightarrow\ell\nu$ decays.

Event selection cutflow for a representative signal sample for the SR-MM med. $m_{CT}$. The masses of next-lightest-neutralinos and LSPs are reported. While the first row of the table reports the total raw MC events produced, all subsequent rows show weighted events. Only statistical uncertainties are shown. Samples are produced with generator filters which selects $h\rightarrow b\bar{b}$ and $W\rightarrow\ell\nu$ decays.

Event selection cutflow for a representative signal sample for the SR-MM high $m_{CT}$. The masses of next-lightest-neutralinos and LSPs are reported. While the first row of the table reports the total raw MC events produced, all subsequent rows show weighted events. Only statistical uncertainties are shown. Samples are produced with generator filters which selects $h\rightarrow b\bar{b}$ and $W\rightarrow\ell\nu$ decays.

Event selection cutflow for a representative signal sample for the SR-HM low $m_{CT}$. The masses of next-lightest-neutralinos and LSPs are reported. While the first row of the table reports the total raw MC events produced, all subsequent rows show weighted events. Only statistical uncertainties are shown. Samples are produced with generator filters which selects $h\rightarrow b\bar{b}$ and $W\rightarrow\ell\nu$ decays.

Event selection cutflow for a representative signal sample for the SR-HM med. $m_{CT}$. The masses of next-lightest-neutralinos and LSPs are reported. While the first row of the table reports the total raw MC events produced, all subsequent rows show weighted events. Only statistical uncertainties are shown. Samples are produced with generator filters which selects $h\rightarrow b\bar{b}$ and $W\rightarrow\ell\nu$ decays.

Event selection cutflow for a representative signal sample for the SR-HM high $m_{CT}$. The masses of next-lightest-neutralinos and LSPs are reported. While the first row of the table reports the total raw MC events produced, all subsequent rows show weighted events. Only statistical uncertainties are shown. Samples are produced with generator filters which selects $h\rightarrow b\bar{b}$ and $W\rightarrow\ell\nu$ decays.

Event selection cutflow for a representative signal sample for the discovery SR-LM. The masses of next-lightest-neutralinos and LSPs are reported. While the first row of the table reports the total raw MC events produced, all subsequent rows show weighted events. Only statistical uncertainties are shown. Samples are produced with generator filters which selects $h\rightarrow b\bar{b}$ and $W\rightarrow\ell\nu$ decays.

Event selection cutflow for a representative signal sample for the discovery SR-MM. The masses of next-lightest-neutralinos and LSPs are reported. While the first row of the table reports the total raw MC events produced, all subsequent rows show weighted events. Only statistical uncertainties are shown. Samples are produced with generator filters which selects $h\rightarrow b\bar{b}$ and $W\rightarrow\ell\nu$ decays.

Event selection cutflow for a representative signal sample for the discovery SR-HM. The masses of next-lightest-neutralinos and LSPs are reported. While the first row of the table reports the total raw MC events produced, all subsequent rows show weighted events. Only statistical uncertainties are shown. Samples are produced with generator filters which selects $h\rightarrow b\bar{b}$ and $W\rightarrow\ell\nu$ decays.