The observed and expected 95% CL upper limits on the production cross section times branching ratio for SM non-resonant Higgs pair production.

Bin-to-bin correlation in the measured cross section for Z+jets as a function of inclusive jet multiplicity, $N_{\text{jets}}^{\text{min}}$.

Measured cross section for Z+jets as a function of the transverse momementum of the Z boson, $p_{\text{T}}(\text{Z})$, for events with at least one jet and breakdown of the relative uncertainty.

Bin-to-bin correlation in the measured cross section for Z+jets as a function of the transverse momementum of the Z boson, $p_{\text{T}}(\text{Z})$, for events with at least one jet.

Measured cross section for Z+jets as a function of the transverse momenum of the first jet, $p_{\text{T}}(\text{j}_1)$, and breakdown of the relative uncertainty.

Bin-to-bin correlation in the measured cross section for Z+jets as a function of the transverse momenum of the first jet, $p_{\text{T}}(\text{j}_1)$.

Measured cross section for Z+jets as a function of the transverse momenum of the second jet, $p_{\text{T}}(\text{j}_2)$, and breakdown of the relative uncertainty.

Bin-to-bin correlation in the measured cross section for Z+jets as a function of the transverse momenum of the second jet, $p_{\text{T}}(\text{j}_2)$.

Measured cross section for Z+jets as a function of the transverse momenum of the third jet, $p_{\text{T}}(\text{j}_3)$, and breakdown of the relative uncertainty.

Bin-to-bin correlation in the measured cross section for Z+jets as a function of the transverse momenum of the third jet, $p_{\text{T}}(\text{j}_3)$.

Measured cross section for Z+jets as a function of the absolute rapidity of the first jet, $|y(\text{j}_1)|$, and breakdown of the relative uncertainty.

Bin-to-bin correlation in the measured cross section for Z+jets as a function of the absolute rapidity of the first jet, $|y(\text{j}_1)|$.

Measured cross section for Z+jets as a function of the absolute rapidity of the second jet, $|y(\text{j}_2)|$, and breakdown of the relative uncertainty.

Bin-to-bin correlation in the measured cross section for Z+jets as a function of the absolute rapidity of the second jet, $|y(\text{j}_2)|$.

Measured cross section for Z+jets as a function of the absolute rapidity of the third jet, $|y(\text{j}_3)|$, and breakdown of the relative uncertainty.

Bin-to-bin correlation in the measured cross section for Z+jets as a function of the absolute rapidity of the third jet, $|y(\text{j}_3)|$.

Measured cross section for Z+jets as a function of the $H_{\text{T}}$ observable for event with at least one jet and breakdown of the relative uncertainty.

Bin-to-bin correlation in the measured cross section for Z+jets as a function of the $H_{\text{T}}$ observable for event with at least one jet.

Measured cross section for Z+jets as a function of the $H_{\text{T}}$ observable for event with at least two jets and breakdown of the relative uncertainty.

Bin-to-bin correlation in the measured cross section for Z+jets as a function of the $H_{\text{T}}$ observable for event with at least two jets.

Measured cross section for Z+jets as a function of the $H_{\text{T}}$ observable for event with at least three jets and breakdown of the relative uncertainty.

Bin-to-bin correlation in the measured cross section for Z+jets as a function of the $H_{\text{T}}$ observable for event with at least three jets.

Measured cross section for Z+jets as a function of the transverse momentum balance between the Z boson and the accompanying jets, $p_{\text{T}}^{\text{bal}}$, for events with at least one jet and breakdown of the relative uncertainty.

Bin-to-bin correlation in the measured cross section for Z+jets as a function of the transverse momentum balance between the Z boson and the accompanying jets, $p_{\text{T}}^{\text{bal}}$, for events with at least one jet.

Measured cross section for Z+jets as a function of the transverse momentum balance between the Z boson and the accompanying jets, $p_{\text{T}}^{\text{bal}}$, for events with at least two jets and breakdown of the relative uncertainty.

Bin-to-bin correlation in the measured cross section for Z+jets as a function of the transverse momentum balance between the Z boson and the accompanying jets, $p_{\text{T}}^{\text{bal}}$, for events with at least two jets.

Measured cross section for Z+jets as a function of the transverse momentum balance between the Z boson and the accompanying jets, $p_{\text{T}}^{\text{bal}}$, for events with at least three jets and breakdown of the relative uncertainty.

Bin-to-bin correlation in the measured cross section for Z+jets as a function of the transverse momentum balance between the Z boson and the accompanying jets, $p_{\text{T}}^{\text{bal}}$, for events with at least three jets.

Measured cross section for Z+jets as a function of the JZB variable with no restriction on $p_{\text{T}}(\text{Z})$ and breakdown of the relative uncertainty.

Bin-to-bin correlation in the measured cross section for Z+jets as a function of the JZB variable with no restriction on $p_{\text{T}}(\text{Z})$.

Measured cross section for Z+jets as a function of the JZB variable for $p_{\text{T}}(\text{Z})\leq50\,$GeV) and breakdown of the relative uncertainty.

Bin-to-bin correlation in the measured cross section for Z+jets as a function of the JZB variable for $p_{\text{T}}(\text{Z})\leq50\,$GeV).

Measured cross section for Z+jets as a function of the JZB variable for $p_{\text{T}}(\text{Z})>50\,$GeV and breakdown of the relative uncertainty.

Bin-to-bin correlation in the measured cross section for Z+jets as a function of the JZB variable for $p_{\text{T}}(\text{Z})>50\,$GeV.

Expected and observed upper limits at the 95% CL on the pp --> X --> ZZ cross section as a function of $m_X$ with $\Gamma_X$=10 GeV in VBF production mode.

Expected and observed upper limits at the 95% CL on the pp --> X --> ZZ cross section as a function of $m_X$ with $\Gamma_X$=100 GeV with VBF fraction profiled.

Expected and observed upper limits at the 95% CL on the pp --> X --> ZZ cross section as a function of $m_X$ with $\Gamma_X$=100 GeV in VBF production mode.

Observed upper limits at the 95% CL on the pp --> X --> ZZ cross section as a function of $m_X$ and $\Gamma_X$ values when VBF fraction is profiled.

Observed upper limits at the 95% CL on the pp --> X --> ZZ cross section as a function of $m_X$ and $\Gamma_X$ values in VBF production mode.

Distribution of $m_{jj}$ in data and simulation for events in the signal regions of the low-mass search. The points show the data while the filled histograms show the background contributions. The gray band shows the statistical and systematic uncertainties in the background, while the dashed vertical region ("Higgs") shows the expected SM Higgs boson mass range, which is excluded from this analysis. A 600 GeV bulk graviton signal prediction is represented by the black dashed histogram; for visibility, the signal cross-section is increased by a factor of 50. The background normalization is derived from the final fit to the $m_{VZ}$ observable in data.

The m$_{j}$ distributions of the events in data, compared to the expected background shape, for the high-mass analysis in the electron channel, for the high-purity category. The expected background shape is extracted from a fit to the data sidebands (Z+jets) or derived from simulation ("top quark" and "VV"). The dashed region around the background sum represents the uncertainty in the Z+jets distribution, while the dashed vertical region ("Higgs") shows the expected SM Higgs boson mass range, excluded from the analysis.

The m$_{j}$ distributions of the events in data, compared to the expected background shape, for the high-mass analysis in the electron channel, for the low-purity category. The expected background shape is extracted from a fit to the data sidebands (Z+jets) or derived from simulation ("top quark" and "VV"). The dashed region around the background sum represents the uncertainty in the Z+jets distribution, while the dashed vertical region ("Higgs") shows the expected SM Higgs boson mass range, excluded from the analysis.

The m$_{j}$ distributions of the events in data, compared to the expected background shape, for the high-mass analysis in the muon channel, for the high-purity category. The expected background shape is extracted from a fit to the data sidebands (Z+jets) or derived from simulation ("top quark" and "VV"). The dashed region around the background sum represents the uncertainty in the Z+jets distribution, while the dashed vertical region ("Higgs") shows the expected SM Higgs boson mass range, excluded from the analysis.

The m$_{j}$ distributions of the events in data, compared to the expected background shape, for the high-mass analysis in the muon channel, for the low-purity category. The expected background shape is extracted from a fit to the data sidebands (Z+jets) or derived from simulation ("top quark" and "VV"). The dashed region around the background sum represents the uncertainty in the Z+jets distribution, while the dashed vertical region ("Higgs") shows the expected SM Higgs boson mass range, excluded from the analysis.

Expected and observed distributions of the resonance candidate mass $m_{\text{VZ}}$ in the high-mass analysis, in the electron channel, high-purity category. The shaded area represents the post-fit uncertainty in the background. The expected contribution from W' signal candidates with mass $m_{\text{X}} = 2000$ GeV, normalized to a cross section of 100 fb$^{-1}$, is also shown.

Expected and observed distributions of the resonance candidate mass $m_{\text{VZ}}$ in the high-mass analysis, in the electron channel, low-purity category. The shaded area represents the post-fit uncertainty in the background. The expected contribution from W' signal candidates with mass $m_{\text{X}} = 2000$ GeV, normalized to a cross section of 100 fb$^{-1}$, is also shown.

Expected and observed distributions of the resonance candidate mass $m_{\text{VZ}}$ in the high-mass analysis, in the muon channel, high-purity category. The shaded area represents the post-fit uncertainty in the background. The expected contribution from W' signal candidates with mass $m_{\text{X}} = 2000$ GeV, normalized to a cross section of 100 fb$^{-1}$, is also shown.

Expected and observed distributions of the resonance candidate mass $m_{\text{VZ}}$ in the high-mass analysis, in the muon channel, low-purity category. The shaded area represents the post-fit uncertainty in the background. The expected contribution from W' signal candidates with mass $m_{\text{X}} = 2000$ GeV, normalized to a cross section of 100 fb$^{-1}$, is also shown.

Sideband $m_{ZV}$ distribution for the low-mass search in the merged V, untagged category, after fitting the sideband data alone. The points show the data while the filled histograms show the background contributions. Electron and muon categories are combined. The gray band indicates the statistical and post-fit systematic uncertainties in the normalization and shape of the background. Larger bin widths are used at higher values of $m_{ZV}$; the bin widths are indicated by the horizontal error bars.

Sideband $m_{ZV}$ distribution for the low-mass search in the merged V, tagged category, after fitting the sideband data alone. The points show the data while the filled histograms show the background contributions. Electron and muon categories are combined. The gray band indicates the statistical and post-fit systematic uncertainties in the normalization and shape of the background. Larger bin widths are used at higher values of $m_{ZV}$; the bin widths are indicated by the horizontal error bars.

Sideband $m_{ZV}$ distribution for the low-mass search in the resolved V, untagged category, after fitting the sideband data alone. The points show the data while the filled histograms show the background contributions. Electron and muon categories are combined. The gray band indicates the statistical and post-fit systematic uncertainties in the normalization and shape of the background. Larger bin widths are used at higher values of $m_{ZV}$; the bin widths are indicated by the horizontal error bars.

Sideband $m_{ZV}$ distribution for the low-mass search in the resolved V, tagged category, after fitting the sideband data alone. The points show the data while the filled histograms show the background contributions. Electron and muon categories are combined. The gray band indicates the statistical and post-fit systematic uncertainties in the normalization and shape of the background. Larger bin widths are used at higher values of $m_{ZV}$; the bin widths are indicated by the horizontal error bars.

The signal region $m_{VZ}$ distribution for the low-mass search, in the the merged V, untagged category, after fitting the signal and sideband regions. Electron and muon categories are combined. A 600 GeV bulk graviton signal prediction is represented by the black dashed histogram. The gray band indicates the statistical and post-fit systematic uncertainties in the normalization and shape of the background. Larger bin widths are used at higher values of $m_{VZ}$; the bin widths are indicated by the horizontal error bars.

The signal region $m_{VZ}$ distribution for the low-mass search, in the the merged V, tagged category, after fitting the signal and sideband regions. Electron and muon categories are combined. A 600 GeV bulk graviton signal prediction is represented by the black dashed histogram. The gray band indicates the statistical and post-fit systematic uncertainties in the normalization and shape of the background. Larger bin widths are used at higher values of $m_{VZ}$; the bin widths are indicated by the horizontal error bars.

The signal region $m_{VZ}$ distribution for the low-mass search, in the the resolved V, untagged category, after fitting the signal and sideband regions. Electron and muon categories are combined. A 600 GeV bulk graviton signal prediction is represented by the black dashed histogram. The gray band indicates the statistical and post-fit systematic uncertainties in the normalization and shape of the background. Larger bin widths are used at higher values of $m_{VZ}$; the bin widths are indicated by the horizontal error bars.

The signal region $m_{VZ}$ distribution for the low-mass search, in the the resolved V, tagged category, after fitting the signal and sideband regions. Electron and muon categories are combined. A 600 GeV bulk graviton signal prediction is represented by the black dashed histogram. The gray band indicates the statistical and post-fit systematic uncertainties in the normalization and shape of the background. Larger bin widths are used at higher values of $m_{VZ}$; the bin widths are indicated by the horizontal error bars.

Observed and expected 95% CL upper limit on $\sigma_{W'} Br(W' \to ZW)$ as a function of the resonance mass, taking into account all statistical and systematic uncertainties. The electron and muon channels and the various categories used in the analysis are combined together. The green (inner) and yellow (outer) bands represent the 68% and 95% coverage of the expected limit in the background-only hypothesis.

Observed and expected 95% CL upper limit on $\sigma_{G} Br(G \to ZZ)$ as a function of the resonance mass, taking into account all statistical and systematic uncertainties. The electron and muon channels and the various categories used in the analysis are combined together. The green (inner) and yellow (outer) bands represent the 68% and 95% coverage of the expected limit in the background-only hypothesis.

Observed local p-values for W' narrow resonances as a function of the resonance mass.

Observed local p-values for G narrow resonances as a function of the resonance mass.

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.

Reconstructed mtb distribution in data and for the background after the fit to the data in the signal region SR2. The statistical uncertainty on data points is calculated using assymetric Poisson confidence intervals.

Reconstructed mtb distribution in data and for the background after the fit to the data in the signal region SR3. The statistical uncertainty on data points is calculated using assymetric Poisson confidence intervals.

Reconstructed mtb distribution in data and for the background after the fit to the data in the validation region VR. The statistical uncertainty on data points is calculated using assymetric Poisson confidence intervals.