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

Measured fiducial cross section in fb as a function of Yll.

Measured fiducial cross section in fb/GeV as a function of the leading jet pT, PT(JET1).

Measured normalised fiducial cross section as a function of Njet. Jet PT>25 GeV for |eta|<2.4 and PT>30 GeV for 2.4<|eta|<4.5.

Measured normalised fiducial cross section in 1/GeV as a function of pTH.

Measured normalised fiducial cross section as a function of Yll.

Measured normalised fiducial cross section in 1/GeV as a function of the leading jet pT, PT(JET1).

Jet veto efficiency as a function of PT(JET1) pT threshold.

Full correlation matrix for the measured differential cross section as a function of Njet, taking into account statistical and systematic uncertainties.

Full correlation matrix for the measured differential cross section as a function of pTH, taking into account statistical and systematic uncertainties.

Full correlation matrix for the measured differential cross section as a function of Yll, taking into account statistical and systematic uncertainties.

Full correlation matrix for the measured differential cross section as a function of the leading jet pT, taking into account statistical and systematic uncertainties.

Full correlation matrix for the measured normalised differential cross section as a function of Njet, taking into account statistical and systematic uncertainties.

Full correlation matrix for the measured normalised differential cross section as a function of pTH, taking into account statistical and systematic uncertainties.

Full correlation matrix for the measured normalised differential cross section as a function of Yll, taking into account statistical and systematic uncertainties.

Full correlation matrix for the measured normalised differential cross section as a function of the leading jet pT, taking into account statistical and systematic uncertainties.

Fiducial cross section of H+0-jet events in femtobarn for different values of PT(JET1) pT threshold with the uncertainties for each bin given in absolute values and in percent. The asterisk for the 25 GeV column header indicates that the results are for a mixed jet pT threshold, which is raised from 25 GeV to 30 GeV for jets with 2.5 < |eta| < 4.5.

Correction factors from inclusive parton level to fiducial particle level for bins of Njet derived with POWHEG NNLOPS+Pythia8. Uncertainties due to the non-perturbative correction are dominant for the correction factors.

Correction factors from inclusive parton level to fiducial particle level for bins of the leading jet pT derived with POWHEG NNLOPS+Pythia8. Uncertainties due to the non-perturbative correction are dominant for the correction factors.

Correction factors from inclusive parton level to fiducial particle level for the jet-veto efficiency with different jet pT thresholds derived with POWHEG NNLOPS+Pythia8. The asterisk on the 25 GeV bin label indicates that the results are for a mixed pT threshold, which is raised from 25 GeV to 30 GeV for jets with 2.5 < |eta| < 4.5, corresponding to the selection used to define the signal regions for the analysis.

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.

Measured differential cross section with associated uncertainties as a function of absolute rapidity of diphoton system. 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 multiplicity of jets with transverse momentum pT(jet) > 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 multiplicity of jets with transverse momentum pT(jet) > 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 multiplicity of jets with transverse momentum pT(jet) > 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 multiplicity of jets with transverse momentum pT(jet) > 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 transverse momentum of the leading jet. 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 transverse momentum of the leading jet. 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 absolute rapidity of the leading jet. 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 absolute rapidity of the leading jet. 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 scalar transverse momentum sum of all jets. 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 scalar transverse momentum sum of all jets. 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 transverse momentum of second leading jet. 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 transverse momentum of second leading jet. 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 rapidity separation of leading two jets. 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 rapidity separation of leading two jets. 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 azimuthal angle between diphoton and dijet systems. 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 azimuthal angle between diphoton and dijet systems. 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 cosine of the decay angle in the Collins-Soper frame. 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 cosine of the decay angle in the Collins-Soper frame. 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 cosine of the decay angle in the Collins-Soper frame. 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 cosine of the decay angle in the Collins-Soper frame. 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 azimuthal angle between the two leading jets. 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 azimuthal angle between the two leading jets. 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 diphoton thrust pT. 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 diphoton thrust pT. 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 rapidity separation of the two photons. 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 rapidity separation of the two photons. 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 tau of highest-tau jet (see paper for description). 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 tau of highest-tau jet (see paper for description). 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 scalar sum of tau for all jets (see paper for description). 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 scalar sum of tau for all jets (see paper for description). 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 absolute rapidity of second leading jet. 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 absolute rapidity of second leading jet. 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 transverse momentum of third leading jet. 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 transverse momentum of third leading jet. 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 dijet invariant mass. 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 dijet invariant mass. 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 transverse momentum of the combined diphoton and dijet system. 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 transverse momentum of the combined diphoton and dijet system. 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 cosine of the decay angle in the Collins-Soper frame in bins of diphoton transverse momentum. 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 cosine of the decay angle in the Collins-Soper frame in bins of diphoton transverse momentum. 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 diphoton transverse momentum in jet multiplicity bins. 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 diphoton transverse momentum in jet multiplicity bins. 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 leading jet transverse momentum for exclusive one-jet events. 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 leading jet transverse momentum for exclusive one-jet events. 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.

Diphoton kinematic acceptances in percent for gluon fusion for each fiducial region/variable bin studied in this paper, defined as the probability to fulfil the diphoton kinematic criteria: pT(2gam) < 0.35 (0.25) for the leading (subleading) photon and |eta(gam)|<2.37. The factors are evaluated using the POWHEG event generator with MPI modelling and hadronisation turned off. Consistent results for the diphoton variables are obtained by HRes 2.2. Uncertainties are taken from PDF variations. QCD scale varaitions have a negligible impact on these factors.

Diphoton kinematic acceptances in percent for gluon fusion for each fiducial region/variable bin studied in this paper, defined as the probability to fulfil the diphoton kinematic criteria: pT(2gam) < 0.35 (0.25) for the leading (subleading) photon and |eta(gam)|<2.37. The factors are evaluated using the POWHEG event generator with MPI modelling and hadronisation turned off. Consistent results for the diphoton variables are obtained by HRes 2.2. Uncertainties are taken from PDF variations. QCD scale varaitions have a negligible impact on these factors.

Isolation efficiencies in percent for gluon fusion H -> 2gam 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 (ETiso < 14 GeV as described in the text) for events that pass the diphoton kinematic criteria. 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~ ef{tab:fid_acceptance}) to extrapolate an inclusive gluon fusion Higgs prediction to the fiducial volume used in this analysis.

Isolation efficiencies in percent for gluon fusion H -> 2gam 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 (ETiso < 14 GeV as described in the text) for events that pass the diphoton kinematic criteria. 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~ ef{tab:fid_acceptance}) to extrapolate an inclusive gluon fusion Higgs prediction to the fiducial volume used in this analysis.

Non-perturbative correction factors in percent accounting for the impact of hadronisation and the underlying event activity for all measured variables and fiducial regions. Uncertainties are evaluated by deriving these factors using different generators and tunes as described in the text.

Non-perturbative correction factors in percent accounting for the impact of hadronisation and the underlying event activity for all measured variables and fiducial regions. Uncertainties are evaluated by deriving these factors using different generators and tunes as described in the text.

Fiducial cross sections in fb from WH, ZH and ttH combined, in each variable bin and fiducial region.

Fiducial cross sections in fb from WH, ZH and ttH combined, in each variable bin and fiducial region.

Statistical bin-to-bin correlations between the 24 bins used in the analysis described in Phys.Lett. B753 (2016) 69-85 (arXiv:1508.02507). See Figure 3 of this paper for the bins used.

Measured cross section in bins of $p_{\rm{T}}^{\rm{j}1}$.

Measured fractions in bins of $p_{\rm{T}}^{\rm{H}}$.

Measured fractions in bins of $|y^{\rm{H}}|$.

Measured fractions in bins of $N_{\rm{jets}}$.

Measured fractions in bins of $p_{\rm{T}}^{\rm{j1}}$.

Correlation matrix for the total uncertainty of the differential cross-section measurement in bins of $p_{\rm{T}}^{\rm{H}}$.

Correlation matrix for the total uncertainty of the differential cross-section measurement in bins of $|y^{\rm{H}}|$.

Correlation matrix for the total uncertainty of the differential cross-section measurement in bins of $N_{\rm{jets}}$.

Correlation matrix for the total uncertainty of the differential cross-section measurement in bins of $p_{\rm{T}}^{\rm{j1}}$.

Correlation matrix for the total uncertainty of the differential shape measurement in bins of $p_{\rm{T}}^{\rm{H}}$.

Correlation matrix for the total uncertainty of the differential shape measurement in bins of $|y^{\rm{H}}|$.

Correlation matrix for the total uncertainty of the differential shape measurement in bins of $N_{\rm{jets}}$.

Correlation matrix for the total uncertainty of the differential shape measurement in bins of $p_{\rm{T}}^{\rm{j1}}$.

The total $pp\to H$ cross section in the $H\to\gamma\gamma$ and $H\to ZZ^*\to 4\ell$ channels and their combination.

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).

Summary of the 95% C.L. upper limits on $\sigma_{ZH} BR(H\rightarrow aa \rightarrow (b\bar{b})(b\bar{b}))$ for the dilepton channel in the resolved analysis. The observed limits are shown, together with the expected limits (dotted black lines). In the case of the expected limits, one- and two-standard-deviation uncertainty bands are also displayed. The data was published in JHEP 10 (2018) 031.

Efficiency and acceptance for simulated $ZH(\rightarrow aa\rightarrow (b\bar{b})(b\bar{b}))$ samples in two signal regions (SR) of the analysis, one with two $a\to b\bar{b}$ candidates in the High Purity Category (HPC), and the other with one $a\to b\bar{b}$ candidate in the High Purity Category (HPC) and one in the Low Purity Category (LPC).

Efficiency and acceptance for simulated $ZH(\rightarrow aa\rightarrow (b\bar{b})(b\bar{b}))$ samples in two signal regions (SR) of the analysis, one with two $a\to b\bar{b}$ candidates in the High Purity Category (HPC), and the other with one $a\to b\bar{b}$ candidate in the High Purity Category (HPC) and one in the Low Purity Category (LPC).

Event yields for a simulated $ZH(\rightarrow aa\rightarrow (b\bar{b})(b\bar{b}))$ sample with $m_a = 17.5\,\text{GeV}$. The signal sample is produced with cross section equals to the standard model $pp\to ZH$, i.e. $0.88\,\text{pb}$. Cut 0 corresponds to the initial number of events. Cut 1 requires the single lepton trigger. Cut 2 requires 2 identified leptons. Cut 3 requires the Z-boson mass window. Cut 4 requires 2 reconstructed $a\to b\bar{b}$ candidates. Cut 5a requires 2 identified $a\to b\bar{b}$ candidates in the 1HPC1LPC region. Cut 6a requires the 2 $a\to b\bar{b}$ candidates in the 1HPC1LPC region to be inside the Higgs mass window. Cut 5b requires 2 identified $a\to b\bar{b}$ candidates in the 2HPC region. Cut 6b requires the 2 $a\to b\bar{b}$ candidates in the 2HPC region to be inside the Higgs mass window.

Event yields for a simulated $ZH(\rightarrow aa\rightarrow (b\bar{b})(b\bar{b}))$ sample with $m_a = 17.5\,\text{GeV}$. The signal sample is produced with cross section equals to the standard model $pp\to ZH$, i.e. $0.88\,\text{pb}$. Cut 0 corresponds to the initial number of events. Cut 1 requires the single lepton trigger. Cut 2 requires 2 identified leptons. Cut 3 requires the Z-boson mass window. Cut 4 requires 2 reconstructed $a\to b\bar{b}$ candidates. Cut 5a requires 2 identified $a\to b\bar{b}$ candidates in the 1HPC1LPC region. Cut 6a requires the 2 $a\to b\bar{b}$ candidates in the 1HPC1LPC region to be inside the Higgs mass window. Cut 5b requires 2 identified $a\to b\bar{b}$ candidates in the 2HPC region. Cut 6b requires the 2 $a\to b\bar{b}$ candidates in the 2HPC region to be inside the Higgs mass window.

Background yield table for Z+jets, $t\bar{t}$, and rare sources. Observed data yield. Signal $ZH(\rightarrow aa\rightarrow (b\bar{b})(b\bar{b}))$ yield with $m_a = 20\,\text{GeV}$. The signal sample is produced with cross section equals to the standard model $pp\to ZH$, i.e. $0.88\,\text{pb}$, with a branching ratio set to 1 for the $H \rightarrow aa$ decay, whereas the ATLAS figure attached to this entry instead uses the upper-limit branching ratio (smaller than 1). The table includes the yields in two signal regions with leptons consistent with an on-shell Z-boson decay, one with 2 $a\to b\bar{b}$ candidates in the 2HPC region and one with 2 $a\to b\bar{b}$ candidates in the 1HPC1LPC region. The table also includes the yields in four control regions, one with leptons consistent with an on-shell Z-boson decay and 2 $a\to b\bar{b}$ candidates in the Low Purity Category (LPC), and three others where the leptons are not consistent an on-shell Z-boson decay.

Background yield table for Z+jets, $t\bar{t}$, and rare sources. Observed data yield. Signal $ZH(\rightarrow aa\rightarrow (b\bar{b})(b\bar{b}))$ yield with $m_a = 20\,\text{GeV}$. The signal sample is produced with cross section equals to the standard model $pp\to ZH$, i.e. $0.88\,\text{pb}$. The table includes the yields in two signal regions with leptons consistent with an on-shell Z-boson decay, one with 2 $a\to b\bar{b}$ candidates in the 2HPC region and one with 2 $a\to b\bar{b}$ candidates in the 1HPC1LPC region. The table also includes the yields in four control regions, one with leptons consistent with an on-shell Z-boson decay and 2 $a\to b\bar{b}$ candidates in the Low Purity Category (LPC), and three others where the leptons are not consistent an on-shell Z-boson decay.

The observed 95\% C.L. limit for $pp\rightarrow G^{*}_{KK}\rightarrow hh\rightarrow b\bar{b}b\bar{b}$ in the bulk RS model with $k/\bar{M}_{Pl} = 2$, as a function of resonance mass.

The observed 95\% C.L. limit for $pp\rightarrow H\rightarrow hh\rightarrow b\bar{b}b\bar{b}$ with fixed $\Gamma_{H} = 1$ GeV, as a function of resonance mass.

Observed and predicted mTtot distribution in the b-tag category of the 2tau_h channel. Despite listing this as an exclusive final state (as there must be at least one b-jets), there is no explicit selection on the presence of additional light-flavour jets. Please note that the bin content is divided by the bin width in the paper figure, but not in the HepData table. The last bin includes overflows. The combined prediction for A and H bosons with masses of 300, 500 and 800 GeV and $\tan\beta$ = 10 in the hMSSM scenario are also provided.

Observed and predicted mTtot distribution for the b-inclusive selection in the 1l1tau_h channel. Please note that the bin content is divided by the bin width in the paper figure, but not in the HepData table. In the paper, the first bin is cut off at 60 GeV for aesthetics but contains underflows down to 50 GeV as in the HepData table. The last bin includes overflows. The prediction for a SSM Zprime with masses of 1500, 2000 and 2500 GeV are also provided.

Observed and predicted mTtot distribution for the b-inclusive selection in the 2tau_h channel. Please note that the bin content is divided by the bin width in the paper figure, but not in the HepData table. The last bin includes overflows. The prediction for a SSM Zprime with masses of 1500, 2000 and 2500 GeV are also provided.

Observed and expected 95% CL upper limits on the b-associated Higgs boson production cross section times ditau branching fraction as a function of the boson mass.

Observed and expected 95% CL upper limits on the Drell Yan production cross section times ditau branching fraction as a function of the Zprime boson mass.

Observed and expected 95% CL upper limits on the Higgs boson production cross section times ditau branching fraction as a function of the boson mass and the relative strength of the b-associated production.

Ratio of the 95% CL upper limits on the production cross section times branching fraction for alternate Zprime models with respect to the SSM, both observed and expected are shown.

Acceptance, acceptance times efficiency and b-tag category fraction for a scalar boson produced by gluon-gluon fusion as a function of the scalar boson mass.

Acceptance, acceptance times efficiency and b-tag category fraction for a scalar boson produced by b-associated production as a function of the scalar boson mass.

Acceptance and acceptance times efficiency for a heavy gauge boson produced by Drell Yan as a function of the gauge boson mass.

Two dimensional likelihood scan of the gluon-gluon fusion cross section times braching fraction, $\sigma(gg\phi)\times B(\phi\to\tau\tau)$, vs the b-associated production times branching fraction, $\sigma(bb\phi)\times B(\phi\to\tau\tau)$ for the Higgs boson mass ($m_\phi$) indicated in the table. For each mass, 10000 points are scanned. At each point $\Delta(\mathrm{NLL})$ is calculated, defined as the negative-log-likelihood (NLL) of the conditional fit with $\sigma(gg\phi)$ and $\sigma(bb\phi)$ fixed to their values at the point and with the minimum NLL value at any point subtracted. Vaules are provided for the fit to the observed data and to the expected data, which is the sum of Standard Model contributions not including the SM Higgs boson. The best-fit point and the preferred 68% and 95% boundaries are found at $2\Delta(\mathrm{NLL})$ values of 0.0, 2.30 and 5.90, respectively.

Two dimensional likelihood scan of the gluon-gluon fusion cross section times braching fraction, $\sigma(gg\phi)\times B(\phi\to\tau\tau)$, vs the b-associated production times branching fraction, $\sigma(bb\phi)\times B(\phi\to\tau\tau)$ for the Higgs boson mass ($m_\phi$) indicated in the table. For each mass, 10000 points are scanned. At each point $\Delta(\mathrm{NLL})$ is calculated, defined as the negative-log-likelihood (NLL) of the conditional fit with $\sigma(gg\phi)$ and $\sigma(bb\phi)$ fixed to their values at the point and with the minimum NLL value at any point subtracted. Vaules are provided for the fit to the observed data and to the expected data, which is the sum of Standard Model contributions not including the SM Higgs boson. The best-fit point and the preferred 68% and 95% boundaries are found at $2\Delta(\mathrm{NLL})$ values of 0.0, 2.30 and 5.90, respectively.

Two dimensional likelihood scan of the gluon-gluon fusion cross section times braching fraction, $\sigma(gg\phi)\times B(\phi\to\tau\tau)$, vs the b-associated production times branching fraction, $\sigma(bb\phi)\times B(\phi\to\tau\tau)$ for the Higgs boson mass ($m_\phi$) indicated in the table. For each mass, 10000 points are scanned. At each point $\Delta(\mathrm{NLL})$ is calculated, defined as the negative-log-likelihood (NLL) of the conditional fit with $\sigma(gg\phi)$ and $\sigma(bb\phi)$ fixed to their values at the point and with the minimum NLL value at any point subtracted. Vaules are provided for the fit to the observed data and to the expected data, which is the sum of Standard Model contributions not including the SM Higgs boson. The best-fit point and the preferred 68% and 95% boundaries are found at $2\Delta(\mathrm{NLL})$ values of 0.0, 2.30 and 5.90, respectively.

Two dimensional likelihood scan of the gluon-gluon fusion cross section times braching fraction, $\sigma(gg\phi)\times B(\phi\to\tau\tau)$, vs the b-associated production times branching fraction, $\sigma(bb\phi)\times B(\phi\to\tau\tau)$ for the Higgs boson mass ($m_\phi$) indicated in the table. For each mass, 10000 points are scanned. At each point $\Delta(\mathrm{NLL})$ is calculated, defined as the negative-log-likelihood (NLL) of the conditional fit with $\sigma(gg\phi)$ and $\sigma(bb\phi)$ fixed to their values at the point and with the minimum NLL value at any point subtracted. Vaules are provided for the fit to the observed data and to the expected data, which is the sum of Standard Model contributions not including the SM Higgs boson. The best-fit point and the preferred 68% and 95% boundaries are found at $2\Delta(\mathrm{NLL})$ values of 0.0, 2.30 and 5.90, respectively.

Two dimensional likelihood scan of the gluon-gluon fusion cross section times braching fraction, $\sigma(gg\phi)\times B(\phi\to\tau\tau)$, vs the b-associated production times branching fraction, $\sigma(bb\phi)\times B(\phi\to\tau\tau)$ for the Higgs boson mass ($m_\phi$) indicated in the table. For each mass, 10000 points are scanned. At each point $\Delta(\mathrm{NLL})$ is calculated, defined as the negative-log-likelihood (NLL) of the conditional fit with $\sigma(gg\phi)$ and $\sigma(bb\phi)$ fixed to their values at the point and with the minimum NLL value at any point subtracted. Vaules are provided for the fit to the observed data and to the expected data, which is the sum of Standard Model contributions not including the SM Higgs boson. The best-fit point and the preferred 68% and 95% boundaries are found at $2\Delta(\mathrm{NLL})$ values of 0.0, 2.30 and 5.90, respectively.

Two dimensional likelihood scan of the gluon-gluon fusion cross section times braching fraction, $\sigma(gg\phi)\times B(\phi\to\tau\tau)$, vs the b-associated production times branching fraction, $\sigma(bb\phi)\times B(\phi\to\tau\tau)$ for the Higgs boson mass ($m_\phi$) indicated in the table. For each mass, 10000 points are scanned. At each point $\Delta(\mathrm{NLL})$ is calculated, defined as the negative-log-likelihood (NLL) of the conditional fit with $\sigma(gg\phi)$ and $\sigma(bb\phi)$ fixed to their values at the point and with the minimum NLL value at any point subtracted. Vaules are provided for the fit to the observed data and to the expected data, which is the sum of Standard Model contributions not including the SM Higgs boson. The best-fit point and the preferred 68% and 95% boundaries are found at $2\Delta(\mathrm{NLL})$ values of 0.0, 2.30 and 5.90, respectively.

Two dimensional likelihood scan of the gluon-gluon fusion cross section times braching fraction, $\sigma(gg\phi)\times B(\phi\to\tau\tau)$, vs the b-associated production times branching fraction, $\sigma(bb\phi)\times B(\phi\to\tau\tau)$ for the Higgs boson mass ($m_\phi$) indicated in the table. For each mass, 10000 points are scanned. At each point $\Delta(\mathrm{NLL})$ is calculated, defined as the negative-log-likelihood (NLL) of the conditional fit with $\sigma(gg\phi)$ and $\sigma(bb\phi)$ fixed to their values at the point and with the minimum NLL value at any point subtracted. Vaules are provided for the fit to the observed data and to the expected data, which is the sum of Standard Model contributions not including the SM Higgs boson. The best-fit point and the preferred 68% and 95% boundaries are found at $2\Delta(\mathrm{NLL})$ values of 0.0, 2.30 and 5.90, respectively.

Two dimensional likelihood scan of the gluon-gluon fusion cross section times braching fraction, $\sigma(gg\phi)\times B(\phi\to\tau\tau)$, vs the b-associated production times branching fraction, $\sigma(bb\phi)\times B(\phi\to\tau\tau)$ for the Higgs boson mass ($m_\phi$) indicated in the table. For each mass, 10000 points are scanned. At each point $\Delta(\mathrm{NLL})$ is calculated, defined as the negative-log-likelihood (NLL) of the conditional fit with $\sigma(gg\phi)$ and $\sigma(bb\phi)$ fixed to their values at the point and with the minimum NLL value at any point subtracted. Vaules are provided for the fit to the observed data and to the expected data, which is the sum of Standard Model contributions not including the SM Higgs boson. The best-fit point and the preferred 68% and 95% boundaries are found at $2\Delta(\mathrm{NLL})$ values of 0.0, 2.30 and 5.90, respectively.

Two dimensional likelihood scan of the gluon-gluon fusion cross section times braching fraction, $\sigma(gg\phi)\times B(\phi\to\tau\tau)$, vs the b-associated production times branching fraction, $\sigma(bb\phi)\times B(\phi\to\tau\tau)$ for the Higgs boson mass ($m_\phi$) indicated in the table. For each mass, 10000 points are scanned. At each point $\Delta(\mathrm{NLL})$ is calculated, defined as the negative-log-likelihood (NLL) of the conditional fit with $\sigma(gg\phi)$ and $\sigma(bb\phi)$ fixed to their values at the point and with the minimum NLL value at any point subtracted. Vaules are provided for the fit to the observed data and to the expected data, which is the sum of Standard Model contributions not including the SM Higgs boson. The best-fit point and the preferred 68% and 95% boundaries are found at $2\Delta(\mathrm{NLL})$ values of 0.0, 2.30 and 5.90, respectively.

Two dimensional likelihood scan of the gluon-gluon fusion cross section times braching fraction, $\sigma(gg\phi)\times B(\phi\to\tau\tau)$, vs the b-associated production times branching fraction, $\sigma(bb\phi)\times B(\phi\to\tau\tau)$ for the Higgs boson mass ($m_\phi$) indicated in the table. For each mass, 10000 points are scanned. At each point $\Delta(\mathrm{NLL})$ is calculated, defined as the negative-log-likelihood (NLL) of the conditional fit with $\sigma(gg\phi)$ and $\sigma(bb\phi)$ fixed to their values at the point and with the minimum NLL value at any point subtracted. Vaules are provided for the fit to the observed data and to the expected data, which is the sum of Standard Model contributions not including the SM Higgs boson. The best-fit point and the preferred 68% and 95% boundaries are found at $2\Delta(\mathrm{NLL})$ values of 0.0, 2.30 and 5.90, respectively.

Two dimensional likelihood scan of the gluon-gluon fusion cross section times braching fraction, $\sigma(gg\phi)\times B(\phi\to\tau\tau)$, vs the b-associated production times branching fraction, $\sigma(bb\phi)\times B(\phi\to\tau\tau)$ for the Higgs boson mass ($m_\phi$) indicated in the table. For each mass, 10000 points are scanned. At each point $\Delta(\mathrm{NLL})$ is calculated, defined as the negative-log-likelihood (NLL) of the conditional fit with $\sigma(gg\phi)$ and $\sigma(bb\phi)$ fixed to their values at the point and with the minimum NLL value at any point subtracted. Vaules are provided for the fit to the observed data and to the expected data, which is the sum of Standard Model contributions not including the SM Higgs boson. The best-fit point and the preferred 68% and 95% boundaries are found at $2\Delta(\mathrm{NLL})$ values of 0.0, 2.30 and 5.90, respectively.

Two dimensional likelihood scan of the gluon-gluon fusion cross section times braching fraction, $\sigma(gg\phi)\times B(\phi\to\tau\tau)$, vs the b-associated production times branching fraction, $\sigma(bb\phi)\times B(\phi\to\tau\tau)$ for the Higgs boson mass ($m_\phi$) indicated in the table. For each mass, 10000 points are scanned. At each point $\Delta(\mathrm{NLL})$ is calculated, defined as the negative-log-likelihood (NLL) of the conditional fit with $\sigma(gg\phi)$ and $\sigma(bb\phi)$ fixed to their values at the point and with the minimum NLL value at any point subtracted. Vaules are provided for the fit to the observed data and to the expected data, which is the sum of Standard Model contributions not including the SM Higgs boson. The best-fit point and the preferred 68% and 95% boundaries are found at $2\Delta(\mathrm{NLL})$ values of 0.0, 2.30 and 5.90, respectively.

Two dimensional likelihood scan of the gluon-gluon fusion cross section times braching fraction, $\sigma(gg\phi)\times B(\phi\to\tau\tau)$, vs the b-associated production times branching fraction, $\sigma(bb\phi)\times B(\phi\to\tau\tau)$ for the Higgs boson mass ($m_\phi$) indicated in the table. For each mass, 10000 points are scanned. At each point $\Delta(\mathrm{NLL})$ is calculated, defined as the negative-log-likelihood (NLL) of the conditional fit with $\sigma(gg\phi)$ and $\sigma(bb\phi)$ fixed to their values at the point and with the minimum NLL value at any point subtracted. Vaules are provided for the fit to the observed data and to the expected data, which is the sum of Standard Model contributions not including the SM Higgs boson. The best-fit point and the preferred 68% and 95% boundaries are found at $2\Delta(\mathrm{NLL})$ values of 0.0, 2.30 and 5.90, respectively.

Two dimensional likelihood scan of the gluon-gluon fusion cross section times braching fraction, $\sigma(gg\phi)\times B(\phi\to\tau\tau)$, vs the b-associated production times branching fraction, $\sigma(bb\phi)\times B(\phi\to\tau\tau)$ for the Higgs boson mass ($m_\phi$) indicated in the table. For each mass, 10000 points are scanned. At each point $\Delta(\mathrm{NLL})$ is calculated, defined as the negative-log-likelihood (NLL) of the conditional fit with $\sigma(gg\phi)$ and $\sigma(bb\phi)$ fixed to their values at the point and with the minimum NLL value at any point subtracted. Vaules are provided for the fit to the observed data and to the expected data, which is the sum of Standard Model contributions not including the SM Higgs boson. The best-fit point and the preferred 68% and 95% boundaries are found at $2\Delta(\mathrm{NLL})$ values of 0.0, 2.30 and 5.90, respectively.

Two dimensional likelihood scan of the gluon-gluon fusion cross section times braching fraction, $\sigma(gg\phi)\times B(\phi\to\tau\tau)$, vs the b-associated production times branching fraction, $\sigma(bb\phi)\times B(\phi\to\tau\tau)$ for the Higgs boson mass ($m_\phi$) indicated in the table. For each mass, 10000 points are scanned. At each point $\Delta(\mathrm{NLL})$ is calculated, defined as the negative-log-likelihood (NLL) of the conditional fit with $\sigma(gg\phi)$ and $\sigma(bb\phi)$ fixed to their values at the point and with the minimum NLL value at any point subtracted. Vaules are provided for the fit to the observed data and to the expected data, which is the sum of Standard Model contributions not including the SM Higgs boson. The best-fit point and the preferred 68% and 95% boundaries are found at $2\Delta(\mathrm{NLL})$ values of 0.0, 2.30 and 5.90, respectively.

Observed and expected 95% CL upper limits on the gluon-gluon fusion Higgs boson production cross section times ditau branching fraction as a function of the boson mass.

DM exclusion limit in the two-dimensional phase space of WIMP mass m_&chi; vs mediator mass m_med determined using the combined ee+&mu;&mu; channel. Both the observed and expected limits are presented, and the 1&sigma; uncertainty band for the expected limits is also provided. Regions bounded by the limit curves are excluded at the 95% CL. The grey line labelled with "m_med = 2m_&chi;'' indicates the kinematic threshold where the mediator can decay on-shell into WIMPs, and the other grey line gives the perturbative limit (arXiv 1603.04156). The relic density line (arXiv 1603.04156) illustrates the combination of m_&chi; and m_med that would explain the observed DM relic density.

DM exclusion limit in the two-dimensional phase space of WIMP mass m_&chi; vs mediator mass m_med determined using the combined ee+&mu;&mu; channel. Both the observed and expected limits are presented, and the 1&sigma; uncertainty band for the expected limits is also provided. Regions bounded by the limit curves are excluded at the 95% CL. The grey line labelled with "m_med = 2m_&chi;'' indicates the kinematic threshold where the mediator can decay on-shell into WIMPs, and the other grey line gives the perturbative limit (arXiv 1603.04156). The relic density line (arXiv 1603.04156) illustrates the combination of m_&chi; and m_med that would explain the observed DM relic density.

DM exclusion limit in the two-dimensional phase space of WIMP mass m_&chi; vs mediator mass m_med determined using the combined ee+&mu;&mu; channel. Both the observed and expected limits are presented, and the 1&sigma; uncertainty band for the expected limits is also provided. Regions bounded by the limit curves are excluded at the 95% CL. The grey line labelled with "m_med = 2m_&chi;'' indicates the kinematic threshold where the mediator can decay on-shell into WIMPs, and the other grey line gives the perturbative limit (arXiv 1603.04156). The relic density line (arXiv 1603.04156) illustrates the combination of m_&chi; and m_med that would explain the observed DM relic density.

The 95% CL upper limits on the production cross-section of the ZH process with the prompt Z&rarr; ee and Z&rarr; &mu;&mu; decays and the invisible Higgs boson decays as a function of m_H, obtained from the combined ee+&mu;&mu; channel. The observed and expected limits are given, as well as the &pm;1&sigma; and &pm;2&sigma; error bands on the expected limits. The signal process is only modelled for the qq initial state, and the theory uncertainties on the signal prediction are not considered.

DM exclusion limit in the two-dimensional phase space of m_&chi; vs m_med set using the combined ee+&mu;&mu; channel. The DM model assumes a vector mediator, a fermionic WIMP, and the coupling parameters g_q = 0.25 and g_&chi; = 1. Both the observed and expected limits are presented, and the 1&sigma; uncertainty band for the expected limits is also provided. Regions bounded by the limit curves are excluded at the 95% CL. The grey line labelled with "m_med = 2m_&chi;'' indicates the kinematic threshold where the mediator can decay on-shell into WIMPs. The relic density line [arXiv 1603.04156] illustrates the combination of m_&chi; and m_med that would explain the observed DM relic density.

DM exclusion limit in the two-dimensional phase space of m_&chi; vs m_med set using the combined ee+&mu;&mu; channel. The DM model assumes a vector mediator, a fermionic WIMP, and the coupling parameters g_q = 0.25 and g_&chi; = 1. Both the observed and expected limits are presented, and the 1&sigma; uncertainty band for the expected limits is also provided. Regions bounded by the limit curves are excluded at the 95% CL. The grey line labelled with "m_med = 2m_&chi;'' indicates the kinematic threshold where the mediator can decay on-shell into WIMPs. The relic density line [arXiv 1603.04156] illustrates the combination of m_&chi; and m_med that would explain the observed DM relic density.

DM exclusion limit in the two-dimensional phase space of m_&chi; vs m_med set using the combined ee+&mu;&mu; channel. The DM model assumes a vector mediator, a fermionic WIMP, and the coupling parameters g_q = 0.25 and g_&chi; = 1. Both the observed and expected limits are presented, and the 1&sigma; uncertainty band for the expected limits is also provided. Regions bounded by the limit curves are excluded at the 95% CL. The grey line labelled with "m_med = 2m_&chi;'' indicates the kinematic threshold where the mediator can decay on-shell into WIMPs. The relic density line [arXiv 1603.04156] illustrates the combination of m_&chi; and m_med that would explain the observed DM relic density.

DM exclusion limit in the two-dimensional phase space of m_&chi; vs m_med set using the combined ee+&mu;&mu; channel. The DM model assumes a vector mediator, a fermionic WIMP, and the coupling parameters g_q = 0.25 and g_&chi; = 1. Both the observed and expected limits are presented, and the 1&sigma; uncertainty band for the expected limits is also provided. Regions bounded by the limit curves are excluded at the 95% CL. The grey line labelled with "m_med = 2m_&chi;'' indicates the kinematic threshold where the mediator can decay on-shell into WIMPs. The relic density line [arXiv 1603.04156] illustrates the combination of m_&chi; and m_med that would explain the observed DM relic density.

Limits on the WIMP and proton scattering cross section as a function of WIMP mass, obtained for simplified dark matter models with an axial-vector mediator (left) and a vector mediator (right). The solid black line shows the observed limit at the 90% confidence level from this search. The left plot incorporates the latest results from these competitive direct-search experiments, LUX [arXiv 1705.03380], PICO-2L [Phys. Rev. Lett. 114 (2015) 231302], and PICO-60 [arXiv 1702.07666]. The right plot incorporates the latest results from the following experiments, CRESST-II [Eur. Phys. J. C76 (2016) 25], CDMSlite [Phys. Rev. Lett. 116 (2016) 071301], PandaX-II [Phys. Rev. Lett. 117 (2016) 121303], LUX [Phys. Rev. Lett. 118 (2017) 021303], and XENON1T [arXiv 1705.06655].

Limits on the WIMP and proton scattering cross section as a function of WIMP mass, obtained for simplified dark matter models with an axial-vector mediator (left) and a vector mediator (right). The solid black line shows the observed limit at the 90% confidence level from this search. The left plot incorporates the latest results from these competitive direct-search experiments, LUX [arXiv 1705.03380], PICO-2L [Phys. Rev. Lett. 114 (2015) 231302], and PICO-60 [arXiv 1702.07666]. The right plot incorporates the latest results from the following experiments, CRESST-II [Eur. Phys. J. C76 (2016) 25], CDMSlite [Phys. Rev. Lett. 116 (2016) 071301], PandaX-II [Phys. Rev. Lett. 117 (2016) 121303], LUX [Phys. Rev. Lett. 118 (2017) 021303], and XENON1T [arXiv 1705.06655].

Lower limits at 95% CL on the mass of $H_{\mathrm{L}}^{\pm \pm}$ and $H_{\mathrm{R}}^{\pm \pm}$ bosons, assuming a 100% branching ratio to $e^{\pm}e^{\pm}$, $e^{\pm}/\mu^{\pm}$ and $\mu^{\pm}\mu^{\pm}$ pairs. The $1\sigma$ variations are also shown for the expected limits.

Expected CLs values in the High mA SR, mA vs ma signal grid.

Observed CLs values in the High mA SR, mA vs ma signal grid.

Expected CLs values in the Low mA SR, mA vs ma signal grid.

Observed CLs values in the High mA SR, mA vs ma signal grid.

CLs+b values in the Low mA SR, mA vs tanB signal grid.

CLs+b values in the High mA SR, mA vs ma signal grid.

CLs+b values in the Low mA SR, mA vs ma signal grid.

Cut flow of the 2HDM+a signal points, gluon–gluon fusion production, Low mA SR. tanB = 1, $sin\theta$ = 0.35. The first two entries in the tables are number of raw MC events, third entry is theoretical prediction, and all other lines include the correct weights. Note that during the generation Higgs boson’s branching ratio to taus has been set to 1. An additional factor of 0.0627 is used to account for that, starting from the ‘Initial’ entry.

Cut flow of the 2HDM+a signal points, gluon–gluon fusion production, High mA SR. tanB = 1, $sin\theta$ = 0.35. The first two entries in the tables are number of raw MC events, third entry is theoretical prediction, and all other lines include the correct weights. Note that during the generation Higgs boson’s branching ratio to taus has been set to 1. An additional factor of 0.0627 is used to account for that, starting from the ‘Initial’ entry.

Cut flow of the 2HDM+a signal points, bb annihilation production, Low mA SR. tanB = 1, $sin\theta$ = 0.35. The first two entries in the tables are number of raw MC events, third entry is theoretical prediction, and all other lines include the correct weights. Note that during the generation Higgs boson’s branching ratio to taus has been set to 1. An additional factor of 0.0627 is used to account for that, starting from the ‘Initial’ entry.

Cut flow of the 2HDM+a signal points, bb annihilation production, High mA SR. tanB = 1, $sin\theta$ = 0.35. The first two entries in the tables are number of raw MC events, third entry is theoretical prediction, and all other lines include the correct weights. Note that during the generation Higgs boson’s branching ratio to taus has been set to 1. An additional factor of 0.0627 is used to account for that, starting from the ‘Initial’ entry.

A comparison of the observed and expected yields in the four bins of the Low mA SR.

A comparison of the observed and expected yields in the two bins of the High mA SR.

Expected exclusion contours at 95% CL as a function of mA and ma.

Observed exclusion contours at 95% CL as a function of mA and ma.

Expected +- 1sigma exclusion contours at 95% CL as a function of mA and ma.

Expected exclusion contours at 95% CL as a function of mA and ma.

Observed exclusion contours at 95% CL as a function of mA and ma.

Expected +- 1sigma exclusion contours at 95% CL as a function of mA and ma.

Acceptance, High mA SR, mA vs tanB grid, 400-750 GeV, bb prod

Acceptance, High mA SR, mA vs tanB grid, >750 GeV, bb prod

Acceptance, Low mA SR, mA vs tanB grid, 100-250 GeV, bb prod

Acceptance, Low mA SR, mA vs tanB grid, 250-400 GeV, bb prod

Acceptance, Low mA SR, mA vs tanB grid, 400-550 GeV, bb prod

Acceptance, Low mA SR, mA vs tanB grid, >550 GeV, bb prod

Acceptance, High mA SR, mA vs ma grid, 400-750 GeV, bb prod

Acceptance, High mA SR, mA vs ma grid, >750 GeV, bb prod

Acceptance, Low mA SR, mA vs ma grid, 100-250 GeV, bb prod

Acceptance, Low mA SR, mA vs ma grid, 250-400 GeV, bb prod

Acceptance, Low mA SR, mA vs ma grid, 400-550 GeV, bb prod

Acceptance, Low mA SR, mA vs ma grid, >550 GeV, bb prod

Acceptance, High mA SR, mA vs tanB grid, 400-750 GeV, ggF prod

Acceptance, High mA SR, mA vs tanB grid, >750 GeV, ggF prod

Acceptance, Low mA SR, mA vs tanB grid, 100-250 GeV, ggF prod

Acceptance, Low mA SR, mA vs tanB grid, 250-400 GeV, ggF prod

Acceptance, Low mA SR, mA vs tanB grid, 400-550 GeV, ggF prod

Acceptance, Low mA SR, mA vs tanB grid, >550 GeV, ggF prod

Acceptance, High mA SR, mA vs ma grid, 400-750 GeV, ggF prod

Acceptance, High mA SR, mA vs ma grid, >750 GeV, ggF prod

Acceptance, Low mA SR, mA vs ma grid, 100-250 GeV, ggF prod

Acceptance, Low mA SR, mA vs ma grid, 250-400 GeV, ggF prod

Acceptance, Low mA SR, mA vs ma grid, 400-550 GeV, ggF prod

Acceptance, Low mA SR, mA vs ma grid, >550 GeV, ggF prod

Efficiency, High mA SR, mA vs tanB grid, 400-750 GeV, bb prod

Efficiency, High mA SR, mA vs tanB grid, >750 GeV, bb prod

Efficiency, Low mA SR, mA vs tanB grid, 100-250 GeV, bb prod

Efficiency, Low mA SR, mA vs tanB grid, 250-400 GeV, bb prod

Efficiency, Low mA SR, mA vs tanB grid, 400-550 GeV, bb prod

Efficiency, Low mA SR, mA vs tanB grid, >550 GeV, bb prod

Efficiency, High mA SR, mA vs ma grid, 400-750 GeV, bb prod

Efficiency, High mA SR, mA vs ma grid, >750 GeV, bb prod

Efficiency, Low mA SR, mA vs ma grid, 100-250 GeV, bb prod

Efficiency, Low mA SR, mA vs ma grid, 250-400 GeV, bb prod

Efficiency, Low mA SR, mA vs ma grid, 400-550 GeV, bb prod

Efficiency, Low mA SR, mA vs ma grid, >550 GeV, bb prod

Efficiency, High mA SR, mA vs tanB grid, 400-750 GeV, ggF prod

Efficiency, High mA SR, mA vs tanB grid, >750 GeV, ggF prod

Efficiency, Low mA SR, mA vs tanB grid, 100-250 GeV, ggF prod

Efficiency, Low mA SR, mA vs tanB grid, 250-400 GeV, ggF prod

Efficiency, Low mA SR, mA vs tanB grid, 400-550 GeV, ggF prod

Efficiency, Low mA SR, mA vs tanB grid, >550 GeV, ggF prod

Efficiency, High mA SR, mA vs ma grid, 400-750 GeV, ggF prod

Efficiency, High mA SR, mA vs ma grid, >750 GeV, ggF prod

Efficiency, Low mA SR, mA vs ma grid, 100-250 GeV, ggF prod

Efficiency, Low mA SR, mA vs ma grid, 250-400 GeV, ggF prod

Efficiency, Low mA SR, mA vs ma grid, 400-550 GeV, ggF prod

Efficiency, Low mA SR, mA vs ma grid, >550 GeV, ggF prod

Expected +- 1 sigma 95% CL exclusion limit for the Zprime-2HDM model.

Expected +- 2 sigma 95% CL exclusion limit for the Zprime-2HDM model.

Observed 95% CL exclusion limit for the 2HDM+a model ggF production.

Expected 95% CL exclusion limit for the 2HDM+a model ggF production.

Expected +- 1 sigma 95% CL exclusion limit for the 2HDM+a model ggF production.

Expected +- 2 sigma 95% CL exclusion limit for the 2HDM+a model ggF production.

Observed 95% CL exclusion limit for the 2HDM+a model bbA production.

Expected 95% CL exclusion limit for the 2HDM+a model bbA production.

Expected +- 1 sigma 95% CL exclusion limit for the 2HDM+a model bbA production.

Expected +- 2 sigma 95% CL exclusion limit for the 2HDM+a model bbA production.

Observed 95% CL exclusion limit for the Zprime-2HDM model with the benchmark used in arXiv:1707.01302.

Expected 95% CL exclusion limit for the Zprime-2HDM model with the benchmark used in arXiv:1707.01302.

Expected +- 1 sigma 95% CL exclusion limit for the Zprime-2HDM model with the benchmark used in arXiv:1707.01302.

Expected +- 2 sigma 95% CL exclusion limit for the Zprime-2HDM model with the benchmark used in arXiv:1707.01302.

Expected and observed upper limits at 95% CL on cross-section for Zprime-2HDM model.

Expected and observed upper limits at 95% CL on cross-section for ggF producton in the 2HDM+a model.

Expected and observed upper limits at 95% CL on cross-section for bbA producton in the 2HDM+a model.

Model-independent upper limits on the visible cross-section $σ_{vis, $h(\bar{b})+DM} ≡ σ_{h+DM} \times B(h \to b\bar{b}) \times \mathcal{A} \times \epsilon$ in the different signal regions.

Theory cross-section for Zprime-2HDM model.

Theory cross-section for bbA production in the 2HDM+a model.

Theory cross-section for ggF production in the 2HDM+a model.

Distribution of Higgs boson candidate mass in 2b region with MET=150-200 GeV.

Distribution of Higgs boson candidate mass in 2b region with MET=200-350 GeV.

Distribution of Higgs boson candidate mass in 2b region with MET=350-500 GeV.

Distribution of Higgs boson candidate mass in 2b region with MET=500-750 GeV.

Distribution of Higgs boson candidate mass in 2b region with MET > 750 GeV.

Distribution of Higgs boson candidate mass in 3b region with MET=150-200 GeV.

Distribution of Higgs boson candidate mass in 3b region with MET=200-350 GeV.

Distribution of Higgs boson candidate mass in 3b region with MET=350-500 GeV.

Distribution of Higgs boson candidate mass in 3b region with MET > 500 GeV.

Yields in 1-lepton control region.

Yields in 2-lepton control region.

MET distribution in 2b region of the 0-lepton channel.

MET distribution in 3b region of the 0-lepton channel.

Expected signal yields after certain selection cuts in 2b region with MET=150-200 GeV.

Expected signal yields after certain selection cuts in 2b region with MET=200-350 GeV.

Expected signal yields after certain selection cuts in 2b region with MET=350-500 GeV.

Expected signal yields after certain selection cuts in 2b region with MET=500-750 GeV.

Expected signal yields after certain selection cuts in 2b region with MET > 750 GeV.

Expected signal yields after certain selection cuts in 3b region with MET=150-200 GeV.

Expected signal yields after certain selection cuts in 3b region with MET=200-350 GeV.

Expected signal yields after certain selection cuts in 3b region with MET=350-500 GeV.

Expected signal yields after certain selection cuts in 3b region with MET > 500 GeV.

Acceptance times efficiency for bbA production in the 2HDM+a model - 2b region with MET=150-200 GeV.

Acceptance times efficiency for bbA production in the 2HDM+a model - 2b region with MET=200-350 GeV.

Acceptance times efficiency for bbA production in the 2HDM+a model - 2b region with MET=350-500 GeV.

Acceptance times efficiency for bbA production in the 2HDM+a model - 2b region with MET=500-750 GeV.

Acceptance times efficiency for bbA production in the 2HDM+a model - 2b region with MET > 750 GeV.

Acceptance times efficiency for bbA production in the 2HDM+a model - 3b region with MET=150-200 GeV.

Acceptance times efficiency for bbA production in the 2HDM+a model - 3b region with MET=200-350 GeV.

Acceptance times efficiency for bbA production in the 2HDM+a model - 3b region with MET=350-500 GeV.

Acceptance times efficiency for bbA production in the 2HDM+a model - 3b region with MET>500 GeV.

Acceptance times efficiency for ggF production in the 2HDM+a model - 2b region with MET=150-200 GeV.

Acceptance times efficiency for ggF production in the 2HDM+a model - 2b region with MET=200-350 GeV.

Acceptance times efficiency for ggF production in the 2HDM+a model - 2b region with MET=350-500 GeV.

Acceptance times efficiency for ggF production in the 2HDM+a model - 2b region with MET=500-750 GeV.

Acceptance times efficiency for ggF production in the 2HDM+a model - 2b region with MET > 750 GeV.

Acceptance times efficiency for ggF production in the 2HDM+a model - 3b region with MET=150-200 GeV.

Acceptance times efficiency for ggF production in the 2HDM+a model - 3b region with MET=200-350 GeV.

Acceptance times efficiency for ggF production in the 2HDM+a model - 3b region with MET=350-500 GeV.

Acceptance times efficiency for ggF production in the 2HDM+a model - 3b region with MET > 500 GeV.

Acceptance times efficiency for ggF production in the Zprime-2HDM model - 2b region with MET=150-200 GeV.

Acceptance times efficiency for ggF production in the Zprime-2HDM model - 2b region with MET=200-350 GeV.

Acceptance times efficiency for ggF production in the Zprime-2HDM model - 2b region with MET=350-500 GeV.

Acceptance times efficiency for ggF production in the Zprime-2HDM model - 2b region with MET=500-750 GeV.

Acceptance times efficiency for ggF production in the Zprime-2HDM model - 2b region with MET > 750 GeV.

Acceptance times efficiency for ggF production in the Zprime-2HDM model - 3b region with MET=150-200 GeV.

Acceptance times efficiency for ggF production in the Zprime-2HDM model - 3b region with MET=200-350 GeV.

Acceptance times efficiency for ggF production in the Zprime-2HDM model - 3b region with MET=350-500 GeV.

Acceptance times efficiency for ggF production in the Zprime-2HDM model - 3b region with MET > 500 GeV.

Distribution of the transverse mass of the $W$-candidate $m_{\rm T}^{W}$ for the one lepton and two $b$-jets channel in SRlbb-1 and SRlbb-2, $m_{bb}$ sidebands. The background histograms are obtained from the background-only fit, and their uncertainty represents the total background uncertainty after the fit. The last bin includes overflow.

Distribution of the number of $b$-jets for the one lepton and two $b$-jets channel in SRlbb-1 and SRlbb-2 without the b-jet multiplicity requirement, central $m_{bb}$ bin. The background histograms are obtained from the background-only fit, and their uncertainty represents the total background uncertainty after the fit.

Distribution of the invariant mass of the $b$-jets $m_{bb}$ for the one lepton and two $b$-jets channel in the SRlbb-1 and SRlbb-2. The background histograms are obtained from the background-only fit, and their uncertainty represents the total background uncertainty after the fit.

Distribution of missing transverse momentum $E_{\rm T}^{\rm miss}$ in the one lepton and two photons signal regions SR$\ell\gamma\gamma$-1 and SR$\ell\gamma\gamma$-2 for the Higgs-mass window ($120\lt m_{\gamma\gamma} \lt 130$ GeV). The final column (background) is a simulation-based cross check. The contributions from non-Higgs backgrounds are scaled by 10 GeV / 50 GeV = 0.2 from the $m_{\gamma\gamma}$ sideband ($100 \lt m_{\gamma\gamma} \lt 120$ GeV and $130 \lt m_{\gamma\gamma} \lt 160$ GeV) into the Higgs-mass window. The last bin includes overflow. The distributions of a signal hypothesis are also shown for $m(\tilde{\chi}^{\pm}_{1}\tilde{\chi}^{0}_{2},\tilde{\chi}^{0}_{1})=(165,35)$ GeV.

Distribution of the azimuth difference between the $W$ and Higgs boson candidates $\Delta\phi(W,h)$ in the one lepton and two photons signal regions SR$\ell\gamma\gamma$-1 and SR$\ell\gamma\gamma$-2 for the Higgs-mass window ($120\lt m_{\gamma\gamma} \lt 130$ GeV). The final column (background) is a simulation-based cross check. The contributions from non-Higgs backgrounds are scaled by 10 GeV / 50 GeV = 0.2 from the $m_{\gamma\gamma}$ sideband ($100 \lt m_{\gamma\gamma} \lt 120$ GeV and $130 \lt m_{\gamma\gamma} \lt 160$ GeV) into the Higgs-mass window. The distributions of a signal hypothesis are also shown for $m(\tilde{\chi}^{\pm}_{1}\tilde{\chi}^{0}_{2},\tilde{\chi}^{0}_{1})=(165,35)$ GeV.

Distribution of the transverse mass of the $W$ and photon system $m_{\rm{T}}^{W\gamma_1}$ in the one lepton and two photons signal regions SR$\ell\gamma\gamma$-1 and SR$\ell\gamma\gamma$-2 for the Higgs-mass window ($120\lt m_{\gamma\gamma} \lt 130$ GeV). The final column (background) is a simulation-based cross check. The contributions from non-Higgs backgrounds are scaled by 10 GeV / 50 GeV = 0.2 from the $m_{\gamma\gamma}$ sideband ($100 \lt m_{\gamma\gamma} \lt 120$ GeV and $130 \lt m_{\gamma\gamma} \lt 160$ GeV) into the Higgs-mass window. The last bin includes overflow. The distributions of a signal hypothesis are also shown for $m(\tilde{\chi}^{\pm}_{1}\tilde{\chi}^{0}_{2},\tilde{\chi}^{0}_{1})=(165,35)$ GeV.

Distribution of the transverse mass of the $W$ and photon system $m_{\rm{T}}^{W\gamma_2}$ in the one lepton and two photons signal regions SR$\ell\gamma\gamma$-1 and SR$\ell\gamma\gamma$-2 for the Higgs-mass window ($120\lt m_{\gamma\gamma} \lt 130$ GeV). The final column (background) is a simulation-based cross check. The contributions from non-Higgs backgrounds are scaled by 10 GeV / 50 GeV = 0.2 from the $m_{\gamma\gamma}$ sideband ($100 \lt m_{\gamma\gamma} \lt 120$ GeV and $130 \lt m_{\gamma\gamma} \lt 160$ GeV) into the Higgs-mass window. The last bin includes overflow. The distributions of a signal hypothesis are also shown for $m(\tilde{\chi}^{\pm}_{1}\tilde{\chi}^{0}_{2},\tilde{\chi}^{0}_{1})=(165,35)$ GeV.

Results of the background-only fit to the diphoton invariant mass, $m_{\gamma\gamma}$, distribution in the one lepton and two photons signal region SR$l\gamma\gamma$-1. The contributions from SM Higgs boson production are constrained to the MC prediction and associated systematic uncertainties. The fit is performed on events with 100 GeV $ \lt m_{\gamma\gamma} \lt $ 160 GeV, with events in SR$l\gamma\gamma$-1 or SR$l\gamma\gamma$-2 in the Higgs-mass window (120 GeV $\le m_{\gamma\gamma} \le$ 130 GeV) excluded from the fit.

Results of the background-only fit to the diphoton invariant mass, $m_{\gamma\gamma}$, distribution in the one lepton and two photons signal region SR$l\gamma\gamma$-2. The contributions from SM Higgs boson production are constrained to the MC prediction and associated systematic uncertainties. The fit is performed on events with 100 GeV $ \lt m_{\gamma\gamma} \lt $ 160 GeV, with events in SR$l\gamma\gamma$-1 or SR$l\gamma\gamma$-2 in the Higgs-mass window (120 GeV $\le m_{\gamma\gamma} \le$ 130 GeV) excluded from the fit.

Results of the background-only fit to the diphoton invariant mass, $m_{\gamma\gamma}$, distribution in the one lepton and two photons validation region VR$l\gamma\gamma$-1. The contributions from SM Higgs boson production are constrained to the MC prediction and associated systematic uncertainties. The fit is performed on events with 100 GeV $ \lt m_{\gamma\gamma} \lt $ 160 GeV, with events in SR$l\gamma\gamma$-1 or SR$l\gamma\gamma$-2 in the Higgs-mass window (120 GeV $\le m_{\gamma\gamma} \le$ 130 GeV) excluded from the fit.

Results of the background-only fit to the diphoton invariant mass, $m_{\gamma\gamma}$, distribution in the one lepton and two photons validation region VR$l\gamma\gamma$-2. The contributions from SM Higgs boson production are constrained to the MC prediction and associated systematic uncertainties. The fit is performed on events with 100 GeV $ \lt m_{\gamma\gamma} \lt $ 160 GeV, with events in SR$l\gamma\gamma$-1 or SR$l\gamma\gamma$-2 in the Higgs-mass window (120 GeV $\le m_{\gamma\gamma} \le$ 130 GeV) excluded from the fit.

Distribution of effective mass $m_{\rm eff}$ in the validation region of the same-sign $e\mu$ channel. This validation region is defined by requiring one, two, or three jets, and reversing the $m_{lj}$, $m_{ljj}$ criteria. The distribution of a signal hypothesis is also shown.

Distribution of effective mass $m_{\rm eff}$ for the same-sign dilepton channel in the signal regions with one jet SR$ll$-1. SR$ll$-1 is the sum of SR$ee$-1, SR$e\mu$-1, and SR$\mu\mu$-1. All selection criteria are applied, except for the one on $m_{\rm eff}$. The distributions of a signal hypothesis are also shown for $m(\tilde{\chi}^{\pm}_{1}\tilde{\chi}^{0}_{2},\tilde{\chi}^{0}_{1})=(130,0)$ GeV. The last bin includes overflow.

Distribution of effective mass $m_{\rm eff}$ for the same-sign dilepton channel in the signal regions with two or three jets SR$ll$-2. SR$ll$-2 is the sum of SR$ee$-2, SR$e\mu$-2, and SR$\mu\mu$-2. All selection criteria are applied, except for the one on $m_{\rm eff}$. The distributions of a signal hypothesis are also shown for $m(\tilde{\chi}^{\pm}_{1}\tilde{\chi}^{0}_{2},\tilde{\chi}^{0}_{1})=(130,0)$ GeV. The last bin includes overflow.

Distribution of largest transverse mass $m_{\rm T}^{\rm max}$ for the same-sign dilepton channel in the signal regions with one jet SR$ll$-1. SR$ll$-1 is the sum of SR$ee$-1, SR$e\mu$-1, and SR$\mu\mu$-1. All selection criteria are applied, except for the one on $m_{\rm T}^{\rm max}$. The distributions of a signal hypothesis are also shown for $m(\tilde{\chi}^{\pm}_{1}\tilde{\chi}^{0}_{2},\tilde{\chi}^{0}_{1})=(130,0)$ GeV. The last bin includes overflow.

Distribution of largest transverse mass $m_{\rm T}^{\rm max}$ for the same-sign dilepton channel in the signal regions with two or three jets SR$ll$-2. SR$ll$-2 is the sum of SR$ee$-2, SR$e\mu$-2, and SR$\mu\mu$-2. All selection criteria are applied, except for the one on $m_{\rm T}^{\rm max}$. The distributions of a signal hypothesis are also shown for $m(\tilde{\chi}^{\pm}_{1}\tilde{\chi}^{0}_{2},\tilde{\chi}^{0}_{1})=(130,0)$ GeV. The last bin includes overflow.

Distribution of invariant mass of lepton and jet $m_{lj}$ for the same-sign dilepton channel in the signal regions with one jet SR$ll$-1. SR$ll$-1 is the sum of SR$ee$-1, SR$e\mu$-1, and SR$\mu\mu$-1. All selection criteria are applied, except for the one on $m_{lj}$. The distributions of a signal hypothesis are also shown for $m(\tilde{\chi}^{\pm}_{1}\tilde{\chi}^{0}_{2},\tilde{\chi}^{0}_{1})=(130,0)$ GeV. The last bin includes overflow.

Distribution of invariant mass of lepton and di-jet $m_{ljj}$ for the same-sign dilepton channel in the signal regions with two or three jets SR$ll$-2. SR$ll$-2 is the sum of SR$ee$-2, SR$e\mu$-2, and SR$\mu\mu$-2. All selection criteria are applied, except for the one on $m_{ljj}$. The distributions of a signal hypothesis are also shown for $m(\tilde{\chi}^{\pm}_{1}\tilde{\chi}^{0}_{2},\tilde{\chi}^{0}_{1})=(130,0)$ GeV. The last bin includes overflow.

One lepton and two $b$-jets channel: observed and expected 95% CL upper limits on the cross section normalised by the simplified model prediction as a function of the common mass $m_{\tilde{\chi}_1^\pm \tilde{\chi}^0_2}$ for $m_{\tilde{\chi}^0_1}=0$.

One lepton and two photons channel: observed and expected 95% CL upper limits on the cross section normalised by the simplified model prediction as a function of the common mass $m_{\tilde{\chi}_1^\pm \tilde{\chi}^0_2}$ for $m_{\tilde{\chi}^0_1}=0$.

Same-sign dilepton channel: observed and expected 95% CL upper limits on the cross section normalised by the simplified model prediction as a function of the common mass $m_{\tilde{\chi}_1^\pm \tilde{\chi}^0_2}$ for $m_{\tilde{\chi}^0_1}=0$.

Combination: observed and expected 95% CL upper limits on the cross section normalised by the simplified model prediction as a function of the common mass $m_{\tilde{\chi}_1^\pm \tilde{\chi}^0_2}$ for $m_{\tilde{\chi}^0_1}=0$. This combination is obtained using the result from the ATLAS three-lepton search, J. High Energy Phys. 04 (2014) 169, in addition to the three channels reported in this paper.

One lepton and two b-jets channel: Expected 95% CL exclusion contour for chargino neutralino production via Wh.

One lepton and two b-jets channel: Observed 95% CL exclusion contour for chargino neutralino production via Wh.

One lepton and two photons channel: Expected 95% CL exclusion contour for chargino neutralino production via Wh.

One lepton and two photons channel: Observed 95% CL exclusion contour for chargino neutralino production via Wh.

Same-sign dilepton channel: Expected 95% CL exclusion contour for chargino neutralino production via Wh.

Same-sign dilepton channel: Observed 95% CL exclusion contour for chargino neutralino production via Wh.

Combination: Expected 95% CL exclusion contour for chargino neutralino production via Wh.

Combination: Observed 95% CL exclusion contour for chargino neutralino production via Wh.

Combination: excluded model cross-section at 95% CL IN PB.

Acceptance SR$\ell\gamma\gamma$-1.

Acceptance SR$\ell\gamma\gamma$-2.

Efficiency SR$\ell\gamma\gamma$-1.

Efficiency SR$\ell\gamma\gamma$-2.

$m_{\gamma\gamma}$ distribution in the SR$l\gamma\gamma$-1 region for the full $m_{\gamma\gamma}$ window. The last column (background) is from a simulation-based cross check. $Z\gamma$ events, with $Z\rightarrow{ee}$ and one electron mis-identified as a photon, are included in the $Z\gamma(\gamma)$ background. The distributions of a signal hypothesis are also shown for $m(\tilde{\chi}^{\pm}_{1}\tilde{\chi}^{0}_{2},\tilde{\chi}^{0}_{1})=(165,35)$ GeV.

$m_{\gamma\gamma}$ distribution in the SR$l\gamma\gamma$-2 region for the full $m_{\gamma\gamma}$ window. The last column (background) is from a simulation-based cross check. $Z\gamma$ events, with $Z\rightarrow{ee}$ and one electron mis-identified as a photon, are included in the $Z\gamma(\gamma)$ background. The distributions of a signal hypothesis are also shown for $m(\tilde{\chi}^{\pm}_{1}\tilde{\chi}^{0}_{2},\tilde{\chi}^{0}_{1})=(165,35)$ GeV.

$m_{\gamma\gamma}$ distribution in the VR$l\gamma\gamma$-1 region for the full $m_{\gamma\gamma}$ window. The last column (background) is from a simulation-based cross check. $Z\gamma$ events, with $Z\rightarrow{ee}$ and one electron mis-identified as a photon, are included in the $Z\gamma(\gamma)$ background. The distributions of a signal hypothesis are also shown for $m(\tilde{\chi}^{\pm}_{1}\tilde{\chi}^{0}_{2},\tilde{\chi}^{0}_{1})=(165,35)$ GeV.

$m_{\gamma\gamma}$ distribution in the VR$l\gamma\gamma$-2 region for the full $m_{\gamma\gamma}$ window. The last column (background) is from a simulation-based cross check. $Z\gamma$ events, with $Z\rightarrow{ee}$ and one electron mis-identified as a photon, are included in the $Z\gamma(\gamma)$ background. The distributions of a signal hypothesis are also shown for $m(\tilde{\chi}^{\pm}_{1}\tilde{\chi}^{0}_{2},\tilde{\chi}^{0}_{1})=(165,35)$ GeV.

Acceptance for SRlbb-1, central $m_{bb}$ bin.

Acceptance for SRlbb-1, $m_{bb}$ sideband.

Acceptance for SRlbb-2, central $m_{bb}$ bin.

Acceptance for SRlbb-2, $m_{bb}$ sideband.

Efficiency for SRlbb-1, central $m_{bb}$ bin.

Efficiency for SRlbb-1, $m_{bb}$ sideband.

Efficiency for SRlbb-2, central $m_{bb}$ bin.

Efficiency for SRlbb-2, $m_{bb}$ sideband.

Acceptance for the same-sign $ee$ channel with one jet.

Acceptance for the same-sign $ee$ channel with two or three jets.

Acceptance for the same-sign $e\mu$ channel with one jet.

Acceptance for the same-sign $e\mu$ channel with two or three jets.

Acceptance for the same-sign $\mu\mu$ channel with one jet.

Acceptance for the same-sign $\mu\mu$ channel with two or three jets.

Efficiency for the same-sign $ee$ channel with one jet.

Efficiency for the same-sign $ee$ channel with two or three jets.

Efficiency for the same-sign $e\mu$ channel with one jet.

Efficiency for the same-sign $e\mu$ channel with two or three jets.

Efficiency for the same-sign $\mu\mu$ channel with one jet.

Efficiency for the same-sign $\mu\mu$ channel with two or three jets.