Observed and expected 95% CL upper limits on the cross section in the ZH channel.

Signal acceptance times efficiency of HVT WH(qqbb) events as a function of the resonance mass at different cut stages. Auxiliary table attached for 2 TeV mass point.

Signal acceptance times efficiency of HVT WH(qqbb) events as a function of the resonance mass at different cut stages. Auxiliary table attached for 2 TeV mass point.

Signal acceptance times efficiency of HVT ZH(qqbb) events as a function of the resonance mass at different cut stages. Auxiliary table attached for 4 TeV mass point.

Signal acceptance times efficiency of HVT ZH(qqbb) events as a function of the resonance mass at different cut stages. Auxiliary table attached for 4 TeV mass point.

Dijet mass distributions in the WH 1-tag signal region. Uncertainties in the background and signal histograms are set to zero. The signal corresponds to that expected for HVT model B with resonance mass of 2 TeV.

Dijet mass distributions in the WH 2-tag signal region. Uncertainties in the background and signal histograms are set to zero. The signal corresponds to that expected for HVT model B with resonance mass of 2 TeV.

Dijet mass distributions in the ZH 1-tag signal region. Uncertainties in the background and signal histograms are set to zero. The signal corresponds to that expected for HVT model B with resonance mass of 2 TeV.

Dijet mass distributions in the ZH 2-tag signal region. Uncertainties in the background and signal histograms are set to zero. The signal corresponds to that expected for HVT model B with resonance mass of 2 TeV.

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

Distribution of $m_{\mathrm{jjj}}$ for preselected events with $\mathrm{N}_{j}$ = 3

Distribution of $m_{\mathrm{j}}$ for preselected events with $\mathrm{N}_{j}$ = 3

Distribution of the deep-WH value of the highest-mass jet with 60 < $m_{\mathrm{j}}$ < 100 GeV for preselected events with $\mathrm{N}_{j}$ = 3

scale factors (SFs) for W, $t^{2}$, and q/g matched jets in the low-$m_{\mathrm{j}}$ and low-$p_{\mathrm{T}}$ (LL) bin, as functions of the deep-W discriminant value.

scale factors (SFs) for W, $t^{2}$, and q/g matched jets in the low-$m_{\mathrm{j}}$ and high-$p_{\mathrm{T}}$ (LH) bin, as functions of the deep-W discriminant value.

scale factors (SFs) for $t^{2}$, $t^{3,4}$, and q/g matched jets in the high-$m_{\mathrm{j}}$ and low-$p_{\mathrm{T}}$ (HL) bin, as functions of the deep-WH discriminant value.

scale factors (SFs) for $t^{2}$, $t^{3,4}$, and q/g matched jets in the high-$m_{\mathrm{j}}$ and high-$p_{\mathrm{T}}$ (HH) bin, as functions of the deep-WH discriminant value.

The deep-W discriminant of the jet with highest mass in the single-lepton sideband for LL samples.

The deep-W discriminant of the jet with highest mass in the single-lepton sideband for LH samples.

The deep-WH discriminant of the jet with highest mass in the single-lepton sideband for HL samples.

The deep-WH discriminant of the jet with highest mass in the single-lepton sideband for HH samples.

Comparison of the distribution of data and simulated backgrounds, as a function of the deep-W discriminant value for the highest-mass jet in CR1, after the scale factors have been applied.

Comparison of the distribution of data and simulated backgrounds, as a function of the deep-WH discriminant value for the highest-mass jet in CR2, after the scale factors have been applied.

Comparison of the distribution of data and simulated backgrounds, as a function of the deep-WH discriminant value for the highest-mass jet in CR3, after the scale factors have been applied.

Comparison of the distribution of data and simulated backgrounds, as a function of the deep-W discriminant value for the highest-mass jet in CR45, after the scale factors have been applied.

Comparison of the distribution of data and simulated backgrounds, as a function of the deep-W discriminant value for the highest-mass jet in CR6, after the scale factors have been applied.

The $m_{\mathrm{j}^{\mathrm{max}}}$ distributions for different jet types for SR1–3 events of the signal with $m_{\mathrm{W}_{\mathrm{KK}}}$ = 2.5 TeV, $m_{\mathrm{R}}$ = 0.2 TeV$ without deep-W (WH) constraints.

The deep-W distribution normalized to unity for the shown components of of the signal with $m_{\mathrm{W}_{\mathrm{KK}}}$ = 2.5 TeV, $m_{\mathrm{R}}$ = 0.2 TeV. The $t^{3,4}$ jets from the preselected sample, normalized to unity, are superimposed to compare shapes with the $R^{3q}$ and $R^{4q}$ distributions.

The deep-WH distribution normalized to unity for the shown components of of the signal with $m_{\mathrm{W}_{\mathrm{KK}}}$ = 2.5 TeV, $m_{\mathrm{R}}$ = 0.2 TeV. The $t^{3,4}$ jets from the preselected sample, normalized to unity, are superimposed to compare shapes with the $R^{3q}$ and $R^{4q}$ distributions.

The $m_{\mathrm{jj}}$ distribution for CR1 for data and simulation.

The $m_{\mathrm{jj}}$ distribution for CR2 for data and simulation.

The $m_{\mathrm{jj}}$ distribution for CR3 for data and simulation.

The $m_{\mathrm{jjj}}$ distribution for CR45 for data and simulation.

The $m_{\mathrm{jjj}}$ distribution for CR6 for data and simulation.

Post-fit distributions of the reconstructed triboson system ($m_{\mathrm{jj}}$) in data and simulation for SR1.

Post-fit distributions of the reconstructed triboson system ($m_{\mathrm{jj}}$) in data and simulation for SR2.

Post-fit distributions of the reconstructed triboson system ($m_{\mathrm{jj}}$) in data and simulation for SR3.

Post-fit distributions of the reconstructed triboson system ($m_{\mathrm{jjj}}$) in data and simulation for SR4.

Post-fit distributions of the reconstructed triboson system ($m_{\mathrm{jjj}}$) in data and simulation for SR5.

Post-fit distributions of the reconstructed triboson system ($m_{\mathrm{jjj}}$) in data and simulation for SR6.

Observed upper limits at $95\%$ CL on the signal cross section $\times$ branching fraction as functions of the $m_{\mathrm{W}_{\mathrm{KK}}}$ and $m_{\mathrm{R}}$ resonance masses.

Expected median lower limit contour on the $m_{\mathrm{W}_{\mathrm{KK}}}$ and $m_{\mathrm{R}}$ plane.

Expected $+ 1$ s.d. lower limit contour on the $m_{\mathrm{W}_{\mathrm{KK}}}$ and $m_{\mathrm{R}}$ plane.

Expected - 1 s.d. lower limit contour on the $m_{\mathrm{W}_{\mathrm{KK}}}$ and $m_{\mathrm{R}}$ plane.

Observed lower limit contour on the $m_{\mathrm{W}_{\mathrm{KK}}}$ and $m_{\mathrm{R}}$ plane.

The signal cross-sections values of the signal QBH RS1 with n=1 extra dimensions, as a function of threshold mass

The product of acceptance and selection efficiency for QBH ADD model, as a function of threshold mass of the signal.

The product of acceptance and selection efficiency for QBH RS1 models, as a function of threshold mass of the signal.

The product of acceptance and selection efficiency for q* signal, as a function of resonance mass

The product of acceptance and selection efficiency for b* signal, as a function of resonance mass of the signal.

Distributions of the γ+jet invariant mass generated from various signal samples for reconstructed mass 4 TeV

Fitting distribution of data with expected q*, and QBH signal events, and fit parameters p0= 3.05674e+00, p1 = 10.5343, p3 = 4.15163 and p4 = 0.107405

Fitting distribution of data with expected b* events, and fit parameters p0= 1.45055e+03, p1 = 2.04830e+01, p3 = -4.16760e-01 and p4 = -5.07768e-01

Fitting distribution of data with expected b* events for 0btag category, and fit parameters p0= 2.44120e+00, p1 = 1.08886e+01, p3 = 4.23198e+00 and p4 = 1.16300e-01

Exclusion limits on the product of cross section and branching fraction for q* signal coupling f = 1.0, as a function of the resonance mass hypothesis.

Theoretical cross-section of q* signal for coupling multiplier f =1.0

Exclusion limits on the product of cross section and branching fraction for the q* signal coupling f = 0.5, as a function of the resonance mass hypothesis.

Theoretical cross-section of q* signal for coupling multiplier f =0.5

Exclusion limits on the product of cross section and branching fraction for q* coupling f = 0.1, as a function of the resonance mass hypothesis.

Theoretical cross-section of q* signal for coupling multiplier f =0.1

Exclusion limits on the product of cross section and branching fraction for b* signal coupling f = 1.0, as a function of the resonance mass hypothesis.

Theoretical cross-section of b* signal for coupling multiplier f =1.0

Exclusion limits on the product of cross section and branching fraction for b* signal coupling f = 0.5, as a function of the resonance mass hypothesis.

Theoretical cross-section of b* signal for coupling multiplier f =0.5

Exclusion limits on the product of cross section and branching fraction for b* signal coupling f = 0.1, as a function of the resonance mass hypothesis.

Theoretical cross-section of b* signal for coupling multiplier f =0.1

Exclusion limits on the product of cross section and branching fraction for QBH ADD signal, as a function of the mass hypothesis.

Theoretical cross-section of ADD signal

Exclusion limits on the product of cross section and branching fraction for QBH RS1 signal, as a function of the mass hypothesis.

Theoretical cross-section of RS1 signal

The expected and observed exclusion mass limits for q* and b* signals, as a function coupling multiplier

The observed exclusion mass limits for b* signal, as a function of threshold mass of the signal.

Dijet mass spectra after the background only fit with the background prediction in the high-mass region with two b-tags.

Dijet mass spectra after the background only fit with the background prediction in the low-mass region with two b-tags.

The online b-tagging efficiency with respect to the offline b-tagging efficiency as a function of pT. The b-tagging online and offline working points correspond to an efficiency of 60% and 70%, respectively.

Observed and expected 95% credibility-level upper limits on the cross-section for the b* model in the high-mass region with inclusive b-jet selection.

Observed and expected 95% credibility-level upper limits on the cross-section times branching ratio for the SSM and leptophobic Z' models in the low- and high-mass region with two b-tags selection.

Observed and expected 95% credibility-level upper limits on the cross-section for DM Z' models in the low-mass region with two b-tags selection. The Z' is expected to decay to all five quark flavors other than the top quark and the mediator to SM quark coupling (gSM) equal to 0.1 is assumed.

Observed and expected 95% credibility-level upper limits on the cross-section times branching ratio for DM Z'->bb models in the high-mass region with two b-tags selection. The Z' is expected to decay to bb only and the mediator to SM quark coupling (gSM) equal to 0.25 is assumed.

Observed and expected 95% credibility-level upper limits on cross section times acceptance times branching ratio of X --> bb, including kinematic acceptance and b-tagging efficiencies, for resonances with intrinsic width smaller than the detector resolution. The width of the Gaussian reconstructed shape is dominated by the dijet mass resolution. The table shows the limits obtained from the high-mass inclusive one b-tag selection.

Observed and expected 95% credibility-level upper limits on cross section times acceptance times branching ratio of X --> bb, including kinematic acceptance and b-tagging efficiencies, for resonances with intrinsic width smaller than the detector resolution. The width of the Gaussian reconstructed shape is dominated by the dijet mass resolution. The table shows the limits obtained from the combined low- and high-mass two b-tags selection.

The mass distributions for the inclusive one b-tagged selection and two b-tagged selection using an integrated luminosity of 36.1 $fb^{-1}$. The inclusive one b-tagged Pythia8 MC distribution is normalized to the inclusive one b-tagged data. The two b-tagged Pythia8 MC distribution is normalized to the two b-tagged data. The systematic uncertainty band is dominated by the b-tagging scale factor and the b-jet energy scale.

Signal acceptance times efficiency for inclusive 1 b-tag and 2 b-tag categories as a function of the simulated b* and the Z' masses.

Signal acceptance times efficiency for 2 b-tag categories as a function of the simulated Z' masses.

The flavor composition of the simulated dijet background as a function of dijet mass before tagging. The fraction of the six combinations of the b-jet , c-jet and light-flavor jet are shown. All offline selections are applied.

The flavor composition of the simulated dijet background as a function of dijet mass with inclusive one b-tag. The fraction of the six combinations of the b-jet , c-jet and light-flavor jet are shown. All offline selections are applied.

The flavor composition of the simulated dijet background as a function of dijet mass with two b-tags. The fraction of the six combinations of the b-jet , c-jet and light-flavor jet are shown. All offline selections are applied.

Observed and expected 95% credibility-level upper limits on cross section times acceptance times branching ratio of X --> bb, including kinematic acceptance and b-tagging efficiencies, for resonances exhibiting a generic Gaussian shape at particle level. The table shows the limits obtained from the inclusive b-jet selection. The limits corresponding to Gaussian-shaped resonances with width of Γ(X)/m(X) = 3%.

Observed and expected 95% credibility-level upper limits on cross section times acceptance times branching ratio of X --> bb, including kinematic acceptance and b-tagging efficiencies, for resonances exhibiting a generic Gaussian shape at particle level. The table shows the limits obtained from the inclusive b-jet selection. The limits corresponding to Gaussian-shaped resonances with width of Γ(X)/m(X) = 7%.

Observed and expected 95% credibility-level upper limits on cross section times acceptance times branching ratio of X --> bb, including kinematic acceptance and b-tagging efficiencies, for resonances exhibiting a generic Gaussian shape at particle level. The table shows the limits obtained from the inclusive b-jet selection. The limits corresponding to Gaussian-shaped resonances with width of Γ(X)/m(X) = 10%.

Observed and expected 95% credibility-level upper limits on cross section times acceptance times branching ratio of X --> bb, including kinematic acceptance and b-tagging efficiencies, for resonances exhibiting a generic Gaussian shape at particle level. The table shows the limits obtained from the inclusive b-jet selection. The limits corresponding to Gaussian-shaped resonances with width of Γ(X)/m(X) = 15%.

Observed and expected 95% credibility-level upper limits on cross section times acceptance times branching ratio of X --> bb, including kinematic acceptance and b-tagging efficiencies, for resonances exhibiting a generic Gaussian shape at particle level. The table shows the limits obtained from the combined low- and high-mass two b-tags selection. The limits corresponding to Gaussian-shaped resonances with width of Γ(X)/m(X) = 3%.

Observed and expected 95% credibility-level upper limits on cross section times acceptance times branching ratio of X --> bb, including kinematic acceptance and b-tagging efficiencies, for resonances exhibiting a generic Gaussian shape at particle level. The table shows the limits obtained from the combined low- and high-mass two b-tags selection. The limits corresponding to Gaussian-shaped resonances with width of Γ(X)/m(X) = 7%.

Observed and expected 95% credibility-level upper limits on cross section times acceptance times branching ratio of X --> bb, including kinematic acceptance and b-tagging efficiencies, for resonances exhibiting a generic Gaussian shape at particle level. The table shows the limits obtained from the combined low- and high-mass two b-tags selection. The limits corresponding to Gaussian-shaped resonances with width of Γ(X)/m(X) = 10%.

Observed and expected 95% credibility-level upper limits on cross section times acceptance times branching ratio of X --> bb, including kinematic acceptance and b-tagging efficiencies, for resonances exhibiting a generic Gaussian shape at particle level. The table shows the limits obtained from the combined low- and high-mass two b-tags selection. The limits corresponding to Gaussian-shaped resonances with width of Γ(X)/m(X) = 15%.

The signal selection acceptance times efficiency (A$\times \epsilon$), defined as the ratio of the number of MC signal events in the category to the number of generated signal events, is presented as a function of the HVT resonance mass in the VBF category. The A$\times \epsilon$ is shown for the combination of all decay channels. For the HVT VBF samples, generator cuts are: m$_{jj} >$ 150 GeV. The decay of $W$ and $Z$ are to all flavors of leptons. The uncertainty shown represents the total statistical and systematic uncertainties.

Signal selection acceptance times efficiency (A$\times \epsilon$), defined as the ratio of the number of MC signal events in the category to the number of generated signal events, as a function of the HVT resonance mass for the quark-antiquark production. The $A\times\epsilon$ is shown for each decay channel and the inclusive one. The decay of $W$ and $Z$ are to all flavors of leptons. The uncertainty shown represents the total statistical and systematic uncertainties.

Signal selection acceptance times efficiency (A$\times \epsilon$), defined as the ratio of the number of MC signal events in the category to the number of generated signal events, as a function of the HVT resonance mass for the quark-antiquark production. The $A\times\epsilon$ is shown for each decay channel and the inclusive one. The decay of $W$ and $Z$ are to all flavors of leptons. The uncertainty shown represents the total statistical and systematic uncertainties.

Observed and expected distributions of the $WZ$ invariant mass in the q$\bar{\mathrm q}$ category after applying all selection cuts. Background contributions are obtained from background-only likelihood fit to the data.

Observed and expected distributions of the $WZ$ invariant mass in the q$\bar{\mathrm q}$ category after applying all selection cuts. Background contributions are obtained from background-only likelihood fit to the data.

Observed and expected distributions of the $WZ$ invariant mass in the $VBF$ category after applying all selection cuts. Background contributions are obtained from background-only likelihood fit to the data.

Observed and expected distributions of the $WZ$ invariant mass in the $VBF$ category after applying all selection cuts. Background contributions are obtained from background-only likelihood fit to the data.

Observed and expected 95% C.L. upper limits on $\sigma\times$BR($W'\to W^\pm Z$) for the q$\bar{\mathrm q}$ production of a $W'$ boson in the HVT model as a function of its mass. The $\pm1\sigma$ and $\pm2\sigma$ uncertainty in the expected limits are given. Limits are calculated in the asymptotic approximation below 900 GeV and are obtained from pseudo-experiments above that.

Observed and expected 95% C.L. upper limits on $\sigma\times$BR($W'\to W^\pm Z$) for the q$\bar{\mathrm q}$ production of a $W'$ boson in the HVT model as a function of its mass. The $\pm1\sigma$ and $\pm2\sigma$ uncertainty in the expected limits are given. Limits are calculated in the asymptotic approximation below 900 GeV and are obtained from pseudo-experiments above that.

Observed and expected 95% C.L. upper limits on $\sigma\times$BR($W'\to W^\pm Z$) for the $VBF$ production of a $W'$ boson in the HVT model as a function of its mass. The $\pm1\sigma$ and $\pm2\sigma$ uncertainty are given in the expected limits. Limits are calculated in the asymptotic approximation below 700 GeV and are obtained from pseudo-experiments above that.

Observed and expected 95% C.L. upper limits on $\sigma\times$BR($W'\to W^\pm Z$) for the $VBF$ production of a $W'$ boson in the HVT model as a function of its mass. The $\pm1\sigma$ and $\pm2\sigma$ uncertainty are given in the expected limits. Limits are calculated in the asymptotic approximation below 700 GeV and are obtained from pseudo-experiments above that.

Observed and expected 95% C.L. upper limits on $\sigma\times$BR($H_5^\pm\to W^\pm Z$) of the Georgi-Machacek Model as a function of $m_{H_5^\pm}$. The $\pm1\sigma$ and $\pm2\sigma$ uncertainty in the expected limits are given. Limits are calculated in the asymptotic approximation below 700 GeV and are obtained from pseudo-experiments above that.

Observed and expected 95% C.L. upper limits on $\sigma\times$BR($H_5^\pm\to W^\pm Z$) of the Georgi-Machacek Model as a function of $m_{H_5^\pm}$. The $\pm1\sigma$ and $\pm2\sigma$ uncertainty in the expected limits are given. Limits are calculated in the asymptotic approximation below 700 GeV and are obtained from pseudo-experiments above that.

Observed and expected 95% C.L. upper limits on the parameter $\sin(\theta_H)$ of the Georgi-Machacek Model as a function of $m_{H_5^\pm}$. The $\pm1\sigma$ and $\pm2\sigma$ uncertainty are given in the expected limits. Limits are calculated in the asymptotic approximation below 700 GeV and are obtained from pseudo-experiments above that.

Observed and expected 95% C.L. upper limits on the parameter $\sin(\theta_H)$ of the Georgi-Machacek Model as a function of $m_{H_5^\pm}$. The $\pm1\sigma$ and $\pm2\sigma$ uncertainty are given in the expected limits. Limits are calculated in the asymptotic approximation below 700 GeV and are obtained from pseudo-experiments above that.

Cumulative efficiency for simulated signal events to pass each stage of the event selection in the hadhad category. The efficiencies are calculated with respect to HH to bbtautau decays in which both tau-leptons decay hadronically. The ‘Pre-selection’ consists of basic requirements, including that at least two tau-had-vis pass loose kinematic requirements, at least one of the tau-had-vis candidate passes identification requirement, and that events do not contain an electron or muon. The ‘Object selections’ require exactly two tau-had-vis candidates, and at least two jets with pT > 25 GeV and abs(eta) < 2.5. The ‘Loose tau ID’ requires that both tau-had-vis candidates pass identification requirements. The ‘DTT offline jet cuts’ are cuts placed on the pT of the reconstructed jet or tau-had-vis that are geometrically matched to the HLT objects, to ensure the efficiencies of the HLT objects reach the plateau region.

Cumulative efficiency for simulated signal events to pass each stage of the event selection in the lephad SLT category. The efficiencies are calculated with respect to HH to bbtautau decays in which one tau-lepton decays hadronically and one decays leptonically. The ‘Pre-selection’ consists of basic requirements, including that at least one tau-had-vis candidate and one lepton pass loose kinematic requirements. The ‘Object selections’ require exactly one tau-had-vis candidate, and at least two jets with pT > 25 GeV and abs(eta) < 2.5. The ‘Trigger specific offline pT cuts’ are cuts placed on the pT of the reconstructed jet or tau-had-vis that are geometrically matched to the HLT objects, to ensure the efficiencies of the HLT objects reach the plateau region.

Cumulative efficiency for simulated signal events to pass each stage of the event selection in the lephad LTT category. The efficiencies are calculated with respect to HH to bbtautau decays in which one tau-lepton decays hadronically and one decays leptonically. The ‘Pre-selection’ consists of basic requirements, including that at least one tau-had-vis candidate and one lepton pass loose kinematic requirements. The ‘Object selections’ require exactly one tau-had-vis candidate, and at least two jets with pT > 25 GeV and abs(eta) < 2.5. The ‘Trigger specific offline pT cuts’ are cuts placed on the pT of the reconstructed jet or tau-had-vis that are geometrically matched to the HLT objects, to ensure the efficiencies of the HLT objects reach the plateau region.

Post-fit expected number of signal and background events and observed number of data events after applying the selection criteria and requiring exactly 2 b-tagged jets and assuming a background-only hypothesis.

Observed and expected limits at 95% CL on the cross-section of HH production, for the non-resonant ggF+VBF HH search, and the resonant HH search for four values of the resonance mass mX.

Acceptance times efficiency for the full analysis selections as a function of the resonance mass mX in the hadhad, lephad SLT and lephad LTT trigger categories, and the combined lephad. The acceptance times efficiency is evaluated for HH to bbtautau decays, with respect to the targeted tau-lepton decay modes (lephad or hadhad).

Post-fit distribution of mHH in the hadhad channel.

Post-fit distribution of mHH in the lephad SLT channel.

Post-fit distribution of mHH in the lephad LTT channel.

Post-fit distribution of mtautau MMC in the hadhad channel.

Post-fit distribution of mtautau MMC in the hadhad channel.

Post-fit distribution of mtautau MMC in the hadhad channel.

Post-fit distribution of of di-b-jet mass in the hadhad channel.

Post-fit distribution of di-b-jet mass in the LTT lephad channel.

Post-fit distribution of di-b-jet mass in the LTT lephad channel.

BDT for SM HH in the hadhad channel

NN for SM HH in the lephad SLT channel

NN for SM HH in the lephad LTT channel

PNN for mX = 500 GeV resonant HH in the hadhad channel

PNN for mX = 500 GeV resonant HH in the lephad SLT channel

PNN for mX = 500 GeV resonant HH in the lephad LTT channel

PNN for mX = 1000 GeV resonant HH in the hadhad channel

PNN for mX = 1000 GeV resonant HH in the lephad SLT channel

PNN for mX = 1000 GeV resonant HH in the lephad LTT channel

Event yields as a function of log10(S/B) for data, background and non-resonant HH signal. Final discriminant bins from the hadhad, lephad SLT and lephad LTT categories are combined into bins of log10(S/B). The B is the fitted background yield assuming background-only hypothesis, and the signal S is scaled to the SM expected cross-section.

Observed and expected limits at 95% CL on the cross-section of the resonant HH production as a function of the scalar resonance mass mX.

Post-fit distribution of delta-R between the taus in the hadhad channel.

Post-fit distribution of delta-R between the b-tagged jets in the hadhad channel.

Post-fit distribution of delta-R between the taus in the lephad SLT channel.

Post-fit distribution of delta-R between the b-tagged jets in the lephad SLT channel.

Post-fit distribution of delta-pT between the tau and lepton in the lephad SLT channel.

Post-fit distribution of pT of the subleading b-tagged jet in the lephad SLT channel.

Post-fit distribution of MTW in the lephad SLT channel.

Post-fit distribution of missing transverse momentum in the lephad SLT channel.

Post-fit distribution of missing transverse momentum centrality in the lephad SLT channel.

Post-fit distribution of delta-phi between the Higgs boson candidates in the lephad SLT channel.

Post-fit distribution of delta-pT between the tau and lepton in the lephad LTT channel.

Post-fit distribution of delta-R between the taus in the lephad LTT channel.

Post-fit distribution of delta-phi between the lepton and the missing transverse momentum in the lephad LTT channel.

Post-fit distribution of delta-phi between the Higgs boson candidates in the lephad LTT channel.

Post-fit distribution of the total transverse momentum s in the lephad LTT channel.

PNN for mX = 300 GeV resonant HH in the hadhad channel

PNN for mX = 300 GeV resonant HH in the lephad SLT channel

PNN for mX = 300 GeV resonant HH in the lephad LTT channel

PNN for mX = 1600 GeV resonant HH in the hadhad channel

PNN for mX = 1600 GeV resonant HH in the lephad SLT channel

PNN for mX = 1600 GeV resonant HH in the lephad LTT channel

Local p-value of the background-only hypothesis as a function of the resonance mass.

The observed 95% CL upper limits on the product of the cross section, branching fraction, and acceptance for resonant production of paired dijet resonances decaying to a quark-gluon pair, with $M(X)/M(Y) = 0.17$. The corresponding expected limits and their variations at the 1 and 2 standard deviation levels are also shown. Limits are compared to predictions for a scalar diquark with couplings to pairs of up quarks, $y_{uu}$ = 0.4, and to pairs of vector-like quarks, $y_{χ}$ = 0.6.

The observed 95% CL upper limits on the product of the cross section, branching fraction, and acceptance for resonant production of paired dijet resonances decaying to a quark-gluon pair, with $M(X)/M(Y) = 0.19$. The corresponding expected limits and their variations at the 1 and 2 standard deviation levels are also shown. Limits are compared to predictions for a scalar diquark with couplings to pairs of up quarks, $y_{uu}$ = 0.4, and to pairs of vector-like quarks, $y_{χ}$ = 0.6.

The observed 95% CL upper limits on the product of the cross section, branching fraction, and acceptance for resonant production of paired dijet resonances decaying to a quark-gluon pair, with $M(X)/M(Y) = 0.21$. The corresponding expected limits and their variations at the 1 and 2 standard deviation levels are also shown. Limits are compared to predictions for a scalar diquark with couplings to pairs of up quarks, $y_{uu}$ = 0.4, and to pairs of vector-like quarks, $y_{χ}$ = 0.6.

The observed 95% CL upper limits on the product of the cross section, branching fraction, and acceptance for resonant production of paired dijet resonances decaying to a quark-gluon pair, with $M(X)/M(Y) = 0.23$. The corresponding expected limits and their variations at the 1 and 2 standard deviation levels are also shown. Limits are compared to predictions for a scalar diquark with couplings to pairs of up quarks, $y_{uu}$ = 0.4, and to pairs of vector-like quarks, $y_{χ}$ = 0.6.

The observed 95% CL upper limits on the product of the cross section, branching fraction, and acceptance for resonant production of paired dijet resonances decaying to a quark-gluon pair, with $M(X)/M(Y) = 0.25$. The corresponding expected limits and their variations at the 1 and 2 standard deviation levels are also shown. Limits are compared to predictions for a scalar diquark with couplings to pairs of up quarks, $y_{uu}$ = 0.4, and to pairs of vector-like quarks, $y_{χ}$ = 0.6.

The observed 95% CL upper limits on the product of the cross section, branching fraction, and acceptance for resonant production of paired dijet resonances decaying to a quark-gluon pair, with $M(X)/M(Y) = 0.27$. The corresponding expected limits and their variations at the 1 and 2 standard deviation levels are also shown. Limits are compared to predictions for a scalar diquark with couplings to pairs of up quarks, $y_{uu}$ = 0.4, and to pairs of vector-like quarks, $y_{χ}$ = 0.6.

The observed 95% CL upper limits on the product of the cross section, branching fraction, and acceptance for resonant production of paired dijet resonances decaying to a quark-gluon pair, with $M(X)/M(Y) = 0.29$. The corresponding expected limits and their variations at the 1 and 2 standard deviation levels are also shown. Limits are compared to predictions for a scalar diquark with couplings to pairs of up quarks, $y_{uu}$ = 0.4, and to pairs of vector-like quarks, $y_{χ}$ = 0.6.

The observed 95% CL upper limits on the product of the cross section, branching fraction, and acceptance for resonant production of paired dijet resonances decaying to a quark-gluon pair, with $M(X)/M(Y) = 0.31$. The corresponding expected limits and their variations at the 1 and 2 standard deviation levels are also shown. Limits are compared to predictions for a scalar diquark with couplings to pairs of up quarks, $y_{uu}$ = 0.4, and to pairs of vector-like quarks, $y_{χ}$ = 0.6.

The observed 95% CL upper limits on the product of the cross section, branching fraction, and acceptance for resonant production of paired dijet resonances decaying to a quark-gluon pair, with $M(X)/M(Y) = 0.33$. The corresponding expected limits and their variations at the 1 and 2 standard deviation levels are also shown. Limits are compared to predictions for a scalar diquark with couplings to pairs of up quarks, $y_{uu}$ = 0.4, and to pairs of vector-like quarks, $y_{χ}$ = 0.6.

The observed 95% CL upper limits on the product of the cross section, branching fraction, and acceptance for resonant production of paired dijet resonances decaying to a quark-gluon pair, with $M(X)/M(Y) = 0.42$. The corresponding expected limits and their variations at the 1 and 2 standard deviation levels are also shown. Limits are compared to predictions for a scalar diquark with couplings to pairs of up quarks, $y_{uu}$ = 0.4, and to pairs of vector-like quarks, $y_{χ}$ = 0.6.

The observed 95% CL upper limits on the product of the cross section, branching fraction, and acceptance for the non-resonant production of top squark pairs in the RPV SUSY decay scenario. The corresponding expected limits and their variations at the 1 and 2 standard deviation levels are also shown. Limits are compared to the top squark model cross section.

Observed differential four-jet mass spectrum for 0.22 < $\\a$ < 0.24. The cross-section is calculated by dividing the event yield by the bin width and luminosity.

Observed differential four-jet mass spectrum for 0.24 < $\\a$ < 0.26. The cross-section is calculated by dividing the event yield by the bin width and luminosity.

Observed differential four-jet mass spectrum for 0.26 < $\\a$ < 0.28. The cross-section is calculated by dividing the event yield by the bin width and luminosity.

Observed differential four-jet mass spectrum for all $\\a$ bins together. The cross-section is calculated by dividing the event yield by the bin width and luminosity.

Observed differential average dijet mass spectrum within three $\\a$ bins of the non-resonant search. The cross-section is calculated by dividing the event yield by the bin width and luminosity.

Corrected $m(HH)$ distribution in the resolved $4b$ signal region (dots), after the fit under the background-only hypothesis. The error bars on the $4b$ points represent the Poisson uncertainties corresponding to their event yields. The background model (teal histogram) is shown with its total post-fit uncertainty (gray band). The final bin includes overflow. Representative spin-0 signal hypotheses (dashed, dotted, and dashed-dotted lines) are overlaid, normalized to the overall expected limits on their cross-sections. The bottom panel shows the difference between the $4b$ distribution and the background model, relative to the background model. No significant excess of data relative to the SM background is observed.

Corrected $m(HH)$ distribution in the resolved $4b$ signal region (dots), after the fit under the background-only hypothesis. The error bars on the $4b$ points represent the Poisson uncertainties corresponding to their event yields. The background model (teal histogram) is shown with its total post-fit uncertainty (gray band). The final bin includes overflow. Representative spin-2 signal hypotheses (dashed, dotted, and dashed-dotted lines) are overlaid, normalized to the overall expected limits on their cross-sections. The bottom panel shows the difference between the $4b$ distribution and the background model, relative to the background model. No significant excess of data relative to the SM background is observed.

Cumulative signal acceptance times efficiency as a function of the resonance mass for various selection steps in the boosted channel. The steps up to the $b$-tag categorization are shown for the spin-0.

Cumulative signal acceptance times efficiency as a function of the resonance mass for various selection steps in the boosted channel. The steps up to the $b$-tag categorization are shown for the spin-2.

Cumulative signal acceptance times efficiency as a function of the resonance mass for various selection steps in the boosted channel. The efficiencies of the three b-tag categories are shown for the spin-0 scenario; this efficiency is obtained after the other selection steps including the SR definition. The signal efficiency in the 4b region has a maximum around 1.5 TeV. Above that value the track-jets start to merge together, and for the highest resonance masses the 2b category becomes the most efficient.

Cumulative signal acceptance times efficiency as a function of the resonance mass for various selection steps in the boosted channel. The efficiencies of the three b-tag categories are shown for the spin-2 scenario; this efficiency is obtained after the other selection steps including the SR definition. The signal efficiency in the 4b region has a maximum around 1.5 TeV. Above that value the track-jets start to merge together, and for the highest resonance masses the 2b category becomes the most efficient.

Comparison of the background model (stacked histograms) with data (dots) in the $2b$ validation region. The error bars on the data points represent the Poisson uncertainties corresponding to their event yields. The background uncertainty (gray band) is computed by adding all individual components in quadrature and is not allowed to extend below zero.

Comparison of the background model (stacked histograms) with data (dots) in the $3b$ validation region. The error bars on the data points represent the Poisson uncertainties corresponding to their event yields. The background uncertainty (gray band) is computed by adding all individual components in quadrature and is not allowed to extend below zero.

Comparison of the background model (stacked histograms) with data (dots) in the $4b$ validation region. The error bars on the data points represent the Poisson uncertainties corresponding to their event yields. The background uncertainty (gray band) is computed by adding all individual components in quadrature and is not allowed to extend below zero.

The $m(HH)$ distributions in the boosted $2b$ signal regions (dots), after the fit under the background-only hypothesis. The error bars on the data points represent the Poisson uncertainties corresponding to their event yields. The background model (stacked histogram) is shown with its total post-fit uncertainty (gray band). The uncertainty bands are defined using an ensemble of curves constructed by sampling a multivariate Gaussian probability density function built from the covariance matrix of the fit. Representative spin-0 signal hypotheses (dashed and dashed-dotted lines) are overlaid, normalized to the overall expected limits on their cross-sections. No significant excess of data relative to the SM background is observed.

The $m(HH)$ distributions in the boosted $2b$ signal regions (dots), after the fit under the background-only hypothesis. The error bars on the data points represent the Poisson uncertainties corresponding to their event yields. The background model (stacked histogram) is shown with its total post-fit uncertainty (gray band). The uncertainty bands are defined using an ensemble of curves constructed by sampling a multivariate Gaussian probability density function built from the covariance matrix of the fit. Representative spin-2 signal hypotheses (dashed and dashed-dotted lines) are overlaid, normalized to the overall expected limits on their cross-sections. No significant excess of data relative to the SM background is observed.

The $m(HH)$ distributions in the boosted $3b$ signal regions (dots), after the fit under the background-only hypothesis. The error bars on the data points represent the Poisson uncertainties corresponding to their event yields. The background model (stacked histogram) is shown with its total post-fit uncertainty (gray band). The uncertainty bands are defined using an ensemble of curves constructed by sampling a multivariate Gaussian probability density function built from the covariance matrix of the fit. Representative spin-0 signal hypotheses (dashed and dashed-dotted lines) are overlaid, normalized to the overall expected limits on their cross-sections. No significant excess of data relative to the SM background is observed.

The $m(HH)$ distributions in the boosted $3b$ signal regions (dots), after the fit under the background-only hypothesis. The error bars on the data points represent the Poisson uncertainties corresponding to their event yields. The background model (stacked histogram) is shown with its total post-fit uncertainty (gray band). The uncertainty bands are defined using an ensemble of curves constructed by sampling a multivariate Gaussian probability density function built from the covariance matrix of the fit. Representative spin-2 signal hypotheses (dashed and dashed-dotted lines) are overlaid, normalized to the overall expected limits on their cross-sections. No significant excess of data relative to the SM background is observed.

The $m(HH)$ distributions in the boosted $4b$ signal regions (dots), after the fit under the background-only hypothesis. The error bars on the data points represent the Poisson uncertainties corresponding to their event yields. The background model (stacked histogram) is shown with its total post-fit uncertainty (gray band). The uncertainty bands are defined using an ensemble of curves constructed by sampling a multivariate Gaussian probability density function built from the covariance matrix of the fit. Representative spin-0 signal hypotheses (dashed and dashed-dotted lines) are overlaid, normalized to the overall expected limits on their cross-sections. No significant excess of data relative to the SM background is observed.

The $m(HH)$ distributions in the boosted $4b$ signal regions (dots), after the fit under the background-only hypothesis. The error bars on the data points represent the Poisson uncertainties corresponding to their event yields. The background model (stacked histogram) is shown with its total post-fit uncertainty (gray band). The uncertainty bands are defined using an ensemble of curves constructed by sampling a multivariate Gaussian probability density function built from the covariance matrix of the fit. Representative spin-2 signal hypotheses (dashed and dashed-dotted lines) are overlaid, normalized to the overall expected limits on their cross-sections. No significant excess of data relative to the SM background is observed.

Expected (dashed black lines) and observed (solid black lines) 95% CL upper limits on the cross-section of resonant $HH$ production in the spin-0 signal models. The $\pm 1 \sigma$ and $\pm 2 \sigma$ uncertainty ranges for the expected limits (colored bands) are shown. Expected limits using each of the resolved and boosted channels individually (dashed colored lines) are shown. The nominal $H\rightarrow b\bar{b}$ branching ratio is taken as 0.582.

Expected (dashed black lines) and observed (solid black lines) 95% CL upper limits on the cross-section of resonant $HH$ production in the spin-2 signal models. The $\pm 1 \sigma$ and $\pm 2 \sigma$ uncertainty ranges for the expected limits (colored bands) are shown. Expected limits using each of the resolved and boosted channels individually (dashed colored lines) are shown. The theoretical prediction for the bulk RS model with $k/\bar{M}_{\text{Pl}} = 1$ (solid red line) is shown; the decrease below 350 GeV is due to a sharp reduction in the $G^{*}_{\text{KK}} \rightarrow HH$ branching ratio. The nominal $H\rightarrow b\bar{b}$ branching ratio is taken as 0.582.