Relative systematic uncertainties in the total cross section measurement.

Measured fraction of events for $N_j = 0, 1, \geq 2$ jets. The total uncertainty is reported.

Measured fraction of events for $N_j = 0, 1, \geq 2$ jets. The total uncertainty is reported.

Correlation matrix from inclusive and normalized cross section measurements

Correlation matrix from inclusive and normalized cross section measurements

Fiducial cross sections for $N_j = 0, 1, \geq 2$ jets. The total uncertainty is reported.

Fiducial cross sections for $N_j = 0, 1, \geq 2$ jets. The total uncertainty is reported.

Covariance Matrix from fiducial and normalized cross sections measurements (total, and one and two jet bin fractions).

Covariance Matrix from fiducial and normalized cross sections measurements (total, and one and two jet bin fractions).

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.

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.

A summary of the uncertainties in the background estimates for SR-Gtt-0L-C. The individual experimental and theoretical uncertainties are assumed to be uncorrelated and are combined by adding in quadrature.

A summary of the uncertainties in the background estimates for SR-Gtt-1L-B. The individual experimental and theoretical uncertainties are assumed to be uncorrelated and are combined by adding in quadrature.

A summary of the uncertainties in the background estimates for SR-Gtt-1L-M1. The individual experimental and theoretical uncertainties are assumed to be uncorrelated and are combined by adding in quadrature.

A summary of the uncertainties in the background estimates for SR-Gtt-1L-M2. The individual experimental and theoretical uncertainties are assumed to be uncorrelated and are combined by adding in quadrature.

A summary of the uncertainties in the background estimates for SR-Gtt-1L-C. The individual experimental and theoretical uncertainties are assumed to be uncorrelated and are combined by adding in quadrature.

A summary of the uncertainties in the background estimates for SR-Gbb-B. The individual experimental and theoretical uncertainties are assumed to be uncorrelated and are combined by adding in quadrature.

A summary of the uncertainties in the background estimates for SR-Gbb-M. The individual experimental and theoretical uncertainties are assumed to be uncorrelated and are combined by adding in quadrature.

A summary of the uncertainties in the background estimates for SR-Gbb-C. The individual experimental and theoretical uncertainties are assumed to be uncorrelated and are combined by adding in quadrature.

A summary of the uncertainties in the background estimates for SR-Gtb-B. The individual experimental and theoretical uncertainties are assumed to be uncorrelated and are combined by adding in quadrature.

A summary of the uncertainties in the background estimates for SR-Gtb-M. The individual experimental and theoretical uncertainties are assumed to be uncorrelated and are combined by adding in quadrature.

A summary of the uncertainties in the background estimates for SR-Gtb-C. The individual experimental and theoretical uncertainties are assumed to be uncorrelated and are combined by adding in quadrature.

A summary of the uncertainties in the background estimates for SR-Gtt-2100-1. The individual experimental and theoretical uncertainties are assumed to be uncorrelated and are combined by adding in quadrature.

A summary of the uncertainties in the background estimates for SR-Gtt-1800-1. The individual experimental and theoretical uncertainties are assumed to be uncorrelated and are combined by adding in quadrature.

A summary of the uncertainties in the background estimates for SR-Gtt-2300-1200. The individual experimental and theoretical uncertainties are assumed to be uncorrelated and are combined by adding in quadrature.

A summary of the uncertainties in the background estimates for SR-Gtt-1900-1400. The individual experimental and theoretical uncertainties are assumed to be uncorrelated and are combined by adding in quadrature.

A summary of the uncertainties in the background estimates for SR-Gbb-2800-1400. The individual experimental and theoretical uncertainties are assumed to be uncorrelated and are combined by adding in quadrature.

A summary of the uncertainties in the background estimates for SR-Gbb-2300-1000. The individual experimental and theoretical uncertainties are assumed to be uncorrelated and are combined by adding in quadrature.

A summary of the uncertainties in the background estimates for SR-Gbb-2100-1600. The individual experimental and theoretical uncertainties are assumed to be uncorrelated and are combined by adding in quadrature.

A summary of the uncertainties in the background estimates for SR-Gbb-2000-1800. The individual experimental and theoretical uncertainties are assumed to be uncorrelated and are combined by adding in quadrature.

Results of the background-only fit extrapolated to SR_Gtt_0L_B in the CC analysis, for both the total expected background yields and the main contributing background processes. The quoted uncertainties include both experimental and theoretical systematics. The data in the SRs are not included in the fit. The background category $t\bar{t}+X$ includes $t\bar{t} W/Z$, $t\bar{t} H$ and $t\bar{t} t\bar{t}$ events. The row ``Pre-fit background'' provides the total background prediction when the $t\bar{t}$ and $Z+$jets normalisations are obtained from theoretical calculation, taking into account the kinematic weights described in Section 5.

Results of the background-only fit extrapolated to SR_Gtt_0L_M1 in the CC analysis, for both the total expected background yields and the main contributing background processes. The quoted uncertainties include both experimental and theoretical systematics. The data in the SRs are not included in the fit. The background category $t\bar{t}+X$ includes $t\bar{t} W/Z$, $t\bar{t} H$ and $t\bar{t} t\bar{t}$ events. The row ``Pre-fit background'' provides the total background prediction when the $t\bar{t}$ and $Z+$jets normalisations are obtained from theoretical calculation, taking into account the kinematic weights described in Section 5.

Results of the background-only fit extrapolated to SR_Gtt_0L_M2 in the CC analysis, for both the total expected background yields and the main contributing background processes. The quoted uncertainties include both experimental and theoretical systematics. The data in the SRs are not included in the fit. The background category $t\bar{t}+X$ includes $t\bar{t} W/Z$, $t\bar{t} H$ and $t\bar{t} t\bar{t}$ events. The row ``Pre-fit background'' provides the total background prediction when the $t\bar{t}$ and $Z+$jets normalisations are obtained from theoretical calculation, taking into account the kinematic weights described in Section 5.

Results of the background-only fit extrapolated to SR_Gtt_0L_C in the CC analysis, for both the total expected background yields and the main contributing background processes. The quoted uncertainties include both experimental and theoretical systematics. The data in the SRs are not included in the fit. The background category $t\bar{t}+X$ includes $t\bar{t} W/Z$, $t\bar{t} H$ and $t\bar{t} t\bar{t}$ events. The row ``Pre-fit background'' provides the total background prediction when the $t\bar{t}$ and $Z+$jets normalisations are obtained from theoretical calculation, taking into account the kinematic weights described in Section 5.

Results of the background-only fit extrapolated to SR_Gtt_1L_B in the CC analysis, for both the total expected background yields and the main contributing background processes. The quoted uncertainties include both experimental and theoretical systematics. The data in the SRs are not included in the fit. The background category $t\bar{t}+X$ includes $t\bar{t} W/Z$, $t\bar{t} H$ and $t\bar{t} t\bar{t}$ events. The row ``Pre-fit background'' provides the total background prediction when the $t\bar{t}$ and $Z+$jets normalisations are obtained from theoretical calculation, taking into account the kinematic weights described in Section 5.

Results of the background-only fit extrapolated to SR_Gtt_1L_M1 in the CC analysis, for both the total expected background yields and the main contributing background processes. The quoted uncertainties include both experimental and theoretical systematics. The data in the SRs are not included in the fit. The background category $t\bar{t}+X$ includes $t\bar{t} W/Z$, $t\bar{t} H$ and $t\bar{t} t\bar{t}$ events. The row ``Pre-fit background'' provides the total background prediction when the $t\bar{t}$ and $Z+$jets normalisations are obtained from theoretical calculation, taking into account the kinematic weights described in Section 5.

Results of the background-only fit extrapolated to SR_Gtt_1L_M2 in the CC analysis, for both the total expected background yields and the main contributing background processes. The quoted uncertainties include both experimental and theoretical systematics. The data in the SRs are not included in the fit. The background category $t\bar{t}+X$ includes $t\bar{t} W/Z$, $t\bar{t} H$ and $t\bar{t} t\bar{t}$ events. The row ``Pre-fit background'' provides the total background prediction when the $t\bar{t}$ and $Z+$jets normalisations are obtained from theoretical calculation, taking into account the kinematic weights described in Section 5.

Results of the background-only fit extrapolated to SR_Gtt_1L_C in the CC analysis, for both the total expected background yields and the main contributing background processes. The quoted uncertainties include both experimental and theoretical systematics. The data in the SRs are not included in the fit. The background category $t\bar{t}+X$ includes $t\bar{t} W/Z$, $t\bar{t} H$ and $t\bar{t} t\bar{t}$ events. The row ``Pre-fit background'' provides the total background prediction when the $t\bar{t}$ and $Z+$jets normalisations are obtained from theoretical calculation, taking into account the kinematic weights described in Section 5.

Results of the background-only fit extrapolated to SR_Gbb_B in the CC analysis, for both the total expected background yields and the main contributing background processes. The quoted uncertainties include both experimental and theoretical systematics. The data in the SRs are not included in the fit. The background category $t\bar{t}+X$ includes $t\bar{t} W/Z$, $t\bar{t} H$ and $t\bar{t} t\bar{t}$ events. The row ``Pre-fit background'' provides the total background prediction when the $t\bar{t}$ and $Z+$jets normalisations are obtained from theoretical calculation, taking into account the kinematic weights described in Section 5.

Results of the background-only fit extrapolated to SR_Gbb_M in the CC analysis, for both the total expected background yields and the main contributing background processes. The quoted uncertainties include both experimental and theoretical systematics. The data in the SRs are not included in the fit. The background category $t\bar{t}+X$ includes $t\bar{t} W/Z$, $t\bar{t} H$ and $t\bar{t} t\bar{t}$ events. The row ``Pre-fit background'' provides the total background prediction when the $t\bar{t}$ and $Z+$jets normalisations are obtained from theoretical calculation, taking into account the kinematic weights described in Section 5.

Results of the background-only fit extrapolated to SR_Gbb_C in the CC analysis, for both the total expected background yields and the main contributing background processes. The quoted uncertainties include both experimental and theoretical systematics. The data in the SRs are not included in the fit. The background category $t\bar{t}+X$ includes $t\bar{t} W/Z$, $t\bar{t} H$ and $t\bar{t} t\bar{t}$ events. The row ``Pre-fit background'' provides the total background prediction when the $t\bar{t}$ and $Z+$jets normalisations are obtained from theoretical calculation, taking into account the kinematic weights described in Section 5.

Results of the background-only fit extrapolated to SR_Gtb_B in the CC analysis, for both the total expected background yields and the main contributing background processes. The quoted uncertainties include both experimental and theoretical systematics. The data in the SRs are not included in the fit. The background category $t\bar{t}+X$ includes $t\bar{t} W/Z$, $t\bar{t} H$ and $t\bar{t} t\bar{t}$ events. The row ``Pre-fit background'' provides the total background prediction when the $t\bar{t}$ and $Z+$jets normalisations are obtained from theoretical calculation, taking into account the kinematic weights described in Section 5.

Results of the background-only fit extrapolated to SR_Gtb_M in the CC analysis, for both the total expected background yields and the main contributing background processes. The quoted uncertainties include both experimental and theoretical systematics. The data in the SRs are not included in the fit. The background category $t\bar{t}+X$ includes $t\bar{t} W/Z$, $t\bar{t} H$ and $t\bar{t} t\bar{t}$ events. The row ``Pre-fit background'' provides the total background prediction when the $t\bar{t}$ and $Z+$jets normalisations are obtained from theoretical calculation, taking into account the kinematic weights described in Section 5.

Results of the background-only fit extrapolated to SR_Gtb_C in the CC analysis, for both the total expected background yields and the main contributing background processes. The quoted uncertainties include both experimental and theoretical systematics. The data in the SRs are not included in the fit. The background category $t\bar{t}+X$ includes $t\bar{t} W/Z$, $t\bar{t} H$ and $t\bar{t} t\bar{t}$ events. The row ``Pre-fit background'' provides the total background prediction when the $t\bar{t}$ and $Z+$jets normalisations are obtained from theoretical calculation, taking into account the kinematic weights described in Section 5.

Results of the background-only fit extrapolated to SR_Gtt_2100_1 in the NN analysis, for both the total expected background yields and the main contributing background processes. The quoted uncertainties include both experimental and theoretical systematics. The data in the SRs are not included in the fit. The background category $t\bar{t}+X$ includes $t\bar{t} W/Z$, $t\bar{t} H$ and $t\bar{t} t\bar{t}$ events. The row ``Pre-fit background'' provides the total background prediction when the $t\bar{t}$ and $Z+$jets normalisations are obtained from theoretical calculation, taking into account the kinematic weights described in Section 5.

Results of the background-only fit extrapolated to SR_Gtt_1800_1 in the NN analysis, for both the total expected background yields and the main contributing background processes. The quoted uncertainties include both experimental and theoretical systematics. The data in the SRs are not included in the fit. The background category $t\bar{t}+X$ includes $t\bar{t} W/Z$, $t\bar{t} H$ and $t\bar{t} t\bar{t}$ events. The row ``Pre-fit background'' provides the total background prediction when the $t\bar{t}$ and $Z+$jets normalisations are obtained from theoretical calculation, taking into account the kinematic weights described in Section 5.

Results of the background-only fit extrapolated to SR_Gtt_2300_1200 in the NN analysis, for both the total expected background yields and the main contributing background processes. The quoted uncertainties include both experimental and theoretical systematics. The data in the SRs are not included in the fit. The background category $t\bar{t}+X$ includes $t\bar{t} W/Z$, $t\bar{t} H$ and $t\bar{t} t\bar{t}$ events. The row ``Pre-fit background'' provides the total background prediction when the $t\bar{t}$ and $Z+$jets normalisations are obtained from theoretical calculation, taking into account the kinematic weights described in Section 5.

Results of the background-only fit extrapolated to SR_Gtt_1900_1400 in the NN analysis, for both the total expected background yields and the main contributing background processes. The quoted uncertainties include both experimental and theoretical systematics. The data in the SRs are not included in the fit. The background category $t\bar{t}+X$ includes $t\bar{t} W/Z$, $t\bar{t} H$ and $t\bar{t} t\bar{t}$ events. The row ``Pre-fit background'' provides the total background prediction when the $t\bar{t}$ and $Z+$jets normalisations are obtained from theoretical calculation, taking into account the kinematic weights described in Section 5.

Results of the background-only fit extrapolated to SR_Gbb_2800_1400 in the NN analysis, for both the total expected background yields and the main contributing background processes. The quoted uncertainties include both experimental and theoretical systematics. The data in the SRs are not included in the fit. The background category $t\bar{t}+X$ includes $t\bar{t} W/Z$, $t\bar{t} H$ and $t\bar{t} t\bar{t}$ events. The row ``Pre-fit background'' provides the total background prediction when the $t\bar{t}$ and $Z+$jets normalisations are obtained from theoretical calculation, taking into account the kinematic weights described in Section 5.

Results of the background-only fit extrapolated to SR_Gbb_2300_1000 in the NN analysis, for both the total expected background yields and the main contributing background processes. The quoted uncertainties include both experimental and theoretical systematics. The data in the SRs are not included in the fit. The background category $t\bar{t}+X$ includes $t\bar{t} W/Z$, $t\bar{t} H$ and $t\bar{t} t\bar{t}$ events. The row ``Pre-fit background'' provides the total background prediction when the $t\bar{t}$ and $Z+$jets normalisations are obtained from theoretical calculation, taking into account the kinematic weights described in Section 5.

Results of the background-only fit extrapolated to SR_Gbb_2100_1600 in the NN analysis, for both the total expected background yields and the main contributing background processes. The quoted uncertainties include both experimental and theoretical systematics. The data in the SRs are not included in the fit. The background category $t\bar{t}+X$ includes $t\bar{t} W/Z$, $t\bar{t} H$ and $t\bar{t} t\bar{t}$ events. The row ``Pre-fit background'' provides the total background prediction when the $t\bar{t}$ and $Z+$jets normalisations are obtained from theoretical calculation, taking into account the kinematic weights described in Section 5.

Results of the background-only fit extrapolated to SR_Gbb_2000_1800 in the NN analysis, for both the total expected background yields and the main contributing background processes. The quoted uncertainties include both experimental and theoretical systematics. The data in the SRs are not included in the fit. The background category $t\bar{t}+X$ includes $t\bar{t} W/Z$, $t\bar{t} H$ and $t\bar{t} t\bar{t}$ events. The row ``Pre-fit background'' provides the total background prediction when the $t\bar{t}$ and $Z+$jets normalisations are obtained from theoretical calculation, taking into account the kinematic weights described in Section 5.

Observed exclusion limit in the $\tilde{g}$--$\tilde\chi^0_1$ mass plane for the Gtt models obtained from the NN analysis. The dashed and solid bold lines show the 95\% CL expected and observed limits, respectively. The shaded bands around the expected limits show the impact of the experimental and background uncertainties. The dotted lines show the impact on the observed limit of the variation of the nominal signal cross-section by $\pm 1 \sigma$ of its theoretical uncertainty.

Observed exclusion limit in the $\tilde{g}$--$\tilde\chi^0_1$ mass plane for the Gtt models obtained from the NN analysis. The dashed and solid bold lines show the 95\% CL expected and observed limits, respectively. The shaded bands around the expected limits show the impact of the experimental and background uncertainties. The dotted lines show the impact on the observed limit of the variation of the nominal signal cross-section by $\pm 1 \sigma$ of its theoretical uncertainty.

Observed exclusion limit in the $\tilde{g}$--$\tilde\chi^0_1$ mass plane for the Gtt models obtained from the NN analysis. The dashed and solid bold lines show the 95\% CL expected and observed limits, respectively. The shaded bands around the expected limits show the impact of the experimental and background uncertainties. The dotted lines show the impact on the observed limit of the variation of the nominal signal cross-section by $\pm 1 \sigma$ of its theoretical uncertainty.

Expected exclusion limit in the $\tilde{g}$--$\tilde\chi^0_1$ mass plane for the Gtt models obtained from the NN analysis. The dashed and solid bold lines show the 95\% CL expected and observed limits, respectively. The shaded bands around the expected limits show the impact of the experimental and background uncertainties. The dotted lines show the impact on the observed limit of the variation of the nominal signal cross-section by $\pm 1 \sigma$ of its theoretical uncertainty.

Expected exclusion limit in the $\tilde{g}$--$\tilde\chi^0_1$ mass plane for the Gtt models obtained from the NN analysis. The dashed and solid bold lines show the 95\% CL expected and observed limits, respectively. The shaded bands around the expected limits show the impact of the experimental and background uncertainties. The dotted lines show the impact on the observed limit of the variation of the nominal signal cross-section by $\pm 1 \sigma$ of its theoretical uncertainty.

Expected exclusion limit in the $\tilde{g}$--$\tilde\chi^0_1$ mass plane for the Gtt models obtained from the NN analysis. The dashed and solid bold lines show the 95\% CL expected and observed limits, respectively. The shaded bands around the expected limits show the impact of the experimental and background uncertainties. The dotted lines show the impact on the observed limit of the variation of the nominal signal cross-section by $\pm 1 \sigma$ of its theoretical uncertainty.

Observed exclusion limit in the $\tilde{g}$--$\tilde\chi^0_1$ mass plane for the Gbb models obtained from the NN analysis. The dashed and solid bold lines show the 95\% CL expected and observed limits, respectively. The shaded bands around the expected limits show the impact of the experimental and background uncertainties. The dotted lines show the impact on the observed limit of the variation of the nominal signal cross-section by $\pm 1 \sigma$ of its theoretical uncertainty.

Observed exclusion limit in the $\tilde{g}$--$\tilde\chi^0_1$ mass plane for the Gbb models obtained from the NN analysis. The dashed and solid bold lines show the 95\% CL expected and observed limits, respectively. The shaded bands around the expected limits show the impact of the experimental and background uncertainties. The dotted lines show the impact on the observed limit of the variation of the nominal signal cross-section by $\pm 1 \sigma$ of its theoretical uncertainty.

Observed exclusion limit in the $\tilde{g}$--$\tilde\chi^0_1$ mass plane for the Gbb models obtained from the NN analysis. The dashed and solid bold lines show the 95\% CL expected and observed limits, respectively. The shaded bands around the expected limits show the impact of the experimental and background uncertainties. The dotted lines show the impact on the observed limit of the variation of the nominal signal cross-section by $\pm 1 \sigma$ of its theoretical uncertainty.

Expected exclusion limit in the $\tilde{g}$--$\tilde\chi^0_1$ mass plane for the Gbb models obtained from the NN analysis. The dashed and solid bold lines show the 95\% CL expected and observed limits, respectively. The shaded bands around the expected limits show the impact of the experimental and background uncertainties. The dotted lines show the impact on the observed limit of the variation of the nominal signal cross-section by $\pm 1 \sigma$ of its theoretical uncertainty.

Expected exclusion limit in the $\tilde{g}$--$\tilde\chi^0_1$ mass plane for the Gbb models obtained from the NN analysis. The dashed and solid bold lines show the 95\% CL expected and observed limits, respectively. The shaded bands around the expected limits show the impact of the experimental and background uncertainties. The dotted lines show the impact on the observed limit of the variation of the nominal signal cross-section by $\pm 1 \sigma$ of its theoretical uncertainty.

Expected exclusion limit in the $\tilde{g}$--$\tilde\chi^0_1$ mass plane for the Gbb models obtained from the NN analysis. The dashed and solid bold lines show the 95\% CL expected and observed limits, respectively. The shaded bands around the expected limits show the impact of the experimental and background uncertainties. The dotted lines show the impact on the observed limit of the variation of the nominal signal cross-section by $\pm 1 \sigma$ of its theoretical uncertainty.

Observed (left) 95\%~CL exclusion limits on the gluino mass as a function of BR$(\tilde{g} \to b\bar{b}\tilde\chi^{0}_{1}$) (vertical) and BR$(\tilde{g} \to t\bar{t}\tilde\chi^{0}_{1}$) (horizontal) for Gtb models with $m(\tilde\chi^{0}_{1}) = 1$~GeV, obtained from the CC analysis.

Expected (right) 95\%~CL exclusion limits on the gluino mass as a function of BR$(\tilde{g} \to b\bar{b}\tilde\chi^{0}_{1}$) (vertical) and BR$(\tilde{g} \to t\bar{t}\tilde\chi^{0}_{1}$) (horizontal) for Gtb models with $m(\tilde\chi^{0}_{1}) = 1$~GeV, obtained from the CC analysis.

Observed (left) 95\%~CL exclusion limits on the gluino mass as a function of BR$(\tilde{g} \to b\bar{b}\tilde\chi^{0}_{1}$) (vertical) and BR$(\tilde{g} \to t\bar{t}\tilde\chi^{0}_{1}$) (horizontal) for Gtb models with $m(\tilde\chi^{0}_{1}) = 600$~GeV, obtained from the CC analysis.

Expected (right) 95\%~CL exclusion limits on the gluino mass as a function of BR$(\tilde{g} \to b\bar{b}\tilde\chi^{0}_{1}$) (vertical) and BR$(\tilde{g} \to t\bar{t}\tilde\chi^{0}_{1}$) (horizontal) for Gtb models with $m(\tilde\chi^{0}_{1}) = 600$~GeV, obtained from the CC analysis.

Observed (left) 95\%~CL exclusion limits on the gluino mass as a function of BR$(\tilde{g} \to b\bar{b}\tilde\chi^{0}_{1}$) (vertical) and BR$(\tilde{g} \to t\bar{t}\tilde\chi^{0}_{1}$) (horizontal) for Gtb models with $m(\tilde\chi^{0}_{1}) = 1$~TeV, obtained from the CC analysis.

Expected (right) 95\%~CL exclusion limits on the gluino mass as a function of BR$(\tilde{g} \to b\bar{b}\tilde\chi^{0}_{1}$) (vertical) and BR$(\tilde{g} \to t\bar{t}\tilde\chi^{0}_{1}$) (horizontal) for Gtb models with $m(\tilde\chi^{0}_{1}) = 1$~TeV, obtained from the CC analysis.

Observed exclusion limit in the $\tilde{g}$--$\tilde\chi^0_1$ mass plane for the Gtt models obtained from the CC analysis. The dashed and solid bold lines show the 95\% CL expected and observed limits, respectively. The shaded bands around the expected limits show the impact of the experimental and background uncertainties. The dotted lines show the impact on the observed limit of the variation of the nominal signal cross-section by $\pm 1 \sigma$ of its theoretical uncertainty.

Observed exclusion limit in the $\tilde{g}$--$\tilde\chi^0_1$ mass plane for the Gtt models obtained from the CC analysis. The dashed and solid bold lines show the 95\% CL expected and observed limits, respectively. The shaded bands around the expected limits show the impact of the experimental and background uncertainties. The dotted lines show the impact on the observed limit of the variation of the nominal signal cross-section by $\pm 1 \sigma$ of its theoretical uncertainty.

Observed exclusion limit in the $\tilde{g}$--$\tilde\chi^0_1$ mass plane for the Gtt models obtained from the CC analysis. The dashed and solid bold lines show the 95\% CL expected and observed limits, respectively. The shaded bands around the expected limits show the impact of the experimental and background uncertainties. The dotted lines show the impact on the observed limit of the variation of the nominal signal cross-section by $\pm 1 \sigma$ of its theoretical uncertainty.

Expected exclusion limit in the $\tilde{g}$--$\tilde\chi^0_1$ mass plane for the Gtt models obtained from the CC analysis. The dashed and solid bold lines show the 95\% CL expected and observed limits, respectively. The shaded bands around the expected limits show the impact of the experimental and background uncertainties. The dotted lines show the impact on the observed limit of the variation of the nominal signal cross-section by $\pm 1 \sigma$ of its theoretical uncertainty.

Expected exclusion limit in the $\tilde{g}$--$\tilde\chi^0_1$ mass plane for the Gtt models obtained from the CC analysis. The dashed and solid bold lines show the 95\% CL expected and observed limits, respectively. The shaded bands around the expected limits show the impact of the experimental and background uncertainties. The dotted lines show the impact on the observed limit of the variation of the nominal signal cross-section by $\pm 1 \sigma$ of its theoretical uncertainty.

Expected exclusion limit in the $\tilde{g}$--$\tilde\chi^0_1$ mass plane for the Gtt models obtained from the CC analysis. The dashed and solid bold lines show the 95\% CL expected and observed limits, respectively. The shaded bands around the expected limits show the impact of the experimental and background uncertainties. The dotted lines show the impact on the observed limit of the variation of the nominal signal cross-section by $\pm 1 \sigma$ of its theoretical uncertainty.

Observed exclusion limit in the $\tilde{g}$--$\tilde\chi^0_1$ mass plane for the Gbb models obtained from the CC analysis. The dashed and solid bold lines show the 95\% CL expected and observed limits, respectively. The shaded bands around the expected limits show the impact of the experimental and background uncertainties. The dotted lines show the impact on the observed limit of the variation of the nominal signal cross-section by $\pm 1 \sigma$ of its theoretical uncertainty.

Observed exclusion limit in the $\tilde{g}$--$\tilde\chi^0_1$ mass plane for the Gbb models obtained from the CC analysis. The dashed and solid bold lines show the 95\% CL expected and observed limits, respectively. The shaded bands around the expected limits show the impact of the experimental and background uncertainties. The dotted lines show the impact on the observed limit of the variation of the nominal signal cross-section by $\pm 1 \sigma$ of its theoretical uncertainty.

Observed exclusion limit in the $\tilde{g}$--$\tilde\chi^0_1$ mass plane for the Gbb models obtained from the CC analysis. The dashed and solid bold lines show the 95\% CL expected and observed limits, respectively. The shaded bands around the expected limits show the impact of the experimental and background uncertainties. The dotted lines show the impact on the observed limit of the variation of the nominal signal cross-section by $\pm 1 \sigma$ of its theoretical uncertainty.

Expected exclusion limit in the $\tilde{g}$--$\tilde\chi^0_1$ mass plane for the Gbb models obtained from the CC analysis. The dashed and solid bold lines show the 95\% CL expected and observed limits, respectively. The shaded bands around the expected limits show the impact of the experimental and background uncertainties. The dotted lines show the impact on the observed limit of the variation of the nominal signal cross-section by $\pm 1 \sigma$ of its theoretical uncertainty.

Expected exclusion limit in the $\tilde{g}$--$\tilde\chi^0_1$ mass plane for the Gbb models obtained from the CC analysis. The dashed and solid bold lines show the 95\% CL expected and observed limits, respectively. The shaded bands around the expected limits show the impact of the experimental and background uncertainties. The dotted lines show the impact on the observed limit of the variation of the nominal signal cross-section by $\pm 1 \sigma$ of its theoretical uncertainty.

Expected exclusion limit in the $\tilde{g}$--$\tilde\chi^0_1$ mass plane for the Gbb models obtained from the CC analysis. The dashed and solid bold lines show the 95\% CL expected and observed limits, respectively. The shaded bands around the expected limits show the impact of the experimental and background uncertainties. The dotted lines show the impact on the observed limit of the variation of the nominal signal cross-section by $\pm 1 \sigma$ of its theoretical uncertainty.

Upper limit at 95\% CL on the cross-section times branching ratio (fb) in the $\tilde{g}$--$\tilde\chi^0_1$ mass plane for the Gbb (right) models obtained from the CC analysis. The numbers give the observed 95\% CL upper limit on the cross section in fb, with the label colour matching the associated best-expected region. Only a lower limit on the excluded cross section (>0.7 fb) is given at some points due to the very small number events expected and observed in the chosen SR. The dashed and solid bold lines show the 95\% CL expected and observed limits, respectively. The shaded bands around the expected limits show the impact of the experimental and background theoretical uncertainties. The dotted lines show the impact on the observed limit of the variation of the nominal signal cross-section by $\pm1\sigma$ of its theoretical uncertainty.

Upper limit at 95\% CL on the cross-section times branching ratio (fb) in the $\tilde{g}$--$\tilde\chi^0_1$ mass plane for the Gbb (right) models obtained from the NN analysis. The numbers give the observed 95\% CL upper limit on the cross section in fb, with the label colour matching the associated best-expected region. Only a lower limit on the excluded cross section (>0.7 fb) is given at some points due to the very small number events expected and observed in the chosen SR. The dashed and solid bold lines show the 95\% CL expected and observed limits, respectively. The shaded bands around the expected limits show the impact of the experimental and background theoretical uncertainties. The dotted lines show the impact on the observed limit of the variation of the nominal signal cross-section by $\pm1\sigma$ of its theoretical uncertainty.

Upper limit at 95\% CL on the cross-section times branching ratio (fb) in the $\tilde{g}$--$\tilde\chi^0_1$ mass plane for the Gtt (left) models obtained from the CC analysis. The numbers give the observed 95\% CL upper limit on the cross section in fb, with the label colour matching the associated best-expected region. Only a lower limit on the excluded cross section (>0.7 fb) is given at some points due to the very small number events expected and observed in the chosen SR. The dashed and solid bold lines show the 95\% CL expected and observed limits, respectively. The shaded bands around the expected limits show the impact of the experimental and background theoretical uncertainties. The dotted lines show the impact on the observed limit of the variation of the nominal signal cross-section by $\pm1\sigma$ of its theoretical uncertainty.

Upper limit at 95\% CL on the cross-section times branching ratio (fb) in the $\tilde{g}$--$\tilde\chi^0_1$ mass plane for the Gtt (left) models obtained from the NN analysis. The numbers give the observed 95\% CL upper limit on the cross section in fb, with the label colour matching the associated best-expected region. Only a lower limit on the excluded cross section (>0.7 fb) is given at some points due to the very small number events expected and observed in the chosen SR. The dashed and solid bold lines show the 95\% CL expected and observed limits, respectively. The shaded bands around the expected limits show the impact of the experimental and background theoretical uncertainties. The dotted lines show the impact on the observed limit of the variation of the nominal signal cross-section by $\pm1\sigma$ of its theoretical uncertainty.

Acceptance for SR-Gtt-0L-B and the $\tilde{g}\rightarrow t\bar{t}\tilde\chi^0_1$ signal process.

Efficiency for SR-Gtt-0L-B and the $\tilde{g}\rightarrow t\bar{t}\tilde\chi^0_1$ signal process.

Acceptance for SR-Gtt-0L-M1 and the $\tilde{g}\rightarrow t\bar{t}\tilde\chi^0_1$ signal process.

Efficiency for SR-Gtt-0L-M1 and the $\tilde{g}\rightarrow t\bar{t}\tilde\chi^0_1$ signal process.

Acceptance for SR-Gtt-0L-M2 and the $\tilde{g}\rightarrow t\bar{t}\tilde\chi^0_1$ signal process.

Efficiency for SR-Gtt-0L-M2 and the $\tilde{g}\rightarrow t\bar{t}\tilde\chi^0_1$ signal process.

Acceptance for SR-Gtt-0L-C and the $\tilde{g}\rightarrow t\bar{t}\tilde\chi^0_1$ signal process.

Efficiency for SR-Gtt-0L-C and the $\tilde{g}\rightarrow t\bar{t}\tilde\chi^0_1$ signal process.

Acceptance for SR-Gtt-1L-B and the $\tilde{g}\rightarrow t\bar{t}\tilde\chi^0_1$ signal process.

Efficiency for SR-Gtt-1L-B and the $\tilde{g}\rightarrow t\bar{t}\tilde\chi^0_1$ signal process.

Acceptance for SR-Gtt-1L-M1 and the $\tilde{g}\rightarrow t\bar{t}\tilde\chi^0_1$ signal process.

Efficiency for SR-Gtt-1L-M1 and the $\tilde{g}\rightarrow t\bar{t}\tilde\chi^0_1$ signal process.

Acceptance for SR-Gtt-1L-M2 and the $\tilde{g}\rightarrow t\bar{t}\tilde\chi^0_1$ signal process.

Efficiency for SR-Gtt-1L-M2 and the $\tilde{g}\rightarrow t\bar{t}\tilde\chi^0_1$ signal process.

Acceptance for SR-Gtt-1L-C and the $\tilde{g}\rightarrow t\bar{t}\tilde\chi^0_1$ signal process.

Efficiency for SR-Gtt-1L-C and the $\tilde{g}\rightarrow t\bar{t}\tilde\chi^0_1$ signal process.

Acceptance for SR-Gbb-B and the $\tilde{g}\rightarrow b\bar{b}\tilde\chi^0_1$ signal process.

Efficiency for SR-Gbb-B and the $\tilde{g}\rightarrow b\bar{b}\tilde\chi^0_1$ signal process.

Acceptance for SR-Gbb-M and the $\tilde{g}\rightarrow b\bar{b}\tilde\chi^0_1$ signal process.

Efficiency for SR-Gbb-M and the $\tilde{g}\rightarrow b\bar{b}\tilde\chi^0_1$ signal process.

Acceptance for SR-Gbb-C and the $\tilde{g}\rightarrow b\bar{b}\tilde\chi^0_1$ signal process.

Efficiency for SR-Gbb-C and the $\tilde{g}\rightarrow b\bar{b}\tilde\chi^0_1$ signal process.

Acceptance for SR-Gtt-2100-1 and the $\tilde{g}\rightarrow t\bar{t}\tilde\chi^0_1$ signal process.

Efficiency for SR-Gtt-2100-1 and the $\tilde{g}\rightarrow t\bar{t}\tilde\chi^0_1$ signal process.

Acceptance for SR-Gtt-1800-1 and the $\tilde{g}\rightarrow t\bar{t}\tilde\chi^0_1$ signal process.

Efficiency for SR-Gtt-1800-1 and the $\tilde{g}\rightarrow t\bar{t}\tilde\chi^0_1$ signal process.

Acceptance for SR-Gtt-2300-1200 and the $\tilde{g}\rightarrow t\bar{t}\tilde\chi^0_1$ signal process.

Efficiency for SR-Gtt-2300-1200 and the $\tilde{g}\rightarrow t\bar{t}\tilde\chi^0_1$ signal process.

Acceptance for SR-Gtt-1900-1400 and the $\tilde{g}\rightarrow t\bar{t}\tilde\chi^0_1$ signal process.

Efficiency for SR-Gtt-1900-1400 and the $\tilde{g}\rightarrow t\bar{t}\tilde\chi^0_1$ signal process.

Acceptance for SR-Gbb-2800-1400 and the $\tilde{g}\rightarrow b\bar{b}\tilde\chi^0_1$ signal process.

Efficiency for SR-Gbb-2800-1400 and the $\tilde{g}\rightarrow b\bar{b}\tilde\chi^0_1$ signal process.

Acceptance for SR-Gbb-2300-1000 and the $\tilde{g}\rightarrow b\bar{b}\tilde\chi^0_1$ signal process.

Efficiency for SR-Gbb-2300-1000 and the $\tilde{g}\rightarrow b\bar{b}\tilde\chi^0_1$ signal process.

Acceptance for SR-Gbb-2100-1600 and the $\tilde{g}\rightarrow b\bar{b}\tilde\chi^0_1$ signal process.

Efficiency for SR-Gbb-2100-1600 and the $\tilde{g}\rightarrow b\bar{b}\tilde\chi^0_1$ signal process.

Acceptance for SR-Gbb-2000-1800 and the $\tilde{g}\rightarrow b\bar{b}\tilde\chi^0_1$ signal process.

Efficiency for SR-Gbb-2000-1800 and the $\tilde{g}\rightarrow b\bar{b}\tilde\chi^0_1$ signal process.

Cutflow for the SR-Gtt-0L-B for a representative Gtt signal. Signal was generated with 30000 events. Expected yields are normalised to a luminosity of 139~fb$^{-1}$.

Cutflow for the SR-Gtt-0L-M1 for a representative Gtt signal. Signal was generated with 30000 events. Expected yields are normalised to a luminosity of 139~fb$^{-1}$.

Cutflow for the SR-Gtt-0L-M2 for a representative Gtt signal. Signal was generated with 30000 events. Expected yields are normalised to a luminosity of 139~fb$^{-1}$.

Cutflow for the SR-Gtt-0L-C for a representative Gtt signal. Signal was generated with 30000 events. Expected yields are normalised to a luminosity of 139~fb$^{-1}$.

Cutflow for the SR-Gtt-1L-B for a representative Gtt signal. Signal was generated with 30000 events. Expected yields are normalised to a luminosity of 139~fb$^{-1}$.

Cutflow for the SR-Gtt-1L-M1 for a representative Gtt signal. Signal was generated with 30000 events. Expected yields are normalised to a luminosity of 139~fb$^{-1}$.

Cutflow for the SR-Gtt-1L-M2 for a representative Gtt signal. Signal was generated with 30000 events. Expected yields are normalised to a luminosity of 139~fb$^{-1}$.

Cutflow for the SR-Gtt-1L-C for a representative Gtt signal. Signal was generated with 30000 events. Expected yields are normalised to a luminosity of 139~fb$^{-1}$.

Cutflow for the SR-Gbb-B for a representative Gbb signal. Signal was generated with 30000 events. Expected yields are normalised to a luminosity of 139~fb$^{-1}$.

Cutflow for the SR-Gbb-M for a representative Gbb signal. Signal was generated with 30000 events. Expected yields are normalised to a luminosity of 139~fb$^{-1}$.

Cutflow for the SR-Gbb-C for a representative Gbb signal. Signal was generated with 30000 events. Expected yields are normalised to a luminosity of 139~fb$^{-1}$.

Cutflow for the SR-Gtb-B for a representative Gtb signal. Signal was generated with 30000 events. Expected yields are normalised to a luminosity of 139~fb$^{-1}$.

Cutflow for the SR-Gtb-M for a representative Gtb signal. Signal was generated with 30000 events. Expected yields are normalised to a luminosity of 139~fb$^{-1}$.

Cutflow for the SR-Gtb-C for a representative Gtb signal. Signal was generated with 30000 events. Expected yields are normalised to a luminosity of 139~fb$^{-1}$.

Cutflow for the SR-Gtt-2100-1 for a representative Gtt signal. Signal was generated with 30000 events. Expected yields are normalised to a luminosity of 139~fb$^{-1}$.

Cutflow for the SR-Gtt-1800-1 for a representative Gtt signal. Signal was generated with 30000 events. Expected yields are normalised to a luminosity of 139~fb$^{-1}$.

Cutflow for the SR-Gtt-2300-1200 for a representative Gtt signal. Signal was generated with 30000 events. Expected yields are normalised to a luminosity of 139~fb$^{-1}$.

Cutflow for the SR-Gtt-1900-1400 for a representative Gtt signal. Signal was generated with 30000 events. Expected yields are normalised to a luminosity of 139~fb$^{-1}$.

Cutflow for the SR-Gbb-2800-1400 for a representative Gbb signal. Signal was generated with 30000 events. Expected yields are normalised to a luminosity of 139~fb$^{-1}$.

Cutflow for the SR-Gbb-2300-1000 for a representative Gbb signal. Signal was generated with 30000 events. Expected yields are normalised to a luminosity of 139~fb$^{-1}$.

Cutflow for the SR-Gbb-2100-1600 for a representative Gbb signal. Signal was generated with 30000 events. Expected yields are normalised to a luminosity of 139~fb$^{-1}$.

Cutflow for the SR-Gbb-2000-1800 for a representative Gbb signal. Signal was generated with 30000 events. Expected yields are normalised to a luminosity of 139~fb$^{-1}$.

Acceptance times efficiency as a function of pseudoscalar resonant mass for each event selection step in the search for a neutral heavy narrow-width pseudoscalar resonance that decays to a 260 GeV scalar resonance and a Z boson, which decays to neutrinos.

Acceptance times efficiency as a function of pseudoscalar resonant mass for each event selection step in the search for a neutral heavy narrow-width pseudoscalar resonance that decays to a 300 GeV scalar resonance and a Z boson, which decays to neutrinos.

Acceptance times efficiency as a function of pseudoscalar resonant mass for each event selection step in the search for a neutral heavy narrow-width pseudoscalar resonance that decays to a 400 GeV scalar resonance and a Z boson, which decays to neutrinos.

Acceptance times efficiency as a function of pseudoscalar resonant mass for each event selection step in the search for a neutral heavy large-width pseudoscalar resonance that decays to a 260 GeV scalar resonance and a Z boson, which decays to neutrinos.

Acceptance times efficiency as a function of pseudoscalar resonant mass for each event selection step in the search for a neutral heavy large-width pseudoscalar resonance that decays to a 300 GeV scalar resonance and a Z boson, which decays to neutrinos.

Acceptance times efficiency as a function of pseudoscalar resonant mass for each event selection step in the search for a neutral heavy large-width pseudoscalar resonance that decays to a 400 GeV scalar resonance and a Z boson, which decays to neutrinos.

Acceptance times efficiency as a function of pseudoscalar resonant mass for each event selection step in the search for a neutral heavy narrow-width pseudoscalar resonance that decays to a 260 GeV scalar resonance and a Z boson, which decays to charged leptons.

Acceptance times efficiency as a function of pseudoscalar resonant mass for each event selection step in the search for a neutral heavy narrow-width pseudoscalar resonance that decays to a 300 GeV scalar resonance and a Z boson, which decays to charged leptons.

Acceptance times efficiency as a function of pseudoscalar resonant mass for each event selection step in the search for a neutral heavy narrow-width pseudoscalar resonance that decays to a 400 GeV scalar resonance and a Z boson, which decays to charged leptons.

Acceptance times efficiency as a function of pseudoscalar resonant mass for each event selection step in the search for a neutral heavy large-width pseudoscalar resonance that decays to a 260 GeV scalar resonance and a Z boson, which decays to charged leptons.

Acceptance times efficiency as a function of pseudoscalar resonant mass for each event selection step in the search for a neutral heavy large-width pseudoscalar resonance that decays to a 300 GeV scalar resonance and a Z boson, which decays to charged leptons.

Acceptance times efficiency as a function of pseudoscalar resonant mass for each event selection step in the search for a neutral heavy large-width pseudoscalar resonance that decays to a 400 GeV scalar resonance and a Z boson, which decays to charged leptons.

Expected and observed 95% CL upper limits on the cross-section of resonant $H\to 4b$ production in association with a W boson. The $\pm 1 \sigma$ and $\pm 2 \sigma$ uncertainty ranges for the expected limits are shown.

Expected and observed 95% CL upper limits on the cross-section of resonant $H\to 4b$ production in association with a Z boson. The $\pm 1 \sigma$ and $\pm 2 \sigma$ uncertainty ranges for the expected limits are shown.

Expected 95% CL upper limits on the cross-section of a heavy narrow-width pseudoscalar resonance decaying to a Z boson and a heavy scalar resonance decaying to $H\to 4b$. The $\pm 1 \sigma$ and $\pm 2 \sigma$ uncertainty ranges for the expected limits are shown.

Observed 95% CL upper limits on the cross-section of a heavy narrow-width pseudoscalar resonance decaying to a Z boson and a heavy scalar resonance decaying to $H\to 4b$. The $\pm 1 \sigma$ and $\pm 2 \sigma$ uncertainty ranges for the expected limits are shown.

Expected 95% CL upper limits on the cross-section of a heavy large-width pseudoscalar resonance decaying to a Z boson and a heavy scalar resonance decaying to $H\to 4b$. The $\pm 1 \sigma$ and $\pm 2 \sigma$ uncertainty ranges for the expected limits are shown.

Observed 95% CL upper limits on the cross-section of a heavy large-width pseudoscalar resonance decaying to a Z boson and a heavy scalar resonance decaying to $H\to 4b$. The $\pm 1 \sigma$ and $\pm 2 \sigma$ uncertainty ranges for the expected limits are shown.

Data and post-fit signal and background from S+B fit for 315 GeV resonant $H\to 4b$ production in association with a W boson.

Data and post-fit signal and background from S+B fit for 400 GeV resonant $H\to 4b$ production in association with a W boson.

Data and post-fit signal and background from S+B fit for 550 GeV resonant $H\to 4b$ production in association with a Z boson.

Data and post-fit signal and background from S+B fit for 400 GeV resonant $H\to 4b$ production in association with a Z boson.

Data and post-fit signal and background from S+B fit for a 790 GeV narrow-width pseudoscalar resonance decaying to a Z boson and a 300 GeV scalar resonance decaying to $H\to 4b$.

Data and post-fit signal and background from S+B fit for a 420 GeV large-width pseudoscalar resonance decaying to a Z boson and a 320 GeV scalar resonance decaying to $H\to 4b$.

Data and post-fit signal and background from S+B fit for a 700 GeV large-width pseudoscalar resonance decaying to a Z boson and a 380 GeV scalar resonance decaying to $H\to 4b$.

Data and post-fit signal and background from S+B fit for SM VHH production, with each Higgs boson decaying to $2b$.

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$