Best fit values of $\sigma(gg\to H\to ZZ)$, $\sigma_i/\sigma_{gg\mathrm{F}}$, and $\mathrm{B}^f/\mathrm{B}^{ZZ}$ relative to their SM prediction from the combined analysis of the $\sqrt{s}$=7 and 8 TeV data. The results are shown for the combination of ATLAS and CMS, and also separately for each experiment, together with their total uncertainties and their breakdown into the four components described in the text. The expected uncertainties in the measurements are also shown.

Best fit values of $\sigma(gg\to H\to WW)$, $\sigma_i/\sigma_{gg\mathrm{F}}$, and $\mathrm{B}^f/\mathrm{B}^{WW}$ from the combined analysis of the $\sqrt{s}$=7 and 8 TeV data. The values involving cross sections are given for $\sqrt{s}$=8 TeV, assuming the SM values for $\sigma_i(7 \mathrm{TeV})/\sigma_i(8 \mathrm{TeV})$. The results are shown for the combination of ATLAS and CMS, and also separately for each experiment, together with their total uncertainties and their breakdown into the four components described in the text. The expected uncertainties in the measurements are also shown.

Best fit values of $\sigma(gg\to H\to WW)$, $\sigma_i/\sigma_{gg\mathrm{F}}$, and $\mathrm{B}^f/\mathrm{B}^{WW}$ relative to their SM prediction from the combined analysis of the $\sqrt{s}$=7 and 8 TeV data. The results are shown for the combination of ATLAS and CMS, and also separately for each experiment, together with their total uncertainties and their breakdown into the four components described in the text. The expected uncertainties in the measurements are also shown.

Best fit values of $\kappa_{gZ}= \kappa_{g}\cdot\kappa_{Z} / \kappa_{H}$ and of the ratios of coupling modifiers, as defined in the most generic parameterisation described in the context of the $\kappa$ framework, from the combined analysis of the $\sqrt{s}$=7 and 8 TeV data. The results are shown for the combination of ATLAS and CMS and also separately for each experiment, together with their total uncertainties and their breakdown into the four components described in the text. The uncertainties in $\lambda_{tg}$ and $\lambda_{WZ}$, for which a negative solution is allowed, are calculated around the overall best fit value. The expected total uncertainties in the measurements are also shown. For those parameters with no sensitivity to the sign, only the absolute values are shown.

Measured global signal strength $\mu$ and its total uncertainty, together with the breakdown of the uncertainty into its four components as defined in the text. The results are shown for the combination of ATLAS and CMS, and separately for each experiment. The expected uncertainty, with its breakdown, is also shown.

Measured signal strengths $\mu$ and their total uncertainties for different Higgs boson production processes. The results are shown for the combination of ATLAS and CMS, and separately for each experiment, for the combined $\sqrt{s}$=7 and 8 TeV data. The expected uncertainties in the measurements are also displayed. These results are obtained assuming that the Higgs boson branching fractions are the same as in the SM.

Measured signal strengths $\mu$ and their total uncertainties for different Higgs boson decay channels. The results are shown for the combination of ATLAS and CMS, and separately for each experiment, for the combined $\sqrt{s}$=7 and 8 TeV data. The expected uncertainties in the measurements are also displayed. These results are obtained assuming that the Higgs boson production process cross sections at $\sqrt{s}$ = 7 and 8 TeV are the same as in the SM.

Results of the ten-parameter fit of $\mu_F^f = \mu_{gg\mathrm{F}+ttH}^f$ and $\mu_V^f = \mu_{\mathrm{VBF}+VH}^f$ for each of the five decay channels. The results are shown for the combination of ATLAS and CMS, together with their measured and expected uncertainties. The measured results are also shown separately for each experiment.

Results of the six-parameter fit of the global ratio $\mu_V/\mu_F = \mu_{\mathrm{VBF}+VH}/\mu_{gg\mathrm{F}+ttH}$ together with $\mu_F^f$ for each of the five decay channels. The results are shown for the combination of ATLAS and CMS, together with their measured and expected uncertainties. The measured results are also shown separately for each experiment.

Fit results allowing BSM loop couplings, assuming that $|\kappa_{V}| \le$ 1, where $\kappa_{V}$ denotes $\kappa_{Z}$ or $\kappa_{W}$, and that $\mathrm{B_{BSM}} \ge$ 0. The results for the combination of ATLAS and CMS are reported with their measured and expected uncertainties. Also shown are the results from each experiment. For the parameters with both signs allowed, the other 1$\sigma$ interval is shown on a second line. For those parameters with no sensitivity to the sign, only the absolute values are shown.

Fit results allowing BSM loop couplings, assuming that there are no additional BSM contributions to the Higgs boson width, i.e. $\mathrm{B_{BSM}}=0$. The results for the combination of ATLAS and CMS are reported with their measured and expected uncertainties. Also shown are the results from each experiment. For the parameters with both signs allowed, the other 1$\sigma$ interval is shown on a second line. When a parameter is constrained and reaches a boundary, namely $\mathrm{B_{BSM}}$ = 0, the uncertainty is not indicated. For those parameters with no sensitivity to the sign, only the absolute values are shown.

Fit results for the parameterisation assuming the absence of BSM particles in the loops ($\mathrm{B_{BSM}}$=0). The results with their measured and expected uncertainties are reported for the combination of ATLAS and CMS, together with the individual results from each experiment. For the parameters with both signs allowed, the other 1$\sigma$ CL interval is shown on a second line. When a parameter is constrained and reaches a boundary, namely $|\kappa_\mu|$ = 0, the uncertainty is not indicated. For those parameters with no sensitivity to the sign, only the absolute values are shown.

Summary of fit results for a parameterisation probing the ratios of coupling modifiers for up-type versus down-type fermions. The results for the combination of ATLAS and CMS are reported together with their measured and expected uncertainties. Also shown are the results from each experiment. The parameter $\kappa_{uu}$ is positive definite since $\kappa_H$ is always assumed to be positive. For the parameter $\lambda_{du}$, for which both signs are allowed, the other 1$\sigma$ CL interval is shown on a second line. Negative values for the parameter $\lambda_{Vu}$ are excluded by more than 4$\sigma$.

Summary of fit results for a parameterisation probing the ratios of coupling modifiers for leptons versus quarks. The results for the combination of ATLAS and CMS are reported together with their measured and expected uncertainties. Also shown are the results from each experiment. The parameter $\kappa_{qq}$ is positive definite since $\kappa_H$ is always assumed to be positive. For the parameter $\lambda_{lq}$, for which there is no sensitivity to the sign, only the absolute values are shown. Negative values for the parameter $\lambda_{Vq}$ are excluded by more than 4$\sigma$.

Correlation matrix obtained from the fit combining the ATLAS and CMS data using the generic parameterisation with 23 parameters, $\sigma_i\cdot\mathrm{B}^f$ for each specific channel $i\to H\to f$. Only 20 parameters are shown because the other three, corresponding to the $H\to ZZ$ decay channel for the $WH$, $ZH$, and $ttH$ production processes, are not measured with a meaningful precision.

Correlation matrix obtained from the fit combining the ATLAS and CMS data using the generic parameterisation with nine parameters, $\sigma(gg\to H\to ZZ)$, $\sigma_i/\sigma_{gg\mathrm{F}}$, and $\mathrm{B}^f/\mathrm{B}^{ZZ}$.

Correlation matrix obtained from the fit combining the ATLAS and CMS data using the generic parameterisation with seven parameters, $\kappa_{gZ}=\kappa_{g}\cdot\kappa_{Z} / \kappa_{H}$ and the ratios of coupling modifiers.

Expected correlation matrix obtained from the fit combining ATLAS and CMS pre-fit Asimov data sets using the generic parameterisation with 23 parameters, $\sigma_i\cdot\mathrm{B}^f$ for each specific channel $i\to H\to f$. Only 20 parameters are shown because the other three, corresponding to the $H\to ZZ$ decay channel for the $WH$, $ZH$, and $ttH$ production processes, are not measured with a meaningful precision.

Expected correlation matrix obtained from the fit combining ATLAS and CMS pre-fit Asimov data sets using the generic parameterisation with nine parameters, $\sigma(gg\to H\to ZZ)$, $\sigma_i/\sigma_{gg\mathrm{F}}$, and $\mathrm{B}^f/\mathrm{B}^{ZZ}$.

Expected correlation matrix obtained from the fit combining ATLAS and CMS pre-fit Asimov data sets using the generic parameterisation with seven parameters, $\kappa_{gZ}=\kappa_{g}\cdot\kappa_{Z} / \kappa_{H}$ and the ratios of coupling modifiers.

ATLAS+CMS combined best-fit and 68% and 95% CL negative log-likelihood contours in the ($\kappa_{g}$,$\kappa_{\gamma}$) plane.

ATLAS best-fit and 68% and 95% CL negative log-likelihood contours in the ($\kappa_{g}$,$\kappa_{\gamma}$) plane.

CMS best-fit and 68% and 95% CL negative log-likelihood contours in the ($\kappa_{g}$,$\kappa_{\gamma}$) plane.

Best-fit and 68% and 95% CL negative log-likelihood contours in the ($\kappa_{F}$,$\kappa_{V}$) plane.

Best-fit and 68% CL negative log-likelihood contour in the ($\kappa_{F}$,$\kappa_{V}$) plane for the $H\rightarrow\gamma\gamma$ channel.

Best-fit and 68% CL negative log-likelihood contour in the ($\kappa_{F}$,$\kappa_{V}$) plane for the $H\rightarrow ZZ$ channel.

Best-fit and 68% CL negative log-likelihood contour in the ($\kappa_{F}$,$\kappa_{V}$) plane for the $H\rightarrow WW$ channel.

Best-fit and 68% CL negative log-likelihood contour in the ($\kappa_{F}$,$\kappa_{V}$) plane for the $H\rightarrow\tau\tau$ channel.

Best-fit and 68% CL negative log-likelihood contour in the ($\kappa_{F}$,$\kappa_{V}$) plane for the $H\rightarrow bb$ channel.

ATLAS best-fit and 68% and 95% CL negative log-likelihood contours in the ($\kappa_{F}$,$\kappa_{V}$) plane.

CMS best-fit and 68% and 95% CL negative log-likelihood contours in the ($\kappa_{F}$,$\kappa_{V}$) plane.

Best-fit and 68% CL negative log-likelihood contour in the ($\kappa_{F}$,$\kappa_{V}$) plane for the $H\rightarrow\gamma\gamma$ channel, assuming $\kappa_{F}>0$.

Best-fit and 68% CL negative log-likelihood contour in the ($\kappa_{F}$,$\kappa_{V}$) plane for the $H\rightarrow ZZ$ channel, assuming $\kappa_{F}>0$.

Best-fit and 68% CL negative log-likelihood contour in the ($\kappa_{F}$,$\kappa_{V}$) plane for the $H\rightarrow WW$ channel, assuming $\kappa_{F}>0$.

Best-fit and 68% CL negative log-likelihood contour in the ($\kappa_{F}$,$\kappa_{V}$) plane for the $H\rightarrow\tau\tau$ channel, assuming $\kappa_{F}>0$.

Best-fit and 68% CL negative log-likelihood contour in the ($\kappa_{F}$,$\kappa_{V}$) plane for the $H\rightarrow bb$ channel, assuming $\kappa_{F}>0$.

Best-fit and 68% and 95% CL negative log-likelihood contours in the ($\mu^{f}_{gg\mathrm{F}+ttH}$,$\mu^{f}_{\mathrm{VBF}+VH}$) plane for the $H\rightarrow\gamma\gamma$ channel.

Best-fit and 68% and 95% CL negative log-likelihood contours in the ($\mu^{f}_{gg\mathrm{F}+ttH}$,$\mu^{f}_{\mathrm{VBF}+VH}$) plane for the $H\rightarrow ZZ$ channel.

Best-fit and 68% and 95% CL negative log-likelihood contours in the ($\mu^{f}_{gg\mathrm{F}+ttH}$,$\mu^{f}_{\mathrm{VBF}+VH}$) plane for the $H\rightarrow WW$ channel.

Best-fit and 68% and 95% CL negative log-likelihood contours in the ($\mu^{f}_{gg\mathrm{F}+ttH}$,$\mu^{f}_{\mathrm{VBF}+VH}$) plane for the $H\rightarrow\tau\tau$ channel.

Best-fit and 68% and 95% CL negative log-likelihood contours in the ($\mu^{f}_{gg\mathrm{F}+ttH}$,$\mu^{f}_{\mathrm{VBF}+VH}$) plane for the $H\rightarrow bb$ channel.

The expected and observed upper limits at 95\% CL on the production cross-section times branching ratio to two photons of the lightest KK graviton as a function of its mass for k/Mpl=0.10. For masses greater than 1000 GeV, pseudo-experiments are used to verify the expected and observed limits, and used in place of the asymptotic limit when differences are observed.

Expected and observed limits computed using asymptotic formulas as a function of the signal mass m_{X} and the relative width $\Gamma_{X}/m_{X}$ for the spin-0 resonance search.

Expected and observed limits computed using asymptotic formulas as a function of the signal mass m_{X} and the relative width $\Gamma_{X}/m_{X}$ for the spin-0 resonance search.

Expected and observed limits computed using asymptotic formulas as a function of the signal mass m_{G*} and k/Mpl for the spin-2 resonance search.

Expected and observed limits computed using asymptotic formulas as a function of the signal mass m_{G*} and k/Mpl for the spin-2 resonance search.

The expected and observed upper limits at 95\% CL on the fiducial cross-section times branching ratio to two photons of a spin-0 resonance as a function of its mass m_X. The decay width of the resonance is equal to 2% of m_X. For masses greater than 1000 GeV, pseudo-experiments are used to verify the expected and observed limits, and used in place of the asymptotic limit when differences are observed.

The expected and observed upper limits at 95\% CL on the fiducial cross-section times branching ratio to two photons of a spin-0 resonance as a function of its mass m_X. The decay width of the resonance is equal to 2% of m_X. For masses greater than 1000 GeV, pseudo-experiments are used to verify the expected and observed limits, and used in place of the asymptotic limit when differences are observed.

The expected and observed upper limits at 95\% CL on the fiducial cross-section times branching ratio to two photons of a spin-0 resonance as a function of its mass m_X. The decay width of the resonance is equal to 6% of m_X. For masses greater than 1000 GeV, pseudo-experiments are used to verify the expected and observed limits, and used in place of the asymptotic limit when differences are observed.

The expected and observed upper limits at 95\% CL on the fiducial cross-section times branching ratio to two photons of a spin-0 resonance as a function of its mass m_X. The decay width of the resonance is equal to 6% of m_X. For masses greater than 1000 GeV, pseudo-experiments are used to verify the expected and observed limits, and used in place of the asymptotic limit when differences are observed.

The expected and observed upper limits at 95\% CL on the fiducial cross-section times branching ratio to two photons of a spin-0 resonance as a function of its mass m_X. The decay width of the resonance is equal to 10% of m_X. For masses greater than 1000 GeV, pseudo-experiments are used to verify the expected and observed limits, and used in place of the asymptotic limit when differences are observed.

The expected and observed upper limits at 95\% CL on the fiducial cross-section times branching ratio to two photons of a spin-0 resonance as a function of its mass m_X. The decay width of the resonance is equal to 10% of m_X. For masses greater than 1000 GeV, pseudo-experiments are used to verify the expected and observed limits, and used in place of the asymptotic limit when differences are observed.

The expected and observed upper limits at 95\% CL on the production cross-section times branching ratio to two photons of the lightest KK graviton as a function of its mass for k/Mpl=0.01. For masses greater than 1000 GeV, pseudo-experiments are used to verify the expected and observed limits, and used in place of the asymptotic limit when differences are observed.

The expected and observed upper limits at 95\% CL on the production cross-section times branching ratio to two photons of the lightest KK graviton as a function of its mass for k/Mpl=0.01. For masses greater than 1000 GeV, pseudo-experiments are used to verify the expected and observed limits, and used in place of the asymptotic limit when differences are observed.

The expected and observed upper limits at 95\% CL on the production cross-section times branching ratio to two photons of the lightest KK graviton as a function of its mass for k/Mpl=0.05. For masses greater than 1000 GeV, pseudo-experiments are used to verify the expected and observed limits, and used in place of the asymptotic limit when differences are observed.

The expected and observed upper limits at 95\% CL on the production cross-section times branching ratio to two photons of the lightest KK graviton as a function of its mass for k/Mpl=0.05. For masses greater than 1000 GeV, pseudo-experiments are used to verify the expected and observed limits, and used in place of the asymptotic limit when differences are observed.

Diphoton invariant mass in the signal region using a binning based on the invariant mass resolution.

Diphoton invariant mass in the signal region using a binning based on the invariant mass resolution.

Diphoton invariant mass in the signal region binned in 1 GeV bins..

Diphoton invariant mass in the signal region binned in 1 GeV bins..

Effect of event selections on a NWA spin-0 and a spin-2 ($k/Mpl=0.01$) MC samples generated for m_X = 1 TeV and on the data. For the MC samples, the event yield is shown after applying event weights; for the data, the absolute yields are shown. The initial yields for data include a trigger preselection that is the OR of a list of single photon and diphoton triggers. The "2 $loose$ photons" step includes the kinematic acceptance cuts. The cross-section times branching ratio for the spin-0 sample is normalized to the MadGraph NLO prediction of 112~fb and the cross-section times branching ratio for the graviton corresponds to the Pythia8 prediction of 12~fb.

Effect of event selections on a NWA spin-0 and a spin-2 ($k/Mpl=0.01$) MC samples generated for m_X = 1 TeV and on the data. For the MC samples, the event yield is shown after applying event weights; for the data, the absolute yields are shown. The initial yields for data include a trigger preselection that is the OR of a list of single photon and diphoton triggers. The "2 $loose$ photons" step includes the kinematic acceptance cuts. The cross-section times branching ratio for the spin-0 sample is normalized to the MadGraph NLO prediction of 112~fb and the cross-section times branching ratio for the graviton corresponds to the Pythia8 prediction of 12~fb.

Parameterization for the CX and AX factors as function of mX for a scalar particle, as obtained from the narrow-width MC signal samples.

Parameterization for the CX and AX factors as function of mX for a scalar particle, as obtained from the narrow-width MC signal samples.

Parameterization for the $C_{X}\cdot A_{X}$ factor as function of $m_{G^*}$ for the RS graviton, as obtained from k/MPl = 0.05 MC signal samples.

Parameterization for the $C_{X}\cdot A_{X}$ factor as function of $m_{G^*}$ for the RS graviton, as obtained from k/MPl = 0.05 MC signal samples.

Parameterization of the Double Sided Crystal Ball function parameters describing the scalar mass resolution model as a function of m_{X} [GeV].

Parameterization of the Double Sided Crystal Ball function parameters describing the scalar mass resolution model as a function of m_{X} [GeV].

Parameterization of the Double Sided Crystal Ball function parameters describing the graviton mass resolution model as a function of m_{X} [GeV].

Parameterization of the Double Sided Crystal Ball function parameters describing the graviton mass resolution model as a function of m_{X} [GeV].

Summary of the observed and expected 95% CL upper limits on the cross section for $\mathrm{LQ}_{3}^{\mathrm{d}}$ pair production as a function of $m_{\mathrm{LQ}_{3}^{\mathrm{d}}}$ under the assumptions of B=1.

Summary of the observed and expected 95% CL upper limits on B as a function of $m_{\mathrm{LQ}_{3}^{\mathrm{d}}}$.

Summary of the observed and expected 95% CL upper limits on B as a function of $m_{\mathrm{LQ}_{3}^{\mathrm{d}}}$.

Cutflow of the preselection requirements (see Section 5) for $\mathrm{LQ}_{3}^{\mathrm{d}}$ signals with $m_{\mathrm{LQ}_{3}^{\mathrm{d}}}$=0.9, 1.1, and 1.3 TeV, assuming B=1. The yields correspond to an integrated luminosity of 139 fb$^{-1}$.

Cutflow of the preselection requirements (see Section 5) for $\mathrm{LQ}_{3}^{\mathrm{d}}$ signals with $m_{\mathrm{LQ}_{3}^{\mathrm{d}}}$=0.9, 1.1, and 1.3 TeV, assuming B=1. The yields correspond to an integrated luminosity of 139 fb$^{-1}$.

Cutflow of the signal region requirements in 1$\ell+\geq 1\tau$ channel (see Table 3) for $\mathrm{LQ}_{3}^{\mathrm{d}}$ signals with $m_{\mathrm{LQ}_{3}^{\mathrm{d}}}$=0.9, 1.1, and 1.3 TeV, assuming B=1. Events that satisfy the preselection requirements are considered. The yields correspond to an integrated luminosity of 139 fb$^{-1}$.

Cutflow of the signal region requirements in 1$\ell+\geq 1\tau$ channel (see Table 3) for $\mathrm{LQ}_{3}^{\mathrm{d}}$ signals with $m_{\mathrm{LQ}_{3}^{\mathrm{d}}}$=0.9, 1.1, and 1.3 TeV, assuming B=1. Events that satisfy the preselection requirements are considered. The yields correspond to an integrated luminosity of 139 fb$^{-1}$.

Cutflow of the signal region requirements in the 2$\ell$OS+$\geq 1\tau$ channel (see Table 4) for $\mathrm{LQ}_{3}^{\mathrm{d}}$ signals with $m_{\mathrm{LQ}_{3}^{\mathrm{d}}}$=0.9, 1.1, and 1.3 TeV, assuming B=1. Events that satisfy the preselection requirements are considered. The yields correspond to an integrated luminosity of 139 fb$^{-1}$.

Cutflow of the signal region requirements in the 2$\ell$OS+$\geq 1\tau$ channel (see Table 4) for $\mathrm{LQ}_{3}^{\mathrm{d}}$ signals with $m_{\mathrm{LQ}_{3}^{\mathrm{d}}}$=0.9, 1.1, and 1.3 TeV, assuming B=1. Events that satisfy the preselection requirements are considered. The yields correspond to an integrated luminosity of 139 fb$^{-1}$.

Cutflow of the signal region requirements in the $2\ell$SS/$3\ell+\geq 1\tau$ channel (see Table 5) for $\mathrm{LQ}_{3}^{\mathrm{d}}$ signals with $m_{\mathrm{LQ}_{3}^{\mathrm{d}}}$=0.9, 1.1, and 1.3 TeV, assuming B=1. Events that satisfy the preselection requirements are considered. In this channel, two signal regions (SR-L and SR-H) are defined based on $p_{\mathrm{T}, 1}^{\tau}$, with SR-L and SR-H requiring $125< p_{\mathrm{T}, 1}^{\tau} < 225$ GeV and $p_{\mathrm{T}, 1}^{\tau}>225$ GeV, respectively. The yields correspond to an integrated luminosity of 139 fb$^{-1}$.

Cutflow of the signal region requirements in the $2\ell$SS/$3\ell+\geq 1\tau$ channel (see Table 5) for $\mathrm{LQ}_{3}^{\mathrm{d}}$ signals with $m_{\mathrm{LQ}_{3}^{\mathrm{d}}}$=0.9, 1.1, and 1.3 TeV, assuming B=1. Events that satisfy the preselection requirements are considered. In this channel, two signal regions (SR-L and SR-H) are defined based on $p_{\mathrm{T}, 1}^{\tau}$, with SR-L and SR-H requiring $125< p_{\mathrm{T}, 1}^{\tau} < 225$ GeV and $p_{\mathrm{T}, 1}^{\tau}>225$ GeV, respectively. The yields correspond to an integrated luminosity of 139 fb$^{-1}$.

Inclusive differential cross section for two electrons and two muons (Max = 1710~GeV).

Mass of the first dilepton pair in the $Z$ mass region (Max = 88.1~GeV).

Mass of the first dilepton pair in the $H$ mass region (Max = 107~GeV).

Mass of the first dilepton pair in the offshell region (Max = 146~GeV).

Mass of the first dilepton pair in the $ZZ$ mass region (Max = 232~GeV).

Mass of the second dilepton pair in the $Z$ mass region (Max = 44.2~GeV).

Mass of the second dilepton pair in the $H$ mass region (Max = 59.6~GeV).

Mass of the second dilepton pair in the offhell region (Max = 88.3~GeV).

Mass of the second dilepton pair in the $ZZ$ mass region (Max = 1000~GeV).

$p_T$ of the first dilepton pair in the $Z$ mass region (Max = 291~GeV).

$p_T$ of the first dilepton pair in the $H$ mass region (Max = 245~GeV).

$p_T$ of the first dilepton pair in the offshell region (Max = 216~GeV).

$p_T$ of the first dilepton pair in the $ZZ$ mass region (Max = 555~GeV).

$p_T$ of the second dilepton pair in the $Z$ mass region (Max = 117~GeV).

$p_T$ of the second dilepton pair in the $H$ mass region (Max = 124~GeV).

$p_T$ of the second dilepton pair in the offhell region (Max = 192~GeV).

$p_T$ of the second dilepton pair in the $ZZ$ mass region (Max = 555~GeV).

$\cos\theta^*$ between first dilepton pair in the $Z$ mass region.

$\cos\theta^*$ between first dilepton pair in the offshell region.

$\cos\theta^*$ between first dilepton pair in the $H$ mass region.

$\cos\theta^*$ between first dilepton pair in the $ZZ$ mass region.

$\cos\theta^*$ between second dilepton pair in the $Z$ mass region.

$\cos\theta^*$ between second dilepton pair in the $H$ mass region.

$\cos\theta^*$ between second dilepton pair in the offshell mass region.

$\cos\theta^*$ between second dilepton pair in the $ZZ$ mass region.

$\Delta y$ between the dilepton pairs in the $Z$ mass region.

$\Delta y$ between the dilepton pairs in the $H$ mass region.

$\Delta y$ between the dilepton pairs in the offshell region.

$\Delta y$ between the dilepton pairs in the $ZZ$ mass region.

$\Delta \phi$ between the dilepton pairs in the $Z$ mass region.

$\Delta \phi$ between the dilepton pairs in the $H$ mass region.

$\Delta \phi$ between the dilepton pairs in the offshell region.

$\Delta \phi$ between the dilepton pairs in the $ZZ$ mass region.

$\Delta \phi$ between the leading lepton pairs in the $Z$ mass region.

$\Delta \phi$ between the leading lepton pairs in the $H$ mass region.

$\Delta \phi$ between the leding lepton pairs in the offshell region.

$\Delta \phi$ between the leading lepton pairs in the $ZZ$ mass region.

Four lepton mass in first $p_T^{4l}$ slice (Max = 726~GeV)

Four lepton mass in second $p_T^{4l}$ slice (Max = 1079~GeV)

Four lepton mass in third $p_T^{4l}$ slice (Max = 1153~GeV)

Four lepton mass in fourth $p_T^{4l}$ slice (Max = 1320~GeV)

Four lepton mass in fifth $p_T^{4l}$ slice (Max = 1413~GeV)

Four lepton mass in first $y^{4l}$ slice (Max = 1710~GeV)

Four lepton mass in second $y^{4l}$ slice (Max = 1116~GeV)

Four lepton mass in third $y^{4l}$ slice (Max = 1413~GeV)

Four lepton mass in fourth $y^{4l}$ slice (Max = 1044~GeV)

Four lepton mass in fifth $y^{4l}$ slice (Max = 918~GeV)

Covariance matrix for M4l

Covariance matrix for M4lvChannel

Covariance matrix for m12_vs_M4l

Covariance matrix for m34_vs_M4l

Covariance matrix for pt12_vs_M4l

Covariance matrix for pt34_vs_M4l

Covariance matrix for CTS12_vs_M4l

Covariance matrix for CTS34_vs_M4l

Covariance matrix for dYpairs_vs_M4l

Covariance matrix for dPhiPairs_vs_M4l

Covariance matrix for dPhiLeadLep_vs_M4l

Covariance matrix for M4lvPt4lbin

Covariance matrix for M4lvRapiditybin

EFT results (linear)

EFT results (full)

Observed and expected limits on $\tan\beta$ as a function of $m_{H^+}$ in the $M_h^{125}(\tilde{\chi})$ scenario. Limits are shown for $\tan\beta$ values in the range of 0.5-60 due to the availability of the model prediction. The bands surrounding the expected limits show the 68% and 95% confidence intervals. Uncertainties in the predicted $H^+$ cross-sections or branching ratios are not considered.

Observed and expected limits on $\tan\beta$ as a function of $m_{H^+}$ in the $M_h^{125}(\tilde{\tau})$ scenario. Limits are shown for $\tan\beta$ values in the range of 0.5-60 due to the availability of the model prediction. The bands surrounding the expected limits show the 68% and 95% confidence intervals. Uncertainties in the predicted $H^+$ cross-sections or branching ratios are not considered.

Observed and expected limits on $\tan\beta$ as a function of $m_{H^+}$ in the $M_h^{125}$(alignment) scenario. Limits are shown for $\tan\beta$ values in the range of 0.5-60 due to the availability of the model prediction. The bands surrounding the expected limits show the 68% and 95% confidence intervals. Uncertainties in the predicted $H^+$ cross-sections or branching ratios are not considered.

Observed and expected limits on $\tan\beta$ as a function of $m_{H^+}$ in the $M_{h_1}^{125}$(CPV) scenario. Limits are shown for $\tan\beta$ values in the range of 0.5-60 due to the availability of the model prediction. The bands surrounding the expected limits show the 68% and 95% confidence intervals. Uncertainties in the predicted $H^+$ cross-sections or branching ratios are not considered.

Event yields of the SM background processes and the 800 GeV $H^{+}$ sample in the four analysis regions before the fit to the data. Uncertainties include both statistical and systematic uncertainties. The yields of the $H^{+}$ signal sample correspond to a cross-section times branching fraction of 10 pb.

Event acceptance for the different $H^+$ mass signal samples.

The \(m_{bb}\) distributions in the MediumBDT categories, overlaid with contributions from the \(H\gamma jj\) signal as well as the resonant \(Z\gamma jj\) and non-resonant \(b\bar{b} \gamma jj\) background fits. The combined \(\chi^2\) per degree of freedom is \(45.2/45\). The bottom panel in each plot presents the significance of the Higgs boson signal relative to the non-resonant \(b\bar{b} \gamma jj\) background in each bin.

The \(m_{bb}\) distributions in the MediumBDT categories, overlaid with contributions from the \(H\gamma jj\) signal as well as the resonant \(Z\gamma jj\) and non-resonant \(b\bar{b} \gamma jj\) background fits. The combined \(\chi^2\) per degree of freedom is \(45.2/45\). The bottom panel in each plot presents the significance of the Higgs boson signal relative to the non-resonant \(b\bar{b} \gamma jj\) background in each bin.

Comparisons of data and simulated event distributions of the BDT input variables \(\Delta R (b1, \gamma)\) in the two \(m_{bb}\) sidebands after kinematic reweighting of the non-resonant \(b\bar{b}\gamma jj\) background. The data are shown as black points, and the background contributions are stacked in coloured histograms. The Higgs boson signal contribution is scaled up and represented by the dashed red line. The bottom panel in each plot shows the ratio between the data and the SM prediction, where the uncertainty band corresponds to the statistical uncertainty only.

Comparisons of data and simulated event distributions of the BDT input variables \(m_{jj}\) in the two \(m_{bb}\) sidebands after kinematic reweighting of the non-resonant \(b\bar{b}\gamma jj\) background. The data are shown as black points, and the background contributions are stacked in coloured histograms. The Higgs boson signal contribution is scaled up and represented by the dashed red line. The bottom panel in each plot shows the ratio between the data and the SM prediction, where the uncertainty band corresponds to the statistical uncertainty only.

Comparisons of data and simulated event distributions of the BDT input variables \(\Delta R (b2, \gamma)\) in the two \(m_{bb}\) sidebands after kinematic reweighting of the non-resonant \(b\bar{b}\gamma jj\) background. The data are shown as black points, and the background contributions are stacked in coloured histograms. The Higgs boson signal contribution is scaled up and represented by the dashed red line. The bottom panel in each plot shows the ratio between the data and the SM prediction, where the uncertainty band corresponds to the statistical uncertainty only.

Comparisons of data and simulated event distributions of the BDT input variables \(\text{centrality}(\gamma, j1, j2)\) in the two \(m_{bb}\) sidebands after kinematic reweighting of the non-resonant \(b\bar{b}\gamma jj\) background. The data are shown as black points, and the background contributions are stacked in coloured histograms. The Higgs boson signal contribution is scaled up and represented by the dashed red line. The bottom panel in each plot shows the ratio between the data and the SM prediction, where the uncertainty band corresponds to the statistical uncertainty only.

Comparisons of data and simulated event distributions of the BDT input variables \(\Delta \phi(b\bar{b}, jj)\) in the two \(m_{bb}\) sidebands after kinematic reweighting of the non-resonant \(b\bar{b}\gamma jj\) background. The data are shown as black points, and the background contributions are stacked in coloured histograms. The Higgs boson signal contribution is scaled up and represented by the dashed red line. The bottom panel in each plot shows the ratio between the data and the SM prediction, where the uncertainty band corresponds to the statistical uncertainty only.

Comparisons of data and simulated event distributions of the BDT input variables \(p_{\mathrm{T}}^{jj}\) in the two \(m_{bb}\) sidebands after kinematic reweighting of the non-resonant \(b\bar{b}\gamma jj\) background. The data are shown as black points, and the background contributions are stacked in coloured histograms. The Higgs boson signal contribution is scaled up and represented by the dashed red line. The bottom panel in each plot shows the ratio between the data and the SM prediction, where the uncertainty band corresponds to the statistical uncertainty only.

Comparisons of data and simulated event distributions of the BDT input variables \(\cos \theta\) in the two \(m_{bb}\) sidebands after kinematic reweighting of the non-resonant \(b\bar{b}\gamma jj\) background. The data are shown as black points, and the background contributions are stacked in coloured histograms. The Higgs boson signal contribution is scaled up and represented by the dashed red line. The bottom panel in each plot shows the ratio between the data and the SM prediction, where the uncertainty band corresponds to the statistical uncertainty only.

Comparisons of data and simulated event distributions of the BDT input variables \(\Delta R(b1,j1)\) in the two \(m_{bb}\) sidebands after kinematic reweighting of the non-resonant \(b\bar{b}\gamma jj\) background. The data are shown as black points, and the background contributions are stacked in coloured histograms. The Higgs boson signal contribution is scaled up and represented by the dashed red line. The bottom panel in each plot shows the ratio between the data and the SM prediction, where the uncertainty band corresponds to the statistical uncertainty only.

Distribution of the four-lepton invariant mass in the four-lepton final state for the ggF-MVA-low category.

Distribution of the four-lepton invariant mass in the four-lepton final state for the VBF-MVA-enriched category.

Distribution of the transverse mass in the llvv final state for the ggF-enriched eevv category.

Distribution of the transverse mass in the llvv final state for the ggF-enriched mumuvv category.

Distribution of the transverse mass in the llvv final state for the VBF-enriched eevv category.

Distribution of the transverse mass in the llvv final state for the VBF-enriched mumuvv category.

The upper limits at 95% CL on the cross section times branching ratio as a function of the heavy resonance mass for the ggF production mode

The upper limits at 95% CL on the cross section times branching ratio as a function of the heavy resonance mass for the VBF production mode

The upper limits at 95% CL on the cross section for the ggF production mode times branching ratio as a function of the heavy resonance mass for an additional heavy scalar assuming a width of 1% of its mass

The upper limits at 95% CL on the cross section for the ggF production mode times branching ratio as a function of the heavy resonance mass for an additional heavy scalar assuming a width of 5% of its mass

The upper limits at 95% CL on the cross section for the ggF production mode times branching ratio as a function of the heavy resonance mass for an additional heavy scalar assuming a width of 10% of its mass

The upper limits at 95% CL on the cross section for the ggF production mode times branching ratio as a function of the heavy resonance mass for an additional heavy scalar assuming a width of 15% of its mass

The upper limits at 95% CL on the cross section times branching ratio for a KK graviton produced with k/M_{PI} = 1

The +1 $\sigma$ band of the observed exclusion contor for the $Z^{\prime}_{B}$ model in the $m_{\chi}$-$m_{Z^{\prime}_{B}}$ plane.

The -1 $\sigma$ band of the observed exclusion contor for the $Z^{\prime}_{B}$ model in the $m_{\chi}$-$m_{Z^{\prime}_{B}}$ plane.

A comparison of the inferred limits to the constraints from direct detection experiments on the spin-independent DM--nucleon cross section in the context of the $Z'_B$ simplified model with vector couplings. Limits are shown at 90% CL.

The observed exclusion contor for the $Z^{\prime}$-2HDM model in the $m_{A}$-$m_{Z^{\prime}}$ plane.

The expected exclusion contor for the $Z^{\prime}$-2HDM model in the $m_{A}$-$m_{Z^{\prime}}$ plane.

The +1 $\sigma$ band of the observed exclusion contor for the $m_{A}$-$Z^{\prime}$-2HDM model in the $m_{Z^{\prime}}$ plane.

The -1 $\sigma$ band of the observed exclusion contor for the $m_{A}$-$Z^{\prime}$-2HDM model in the $m_{Z^{\prime}}$ plane.

The observed exclusion contor for the 2HDM-a model in the $m_{A}$-$m_{a}$ plane.

The expected exclusion contor for the 2HDM-a model in the $m_{A}$-$m_{a}$ plane.

The +1 $\sigma$ band of the observed exclusion contor for the 2HDM-a model in the $m_{A}$-$m_{a}$ plane.

The -1 $\sigma$ band of the observed exclusion contor for the 2HDM-a model in the $m_{A}$-$m_{a}$ plane.

The observed exclusion contor for the 2HDM-a model in the $tan\beta$-$m_{a}$ plane.

The expected exclusion contor for the 2HDM-a model in the $tan\beta$-$m_{a}$ plane.

The +1 $\sigma$ band of the observed exclusion contor for the 2HDM-a model in the $tan\beta$-$m_{a}$ plane.

The -1 $\sigma$ band of the observed exclusion contor for the 2HDM-a model in the $tan\beta$-$m_{a}$ plane.

The exclusion limits at 95% CL for the 2HDM+a model as a function of $\sin \theta$ for $m_{A,H^{\pm},H}$= 600GeV, $m_a$ = 200GeV, $\tan \beta$ = 1.0$.

The exclusion limits at 95% CL for the 2HDM+a model as a function of $\sin \theta$ for $m_{A,H^{\pm},H}$= 1000GeV, $m_a$ = 350GeV, $\tan \beta$ = 1.0$.

Breakdown of the dominant systematic uncertainties.

Event yields in the range of 120 $<m_{\gamma\gamma}<$ 130 GeV for data, signal models, the SM Higgs boson background and non-resonant background in each analysis category, for an integrated luminosity of $139$fb$^{-1}$.

Detailed background contributions from the SM Higgs boson and continuum background for each cut

Detailed contributions from the signals for each cut.

Acceptance times efficiency for several signals in each category.

Efficiency in the Single-bin SR as a function of the M(Sbottom) vs. M(N2) with the $\Delta$ M(N2,N1) $= 130$ GeV. Keep in mind that the efficiency is given in units of $10^{-2}$.

Acceptance in the Multi-bin SR, $\min_{\Theta} < 0.5$ bin as a function of the M(Sbottom) vs. M(N2) with the $\Delta M$(N2,N1) = 130 GeV. Keep in mind that the acceptance is given in units of $10^{-4}$.

Efficiency in the Multi-bin SR, $\min_{\Theta} < 0.5$ bin as a function of the M(Sbottom) vs. M(N2) with the $\Delta M$(N2,N1) = 130 GeV. Keep in mind that the efficiency is given in units of $10^{-2}$.

Acceptance in the Multi-bin SR, $0.5 < \min_{\Theta} < 1.0$ bin as a function of the M(Sbottom) vs. M(N2) with the $\Delta M$(N2,N1) = 130 GeV. Keep in mind that the acceptance is given in units of $10^{-4}$.

Efficiency in the Multi-bin SR, $0.5 < \min_{\Theta} < 1.0$ bin as a function of the M(Sbottom) vs. M(N2) with the $\Delta M$(N2,N1) = 130 GeV. Keep in mind that the efficiency is given in units of $10^{-2}$.

Acceptance in the Multi-bin SR, $\min_{\Theta} > 1.0$ bin as a function of the M(Sbottom) vs. M(N2) with the $\Delta M$(N2,N1) = 130 GeV. Keep in mind that the acceptance is given in units of $10^{-4}$.

Efficiency in the Multi-bin SR, $\min_{\Theta} > 1.0$ bin as a function of the M(Sbottom) vs. M(N2) with the $\Delta M$(N2,N1) = 130 GeV. Keep in mind that the efficiency is given in units of $10^{-2}$.

Observed upper limits on the signal cross section as a function of the M(Sbottom) vs. M(N2) with the $\Delta M$(N2,N1) = 130 GeV.

Expected upper limits on the signal cross section as a function of the M(Sbottom) vs. M(N2) with the $\Delta M$(N2,N1) = 130 GeV.

Cutflows for the bechmarl signal point M(Sbottom) = 800 GeV, M(N2) = 180 GeV. Weighted event yields are reported starting with the "Preselection" line, normalized to an integrated luminosity of $139$ fb$^{−1}$.

Comparison of the expected and observed event yields in the signal regions. The top-quark and Z(mumu) background contributions are scaled with the normalization factors obtained from the background-only fit. The other contribution includes all the backgrounds not explicitly listed in the legend (V+jets except Z(mumu)+jets, di-/triboson, multijet). The hatched band indicates the total statistical and systematic uncertainties in the SM background. The contributions from three signal models to the signal regions are also displayed, where the masses M(Sbottom) and M(N2) are given in GeV in the legend. The lower panel shows the significance of the deviation of the observed yield from the expected background yield.

Dominant systematic uncertainties in the background prediction for the signal regions after the fit to the control regions. “Other” includes the uncertainties arising from muons, jet-vertex tagging, modeling of pile-up, the $E_{T}^{miss}$ computation, multijet background, and luminosity. The individual uncertainties can be correlated and do not necessarily add up quadratically to the total uncertainty.