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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

Theory cross-section for Zprime-2HDM model.

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

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

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

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

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

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

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

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

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

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

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

Yields in 1-lepton control region.

Yields in 2-lepton control region.

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

Correlation coefficients $\rho_{i,j}$ for the statistical and systematic uncertainties between the $i$-th and $j$-th bin of the differential $A_E$ measurement as a function of the jet scattering angle $\theta_j$

The effects on the energy asymmetry of $1\sigma$ variations in its influencing nuisance parameters for the three $\theta_j$ bins. These are extracted from the samples of the posterior distribution with $\sigma_i^{(j)} = c_{ij}/\sqrt{c_{jj}}$ being the estimated shift of bin $i$ in conjunction with a shift $\Delta\theta_j$ of nuisance parameter $j$. The data statistical (Data stat.) uncertainty is obtained from running the unfolding with all nuisance parameters being fixed to their post-marginalised values, the MC statistical uncertainty on the response matrix ($t\bar{t}$ response MC stat.) is evaluated using a bootstrapping method from the covariance matrix of the ensemble of repeated unfolding results with varied response matrices. The $\gamma$ variations denote the statistical uncertainties of the background predictions in the corresponding bin of the $\Delta E$ vs $\theta_{j}$ distribution. The numbers appended to the $W$+jets PDF variations denote the corresponding NNPDF3.0 PDF sets.

The effects on the energy asymmetry of $1\sigma$ variations in its influencing nuisance parameters for the three $\theta_j$ bins. These are extracted from the samples of the posterior distribution with $\sigma_i^{(j)} = c_{ij}/\sqrt{c_{jj}}$ being the estimated shift of bin $i$ in conjunction with a shift $\Delta\theta_j$ of nuisance parameter $j$. The data statistical (Data stat.) uncertainty is obtained from running the unfolding with all nuisance parameters being fixed to their post-marginalised values, the MC statistical uncertainty on the response matrix ($t\bar{t}$ response MC stat.) is evaluated using a bootstrapping method from the covariance matrix of the ensemble of repeated unfolding results with varied response matrices. The $\gamma$ variations denote the statistical uncertainties of the background predictions in the corresponding bin of the $\Delta E$ vs $\theta_{j}$ distribution. The numbers appended to the $W$+jets PDF variations denote the corresponding NNPDF3.0 PDF sets.

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

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

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.

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.

The BDT distribution of the di-Higgs ggF signal for two different values of $\kappa_{\lambda}$ and the main backgrounds in the low 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.881 in the low mass region are discarded.

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 low 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.881 in the low mass region are discarded. The range of BDT scores is from 0.8 to 1.

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 low 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.881 in the low mass region are discarded. The range of BDT scores is from 0.8 to 1.

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 low 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.881 in the low mass region are discarded. The range of BDT scores is from 0.8 to 1.

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 low 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.881 in the low mass region are discarded. The range of BDT scores is from 0.8 to 1.

The BDT distribution 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 BDT distribution 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 BDT distribution 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 BDT distribution 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 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 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 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 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}$ = 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}$ = 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}$ = 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.

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.

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.

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

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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 nonresonant di-Higgs ggF signal sample, yields are normalized to 139 $fb^{-1}$.

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

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}$ = 300 GeV, 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}$ = 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}$.

Cutflow for resonant signal sample, with $m_{X}$ = 500 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}$.

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.

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.

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.

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.

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.

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.

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

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.

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.

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

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.

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.

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

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.

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.

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.

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.

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.

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.

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.

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.

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.

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

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

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

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

Lower mass requirement for signal regions.

Lower mass requirement for signal regions.

Lower mass requirement for signal regions.

Expected and observed events in the 16 discovery regions along with the according control regions.

Expected and observed events in the 16 discovery regions along with the according control regions.

Expected signal yield and acceptance x efficiency, estimated background and observed number of events in data for the full range of simulated masses in the MS-agnostic R-hadron search.

Expected signal yield and acceptance x efficiency, estimated background and observed number of events in data for the full range of simulated masses in the MS-agnostic R-hadron search.

Expected signal yield and acceptance x efficiency, estimated background and observed number of events in data for the full range of simulated masses in the full-detector R-hadron search.

Expected signal yield and acceptance x efficiency, estimated background and observed number of events in data for the full range of simulated masses in the full-detector R-hadron search.

p0-values and model-independent upper limits on cross-section x acceptance x efficiency for the 16 discovery regions.

p0-values and model-independent upper limits on cross-section x acceptance x efficiency for the 16 discovery regions.

Expected signal yield and acceptance x efficiency, estimated background and observed number of events in data for the full range of simulated masses in the MS-agnostic search for metastable gluino R-hadrons.

Expected signal yield and acceptance x efficiency, estimated background and observed number of events in data for the full range of simulated masses in the MS-agnostic search for metastable gluino R-hadrons.

Expected signal yield and acceptance x efficiency, estimated background and observed number of events in data for the full range of simulated masses in the full-detector direct-stau search.

Expected signal yield and acceptance x efficiency, estimated background and observed number of events in data for the full range of simulated masses in the full-detector direct-stau search.

Expected signal yield and acceptance x efficiency, estimated background and observed number of events in data for the full range of simulated masses in the full-detector chargino search.

Expected signal yield and acceptance x efficiency, estimated background and observed number of events in data for the full range of simulated masses in the full-detector chargino search.

Upper cross-section limit in gluino R-hadron search.

Upper cross-section limit in gluino R-hadron search.

Upper cross-section limit in sbottom R-hadron search.

Upper cross-section limit in sbottom R-hadron search.

Upper cross-section limit in stop R-hadron search.

Upper cross-section limit in stop R-hadron search.

Upper cross-section limit in stau search.

Upper cross-section limit in stau search.

Upper cross-section limit in chargino search.

Upper cross-section limit in chargino search.

Lower mass limit as function of gluino lifetime.

Lower mass limit as function of gluino lifetime.

Acceptance x efficiency, acceptance and efficiency for the full range of simulated masses in the MS-agnostic R-hadron search.

Acceptance x efficiency, acceptance and efficiency for the full range of simulated masses in the MS-agnostic R-hadron search.

Upper cross-section limit in meta-stable gluino R-hadron search.

Upper cross-section limit in meta-stable gluino R-hadron search.

Flavor composition of 800 GeV stop R-hadrons simulated using the generic model as a function of radial distance from the interaction point.

Flavor composition of 800 GeV stop R-hadrons simulated using the generic model as a function of radial distance from the interaction point.

Flavor composition of 800 GeV anti-stop R-hadrons simulated using the generic model as a function of radial distance from the interaction point.

Flavor composition of 800 GeV anti-stop R-hadrons simulated using the generic model as a function of radial distance from the interaction point.

Flavor composition of 800 GeV stop R-hadrons simulated using the Regge model as a function of radial distance from the interaction point.

Flavor composition of 800 GeV stop R-hadrons simulated using the Regge model as a function of radial distance from the interaction point.

Flavor composition of 800 GeV anti-stop R-hadrons simulated using the Regge model as a function of radial distance from the interaction point.

Flavor composition of 800 GeV anti-stop R-hadrons simulated using the Regge model as a function of radial distance from the interaction point.

ETmiss trigger efficiency as function of true ETmiss (EtmissTurnOn).

ETmiss trigger efficiency as function of true ETmiss (EtmissTurnOn).

Single-muon trigger efficiency as function of $|\eta|$ and $\beta$ (SingleMuTurnOn).

Single-muon trigger efficiency as function of $|\eta|$ and $\beta$ (SingleMuTurnOn).

Candidate reconstruction efficiency for ID+Calo selection (IDCaloEff).

Candidate reconstruction efficiency for ID+Calo selection (IDCaloEff).

Candidate reconstruction efficiency for loose selection (LooseEff).

Candidate reconstruction efficiency for loose selection (LooseEff).

Efficiency for a loose candidate to be promoted to a tight candidate (TightPromotionEff).

Efficiency for a loose candidate to be promoted to a tight candidate (TightPromotionEff).

Resolution and average of reconstructed dE/dx mass for a given simulated mass for ID+calo candidates.

Resolution and average of reconstructed dE/dx mass for a given simulated mass for ID+calo candidates.

Resolution and average of reconstructed ToF mass for a given simulated mass for ID+calo candidates.

Resolution and average of reconstructed ToF mass for a given simulated mass for ID+calo candidates.

Resolution and average of reconstructed ToF mass for a given simulated mass for FullDet candidates.

Resolution and average of reconstructed ToF mass for a given simulated mass for FullDet candidates.

List of relative uncertainties in the measured cross sections of the $t\bar{t}Z$ and $t\bar{t}W$ processes from the fit, grouped in categories. All uncertainties are symmetrized. The sum in quadrature may not be equal to the total due to correlations between uncertainties introduced by the fit.

The expected and observed 68% and 95% confidence intervals, which include the value 0, for $\mathcal{C}_{i}/\Lambda^{2}$ for the EFT coefficients $\mathcal{C}_{\phi Q}^{(3)}$, $\mathcal{C}_{\phi t}$, $\mathcal{C}_{tB}$ and $\mathcal{C}_{tW}$. The intervals for $\mathcal{C}_{\phi Q}^{(3)}$ are derived setting $\mathcal{C}_{\phi Q}^{(1)}$ to zero; the measurement is sensitive to the difference $\mathcal{C}_{\phi Q}^{(3)}-\mathcal{C}_{\phi Q}^{(1)}$. All results are given in units of 1/TeV$^{2}$. Limits from fits for the EFT coefficients with only the linear term are also shown.

Expected 95% CL exclusion contour for the Gbb signal.

Observed 95% CL exclusion contour for the Gbb signal.

Expected 95% CL exclusion contour for the Gtt combination.

Observed 95% CL exclusion contour for the Gtt combination.

Acceptances for the Gbb model in SR-Gbb-A. Acceptance is evaluated at truth level, with only leptons from heavy bosons and taus considered, and no further quality or isolation criteria applied in their selection.

Acceptances for the Gbb model in SR-Gbb-B. Acceptance is evaluated at truth level, with only leptons from heavy bosons and taus considered, and no further quality or isolation criteria applied in their selection.

Acceptances for the Gbb model in SR-Gbb-C. Acceptance is evaluated at truth level, with only leptons from heavy bosons and taus considered, and no further quality or isolation criteria applied in their selection.

Acceptances for the Gtt model in SR-Gtt-0L-A. Acceptance is evaluated at truth level, with only leptons from heavy bosons and taus considered, and no further quality or isolation criteria applied in their selection.

Acceptances for the Gtt model in SR-Gtt-0L-B. Acceptance is evaluated at truth level, with only leptons from heavy bosons and taus considered, and no further quality or isolation criteria applied in their selection.

Acceptances for the Gtt model in SR-Gtt-0L-C. Acceptance is evaluated at truth level, with only leptons from heavy bosons and taus considered, and no further quality or isolation criteria applied in their selection.

Acceptances for the Gtt model in SR-Gtt-1L-A. Acceptance is evaluated at truth level, with only leptons from heavy bosons and taus considered, and no further quality or isolation criteria applied in their selection.

Acceptances for the Gtt model in SR-Gtt-1L-B. Acceptance is evaluated at truth level, with only leptons from heavy bosons and taus considered, and no further quality or isolation criteria applied in their selection.

Acceptance times efficiency for the Gbb model in SR-Gbb-A.

Acceptance times efficiency for the Gbb model in SR-Gbb-B.

Acceptance times efficiency for the Gbb model in SR-Gbb-C.

Acceptance times efficiency for the Gtt model in SR-Gtt-0L-A.

Acceptance times efficiency for the Gtt model in SR-Gtt-0L-B.

Acceptance times efficiency for the Gtt model in SR-Gtt-0L-C.

Acceptance times efficiency for the Gtt model in SR-Gtt-1L-A.

Acceptance times efficiency for the Gtt model in SR-Gtt-1L-B.

95% CL upper limit on the cross-section times branching ratio (in fb) for the Gbb model in SR-Gbb-A.

95% CL upper limit on the cross-section times branching ratio (in fb) for the Gbb model in SR-Gbb-B.

95% CL upper limit on the cross-section times branching ratio (in fb) for the Gbb model in SR-Gbb-C.

95% CL upper limit on the cross-section times branching ratio (in fb) for the Gtt model in SR-Gtt-0L-A.

95% CL upper limit on the cross-section times branching ratio (in fb) for the Gtt model in SR-Gtt-0L-B.

95% CL upper limit on the cross-section times branching ratio (in fb) for the Gtt model in SR-Gtt-0L-C.

95% CL upper limit on the cross-section times branching ratio (in fb) for the Gtt model in SR-Gtt-1L-A.

95% CL upper limit on the cross-section times branching ratio (in fb) for the Gtt model in SR-Gtt-1L-B.

Signal region yielding the best expected sensitivity for each point of the parameter space in the Gbb model.

Signal region yielding the best expected sensitivity for each point of the parameter space in the Gtt model for the 0-lepton channel.

Signal region yielding the best expected sensitivity for each point of the parameter space in the Gtt model for the 1-lepton channel.

Combination of two 0-lepton and 1-lepton signal regions yielding the best expected sensitivity for each point of the parameter space in the Gtt model.

$E_{\mathrm{T}}^{\mathrm{miss}} / \sqrt{H_{\mathrm{T}}}$ distribution in signal region 10j50 2b. Two benchmark signal models are overlaid on the plot for comparison. Labelled `pMSSM' and `2-step', they show signal distributions from the example SUSY models (as described in the paper): a pMSSM slice model with ($m \tilde{g}$, $m \tilde{\chi_{1}^{\pm}}$) = (1300, 200) GeV and a cascade decay model with ($m \tilde{g}$, $m \tilde{\chi_{1}^{0}}$) = (1300, 200) GeV.

$E_{\mathrm{T}}^{\mathrm{miss}} / \sqrt{H_{\mathrm{T}}}$ distribution in signal region 8j80 0b. Two benchmark signal models are overlaid on the plot for comparison. Labelled `pMSSM' and `2-step', they show signal distributions from the example SUSY models (as described in the paper): a pMSSM slice model with ($m \tilde{g}$, $m \tilde{\chi_{1}^{\pm}}$) = (1300, 200) GeV and a cascade decay model with ($m \tilde{g}$, $m \tilde{\chi_{1}^{0}}$) = (1300, 200) GeV.

$E_{\mathrm{T}}^{\mathrm{miss}} / \sqrt{H_{\mathrm{T}}}$ distribution in signal region 8j80 2b. Two benchmark signal models are overlaid on the plot for comparison. Labelled `pMSSM' and `2-step', they show signal distributions from the example SUSY models (as described in the paper): a pMSSM slice model with ($m \tilde{g}$, $m \tilde{\chi_{1}^{\pm}}$) = (1300, 200) GeV and a cascade decay model with ($m \tilde{g}$, $m \tilde{\chi_{1}^{0}}$) = (1300, 200) GeV.

Observed 95% CL limit for the pMSSM grid.

Observed 95% CL limit for the pMSSM grid when the signal cross section is increased by one standard deviation.

Observed 95% CL limit for the pMSSM grid when the signal cross section is decreased by one standard deviation.

Expected 95% CL limit for the pMSSM grid.

+1 sigma excursion of the expected 95% CL limit for the pMSSM grid.

-1 sigma excursion of the expected 95% CL limit for the pMSSM grid.

Observed 95% CL limit for the 2Step grid.

Observed 95% CL limit for the 2Step grid when the signal cross section is increased by one standard deviation.

Observed 95% CL limit for the 2Step grid when the signal cross section is decreased by one standard deviation.

Expected 95% CL limit for the 2Step grid.

+1 sigma excursion of the expected 95% CL limit for the 2Step grid.

-1 sigma excursion of the expected 95% CL limit for the 2Step grid.

$E_{\mathrm{T}}^{\mathrm{miss}} / \sqrt{H_{\mathrm{T}}}$ distribution in signal region 8j50 0b. Two benchmark signal models are overlaid on the plot for comparison. Labelled `pMSSM' and `2-step', they show signal distributions from the example SUSY models (as described in the paper): a pMSSM slice model with ($m \tilde{g}$, $m \tilde{\chi_{1}^{\pm}}$) = (1300, 200) GeV and a cascade decay model with ($m \tilde{g}$, $m \tilde{\chi_{1}^{0}}$) = (1300, 200) GeV.

$E_{\mathrm{T}}^{\mathrm{miss}} / \sqrt{H_{\mathrm{T}}}$ distribution in signal region 8j50 1b. Two benchmark signal models are overlaid on the plot for comparison. Labelled `pMSSM' and `2-step', they show signal distributions from the example SUSY models (as described in the paper): a pMSSM slice model with ($m \tilde{g}$, $m \tilde{\chi_{1}^{\pm}}$) = (1300, 200) GeV and a cascade decay model with ($m \tilde{g}$, $m \tilde{\chi_{1}^{0}}$) = (1300, 200) GeV.

$E_{\mathrm{T}}^{\mathrm{miss}} / \sqrt{H_{\mathrm{T}}}$ distribution in signal region 8j50 2b. Two benchmark signal models are overlaid on the plot for comparison. Labelled `pMSSM' and `2-step', they show signal distributions from the example SUSY models (as described in the paper): a pMSSM slice model with ($m \tilde{g}$, $m \tilde{\chi_{1}^{\pm}}$) = (1300, 200) GeV and a cascade decay model with ($m \tilde{g}$, $m \tilde{\chi_{1}^{0}}$) = (1300, 200) GeV.

$E_{\mathrm{T}}^{\mathrm{miss}} / \sqrt{H_{\mathrm{T}}}$ distribution in signal region 9j50 0b. Two benchmark signal models are overlaid on the plot for comparison. Labelled `pMSSM' and `2-step', they show signal distributions from the example SUSY models (as described in the paper): a pMSSM slice model with ($m \tilde{g}$, $m \tilde{\chi_{1}^{\pm}}$) = (1300, 200) GeV and a cascade decay model with ($m \tilde{g}$, $m \tilde{\chi_{1}^{0}}$) = (1300, 200) GeV.

$E_{\mathrm{T}}^{\mathrm{miss}} / \sqrt{H_{\mathrm{T}}}$ distribution in signal region 9j50 1b. Two benchmark signal models are overlaid on the plot for comparison. Labelled `pMSSM' and `2-step', they show signal distributions from the example SUSY models (as described in the paper): a pMSSM slice model with ($m \tilde{g}$, $m \tilde{\chi_{1}^{\pm}}$) = (1300, 200) GeV and a cascade decay model with ($m \tilde{g}$, $m \tilde{\chi_{1}^{0}}$) = (1300, 200) GeV.

$E_{\mathrm{T}}^{\mathrm{miss}} / \sqrt{H_{\mathrm{T}}}$ distribution in signal region 9j50 2b. Two benchmark signal models are overlaid on the plot for comparison. Labelled `pMSSM' and `2-step', they show signal distributions from the example SUSY models (as described in the paper): a pMSSM slice model with ($m \tilde{g}$, $m \tilde{\chi_{1}^{\pm}}$) = (1300, 200) GeV and a cascade decay model with ($m \tilde{g}$, $m \tilde{\chi_{1}^{0}}$) = (1300, 200) GeV.

$E_{\mathrm{T}}^{\mathrm{miss}} / \sqrt{H_{\mathrm{T}}}$ distribution in signal region 10j50 0b. Two benchmark signal models are overlaid on the plot for comparison. Labelled `pMSSM' and `2-step', they show signal distributions from the example SUSY models (as described in the paper): a pMSSM slice model with ($m \tilde{g}$, $m \tilde{\chi_{1}^{\pm}}$) = (1300, 200) GeV and a cascade decay model with ($m \tilde{g}$, $m \tilde{\chi_{1}^{0}}$) = (1300, 200) GeV.

$E_{\mathrm{T}}^{\mathrm{miss}} / \sqrt{H_{\mathrm{T}}}$ distribution in signal region 10j50 1b. Two benchmark signal models are overlaid on the plot for comparison. Labelled `pMSSM' and `2-step', they show signal distributions from the example SUSY models (as described in the paper): a pMSSM slice model with ($m \tilde{g}$, $m \tilde{\chi_{1}^{\pm}}$) = (1300, 200) GeV and a cascade decay model with ($m \tilde{g}$, $m \tilde{\chi_{1}^{0}}$) = (1300, 200) GeV.

$E_{\mathrm{T}}^{\mathrm{miss}} / \sqrt{H_{\mathrm{T}}}$ distribution in signal region 10j50 2b. Two benchmark signal models are overlaid on the plot for comparison. Labelled `pMSSM' and `2-step', they show signal distributions from the example SUSY models (as described in the paper): a pMSSM slice model with ($m \tilde{g}$, $m \tilde{\chi_{1}^{\pm}}$) = (1300, 200) GeV and a cascade decay model with ($m \tilde{g}$, $m \tilde{\chi_{1}^{0}}$) = (1300, 200) GeV.

$E_{\mathrm{T}}^{\mathrm{miss}} / \sqrt{H_{\mathrm{T}}}$ distribution in signal region 7j80 0b. Two benchmark signal models are overlaid on the plot for comparison. Labelled `pMSSM' and `2-step', they show signal distributions from the example SUSY models (as described in the paper): a pMSSM slice model with ($m \tilde{g}$, $m \tilde{\chi_{1}^{\pm}}$) = (1300, 200) GeV and a cascade decay model with ($m \tilde{g}$, $m \tilde{\chi_{1}^{0}}$) = (1300, 200) GeV.

$E_{\mathrm{T}}^{\mathrm{miss}} / \sqrt{H_{\mathrm{T}}}$ distribution in signal region 7j80 1b. Two benchmark signal models are overlaid on the plot for comparison. Labelled `pMSSM' and `2-step', they show signal distributions from the example SUSY models (as described in the paper): a pMSSM slice model with ($m \tilde{g}$, $m \tilde{\chi_{1}^{\pm}}$) = (1300, 200) GeV and a cascade decay model with ($m \tilde{g}$, $m \tilde{\chi_{1}^{0}}$) = (1300, 200) GeV.

$E_{\mathrm{T}}^{\mathrm{miss}} / \sqrt{H_{\mathrm{T}}}$ distribution in signal region 7j80 2b. Two benchmark signal models are overlaid on the plot for comparison. Labelled `pMSSM' and `2-step', they show signal distributions from the example SUSY models (as described in the paper): a pMSSM slice model with ($m \tilde{g}$, $m \tilde{\chi_{1}^{\pm}}$) = (1300, 200) GeV and a cascade decay model with ($m \tilde{g}$, $m \tilde{\chi_{1}^{0}}$) = (1300, 200) GeV.

$E_{\mathrm{T}}^{\mathrm{miss}} / \sqrt{H_{\mathrm{T}}}$ distribution in signal region 8j80 0b. Two benchmark signal models are overlaid on the plot for comparison. Labelled `pMSSM' and `2-step', they show signal distributions from the example SUSY models (as described in the paper): a pMSSM slice model with ($m \tilde{g}$, $m \tilde{\chi_{1}^{\pm}}$) = (1300, 200) GeV and a cascade decay model with ($m \tilde{g}$, $m \tilde{\chi_{1}^{0}}$) = (1300, 200) GeV.

$E_{\mathrm{T}}^{\mathrm{miss}} / \sqrt{H_{\mathrm{T}}}$ distribution in signal region 8j80 1b. Two benchmark signal models are overlaid on the plot for comparison. Labelled `pMSSM' and `2-step', they show signal distributions from the example SUSY models (as described in the paper): a pMSSM slice model with ($m \tilde{g}$, $m \tilde{\chi_{1}^{\pm}}$) = (1300, 200) GeV and a cascade decay model with ($m \tilde{g}$, $m \tilde{\chi_{1}^{0}}$) = (1300, 200) GeV.

$E_{\mathrm{T}}^{\mathrm{miss}} / \sqrt{H_{\mathrm{T}}}$ distribution in signal region 8j80 2b. Two benchmark signal models are overlaid on the plot for comparison. Labelled `pMSSM' and `2-step', they show signal distributions from the example SUSY models (as described in the paper): a pMSSM slice model with ($m \tilde{g}$, $m \tilde{\chi_{1}^{\pm}}$) = (1300, 200) GeV and a cascade decay model with ($m \tilde{g}$, $m \tilde{\chi_{1}^{0}}$) = (1300, 200) GeV.

Degree of multijet closure for signal and vaidation regions with at no b-jet requirement. The solid lines are the pre-fit predicted numbers of events and the points are the observed numbers. The blue hatched band shows only the statistical (MC and data) uncertainty on the background estimate. The bins labelled in bold are signal regions, while the others are validation regions. The template closure uncertainty for each SR bin is given by the maximal deviation of data from prediction in any non-SR bin to its left on this plot (although those for 80 GeV regions are independent of deviations in 50 GeV regions).

Degree of multijet closure for signal and vaidation regions with at least 1 b-jet. The solid lines are the pre-fit predicted numbers of events and the points are the observed numbers. The blue hatched band shows only the statistical (MC and data) uncertainty on the background estimate. The bins labelled in bold are signal regions, while the others are validation regions. The template closure uncertainty for each SR bin is given by the maximal deviation of data from prediction in any non-SR bin to its left on this plot (although those for 80 GeV regions are independent of deviations in 50 GeV regions).

Degree of multijet closure for signal and vaidation regions with at least 2 b-jets. The solid lines are the pre-fit predicted numbers of events and the points are the observed numbers. The blue hatched band shows only the statistical (MC and data) uncertainty on the background estimate. The bins labelled in bold are signal regions, while the others are validation regions. The template closure uncertainty for each SR bin is given by the maximal deviation of data from prediction in any non-SR bin to its left on this plot (although those for 80 GeV regions are independent of deviations in 50 GeV regions).

Summary of all 15 signal regions (post-fit).

Signal region yielding the best-expected CLs value, the best expected CLs value, and the corresponding observed CLs value for the 2Step grid.

Signal region yielding the best-expected CLs value, the best expected CLs value, and the corresponding observed CLs value for the pMSSM grid.

95% CLs observed upper limit on model cross-section for 2-step signal points for the best-expected signal region.

Performance of the 8j50-0b selection for the pMSSM grid: number of generated signal events; total signal cross-section; acceptance; efficiency (fractional); observed CL using this region alone; expected CL using this region alone.

Performance of the 8j50-1b selection for the pMSSM grid: number of generated signal events; total signal cross-section; acceptance; efficiency (fractional); observed CL using this region alone; expected CL using this region alone.

Performance of the 8j50-2b selection for the pMSSM grid: number of generated signal events; total signal cross-section; acceptance; efficiency (fractional); observed CL using this region alone; expected CL using this region alone.

Performance of the 9j50-0b selection for the pMSSM grid: number of generated signal events; total signal cross-section; acceptance; efficiency (fractional); observed CL using this region alone; expected CL using this region alone.

Performance of the 9j50-1b selection for the pMSSM grid: number of generated signal events; total signal cross-section; acceptance; efficiency (fractional); observed CL using this region alone; expected CL using this region alone.

Performance of the 9j50-2b selection for the pMSSM grid: number of generated signal events; total signal cross-section; acceptance; efficiency (fractional); observed CL using this region alone; expected CL using this region alone.

Performance of the 10j50-0b selection for the pMSSM grid: number of generated signal events; total signal cross-section; acceptance; efficiency (fractional); observed CL using this region alone; expected CL using this region alone.

Performance of the 10j50-1b selection for the pMSSM grid: number of generated signal events; total signal cross-section; acceptance; efficiency (fractional); observed CL using this region alone; expected CL using this region alone.

Performance of the 10j50-2b selection for the pMSSM grid: number of generated signal events; total signal cross-section; acceptance; efficiency (fractional); observed CL using this region alone; expected CL using this region alone.

Performance of the 7j80-0b selection for the pMSSM grid: number of generated signal events; total signal cross-section; acceptance; efficiency (fractional); observed CL using this region alone; expected CL using this region alone.

Performance of the 7j80-1b selection for the pMSSM grid: number of generated signal events; total signal cross-section; acceptance; efficiency (fractional); observed CL using this region alone; expected CL using this region alone.

Performance of the 7j80-2b selection for the pMSSM grid: number of generated signal events; total signal cross-section; acceptance; efficiency (fractional); observed CL using this region alone; expected CL using this region alone.

Performance of the 8j80-0b selection for the pMSSM grid: number of generated signal events; total signal cross-section; acceptance; efficiency (fractional); observed CL using this region alone; expected CL using this region alone.

Performance of the 8j80-1b selection for the pMSSM grid: number of generated signal events; total signal cross-section; acceptance; efficiency (fractional); observed CL using this region alone; expected CL using this region alone.

Performance of the 8j80-2b selection for the pMSSM grid: number of generated signal events; total signal cross-section; acceptance; efficiency (fractional); observed CL using this region alone; expected CL using this region alone.

Performance of the 8j50-0b selection for the 2Step grid: number of generated signal events; total signal cross-section; acceptance; efficiency (fractional); observed CL using this region alone; expected CL using this region alone.

Performance of the 8j50-1b selection for the 2Step grid: number of generated signal events; total signal cross-section; acceptance; efficiency (fractional); observed CL using this region alone; expected CL using this region alone.

Performance of the 8j50-2b selection for the 2Step grid: number of generated signal events; total signal cross-section; acceptance; efficiency (fractional); observed CL using this region alone; expected CL using this region alone.

Performance of the 9j50-0b selection for the 2Step grid: number of generated signal events; total signal cross-section; acceptance; efficiency (fractional); observed CL using this region alone; expected CL using this region alone.

Performance of the 9j50-1b selection for the 2Step grid: number of generated signal events; total signal cross-section; acceptance; efficiency (fractional); observed CL using this region alone; expected CL using this region alone.

Performance of the 9j50-2b selection for the 2Step grid: number of generated signal events; total signal cross-section; acceptance; efficiency (fractional); observed CL using this region alone; expected CL using this region alone.

Performance of the 10j50-0b selection for the 2Step grid: number of generated signal events; total signal cross-section; acceptance; efficiency (fractional); observed CL using this region alone; expected CL using this region alone.

Performance of the 10j50-1b selection for the 2Step grid: number of generated signal events; total signal cross-section; acceptance; efficiency (fractional); observed CL using this region alone; expected CL using this region alone.

Performance of the 10j50-2b selection for the 2Step grid: number of generated signal events; total signal cross-section; acceptance; efficiency (fractional); observed CL using this region alone; expected CL using this region alone.

Performance of the 7j80-0b selection for the 2Step grid: number of generated signal events; total signal cross-section; acceptance; efficiency (fractional); observed CL using this region alone; expected CL using this region alone.

Performance of the 7j80-1b selection for the 2Step grid: number of generated signal events; total signal cross-section; acceptance; efficiency (fractional); observed CL using this region alone; expected CL using this region alone.

Performance of the 7j80-2b selection for the 2Step grid: number of generated signal events; total signal cross-section; acceptance; efficiency (fractional); observed CL using this region alone; expected CL using this region alone.

Performance of the 8j80-0b selection for the 2Step grid: number of generated signal events; total signal cross-section; acceptance; efficiency (fractional); observed CL using this region alone; expected CL using this region alone.

Performance of the 8j80-1b selection for the 2Step grid: number of generated signal events; total signal cross-section; acceptance; efficiency (fractional); observed CL using this region alone; expected CL using this region alone.

Performance of the 8j80-2b selection for the 2Step grid: number of generated signal events; total signal cross-section; acceptance; efficiency (fractional); observed CL using this region alone; expected CL using this region alone.

Observed and expected background and signal effective mass distributions for SR2jl. For signal, a squark direct decay model with $m(\tilde q)=800$ GeV and $m(\tilde\chi^0_1)=400$ GeV is shown.

Observed and expected background and signal effective mass distributions for SR2jm. For signal, a gluino direct decay model with $m(\tilde g)=750$ GeV and $m(\tilde\chi^0_1)=660$ GeV is shown.

Observed and expected background and signal effective mass distributions for SR2jt. For signal, a squark direct decay model with $m(\tilde q)=1200$ GeV and $m(\tilde\chi^0_1)=0$ GeV is shown.

Observed and expected background and signal effective mass distributions for SR4jt. For signal, a gluino direct decay model with $m(\tilde g)=1400$ GeV and $m(\tilde\chi^0_1)=0$ GeV is shown.

Observed and expected background and signal effective mass distributions for SR5j. For signal, a gluino one-step decay model with $m(\tilde g)=1265$ GeV, $m(\tilde\chi^\pm_1)=945$ GeV and $m(\tilde\chi^0_1)=625$ GeV is shown.

Observed and expected background and signal effective mass distributions for SR6jm. For signal, a gluino one-step decay model with $m(\tilde g)=1265$ GeV, $m(\tilde\chi^\pm_1)=945$ GeV and $m(\tilde\chi^0_1)=625$ GeV is shown.

Observed and expected background and signal effective mass distributions for SR6jt. For signal, a gluino one-step decay model with $m(\tilde g)=1385$ GeV, $m(\tilde\chi^\pm_1)=705$ GeV and $m(\tilde\chi^0_1)=25$ GeV is shown.

Expected limit at 95% CL for squark direct decay model grid.

Expected limits at 95% CL +1 sigma excursion due to experimental and background-only theoretical uncertainties for squark direct decay model grid.

Expected limits at 95% CL -1 sigma excursion due to experimental and background-only theoretical uncertainties for squark direct decay model grid.

Observed limits at 95% CL for squark direct decay model grid.

Observed limits at 95% CL +1 sigma excursion due to the signal cross-section uncertainty for squark direct decay model grid.

Observed limits at 95% CL -1 sigma excursion due to the signal cross-section uncertainty for squark direct decay model grid.

Expected limit at 95% CL for gluino direct decay model grid.

Expected limits at 95% CL +1 sigma excursion due to experimental and background-only theoretical uncertainties for gluino direct decay model grid.

Expected limits at 95% CL -1 sigma excursion due to experimental and background-only theoretical uncertainties for gluino direct decay model grid.

Observed limits at 95% CL for gluino direct decay model grid.

Observed limits at 95% CL +1 sigma excursion due to the signal cross-section uncertainty for gluino direct decay model grid.

Observed limits at 95% CL -1 sigma excursion due to the signal cross-section uncertainty for gluino direct decay model grid.

Expected limit at 95% CL for gluino one-step decay model grid.

Expected limits at 95% CL +1 sigma excursion due to experimental and background-only theoretical uncertainties for gluino one-step decay model grid.

Expected limits at 95% CL -1 sigma excursion due to experimental and background-only theoretical uncertainties for gluino one-step decay model grid.

Observed limits at 95% CL for gluino one-step decay model grid.

Observed limits at 95% CL +1 sigma excursion due to the signal cross-section uncertainty for gluino one-step decay model grid.

Observed limits at 95% CL -1 sigma excursion due to the signal cross-section uncertainty for gluino one-step decay model grid.

Observed and expected background effective mass distributions in control region CRgamma for SR2jl.

Observed and expected background effective mass distributions in validation region VRZ for SR2jl.

Observed and expected background effective mass distributions in control region CRW for SR2jl.

Observed and expected background effective mass distributions in control region CRT for SR2jl.

Observed and expected background effective mass distributions in control region CRgamma for SR2jm.

Observed and expected background effective mass distributions in validation region VRZ for SR2jm.

Observed and expected background effective mass distributions in control region CRW for SR2jm.

Observed and expected background effective mass distributions in control region CRT for SR2jm.

Observed and expected background effective mass distributions in control region CRgamma for SR2jt.

Observed and expected background effective mass distributions in validation region VRZ for SR2jt.

Observed and expected background effective mass distributions in control region CRW for SR2jt.

Observed and expected background effective mass distributions in control region CRT for SR2jt.

Observed and expected background effective mass distributions in control region CRgamma for SR4jt.

Observed and expected background effective mass distributions in validation region VRZ for SR4jt.

Observed and expected background effective mass distributions in control region CRW for SR4jt.

Observed and expected background effective mass distributions in control region CRT for SR4jt.

Observed and expected background effective mass distributions in control region CRgamma for SR5j.

Observed and expected background effective mass distributions in validation region VRZ for SR5j.

Observed and expected background effective mass distributions in control region CRW for SR5j.

Observed and expected background effective mass distributions in control region CRT for SR5j.

Observed and expected background effective mass distributions in control region CRgamma for SR6jm.

Observed and expected background effective mass distributions in validation region VRZ for SR6jm.

Observed and expected background effective mass distributions in control region CRW for SR6jm.

Observed and expected background effective mass distributions in control region CRT for SR6jm.

Observed and expected background effective mass distributions in control region CRgamma for SR6jt.

Observed and expected background effective mass distributions in validation region VRZ for SR6jt.

Observed and expected background effective mass distributions in control region CRW for SR6jt.

Observed and expected background effective mass distributions in control region CRT for SR6jt.

Observed and expected event yields in VRZ as a function of signal region.

Observed and expected event yields in VRW as a function of signal region.

Observed and expected event yields in VRWv as a function of signal region.

Observed and expected event yields in VRT as a function of signal region.

Observed and expected event yields in VRTv as a function of signal region.

Observed and expected event yields in VRQa as a function of signal region.

Observed and expected event yields in VRQb as a function of signal region.

Signal acceptance for SR2jl in squark direct decay model grid.

Signal acceptance times efficiency for SR2jl in squark direct decay model grid.

Signal acceptance for SR2jm in squark direct decay model grid.

Signal acceptance times efficiency for SR2jm in squark direct decay model grid.

Signal acceptance for SR2jt in squark direct decay model grid.

Signal acceptance times efficiency for SR2jt in squark direct decay model grid.

Signal acceptance for SR4jt in squark direct decay model grid.

Signal acceptance times efficiency for SR4jt in squark direct decay model grid.

Signal acceptance for SR5j in squark direct decay model grid.

Signal acceptance times efficiency for SR5j in squark direct decay model grid.

Signal acceptance for SR6jm in squark direct decay model grid.

Signal acceptance times efficiency for SR6jm in squark direct decay model grid.

Signal acceptance for SR6jt in squark direct decay model grid.

Signal acceptance times efficiency for SR6jt in squark direct decay model grid.

Signal acceptance for SR2jl in gluino direct decay model grid.

Signal acceptance times efficiency for SR2jl in gluino direct decay model grid.

Signal acceptance for SR2jm in gluino direct decay model grid.

Signal acceptance times efficiency for SR2jm in gluino direct decay model grid.

Signal acceptance for SR2jt in gluino direct decay model grid.

Signal acceptance times efficiency for SR2jt in gluino direct decay model grid.

Signal acceptance for SR4jt in gluino direct decay model grid.

Signal acceptance times efficiency for SR4jt in gluino direct decay model grid.

Signal acceptance for SR5j in gluino direct decay model grid.

Signal acceptance times efficiency for SR5j in gluino direct decay model grid.

Signal acceptance for SR6jm in gluino direct decay model grid.

Signal acceptance times efficiency for SR6jm in gluino direct decay model grid.

Signal acceptance for SR6jt in gluino direct decay model grid.

Signal acceptance times efficiency for SR6jt in gluino direct decay model grid.

Signal acceptance for SR2jl in gluino one-step decay model grid.

Signal acceptance times efficiency for SR2jl in gluino one-step decay model grid.

Signal acceptance for SR2jm in gluino one-step decay model grid.

Signal acceptance times efficiency for SR2jm in gluino one-step decay model grid.

Signal acceptance for SR2j5 in gluino one-step decay model grid.

Signal acceptance times efficiency for SR2jt in gluino one-step decay model grid.

Signal acceptance for SR4jt in gluino one-step decay model grid.

Signal acceptance times efficiency for SR4jt in gluino one-step decay model grid.

Signal acceptance for SR5j in gluino one-step decay model grid.

Signal acceptance times efficiency for SR5j in gluino one-step decay model grid.

Signal acceptance for SR6jm in gluino one-step decay model grid.

Signal acceptance times efficiency for SR6jm in gluino one-step decay model grid.

Signal acceptance for SR6jt in gluino one-step decay model grid.

Signal acceptance times efficiency for SR6jt in gluino one-step decay model grid.

The distribution of the missing transverse momentum is shown in soft-lepton 2-jet W+jets control regions after normalising the ttbar and W+jets background processes in the simultaneous fit.

Expected background yields as obtained in the background-only fits in all hard-lepton and soft-lepton validation together with observed data are given. Uncertainties in the fitted background estimates combine statistical (in the simulated event yields) and systematic uncertainties.

Expected background yields as obtained in the background-only fits in all hard-lepton and soft-lepton signal together with observed data are given. Uncertainties in the fitted background estimates combine statistical (in the simulated event yields) and systematic uncertainties.

Distributions of mt for the hard-lepton 4-jet low-x signal region. The requirement on the variable plotted is removed from the definitions of the signal regions, where the arrow indicates the position of the cut in the signal region. The lower panels of the plots show the ratio of the observed data to the total background prediction as derived in the background-only fit. The uncertainty bands plotted include all statistical and systematic uncertainties as discussed in Section 7. The component `Others' is the sum of Z+jets and ttbar+V. The last bin includes the overflow.

Distributions of met/meff for the 4-jet high-x signal region. The requirement on the variable plotted is removed from the definitions of the signal regions, where the arrow indicates the position of the cut in the signal region. The lower panels of the plots show the ratio of the observed data to the total background prediction as derived in the background-only fit. The uncertainty bands plotted include all statistical and systematic uncertainties as discussed in Section 7. The component `Others' is the sum of Z+jets and ttbar+V. The last bin includes the overflow.

Distributions of mt for the hard-lepton 5-jet signal region. The requirement on the variable plotted is removed from the definitions of the signal regions, where the arrow indicates the position of the cut in the signal region. The lower panels of the plots show the ratio of the observed data to the total background prediction as derived in the background-only fit. The uncertainty bands plotted include all statistical and systematic uncertainties as discussed in Section 7. The component `Others' is the sum of Z+jets and ttbar+V. The last bin includes the overflow.

Distributions of mt for the hard-lepton 6-jet signal region. The requirement on the variable plotted is removed from the definitions of the signal regions, where the arrow indicates the position of the cut in the signal region. The lower panels of the plots show the ratio of the observed data to the total background prediction as derived in the background-only fit. The uncertainty bands plotted include all statistical and systematic uncertainties as discussed in Section 7. The component `Others' is the sum of Z+jets and ttbar+V. The last bin includes the overflow.

Distributions of met for the soft-lepton 2-jet signal region. The requirement on the variable plotted is removed from the definitions of the signal regions, where the arrow indicates the position of the cut in the signal region. The lower panels of the plots show the ratio of the observed data to the total background prediction as derived in the background-only fit. The uncertainty bands plotted include all statistical and systematic uncertainties as discussed in Section 7. The component `Others' is the sum of Z+jets and ttbar+V. The last bin includes the overflow.

Distributions of met for the soft-lepton 5-jet signal region. The requirement on the variable plotted is removed from the definitions of the signal regions, where the arrow indicates the position of the cut in the signal region. The lower panels of the plots show the ratio of the observed data to the total background prediction as derived in the background-only fit. The uncertainty bands plotted include all statistical and systematic uncertainties as discussed in Section 7. The component `Others' is the sum of Z+jets and ttbar+V. The last bin includes the overflow.

The observed combined 95% CL exclusion limits in the the gluino simplified models using for each model point the signal region with the best expected sensitivity. The limits are presented in the (gluino, chargino) mass plane for the scenario where the mass of the chargino is fixed to $x=(m(\tilde\chi^\pm_1)-m(\tilde\chi^0_1))/(m(\tilde g) - m(\tilde\chi^0_1)) = 1/2$ models.

The expected combined 95% CL exclusion limits in the the gluino simplified models using for each model point the signal region with the best expected sensitivity. The limits are presented in the (gluino, chargino) mass plane for the scenario where the mass of the chargino is fixed to $x=(m(\tilde\chi^\pm_1)-m(\tilde\chi^0_1))/(m(\tilde g) - m(\tilde\chi^0_1)) = 1/2$ models.

The yellow band ($+ 1 \sigma$) of the combined 95% CL exclusion limits in the the gluino simplified models using for each model point the signal region with the best expected sensitivity. The limits are presented in the (gluino, chargino) mass plane for the scenario where the mass of the chargino is fixed to $x=(m(\tilde\chi^\pm_1)-m(\tilde\chi^0_1))/(m(\tilde g) - m(\tilde\chi^0_1)) = 1/2$ models. The yellow band represents the $\pm 1 \sigma$ variation of the median expected limit due to the experimental and theoretical uncertainties.

The yellow band ($- 1 \sigma$) of the combined 95% CL exclusion limits in the the gluino simplified models using for each model point the signal region with the best expected sensitivity. The limits are presented in the (gluino, chargino) mass plane for the scenario where the mass of the chargino is fixed to $x=(m(\tilde\chi^\pm_1)-m(\tilde\chi^0_1))/(m(\tilde g) - m(\tilde\chi^0_1)) = 1/2$ models. The yellow band represents the $\pm 1 \sigma$ variation of the median expected limit due to the experimental and theoretical uncertainties.

The observed combined 95% CL exclusion limits in the the gluino simplified models using for each model point the signal region with the best expected sensitivity. The limits are presented in the (gluino, x) plane for the chargino = 60 GeV models where $x=(m(\tilde\chi^\pm_1)-m(\tilde\chi^0_1))/(m(\tilde g) - m(\tilde\chi^0_1))$.

The expected combined 95% CL exclusion limits in the the gluino simplified models using for each model point the signal region with the best expected sensitivity. The limits are presented in the (gluino, x) plane for the chargino = 60 GeV models where $x=(m(\tilde\chi^\pm_1)-m(\tilde\chi^0_1))/(m(\tilde g) - m(\tilde\chi^0_1))$.

The yellow band ($+ 1 \sigma$) of the combined 95% CL exclusion limits in the the gluino simplified models using for each model point the signal region with the best expected sensitivity. The limits are presented in the (gluino, x) plane for the chargino = 60 GeV models where $x=(m(\tilde\chi^\pm_1)-m(\tilde\chi^0_1))/(m(\tilde g) - m(\tilde\chi^0_1))$. The yellow band represents the $\pm 1 \sigma$ variation of the median expected limit due to the experimental and theoretical uncertainties.

The yellow band ($- 1 \sigma$) of the combined 95% CL exclusion limits in the the gluino simplified models using for each model point the signal region with the best expected sensitivity. The limits are presented in the (gluino, x) plane for the chargino = 60 GeV models where $x=(m(\tilde\chi^\pm_1)-m(\tilde\chi^0_1))/(m(\tilde g) - m(\tilde\chi^0_1))$. The yellow band represents the $\pm 1 \sigma$ variation of the median expected limit due to the experimental and theoretical uncertainties.

The observed limits for the soft-lepton 2-jet signal region. The limits are presented in the (gluino, chargino) mass plane for the scenario where the mass of the chargino is fixed to $x=(m(\tilde\chi^\pm_1)-m(\tilde\chi^0_1))/(m(\tilde g) - m(\tilde\chi^0_1)) = 1/2$ models.

The expected limits for the soft-lepton 2-jet signal region. The limits are presented in the (gluino, chargino) mass plane for the scenario where the mass of the chargino is fixed to $x=(m(\tilde\chi^\pm_1)-m(\tilde\chi^0_1))/(m(\tilde g) - m(\tilde\chi^0_1)) = 1/2$ models.

The observed limits for the hard-lepton 4-jet low-x signal region. The limits are presented in the (gluino, chargino) mass plane for the scenario where the mass of the chargino is fixed to $x=(m(\tilde\chi^\pm_1)-m(\tilde\chi^0_1))/(m(\tilde g) - m(\tilde\chi^0_1)) = 1/2$ models.

The expected limits for the hard-lepton 4-jet low-x signal region. The limits are presented in the (gluino, chargino) mass plane for the scenario where the mass of the chargino is fixed to $x=(m(\tilde\chi^\pm_1)-m(\tilde\chi^0_1))/(m(\tilde g) - m(\tilde\chi^0_1)) = 1/2$ models.

The observed limits for the hard-lepton 5-jet signal region. The limits are presented in the (gluino, chargino) mass plane for the scenario where the mass of the chargino is fixed to $x=(m(\tilde\chi^\pm_1)-m(\tilde\chi^0_1))/(m(\tilde g) - m(\tilde\chi^0_1)) = 1/2$ models.

The expected limits for the hard-lepton 5-jet signal region. The limits are presented in the (gluino, chargino) mass plane for the scenario where the mass of the chargino is fixed to $x=(m(\tilde\chi^\pm_1)-m(\tilde\chi^0_1))/(m(\tilde g) - m(\tilde\chi^0_1)) = 1/2$ models.

The observed limits for the hard-lepton 6-jet signal region. The limits are presented in the (gluino, chargino) mass plane for the scenario where the mass of the chargino is fixed to $x=(m(\tilde\chi^\pm_1)-m(\tilde\chi^0_1))/(m(\tilde g) - m(\tilde\chi^0_1)) = 1/2$ models.

The expected limits for the hard-lepton 6-jet signal region. The limits are presented in the (gluino, chargino) mass plane for the scenario where the mass of the chargino is fixed to $x=(m(\tilde\chi^\pm_1)-m(\tilde\chi^0_1))/(m(\tilde g) - m(\tilde\chi^0_1)) = 1/2$ models.

The observed limits for the soft-lepton 5-jet signal region. The limits are presented in the (gluino, x) plane for the chargino = 60 GeV models where $x=(m(\tilde\chi^\pm_1)-m(\tilde\chi^0_1))/(m(\tilde g) - m(\tilde\chi^0_1))$.

The expected limits for the soft-lepton 5-jet signal region. The limits are presented in the (gluino, x) plane for the chargino = 60 GeV models where $x=(m(\tilde\chi^\pm_1)-m(\tilde\chi^0_1))/(m(\tilde g) - m(\tilde\chi^0_1))$.

The observed limits for the hard-lepton 4-jet low-x signal region. The limits are presented in the (gluino, x) plane for the chargino = 60 GeV models where $x=(m(\tilde\chi^\pm_1)-m(\tilde\chi^0_1))/(m(\tilde g) - m(\tilde\chi^0_1))$.

The expected limits for the hard-lepton 4-jet low-x signal region. The limits are presented in the (gluino, x) plane for the chargino = 60 GeV models where $x=(m(\tilde\chi^\pm_1)-m(\tilde\chi^0_1))/(m(\tilde g) - m(\tilde\chi^0_1))$.

The observed limits for the hard-lepton 4-jet high-x signal region. The limits are presented in the (gluino, x) plane for the chargino = 60 GeV models where $x=(m(\tilde\chi^\pm_1)-m(\tilde\chi^0_1))/(m(\tilde g) - m(\tilde\chi^0_1))$.

The expected limits for the hard-lepton 4-jet high-x signal region. The limits are presented in the (gluino, x) plane for the chargino = 60 GeV models where $x=(m(\tilde\chi^\pm_1)-m(\tilde\chi^0_1))/(m(\tilde g) - m(\tilde\chi^0_1))$.

The observed limits for the hard-lepton 5-jet signal region. The limits are presented in the (gluino, x) plane for the chargino = 60 GeV models where $x=(m(\tilde\chi^\pm_1)-m(\tilde\chi^0_1))/(m(\tilde g) - m(\tilde\chi^0_1))$.

The expected limits for the hard-lepton 5-jet signal region. The limits are presented in the (gluino, x) plane for the chargino = 60 GeV models where $x=(m(\tilde\chi^\pm_1)-m(\tilde\chi^0_1))/(m(\tilde g) - m(\tilde\chi^0_1))$.

The observed limits for the hard-lepton 6-jet signal region. The limits are presented in the (gluino, x) plane for the chargino = 60 GeV models where $x=(m(\tilde\chi^\pm_1)-m(\tilde\chi^0_1))/(m(\tilde g) - m(\tilde\chi^0_1))$.

The expected limits for the hard-lepton 6-jet signal region. The limits are presented in the (gluino, x) plane for the chargino = 60 GeV models where $x=(m(\tilde\chi^\pm_1)-m(\tilde\chi^0_1))/(m(\tilde g) - m(\tilde\chi^0_1))$.

Number of generated events in the (gluino, chargino) mass plane for the scenario where the mass of the chargino is fixed to $x=(m(\tilde\chi^\pm_1)-m(\tilde\chi^0_1))/(m(\tilde g) - m(\tilde\chi^0_1)) = 1/2$.

Number of generated events in the (gluino, x) plane for the chargino = 60 GeV models.

Production cross-section in the (gluino, chargino) mass plane for the scenario where the mass of the chargino is fixed to $x=(m(\tilde\chi^\pm_1)-m(\tilde\chi^0_1))/(m(\tilde g) - m(\tilde\chi^0_1)) = 1/2$.

Production cross-section in the (gluino, x) plane for the chargino = 60 GeV models.

Acceptance times efficiency obtained in the different signal regions in the (gluino, chargino) mass plane for the scenario where the mass of the chargino is fixed to $x=(m(\tilde\chi^\pm_1)-m(\tilde\chi^0_1))/(m(\tilde g) - m(\tilde\chi^0_1)) = 1/2$ (hard-lepton 4-jet low-x).

Acceptance times efficiency in the different signal regions in the (gluino, chargino) mass plane for the scenario where the mass of the chargino is fixed to $x=(m(\tilde\chi^\pm_1)-m(\tilde\chi^0_1))/(m(\tilde g) - m(\tilde\chi^0_1)) = 1/2$ (hard-lepton 4-jet high-x).

Acceptance times efficiency in the different signal regions in the (gluino, chargino) mass plane for the scenario where the mass of the chargino is fixed to $x=(m(\tilde\chi^\pm_1)-m(\tilde\chi^0_1))/(m(\tilde g) - m(\tilde\chi^0_1)) = 1/2$ (hard-lepton 5-jet).

Acceptance times efficiency in the different signal regions in the (gluino, chargino) mass plane for the scenario where the mass of the chargino is fixed to $x=(m(\tilde\chi^\pm_1)-m(\tilde\chi^0_1))/(m(\tilde g) - m(\tilde\chi^0_1)) = 1/2$ (hard-lepton 6-jet).

Acceptance times efficiency in the different signal regions in the (gluino, chargino) mass plane for the scenario where the mass of the chargino is fixed to $x=(m(\tilde\chi^\pm_1)-m(\tilde\chi^0_1))/(m(\tilde g) - m(\tilde\chi^0_1)) = 1/2$ (soft-lepton 2-jet).

Acceptance times efficiency obtained in the different signal regions in the (gluino, x) plane for the chargino = 60 GeV models (hard-lepton 4-jet low-x).

Acceptance times efficiency in the different signal regions in the (gluino, x) plane for the chargino = 60 GeV models (hard-lepton 4-jet high-x).

Acceptance times efficiency in the different signal regions in the (gluino, x) plane for the chargino = 60 GeV models (hard-lepton 5-jet).

Acceptance times efficiency in the different signal regions in the (gluino, x) plane for the chargino = 60 GeV models (hard-lepton 6-jet).

Acceptance times efficiency in the different signal regions in the (gluino, x) plane for the chargino = 60 GeV models (soft-lepton 5-jet).

The observed CLs values as obtained in the different signal regions in the (gluino, chargino) mass plane for the scenario where the mass of the chargino is fixed to $x=(m(\tilde\chi^\pm_1)-m(\tilde\chi^0_1))/(m(\tilde g) - m(\tilde\chi^0_1)) = 1/2$ (hard-lepton 4-jet low-x).

The observed CLs values as obtained in the different signal regions in the (gluino, chargino) mass plane for the scenario where the mass of the chargino is fixed to $x=(m(\tilde\chi^\pm_1)-m(\tilde\chi^0_1))/(m(\tilde g) - m(\tilde\chi^0_1)) = 1/2$ (hard-lepton 4-jet high-x).

The observed CLs values as obtained in the different signal regions in the (gluino, chargino) mass plane for the scenario where the mass of the chargino is fixed to $x=(m(\tilde\chi^\pm_1)-m(\tilde\chi^0_1))/(m(\tilde g) - m(\tilde\chi^0_1)) = 1/2$ (hard-lepton 5-jet).

The observed CLs values as obtained in the different signal regions in the (gluino, chargino) mass plane for the scenario where the mass of the chargino is fixed to $x=(m(\tilde\chi^\pm_1)-m(\tilde\chi^0_1))/(m(\tilde g) - m(\tilde\chi^0_1)) = 1/2$ (hard-lepton 6-jet).

The observed CLs values as obtained in the different signal regions in the (gluino, chargino) mass plane for the scenario where the mass of the chargino is fixed to $x=(m(\tilde\chi^\pm_1)-m(\tilde\chi^0_1))/(m(\tilde g) - m(\tilde\chi^0_1)) = 1/2$ (soft-lepton 2-jet).

The expected CLs values as obtained in the different signal regions in the (gluino, chargino) mass plane for the scenario where the mass of the chargino is fixed to $x=(m(\tilde\chi^\pm_1)-m(\tilde\chi^0_1))/(m(\tilde g) - m(\tilde\chi^0_1)) = 1/2$ (hard-lepton 4-jet low-x).

The expected CLs values as obtained in the different signal regions in the (gluino, chargino) mass plane for the scenario where the mass of the chargino is fixed to $x=(m(\tilde\chi^\pm_1)-m(\tilde\chi^0_1))/(m(\tilde g) - m(\tilde\chi^0_1)) = 1/2$ (hard-lepton 4-jet high-x).

The expected CLs values as obtained in the different signal regions in the (gluino, chargino) mass plane for the scenario where the mass of the chargino is fixed to $x=(m(\tilde\chi^\pm_1)-m(\tilde\chi^0_1))/(m(\tilde g) - m(\tilde\chi^0_1)) = 1/2$ (hard-lepton 5-jet).

The expected CLs values as obtained in the different signal regions in the (gluino, chargino) mass plane for the scenario where the mass of the chargino is fixed to $x=(m(\tilde\chi^\pm_1)-m(\tilde\chi^0_1))/(m(\tilde g) - m(\tilde\chi^0_1)) = 1/2$ (hard-lepton 6-jet).

The expected CLs values as obtained in the different signal regions in the (gluino, chargino) mass plane for the scenario where the mass of the chargino is fixed to $x=(m(\tilde\chi^\pm_1)-m(\tilde\chi^0_1))/(m(\tilde g) - m(\tilde\chi^0_1)) = 1/2$ (soft-lepton 2-jet).

The observed CLs values as obtained in the different signal regions in the (gluino, x) plane for the chargino = 60 GeV models (hard-lepton 4-jet low-x).

The observed CLs values as obtained in the different signal regions in the (gluino, x) plane for the chargino = 60 GeV models (hard-lepton 4-jet high-x).

The observed CLs values as obtained in the different signal regions in the (gluino, x) plane for the chargino = 60 GeV models (hard-lepton 5-jet).

The observed CLs values as obtained in the different signal regions in the (gluino, x) plane for the chargino = 60 GeV models (hard-lepton 6-jet).

The observed CLs values as obtained in the different signal regions in the (gluino, x) plane for the chargino = 60 GeV models (soft-lepton 5-jet).

The expected CLs values as obtained in the different signal regions in the (gluino, x) plane for the chargino = 60 GeV models (hard-lepton 4-jet low-x).

The expected CLs values as obtained in the different signal regions in the (gluino, x) plane for the chargino = 60 GeV models (hard-lepton 4-jet high-x).

The expected CLs values as obtained in the different signal regions in the (gluino, x) plane for the chargino = 60 GeV models (hard-lepton 5-jet).

The expected CLs values as obtained in the different signal regions in the (gluino, x) plane for the chargino = 60 GeV models (hard-lepton 6-jet).

The expected CLs values as obtained in the different signal regions in the (gluino, x) plane for the chargino = 60 GeV models (soft-lepton 5-jet).

The signal region yielding in the best expected limit is indicated for every signal point used in the the gluino simplified models for the (gluino, chargino) mass plane for the scenario where the mass of the chargino is fixed to $x=(m(\tilde\chi^\pm_1)-m(\tilde\chi^0_1))/(m(\tilde g) - m(\tilde\chi^0_1)) = 1/2$.

The signal region yielding in the best expected limit is indicated for every signal point used in the the gluino simplified models for the (gluino, x) mass plane where for the chargino = 60 GeV and $x=(m(\tilde\chi^\pm_1)-m(\tilde\chi^0_1))/(m(\tilde g) - m(\tilde\chi^0_1))$.

Model-dependent 95% CL upper limits on the visible cross-section in addition to the observed and expected exclusion limits for the (gluino, chargino) mass plane for the scenario where the mass of the chargino is fixed to $x=(m(\tilde\chi^\pm_1)-m(\tilde\chi^0_1))/(m(\tilde g) - m(\tilde\chi^0_1)) = 1/2$.

Model-dependent 95% CL upper limits on the visible cross-section in addition to the observed and expected exclusion limits for the (gluino, x) mass plane where for the chargino = 60 GeV and $x=(m(\tilde\chi^\pm_1)-m(\tilde\chi^0_1))/(m(\tilde g) - m(\tilde\chi^0_1))$.

Simulated background event samples: the corresponding generator, parton shower, cross-section normalisation, PDF set and underlying-event tune are shown.

Overview of the selection criteria for the soft-lepton signal regions. The symbol $p_{T}^{l}$ refers to signal leptons.

Overview of the selection criteria for the hard-lepton signal regions. The symbol $p_{T}^{l}$ refers to signal leptons.

Background fit results for the hard-lepton and soft-lepton signal regions, for an integrated luminosity of 3.2 fb-1. Uncertainties in the fitted background estimates combine statistical (in the simulated event yields) and systematic uncertainties. The uncertainties in this table are symmetrised for propagation purposes but truncated at zero to remain within the physical boundaries.

Breakdown of upper limits. The columns show from left to right: the name of the respective signal region; the 95% confidence level (CL) upper limits on the visible cross-section and on the number of signal events the 95% CL upper limit on the number of signal events, given the expected number (and $\pm 1 \sigma$ variations on the expectation) of background events; the two-sided CLb value, i.e. the confidence level observed for the background-only hypothesis and the one-sided discovery p-value (p(s = 0)). The discovery p-values are capped to 0.5 in the case of observing less events than the fitted background estimates.

Table shows the data, fitted background and expected signal event counts for a benchmark signal point in each bin of the mt distribution shown in figure 5 (top left). The fit results are shown for an integrated luminosity of 3.2 fb-1.

Table shows the data, fitted background and expected signal event counts for a benchmark signal point in each bin of the met/meff distribution shown in figure 5 (top right). The fit results are shown for an integrated luminosity of 3.2 fb-1.

Table shows the data, fitted background and expected signal event counts for a benchmark signal point in each bin of the mt distribution shown in figure 5 (middle left). The fit results are shown for an integrated luminosity of 3.2 fb-1.

Table shows the data, fitted background and expected signal event counts for a benchmark signal point in each bin of the mt distribution shown in figure 5 (middle right). The fit results are shown for an integrated luminosity of 3.2 fb-1.

Table shows the data, fitted background and expected signal event counts for a benchmark signal point in each bin of the met distribution shown in figure 5 (bottom left). The fit results are shown for an integrated luminosity of 3.2 fb-1.

Table shows the data, fitted background and expected signal event counts for a benchmark signal point in each bin of the met distribution shown in figure 5 (bottom right). The fit results are shown for an integrated luminosity of 3.2 fb-1.

Cutflow table for the hard-lepton signal regions with representative target signal models. The weighted numbers are normalized to 3.2 fb-1 and rounded to the statistical error.

Cutflow table for the hard-lepton signal regions with representative target signal models. The weighted numbers are normalized to 3.2 fb-1 and rounded to the statistical error.