Theory cross sections vs m($\widetilde{\chi}^0_1$) for SMS model TChiHH-G, $\widetilde{\chi}^0_1, \widetilde{\chi}^0_2$ NLSPs only.

The 95% CL observed and expected upper limits on the production cross sections of the TChiHH signal model as a function of the $\widetilde{\chi}^0_2$ and $\widetilde{\chi}^0_1$ masses.

Exclusion limits assuming the approximate-NNLO+NNLL cross sections

Exclusion limits assuming the approximate-NNLO+NNLL cross sections

Exclusion limits assuming the approximate-NNLO+NNLL cross sections

Exclusion limits assuming the approximate-NNLO+NNLL cross sections

Exclusion limits assuming the approximate-NNLO+NNLL cross sections

Exclusion limits assuming the approximate-NNLO+NNLL cross sections

Cross section, 95% CL upper limit plus theory vs m($\widetilde{g}$) for SMS model T5HH.

Pre-fit background covariance matrix $\sigma_{xy}$ for the 22 analysis bins, ordered as in Fig. 10.

Pre-fit background correlation matrix $\rho_{xy}$ for the 22 analysis bins, ordered as in Fig. 10.

Significance vs $(m(\widetilde{\chi}^0_2), m(\widetilde{\chi}^0_1))$ for SMS model TChiHH.

Significance vs $m(\widetilde{g})$ for SMS model T5HH.

Efficiency vs $m(\widetilde{\chi}^0_1)$ for SMS model TChiHH-G, $\mathcal{B}(\mathrm{H}$-->$\mathrm{b}\overline{\mathrm{b}})=100\%$.

Efficiency vs $m(\widetilde{g})$ for SMS model T5HH.

Efficiency vs $(m(\widetilde{\chi}^0_2), m(\widetilde{\chi}^0_1))$ for SMS model TChiHH.

Signal and control region populations for the resolved signature.

Signal and control region populations for the boosted signature.

Cut flow table for the resolved topology showing expected event yields at 137 fb$^{-1}$ for selected signal models as successive selections are applied.

Cut flow table for the boosted topology showing expected event yields at 137 fb$^{-1}$ for selected signal model points as successive selections are applied. Here the hadronic baseline is defined by $N_{\mathrm{jet}}\geq$2, \ptmiss$>$300 GeV, $\Delta\phi$ cuts, and the isolated lepton and track vetoes.

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 $m_{\mathrm{jjj}}$ for preselected events with $\mathrm{N}_{j}$ = 3

Distribution of $m_{\mathrm{j}}$ for preselected events with $\mathrm{N}_{j}$ = 3

Distribution of the deep-WH value of the highest-mass jet with 60 < $m_{\mathrm{j}}$ < 100 GeV for preselected events with $\mathrm{N}_{j}$ = 3

scale factors (SFs) for W, $t^{2}$, and q/g matched jets in the low-$m_{\mathrm{j}}$ and low-$p_{\mathrm{T}}$ (LL) bin, as functions of the deep-W discriminant value.

scale factors (SFs) for W, $t^{2}$, and q/g matched jets in the low-$m_{\mathrm{j}}$ and high-$p_{\mathrm{T}}$ (LH) bin, as functions of the deep-W discriminant value.

scale factors (SFs) for $t^{2}$, $t^{3,4}$, and q/g matched jets in the high-$m_{\mathrm{j}}$ and low-$p_{\mathrm{T}}$ (HL) bin, as functions of the deep-WH discriminant value.

scale factors (SFs) for $t^{2}$, $t^{3,4}$, and q/g matched jets in the high-$m_{\mathrm{j}}$ and high-$p_{\mathrm{T}}$ (HH) bin, as functions of the deep-WH discriminant value.

The deep-W discriminant of the jet with highest mass in the single-lepton sideband for LL samples.

The deep-W discriminant of the jet with highest mass in the single-lepton sideband for LH samples.

The deep-WH discriminant of the jet with highest mass in the single-lepton sideband for HL samples.

The deep-WH discriminant of the jet with highest mass in the single-lepton sideband for HH samples.

Comparison of the distribution of data and simulated backgrounds, as a function of the deep-W discriminant value for the highest-mass jet in CR1, after the scale factors have been applied.

Comparison of the distribution of data and simulated backgrounds, as a function of the deep-WH discriminant value for the highest-mass jet in CR2, after the scale factors have been applied.

Comparison of the distribution of data and simulated backgrounds, as a function of the deep-WH discriminant value for the highest-mass jet in CR3, after the scale factors have been applied.

Comparison of the distribution of data and simulated backgrounds, as a function of the deep-W discriminant value for the highest-mass jet in CR45, after the scale factors have been applied.

Comparison of the distribution of data and simulated backgrounds, as a function of the deep-W discriminant value for the highest-mass jet in CR6, after the scale factors have been applied.

The $m_{\mathrm{j}^{\mathrm{max}}}$ distributions for different jet types for SR1–3 events of the signal with $m_{\mathrm{W}_{\mathrm{KK}}}$ = 2.5 TeV, $m_{\mathrm{R}}$ = 0.2 TeV$ without deep-W (WH) constraints.

The deep-W distribution normalized to unity for the shown components of of the signal with $m_{\mathrm{W}_{\mathrm{KK}}}$ = 2.5 TeV, $m_{\mathrm{R}}$ = 0.2 TeV. The $t^{3,4}$ jets from the preselected sample, normalized to unity, are superimposed to compare shapes with the $R^{3q}$ and $R^{4q}$ distributions.

The deep-WH distribution normalized to unity for the shown components of of the signal with $m_{\mathrm{W}_{\mathrm{KK}}}$ = 2.5 TeV, $m_{\mathrm{R}}$ = 0.2 TeV. The $t^{3,4}$ jets from the preselected sample, normalized to unity, are superimposed to compare shapes with the $R^{3q}$ and $R^{4q}$ distributions.

The $m_{\mathrm{jj}}$ distribution for CR1 for data and simulation.

The $m_{\mathrm{jj}}$ distribution for CR2 for data and simulation.

The $m_{\mathrm{jj}}$ distribution for CR3 for data and simulation.

The $m_{\mathrm{jjj}}$ distribution for CR45 for data and simulation.

The $m_{\mathrm{jjj}}$ distribution for CR6 for data and simulation.

Post-fit distributions of the reconstructed triboson system ($m_{\mathrm{jj}}$) in data and simulation for SR1.

Post-fit distributions of the reconstructed triboson system ($m_{\mathrm{jj}}$) in data and simulation for SR2.

Post-fit distributions of the reconstructed triboson system ($m_{\mathrm{jj}}$) in data and simulation for SR3.

Post-fit distributions of the reconstructed triboson system ($m_{\mathrm{jjj}}$) in data and simulation for SR4.

Post-fit distributions of the reconstructed triboson system ($m_{\mathrm{jjj}}$) in data and simulation for SR5.

Post-fit distributions of the reconstructed triboson system ($m_{\mathrm{jjj}}$) in data and simulation for SR6.

Observed upper limits at $95\%$ CL on the signal cross section $\times$ branching fraction as functions of the $m_{\mathrm{W}_{\mathrm{KK}}}$ and $m_{\mathrm{R}}$ resonance masses.

Expected median lower limit contour on the $m_{\mathrm{W}_{\mathrm{KK}}}$ and $m_{\mathrm{R}}$ plane.

Expected $+ 1$ s.d. lower limit contour on the $m_{\mathrm{W}_{\mathrm{KK}}}$ and $m_{\mathrm{R}}$ plane.

Expected - 1 s.d. lower limit contour on the $m_{\mathrm{W}_{\mathrm{KK}}}$ and $m_{\mathrm{R}}$ plane.

Observed lower limit contour on the $m_{\mathrm{W}_{\mathrm{KK}}}$ and $m_{\mathrm{R}}$ plane.

Differential signal yields for various signal hypotheses.

Differential signal yields for various signal hypotheses.

Differential signal yields for various signal hypotheses.

Differential signal yields for various signal hypotheses.

Differential signal yields for various signal hypotheses.

Background and data yields in the control and signal region bins. The prediction before ("prefit") and after the background only fit ("b-only") are given separately.

Background and data yields in the control and signal region bins. The prediction before ("prefit") and after the background only fit ("b-only") are given separately.

Background and data yields in the control and signal region bins. The prediction before ("prefit") and after the background only fit ("b-only") are given separately.

Matrix of covariance coefficients between signal region bins. The coefficients are obtained from the background-only fit to the control regions, and serve as input to the simplified likelihood reinterpretation scheme.

Background prediction and observed data yields in the signal region bins. The background yields are obtained from the background-only fit to the corresponding control regions, and serve as input to the simplified likelihood reinterpretation scheme.

Matrix of covariance coefficients between signal region bins. The coefficients are obtained from the background-only fit to the control regions, and serve as input to the simplified likelihood reinterpretation scheme.

Background prediction and observed data yields in the signal region bins. The background yields are obtained from the background-only fit to the corresponding control regions, and serve as input to the simplified likelihood reinterpretation scheme.

Matrix of covariance coefficients between signal region bins. The coefficients are obtained from the background-only fit to the control regions, and serve as input to the simplified likelihood reinterpretation scheme.

Background prediction and observed data yields in the signal region bins. The background yields are obtained from the background-only fit to the corresponding control regions, and serve as input to the simplified likelihood reinterpretation scheme.

Unweighted signal acceptance times efficiency at every cut stage. The requirements called "HCAL mitigation" refer to the requirements imposed in the 2018 data set in order to mitigate the localized failure of the HCAL detector.

Unweighted signal acceptance times efficiency at every cut stage. The requirements called "HCAL mitigation" refer to the requirements imposed in the 2018 data set in order to mitigate the localized failure of the HCAL detector.

Unweighted signal acceptance times efficiency at every cut stage. The requirements called "HCAL mitigation" refer to the requirements imposed in the 2018 data set in order to mitigate the localized failure of the HCAL detector.

Unweighted signal acceptance times efficiency at every cut stage. The requirements called "HCAL mitigation" refer to the requirements imposed in the 2018 data set in order to mitigate the localized failure of the HCAL detector.

Unweighted signal acceptance times efficiency at every cut stage. The requirements called "HCAL mitigation" refer to the requirements imposed in the 2018 data set in order to mitigate the localized failure of the HCAL detector.

Unweighted signal acceptance times efficiency at every cut stage. The requirements called "HCAL mitigation" refer to the requirements imposed in the 2018 data set in order to mitigate the localized failure of the HCAL detector.

Unweighted signal acceptance times efficiency at every cut stage. The requirements called "HCAL mitigation" refer to the requirements imposed in the 2018 data set in order to mitigate the localized failure of the HCAL detector.

Unweighted signal acceptance times efficiency at every cut stage. The requirements called "HCAL mitigation" refer to the requirements imposed in the 2018 data set in order to mitigate the localized failure of the HCAL detector.

Unweighted signal acceptance times efficiency at every cut stage. The requirements called "HCAL mitigation" refer to the requirements imposed in the 2018 data set in order to mitigate the localized failure of the HCAL detector.

Unweighted signal acceptance times efficiency at every cut stage. The requirements called "HCAL mitigation" refer to the requirements imposed in the 2018 data set in order to mitigate the localized failure of the HCAL detector.

Unweighted signal acceptance times efficiency at every cut stage. The requirements called "HCAL mitigation" refer to the requirements imposed in the 2018 data set in order to mitigate the localized failure of the HCAL detector.

Unweighted signal acceptance times efficiency at every cut stage. The requirements called "HCAL mitigation" refer to the requirements imposed in the 2018 data set in order to mitigate the localized failure of the HCAL detector.

Unweighted signal acceptance times efficiency at every cut stage. The requirements called "HCAL mitigation" refer to the requirements imposed in the 2018 data set in order to mitigate the localized failure of the HCAL detector.

Unweighted signal acceptance times efficiency at every cut stage. The requirements called "HCAL mitigation" refer to the requirements imposed in the 2018 data set in order to mitigate the localized failure of the HCAL detector.

Unweighted signal acceptance times efficiency at every cut stage. The requirements called "HCAL mitigation" refer to the requirements imposed in the 2018 data set in order to mitigate the localized failure of the HCAL detector.

Unweighted signal acceptance times efficiency at every cut stage. The requirements called "HCAL mitigation" refer to the requirements imposed in the 2018 data set in order to mitigate the localized failure of the HCAL detector.

Unweighted signal acceptance times efficiency at every cut stage. The requirements called "HCAL mitigation" refer to the requirements imposed in the 2018 data set in order to mitigate the localized failure of the HCAL detector.

Unweighted signal acceptance times efficiency at every cut stage. The requirements called "HCAL mitigation" refer to the requirements imposed in the 2018 data set in order to mitigate the localized failure of the HCAL detector.

Median Expected exclusion contour in the $m_{med}$-$m_{\chi}$ plane in the simplified model with axial couplings.

Observed exclusion contour in the $m_{med}$-$m_{\chi}$ plane in the simplified model with axial couplings.

Expected plus 1 s.d. exclusion contour in the $m_{med}$-$m_{\chi}$ plane in the simplified model with axial couplings.

Expected minus 1 s.d. exclusion contour in the $m_{med}$-$m_{\chi}$ plane in the simplified model with axial couplings.

Median Expected exclusion contour in the $m_{med}$-$m_{\chi}$ plane in the simplified model with vector couplings.

Observed exclusion contour in the $m_{med}$-$m_{\chi}$ plane in the simplified model with vector couplings.

Expected plus 1 s.d. exclusion contour in the $m_{med}$-$m_{\chi}$ plane in the simplified model with vector couplings.

Expected minus 1 s.d. exclusion contour in the $m_{med}$-$m_{\chi}$ plane in the simplified model with vector couplings.

Upper limits on the coupling $g_{\chi}$ in the simplified model with a axial mediator.

Upper limits on the coupling $g_{q}$ in the simplified model with a axial mediator.

Upper limits on the coupling $g_{\chi}$ in the simplified model with a vector mediator.

Upper limits on the coupling $g_{q}$ in the simplified model with a vector mediator.

Exclusion limits on the signal strength in the simplified model with scalar couplings.

Exclusion limits on the signal strength in the simplified model with pseudoscalar couplings.

Exclusion limits on the fundamental Planck scale $M_{D}$ as a function of the number of extra dimensions $d$.

Median Expected exclusion contour in the $m_{med}$-$m_{\chi}$ plane in the fermion portal model.

Observed exclusion contour in the $m_{med}$-$m_{\chi}$ plane in the fermion portal model.

Expected plus 1 s.d. exclusion contour in the $m_{med}$-$m_{\chi}$ plane in the fermion portal model.

Expected minus 1 s.d. exclusion contour in the $m_{med}$-$m_{\chi}$ plane in the fermion portal model.

Tagging efficiency for AK8 jets. The efficiency includes the effect of the machine-learning based DeepAK8 tagger, as well as the application of the mass window requirement on the jet. The efficiency is split depending on the matching generator-level object. For reinterpretation purposes, this efficiency can directly be applied to any AK8 jet that is matched to a given type of generator-level object, with no prior selection on the jet mass. Other acceptance requirements on jet $\p_{T}$ and $\eta$ should still be applied.

Simultaneous EW and QCD WV signal strength fit.

The BDT distribution (with x-axis zoomed in) of the di-Higgs ggF signal for two different values of $\kappa_{\lambda}$ and the main backgrounds in the high mass region ($m^{*}_{b\bar{b}\gamma\gamma} > 350$ GeV). Distributions are normalized to unit area. The dotted lines denote the category boundaries. Events with a BDT score below 0.857 in the high mass region are discarded. The range of BDT scores is from 0.8 to 1.

The $BDT_{tot}$ score for the benchmark signal $m_{X}$ = 300 GeV and for the main backgrounds. Distributions are normalized to unit area. The dotted line denotes the event selection threshold. Events with a $BDT_{tot}$ score below 0.85 for $m_{X}$ = 300 GeV are discarded.

The $BDT_{tot}$ score for the benchmark signal $m_{X}$ = 500 GeV and for the main backgrounds. Distributions are normalized to unit area. The dotted line denotes the event selection threshold. Events with a $BDT_{tot}$ score below 0.75 for $m_{X}$ = 500 GeV are discarded.

Distributions of $m_{\gamma\gamma}$ in high mass BDT tight category for the nonresonant $HH$ search. The data-derived fractions of nonresonant $\gamma\gamma$, $\gamma$-jet or jet-$\gamma$, and dijet background are applied and the total background is normalized to the data sideband.

Distributions of $m_{\gamma\gamma}$ in high mass BDT loose category for the nonresonant $HH$ search. The data-derived fractions of nonresonant $\gamma\gamma$, $\gamma$-jet or jet-$\gamma$, and dijet background are applied and the total background is normalized to the data sideband.

Distributions of $m_{\gamma\gamma}$ in low mass BDT tight category for the nonresonant $HH$ search. The data-derived fractions of nonresonant $\gamma\gamma$, $\gamma$-jet or jet-$\gamma$, and dijet background are applied and the total background is normalized to the data sideband.

Distributions of $m_{\gamma\gamma}$ in low mass BDT loose category for the nonresonant $HH$ search. The data-derived fractions of nonresonant $\gamma\gamma$, $\gamma$-jet or jet-$\gamma$, and dijet background are applied and the total background is normalized to the data sideband.

Distributions of $m_{\gamma\gamma}$ for the selections used for the resonance mass point $m_{X}$ = 300 GeV for the resonant search. The data-derived fractions of nonresonant $\gamma\gamma$, $\gamma$-jet or jet-$\gamma$, and dijet background are applied and the total background is normalized to the data sideband. The scalar resonance signal is scaled to a total production cross section $\sigma(pp \rightarrow X \rightarrow HH)$ = 370 fb.

Distributions of $m_{\gamma\gamma}$ for the selections used for the resonance mass point $m_{X}$ = 500 GeV for the resonant search. The data-derived fractions of nonresonant $\gamma\gamma$, $\gamma$-jet or jet-$\gamma$, and dijet background are applied and the total background is normalized to the data sideband. The scalar resonance signal is scaled to a total production cross section $\sigma(pp \rightarrow X \rightarrow HH)$ = 67 fb.

The number of data events observed in the 120 GeV < $m_{\gamma\gamma}$ < 130 GeV window, the number of $HH$ signal events expected for $\kappa_{\lambda}$ = 1 and for $\kappa_{\lambda}$ = 10, and events expected for single Higgs boson production (estimated using MC simulation), as well as for continuum background. For the single Higgs boson, Rest includes VBF, $WH$, $tHqb$, and $tHW$. The values are obtained from a fit of the Asimov data set generated under the SM signal-plus-background hypothesis, $\kappa_{\lambda}$ = 1. The continuum background component of the Asimov data set is obtained from the fit of the data sideband. The uncertainties in $HH$ signals and single Higgs boson background include the systematic uncertainties discussed in Section 6. The uncertainty in the continuum background is given by the sum in quadrature of the statistical uncertainty from the fit to the data and the spurious-signal uncertainty.

Observed and expected limits at 95% CL on the cross section of nonresonant Higgs boson pair production as a function of the Higgs boson self-coupling modifier $\kappa_{\lambda}= \lambda_{HHH}/\lambda^{\textrm{SM}}_{HHH}$. The expected constraints on $\kappa_{\lambda}$ are obtained with a background hypothesis excluding $pp \rightarrow HH$ production. The $\pm 1\sigma$ and $\pm 2\sigma$ variations about the expected limit due to statistical and systematic uncertainties are also shown. The theory prediction curve represents the scenario where all parameters and couplings are set to their SM values except for $\kappa_{\lambda}$. The uncertainty band of the theory prediction curve shows the cross-section uncertainty.

Values of the negative log-profile-likelihood ratio ($-2ln\Lambda$) as a function of $\kappa_{\lambda}$ evaluated for the combination of all the categories of the nonresonant search. The coupling of the Higgs boson to fermions and gauge bosons is set to SM values in the profile likelihood calculation. The Asimov data set is generated under the SM signal-plus-background hypothesis, $\kappa_{\lambda}$= 1. All systematic uncertainties, including the theoretical uncertainties in the di-Higgs boson production cross section, are included. The intersections of the solid curves and the horizontal dashed lines indicate the 1$\sigma$ and 2$\sigma$ confidence-level intervals.

The number of events observed in the 120 < $m_{\gamma\gamma}$ < 130 GeV window in data, the number of events expected for scalar resonance signals of masses $m_{X}$ = 300 GeV and $m_{X}$ = 500 GeV assuming a total production cross section $\sigma(pp \rightarrow X \rightarrow HH)$ equal to the observed exclusion limits of Figure 15, and events expected for SM $HH$ and single Higgs boson production (estimated using MC simulation), as well as for continuum background. The values are obtained from a fit of the Asimov data set generated under the signal-plus-background hypothesis. The continuum background component of the Asimov data set is obtained from the fit of the data sideband. The uncertainties in the resonant signals and the SM $HH$ and single-Higgs-boson backgrounds include the systematic uncertainties discussed in Section 6. The uncertainty in the continuum background is given by the sum in quadrature of the statistical uncertainty from the fit to the data and the spurious-signal uncertainty.

Observed and expected limits at 95% CL on the production cross section of a narrow-width scalar resonance $X$ as a function of the mass $m_{X}$ of the hypothetical scalar particle. The black solid line represents the observed upper limits. The dashed line represents the expected upper limits. The $\pm 1\sigma$ and $\pm 2\sigma$ variations about the expected limit due to statistical and systematic uncertainties are also shown.

Breakdown of the dominant systematic uncertainties. The impact of the uncertainties is defined according to the statistical analysis described in Section 7. It corresponds to the relative variation of the expected upper limit on the cross section when re-evaluating the profile likelihood ratio after fixing the nuisance parameter in question to its best-fit value, while all remaining nuisance parameters remain free to float. The impact is shown in %. Only systematic uncertainties with an impact of at least 0.2% are shown. Uncertainties of the "Norm. + Shape" type affect both the normalization and the parameters of the functional form. The rest of the uncertainties affect only the yields.

Cutflow for nonresonant di-Higgs ggF signal sample, yields are normalized to 139 $fb^{-1}$.

Cutflow for resonant signal sample, with $m_{X}$ = 300 GeV, yields are normalized to 139 $fb^{-1}$.

Cutflow for resonant signal sample, with $m_{X}$ = 500 GeV, yields are normalized to 139 $fb^{-1}$.

Comparison of $m_{b\bar{b}}$ distributions when applying the specific b-jet energy calibration and the nominal jet energy calibration. The distributions are fitted using a Bukin function, and the values of the peak position, resolution and the relative improvement are reported in the legend.

Fit results of $m_{b\bar{b}}$ distributions when applying the specific b-jet energy calibration and the nominal jet energy calibration. The distributions are fitted using a Bukin function, and the values of the peak position, resolution and the relative improvement are reported in the legend.

The relative amount (purity) of expected events from SM $HH$ and single Higgs boson production processes for each of the four categories of the nonresonant search. The Higgs boson pair production with $k_{\lambda} = 1$ is considered as signal in (a), while the case with $k_{\lambda} = 10$ is considered as signal in (b).

The expected significance in each of the four categories of the nonresonant search. The Higgs boson pair production with $k_{\lambda} = 1$ is considered as signal in (a), while the case with $k_{\lambda} = 10$ is considered as signal in (b). The single Higgs boson processes and the di-photon continuum spectrum are considered as background.

Spurious signal result for the exponential pdf for the various ggF nonresonant di-Higgs categories. In each category, the spurious signal value ($n_{sp}$) and its ratio to the expected statistical error ($Z_{spur}$) from data are shown.

Spurious signal result for the exponential pdf as function of the resonant di-Higgs signal mass. The spurious signal value and its ratio to the expected statistical error from data are shown.

Distributions of the signal efficiency as a function of $\kappa_{\lambda}$, for the di-Higgs boson ggF nonresonant production mode. The range of $\kappa_{\lambda}$ in the table is from -10 to 10.

Distributions of the signal efficiency as a function of $\kappa_{\lambda}$, for the di-Higgs boson VBF nonresonant production mode. The range of $\kappa_{\lambda}$ in the table is from -10 to 10.

Values of the negative log-profile-likelihood ($-2ln\Lambda$) as a function of $\kappa_{\lambda}$ evaluated for the combination of all the categories of the nonresonant search. The coupling of the Higgs boson to fermions and gauge bosons is set to SM values in the profile likelihood calculation. The Asimov data set is generated under the SM signal-plus-background hypothesis, $\kappa_{\lambda}$= 1. All systematic uncertainties, including the theoretical uncertainties on the di-Higgs boson production cross section, are included. The intersections of the solid curves and the horizontal dashed lines indicate the 1$\sigma$ and 2$\sigma$ confidence level intervals. The best fit value corresponds to $\kappa_{\lambda}$ = $2.8^{+2.0}_{-2.2}(^{+3.8}_{-4.3})$ for the 1$\sigma$(2$\sigma$) confidence interval. The expected value corresponds to $\kappa_{\lambda}$ = $1.0^{+5.5}_{-2.4}(^{+7.3}_{-4.2})$ for the 1$\sigma$(2$\sigma$) confidence interval. The dashed curves represent values of the negative log-profile-likelihood where the Higgs boson branching fractions and the cross section of the production modes are varied as a function of $\kappa_{\lambda}$. In this case,the best fit value corresponds to $\kappa_{\lambda}$ = $2.7^{+2.0}_{-2.2}(^{+3.8}_{-4.3})$ and the expected value corresponds to $\kappa_{\lambda}$ = $1.0^{+5.4}_{-2.5}(^{+7.3}_{-4.3})$ for the 1$\sigma$(2$\sigma$) confidence interval.

Minimum BDT value of the events passing the selection criteria of the resonant search. The combined BDT score is formed using as coefficients $C_{1}$ = 0.65 and $C_{2}$ = 1 − $C_{1}$. The selection efficiency for the resonant $X \rightarrow HH$ signal is also shown.

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.

Normalized differential cross section distribution as a function of the leading b jet transverse momentum for the Z + >= 1 b jet events

Differential cross section distribution as a function of the leading b jet absolute pseudorapidity for the Z + >= 1 b jet events

Normalized differential cross section distribution as a function of the leading b jet absolute pseudorapidity for the Z + >= 1 b jet events

Differential cross section distribution as a function of the leading b jet transverse momentum for Z +>= 1 b jet events

Normalized differential cross section distribution as a function of azimuthal difference between Z boson and the leading b jet for the Z + >= 1 b jet events

Differential cross section distribution as a function of the rapidity difference between Z boson and the leading b jet for the Z + >= 1 b jet events

Normalized differential cross section distribution as a function of the rapidity difference between Z boson and the leading b jet for the Z + >= 1 b jet events

Differential cross section distribution as a function of the angular separation between the Z boson and the leading b jet for the Z + >= 1 b jet events

Normalized differential cross section distribution as a function of the angular separation between the Z boson and the leading b jet for the Z + >= 1 b jet events

Differential cross section distribution as a function of the leading b jet transverse momentum for the Z + >= 2 b jet events

Normalized differential cross section distribution as a function of the leading b jet transverse momentum for the Z + >= 2 b jet events

Differential cross section distribution as a function of the leading b jet absolute pseudorapidity

Normalized differential cross section distribution as a function of the leading b jet absolute pseudorapidity

Differential cross section distribution as a function of the subleading b jet transverse momentum for the Z + >= 2 b jet events

Normalized differential cross section distribution as a function of the subleading b jet transverse momentum for the Z + >= 2 b jet events

Differential cross section distribution as a function of the Z boson transverse momentum for the Z + >= 2 b jet events

Normalized differential cross section distribution as a function of the Z boson transverse momentum for the Z + >= 2 b jet events

Differential cross section as a function of the angular separation between two b jets for the Z + >= 2 b jet events

Normalized differential cross section as a function of the angular separtion between two b jets for the Z + >= 2 b jet events

Differential cross section as a function of the minimum angular separation between the Z boson and two b jets for the Z + >= 2 b jet events

Normalized differential cross section as a function of the minimum angular separation between the Z boson and two b jets for the Z + >= 2 b jet events

Differential cross section as a function of the asymmetry of the Z + >= 2 b jets system

Normalized differential cross section as a function of the asymmetry of the Z + >= 2 b jets system

Differential cross section as a function of the invariant mass of two b jets for the Z + >= 2 b jet events

Normalized differential cross section as a function of the invariant mass of two b jets for the Z + >= 2 b jet events

Differential cross section as a function of the invariant mass of the Z boson and two b jets for the Z + >= 2 b jet events

Normalized differential cross section as a function of invariant mass of the Z boson and two b jets for the Z + >= 2 b jet events

Distributions of the cross section ratios as a function of the leading b jet transverse momentum

Distributions of the cross section ratios as a function of the leading b jet absolute pseudorapidity

The inclusive jet production cross section as a function of the jet transverse momentum~$p_\mathrm{T}$ measured in $1.5 < |y| < 2.0$ for jets clustered using the anti-$k_\mathrm{t}$ algorithm with $R=0.4$.

The correlation matrix at the particle level of inclusive jet production cross section as a function of the jet transverse momentum~$p_\mathrm{T}$ for jets in $|y| < 0.5$ clustered using the anti-$k_\mathrm{t}$ algorithm with $R=0.4$. It contains contributions from the data and from the \textsc{PYTHIA}~8 sample used to perform the unfolding.

The correlation matrix at the particle level of inclusive jet production cross section as a function of the jet transverse momentum~$p_\mathrm{T}$ for jets in $|y| < 0.5$and $0.5 < |y| < 1.0$ clustered using the anti-$k_\mathrm{t}$ algorithm with $R=0.4$. It contains contributions from the data and from the \textsc{PYTHIA}~8 sample used to perform the unfolding.

The correlation matrix at the particle level of inclusive jet production cross section as a function of the jet transverse momentum~$p_\mathrm{T}$ for jets in $|y| < 0.5$and $1.0 < |y| < 1.5$ clustered using the anti-$k_\mathrm{t}$ algorithm with $R=0.4$. It contains contributions from the data and from the \textsc{PYTHIA}~8 sample used to perform the unfolding.

The correlation matrix at the particle level of inclusive jet production cross section as a function of the jet transverse momentum~$p_\mathrm{T}$ for jets in $|y| < 0.5$and $1.5 < |y| < 2.0$ clustered using the anti-$k_\mathrm{t}$ algorithm with $R=0.4$. It contains contributions from the data and from the \textsc{PYTHIA}~8 sample used to perform the unfolding.

The correlation matrix at the particle level of inclusive jet production cross section as a function of the jet transverse momentum~$p_\mathrm{T}$ for jets in $0.5 < |y| < 1.0$and $|y| < 0.5$ clustered using the anti-$k_\mathrm{t}$ algorithm with $R=0.4$. It contains contributions from the data and from the \textsc{PYTHIA}~8 sample used to perform the unfolding.

The correlation matrix at the particle level of inclusive jet production cross section as a function of the jet transverse momentum~$p_\mathrm{T}$ for jets in $0.5 < |y| < 1.0$ clustered using the anti-$k_\mathrm{t}$ algorithm with $R=0.4$. It contains contributions from the data and from the \textsc{PYTHIA}~8 sample used to perform the unfolding.

The correlation matrix at the particle level of inclusive jet production cross section as a function of the jet transverse momentum~$p_\mathrm{T}$ for jets in $0.5 < |y| < 1.0$and $1.0 < |y| < 1.5$ clustered using the anti-$k_\mathrm{t}$ algorithm with $R=0.4$. It contains contributions from the data and from the \textsc{PYTHIA}~8 sample used to perform the unfolding.

The correlation matrix at the particle level of inclusive jet production cross section as a function of the jet transverse momentum~$p_\mathrm{T}$ for jets in $0.5 < |y| < 1.0$and $1.5 < |y| < 2.0$ clustered using the anti-$k_\mathrm{t}$ algorithm with $R=0.4$. It contains contributions from the data and from the \textsc{PYTHIA}~8 sample used to perform the unfolding.

The correlation matrix at the particle level of inclusive jet production cross section as a function of the jet transverse momentum~$p_\mathrm{T}$ for jets in $1.0 < |y| < 1.5$and $|y| < 0.5$ clustered using the anti-$k_\mathrm{t}$ algorithm with $R=0.4$. It contains contributions from the data and from the \textsc{PYTHIA}~8 sample used to perform the unfolding.

The correlation matrix at the particle level of inclusive jet production cross section as a function of the jet transverse momentum~$p_\mathrm{T}$ for jets in $1.0 < |y| < 1.5$and $0.5 < |y| < 1.0$ clustered using the anti-$k_\mathrm{t}$ algorithm with $R=0.4$. It contains contributions from the data and from the \textsc{PYTHIA}~8 sample used to perform the unfolding.

The correlation matrix at the particle level of inclusive jet production cross section as a function of the jet transverse momentum~$p_\mathrm{T}$ for jets in $1.0 < |y| < 1.5$ clustered using the anti-$k_\mathrm{t}$ algorithm with $R=0.4$. It contains contributions from the data and from the \textsc{PYTHIA}~8 sample used to perform the unfolding.

The correlation matrix at the particle level of inclusive jet production cross section as a function of the jet transverse momentum~$p_\mathrm{T}$ for jets in $1.0 < |y| < 1.5$and $1.5 < |y| < 2.0$ clustered using the anti-$k_\mathrm{t}$ algorithm with $R=0.4$. It contains contributions from the data and from the \textsc{PYTHIA}~8 sample used to perform the unfolding.

The correlation matrix at the particle level of inclusive jet production cross section as a function of the jet transverse momentum~$p_\mathrm{T}$ for jets in $1.5 < |y| < 2.0$and $|y| < 0.5$ clustered using the anti-$k_\mathrm{t}$ algorithm with $R=0.4$. It contains contributions from the data and from the \textsc{PYTHIA}~8 sample used to perform the unfolding.

The correlation matrix at the particle level of inclusive jet production cross section as a function of the jet transverse momentum~$p_\mathrm{T}$ for jets in $1.5 < |y| < 2.0$and $0.5 < |y| < 1.0$ clustered using the anti-$k_\mathrm{t}$ algorithm with $R=0.4$. It contains contributions from the data and from the \textsc{PYTHIA}~8 sample used to perform the unfolding.

The correlation matrix at the particle level of inclusive jet production cross section as a function of the jet transverse momentum~$p_\mathrm{T}$ for jets in $1.5 < |y| < 2.0$and $1.0 < |y| < 1.5$ clustered using the anti-$k_\mathrm{t}$ algorithm with $R=0.4$. It contains contributions from the data and from the \textsc{PYTHIA}~8 sample used to perform the unfolding.

The correlation matrix at the particle level of inclusive jet production cross section as a function of the jet transverse momentum~$p_\mathrm{T}$ for jets in $1.5 < |y| < 2.0$ clustered using the anti-$k_\mathrm{t}$ algorithm with $R=0.4$. It contains contributions from the data and from the \textsc{PYTHIA}~8 sample used to perform the unfolding.

The inclusive jet production cross section as a function of the jet transverse momentum~$p_\mathrm{T}$ measured in $|y| < 0.5$ for jets clustered using the anti-$k_\mathrm{t}$ algorithm with $R=0.7$.

The inclusive jet production cross section as a function of the jet transverse momentum~$p_\mathrm{T}$ measured in $0.5 < |y| < 1.0$ for jets clustered using the anti-$k_\mathrm{t}$ algorithm with $R=0.7$.

The inclusive jet production cross section as a function of the jet transverse momentum~$p_\mathrm{T}$ measured in $1.0 < |y| < 1.5$ for jets clustered using the anti-$k_\mathrm{t}$ algorithm with $R=0.7$.

The inclusive jet production cross section as a function of the jet transverse momentum~$p_\mathrm{T}$ measured in $1.5 < |y| < 2.0$ for jets clustered using the anti-$k_\mathrm{t}$ algorithm with $R=0.7$.

The correlation matrix at the particle level of inclusive jet production cross section as a function of the jet transverse momentum~$p_\mathrm{T}$ for jets in $|y| < 0.5$ clustered using the anti-$k_\mathrm{t}$ algorithm with $R=0.7$. It contains contributions from the data and from the \textsc{PYTHIA}~8 sample used to perform the unfolding.

The correlation matrix at the particle level of inclusive jet production cross section as a function of the jet transverse momentum~$p_\mathrm{T}$ for jets in $|y| < 0.5$and $0.5 < |y| < 1.0$ clustered using the anti-$k_\mathrm{t}$ algorithm with $R=0.7$. It contains contributions from the data and from the \textsc{PYTHIA}~8 sample used to perform the unfolding.

The correlation matrix at the particle level of inclusive jet production cross section as a function of the jet transverse momentum~$p_\mathrm{T}$ for jets in $|y| < 0.5$and $1.0 < |y| < 1.5$ clustered using the anti-$k_\mathrm{t}$ algorithm with $R=0.7$. It contains contributions from the data and from the \textsc{PYTHIA}~8 sample used to perform the unfolding.

The correlation matrix at the particle level of inclusive jet production cross section as a function of the jet transverse momentum~$p_\mathrm{T}$ for jets in $|y| < 0.5$and $1.5 < |y| < 2.0$ clustered using the anti-$k_\mathrm{t}$ algorithm with $R=0.7$. It contains contributions from the data and from the \textsc{PYTHIA}~8 sample used to perform the unfolding.

The correlation matrix at the particle level of inclusive jet production cross section as a function of the jet transverse momentum~$p_\mathrm{T}$ for jets in $0.5 < |y| < 1.0$and $|y| < 0.5$ clustered using the anti-$k_\mathrm{t}$ algorithm with $R=0.7$. It contains contributions from the data and from the \textsc{PYTHIA}~8 sample used to perform the unfolding.

The correlation matrix at the particle level of inclusive jet production cross section as a function of the jet transverse momentum~$p_\mathrm{T}$ for jets in $0.5 < |y| < 1.0$ clustered using the anti-$k_\mathrm{t}$ algorithm with $R=0.7$. It contains contributions from the data and from the \textsc{PYTHIA}~8 sample used to perform the unfolding.

The correlation matrix at the particle level of inclusive jet production cross section as a function of the jet transverse momentum~$p_\mathrm{T}$ for jets in $0.5 < |y| < 1.0$and $1.0 < |y| < 1.5$ clustered using the anti-$k_\mathrm{t}$ algorithm with $R=0.7$. It contains contributions from the data and from the \textsc{PYTHIA}~8 sample used to perform the unfolding.

The correlation matrix at the particle level of inclusive jet production cross section as a function of the jet transverse momentum~$p_\mathrm{T}$ for jets in $0.5 < |y| < 1.0$and $1.5 < |y| < 2.0$ clustered using the anti-$k_\mathrm{t}$ algorithm with $R=0.7$. It contains contributions from the data and from the \textsc{PYTHIA}~8 sample used to perform the unfolding.

The correlation matrix at the particle level of inclusive jet production cross section as a function of the jet transverse momentum~$p_\mathrm{T}$ for jets in $1.0 < |y| < 1.5$and $|y| < 0.5$ clustered using the anti-$k_\mathrm{t}$ algorithm with $R=0.7$. It contains contributions from the data and from the \textsc{PYTHIA}~8 sample used to perform the unfolding.

The correlation matrix at the particle level of inclusive jet production cross section as a function of the jet transverse momentum~$p_\mathrm{T}$ for jets in $1.0 < |y| < 1.5$and $0.5 < |y| < 1.0$ clustered using the anti-$k_\mathrm{t}$ algorithm with $R=0.7$. It contains contributions from the data and from the \textsc{PYTHIA}~8 sample used to perform the unfolding.

The correlation matrix at the particle level of inclusive jet production cross section as a function of the jet transverse momentum~$p_\mathrm{T}$ for jets in $1.0 < |y| < 1.5$ clustered using the anti-$k_\mathrm{t}$ algorithm with $R=0.7$. It contains contributions from the data and from the \textsc{PYTHIA}~8 sample used to perform the unfolding.

The correlation matrix at the particle level of inclusive jet production cross section as a function of the jet transverse momentum~$p_\mathrm{T}$ for jets in $1.0 < |y| < 1.5$and $1.5 < |y| < 2.0$ clustered using the anti-$k_\mathrm{t}$ algorithm with $R=0.7$. It contains contributions from the data and from the \textsc{PYTHIA}~8 sample used to perform the unfolding.

The correlation matrix at the particle level of inclusive jet production cross section as a function of the jet transverse momentum~$p_\mathrm{T}$ for jets in $1.5 < |y| < 2.0$and $|y| < 0.5$ clustered using the anti-$k_\mathrm{t}$ algorithm with $R=0.7$. It contains contributions from the data and from the \textsc{PYTHIA}~8 sample used to perform the unfolding.

The correlation matrix at the particle level of inclusive jet production cross section as a function of the jet transverse momentum~$p_\mathrm{T}$ for jets in $1.5 < |y| < 2.0$and $0.5 < |y| < 1.0$ clustered using the anti-$k_\mathrm{t}$ algorithm with $R=0.7$. It contains contributions from the data and from the \textsc{PYTHIA}~8 sample used to perform the unfolding.

The correlation matrix at the particle level of inclusive jet production cross section as a function of the jet transverse momentum~$p_\mathrm{T}$ for jets in $1.5 < |y| < 2.0$and $1.0 < |y| < 1.5$ clustered using the anti-$k_\mathrm{t}$ algorithm with $R=0.7$. It contains contributions from the data and from the \textsc{PYTHIA}~8 sample used to perform the unfolding.

The correlation matrix at the particle level of inclusive jet production cross section as a function of the jet transverse momentum~$p_\mathrm{T}$ for jets in $1.5 < |y| < 2.0$ clustered using the anti-$k_\mathrm{t}$ algorithm with $R=0.7$. It contains contributions from the data and from the \textsc{PYTHIA}~8 sample used to perform the unfolding.

The corrections were calculated by Dittmaier, Huss, Speckner.

The corrections were calculated by Dittmaier, Huss, Speckner.

The corrections were calculated by Dittmaier, Huss, Speckner.

The corrections were calculated by Dittmaier, Huss, Speckner.

The corrections were calculated by Dittmaier, Huss, Speckner.

The corrections were calculated by Dittmaier, Huss, Speckner.

The corrections were calculated by Dittmaier, Huss, Speckner.

The corrections were calculated by Dittmaier, Huss, Speckner.

The values correspond to the average of the corrections obtained with \textsc{{PYTHIA}}~8 and with \textsc{{HERWIG}}++.

The values correspond to the average of the corrections obtained with \textsc{{PYTHIA}}~8 and with \textsc{{HERWIG}}++.

The values correspond to the average of the corrections obtained with \textsc{{PYTHIA}}~8 and with \textsc{{HERWIG}}++.

The values correspond to the average of the corrections obtained with \textsc{{PYTHIA}}~8 and with \textsc{{HERWIG}}++.

The values correspond to the average of the corrections obtained with \textsc{{PYTHIA}}~8 and with \textsc{{HERWIG}}++.

The values correspond to the average of the corrections obtained with \textsc{{PYTHIA}}~8 and with \textsc{{HERWIG}}++.

The values correspond to the average of the corrections obtained with \textsc{{PYTHIA}}~8 and with \textsc{{HERWIG}}++.

The values correspond to the average of the corrections obtained with \textsc{{PYTHIA}}~8 and with \textsc{{HERWIG}}++.

Pre-fit distribution of the reconstructed Higgs boson candidate $p_T^H$ for the dilepton $SR^{\geq 4j}_{\geq 4b}$ signal region. The dashed line shows the ${t\bar {t}H}$ signal distribution normalised to the total background prediction. The uncertainty band includes all uncertainties and their correlations, except for the uncertainty in the $k({t\bar {t}+{\geq }1b})$ normalisation factor which is not defined pre-fit. The last bin includes the overflow.

Pre-fit distribution of the reconstructed Higgs boson candidate $p_T^H$ for the single-lepton resolved $SR^{\geq 6j}_{\geq 4b}$ signal region. The dashed line shows the ${t\bar {t}H}$ signal distribution normalised to the total background prediction. The uncertainty band includes all uncertainties and their correlations, except for the uncertainty in the $k({t\bar {t}+{\geq }1b})$ normalisation factor which is not defined pre-fit. The last bin includes the overflow.

Pre-fit distribution of the reconstructed Higgs boson candidate $p_T^H$ for the single-lepton boosted ${{SR}_{{boosted}}}$ signal region. The dashed line shows the ${t\bar {t}H}$ signal distribution normalised to the total background prediction. The uncertainty band includes all uncertainties and their correlations, except for the uncertainty in the $k({t\bar {t}+{\geq }1b})$ normalisation factor which is not defined pre-fit. The last bin includes the overflow.

Comparison of predicted and observed event yields in each of the control and signal regions in the dilepton channel after the fit to the data. The uncertainty band includes all uncertainties and their correlations.

Comparison of predicted and observed event yields in each of the control and signal regions in the single-lepton channels after the fit to the data. The uncertainty band includes all uncertainties and their correlations.

Comparison between data and prediction for the BDT discriminant in the dilepton SRs after the inclusive fit to the data for $0\le p_T^H<120$ GeV. The ${t\bar {t}H}$ signal yield (solid red) is normalised to the fitted $\mu $ value from the inclusive fit. The dashed line shows the ${t\bar {t}H}$ signal distribution normalised to the total background prediction. The uncertainty band includes all uncertainties and their correlations.

Comparison between data and prediction for the BDT discriminant in the dilepton SRs after the inclusive fit to the data for $120\le p_T^H<200$ GeV. The ${t\bar {t}H}$ signal yield (solid red) is normalised to the fitted $\mu $ value from the inclusive fit. The dashed line shows the ${t\bar {t}H}$ signal distribution normalised to the total background prediction. The uncertainty band includes all uncertainties and their correlations.

Comparison between data and prediction for the BDT discriminant in the dilepton SRs after the inclusive fit to the data for $200\le p_T^H<300$ GeV. The ${t\bar {t}H}$ signal yield (solid red) is normalised to the fitted $\mu $ value from the inclusive fit. The dashed line shows the ${t\bar {t}H}$ signal distribution normalised to the total background prediction. The uncertainty band includes all uncertainties and their correlations.

Comparison between data and prediction for the BDT discriminant in the dilepton SRs after the inclusive fit to the data for $p_{{T}}^{H}\ge 300$ GeV. The ${t\bar {t}H}$ signal yield (solid red) is normalised to the fitted $\mu $ value from the inclusive fit. The dashed line shows the ${t\bar {t}H}$ signal distribution normalised to the total background prediction. The uncertainty band includes all uncertainties and their correlations.

Comparison between data and prediction for the BDT discriminant in the single-lepton resolved SRs after the inclusive fit to the data for $0\le p_T^H<120$ GeV. The ${t\bar {t}H}$ signal yield (solid red) is normalised to the fitted $\mu $ value from the inclusive fit. The dashed line shows the ${t\bar {t}H}$ signal distribution normalised to the total background prediction. The uncertainty band includes all uncertainties and their correlations.

Comparison between data and prediction for the BDT discriminant in the single-lepton resolved SRs after the inclusive fit to the data for $120\le p_T^H<200$ GeV. The ${t\bar {t}H}$ signal yield (solid red) is normalised to the fitted $\mu $ value from the inclusive fit. The dashed line shows the ${t\bar {t}H}$ signal distribution normalised to the total background prediction. The uncertainty band includes all uncertainties and their correlations.

Comparison between data and prediction for the BDT discriminant in the single-lepton resolved SRs after the inclusive fit to the data for $200\le p_T^H<300$ GeV. The ${t\bar {t}H}$ signal yield (solid red) is normalised to the fitted $\mu $ value from the inclusive fit. The dashed line shows the ${t\bar {t}H}$ signal distribution normalised to the total background prediction. The uncertainty band includes all uncertainties and their correlations.

Comparison between data and prediction for the BDT discriminant in the single-lepton resolved SRs after the inclusive fit to the data for $300\le p_T^H<450$ GeV. The ${t\bar {t}H}$ signal yield (solid red) is normalised to the fitted $\mu $ value from the inclusive fit. The dashed line shows the ${t\bar {t}H}$ signal distribution normalised to the total background prediction. The uncertainty band includes all uncertainties and their correlations.

Comparison between data and prediction for the BDT discriminant in the single-lepton resolved SRs after the inclusive fit to the data for $p_{{T}}^{H}\ge 450$ GeV (yield only). The ${t\bar {t}H}$ signal yield (solid red) is normalised to the fitted $\mu $ value from the inclusive fit. The dashed line shows the ${t\bar {t}H}$ signal distribution normalised to the total background prediction. The uncertainty band includes all uncertainties and their correlations.

Comparison between data and prediction for the BDT discriminant in the single-lepton boosted SRs after the inclusive fit to the data for $300\le p_T^H<450$ GeV. The ${t\bar {t}H}$ signal yield (solid red) is normalised to the fitted $\mu $ value from the inclusive fit. The dashed line shows the ${t\bar {t}H}$ signal distribution normalised to the total background prediction. The uncertainty band includes all uncertainties and their correlations.

Comparison between data and prediction for the BDT discriminant in the single-lepton boosted SRs after the inclusive fit to the data for $p_{{T}}^{H}\ge 450$ GeV. The ${t\bar {t}H}$ signal yield (solid red) is normalised to the fitted $\mu $ value from the inclusive fit. The dashed line shows the ${t\bar {t}H}$ signal distribution normalised to the total background prediction. The uncertainty band includes all uncertainties and their correlations.

Comparison between data and prediction for ${\Delta R^{{avg}}_{bb}}$ after the inclusive fit to the data in the single-lepton $CR^{5j}_{{\geq}4b\ lo}$ control region. The ${t\bar {t}H}$ signal yield (solid red) is normalised to the fitted $\mu $ value from the inclusive fit. The dashed line shows the ${t\bar {t}H}$ signal distribution normalised to the total background prediction. The uncertainty band includes all uncertainties and their correlations. The first (last) bin includes the underflow (overflow).

Comparison between data and prediction for ${\Delta R^{{avg}}_{bb}}$ after the inclusive fit to the data in the single-lepton $CR^{5j}_{{\geq}4b\ hi}$ control region. The ${t\bar {t}H}$ signal yield (solid red) is normalised to the fitted $\mu $ value from the inclusive fit. The dashed line shows the ${t\bar {t}H}$ signal distribution normalised to the total background prediction. The uncertainty band includes all uncertainties and their correlations. The first (last) bin includes the underflow (overflow).

Post-fit yields of signal ($S$) and total background ($B$) as a function of $\log (S/B)$, compared with data. Final-discriminant bins in all dilepton and single-lepton analysis regions are combined into bins of $\log (S/B)$, with the signal normalised to the SM prediction used for the computation of $\log (S/B)$. The signal is then shown normalised to the best-fit value and the SM prediction. The lower frame reports the ratio of data to background, and this is compared with the expected ${t\bar {t}H}$-signal-plus-background yield divided by the background-only yield for the best-fit signal strength (solid red line) and the SM prediction (dashed orange line).

Comparison between data and prediction for the reconstruction BDT score for the Higgs boson candidate identified using Higgs boson information, after the inclusive fit to the data in the dilepton resolved channel for $0\le p_T^H<120$ GeV. The ${t\bar {t}H}$ signal yield (solid red) is normalised to the fitted $\mu $ value from the inclusive fit. The dashed line shows the ${t\bar {t}H}$ signal distribution normalised to the total background prediction. The uncertainty band includes all uncertainties and their correlations.

Comparison between data and prediction for the average $\Delta \eta $ between $b$-tagged jets, after the inclusive fit to the data in the dilepton resolved channel for $0\le p_T^H<120$ GeV. The ${t\bar {t}H}$ signal yield (solid red) is normalised to the fitted $\mu $ value from the inclusive fit. The dashed line shows the ${t\bar {t}H}$ signal distribution normalised to the total background prediction. The uncertainty band includes all uncertainties and their correlations.

Comparison between data and prediction for the likelihood discriminant, after the inclusive fit to the data in the single-lepton resolved channel for $0\le p_T^H<120$ GeV. The ${t\bar {t}H}$ signal yield (solid red) is normalised to the fitted $\mu $ value from the inclusive fit. The dashed line shows the ${t\bar {t}H}$ signal distribution normalised to the total background prediction. The uncertainty band includes all uncertainties and their correlations.

Comparison between data and prediction for the average $\Delta R$ for all possible combinations of $b$-tagged jet pairs, after the inclusive fit to the data in the single-lepton resolved channel for $0\le p_T^H<120$ GeV. The ${t\bar {t}H}$ signal yield (solid red) is normalised to the fitted $\mu $ value from the inclusive fit. The dashed line shows the ${t\bar {t}H}$ signal distribution normalised to the total background prediction. The uncertainty band includes all uncertainties and their correlations.

Comparison between data and prediction for the DNN $P(H)$ output for the Higgs boson candidate after the inclusive fit to the data in the single-lepton boosted channel for $300\le p_T^H<450$ GeV. The ${t\bar {t}H}$ signal yield (solid red) is normalised to the fitted $\mu $ value from the inclusive fit. The dashed line shows the ${t\bar {t}H}$ signal distribution normalised to the total background prediction. The uncertainty band includes all uncertainties and their correlations.

Comparison between data and prediction for the DNN $P(H)$ output for the Higgs boson candidate after the inclusive fit to the data in the single-lepton boosted channel for $p_{{T}}^{H}\ge 450$ GeV. The ${t\bar {t}H}$ signal yield (solid red) is normalised to the fitted $\mu $ value from the inclusive fit. The dashed line shows the ${t\bar {t}H}$ signal distribution normalised to the total background prediction. The uncertainty band includes all uncertainties and their correlations.

Post-fit distribution of the reconstructed Higgs boson candidate mass for the dilepton $SR^{\geq 4j}_{\geq 4b}$ signal region. The ${t\bar {t}H}$ signal yield (solid red) is normalised to the fitted $\mu $ value from the inclusive fit. The dashed line shows the ${t\bar {t}H}$ signal distribution normalised to the total background prediction. The uncertainty band includes all uncertainties and their correlations. The first (last) bin includes the underflow (overflow).

Post-fit distribution of the reconstructed Higgs boson candidate mass for the single-lepton resolved $SR^{\geq 6j}_{\geq 4b}$ signal region. The ${t\bar {t}H}$ signal yield (solid red) is normalised to the fitted $\mu $ value from the inclusive fit. The dashed line shows the ${t\bar {t}H}$ signal distribution normalised to the total background prediction. The uncertainty band includes all uncertainties and their correlations. The first (last) bin includes the underflow (overflow).

Post-fit distribution of the reconstructed Higgs boson candidate mass for the single-lepton boosted ${{SR}_{{boosted}}}$ signal region. The ${t\bar {t}H}$ signal yield (solid red) is normalised to the fitted $\mu $ value from the inclusive fit. The dashed line shows the ${t\bar {t}H}$ signal distribution normalised to the total background prediction. The uncertainty band includes all uncertainties and their correlations. The first (last) bin includes the underflow (overflow).

Fitted values of the ${t\bar {t}H}$ signal strength parameter in the individual channels and in the inclusive signal-strength measurement.

Ranking of the 20 nuisance parameters with the largest post-fit impact on $\mu $ in the fit. Nuisance parameters corresponding to statistical uncertainties in the simulated event samples are not included. The empty blue rectangles correspond to the pre-fit impact on $\mu $ and the filled blue ones to the post-fit impact on $\mu $, both referring to the upper scale. The impact of each nuisance parameter, $\Delta \mu $, is computed by comparing the nominal best-fit value of $\mu $ with the result of the fit when fixing the considered nuisance parameter to its best-fit value, $\hat{\theta }$, shifted by its pre-fit (post-fit) uncertainties $\pm \Delta \theta $ ($\pm \Delta \hat{\theta }$). The black points show the pulls of the nuisance parameters relative to their nominal values, $\theta _0$. These pulls and their relative post-fit errors, $\Delta \hat{\theta }/\Delta \theta $, refer to the lower scale. The `ljets' (`dilep') label refers to the single-lepton (dilepton) channel.

Pre-fit distribution of the number of jets in the dilepton $SR^{\geq 4j}_{\geq 4b}$ signal region. The ${t\bar {t}H}$ signal yield (solid red) is normalised to the Standard Model expectation. The dashed line shows the ${t\bar {t}H}$ signal distribution normalised to the total background prediction. The uncertainty band includes all uncertainties and their correlations, except the uncertainty in the $k({t\bar {t}+{\geq }1b})$ normalisation factor that is not defined pre-fit.

Pre-fit distribution of the number of jets in the single-lepton resolved $SR^{\geq 6j}_{\geq 4b}$ signal region. The ${t\bar {t}H}$ signal yield (solid red) is normalised to the Standard Model expectation. The dashed line shows the ${t\bar {t}H}$ signal distribution normalised to the total background prediction. The uncertainty band includes all uncertainties and their correlations, except the uncertainty in the $k({t\bar {t}+{\geq }1b})$ normalisation factor that is not defined pre-fit.

Pre-fit distribution of the number of jets in the single-lepton boosted ${{SR}_{{boosted}}}$ signal region. The ${t\bar {t}H}$ signal yield (solid red) is normalised to the Standard Model expectation. The dashed line shows the ${t\bar {t}H}$ signal distribution normalised to the total background prediction. The uncertainty band includes all uncertainties and their correlations, except the uncertainty in the $k({t\bar {t}+{\geq }1b})$ normalisation factor that is not defined pre-fit.

Post-fit distribution of the number of jets in the dilepton $SR^{\geq 4j}_{\geq 4b}$ signal region. The ${t\bar {t}H}$ signal yield (solid red) is normalised to the fitted $\mu $ value from the inclusive fit. The dashed line shows the ${t\bar {t}H}$ signal distribution normalised to the total background prediction. The uncertainty band includes all uncertainties and their correlations.

Post-fit distribution of the number of jets in the single-lepton resolved $SR^{\geq 6j}_{\geq 4b}$ signal region. The ${t\bar {t}H}$ signal yield (solid red) is normalised to the fitted $\mu $ value from the inclusive fit. The dashed line shows the ${t\bar {t}H}$ signal distribution normalised to the total background prediction. The uncertainty band includes all uncertainties and their correlations.

Post-fit distribution of the number of jets in the single-lepton boosted ${{SR}_{{boosted}}}$ signal region. The ${t\bar {t}H}$ signal yield (solid red) is normalised to the fitted $\mu $ value from the inclusive fit. The dashed line shows the ${t\bar {t}H}$ signal distribution normalised to the total background prediction. The uncertainty band includes all uncertainties and their correlations.

Post-fit distribution of the reconstructed Higgs boson candidate $p_T^H$ for the dilepton $SR^{\geq 4j}_{\geq 4b}$ signal region. The ${t\bar {t}H}$ signal yield (solid red) is normalised to the fitted $\mu $ value from the inclusive fit. The dashed line shows the ${t\bar {t}H}$ signal distribution normalised to the total background prediction. The uncertainty band includes all uncertainties and their correlations. The last bin includes the overflow.

Post-fit distribution of the reconstructed Higgs boson candidate $p_T^H$ for the single-lepton resolved $SR^{\geq 6j}_{\geq 4b}$ signal region. The ${t\bar {t}H}$ signal yield (solid red) is normalised to the fitted $\mu $ value from the inclusive fit. The dashed line shows the ${t\bar {t}H}$ signal distribution normalised to the total background prediction. The uncertainty band includes all uncertainties and their correlations. The last bin includes the overflow.

Post-fit distribution of the reconstructed Higgs boson candidate $p_T^H$ for the single-lepton boosted ${{SR}_{{boosted}}}$ signal region. The ${t\bar {t}H}$ signal yield (solid red) is normalised to the fitted $\mu $ value from the inclusive fit. The dashed line shows the ${t\bar {t}H}$ signal distribution normalised to the total background prediction. The uncertainty band includes all uncertainties and their correlations. The last bin includes the overflow.

Signal-strength measurements in the individual STXS ${\hat{p}_{{T}}^{H}}$ bins, as well as the inclusive signal strength.

95% CL simplified template cross-section upper limits in the individual STXS ${\hat{p}_{{T}}^{H}}$ bins, as well as the inclusive limit. The observed limits are shown (solid black lines), together with the expected limits both in the background-only hypothesis (dotted black lines) and in the SM hypothesis (dotted red lines). In the case of the expected limits in the background-only hypothesis, one- and two-standard-deviation uncertainty bands are also shown. The hatched uncertainty bands correspond to the theory uncertainty in the fiducial cross-section prediction in each bin.

The ratios $S/B$ (black solid line, referring to the vertical axis on the left) and $S/\sqrt{B}$ (red dashed line, referring to the vertical axis on the right) for each category in the inclusive analysis in the dilepton channel (left) and in the single-lepton channels (right), where $S$ ($B$) is the number of selected signal (background) events predicted by the simulation and normalised to a luminosity of 139 fb$^{-1}$ .

Comparison between data and prediction for the $\Delta R$ between the Higgs candidate and the ${t\bar {t}}$ candidate system, after the inclusive fit to the data in the dilepton resolved channel for $0\le p_T^H<120$ GeV. The ${t\bar {t}H}$ signal yield (solid red) is normalised to the fitted $\mu $ value from the inclusive fit. The dashed line shows the ${t\bar {t}H}$ signal distribution normalised to the total background prediction. The uncertainty band includes all uncertainties and their correlations.

Comparison between data and prediction for the number of $b$-tagged jet pairs with an invariant mass within 30 GeV of 125 GeV, after the inclusive fit to the data in the dilepton resolved channel for $0\le p_T^H<120$ GeV. The ${t\bar {t}H}$ signal yield (solid red) is normalised to the fitted $\mu $ value from the inclusive fit. The dashed line shows the ${t\bar {t}H}$ signal distribution normalised to the total background prediction. The uncertainty band includes all uncertainties and their correlations.

Comparison between data and prediction for the reconstruction BDT score for the Higgs boson candidate identified using Higgs boson information, after the inclusive fit to the data in the single-lepton resolved channel for $0\le p_T^H<120$ GeV. The ${t\bar {t}H}$ signal yield (solid red) is normalised to the fitted $\mu $ value from the inclusive fit. The dashed line shows the ${t\bar {t}H}$ signal distribution normalised to the total background prediction. The uncertainty band includes all uncertainties and their correlations.

Comparison between data and prediction for the $\Delta R$ between the two highest ${p_{{T}}}$ $b$-tagged jets, after the inclusive fit to the data in the single-lepton resolved channel for $0\le p_T^H<120$ GeV. The ${t\bar {t}H}$ signal yield (solid red) is normalised to the fitted $\mu $ value from the inclusive fit. The dashed line shows the ${t\bar {t}H}$ signal distribution normalised to the total background prediction. The uncertainty band includes all uncertainties and their correlations.

Comparison between data and prediction for the sum of $b$-tagging discriminants of jets from Higgs, hadronic top and leptonic top candidates, after the inclusive fit to the data in the single-lepton boosted channel for $300\le p_T^H<450$ GeV. The ${t\bar {t}H}$ signal yield (solid red) is normalised to the fitted $\mu $ value from the inclusive fit. The dashed line shows the ${t\bar {t}H}$ signal distribution normalised to the total background prediction. The uncertainty band includes all uncertainties and their correlations. The first (last) bin includes the underflow (overflow).

Comparison between data and prediction for the sum of $b$-tagging discriminants of jets from Higgs, hadronic top and leptonic top candidates, after the inclusive fit to the data in the single-lepton boosted channel for $p_{{T}}^{H}\ge 450$ GeV. The ${t\bar {t}H}$ signal yield (solid red) is normalised to the fitted $\mu $ value from the inclusive fit. The dashed line shows the ${t\bar {t}H}$ signal distribution normalised to the total background prediction. The uncertainty band includes all uncertainties and their correlations. The first (last) bin includes the underflow (overflow).

Comparison between data and prediction for the hadronic top candidate invariant mass, after the inclusive fit to the data in the single-lepton boosted channel for $300\le p_T^H<450$ GeV. The ${t\bar {t}H}$ signal yield (solid red) is normalised to the fitted $\mu $ value from the inclusive fit. The dashed line shows the ${t\bar {t}H}$ signal distribution normalised to the total background prediction. The uncertainty band includes all uncertainties and their correlations. The first (last) bin includes the underflow (overflow).

Comparison between data and prediction for the hadronic top candidate invariant mass, after the inclusive fit to the data in the single-lepton boosted channel for $p_{{T}}^{H}\ge 450$ GeV. The ${t\bar {t}H}$ signal yield (solid red) is normalised to the fitted $\mu $ value from the inclusive fit. The dashed line shows the ${t\bar {t}H}$ signal distribution normalised to the total background prediction. The uncertainty band includes all uncertainties and their correlations. The first (last) bin includes the underflow (overflow).

Comparison between data and prediction for the fraction of the sum of $b$-tagging discriminants of all jets not associated to the Higgs or hadronic top candidates, after the inclusive fit to the data in the single-lepton boosted channel for $300\le p_T^H<450$ GeV. The ${t\bar {t}H}$ signal yield (solid red) is normalised to the fitted $\mu $ value from the inclusive fit. The dashed line shows the ${t\bar {t}H}$ signal distribution normalised to the total background prediction. The uncertainty band includes all uncertainties and their correlations. The first (last) bin includes the underflow (overflow).

Comparison between data and prediction for the fraction of the sum of $b$-tagging discriminants of all jets not associated to the Higgs or hadronic top candidates, after the inclusive fit to the data in the single-lepton boosted channel for $p_{{T}}^{H}\ge 450$ GeV. The ${t\bar {t}H}$ signal yield (solid red) is normalised to the fitted $\mu $ value from the inclusive fit. The dashed line shows the ${t\bar {t}H}$ signal distribution normalised to the total background prediction. The uncertainty band includes all uncertainties and their correlations. The first (last) bin includes the underflow (overflow).

Ranking of the 20 nuisance parameters with the largest post-fit impact on $\mu $ in the STXS fit for $0\le {\hat{p}_{{T}}^{H}}<120$ GeV. Nuisance parameters corresponding to statistical uncertainties in the simulated event samples are not included. The empty blue rectangles correspond to the pre-fit impact on $\mu $ and the filled blue ones to the post-fit impact on $\mu $, both referring to the upper scale. The impact of each nuisance parameter, $\Delta \mu $, is computed by comparing the nominal best-fit value of $\mu $ with the result of the fit when fixing the considered nuisance parameter to its best-fit value, $\hat{\theta }$, shifted by its pre-fit (post-fit) uncertainties $\pm \Delta \theta $ ($\pm \Delta \hat{\theta }$). The black points show the pulls of the nuisance parameters relative to their nominal values, $\theta _0$. These pulls and their relative post-fit errors, $\Delta \hat{\theta }/\Delta \theta $, refer to the lower scale. For experimental uncertainties that are decomposed into several independent sources, NP X corresponds to the X$^{th}$ nuisance parameter, ordered by their impact on $\mu $. The `ljets' (`dilep') label refers to the single-lepton (dilepton) channel.

Ranking of the 20 nuisance parameters with the largest post-fit impact on $\mu $ in the STXS fit for $120\le {\hat{p}_{{T}}^{H}}<200$ GeV. Nuisance parameters corresponding to statistical uncertainties in the simulated event samples are not included. The empty blue rectangles correspond to the pre-fit impact on $\mu $ and the filled blue ones to the post-fit impact on $\mu $, both referring to the upper scale. The impact of each nuisance parameter, $\Delta \mu $, is computed by comparing the nominal best-fit value of $\mu $ with the result of the fit when fixing the considered nuisance parameter to its best-fit value, $\hat{\theta }$, shifted by its pre-fit (post-fit) uncertainties $\pm \Delta \theta $ ($\pm \Delta \hat{\theta }$). The black points show the pulls of the nuisance parameters relative to their nominal values, $\theta _0$. These pulls and their relative post-fit errors, $\Delta \hat{\theta }/\Delta \theta $, refer to the lower scale. For experimental uncertainties that are decomposed into several independent sources, NP X corresponds to the X$^{th}$ nuisance parameter, ordered by their impact on $\mu $. The `ljets' (`dilep') label refers to the single-lepton (dilepton) channel.

Ranking of the 20 nuisance parameters with the largest post-fit impact on $\mu $ in the STXS fit for $200\le {\hat{p}_{{T}}^{H}}<300$ GeV. Nuisance parameters corresponding to statistical uncertainties in the simulated event samples are not included. The empty blue rectangles correspond to the pre-fit impact on $\mu $ and the filled blue ones to the post-fit impact on $\mu $, both referring to the upper scale. The impact of each nuisance parameter, $\Delta \mu $, is computed by comparing the nominal best-fit value of $\mu $ with the result of the fit when fixing the considered nuisance parameter to its best-fit value, $\hat{\theta }$, shifted by its pre-fit (post-fit) uncertainties $\pm \Delta \theta $ ($\pm \Delta \hat{\theta }$). The black points show the pulls of the nuisance parameters relative to their nominal values, $\theta _0$. These pulls and their relative post-fit errors, $\Delta \hat{\theta }/\Delta \theta $, refer to the lower scale. For experimental uncertainties that are decomposed into several independent sources, NP X corresponds to the X$^{th}$ nuisance parameter, ordered by their impact on $\mu $. The `ljets' (`dilep') label refers to the single-lepton (dilepton) channel.

Ranking of the 20 nuisance parameters with the largest post-fit impact on $\mu $ in the STXS fit for $300\le {\hat{p}_{{T}}^{H}}<450$ GeV. Nuisance parameters corresponding to statistical uncertainties in the simulated event samples are not included. The empty blue rectangles correspond to the pre-fit impact on $\mu $ and the filled blue ones to the post-fit impact on $\mu $, both referring to the upper scale. The impact of each nuisance parameter, $\Delta \mu $, is computed by comparing the nominal best-fit value of $\mu $ with the result of the fit when fixing the considered nuisance parameter to its best-fit value, $\hat{\theta }$, shifted by its pre-fit (post-fit) uncertainties $\pm \Delta \theta $ ($\pm \Delta \hat{\theta }$). The black points show the pulls of the nuisance parameters relative to their nominal values, $\theta _0$. These pulls and their relative post-fit errors, $\Delta \hat{\theta }/\Delta \theta $, refer to the lower scale. For experimental uncertainties that are decomposed into several independent sources, NP X corresponds to the X$^{th}$ nuisance parameter, ordered by their impact on $\mu $. The `ljets' (`dilep') label refers to the single-lepton (dilepton) channel.

Ranking of the 20 nuisance parameters with the largest post-fit impact on $\mu $ in the STXS fit for ${\hat{p}_{{T}}^{H}}\ge 450$ GeV. Nuisance parameters corresponding to statistical uncertainties in the simulated event samples are not included. The empty blue rectangles correspond to the pre-fit impact on $\mu $ and the filled blue ones to the post-fit impact on $\mu $, both referring to the upper scale. The impact of each nuisance parameter, $\Delta \mu $, is computed by comparing the nominal best-fit value of $\mu $ with the result of the fit when fixing the considered nuisance parameter to its best-fit value, $\hat{\theta }$, shifted by its pre-fit (post-fit) uncertainties $\pm \Delta \theta $ ($\pm \Delta \hat{\theta }$). The black points show the pulls of the nuisance parameters relative to their nominal values, $\theta _0$. These pulls and their relative post-fit errors, $\Delta \hat{\theta }/\Delta \theta $, refer to the lower scale. For experimental uncertainties that are decomposed into several independent sources, NP X corresponds to the X$^{th}$ nuisance parameter, ordered by their impact on $\mu $. The `ljets' (`dilep') label refers to the single-lepton (dilepton) channel.

95% confidence level upper limits on signal-strength measurements in the individual STXS ${\hat{p}_{{T}}^{H}}$ bins, as well as the inclusive signal-strength limit, after the fit used to extract multiple signal-strength parameters. The observed limits are shown (solid black lines), together with the expected limits both in the background-only hypothesis (dotted black lines) and in the SM hypothesis (dotted red lines). In the case of the expected limits in the background-only hypothesis, one- and two-standard-deviation uncertainty bands are also shown.

Post-fit correlation matrix (in percentages) between the $\mu $ values obtained in the STXS bins.

Performance of the Higgs boson reconstruction algorithms. For each row of `truth' ${\hat{p}_{{T}}^{H}}$, the matrix shows (in percentages) the fraction of Higgs boson candidates which are truth-matched to ${b\bar {b}}$ decays, with reconstructed $p_T^H$ in the various bins of the dilepton (left), single lepton resolved (middle) and boosted (right) channels.

Pre-fit event yields in the dilepton signal regions and control regions. All uncertainties are included except the $k({t\bar {t}+{\geq }1b})$ uncertainty that is not defined pre-fit. For the ${t\bar {t}H}$ signal, the pre-fit yield values correspond to the theoretical prediction and corresponding uncertainties. `Other sources' refers to s-channel, t-channel, $tW$, $tWZ$, $tZq$, $Z+$ jets and diboson events.

Post-fit event yields in the dilepton signal regions and control regions, after the inclusive fit in all channels. All uncertainties are included, taking into account correlations. For the ${t\bar {t}H}$ signal, the post-fit yield and uncertainties correspond to those in the inclusive signal-strength measurement. `Other sources' refers to s-channel, t-channel, $tW$, $tWZ$, $tZq$, $Z+$ jets and diboson events.

Pre-fit event yields in the single-lepton resolved and boosted signal regions and control regions. All uncertainties are included except the $k({t\bar {t}+{\geq }1b})$ uncertainty that is not defined pre-fit. For the ${t\bar {t}H}$ signal, the pre-fit yield values correspond to the theoretical prediction and corresponding uncertainties. `Other top sources' refers to s-channel, t-channel, $tWZ$ and $tZq$ events.

Post-fit event yields in the single-lepton resolved and boosted signal regions and control regions, after the inclusive fit in all channels. All uncertainties are included, taking into account correlations. For the ${t\bar {t}H}$ signal, the post-fit yield and uncertainties correspond to those in the inclusive signal-strength measurement. `Other top sources' refers to s-channel, t-channel, $tWZ$ and $tZq$ events.

Breakdown of the contributions to the uncertainties in $\mu$. The contributions from the different sources of uncertainty are evaluated after the fit. The $\Delta \mu $ values are obtained by repeating the fit after having fixed a certain set of nuisance parameters corresponding to a group of systematic uncertainties, and then evaluating $(\Delta \mu)^2$ by subtracting the resulting squared uncertainty of $\mu $ from its squared uncertainty found in the full fit. The same procedure is followed when quoting the effect of the ${t\bar {t}+{\geq }1b}$ normalisation. The total uncertainty is different from the sum in quadrature of the different components due to correlations between nuisance parameters existing in the fit.

Fraction (in percentages) of signal events, after SR and CR selections, originating from $b\bar {b}$, $WW$ and other remaining Higgs boson decay modes in the dilepton channel.

Fraction (in percentages) of signal events, after SR and CR selections, originating from $b\bar {b}$, $WW$ and other remaining Higgs boson decay modes in the single-lepton channels.

Predicted SM ${t\bar {t}H}$ cross-section in each of the five STXS ${\hat{p}_{{T}}^{H}}$ bins and signal acceptance times efficiency (including all event selection criteria) in each STXS bin as well as for the inclusive ${\hat{p}_{{T}}^{H}}$ range.

Number of expected signal events before the fit, after each selection requirement applied to enter the dilepton channel $SR^{\geq 4j}_{\geq 4b}$ region. All ${t\bar {t}H}$ signal events are included, regardless of the $H$ or ${t\bar {t}H}$ decay mode. All object corrections are applied, except for the initial number of events which is calculated using the NLO QCD+EW theoretical prediction. All quoted numbers are rounded to unity. More details on the selection criteria can be found in the text.

Number of expected signal events before the fit, after each selection requirement applied to enter the single-lepton channel resolved $SR^{\geq 6j}_{\geq 4b}$ region. All ${t\bar {t}H}$ signal events are included, regardless of the $H$ or ${t\bar {t}H}$ decay mode. All object corrections are applied, except for the initial number of events which is calculated using the NLO QCD+EW theoretical prediction. All quoted numbers are rounded to unity. More details on the selection criteria can be found in the text.

Number of expected signal events before the fit, after each selection requirement applied to enter the single-lepton channel boosted $SR_{boosted}$ region. All ${t\bar {t}H}$ signal events are included, regardless of the $H$ or ${t\bar {t}H}$ decay mode. All object corrections are applied, except for the initial number of events which is calculated using the NLO QCD+EW theoretical prediction. All quoted numbers are rounded to unity. More details on the selection criteria can be found in the text.