Distribution of the missing transverse energy in CRTop. The Z+HF and $t\bar{t}$ scale factors, described in the text, have been applied to the simulated samples. The distribution labeled "Signal" is for the model with ($m_a$, $m_{\tilde{\chi}_{1}^{0}}$, $m_{\tilde{\chi}_{2}^{0}}$) = (45 GeV, 10 GeV, 80 GeV).

Distribution of the dijet invariant mass in the signal region, shown together with the parameterized background model (labelled "Bkg Model"). For reference, the MC prediction for the SM background is also shown (labelled "SM MC"). The Z+HF and $t\bar{t}$ scale factors, described in the text, have been applied to the simulated samples. The signal region is defined to have dijet invariant mass > 20 GeV. The distribution labeled "Signal" is for the model with ($m_a$, $m_{\tilde{\chi}_{1}^{0}}$, $m_{\tilde{\chi}_{2}^{0}}$) = (45 GeV, 10 GeV, 80 GeV).

Upper limits at 95% CL on the $pp \rightarrow ZH$ cross section times the branching ratio for $Z \rightarrow \ell^{+}\ell^{-}$ (where $\ell = e, \mu \;\mathrm{or}\; \tau$) and $H \rightarrow \tilde{\chi}_{2}^{0}\tilde{\chi}_{1}^{0} \rightarrow a \tilde{\chi}_{1}^{0}\tilde{\chi}_{1}^{0} \rightarrow b\bar{b} \tilde{\chi}_{1}^{0}\tilde{\chi}_{1}^{0}$ as a function of $m_{a}$ for $m_{\tilde{\chi}_{1}^{0}} = 10$ GeV and $m_{\tilde{\chi}_{2}^{0}} = 65$ GeV in the NMSSM scenario described in the text. All branching ratios in the Higgs boson decay chain after the decay $H \rightarrow \tilde{\chi}_{2}^{0}\tilde{\chi}_{1}^{0}$ are set to 100%.

Upper limits at 95% CL on the $pp \rightarrow ZH$ cross section times the branching ratio for $Z \rightarrow \ell^{+}\ell^{-}$ (where $\ell = e, \mu \;\mathrm{or}\; \tau$) and $H \rightarrow \tilde{\chi}_{2}^{0}\tilde{\chi}_{1}^{0} \rightarrow a \tilde{\chi}_{1}^{0}\tilde{\chi}_{1}^{0} \rightarrow b\bar{b} \tilde{\chi}_{1}^{0}\tilde{\chi}_{1}^{0}$ as a function of $m_{a}$ for $m_{\tilde{\chi}_{1}^{0}} = 10$ GeV and $m_{\tilde{\chi}_{2}^{0}} = 80$ GeV in the NMSSM scenario described in the text. All branching ratios in the Higgs boson decay chain after the decay $H \rightarrow \tilde{\chi}_{2}^{0}\tilde{\chi}_{1}^{0}$ are set to 100%.

Upper limits at 95% CL on the $pp \rightarrow ZH$ cross section times the branching ratio for $Z \rightarrow \ell^{+}\ell^{-}$ (where $\ell = e, \mu\; \mathrm{or}\; \tau$) and $H \rightarrow \tilde{\chi}_{2}^{0}\tilde{\chi}_{1}^{0} \rightarrow a \tilde{\chi}_{1}^{0}\tilde{\chi}_{1}^{0} \rightarrow b\bar{b} \tilde{\chi}_{1}^{0}\tilde{\chi}_{1}^{0}$ as a function of $m_{a}$ for $m_{\tilde{\chi}_{1}^{0}} = 10$ GeV and $m_{\tilde{\chi}_{2}^{0}} = 95$ GeV in the NMSSM scenario described in the text. All branching ratios in the Higgs boson decay chain after the decay $H \rightarrow \tilde{\chi}_{2}^{0}\tilde{\chi}_{1}^{0}$ are set to 100%.

Upper limits at 95% CL on the $pp \rightarrow ZH$ cross section times the branching ratio for $Z \rightarrow \ell^{+}\ell^{-}$ (where $\ell = e, \mu\; \mathrm{or}\; \tau$) and $H \rightarrow \tilde{\chi}_{2}^{0}\tilde{\chi}_{1}^{0} \rightarrow a \tilde{\chi}_{1}^{0}\tilde{\chi}_{1}^{0} \rightarrow b\bar{b} \tilde{\chi}_{1}^{0}\tilde{\chi}_{1}^{0}$ as a function of $m_{a}$ for $m_{\tilde{\chi}_{1}^{0}} = 10$ GeV and $m_{\tilde{\chi}_{2}^{0}} = 110$ GeV in the NMSSM scenario described in the text. All branching ratios in the Higgs boson decay chain after the decay $H \rightarrow \tilde{\chi}_{2}^{0}\tilde{\chi}_{1}^{0}$ are set to 100%.

Upper limits at 95% CL on the $pp \rightarrow ZH$ cross section times the branching ratio for $Z \rightarrow \ell^{+}\ell^{-}$ (where $\ell = e, \mu\; \mathrm{or}\; \tau$) and $H \rightarrow \tilde{\chi}_{2}^{0}\tilde{\chi}_{1}^{0} \rightarrow a \tilde{\chi}_{1}^{0}\tilde{\chi}_{1}^{0} \rightarrow b\bar{b} \tilde{\chi}_{1}^{0}\tilde{\chi}_{1}^{0}$ as a function of $m_{a}$ for $m_{\tilde{\chi}_{1}^{0}} = 20$ GeV and $m_{\tilde{\chi}_{2}^{0}} = 80$ GeV in the NMSSM scenario described in the text. All branching ratios in the Higgs boson decay chain after the decay $H \rightarrow \tilde{\chi}_{2}^{0}\tilde{\chi}_{1}^{0}$ are set to 100%.

Upper limits at 95% CL on the $pp \rightarrow ZH$ cross section times the branching ratio for $Z \rightarrow \ell^{+}\ell^{-}$ (where $\ell = e, \mu \;\mathrm{or}\; \tau$) and $H \rightarrow \tilde{\chi}_{2}^{0}\tilde{\chi}_{1}^{0} \rightarrow a \tilde{\chi}_{1}^{0}\tilde{\chi}_{1}^{0} \rightarrow b\bar{b} \tilde{\chi}_{1}^{0}\tilde{\chi}_{1}^{0}$ as a function of $m_{a}$ for $m_{\tilde{\chi}_{1}^{0}} = 30$ GeV and $m_{\tilde{\chi}_{2}^{0}} = 80$ GeV in the NMSSM scenario described in the text. All branching ratios in the Higgs boson decay chain after the decay $H \rightarrow \tilde{\chi}_{2}^{0}\tilde{\chi}_{1}^{0}$ are set to 100%.

Upper limits at 95% CL on the branching ratio $H \rightarrow \tilde{\chi}_{2}^{0}\tilde{\chi}_{1}^{0} \rightarrow a \tilde{\chi}_{1}^{0}\tilde{\chi}_{1}^{0} \rightarrow b\bar{b} \tilde{\chi}_{1}^{0}\tilde{\chi}_{1}^{0}$ as a function of $m_{a}$ for $m_{\tilde{\chi}_{1}^{0}} = 10$ GeV and $m_{\tilde{\chi}_{2}^{0}} = 65$ GeV in the NMSSM scenario described in the text. The SM $ZH$ cross section is assumed.

Upper limits at 95% CL on the branching ratio $H \rightarrow \tilde{\chi}_{2}^{0}\tilde{\chi}_{1}^{0} \rightarrow a \tilde{\chi}_{1}^{0}\tilde{\chi}_{1}^{0} \rightarrow b\bar{b} \tilde{\chi}_{1}^{0}\tilde{\chi}_{1}^{0}$ as a function of $m_{a}$ for $m_{\tilde{\chi}_{1}^{0}} = 10$ GeV and $m_{\tilde{\chi}_{2}^{0}} = 80$ GeV in the NMSSM scenario described in the text. The SM $ZH$ cross section is assumed.

Upper limits at 95% CL on the branching ratio $H \rightarrow \tilde{\chi}_{2}^{0}\tilde{\chi}_{1}^{0} \rightarrow a \tilde{\chi}_{1}^{0}\tilde{\chi}_{1}^{0} \rightarrow b\bar{b} \tilde{\chi}_{1}^{0}\tilde{\chi}_{1}^{0}$ as a function of $m_{a}$ for $m_{\tilde{\chi}_{1}^{0}} = 10$ GeV and $m_{\tilde{\chi}_{2}^{0}} = 95$ GeV in the NMSSM scenario described in the text. The SM $ZH$ cross section is assumed.

Upper limits at 95% CL on the branching ratio $H \rightarrow \tilde{\chi}_{2}^{0}\tilde{\chi}_{1}^{0} \rightarrow a \tilde{\chi}_{1}^{0}\tilde{\chi}_{1}^{0} \rightarrow b\bar{b} \tilde{\chi}_{1}^{0}\tilde{\chi}_{1}^{0}$ as a function of $m_{a}$ for $m_{\tilde{\chi}_{1}^{0}} = 10$ GeV and $m_{\tilde{\chi}_{2}^{0}} = 110$ GeV in the NMSSM scenario described in the text. The SM $ZH$ cross section is assumed.

Upper limits at 95% CL on the branching ratio $H \rightarrow \tilde{\chi}_{2}^{0}\tilde{\chi}_{1}^{0} \rightarrow a \tilde{\chi}_{1}^{0}\tilde{\chi}_{1}^{0} \rightarrow b\bar{b} \tilde{\chi}_{1}^{0}\tilde{\chi}_{1}^{0}$ as a function of $m_{a}$ for $m_{\tilde{\chi}_{1}^{0}} = 20$ GeV and $m_{\tilde{\chi}_{2}^{0}} = 80$ GeV in the NMSSM scenario described in the text. The SM $ZH$ cross section is assumed.

Upper limits at 95% CL on the branching ratio $H \rightarrow \tilde{\chi}_{2}^{0}\tilde{\chi}_{1}^{0} \rightarrow a \tilde{\chi}_{1}^{0}\tilde{\chi}_{1}^{0} \rightarrow b\bar{b} \tilde{\chi}_{1}^{0}\tilde{\chi}_{1}^{0}$ as a function of $m_{a}$ for $m_{\tilde{\chi}_{1}^{0}} = 30$ GeV and $m_{\tilde{\chi}_{2}^{0}} = 80$ GeV in the NMSSM scenario described in the text. The SM $ZH$ cross section is assumed.

Unweighted and weighted number of events after each stage of selection for the NMSSM scenario with $pp \rightarrow ZH$, $Z \rightarrow \ell^{+}\ell^{-}$, $H \rightarrow \tilde{\chi}_{2}^{0}\tilde{\chi}_{1}^{0} \rightarrow a \tilde{\chi}_{1}^{0}\tilde{\chi}_{1}^{0} \rightarrow b\bar{b} \tilde{\chi}_{1}^{0}\tilde{\chi}_{1}^{0}$ where $(m_{a}, m_{\tilde{\chi}_{1}^{0}}, m_{\tilde{\chi}_{2}^{0}}) = (45,10,80)$ GeV and the Z boson decaying to $e^{+}e^{-}, \;\mu^{+}\mu^{-} \;\mathrm{or}\; \tau^{+}\tau^{-}$. All branching ratios in the Higgs boson decay chain are set to 100%. The weighted number of events corresponds to an integrated luminosity of $139 \;\mathrm{fb}^{-1}$. The "skimming selection" required either a single electron or muon with $p_{T} > 100\;$ GeV or a pair of electrons or muons each with $p_{T} > 20\;$ GeV. The preselection requirements include the trigger, absence of a bad jet or bad muon, and two leptons, where a lepton is either an electron or a muon.

Acceptance and efficiency of this analysis for the signal models considered in this paper. The signal is an NMSSM scenario with $pp \rightarrow ZH$, $Z \rightarrow \ell^{+}\ell^{-}$, $H \rightarrow \tilde{\chi}_{2}^{0}\tilde{\chi}_{1}^{0} \rightarrow a \tilde{\chi}_{1}^{0}\tilde{\chi}_{1}^{0} \rightarrow b\bar{b} \tilde{\chi}_{1}^{0}\tilde{\chi}_{1}^{0}$ where the values of $(m_{a}, m_{\tilde{\chi}_{1}^{0}}, m_{\tilde{\chi}_{2}^{0}})$ are varied. The Z boson decays to $e^{+}e^{-}, \mu^{+}\mu^{-} \;\mathrm{or}\; \tau^{+}\tau^{-}$. All branching ratios in the Higgs boson decay chain are set to 100%. The product of acceptance times efficiency is defined by the fraction of simulated events that pass all the selection criteria of this analysis. The acceptance is defined by the fraction of simulated events that pass all the selection criteria as applied to Monte Carlo truth-level quantities. The efficiency is then defined as the acceptance times efficiency, divided by the acceptance.

Acceptance and efficiency of this analysis for the signal models considered in this paper. The signal is an NMSSM scenario with $pp \rightarrow ZH$, $Z \rightarrow \ell^{+}\ell^{-}$, $H \rightarrow \tilde{\chi}_{2}^{0}\tilde{\chi}_{1}^{0} \rightarrow a \tilde{\chi}_{1}^{0}\tilde{\chi}_{1}^{0} \rightarrow b\bar{b} \tilde{\chi}_{1}^{0}\tilde{\chi}_{1}^{0}$ where the values of $(m_{a}, m_{\tilde{\chi}_{1}^{0}}, m_{\tilde{\chi}_{2}^{0}})$ are varied. The Z boson decays to $e^{+}e^{-}, \mu^{+}\mu^{-} \;\mathrm{or}\; \tau^{+}\tau^{-}$. All branching ratios in the Higgs boson decay chain are set to 100%. The product of acceptance times efficiency is defined by the fraction of simulated events that pass all the selection criteria of this analysis. The acceptance is defined by the fraction of simulated events that pass all the selection criteria as applied to Monte Carlo truth-level quantities. The efficiency is then defined as the acceptance times efficiency, divided by the acceptance.

Post-fit $m_\text{T}(\gamma, E_\text{T}^\text{miss})$ distribution in the inclusive signal region for the dark-photon search with the 125 GeV mass $\mathcal{B}(H \to \gamma \gamma_\text{d})$ signal normalization set to zero. A $H \to \gamma \gamma_\text{d}$ decay signal is shown for two different mass hypotheses, 125 GeV and 500 GeV, and scaled to a $\mathcal{B}(H \to \gamma \gamma_\text{d})$ of 2% and 1%, respectively. Events with $m_\text{T}(\gamma, E_\text{T}^\text{miss})$ larger than the rightmost bin boundary are added to that bin.

The 95% CL upper limit on the Higgs boson production cross-section times branching ratio to $\gamma \gamma_\text{d}$ is shown for different VBF-produced scalar-mediator-mass hypotheses in the NWA. The theoretically predicted cross-section of a Higgs boson produced via VBF and with the $\mathcal{B}(H \to \gamma \gamma_\text{d}) =$ 5% is superimposed on the $\pm 1\sigma$ and $\pm 2\sigma$ NNLO QCD + NLO EW uncertainty band of the expected production cross-section limit.

Post-fit $m_\text{jj}$ distribution in the inclusive signal region. The Higgs boson invisible decay signal is scaled to a $\mathcal{B}_\text{inv}$ of 37%. Events with $m_\text{jj}$ larger than the rightmost bin boundary are added to that bin.

Post-fit $m_\text{jj}$ distribution in the one-lepton control region $W_{\ell \nu}^\gamma$ CR. Events with $m_\text{jj}$ larger than the rightmost bin boundary are added to that bin.

Post-fit $m_\text{T}$ distribution in the one lepton control region. Events with $m_\text{T}$ larger than the rightmost bin boundary are added to that bin.

Post-fit photon centrality distribution in the zero lepton signal plus control region with the $\mathcal{B}_\text{inv}$ signal normalization set to zero in the fit.

Post-fit photon $E_\text{T}$ distribution in the zero lepton signal region with the $\mathcal{B}_\text{inv}$ signal normalization set to zero in the fit.

Post-fit photon centrality distribution in the zero lepton signal plus control region resulting from the fit to the $m_\text{jj}$ distribution for EW $Z \gamma + \text{jets}$. The post-fit uncertainties include statistical, experimental, and theory contributions.

Post-fit photon $E_\text{T}$ distribution in the zero lepton signal region resulting from the fit to the $m_\text{jj}$ distribution for EW $Z \gamma + \text{jets}$. The post-fit uncertainties include statistical, experimental, and theory contributions.

Post-fit DNN output score distribution in the one lepton control region.

Yields for the EW $Z \gamma + \text{jets}$ process are shown after each selection along with relative and absolute signal acceptance efficiencies.

Yields for the 125 GeV Higgs boson with $\mathcal{B}_\text{inv.} =$ 1 signal produced by the vector boson fusion process in association with a final state photon are shown after each selection along with relative and absolute signal acceptance efficiencies.

Yields for the 125 GeV Higgs boson with $\mathcal{B}(H \to \gamma \gamma_\text{d}) =$ 1 signal produced by the vector boson fusion process are shown after each selection along with relative and absolute signal acceptance efficiencies.

Expected +- 2 sigma 95% CL exclusion limit for the Zprime-2HDM model.

Observed 95% CL exclusion limit for the 2HDM+a model ggF production.

Expected 95% CL exclusion limit for the 2HDM+a model ggF production.

Expected +- 1 sigma 95% CL exclusion limit for the 2HDM+a model ggF production.

Expected +- 2 sigma 95% CL exclusion limit for the 2HDM+a model ggF production.

Observed 95% CL exclusion limit for the 2HDM+a model bbA production.

Expected 95% CL exclusion limit for the 2HDM+a model bbA production.

Expected +- 1 sigma 95% CL exclusion limit for the 2HDM+a model bbA production.

Expected +- 2 sigma 95% CL exclusion limit for the 2HDM+a model bbA production.

Observed 95% CL exclusion limit for the Zprime-2HDM model with the benchmark used in arXiv:1707.01302.

Expected 95% CL exclusion limit for the Zprime-2HDM model with the benchmark used in arXiv:1707.01302.

Expected +- 1 sigma 95% CL exclusion limit for the Zprime-2HDM model with the benchmark used in arXiv:1707.01302.

Expected +- 2 sigma 95% CL exclusion limit for the Zprime-2HDM model with the benchmark used in arXiv:1707.01302.

Expected and observed upper limits at 95% CL on cross-section for Zprime-2HDM model.

Expected and observed upper limits at 95% CL on cross-section for ggF producton in the 2HDM+a model.

Expected and observed upper limits at 95% CL on cross-section for bbA producton in the 2HDM+a model.

Model-independent upper limits on the visible cross-section $σ_{vis, $h(\bar{b})+DM} ≡ σ_{h+DM} \times B(h \to b\bar{b}) \times \mathcal{A} \times \epsilon$ in the different signal regions.

Theory cross-section for Zprime-2HDM model.

Theory cross-section for bbA production in the 2HDM+a model.

Theory cross-section for ggF production in the 2HDM+a model.

Distribution of Higgs boson candidate mass in 2b region with MET=150-200 GeV.

Distribution of Higgs boson candidate mass in 2b region with MET=200-350 GeV.

Distribution of Higgs boson candidate mass in 2b region with MET=350-500 GeV.

Distribution of Higgs boson candidate mass in 2b region with MET=500-750 GeV.

Distribution of Higgs boson candidate mass in 2b region with MET > 750 GeV.

Distribution of Higgs boson candidate mass in 3b region with MET=150-200 GeV.

Distribution of Higgs boson candidate mass in 3b region with MET=200-350 GeV.

Distribution of Higgs boson candidate mass in 3b region with MET=350-500 GeV.

Distribution of Higgs boson candidate mass in 3b region with MET > 500 GeV.

Yields in 1-lepton control region.

Yields in 2-lepton control region.

MET distribution in 2b region of the 0-lepton channel.

MET distribution in 3b region of the 0-lepton channel.

Expected signal yields after certain selection cuts in 2b region with MET=150-200 GeV.

Expected signal yields after certain selection cuts in 2b region with MET=200-350 GeV.

Expected signal yields after certain selection cuts in 2b region with MET=350-500 GeV.

Expected signal yields after certain selection cuts in 2b region with MET=500-750 GeV.

Expected signal yields after certain selection cuts in 2b region with MET > 750 GeV.

Expected signal yields after certain selection cuts in 3b region with MET=150-200 GeV.

Expected signal yields after certain selection cuts in 3b region with MET=200-350 GeV.

Expected signal yields after certain selection cuts in 3b region with MET=350-500 GeV.

Expected signal yields after certain selection cuts in 3b region with MET > 500 GeV.

Acceptance times efficiency for bbA production in the 2HDM+a model - 2b region with MET=150-200 GeV.

Acceptance times efficiency for bbA production in the 2HDM+a model - 2b region with MET=200-350 GeV.

Acceptance times efficiency for bbA production in the 2HDM+a model - 2b region with MET=350-500 GeV.

Acceptance times efficiency for bbA production in the 2HDM+a model - 2b region with MET=500-750 GeV.

Acceptance times efficiency for bbA production in the 2HDM+a model - 2b region with MET > 750 GeV.

Acceptance times efficiency for bbA production in the 2HDM+a model - 3b region with MET=150-200 GeV.

Acceptance times efficiency for bbA production in the 2HDM+a model - 3b region with MET=200-350 GeV.

Acceptance times efficiency for bbA production in the 2HDM+a model - 3b region with MET=350-500 GeV.

Acceptance times efficiency for bbA production in the 2HDM+a model - 3b region with MET>500 GeV.

Acceptance times efficiency for ggF production in the 2HDM+a model - 2b region with MET=150-200 GeV.

Acceptance times efficiency for ggF production in the 2HDM+a model - 2b region with MET=200-350 GeV.

Acceptance times efficiency for ggF production in the 2HDM+a model - 2b region with MET=350-500 GeV.

Acceptance times efficiency for ggF production in the 2HDM+a model - 2b region with MET=500-750 GeV.

Acceptance times efficiency for ggF production in the 2HDM+a model - 2b region with MET > 750 GeV.

Acceptance times efficiency for ggF production in the 2HDM+a model - 3b region with MET=150-200 GeV.

Acceptance times efficiency for ggF production in the 2HDM+a model - 3b region with MET=200-350 GeV.

Acceptance times efficiency for ggF production in the 2HDM+a model - 3b region with MET=350-500 GeV.

Acceptance times efficiency for ggF production in the 2HDM+a model - 3b region with MET > 500 GeV.

Acceptance times efficiency for ggF production in the Zprime-2HDM model - 2b region with MET=150-200 GeV.

Acceptance times efficiency for ggF production in the Zprime-2HDM model - 2b region with MET=200-350 GeV.

Acceptance times efficiency for ggF production in the Zprime-2HDM model - 2b region with MET=350-500 GeV.

Acceptance times efficiency for ggF production in the Zprime-2HDM model - 2b region with MET=500-750 GeV.

Acceptance times efficiency for ggF production in the Zprime-2HDM model - 2b region with MET > 750 GeV.

Acceptance times efficiency for ggF production in the Zprime-2HDM model - 3b region with MET=150-200 GeV.

Acceptance times efficiency for ggF production in the Zprime-2HDM model - 3b region with MET=200-350 GeV.

Acceptance times efficiency for ggF production in the Zprime-2HDM model - 3b region with MET=350-500 GeV.

Acceptance times efficiency for ggF production in the Zprime-2HDM model - 3b region with MET > 500 GeV.

The $m_{\mathrm{eff}}$ distribution in SR0$_{\mathrm{bveto}}^{\mathrm{loose}}$ and SR0$_{\mathrm{bveto}}^{\mathrm{tight}}$ for events passing the signal region requirements except the $m_{\mathrm{eff}}$ requirement. Distributions for data, the estimated SM backgrounds after the background-only fit, and an example SUSY scenario are shown. "Other" is the sum of the $tWZ$, $t\bar{t}WW$, $t\bar{t} ZZ$, $t\bar{t} WH$, $t\bar{t} HH$, $t\bar{t} tW$, and $t\bar{t}t\bar{t}$ backgrounds. The last bin captures the overflow events. The lower panel shows the ratio of the observed data to the expected SM background yield in each bin. Both the statistical and systematic uncertainties in the SM background are included in the shaded band. The red arrows indicate the $m_{\mathrm{eff}}$ selections in the signal regions.

The $m_{\mathrm{eff}}$ distribution in SR1$_{\mathrm{bveto}}^{\mathrm{loose}}$ and SR1$_{\mathrm{bveto}}^{\mathrm{tight}}$ for events passing the signal region requirements except the $m_{\mathrm{eff}}$ requirement. Distributions for data, the estimated SM backgrounds after the background-only fit, and an example SUSY scenario are shown. "Other" is the sum of the $tWZ$, $t\bar{t}WW$, $t\bar{t} ZZ$, $t\bar{t} WH$, $t\bar{t} HH$, $t\bar{t} tW$, and $t\bar{t}t\bar{t}$ backgrounds. The last bin captures the overflow events. The lower panel shows the ratio of the observed data to the expected SM background yield in each bin. Both the statistical and systematic uncertainties in the SM background are included in the shaded band. The red arrows indicate the $m_{\mathrm{eff}}$ selections in the signal regions.

The $m_{\mathrm{eff}}$ distribution in SR2$_{\mathrm{bveto}}^{\mathrm{loose}}$ and SR2$_{\mathrm{bveto}}^{\mathrm{tight}}$ for events passing the signal region requirements except the $m_{\mathrm{eff}}$ requirement. Distributions for data, the estimated SM backgrounds after the background-only fit, and an example SUSY scenario are shown. "Other" is the sum of the $tWZ$, $t\bar{t}WW$, $t\bar{t} ZZ$, $t\bar{t} WH$, $t\bar{t} HH$, $t\bar{t} tW$, and $t\bar{t}t\bar{t}$ backgrounds. The last bin captures the overflow events. The lower panel shows the ratio of the observed data to the expected SM background yield in each bin. Both the statistical and systematic uncertainties in the SM background are included in the shaded band. The red arrows indicate the $m_{\mathrm{eff}}$ selections in the signal regions.

The $m_{\mathrm{eff}}$ distribution in SR0$_{\mathrm{breq}}$ for events passing the signal region requirements except the $m_{\mathrm{eff}}$ requirement. Distributions for data, the estimated SM backgrounds after the background-only fit, and an example SUSY scenario are shown. "Other" is the sum of the $tWZ$, $t\bar{t}WW$, $t\bar{t} ZZ$, $t\bar{t} WH$, $t\bar{t} HH$, $t\bar{t} tW$, and $t\bar{t}t\bar{t}$ backgrounds. The last bin captures the overflow events. The lower panel shows the ratio of the observed data to the expected SM background yield in each bin. Both the statistical and systematic uncertainties in the SM background are included in the shaded band. The red arrows indicate the $m_{\mathrm{eff}}$ selections in the signal regions.

The $m_{\mathrm{eff}}$ distribution in SR1$_{\mathrm{breq}}$ for events passing the signal region requirements except the $m_{\mathrm{eff}}$ requirement. Distributions for data, the estimated SM backgrounds after the background-only fit, and an example SUSY scenario are shown. "Other" is the sum of the $tWZ$, $t\bar{t}WW$, $t\bar{t} ZZ$, $t\bar{t} WH$, $t\bar{t} HH$, $t\bar{t} tW$, and $t\bar{t}t\bar{t}$ backgrounds. The last bin captures the overflow events. The lower panel shows the ratio of the observed data to the expected SM background yield in each bin. Both the statistical and systematic uncertainties in the SM background are included in the shaded band. The red arrows indicate the $m_{\mathrm{eff}}$ selections in the signal regions.

The $m_{\mathrm{eff}}$ distribution in SR2$_{\mathrm{breq}}$ for events passing the signal region requirements except the $m_{\mathrm{eff}}$ requirement. Distributions for data, the estimated SM backgrounds after the background-only fit, and an example SUSY scenario are shown. "Other" is the sum of the $tWZ$, $t\bar{t}WW$, $t\bar{t} ZZ$, $t\bar{t} WH$, $t\bar{t} HH$, $t\bar{t} tW$, and $t\bar{t}t\bar{t}$ backgrounds. The last bin captures the overflow events. The lower panel shows the ratio of the observed data to the expected SM background yield in each bin. Both the statistical and systematic uncertainties in the SM background are included in the shaded band. The red arrows indicate the $m_{\mathrm{eff}}$ selections in the signal regions.

Expected 95% CL exclusion limits on the higgsino GGM models. The limits are set using the statistical combination of disjoint signal regions. Where two (or more) signal regions overlap, the signal region contributing its observed $\mathrm{CL}_{\mathrm{s}}$ value to the combination is the one with the better (best) expected $\mathrm{CL}_{\mathrm{s}}$ value.

$+1\sigma$ expected 95% CL exclusion limits on the higgsino GGM models. The limits are set using the statistical combination of disjoint signal regions. Where two (or more) signal regions overlap, the signal region contributing its observed $\mathrm{CL}_{\mathrm{s}}$ value to the combination is the one with the better (best) expected $\mathrm{CL}_{\mathrm{s}}$ value.

$-1\sigma$ expected 95% CL exclusion limits on the higgsino GGM models. The limits are set using the statistical combination of disjoint signal regions. Where two (or more) signal regions overlap, the signal region contributing its observed $\mathrm{CL}_{\mathrm{s}}$ value to the combination is the one with the better (best) expected $\mathrm{CL}_{\mathrm{s}}$ value.

Observed 95% CL exclusion limits on the higgsino GGM models. The limits are set using the statistical combination of disjoint signal regions. Where two (or more) signal regions overlap, the signal region contributing its observed $\mathrm{CL}_{\mathrm{s}}$ value to the combination is the one with the better (best) expected $\mathrm{CL}_{\mathrm{s}}$ value.

$+1\sigma$ observed 95% CL exclusion limits on the higgsino GGM models. The limits are set using the statistical combination of disjoint signal regions. Where two (or more) signal regions overlap, the signal region contributing its observed $\mathrm{CL}_{\mathrm{s}}$ value to the combination is the one with the better (best) expected $\mathrm{CL}_{\mathrm{s}}$ value.

$-1\sigma$ observed 95% CL exclusion limits on the higgsino GGM models. The limits are set using the statistical combination of disjoint signal regions. Where two (or more) signal regions overlap, the signal region contributing its observed $\mathrm{CL}_{\mathrm{s}}$ value to the combination is the one with the better (best) expected $\mathrm{CL}_{\mathrm{s}}$ value.

Expected 95% CL exclusion limits on wino NLSP pair production with RPV LSP decays via $\lambda_{12k}$, where $k \in{1,2}$. The limits are set using the statistical combination of disjoint signal regions. Where two (or more) signal regions overlap, the signal region contributing its observed $\mathrm{CL}_{\mathrm{s}}$ value to the combination is the one with the better (best) expected $\mathrm{CL}_{\mathrm{s}}$ value.

$+1\sigma$ expected 95% CL exclusion limits on wino NLSP pair production with RPV LSP decays via $\lambda_{12k}$, where $k \in{1,2}$. The limits are set using the statistical combination of disjoint signal regions. Where two (or more) signal regions overlap, the signal region contributing its observed $\mathrm{CL}_{\mathrm{s}}$ value to the combination is the one with the better (best) expected $\mathrm{CL}_{\mathrm{s}}$ value.

$-1\sigma$ expected 95% CL exclusion limits on wino NLSP pair production with RPV LSP decays via $\lambda_{12k}$, where $k \in{1,2}$. The limits are set using the statistical combination of disjoint signal regions. Where two (or more) signal regions overlap, the signal region contributing its observed $\mathrm{CL}_{\mathrm{s}}$ value to the combination is the one with the better (best) expected $\mathrm{CL}_{\mathrm{s}}$ value.

Observed 95% CL exclusion limits on wino NLSP pair production with RPV LSP decays via $\lambda_{12k}$, where $k \in{1,2}$. The limits are set using the statistical combination of disjoint signal regions. Where two (or more) signal regions overlap, the signal region contributing its observed $\mathrm{CL}_{\mathrm{s}}$ value to the combination is the one with the better (best) expected $\mathrm{CL}_{\mathrm{s}}$ value.

$+1\sigma$ bserved 95% CL exclusion limits on wino NLSP pair production with RPV LSP decays via $\lambda_{12k}$, where $k \in{1,2}$. The limits are set using the statistical combination of disjoint signal regions. Where two (or more) signal regions overlap, the signal region contributing its observed $\mathrm{CL}_{\mathrm{s}}$ value to the combination is the one with the better (best) expected $\mathrm{CL}_{\mathrm{s}}$ value.

$-1\sigma$ observed 95% CL exclusion limits on wino NLSP pair production with RPV LSP decays via $\lambda_{12k}$, where $k \in{1,2}$. The limits are set using the statistical combination of disjoint signal regions. Where two (or more) signal regions overlap, the signal region contributing its observed $\mathrm{CL}_{\mathrm{s}}$ value to the combination is the one with the better (best) expected $\mathrm{CL}_{\mathrm{s}}$ value.

Expected 95% CL exclusion limits on wino NLSP pair production with RPV LSP decays via $\lambda_{i33}$, where $i \in{1,2}$. The limits are set using the statistical combination of disjoint signal regions. Where two (or more) signal regions overlap, the signal region contributing its observed $\mathrm{CL}_{\mathrm{s}}$ value to the combination is the one with the better (best) expected $\mathrm{CL}_{\mathrm{s}}$ value.

$+1\sigma$ expected 95% CL exclusion limits on wino NLSP pair production with RPV LSP decays via $\lambda_{i33}$, where $i \in{1,2}$. The limits are set using the statistical combination of disjoint signal regions. Where two (or more) signal regions overlap, the signal region contributing its observed $\mathrm{CL}_{\mathrm{s}}$ value to the combination is the one with the better (best) expected $\mathrm{CL}_{\mathrm{s}}$ value.

$-1\sigma$ expected 95% CL exclusion limits on wino NLSP pair production with RPV LSP decays via $\lambda_{i33}$, where $i \in{1,2}$. The limits are set using the statistical combination of disjoint signal regions. Where two (or more) signal regions overlap, the signal region contributing its observed $\mathrm{CL}_{\mathrm{s}}$ value to the combination is the one with the better (best) expected $\mathrm{CL}_{\mathrm{s}}$ value.

Observed 95% CL exclusion limits on wino NLSP pair production with RPV LSP decays via $\lambda_{i33}$, where $i \in{1,2}$. The limits are set using the statistical combination of disjoint signal regions. Where two (or more) signal regions overlap, the signal region contributing its observed $\mathrm{CL}_{\mathrm{s}}$ value to the combination is the one with the better (best) expected $\mathrm{CL}_{\mathrm{s}}$ value.

$+1\sigma$ observed 95% CL exclusion limits on wino NLSP pair production with RPV LSP decays via $\lambda_{i33}$, where $i \in{1,2}$. The limits are set using the statistical combination of disjoint signal regions. Where two (or more) signal regions overlap, the signal region contributing its observed $\mathrm{CL}_{\mathrm{s}}$ value to the combination is the one with the better (best) expected $\mathrm{CL}_{\mathrm{s}}$ value.

$-1\sigma$ observed 95% CL exclusion limits on wino NLSP pair production with RPV LSP decays via $\lambda_{i33}$, where $i \in{1,2}$. The limits are set using the statistical combination of disjoint signal regions. Where two (or more) signal regions overlap, the signal region contributing its observed $\mathrm{CL}_{\mathrm{s}}$ value to the combination is the one with the better (best) expected $\mathrm{CL}_{\mathrm{s}}$ value.

Expected 95% CL exclusion limits on slepton/sneutrino NLSP pair production with RPV LSP decays via $\lambda_{12k}$, where $k \in{1,2}$. The limits are set using the statistical combination of disjoint signal regions. Where two (or more) signal regions overlap, the signal region contributing its observed $\mathrm{CL}_{\mathrm{s}}$ value to the combination is the one with the better (best) expected $\mathrm{CL}_{\mathrm{s}}$ value.

$+1\sigma$ expected 95% CL exclusion limits on slepton/sneutrino NLSP pair production with RPV LSP decays via $\lambda_{12k}$, where $k \in{1,2}$. The limits are set using the statistical combination of disjoint signal regions. Where two (or more) signal regions overlap, the signal region contributing its observed $\mathrm{CL}_{\mathrm{s}}$ value to the combination is the one with the better (best) expected $\mathrm{CL}_{\mathrm{s}}$ value.

$-1\sigma$ expected 95% CL exclusion limits on slepton/sneutrino NLSP pair production with RPV LSP decays via $\lambda_{12k}$, where $k \in{1,2}$. The limits are set using the statistical combination of disjoint signal regions. Where two (or more) signal regions overlap, the signal region contributing its observed $\mathrm{CL}_{\mathrm{s}}$ value to the combination is the one with the better (best) expected $\mathrm{CL}_{\mathrm{s}}$ value.

Observed 95% CL exclusion limits on slepton/sneutrino NLSP pair production with RPV LSP decays via $\lambda_{12k}$, where $k \in{1,2}$. The limits are set using the statistical combination of disjoint signal regions. Where two (or more) signal regions overlap, the signal region contributing its observed $\mathrm{CL}_{\mathrm{s}}$ value to the combination is the one with the better (best) expected $\mathrm{CL}_{\mathrm{s}}$ value.

$+1\sigma$ observed 95% CL exclusion limits on slepton/sneutrino NLSP pair production with RPV LSP decays via $\lambda_{12k}$, where $k \in{1,2}$. The limits are set using the statistical combination of disjoint signal regions. Where two (or more) signal regions overlap, the signal region contributing its observed $\mathrm{CL}_{\mathrm{s}}$ value to the combination is the one with the better (best) expected $\mathrm{CL}_{\mathrm{s}}$ value.

$-1\sigma$ observed 95% CL exclusion limits on slepton/sneutrino NLSP pair production with RPV LSP decays via $\lambda_{12k}$, where $k \in{1,2}$. The limits are set using the statistical combination of disjoint signal regions. Where two (or more) signal regions overlap, the signal region contributing its observed $\mathrm{CL}_{\mathrm{s}}$ value to the combination is the one with the better (best) expected $\mathrm{CL}_{\mathrm{s}}$ value.

Expected 95% CL exclusion limits on slepton/sneutrino NLSP pair production with RPV LSP decays via $\lambda_{i33}$, where $i \in{1,2}$. The limits are set using the statistical combination of disjoint signal regions. Where two (or more) signal regions overlap, the signal region contributing its observed $\mathrm{CL}_{\mathrm{s}}$ value to the combination is the one with the better (best) expected $\mathrm{CL}_{\mathrm{s}}$ value.

$+1\sigma$ expected 95% CL exclusion limits on slepton/sneutrino NLSP pair production with RPV LSP decays via $\lambda_{i33}$, where $i \in{1,2}$. The limits are set using the statistical combination of disjoint signal regions. Where two (or more) signal regions overlap, the signal region contributing its observed $\mathrm{CL}_{\mathrm{s}}$ value to the combination is the one with the better (best) expected $\mathrm{CL}_{\mathrm{s}}$ value.

$-1\sigma$ expected 95% CL exclusion limits on slepton/sneutrino NLSP pair production with RPV LSP decays via $\lambda_{i33}$, where $i \in{1,2}$. The limits are set using the statistical combination of disjoint signal regions. Where two (or more) signal regions overlap, the signal region contributing its observed $\mathrm{CL}_{\mathrm{s}}$ value to the combination is the one with the better (best) expected $\mathrm{CL}_{\mathrm{s}}$ value.

Observed 95% CL exclusion limits on slepton/sneutrino NLSP pair production with RPV LSP decays via $\lambda_{i33}$, where $i \in{1,2}$. The limits are set using the statistical combination of disjoint signal regions. Where two (or more) signal regions overlap, the signal region contributing its observed $\mathrm{CL}_{\mathrm{s}}$ value to the combination is the one with the better (best) expected $\mathrm{CL}_{\mathrm{s}}$ value.

$+1\sigma$ observed 95% CL exclusion limits on slepton/sneutrino NLSP pair production with RPV LSP decays via $\lambda_{i33}$, where $i \in{1,2}$. The limits are set using the statistical combination of disjoint signal regions. Where two (or more) signal regions overlap, the signal region contributing its observed $\mathrm{CL}_{\mathrm{s}}$ value to the combination is the one with the better (best) expected $\mathrm{CL}_{\mathrm{s}}$ value.

$-1\sigma$ observed 95% CL exclusion limits on slepton/sneutrino NLSP pair production with RPV LSP decays via $\lambda_{i33}$, where $i \in{1,2}$. The limits are set using the statistical combination of disjoint signal regions. Where two (or more) signal regions overlap, the signal region contributing its observed $\mathrm{CL}_{\mathrm{s}}$ value to the combination is the one with the better (best) expected $\mathrm{CL}_{\mathrm{s}}$ value.

Expected 95% CL exclusion limits on gluino NLSP pair production with RPV LSP decays via $\lambda_{12k}$, where $k \in{1,2}$. The limits are set using the statistical combination of disjoint signal regions. Where two (or more) signal regions overlap, the signal region contributing its observed $\mathrm{CL}_{\mathrm{s}}$ value to the combination is the one with the better (best) expected $\mathrm{CL}_{\mathrm{s}}$ value.

$+1\sigma$ expected 95% CL exclusion limits on gluino NLSP pair production with RPV LSP decays via $\lambda_{12k}$, where $k \in{1,2}$. The limits are set using the statistical combination of disjoint signal regions. Where two (or more) signal regions overlap, the signal region contributing its observed $\mathrm{CL}_{\mathrm{s}}$ value to the combination is the one with the better (best) expected $\mathrm{CL}_{\mathrm{s}}$ value.

$-1\sigma$ expected 95% CL exclusion limits on gluino NLSP pair production with RPV LSP decays via $\lambda_{12k}$, where $k \in{1,2}$. The limits are set using the statistical combination of disjoint signal regions. Where two (or more) signal regions overlap, the signal region contributing its observed $\mathrm{CL}_{\mathrm{s}}$ value to the combination is the one with the better (best) expected $\mathrm{CL}_{\mathrm{s}}$ value.

Observed 95% CL exclusion limits on gluino NLSP pair production with RPV LSP decays via $\lambda_{12k}$, where $k \in{1,2}$. The limits are set using the statistical combination of disjoint signal regions. Where two (or more) signal regions overlap, the signal region contributing its observed $\mathrm{CL}_{\mathrm{s}}$ value to the combination is the one with the better (best) expected $\mathrm{CL}_{\mathrm{s}}$ value.

$+1\sigma$ observed 95% CL exclusion limits on gluino NLSP pair production with RPV LSP decays via $\lambda_{12k}$, where $k \in{1,2}$. The limits are set using the statistical combination of disjoint signal regions. Where two (or more) signal regions overlap, the signal region contributing its observed $\mathrm{CL}_{\mathrm{s}}$ value to the combination is the one with the better (best) expected $\mathrm{CL}_{\mathrm{s}}$ value.

$-1\sigma$ observed 95% CL exclusion limits on gluino NLSP pair production with RPV LSP decays via $\lambda_{12k}$, where $k \in{1,2}$. The limits are set using the statistical combination of disjoint signal regions. Where two (or more) signal regions overlap, the signal region contributing its observed $\mathrm{CL}_{\mathrm{s}}$ value to the combination is the one with the better (best) expected $\mathrm{CL}_{\mathrm{s}}$ value.

Expected 95% CL exclusion limits on gluino NLSP pair production with RPV LSP decays via $\lambda_{i33}$, where $i \in{1,2}$. The limits are set using the statistical combination of disjoint signal regions. Where two (or more) signal regions overlap, the signal region contributing its observed $\mathrm{CL}_{\mathrm{s}}$ value to the combination is the one with the better (best) expected $\mathrm{CL}_{\mathrm{s}}$ value.

$+1\sigma$ expected 95% CL exclusion limits on gluino NLSP pair production with RPV LSP decays via $\lambda_{i33}$, where $i \in{1,2}$. The limits are set using the statistical combination of disjoint signal regions. Where two (or more) signal regions overlap, the signal region contributing its observed $\mathrm{CL}_{\mathrm{s}}$ value to the combination is the one with the better (best) expected $\mathrm{CL}_{\mathrm{s}}$ value.

$-1\sigma$ expected 95% CL exclusion limits on gluino NLSP pair production with RPV LSP decays via $\lambda_{i33}$, where $i \in{1,2}$. The limits are set using the statistical combination of disjoint signal regions. Where two (or more) signal regions overlap, the signal region contributing its observed $\mathrm{CL}_{\mathrm{s}}$ value to the combination is the one with the better (best) expected $\mathrm{CL}_{\mathrm{s}}$ value.

Observed 95% CL exclusion limits on gluino NLSP pair production with RPV LSP decays via $\lambda_{i33}$, where $i \in{1,2}$. The limits are set using the statistical combination of disjoint signal regions. Where two (or more) signal regions overlap, the signal region contributing its observed $\mathrm{CL}_{\mathrm{s}}$ value to the combination is the one with the better (best) expected $\mathrm{CL}_{\mathrm{s}}$ value.

$+1\sigma$ observed 95% CL exclusion limits on gluino NLSP pair production with RPV LSP decays via $\lambda_{i33}$, where $i \in{1,2}$. The limits are set using the statistical combination of disjoint signal regions. Where two (or more) signal regions overlap, the signal region contributing its observed $\mathrm{CL}_{\mathrm{s}}$ value to the combination is the one with the better (best) expected $\mathrm{CL}_{\mathrm{s}}$ value.

$-1\sigma$ observed 95% CL exclusion limits on gluino NLSP pair production with RPV LSP decays via $\lambda_{i33}$, where $i \in{1,2}$. The limits are set using the statistical combination of disjoint signal regions. Where two (or more) signal regions overlap, the signal region contributing its observed $\mathrm{CL}_{\mathrm{s}}$ value to the combination is the one with the better (best) expected $\mathrm{CL}_{\mathrm{s}}$ value.

Observed upper limit on the signal cross section in fb for the wino NLSP models with RPV LSP decays via $\lambda_{12k}$ where $k \in{1,2}$. The interpolation between signal scenarios studied is included for illustration purposes only and may be subject to interpolation effects in sparsely populated areas.

Observed upper limit on the signal cross section in fb for the wino NLSP models with RPV LSP decays via $\lambda_{i33}$ where $i \in{1,2}$. The interpolation between signal scenarios studied is included for illustration purposes only and may be subject to interpolation effects in sparsely populated areas.

Observed upper limit on the signal cross section in fb for the slepton/sneutrino NLSP models with RPV LSP decays via $\lambda_{12k}$ where $k \in{1,2}$. The interpolation between signal scenarios studied is included for illustration purposes only and may be subject to interpolation effects in sparsely populated areas.

Observed upper limit on the signal cross section in fb for the slepton/sneutrino NLSP models with RPV LSP decays via $\lambda_{i33}$ where $i \in{1,2}$. The interpolation between signal scenarios studied is included for illustration purposes only and may be subject to interpolation effects in sparsely populated areas.

Observed upper limit on the signal cross section in fb for the gluino NLSP models with RPV LSP decays via $\lambda_{12k}$ where $k \in{1,2}$. The interpolation between signal scenarios studied is included for illustration purposes only and may be subject to interpolation effects in sparsely populated areas.

Observed upper limit on the signal cross section in fb for the gluino NLSP models with RPV LSP decays via $\lambda_{i33}$ where $i \in{1,2}$. The interpolation between signal scenarios studied is included for illustration purposes only and may be subject to interpolation effects in sparsely populated areas.

Observed upper limit on the signal cross section in fb for the higgsino GGM models. The interpolation between signal scenarios studied is included for illustration purposes only and may be subject to interpolation effects in sparsely populated areas.

Best expected SR for the wino NLSP models with RPV LSP decays via $\lambda_{12k}$ where $k \in{1,2}$. A value of 1 corresponds to SR0$_{\mathrm{bveto}}^{\mathrm{loose}}$SR1$_{\mathrm{bveto}}^{\mathrm{loose}}$SR2$_{\mathrm{bveto}}^{\mathrm{loose}}$, 2 corresponds to SR0$_{\mathrm{bveto}}^{\mathrm{loose}}$SR1$_{\mathrm{bveto}}^{\mathrm{loose}}$SR2$_{\mathrm{bveto}}^{\mathrm{tight}}$, 3 corresponds to SR0$_{\mathrm{bveto}}^{\mathrm{loose}}$SR1$_{\mathrm{bveto}}^{\mathrm{tight}}$SR2$_{\mathrm{bveto}}^{\mathrm{loose}}$, 4 corresponds to SR0$_{\mathrm{bveto}}^{\mathrm{loose}}$SR1$_{\mathrm{bveto}}^{\mathrm{tight}}$SR2$_{\mathrm{bveto}}^{\mathrm{tight}}$, and 5 corresponds to SR0$_{\mathrm{bveto}}^{\mathrm{tight}}$SR1$_{\mathrm{bveto}}^{\mathrm{tight}}$SR2$_{\mathrm{bveto}}^{\mathrm{tight}}$.

Best expected SR for the wino NLSP models with RPV LSP decays via $\lambda_{i33}$ where $i \in{1,2}$. A value of 1 corresponds to SR0$_{\mathrm{bveto}}^{\mathrm{loose}}$SR1$_{\mathrm{bveto}}^{\mathrm{loose}}$SR2$_{\mathrm{bveto}}^{\mathrm{loose}}$, 2 corresponds to SR0$_{\mathrm{bveto}}^{\mathrm{loose}}$SR1$_{\mathrm{bveto}}^{\mathrm{loose}}$SR2$_{\mathrm{bveto}}^{\mathrm{tight}}$, 3 corresponds to SR0$_{\mathrm{bveto}}^{\mathrm{loose}}$SR1$_{\mathrm{bveto}}^{\mathrm{tight}}$SR2$_{\mathrm{bveto}}^{\mathrm{loose}}$, 4 corresponds to SR0$_{\mathrm{bveto}}^{\mathrm{loose}}$SR1$_{\mathrm{bveto}}^{\mathrm{tight}}$SR2$_{\mathrm{bveto}}^{\mathrm{tight}}$, and 5 corresponds to SR0$_{\mathrm{bveto}}^{\mathrm{tight}}$SR1$_{\mathrm{bveto}}^{\mathrm{tight}}$SR2$_{\mathrm{bveto}}^{\mathrm{tight}}$.

Best expected SR for the slepton/sneutrino NLSP models with RPV LSP decays via $\lambda_{12k}$ where $k \in{1,2}$. A value of 1 corresponds to SR0$_{\mathrm{bveto}}^{\mathrm{loose}}$SR1$_{\mathrm{bveto}}^{\mathrm{loose}}$SR2$_{\mathrm{bveto}}^{\mathrm{loose}}$, 2 corresponds to SR0$_{\mathrm{bveto}}^{\mathrm{loose}}$SR1$_{\mathrm{bveto}}^{\mathrm{loose}}$SR2$_{\mathrm{bveto}}^{\mathrm{tight}}$, 3 corresponds to SR0$_{\mathrm{bveto}}^{\mathrm{loose}}$SR1$_{\mathrm{bveto}}^{\mathrm{tight}}$SR2$_{\mathrm{bveto}}^{\mathrm{loose}}$, 4 corresponds to SR0$_{\mathrm{bveto}}^{\mathrm{loose}}$SR1$_{\mathrm{bveto}}^{\mathrm{tight}}$SR2$_{\mathrm{bveto}}^{\mathrm{tight}}$, and 5 corresponds to SR0$_{\mathrm{bveto}}^{\mathrm{tight}}$SR1$_{\mathrm{bveto}}^{\mathrm{tight}}$SR2$_{\mathrm{bveto}}^{\mathrm{tight}}$.

Best expected SR for the slepton/sneutrino NLSP models with RPV LSP decays via $\lambda_{i33}$ where $i \in{1,2}$. A value of 1 corresponds to SR0$_{\mathrm{bveto}}^{\mathrm{loose}}$SR1$_{\mathrm{bveto}}^{\mathrm{loose}}$SR2$_{\mathrm{bveto}}^{\mathrm{loose}}$, 2 corresponds to SR0$_{\mathrm{bveto}}^{\mathrm{loose}}$SR1$_{\mathrm{bveto}}^{\mathrm{loose}}$SR2$_{\mathrm{bveto}}^{\mathrm{tight}}$, 3 corresponds to SR0$_{\mathrm{bveto}}^{\mathrm{loose}}$SR1$_{\mathrm{bveto}}^{\mathrm{tight}}$SR2$_{\mathrm{bveto}}^{\mathrm{loose}}$, 4 corresponds to SR0$_{\mathrm{bveto}}^{\mathrm{loose}}$SR1$_{\mathrm{bveto}}^{\mathrm{tight}}$SR2$_{\mathrm{bveto}}^{\mathrm{tight}}$, and 5 corresponds to SR0$_{\mathrm{bveto}}^{\mathrm{tight}}$SR1$_{\mathrm{bveto}}^{\mathrm{tight}}$SR2$_{\mathrm{bveto}}^{\mathrm{tight}}$.

Best expected SR for the gluino NLSP models with RPV LSP decays via $\lambda_{12k}$ where $k \in{1,2}$. A value of 1 corresponds to SR0$_{\mathrm{bveto}}^{\mathrm{loose}}$SR1$_{\mathrm{bveto}}^{\mathrm{loose}}$SR2$_{\mathrm{bveto}}^{\mathrm{loose}}$, 2 corresponds to SR0$_{\mathrm{bveto}}^{\mathrm{loose}}$SR1$_{\mathrm{bveto}}^{\mathrm{loose}}$SR2$_{\mathrm{bveto}}^{\mathrm{tight}}$, 3 corresponds to SR0$_{\mathrm{bveto}}^{\mathrm{loose}}$SR1$_{\mathrm{bveto}}^{\mathrm{tight}}$SR2$_{\mathrm{bveto}}^{\mathrm{loose}}$, 4 corresponds to SR0$_{\mathrm{bveto}}^{\mathrm{loose}}$SR1$_{\mathrm{bveto}}^{\mathrm{tight}}$SR2$_{\mathrm{bveto}}^{\mathrm{tight}}$, and 5 corresponds to SR0$_{\mathrm{bveto}}^{\mathrm{tight}}$SR1$_{\mathrm{bveto}}^{\mathrm{tight}}$SR2$_{\mathrm{bveto}}^{\mathrm{tight}}$.

Best expected SR for the gluino NLSP models with RPV LSP decays via $\lambda_{i33}$ where $i \in{1,2}$. A value of 1 corresponds to SR0$_{\mathrm{bveto}}^{\mathrm{loose}}$SR1$_{\mathrm{bveto}}^{\mathrm{loose}}$SR2$_{\mathrm{bveto}}^{\mathrm{loose}}$, 2 corresponds to SR0$_{\mathrm{bveto}}^{\mathrm{loose}}$SR1$_{\mathrm{bveto}}^{\mathrm{loose}}$SR2$_{\mathrm{bveto}}^{\mathrm{tight}}$, 3 corresponds to SR0$_{\mathrm{bveto}}^{\mathrm{loose}}$SR1$_{\mathrm{bveto}}^{\mathrm{tight}}$SR2$_{\mathrm{bveto}}^{\mathrm{loose}}$, 4 corresponds to SR0$_{\mathrm{bveto}}^{\mathrm{loose}}$SR1$_{\mathrm{bveto}}^{\mathrm{tight}}$SR2$_{\mathrm{bveto}}^{\mathrm{tight}}$, and 5 corresponds to SR0$_{\mathrm{bveto}}^{\mathrm{tight}}$SR1$_{\mathrm{bveto}}^{\mathrm{tight}}$SR2$_{\mathrm{bveto}}^{\mathrm{tight}}$.

Best expected SR for the higgsino GGM models. A value of 6 corresponds to SR0-ZZ$^{\mathrm{loose}}$, 7 corresponds to SR0-ZZ$^{\mathrm{tight}}$, 8 corresponds to SR0-ZZ$^{\mathrm{loose}}_{\mathrm{bveto}}$, and 9 corresponds to SR0-ZZ$^{\mathrm{tight}}_{\mathrm{bveto}}$.

Acceptance across the wino NLSP $\lambda_{12k}\neq 0$ models for SR0$_{\mathrm{bveto}}^{\mathrm{loose}}$. The interpolation between signal scenarios studied is included for illustration purposes only and may be subject to interpolation effects in sparsely populated areas.

Efficiency across the wino NLSP $\lambda_{12k}\neq 0$ models for SR0$_{\mathrm{bveto}}^{\mathrm{loose}}$. The interpolation between signal scenarios studied is included for illustration purposes only and may be subject to interpolation effects in sparsely populated areas.

Acceptance across the wino NLSP $\lambda_{12k}\neq 0$ models for SR0$_{\mathrm{bveto}}^{\mathrm{tight}}$. The interpolation between signal scenarios studied is included for illustration purposes only and may be subject to interpolation effects in sparsely populated areas.

Efficiency across the wino NLSP $\lambda_{12k}\neq 0$ models for SR0$_{\mathrm{bveto}}^{\mathrm{tight}}$. The interpolation between signal scenarios studied is included for illustration purposes only and may be subject to interpolation effects in sparsely populated areas.

Acceptance across the wino NLSP $\lambda_{12k}\neq 0$ models for SR0$_{\mathrm{breq}}$. The interpolation between signal scenarios studied is included for illustration purposes only and may be subject to interpolation effects in sparsely populated areas.

Efficiency across the wino NLSP $\lambda_{12k}\neq 0$ models for SR0$_{\mathrm{breq}}$. The interpolation between signal scenarios studied is included for illustration purposes only and may be subject to interpolation effects in sparsely populated areas.

Acceptance across the wino NLSP $\lambda_{i33}\neq 0$ models for SR1$_{\mathrm{bveto}}^{\mathrm{loose}}$. The interpolation between signal scenarios studied is included for illustration purposes only and may be subject to interpolation effects in sparsely populated areas.

Efficiency across the wino NLSP $\lambda_{i33}\neq 0$ models for SR1$_{\mathrm{bveto}}^{\mathrm{loose}}$. The interpolation between signal scenarios studied is included for illustration purposes only and may be subject to interpolation effects in sparsely populated areas.

Acceptance across the wino NLSP $\lambda_{i33}\neq 0$ models for SR1$_{\mathrm{bveto}}^{\mathrm{tight}}$. The interpolation between signal scenarios studied is included for illustration purposes only and may be subject to interpolation effects in sparsely populated areas.

Efficiency across the wino NLSP $\lambda_{i33}\neq 0$ models for SR1$_{\mathrm{bveto}}^{\mathrm{tight}}$. The interpolation between signal scenarios studied is included for illustration purposes only and may be subject to interpolation effects in sparsely populated areas.

Acceptance across the wino NLSP $\lambda_{i33}\neq 0$ models for SR1$_{\mathrm{breq}}$. The interpolation between signal scenarios studied is included for illustration purposes only and may be subject to interpolation effects in sparsely populated areas.

Efficiency across the wino NLSP $\lambda_{i33}\neq 0$ models for SR1$_{\mathrm{breq}}$. The interpolation between signal scenarios studied is included for illustration purposes only and may be subject to interpolation effects in sparsely populated areas.

Acceptance across the wino NLSP $\lambda_{i33}\neq 0$ models for SR2$_{\mathrm{bveto}}^{\mathrm{loose}}$. The interpolation between signal scenarios studied is included for illustration purposes only and may be subject to interpolation effects in sparsely populated areas.

Efficiency across the wino NLSP $\lambda_{i33}\neq 0$ models for SR2$_{\mathrm{bveto}}^{\mathrm{loose}}$. The interpolation between signal scenarios studied is included for illustration purposes only and may be subject to interpolation effects in sparsely populated areas.

Acceptance across the wino NLSP $\lambda_{i33}\neq 0$ models for SR2$_{\mathrm{bveto}}^{\mathrm{tight}}$. The interpolation between signal scenarios studied is included for illustration purposes only and may be subject to interpolation effects in sparsely populated areas.

Efficiency across the wino NLSP $\lambda_{i33}\neq 0$ models for SR2$_{\mathrm{bveto}}^{\mathrm{tight}}$. The interpolation between signal scenarios studied is included for illustration purposes only and may be subject to interpolation effects in sparsely populated areas.

Acceptance across the wino NLSP $\lambda_{i33}\neq 0$ models for SR2$_{\mathrm{breq}}$. The interpolation between signal scenarios studied is included for illustration purposes only and may be subject to interpolation effects in sparsely populated areas.

Efficiency across the wino NLSP $\lambda_{i33}\neq 0$ models for SR2$_{\mathrm{breq}}$. The interpolation between signal scenarios studied is included for illustration purposes only and may be subject to interpolation effects in sparsely populated areas.

Acceptance across the GGM Higgsino grid for SR0-ZZ$^{\mathrm{loose}}$. The interpolation between signal scenarios studied is included for illustration purposes only and may be subject to interpolation effects in sparsely populated areas.

Efficiency across the GGM Higgsino grid for SR0-ZZ$^{\mathrm{loose}}$. The interpolation between signal scenarios studied is included for illustration purposes only and may be subject to interpolation effects in sparsely populated areas.

Acceptance across the GGM Higgsino grid for SR0-ZZ$^{\mathrm{tight}}$. The interpolation between signal scenarios studied is included for illustration purposes only and may be subject to interpolation effects in sparsely populated areas.

Efficiency across the GGM Higgsino grid for SR0-ZZ$^{\mathrm{tight}}$. The interpolation between signal scenarios studied is included for illustration purposes only and may be subject to interpolation effects in sparsely populated areas.

Acceptance across the GGM Higgsino grid for SR0-ZZ$_{\mathrm{bveto}}^{\mathrm{loose}}$. The interpolation between signal scenarios studied is included for illustration purposes only and may be subject to interpolation effects in sparsely populated areas.

Efficiency across the GGM Higgsino grid for SR0-ZZ$_{\mathrm{bveto}}^{\mathrm{loose}}$. The interpolation between signal scenarios studied is included for illustration purposes only and may be subject to interpolation effects in sparsely populated areas.

Acceptance across the GGM Higgsino grid for SR0-ZZ$_{\mathrm{bveto}}^{\mathrm{tight}}$. The interpolation between signal scenarios studied is included for illustration purposes only and may be subject to interpolation effects in sparsely populated areas.

Efficiency across the GGM Higgsino grid for SR0-ZZ$_{\mathrm{bveto}}^{\mathrm{tight}}$. The interpolation between signal scenarios studied is included for illustration purposes only and may be subject to interpolation effects in sparsely populated areas.

The $p_{\mathrm{T}}$ of the light leptons in distribution in SR0$_{\mathrm{bveto}}^{\mathrm{loose}}$. Distributions for data, the estimated SM backgrounds after the background-only fit, and an example SUSY scenario are shown. "Other" is the sum of the $tWZ$, $t\bar{t}WW$, $t\bar{t} ZZ$, $t\bar{t} WH$, $t\bar{t} HH$, $t\bar{t} tW$, and $t\bar{t}t\bar{t}$ backgrounds. The last bin captures the overflow events. The lower panel shows the ratio of the observed data to the expected SM background yield in each bin. Both the statistical and systematic uncertainties in the SM background are included in the shaded band.

The $p_{\mathrm{T}}$ of the light leptons in distribution in SR0-ZZ$^{\mathrm{loose}}$. Distributions for data, the estimated SM backgrounds after the background-only fit, and an example SUSY scenario are shown. "Other" is the sum of the $tWZ$, $t\bar{t}WW$, $t\bar{t} ZZ$, $t\bar{t} WH$, $t\bar{t} HH$, $t\bar{t} tW$, and $t\bar{t}t\bar{t}$ backgrounds. The last bin captures the overflow events. The lower panel shows the ratio of the observed data to the expected SM background yield in each bin. Both the statistical and systematic uncertainties in the SM background are included in the shaded band.

The $p_{\mathrm{T}}$ of the light leptons in distribution in SR0-ZZ$^{\mathrm{tight}}$. Distributions for data, the estimated SM backgrounds after the background-only fit, and an example SUSY scenario are shown. "Other" is the sum of the $tWZ$, $t\bar{t}WW$, $t\bar{t} ZZ$, $t\bar{t} WH$, $t\bar{t} HH$, $t\bar{t} tW$, and $t\bar{t}t\bar{t}$ backgrounds. The last bin captures the overflow events. The lower panel shows the ratio of the observed data to the expected SM background yield in each bin. Both the statistical and systematic uncertainties in the SM background are included in the shaded band.

The $p_{\mathrm{T}}$ of the light leptons in distribution in SR0-ZZ$_{\mathrm{bveto}}^{\mathrm{loose}}$. Distributions for data, the estimated SM backgrounds after the background-only fit, and an example SUSY scenario are shown. "Other" is the sum of the $tWZ$, $t\bar{t}WW$, $t\bar{t} ZZ$, $t\bar{t} WH$, $t\bar{t} HH$, $t\bar{t} tW$, and $t\bar{t}t\bar{t}$ backgrounds. The last bin captures the overflow events. The lower panel shows the ratio of the observed data to the expected SM background yield in each bin. Both the statistical and systematic uncertainties in the SM background are included in the shaded band.

The $p_{\mathrm{T}}$ of the light leptons in distribution in SR5L. Distributions for data, the estimated SM backgrounds after the background-only fit, and an example SUSY scenario are shown. "Other" is the sum of the $tWZ$, $t\bar{t}WW$, $t\bar{t} ZZ$, $t\bar{t} WH$, $t\bar{t} HH$, $t\bar{t} tW$, and $t\bar{t}t\bar{t}$ backgrounds. The last bin captures the overflow events. The lower panel shows the ratio of the observed data to the expected SM background yield in each bin. Both the statistical and systematic uncertainties in the SM background are included in the shaded band.

The $p_{\mathrm{T}}$ of the light leptons in distribution in SR1$_{\mathrm{bveto}}^{\mathrm{loose}}$. Distributions for data, the estimated SM backgrounds after the background-only fit, and an example SUSY scenario are shown. "Other" is the sum of the $tWZ$, $t\bar{t}WW$, $t\bar{t} ZZ$, $t\bar{t} WH$, $t\bar{t} HH$, $t\bar{t} tW$, and $t\bar{t}t\bar{t}$ backgrounds. The last bin captures the overflow events. The lower panel shows the ratio of the observed data to the expected SM background yield in each bin. Both the statistical and systematic uncertainties in the SM background are included in the shaded band.

The $p_{\mathrm{T}}$ of the taus leptons in distribution in SR1$_{\mathrm{bveto}}^{\mathrm{loose}}$. Distributions for data, the estimated SM backgrounds after the background-only fit, and an example SUSY scenario are shown. "Other" is the sum of the $tWZ$, $t\bar{t}WW$, $t\bar{t} ZZ$, $t\bar{t} WH$, $t\bar{t} HH$, $t\bar{t} tW$, and $t\bar{t}t\bar{t}$ backgrounds. The last bin captures the overflow events. The lower panel shows the ratio of the observed data to the expected SM background yield in each bin. Both the statistical and systematic uncertainties in the SM background are included in the shaded band.

The $p_{\mathrm{T}}$ of the light taus in distribution in SR2$_{\mathrm{bveto}}^{\mathrm{loose}}$. Distributions for data, the estimated SM backgrounds after the background-only fit, and an example SUSY scenario are shown. "Other" is the sum of the $tWZ$, $t\bar{t}WW$, $t\bar{t} ZZ$, $t\bar{t} WH$, $t\bar{t} HH$, $t\bar{t} tW$, and $t\bar{t}t\bar{t}$ backgrounds. The last bin captures the overflow events. The lower panel shows the ratio of the observed data to the expected SM background yield in each bin. Both the statistical and systematic uncertainties in the SM background are included in the shaded band.

The lepton flavour and multiplicities in events with four light leptons and a Z veto. Distributions for data, the estimated SM backgrounds after the background-only fit, and an example SUSY scenario are shown. "Other" is the sum of the $tWZ$, $t\bar{t}WW$, $t\bar{t} ZZ$, $t\bar{t} WH$, $t\bar{t} HH$, $t\bar{t} tW$, and $t\bar{t}t\bar{t}$ backgrounds. The last bin captures the overflow events. The lower panel shows the ratio of the observed data to the expected SM background yield in each bin. Both the statistical and systematic uncertainties in the SM background are included in the shaded band.

The lepton flavour and multiplicities in events with four light leptons and one Z candidate. Distributions for data, the estimated SM backgrounds after the background-only fit, and an example SUSY scenario are shown. "Other" is the sum of the $tWZ$, $t\bar{t}WW$, $t\bar{t} ZZ$, $t\bar{t} WH$, $t\bar{t} HH$, $t\bar{t} tW$, and $t\bar{t}t\bar{t}$ backgrounds. The last bin captures the overflow events. The lower panel shows the ratio of the observed data to the expected SM background yield in each bin. Both the statistical and systematic uncertainties in the SM background are included in the shaded band.

The lepton flavour and multiplicities in events with four light leptons and two Z candidates. Distributions for data, the estimated SM backgrounds after the background-only fit, and an example SUSY scenario are shown. "Other" is the sum of the $tWZ$, $t\bar{t}WW$, $t\bar{t} ZZ$, $t\bar{t} WH$, $t\bar{t} HH$, $t\bar{t} tW$, and $t\bar{t}t\bar{t}$ backgrounds. The last bin captures the overflow events. The lower panel shows the ratio of the observed data to the expected SM background yield in each bin. Both the statistical and systematic uncertainties in the SM background are included in the shaded band.

The lepton flavour and multiplicities in events with exactly five light leptons. Distributions for data, the estimated SM backgrounds after the background-only fit, and an example SUSY scenario are shown. "Other" is the sum of the $tWZ$, $t\bar{t}WW$, $t\bar{t} ZZ$, $t\bar{t} WH$, $t\bar{t} HH$, $t\bar{t} tW$, and $t\bar{t}t\bar{t}$ backgrounds. The last bin captures the overflow events. The lower panel shows the ratio of the observed data to the expected SM background yield in each bin. Both the statistical and systematic uncertainties in the SM background are included in the shaded band.

The lepton flavour and multiplicities in events with three light leptons and one tau and a Z veto. Distributions for data, the estimated SM backgrounds after the background-only fit, and an example SUSY scenario are shown. "Other" is the sum of the $tWZ$, $t\bar{t}WW$, $t\bar{t} ZZ$, $t\bar{t} WH$, $t\bar{t} HH$, $t\bar{t} tW$, and $t\bar{t}t\bar{t}$ backgrounds. The last bin captures the overflow events. The lower panel shows the ratio of the observed data to the expected SM background yield in each bin. Both the statistical and systematic uncertainties in the SM background are included in the shaded band.

The lepton flavour and multiplicities in events with three light leptons and one tau and one Z candidate. Distributions for data, the estimated SM backgrounds after the background-only fit, and an example SUSY scenario are shown. "Other" is the sum of the $tWZ$, $t\bar{t}WW$, $t\bar{t} ZZ$, $t\bar{t} WH$, $t\bar{t} HH$, $t\bar{t} tW$, and $t\bar{t}t\bar{t}$ backgrounds. The last bin captures the overflow events. The lower panel shows the ratio of the observed data to the expected SM background yield in each bin. Both the statistical and systematic uncertainties in the SM background are included in the shaded band.

The lepton flavour and multiplicities in events with two light leptons and two taus and a Z veto. Distributions for data, the estimated SM backgrounds after the background-only fit, and an example SUSY scenario are shown. "Other" is the sum of the $tWZ$, $t\bar{t}WW$, $t\bar{t} ZZ$, $t\bar{t} WH$, $t\bar{t} HH$, $t\bar{t} tW$, and $t\bar{t}t\bar{t}$ backgrounds. The last bin captures the overflow events. The lower panel shows the ratio of the observed data to the expected SM background yield in each bin. Both the statistical and systematic uncertainties in the SM background are included in the shaded band.

The lepton flavour and multiplicities in events with two light leptons and two taus and one Z candidate. Distributions for data, the estimated SM backgrounds after the background-only fit, and an example SUSY scenario are shown. "Other" is the sum of the $tWZ$, $t\bar{t}WW$, $t\bar{t} ZZ$, $t\bar{t} WH$, $t\bar{t} HH$, $t\bar{t} tW$, and $t\bar{t}t\bar{t}$ backgrounds. The last bin captures the overflow events. The lower panel shows the ratio of the observed data to the expected SM background yield in each bin. Both the statistical and systematic uncertainties in the SM background are included in the shaded band.

Cutflow event yields in regions SR0$_{\mathrm{bveto}}^{\mathrm{loose}}$, SR0$_{\mathrm{bveto}}^{\mathrm{tight}}$, SR0$_{\mathrm{breq}}$, and SR5L for RPV models with the $\lambda_{12k}\neq 0$ coupling. All yields correspond to weighted events, so that effects from lepton reconstruction efficiencies, trigger corrections, pileup reweighting, etc., are included. They are normalized to the integrated luminosity of the data sample, $\int L dt = 139\,\mbox{fb\(^{-1}\)}$. The preliminary event reduction is a centralized stage where at least two electrons/muons with uncalibrated $p_{\mathrm{T}} >$ 9 GeV are required.

Cutflow event yields in regions SR1$_{\mathrm{bveto}}^{\mathrm{loose}}$, SR1$_{\mathrm{bveto}}^{\mathrm{tight}}$, and SR1$_{\mathrm{breq}}$ for RPV models with the $\lambda_{i33}\neq 0$ coupling. All yields correspond to weighted events, so that effects from lepton reconstruction efficiencies, trigger corrections, pileup reweighting, etc., are included. They are normalized to the integrated luminosity of the data sample, $\int L dt = 139\,\mbox{fb\(^{-1}\)}$. The preliminary event reduction is a centralized stage where at least two electrons/muons with uncalibrated $p_{\mathrm{T}} >$ 9 GeV are required.

Cutflow event yields in regions SR2$_{\mathrm{bveto}}^{\mathrm{loose}}$, SR2$_{\mathrm{bveto}}^{\mathrm{tight}}$, and SR2$_{\mathrm{breq}}$ for RPV models with the $\lambda_{i33}\neq 0$ coupling. All yields correspond to weighted events, so that effects from lepton reconstruction efficiencies, trigger corrections, pileup reweighting, etc., are included. They are normalized to the integrated luminosity of the data sample, $\int L dt = 139\,\mbox{fb\(^{-1}\)}$. The preliminary event reduction is a centralized stage where at least two electrons/muons with uncalibrated $p_{\mathrm{T}} >$ 9 GeV are required.

Cutflow event yields in regions SR0-ZZ$^{\mathrm{loose}}$, SR0-ZZ$^{\mathrm{tight}}$, SR0-ZZ$_{\mathrm{bveto}}^{\mathrm{loose}}$, SR0-ZZ$_{\mathrm{bveto}}^{\mathrm{tight}}$, and SR5L the higgsino GGM RPC model with BR($\tilde{\chi}^{0}_1 \rightarrow Z \tilde{G}$) = 50% and higgsino masses of 200 GeV, or BR($\tilde{\chi}^{0}_1 \rightarrow Z \tilde{G}$) = 100% and higgsino masses of 300 GeV. All yields correspond to weighted events, so that effects from lepton reconstruction efficiencies, trigger corrections, pileup reweighting, etc., are included. They are normalized to the integrated luminosity of the data sample, $\int L dt = 139\,\mbox{fb\(^{-1}\)}$. The generator filter is a selection of $\geq$4e/$\mu$/$\tau_{\mathrm{had-vis}}$ leptons with $p_{\mathrm{T}}(e,\mu)>4$GeV, $p_{\mathrm{T}}(\tau_{\mathrm{had-vis}})>15$GeV and $|\eta|<2.8$ and is applied during the MC generation of the simulated events. The preliminary event reduction is a centralized stage where at least two electrons/muons with uncalibrated $p_{\mathrm{T}} > 9$ GeV are required.

The definitions of the four tetralepton signal regions. The regions are defined to target different $b$-jet multiplicities and flavour combinations of the non-Z leptons.

The absolute particle-level differential cross-section measured in the fiducial phase-space as a function of the $p_{\textrm{T}}$ of the $Z$ boson in the 3$\ell$+4$\ell$ channel. The uncertainty is decomposed into four components which are the signal modelling uncertainty, the background modelling uncertainty, the experimental uncertainty, and the data statistical uncertainty.

The absolute parton-level differential cross-section measured in the fiducial phase-space as a function of the $p_{\textrm{T}}$ of the $Z$ boson in the 3$\ell$+4$\ell$ channel. The uncertainty is decomposed into four components which are the signal modelling uncertainty, the background modelling uncertainty, the experimental uncertainty, and the data statistical uncertainty.

The normalised particle-level differential cross-section measured in the fiducial phase-space as a function of the $p_{\textrm{T}}$ of the $Z$ boson in the 3$\ell$+4$\ell$ channel. The uncertainty is decomposed into four components which are the signal modelling uncertainty, the background modelling uncertainty, the experimental uncertainty, and the data statistical uncertainty.

The normalised parton-level differential cross-section measured in the fiducial phase-space as a function of the $p_{\textrm{T}}$ of the $Z$ boson in the 3$\ell$+4$\ell$ channel. The uncertainty is decomposed into four components which are the signal modelling uncertainty, the background modelling uncertainty, the experimental uncertainty, and the data statistical uncertainty.

The absolute parton-level differential cross-section measured in the fiducial phase-space as a function of the absolute value of rapidity of the $Z$ boson in the 3$\ell$+4$\ell$ channel. The uncertainty is decomposed into four components which are the signal modelling uncertainty, the background modelling uncertainty, the experimental uncertainty, and the data statistical uncertainty.

The absolute particle-level differential cross-section measured in the fiducial phase-space as a function of the number of jets in the 3$\ell$ channel. The uncertainty is decomposed into four components which are the signal modelling uncertainty, the background modelling uncertainty, the experimental uncertainty, and the data statistical uncertainty.

The absolute particle-level differential cross-section measured in the fiducial phase-space as a function of the number of jets in the 4$\ell$ channel. The uncertainty is decomposed into four components which are the signal modelling uncertainty, the background modelling uncertainty, the experimental uncertainty, and the data statistical uncertainty.

The absolute parton-level differential cross-section measured in the fiducial phase-space as a function of the $p_{\textrm{T}}^{l \textrm{non-}Z}$ in the 3$\ell$ channel. The uncertainty is decomposed into four components which are the signal modelling uncertainty, the background modelling uncertainty, the experimental uncertainty, and the data statistical uncertainty.

The absolute parton-level differential cross-section measured in the fiducial phase-space as a function of the $|\Delta \phi (Z, t_{\textrm{lep}})|/\pi$ in the 3$\ell$ channel. The uncertainty is decomposed into four components which are the signal modelling uncertainty, the background modelling uncertainty, the experimental uncertainty, and the data statistical uncertainty.

The absolute parton-level differential cross-section measured in the fiducial phase-space as a function of the $|\Delta y (Z, t_{\textrm{lep}})|$ in the 3$\ell$ channel. The uncertainty is decomposed into four components which are the signal modelling uncertainty, the background modelling uncertainty, the experimental uncertainty, and the data statistical uncertainty.

The absolute parton-level differential cross-section measured in the fiducial phase-space as a function of the $|\Delta \phi (l_{t}^{+}, l_{\bar{t}}^{-})|/\pi$ in the 4$\ell$ channel. The uncertainty is decomposed into four components which are the signal modelling uncertainty, the background modelling uncertainty, the experimental uncertainty, and the data statistical uncertainty.

The absolute parton-level differential cross-section measured in the fiducial phase-space as a function of the $|\Delta \phi (t\bar{t}, Z)|/\pi$ in the 4$\ell$ channel. The uncertainty is decomposed into four components which are the signal modelling uncertainty, the background modelling uncertainty, the experimental uncertainty, and the data statistical uncertainty.

The absolute parton-level differential cross-section measured in the fiducial phase-space as a function of the $p_{\textrm{T}}^{t\bar{t}}$ in the 4$\ell$ channel. The uncertainty is decomposed into four components which are the signal modelling uncertainty, the background modelling uncertainty, the experimental uncertainty, and the data statistical uncertainty.

The normalised parton-level differential cross-section measured in the fiducial phase-space as a function of the $p_{\textrm{T}}^{l \textrm{non-}Z}$ in the 3$\ell$ channel. The uncertainty is decomposed into four components which are the signal modelling uncertainty, the background modelling uncertainty, the experimental uncertainty, and the data statistical uncertainty.

The normalised parton-level differential cross-section measured in the fiducial phase-space as a function of the $|\Delta \phi (Z, t_{\textrm{lep}})|/\pi$ in the 3$\ell$ channel. The uncertainty is decomposed into four components which are the signal modelling uncertainty, the background modelling uncertainty, the experimental uncertainty, and the data statistical uncertainty.

The normalised parton-level differential cross-section measured in the fiducial phase-space as a function of the $|\Delta y (Z, t_{\textrm{lep}})|$ in the 3$\ell$ channel. The uncertainty is decomposed into four components which are the signal modelling uncertainty, the background modelling uncertainty, the experimental uncertainty, and the data statistical uncertainty.

The normalised parton-level differential cross-section measured in the fiducial phase-space as a function of the $|\Delta \phi (l_{t}^{+}, l_{\bar{t}}^{-})|/\pi$ in the 4$\ell$ channel. The uncertainty is decomposed into four components which are the signal modelling uncertainty, the background modelling uncertainty, the experimental uncertainty, and the data statistical uncertainty.

The normalised parton-level differential cross-section measured in the fiducial phase-space as a function of the $|\Delta \phi (t\bar{t}, Z)|/\pi$ in the 4$\ell$ channel. The uncertainty is decomposed into four components which are the signal modelling uncertainty, the background modelling uncertainty, the experimental uncertainty, and the data statistical uncertainty.

The normalised parton-level differential cross-section measured in the fiducial phase-space as a function of the $p_{\textrm{T}}^{t\bar{t}}$ in the 4$\ell$ channel. The uncertainty is decomposed into four components which are the signal modelling uncertainty, the background modelling uncertainty, the experimental uncertainty, and the data statistical uncertainty.

The normalised parton-level differential cross-section measured in the fiducial phase-space as a function of the absolute value of rapidity of the $Z$ boson in the 3$\ell$+4$\ell$ channel. The uncertainty is decomposed into four components which are the signal modelling uncertainty, the background modelling uncertainty, the experimental uncertainty, and the data statistical uncertainty.

The absolute particle-level differential cross-section measured in the fiducial phase-space as a function of the $p_{\textrm{T}}^{l \textrm{non-}Z}$ in the 3$\ell$ channel. The uncertainty is decomposed into four components which are the signal modelling uncertainty, the background modelling uncertainty, the experimental uncertainty, and the data statistical uncertainty.

The absolute particle-level differential cross-section measured in the fiducial phase-space as a function of the $|\Delta \phi (Z, t_{\textrm{lep}})|/\pi$ in the 3$\ell$ channel. The uncertainty is decomposed into four components which are the signal modelling uncertainty, the background modelling uncertainty, the experimental uncertainty, and the data statistical uncertainty.

The absolute particle-level differential cross-section measured in the fiducial phase-space as a function of the $|\Delta y (Z, t_{\textrm{lep}})|$ in the 3$\ell$ channel. The uncertainty is decomposed into four components which are the signal modelling uncertainty, the background modelling uncertainty, the experimental uncertainty, and the data statistical uncertainty.

The absolute particle-level differential cross-section measured in the fiducial phase-space as a function of the $|\Delta \phi (l_{t}^{+}, l_{\bar{t}}^{-})|/\pi$ in the 4$\ell$ channel. The uncertainty is decomposed into four components which are the signal modelling uncertainty, the background modelling uncertainty, the experimental uncertainty, and the data statistical uncertainty.

The absolute particle-level differential cross-section measured in the fiducial phase-space as a function of the $|\Delta \phi (t\bar{t}, Z)|/\pi$ in the 4$\ell$ channel. The uncertainty is decomposed into four components which are the signal modelling uncertainty, the background modelling uncertainty, the experimental uncertainty, and the data statistical uncertainty.

The absolute particle-level differential cross-section measured in the fiducial phase-space as a function of the $p_{\textrm{T}}^{t\bar{t}}$ in the 4$\ell$ channel. The uncertainty is decomposed into four components which are the signal modelling uncertainty, the background modelling uncertainty, the experimental uncertainty, and the data statistical uncertainty.

The absolute particle-level differential cross-section measured in the fiducial phase-space as a function of the absolute value of rapidity of the $Z$ boson in the 3$\ell$+4$\ell$ channel. The uncertainty is decomposed into four components which are the signal modelling uncertainty, the background modelling uncertainty, the experimental uncertainty, and the data statistical uncertainty.

The normalised particle-level differential cross-section measured in the fiducial phase-space as a function of the $p_{\textrm{T}}^{l \textrm{non-}Z}$ in the 3$\ell$ channel. The uncertainty is decomposed into four components which are the signal modelling uncertainty, the background modelling uncertainty, the experimental uncertainty, and the data statistical uncertainty.

The normalised particle-level differential cross-section measured in the fiducial phase-space as a function of the $|\Delta \phi (Z, t_{\textrm{lep}})|/\pi$ in the 3$\ell$ channel. The uncertainty is decomposed into four components which are the signal modelling uncertainty, the background modelling uncertainty, the experimental uncertainty, and the data statistical uncertainty.

The normalised particle-level differential cross-section measured in the fiducial phase-space as a function of the $|\Delta y (Z, t_{\textrm{lep}})|$ in the 3$\ell$ channel. The uncertainty is decomposed into four components which are the signal modelling uncertainty, the background modelling uncertainty, the experimental uncertainty, and the data statistical uncertainty.

The normalised particle-level differential cross-section measured in the fiducial phase-space as a function of the $|\Delta \phi (l_{t}^{+}, l_{\bar{t}}^{-})|/\pi$ in the 4$\ell$ channel. The uncertainty is decomposed into four components which are the signal modelling uncertainty, the background modelling uncertainty, the experimental uncertainty, and the data statistical uncertainty.

The normalised particle-level differential cross-section measured in the fiducial phase-space as a function of the $|\Delta \phi (t\bar{t}, Z)|/\pi$ in the 4$\ell$ channel. The uncertainty is decomposed into four components which are the signal modelling uncertainty, the background modelling uncertainty, the experimental uncertainty, and the data statistical uncertainty.

The normalised particle-level differential cross-section measured in the fiducial phase-space as a function of the $p_{\textrm{T}}^{t\bar{t}}$ in the 4$\ell$ channel. The uncertainty is decomposed into four components which are the signal modelling uncertainty, the background modelling uncertainty, the experimental uncertainty, and the data statistical uncertainty.

The normalised particle-level differential cross-section measured in the fiducial phase-space as a function of the absolute value of rapidity of the $Z$ boson in the 3$\ell$+4$\ell$ channel. The uncertainty is decomposed into four components which are the signal modelling uncertainty, the background modelling uncertainty, the experimental uncertainty, and the data statistical uncertainty.

The normalised particle-level differential cross-section measured in the fiducial phase-space as a function of the number of jets in the 3$\ell$ channel. The uncertainty is decomposed into four components which are the signal modelling uncertainty, the background modelling uncertainty, the experimental uncertainty, and the data statistical uncertainty.

The normalised particle-level differential cross-section measured in the fiducial phase-space as a function of the number of jets in the 4$\ell$ channel. The uncertainty is decomposed into four components which are the signal modelling uncertainty, the background modelling uncertainty, the experimental uncertainty, and the data statistical uncertainty.

The total correlation matrix of the absolute parton-level differential cross-section measured in the fiducial phase-space as a function of the $p_{\textrm{T}}$ of the $Z$ boson in the 3$\ell$+4$\ell$ channel.

The total correlation matrix of the absolute parton-level differential cross-section measured in the fiducial phase-space as a function of the absolute value of rapidity of the $Z$ boson in the 3$\ell$+4$\ell$ channel.

The total correlation matrix of the absolute parton-level differential cross-section measured in the fiducial phase-space as a function of the $p_{\textrm{T}}^{l \textrm{non-}Z}$ in the 3$\ell$ channel.

The total correlation matrix of the absolute parton-level differential cross-section measured in the fiducial phase-space as a function of the $|\Delta \phi (Z, t_{\textrm{lep}})|/\pi$ in the 3$\ell$ channel.

The total correlation matrix of the absolute parton-level differential cross-section measured in the fiducial phase-space as a function of the $|\Delta y (Z, t_{\textrm{lep}})|/\pi$ in the 3$\ell$ channel.

The total correlation matrix of the absolute parton-level differential cross-section measured in the fiducial phase-space as a function of the $|\Delta \phi (l_{t}^{+}, l_{\bar{t}}^{-})|/\pi$ in the 4$\ell$ channel.

The total correlation matrix of the absolute parton-level differential cross-section measured in the fiducial phase-space as a function of the $|\Delta \phi (t\bar{t}, Z)|/\pi$ in the 4$\ell$ channel.

The total correlation matrix of the absolute parton-level differential cross-section measured in the fiducial phase-space as a function of the $p_{\textrm{T}}^{t\bar{t}}$ in the 4$\ell$ channel.

The total correlation matrix of the normalised parton-level differential cross-section measured in the fiducial phase-space as a function of the $p_{\textrm{T}}$ of the $Z$ boson in the 3$\ell$+4$\ell$ channel.

The total correlation matrix of the normalised parton-level differential cross-section measured in the fiducial phase-space as a function of the absolute value of rapidity of the $Z$ boson in the 3$\ell$+4$\ell$ channel.

The total correlation matrix of the normalised parton-level differential cross-section measured in the fiducial phase-space as a function of the $p_{\textrm{T}}^{l \textrm{non-}Z}$ in the 3$\ell$ channel.

The total correlation matrix of the normalised parton-level differential cross-section measured in the fiducial phase-space as a function of the $|\Delta \phi (Z, t_{\textrm{lep}})|/\pi$ in the 3$\ell$ channel.

The total correlation matrix of the normalised parton-level differential cross-section measured in the fiducial phase-space as a function of the $|\Delta y (Z, t_{\textrm{lep}})|/\pi$ in the 3$\ell$ channel.

The total correlation matrix of the normalised parton-level differential cross-section measured in the fiducial phase-space as a function of the $|\Delta \phi (l_{t}^{+}, l_{\bar{t}}^{-})|/\pi$ in the 4$\ell$ channel.

The total correlation matrix of the normalised parton-level differential cross-section measured in the fiducial phase-space as a function of the $|\Delta \phi (t\bar{t}, Z)|/\pi$ in the 4$\ell$ channel.

The total correlation matrix of the normalised parton-level differential cross-section measured in the fiducial phase-space as a function of the $p_{\textrm{T}}^{t\bar{t}}$ in the 4$\ell$ channel.

The total correlation matrix of the absolute particle-level differential cross-section measured in the fiducial phase-space as a function of the $p_{\textrm{T}}$ of the $Z$ boson in the 3$\ell$+4$\ell$ channel.

The total correlation matrix of the absolute particle-level differential cross-section measured in the fiducial phase-space as a function of the absolute value of rapidity of the $Z$ boson in the 3$\ell$+4$\ell$ channel.

The total correlation matrix of the absolute particle-level differential cross-section measured in the fiducial phase-space as a function of the $p_{\textrm{T}}^{l \textrm{non-}Z}$ in the 3$\ell$ channel.

The total correlation matrix of the absolute particle-level differential cross-section measured in the fiducial phase-space as a function of the $|\Delta \phi (Z, t_{\textrm{lep}})|/\pi$ in the 3$\ell$ channel.

The total correlation matrix of the absolute particle-level differential cross-section measured in the fiducial phase-space as a function of the $|\Delta y (Z, t_{\textrm{lep}})|/\pi$ in the 3$\ell$ channel.

The total correlation matrix of the absolute particle-level differential cross-section measured in the fiducial phase-space as a function of the number of jets in the 3$\ell$ channel.

The total correlation matrix of the absolute particle-level differential cross-section measured in the fiducial phase-space as a function of the $|\Delta \phi (l_{t}^{+}, l_{\bar{t}}^{-})|/\pi$ in the 4$\ell$ channel.

The total correlation matrix of the absolute particle-level differential cross-section measured in the fiducial phase-space as a function of the $|\Delta \phi (t\bar{t}, Z)|/\pi$ in the 4$\ell$ channel.

The total correlation matrix of the absolute particle-level differential cross-section measured in the fiducial phase-space as a function of the $p_{\textrm{T}}^{t\bar{t}}$ in the 4$\ell$ channel.

The total correlation matrix of the absolute particle-level differential cross-section measured in the fiducial phase-space as a function of the number of jets in the 4$\ell$ channel.

The total correlation matrix of the normalised particle-level differential cross-section measured in the fiducial phase-space as a function of the $p_{\textrm{T}}$ of the $Z$ boson in the 3$\ell$+4$\ell$ channel.

The total correlation matrix of the normalised particle-level differential cross-section measured in the fiducial phase-space as a function of the absolute value of rapidity of the $Z$ boson in the 3$\ell$+4$\ell$ channel.

The total correlation matrix of the normalised particle-level differential cross-section measured in the fiducial phase-space as a function of the $p_{\textrm{T}}^{l \textrm{non-}Z}$ in the 3$\ell$ channel.

The total correlation matrix of the normalised particle-level differential cross-section measured in the fiducial phase-space as a function of the $|\Delta \phi (Z, t_{\textrm{lep}})|/\pi$ in the 3$\ell$ channel.

The total correlation matrix of the normalised particle-level differential cross-section measured in the fiducial phase-space as a function of the $|\Delta y (Z, t_{\textrm{lep}})|/\pi$ in the 3$\ell$ channel.

The total correlation matrix of the normalised particle-level differential cross-section measured in the fiducial phase-space as a function of the number of jets in the 3$\ell$ channel.

The total correlation matrix of the normalised particle-level differential cross-section measured in the fiducial phase-space as a function of the $|\Delta \phi (l_{t}^{+}, l_{\bar{t}}^{-})|/\pi$ in the 4$\ell$ channel.

The total correlation matrix of the normalised particle-level differential cross-section measured in the fiducial phase-space as a function of the $|\Delta \phi (t\bar{t}, Z)|/\pi$ in the 4$\ell$ channel.

The total correlation matrix of the normalised particle-level differential cross-section measured in the fiducial phase-space as a function of the $p_{\textrm{T}}^{t\bar{t}}$ in the 4$\ell$ channel.

The total correlation matrix of the normalised particle-level differential cross-section measured in the fiducial phase-space as a function of the number of jets in the 4$\ell$ channel.

The contribution from different systematic uncertainties to the measured $t\bar{t}t\bar{t}$ production cross section, grouped in categories.

The results of the fitted signal strength $\mu$ in the 1L/2LOS and 2LSS/3L combined channel.

The results of fitted inclusive ${t\bar{t}t\bar{t}}$ cross-section in the 1L/2LOS and 2LSS/3L combined channel.

Comparison between data and prediction for the distribution of the sum of the pseudo-continuous b-tagging score over the six jets with the highest score in the 1L,$\geq$9j,$\geq$3b region before the fit.

Comparison between data and prediction for the distribution of the sum of the pseudo-continuous b-tagging score over the six jets with the highest score in the 1L,$\geq$9j,$\geq$3b region after the fit.

Comparison between data and prediction for the distribution of the sum of the pseudo-continuous b-tagging score over the six jets with the highest score in the 2LOS,$\geq$7j,$\geq$3b region before the fit.

Comparison between data and prediction for the distribution of the sum of the pseudo-continuous b-tagging score over the six jets with the highest score in the 2LOS,$\geq$7j,$\geq$3b region after the fit.

Comparison between data and prediction for the distribution of the scalar sum of all jet and lepton pT in the event in the 1L,$\geq$8j,$\geq$3b region before the fit.

Comparison between data and prediction for the distribution of the scalar sum of all jet and lepton pT in the event in the 1L,$\geq$8j,$\geq$3b region after the fit.

Comparison between data and prediction for the distribution of the scalar sum of all jet and lepton pT in the event in the 2LOS,$\geq$6j,$\geq$3b region before the fit.

Comparison between data and prediction for the distribution of the scalar sum of all jet and lepton pT in the event in the 2LOS,$\geq$6j,$\geq$3b region after the fit.

Comparison between data and prediction for the distribution of number of jets in the 1L,$\geq$8j,$\geq$3b region before the fit.

Comparison between data and prediction for the distribution of number of jets in the 1L,$\geq$8j,$\geq$3b region after the fit.

Comparison between data and prediction for the distribution of number of jets in the 2LOS,$\geq$6j,$\geq$3b region before the fit.

Comparison between data and prediction for the distribution of number of jets in the 2LOS,$\geq$6j,$\geq$3b region after the fit.

Comparison between data and prediction for the distribution of b-jets multiplicity in the 1L,$\geq$8j,$\geq$3b region before the fit.

Comparison between data and prediction for the distribution of b-jets multiplicity in the 1L,$\geq$8j,$\geq$3b region after the fit.

Comparison between data and prediction for the distribution of b-jets multiplicity in the 2LOS,$\geq$6j,$\geq$3b region before the fit.

Comparison between data and prediction for the distribution of b-jets multiplicity in the 2LOS,$\geq$6j,$\geq$3b region after the fit.

Comparison between data and prediction for the distribution of the BDT score in the 1L,9j,4b signal region before the fit.

Comparison between data and prediction for the distribution of the BDT score in the 1L,9j,4b signal region after the fit.

Comparison between data and prediction for the distribution of the BDT score in the 1L,9j,$\geq$5b signal region before the fit.

Comparison between data and prediction for the distribution of the BDT score in the 1L,9j,$\geq$5b signal region after the fit.

Comparison between data and prediction for the distribution of the BDT score in the 1L,$\geq$10j,3bL signal region before the fit.

Comparison between data and prediction for the distribution of the BDT score in the 1L,$\geq$10j,3bL signal region after the fit.

Comparison between data and prediction for the distribution of the BDT score in the 1L,$\geq$10j,3bH signal region before the fit.

Comparison between data and prediction for the distribution of the BDT score in the 1L,$\geq$10j,3bH signal region after the fit.

Comparison between data and prediction for the distribution of the BDT score in the 1L,$\geq$10j,4b signal region before the fit.

Comparison between data and prediction for the distribution of the BDT score in the 1L,$\geq$10j,4b signal region after the fit.

Comparison between data and prediction for the distribution of the BDT score in the 1L,$\geq$10j,$\geq$5b signal region before the fit.

Comparison between data and prediction for the distribution of the BDT score in the 1L,$\geq$10j,$\geq$5b signal region after the fit.

Comparison between data and prediction for the distribution of the BDT score in the 2LOS,7j,$\geq$4b signal region before the fit.

Comparison between data and prediction for the distribution of the BDT score in the 2LOS,7j,$\geq$4b signal region after the fit.

Comparison between data and prediction for the distribution of the BDT score in the 2LOS,$\geq$8j,3bL signal region before the fit.

Comparison between data and prediction for the distribution of the BDT score in the 2LOS,$\geq$8j,3bL signal region after the fit.

Comparison between data and prediction for the distribution of the BDT score in the 2LOS,$\geq$8j,3bH signal region before the fit.

Comparison between data and prediction for the distribution of the BDT score in the 2LOS,$\geq$8j,3bH signal region after the fit.

Comparison between data and prediction for the distribution of the BDT score in the 2LOS,$\geq$8j,$\geq$4b signal region before the fit.

Comparison between data and prediction for the distribution of the BDT score in the 2LOS,$\geq$8j,$\geq$4b signal region after the fit.

Comparison between data and prediction for the distribution of the scalar sum of all jet and lepton pT in the event in the 1L,8j,3bV validation region before the fit.

Comparison between data and prediction for the distribution of the scalar sum of all jet and lepton pT in the event in the 1L,8j,3bV validation region after the fit.

Comparison between data and prediction for the distribution of the BDT score in the 1L,9j,3bV validation region before the fit.

Comparison between data and prediction for the distribution of the BDT score in the 1L,9j,3bV validation region after the fit.

Comparison between data and prediction for the distribution of the BDT score in the 1L,$\geq$10j,3bV validation region before the fit.

Comparison between data and prediction for the distribution of the BDT score in the 1L,$\geq$10j,3bV validation region after the fit.

Comparison between data and prediction for the distribution of the scalar sum of all jet and lepton pT in the event in the 2LOS,6j,3bV validation region before the fit.

Comparison between data and prediction for the distribution of the scalar sum of all jet and lepton pT in the event in the 2LOS,6j,3bV validation region after the fit.

Comparison between data and prediction for the distribution of the BDT score in the 2LOS,7j,3bV validation region before the fit.

Comparison between data and prediction for the distribution of the BDT score in the 2LOS,7j,3bV validation region after the fit.

Comparison between data and prediction for the distribution of the BDT score in the 2LOS,$\geq$8j,3bV validation region before the fit.

Comparison between data and prediction for the distribution of the BDT score in the 2LOS,$\geq$8j,3bV validation region after the fit.

Comparison between data and prediction for the distribution of the scalar sum of all jet and lepton pT in the event in the 1L,8j,3bL control region before the fit.

Comparison between data and prediction for the distribution of the scalar sum of all jet and lepton pT in the event in the 1L,8j,3bL control region after the fit.

Comparison between data and prediction for the distribution of the scalar sum of all jet and lepton pT in the event in the 1L,8j,3bH control region before the fit.

Comparison between data and prediction for the distribution of the scalar sum of all jet and lepton pT in the event in the 1L,8j,3bH control region after the fit.

Comparison between data and prediction for the distribution of the scalar sum of all jet and lepton pT in the event in the 1L,9j,3bL control region before the fit.

Comparison between data and prediction for the distribution of the scalar sum of all jet and lepton pT in the event in the 1L,9j,3bL control region after the fit.

Comparison between data and prediction for the distribution of the scalar sum of all jet and lepton pT in the event in the 1L,9j,3bH control region before the fit.

Comparison between data and prediction for the distribution of the scalar sum of all jet and lepton pT in the event in the 1L,9j,3bH control region after the fit.

Comparison between data and prediction for the distribution of the scalar sum of all jet and lepton pT in the event in the 1L,8j,4b control region before the fit.

Comparison between data and prediction for the distribution of the scalar sum of all jet and lepton pT in the event in the 1L,8j,4b control region after the fit.

Comparison between data and prediction for the distribution of the scalar sum of all jet and lepton pT in the event in the 1L,8j,$\geq$5b control region before the fit.

Comparison between data and prediction for the distribution of the scalar sum of all jet and lepton pT in the event in the 1L,8j,$\geq$5b control region after the fit.

Comparison between data and prediction for the distribution of the scalar sum of all jet and lepton pT in the event in the 2LOS,6j,3bL control region before the fit.

Comparison between data and prediction for the distribution of the scalar sum of all jet and lepton pT in the event in the 2LOS,6j,3bL control region after the fit.

Comparison between data and prediction for the distribution of the scalar sum of all jet and lepton pT in the event in the 2LOS,6j,3bH control region before the fit.

Comparison between data and prediction for the distribution of the scalar sum of all jet and lepton pT in the event in the 2LOS,6j,3bH control region after the fit.

Comparison between data and prediction for the distribution of the scalar sum of all jet and lepton pT in the event in the 2LOS,7j,3bL control region before the fit.

Comparison between data and prediction for the distribution of the scalar sum of all jet and lepton pT in the event in the 2LOS,7j,3bL control region after the fit.

Comparison between data and prediction for the distribution of the scalar sum of all jet and lepton pT in the event in the 2LOS,7j,3bH control region before the fit.

Comparison between data and prediction for the distribution of the scalar sum of all jet and lepton pT in the event in the 2LOS,7j,3bH control region after the fit.

Comparison between data and prediction for the distribution of the scalar sum of all jet and lepton pT in the event in the 2LOS,6j,$\geq$4b control region before the fit.

Comparison between data and prediction for the distribution of the scalar sum of all jet and lepton pT in the event in the 2LOS,6j,$\geq$4b control region after the fit.

95% CL exclusion limits on $\mathcal{B}(H\rightarrow aa \rightarrow b\bar{b}b\bar{b})$ for $m_a = 45$ GeV.

95% CL exclusion limits on $\mathcal{B}(H\rightarrow aa \rightarrow b\bar{b}b\bar{b})$ for $m_a = 55$ GeV.

The fraction of $a$ boson decays matched to reconstructed displaced vertices passing all vertex selections in signal MC.

The extrapolated signal selection efficiency as a function of $c\tau_{a}$.

Post-fit $m_{T}$ distribution in the SR 4J b-tag N-1 region. N-1 refers to all cuts except for the requirement on $m_T$ being applied. Uncertainties include statistical and systematic uncertainties. The value 9999 is used as a placeholder for infinity.

Post-fit $m_{T}$ distribution in the SR 6J b-veto N-1 region. N-1 refers to all cuts except for the requirement on $m_T$ being applied. Uncertainties include statistical and systematic uncertainties. The value 9999 is used as a placeholder for infinity.

Post-fit $m_{T}$ distribution in the SR 6J b-tag N-1 region. N-1 refers to all cuts except for the requirement on $m_T$ being applied. Uncertainties include statistical and systematic uncertainties. The value 9999 is used as a placeholder for infinity.

Post-fit $m_{eff}$ distribution in the 2J b-tag signal region. Uncertainties include statistical and systematic uncertainties. Including exemplary signal points. The value 9999 is used as a placeholder for infinity.

Post-fit $m_{eff}$ distribution in the 2J b-veto signal region. Uncertainties include statistical and systematic uncertainties. Including exemplary signal points. The value 9999 is used as a placeholder for infinity.

Post-fit $m_{eff}$ distribution in the 4J low-x b-tag signal region. Uncertainties include statistical and systematic uncertainties. Including exemplary signal points. The value 9999 is used as a placeholder for infinity.

Post-fit $m_{eff}$ distribution in the 4J low-x b-veto signal region. Uncertainties include statistical and systematic uncertainties. Including exemplary signal points. The value 9999 is used as a placeholder for infinity.

Post-fit $m_{eff}$ distribution in the 4J high-x b-tag signal region. Uncertainties include statistical and systematic uncertainties. Including exemplary signal points. The value 9999 is used as a placeholder for infinity.

Post-fit $m_{eff}$ distribution in the 4J high-x b-veto signal region. Uncertainties include statistical and systematic uncertainties. Including exemplary signal points. The value 9999 is used as a placeholder for infinity.

Post-fit $m_{eff}$ distribution in the 6J b-tag signal region. Uncertainties include statistical and systematic uncertainties. Including exemplary signal points. The value 9999 is used as a placeholder for infinity.

Post-fit $m_{eff}$ distribution in the 6J b-veto signal region. Uncertainties include statistical and systematic uncertainties. Including exemplary signal points. The value 9999 is used as a placeholder for infinity.

Observed 95% CL exclusion contours for the gluino one-step x = 1/2 model.

Expected 95% CL exclusion contours for the gluino one-step x = 1/2 model. space.

Observed 95% CL exclusion contours for the gluino one-step variable-x

Expected 95% CL exclusion contours for the gluino one-step variable-x

Observed 95% CL exclusion contours for the squark one-step x = 1/2 model.

Observed 95% CL exclusion contours for the squark one-step x = 1/2 model.

Observed 95% CL exclusion contours for one-flavour schemes in one-step x = 1/2 model.

Expected 95% CL exclusion contours for the squark one-step variable-x

Expected 95% CL exclusion contours for the squark one-step variable-x

Expected 95% CL exclusion contours for the squark one-flavour schemes in variable-x

Expected 95% CL exclusion contours for the squark one-flavour schemes in variable-x

Post-fit $m_{eff}$ distribution in the 2J b-tag validation region. Uncertainties include statistical and systematic uncertainties. The value 9999 is used as a placeholder for infinity.

Post-fit $m_{eff}$ distribution in the 2J b-veto validation region. Uncertainties include statistical and systematic uncertainties. The value 9999 is used as a placeholder for infinity.

Post-fit $m_{eff}$ distribution in the 4J b-tag validation region. Uncertainties include statistical and systematic uncertainties. The value 9999 is used as a placeholder for infinity.

Post-fit $m_{eff}$ distribution in the 4J b-veto validation region. Uncertainties include statistical and systematic uncertainties. The value 9999 is used as a placeholder for infinity.

Post-fit $m_{eff}$ distribution in the 6J b-tag validation region. Uncertainties include statistical and systematic uncertainties.

Post-fit $m_{eff}$ distribution in the 6J b-veto validation region. Uncertainties include statistical and systematic uncertainties.

Event selection cutflow for two representative signal samples for the SR2JBT. The gluino, squark, chargino and neutralino masses are reported. Weighted events including statistical uncertainties are shown.

Event selection cutflow for two representative signal samples for the SR2JBV. The gluino, squark, chargino and neutralino masses are reported. Weighted events including statistical uncertainties are shown.

Event selection cutflow for two representative signal samples for the SR4JBT. The gluino, squark, chargino and neutralino masses are reported. Weighted events including statistical uncertainties are shown.

Event selection cutflow for two representative signal samples for the SR4JBV. The gluino, squark, chargino and neutralino masses are reported. Weighted events including statistical uncertainties are shown.

Event selection cutflow for two representative signal samples for the SR6JBT. The gluino, squark, chargino and neutralino masses are reported. Weighted events including statistical uncertainties are shown.

Event selection cutflow for two representative signal samples for the SR6JBV. The gluino, squark, chargino and neutralino masses are reported. Weighted events including statistical uncertainties are shown.

Signal acceptance in SR2J b-Tag bin1 region for gluino production one-step x = 1/2 simplified models

Signal acceptance in SR2J b-Tag bin2 region for gluino production one-step x = 1/2 simplified models

Signal acceptance in SR2J b-Tag bin3 region for gluino production one-step x = 1/2 simplified models

Signal acceptance in SR2J b-Veto bin1 region for gluino production one-step x = 1/2 simplified models

Signal acceptance in SR2J b-Veto bin2 region for gluino production one-step x = 1/2 simplified models

Signal acceptance in SR2J b-Veto bin3 region for gluino production one-step x = 1/2 simplified models

Signal acceptance in SR2J discovery high region for gluino production one-step x = 1/2 simplified models

Signal acceptance in SR2J discovery low region for gluino production one-step x = 1/2 simplified models

Signal acceptance in SR4Jhx discovery region for gluino production one-step x = 1/2 simplified models

Signal acceptance in SR4Jhx b-Tag bin1 region for gluino production one-step x = 1/2 simplified models

Signal acceptance in SR4Jhx b-Tag bin2 region for gluino production one-step x = 1/2 simplified models

Signal acceptance in SR4Jhx b-Tag bin3 region for gluino production one-step x = 1/2 simplified models

Signal acceptance in SR4Jhx b-Veto bin1 region for gluino production one-step x = 1/2 simplified models

Signal acceptance in SR4Jhx b-Veto bin2 region for gluino production one-step x = 1/2 simplified models

Signal acceptance in SR4Jhx b-Veto bin3 region for gluino production one-step x = 1/2 simplified models

Signal acceptance in SR4Jlx discovery region for gluino production one-step x = 1/2 simplified models

Signal acceptance in SR4Jlx b-Tag bin1 region for gluino production one-step x = 1/2 simplified models

Signal acceptance in SR4Jlx b-Tag bin2 region for gluino production one-step x = 1/2 simplified models

Signal acceptance in SR4Jlx b-Tag bin3 region for gluino production one-step x = 1/2 simplified models

Signal acceptance in SR4Jlx b-Veto bin1 region for gluino production one-step x = 1/2 simplified models

Signal acceptance in SR4Jlx b-Veto bin2 region for gluino production one-step x = 1/2 simplified models

Signal acceptance in SR4Jlx b-Veto bin3 region for gluino production one-step x = 1/2 simplified models

Signal acceptance in SR6J b-Tag bin1 region for gluino production one-step x = 1/2 simplified models

Signal acceptance in SR6J b-Tag bin2 region for gluino production one-step x = 1/2 simplified models

Signal acceptance in SR6J b-Tag bin3 region for gluino production one-step x = 1/2 simplified models

Signal acceptance in SR6J b-Tag bin4 region for gluino production one-step x = 1/2 simplified models

Signal acceptance in SR6J b-Veto bin1 region for gluino production one-step x = 1/2 simplified models

Signal acceptance in SR6J b-Veto bin2 region for gluino production one-step x = 1/2 simplified models

Signal acceptance in SR6J b-Veto bin3 region for gluino production one-step x = 1/2 simplified models

Signal acceptance in SR6J b-Veto bin4 region for gluino production one-step x = 1/2 simplified models

Signal acceptance in SR6J discovery high region for gluino production one-step x = 1/2 simplified models

Signal acceptance in SR6J discovery low region for gluino production one-step x = 1/2 simplified models

Signal acceptance in SR2J b-Tag bin1 region for gluino production one-step variable-x simplified models

Signal acceptance in SR2J b-Tag bin2 region for gluino production one-step variable-x simplified models

Signal acceptance in SR2J b-Tag bin3 region for gluino production one-step variable-x simplified models

Signal acceptance in SR2J b-Veto bin1 region for gluino production one-step variable-x simplified models

Signal acceptance in SR2J b-Veto bin2 region for gluino production one-step variable-x simplified models

Signal acceptance in SR2J b-Veto bin3 region for gluino production one-step variable-x simplified models

Signal acceptance in SR2J discovery high region for gluino production one-step variable-x simplified models

Signal acceptance in SR2J discovery low region for gluino production one-step variable-x simplified models

Signal acceptance in SR4Jhx discovery region for gluino production one-step variable-x simplified models

Signal acceptance in SR4Jhx b-Tag bin1 region for gluino production one-step variable-x simplified models

Signal acceptance in SR4Jhx b-Tag bin2 region for gluino production one-step variable-x simplified models

Signal acceptance in SR4Jhx b-Tag bin3 region for gluino production one-step variable-x simplified models

Signal acceptance in SR4Jhx b-Veto bin1 region for gluino production one-step variable-x simplified models

Signal acceptance in SR4Jhx b-Veto bin2 region for gluino production one-step variable-x simplified models

Signal acceptance in SR4Jhx b-Veto bin3 region for gluino production one-step variable-x simplified models

Signal acceptance in SR4Jlx discovery region for gluino production one-step variable-x simplified models

Signal acceptance in SR4Jlx b-Tag bin1 region for gluino production one-step variable-x simplified models

Signal acceptance in SR4Jlx b-Tag bin2 region for gluino production one-step variable-x simplified models

Signal acceptance in SR4Jlx b-Tag bin3 region for gluino production one-step variable-x simplified models

Signal acceptance in SR4Jlx b-Veto bin1 region for gluino production one-step variable-x simplified models

Signal acceptance in SR4Jlx b-Veto bin2 region for gluino production one-step variable-x simplified models

Signal acceptance in SR4Jlx b-Veto bin3 region for gluino production one-step variable-x simplified models

Signal acceptance in SR6J b-Tag bin1 region for gluino production one-step variable-x simplified models

Signal acceptance in SR6J b-Tag bin2 region for gluino production one-step variable-x simplified models

Signal acceptance in SR6J b-Tag bin3 region for gluino production one-step variable-x simplified models

Signal acceptance in SR6J b-Tag bin4 region for gluino production one-step variable-x simplified models

Signal acceptance in SR6J b-Veto bin1 region for gluino production one-step variable-x simplified models

Signal acceptance in SR6J b-Veto bin2 region for gluino production one-step variable-x simplified models

Signal acceptance in SR6J b-Veto bin3 region for gluino production one-step variable-x simplified models

Signal acceptance in SR6J b-Veto bin4 region for gluino production one-step variable-x simplified models

Signal acceptance in SR6J discovery high region for gluino production one-step variable-x simplified models

Signal acceptance in SR6J discovery low region for gluino production one-step variable-x simplified models

Signal acceptance in SR2J b-Tag bin1 region for squark production one-step x = 1/2 simplified models

Signal acceptance in SR2J b-Tag bin2 region for squark production one-step x = 1/2 simplified models

Signal acceptance in SR2J b-Tag bin3 region for squark production one-step x = 1/2 simplified models

Signal acceptance in SR2J b-Veto bin1 region for squark production one-step x = 1/2 simplified models

Signal acceptance in SR2J b-Veto bin2 region for squark production one-step x = 1/2 simplified models

Signal acceptance in SR2J b-Veto bin3 region for squark production one-step x = 1/2 simplified models

Signal acceptance in SR2J discovery high region for squark production one-step x = 1/2 simplified models

Signal acceptance in SR2J discovery low region for squark production one-step x = 1/2 simplified models

Signal acceptance in SR4Jhx discovery region for squark production one-step x = 1/2 simplified models

Signal acceptance in SR4Jhx b-Tag bin1 region for squark production one-step x = 1/2 simplified models

Signal acceptance in SR4Jhx b-Tag bin2 region for squark production one-step x = 1/2 simplified models

Signal acceptance in SR4Jhx b-Tag bin3 region for squark production one-step x = 1/2 simplified models

Signal acceptance in SR4Jhx b-Veto bin1 region for squark production one-step x = 1/2 simplified models

Signal acceptance in SR4Jhx b-Veto bin2 region for squark production one-step x = 1/2 simplified models

Signal acceptance in SR4Jhx b-Veto bin3 region for squark production one-step x = 1/2 simplified models

Signal acceptance in SR4Jlx discovery region for squark production one-step x = 1/2 simplified models

Signal acceptance in SR4Jlx b-Tag bin1 region for squark production one-step x = 1/2 simplified models

Signal acceptance in SR4Jlx b-Tag bin2 region for squark production one-step x = 1/2 simplified models

Signal acceptance in SR4Jlx b-Tag bin3 region for squark production one-step x = 1/2 simplified models

Signal acceptance in SR4Jlx b-Veto bin1 region for squark production one-step x = 1/2 simplified models

Signal acceptance in SR4Jlx b-Veto bin2 region for squark production one-step x = 1/2 simplified models

Signal acceptance in SR4Jlx b-Veto bin3 region for squark production one-step x = 1/2 simplified models

Signal acceptance in SR6J b-Tag bin1 region for squark production one-step x = 1/2 simplified models

Signal acceptance in SR6J b-Tag bin2 region for squark production one-step x = 1/2 simplified models

Signal acceptance in SR6J b-Tag bin3 region for squark production one-step x = 1/2 simplified models

Signal acceptance in SR6J b-Tag bin4 region for squark production one-step x = 1/2 simplified models

Signal acceptance in SR6J b-Veto bin1 region for squark production one-step x = 1/2 simplified models

Signal acceptance in SR6J b-Veto bin2 region for squark production one-step x = 1/2 simplified models

Signal acceptance in SR6J b-Veto bin3 region for squark production one-step x = 1/2 simplified models

Signal acceptance in SR6J b-Veto bin4 region for squark production one-step x = 1/2 simplified models

Signal acceptance in SR6J discovery high region for squark production one-step x = 1/2 simplified models

Signal acceptance in SR6J discovery low region for squark production one-step x = 1/2 simplified models

Signal acceptance in SR2J b-Tag bin1 region for squark production one-step variable-x simplified models

Signal acceptance in SR2J b-Tag bin2 region for squark production one-step variable-x simplified models

Signal acceptance in SR2J b-Tag bin3 region for squark production one-step variable-x simplified models

Signal acceptance in SR2J b-Veto bin1 region for squark production one-step variable-x simplified models

Signal acceptance in SR2J b-Veto bin2 region for squark production one-step variable-x simplified models

Signal acceptance in SR2J b-Veto bin3 region for squark production one-step variable-x simplified models

Signal acceptance in SR2J discovery high region for squark production one-step variable-x simplified models

Signal acceptance in SR2J discovery low region for squark production one-step variable-x simplified models

Signal acceptance in SR4Jhx discovery region for squark production one-step variable-x simplified models

Signal acceptance in SR4Jhx b-Tag bin1 region for squark production one-step variable-x simplified models

Signal acceptance in SR4Jhx b-Tag bin2 region for squark production one-step variable-x simplified models

Signal acceptance in SR4Jhx b-Tag bin3 region for squark production one-step variable-x simplified models

Signal acceptance in SR4Jhx b-Veto bin1 region for squark production one-step variable-x simplified models

Signal acceptance in SR4Jhx b-Veto bin2 region for squark production one-step variable-x simplified models

Signal acceptance in SR4Jhx b-Veto bin3 region for squark production one-step variable-x simplified models

Signal acceptance in SR4Jlx discovery region for squark production one-step variable-x simplified models

Signal acceptance in SR4Jlx b-Tag bin1 region for squark production one-step variable-x simplified models

Signal acceptance in SR4Jlx b-Tag bin2 region for squark production one-step variable-x simplified models

Signal acceptance in SR4Jlx b-Tag bin3 region for squark production one-step variable-x simplified models

Signal acceptance in SR4Jlx b-Veto bin1 region for squark production one-step variable-x simplified models

Signal acceptance in SR4Jlx b-Veto bin2 region for squark production one-step variable-x simplified models

Signal acceptance in SR4Jlx b-Veto bin3 region for squark production one-step variable-x simplified models

Signal acceptance in SR6J b-Tag bin1 region for squark production one-step variable-x simplified models

Signal acceptance in SR6J b-Tag bin2 region for squark production one-step variable-x simplified models

Signal acceptance in SR6J b-Tag bin3 region for squark production one-step variable-x simplified models

Signal acceptance in SR6J b-Tag bin4 region for squark production one-step variable-x simplified models

Signal acceptance in SR6J b-Veto bin1 region for squark production one-step variable-x simplified models

Signal acceptance in SR6J b-Veto bin2 region for squark production one-step variable-x simplified models

Signal acceptance in SR6J b-Veto bin3 region for squark production one-step variable-x simplified models

Signal acceptance in SR6J b-Veto bin4 region for squark production one-step variable-x simplified models

Signal acceptance in SR6J discovery high region for squark production one-step variable-x simplified models

Signal acceptance in SR6J discovery low region for squark production one-step variable-x simplified models

Signal efficiency in SR2J b-Tag bin1 region for gluino production one-step x = 1/2 simplified models. The -1 value indicates the truth yields for this point is 0 but the reco yields is not 0

Signal efficiency in SR2J b-Tag bin2 region for gluino production one-step x = 1/2 simplified models. The -1 value indicates the truth yields for this point is 0 but the reco yields is not 0

Signal efficiency in SR2J b-Tag bin3 region for gluino production one-step x = 1/2 simplified models. The -1 value indicates the truth yields for this point is 0 but the reco yields is not 0

Signal efficiency in SR2J b-Veto bin1 region for gluino production one-step x = 1/2 simplified models. The -1 value indicates the truth yields for this point is 0 but the reco yields is not 0

Signal efficiency in SR2J b-Veto bin2 region for gluino production one-step x = 1/2 simplified models. The -1 value indicates the truth yields for this point is 0 but the reco yields is not 0

Signal efficiency in SR2J b-Veto bin3 region for gluino production one-step x = 1/2 simplified models. The -1 value indicates the truth yields for this point is 0 but the reco yields is not 0

Signal efficiency in SR2J discovery high region for gluino production one-step x = 1/2 simplified models. The -1 value indicates the truth yields for this point is 0 but the reco yields is not 0

Signal efficiency in SR2J discovery low region for gluino production one-step x = 1/2 simplified models. The -1 value indicates the truth yields for this point is 0 but the reco yields is not 0

Signal efficiency in SR4Jhx discovery region for gluino production one-step x = 1/2 simplified models. The -1 value indicates the truth yields for this point is 0 but the reco yields is not 0

Signal efficiency in SR4Jhx b-Tag bin1 region for gluino production one-step x = 1/2 simplified models. The -1 value indicates the truth yields for this point is 0 but the reco yields is not 0

Signal efficiency in SR4Jhx b-Tag bin2 region for gluino production one-step x = 1/2 simplified models. The -1 value indicates the truth yields for this point is 0 but the reco yields is not 0

Signal efficiency in SR4Jhx b-Tag bin3 region for gluino production one-step x = 1/2 simplified models. The -1 value indicates the truth yields for this point is 0 but the reco yields is not 0

Signal efficiency in SR4Jhx b-Veto bin1 region for gluino production one-step x = 1/2 simplified models. The -1 value indicates the truth yields for this point is 0 but the reco yields is not 0

Signal efficiency in SR4Jhx b-Veto bin2 region for gluino production one-step x = 1/2 simplified models. The -1 value indicates the truth yields for this point is 0 but the reco yields is not 0

Signal efficiency in SR4Jhx b-Veto bin3 region for gluino production one-step x = 1/2 simplified models. The -1 value indicates the truth yields for this point is 0 but the reco yields is not 0

Signal efficiency in SR4Jlx discovery region for gluino production one-step x = 1/2 simplified models. The -1 value indicates the truth yields for this point is 0 but the reco yields is not 0

Signal efficiency in SR4Jlx b-Tag bin1 region for gluino production one-step x = 1/2 simplified models. The -1 value indicates the truth yields for this point is 0 but the reco yields is not 0

Signal efficiency in SR4Jlx b-Tag bin2 region for gluino production one-step x = 1/2 simplified models. The -1 value indicates the truth yields for this point is 0 but the reco yields is not 0

Signal efficiency in SR4Jlx b-Tag bin3 region for gluino production one-step x = 1/2 simplified models. The -1 value indicates the truth yields for this point is 0 but the reco yields is not 0

Signal efficiency in SR4Jlx b-Veto bin1 region for gluino production one-step x = 1/2 simplified models. The -1 value indicates the truth yields for this point is 0 but the reco yields is not 0

Signal efficiency in SR4Jlx b-Veto bin2 region for gluino production one-step x = 1/2 simplified models. The -1 value indicates the truth yields for this point is 0 but the reco yields is not 0

Signal efficiency in SR4Jlx b-Veto bin3 region for gluino production one-step x = 1/2 simplified models. The -1 value indicates the truth yields for this point is 0 but the reco yields is not 0

Signal efficiency in SR6J b-Tag bin1 region for gluino production one-step x = 1/2 simplified models. The -1 value indicates the truth yields for this point is 0 but the reco yields is not 0

Signal efficiency in SR6J b-Tag bin2 region for gluino production one-step x = 1/2 simplified models. The -1 value indicates the truth yields for this point is 0 but the reco yields is not 0

Signal efficiency in SR6J b-Tag bin3 region for gluino production one-step x = 1/2 simplified models. The -1 value indicates the truth yields for this point is 0 but the reco yields is not 0

Signal efficiency in SR6J b-Tag bin4 region for gluino production one-step x = 1/2 simplified models. The -1 value indicates the truth yields for this point is 0 but the reco yields is not 0

Signal efficiency in SR6J b-Veto bin1 region for gluino production one-step x = 1/2 simplified models. The -1 value indicates the truth yields for this point is 0 but the reco yields is not 0

Signal efficiency in SR6J b-Veto bin2 region for gluino production one-step x = 1/2 simplified models. The -1 value indicates the truth yields for this point is 0 but the reco yields is not 0

Signal efficiency in SR6J b-Veto bin3 region for gluino production one-step x = 1/2 simplified models. The -1 value indicates the truth yields for this point is 0 but the reco yields is not 0

Signal efficiency in SR6J b-Veto bin4 region for gluino production one-step x = 1/2 simplified models. The -1 value indicates the truth yields for this point is 0 but the reco yields is not 0

Signal efficiency in SR6J discovery high region for gluino production one-step x = 1/2 simplified models. The -1 value indicates the truth yields for this point is 0 but the reco yields is not 0

Signal efficiency in SR6J discovery low region for gluino production one-step x = 1/2 simplified models. The -1 value indicates the truth yields for this point is 0 but the reco yields is not 0

Signal efficiency in SR2J b-Tag bin1 region for gluino production one-step variable-x simplified models. The -1 value indicates the truth yields for this point is 0 but the reco yields is not 0

Signal efficiency in SR2J b-Tag bin2 region for gluino production one-step variable-x simplified models. The -1 value indicates the truth yields for this point is 0 but the reco yields is not 0

Signal efficiency in SR2J b-Tag bin3 region for gluino production one-step variable-x simplified models. The -1 value indicates the truth yields for this point is 0 but the reco yields is not 0

Signal efficiency in SR2J b-Veto bin1 region for gluino production one-step variable-x simplified models. The -1 value indicates the truth yields for this point is 0 but the reco yields is not 0

Signal efficiency in SR2J b-Veto bin2 region for gluino production one-step variable-x simplified models. The -1 value indicates the truth yields for this point is 0 but the reco yields is not 0

Signal efficiency in SR2J b-Veto bin3 region for gluino production one-step variable-x simplified models. The -1 value indicates the truth yields for this point is 0 but the reco yields is not 0

Signal efficiency in SR2J discovery high region for gluino production one-step variable-x simplified models. The -1 value indicates the truth yields for this point is 0 but the reco yields is not 0

Signal efficiency in SR2J discovery low region for gluino production one-step variable-x simplified models. The -1 value indicates the truth yields for this point is 0 but the reco yields is not 0

Signal efficiency in SR4Jhx discovery region for gluino production one-step variable-x simplified models. The -1 value indicates the truth yields for this point is 0 but the reco yields is not 0

Signal efficiency in SR4Jhx b-Tag bin1 region for gluino production one-step variable-x simplified models. The -1 value indicates the truth yields for this point is 0 but the reco yields is not 0

Signal efficiency in SR4Jhx b-Tag bin2 region for gluino production one-step variable-x simplified models. The -1 value indicates the truth yields for this point is 0 but the reco yields is not 0

Signal efficiency in SR4Jhx b-Tag bin3 region for gluino production one-step variable-x simplified models. The -1 value indicates the truth yields for this point is 0 but the reco yields is not 0

Signal efficiency in SR4Jhx b-Veto bin1 region for gluino production one-step variable-x simplified models. The -1 value indicates the truth yields for this point is 0 but the reco yields is not 0

Signal efficiency in SR4Jhx b-Veto bin2 region for gluino production one-step variable-x simplified models. The -1 value indicates the truth yields for this point is 0 but the reco yields is not 0

Signal efficiency in SR4Jhx b-Veto bin3 region for gluino production one-step variable-x simplified models. The -1 value indicates the truth yields for this point is 0 but the reco yields is not 0

Signal efficiency in SR4Jlx discovery region for gluino production one-step variable-x simplified models. The -1 value indicates the truth yields for this point is 0 but the reco yields is not 0

Signal efficiency in SR4Jlx b-Tag bin1 region for gluino production one-step variable-x simplified models. The -1 value indicates the truth yields for this point is 0 but the reco yields is not 0

Signal efficiency in SR4Jlx b-Tag bin2 region for gluino production one-step variable-x simplified models. The -1 value indicates the truth yields for this point is 0 but the reco yields is not 0

Signal efficiency in SR4Jlx b-Tag bin3 region for gluino production one-step variable-x simplified models. The -1 value indicates the truth yields for this point is 0 but the reco yields is not 0

Signal efficiency in SR4Jlx b-Veto bin1 region for gluino production one-step variable-x simplified models. The -1 value indicates the truth yields for this point is 0 but the reco yields is not 0

Signal efficiency in SR4Jlx b-Veto bin2 region for gluino production one-step variable-x simplified models. The -1 value indicates the truth yields for this point is 0 but the reco yields is not 0

Signal efficiency in SR4Jlx b-Veto bin3 region for gluino production one-step variable-x simplified models. The -1 value indicates the truth yields for this point is 0 but the reco yields is not 0

Signal efficiency in SR6J b-Tag bin1 region for gluino production one-step variable-x simplified models. The -1 value indicates the truth yields for this point is 0 but the reco yields is not 0

Signal efficiency in SR6J b-Tag bin2 region for gluino production one-step variable-x simplified models. The -1 value indicates the truth yields for this point is 0 but the reco yields is not 0

Signal efficiency in SR6J b-Tag bin3 region for gluino production one-step variable-x simplified models. The -1 value indicates the truth yields for this point is 0 but the reco yields is not 0

Signal efficiency in SR6J b-Tag bin4 region for gluino production one-step variable-x simplified models. The -1 value indicates the truth yields for this point is 0 but the reco yields is not 0

Signal efficiency in SR6J b-Veto bin1 region for gluino production one-step variable-x simplified models. The -1 value indicates the truth yields for this point is 0 but the reco yields is not 0

Signal efficiency in SR6J b-Veto bin2 region for gluino production one-step variable-x simplified models. The -1 value indicates the truth yields for this point is 0 but the reco yields is not 0

Signal efficiency in SR6J b-Veto bin3 region for gluino production one-step variable-x simplified models. The -1 value indicates the truth yields for this point is 0 but the reco yields is not 0

Signal efficiency in SR6J b-Veto bin4 region for gluino production one-step variable-x simplified models. The -1 value indicates the truth yields for this point is 0 but the reco yields is not 0

Signal efficiency in SR6J discovery high region for gluino production one-step variable-x simplified models. The -1 value indicates the truth yields for this point is 0 but the reco yields is not 0

Signal efficiency in SR6J discovery low region for gluino production one-step variable-x simplified models. The -1 value indicates the truth yields for this point is 0 but the reco yields is not 0

Signal efficiency in SR2J b-Tag bin1 region for squark production one-step x = 1/2 simplified models. The -1 value indicates the truth yields for this point is 0 but the reco yields is not 0

Signal efficiency in SR2J b-Tag bin2 region for squark production one-step x = 1/2 simplified models. The -1 value indicates the truth yields for this point is 0 but the reco yields is not 0

Signal efficiency in SR2J b-Tag bin3 region for squark production one-step x = 1/2 simplified models. The -1 value indicates the truth yields for this point is 0 but the reco yields is not 0

Signal efficiency in SR2J b-Veto bin1 region for squark production one-step x = 1/2 simplified models. The -1 value indicates the truth yields for this point is 0 but the reco yields is not 0

Signal efficiency in SR2J b-Veto bin2 region for squark production one-step x = 1/2 simplified models. The -1 value indicates the truth yields for this point is 0 but the reco yields is not 0

Signal efficiency in SR2J b-Veto bin3 region for squark production one-step x = 1/2 simplified models. The -1 value indicates the truth yields for this point is 0 but the reco yields is not 0

Signal efficiency in SR2J discovery high region for squark production one-step x = 1/2 simplified models. The -1 value indicates the truth yields for this point is 0 but the reco yields is not 0

Signal efficiency in SR2J discovery low region for squark production one-step x = 1/2 simplified models. The -1 value indicates the truth yields for this point is 0 but the reco yields is not 0

Signal efficiency in SR4Jhx discovery region for squark production one-step x = 1/2 simplified models. The -1 value indicates the truth yields for this point is 0 but the reco yields is not 0

Signal efficiency in SR4Jhx b-Tag bin1 region for squark production one-step x = 1/2 simplified models. The -1 value indicates the truth yields for this point is 0 but the reco yields is not 0

Signal efficiency in SR4Jhx b-Tag bin2 region for squark production one-step x = 1/2 simplified models. The -1 value indicates the truth yields for this point is 0 but the reco yields is not 0

Signal efficiency in SR4Jhx b-Tag bin3 region for squark production one-step x = 1/2 simplified models. The -1 value indicates the truth yields for this point is 0 but the reco yields is not 0

Signal efficiency in SR4Jhx b-Veto bin1 region for squark production one-step x = 1/2 simplified models. The -1 value indicates the truth yields for this point is 0 but the reco yields is not 0

Signal efficiency in SR4Jhx b-Veto bin2 region for squark production one-step x = 1/2 simplified models. The -1 value indicates the truth yields for this point is 0 but the reco yields is not 0

Signal efficiency in SR4Jhx b-Veto bin3 region for squark production one-step x = 1/2 simplified models. The -1 value indicates the truth yields for this point is 0 but the reco yields is not 0

Signal efficiency in SR4Jlx discovery region for squark production one-step x = 1/2 simplified models. The -1 value indicates the truth yields for this point is 0 but the reco yields is not 0

Signal efficiency in SR4Jlx b-Tag bin1 region for squark production one-step x = 1/2 simplified models. The -1 value indicates the truth yields for this point is 0 but the reco yields is not 0

Signal efficiency in SR4Jlx b-Tag bin2 region for squark production one-step x = 1/2 simplified models. The -1 value indicates the truth yields for this point is 0 but the reco yields is not 0

Signal efficiency in SR4Jlx b-Tag bin3 region for squark production one-step x = 1/2 simplified models. The -1 value indicates the truth yields for this point is 0 but the reco yields is not 0

Signal efficiency in SR4Jlx b-Veto bin1 region for squark production one-step x = 1/2 simplified models. The -1 value indicates the truth yields for this point is 0 but the reco yields is not 0

Signal efficiency in SR4Jlx b-Veto bin2 region for squark production one-step x = 1/2 simplified models. The -1 value indicates the truth yields for this point is 0 but the reco yields is not 0

Signal efficiency in SR4Jlx b-Veto bin3 region for squark production one-step x = 1/2 simplified models. The -1 value indicates the truth yields for this point is 0 but the reco yields is not 0

Signal efficiency in SR6J b-Tag bin1 region for squark production one-step x = 1/2 simplified models. The -1 value indicates the truth yields for this point is 0 but the reco yields is not 0

Signal efficiency in SR6J b-Tag bin2 region for squark production one-step x = 1/2 simplified models. The -1 value indicates the truth yields for this point is 0 but the reco yields is not 0

Signal efficiency in SR6J b-Tag bin3 region for squark production one-step x = 1/2 simplified models. The -1 value indicates the truth yields for this point is 0 but the reco yields is not 0

Signal efficiency in SR6J b-Tag bin4 region for squark production one-step x = 1/2 simplified models. The -1 value indicates the truth yields for this point is 0 but the reco yields is not 0

Signal efficiency in SR6J b-Veto bin1 region for squark production one-step x = 1/2 simplified models. The -1 value indicates the truth yields for this point is 0 but the reco yields is not 0

Signal efficiency in SR6J b-Veto bin2 region for squark production one-step x = 1/2 simplified models. The -1 value indicates the truth yields for this point is 0 but the reco yields is not 0

Signal efficiency in SR6J b-Veto bin3 region for squark production one-step x = 1/2 simplified models. The -1 value indicates the truth yields for this point is 0 but the reco yields is not 0

Signal efficiency in SR6J b-Veto bin4 region for squark production one-step x = 1/2 simplified models. The -1 value indicates the truth yields for this point is 0 but the reco yields is not 0

Signal efficiency in SR6J discovery high region for squark production one-step x = 1/2 simplified models. The -1 value indicates the truth yields for this point is 0 but the reco yields is not 0

Signal efficiency in SR6J discovery low region for squark production one-step x = 1/2 simplified models. The -1 value indicates the truth yields for this point is 0 but the reco yields is not 0

Signal efficiency in SR2J b-Tag bin1 region for squark production one-step variable-x simplified models. The -1 value indicates the truth yields for this point is 0 but the reco yields is not 0

Signal efficiency in SR2J b-Tag bin2 region for squark production one-step variable-x simplified models. The -1 value indicates the truth yields for this point is 0 but the reco yields is not 0

Signal efficiency in SR2J b-Tag bin3 region for squark production one-step variable-x simplified models. The -1 value indicates the truth yields for this point is 0 but the reco yields is not 0

Signal efficiency in SR2J b-Veto bin1 region for squark production one-step variable-x simplified models. The -1 value indicates the truth yields for this point is 0 but the reco yields is not 0

Signal efficiency in SR2J b-Veto bin2 region for squark production one-step variable-x simplified models. The -1 value indicates the truth yields for this point is 0 but the reco yields is not 0

Signal efficiency in SR2J b-Veto bin3 region for squark production one-step variable-x simplified models. The -1 value indicates the truth yields for this point is 0 but the reco yields is not 0

Signal efficiency in SR2J discovery high region for squark production one-step variable-x simplified models. The -1 value indicates the truth yields for this point is 0 but the reco yields is not 0

Signal efficiency in SR2J discovery low region for squark production one-step variable-x simplified models. The -1 value indicates the truth yields for this point is 0 but the reco yields is not 0

Signal efficiency in SR4Jhx discovery region for squark production one-step variable-x simplified models. The -1 value indicates the truth yields for this point is 0 but the reco yields is not 0

Signal efficiency in SR4Jhx b-Tag bin1 region for squark production one-step variable-x simplified models. The -1 value indicates the truth yields for this point is 0 but the reco yields is not 0

Signal efficiency in SR4Jhx b-Tag bin2 region for squark production one-step variable-x simplified models. The -1 value indicates the truth yields for this point is 0 but the reco yields is not 0

Signal efficiency in SR4Jhx b-Tag bin3 region for squark production one-step variable-x simplified models. The -1 value indicates the truth yields for this point is 0 but the reco yields is not 0

Signal efficiency in SR4Jhx b-Veto bin1 region for squark production one-step variable-x simplified models. The -1 value indicates the truth yields for this point is 0 but the reco yields is not 0

Signal efficiency in SR4Jhx b-Veto bin2 region for squark production one-step variable-x simplified models. The -1 value indicates the truth yields for this point is 0 but the reco yields is not 0

Signal efficiency in SR4Jhx b-Veto bin3 region for squark production one-step variable-x simplified models. The -1 value indicates the truth yields for this point is 0 but the reco yields is not 0

Signal efficiency in SR4Jlx discovery region for squark production one-step variable-x simplified models. The -1 value indicates the truth yields for this point is 0 but the reco yields is not 0

Signal efficiency in SR4Jlx b-Tag bin1 region for squark production one-step variable-x simplified models. The -1 value indicates the truth yields for this point is 0 but the reco yields is not 0

Signal efficiency in SR4Jlx b-Tag bin2 region for squark production one-step variable-x simplified models. The -1 value indicates the truth yields for this point is 0 but the reco yields is not 0

Signal efficiency in SR4Jlx b-Tag bin3 region for squark production one-step variable-x simplified models. The -1 value indicates the truth yields for this point is 0 but the reco yields is not 0

Signal efficiency in SR4Jlx b-Veto bin1 region for squark production one-step variable-x simplified models. The -1 value indicates the truth yields for this point is 0 but the reco yields is not 0

Signal efficiency in SR4Jlx b-Veto bin2 region for squark production one-step variable-x simplified models. The -1 value indicates the truth yields for this point is 0 but the reco yields is not 0

Signal efficiency in SR4Jlx b-Veto bin3 region for squark production one-step variable-x simplified models. The -1 value indicates the truth yields for this point is 0 but the reco yields is not 0

Signal efficiency in SR6J b-Tag bin1 region for squark production one-step variable-x simplified models. The -1 value indicates the truth yields for this point is 0 but the reco yields is not 0

Signal efficiency in SR6J b-Tag bin2 region for squark production one-step variable-x simplified models. The -1 value indicates the truth yields for this point is 0 but the reco yields is not 0

Signal efficiency in SR6J b-Tag bin3 region for squark production one-step variable-x simplified models. The -1 value indicates the truth yields for this point is 0 but the reco yields is not 0

Signal efficiency in SR6J b-Tag bin4 region for squark production one-step variable-x simplified models. The -1 value indicates the truth yields for this point is 0 but the reco yields is not 0

Signal efficiency in SR6J b-Veto bin1 region for squark production one-step variable-x simplified models. The -1 value indicates the truth yields for this point is 0 but the reco yields is not 0

Signal efficiency in SR6J b-Veto bin2 region for squark production one-step variable-x simplified models. The -1 value indicates the truth yields for this point is 0 but the reco yields is not 0

Signal efficiency in SR6J b-Veto bin3 region for squark production one-step variable-x simplified models. The -1 value indicates the truth yields for this point is 0 but the reco yields is not 0

Signal efficiency in SR6J b-Veto bin4 region for squark production one-step variable-x simplified models. The -1 value indicates the truth yields for this point is 0 but the reco yields is not 0

Signal efficiency in SR6J discovery high region for squark production one-step variable-x simplified models. The -1 value indicates the truth yields for this point is 0 but the reco yields is not 0

Signal efficiency in SR6J discovery low region for squark production one-step variable-x simplified models. The -1 value indicates the truth yields for this point is 0 but the reco yields is not 0

Covariances $c_{ij}$ for the highest ranked systematic uncertainties in the energy asymmetry measurement differential in $\theta_j$.

Definition of the fiducial phase space with lepton candidate ($\ell$), hadronic top candidate ($j_h$), leptonic top candidate ($j_l$) and associated jet candidate ($j_a$).

Bounds on individual Wilson coefficients $C($TeV$/\Lambda)^2$ from the energy asymmetry, obtained from a combined fit to its values in all three $\theta_j$ bins.

Distribution of $S$, which is one of the discriminating variables used to define the $2\ell^{sc}$ SRs. The events are selected with the preselection requirements listed in Table 4 in the paper. The data (dots) are compared with the expected contributions from the relevant background sources (histograms). The expected signal distributions for $m_{H^{\pm\pm}} = 300~GeV$ are also shown, scaled to the observed number of events. The last bin includes overflows.

Distribution of $E_{T}^{miss}$, which is one of the discriminating variables used to define the $3\ell$ SRs. The events are selected with the preselection requirements listed in Table 4 in the paper. The data (dots) are compared with the expected contributions from the relevant background sources (histograms). The expected signal distributions for $m_{H^{\pm\pm}} = 300~GeV$ are also shown, scaled to the observed number of events. The last bin includes overflows.

Distribution of $\Delta R_{\ell^{\pm}\ell^{\pm}}$, which is one of the discriminating variables used to define the $3\ell$ SRs. The events are selected with the preselection requirements listed in Table 4 in the paper. The data (dots) are compared with the expected contributions from the relevant background sources (histograms). The expected signal distributions for $m_{H^{\pm\pm}} = 300~GeV$ are also shown, scaled to the observed number of events. The last bin includes overflows.

Distribution of $m_{x\ell}$ ($x$=3), which is one of the discriminating variables used to define the $3\ell$ SRs. The events are selected with the preselection requirements listed in Table 4 in the paper. The data (dots) are compared with the expected contributions from the relevant background sources (histograms). The expected signal distributions for $m_{H^{\pm\pm}} = 300~GeV$ are also shown, scaled to the observed number of events. The last bin includes overflows.

Distribution of $p_{T}^{leading jet}$, which is one of the discriminating variables used to define the $3\ell$ SRs. The events are selected with the preselection requirements listed in Table 4 in the paper. The data (dots) are compared with the expected contributions from the relevant background sources (histograms). The expected signal distributions for $m_{H^{\pm\pm}} = 300~GeV$ are also shown, scaled to the observed number of events. The last bin includes overflows.

Distribution of $E_{T}^{miss}$, which is one of the discriminating variables used to define the $4\ell$ SRs. The events are selected with the preselection requirements listed in Table 4 in the paper. The data (dots) are compared with the expected contributions from the relevant background sources (histograms). The expected signal distributions for $m_{H^{\pm\pm}} = 300~GeV$ are also shown, scaled to the observed number of events. The last bin includes overflows.

Distribution of $\Delta R_{\ell^{\pm}\ell^{\pm}}^{min}$, which is one of the discriminating variables used to define the $4\ell$ SRs. The events are selected with the preselection requirements listed in Table 4 in the paper. The data (dots) are compared with the expected contributions from the relevant background sources (histograms). The expected signal distributions for $m_{H^{\pm\pm}} = 300~GeV$ are also shown, scaled to the observed number of events. The last bin includes overflows.

Distribution of $m_{x\ell}$ ($x$=4), which is one of the discriminating variables used to define the $4\ell$ SRs. The events are selected with the preselection requirements listed in Table 4 in the paper. The data (dots) are compared with the expected contributions from the relevant background sources (histograms). The expected signal distributions for $m_{H^{\pm\pm}} = 300~GeV$ are also shown, scaled to the observed number of events. The last bin includes overflows.

Distribution of $p_{T}^{\ell_{1}}$, which is one of the discriminating variables used to define the $4\ell$ SRs. The events are selected with the preselection requirements listed in Table 4 in the paper. The data (dots) are compared with the expected contributions from the relevant background sources (histograms). The expected signal distributions for $m_{H^{\pm\pm}} = 300~GeV$ are also shown, scaled to the observed number of events. The last bin includes overflows.

The $E_{T}^{miss}$ distribution for the SRs of the $m_{H^{\pm\pm}} = 300~GeV$ signal mass hypothesis, where the selection requirement on $E_{T}^{miss}$ has been removed. In the attached plot, the signals are stacked on top of the backgrounds while individuals contributions of the $2\ell^{sc}$ channel are shown in HEPData. The last bin, isolated by a vertical red dashed line, is inclusive and corresponds to the SR.

The $E_{T}^{miss}$ distribution for the SRs of the $m_{H^{\pm\pm}} = 300~GeV$ signal mass hypothesis, where the selection requirement on $E_{T}^{miss}$ has been removed. In the attached plot, the signals are stacked on top of the backgrounds while individuals contributions of the $3\ell$ channel are shown in HEPData. The last bin, isolated by a vertical red dashed line, is inclusive and corresponds to the SR.

The $E_{T}^{miss}$ distribution for the SRs of the $m_{H^{\pm\pm}} = 300~GeV$ signal mass hypothesis, where the selection requirement on $E_{T}^{miss}$ has been removed. In the attached plot, the signals are stacked on top of the backgrounds while individuals contributions of the $4\ell$ channel are shown in HEPData. The last bin, isolated by a vertical red dashed line, is inclusive and corresponds to the SR.