Values of the test statistic $t_{\kappa_{g,\mathrm{off-shell}}}$ as a function of $\kappa_{g,\mathrm{off-shell}}$ while profiling $\kappa_{V,\mathrm{off-shell}}$ obtained with an Asimov dataset and with data. The values with all nuisance parameters fixed at their best-fit values (stat-only) are added for comparison. The 68% and 95% confidence intervals obtained from the Neyman construction are also added.

Values of the test statistic $t_{\kappa_{V,\mathrm{off-shell}}}$ as a function of $\kappa_{V,\mathrm{off-shell}}$ while profiling $\kappa_{g,\mathrm{off-shell}}$ obtained with an Asimov dataset and with data. The values with all nuisance parameters fixed at their best-fit values (stat-only) are added for comparison. The 68% and 95% confidence intervals obtained from the Neyman construction are also added.

Values of the test statistic $t_{\kappa_{H}}$ as a function of $\kappa_{H}=\Gamma_H/\Gamma_H^{\mathrm{SM}}$ obtained with an Asimov dataset and with data. The values with all nuisance parameters fixed at their best-fit values (stat-only) are added for comparison. The 68% and 95% confidence intervals obtained from the Neyman construction are also added.

Values of the test statistic $t_{R_{gg}}$ as a function of $R_{gg}=\kappa_{g,\mathrm{on-shell}}^2/\kappa_{g,\mathrm{off-shell}}^2$ while profiling $R_{VV}=\kappa_{V,\mathrm{on-shell}}^2/\kappa_{V,\mathrm{off-shell}}^2$ and fixing $\kappa_H=1$ obtained with an Asimov dataset and with data. The values with all nuisance parameters fixed at their best-fit values (stat-only) are added for comparison. The 68% and 95% confidence intervals obtained from the Neyman construction are also added.

Values of the test statistic $t_{R_{VV}}$ as a function of $R_{VV}=\kappa_{V,\mathrm{on-shell}}^2/\kappa_{V,\mathrm{off-shell}}^2$ while profiling $R_{gg}=\kappa_{g,\mathrm{on-shell}}^2/\kappa_{g,\mathrm{off-shell}}^2$ and fixing $\kappa_H=1$ obtained with an Asimov dataset and with data. The values with all nuisance parameters fixed at their best-fit values (stat-only) are added for comparison. The 68% and 95% confidence intervals obtained from the Neyman construction are also added.

Values of the test statistic $t_{\mu_{\mathrm{off-shell}}}$ assuming a single parameter of interest $\mu_{\mathrm{off-shell}}$ obtained with an Asimov dataset and with data when combining the $H^*\rightarrow ZZ\rightarrow 4\ell$ and $H^*\rightarrow ZZ\rightarrow 2\ell 2\nu$ decay channels. The values from the histogram-based analysis (Phys. Lett. B 846 (2023) 138223) are added for comparison. The 68% and 95% confidence intervals obtained from the Neyman construction are also added.

Values of the test statistic $t_{\kappa_{H}}$ as a function of $\kappa_{H}=\Gamma_H/\Gamma_H^{\mathrm{SM}}$ obtained with an Asimov dataset and with data in the $H^*\rightarrow ZZ\rightarrow 4\ell$ decay channel. The values from the histogram-based analysis (Phys. Lett. B 846 (2023) 138223) are added for comparison. The 68% and 95% confidence intervals obtained from the Neyman construction are also added.

Values of the test statistic $t_{\kappa_{H}}$ as a function of $\kappa_{H}=\Gamma_H/\Gamma_H^{\mathrm{SM}}$ obtained with an Asimov dataset and with data when combining the $H^*\rightarrow ZZ\rightarrow 4\ell$ and $H^*\rightarrow ZZ\rightarrow 2\ell 2\nu$ decay channels. The values from the histogram-based analysis (Phys. Lett. B 846 (2023) 138223) are added for comparison. The 68% and 95% confidence intervals obtained from the Neyman construction are also added.

Values of the test statistic $t_{R_{gg}}$ as a function of $R_{gg}$ while profiling $R_{VV}$ and fixing $\kappa_H=1$ obtained with an Asimov dataset and with data. The values from the histogram-based analysis (Phys. Lett. B 846 (2023) 138223) are added for comparison. The 68% and 95% confidence intervals obtained from the Neyman construction are also added.

Values of the test statistic $t_{R_{VV}}$ as a function of $R_{VV}$ while profiling $R_{gg}$ and fixing $\kappa_H=1$ obtained with an Asimov dataset and with data. The values from the histogram-based analysis (Phys. Lett. B 846 (2023) 138223) are added for comparison. The 68% and 95% confidence intervals obtained from the Neyman construction are also added.

The pNN input variable visible transverse momentum $p_{\mathrm{T}}^{\mathrm{vis}}(\mu\tau_{\mathrm{had}})$ is shown in the SR with no cut on the pNN discriminant. The signal shape is normalized to the same integral as the total background prediction. Overflow events are included in the last bins.

The pNN input variable transverse mass $m_{\mathrm{T}}(\mu,E_{\mathrm{T}}^{\mathrm{miss}})$ is shown in the SR with no cut on the pNN discriminant. The signal shape is normalized to the same integral as the total background prediction. Overflow events are included in the last bins.

The pNN input variable MMC mass $m^{\mathrm{MMC}}(bb\tau\tau)$ is shown in the SR with no cut on the pNN discriminant. The signal shape is normalized to the same integral as the total background prediction. Overflow events are included in the last bins.

Mass distribution of the $\tau^+\tau^-$ system for each of the generated MC signal mass points for the $\mu\tau_{\mathrm{h}}$ category. All distributions are normalized assuming $\mathcal{B}(H\rightarrow aa\rightarrow b\bar{b}\tau^+\tau^-)=10\%$.

Mass distribution of the $\tau^+\tau^-$ system for each of the generated MC signal mass points for the $e\mu$ category. All distributions are normalized assuming $\mathcal{B}(H\rightarrow aa\rightarrow b\bar{b}\tau^+\tau^-)=10\%$.

The pNN spectrum in the three bins are shown separately for each mass hypothesis in the ($e\mu$, $1b$) category. The signal shape (normalized to $\mathcal{B}(H\rightarrow aa\rightarrow b\bar{b}\tau^+\tau^-)=1$) for each corresponding mass hypothesis is overlaid on top of the SM prediction. Bins used for different $m_a$ hypotheses are not statistically independent.

The pNN spectrum in the three bins are shown separately for each mass hypothesis in the ($\mu\tau_{\mathrm{had}}$,$2b$) category. The signal shape (normalized to $\mathcal{B}(H\rightarrow aa\rightarrow b\bar{b}\tau^+\tau^-)=1$) for each corresponding mass hypothesis is overlaid on top of the SM prediction. Bins used for different $m_a$ hypotheses are not statistically independent.

The pNN spectrum in the three bins are shown separately for each mass hypothesis in the ($e\mu$,$1B$) category. The signal shape (normalized to $\mathcal{B}(H\rightarrow aa\rightarrow b\bar{b}\tau^+\tau^-)=1$) for each corresponding mass hypothesis is overlaid on top of the SM prediction. Bins used for different $m_a$ hypotheses are not statistically independent.

The pNN spectrum in the three bins are shown separately for each mass hypothesis in the ($e\tau_{\mathrm{had}}$,$2b$) category. The signal shape (normalized to $\mathcal{B}(H\rightarrow aa\rightarrow b\bar{b}\tau^+\tau^-)=1$) for each corresponding mass hypothesis is overlaid on top of the SM prediction. Bins used for different $m_a$ hypotheses are not statistically independent.

The observed (solid) 95% C.L. upper limits on $(\sigma(H)/\sigma_{\mathrm{SM}}(H)) \mathcal{B}(H\rightarrow aa\rightarrow b\bar{b}\tau^+\tau^-)$ as a function of $m_a$ and the expected (dashed) limits under the background-only hypothesis when combining all categories. The inner green and outer yellow shaded bands show the $\pm 1\sigma$ and $\pm 2\sigma$ uncertainties of the expected limits.

The observed (solid) 95% C.L. upper limits on $(\sigma(H)/\sigma_{\mathrm{SM}}(H)) \mathcal{B}(H\rightarrow aa\rightarrow b\bar{b}\tau^+\tau^-)$ as a function of $m_a$ and the expected (dashed) limits under the background-only hypothesis when considering different categories based on the heavy-flavor objects separately. When the category is not used for a mass hypothesis, the limit is listed as zero.

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.

Measured total $t\bar{t}\gamma$ production and decay fiducial inclusive cross-sections in both decay channels and in the combination.

Summary of the impact of the systematic uncertainties on the $t\bar{t}\gamma$ fiducial inclusive cross-section in the single-lepton and dilepton channels and their combination grouped into different categories. The quoted relative uncertainties are obtained by repeating the fit, fixing a set of nuisance parameters of the sources corresponding to each category to their post-fit values, and subtracting in quadrature the resulting uncertainty from the total uncertainty of the nominal fit. The total uncertainty is different from the sum in quadrature of the components due to correlations among nuisance parameters.

Absolute $t\bar{t}\gamma$ production differential cross-sections measured in fiducial phase space in the single-lepton channel as a function of the $\Delta R(\gamma,b)_{min}$. The last bin of the distribution includes the overflow events.

Normalised $t\bar{t}\gamma$ production differential cross-sections measured in fiducial phase space in the single-lepton channel as a function of the $\Delta R(\gamma,b)_{min}$. The last bin of the distribution includes the overflow events.

Absolute $t\bar{t}\gamma$ production differential cross-sections measured in fiducial phase space in the single-lepton channel as a function of the $\Delta R(\gamma,l)$. The last bin of the distribution includes the overflow events.

Normalised $t\bar{t}\gamma$ production differential cross-sections measured in fiducial phase space in the single-lepton channel as a function of the $\Delta R(\gamma,l)$. The last bin of the distribution includes the overflow events.

Absolute $t\bar{t}\gamma$ production differential cross-sections measured in fiducial phase space in the single-lepton channel as a function of the $\Delta R(l,j)_{min}$. The last bin of the distribution includes the overflow events.

Normalised $t\bar{t}\gamma$ production differential cross-sections measured in fiducial phase space in the single-lepton channel as a function of the $\Delta R(l,j)_{min}$. The last bin of the distribution includes the overflow events.

Absolute $t\bar{t}\gamma$ production differential cross-sections measured in fiducial phase space in the dilepton channel as a function of the $\Delta R(\gamma,b)_{min}$. The last bin of the distribution includes the overflow events.

Normalised $t\bar{t}\gamma$ production differential cross-sections measured in fiducial phase space in the dilepton channel as a function of the $\Delta R(\gamma,b)_{min}$. The last bin of the distribution includes the overflow events.

Absolute $t\bar{t}\gamma$ production differential cross-sections measured in fiducial phase space in the dilepton channel as a function of the $\Delta R(\gamma,l)_{min}$. The last bin of the distribution includes the overflow events.

Normalised $t\bar{t}\gamma$ production differential cross-sections measured in fiducial phase space in the dilepton channel as a function of the $\Delta R(\gamma,l)_{min}$. The last bin of the distribution includes the overflow events.

Absolute $t\bar{t}\gamma$ production differential cross-sections measured in fiducial phase space in the dilepton channel as a function of the $\Delta R(l,j)_{min}$. The last bin of the distribution includes the overflow events.

Normalised $t\bar{t}\gamma$ production differential cross-sections measured in fiducial phase space in the dilepton channel as a function of the $\Delta R(l,j)_{min}$. The last bin of the distribution includes the overflow events.

Absolute $t\bar{t}\gamma$ production differential cross-sections measured in the combined fiducial phase space of the single-lepton and dilepton channels as a function of the photon $p_T$. The last bin of the distributions includes the overflow events.

Normalised $t\bar{t}\gamma$ production differential cross-sections measured in the combined fiducial phase space of the single-lepton and dilepton channels as a function of the photon $p_T$. The last bin of the distributions includes the overflow events.

Absolute $t\bar{t}\gamma$ production differential cross-sections measured in the combined fiducial phase space of the single-lepton and dilepton channels as a function of the photon $|\eta|$.

Normalised $t\bar{t}\gamma$ production differential cross-sections measured in the combined fiducial phase space of the single-lepton and dilepton channels as a function of the photon $|\eta|$.

Absolute differential cross-section of the total $t\bar{t}\gamma$ production and decay measured in fiducial phase space in the single-lepton channel as a function of the $p_T(\gamma)$. The last bin of the distribution includes the overflow events.

Normalised differential cross-section of the total $t\bar{t}\gamma$ production and decay measured in fiducial phase space in the single-lepton channel as a function of the $p_T(\gamma)$. The last bin of the distribution includes the overflow events.

Absolute differential cross-section of the total $t\bar{t}\gamma$ production and decay measured in fiducial phase space in the single-lepton channel as a function of the $|\eta(\gamma)|$.

Normalised differential cross-section of the total $t\bar{t}\gamma$ production and decay measured in fiducial phase space in the single-lepton channel as a function of the $|\eta(\gamma)|$.

Absolute differential cross-section of the total $t\bar{t}\gamma$ production and decay measured in fiducial phase space in the single-lepton channel as a function of the $p_T(j_1)$. The last bin of the distribution includes the overflow events.

Normalised differential cross-section of the total $t\bar{t}\gamma$ production and decay measured in fiducial phase space in the single-lepton channel as a function of the $p_T(j_1)$. The last bin of the distribution includes the overflow events.

Absolute differential cross-section of the total $t\bar{t}\gamma$ production and decay measured in fiducial phase space in the single-lepton channel as a function of the $\Delta R(\gamma,b)_{min}$. The last bin of the distribution includes the overflow events.

Normalised differential cross-section of the total $t\bar{t}\gamma$ production and decay measured in fiducial phase space in the single-lepton channel as a function of the $\Delta R(\gamma,b)_{min}$. The last bin of the distribution includes the overflow events.

Absolute differential cross-section of the total $t\bar{t}\gamma$ production and decay measured in fiducial phase space in the single-lepton channel as a function of the $\Delta R(\gamma,l)$. The last bin of the distribution includes the overflow events.

Normalised differential cross-section of the total $t\bar{t}\gamma$ production and decay measured in fiducial phase space in the single-lepton channel as a function of the $\Delta R(\gamma,l)$. The last bin of the distribution includes the overflow events.

Absolute differential cross-section of the total $t\bar{t}\gamma$ production and decay measured in fiducial phase space in the single-lepton channel as a function of the $\Delta R(l,j)_{min}$. The last bin of the distribution includes the overflow events.

Normalised differential cross-section of the total $t\bar{t}\gamma$ production and decay measured in fiducial phase space in the single-lepton channel as a function of the $\Delta R(l,j)_{min}$. The last bin of the distribution includes the overflow events.

Absolute differential cross-section of the total $t\bar{t}\gamma$ production and decay measured in fiducial phase space in the dilepton channel as a function of the $|\Delta \eta(l,l)|$ . The last bin of the distribution includes the overflow events.

Normalised differential cross-section of the total $t\bar{t}\gamma$ production and decay measured in fiducial phase space in the dilepton channel as a function of the $|\Delta \eta(l,l)|$ . The last bin of the distribution includes the overflow events.

Absolute differential cross-section of the total $t\bar{t}\gamma$ production and decay measured in fiducial phase space in the dilepton channel as a function of the $\Delta \phi(l,l)$ . The last bin of the distribution includes the overflow events.

Normalised differential cross-section of the total $t\bar{t}\gamma$ production and decay measured in fiducial phase space in the dilepton channel as a function of the $\Delta \phi(l,l)$ . The last bin of the distribution includes the overflow events.

Absolute differential cross-section of the total $t\bar{t}\gamma$ production and decay measured in fiducial phase space in the dilepton channel as a function of the $p_T(l,l)$. The last bin of the distribution includes the overflow events.

Normalised differential cross-section of the total $t\bar{t}\gamma$ production and decay measured in fiducial phase space in the dilepton channel as a function of the $p_T(l,l)$. The last bin of the distribution includes the overflow events.

Absolute differential cross-section of the total $t\bar{t}\gamma$ production and decay measured in fiducial phase space in the dilepton channel as a function of the $\Delta R(\gamma,l_1)$. The last bin of the distribution includes the overflow events.

Normalised differential cross-section of the total $t\bar{t}\gamma$ production and decay measured in fiducial phase space in the dilepton channel as a function of the $\Delta R(\gamma,l_1)$. The last bin of the distribution includes the overflow events.

Absolute differential cross-section of the total $t\bar{t}\gamma$ production and decay measured in fiducial phase space in the dilepton channel as a function of the $\Delta R(\gamma,l_2)$. The last bin of the distribution includes the overflow events.

Normalised differential cross-section of the total $t\bar{t}\gamma$ production and decay measured in fiducial phase space in the dilepton channel as a function of the $\Delta R(\gamma,l_2)$. The last bin of the distribution includes the overflow events.

Absolute $t\bar{t}\gamma$ production differential cross-sections measured in fiducial phase space in the single-lepton channel as a function of the $p_T(\gamma)$. The last bin of the distribution includes the overflow events.

Normalised $t\bar{t}\gamma$ production differential cross-sections measured in fiducial phase space in the single-lepton channel as a function of the $p_T(\gamma)$. The last bin of the distribution includes the overflow events.

Absolute $t\bar{t}\gamma$ production differential cross-sections measured in fiducial phase space in the single-lepton channel as a function of the $|\eta(\gamma)|$.

Normalised $t\bar{t}\gamma$ production differential cross-sections measured in fiducial phase space in the single-lepton channel as a function of the $|\eta(\gamma)|$.

Absolute $t\bar{t}\gamma$ production differential cross-sections measured in fiducial phase space in the single-lepton channel as a function of the $p_T(j_1)$. The last bin of the distribution includes the overflow events.

Normalised $t\bar{t}\gamma$ production differential cross-sections measured in fiducial phase space in the single-lepton channel as a function of the $p_T(j_1)$. The last bin of the distribution includes the overflow events.

Absolute $t\bar{t}\gamma$ production differential cross-sections measured in fiducial phase space in the dilepton channel as a function of the $p_T(\gamma)$. The last bin of the distribution includes the overflow events.

Normalised $t\bar{t}\gamma$ production differential cross-sections measured in fiducial phase space in the dilepton channel as a function of the $p_T(\gamma)$. The last bin of the distribution includes the overflow events.

Absolute $t\bar{t}\gamma$ production differential cross-sections measured in fiducial phase space in the dilepton channel as a function of the $|\eta(\gamma)|$.

Normalised $t\bar{t}\gamma$ production differential cross-sections measured in fiducial phase space in the dilepton channel as a function of the $|\eta(\gamma)|$.

Absolute $t\bar{t}\gamma$ production differential cross-sections measured in fiducial phase space in the dilepton channel as a function of the $p_T(j_1)$. The last bin of the distribution includes the overflow events.

Normalised $t\bar{t}\gamma$ production differential cross-sections measured in fiducial phase space in the dilepton channel as a function of the $p_T(j_1)$. The last bin of the distribution includes the overflow events.

Absolute differential cross-section of the total $t\bar{t}\gamma$ production and decay measured in fiducial phase space in the dilepton channel as a function of the $p_T(\gamma)$. The last bin of the distribution includes the overflow events.

Normalised differential cross-section of the total $t\bar{t}\gamma$ production and decay measured in fiducial phase space in the dilepton channel as a function of the $p_T(\gamma)$. The last bin of the distribution includes the overflow events.

Absolute differential cross-section of the total $t\bar{t}\gamma$ production and decay measured in fiducial phase space in the dilepton channel as a function of the $|\eta(\gamma)|$.

Normalised differential cross-section of the total $t\bar{t}\gamma$ production and decay measured in fiducial phase space in the dilepton channel as a function of the $|\eta(\gamma)|$.

Absolute differential cross-section of the total $t\bar{t}\gamma$ production and decay measured in fiducial phase space in the dilepton channel as a function of the $p_T(j_1)$. The last bin of the distribution includes the overflow events.

Normalised differential cross-section of the total $t\bar{t}\gamma$ production and decay measured in fiducial phase space in the single-lepton channel as a function of the $p_T(j_1)$. The last bin of the distribution includes the overflow events.

Absolute differential cross-section of the total $t\bar{t}\gamma$ production and decay measured in fiducial phase space in the dilepton channel as a function of the $\Delta R(\gamma,b)_{min}$. The last bin of the distribution includes the overflow events.

Normalised differential cross-section of the total $t\bar{t}\gamma$ production and decay measured in fiducial phase space in the dilepton channel as a function of the $\Delta R(\gamma,b)_{min}$. The last bin of the distribution includes the overflow events.

Absolute differential cross-section of the total $t\bar{t}\gamma$ production and decay measured in fiducial phase space in the dilepton channel as a function of the $\Delta R(\gamma,l)_{min}$. The last bin of the distribution includes the overflow events.

Normalised differential cross-section of the total $t\bar{t}\gamma$ production and decay measured in fiducial phase space in the dilepton channel as a function of the $\Delta R(\gamma,l)_{min}$. The last bin of the distribution includes the overflow events.

Absolute differential cross-section of the total $t\bar{t}\gamma$ production and decay measured in fiducial phase space in the dilepton channel as a function of the $\Delta R(l,j)_{min}$. The last bin of the distribution includes the overflow events.

Normalised differential cross-section of the total $t\bar{t}\gamma$ production and decay measured in fiducial phase space in the dilepton channel as a function of the $\Delta R(l,j)_{min}$. The last bin of the distribution includes the overflow events.

Observed and expected 68% and 95% credible intervals on $C/\Lambda^2\ [\text{TeV}^{-2} ]$ for the $C_{tW}$ operators and $C_{tB}$ obtained from the $t\bar{t}\gamma$ measurement, comparing the results obtained from the marginalised quadratic fit (‘marg.’) and the independent quadratic fits (‘indep.’). Also shown are the best-fit values (global mode) for each operator.

Observed and expected 68% and 95% credible intervals on $C/\Lambda^2\ [\text{TeV}^{-2} ]$ for the $C_{tW}$ operators and $C_{tB}$ obtained from a combined $t\bar{t}Z$ and $t\bar{t}\gamma$ measurements comparing the results obtained from the marginalised global quadratic fit (‘marg.’) and the independent quadratic fits (‘indep.’). Also shown are the best-fit values (global mode) for each operator.

Contributions of the systematic uncertainties grouped in categories in each bin of the measurement of the absolute $t\bar{t}\gamma$ production differential cross-sections measured in the combined fiducial phase space of the single-lepton and dilepton channels as a function of the photon $p_T$. The relative uncertainties quoted are obtained by repeating the fit, fixing the set of nuisance parameters of the sources corresponding to each category to their post-fit values, and subtracting in quadrature the resulting uncertainty from the total uncertainty of the nominal fit.

Contributions of the systematic uncertainties grouped in categories in each bin of the measurement of the absolute $t\bar{t}\gamma$ production differential cross-sections measured in the combined fiducial phase space of the single-lepton and dilepton channels as a function of the photon $\eta$. The relative uncertainties quoted are obtained by repeating the fit, fixing the set of nuisance parameters of the sources corresponding to each category to their post-fit values, and subtracting in quadrature the resulting uncertainty from the total uncertainty of the nominal fit.

Post-fit SR yields of the high-mass channel. The upper panel shows the observed number of events, as well the post-fit background predictions in each region. The bottom panel shows the ratio of the observed data and the total background prediction. The shaded areas correspond to the total statistical and systematic uncertainties obtained after the fit and described in Section 6.

Post-fit SR yields of the high-mass channel. The upper panel shows the observed number of events, as well the post-fit background predictions in each region. The bottom panel shows the ratio of the observed data and the total background prediction. The shaded areas correspond to the total statistical and systematic uncertainties obtained after the fit and described in Section 6.

Post-fit SR yields of the high-mass channel. The upper panel shows the observed number of events, as well the post-fit background predictions in each region. The bottom panel shows the ratio of the observed data and the total background prediction. The shaded areas correspond to the total statistical and systematic uncertainties obtained after the fit and described in Section 6.

Post-fit SR yields of the high-mass channel. The upper panel shows the observed number of events, as well the post-fit background predictions in each region. The bottom panel shows the ratio of the observed data and the total background prediction. The shaded areas correspond to the total statistical and systematic uncertainties obtained after the fit and described in Section 6.

Post-fit SR yields of the high-mass channel. The upper panel shows the observed number of events, as well the post-fit background predictions in each region. The bottom panel shows the ratio of the observed data and the total background prediction. The shaded areas correspond to the total statistical and systematic uncertainties obtained after the fit and described in Section 6.

Post-fit SR yields of the high-mass channel. The upper panel shows the observed number of events, as well the post-fit background predictions in each region. The bottom panel shows the ratio of the observed data and the total background prediction. The shaded areas correspond to the total statistical and systematic uncertainties obtained after the fit and described in Section 6.

Post-fit SR yields of the high-mass channel. The upper panel shows the observed number of events, as well the post-fit background predictions in each region. The bottom panel shows the ratio of the observed data and the total background prediction. The shaded areas correspond to the total statistical and systematic uncertainties obtained after the fit and described in Section 6.

Pre-fit data and background (reweighted $2b$) predictions for each $4b$ SR $E_\text{T}^\text{miss}$ and $m_\text{eff}$ bin of the low-mass channel for the 2016 data-taking period. The bottom panel shows the significance of any differences between the observed $4b$ data and the background prediction. The $1\sigma$ and $2\sigma$ bands are shown in green and yellow, respectively. All systematics are included except the background normalization, which is 2.3%.

Pre-fit data and background (reweighted $2b$) predictions for each $4b$ SR $E_\text{T}^\text{miss}$ and $m_\text{eff}$ bin of the low-mass channel for the 2017 data-taking period. The bottom panel shows the significance of any differences between the observed $4b$ data and the background prediction. The $1\sigma$ and $2\sigma$ bands are shown in green and yellow, respectively. All systematics are included except the background normalization, which is 3.7%.

Pre-fit data and background (reweighted $2b$) predictions for each $4b$ SR $E_\text{T}^\text{miss}$ and $m_\text{eff}$ bin of the low-mass channel for the 2018 data-taking period. The bottom panel shows the significance of any differences between the observed $4b$ data and the background prediction. The $1\sigma$ and $2\sigma$ bands are shown in green and yellow, respectively. All systematics are included except the background normalization, which is 1.8%.

Exclusion limits of the low-mass and high-mass channels. The low-mass channel is used for $m_{\tilde{H}}<250$ GeV while the high-mass channel is used for $m_{\tilde{H}}\ge250$ GeV. The plot shows the observed (solid) and expected (dashed) 95% CL upper limits on the cross section of higgsino pair production, assuming a higgsino decay branching ratio of $\mathcal{B}(\tilde{H}\rightarrow h + \tilde{G})=100\%$. The theory cross section and its uncertainty are shown by the solid red line and red shading. Results from a previous ATLAS search using 24.3-36.1 fb$^{-1}$ [13] are shown by the solid (observed) and dashed (expected) blue lines. The bottom panel shows the ratio of the limits to the theory cross section. The phase space above the lines is excluded.

Exclusion limits of the low-mass and high-mass channels. The low-mass channel is used for $m_{\tilde{H}}<250$ GeV while the high-mass channel is used for $m_{\tilde{H}}\ge250$ GeV. The plot shows the observed (solid) and expected (dashed) 95% CL upper limits on the cross section of higgsino pair production, assuming a higgsino decay branching ratio of $\mathcal{B}(\tilde{H}\rightarrow h + \tilde{G})=100\%$. The theory cross section and its uncertainty are shown by the solid red line and red shading. Results from a previous ATLAS search using 24.3-36.1 fb$^{-1}$ [13] are shown by the solid (observed) and dashed (expected) blue lines. The bottom panel shows the ratio of the limits to the theory cross section. The phase space above the lines is excluded.

Exclusion limits of the low-mass and high-mass channels. The low-mass channel is used for $m_{\tilde{H}}<250$ GeV while the high-mass channel is used for $m_{\tilde{H}}\ge250$ GeV. The plot shows the observed (solid) and expected (dashed) 95% CL upper limits on the cross section of higgsino pair production, assuming a higgsino decay branching ratio of $\mathcal{B}(\tilde{H}\rightarrow h + \tilde{G})=100\%$. The theory cross section and its uncertainty are shown by the solid red line and red shading. Results from a previous ATLAS search using 24.3-36.1 fb$^{-1}$ [13] are shown by the solid (observed) and dashed (expected) blue lines. The bottom panel shows the ratio of the limits to the theory cross section. The phase space above the lines is excluded.

Exclusion limits of the low-mass and high-mass channels. The low-mass channel is used for $m_{\tilde{H}}<250$ GeV while the high-mass channel is used for $m_{\tilde{H}}\ge250$ GeV. The plot shows the observed (solid) and expected (dashed) 95% CL upper limits on the cross section of higgsino pair production, assuming a higgsino decay branching ratio of $\mathcal{B}(\tilde{H}\rightarrow h + \tilde{G})=100\%$. The theory cross section and its uncertainty are shown by the solid red line and red shading. Results from a previous ATLAS search using 24.3-36.1 fb$^{-1}$ [13] are shown by the solid (observed) and dashed (expected) blue lines. The bottom panel shows the ratio of the limits to the theory cross section. The phase space above the lines is excluded.

Exclusion limits of the low-mass and high-mass channels. The low-mass channel is used for $m_{\tilde{H}}<250$ GeV while the high-mass channel is used for $m_{\tilde{H}}\ge250$ GeV. The plot shows the observed (solid) and expected (dashed) 95% CL upper limits on the cross section of higgsino pair production, assuming a higgsino decay branching ratio of $\mathcal{B}(\tilde{H}\rightarrow h + \tilde{G})=100\%$. The theory cross section and its uncertainty are shown by the solid red line and red shading. Results from a previous ATLAS search using 24.3-36.1 fb$^{-1}$ [13] are shown by the solid (observed) and dashed (expected) blue lines. The bottom panel shows the ratio of the limits to the theory cross section. The phase space above the lines is excluded.

Exclusion limits of the low-mass and high-mass channels. The low-mass channel is used for $m_{\tilde{H}}<250$ GeV while the high-mass channel is used for $m_{\tilde{H}}\ge250$ GeV. The plot shows the observed (solid) and expected (dashed) 95% CL upper limits on the cross section of higgsino pair production, assuming a higgsino decay branching ratio of $\mathcal{B}(\tilde{H}\rightarrow h + \tilde{G})=100\%$. The theory cross section and its uncertainty are shown by the solid red line and red shading. Results from a previous ATLAS search using 24.3-36.1 fb$^{-1}$ [13] are shown by the solid (observed) and dashed (expected) blue lines. The bottom panel shows the ratio of the limits to the theory cross section. The phase space above the lines is excluded.

Exclusion limits of the low-mass and high-mass channels. The low-mass channel is used for $m_{\tilde{H}}<250$ GeV while the high-mass channel is used for $m_{\tilde{H}}\ge250$ GeV. The plot shows the observed (solid) and expected (dashed) 95% CL upper limits on the cross section of higgsino pair production, assuming a higgsino decay branching ratio of $\mathcal{B}(\tilde{H}\rightarrow h + \tilde{G})=100\%$. The theory cross section and its uncertainty are shown by the solid red line and red shading. Results from a previous ATLAS search using 24.3-36.1 fb$^{-1}$ [13] are shown by the solid (observed) and dashed (expected) blue lines. The bottom panel shows the ratio of the limits to the theory cross section. The phase space above the lines is excluded.

Exclusion limits of the low-mass and high-mass channels. The low-mass channel is used for $m_{\tilde{H}}<250$ GeV while the high-mass channel is used for $m_{\tilde{H}}\ge250$ GeV. The plot shows the observed (solid) and expected (dashed) 95% CL upper limits on the cross section of higgsino pair production, assuming a higgsino decay branching ratio of $\mathcal{B}(\tilde{H}\rightarrow h + \tilde{G})=100\%$. The theory cross section and its uncertainty are shown by the solid red line and red shading. Results from a previous ATLAS search using 24.3-36.1 fb$^{-1}$ [13] are shown by the solid (observed) and dashed (expected) blue lines. The bottom panel shows the ratio of the limits to the theory cross section. The phase space above the lines is excluded.

Exclusion limits of the low-mass and high-mass channels. The low-mass channel is used for $m_{\tilde{H}}<250$ GeV while the high-mass channel is used for $m_{\tilde{H}}\ge250$ GeV. The plot shows the observed (solid) and expected (dashed) 95% CL upper limits on the cross section of higgsino pair production, assuming a higgsino decay branching ratio of $\mathcal{B}(\tilde{H}\rightarrow h + \tilde{G})=100\%$. The theory cross section and its uncertainty are shown by the solid red line and red shading. Results from a previous ATLAS search using 24.3-36.1 fb$^{-1}$ [13] are shown by the solid (observed) and dashed (expected) blue lines. The bottom panel shows the ratio of the limits to the theory cross section. The phase space above the lines is excluded.

Exclusion limits of the low-mass and high-mass channels. The low-mass channel is used for $m_{\tilde{H}}<250$ GeV while the high-mass channel is used for $m_{\tilde{H}}\ge250$ GeV. The plot shows the observed (solid) and expected (dashed) 95% CL upper limits on the cross section of higgsino pair production, assuming a higgsino decay branching ratio of $\mathcal{B}(\tilde{H}\rightarrow h + \tilde{G})=100\%$. The theory cross section and its uncertainty are shown by the solid red line and red shading. Results from a previous ATLAS search using 24.3-36.1 fb$^{-1}$ [13] are shown by the solid (observed) and dashed (expected) blue lines. The bottom panel shows the ratio of the limits to the theory cross section. The phase space above the lines is excluded.

Exclusion limits of the low-mass and high-mass channels. The low-mass channel is used for $m_{\tilde{H}}<250$ GeV while the high-mass channel is used for $m_{\tilde{H}}\ge250$ GeV. The plot shows the 95% CL observed (solid) and expected (dashed) upper limits on $\mathcal{B}(\tilde{H}\rightarrow h + \tilde{G})$, assuming the theory cross section for higgsino pair production. The higgsinos are assumed to decay as $\tilde{H}\rightarrow h + \tilde{G}$ or $\tilde{H}\rightarrow Z + \tilde{G}$. The phase space above the lines is excluded.

Exclusion limits of the low-mass and high-mass channels. The low-mass channel is used for $m_{\tilde{H}}<250$ GeV while the high-mass channel is used for $m_{\tilde{H}}\ge250$ GeV. The plot shows the 95% CL observed (solid) and expected (dashed) upper limits on $\mathcal{B}(\tilde{H}\rightarrow h + \tilde{G})$, assuming the theory cross section for higgsino pair production. The higgsinos are assumed to decay as $\tilde{H}\rightarrow h + \tilde{G}$ or $\tilde{H}\rightarrow Z + \tilde{G}$. The phase space above the lines is excluded.

Exclusion limits of the low-mass and high-mass channels. The low-mass channel is used for $m_{\tilde{H}}<250$ GeV while the high-mass channel is used for $m_{\tilde{H}}\ge250$ GeV. The plot shows the 95% CL observed (solid) and expected (dashed) upper limits on $\mathcal{B}(\tilde{H}\rightarrow h + \tilde{G})$, assuming the theory cross section for higgsino pair production. The higgsinos are assumed to decay as $\tilde{H}\rightarrow h + \tilde{G}$ or $\tilde{H}\rightarrow Z + \tilde{G}$. The phase space above the lines is excluded.

Exclusion limits of the low-mass and high-mass channels. The low-mass channel is used for $m_{\tilde{H}}<250$ GeV while the high-mass channel is used for $m_{\tilde{H}}\ge250$ GeV. The plot shows the 95% CL observed (solid) and expected (dashed) upper limits on $\mathcal{B}(\tilde{H}\rightarrow h + \tilde{G})$, assuming the theory cross section for higgsino pair production. The higgsinos are assumed to decay as $\tilde{H}\rightarrow h + \tilde{G}$ or $\tilde{H}\rightarrow Z + \tilde{G}$. The phase space above the lines is excluded.

Exclusion limits of the low-mass and high-mass channels. The low-mass channel is used for $m_{\tilde{H}}<250$ GeV while the high-mass channel is used for $m_{\tilde{H}}\ge250$ GeV. The plot shows the 95% CL observed (solid) and expected (dashed) upper limits on $\mathcal{B}(\tilde{H}\rightarrow h + \tilde{G})$, assuming the theory cross section for higgsino pair production. The higgsinos are assumed to decay as $\tilde{H}\rightarrow h + \tilde{G}$ or $\tilde{H}\rightarrow Z + \tilde{G}$. The phase space above the lines is excluded.

Exclusion limits of the low-mass channel. The plot shows the observed (solid) and expected (dashed) 95% CL upper limits on the cross section of higgsino pair production, assuming a higgsino decay branching ratio of $\mathcal{B}(\tilde{H}\rightarrow h + \tilde{G})=100\%$. The theory cross section and its uncertainty are shown by the solid red line and red shading. The bottom panel shows the ratio of the limits to the theory cross section. The phase space above the lines is excluded.

Exclusion limits of the low-mass channel. The plot shows the observed (solid) and expected (dashed) 95% CL upper limits on the cross section of higgsino pair production, assuming a higgsino decay branching ratio of $\mathcal{B}(\tilde{H}\rightarrow h + \tilde{G})=100\%$. The theory cross section and its uncertainty are shown by the solid red line and red shading. The bottom panel shows the ratio of the limits to the theory cross section. The phase space above the lines is excluded.

Exclusion limits of the high-mass channel. The plot shows the observed (solid) and expected (dashed) 95% CL upper limits on the cross section of higgsino pair production, assuming a higgsino decay branching ratio of $\mathcal{B}(\tilde{H}\rightarrow h + \tilde{G})=100\%$. The theory cross section and its uncertainty are shown by the solid red line and red shading. The bottom panel shows the ratio of the limits to the theory cross section. The phase space above the lines is excluded.

Exclusion limits of the high-mass channel. The plot shows the observed (solid) and expected (dashed) 95% CL upper limits on the cross section of higgsino pair production, assuming a higgsino decay branching ratio of $\mathcal{B}(\tilde{H}\rightarrow h + \tilde{G})=100\%$. The theory cross section and its uncertainty are shown by the solid red line and red shading. The bottom panel shows the ratio of the limits to the theory cross section. The phase space above the lines is excluded.

Results of the background-only fit in the low-mass channel discovery region SR_LM_150. Both pre-fit and post-fit values are shown.

Results of the background-only fit in the low-mass channel discovery region SR_LM_300. Both pre-fit and post-fit values are shown.

The experimental efficiency of the low-mass channel for the exclusion and discovery signal regions as a function of higgsino mass. The experimental efficiency is defined as the number of events passing the detector-level event selections divided by the number of events passing the event selections for a perfect detector. The denominator is obtained by implementing particle-level event selections that emulate the detector-level selections. This treats the lack of availability of $b$-jet triggers as an inefficiency.

The particle-level acceptance for the low-mass exclusion and discovery signal regions, shown as a function of higgsino mass. The acceptance is defined as the fraction of signal events passing the particle-level event selection that emulates the detector-level selection. The acceptance calculation considers only those signal events where both higgsinos decay to Higgs bosons.

The experimental efficiency of the high-mass channel discovery regions as a function of higgsino mass. For each higgsino mass, the efficiency is shown for the SR-1 region corresponding to the mass. For masses above 1100 GeV, SR-1-1100 is used. The experimental efficiency is defined as the number of events passing the detector-level event selections divided by the number of events passing the event selections for a perfect detector. The denominator is obtained by implementing particle-level event selections that emulate the detector-level selections. The efficiency calculation considers only those signal events where both higgsinos decay to Higgs bosons.

The particle-level acceptance for the high-mass signal regions, shown as a function of higgsino mass. For each higgsino mass, the acceptance is shown for the SR-1 region corresponding to the mass. For masses above 1100 GeV, SR-1-1100 is used. The acceptance is defined as the fraction of signal events passing the particle-level event selection that emulates the detector-level selection. The acceptance calculation considers only those signal events where both higgsinos decay to Higgs bosons.

Cutflow for the low-mass channel for a representative 130 GeV signal. The preselection requires 4 or more jets with $p_\text{T}>25$ GeV and 2 or more $b$-jets with $p_\text{T}>25$ GeV. The $b$-jet cut requires 4 or more $b$-jets with $p_\text{T}>40$ GeV. As the samples are generated with $\mathcal{B}(\tilde{H}\rightarrow h\tilde{G})$=50%, $\mathcal{B}(\tilde{H}\rightarrow Z\tilde{G})$=50% to allow for both decays to be studied, the $hh$ events selection is used to select the events where each of the higgsinos decays to a Higgs boson. Expected yields are normalized to a luminosity of 139 fb$^{-1}$, with the availability of $b$-jet triggers lowering the luminosity to 126 fb$^{-1}$. All selections are cumulative, with the exception of the SR cuts, which are each applied separately.

Cutflow for the low-mass channel for a representative 150 GeV signal. The preselection requires 4 or more jets with $p_\text{T}>25$ GeV and 2 or more $b$-jets with $p_\text{T}>25$ GeV. The $b$-jet cut requires 4 or more $b$-jets with $p_\text{T}>40$ GeV. As the samples are generated with $\mathcal{B}(\tilde{H}\rightarrow h\tilde{G})$=50%, $\mathcal{B}(\tilde{H}\rightarrow Z\tilde{G})$=50% to allow for both decays to be studied, the $hh$ events selection is used to select the events where each of the higgsinos decays to a Higgs boson. Expected yields are normalized to a luminosity of 139 fb$^{-1}$, with the availability of $b$-jet triggers lowering the luminosity to 126 fb$^{-1}$. All selections are cumulative, with the exception of the SR cuts, which are each applied separately.

Cutflow for the low-mass channel for a representative 200 GeV signal. The preselection requires 4 or more jets with $p_\text{T}>25$ GeV and 2 or more $b$-jets with $p_\text{T}>25$ GeV. The $b$-jet cut requires 4 or more $b$-jets with $p_\text{T}>40$ GeV. As the samples are generated with $\mathcal{B}(\tilde{H}\rightarrow h\tilde{G})$=50%, $\mathcal{B}(\tilde{H}\rightarrow Z\tilde{G})$=50% to allow for both decays to be studied, the $hh$ events selection is used to select the events where each of the higgsinos decays to a Higgs boson. Expected yields are normalized to a luminosity of 139 fb$^{-1}$, with the availability of $b$-jet triggers lowering the luminosity to 126 fb$^{-1}$. All selections are cumulative, with the exception of the SR cuts, which are each applied separately.

Cutflow for the low-mass channel for a representative 250 GeV signal. The preselection requires 4 or more jets with $p_\text{T}>25$ GeV and 2 or more $b$-jets with $p_\text{T}>25$ GeV. The $b$-jet cut requires 4 or more $b$-jets with $p_\text{T}>40$ GeV. As the samples are generated with $\mathcal{B}(\tilde{H}\rightarrow h\tilde{G})$=50%, $\mathcal{B}(\tilde{H}\rightarrow Z\tilde{G})$=50% to allow for both decays to be studied, the $hh$ events selection is used to select the events where each of the higgsinos decays to a Higgs boson. Expected yields are normalized to a luminosity of 139 fb$^{-1}$, with the availability of $b$-jet triggers lowering the luminosity to 126 fb$^{-1}$. All selections are cumulative, with the exception of the SR cuts, which are each applied separately.

Cutflow for the low-mass channel for a representative 300 GeV signal. The preselection requires 4 or more jets with $p_\text{T}>25$ GeV and 2 or more $b$-jets with $p_\text{T}>25$ GeV. The $b$-jet cut requires 4 or more $b$-jets with $p_\text{T}>40$ GeV. As the samples are generated with $\mathcal{B}(\tilde{H}\rightarrow h\tilde{G})$=50%, $\mathcal{B}(\tilde{H}\rightarrow Z\tilde{G})$=50% to allow for both decays to be studied, the $hh$ events selection is used to select the events where each of the higgsinos decays to a Higgs boson. Expected yields are normalized to a luminosity of 139 fb$^{-1}$, with the availability of $b$-jet triggers lowering the luminosity to 126 fb$^{-1}$. All selections are cumulative, with the exception of the SR cuts, which are each applied separately.

Cutflow for the low-mass channel for a representative 400 GeV signal. The preselection requires 4 or more jets with $p_\text{T}>25$ GeV and 2 or more $b$-jets with $p_\text{T}>25$ GeV. The $b$-jet cut requires 4 or more $b$-jets with $p_\text{T}>40$ GeV. As the samples are generated with $\mathcal{B}(\tilde{H}\rightarrow h\tilde{G})$=50%, $\mathcal{B}(\tilde{H}\rightarrow Z\tilde{G})$=50% to allow for both decays to be studied, the $hh$ events selection is used to select the events where each of the higgsinos decays to a Higgs boson. Expected yields are normalized to a luminosity of 139 fb$^{-1}$, with the availability of $b$-jet triggers lowering the luminosity to 126 fb$^{-1}$. All selections are cumulative, with the exception of the SR cuts, which are each applied separately.

Cutflow for the low-mass channel for a representative 500 GeV signal. The preselection requires 4 or more jets with $p_\text{T}>25$ GeV and 2 or more $b$-jets with $p_\text{T}>25$ GeV. The $b$-jet cut requires 4 or more $b$-jets with $p_\text{T}>40$ GeV. As the samples are generated with $\mathcal{B}(\tilde{H}\rightarrow h\tilde{G})$=50%, $\mathcal{B}(\tilde{H}\rightarrow Z\tilde{G})$=50% to allow for both decays to be studied, the $hh$ events selection is used to select the events where each of the higgsinos decays to a Higgs boson. Expected yields are normalized to a luminosity of 139 fb$^{-1}$, with the availability of $b$-jet triggers lowering the luminosity to 126 fb$^{-1}$. All selections are cumulative, with the exception of the SR cuts, which are each applied separately.

Cutflow for the low-mass channel for a representative 600 GeV signal. The preselection requires 4 or more jets with $p_\text{T}>25$ GeV and 2 or more $b$-jets with $p_\text{T}>25$ GeV. The $b$-jet cut requires 4 or more $b$-jets with $p_\text{T}>40$ GeV. As the samples are generated with $\mathcal{B}(\tilde{H}\rightarrow h\tilde{G})$=50%, $\mathcal{B}(\tilde{H}\rightarrow Z\tilde{G})$=50% to allow for both decays to be studied, the $hh$ events selection is used to select the events where each of the higgsinos decays to a Higgs boson. Expected yields are normalized to a luminosity of 139 fb$^{-1}$, with the availability of $b$-jet triggers lowering the luminosity to 126 fb$^{-1}$. All selections are cumulative, with the exception of the SR cuts, which are each applied separately.

Cutflow for the low-mass channel for a representative 700 GeV signal. The preselection requires 4 or more jets with $p_\text{T}>25$ GeV and 2 or more $b$-jets with $p_\text{T}>25$ GeV. The $b$-jet cut requires 4 or more $b$-jets with $p_\text{T}>40$ GeV. As the samples are generated with $\mathcal{B}(\tilde{H}\rightarrow h\tilde{G})$=50%, $\mathcal{B}(\tilde{H}\rightarrow Z\tilde{G})$=50% to allow for both decays to be studied, the $hh$ events selection is used to select the events where each of the higgsinos decays to a Higgs boson. Expected yields are normalized to a luminosity of 139 fb$^{-1}$, with the availability of $b$-jet triggers lowering the luminosity to 126 fb$^{-1}$. All selections are cumulative, with the exception of the SR cuts, which are each applied separately.

Cutflow for the low-mass channel for a representative 800 GeV signal. The preselection requires 4 or more jets with $p_\text{T}>25$ GeV and 2 or more $b$-jets with $p_\text{T}>25$ GeV. The $b$-jet cut requires 4 or more $b$-jets with $p_\text{T}>40$ GeV. As the samples are generated with $\mathcal{B}(\tilde{H}\rightarrow h\tilde{G})$=50%, $\mathcal{B}(\tilde{H}\rightarrow Z\tilde{G})$=50% to allow for both decays to be studied, the $hh$ events selection is used to select the events where each of the higgsinos decays to a Higgs boson. Expected yields are normalized to a luminosity of 139 fb$^{-1}$, with the availability of $b$-jet triggers lowering the luminosity to 126 fb$^{-1}$. All selections are cumulative, with the exception of the SR cuts, which are each applied separately.

Cutflow for the low-mass channel for a representative 900 GeV signal. The preselection requires 4 or more jets with $p_\text{T}>25$ GeV and 2 or more $b$-jets with $p_\text{T}>25$ GeV. The $b$-jet cut requires 4 or more $b$-jets with $p_\text{T}>40$ GeV. As the samples are generated with $\mathcal{B}(\tilde{H}\rightarrow h\tilde{G})$=50%, $\mathcal{B}(\tilde{H}\rightarrow Z\tilde{G})$=50% to allow for both decays to be studied, the $hh$ events selection is used to select the events where each of the higgsinos decays to a Higgs boson. Expected yields are normalized to a luminosity of 139 fb$^{-1}$, with the availability of $b$-jet triggers lowering the luminosity to 126 fb$^{-1}$. All selections are cumulative, with the exception of the SR cuts, which are each applied separately.

Cutflow for the low-mass channel for a representative 1000 GeV signal. The preselection requires 4 or more jets with $p_\text{T}>25$ GeV and 2 or more $b$-jets with $p_\text{T}>25$ GeV. The $b$-jet cut requires 4 or more $b$-jets with $p_\text{T}>40$ GeV. As the samples are generated with $\mathcal{B}(\tilde{H}\rightarrow h\tilde{G})$=50%, $\mathcal{B}(\tilde{H}\rightarrow Z\tilde{G})$=50% to allow for both decays to be studied, the $hh$ events selection is used to select the events where each of the higgsinos decays to a Higgs boson. Expected yields are normalized to a luminosity of 139 fb$^{-1}$, with the availability of $b$-jet triggers lowering the luminosity to 126 fb$^{-1}$. All selections are cumulative, with the exception of the SR cuts, which are each applied separately.

Cutflow for the low-mass channel for a representative 1100 GeV signal. The preselection requires 4 or more jets with $p_\text{T}>25$ GeV and 2 or more $b$-jets with $p_\text{T}>25$ GeV. The $b$-jet cut requires 4 or more $b$-jets with $p_\text{T}>40$ GeV. As the samples are generated with $\mathcal{B}(\tilde{H}\rightarrow h\tilde{G})$=50%, $\mathcal{B}(\tilde{H}\rightarrow Z\tilde{G})$=50% to allow for both decays to be studied, the $hh$ events selection is used to select the events where each of the higgsinos decays to a Higgs boson. Expected yields are normalized to a luminosity of 139 fb$^{-1}$, with the availability of $b$-jet triggers lowering the luminosity to 126 fb$^{-1}$. All selections are cumulative, with the exception of the SR cuts, which are each applied separately.

Cutflow for the high-mass channel for a representative 200 GeV signal. The preselection requires 4 or more jets with $p_\text{T}>25$ GeV and 2 or more $b$-jets with $p_\text{T}>25$ GeV. As the samples are generated with $\mathcal{B}(\tilde{H}\rightarrow h\tilde{G})$=50%, $\mathcal{B}(\tilde{H}\rightarrow Z\tilde{G})$=50% to allow for both decays to be studied, the $hh$ events selection is used to select the events where each of the higgsinos decays to a Higgs boson. Expected yields are normalized to a luminosity of 139 fb$^{-1}$. All selections are cumulative, with the exception of the SR cuts, which are each applied separately.

Cutflow for the high-mass channel for a representative 250 GeV signal. The preselection requires 4 or more jets with $p_\text{T}>25$ GeV and 2 or more $b$-jets with $p_\text{T}>25$ GeV. As the samples are generated with $\mathcal{B}(\tilde{H}\rightarrow h\tilde{G})$=50%, $\mathcal{B}(\tilde{H}\rightarrow Z\tilde{G})$=50% to allow for both decays to be studied, the $hh$ events selection is used to select the events where each of the higgsinos decays to a Higgs boson. Expected yields are normalized to a luminosity of 139 fb$^{-1}$. All selections are cumulative, with the exception of the SR cuts, which are each applied separately.

Cutflow for the high-mass channel for a representative 300 GeV signal. The preselection requires 4 or more jets with $p_\text{T}>25$ GeV and 2 or more $b$-jets with $p_\text{T}>25$ GeV. As the samples are generated with $\mathcal{B}(\tilde{H}\rightarrow h\tilde{G})$=50%, $\mathcal{B}(\tilde{H}\rightarrow Z\tilde{G})$=50% to allow for both decays to be studied, the $hh$ events selection is used to select the events where each of the higgsinos decays to a Higgs boson. Expected yields are normalized to a luminosity of 139 fb$^{-1}$. All selections are cumulative, with the exception of the SR cuts, which are each applied separately.

Cutflow for the high-mass channel for a representative 400 GeV signal. The preselection requires 4 or more jets with $p_\text{T}>25$ GeV and 2 or more $b$-jets with $p_\text{T}>25$ GeV. As the samples are generated with $\mathcal{B}(\tilde{H}\rightarrow h\tilde{G})$=50%, $\mathcal{B}(\tilde{H}\rightarrow Z\tilde{G})$=50% to allow for both decays to be studied, the $hh$ events selection is used to select the events where each of the higgsinos decays to a Higgs boson. Expected yields are normalized to a luminosity of 139 fb$^{-1}$. All selections are cumulative, with the exception of the SR cuts, which are each applied separately.

Cutflow for the high-mass channel for a representative 500 GeV signal. The preselection requires 4 or more jets with $p_\text{T}>25$ GeV and 2 or more $b$-jets with $p_\text{T}>25$ GeV. As the samples are generated with $\mathcal{B}(\tilde{H}\rightarrow h\tilde{G})$=50%, $\mathcal{B}(\tilde{H}\rightarrow Z\tilde{G})$=50% to allow for both decays to be studied, the $hh$ events selection is used to select the events where each of the higgsinos decays to a Higgs boson. Expected yields are normalized to a luminosity of 139 fb$^{-1}$. All selections are cumulative, with the exception of the SR cuts, which are each applied separately.

Cutflow for the high-mass channel for a representative 600 GeV signal. The preselection requires 4 or more jets with $p_\text{T}>25$ GeV and 2 or more $b$-jets with $p_\text{T}>25$ GeV. As the samples are generated with $\mathcal{B}(\tilde{H}\rightarrow h\tilde{G})$=50%, $\mathcal{B}(\tilde{H}\rightarrow Z\tilde{G})$=50% to allow for both decays to be studied, the $hh$ events selection is used to select the events where each of the higgsinos decays to a Higgs boson. Expected yields are normalized to a luminosity of 139 fb$^{-1}$. All selections are cumulative, with the exception of the SR cuts, which are each applied separately.

Cutflow for the high-mass channel for a representative 700 GeV signal. The preselection requires 4 or more jets with $p_\text{T}>25$ GeV and 2 or more $b$-jets with $p_\text{T}>25$ GeV. As the samples are generated with $\mathcal{B}(\tilde{H}\rightarrow h\tilde{G})$=50%, $\mathcal{B}(\tilde{H}\rightarrow Z\tilde{G})$=50% to allow for both decays to be studied, the $hh$ events selection is used to select the events where each of the higgsinos decays to a Higgs boson. Expected yields are normalized to a luminosity of 139 fb$^{-1}$. All selections are cumulative, with the exception of the SR cuts, which are each applied separately.

Cutflow for the high-mass channel for a representative 800 GeV signal. The preselection requires 4 or more jets with $p_\text{T}>25$ GeV and 2 or more $b$-jets with $p_\text{T}>25$ GeV. As the samples are generated with $\mathcal{B}(\tilde{H}\rightarrow h\tilde{G})$=50%, $\mathcal{B}(\tilde{H}\rightarrow Z\tilde{G})$=50% to allow for both decays to be studied, the $hh$ events selection is used to select the events where each of the higgsinos decays to a Higgs boson. Expected yields are normalized to a luminosity of 139 fb$^{-1}$. All selections are cumulative, with the exception of the SR cuts, which are each applied separately.

Cutflow for the high-mass channel for a representative 900 GeV signal. The preselection requires 4 or more jets with $p_\text{T}>25$ GeV and 2 or more $b$-jets with $p_\text{T}>25$ GeV. As the samples are generated with $\mathcal{B}(\tilde{H}\rightarrow h\tilde{G})$=50%, $\mathcal{B}(\tilde{H}\rightarrow Z\tilde{G})$=50% to allow for both decays to be studied, the $hh$ events selection is used to select the events where each of the higgsinos decays to a Higgs boson. Expected yields are normalized to a luminosity of 139 fb$^{-1}$. All selections are cumulative, with the exception of the SR cuts, which are each applied separately.

Cutflow for the high-mass channel for a representative 1000 GeV signal. The preselection requires 4 or more jets with $p_\text{T}>25$ GeV and 2 or more $b$-jets with $p_\text{T}>25$ GeV. As the samples are generated with $\mathcal{B}(\tilde{H}\rightarrow h\tilde{G})$=50%, $\mathcal{B}(\tilde{H}\rightarrow Z\tilde{G})$=50% to allow for both decays to be studied, the $hh$ events selection is used to select the events where each of the higgsinos decays to a Higgs boson. Expected yields are normalized to a luminosity of 139 fb$^{-1}$. All selections are cumulative, with the exception of the SR cuts, which are each applied separately.

Cutflow for the high-mass channel for a representative 1100 GeV signal. The preselection requires 4 or more jets with $p_\text{T}>25$ GeV and 2 or more $b$-jets with $p_\text{T}>25$ GeV. As the samples are generated with $\mathcal{B}(\tilde{H}\rightarrow h\tilde{G})$=50%, $\mathcal{B}(\tilde{H}\rightarrow Z\tilde{G})$=50% to allow for both decays to be studied, the $hh$ events selection is used to select the events where each of the higgsinos decays to a Higgs boson. Expected yields are normalized to a luminosity of 139 fb$^{-1}$. All selections are cumulative, with the exception of the SR cuts, which are each applied separately.

Cutflow for the high-mass channel for a representative 1200 GeV signal. The preselection requires 4 or more jets with $p_\text{T}>25$ GeV and 2 or more $b$-jets with $p_\text{T}>25$ GeV. As the samples are generated with $\mathcal{B}(\tilde{H}\rightarrow h\tilde{G})$=50%, $\mathcal{B}(\tilde{H}\rightarrow Z\tilde{G})$=50% to allow for both decays to be studied, the $hh$ events selection is used to select the events where each of the higgsinos decays to a Higgs boson. Expected yields are normalized to a luminosity of 139 fb$^{-1}$. All selections are cumulative, with the exception of the SR cuts, which are each applied separately.

Cutflow for the high-mass channel for a representative 1300 GeV signal. The preselection requires 4 or more jets with $p_\text{T}>25$ GeV and 2 or more $b$-jets with $p_\text{T}>25$ GeV. As the samples are generated with $\mathcal{B}(\tilde{H}\rightarrow h\tilde{G})$=50%, $\mathcal{B}(\tilde{H}\rightarrow Z\tilde{G})$=50% to allow for both decays to be studied, the $hh$ events selection is used to select the events where each of the higgsinos decays to a Higgs boson. Expected yields are normalized to a luminosity of 139 fb$^{-1}$. All selections are cumulative, with the exception of the SR cuts, which are each applied separately.

Cutflow for the high-mass channel for a representative 1400 GeV signal. The preselection requires 4 or more jets with $p_\text{T}>25$ GeV and 2 or more $b$-jets with $p_\text{T}>25$ GeV. As the samples are generated with $\mathcal{B}(\tilde{H}\rightarrow h\tilde{G})$=50%, $\mathcal{B}(\tilde{H}\rightarrow Z\tilde{G})$=50% to allow for both decays to be studied, the $hh$ events selection is used to select the events where each of the higgsinos decays to a Higgs boson. Expected yields are normalized to a luminosity of 139 fb$^{-1}$. All selections are cumulative, with the exception of the SR cuts, which are each applied separately.

Cutflow for the high-mass channel for a representative 1500 GeV signal. The preselection requires 4 or more jets with $p_\text{T}>25$ GeV and 2 or more $b$-jets with $p_\text{T}>25$ GeV. As the samples are generated with $\mathcal{B}(\tilde{H}\rightarrow h\tilde{G})$=50%, $\mathcal{B}(\tilde{H}\rightarrow Z\tilde{G})$=50% to allow for both decays to be studied, the $hh$ events selection is used to select the events where each of the higgsinos decays to a Higgs boson. Expected yields are normalized to a luminosity of 139 fb$^{-1}$. All selections are cumulative, with the exception of the SR cuts, which are each applied separately.

Expected (dashed line) value of $-2\ln\Lambda$ as a function of $\kappa_{2V}$, 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 observed contours at 68% CL (solid line) and 95% CL (dashed line) in the $c_{gghh}$ versus $c_{hhh}$ HEFT parameter space, with the remaining coefficient fixed to its SM value. 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) 95% CL (yellow shaded region) in the $c_{gghh}$ versus $c_{hhh}$ HEFT parameter space, with the remaining coefficient fixed to its SM value. In the plot, the SM prediction is indicated by the star.

Likelihood observed contours at 68% CL (solid line) and 95% CL (dashed line) in the $c_{tthh}$ versus $c_{hhh}$ HEFT parameter space, with the remaining coefficient fixed to its SM value. The best-fit value, denoted by a cross in the plot, corresponds to the very first entry in the table.

Likelihood exptected contours at 68% CL (teal shaded region) and 95% CL (yellow shaded region) in the $c_{tthh}$ versus $c_{hhh}$ HEFT parameter space, with the remaining coefficient fixed to its SM value. In the plot, the SM prediction is indicated by the star.

The observed (filled circles) and expected (hollow circles) 95% CL upper limits on the $HH$ ggF production cross-section in the Standard Model and for seven HEFT benchmark points. The teal and yellow colored bands indicate the $\pm 1\sigma$ and $\pm 2\sigma$ variations on the expected limit due to statistical and systematic uncertainties. The predicted cross-sections of each of the models under consideration are shown by the red crosses. The contribution from VBF production to the total $HH$ production cross-section is neglected.

Likelihood observed contours at 68% CL (solid line) and 95% CL (dashed line) in the $c_{{H}\boxed{}}$ versus $c_{H}$ SMEFT parameter space, with the remaining coefficient fixed to its SM value. 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) 95% CL (yellow shaded region) in the $c_{{H}\boxed{}}$ versus $c_{H}$ SMEFT parameter space, with the remaining coefficient fixed to its SM value. In the plot, the SM prediction is indicated by the star.

The acceptance times efficiency for the signal ggF $HH$ process as a function of the coupling modifier $\kappa_{\lambda}$ in each analysis category. The other coupling modifiers affecting ggF $HH$ production are fixed to their SM predictions. The bands indicate the simulation statistical uncertainty.

The acceptance times efficiency for the signal VBF $HH$ process as a function of the coupling modifier $\kappa_{\lambda}$ in each analysis category. When varying $\kappa_{\lambda}$, the other coupling modifiers affecting the VBF $HH$ production mode are set to their SM predictions. The bands indicate the simulation statistical uncertainty.

The acceptance times efficiency for the signal VBF $HH$ process as a function of the coupling modifier $\kappa_{2V}$ in each analysis category. When varying $\kappa_{2V}$, the other coupling modifiers affecting the VBF $HH$ production mode are set to their SM predictions. The bands indicate the simulation statistical uncertainty.

Expected yields of the signal $HH$ process as a function of the coupling modifier $\kappa_{\lambda}$ after applying the analysis preselection and in each analysis category after the final selection. The $HH$ yield is obtained considering both the $ggF$ and $VBF$ contributions. The other coupling modifiers affecting $HH$ production are fixed to their SM predictions. The bands indicate the simulation statistical uncertainty.

The acceptance times efficiency for the signal ggF $HH$ process in each analysis category as a function of the $c_{hhh}$ HEFT coefficients. The dashed lines denote values that are excluded at 95% CL. The bottom panels show the efficiency of the sum of all analysis categories.

The acceptance times efficiency for the signal ggF $HH$ process in each analysis category as a function of the $c_{tthh}$ HEFT coefficients. The dashed lines denote values that are excluded at 95% CL. The bottom panels show the efficiency of the sum of all analysis categories.

The acceptance times efficiency for the signal ggF $HH$ process in each analysis category as a function of the $c_{gghh}$ HEFT coefficients. The dashed lines denote values that are excluded at 95% CL. The bottom panels show the efficiency of the sum of all analysis categories.

The acceptance times efficiency for the signal ggF $HH$ process in each analysis category as a function of the $c_{H}$ SMEFT coefficients. The dashed lines denote values that are excluded at 95% CL. The bottom panels show the efficiency of the sum of all analysis categories.

The acceptance times efficiency for the signal ggF $HH$ process in each analysis category as a function of the $c_{{H}\boxed{}}$ SMEFT coefficients. The dashed lines denote values that are excluded at 95% CL. The bottom panels show the efficiency of the sum of all analysis categories.

The yield of the signal ggF $HH$ process in each analysis category as a function of the $c_{hhh}$ HEFT coefficients. The dashed lines denote values that are excluded at 95% CL. The bottom panels show the efficiency of the sum of all analysis categories.

The yield of the signal ggF $HH$ process in each analysis category as a function of the $c_{tthh}$ HEFT coefficients. The dashed lines denote values that are excluded at 95% CL. The bottom panels show the efficiency of the sum of all analysis categories.

The yield of the signal ggF $HH$ process in each analysis category as a function of the $c_{gghh}$ HEFT coefficients. The dashed lines denote values that are excluded at 95% CL. The bottom panels show the efficiency of the sum of all analysis categories.

The yield of the signal ggF $HH$ process in each analysis category as a function of the $c_{H}$ SMEFT coefficients. The dashed lines denote values that are excluded at 95% CL. The bottom panels show the efficiency of the sum of all analysis categories.

The yield of the signal ggF $HH$ process in each analysis category as a function of the $c_{{H}\boxed{}}$ SMEFT coefficients. The dashed lines denote values that are excluded at 95% CL. The bottom panels show the efficiency of the sum of all analysis categories.

Distributions of the diphoton invariant mass for events in data (black dots with error bars) compared to the sum of the expected signal and backgrounds (histograms) in the low mass 1 region.

Distributions of the diphoton invariant mass for events in data (black dots with error bars) compared to the sum of the expected signal and backgrounds (histograms) in the low mass 2 region.

Distributions of the diphoton invariant mass for events in data (black dots with error bars) compared to the sum of the expected signal and backgrounds (histograms) in the low mass 3 region.

Distributions of the diphoton invariant mass for events in data (black dots with error bars) compared to the sum of the expected signal and backgrounds (histograms) in the low mass 4 region.

Distributions of the diphoton invariant mass for events in data (black dots with error bars) compared to the sum of the expected signal and backgrounds (histograms) in the high mass 1 region.

Distributions of the diphoton invariant mass for events in data (black dots with error bars) compared to the sum of the expected signal and backgrounds (histograms) in the high mass 2 region.

Distributions of the diphoton invariant mass for events in data (black dots with error bars) compared to the sum of the expected signal and backgrounds (histograms) in the high mass 3 region.

Observed (solid line) value of $-2\ln\Lambda$ as a function of the $c_{hhh}$ HEFT coefficients when all other coupling modifiers are fixed to the SM value

Expected (dashed line) value of $-2\ln\Lambda$ as a function of the $c_{hhh}$ HEFT coefficients when all other coupling modifiers are fixed to the SM value

Observed (solid line) value of $-2\ln\Lambda$ as a function of the $c_{gghh}$ HEFT coefficients when all other coupling modifiers are fixed to the SM value

Expected (dashed line) value of $-2\ln\Lambda$ as a function of the $c_{gghh}$ HEFT coefficients when all other coupling modifiers are fixed to the SM value

Observed (solid line) value of $-2\ln\Lambda$ as a function of the $c_{tthh}$ HEFT coefficients when all other coupling modifiers are fixed to the SM value

Expected (dashed line) value of $-2\ln\Lambda$ as a function of the $c_{tthh}$ HEFT coefficients when all other coupling modifiers are fixed to the SM value

Observed (solid line) value of $-2\ln\Lambda$ as a function of the $c_{H}$ SMEFT coefficients when all other coupling modifiers are fixed to the SM value.

Expected (dashed line) value of $-2\ln\Lambda$ as a function of the $c_{H}$ SMEFT coefficients when all other coupling modifiers are fixed to the SM value.

Observed (solid line) value of $-2\ln\Lambda$ as a function of the $c_{{H}\boxed{}}$ SMEFT coefficients when all other coupling modifiers are fixed to the SM value.

Expected (dashed line) value of $-2\ln\Lambda$ as a function of the $c_{{H}\boxed{}}$ SMEFT coefficients when all other coupling modifiers are fixed to the SM value.

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.

Likelihood profiles of the $\kappa_\lambda$ parameter.

Likelihood profiles of the $\kappa_{2V}$ parameter.

The 95% CL upper limits on the $HH$ cross-section for different values of $\kappa_\lambda$. The expected limits assume no $HH$ production. The red line shows the theory prediction for the combined ggF $HH$ and VBF $HH$ cross-section as a function of $\kappa_\lambda$ where all parameters and couplings are set to their SM values except $\kappa_\lambda$. The bands surrounding the red cross-section line indicate the theoretical uncertainty on the predicted cross-section. The $\kappa_\lambda$ values in the interval $[-6.7, 15.3]$ cannot be excluded by this analysis.

The 95% CL upper limits on the VBF $HH$ cross section for different values of $\kappa_{2V}$. The expected limits assume no VBF $HH$ production. SM ggF $HH$ production is considered as background. The red line shows the theory prediction for the VBF $HH$ cross-section as a function of $\kappa_{2V}$ where all parameters and couplings are set to their SM values except $\kappa_{2V}$. The uncertainty band is smaller than the width of the plotted line. The $\kappa_{2V}$ values in the interval $[-0.3, 2.4]$ cannot be excluded by this analysis.

Cross sections (fb) at large Z $p_{T}$ values in the Z -> $\ell\ell$ and Z -> $\nu\nu$ channels, and their combination. The theoretical predictions from Madgraph at NLO in QCD and corrected to NLO in electroweak using the NNPDF 3.0 are also reported. With the exception of the largest Z $p_{T}$ bin, the statistical uncertainties in the measurements are much smaller than the systematic uncertainties. Both measurements and predictions correspond to $\sigma \mathcal{B}(Z -> \ell\ell)$, where $\sigma$ is the total fiducial cross section, $\mathcal{B}$ is the branching fraction, and $\ell$ is a charged lepton. The $Z -> \nu\nu$ measurement corresponds to $\sigma \mathcal{B}(Z -> \ell\ell)/\mathcal{B}(Z -> \nu\nu)$.

Cross sections normalized to the total cross section measurements at high Z $p_{T}$ values in the Z -> $\ell\ell$ and Z -> $\nu\nu$ channels, and in their combination. The uncertainty includes both statistical and systematic uncertainties.

The relative statistical and systematic uncertainties from various sources for the absolute cross section measurements in bins of Z $p_{T}$ on charged leptons.

The relative statistical and systematic uncertainties from various sources for the relative cross section measurements in bins of Z $p_{T}$ on charged leptons.

The relative statistical and systematic uncertainties from various sources for the absolute cross section measurements in bins of Z $p_{T}$ on neutrinos.

The relative statistical and systematic uncertainties from various sources for the relative cross section measurements in bins of Z $p_{T}$ on neutrinos.

The relative statistical and systematic uncertainties from various sources for the absolute cross section measurements in bins of Z $p_{T}$ on neutrinos and charged leptons.

The relative statistical and systematic uncertainties from various sources for the relative cross section measurements in bins of Z $p_{T}$ on neutrinos and charged leptons.

Cross section measurements in bins of Z $p_{T}$ on neutrinos and charged leptons.

Cross section measurements in bins of Z $p_{T}$ on neutrinos and charged leptons.

Higgs boson transverse momentum spectrum, $H\rightarrow b\overline{b}$

Combined Higgs boson transverse momentum spectrum, ggH only

Combined Higgs boson jet multiplicity spectrum

Higgs boson jet multiplicity spectrum, $H\rightarrow\gamma\gamma$

Higgs boson jet multiplicity spectrum, $H\rightarrow ZZ$

Combined Higgs boson rapidity spectrum

Higgs boson rapidity spectrum, $H\rightarrow\gamma\gamma$

Higgs boson rapidity spectrum, $H\rightarrow ZZ$

Combined leading jet transverse momentum spectrum

Leading jet transverse momentum spectrum, $H\rightarrow\gamma\gamma$

Leading jet transverse momentum spectrum, $H\rightarrow ZZ$