The observed 95% CL upper limit on $\sigma( p p \rightarrow H) \cdot (H\rightarrow Z \gamma)$

Expected 2-D likelihood scan of the signal strength $\mu_{\tau\tau}$ versus the $CP$-mixing angle $\phi_{\tau}$.

Free-floating parameters in the measurement.

Measured and expected differential cross-section $\text{d}\sigma / \text{d} m_{4l}$ as a function of $m_{4l}$ in bin of 50$< p_{T}^{4l} <$100 GeV

Measured and expected differential cross-section $\text{d}\sigma / \text{d} m_{4l}$ as a function of $m_{4l}$ in bin of 100$< p_{T}^{4l} <$600 GeV

Measured and expected double differential cross-section $\text{d}\sigma / \text{d} m_{4l}$ as a function of $m_{4l}$ in bin of 0.0$< y_{4l} <$0.4

Measured and expected double differential cross-section $\text{d}\sigma / \text{d} m_{4l}$ as a function of $m_{4l}$ in bin of 0.4$< y_{4l} <$0.8

Measured and expected double differential cross-section $\text{d}\sigma / \text{d} m_{4l}$ as a function of $m_{4l}$ in bin of 0.8$< y_{4l} <$1.2

Measured and expected double differential cross-section $\text{d}\sigma / \text{d} m_{4l}$ as a function of $m_{4l}$ in bin of 1.2$< y_{4l} <$2.5

Measured and expected double differential cross-section $\text{d}\sigma / \text{d} m_{4l}$ as a function of $m_{4l}$ for MELA$<$1.4 bin

Measured and expected double differential cross-section $\text{d}\sigma / \text{d} m_{4l}$ as a function of $m_{4l}$ for MELA$>$1.4 bin

Measured and expected double differential cross-section $\text{d}\sigma / \text{d} m_{4l}$ as a function of $m_{4l}$ for $\mu\mu\mu\mu$ events

Measured and expected double differential cross-section $\text{d}\sigma / \text{d} m_{4l}$ as a function of $m_{4l}$ for $eeee$ events

Measured and expected double differential cross-section $\text{d}\sigma / \text{d} m_{4l}$ as a function of $m_{4l}$ for $ee\mu\mu$ events

Systematic covariance matrix for the differential $m_{4l}$ distribution.

Statistical covariance matrix for the differential $m_{4l}$ distribution.

Background covariance matrix for the differential $m_{4l}$ distribution.

Systematic covariance matrix for the differential $m_{4l}$-$p_{T}^{4l}$ distribution.<br><br> Bins labelled 1-9 correspond to the 0$< p_{T}^{4l} < $20 GeV bin with $m_{4l}$ values as listed in Table 2.<br> Bins labelled 10-16 correspond to the 20$< p_{T}^{4l} <$50 GeV bin with $m_{4l}$ values as listed in Table 3.<br> Bins labelled 17-23 correspond to the 50$< p_{T}^{4l} <$100 GeV bin with $m_{4l}$ values as listed in Table 4.<br> Bins labelled 24-30 correspond to the 100$< p_{T}^{4l} <$600 GeV bin with $m_{4l}$ values as listed in Table 5.

Statistical covariance matrix for the differential $m_{4l}$-$p_{T}^{4l}$ distribution. <br><br> Bins labelled 1-9 correspond to the 0$< p_{T}^{4l} < $20 GeV bin with $m_{4l}$ values as listed in Table 2.<br> Bins labelled 10-16 correspond to the 20$< p_{T}^{4l} <$50 GeV bin with $m_{4l}$ values as listed in Table 3.<br> Bins labelled 17-23 correspond to the 50$< p_{T}^{4l} <$100 GeV bin with $m_{4l}$ values as listed in Table 4.<br> Bins labelled 24-30 correspond to the 100$< p_{T}^{4l} <$600 GeV bin with $m_{4l}$ values as listed in Table 5.

Background covariance matrix for the differential $m_{4l}$-$p_{T}^{4l}$ distribution. <br><br> Bins labelled 1-9 correspond to the 0$< p_{T}^{4l} < $20 GeV bin with $m_{4l}$ values as listed in Table 2.<br> Bins labelled 10-16 correspond to the 20$< p_{T}^{4l} <$50 GeV bin with $m_{4l}$ values as listed in Table 3.<br> Bins labelled 17-23 correspond to the 50$< p_{T}^{4l} <$100 GeV bin with $m_{4l}$ values as listed in Table 4.<br> Bins labelled 24-30 correspond to the 100$< p_{T}^{4l} <$600 GeV bin with $m_{4l}$ values as listed in Table 5.

Systematic covariance matrix for the differential $m_{4l}$-$y_{4l}$ distribution. <br><br> Bins labelled 1-9 correspond to the 0.0$< y_{4l} < $0.4 bin with $m_{4l}$ values as listed in Table 6.<br> Bins labelled 10-18 correspond to the 0.4$< y_{4l} <$0.8 bin with $m_{4l}$ values as listed in Table 7.<br> Bins labelled 19-26 correspond to the 0.8$< y_{4l} <$1.2 bin with $m_{4l}$ values as listed in Table 8.<br> Bins labelled 27-34 correspond to the 1.2$< y_{4l} <$2.5 bin with $m_{4l}$ values as listed in Table 9.

Statistical covariance matrix for the differential $m_{4l}$-$y_{4l}$ distribution. <br><br> Bins labelled 1-9 correspond to the 0.0$< y_{4l} < $0.4 bin with $m_{4l}$ values as listed in Table 6.<br> Bins labelled 10-18 correspond to the 0.4$< y_{4l} <$0.8 bin with $m_{4l}$ values as listed in Table 7.<br> Bins labelled 19-26 correspond to the 0.8$< y_{4l} <$1.2 bin with $m_{4l}$ values as listed in Table 8.<br> Bins labelled 27-34 correspond to the 1.2$< y_{4l} <$2.5 bin with $m_{4l}$ values as listed in Table 9.

Background covariance matrix for the differential $m_{4l}$-$y_{4l}$ distribution. <br><br> Bins labelled 1-9 correspond to the 0.0$< y_{4l} < $0.4 bin with $m_{4l}$ values as listed in Table 6.<br> Bins labelled 10-18 correspond to the 0.4$< y_{4l} <$0.8 bin with $m_{4l}$ values as listed in Table 7.<br> Bins labelled 19-26 correspond to the 0.8$< y_{4l} <$1.2 bin with $m_{4l}$ values as listed in Table 8.<br> Bins labelled 27-34 correspond to the 1.2$< y_{4l} <$2.5 bin with $m_{4l}$ values as listed in Table 9.

Systematic covariance matrix for the differential $m_{4l}$-MELA distribution. <br><br> Bins labelled 1-7 correspond to the MELA$<$ 0.4 bin with $m_{4l}$ values as listed in Table 10.<br> Bins labelled 8-12 correspond to the MELA$>$ 0.4 bin with $m_{4l}$ values as listed in Table 11.

Statistical covariance matrix for the differential $m_{4l}$-MELA distribution. <br><br> Bins labelled 1-7 correspond to the MELA$<$ 0.4 bin with $m_{4l}$ values as listed in Table 10.<br> Bins labelled 8-12 correspond to the MELA$>$ 0.4 bin with $m_{4l}$ values as listed in Table 11.

Background covariance matrix for the differential $m_{4l}$-MELA distribution. <br><br> Bins labelled 1-7 correspond to the MELA$<$ 0.4 bin with $m_{4l}$ values as listed in Table 10.<br> Bins labelled 8-12 correspond to the MELA$>$ 0.4 bin with $m_{4l}$ values as listed in Table 11.

Systematic covariance matrix for the differential $m_{4l}$-lepton flavour distribution. <br><br> Bins labelled 1-9 correspond to $\mu\mu\mu\mu$ events with $m_{4l}$ values as listed in Table 12.<br> Bins labelled 10-18 correspond to $eeee$ events with $m_{4l}$ values as listed in Table 13.<br> Bins labelled 19-27 correspond to $\mu\mu\mu\mu$ events with $m_{4l}$ values as listed in Table 14.<br>

Statistical covariance matrix for the differential $m_{4l}$-lepton flavour distribution. <br><br> Bins labelled 1-9 correspond to $\mu\mu\mu\mu$ events with $m_{4l}$ values as listed in Table 12.<br> Bins labelled 10-18 correspond to $eeee$ events with $m_{4l}$ values as listed in Table 13.<br> Bins labelled 19-27 correspond to $\mu\mu\mu\mu$ events with $m_{4l}$ values as listed in Table 14.<br>

Background covariance matrix for the differential $m_{4l}$-lepton flavour distribution. <br><br> Bins labelled 1-9 correspond to $\mu\mu\mu\mu$ events with $m_{4l}$ values as listed in Table 12.<br> Bins labelled 10-18 correspond to $eeee$ events with $m_{4l}$ values as listed in Table 13.<br> Bins labelled 19-27 correspond to $\mu\mu\mu\mu$ events with $m_{4l}$ values as listed in Table 14.<br>

Observed correlations between the measured reduced stage-1.2 simplified template $VH, H \rightarrow b\bar{b}$ cross-sections (STXS), including both the statistical and systematic uncertainties. All leptonic decays of the $V$ bosons (including those to $\tau$ leptons, $\ell = e, \mu, \tau$) are considered.

Linear combination of Wilson coefficients corresponding to SMEFT operators in the Warsaw basis for which this analysis provides orthogonal constraints (eigenvectors). Eigenvalues are shown for each eigenvector, which provides a measure of the experimental sensitivity to that linear combination. The modifications to the $qq\rightarrow ZH$ and $qq\rightarrow WH$ processes due to SMEFT operators are computed at LO, and changes to the $gg\rightarrow ZH$ process are neglected.

Linear combination of Wilson coefficients corresponding to SMEFT operators in the Warsaw basis for which this analysis provides orthogonal constraints (eigenvectors). Eigenvalues are shown for each eigenvector, which provides a measure of the experimental sensitivity to that linear combination. The modifications to the $qq\rightarrow ZH$ and $qq\rightarrow WH$ processes due to SMEFT operators are computed at LO, and changes to the $gg\rightarrow ZH$ process are neglected.

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

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

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

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

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_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 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_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 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_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 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_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 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_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_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_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_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_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 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_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 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_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 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_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_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_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_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_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 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_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 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_0 signal region category of the ll channel.

Observed yields in the VBF_1 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.

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

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.

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.

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.

Measured differential cross section with associated uncertainties as a function of COS(THETA*). Each systematic uncertainty sources 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 DELTAYRAP(2GAMMA). Each systematic uncertainty sources 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 MULT(JET,PT>30 GEV). Each systematic uncertainty sources 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 MULT(JET,PT>50 GEV). Each systematic uncertainty sources 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 PT(JET1). Each systematic uncertainty sources 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 PT(JET2). Each systematic uncertainty sources 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 HT. Each systematic uncertainty sources 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 YRAP(JET1). Each systematic uncertainty sources 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 YRAP(JET2). Each systematic uncertainty sources 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(2JET). Each systematic uncertainty sources 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 DELTAYRAP(2JET). Each systematic uncertainty sources 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 ABSDPHI(2JET). Each systematic uncertainty sources 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 DPHI(2JET). Each systematic uncertainty sources 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 PT(2GAMMA2JET). Each systematic uncertainty sources 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 DPHI(2GAMMA,2JET). Each systematic uncertainty sources 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 TAUJET. Each systematic uncertainty sources 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 SUM(TAUJET). Each systematic uncertainty sources 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 PT(2GAMMA) [NJET=0,PT>30 GEV]. Each systematic uncertainty sources 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 PT(2GAMMA) [NJET=1,PT>30 GEV]. Each systematic uncertainty sources 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 PT(2GAMMA) [NJET=2,PT>30 GEV]. Each systematic uncertainty sources 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 PT(2GAMMA) [NJET>=3,PT>30 GEV]. Each systematic uncertainty sources 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.

The measured cross sections or cross section limits of the diphoton, VBF-enhanced, Nlepton $\geq$ 1, high $E_{T}^{miss}$, and ttH-enhanced fiducial regions are shown.

Measured differential cross section with associated uncertainties as a function of diphoton transverse momentum in bins of ABS(COS(THETA*)). Each systematic uncertainty sources 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.Each systematic uncertainty sources is fully uncorrelated with the other sources.

Non-perturbative correction factors in percent accounting for the impact of hadronisation and the underlying event activity for all measured variables and fiducial regions. Regions of phase space where no reliable estimate could be obtained are listed as 100 without uncertainties. Uncertainties are evaluated by deriving these factors using different generators and tunes as described in the text. No factor are given for the Nlepton $\geq$ 1 and High-$E_{T}^{miss}$ fiducial regions as the gluon fusion contamination in both is negligible.

Isolation efficiencies in percent for gluon fusion $H\rightarrow\gamma\gamma$ for each fiducial region/variable bin measured in this analysis. The isolation efficiency is defined as the probability for both photons to fulfil the isolation criteria (as described in Section 9.1) for events that pass the diphoton kinematic criteria. Regions of phase space where no reliable estimate could be obtained are listed as 100 without uncertainties. Uncertainties are assigned in the same way as for the non-perturbative correction factors: by varying the fragmentation and underlying event modelling. These factors can be multiplied by the kinematic acceptance factors (see Table 29) to extrapolate an inclusive gluon fusion Higgs prediction to the fiducial volume used in this analysis. No factors for the Nlepton $\geq$ 1 and High $E_{T}^{miss}$ fiducial regions are provided as the gluon fusion contamination is negligible.

Combined non-perturbative (Table 27) and particle-level isolation correction factors (Table 28) in percent accounting for the impact of hadronisation and the underlying event activity for all measured variables and fiducial regions. Regions of phase space where no reliable estimate could be obtained are listed as 100 without uncertainties. The uncertainties on the combined values properly take into account the correlations between both multiplicative factors.

Diphoton kinematic acceptances in percent for gluon-gluon fusion for the diphoton fiducial region and all differential variable bins studied in this paper, defined as the probability to fulfill the diphoton kinematic criteria: $p_{T}$/$m_{\gamma\gamma}$ < 0.35 (0.25) for the leading (subleading) photon and $|\eta_{\gamma\gamma}|$ < 2.37. The factors are evaluated using the Powheg NNLOPSevent generator. Uncertainties are taken from PDF variations. QCD scale variations have a negligible impact on these factors. The range of each bin is given in Table 26.

observed statistical correlations between pTyy, Njets, mjj, |DeltaPhijj| and pTj1

ggH default MC + XH predictions

XH ( = VBF + VH + ttH + bbH ) MC predictions

Best-fit values and uncertainties of the production-mode cross sections times branching ratio.

Best-fit values and uncertainties of the simplified template cross sections times branching ratio.

Observed correlations between the measured simplified template cross sections, including both the statistical and systematic uncertainties.

Best-fit values and uncertainties of the simplified template cross sections times branching ratio.

Observed correlations between the measured simplified template cross sections, including both the statistical and systematic uncertainties.

Observed correlations between the measured simplified template cross sections, including both the statistical and systematic uncertainties.

Post-fit $m_\text{T}(\gamma, E_\text{T}^\text{miss})$ distribution in the inclusive signal region for the dark-photon search with the 125 GeV mass $\mathcal{B}(H \to \gamma \gamma_\text{d})$ signal normalization set to zero. A $H \to \gamma \gamma_\text{d}$ decay signal is shown for two different mass hypotheses, 125 GeV and 500 GeV, and scaled to a $\mathcal{B}(H \to \gamma \gamma_\text{d})$ of 2% and 1%, respectively. Events with $m_\text{T}(\gamma, E_\text{T}^\text{miss})$ larger than the rightmost bin boundary are added to that bin.

The 95% CL upper limit on the Higgs boson production cross-section times branching ratio to $\gamma \gamma_\text{d}$ is shown for different VBF-produced scalar-mediator-mass hypotheses in the NWA. The theoretically predicted cross-section of a Higgs boson produced via VBF and with the $\mathcal{B}(H \to \gamma \gamma_\text{d}) =$ 5% is superimposed on the $\pm 1\sigma$ and $\pm 2\sigma$ NNLO QCD + NLO EW uncertainty band of the expected production cross-section limit.

Post-fit $m_\text{jj}$ distribution in the inclusive signal region. The Higgs boson invisible decay signal is scaled to a $\mathcal{B}_\text{inv}$ of 37%. Events with $m_\text{jj}$ larger than the rightmost bin boundary are added to that bin.

Post-fit $m_\text{jj}$ distribution in the one-lepton control region $W_{\ell \nu}^\gamma$ CR. Events with $m_\text{jj}$ larger than the rightmost bin boundary are added to that bin.

Post-fit $m_\text{T}$ distribution in the one lepton control region. Events with $m_\text{T}$ larger than the rightmost bin boundary are added to that bin.

Post-fit photon centrality distribution in the zero lepton signal plus control region with the $\mathcal{B}_\text{inv}$ signal normalization set to zero in the fit.

Post-fit photon $E_\text{T}$ distribution in the zero lepton signal region with the $\mathcal{B}_\text{inv}$ signal normalization set to zero in the fit.

Post-fit photon centrality distribution in the zero lepton signal plus control region resulting from the fit to the $m_\text{jj}$ distribution for EW $Z \gamma + \text{jets}$. The post-fit uncertainties include statistical, experimental, and theory contributions.

Post-fit photon $E_\text{T}$ distribution in the zero lepton signal region resulting from the fit to the $m_\text{jj}$ distribution for EW $Z \gamma + \text{jets}$. The post-fit uncertainties include statistical, experimental, and theory contributions.

Post-fit DNN output score distribution in the one lepton control region.

Yields for the EW $Z \gamma + \text{jets}$ process are shown after each selection along with relative and absolute signal acceptance efficiencies.

Yields for the 125 GeV Higgs boson with $\mathcal{B}_\text{inv.} =$ 1 signal produced by the vector boson fusion process in association with a final state photon are shown after each selection along with relative and absolute signal acceptance efficiencies.

Yields for the 125 GeV Higgs boson with $\mathcal{B}(H \to \gamma \gamma_\text{d}) =$ 1 signal produced by the vector boson fusion process are shown after each selection along with relative and absolute signal acceptance efficiencies.

Distribution of $m^{vis}_{HH}$ after applying all the event selection that define the $HH$ SR, except the requirement on $m^{vis}_{HH}$. The background labelled as "Others" contains $W$+jets, diboson, $t\bar{t}$ and single-top-quark processes. The $X\rightarrow HH \rightarrow b\bar{b}\tau^{+}\tau^{-}$ signal is overlaid for two resonance mass hypotheses with a cross-section set to the expected limit, while all backgrounds are pre-fit. The first and the last bins contains the under-flow and over-flow bin entries, respectively. The hatched bands represent combined statistical and systematic uncertainties.

Event yields of the various estimated backgrounds and data, computed in the signal region of the search for $X\rightarrow HH \rightarrow b\bar{b}\tau^{+}\tau^{-}$. The background labelled as "Others" contains $W$+jets, diboson, $t\bar{t}$ and single-top-quark processes. Statistical and systematic uncertainties are quoted. The background yields and uncertainties are pre-fit and are found to be similar to those post-fit.

Event yields of the various estimated backgrounds and data, computed in the signal region of the search for $X\rightarrow HH \rightarrow b\bar{b}\tau^{+}\tau^{-}$. The background labelled as "Others" contains $W$+jets, diboson, $t\bar{t}$ and single-top-quark processes. Statistical and systematic uncertainties are quoted. The background yields and uncertainties are pre-fit and are found to be similar to those post-fit.

Expected and observed 95% CL upper limits on the production of a heavy, narrow-width, scalar resonance decaying to a pair of Higgs bosons ($X\rightarrow HH$). The final state used in the search consists of a boosted $b\bar{b}$ pair and a boosted hadronically decaying $\tau^{+}\tau^{-}$ pair, and the SM braching ratio of the Higgs boson are assumed. The $\pm 1\sigma$ and $\pm 2\sigma$ variations about the expected limit are indicated by the error bands. Two different requirements are applied on the visible mass of the two boosted Higgs boson candidates for the resonance mass hypotheses of 1.6 TeV and 2.5 TeV, leading to discontinuities in the limits (at 1.6 TeV, the difference between imposing no requirement and $m^{vis}_{HH}>900$ GeV is less than 1% though).

Expected and observed 95% CL upper limits on the production of a heavy, narrow-width, scalar resonance decaying to a pair of Higgs bosons ($X\rightarrow HH$). The final state used in the search consists of a boosted $b\bar{b}$ pair and a boosted hadronically decaying $\tau^{+}\tau^{-}$ pair, and the SM braching ratio of the Higgs boson are assumed. The $\pm 1\sigma$ and $\pm 2\sigma$ variations about the expected limit are indicated by the error bands. Two different requirements are applied on the visible mass of the two boosted Higgs boson candidates for the resonance mass hypotheses of 1.6 TeV and 2.5 TeV, leading to discontinuities in the limits (at 1.6 TeV, the difference between imposing no requirement and $m^{vis}_{HH}>900$ GeV is less than 1% though).