Measured differential cross section with associated uncertainties as a function of $N$($b$-jets). Each systematic uncertainty source is fully uncorrelated with the other sources and fully correlated across bins, except for the background modelling systematics for which an uncorrelated treatment across bins is more appropriate. For the meaning of the bins please refer to the figure.

Measured differential cross section with associated uncertainties as a function of $p_{T}^{j1}$. Each systematic uncertainty source is fully uncorrelated with the other sources and fully correlated across bins, except for the background modelling systematics for which an uncorrelated treatment across bins is more appropriate.

Measured differential cross section with associated uncertainties as a function of $p_{T}^{\gamma\gamma} (N_\mathrm{jets}^{p_T>30\ \mathrm{GeV}}=0)$. Each systematic uncertainty source is fully uncorrelated with the other sources and fully correlated across bins, except for the background modelling systematics for which an uncorrelated treatment across bins is more appropriate.

Measured differential cross section with associated uncertainties as a function of $m_{jj}$. Each systematic uncertainty source is fully uncorrelated with the other sources and fully correlated across bins, except for the background modelling systematics for which an uncorrelated treatment across bins is more appropriate.

Measured differential cross section with associated uncertainties as a function of $\Delta\phi_{jj}$. Each systematic uncertainty source is fully uncorrelated with the other sources and fully correlated across bins, except for the background modelling systematics for which an uncorrelated treatment across bins is more appropriate.

Measured differential cross section with associated uncertainties as a function of $p_{T}^{\gamma\gamma}$ in bins of $|y_{\gamma\gamma}|$. Each systematic uncertainty source is fully uncorrelated with the other sources and fully correlated across bins, except for the background modelling systematics for which an uncorrelated treatment across bins is more appropriate.

Measured differential cross section with associated uncertainties as a function of $\Delta\phi_{jj}$ in the VBF fiducial region. Each systematic uncertainty source is fully uncorrelated with the other sources and fully correlated across bins, except for the background modelling systematics for which an uncorrelated treatment across bins is more appropriate.

Observed statistical correlations, evaluated with a bootstrapping technique, between $p_\mathrm{T}^{\gamma\gamma}$, $N_\mathrm{jets}$, $m_{jj}$, $\Delta\phi_{jj}$, and $p_\mathrm{T}^{j1}$

Post-fit correlation matrix for the differential cross section measured as a function of $p_{T}^{\gamma\gamma}$.

Post-fit correlation matrix for the differential cross section measured as a function of $N_\mathrm{jets}$.

Post-fit correlation matrix for the differential cross section measured as a function of $N$($b$-jets).

Post-fit correlation matrix for the differential cross section measured as a function of $p_{T}^{j1}$.

Post-fit correlation matrix for the differential cross section measured as a function of $p_{T}^{\gamma\gamma} (N_\mathrm{jets}^{p_T>30\ \mathrm{GeV}}=0)$.

Post-fit correlation matrix for the differential cross section measured as a function of $m_{jj}$.

Post-fit correlation matrix for the differential cross section measured as a function of $\Delta\phi_{jj}$.

Post-fit correlation matrix for the differential cross section measured as a function of $\Delta\phi_{jj}$ in the VBF fiducial region.

Post-fit correlation matrix for the differential cross section measured as a function of $p_{T}^{\gamma\gamma}$ in bins of $|y_{\gamma\gamma}|$. The bins are the same as in the corresponding plot of the differential cross section.

Measured differential cross section with associated uncertainties as a function of $p_{T}^{\gamma1}/m_{\gamma\gamma}$. Each systematic uncertainty source is fully uncorrelated with the other sources and fully correlated across bins, except for the background modelling systematics for which an uncorrelated treatment across bins is more appropriate.

Post-fit correlation matrix for the differential cross section measured as a function of $p_{T}^{\gamma1}/m_{\gamma\gamma}$.

Measured differential cross section with associated uncertainties as a function of $p_{T}^{\gamma2}/m_{\gamma\gamma}$. Each systematic uncertainty source is fully uncorrelated with the other sources and fully correlated across bins, except for the background modelling systematics for which an uncorrelated treatment across bins is more appropriate.

Post-fit correlation matrix for the differential cross section measured as a function of $p_{T}^{\gamma2}/m_{\gamma\gamma}$.

Measured differential cross section with associated uncertainties as a function of $|y_{\gamma\gamma}|$. Each systematic uncertainty source is fully uncorrelated with the other sources and fully correlated across bins, except for the background modelling systematics for which an uncorrelated treatment across bins is more appropriate.

Post-fit correlation matrix for the differential cross section measured as a function of $|y_{\gamma\gamma}|$.

Measured differential cross section with associated uncertainties as a function of $m_{\gamma\gamma j}$. Each systematic uncertainty source is fully uncorrelated with the other sources and fully correlated across bins, except for the background modelling systematics for which an uncorrelated treatment across bins is more appropriate.

Post-fit correlation matrix for the differential cross section measured as a function of $m_{\gamma\gamma j}$.

Measured differential cross section with associated uncertainties as a function of $p_{T}^{\gamma\gamma j}$. Each systematic uncertainty source is fully uncorrelated with the other sources and fully correlated across bins, except for the background modelling systematics for which an uncorrelated treatment across bins is more appropriate.

Post-fit correlation matrix for the differential cross section measured as a function of $p_{T}^{\gamma\gamma j}$.

Measured differential cross section with associated uncertainties as a function of $H_{T}$. Each systematic uncertainty source is fully uncorrelated with the other sources and fully correlated across bins, except for the background modelling systematics for which an uncorrelated treatment across bins is more appropriate.

Post-fit correlation matrix for the differential cross section measured as a function of $H_{T}$.

Measured differential cross section with associated uncertainties as a function of $\tau_{C,j1}$. Each systematic uncertainty source is fully uncorrelated with the other sources and fully correlated across bins, except for the background modelling systematics for which an uncorrelated treatment across bins is more appropriate.

Post-fit correlation matrix for the differential cross section measured as a function of $\tau_{C,j1}$.

Measured differential cross section with associated uncertainties as a function of $\sum\tau_{C,j}$. Each systematic uncertainty source is fully uncorrelated with the other sources and fully correlated across bins, except for the background modelling systematics for which an uncorrelated treatment across bins is more appropriate.

Post-fit correlation matrix for the differential cross section measured as a function of $\sum\tau_{C,j}$.

Measured differential cross section with associated uncertainties as a function of $p_{T}^{\gamma\gamma} (N_\mathrm{jets}^{p_T>40\ \mathrm{GeV}}=0)$. Each systematic uncertainty source is fully uncorrelated with the other sources and fully correlated across bins, except for the background modelling systematics for which an uncorrelated treatment across bins is more appropriate.

Post-fit correlation matrix for the differential cross section measured as a function of $p_{T}^{\gamma\gamma} (N_\mathrm{jets}^{p_T>40\ \mathrm{GeV}}=0)$.

Measured differential cross section with associated uncertainties as a function of $p_{T}^{\gamma\gamma} (N_\mathrm{jets}^{p_T>50\ \mathrm{GeV}}=0)$. Each systematic uncertainty source is fully uncorrelated with the other sources and fully correlated across bins, except for the background modelling systematics for which an uncorrelated treatment across bins is more appropriate.

Post-fit correlation matrix for the differential cross section measured as a function of $p_{T}^{\gamma\gamma} (N_\mathrm{jets}^{p_T>50\ \mathrm{GeV}}=0)$.

Measured differential cross section with associated uncertainties as a function of $p_{T}^{\gamma\gamma} (N_\mathrm{jets}^{p_T>60\ \mathrm{GeV}}=0)$. Each systematic uncertainty source is fully uncorrelated with the other sources and fully correlated across bins, except for the background modelling systematics for which an uncorrelated treatment across bins is more appropriate.

Post-fit correlation matrix for the differential cross section measured as a function of $p_{T}^{\gamma\gamma} (N_\mathrm{jets}^{p_T>60\ \mathrm{GeV}}=0)$.

Measured differential cross section with associated uncertainties as a function of $\pi - |\Delta\phi_{\gamma\gamma,jj}|$. Each systematic uncertainty source is fully uncorrelated with the other sources and fully correlated across bins, except for the background modelling systematics for which an uncorrelated treatment across bins is more appropriate.

Post-fit correlation matrix for the differential cross section measured as a function of $\pi - |\Delta\phi_{\gamma\gamma,jj}|$.

Measured differential cross section with associated uncertainties as a function of $p_{T}^{\gamma\gamma jj}$. Each systematic uncertainty source is fully uncorrelated with the other sources and fully correlated across bins, except for the background modelling systematics for which an uncorrelated treatment across bins is more appropriate.

Post-fit correlation matrix for the differential cross section measured as a function of $p_{T}^{\gamma\gamma jj}$.

Measured differential cross section with associated uncertainties as a function of $p_{T}^{\gamma\gamma}$ in bins of $p_{T}^{\gamma\gamma j}$. Each systematic uncertainty source is fully uncorrelated with the other sources and fully correlated across bins, except for the background modelling systematics for which an uncorrelated treatment across bins is more appropriate.

Post-fit correlation matrix for the differential cross section measured as a function of $p_{T}^{\gamma\gamma}$ in bins of $p_{T}^{\gamma\gamma j}$. The bins are the same as in the corresponding plot of the differential cross section.

Measured differential cross section with associated uncertainties as a function of $p_{T}^{\gamma\gamma}$ in bins of $\tau_{C,j1}$. Each systematic uncertainty source is fully uncorrelated with the other sources and fully correlated across bins, except for the background modelling systematics for which an uncorrelated treatment across bins is more appropriate.

Post-fit correlation matrix for the differential cross section measured as a function of $p_{T}^{\gamma\gamma}$ in bins of $\tau_{C,j1}$. The bins are the same as in the corresponding plot of the differential cross section.

Measured differential cross section with associated uncertainties as a function of $(p_{T}^{\gamma1}-p_{T}^{\gamma2})/m_{\gamma\gamma}$ in bins of $(p_{T}^{\gamma1}+p_{T}^{\gamma2})/m_{\gamma\gamma}$. Each systematic uncertainty source is fully uncorrelated with the other sources and fully correlated across bins, except for the background modelling systematics for which an uncorrelated treatment across bins is more appropriate.

Post-fit correlation matrix for the differential cross section measured as a function of $(p_{T}^{\gamma1}-p_{T}^{\gamma2})/m_{\gamma\gamma}$ in bins of $(p_{T}^{\gamma1}+p_{T}^{\gamma2})/m_{\gamma\gamma}$. The bins are the same as in the corresponding plot of the differential cross section.

Measured differential cross section with associated uncertainties as a function of $|\eta^*|$ in the VBF fiducial region. Each systematic uncertainty source is fully uncorrelated with the other sources and fully correlated across bins, except for the background modelling systematics for which an uncorrelated treatment across bins is more appropriate.

Post-fit correlation matrix for the differential cross section measured as a function of $|\eta^*|$ in the VBF fiducial region.

Measured differential cross section with associated uncertainties as a function of $p_{T}^{\gamma\gamma jj}$ in the VBF fiducial region. Each systematic uncertainty source is fully uncorrelated with the other sources and fully correlated across bins, except for the background modelling systematics for which an uncorrelated treatment across bins is more appropriate.

Post-fit correlation matrix for the differential cross section measured as a function of $p_{T}^{\gamma\gamma jj}$ in the VBF fiducial region.

Measured differential cross section with associated uncertainties as a function of $p_{T}^{j1}$ in the VBF fiducial region. Each systematic uncertainty source is fully uncorrelated with the other sources and fully correlated across bins, except for the background modelling systematics for which an uncorrelated treatment across bins is more appropriate.

Post-fit correlation matrix for the differential cross section measured as a function of $p_{T}^{j1}$ in the VBF fiducial region.

Measured differential cross section with associated uncertainties as a function of $p_{T}^{j1}$ in bins of $\Delta\phi_{jj}$ in the VBF fiducial region. Each systematic uncertainty source is fully uncorrelated with the other sources and fully correlated across bins, except for the background modelling systematics for which an uncorrelated treatment across bins is more appropriate.

Post-fit correlation matrix for the differential cross section measured as a function of $p_{T}^{j1}$ in bins of $\Delta\phi_{jj}$ in the VBF fiducial region. The bins are the same as in the corresponding plot of the differential cross section.

The Au+Au 200 GeV measurements of the two-, four- and six-particle elliptic flow and the elliptic flow fluctuations.

The Au+Au 54.4 GeV measurements of the two-, four- and six-particle elliptic flow and the elliptic flow fluctuations.

The Au+Au 39 GeV measurements of the two-, four- and six-particle elliptic flow and the elliptic flow fluctuations.

The Au+Au 27 GeV measurements of the two-, four- and six-particle elliptic flow and the elliptic flow fluctuations.

The Au+Au 19.6 GeV measurements of the two- and four-particle elliptic flow and the elliptic flow fluctuations.

The Au+Au 11.5 GeV measurements of the two- and four-particle elliptic flow and the elliptic flow fluctuations.

The U+U 193 GeV measurements of the two-, four- and six-particle elliptic flow and the elliptic flow fluctuations.

The Cu+Au 200 GeV measurements of the two-, four- and six-particle elliptic flow and the elliptic flow fluctuations.

Observed yields in the boost_1_ge2J signal region category of the hh channel.

Observed yields in the boost_2 signal region category of the hh channel.

Observed yields in the boost_3 signal region category of the hh channel.

Observed yields in the ttH_0 signal region category of the hh channel.

Observed yields in the ttH_1 signal region category of the hh channel.

Observed yields in the VH_0 signal region category of the hh channel.

Observed yields in the VH_1 signal region category of the hh channel.

Observed yields in the VBF_0 signal region category of the hh channel.

Observed yields in the VBF_1 signal region category of the hh channel.

Observed yields in the boost_0_1J signal region category of the lh channel.

Observed yields in the boost_0_ge2J signal region category of the lh channel.

Observed yields in the boost_1_1J signal region category of the lh channel.

Observed yields in the boost_1_ge2J signal region category of the lh channel.

Observed yields in the boost_2 signal region category of the lh channel.

Observed yields in the boost_3 signal region category of the lh channel.

Observed yields in the VH_0 signal region category of the lh channel.

Observed yields in the VH_1 signal region category of the lh channel.

Observed yields in the VBF_0 signal region category of the lh channel.

Observed yields in the VBF_1 signal region category of the lh channel.

Observed yields in the boost_0_1J signal region category of the ll channel.

Observed yields in the boost_0_ge2J signal region category of the ll channel.

Observed yields in the boost_1_1J signal region category of the ll channel.

Observed yields in the boost_1_ge2J signal region category of the ll channel.

Observed yields in the boost_2 signal region category of the ll channel.

Observed yields in the boost_3 signal region category of the ll channel.

Observed yields in the VH_0 signal region category of the ll channel.

Observed yields in the VH_1 signal region category of the ll channel.

Observed yields in the VBF_0 signal region category of the ll channel.

Observed yields in the VBF_1 signal region category of the ll channel.

Best-fit values and uncertainties for the pp $\rightarrow H\rightarrow\tau\tau$ cross-section measurement and the measurement in the four dominant production modes. All measurements include the branching ratio of $H\rightarrow\tau\tau$ and are performed with true Higgs boson rapidity $|y_{H}|<2.5$. The SM predictions for each region, computed using the inclusive cross-section calculations and the simulated event samples are also shown. The contributions to the total uncertainty in the measurements from statistical (Stat. unc.) or systematic uncertainties (Syst. unc.) are given separately.

The measured correlations between each parameter of interest in the measurement of the cross-sections per production mode

Best-fit values and uncertainties for the $H\rightarrow\tau\tau$ cross-sections, in the reduced stage 1.2 STXS scheme. The EW production mode includes vector-boson fusion and qq$\rightarrow$V($\rightarrow$qq)H processes. All measurements include the branching ratio of $H\rightarrow\tau\tau$ and are performed with true Higgs boson rapidity $|y_{H}|<2.5$. The SM predictions for each region, computed using the inclusive cross-section calculations and the simulated event samples are also shown. The contributions to the total uncertainty in the measurements from statistical (Stat. unc.) or systematic uncertainties (Syst. unc.) are given separately. The asterisk symbol (*) indicates that the criteria for m$_{jj}$ only apply to events with at least two reconstructed jets.

The measured correlations between each pair of parameters of interest in the STXS measurement. The asterisk symbol (*) indicates that the criteria for mjj only apply to events with at least two reconstructed jets.

Observed and expected 95% CL lower limits on the $T$-quark mass as a function of the $T$-quark width-to-mass ratio and the branching fraction of the $T \rightarrow Ht$ decay ($\Gamma_{T}$ is the $T$-quark width).

Cutflow table listing the number of events passing each criterion for a $T$-quark hypothesis with a mass of 1.6 TeV and $\kappa_{T} = 0.5$. The initial signal event yield is the predicted number of $T$-quark events inclusive in the Higgs-boson and top-quark decays for 139 fb$^{-1}$.

Observed 95% CL upper limits on the single $T$-quark production cross-section as a function of the $T$-quark coupling $\kappa_{T}$ and $m_{T}$.

Expected 95% CL upper limits on the single $T$-quark production cross-section as a function of the $T$-quark coupling $\kappa_{T}$ and $m_{T}$.

Observed and expected 95% CL lower limits on the $T$-quark mass as a function of the $T$-quark width-to-mass ratio and the branching fraction of the $T \rightarrow Wb$ decay ($\Gamma_{T}$ is the $T$-quark width).

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

The $1\sigma_{\text{exp}}$ variation of the expected exclusion contour at 95% CL as a function of the smuon and neutralino masses, for $\lambda_{231}^{'}=1.0$.

The observed exclusion contour at 95% CL as a function of the smuon and neutralino masses, for $\lambda_{231}^{'}=0.1$.

The expected exclusion contour at 95% CL as a function of the smuon and neutralino masses, for $\lambda_{231}^{'}=0.1$.

The observed exclusion contour at 95% CL as a function of the smuon and neutralino masses, for $\lambda_{231}^{'}=0.15$.

The expected exclusion contour at 95% CL as a function of the smuon and neutralino masses, for $\lambda_{231}^{'}=0.15$.

The observed exclusion contour at 95% CL as a function of the smuon and neutralino masses, for $\lambda_{231}^{'}=0.2$.

The expected exclusion contour at 95% CL as a function of the smuon and neutralino masses, for $\lambda_{231}^{'}=0.2$.

The observed exclusion contour at 95% CL as a function of the smuon and neutralino masses, for $\lambda_{231}^{'}=0.4$.

The expected exclusion contour at 95% CL as a function of the smuon and neutralino masses, for $\lambda_{231}^{'}=0.4$.

The observed exclusion contour at 95% CL as a function of the smuon and neutralino masses, for $\lambda_{231}^{'}=0.6$.

The expected exclusion contour at 95% CL as a function of the smuon and neutralino masses, for $\lambda_{231}^{'}=0p6$.

The observed exclusion contour at 95% CL as a function of the smuon and neutralino masses, for $\lambda_{231}^{'}=1.5$.

The expected exclusion contour at 95% CL as a function of the smuon and neutralino masses, for $\lambda_{231}^{'}=1.5$.

The observed exclusion contour at 95% CL as a function of the leptoquark mass and coupling strength.

The expected exclusion contour at 95% CL as a function of the leptoquark mass and coupling strength.

The minus $1\sigma_{\text{theory}}$ variation of the observed exclusion contour at 95% CL as a function of the leptoquark mass and coupling strength.

The plus $1\sigma_{\text{theory}}$ variation of the observed exclusion contour at 95% CL as a function of the leptoquark mass and coupling strength.

The $1\sigma_{\text{exp}}$ variation of the expected exclusion contour at 95% CL as a function of the leptoquark mass and coupling strength.

Observed yields, and fake lepton background yields in the $e^{+}\mu^{-}$ and $e^{-}\mu^{+}$ channels of SR-MET, along with the results of the $e^{+}\mu^{-}/e^{-}\mu^{+}$ ratio measurement and 1-sided p-value in SR-MET, binned in $M_{T2}$.

Observed yields, and fake lepton background yields in the $e^{+}\mu^{-}$ and $e^{-}\mu^{+}$ channels of SR-JET, along with the results of the $e^{+}\mu^{-}/e^{-}\mu^{+}$ ratio measurement and 1-sided p-value in SR-JET, binned in $H_{\text{P}}$.

Observed and expected 95% CL upper limits on the total number of signal events entering the $e^{+}\mu^{-}$ and $e^{-}\mu^{+}$ channels of each bin of SR-MET. The regions are binned in the same way as the ratio $\rho$ measurement. The limits are shown for a selection of 'z' values, where 'z' is the fraction of the total signal events entering the $e^{+}\mu^{-}$ channel.

Observed and expected 95% CL upper limits on the total number of signal events entering the $e^{+}\mu^{-}$ and $e^{-}\mu^{+}$ channels of each bin of SR-JET. The regions are binned in the same way as the ratio $\rho$ measurement. The limits are shown for a selection of 'z' values, where 'z' is the fraction of the total signal events entering the $e^{+}\mu^{-}$ channel.

Signal yields following each cut in the analysis, for representative $R$-parity-violating supersymmetry and leptoquark signals. All yields are MC generator-weighted and normalised to $139~\text{fb}^{-1}$. The cut labelled `Preselection' includes trigger requirements, and requires exactly one Baseline electron and one Baseline muon. At this point, the muon charge-bias correction weights are also applied. The $R$-parity-violating supersymmetry models were generated by specifying a top-quark in the final state and applying a two-lepton filter, hence the first row also includes events where the top quark decays to an electron.

The efficiency for simulated VBF events to pass each analysis cut.

The efficiency for simulated VH events to pass each analysis cut.

The efficiency for simulated ttH events to pass each analysis cut.

The fractional contribution of each production mode to a given analysis region around the Higgs boson peak, defined as 105 < mJ < 140 GeV. The fraction is given with respect to the total signal yield in the analysis region in question.

Signal acceptance times efficiency within the fiducial volume used in the fiducial region.

Signal acceptance times efficiency for the STXS volumes in the differential measurement. Along with the pTH requirements shown, |yH | < 2 is required. For events with pTH < 300 GeV, the acceptance times efficiency is less than 0.1 × 10−2.

Event yields and associated uncertainties after the global likelihood fit in the pT(H) > 450 GeV fiducial region. The total background yield can differ from the sum of the individual components due to rounding.

Event yields and associated uncertainties after the global likelihood fit of the differential regions. The total background yield can differ from the sum of the individual components due to rounding. The five STXS volumes are labeled pT0–pT4 corresponding to pTH < 300 GeV, 300 – 450 GeV, 450 – 650 GeV, 650 – 1000 GeV, and > 1000 GeV, respectively. The pT0 event yield is constrained to its SM value within the theoretical and experimental uncertainties and free parameters act independently on the remaining four volumes.

The observed exclusion contours as a function of (m(med), m($\chi$)), with Vector mediator)

The expected exclusion contours as a function of (m(a), tan($\beta$)), with sin($\theta$) = 0.35)

The observed exclusion contours as a function of (m(a), tan($\beta$)), with sin($\theta$) = 0.35)

The expected exclusion contours as a function of (m(a), tan($\beta$)), with sin($\theta$) = 0.7)

The observed exclusion contours as a function of (m(a), tan($\beta$)), with sin($\theta$) = 0.7)

The expected exclusion contours as a function of (m(H), tan($\beta$)), with sin($\theta$) = 0.35)

The observed exclusion contours as a function of (m(H), tan($\beta$)), with sin($\theta$) = 0.35)

The expected exclusion contours as a function of (m(H), tan($\beta$)), with sin($\theta$) = 0.7)

The observed exclusion contours as a function of (m(H), tan($\beta$)), with sin($\theta$) = 0.7)

The expected exclusion contours as a function of (m(a), m(H)), with sin($\theta$) = 0.35)

The observed exclusion contours as a function of (m(a), m(H)), with sin($\theta$) = 0.35)

The expected exclusion contours as a function of (m(a), m(H)), with sin($\theta$) = 0.7)

The observed exclusion contours as a function of (m(a), m(H)), with sin($\theta$) = 0.7)

Expected lower limit on signal strength at 95% CL as a function of sin($\theta$), with m(a) = 200 GeV, m(H) = 600 GeV.

Observed lower limit on signal strength at 95% CL as a function of sin($\theta$), with m(a) = 200 GeV, m(H) = 600 GeV.

Expected lower limit on signal strength at 95% CL as a function of sin($\theta$), with m(a) = 350 GeV, m(H) = 1000 GeV.

Observed lower limit on signal strength at 95% CL as a function of sin($\theta$), with m(a) = 350 GeV, m(H) = 1000 GeV.

Observed lower limit on WIMP-nucleon cross section at 90% CL as a function of m(WIMP), assuming Higgs-portal scenario with Scalar WIMP.

Observed lower limit on WIMP-nucleon cross section at 90% CL as a function of m(WIMP), assuming Higgs-portal scenario with Majorana WIMP.

Observed lower limit on the spin-dependent WIMP–proton scattering cross-section.

Observed lower limit on the spin-independent WIMP–nucleon scattering cross-section.

Cutflow of unweighted and weighted events of $ZH$ signals in the electron channel.

Cutflow of unweighted and weighted events of $ZH$ signals in the muon channel.

Cutflow of unweighted and weighted events of DM signals (simplified DM axial-vector with m(med) = 900 GeV and m($\chi$) = 40 GeV, 2HDM+$a$ signal with m(A) = 600 GeV, m(a) = 400 GeV, tan($\beta$) = 1.0 and m($\chi$) = 10 GeV) in the electron channel.

Cutflow of unweighted and weighted events of DM signals (simplified DM axial-vector with m(med) = 900 GeV and m($\chi$) = 40 GeV, 2HDM+$a$ signal with m(A) = 600 GeV, m(a) = 400 GeV, tan($\beta$) = 1.0 and m($\chi$) = 10 GeV) in the muon channel.

The JER for R = 0.2 jets in Pb+Pb collisions as a function of $2|\Psi_2-\phi^{reco}|$ for centrality selections of 0-5%, 5-10%, 10-20%, 20-40% and 40-60%.

The systematic uncertainties in v2 for 20-40% centrality Pb+Pb collisions as a function of pT.

The systematic uncertainties in v2 for 5-10% centrality Pb+Pb collisions as a function of pT.

The systematic uncertainties in v3 for 20-40% centrality Pb+Pb collisions as a function of pT.

The systematic uncertainties in v3 for 5-10% centrality Pb+Pb collisions as a function of pT.

The systematic uncertainties in v4 for 20-40% centrality Pb+Pb collisions as a function of pT.

The systematic uncertainties in v4 for 5-10% centrality Pb+Pb collisions as a function of pT.

The systematic uncertainties in v2 for pT = 71--398 GeV jets as a function of centrality.

The systematic uncertainties in v3 for pT = 71--398 GeV jets as a function of centrality.

The systematic uncertainties in v4 for pT = 71--398 GeV jets as a function of centrality.

Angular distribution of jets as a function of the observed psi2 plane for jets with 71 < pT < 79 GeV in the 10-20% centrality bin.

Angular distribution of jets as a function of the observed psi3 plane for jets with 71 < pT < 79 GeV in the 10-20% centrality bin.

Angular distribution of jets as a function of the observed psi4 plane for jets with 71 < pT < 79 GeV in the 10-20% centrality bin.

The v2 values for R = 0.2 jets as a function of centrality for jets in several pT ranges.

The v2 values for R = 0.2 jets as a function of pT for 0-5%, 5-10%, and 20-40% centrality collisions.

The v2, v3, and v4 as a function of centrality for jets with pT = 71-398 GeV.

The v3 values for R = 0.2 jets as a function of centrality for jets in several pT ranges.

The v4 values for R = 0.2 jets as a function of centrality for jets in several pT ranges.

The v2 as a function of pT for jets in 10-20% centrality collisions.

The v3 as a function of pT for jets in 10-20% centrality collisions.

The v2 as a function of pT for jets in 20-40% centrality collisions.

The v3 as a function of pT for jets in 20-40% centrality collisions.

The v2 for jets in 10-20% centrality collisions.

The v3 for jets in 10-20% centrality collisions.

R2max as a function of pT (filled circles). Also shown is 1 - 4v2/(1+2v2).

R3max as a function of pT (filled circles). Also shown is 1 - 4v3/(1+2v3).

The systematic uncertainties in v2 for 40-60% centrality Pb+Pb collisions as a function of pT.

The systematic uncertainties in v3 for 40-60% centrality Pb+Pb collisions as a function of pT.

The systematic uncertainties in v4 for 40-60% centrality Pb+Pb collisions as a function of pT.

The systematic uncertainties in v2 for 10-20% centrality Pb+Pb collisions as a function of pT.

The systematic uncertainties in v3 for 10-20% centrality Pb+Pb collisions as a function of pT.

The systematic uncertainties in v4 for 10-20% centrality Pb+Pb collisions as a function of pT.

The systematic uncertainties in v2 for 0-5% centrality Pb+Pb collisions as a function of pT.

The systematic uncertainties in v3 for 0-5% centrality Pb+Pb collisions as a function of pT.

The systematic uncertainties in v4 for 0-5% centrality Pb+Pb collisions as a function of pT.

The v2 as a function of pT for jets in 40-60% centrality collisions.

The v2 as a function of pT for jets in 5-10% centrality collisions.

The v2 as a function of pT for jets in 0-5% centrality collisions.

The v3 as a function of pT for jets in 40-60% centrality collisions.

The v3 as a function of pT for jets in 5-10% centrality collisions.

The v3 as a function of pT for jets in 0-5% centrality collisions.

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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