Correlation matrix for the individual sources of systematic uncertainty considered in the analysis, for the combined opposite sign (OS) and same sign (SS) binned-template profile likelihood fit to data. The OS or SS refers to the charge signs of the primary lepton and the soft muon. The gamma parameters are NPs used to describe the effect of the limited statistics of the samples.

Observed 95% CL limits on the production cross-section for scalar up-type LQs decaying into a top quark and a neutrino or a bottom quark and an electron.

Expected 95% CL limits on the production cross-section for scalar up-type LQs decaying into a top quark and a neutrino or a bottom quark and an electron.

Observed 95% CL limits on the production cross-section for scalar down-type LQs decaying into a bottom quark and a neutrino or a top quark and a muon.

Expected 95% CL limits on the production cross-section for scalar down-type LQs decaying into a bottom quark and a neutrino or a top quark and a muon.

Observed 95% CL limits on the production cross-section for scalar down-type LQs decaying into a bottom quark and a neutrino or a top quark and an electron.

Expected 95% CL limits on the production cross-section for scalar down-type LQs decaying into a bottom quark and a neutrino or a top quark and an electron.

Observed 95% CL limits on the production cross-section for vector up-type LQs in the Yang-Mills coupling scenario decaying into a top quark and a neutrino or a bottom quark and a muon.

Expected 95% CL limits on the production cross-section for vector up-type LQs in the Yang-Mills coupling scenario decaying into a top quark and a neutrino or a bottom quark and a muon.

Observed 95% CL limits on the production cross-section for vector up-type LQs in the Yang-Mills coupling scenario decaying into a top quark and a neutrino or a bottom quark and an electron.

Expected 95% CL limits on the production cross-section for vector up-type LQs in the Yang-Mills coupling scenario decaying into a top quark and a neutrino or a bottom quark and an electron.

Observed 95% CL limits on the production cross-section for vector up-type LQs in the minimal coupling scenario decaying into a top quark and a neutrino or a bottom quark and a muon.

Expected 95% CL limits on the production cross-section for vector up-type LQs in the minimal coupling scenario decaying into a top quark and a neutrino or a bottom quark and a muon.

Observed 95% CL limits on the production cross-section for vector up-type LQs in the minimal coupling scenario decaying into a top quark and a neutrino or a bottom quark and an electron.

Expected 95% CL limits on the production cross-section for vector up-type LQs in the minimal coupling scenario decaying into a top quark and a neutrino or a bottom quark and an electron.

Product of signal acceptance and efficiency in the training region for scalar up-type LQs decaying into top quarks and neutrinos or bottom quarks and muons depending on the signal mass, $m_{LQ}$, and $\mathcal{B}(LQ_{mix}^u\rightarrow b\mu)$.

Product of signal acceptance and efficiency in the training region for scalar up-type LQs decaying into top quarks and neutrinos or bottom quarks and electrons depending on the signal mass, $m_{LQ}$, and $\mathcal{B}(LQ_{mix}^u\rightarrow be)$.

Product of signal acceptance and efficiency in the training region for scalar down-type LQs decaying into top quarks and muons or bottom quarks and neutrinos depending on the signal mass, $m_{LQ}$, and $\mathcal{B}(LQ_{mix}^d\rightarrow t\mu)$.

Product of signal acceptance and efficiency in the training region for scalar down-type LQs decaying into top quarks and electrons or bottom quarks and neutrinos depending on the signal mass, $m_{LQ}$, and $\mathcal{B}(LQ_{mix}^d\rightarrow te)$.

Product of signal acceptance and efficiency in the training region for vector LQs in the Yang-Mills coupling scenario decaying into top quarks and neutrinos or bottom quarks and muons depending on the signal mass, $m_{LQ}$, and $\mathcal{B}(vLQ_{mix}^{YM}\rightarrow b\mu)$.

Product of signal acceptance and efficiency in the training region for vector LQs in the Yang-Mills coupling scenario decaying into top quarks and neutrinos or bottom quarks and electrons depending on the signal mass, $m_{LQ}$, and $\mathcal{B}(vLQ_{mix}^{YM}\rightarrow be)$.

Product of signal acceptance and efficiency in the training region for vector LQs in the minimal coupling scenario decaying into top quarks and neutrinos or bottom quarks and muons depending on the signal mass, $m_{LQ}$, and $\mathcal{B}(vLQ_{mix}^{min}\rightarrow b\mu)$.

Product of signal acceptance and efficiency in the training region for vector LQs in the minimal coupling scenario decaying into top quarks and neutrinos or bottom quarks and electrons depending on the signal mass, $m_{LQ}$, and $\mathcal{B}(vLQ_{mix}^{min}\rightarrow be)$.

Cut-flow for the four models containing scalar LQs at a mass of $m_{LQ}=1300$ GeV and $\mathcal{B}(LQ_{mix}\rightarrow q_3\ell)=0.5$.

Cut-flow for the four models containing vector LQs at a mass of $m_{LQ}=1300$ GeV and $\mathcal{B}(LQ_{mix}\rightarrow q_3\ell)=0.5$.

Expected and observed 95% CL for the 1SFH mu Majorana model.

Expected and observed 95% CL for the 2QDH NH Dirac model.

Expected and observed 95% CL for the 2QDH NH Majorana model.

Expected and observed 95% CL for the 2QDH IH Dirac model.

Expected and observed 95% CL for the 2QDH IH Majorana model.

Cutflow for six simulated signal channels showing the weighted number of expected events based on the single-flavour mixing model in the Majorana limit. Each column uses the generated signal sample with the mass hypothesis $m_N = 10$ GeV and proper decay length $c\tau_N = 10$ mm.

Cutflow for the six channels in data showing the number of events passing each successive signal selection for Majorana HNLs.

The event selection efficiency for each mass-lifetime point in all six studied channels. Shown is the fraction of the produced MC simulation events that pass all signal region selections. An entry of 0 indicates no events were selected.

The dominant signal uncertainty is due to differences in reconstruction of displaced vertices and tracks between data and MC. This is evaluated by comparing $K^{0}_{S} \rightarrow \pi^+\pi^-$ event yields in the VR and in MC produced with Pythia8.186 in bins of $p_\mathrm{T}$ and $r_\mathrm{DV}$. The data/MC ratio is normalized to the bin nearest the IP where the tracking and vertexing reconstruction algorithms are expected to be most robust. The symmetrized difference from 1.0 is applied to each signal vertex as a per-event systematic variation.

Expected and observed yields in the different analysis regions (prefit) for the 1SFH e Dirac model (10 GeV, 10mm).

Expected and observed yields in the different analysis regions (postfit) for the 1SFH e Dirac model (10 GeV, 10mm).

Expected and observed yields in the different analysis regions (prefit) for the 1SFH e Majorana model (10 GeV, 10mm).

Expected and observed yields in the different analysis regions (postfit) for the 1SFH e Majorana model (10 GeV, 10mm).

Expected and observed yields in the different analysis regions (prefit) for the 1SFH u Dirac model (10 GeV, 10mm).

Expected and observed yields in the different analysis regions (postfit) for the 1SFH u Dirac model (10 GeV, 10mm).

Expected and observed yields in the different analysis regions (prefit) for the 1SFH u Majorana model (10 GeV, 10mm).

Expected and observed yields in the different analysis regions (postfit) for the 1SFH u Majorana model (10 GeV, 10mm).

Expected and observed yields in the different analysis regions (prefit) for the 2QDH (NH) Dirac model (10 GeV, 10mm).

Expected and observed yields in the different analysis regions (postfit) for the 2QDH (NH) Dirac model (10 GeV, 10mm).

Expected and observed yields in the different analysis regions (prefit) for the 2QDH (NH) Majorana model (10 GeV, 10mm).

Expected and observed yields in the different analysis regions (postfit) for the 2QDH (NH) Majorana model (10 GeV, 10mm).

Expected and observed yields in the different analysis regions (prefit) for the 2QDH (IH) Dirac model (10 GeV, 10mm).

Expected and observed yields in the different analysis regions (postfit) for the 2QDH (IH) Dirac model (10 GeV, 10mm).

Expected and observed yields in the different analysis regions (prefit) for the 2QDH (IH) Majorana model (10 GeV, 10mm).

Expected and observed yields in the different analysis regions (postfit) for the 2QDH (IH) Majorana model (10 GeV, 10mm).

The total displaced vertexing efficiency as a function of $r_{DV}$ for the custom configuration used in this analysis. The definition of the secondary vertex efficiency can be found in defined in \cite{ATL-PHYS-PUB-2019-013}. The efficiency is shown for $\mu-\mu\mu$, $\mu-\mu e$ and $\mu-ee$ signals with $m_N=10$~GeV and $c\tau_N=10$~mm.

Acceptance * reconstruction efficiency for the P P --> bbA --> Zh --> llbb signals in the 2-lepton channel.

Acceptance * reconstruction efficiency for the P P --> Wprime --> Zh --> lvbb/cc signals in the 0-lepton channel.

Acceptance * reconstruction efficiency for the P P --> Wprime --> Zh --> lvbb/cc signals in the 1-lepton channel.

Event distributions of $m_{T,Vh}$ for the 0-lepton channel in the resolved 1 b-tag signal region. The background prediction is shown after a background-only maximum-likelihood Z' fit to the data. In the plot, the last bin contains the overflow.

Event distributions of $m_{T,Vh}$ for the 0-lepton channel in the resolved 2 b-tag signal region. The background prediction is shown after a background-only maximum-likelihood Z' fit to the data. In the plot, the last bin contains the overflow.

Event distributions of $m_{Vh}$ for the 2-lepton channel in the resolved 1 b-tag signal region. The background prediction is shown after a background-only maximum-likelihood Z' fit to the data. In the plot, the last bin contains the overflow.

Event distributions of $m_{Vh}$ for the 2-lepton channel in the resolved 2 b-tag signal region. The background prediction is shown after a background-only maximum-likelihood Z' fit to the data. In the plot, the last bin contains the overflow.

Event distributions of $m_{T,Vh}$ for the 0-lepton channel in the merged 1 b-tag signal region. The background prediction is shown after a background-only maximum-likelihood Z' fit to the data. In the plot, the last bin contains the overflow.

Event distributions of $m_{T,Vh}$ for the 0-lepton channel in the merged 2 b-tag signal region. The background prediction is shown after a background-only maximum-likelihood Z' fit to the data. In the plot, the last bin contains the overflow.

Event distributions of $m_{Vh}$ for the 2-lepton channel in the merged 1 b-tag signal region. The background prediction is shown after a background-only maximum-likelihood Z' fit to the data. In the plot, the last bin contains the overflow.

Event distributions of $m_{Vh}$ for the 2-lepton channel in the merged 2 b-tag signal region. The background prediction is shown after a background-only maximum-likelihood Z' fit to the data. In the plot, the last bin contains the overflow.

Event distributions of $m_{T,Vh}$ for the 1-lepton channel in the resolved 1 b-tag signal region. The background prediction is shown after a background-only maximum-likelihood W' fit to the data. In the plot, the last bin contains the overflow.

Event distributions of $m_{T,Vh}$ for the 1-lepton channel in the resolved 2 b-tag signal region. The background prediction is shown after a background-only maximum-likelihood W' fit to the data. In the plot, the last bin contains the overflow.

Event distributions of $m_{T,Vh}$ for the 1-lepton channel in the merged 1 b-tag signal region. The background prediction is shown after a background-only maximum-likelihood W' fit to the data. In the plot, the last bin contains the overflow.

Event distributions of $m_{T,Vh}$ for the 1-lepton channel in the merged 2 b-tag signal region. The background prediction is shown after a background-only maximum-likelihood W' fit to the data. In the plot, the last bin contains the overflow.

Event distributions of $m_{T,Vh}$ for the 0-lepton channel in the resolved 3+ b-tag signal region. The background prediction is shown after a background-only maximum-likelihood bbA fit to the data. In the plot, the last bin contains the overflow.

Event distributions of $m_{Vh}$ for the 2-lepton channel in the resolved 3+ b-tag signal region. The background prediction is shown after a background-only maximum-likelihood bbA fit to the data. In the plot, the last bin contains the overflow.

Event distributions of $m_{T,Vh}$ for the 0-lepton channel in the merged 2 b-tag signal region with additional b-tagged track jets not associated with the large-R jet. The background prediction is shown after a background-only maximum-likelihood bbA fit to the data. In the plot, the last bin contains the overflow.

Event distributions of $m_{Vh}$ for the 2-lepton channel in the merged 1+2 b-tag signal region with additional b-tagged track jets not associated with the large-R jet. The background prediction is shown after a background-only maximum-likelihood bbA fit to the data. In the plot, the last bin contains the overflow.

Event distributions of $m_{T,Vh}$ for the 0-lepton channel in the resolved 1 b-tag sideband control region. The background prediction is shown after a background-only maximum-likelihood W' fit to the data. In the plot, the last bin contains the overflow.

Event distributions of $m_{T,Vh}$ for the 0-lepton channel in the resolved 2 b-tag sideband control region. The background prediction is shown after a background-only maximum-likelihood W' fit to the data. In the plot, the last bin contains the overflow.

Event distributions of $m_{T,Vh}$ for the 1-lepton channel in the resolved 1 b-tag sideband control region. The background prediction is shown after a background-only maximum-likelihood W' fit to the data. In the plot, the last bin contains the overflow.

Event distributions of $m_{T,Vh}$ for the 1-lepton channel in the resolved 2 b-tag sideband control region. The background prediction is shown after a background-only maximum-likelihood W' fit to the data. In the plot, the last bin contains the overflow.

Event distributions of $m_{T,Vh}$ for the 0-lepton channel in the merged 1 b-tag sideband control region. The background prediction is shown after a background-only maximum-likelihood W' fit to the data. In the plot, the last bin contains the overflow.

Event distributions of $m_{T,Vh}$ for the 0-lepton channel in the merged 2 b-tag sideband control region. The background prediction is shown after a background-only maximum-likelihood W' fit to the data. In the plot, the last bin contains the overflow.

Event distributions of $m_{T,Vh}$ for the 1-lepton channel in the merged 1 b-tag sideband control region. The background prediction is shown after a background-only maximum-likelihood W' fit to the data. In the plot, the last bin contains the overflow.

Event distributions of $m_{T,Vh}$ for the 1-lepton channel in the merged 2 b-tag sideband control region. The background prediction is shown after a background-only maximum-likelihood W' fit to the data. In the plot, the last bin contains the overflow.

Event distributions of $m_{T,Vh}$ for the 0-lepton channel in the resolved 1 b-tag sideband control region. The background prediction is shown after a background-only maximum-likelihood Z' fit to the data. In the plot, the last bin contains the overflow.

Event distributions of $m_{T,Vh}$ for the 0-lepton channel in the resolved 2 b-tag sideband control region. The background prediction is shown after a background-only maximum-likelihood Z' fit to the data. In the plot, the last bin contains the overflow.

Event distributions of $m_{Vh}$ for the 2-lepton channel in the resolved top control region. The background prediction is shown after a background-only maximum-likelihood Z' fit to the data. In the plot, the last bin contains the overflow.

Event distributions of $m_{T,Vh}$ for the 0-lepton channel in the merged 1 b-tag sideband control region. The background prediction is shown after a background-only maximum-likelihood Z' fit to the data. In the plot, the last bin contains the overflow.

Event distributions of $m_{T,Vh}$ for the 0-lepton channel in the merged 2 b-tag sideband control region. The background prediction is shown after a background-only maximum-likelihood Z' fit to the data. In the plot, the last bin contains the overflow.

Upper limits on Zprime to Z h production cross section times branching fraction in pb.

Upper limits on Wprime to W h production cross section times branching fraction in pb.

Upper limits on ggA to Z h production cross section times branching fraction in pb.

Upper limits on bbA to Z h production cross section times branching fraction in pb.

Expected and observed two-dimensional likelihood scans of the b-associated production cross section times branching fraction vs the gluon-fusion production cross section times branching fraction at $m_{A}$ = 220 GeV.

Expected and observed two-dimensional likelihood scans of the b-associated production cross section times branching fraction vs the gluon-fusion production cross section times branching fraction at $m_{A}$ = 260 GeV.

Expected and observed two-dimensional likelihood scans of the b-associated production cross section times branching fraction vs the gluon-fusion production cross section times branching fraction at $m_{A}$ = 300 GeV.

Expected and observed two-dimensional likelihood scans of the b-associated production cross section times branching fraction vs the gluon-fusion production cross section times branching fraction at $m_{A}$ = 340 GeV.

Expected and observed two-dimensional likelihood scans of the b-associated production cross section times branching fraction vs the gluon-fusion production cross section times branching fraction at $m_{A}$ = 380 GeV.

Expected and observed two-dimensional likelihood scans of the b-associated production cross section times branching fraction vs the gluon-fusion production cross section times branching fraction at $m_{A}$ = 400 GeV.

Expected and observed two-dimensional likelihood scans of the b-associated production cross section times branching fraction vs the gluon-fusion production cross section times branching fraction at $m_{A}$ = 420 GeV.

Expected and observed two-dimensional likelihood scans of the b-associated production cross section times branching fraction vs the gluon-fusion production cross section times branching fraction at $m_{A}$ = 440 GeV.

Expected and observed two-dimensional likelihood scans of the b-associated production cross section times branching fraction vs the gluon-fusion production cross section times branching fraction at $m_{A}$ = 460 GeV.

Expected and observed two-dimensional likelihood scans of the b-associated production cross section times branching fraction vs the gluon-fusion production cross section times branching fraction at $m_{A}$ = 500 GeV.

Expected and observed two-dimensional likelihood scans of the b-associated production cross section times branching fraction vs the gluon-fusion production cross section times branching fraction at $m_{A}$ = 600 GeV.

Expected and observed two-dimensional likelihood scans of the b-associated production cross section times branching fraction vs the gluon-fusion production cross section times branching fraction at $m_{A}$ = 700 GeV.

Expected and observed two-dimensional likelihood scans of the b-associated production cross section times branching fraction vs the gluon-fusion production cross section times branching fraction at $m_{A}$ = 800 GeV.

Expected and observed two-dimensional likelihood scans of the b-associated production cross section times branching fraction vs the gluon-fusion production cross section times branching fraction at $m_{A}$ = 900 GeV.

Expected and observed two-dimensional likelihood scans of the b-associated production cross section times branching fraction vs the gluon-fusion production cross section times branching fraction at $m_{A}$ = 1000 GeV.

Expected and observed two-dimensional likelihood scans of the b-associated production cross section times branching fraction vs the gluon-fusion production cross section times branching fraction at $m_{A}$ = 1200 GeV.

Expected and observed two-dimensional likelihood scans of the b-associated production cross section times branching fraction vs the gluon-fusion production cross section times branching fraction at $m_{A}$ = 1400 GeV.

Expected and observed two-dimensional likelihood scans of the b-associated production cross section times branching fraction vs the gluon-fusion production cross section times branching fraction at $m_{A}$ = 1600 GeV.

Expected and observed two-dimensional likelihood scans of the b-associated production cross section times branching fraction vs the gluon-fusion production cross section times branching fraction at $m_{A}$ = 2000 GeV.

Acceptance * reconstruction efficiency for the P P --> A --> Zh --> vvbb signal in the 0-lepton channel.

Acceptance * reconstruction efficiency for the P P --> A --> Zh --> llbb signal in the 2-lepton channel.

Event distributions of $m_{T,Vh}$ for the 0-lepton channel in the resolved 1 b-tag signal region. The background prediction is shown after a background-only maximum-likelihood W' fit to the data. In the plot, the last bin contains the overflow.

Event distributions of $m_{T,Vh}$ for the 0-lepton channel in the resolved 2 b-tag signal region. The background prediction is shown after a background-only maximum-likelihood W' fit to the data. In the plot, the last bin contains the overflow.

Event distributions of $m_{T,Vh}$ for the 0-lepton channel in the merged 1 b-tag signal region. The background prediction is shown after a background-only maximum-likelihood W' fit to the data. In the plot, the last bin contains the overflow.

Event distributions of $m_{T,Vh}$ for the 0-lepton channel in the merged 2 b-tag signal region. The background prediction is shown after a background-only maximum-likelihood W' fit to the data. In the plot, the last bin contains the overflow.

Distributions of expected upper limits at 95% confidence level on the cross section of P P --> A --> Zh as a function of bbA fraction an signal mass.

Distributions of observed upper limits at 95% confidence level on the cross section of P P --> A --> Zh as a function of bbA fraction an signal mass.

Best fit values of $\sigma(gg\to H\to ZZ)$, $\sigma_i/\sigma_{gg\mathrm{F}}$, and $\mathrm{B}^f/\mathrm{B}^{ZZ}$ relative to their SM prediction from the combined analysis of the $\sqrt{s}$=7 and 8 TeV data. The results are shown for the combination of ATLAS and CMS, and also separately for each experiment, together with their total uncertainties and their breakdown into the four components described in the text. The expected uncertainties in the measurements are also shown.

Best fit values of $\sigma(gg\to H\to WW)$, $\sigma_i/\sigma_{gg\mathrm{F}}$, and $\mathrm{B}^f/\mathrm{B}^{WW}$ from the combined analysis of the $\sqrt{s}$=7 and 8 TeV data. The values involving cross sections are given for $\sqrt{s}$=8 TeV, assuming the SM values for $\sigma_i(7 \mathrm{TeV})/\sigma_i(8 \mathrm{TeV})$. The results are shown for the combination of ATLAS and CMS, and also separately for each experiment, together with their total uncertainties and their breakdown into the four components described in the text. The expected uncertainties in the measurements are also shown.

Best fit values of $\sigma(gg\to H\to WW)$, $\sigma_i/\sigma_{gg\mathrm{F}}$, and $\mathrm{B}^f/\mathrm{B}^{WW}$ relative to their SM prediction from the combined analysis of the $\sqrt{s}$=7 and 8 TeV data. The results are shown for the combination of ATLAS and CMS, and also separately for each experiment, together with their total uncertainties and their breakdown into the four components described in the text. The expected uncertainties in the measurements are also shown.

Best fit values of $\kappa_{gZ}= \kappa_{g}\cdot\kappa_{Z} / \kappa_{H}$ and of the ratios of coupling modifiers, as defined in the most generic parameterisation described in the context of the $\kappa$ framework, from the combined analysis of the $\sqrt{s}$=7 and 8 TeV data. The results are shown for the combination of ATLAS and CMS and also separately for each experiment, together with their total uncertainties and their breakdown into the four components described in the text. The uncertainties in $\lambda_{tg}$ and $\lambda_{WZ}$, for which a negative solution is allowed, are calculated around the overall best fit value. The expected total uncertainties in the measurements are also shown. For those parameters with no sensitivity to the sign, only the absolute values are shown.

Measured global signal strength $\mu$ and its total uncertainty, together with the breakdown of the uncertainty into its four components as defined in the text. The results are shown for the combination of ATLAS and CMS, and separately for each experiment. The expected uncertainty, with its breakdown, is also shown.

Measured signal strengths $\mu$ and their total uncertainties for different Higgs boson production processes. The results are shown for the combination of ATLAS and CMS, and separately for each experiment, for the combined $\sqrt{s}$=7 and 8 TeV data. The expected uncertainties in the measurements are also displayed. These results are obtained assuming that the Higgs boson branching fractions are the same as in the SM.

Measured signal strengths $\mu$ and their total uncertainties for different Higgs boson decay channels. The results are shown for the combination of ATLAS and CMS, and separately for each experiment, for the combined $\sqrt{s}$=7 and 8 TeV data. The expected uncertainties in the measurements are also displayed. These results are obtained assuming that the Higgs boson production process cross sections at $\sqrt{s}$ = 7 and 8 TeV are the same as in the SM.

Results of the ten-parameter fit of $\mu_F^f = \mu_{gg\mathrm{F}+ttH}^f$ and $\mu_V^f = \mu_{\mathrm{VBF}+VH}^f$ for each of the five decay channels. The results are shown for the combination of ATLAS and CMS, together with their measured and expected uncertainties. The measured results are also shown separately for each experiment.

Results of the six-parameter fit of the global ratio $\mu_V/\mu_F = \mu_{\mathrm{VBF}+VH}/\mu_{gg\mathrm{F}+ttH}$ together with $\mu_F^f$ for each of the five decay channels. The results are shown for the combination of ATLAS and CMS, together with their measured and expected uncertainties. The measured results are also shown separately for each experiment.

Fit results allowing BSM loop couplings, assuming that $|\kappa_{V}| \le$ 1, where $\kappa_{V}$ denotes $\kappa_{Z}$ or $\kappa_{W}$, and that $\mathrm{B_{BSM}} \ge$ 0. The results for the combination of ATLAS and CMS are reported with their measured and expected uncertainties. Also shown are the results from each experiment. For the parameters with both signs allowed, the other 1$\sigma$ interval is shown on a second line. For those parameters with no sensitivity to the sign, only the absolute values are shown.

Fit results allowing BSM loop couplings, assuming that there are no additional BSM contributions to the Higgs boson width, i.e. $\mathrm{B_{BSM}}=0$. The results for the combination of ATLAS and CMS are reported with their measured and expected uncertainties. Also shown are the results from each experiment. For the parameters with both signs allowed, the other 1$\sigma$ interval is shown on a second line. When a parameter is constrained and reaches a boundary, namely $\mathrm{B_{BSM}}$ = 0, the uncertainty is not indicated. For those parameters with no sensitivity to the sign, only the absolute values are shown.

Fit results for the parameterisation assuming the absence of BSM particles in the loops ($\mathrm{B_{BSM}}$=0). The results with their measured and expected uncertainties are reported for the combination of ATLAS and CMS, together with the individual results from each experiment. For the parameters with both signs allowed, the other 1$\sigma$ CL interval is shown on a second line. When a parameter is constrained and reaches a boundary, namely $|\kappa_\mu|$ = 0, the uncertainty is not indicated. For those parameters with no sensitivity to the sign, only the absolute values are shown.

Summary of fit results for a parameterisation probing the ratios of coupling modifiers for up-type versus down-type fermions. The results for the combination of ATLAS and CMS are reported together with their measured and expected uncertainties. Also shown are the results from each experiment. The parameter $\kappa_{uu}$ is positive definite since $\kappa_H$ is always assumed to be positive. For the parameter $\lambda_{du}$, for which both signs are allowed, the other 1$\sigma$ CL interval is shown on a second line. Negative values for the parameter $\lambda_{Vu}$ are excluded by more than 4$\sigma$.

Summary of fit results for a parameterisation probing the ratios of coupling modifiers for leptons versus quarks. The results for the combination of ATLAS and CMS are reported together with their measured and expected uncertainties. Also shown are the results from each experiment. The parameter $\kappa_{qq}$ is positive definite since $\kappa_H$ is always assumed to be positive. For the parameter $\lambda_{lq}$, for which there is no sensitivity to the sign, only the absolute values are shown. Negative values for the parameter $\lambda_{Vq}$ are excluded by more than 4$\sigma$.

Correlation matrix obtained from the fit combining the ATLAS and CMS data using the generic parameterisation with 23 parameters, $\sigma_i\cdot\mathrm{B}^f$ for each specific channel $i\to H\to f$. Only 20 parameters are shown because the other three, corresponding to the $H\to ZZ$ decay channel for the $WH$, $ZH$, and $ttH$ production processes, are not measured with a meaningful precision.

Correlation matrix obtained from the fit combining the ATLAS and CMS data using the generic parameterisation with nine parameters, $\sigma(gg\to H\to ZZ)$, $\sigma_i/\sigma_{gg\mathrm{F}}$, and $\mathrm{B}^f/\mathrm{B}^{ZZ}$.

Correlation matrix obtained from the fit combining the ATLAS and CMS data using the generic parameterisation with seven parameters, $\kappa_{gZ}=\kappa_{g}\cdot\kappa_{Z} / \kappa_{H}$ and the ratios of coupling modifiers.

Expected correlation matrix obtained from the fit combining ATLAS and CMS pre-fit Asimov data sets using the generic parameterisation with 23 parameters, $\sigma_i\cdot\mathrm{B}^f$ for each specific channel $i\to H\to f$. Only 20 parameters are shown because the other three, corresponding to the $H\to ZZ$ decay channel for the $WH$, $ZH$, and $ttH$ production processes, are not measured with a meaningful precision.

Expected correlation matrix obtained from the fit combining ATLAS and CMS pre-fit Asimov data sets using the generic parameterisation with nine parameters, $\sigma(gg\to H\to ZZ)$, $\sigma_i/\sigma_{gg\mathrm{F}}$, and $\mathrm{B}^f/\mathrm{B}^{ZZ}$.

Expected correlation matrix obtained from the fit combining ATLAS and CMS pre-fit Asimov data sets using the generic parameterisation with seven parameters, $\kappa_{gZ}=\kappa_{g}\cdot\kappa_{Z} / \kappa_{H}$ and the ratios of coupling modifiers.

ATLAS+CMS combined best-fit and 68% and 95% CL negative log-likelihood contours in the ($\kappa_{g}$,$\kappa_{\gamma}$) plane.

ATLAS best-fit and 68% and 95% CL negative log-likelihood contours in the ($\kappa_{g}$,$\kappa_{\gamma}$) plane.

CMS best-fit and 68% and 95% CL negative log-likelihood contours in the ($\kappa_{g}$,$\kappa_{\gamma}$) plane.

Best-fit and 68% and 95% CL negative log-likelihood contours in the ($\kappa_{F}$,$\kappa_{V}$) plane.

Best-fit and 68% CL negative log-likelihood contour in the ($\kappa_{F}$,$\kappa_{V}$) plane for the $H\rightarrow\gamma\gamma$ channel.

Best-fit and 68% CL negative log-likelihood contour in the ($\kappa_{F}$,$\kappa_{V}$) plane for the $H\rightarrow ZZ$ channel.

Best-fit and 68% CL negative log-likelihood contour in the ($\kappa_{F}$,$\kappa_{V}$) plane for the $H\rightarrow WW$ channel.

Best-fit and 68% CL negative log-likelihood contour in the ($\kappa_{F}$,$\kappa_{V}$) plane for the $H\rightarrow\tau\tau$ channel.

Best-fit and 68% CL negative log-likelihood contour in the ($\kappa_{F}$,$\kappa_{V}$) plane for the $H\rightarrow bb$ channel.

ATLAS best-fit and 68% and 95% CL negative log-likelihood contours in the ($\kappa_{F}$,$\kappa_{V}$) plane.

CMS best-fit and 68% and 95% CL negative log-likelihood contours in the ($\kappa_{F}$,$\kappa_{V}$) plane.

Best-fit and 68% CL negative log-likelihood contour in the ($\kappa_{F}$,$\kappa_{V}$) plane for the $H\rightarrow\gamma\gamma$ channel, assuming $\kappa_{F}>0$.

Best-fit and 68% CL negative log-likelihood contour in the ($\kappa_{F}$,$\kappa_{V}$) plane for the $H\rightarrow ZZ$ channel, assuming $\kappa_{F}>0$.

Best-fit and 68% CL negative log-likelihood contour in the ($\kappa_{F}$,$\kappa_{V}$) plane for the $H\rightarrow WW$ channel, assuming $\kappa_{F}>0$.

Best-fit and 68% CL negative log-likelihood contour in the ($\kappa_{F}$,$\kappa_{V}$) plane for the $H\rightarrow\tau\tau$ channel, assuming $\kappa_{F}>0$.

Best-fit and 68% CL negative log-likelihood contour in the ($\kappa_{F}$,$\kappa_{V}$) plane for the $H\rightarrow bb$ channel, assuming $\kappa_{F}>0$.

Best-fit and 68% and 95% CL negative log-likelihood contours in the ($\mu^{f}_{gg\mathrm{F}+ttH}$,$\mu^{f}_{\mathrm{VBF}+VH}$) plane for the $H\rightarrow\gamma\gamma$ channel.

Best-fit and 68% and 95% CL negative log-likelihood contours in the ($\mu^{f}_{gg\mathrm{F}+ttH}$,$\mu^{f}_{\mathrm{VBF}+VH}$) plane for the $H\rightarrow ZZ$ channel.

Best-fit and 68% and 95% CL negative log-likelihood contours in the ($\mu^{f}_{gg\mathrm{F}+ttH}$,$\mu^{f}_{\mathrm{VBF}+VH}$) plane for the $H\rightarrow WW$ channel.

Best-fit and 68% and 95% CL negative log-likelihood contours in the ($\mu^{f}_{gg\mathrm{F}+ttH}$,$\mu^{f}_{\mathrm{VBF}+VH}$) plane for the $H\rightarrow\tau\tau$ channel.

Best-fit and 68% and 95% CL negative log-likelihood contours in the ($\mu^{f}_{gg\mathrm{F}+ttH}$,$\mu^{f}_{\mathrm{VBF}+VH}$) plane for the $H\rightarrow bb$ channel.

Expected 95% CL upper limits with 1-sigma uncertainty on the production of a heavy Higgs boson as functions of fw and fww with mass equal to 300 GeV.

Expected 95% CL upper limits with 2-sigma uncertainty on the production of a heavy Higgs boson as functions of fw and fww with mass equal to 300 GeV.

Observed 95% CL upper limits on the production of a heavy Higgs boson as functions of fw and fww with mass equal to 300 GeV.

Expected 95% CL upper limits on the production of a heavy Higgs boson as functions of fw and fww with mass equal to 600 GeV.

Expected 95% CL upper limits with 1-sigma uncertainty on the production of a heavy Higgs boson as functions of fw and fww with mass equal to 600 GeV.

Expected 95% CL upper limits with 2-sigma uncertainty on the production of a heavy Higgs boson as functions of fw and fww with mass equal to 600 GeV.

Observed 95% CL upper limits on the production of a heavy Higgs boson as functions of fw and fww with mass equal to 600 GeV.

Expected 95% CL upper limits on the production of a heavy Higgs boson as functions of fw and fww with mass equal to 900 GeV.

Expected 95% CL upper limits with 1-sigma uncertainty on the production of a heavy Higgs boson as functions of fw and fww with mass equal to 900 GeV.

Expected 95% CL upper limits with 2-sigma uncertainty on the production of a heavy Higgs boson as functions of fw and fww with mass equal to 900 GeV.

Observed 95% CL upper limits on the production of a heavy Higgs boson as functions of fw and fww with mass equal to 900 GeV.

Observed (black solid curve) and expected (black dashed curve) 95% CL upper limits on the production of a heavy Higgs boson as functions of its mass with ( fW, fWW) fixed at (0, 6200). The green (inner) and yellow (outer) bands represent 1 and 2 sigma uncertainty in the expected limits. Theoretical predictions (red solid curve) as functions of the heavy Higgs boson mass are overlaid.

Observed (black solid curve) and expected (black dashed curve) 95% CL upper limits on the production of a heavy Higgs boson as functions of its mass with ( fW, fWW) fixed at (1350, 0). The green (inner) and yellow (outer) bands represent 1 and 2 sigma uncertainty in the expected limits. Theoretical predictions (red solid curve) as functions of the heavy Higgs boson mass are overlaid.

Event yields of the signal and SM background processes in the six analysis regions after the fit to the data under the $t \rightarrow uX$ hypothesis assuming $m_X = 30$ GeV. Total includes signal and background.The quoted uncertainties take into account correlations and constraints of the nuisance parameters and include both the statistical and systematic uncertainties. Negative correlations between the $t\bar{t} +$ light, $t\bar{t} + \geq1b$ and $t\bar{t} + \geq1c$ modelling uncertainties can make the uncertainty in the total yields smaller than in the individual components.

Event yields of the signal and SM background processes in the six analysis regions after the fit to the data under the $t \rightarrow cX$ hypothesis assuming $m_X = 30$ GeV. Total includes signal and background.The quoted uncertainties take into account correlations and constraints of the nuisance parameters and include both the statistical and systematic uncertainties. Negative correlations between the $t\bar{t} +$ light, $t\bar{t} + \geq1b$ and $t\bar{t} + \geq1c$ modelling uncertainties can make the uncertainty in the total yields smaller than in the individual components.

Event acceptance times efficiency in percent for every $t\rightarrow uX$ and $t\rightarrow cX$ mass signal sample.

Cut flow for the scalar signal in the $t\rightarrow uX$ decay combining both quark and anti-quark samples. Shown for each signal are the corresponding mass, the number of generated events, the number of reconstructed ("Reco") events, the events that pass the lepton triggers, the events that have only one electron ("el") or only one muon ("mu") with p$_{\text{T}}$ larger than 27 GeV, the number of events with at least four jets with p$_{\text{T}}$ larger than 25 GeV, the number of events with at least three $b$-tagged jets at the 70% efficiency working point, and the number of events with at least two $b$-tagged jets at the 60% efficiency working point and at least another one at the 70%. The quoted yields do not include reweighting.

Cut flow for the scalar signal in the $t\rightarrow cX$ decay combining both quark and anti-quark samples. Shown for each signal are the corresponding mass, the number of generated events, the number of reconstructed ("Reco") events, the events that pass the lepton triggers, the events that have only one electron ("el") or only one muon ("mu") with p$_{\text{T}}$ larger than 27 GeV, the number of events with at least four jets with p$_{\text{T}}$ larger than 25 GeV, the number of events with at least three $b$-tagged jets at the 70% efficiency working point, and the number of events with at least two $b$-tagged jets at the 60% efficiency working point and at least another one at the 70%. The quoted yields do not include reweighting.

Probability of finding at least one TAR jet, where the p_T-leading TAR jet passes the m_Wcand and D₂^β=1 requirements, as a function of m_s. The probability is determined in a sample of signal events with m_Z'=2100 GeV, with the preselections applied.

Probability of finding a W_had candidate reconstructed as a pair of R=0.4 PFlow jets, as a function of m_s. The probability is determined in a sample of signal events with m_Z'=500 GeV, with the preselections applied that do not pass the requirements of the merged category.

Probability of finding a W_had candidate reconstructed as a pair of R=0.4 PFlow jets, as a function of m_s. The probability is determined in a sample of signal events with m_Z'=1000 GeV, with the preselections applied that do not pass the requirements of the merged category.

Probability of finding a W_had candidate reconstructed as a pair of R=0.4 PFlow jets, as a function of m_s. The probability is determined in a sample of signal events with m_Z'=1700 GeV, with the preselections applied that do not pass the requirements of the merged category.

Probability of finding a W_had candidate reconstructed as a pair of R=0.4 PFlow jets, as a function of m_s. The probability is determined in a sample of signal events with m_Z'=2100 GeV, with the preselections applied that do not pass the requirements of the merged category.

Observed exclusion contour at 95% C.L. for the dark Higgs model in the (m_Z', m_s) plane for g_q=0.25, g_χ=1, m_χ=200 GeV, and sinθ=0.01.

Expected exclusion contour at 95% C.L. for the dark Higgs model in the (m_Z', m_s) plane for g_q=0.25, g_χ=1, m_χ=200 GeV, and sinθ=0.01.

Expected+1σ exclusion contour at 95% C.L. for the dark Higgs model in the (m_Z', m_s) plane for g_q=0.25, g_χ=1, m_χ=200 GeV, and sinθ=0.01.

Expected-1σ exclusion contour at 95% C.L. for the dark Higgs model in the (m_Z', m_s) plane for g_q=0.25, g_χ=1, m_χ=200 GeV, and sinθ=0.01.

Expected+2σ exclusion contour at 95% C.L. for the dark Higgs model in the (m_Z', m_s) plane for g_q=0.25, g_χ=1, m_χ=200 GeV, and sinθ=0.01.

Expected-2σ exclusion contour at 95% C.L. for the dark Higgs model in the (m_Z', m_s) plane for g_q=0.25, g_χ=1, m_χ=200 GeV, and sinθ=0.01.

Observed upper limits at 95% C.L. on σ(pp → s χχ) × B(s → W^± W^∓) for m_Z'=0.5 TeV as a function of m_s. The expected limits, varied up and down by one and two standard deviations, are shown as green and yellow bands, respectively. The observed and expected limits are compared to the theoretical LO cross section for the σ(pp → s χχ) × B(s → W^± W^∓) process for m_Z'=0.5 TeV, shown in dashed blue.

Observed upper limits at 95% C.L. on σ(pp → s χχ) × B(s → W^± W^∓) for m_Z'=1 TeV as a function of m_s. The expected limits, varied up and down by one and two standard deviations, are shown as green and yellow bands, respectively. The observed and expected limits are compared to the theoretical LO cross section for the σ(pp → s χχ) × B(s → W^± W^∓) process for m_Z'=1 TeV, shown in dashed blue.

Observed upper limits at 95% C.L. on σ(pp → s χχ) × B(s → W^± W^∓) for m_Z'=1.7 TeV as a function of m_s. The expected limits, varied up and down by one and two standard deviations, are shown as green and yellow bands, respectively. The observed and expected limits are compared to the theoretical LO cross section for the σ(pp → s χχ) × B(s → W^± W^∓) process for m_Z'=1.7 TeV, shown in dashed blue.

Observed upper limits at 95% C.L. on σ(pp → s χχ) × B(s → W^± W^∓) for m_Z'=2.1 TeV as a function of m_s. The expected limits, varied up and down by one and two standard deviations, are shown as green and yellow bands, respectively. The observed and expected limits are compared to the theoretical LO cross section for the σ(pp → s χχ) × B(s → W^± W^∓) process for m_Z'=2.1 TeV, shown in dashed blue.

Data overlaid on SM background yields stacked in each SR and CR category after the fit to data ('Post-fit'). The yields in the SR are broken down into their contributions to the individual bins. The maximum-likelihood estimators are set to the conditional values of the CR-only fit, and propagated to SR and CRs.

Dominant sources of uncertainty for three dark Higgs scenarios after the fit to data. The uncertainties are quantified in terms of their contribution to the fitted signal uncertainty that is expressed relative to the theory prediction. Three representative dark Higgs signal scenarios with g_q=0.25, g_χ=1.0, sinθ=0.01 and m_χ=200 GeV are considered, which are indicated using the (m_Z', m_s) format in units of GeV in the table columns.

Cumulative efficiencies in the merged category for three representative dark Higgs signal scenarios with g_q=0.25, g_χ=1.0, sinθ=0.01, m_Z' = 1 TeV, and m_χ=200 GeV considering s→W(ℓν)W(qq) decays only.

Cumulative efficiencies in the resolved category for three representative dark Higgs signal scenarios with g_q=0.25, g_χ=1.0, sinθ=0.01, m_Z' = 1 TeV, and m_χ=200 GeV considering s→W(ℓν)W(qq) decays only.

Theoretical cross section for &sigma;(pp &rarr; s&chi;&chi;) &times; B(s &rarr; W^&plusmn;W^&#8723;) for each of the dark Higgs signal points at m_Z′ ={300, 350, 400, 500, 750} GeV, with g_q = 0.25, g<sub>&chi; = 1.0, sin&theta; = 0.01, m_Z′ = 1 TeV , and m_&chi; = 200 GeV. Also shown are experimentally excluded cross sections of &sigma;(pp &rarr; s&chi;&chi;) &times; B(s &rarr; W^&plusmn;W^&#8723;) (Obs.) together with the expectations (Exp.) varied up and down by one standard deviation (&plusmn;1&sigma;).

Theoretical cross section for &sigma;(pp &rarr; s&chi;&chi;) &times; B(s &rarr; W^&plusmn;W^&#8723;) for each of the dark Higgs signal points at m_Z′ ={1000, 1700} GeV, with g_q = 0.25, g<sub>&chi; = 1.0, sin&theta; = 0.01, m_Z′ = 1 TeV , and m_&chi; = 200 GeV. Also shown are experimentally excluded cross sections of &sigma;(pp &rarr; s&chi;&chi;) &times; B(s &rarr; W^&plusmn;W^&#8723;) (Obs.) together with the expectations (Exp.) varied up and down by one standard deviation (&plusmn;1&sigma;).

Theoretical cross section for &sigma;(pp &rarr; s&chi;&chi;) &times; B(s &rarr; W^&plusmn;W^&#8723;) for each of the dark Higgs signal points at m_Z′ ={2100, 2500, 2900, 3300} GeV, with g_q = 0.25, g<sub>&chi; = 1.0, sin&theta; = 0.01, m_Z′ = 1 TeV , and m_&chi; = 200 GeV. Also shown are experimentally excluded cross sections of &sigma;(pp &rarr; s&chi;&chi;) &times; B(s &rarr; W^&plusmn;W^&#8723;) (Obs.) together with the expectations (Exp.) varied up and down by one standard deviation (&plusmn;1&sigma;).

The $m_{jj}$ distributions for the CRs and the BDT classifier response distribution for the SR after the fit in all regions. The dashed line shows the total background distribution before the fit. The vertical error bars on the data points correspond to the data's statistical uncertainty. Overflows are included in the last bin. The uncertainty band corresponds to the combination of the MC statistical uncertainty and systematic uncertainties obtained in the fit.

The $m_{jj}$ distributions for the CRs and the BDT classifier response distribution for the SR after the fit in all regions. The dashed line shows the total background distribution before the fit. The vertical error bars on the data points correspond to the data's statistical uncertainty. Overflows are included in the last bin. The uncertainty band corresponds to the combination of the MC statistical uncertainty and systematic uncertainties obtained in the fit.

The $m_{jj}$ distributions for the CRs and the BDT classifier response distribution for the SR after the fit in all regions. The dashed line shows the total background distribution before the fit. The vertical error bars on the data points correspond to the data's statistical uncertainty. Overflows are included in the last bin. The uncertainty band corresponds to the combination of the MC statistical uncertainty and systematic uncertainties obtained in the fit.

Observed and expected event yields for the signal and all of the background processes considered in this analysis after the fit to the data in all of the regions. The uncertainty on the expected yield is the combination of statistical and systematic uncertainties obtained in the fit. The individual uncertainties can be correlated and do not necessarily add in quadrature to equal the total background uncertainty.

Summary of the event yield for processes in all regions, after the fit over all regions. The dashed line shows the total background distribution before the fit. The vertical error bars on the data points correspond to the data's statistical uncertainty. The uncertainty band corresponds to the combination of the MC statistical uncertainty and systematic uncertainties obtained in the fit.

Impact of different components of systematic uncertainty on the measured cross section, without taking into account the correlations. The impact is calculated by fixing the value of the corresponding nuisance parameters to the values obtained in the fit used to measure the cross section, performing the fit, estimating the signal strength uncertainty and subtracting it from the nominal uncertainty in quadrature.

Fitted POI values for this analysis, the previous ATLAS analysis, and their combination. The first and second columns present the values obtained in the individual analyses. The third column presents the values obtained in the combination.

Observed and expected one-dimensional limits on dimension 8 aQGC parameters. Limits are obtained by setting all aQGCs parameters except one to zero. Unitarity is not preserved.

Observed and expected one-dimensional limits on dimension 8 aQGC parameters in the region, where unitarity is preserved. Cutoff scales in MC, $E_{c}$, are given for each parameter. Limits are obtained by setting all aQGCs parameters except one to zero.

The $E_{T}^{\gamma}$ distribution in the SR after the fit in the control regions. The red (green) line shows the expected number of events in the case of non-zero EFT coefficient $f_{T0}/\Lambda^4$ ($f_{M0}/\Lambda^4$) with the value shown in the legend. The vertical error bars on the data points correspond to the data statistical uncertainty. Overflows are included in the last bin. The uncertainty band corresponds to the combination of the MC statistical uncertainty and systematic uncertainties obtained in the fit.

Evolution of the expected (red line) and observed (blue line) limits versus $E_{c}$ values for $f_{T0}/\Lambda^4$. The unitarity bound is shown by the black line. The $E_{c}<4$ TeV regime was obtained with $E_{T}^{\gamma} > 600$ GeV. The $E_{c} \geq 4$ TeV regime was obtained with $E_{T}^{\gamma} > 900$ GeV.

Evolution of the expected (red line) and observed (blue line) limits versus $E_{c}$ values for $f_{T5}/\Lambda^4$. The unitarity bound is shown by the black line. The $E_{c}<4$ TeV regime was obtained with $E_{T}^{\gamma} > 600$ GeV. The $E_{c} \geq 4$ TeV regime was obtained with $E_{T}^{\gamma} > 900$ GeV.

Evolution of the expected (red line) and observed (blue line) limits versus $E_{c}$ values for $f_{T8}/\Lambda^4$. The unitarity bound is shown by the black line. The $E_{c}<4$ TeV regime was obtained with $E_{T}^{\gamma} > 600$ GeV. The $E_{c} \geq 4$ TeV regime was obtained with $E_{T}^{\gamma} > 900$ GeV.

Evolution of the expected (red line) and observed (blue line) limits versus $E_{c}$ values for $f_{T9}/\Lambda^4$. The unitarity bound is shown by the black line. The $E_{c}<4$ TeV regime was obtained with $E_{T}^{\gamma} > 600$ GeV. The $E_{c} \geq 4$ TeV regime was obtained with $E_{T}^{\gamma} > 900$ GeV.

Evolution of the expected (red line) and observed (blue line) limits versus $E_{c}$ values for $f_{M0}/\Lambda^4$. The unitarity bound is shown by the black line. The $E_{c}<4$ TeV regime was obtained with $E_{T}^{\gamma} > 400$ GeV. The $E_{c} \geq 4$ TeV regime was obtained with $E_{T}^{\gamma} > 900$ GeV.

Evolution of the expected (red line) and observed (blue line) limits versus $E_{c}$ values for $f_{M1}/\Lambda^4$. The unitarity bound is shown by the black line. The $E_{c}<4$ TeV regime was obtained with $E_{T}^{\gamma} > 400$ GeV. The $E_{c} \geq 4$ TeV regime was obtained with $E_{T}^{\gamma} > 900$ GeV.

Evolution of the expected (red line) and observed (blue line) limits versus $E_{c}$ values for $f_{M2}/\Lambda^4$. The unitarity bound is shown by the black line. The $E_{c}<4$ TeV regime was obtained with $E_{T}^{\gamma} > 400$ GeV. The $E_{c} \geq 4$ TeV regime was obtained with $E_{T}^{\gamma} > 900$ GeV.

Correlation matrix of the parameters for the fit in all regions. Only parameters with absolute value of the correlation coefficient larger than 10 are shown.

Results of applying the BH algorithm to the mass spectra in the leading photon category, using the fitted background estimations from the initial step

Results of applying the BH algorithm to the mass spectra in the leading electron category, using the fitted background estimations from the initial step

Results of applying the BH algorithm to the mass spectra in the leading muon category, using the fitted background estimations from the initial step

Results of applying the BH algorithm to the mass spectra in the leading small-R jet category, using the fitted background estimations from the initial step

Results of applying the BH algorithm to the mass spectra in the leading bjet category, using the fitted background estimations from the initial step

Results of applying the BH algorithm to the mass spectra in the leading large-R jet category, using the fitted background estimations from the initial step

Results of applying the BH algorithm to the mass spectra in the leading photon category, using the fitted background estimations from the initial step

Results of applying the BH algorithm to the mass spectra in the leading electron category, using the fitted background estimations from the initial step

Results of applying the BH algorithm to the mass spectra in the leading muon category, using the fitted background estimations from the initial step

Results of applying the BH algorithm to the mass spectra in the inclusive category, using the fitted background estimations from the initial step

Results of applying the BH algorithm to the mass spectra in the inclusive category, using the fitted background estimations from the initial step

Upper limits at 95% CL on the cross section times branching fraction times acceptance for a Gaussian-shaped signal with a relative width value of 3% as a function of mass in the leading small-R jet category

Upper limits at 95% CL on the cross section times branching fraction times acceptance for a Gaussian-shaped signal with a relative width value of 3% as a function of mass in the leading bjet category

Upper limits at 95% CL on the cross section times branching fraction times acceptance for a Gaussian-shaped signal with a relative width value of 3% as a function of mass in the leading large-R jet category

Upper limits at 95% CL on the cross section times branching fraction times acceptance for a Gaussian-shaped signal with a relative width value of 3% as a function of mass in the leading photon category

Upper limits at 95% CL on the cross section times branching fraction times acceptance for a Gaussian-shaped signal with a relative width value of 3% as a function of mass in the leading electron category

Upper limits at 95% CL on the cross section times branching fraction times acceptance for a Gaussian-shaped signal with a relative width value of 3% as a function of mass in the leading muon category

Upper limits at 95% CL on the cross section times branching fraction times acceptance for a Gaussian-shaped signal with a relative width value of 3% as a function of mass in the leading small-R jet category

Upper limits at 95% CL on the cross section times branching fraction times acceptance for a Gaussian-shaped signal with a relative width value of 3% as a function of mass in the leading bjet category

Upper limits at 95% CL on the cross section times branching fraction times acceptance for a Gaussian-shaped signal with a relative width value of 3% as a function of mass in the leading large-R jet category

Upper limits at 95% CL on the cross section times branching fraction times acceptance for a Gaussian-shaped signal with a relative width value of 3% as a function of mass in the leading photon category

Upper limits at 95% CL on the cross section times branching fraction times acceptance for a Gaussian-shaped signal with a relative width value of 3% as a function of mass in the leading electron category

Upper limits at 95% CL on the cross section times branching fraction times acceptance for a Gaussian-shaped signal with a relative width value of 3% as a function of mass in the leading muon category

Upper limits at 95% CL on the cross section times branching fraction times acceptance for a Gaussian-shaped signal with a relative width value of 3% as a function of mass in the inclusive category

Upper limits at 95% CL on the cross section times branching fraction times acceptance for a Gaussian-shaped signal with a relative width value of 3% as a function of mass in the inclusive category

Comparison of observed upper limits at 95% CL on the cross section times branching fraction times acceptance for a Gaussian-shaped signals with relative width values as a function of mass in the leading small-R jet category

Comparison of observed upper limits at 95% CL on the cross section times branching fraction times acceptance for a Gaussian-shaped signals with relative width values as a function of mass in the leading bjet category

Comparison of observed upper limits at 95% CL on the cross section times branching fraction times acceptance for a Gaussian-shaped signals with relative width values as a function of mass in the leading large-R jet category

Comparison of observed upper limits at 95% CL on the cross section times branching fraction times acceptance for a Gaussian-shaped signals with relative width values as a function of mass in the leading photon category

Comparison of observed upper limits at 95% CL on the cross section times branching fraction times acceptance for a Gaussian-shaped signals with relative width values as a function of mass in the leading electron category

Comparison of observed upper limits at 95% CL on the cross section times branching fraction times acceptance for a Gaussian-shaped signals with relative width values as a function of mass in the leading muon category

Comparison of observed upper limits at 95% CL on the cross section times branching fraction times acceptance for a Gaussian-shaped signals with relative width values as a function of mass in the leading small-R jet category

Comparison of observed upper limits at 95% CL on the cross section times branching fraction times acceptance for a Gaussian-shaped signals with relative width values as a function of mass in the leading bjet category

Comparison of observed upper limits at 95% CL on the cross section times branching fraction times acceptance for a Gaussian-shaped signals with relative width values as a function of mass in the leading large-R jet category

Comparison of observed upper limits at 95% CL on the cross section times branching fraction times acceptance for a Gaussian-shaped signals with relative width values as a function of mass in the leading photon category

Comparison of observed upper limits at 95% CL on the cross section times branching fraction times acceptance for a Gaussian-shaped signals with relative width values as a function of mass in the leading electron category

Comparison of observed upper limits at 95% CL on the cross section times branching fraction times acceptance for a Gaussian-shaped signals with relative width values as a function of mass in the leading muon category

Comparison of observed upper limits at 95% CL on the cross section times branching fraction times acceptance for a Gaussian-shaped signals with relative width values as a function of mass in the inclusive category

Comparison of observed upper limits at 95% CL on the cross section times branching fraction times acceptance for a Gaussian-shaped signals with relative width values as a function of mass in the inclusive category

Acceptance times efficiency in an HVT model as a function of mass in the leading large-R jet category

Upper limits at 95% CL on the cross section times the branching fraction($W^{\prime} \to ZW$) for the HVT signal as functions of mass($m_{ZX}$) in the leading large-R jet category.The full(dashed) red curves correspond to the theoretical predictions of HVT model A(B)

Comparison of the identification efficiencies using standard and merged-ee reconstruction as a function of true PT(Z)

Background rejection factor as a function of signal efficiency. The red curve shows the BDT performance whereas the blue curve corresponds to that of a cut-based analysis relying only on the jet energyfraction deposited in the EM calorimeter

BH p-values of the 100 pseudo-experiments as a function of in the leading small-R jet category for an injected Gaussian-shaped signal with a relative width value of 3%. The fractions of the pseudo-experiments that have the correctly identified interval and BH p-values below the threshold of 0.01 are indicated. The background is derived by the background-only fit in the full fitting range.

BH p-values of the 100 pseudo-experiments as a function of in the leading small-R jet category for an injected Gaussian-shaped signal with a relative width value of 3%. The fractions of the pseudo-experiments that have the correctly identified interval and BH p-values below the threshold of 0.01 are indicated. The background is derived by the background-only fit in the range excluding the BH interval.

Fractions of pseudo-experiments in which the detected BH interval agrees with the injected mass point and the BH p-value is below 0.01 as a function of mass in the leading small-R jet category for Gaussian-shaped signal with a relative width of 3%.

Fractions of pseudo-experiments in which the detected BH interval agrees with the injected mass point and the BH p-value is below 0.01 as a function of mass in the leading small-R jet category for Gaussian-shaped signal with a relative width of 3%.

Distribution of exclusion upper limits on signal event yields at 95% CL from 1000 pseudo-experiments for Gaussian-shaped signals with relative width values of 3% at the low boundary of the limit-sensitive mass range for the ZX spectrum of the leading small-R jet category. The vertical line corresponds to the expected nominal exclusion limit at 95% CL.

Distribution of exclusion upper limits on signal event yields at 95% CL from 1000 pseudo-experiments for Gaussian-shaped signals with relative width values of 3% at the high boundary of the limit-sensitive mass range for the ZX spectrum of the leading small-R jet category. The vertical line corresponds to the expected nominal exclusion limit at 95% CL.

Distribution of exclusion upper limits on signal event yields at 95% CL from 1000 pseudo-experiments for Gaussian-shaped signals with relative width values of 5% at the low boundary of the limit-sensitive mass range for the ZX spectrum of the leading small-R jet category. The vertical line corresponds to the expected nominal exclusion limit at 95% CL.

Distribution of exclusion upper limits on signal event yields at 95% CL from 1000 pseudo-experiments for Gaussian-shaped signals with relative width values of 5% at the high boundary of the limit-sensitive mass range for the ZX spectrum of the leading small-R jet category. The vertical line corresponds to the expected nominal exclusion limit at 95% CL.

Distribution of exclusion upper limits on signal event yields at 95% CL from 1000 pseudo-experiments for Gaussian-shaped signals with relative width values of 10% at the low boundary of the limit-sensitive mass range for the ZX spectrum of the leading small-R jet category. The vertical line corresponds to the expected nominal exclusion limit at 95% CL.

Distribution of exclusion upper limits on signal event yields at 95% CL from 1000 pseudo-experiments for Gaussian-shaped signals with relative width values of 10% at the high boundary of the limit-sensitive mass range for the ZX spectrum of the leading small-R jet category. The vertical line corresponds to the expected nominal exclusion limit at 95% CL.

Data yields of the six event categories in the $Z\to e^+e^-$ and $\mu^+\mu^-$ decay channels. The merged-$e^+e^-$ events are included in the $e^+e^-$ channel, increasing the event yield, mainly in the leading large-$R$-jet category, by 0.6 %.

A list of mass spectra, event categories and their corresponding fit ranges, functional forms, numbers of free parameters and global $\chi^2$ $p$-values from background-only fits. The fit range values are rounded to the nearest 5 GeV. Here $f_1(x)=p_0\left(\mathrm{e}^{-p_1x}+p_2\mathrm{e}^{-(p_1+p_3)x}+p_4\mathrm{e}^{-(p_1+p_3+p_5)x}+\cdots\right)$ and $f_2(x)=p_0(1-x)^{p_1}x^{p_2+p_3\ln\!x+p_4\ln^2\!x+\cdots}$

A list of mass spectra, event categories and their corresponding signal-sensitive mass ranges. The initial BH $p$-value is obtained by using the background derived from the background-only fit in the full fit range, whereas the new BH $p$-value uses the background derived from the background-only fit in the range excluding the initial BH interval.

A list of mass spectra, event categories, relative width values of Gaussian-shaped signals and limit-sensitive mass ranges and fractions, corresponding to the mass values in the previous column, of pseudo-experiments having exclusion upper limits higher than the nominal exclusion limit at 95% CL

Cutflow of HVT model $A$ signals ($W^\prime \to ZW \to \ell\ell qq$) with $m_{W^\prime} = 1$ TeV and $m_{W^\prime} = 4$ TeV based on MC simulations.

Acceptance times efficiency ($\mathcal{A} \times \epsilon$) values in % in the dominant event category for $p^Z_{\mathrm{T}} > 100$ GeV in the $Z \to \ell^{+}\ell^{-}$ decay channel.