The expected upper $68\%$ exclusions for the nominal $\sigma$ cross section in the plane of $m_{\phi}$ and $T_D$, for various $m_S$ values, for the case $m_{A'}=1.0$ GeV ($A' \rightarrow \pi^+\pi^-$ with $\mathcal{BR}=100\%$).

The expected lower $68\%$ exclusions for the nominal $\sigma$ cross section in the plane of $m_{\phi}$ and $T_D$, for various $m_S$ values, for the case $m_{A'}=1.0$ GeV ($A' \rightarrow \pi^+\pi^-$ with $\mathcal{BR}=100\%$).

The $95\%$ CL observed upper limits on the cross section shown as a functions of $T_D$ and the pseudoscalar mass $m_{\phi}$ for $m_S = 300$ GeV, for the case $m_{A'}=1.0$ GeV ($A' \rightarrow \pi^+\pi^-$ with $\mathcal{BR}=100\%$).

The $95\%$ CL observed upper limits on the cross section shown as a functions of $T_D$ and the pseudoscalar mass $m_{\phi}$ for $m_S = 600$ GeV, for the case $m_{A'}=1.0$ GeV ($A' \rightarrow \pi^+\pi^-$ with $\mathcal{BR}=100\%$).

The $95\%$ CL observed upper limits on the cross section shown as a functions of $T_D$ and the pseudoscalar mass $m_{\phi}$ for $m_S = 1000$ GeV, for the case $m_{A'}=1.0$ GeV ($A' \rightarrow \pi^+\pi^-$ with $\mathcal{BR}=100\%$).

The observed exclusions for the nominal $\sigma$ cross section in the plane of $m_{\phi}$ and $T_D$, for various $m_S$ values, for the case $m_{A'} = 0.5$ GeV ($A' \rightarrow e^+e^-, \mu^+\mu^-, \pi^+\pi^-$ with $\mathcal{BR} = 40, 40, 20\%$).

The expected exclusions for the nominal $\sigma$ cross section in the plane of $m_{\phi}$ and $T_D$, for various $m_S$ values, for the case $m_{A'} = 0.5$ GeV ($A' \rightarrow e^+e^-, \mu^+\mu^-, \pi^+\pi^-$ with $\mathcal{BR} = 40, 40, 20\%$).

The expected upper $68\%$ exclusions for the nominal $\sigma$ cross section in the plane of $m_{\phi}$ and $T_D$, for various $m_S$ values, for the case $m_{A'} = 0.5$ GeV ($A' \rightarrow e^+e^-, \mu^+\mu^-, \pi^+\pi^-$ with $\mathcal{BR} = 40, 40, 20\%$).

The expected lower $68\%$ exclusions for the nominal $\sigma$ cross section in the plane of $m_{\phi}$ and $T_D$, for various $m_S$ values, for the case $m_{A'} = 0.5$ GeV ($A' \rightarrow e^+e^-, \mu^+\mu^-, \pi^+\pi^-$ with $\mathcal{BR} = 40, 40, 20\%$).

The observed exclusions for the nominal $\sigma$ cross section in the plane of $m_{\phi}$ and $T_D$, for various $m_S$ values, for the case $m_{A'} = 0.7$ GeV ($A' \rightarrow e^+e^-, \mu^+\mu^-, \pi^+\pi^-$ with $\mathcal{BR} = 15, 15, 70\%$).

The expected exclusions for the nominal $\sigma$ cross section in the plane of $m_{\phi}$ and $T_D$, for various $m_S$ values, for the case $m_{A'} = 0.7$ GeV ($A' \rightarrow e^+e^-, \mu^+\mu^-, \pi^+\pi^-$ with $\mathcal{BR} = 15, 15, 70\%$).

The expected upper $68\%$ exclusions for the nominal $\sigma$ cross section in the plane of $m_{\phi}$ and $T_D$, for various $m_S$ values, for the case $m_{A'} = 0.7$ GeV ($A' \rightarrow e^+e^-, \mu^+\mu^-, \pi^+\pi^-$ with $\mathcal{BR} = 15, 15, 70\%$).

The expected lower $68\%$ exclusions for the nominal $\sigma$ cross section in the plane of $m_{\phi}$ and $T_D$, for various $m_S$ values, for the case $m_{A'} = 0.7$ GeV ($A' \rightarrow e^+e^-, \mu^+\mu^-, \pi^+\pi^-$ with $\mathcal{BR} = 15, 15, 70\%$).

For a selection of signal samples in era 2018 with $m_{\phi} = T_D = 3$ GeV, the relative efficiencies and the total efficiency with respect to the selection on the generator-level $H_T$,$H^{gen}_T > 1000$ GeV.

The $95\%$ CL cross section limit for all values of $m_{A'}$, $T_D$, $m_{\phi}$, and $m_S$ used in the scan.

Slices of 2D distributions of observed events and the post-fit templates in the TT pass region, projected onto the plane of leading jet mass mJ1, including expected radion signal at 1.5 TeV.

Slices of 2D distributions of observed events and the post-fit templates in the TT pass region, projected onto the plane of leading jet mass mJ1, including expected radion signal at 1.5 TeV.

Slices of 2D distributions of observed events and the post-fit templates in the TT pass region, projected onto the plane of leading jet mass mJ1, including expected radion signal at 1.5 TeV.

Slices of 2D distributions of observed events and the post-fit templates in the LL pass region, projected onto the plane of corrected HH mass (mHH) mHH, including expected radion signal at 1.5 TeV.

Slices of 2D distributions of observed events and the post-fit templates in the LL pass region, projected onto the plane of corrected HH mass (mHH) mHH, including expected radion signal at 1.5 TeV.

Slices of 2D distributions of observed events and the post-fit templates in the LL pass region, projected onto the plane of corrected HH mass (mHH) mHH, including expected radion signal at 1.5 TeV.

Slices of 2D distributions of observed events and the post-fit templates in the TT pass region, projected onto the plane of corrected HH mass (mHH) mHH, including expected radion signal at 1.5 TeV.

Slices of 2D distributions of observed events and the post-fit templates in the TT pass region, projected onto the plane of corrected HH mass (mHH) mHH, including expected radion signal at 1.5 TeV.

Slices of 2D distributions of observed events and the post-fit templates in the TT pass region, projected onto the plane of corrected HH mass (mHH) mHH, including expected radion signal at 1.5 TeV.

Slices of 2D distributions of observed events and the post-fit templates in the semi-resolved pass region, projected onto the plane of corrected HH mass (mHH) mHH, including expected radion signal at 1.5 TeV.

Slices of 2D distributions of observed events and the post-fit templates in the semi-resolved pass region, projected onto the plane of corrected HH mass (mHH) mHH, including expected radion signal at 1.5 TeV.

Slices of 2D distributions of observed events and the post-fit templates in the semi-resolved pass region, projected onto the plane of corrected HH mass (mHH) mHH, including expected radion signal at 1.5 TeV.

Slices of 2D distributions of observed events and the post-fit templates in the semi-resolved pass region, projected onto the plane of leading jet mass mJ1, including expected radion signal at 1.5 TeV.

Slices of 2D distributions of observed events and the post-fit templates in the semi-resolved pass region, projected onto the plane of leading jet mass mJ1, including expected radion signal at 1.5 TeV.

Slices of 2D distributions of observed events and the post-fit templates in the semi-resolved pass region, projected onto the plane of leading jet mass mJ1, including expected radion signal at 1.5 TeV.

The observed (solid black line) and expected (dashed black line) upper limits at 95% CL on $\sigma$(pp->X)B(X->HH->4b) for the narrow spin-0 radion model. The green (yellow) bands represent one (two) standard deviations from the expected limit. The predicted theoretical cross sections are also shown.

The observed (solid black line) and expected (dashed black line) upper limits at 95% CL on $\sigma$(pp->X)B(X->HH->4b) for the narrow width spin-2 bulk graviton model. The green (yellow) bands represent one (two) standard deviations from the expected limit. The predicted theoretical cross sections are also shown.

The observed (solid black line) and expected (dashed black line) upper limits at 95% CL on $\sigma$(pp->X)B(X->HH->4b) for the 10%-width spin-0 radion model. The green (yellow) bands represent one (two) standard deviations from the expected limit. The predicted theoretical cross sections are also shown.

The observed (solid black line) and expected (dashed black line) upper limits at 95% CL on $\sigma$(pp->X)B(X->HH->4b) for the 10%-width spin-2 bulk graviton model. The green (yellow) bands represent one (two) standard deviations from the expected limit. The predicted theoretical cross sections are also shown.

Correlations of the ttH signal-strength modifiers in the different Higgs pT bins of the STXS measurement.

Observed and expected upper 95% CL limit on the tH signal strength modifier $\mu_{tH}$ with $\mu_{ttH}=1$ for different channels and years, where the uncertainties are uncorrelated between the channels and years, and in their combination. The one and two standard deviation confidence intervals on the expected limit are also listed (sigma).

Likelihood-ratio test statistic as a function of the ttH and tH signal strength modifiers $\mu_{ttH}$ and $\mu_{tH}$. The observed best fit point is $(\mu_{tH}, \mu_{ttH}) = (-3.83, +0.35)$.

Observed and expected likelihood ratio test statistic as a function of $\kappa_{t}$ and $\kappa_{V}$. The observed best fit point is $(\kappa_{t}, \kappa_{V}) = (+0.59, +1.40)$.

Observed and expected likelihood ratio test statistic as a function of $\kappa_{t}$ for $\kappa_{V}=1$. The observed best fit point is $\kappa_{t} = +0.54$.

Observed and expected likelihood ratio test statistic as a function of $\kappa_{t}$ and $\tilde{\kappa}_{t}$ with $\kappa_{V}=1$. The observed best fit point is $(\kappa_{t}, \tilde{\kappa}_{t}) = (+0.53, 0.00)$.

Observed and expected likelihood ratio test statistic as a function of $f_{CP}$. The overall modifier of the top-Higgs coupling strength is profiled such that the likelihood attains its minimum at each point in the plane, and $\kappa_{V}$ is set to 1. The observed best fit point is $f_{CP} = 0.00$.

Observed and expected likelihood ratio test statistic as a function of $\cos(\alpha)$. The overall modifier of the top-Higgs coupling strength is profiled such that the likelihood attains its minimum at each point in the plane, and $\kappa_{V}$ is set to 1. The observed best fit point is $\cos(\alpha) = +1.00$.

Observed and expected likelihood ratio test statistic as a function of $\kappa_{t}$ and $\tilde{\kappa}_{t}$ with $\kappa_{V}=1$, in the $H \rightarrow bb$, $H\rightarrow\gamma\gamma$, $H\rightarrow ZZ$, $H\rightarrow WW$, and $H\rightarrow\tau\tau$ decay channels and in their combination.

Observed and expected likelihood ratio test statistic as a function of $f_{CP}$, in the $H \rightarrow bb$, $H\rightarrow\gamma\gamma$, $H\rightarrow ZZ$, $H\rightarrow WW$, and $H\rightarrow\tau\tau$ decay channel and in their combination. The overall modifier of the top-Higgs coupling strength is profiled such that the likelihood attains its minimum at each point in the plane, and $\kappa_{V}$ is set to 1.

Observed and expected likelihood ratio test statistic as a function of $\cos(\alpha)$, in the $H \rightarrow bb$, $H\rightarrow\gamma\gamma$, $H\rightarrow ZZ$, $H\rightarrow WW$, and $H\rightarrow\tau\tau$ decay channels and in their combination. The overall modifier of the top-Higgs coupling strength is profiled such that the likelihood attains its minimum at each point in the plane, and $\kappa_{V}$ is set to 1.

non prompt $\psi(2S)$ $\lambda_\theta$

The product of detector acceptance and analysis selection efficiency in the muon channel as functions of the particle X mass. Three analysis requirements are applied consecutively: event reconstruction, HLT, and final signal selection. The product of detector acceptance and analysis selection efficiencies are shown at each stage in red, blue, and orange, respectively.

Background-only fit to data. The lower panel contains the pull distribution, defined as the difference between the data yield and the background prediction divided by their combined uncertainty (electron channel).

Background-only fit to data. The lower panel contains the pull distribution, defined as the difference between the data yield and the background prediction divided by their combined uncertainty (muon channel).

The systematic uncertainties affecting signal normalization (electron channel).

The systematic uncertainties affecting signal normalization (muon channel).

The 95% CL expected and observed limits on P P --> X --> W gamma as a function of the resonant mass with leptonic decays of W bosons (narrow width).

The 95% CL expected and observed limits on P P --> X --> W gamma as a function of the resonant mass with leptonic decays of W bosons (broad width).

The 95% CL expected and observed limits on P P --> X --> W gamma as a function of the resonant mass with both leptonic and hadronic decays of W bosons (narrow width).

The 95% CL expected and observed limits on P P --> X --> W gamma as a function of the resonant mass with both leptonic and hadronic decays of W bosons (broad width).

The observed local p-values for narrow resonance width hypotheses.

The observed local p-values for broad resonance width hypotheses.

The 95% CL limits on $|V_{\mu N}|^2$ as a function of the HNL mass for a Majorana HNL.

The 95% CL limits on $|V_{\tau N}|^2$ as a function of the HNL mass for a Majorana HNL.

The 95% CL limits on $|V_{\tau N}|^2$ as a function of the HNL mass for a Majorana HNL.

The 95% CL limits on $|V_{eN}|^2$ as a function of the HNL mass for a Dirac HNL.

The 95% CL limits on $|V_{eN}|^2$ as a function of the HNL mass for a Dirac HNL.

The 95% CL limits on $|V_{\mu N}|^2$ as a function of the HNL mass for a Dirac HNL.

The 95% CL limits on $|V_{\mu N}|^2$ as a function of the HNL mass for a Dirac HNL.

The 95% CL limits on $|V_{\tau N}|^2$ as a function of the HNL mass for a Dirac HNL.

The 95% CL limits on $|V_{\tau N}|^2$ as a function of the HNL mass for a Dirac HNL.

Comparison of the number of predicted background events and the observed events in search regions La. The predicted background yields are shown with the values of the normalizations and nuisance parameters obtained in the background-only fit applied.

Comparison of the number of predicted background events and the observed events in search regions La. The predicted background yields are shown with the values of the normalizations and nuisance parameters obtained in the background-only fit applied.

Signal yields predicted in search regions La as function of HNL mass and coupling strength for a scenario with exclusive coupling between a Majorana HNL and an electron neutrino.

Signal yields predicted in search regions La as function of HNL mass and coupling strength for a scenario with exclusive coupling between a Majorana HNL and an electron neutrino.

Signal yields predicted in search regions La as function of HNL mass and coupling strength for a scenario with exclusive coupling between a Majorana HNL and a muon neutrino.

Signal yields predicted in search regions La as function of HNL mass and coupling strength for a scenario with exclusive coupling between a Majorana HNL and a muon neutrino.

Signal yields predicted in search regions La as function of HNL mass and coupling strength for a scenario with exclusive coupling between a Majorana HNL and a tau neutrino.

Signal yields predicted in search regions La as function of HNL mass and coupling strength for a scenario with exclusive coupling between a Majorana HNL and a tau neutrino.

Comparison of the number of predicted background events and the observed events in search regions Lb. The predicted background yields are shown with the values of the normalizations and nuisance parameters obtained in the background-only fit applied.

Comparison of the number of predicted background events and the observed events in search regions Lb. The predicted background yields are shown with the values of the normalizations and nuisance parameters obtained in the background-only fit applied.

Signal yields predicted in search regions Lb as function of HNL mass and coupling strength for a scenario with exclusive coupling between a Majorana HNL and an electron neutrino.

Signal yields predicted in search regions Lb as function of HNL mass and coupling strength for a scenario with exclusive coupling between a Majorana HNL and an electron neutrino.

Signal yields predicted in search regions Lb as function of HNL mass and coupling strength for a scenario with exclusive coupling between a Majorana HNL and a muon neutrino.

Signal yields predicted in search regions Lb as function of HNL mass and coupling strength for a scenario with exclusive coupling between a Majorana HNL and a muon neutrino.

Signal yields predicted in search regions Lb as function of HNL mass and coupling strength for a scenario with exclusive coupling between a Majorana HNL and a tau neutrino.

Signal yields predicted in search regions Lb as function of HNL mass and coupling strength for a scenario with exclusive coupling between a Majorana HNL and a tau neutrino.

Comparison of the number of predicted background events and the observed events in search regions Ha. The predicted background yields are shown with the values of the normalizations and nuisance parameters obtained in the background-only fit applied.

Comparison of the number of predicted background events and the observed events in search regions Ha. The predicted background yields are shown with the values of the normalizations and nuisance parameters obtained in the background-only fit applied.

Signal yields predicted in search regions Ha as function of HNL mass and coupling strength for a scenario with exclusive coupling between a Majorana HNL and an electron neutrino.

Signal yields predicted in search regions Ha as function of HNL mass and coupling strength for a scenario with exclusive coupling between a Majorana HNL and an electron neutrino.

Signal yields predicted in search regions Ha as function of HNL mass and coupling strength for a scenario with exclusive coupling between a Majorana HNL and a muon neutrino.

Signal yields predicted in search regions Ha as function of HNL mass and coupling strength for a scenario with exclusive coupling between a Majorana HNL and a muon neutrino.

Signal yields predicted in search regions Ha as function of HNL mass and coupling strength for a scenario with exclusive coupling between a Majorana HNL and a tau neutrino.

Signal yields predicted in search regions Ha as function of HNL mass and coupling strength for a scenario with exclusive coupling between a Majorana HNL and a tau neutrino.

Comparison of the number of predicted background events and the observed events in search regions Hb. The predicted background yields are shown with the values of the normalizations and nuisance parameters obtained in the background-only fit applied. Search regions Hb2 and Hb3 are merged for final states including a hadronically decayed tau lepton. This is represented by the value of -1 in the Hb2 entry.

Comparison of the number of predicted background events and the observed events in search regions Hb. The predicted background yields are shown with the values of the normalizations and nuisance parameters obtained in the background-only fit applied. Search regions Hb2 and Hb3 are merged for final states including a hadronically decayed tau lepton. This is represented by the value of -1 in the Hb2 entry.

Signal yields predicted in search regions Hb as function of HNL mass and coupling strength for a scenario with exclusive coupling between a Majorana HNL and an electron neutrino. Search regions Hb2 and Hb3 are merged for final states including a hadronically decayed tau lepton. This is represented by the value of -1 in the Hb2 entry.

Signal yields predicted in search regions Hb as function of HNL mass and coupling strength for a scenario with exclusive coupling between a Majorana HNL and an electron neutrino. Search regions Hb2 and Hb3 are merged for final states including a hadronically decayed tau lepton. This is represented by the value of -1 in the Hb2 entry.

Signal yields predicted in search regions Hb as function of HNL mass and coupling strength for a scenario with exclusive coupling between a Majorana HNL and a muon neutrino. Search regions Hb2 and Hb3 are merged for final states including a hadronically decayed tau lepton. This is represented by the value of -1 in the Hb2 entry.

Signal yields predicted in search regions Hb as function of HNL mass and coupling strength for a scenario with exclusive coupling between a Majorana HNL and a muon neutrino. Search regions Hb2 and Hb3 are merged for final states including a hadronically decayed tau lepton. This is represented by the value of -1 in the Hb2 entry.

Signal yields predicted in search regions Hb as function of HNL mass and coupling strength for a scenario with exclusive coupling between a Majorana HNL and a tau neutrino. Search regions Hb2 and Hb3 are merged for final states including a hadronically decayed tau lepton. This is represented by the value of -1 in the Hb2 entry.

Signal yields predicted in search regions Hb as function of HNL mass and coupling strength for a scenario with exclusive coupling between a Majorana HNL and a tau neutrino. Search regions Hb2 and Hb3 are merged for final states including a hadronically decayed tau lepton. This is represented by the value of -1 in the Hb2 entry.

Covariance matrix for the postfit background prediction in signal region bins La.

Covariance matrix for the postfit background prediction in signal region bins La.

Covariance matrix for the postfit background prediction in signal region bins Lb.

Covariance matrix for the postfit background prediction in signal region bins Lb.

Covariance matrix for the postfit background prediction in signal region bins Ha.

Covariance matrix for the postfit background prediction in signal region bins Ha.

Covariance matrix for the postfit background prediction in signal region bins Hb.

Covariance matrix for the postfit background prediction in signal region bins Hb.

Simulated signal yields passing each stage of the low mass region selection as function of HNL mass and signal strength. Shown for exclusive coupling between a Majorana HNL and an electron neutrino.

Simulated signal yields passing each stage of the low mass region selection as function of HNL mass and signal strength. Shown for exclusive coupling between a Majorana HNL and an electron neutrino.

Simulated signal yields passing each stage of the low mass region selection as function of HNL mass and signal strength. Shown for exclusive coupling between a Majorana HNL and a muon neutrino.

Simulated signal yields passing each stage of the low mass region selection as function of HNL mass and signal strength. Shown for exclusive coupling between a Majorana HNL and a muon neutrino.

Simulated signal yields passing each stage of the low mass region selection as function of HNL mass and signal strength. Shown for exclusive coupling between a Majorana HNL and a tau neutrino.

Simulated signal yields passing each stage of the low mass region selection as function of HNL mass and signal strength. Shown for exclusive coupling between a Majorana HNL and a tau neutrino.

Simulated signal yields passing each stage of the high mass region selection as function of HNL mass and signal strength. Shown for exclusive coupling between a Majorana HNL and an electron neutrino.

Simulated signal yields passing each stage of the high mass region selection as function of HNL mass and signal strength. Shown for exclusive coupling between a Majorana HNL and an electron neutrino.

Simulated signal yields passing each stage of the high mass region selection as function of HNL mass and signal strength. Shown for exclusive coupling between a Majorana HNL and a muon neutrino.

Simulated signal yields passing each stage of the high mass region selection as function of HNL mass and signal strength. Shown for exclusive coupling between a Majorana HNL and a muon neutrino.

Simulated signal yields passing each stage of the high mass region selection as function of HNL mass and signal strength. Shown for exclusive coupling between a Majorana HNL and a tau neutrino.

Simulated signal yields passing each stage of the high mass region selection as function of HNL mass and signal strength. Shown for exclusive coupling between a Majorana HNL and a tau neutrino.

Histogram of the Hbb AK8 jet softdrop mass for data, simulated signal and background events, and the background yields estimated from data in region A.

Histogram of the Hbb AK8 jet softdrop mass for data and simulated signal and background events.

Histogram of $p_{T}(H) + p_{T}(W)$ at the LHE level.

The signal yield estimated by simulation, background yield predicted from data, and observed yield in the signal region.

Observed and expected BDT output in the signal region used for the extraction of the SM-like signal strength. The yields and uncertainties of the expected contributions are shown after the fit to the data. The rightmost bin includes the overflow.

Limits on the signal strength of SM VBS WH production

Observed and expected significance of the measurement of SM VBS WH production

95% confidence limits on the signal strength in the SM analysis

A scan of $-2\Delta\ln L$, where $L$ is the likelihood, for the SM analysis plotted as a function of the signal strength $\mu$ for the expected (red) and the observed (black) data.

The data, fitted signal, and background distributions for the VBS dijet mass in the $t\overline{t}$ control region of the SM analysis.

The data, fitted signal, and background distributions for the W+jets BDT score in the W+jets control region of the SM analysis.

Expected and observed exclusion limits at 95% CL in the asymptotic approximation, on the product of the production cross section and branching fraction into two photons for an additional SM-like Higgs boson, for the VBF plus VH processes, from the analysis of the combined data from 2016, 2017, and 2018.

Expected and observed exclusion limits at 95% CL in the asymptotic approximation, on the product of the production cross section and branching fraction into two photons for an additional SM-like Higgs boson, assuming 100% production via the VBF processes, from the analysis of the combined data from 2016, 2017, and 2018.

Expected and observed exclusion limits at 95% CL in the asymptotic approximation, on the product of the production cross section and branching fraction into two photons for an additional SM-like Higgs boson, assuming 100% production via the VH processes, from the analysis of the combined data from 2016, 2017, and 2018.