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.

Post-fit event yields in the signal and control regions obtained from a scan over $\epsilon_{VV}$ exploiting both shape and rate information. The quoted uncertainties include the theoretical and experimental systematic sources and those due to sample statistics. The fit constrains the total expected yield to the observed yield. The diboson background is split into $W W$ and non-$W W$ contributions.

Best-fit values and their uncertainties as obtained from the shape-only and shape-plus-rate likelihood fits to the Asimov dataset and to ATLAS data. Results of both shape-only and shape+rate fits for $a_L$ and $a_T$ are shown. Results of fits to one parameter with the other one fixed or profiled are presented.

Best-fit values and their uncertainties as obtained from the shape-only and shape-plus-rate likelihood fits to the Asimov dataset and to ATLAS data. Results of both shape-only and shape+rate fits for $\epsilon_{VV}$ and $\kappa_{VV}$ are shown. Results of fits to one parameter with the other one fixed or profiled are presented.

The contributions of the leading individual systematic uncertainties together with the data statistical uncertainties, in the one dimensional fit for the pseudo-observables $\kappa_{VV}$ (a) and $\epsilon_{VV}$ (b) for electroweak-boson polarisation in the VBF $H\to WW$ channel. Both shape and rate informations are exploited in the fit. The theoretical and experimental uncertainties are subdivided further into categories.

The contributions of the leading individual systematic uncertainties together with the data statistical uncertainties, in the one dimensional fit for the pseudo-observables $\kappa_{VV}$ (a) and $\epsilon_{VV}$ (b) for electroweak-boson polarisation in the VBF $H\to WW$ channel. Both shape and rate informations are exploited in the fit. The theoretical and experimental uncertainties are subdivided further into categories.

Post-fit distribution of the BDT response observable presented in the four $|\Delta \eta jj|$ categories of the ggF +2 jets signal region, with signal and background yields fixed from the fit for $\mu^{\text{ggF+2jets}}$. The distributions of the ggF + 2 jets and VBF processes are overlaid with their respective contributions multiplied by 50.

The weighted $\Delta \Phi_{jj}$ post-fit distribution in the ggF +2 jets signal region, with signal and background yields fixed from the fit to $\tan \alpha$ using shape and rate information.

Expected and observed likelihood curves for scans over $\tan \alpha$ where only the shape is taken into account in the fit, $\mu_{VBF}$ is fixed.

Expected and observed likelihood curves for scans over $\tan \alpha$ where both shape and normalisation are taken into account in the fit, $\mu_{VBF}$ is fixed.

68% and 95% CL two-dimensional likelihood contours of the CP-even and CP-odd coupling parameters $K_{gg} \cos(\alpha)$ and $K_{gg} \sin(\alpha)$. The minima are represented by black stars, while the SM value is shown as a red star.

The weighted $\Delta \Phi jj$ distribution in the VBF signal region, with signal and background yields fixed from the fit for $\epsilon_{VV}$ using shape and rate information.

Likelihood scans over the transversally polarised couplings. The fit is using shape-only information. All relevant experimental and modelling systematic uncertainties are considered in the fit.

Likelihood scans over the transversally polarised couplings. The fit is using shape + rate information. All relevant experimental and modelling systematic uncertainties are considered in the fit.

Likelihood scans over the longitudinally polarised couplings. The fit is using shape + rate information. All relevant experimental and modelling systematic uncertainties are considered in the fit.

Likelihood scans over $\kappa_{VV}$ with the $\epsilon_{VV}$ profiled. The fit is performed using both shape and rate information. All relevant experimental and theoretical systematic uncertainties are considered in the fit.

Likelihood scans over $\epsilon_{VV}$ with the $\kappa_{VV}$ profiled. The fit is performed using both shape and rate information. All relevant experimental and theoretical systematic uncertainties are considered in the fit.

The contributions of the leading individual systematic uncertainties together with the data statistical uncertainties, in the one dimensional fit for electroweak-boson polarisation in the VBF $H\to WW$ channel, using (aL, aT) parametrisation. Both shape and rate informations are exploited in the fit. The theoretical and experimental uncertainties are subdivided further into categories.

The contributions of the leading individual systematic uncertainties together with the data statistical uncertainties, in the one dimensional fit for electroweak-boson polarisation in the VBF $H\to WW$ channel, using (aL, aT) parametrisation.. Both shape and rate informations are exploited in the fit. The theoretical and experimental uncertainties are subdivided further into categories.

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 +- 1 sigma 95% CL exclusion limit for the Zprime-2HDM model.

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.

Transverse profile for 70 GeV < pT < 100 GeV.

Longitudinal profile for pT > 100 GeV.

Transverse profile for pT > 100 GeV.

Average longitudinal profile

Transverse profile for pT > 100 GeV.

Differential cross section as a function of $m_{\gamma\gamma}$. The table contains the values measured in data and theory predictions from SHERPA, DIPHOX and NNLOJET.

Differential cross section as a function of $p_{T,\gamma\gamma}$. The table contains the values measured in data and theory predictions from SHERPA, DIPHOX and NNLOJET.

Differential cross section as a function of $a_{T,\gamma\gamma}$. The table contains the values measured in data and theory predictions from SHERPA, DIPHOX and NNLOJET.

Differential cross section as a function of $\phi_{\eta}*$. The table contains the values measured in data and theory predictions from SHERPA, DIPHOX and NNLOJET.

Differential cross section as a function of $\pi-\Delta\phi_{\gamma\gamma}$. The table contains the values measured in data and theory predictions from SHERPA, DIPHOX and NNLOJET.

Differential cross section as a function of $|cos\theta*|^{(CS)}$. The table contains the values measured in data and theory predictions from SHERPA, DIPHOX and NNLOJET.

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}$.

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.

The observed data event yields and the corresponding estimates for the backgrounds in the different $b$-jet multiplicity bins for the 20 GeV jet $p_{\mathrm{T}}$ threshold regions defined for the EWK analysis in the $1\ell$ category for 7 jets. The background is estimated by including all bins in the fit. All uncertainties, which may be correlated across the bins, are included in the total background uncertainty.

The observed data event yields and the corresponding estimates for the backgrounds in the different $b$-jet multiplicity bins for the 20 GeV jet $p_{\mathrm{T}}$ threshold regions defined for the EWK analysis in the $1\ell$ category for 8 jets. The background is estimated by including all bins in the fit. All uncertainties, which may be correlated across the bins, are included in the total background uncertainty.

The observed data event yields and the corresponding estimates for the backgrounds in the different $b$-jet multiplicity bins for the 20 GeV jet $p_{\mathrm{T}}$ threshold regions defined for the EWK analysis in the $1\ell$ category for 9 jets. The background is estimated by including all bins in the fit. All uncertainties, which may be correlated across the bins, are included in the total background uncertainty.

The observed data event yields and the corresponding estimates for the backgrounds in the different $b$-jet multiplicity bins for the 20 GeV jet $p_{\mathrm{T}}$ threshold regions defined for the EWK analysis in the $1\ell$ category for 10 jets. The background is estimated by including all bins in the fit. All uncertainties, which may be correlated across the bins, are included in the total background uncertainty.

The observed data event yields and the corresponding estimates for the backgrounds in the different $b$-jet multiplicity bins for the 20 GeV jet $p_{\mathrm{T}}$ threshold regions defined for the EWK analysis in the $1\ell$ category for 11 jets. The background is estimated by including all bins in the fit. All uncertainties, which may be correlated across the bins, are included in the total background uncertainty.

The observed data event yields and the corresponding estimates for the backgrounds in the different $b$-jet multiplicity bins for the 20 GeV jet $p_{\mathrm{T}}$ threshold regions defined for the EWK analysis in the $1\ell$ category for 12 jets. The background is estimated by including all bins in the fit. All uncertainties, which may be correlated across the bins, are included in the total background uncertainty.

The observed data event yields and the corresponding estimates for the backgrounds in the different $b$-jet multiplicity bins for the 20 GeV jet $p_{\mathrm{T}}$ threshold regions defined for the EWK analysis in the $1\ell$ category for 13 jets. The background is estimated by including all bins in the fit. All uncertainties, which may be correlated across the bins, are included in the total background uncertainty.

The observed data event yields and the corresponding estimates for the backgrounds in the different $b$-jet multiplicity bins for the 20 GeV jet $p_{\mathrm{T}}$ threshold regions defined for the EWK analysis in the $1\ell$ category for 14 jets. The background is estimated by including all bins in the fit. All uncertainties, which may be correlated across the bins, are included in the total background uncertainty.

The observed data event yields and the corresponding estimates for the backgrounds in the different $b$-jet multiplicity bins for the 20 GeV jet $p_{\mathrm{T}}$ threshold regions defined for the EWK analysis in the $1\ell$ category for at least 15 jets. The background is estimated by including all bins in the fit. All uncertainties, which may be correlated across the bins, are included in the total background uncertainty.

The observed data event yields and the corresponding estimates for the backgrounds in the different $b$-jet multiplicity bins for the 20 GeV jet $p_{\mathrm{T}}$ threshold regions defined for the EWK analysis in the $2\ell^{\mathrm{sc}}$ category for 4 jets. The background is estimated by including all bins in the fit. All uncertainties, which may be correlated across the bins, are included in the total background uncertainty.

The observed data event yields and the corresponding estimates for the backgrounds in the different $b$-jet multiplicity bins for the 20 GeV jet $p_{\mathrm{T}}$ threshold regions defined for the EWK analysis in the $2\ell^{\mathrm{sc}}$ category for 5 jets. The background is estimated by including all bins in the fit. All uncertainties, which may be correlated across the bins, are included in the total background uncertainty.

The observed data event yields and the corresponding estimates for the backgrounds in the different $b$-jet multiplicity bins for the 20 GeV jet $p_{\mathrm{T}}$ threshold regions defined for the EWK analysis in the $2\ell^{\mathrm{sc}}$ category for 6 jets. The background is estimated by including all bins in the fit. All uncertainties, which may be correlated across the bins, are included in the total background uncertainty.

The observed data event yields and the corresponding estimates for the backgrounds in the different $b$-jet multiplicity bins for the 20 GeV jet $p_{\mathrm{T}}$ threshold regions defined for the EWK analysis in the $2\ell^{\mathrm{sc}}$ category for 7 jets. The background is estimated by including all bins in the fit. All uncertainties, which may be correlated across the bins, are included in the total background uncertainty.

The observed data event yields and the corresponding estimates for the backgrounds in the different $b$-jet multiplicity bins for the 20 GeV jet $p_{\mathrm{T}}$ threshold regions defined for the EWK analysis in the $2\ell^{\mathrm{sc}}$ category for 8 jets. The background is estimated by including all bins in the fit. All uncertainties, which may be correlated across the bins, are included in the total background uncertainty.

The observed data event yields and the corresponding estimates for the backgrounds in the different $b$-jet multiplicity bins for the 20 GeV jet $p_{\mathrm{T}}$ threshold regions defined for the EWK analysis in the $2\ell^{\mathrm{sc}}$ category for 9 jets. The background is estimated by including all bins in the fit. All uncertainties, which may be correlated across the bins, are included in the total background uncertainty.

The observed data event yields and the corresponding estimates for the backgrounds in the different $b$-jet multiplicity bins for the 20 GeV jet $p_{\mathrm{T}}$ threshold regions defined for the EWK analysis in the $2\ell^{\mathrm{sc}}$ category for at least 10 jets. The background is estimated by including all bins in the fit. All uncertainties, which may be correlated across the bins, are included in the total background uncertainty.

The observed data event yields and the corresponding estimates for the backgrounds in the different $b$-jet multiplicity bins for the 20 GeV jet $p_{\mathrm{T}}$ threshold region in the $1\ell$ category for at least 15 jets. The background is estimated by including all bins in the fit. All uncertainties, which may be correlated across the bins, are included in the total background uncertainty.

The observed data event yields and the corresponding estimates for the backgrounds in the different $b$-jet multiplicity bins for the 40 GeV jet $p_{\mathrm{T}}$ threshold region in the $1\ell$ category for at least 12 jets. The background is estimated by including all bins in the fit. All uncertainties, which may be correlated across the bins, are included in the total background uncertainty.

The observed data event yields and the corresponding estimates for the backgrounds in the different $b$-jet multiplicity bins for the 60 GeV jet $p_{\mathrm{T}}$ threshold region in the $1\ell$ category for at least 11 jets. The background is estimated by including all bins in the fit. All uncertainties, which may be correlated across the bins, are included in the total background uncertainty.

The observed data event yields and the corresponding estimates for the backgrounds in the different $b$-jet multiplicity bins for the 80 GeV jet $p_{\mathrm{T}}$ threshold region in the $1\ell$ category for at least 10 jets. The background is estimated by including all bins in the fit. All uncertainties, which may be correlated across the bins, are included in the total background uncertainty.

The observed data event yields and the corresponding estimates for the backgrounds in the different $b$-jet multiplicity bins for the 100 GeV jet $p_{\mathrm{T}}$ threshold region in the $1\ell$ category for at least 8 jets. The background is estimated by including all bins in the fit. All uncertainties, which may be correlated across the bins, are included in the total background uncertainty.

The observed data event yields and the corresponding estimates for the backgrounds in the different $b$-jet multiplicity bins for the 20 GeV jet $p_{\mathrm{T}}$ threshold region in the $2\ell^{\mathrm{sc}}$ category for at least 10 jets. The background is estimated by including all bins in the fit. All uncertainties, which may be correlated across the bins, are included in the total background uncertainty.

The observed data event yields and the corresponding estimates for the backgrounds in the different $b$-jet multiplicity bins for the 40 GeV jet $p_{\mathrm{T}}$ threshold region in the $2\ell^{\mathrm{sc}}$ category for at least 8 jets. The background is estimated by including all bins in the fit. All uncertainties, which may be correlated across the bins, are included in the total background uncertainty.

The observed data event yields and the corresponding estimates for the backgrounds in the different $b$-jet multiplicity bins for the 60 GeV jet $p_{\mathrm{T}}$ threshold region in the $2\ell^{\mathrm{sc}}$ category for at least 7 jets. The background is estimated by including all bins in the fit. All uncertainties, which may be correlated across the bins, are included in the total background uncertainty.

The observed data event yields and the corresponding estimates for the backgrounds in the different $b$-jet multiplicity bins for the 80 GeV jet $p_{\mathrm{T}}$ threshold region in the $2\ell^{\mathrm{sc}}$ category for at least 7 jets. The background is estimated by including all bins in the fit. All uncertainties, which may be correlated across the bins, are included in the total background uncertainty.

The observed data event yields and the corresponding estimates for the backgrounds in the different $b$-jet multiplicity bins for the 100 GeV jet $p_{\mathrm{T}}$ threshold region in the $2\ell^{\mathrm{sc}}$ category for at least 6 jets. The background is estimated by including all bins in the fit. All uncertainties, which may be correlated across the bins, are included in the total background uncertainty.

Data event yields compared with the expected contributions from relevant background sources, in the discovery signal regions defined for the $1\ell$ category, as well as the observed and expected 95% CL model-independent upper limits on product of cross-section, acceptance and efficiency (in fb). The parameters of the model are determined in a fit to a reduced set of bins, corresponding to the model-independent fit discussed in the text.

Data event yields compared with the expected contributions from relevant background sources, in the discovery signal regions defined for the $2\ell^{\mathrm{sc}}$ category, as well as the observed and expected 95% CL model-independent upper limits on product of cross-section, acceptance and efficiency (in fb). The parameters of the model are determined in a fit to a reduced set of bins, corresponding to the model-independent fit discussed in the text.

Expected 95% CL exclusion contour for the gluino $\tilde{g}\rightarrow q\bar{q}\tilde{\chi}^0_1\rightarrow q\bar{q}q\bar{q} \ell\nu$ model.

Observed 95% CL exclusion contour for the gluino $\tilde{g}\rightarrow q\bar{q}\tilde{\chi}^0_1\rightarrow q\bar{q}q\bar{q} \ell\nu$ model.

Expected 95% CL exclusion contour for the gluino $\tilde{g}\rightarrow \bar{t}\tilde{t} \rightarrow \bar{t}\bar{b}\bar{s}$ model.

Observed 95% CL exclusion contour for the gluino $\tilde{g}\rightarrow \bar{t}\tilde{t} \rightarrow \bar{t}\bar{b}\bar{s}$ model.

Expected 95% CL exclusion contour for the gluino $\tilde{g}\rightarrow t\bar{t}\tilde{\chi}^0_1 \rightarrow t\bar{t} tbs$ model with bino LSP.

Observed 95% CL exclusion contour for the gluino $\tilde{g}\rightarrow t\bar{t}\tilde{\chi}^0_1 \rightarrow t\bar{t} tbs$ model with bino LSP.

Expected 95% CL exclusion contour for the gluino $\tilde{g}\rightarrow t\bar{t}\tilde{\chi}^0_1 \rightarrow t\bar{t} tbs$ model with higgsino LSP.

Observed 95% CL exclusion contour for the gluino $\tilde{g}\rightarrow t\bar{t}\tilde{\chi}^0_1 \rightarrow t\bar{t} tbs$ model with higgsino LSP.

Expected 95% CL exclusion contour for the gluino $\tilde{g}\rightarrow t\bar{t}\tilde{\chi}^0_1 \rightarrow t\bar{t} tbs$ model with wino LSP.

Observed 95% CL exclusion contour for the gluino $\tilde{g}\rightarrow t\bar{t}\tilde{\chi}^0_1 \rightarrow t\bar{t} tbs$ model with wino LSP.

Expected 95% CL exclusion contour for the stop pair production model with bino LSP.

Observed 95% CL exclusion contour for the stop pair production model with bino LSP.

Expected 95% CL exclusion contour for the stop pair production model with higgsino LSP.

Observed 95% CL exclusion contour for the stop pair production model with higgsino LSP.

Expected 95% CL exclusion contour for the stop pair production model with wino LSP.

Observed 95% CL exclusion contour for the stop pair production model with wino LSP.

Expected 95% CL excluded cross section for the RPV model with electroweakino prodiction with higgsino LSP hypothesis.

Observed 95% CL excluded cross section for the RPV model with electroweakino prodiction with higgsino LSP hypothesis.

Expected 95% CL excluded cross section for the RPV model with electroweakino prodiction with wino LSP hypothesis.

Observed 95% CL excluded cross section for the RPV model with electroweakino prodiction with wino LSP hypothesis.

Data event yields compared with the expected contributions from relevant background sources, in the discovery signal regions defined for the $1\ell$ category with 20 and 40 GeV jet $p_{\mathrm{T}}$ thresholds. The uncertainties across backgrounds can exhibit strong anticorrelations.

Data event yields compared with the expected contributions from relevant background sources, in the discovery signal regions defined for the $1\ell$ category with 60, 80 and 100 GeV jet $p_{\mathrm{T}}$ thresholds. The uncertainties across backgrounds can exhibit strong anticorrelations.

Data event yields compared with the expected contributions from relevant background sources, in the discovery signal regions defined for the $2\ell^{\mathrm{sc}}$ category with 20 and 40 GeV jet $p_{\mathrm{T}}$ thresholds. The uncertainties across backgrounds can exhibit strong anticorrelations.

Data event yields compared with the expected contributions from relevant background sources, in the discovery signal regions defined for the $2\ell^{\mathrm{sc}}$ category with 60, 80 and 100 GeV jet $p_{\mathrm{T}}$ thresholds. The uncertainties across backgrounds can exhibit strong anticorrelations.

Expected 95% CL excluded cross section for the gluino $\tilde{g}\rightarrow q\bar{q}\tilde{\chi}^0_1\rightarrow q\bar{q}q\bar{q} \ell\nu$ model.

Observed 95% CL excluded cross section for the gluino $\tilde{g}\rightarrow q\bar{q}\tilde{\chi}^0_1\rightarrow q\bar{q}q\bar{q} \ell\nu$ model.

Expected 95% CL excluded cross section for the gluino $\tilde{g}\rightarrow \bar{t}\tilde{t} \rightarrow \bar{t}\bar{b}\bar{s}$ model.

Observed 95% CL excluded cross section for the gluino $\tilde{g}\rightarrow \bar{t}\tilde{t} \rightarrow \bar{t}\bar{b}\bar{s}$ model.

Expected 95% CL excluded cross section for the gluino $\tilde{g}\rightarrow t\bar{t}\tilde{\chi}^0_1 \rightarrow t\bar{t} tbs$ model with bino LSP.

Observed 95% CL excluded cross section for the gluino $\tilde{g}\rightarrow t\bar{t}\tilde{\chi}^0_1 \rightarrow t\bar{t} tbs$ model with bino LSP.

Expected 95% CL excluded cross section for the gluino $\tilde{g}\rightarrow t\bar{t}\tilde{\chi}^0_1 \rightarrow t\bar{t} tbs$ model with wino LSP.

Observed 95% CL excluded cross section for the gluino $\tilde{g}\rightarrow t\bar{t}\tilde{\chi}^0_1 \rightarrow t\bar{t} tbs$ model with wino LSP.

Expected 95% CL excluded cross section for the gluino $\tilde{g}\rightarrow t\bar{t}\tilde{\chi}^0_1 \rightarrow t\bar{t} tbs$ model with higgsino LSP.

Observed 95% CL excluded cross section for the gluino $\tilde{g}\rightarrow t\bar{t}\tilde{\chi}^0_1 \rightarrow t\bar{t} tbs$ model with higgsino LSP.

Expected 95% CL excluded cross section for the stop pair production model with bino LSP.

Observed 95% CL excluded cross section for the stop pair production model with bino LSP.

Expected 95% CL excluded cross section for the stop pair production model with wino LSP.

Observed 95% CL excluded cross section for the stop pair production model with wino LSP.

Expected 95% CL excluded cross section for the stop pair production model with higgsino LSP.

Observed 95% CL excluded cross section for the stop pair production model with higgsino LSP.

Acceptance and efficiency for the gluino $\tilde{g}\rightarrow q\bar{q}\tilde{\chi}^0_1\rightarrow q\bar{q}q\bar{q} \ell\nu$ model in the $1\ell$ category.

Acceptance and efficiency for the gluino $\tilde{g}\rightarrow q\bar{q}\tilde{\chi}^0_1\rightarrow q\bar{q}q\bar{q} \ell\nu$ model in the $2\ell^{\mathrm{sc}}$ category.

Acceptance and efficiency for the gluino $\tilde{g}\rightarrow \bar{t}\tilde{t} \rightarrow \bar{t}\bar{b}\bar{s}$ model in the $1\ell$ category.

Acceptance and efficiency for the gluino $\tilde{g}\rightarrow \bar{t}\tilde{t} \rightarrow \bar{t}\bar{b}\bar{s}$ model in the $2\ell^{\mathrm{sc}}$ category.

Acceptance and efficiency for the gluino $\tilde{g}\rightarrow t\bar{t}\tilde{\chi}^0_1 \rightarrow t\bar{t} tbs$ model with higgsino LSP in the $1\ell$ category.

Acceptance and efficiency for the gluino $\tilde{g}\rightarrow t\bar{t}\tilde{\chi}^0_1 \rightarrow t\bar{t} tbs$ model with higgsino LSP in the $2\ell^{\mathrm{sc}}$ category.

Acceptance and efficiency for the stop model with higgsino LSP in the $1\ell$ category.

Acceptance and efficiency for the stop model with higgsino LSP in the $2\ell^{\mathrm{sc}}$ category.

Acceptance and efficiency for the or the electroweakino production model in the EWK analysis discovery SR for the $1\ell$ category with 4 jets, considering $\tilde{\chi}^{\pm}_1 \tilde{\chi}^0_1$ production.

Acceptance and efficiency for the or the electroweakino production model in the EWK analysis discovery SR for the $1\ell$ category with 5 jets, considering $\tilde{\chi}^{\pm}_1 \tilde{\chi}^0_1$ production.

Acceptance and efficiency for the or the electroweakino production model in the EWK analysis discovery SR for the $1\ell$ category with 6 jets, considering $\tilde{\chi}^{\pm}_1 \tilde{\chi}^0_1$ production.

Acceptance and efficiency for the or the electroweakino production model in the EWK analysis discovery SR for the $1\ell$ category with 7 jets, considering $\tilde{\chi}^{\pm}_1 \tilde{\chi}^0_1$ production.

Acceptance and efficiency for the or the electroweakino production model in the EWK analysis discovery SR for the $1\ell$ category with 8 jets, considering $\tilde{\chi}^{\pm}_1 \tilde{\chi}^0_1$ production.

Acceptance and efficiency for the or the electroweakino production model in the EWK analysis discovery SR for the $2\ell^{\mathrm{sc}}$ category, considering $\tilde{\chi}^{\pm}_1 \tilde{\chi}^0_1$ production.

Acceptance and efficiency for the or the electroweakino production model in the EWK analysis discovery SR for the $1\ell$ category with 4 jets, considering $\tilde{\chi}^0_1 \tilde{\chi}^0_2$ production.

Acceptance and efficiency for the or the electroweakino production model in the EWK analysis discovery SR for the $1\ell$ category with 5 jets, considering $\tilde{\chi}^0_1 \tilde{\chi}^0_2$ production.

Acceptance and efficiency for the or the electroweakino production model in the EWK analysis discovery SR for the $1\ell$ category with 6 jets, considering $\tilde{\chi}^0_1 \tilde{\chi}^0_2$ production.

Acceptance and efficiency for the or the electroweakino production model in the EWK analysis discovery SR for the $1\ell$ category with 7 jets, considering $\tilde{\chi}^0_1 \tilde{\chi}^0_2$ production.

Acceptance and efficiency for the or the electroweakino production model in the EWK analysis discovery SR for the $1\ell$ category with 8 jets, considering $\tilde{\chi}^0_1 \tilde{\chi}^0_2$ production.

Acceptance and efficiency for the or the electroweakino production model in the EWK analysis discovery SR for the $2\ell^{\mathrm{sc}}$ category, considering $\tilde{\chi}^0_1 \tilde{\chi}^0_2$ production.

Cut flow for the gluino $\tilde{g}\rightarrow q\bar{q}\tilde{\chi}^0_1\rightarrow q\bar{q}q\bar{q} \ell\nu$ model with $m_{\tilde{g}} = 2$ TeV and $m_{\tilde{\chi}^0_1} = 1$ TeV. The column labelled $\mathrm{N}_{\mathrm{raw}}$ shows the number of generated events, while $\mathrm{N}_{\mathrm{events}}$ shows the expected number of events with a luminosity of 139fb$^{−1}$. The last column shows the cut flow efficiency with respect to all weighted events. The events are skimmed by requiring at least one electron or muon that satisfies very loose identification criteria, where the lepton satisfies $p_{\mathrm{T}} > 25$ GeV. The efficiency of this cut is considered in the quoted efficiency of the lepton trigger requirement. Selections with negligible inefficiencies on the given sample, such as data quality requirements, are not displayed. Selections that have not been evaluated in the analysis or are not applicable are denoted with a dash (--).

Cut flow for the gluino $\tilde{g}\rightarrow \bar{t}\tilde{t} \rightarrow \bar{t}\bar{b}\bar{s}$ model with with $m_{\tilde{g}} = 1.6$ TeV and $m_{\tilde{t}} = 1$ TeV. The column labelled $\mathrm{N}_{\mathrm{raw}}$ shows the number of generated events, while $\mathrm{N}_{\mathrm{events}}$ shows the expected number of events with a luminosity of 139fb$^{−1}$. The last column shows the cut flow efficiency with respect to all weighted events. The events are skimmed by requiring at least one electron or muon that satisfies very loose identification criteria, where the lepton satisfies $p_{\mathrm{T}} > 25$ GeV. The efficiency of this cut is considered in the quoted efficiency of the lepton trigger requirement. Selections with negligible inefficiencies on the given sample, such as data quality requirements, are not displayed. Selections that have not been evaluated in the analysis or are not applicable are denoted with a dash (--).

Cut flow for the gluino $\tilde{g}\rightarrow t\bar{t}\tilde{\chi}^0_1 \rightarrow t\bar{t} tbs$ model with $m_{\tilde{g}} = 2.2$ TeV and $m_{\tilde{\chi}^0_1} = 1.05$ TeV. The column labelled $\mathrm{N}_{\mathrm{raw}}$ shows the number of generated events, while $\mathrm{N}_{\mathrm{events}}$ shows the expected number of events with a luminosity of 139fb$^{−1}$. The last column shows the cut flow efficiency with respect to all weighted events. The events are skimmed by requiring at least one electron or muon that satisfies very loose identification criteria, where the lepton satisfies $p_{\mathrm{T}} > 25$ GeV. The efficiency of this cut is considered in the quoted efficiency of the lepton trigger requirement. Selections with negligible inefficiencies on the given sample, such as data quality requirements, are not displayed. Selections that have not been evaluated in the analysis or are not applicable are denoted with a dash (--).

Cut flow for the stop model with $m_{\tilde{t}} = 1.175$ TeV and $m_{\tilde{\chi}^0_1} = 0.7$ TeV. The column labelled $\mathrm{N}_{\mathrm{raw}}$ shows the number of generated events, while $\mathrm{N}_{\mathrm{events}}$ shows the expected number of events with a luminosity of 139fb$^{−1}$. The last column shows the cut flow efficiency with respect to all weighted events. The events are skimmed by requiring at least one electron or muon that satisfies very loose identification criteria, where the lepton satisfies $p_{\mathrm{T}} > 25$ GeV. The efficiency of this cut is considered in the quoted efficiency of the lepton trigger requirement. Selections with negligible inefficiencies on the given sample, such as data quality requirements, are not displayed. Selections that have not been evaluated in the analysis or are not applicable are denoted with a dash (--).

Cut flow for the electroweakino production model, considering only the production of $\tilde{\chi}^{\pm}_1 \tilde{\chi}^0_1$, with $m(\tilde{\chi}^{\pm}_1,\tilde{\chi}^0_1)= 250$ GeV. The column labelled $\mathrm{N}_{\mathrm{raw}}$ shows the number of generated events, while $\mathrm{N}_{\mathrm{events}}$ shows the expected number of events with a luminosity of 139fb$^{−1}$. The last column shows the cut flow efficiency with respect to all weighted events. The events are skimmed by requiring at least one electron or muon that satisfies very loose identification criteria, where the lepton satisfies $p_{\mathrm{T}} > 25$ GeV. The efficiency of this cut is considered in the quoted efficiency of the lepton trigger requirement. Selections with negligible inefficiencies on the given sample, such as data quality requirements, are not displayed. In the $2\ell^{\mathrm{sc}}$ category no events are expected, as only one lepton is expected to be produced in the decay.Selections that have not been evaluated in the analysis or are not applicable are denoted with a dash (--).

Cut flow for the electroweakino production model, considering only the production of $\tilde{\chi}^0_1 \tilde{\chi}^0_2$, with $m(\tilde{\chi}^0_1,\tilde{\chi}^0_2)= 250$ GeV. The column labelled $\mathrm{N}_{\mathrm{raw}}$ shows the number of generated events, while $\mathrm{N}_{\mathrm{events}}$ shows the expected number of events with a luminosity of 139fb$^{−1}$. The last column shows the cut flow efficiency with respect to all weighted events. The events are skimmed by requiring at least one electron or muon that satisfies very loose identification criteria, where the lepton satisfies $p_{\mathrm{T}} > 25$ GeV. The efficiency of this cut is considered in the quoted efficiency of the lepton trigger requirement. Selections with negligible inefficiencies on the given sample, such as data quality requirements, are not displayed. Selections that have not been evaluated in the analysis or are not applicable are denoted with a dash (--).