Observed likelihood scan in the ($\kappa_\tau$, $\tilde{\kappa}_\tau$) plane.

Observed likelihood scan in the ($\alpha^{\mathrm{H}\tau\tau}$, $\mu_{\text{ggH}}$) plane. The $\mu_{\text{qqH}}$ signal strength modifier is profiled.

Observed likelihood scan in the ($\alpha^{\mathrm{H}\tau\tau}$, $\mu_{\text{qqH}}$) plane. The $\mu_{\text{ggH}}$ signal strength modifier is profiled.

Observed likelihood scan in the ($\mu_{\text{qqH}}$, $\mu_{\text{ggH}}$) plane. The value of $\alpha^{\mathrm{H}\tau\tau}$ is profiled.

Observed profiled likelihood on $f_{\Lambda1}$ using Approach 1. The signal strength modifiers are treated as free parameters. Axis scales are varied to improve the visibility of important features.

Expected profiled likelihood on $f_{a3}$ using Approach 1. The signal strength modifiers are treated as free parameters.

Observed profiled likelihood on $f_{a3}$ using Approach 1. The signal strength modifiers are treated as free parameters.

Expected profiled likelihood on $f_{\Lambda1}^{Z\gamma}$ using Approach 1. The signal strength modifiers are treated as free parameters.

Observed profiled likelihood on $f_{\Lambda1}^{Z\gamma}$ using Approach 1. The signal strength modifiers are treated as free parameters.

Expected profiled likelihood on $f_{a2}$ using Approach 2. The other two anomalous coupling cross section fractions are fixed to zero. The signal strength modifiers are treated as free parameters. Axis scales are varied to improve the visibility of important features.

Observed profiled likelihood on $f_{a2}$ using Approach 2. The other two anomalous coupling cross section fractions are fixed to zero. The signal strength modifiers are treated as free parameters. Axis scales are varied to improve the visibility of important features.

Expected profiled likelihood on $f_{\Lambda1}$ using Approach 2. The other two anomalous coupling cross section fractions are fixed to zero. The signal strength modifiers are treated as free parameters. Axis scales are varied to improve the visibility of important features.

Observed profiled likelihood on $f_{\Lambda1}$ using Approach 2. The other two anomalous coupling cross section fractions are fixed to zero. The signal strength modifiers are treated as free parameters. Axis scales are varied to improve the visibility of important features.

Expected profiled likelihood on $f_{a3}$ using Approach 2. The other two anomalous coupling cross section fractions are fixed to zero. The signal strength modifiers are treated as free parameters.

Observed profiled likelihood on $f_{a3}$ using Approach 2. The other two anomalous coupling cross section fractions are fixed to zero. The signal strength modifiers are treated as free parameters.

Expected profiled likelihood on $f_{a2}$ using Approach 2. The other two anomalous coupling cross section fractions are left floating in the fit. The signal strength modifiers are treated as free parameters. Axis scales are varied to improve the visibility of important features.

Observed profiled likelihood on $f_{a2}$ using Approach 2. The other two anomalous coupling cross section fractions are left floating in the fit. The signal strength modifiers are treated as free parameters. Axis scales are varied to improve the visibility of important features.

Expected profiled likelihood on $f_{\Lambda1}$ using Approach 2. The other two anomalous coupling cross section fractions are left floating in the fit. The signal strength modifiers are treated as free parameters. Axis scales are varied to improve the visibility of important features.

Observed profiled likelihood on $f_{\Lambda1}$ using Approach 2. The other two anomalous coupling cross section fractions are left floating in the fit. The signal strength modifiers are treated as free parameters. Axis scales are varied to improve the visibility of important features.

Expected profiled likelihood on $f_{a3}$ using Approach 2. The other two anomalous coupling cross section fractions are left floating in the fit. The signal strength modifiers are treated as free parameters.

Observed profiled likelihood on $f_{a3}$ using Approach 2. The other two anomalous coupling cross section fractions are left floating in the fit. The signal strength modifiers are treated as free parameters.

Expected profiled likelihood on the $\delta c_{Z}$ coupling of the SMEFT Higgs basis.

Observed profiled likelihood on the $\delta c_{Z}$ coupling of the SMEFT Higgs basis.

Expected profiled likelihood on the $c_{Z\Box}$ coupling of the SMEFT Higgs basis.

Observed profiled likelihood on the $c_{Z\Box}$ coupling of the SMEFT Higgs basis.

Expected profiled likelihood on the $c_{ZZ}$ coupling of the SMEFT Higgs basis.

Observed profiled likelihood on the $c_{ZZ}$ coupling of the SMEFT Higgs basis.

Expected profiled likelihood on the $\tilde{c}_{ZZ}$ coupling of the SMEFT Higgs basis.

Observed profiled likelihood on the $\tilde{c}_{ZZ}$ coupling of the SMEFT Higgs basis.

Expected profiled likelihood on $f_{a3}^{ggH}$. The signal strength modifiers and the CP-odd HVV anomalous coupling cross section fraction are treated as free parameters.

Observed profiled likelihood on $f_{a3}^{ggH}$. The signal strength modifiers and the CP-odd HVV anomalous coupling cross section fraction are treated as free parameters.

Expected likelihood scan of $f_{a3}^{\mathrm{ggH}}$ ($\Delta\phi_{\mathrm{jj}}$ method).

Observed likelihood scan of $\alpha_{\mathrm{Hff}}$ (MELA method).

Expected likelihood scan of $\alpha_{\mathrm{Hff}}$ (MELA method).

Observed likelihood scan of $\alpha_{\mathrm{Hff}}$ ($\Delta\phi_{\mathrm{jj}}$ method).

Expected likelihood scan of $\alpha_{\mathrm{Hff}}$ ($\Delta\phi_{\mathrm{jj}}$ method).

Observed likelihood scan of $f_{a3}$ (Approach 1).

Expected likelihood scan of $f_{a3}$ (Approach 1).

Observed likelihood scan of $f_{a2}$.

Expected likelihood scan of $f_{a2}$.

Observed likelihood scan of $f_{\Lambda 1}$.

Expected likelihood scan of $f_{\Lambda 1}$.

Observed likelihood scan of $f_{\Lambda 1}^{\mathrm{Z}\gamma}$.

Expected likelihood scan of $f_{\Lambda 1}^{\mathrm{Z}\gamma}$.

Observed likelihood scan of $f_{a3}$ (Approach 2).

Expected likelihood scan of $f_{a3}$ (Approach 2).

Observed likelihood scan of $f_{a3}$ (Approach 1). $\mathrm{H}\to\tau\tau$ and $\mathrm{H}\to4\ell$ channel results are combined.

Expected likelihood scan of $f_{a3}$ (Approach 1). $\mathrm{H}\to\tau\tau$ and $\mathrm{H}\to4\ell$ channel results are combined.

Observed likelihood scan of $f_{a2}$. $\mathrm{H}\to\tau\tau$ and $\mathrm{H}\to4\ell$ channel results are combined.

Expected likelihood scan of $f_{a2}$. $\mathrm{H}\to\tau\tau$ and $\mathrm{H}\to4\ell$ channel results are combined.

Observed likelihood scan of $f_{\Lambda 1}$. $\mathrm{H}\to\tau\tau$ and $\mathrm{H}\to4\ell$ channel results are combined.

Expected likelihood scan of $f_{\Lambda 1}$. $\mathrm{H}\to\tau\tau$ and $\mathrm{H}\to4\ell$ channel results are combined.

Observed likelihood scan of $f_{\Lambda 1}^{\mathrm{Z}\gamma}$. $\mathrm{H}\to\tau\tau$ and $\mathrm{H}\to4\ell$ channel results are combined.

Expected likelihood scan of $f_{\Lambda 1}^{\mathrm{Z}\gamma}$. $\mathrm{H}\to\tau\tau$ and $\mathrm{H}\to4\ell$ channel results are combined.

Observed likelihood scan of $f_{a3}$ (Approach 2). $\mathrm{H}\to\tau\tau$ and $\mathrm{H}\to4\ell$ channel results are combined.

Expected likelihood scan of $f_{a3}$ (Approach 2). $\mathrm{H}\to\tau\tau$ and $\mathrm{H}\to4\ell$ channel results are combined.

Observed likelihood scan of $f_{a3}^{\mathrm{ggH}}$. $\mathrm{H}\to\tau\tau$ and $\mathrm{H}\to4\ell$ channel results are combined.

Expected likelihood scan of $f_{a3}^{\mathrm{ggH}}$. $\mathrm{H}\to\tau\tau$ and $\mathrm{H}\to4\ell$ channel results are combined.

Observed likelihood scan of $f_{\textit{CP}}^{\mathrm{Htt}}$. $\mathrm{H}\to\tau\tau$, $\mathrm{H}\to4\ell$, and $\mathrm{H}\to\gamma\gamma$ channel results are combined.

Expected likelihood scan of $f_{\textit{CP}}^{\mathrm{Htt}}$. $\mathrm{H}\to\tau\tau$, $\mathrm{H}\to4\ell$, and $\mathrm{H}\to\gamma\gamma$ channel results are combined.

Observed likelihood scan of $c_{\mathrm{gg}}$ ($\kappa$ and $\widetilde{\kappa}_\mathrm{t}$ profiled). $\mathrm{H}\to\tau\tau$, $\mathrm{H}\to4\ell$, and $\mathrm{H}\to\gamma\gamma$ channel results are combined.

Expected likelihood scan of $c_{\mathrm{gg}}$ ($\kappa$ and $\widetilde{\kappa}_\mathrm{t}$ profiled). $\mathrm{H}\to\tau\tau$, $\mathrm{H}\to4\ell$, and $\mathrm{H}\to\gamma\gamma$ channel results are combined.

Observed likelihood scan of $\widetilde{c}_{\mathrm{gg}}$ ($\kappa$ and $\widetilde{\kappa}_\mathrm{t}$ profiled)

Expected likelihood scan of $\widetilde{c}_{\mathrm{gg}}$ ($\kappa$ and $\widetilde{\kappa}_\mathrm{t}$ profiled). $\mathrm{H}\to\tau\tau$, $\mathrm{H}\to4\ell$, and $\mathrm{H}\to\gamma\gamma$ channel results are combined.

Observed likelihood scan of $c_{\mathrm{gg}}$ ($\kappa = 1$ and $\widetilde{\kappa}_\mathrm{t} = 0$). $\mathrm{H}\to\tau\tau$, $\mathrm{H}\to4\ell$, and $\mathrm{H}\to\gamma\gamma$ channel results are combined.

Expected likelihood scan of $c_{\mathrm{gg}}$ ($\kappa = 1$ and $\widetilde{\kappa}_\mathrm{t} = 0$). $\mathrm{H}\to\tau\tau$, $\mathrm{H}\to4\ell$, and $\mathrm{H}\to\gamma\gamma$ channel results are combined.

Observed likelihood scan of $\widetilde{c}_{\mathrm{gg}}$ ($\kappa = 1$ and $\widetilde{\kappa}_\mathrm{t} = 0$). $\mathrm{H}\to\tau\tau$, $\mathrm{H}\to4\ell$, and $\mathrm{H}\to\gamma\gamma$ channel results are combined.

Expected likelihood scan of $\widetilde{c}_{\mathrm{gg}}$ ($\kappa = 1$ and $\widetilde{\kappa}_\mathrm{t} = 0$). $\mathrm{H}\to\tau\tau$, $\mathrm{H}\to4\ell$, and $\mathrm{H}\to\gamma\gamma$ channel results are combined.

Higgs boson transverse momentum spectrum, $H\rightarrow b\overline{b}$

Combined Higgs boson transverse momentum spectrum, ggH only

Combined Higgs boson jet multiplicity spectrum

Higgs boson jet multiplicity spectrum, $H\rightarrow\gamma\gamma$

Higgs boson jet multiplicity spectrum, $H\rightarrow ZZ$

Combined Higgs boson rapidity spectrum

Higgs boson rapidity spectrum, $H\rightarrow\gamma\gamma$

Higgs boson rapidity spectrum, $H\rightarrow ZZ$

Combined leading jet transverse momentum spectrum

Leading jet transverse momentum spectrum, $H\rightarrow\gamma\gamma$

Leading jet transverse momentum spectrum, $H\rightarrow ZZ$

The correlation matrix for the Higgs pT measurements, both for the unregularized and regularized fits. The last bin is inclusive.

The correlation matrix for the jet multiplicity measurements, both for the unregularized and regularized fits. The last bin is inclusive.

The correlation matrix for the leading jet pT measurements, both for the unregularized and regularized fits. The last bin is inclusive.

The fiducial integrated signal strength and cross section.

The correlation matrix for the njet measurement. The last bin is inclusive.

The fiducial integrated signal strength and cross section The uncertainty breakdown is given in terms of statistical (stat), experimental (exp), theoretical uncertainties on the background (bkg) and on the signal (sig), and the luminosity uncertainty (lumi). The fiducial cross section and its full uncertainty in each bin are also given.

Distribution of an input to the BDT classifier in the $3\ell$ category: The linear discrimminant $p_{T}^{miss,LD}$. The discriminant is defined by the relation $p_{T}^{miss,LD} = 0.6 p_{T}^{miss} + 0.4 H_{T}^{miss}$, where $H_{T}^{miss}$ corresponds to the magnitude of the vector $p_{T}$ sum of all electrons, muons, taus and small radius jets passing the selection critertia, and $p_{T}^{miss}$ represents the missing transverse momentum vector computed as the negative vector $p_{T}$ sum of all the particles reconstructed by the PF algorithm in the event.

Distributions of the transverse mass $m_{T}$ of the lepton that is not originating from the Z boson decay and the missing transverse momentum in the $3\ell$/WZ control region. The distributions expected for the WZ and ZZ as well as for other background processes are shown for the values of nuisance parameters obtained from the ML fit used in the signal extraction.

Distributions of the four lepton invariant mass $m_{4\ell}$ in the $4\ell$/ZZ control region. The distributions expected for the ZZ as well as for other background processes are shown for the values of nuisance parameters obtained from the ML fit used in the signal extraction.

Distributions of the transverse mass $m_{T}$ of the leading lepton and the missing transverse momentum in the $2\ell$(ss) control region. The distributions expected for the misidentified $\ell$ background as well as for other background processes is shown for the values of nuisance parameters obtained from the background-only ML fit, in which the HH signal is constrained to be zero.

Distributions of the reconstructed HH mass $m_{HH}$ in the $2\ell+2 au_{h}$ control region. The distributions expected for the misidentified $\ell/ au_{h}$ background as well as for other background processes is shown for the values of nuisance parameters obtained from the background-only ML fit, in which the HH signal is constrained to be zero.

Distribution in the output of the BDT trained for nonresonant HH production and evaluated for the benchmark scenario JHEP04 BM7 for the $2\ell$(ss) category. The SM HH signal is shown for a cross section amounting to $30$ times the value predicted in the SM. The distributions expected for the background processes are shown for the values of nuisance parameters obtained from the ML fit of the signal+background hypothesis to the data. The binning is chosen to approximate signal quantiles.

Distribution in the output of the BDT trained for nonresonant HH production and evaluated for the benchmark scenario JHEP04 BM7 for the $3\ell$ category. The SM HH signal is shown for a cross section amounting to $30$ times the value predicted in the SM. The distributions expected for the background processes are shown for the values of nuisance parameters obtained from the ML fit of the signal+background hypothesis to the data. The binning is chosen to approximate signal quantiles.

Distribution in the output of the BDT trained for nonresonant HH production and evaluated for the benchmark scenario JHEP04 BM7 for the $4\ell$ category. The SM HH signal is shown for a cross section amounting to $30$ times the value predicted in the SM. The distributions expected for the background processes are shown for the values of nuisance parameters obtained from the ML fit of the signal+background hypothesis to the data. The binning is chosen to approximate signal quantiles.

Distribution in the output of the BDT trained for nonresonant HH production and evaluated for the benchmark scenario JHEP04 BM7 for the $3\ell+1 au_{h}$ category. The SM HH signal is shown for a cross section amounting to $30$ times the value predicted in the SM. The distributions expected for the background processes are shown for the values of nuisance parameters obtained from the ML fit of the signal+background hypothesis to the data. The binning is chosen to approximate background quantiles.

Distribution in the output of the BDT trained for nonresonant HH production and evaluated for the benchmark scenario JHEP04 BM7 for the $2\ell+2 au_{h}$ category. The SM HH signal is shown for a cross section amounting to $30$ times the value predicted in the SM. The distributions expected for the background processes are shown for the values of nuisance parameters obtained from the ML fit of the signal+background hypothesis to the data. The binning is chosen to approximate background quantiles.

Distribution in the output of the BDT trained for nonresonant HH production and evaluated for the benchmark scenario JHEP04 BM7 for the $1\ell+3 au_{h}$ category. The SM HH signal is shown for a cross section amounting to $30$ times the value predicted in the SM. The distributions expected for the background processes are shown for the values of nuisance parameters obtained from the ML fit of the signal+background hypothesis to the data. The binning is chosen to approximate background quantiles.

Distribution in the output of the BDT trained for nonresonant HH production and evaluated for the benchmark scenario JHEP04 BM7 for the $4 au_{h}$ category. The SM HH signal is shown for a cross section amounting to $30$ times the value predicted in the SM. The distributions expected for the background processes are shown for the values of nuisance parameters obtained from the ML fit of the signal+background hypothesis to the data. The binning is chosen to approximate background quantiles.

Observed and expected $95\%$ CL upper limits on the SM HH production cross section, obtained for both individual search categories and from a simultaneous fit of all seven categories combined.

Observed and expected $95\%$ CL upper limits on the HH production cross section as a function of the Higgs boson self-coupling strength modifier $\kappa_\lambda$ for the combination of all seven categories. All Higgs boson couplings other than $\lambda$ are assumed to have the values predicted in the SM.

Observed and expected $95\%$ CL upper limits on the HH production cross section as a function of the Higgs boson self-coupling strength modifier $\kappa_\lambda$ for the seven different categories as well as their combination. All Higgs boson couplings other than $\lambda$ are assumed to have the values predicted in the SM.

Observed and expected $95\%$ CL upper limits on the HH production cross section for the twelve benchmark scenarios from doi:10.1007/JHEP04(2016)126, the additional benchmark scenario 8a from doi:10.1007/JHEP09(2018)057, the seven benchmark scenarios from doi:10.1007/JHEP03(2020)091, and for the SM.

Observed and expected $95\%$ CL upper limits on the HH production cross section for the twelve benchmark scenarios from doi:10.1007/JHEP04(2016)126, the additional benchmark scenario 8a from doi:10.1007/JHEP09(2018)057, the seven benchmark scenarios from doi:10.1007/JHEP03(2020)091, and for the SM. The limits are shown for all seven search categories as well as for their combination.

Observed and expected $95\%$ CL upper limits on the HH production cross section as a function of the effective coupling $c_{2}$ for the combination of all seven categories. All Higgs boson couplings other than $c_{2}$ are assumed to have the values predicted in the SM.

Observed and expected $95\%$ CL upper limits on the HH production cross section as a function of the effective coupling $c_{2}$ and the top Yukawa coupling modifier $\kappa_{t}$ for the combination of all seven categories. All Higgs boson couplings other than $c_{2}$ and $\kappa_{t}$ are assumed to have the values predicted in the SM. The position of the SM in the ($c_{2}-\kappa_{t}$) plane, as well as the best fit value of $(c_{2},\kappa_{t})=(1.05, 1.74)$ together with contours for the theory cross section are shown as well.

Observed and expected $95\%$ CL upper limits on the HH production cross section as a function of the Higgs boson self-coupling modifier $\kappa_{\lambda}$ and the top Yukawa coupling modifier $\kappa_{t}$ for the combination of all seven categories. All Higgs boson couplings other than $\kappa_{\lambda}$ and $\kappa_{t}$ are assumed to have the values predicted in the SM. The position of the SM in the ($\kappa_{\lambda}-\kappa_{t}$) plane, as well as the best fit value of $(\kappa_{\lambda},\kappa_{t})=(-3.54, -1.70)$ together with contours for the theory cross section are shown as well.

Observed and expected $95\%$ CL upper limits on the HH production cross section as a function of the effective coupling $c_{2}$ and the Higgs boson self-coupling modifier $\kappa_{\lambda}$ for the combination of all seven categories. All Higgs boson couplings other than $c_{2}$ and $\kappa_{\lambda}$ are assumed to have the values predicted in the SM. The position of the SM in the ($c_{2}-\kappa_{\lambda}$) plane, as well as the best fit value of $(c_{2},\kappa_{\lambda})=(-0.66, 5.33)$ together with contours for the theory cross section are shown as well.

Observed and expected $95\%$ CL upper limits on the production of new particles X of spin $0$ and mass $m_{X}$ in the range $250$-$1000$ GeV, which decay to Higgs boson pairs.

Observed and expected $95\%$ CL upper limits on the production of new particles X of spin $0$ and mass $m_{X}$ in the range $250$-$1000$ GeV, which decay to Higgs boson pairs. The Limit is shown for the seven different search categories as well as their combination.

Observed and expected $95\%$ CL upper limits on the production of new particles X of spin $2$ and mass $m_{X}$ in the range $250$-$1000$ GeV, which decay to Higgs boson pairs.

Observed and expected $95\%$ CL upper limits on the production of new particles X of spin $2$ and mass $m_{X}$ in the range $250$-$1000$ GeV, which decay to Higgs boson pairs. The Limit is shown for the seven different search categories as well as their combination.

Distribution in the output of the BDT trained for resonances of spin 2 and mass $750\,$GeV production for the $2\ell$(ss) category. The resonant HH signal is shown for a cross section amounting to $1\,$pb. The distributions expected for the background processes are shown for the values of nuisance parameters obtained from the ML fit of the signal+background hypothesis to the data. The binning is chosen to approximate signal quantiles.

Distribution in the output of the BDT trained for resonances of spin 2 and mass $750\,$GeV production for the $3\ell$ category. The resonant HH signal is shown for a cross section amounting to $1\,$pb. The distributions expected for the background processes are shown for the values of nuisance parameters obtained from the ML fit of the signal+background hypothesis to the data. The binning is chosen to approximate signal quantiles.

The product of signal acceptance and efficiency in the 2l categories for the signal produced via vector boson fusion.

$m_{X}^{T}$ distribution in data in the 0l 2b non-VBF category. The data is shown as Events/10 GeV. The distribution is shown up 4000 GeV, which corresponds to the event with the highest $m_{X}^{T}$ observed in the SR.

$m_{X}^{T}$ distribution in data in the 0l $\leq$1b non-VBF category. The data is shown as Events/10 GeV. The distribution is shown up 4000 GeV, which corresponds to the event with the highest $m_{X}^{T}$ observed in the SR.

$m_{X}$ distribution in data in the 2e 2b non-VBF category. The data is shown as Events/10 GeV. The distribution is shown up 4000 GeV, which corresponds to the event with the highest $m_{X}$ observed in the SR.

$m_{X}$ distribution in data in the 2e $\leq$1b non-VBF category. The data is shown as Events/10 GeV. The distribution is shown up 4000 GeV, which corresponds to the event with the highest $m_{X}$ observed in the SR.

$m_{X}$ distribution in data in the 2$\mu$ 2b non-VBF category. The data is shown as Events/10 GeV. The distribution is shown up 4000 GeV, which corresponds to the event with the highest $m_{X}$ observed in the SR.

$m_{X}$ distribution in data in the 2$\mu$ $\leq$1b non-VBF category. The data is shown as Events/10 GeV. The distribution is shown up 4000 GeV, which corresponds to the event with the highest $m_{X}$ observed in the SR.

$m_{X}^{T}$ distribution in data in the 0l 2b VBF category. The data is shown as Events/10 GeV. The distribution is shown up 4000 GeV, which corresponds to the event with the highest $m_{X}^{T}$ observed in the SR.

$m_{X}^{T}$ distribution in data in the 0l $\leq$1b VBF category. The data is shown as Events/10 GeV. The distribution is shown up 4000 GeV, which corresponds to the event with the highest $m_{X}^{T}$ observed in the SR.

$m_{X}$ distribution in data in the 2e 2b VBF category. The data is shown as Events/10 GeV. The distribution is shown up 4000 GeV, which corresponds to the event with the highest $m_{X}$ observed in the SR.

$m_{X}$ distribution in data in the 2e $\leq$1b VBF category. The data is shown as Events/10 GeV. The distribution is shown up 4000 GeV, which corresponds to the event with the highest $m_{X}$ observed in the SR.

$m_{X}$ distribution in data in the 2$\mu$ 2b VBF category. The data is shown as Events/10 GeV. The distribution is shown up 4000 GeV, which corresponds to the event with the highest $m_{X}$ observed in the SR.

$m_{X}$ distribution in data in the 2$\mu$ $\leq$1b VBF category. The data is shown as Events/10 GeV. The distribution is shown up 4000 GeV, which corresponds to the event with the highest $m_{X}$ observed in the SR.

Observed and expected 95% CL upper limit on $\sigma \mathcal{B}$(Z'-> ZH) with all categories combined for the non-VBF signal, including all statistical and systematic uncertainties. The inner green band and the outer yellow band indicate the regions containing 68 and 95%, respectively, of the distribution of expected limits under the background-only hypothesis. The CMS search for a heavy resonance using 2016 data and the same final state [JHEP 11 (2018) 172] is shown as a comparison.

Observed and expected 95% CL upper limit on $\sigma \mathcal{B}$(Z'-> ZH) with all categories combined for the VBF signal, including all statistical and systematic uncertainties. The inner green band and the outer yellow band indicate the regions containing 68 and 95%, respectively, of the distribution of expected limits under the background-only hypothesis.

Observed exclusion limit in the space of the HVT model parameters [$g_{V}c_{H}$, $g^{2}c_{F}/g_{V}$] for mass hypotheses of 2 TeV for the non-VBF signal.

Observed exclusion limit in the space of the HVT model parameters [$g_{V}c_{H}$, $g^{2}c_{F}/g_{V}$] for mass hypotheses of 3 TeV for the non-VBF signal.

Observed exclusion limit in the space of the HVT model parameters [$g_{V}c_{H}$, $g^{2}c_{F}/g_{V}$] for mass hypotheses of 4 TeV for the non-VBF signal.