95% CL observed upper limits on the Yukawa couplings

The correlation matrix for the Higgs pT measurements, both for the unregularized and regularized fits. The last bin is inclusive.

The correlation matrix for the jet multiplicity measurements, both for the unregularized and regularized fits. The last bin is inclusive.

The correlation matrix for the leading jet pT measurements, both for the unregularized and regularized fits. The last bin is inclusive.

The fiducial integrated signal strength and cross section.

Upper limits on the production cross section times branching fraction of the b* RH hypothesis at a 95% CL. Dashed colored lines show the expected limits from the l+jets and all-hadronic channels, where the latter start at resonance masses of 1.4 TeV. The observed and expected limits from the combination are shown as solid and dashed black lines, respectively. The green and yellow bands show the 68 and 95% confidence intervals on the combined expected limits.

Upper limits on the production cross section times branching fraction of the b* VL hypothesis at a 95% CL. Dashed colored lines show the expected limits from the l+jets and all-hadronic channels, where the latter start at resonance masses of 1.4 TeV. The observed and expected limits from the combination are shown as solid and dashed black lines, respectively. The green and yellow bands show the 68 and 95% confidence intervals on the combined expected limits.

Upper limits on the production cross section times branching fraction of the B+b hypothesis at a 95% CL. Dashed colored lines show the expected limits from the l+jets and all-hadronic channels, where the latter start at resonance masses of 1.4 TeV. The observed and expected limits from the combination are shown as solid and dashed black lines, respectively. The green and yellow bands show the 68 and 95% confidence intervals on the combined expected limits.

Upper limits on the production cross section times branching fraction of the B+t hypothesis at a 95% CL. Dashed colored lines show the expected limits from the l+jets and all-hadronic channels, where the latter start at resonance masses of 1.4 TeV. The observed and expected limits from the combination are shown as solid and dashed black lines, respectively. The green and yellow bands show the 68 and 95% confidence intervals on the combined expected limits.

The correlation matrix for the njet measurement. The last bin is inclusive.

The fiducial integrated signal strength and cross section The uncertainty breakdown is given in terms of statistical (stat), experimental (exp), theoretical uncertainties on the background (bkg) and on the signal (sig), and the luminosity uncertainty (lumi). The fiducial cross section and its full uncertainty in each bin are also given.

Summary of the sources of systematic and statistical uncertainties and their impact on the measurement of the ttH and tH signal rates, and the measured value of the unconstrained nuisance parameters. The quantity $\Delta r_{x}/r_{x}$ corresponds to the change in uncertainty when fixing the nuisance parameters associated with that uncertainty in the fit. Under the label "MC and sideband statistical uncertainty" are the uncertainties associated with the limited number of simulated MC events and the amount of data events in the application region of the MP method.

Number of events selected in the 2$\ell$ss+0$\tau_{\mathrm{h}}$, 3$\ell$+0$\tau_{\mathrm{h}}$, and 2$\ell$ss+1$\tau_{\mathrm{h}}$ analysis channels compared to the event yields expected from the $ tH$ and $ H$ signals and from background processes. The expected event yields are computed for the values of nuisance parameters and of the POI obtained from the ML fit. The best fit values of the POI amount to $\hat{r_{ttH}} = 0.92$ and $\hat{r_{tH}} = 5.7$. Quoted uncertainties represent the sum of statistical and systematic components.

Number of events selected in the 2$\ell$ss+0$\tau_{\mathrm{h}}$, 3$\ell$+0$\tau_{\mathrm{h}}$, and 2$\ell$ss+1$\tau_{\mathrm{h}}$ analysis channels compared to the event yields expected from the $ tH$ and $ H$ signals and from background processes. The expected event yields are computed for the values of nuisance parameters and of the POI obtained from the ML fit. The best fit values of the POI amount to $\hat{r_{ttH}} = 0.92$ and $\hat{r_{tH}} = 5.7$. Quoted uncertainties represent the sum of statistical and systematic components.

Number of events selected in the 1$\ell$+1$\tau_{\mathrm{h}}$, 0$\ell$+2$\tau_{\mathrm{h}}$, 2$\ell$os+1$\tau_{\mathrm{h}}$, and 1$\ell$+2$\tau_{\mathrm{h}}$ analysis channels compared to the event yields expected from the $ tH$ and $ H$ signals and from background processes. The expected event yields are computed for the values of nuisance parameters and of the POI obtained from the ML fit. The best fit values of the POI amount to $\hat{r_{ttH}} = 0.92$ and $\hat{r_{tH}} = 5.7$. Quoted uncertainties represent the sum of statistical and systematic components.

Number of events selected in the 1$\ell$+1$\tau_{\mathrm{h}}$, 0$\ell$+2$\tau_{\mathrm{h}}$, 2$\ell$os+1$\tau_{\mathrm{h}}$, and 1$\ell$+2$\tau_{\mathrm{h}}$ analysis channels compared to the event yields expected from the $ tH$ and $ H$ signals and from background processes. The expected event yields are computed for the values of nuisance parameters and of the POI obtained from the ML fit. The best fit values of the POI amount to $\hat{r_{ttH}} = 0.92$ and $\hat{r_{tH}} = 5.7$. Quoted uncertainties represent the sum of statistical and systematic components.

Number of events selected in the 4$\ell$+0$\tau_{\mathrm{h}}$, 3$\ell$+1$\tau_{\mathrm{h}}$, and 2$\ell$+2$\tau_{\mathrm{h}}$ analysis channels compared to the event yields expected from the $ tH$ and $ H$ signals and from background processes. The expected event yields are computed for the values of nuisance parameters and of the POI obtained from the ML fit. The best fit values of the POI amount to $\hat{r_{ttH}} = 0.92$ and $\hat{r_{tH}} = 5.7$. Quoted uncertainties represent the sum of statistical and systematic components.

Number of events selected in the 4$\ell$+0$\tau_{\mathrm{h}}$, 3$\ell$+1$\tau_{\mathrm{h}}$, and 2$\ell$+2$\tau_{\mathrm{h}}$ analysis channels compared to the event yields expected from the $ tH$ and $ H$ signals and from background processes. The expected event yields are computed for the values of nuisance parameters and of the POI obtained from the ML fit. The best fit values of the POI amount to $\hat{r_{ttH}} = 0.92$ and $\hat{r_{tH}} = 5.7$. Quoted uncertainties represent the sum of statistical and systematic components.

Probability for ttH signal events to pass the event selection in each of the analysis categories employed.

Observed and expected 95% CL upper limits on the product of the production cross section ($\sigma$) and the branching fraction, obtained after combining all categories with 138 $\mathrm{fb}^{−1}$ of data at $\sqrt{s}$ = 13 TeV for $\mathrm{V'}$ to $VV$ + $VH$ signal in HVT model C.

Total signal efficiency including all event categories as a function of the resonance mass $M_X$ for all DY/ggF signal models under consideration. The estimates from MC simulation (markers) are interpolated with a functional form (lines).

Total signal efficiency including all event categories as a function of the resonance mass $M_X$ for all VBF signal models under consideration. The estimates from MC simulation (markers) are interpolated with a functional form (lines).

Expected and observed upper limits on the AQGC operators $a^Z_C/\Lambda^2$, with no unitarization. The $y$ axis shows the limit on the ratio of the observed cross section to the cross section predicted for each anomalous coupling value ($\sigma_\mathrm{AQGC}$).

Limits on LEP-like dimension-6 anomalous quartic gauge coupling parameters, with and without unitarization via a clipping procedure.

Conversion of limits on $a^W_0$ to dimension-8 $f_{M,i}$ operators, using the assumption of vanishing $WWZ\gamma$ couplings to eliminate some parameters. When quoting limits on one of the operators, the other is fixed to zero. The results for $|f_{M,0}/\Lambda^{4}|$ and $|f_{M,4}/\Lambda^{4}|$ are shown with and without clipping of the signal model at 1.4 TeV, when the other parameter is fixed to the SM value of zero.

Conversion of limits on $a^W_0$ and $a^W_C$ to dimension-8 $f_{M,i}$ operators, using the assumption that all $f_{M,i}$ except one are equal to zero. The results are shown with and without clipping of the signal model at 1.4 TeV.

{Expected and observed limits in the two-dimensional plane of $a^W_0/\Lambda^2$ vs. $a^W_C/\Lambda^2$. The limits are described by analytical ellipses of equation $(x-x0)^2/a^2 + (y-y0)^2/b^2 = 1$ and rotated counter-clockwise by $\theta$ degrees, where $x$ and $y$ in the equation correspond to the $a_0^W$ and $a_C^W$ couplings, respectively.

{Expected and observed limits in the two-dimensional plane of $a^W_0/\Lambda^2$ vs. $a^W_C/\Lambda^2$ with unitarization imposed by clipping the signal model at 1.4 TeV. The limits are described by analytical ellipses of equation $(x-x0)^2/a^2 + (y-y0)^2/b^2 = 1$ and rotated counter-clockwise by $\theta$ degrees, where $x$ and $y$ in the equation correspond to the $a_0^W$ and $a_C^W$ couplings, respectively.

{Expected and observed limits in the two-dimensional plane of $a^Z_0/\Lambda^2$ vs. $a^Z_C/\Lambda^2$. The limits are described by analytical ellipses of equation $(x-x0)^2/a^2 + (y-y0)^2/b^2 = 1$ and rotated counter-clockwise by $\theta$ degrees, where $x$ and $y$ in the equation correspond to the $a_0^Z$ and $a_C^Z$ couplings, respectively.

Fits to the $m_{\ell^+\ell^-\gamma}$ data distribution in the dijet-1 categories. In the upper panel, the red solid line shows the result of a signal-plus-background fit to the given category. The red dashed line shows the background component of the fit. The green and yellow bands represent the 68 and 95% CL uncertainties in the fit. Also plotted is the expected SM signal, scaled by a factor of 10. In the lower panel, the data minus the background component of the fit is shown.

Fits to the $m_{\ell^+\ell^-\gamma}$ data distribution in the dijet-2 categories. In the upper panel, the red solid line shows the result of a signal-plus-background fit to the given category. The red dashed line shows the background component of the fit. The green and yellow bands represent the 68 and 95% CL uncertainties in the fit. Also plotted is the expected SM signal, scaled by a factor of 10. In the lower panel, the data minus the background component of the fit is shown.

Fits to the $m_{\ell^+\ell^-\gamma}$ data distribution in the dijet-3 categories. In the upper panel, the red solid line shows the result of a signal-plus-background fit to the given category. The red dashed line shows the background component of the fit. The green and yellow bands represent the 68 and 95% CL uncertainties in the fit. Also plotted is the expected SM signal, scaled by a factor of 10. In the lower panel, the data minus the background component of the fit is shown.

Fits to the $m_{\ell^+\ell^-\gamma}$ data distribution in the untagged 1 categories. In the upper panel, the red solid line shows the result of a signal-plus-background fit to the given category. The red dashed line shows the background component of the fit. The green and yellow bands represent the 68 and 95% CL uncertainties in the fit. Also plotted is the expected SM signal, scaled by a factor of 10. In the lower panel, the data minus the background component of the fit is shown.

Fits to the $m_{\ell^+\ell^-\gamma}$ data distribution in the untagged 2 categories. In the upper panel, the red solid line shows the result of a signal-plus-background fit to the given category. The red dashed line shows the background component of the fit. The green and yellow bands represent the 68 and 95% CL uncertainties in the fit. Also plotted is the expected SM signal, scaled by a factor of 10. In the lower panel, the data minus the background component of the fit is shown.

Fits to the $m_{\ell^+\ell^-\gamma}$ data distribution in the untagged 3 categories. In the upper panel, the red solid line shows the result of a signal-plus-background fit to the given category. The red dashed line shows the background component of the fit. The green and yellow bands represent the 68 and 95% CL uncertainties in the fit. Also plotted is the expected SM signal, scaled by a factor of 10. In the lower panel, the data minus the background component of the fit is shown.

Fits to the $m_{\ell^+\ell^-\gamma}$ data distribution in the untagged 4 categories. In the upper panel, the red solid line shows the result of a signal-plus-background fit to the given category. The red dashed line shows the background component of the fit. The green and yellow bands represent the 68 and 95% CL uncertainties in the fit. Also plotted is the expected SM signal, scaled by a factor of 10. In the lower panel, the data minus the background component of the fit is shown.

Observed signal strength ($\mu$) for a SM Higgs boson at 125.38 GeV. The labels untagged combined,dijet combined and combined represent the results obtained from simultaneous fits of the untagged categories, dijet categories, and full set of categories, respectively. The black solid line shows $\mu=1$, and the red dashed line shows the best fit value $\hat{\mu}=2.4$ of all categories combined.

Upper limit (95% CL) on the signal strength ($\mu$) relative to the SM prediction, as a function of the assumed value of the Higgs boson mass used in the fit.

The observed DNN output distribution for data collected in 2016 in the VBF-SB region compared to the post-fit background estimate from SM processes. The predicted backgrounds are obtained from a S+B fit performed across analysis regions and years. The total post-fit and pre-fit uncertainties on the background prediction are also reported.

The observed DNN output distribution for data collected in 2017 in the VBF-SB region compared to the post-fit background estimate from SM processes. The predicted backgrounds are obtained from a S+B fit performed across analysis regions and years. The total post-fit and pre-fit uncertainties on the background prediction are also reported.

The observed DNN output distribution for data collected in 2018 in the VBF-SB region compared to the post-fit background estimate from SM processes. The predicted backgrounds are obtained from a S+B fit performed across analysis regions and years. The total post-fit and pre-fit uncertainties on the background prediction are also reported.

The observed DNN output distribution in the VBF-SB regions for the combination of 2016, 2017, and 2018 data, compared to the post-fit prediction from SM processes. The total post-fit and pre-fit uncertainties on the background prediction are also reported.

The observed DNN output distribution in the VBF-SR region for the combination of 2016, 2017, and 2018 data, compared to the post-fit prediction from SM processes. The post-fit distributions for the Higgs boson signal produced via ggH and VBF modes with mass of 125.38 GeV, along with the total post-fit and pre-fit uncertainties on the background prediction, are also reported.

The observed BDT output distribution compared to the prediction from the simulation of various SM background processes. Dimuon events passing the event selection requirements of the ggH category, with dimuon mass between 110–150 GeV, are considered. The expected distributions for ggH, VBF, and other signal processes are also reported. The background uncertainty is obtained by summing in quadrature the uncertainty due to the limited size of the simulated samples with those arising from experimental and theoretical sources of systematic uncertainty.

The observed BDT output distribution compared to the prediction from the simulation of various SM background processes. Dimuon events passing the event selection requirements of the ttH hadronic category, with dimuon mass between 110–150 GeV, are considered. The expected distributions from different production modes of the Higgs boson are also reported. The background uncertainty is obtained by summing in quadrature the uncertainty due to the limited size of the simulated samples with those arising from experimental and theoretical sources of systematic uncertainty.

The observed BDT output distribution compared to the prediction from the simulation of various SM background processes. Dimuon events passing the event selection requirements of the ttH leptonic category, with dimuon mass between 110–150 GeV, are considered. The expected distributions from different production modes of the Higgs boson are also reported. The background uncertainty is obtained by summing in quadrature the uncertainty due to the limited size of the simulated samples with those arising from experimental and theoretical sources of systematic uncertainty.

The observed BDT output distribution compared to the prediction from the simulation of various SM background processes. Dimuon events passing the event selection requirements of the WH leptonic category, with dimuon mass between 110–150 GeV, are considered. The expected distributions from different production modes of the Higgs boson are also reported. The background uncertainty is obtained by summing in quadrature the uncertainty due to the limited size of the simulated samples with those arising from experimental and theoretical sources of systematic uncertainty.

The observed BDT output distribution compared to the prediction from the simulation of various SM background processes. Dimuon events passing the event selection requirements of the ZH leptonic category, with dimuon mass between 110–150 GeV, are considered. The expected distributions from different production modes of the Higgs boson are also reported. The background uncertainty is obtained by summing in quadrature the uncertainty due to the limited size of the simulated samples with those arising from experimental and theoretical sources of systematic uncertainty.

The mass distribution for the weighted combination of VBF-SB and VBF-SR events. Each event is weighted proportionally to the S/(S+B) ratio, calculated as a function of the mass-decorrelated DNN output. The total post-fit and pre-fit uncertainties on the background prediction, along with the best fit signal contribution, are also reported.

The observed profile likelihood ratio as a function of $\kappa_{\mu}$ for a Higgs boson mass of 125.38 GeV, obtained from a combined fit with Ref.[10] in the resolved $\kappa$-framework. The best fit value for $\kappa_{\mu}$ is 1.07 and the corresponding observed 68% CL interval is $0.85 < \kappa_{\mu} < 1.29$.

The best fit estimates for the reduced coupling modifiers extracted for fermions and weak bosons from the resolved $\kappa$-framework compared to their corresponding prediction from the SM. The error bars represent 68% CL intervals for the measured parameters.

Distribution of the transverse momentum of the diphoton system for the $\mathrm{Z}\gamma\gamma$ muon channel. The predicted yields are shown with their pre-fit normalisations. The observed data, the expected signal contribution and the background estimates are presented with error bars showing the corresponding statistical uncertainties.

Distribution of the transverse momentum of the diphoton system, obtained in the control region enriched in misidentified photons for the $\mathrm{W}\gamma\gamma$ electron channel. The predicted yields are shown with their pre-fit normalisations. The observed data, the expected signal contribution and the background estimates are presented with error bars showing the corresponding statistical uncertainties.

Distribution of the transverse momentum of the diphoton system, obtained in the control region enriched in misidentified photons for the $\mathrm{Z}\gamma\gamma$ electron channel. The predicted yields are shown with their pre-fit normalisations. The observed data, the expected signal contribution and the background estimates are presented with error bars showing the corresponding statistical uncertainties.

The measured $\mathrm{W}\gamma\gamma$ and $\mathrm{Z}\gamma\gamma$ fiducial cross sections and corresponding theoretical predictions from MadGraph5_aMC@NLO. The measured yields in the electron and muon channels are extrapolated to a common fiducial phase space determined from simulated signal events at the generated particle level. Generated particles are considered stable if their mean decay length is larger than $1\,\mathrm{cm}$. Generated leptons are required to have a $\mathrm{p}_\mathrm{t}>15\,\mathrm{GeV}$ and $|\eta|<2.5$. The momenta of photons in a cone of $\Delta R=0.1$, the same cone size as the one applied to reconstructed data, are added to the charged lepton momentum to correct for final-state radiation. Generated photons are required to have $\mathrm{p}_\mathrm{t}>15\,\mathrm{GeV}$ and $|\eta|<2.5$. Additionally, the candidate photons are required to have no selected leptons or photons in a cone of radius $\Delta R=0.4$ and no other stable particles, apart from photons and neutrinos, in a cone of radius $\Delta R=0.1$. Events are then selected in the $\mathrm{W}\gamma\gamma$ channel by requiring exactly one electron (muon) with $\mathrm{p}_\mathrm{t}>30\,\mathrm{GeV}$ and at least two photons with $\mathrm{p}_\mathrm{t}>20\,\mathrm{GeV}$. Events are selected in the $\mathrm{Z}\gamma\gamma$ channel by requiring two electrons (muons), at least one of them with $\mathrm{p}_\mathrm{t}>30\,\mathrm{GeV}$, and not less than two photons, each of them with $\mathrm{p}_\mathrm{t}>20\,\mathrm{GeV}$. Additionally, the invariant mass of the dilepton system is required to be $m_{\ell\ell}>55\,\mathrm{GeV}$.

Expected and observed 95% confidence level intervals for the different anomalous couplings parameters in both the $\ell\nu\gamma\gamma$ and $\ell\ell\gamma\gamma$ channels. All parameters are fixed to their SM values except the one that is fitted. No unitarity regularisation scheme is applied. Intervals with a NAN value are not measured in the corresponding channels.