Expected value of the test statistic (-2ln$\Lambda$), as a function of the $\kappa_\lambda$ parameter for $b\bar{b}\gamma\gamma$.

Expected value of the test statistic (-2ln$\Lambda$), as a function of the $\kappa_\lambda$ parameter for $b\bar{b}\ell\ell+E_\text{T}^{\text{miss}}$.

Expected value of the test statistic (-2ln$\Lambda$), as a function of the $\kappa_\lambda$ parameter for Multilepton.

Expected value of the test statistic (-2ln$\Lambda$), as a function of the $\kappa_\lambda$ parameter for combination.

Expected value of the test statistic (-2ln$\Lambda$), as a function of the $\kappa_{2V}$ parameter for $b\bar{b}b\bar{b}$.

Expected value of the test statistic (-2ln$\Lambda$), as a function of the $\kappa_{2V}$ parameter for $b\bar{b}\tau\tau$.

Expected value of the test statistic (-2ln$\Lambda$), as a function of the $\kappa_{2V}$ parameter for $b\bar{b}\gamma\gamma$.

Expected value of the test statistic (-2ln$\Lambda$), as a function of the $\kappa_{2V}$ parameter for $b\bar{b}\ell\ell+E_\text{T}^{\text{miss}}$.

Expected value of the test statistic (-2ln$\Lambda$), as a function of the $\kappa_{2V}$ parameter for Multilepton.

Expected value of the test statistic (-2ln$\Lambda$), as a function of the $\kappa_{2V}$ parameter for combination.

Observed and expected 95% CL upper limits on the signal strength for the inclusive ggF HH production from the bb̄τ⁺τ^-, bb̄γγ, bb̄bb̄, multilepton and bb̄ℓℓ+E_T^miss decay channels, and their statistical combination. When deriving the limit on μ_ggF^HH (μ_VBF^HH), the VBF (ggF) HH production cross-section is fixed to the SM predicted value for m_H=125 GeV. The expected limit, along with the ±1σ and ±2σ bands, is calculated under the assumption of no HH process and with all NPs profiled to the observed data.

Observed and expected 95% CL upper limits on the signal strength for the VBF HH production from the bb̄τ⁺τ^-, bb̄γγ, bb̄bb̄, multilepton and bb̄ℓℓ+E_T^miss decay channels, and their statistical combination. When deriving the limit on μ_ggF^HH (μ_VBF^HH), the VBF (ggF) HH production cross-section is fixed to the SM predicted value for m_H=125 GeV. The expected limit, along with the ±1σ and ±2σ bands, is calculated under the assumption of no HH process and with all NPs profiled to the observed data.

Observed and expected upper limits at 95% CL on the inclusive ggF and VBF HH production cross-section from the bb̄τ⁺τ^-, bb̄γγ, bb̄bb̄, multilepton and bb̄ℓℓ+E_T^miss decay channels, and their statistical combination. The Higgs boson mass is set to 125 GeV when deriving the predicted SM cross-section (red band). The expected limits, along with the ±1σ and ±2σ bands, are calculated under the assumption of no Higgs boson pair production and with all NPs profiled to the observed data.

Observed value of the test statistic (-2ln$\Lambda$), as a function of the $\kappa_\lambda$ parameter for $b\bar{b}b\bar{b}$.

Observed value of the test statistic (-2ln$\Lambda$), as a function of the $\kappa_\lambda$ parameter for $b\bar{b}\tau\tau$.

Observed value of the test statistic (-2ln$\Lambda$), as a function of the $\kappa_\lambda$ parameter for $b\bar{b}\gamma\gamma$.

Observed value of the test statistic (-2ln$\Lambda$), as a function of the $\kappa_\lambda$ parameter for $b\bar{b}\ell\ell+E_\text{T}^{\text{miss}}$.

Observed value of the test statistic (-2ln$\Lambda$), as a function of the $\kappa_\lambda$ parameter for Multilepton.

Observed value of the test statistic (-2ln$\Lambda$), as a function of the $\kappa_\lambda$ parameter for combination.

Observed value of the test statistic (-2ln$\Lambda$), as a function of the $\kappa_{2V}$ parameter for $b\bar{b}b\bar{b}$.

Observed value of the test statistic (-2ln$\Lambda$), as a function of the $\kappa_{2V}$ parameter for $b\bar{b}\tau\tau$.

Observed value of the test statistic (-2ln$\Lambda$), as a function of the $\kappa_{2V}$ parameter for $b\bar{b}\gamma\gamma$.

Observed value of the test statistic (-2ln$\Lambda$), as a function of the $\kappa_{2V}$ parameter for $b\bar{b}\ell\ell+E_\text{T}^{\text{miss}}$.

Observed value of the test statistic (-2ln$\Lambda$), as a function of the $\kappa_{2V}$ parameter for Multilepton.

Observed value of the test statistic (-2ln$\Lambda$), as a function of the $\kappa_{2V}$ parameter for combination.

Observed (solid lines) and expected (dashed lines) 95% CL exclusion limits on the HH production cross-sections of (a) the inclusive ggF and VBF processes as a function of κ_λ, for the bb̄γγ (purple), bb̄τ⁺τ^- (green), multilepton (cyan), bb̄bb̄ (blue) and bb̄ℓℓ+E_T^miss (brown) decay channels and their combination (black). The expected limits assume no HH production. The red line shows the theory prediction for the ggF and VBF HH production cross-section as a function of κ_λ. The bands surrounding the red cross-section lines indicate the theoretical uncertainties on the predicted cross-section.

Observed (solid lines) and expected (dashed lines) 95% CL exclusion limits on the HH production cross-sections of the VBF process as a function of κ_2V, for the bb̄γγ (purple), bb̄τ⁺τ^- (green), multilepton (cyan), bb̄bb̄ (blue) and bb̄ℓℓ+E_T^miss (brown) decay channels and their combination (black). The expected limits assume no VBF HH production. The ggF HH production cross-section is assumed to be as predicted by the SM. The ggF HH production cross-section is assumed to be as predicted by the SM. The red line shows the predicted VBF HH cross-section as a function of κ_2V. The bands surrounding the red cross-section lines indicate the theoretical uncertainties on the predicted cross-section.

Observed value of the test statistic (-2ln$\Lambda$), as a function of the $c_{gghh}$ parameter for $b\bar{b}b\bar{b}$.

Observed value of the test statistic (-2ln$\Lambda$), as a function of the $c_{gghh}$ parameter for $b\bar{b}\tau\tau$.

Observed value of the test statistic (-2ln$\Lambda$), as a function of the $c_{gghh}$ parameter for $b\bar{b}\gamma\gamma$.

Observed value of the test statistic (-2ln$\Lambda$), as a function of the $c_{gghh}$ parameter for combination.

Expected value of the test statistic (-2ln$\Lambda$), as a function of the $c_{gghh}$ parameter for $b\bar{b}b\bar{b}$.

Expected value of the test statistic (-2ln$\Lambda$), as a function of the $c_{gghh}$ parameter for $b\bar{b}\tau\tau$.

Expected value of the test statistic (-2ln$\Lambda$), as a function of the $c_{gghh}$ parameter for $b\bar{b}\gamma\gamma$.

Expected value of the test statistic (-2ln$\Lambda$), as a function of the $c_{gghh}$ parameter for combination.

Observed value of the test statistic (-2ln$\Lambda$), as a function of the $c_{tthh}$ parameter for $b\bar{b}b\bar{b}$.

Observed value of the test statistic (-2ln$\Lambda$), as a function of the $c_{tthh}$ parameter for $b\bar{b}\tau\tau$.

Observed value of the test statistic (-2ln$\Lambda$), as a function of the $c_{tthh}$ parameter for $b\bar{b}\gamma\gamma$.

Observed value of the test statistic (-2ln$\Lambda$), as a function of the $c_{tthh}$ parameter for combination.

Expected value of the test statistic (-2ln$\Lambda$), as a function of the $c_{tthh}$ parameter for $b\bar{b}b\bar{b}$.

Expected value of the test statistic (-2ln$\Lambda$), as a function of the $c_{tthh}$ parameter for $b\bar{b}\tau\tau$.

Expected value of the test statistic (-2ln$\Lambda$), as a function of the $c_{tthh}$ parameter for $b\bar{b}\gamma\gamma$.

Expected value of the test statistic (-2ln$\Lambda$), as a function of the $c_{tthh}$ parameter for combination.

Observed and expected 95% CL combined upper limits on the cross-section for the SM and seven BSM HEFT benchmarks (arXiv:2304.01968 [hep-ph]) in the ggF process, describing representative signal kinematics and m_HH shape features obtained by varying multiple HEFT coefficients. The expected limits from the bb̄τ⁺τ^-, bb̄γγ and bb̄bb̄ decay channels are presented as well. Theoretical predictions, estimated using specific sets of coefficient values defined in the benchmarks, are shown as red cross dots.

Expected (dashed line) value of $-2\ln\Lambda$ as a function of $\kappa_{\lambda}$ for the combined fit, when all other coupling modifiers are fixed to their SM predictions.

Observed (solid line) value of $-2\ln\Lambda$ as a function of $\kappa_{2V}$ for the combined fit, when all other coupling modifiers are fixed to their SM predictions.

Expected (dashed line) value of $-2\ln\Lambda$ as a function of $\kappa_{2V}$ for the combined fit, when all other coupling modifiers are fixed to their SM predictions.

Likelihood observed contours at 68% CL (solid line) and 95% CL (dashed line) in the ($\kappa_{\lambda}$, $\kappa_{2V}$) parameter space, when all other coupling modifiers are fixed to their SM predictions. The best-fit value, denoted by a cross in the plot, corresponds to the very first entry in the table.

Likelihood expected contours at 68% CL (teal shaded region) and 95% CL (yellow shaded region) in the ($\kappa_{\lambda}$, $\kappa_{2V}$) parameter space, when all other coupling modifiers are fixed to their SM predictions. In the plot, the SM prediction is indicated by the star.

Likelihood contours at 68% (solid line) and 95% (dashed line) CL in the ($c_{hhh}$, $c_{gghh}$) parameter space, when all other Wilson coefficients are fixed to their SM values. The corresponding expected contours are shown by the inner and outer shaded regions. The SM prediction is indicated by the star, while the observed best-fit value is denoted by the black cross. The best-fit value corresponds to the very first entry in the table.

Likelihood contours at 68% (solid line) and 95% (dashed line) CL in the ($c_{hhh}$, $c_{gghh}$) parameter space, when all other Wilson coefficients are fixed to their SM values. The corresponding expected contours are shown by the inner and outer shaded regions. The SM prediction is indicated by the star, while the observed best-fit value is denoted by the black cross.

Likelihood contours at 68% (solid line) and 95% (dashed line) CL in the ($c_{hhh}$, $c_{tthh}$) parameter space, when all other Wilson coefficients are fixed to their SM values. The corresponding expected contours are shown by the inner and outer shaded regions. The SM prediction is indicated by the star, while the observed best-fit value is denoted by the black cross. The best-fit value corresponds to the very first entry in the table.

Likelihood contours at 68% (solid line) and 95% (dashed line) CL in the ($c_{hhh}$, $c_{tthh}$) parameter space, when all other Wilson coefficients are fixed to their SM values. The corresponding expected contours are shown by the inner and outer shaded regions. The SM prediction is indicated by the star, while the observed best-fit value is denoted by the black cross.

Likelihood contours at 68% (solid line) and 95% (dashed line) CL in the ($c_{H}$, $c_{H◻}$) parameter space, when all other Wilson coefficients are fixed to their SM values. The corresponding expected contours are shown by the inner and outer shaded regions. The SM prediction is indicated by the star, while the observed best-fit value is denoted by the black cross. The best-fit value corresponds to the very first entry in the table.

Likelihood contours at 68% (solid line) and 95% (dashed line) CL in the ($c_{H}$, $c_{H◻}$) parameter space, when all other Wilson coefficients are fixed to their SM values. The corresponding expected contours are shown by the inner and outer shaded regions. The SM prediction is indicated by the star, while the observed best-fit value is denoted by the black cross.

The acceptance times efficiency [in %] for the signal ggF $HH$ process as a function of the coupling modifier $\kappa_{\lambda}$ in each analysis category for the $\tau_{\mathrm{had}}\tau_{\mathrm{had}}$ channel. The other coupling modifiers affecting ggF $HH$ production are fixed to their SM predictions.

The acceptance times efficiency [in %] for the signal ggF $HH$ process as a function of the coupling modifier $\kappa_{\lambda}$ in each analysis category for the $\tau_{\mathrm{lep}}\tau_{\mathrm{had}}$ SLT channel. The other coupling modifiers affecting ggF $HH$ production are fixed to their SM predictions.

The acceptance times efficiency [in %] for the signal ggF $HH$ process as a function of the coupling modifier $\kappa_{\lambda}$ in each analysis category for the $\tau_{\mathrm{lep}}\tau_{\mathrm{had}}$ LTT channel. The other coupling modifiers affecting ggF $HH$ production are fixed to their SM predictions.

The acceptance times efficiency [in %] for the signal VBF $HH$ process as a function of the coupling modifier $\kappa_{\lambda}$ in each analysis category for the $\tau_{\mathrm{had}}\tau_{\mathrm{had}}$ channel. The other coupling modifiers affecting VBF $HH$ production are fixed to their SM predictions.

The acceptance times efficiency [in %] for the signal VBF $HH$ process as a function of the coupling modifier $\kappa_{\lambda}$ in each analysis category for the $\tau_{\mathrm{lep}}\tau_{\mathrm{had}}$ SLT channel. The other coupling modifiers affecting VBF $HH$ production are fixed to their SM predictions.

The acceptance times efficiency [in %] for the signal VBF $HH$ process as a function of the coupling modifier $\kappa_{\lambda}$ in each analysis category for the $\tau_{\mathrm{lep}}\tau_{\mathrm{had}}$ LTT channel. The other coupling modifiers affecting VBF $HH$ production are fixed to their SM predictions.

The acceptance times efficiency [in %] for the signal VBF $HH$ process as a function of the coupling modifier $\kappa_{2V}$ in each analysis category for the $\tau_{\mathrm{had}}\tau_{\mathrm{had}}$ channel. The other coupling modifiers affecting VBF $HH$ production are fixed to their SM predictions.

The acceptance times efficiency [in %] for the signal VBF $HH$ process as a function of the coupling modifier $\kappa_{2V}$ in each analysis category for the $\tau_{\mathrm{lep}}\tau_{\mathrm{had}}$ SLT channel. The other coupling modifiers affecting VBF $HH$ production are fixed to their SM predictions.

The acceptance times efficiency [in %] for the signal VBF $HH$ process as a function of the coupling modifier $\kappa_{2V}$ in each analysis category for the $\tau_{\mathrm{lep}}\tau_{\mathrm{had}}$ LTT channel. The other coupling modifiers affecting VBF $HH$ production are fixed to their SM predictions.

Expected (dashed line) value of $-2\ln\Lambda$ as a function of $\kappa_{\lambda}$ for the fit to the $\tau_{\mathrm{had}}\tau_{\mathrm{had}}$ channel, when all other coupling modifiers are fixed to their SM predictions.

Expected (dashed line) value of $-2\ln\Lambda$ as a function of $\kappa_{\lambda}$ for the fit to the $\tau_{\mathrm{lep}}\tau_{\mathrm{had}}$ channel, when all other coupling modifiers are fixed to their SM predictions.

Observed (solid line) value of $-2\ln\Lambda$ as a function of $\kappa_{\lambda}$ for the fit to the $\tau_{\mathrm{had}}\tau_{\mathrm{had}}$ channel, when all other coupling modifiers are fixed to their SM predictions.

Observed (solid line) value of $-2\ln\Lambda$ as a function of $\kappa_{\lambda}$ for the fit to the $\tau_{\mathrm{lep}}\tau_{\mathrm{had}}$ channel, when all other coupling modifiers are fixed to their SM predictions.

Expected (dashed line) value of $-2\ln\Lambda$ as a function of $\kappa_{2V}$ for the fit to the $\tau_{\mathrm{had}}\tau_{\mathrm{had}}$ channel, when all other coupling modifiers are fixed to their SM predictions.

Expected (dashed line) value of $-2\ln\Lambda$ as a function of $\kappa_{2V}$ for the fit to the $\tau_{\mathrm{lep}}\tau_{\mathrm{had}}$ channel, when all other coupling modifiers are fixed to their SM predictions.

Observed (solid line) value of $-2\ln\Lambda$ as a function of $\kappa_{2V}$ for the fit to the $\tau_{\mathrm{had}}\tau_{\mathrm{had}}$ channel, when all other coupling modifiers are fixed to their SM predictions.

Observed (solid line) value of $-2\ln\Lambda$ as a function of $\kappa_{2V}$ for the fit to the $\tau_{\mathrm{lep}}\tau_{\mathrm{had}}$ channel, when all other coupling modifiers are fixed to their SM predictions.

Summary of the expected (blue) and observed (orange) one-dimensional constraints at 68% (solid error bars) and 95% (dashed error bars) CL on the coefficients of selected HEFT (top) and SMEFT (bottom) operators describing Higgs boson interactions affecting $HH$ production through gluon-gluon fusion. The results are obtained from one dimensional scans of the profile log-likelihood assuming that all other Wilson coefficients are fixed to their SM values. The contribution from VBF $HH$ production is neglected.