Distribution of the four-lepton invariant mass in the four-lepton final state for the ggF-MVA-low category.

Distribution of the four-lepton invariant mass in the four-lepton final state for the VBF-MVA-enriched category.

Distribution of the transverse mass in the llvv final state for the ggF-enriched eevv category.

Distribution of the transverse mass in the llvv final state for the ggF-enriched mumuvv category.

Distribution of the transverse mass in the llvv final state for the VBF-enriched eevv category.

Distribution of the transverse mass in the llvv final state for the VBF-enriched mumuvv category.

The upper limits at 95% CL on the cross section times branching ratio as a function of the heavy resonance mass for the ggF production mode

The upper limits at 95% CL on the cross section times branching ratio as a function of the heavy resonance mass for the VBF production mode

The upper limits at 95% CL on the cross section for the ggF production mode times branching ratio as a function of the heavy resonance mass for an additional heavy scalar assuming a width of 1% of its mass

The upper limits at 95% CL on the cross section for the ggF production mode times branching ratio as a function of the heavy resonance mass for an additional heavy scalar assuming a width of 5% of its mass

The upper limits at 95% CL on the cross section for the ggF production mode times branching ratio as a function of the heavy resonance mass for an additional heavy scalar assuming a width of 10% of its mass

The upper limits at 95% CL on the cross section for the ggF production mode times branching ratio as a function of the heavy resonance mass for an additional heavy scalar assuming a width of 15% of its mass

The upper limits at 95% CL on the cross section times branching ratio for a KK graviton produced with k/M_{PI} = 1

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

Summary of the observed and expected 95% CL upper limits on the cross section for $\mathrm{LQ}_{3}^{\mathrm{d}}$ pair production as a function of $m_{\mathrm{LQ}_{3}^{\mathrm{d}}}$ under the assumptions of B=1.

Summary of the observed and expected 95% CL upper limits on B as a function of $m_{\mathrm{LQ}_{3}^{\mathrm{d}}}$.

Summary of the observed and expected 95% CL upper limits on B as a function of $m_{\mathrm{LQ}_{3}^{\mathrm{d}}}$.

Cutflow of the preselection requirements (see Section 5) for $\mathrm{LQ}_{3}^{\mathrm{d}}$ signals with $m_{\mathrm{LQ}_{3}^{\mathrm{d}}}$=0.9, 1.1, and 1.3 TeV, assuming B=1. The yields correspond to an integrated luminosity of 139 fb$^{-1}$.

Cutflow of the preselection requirements (see Section 5) for $\mathrm{LQ}_{3}^{\mathrm{d}}$ signals with $m_{\mathrm{LQ}_{3}^{\mathrm{d}}}$=0.9, 1.1, and 1.3 TeV, assuming B=1. The yields correspond to an integrated luminosity of 139 fb$^{-1}$.

Cutflow of the signal region requirements in 1$\ell+\geq 1\tau$ channel (see Table 3) for $\mathrm{LQ}_{3}^{\mathrm{d}}$ signals with $m_{\mathrm{LQ}_{3}^{\mathrm{d}}}$=0.9, 1.1, and 1.3 TeV, assuming B=1. Events that satisfy the preselection requirements are considered. The yields correspond to an integrated luminosity of 139 fb$^{-1}$.

Cutflow of the signal region requirements in 1$\ell+\geq 1\tau$ channel (see Table 3) for $\mathrm{LQ}_{3}^{\mathrm{d}}$ signals with $m_{\mathrm{LQ}_{3}^{\mathrm{d}}}$=0.9, 1.1, and 1.3 TeV, assuming B=1. Events that satisfy the preselection requirements are considered. The yields correspond to an integrated luminosity of 139 fb$^{-1}$.

Cutflow of the signal region requirements in the 2$\ell$OS+$\geq 1\tau$ channel (see Table 4) for $\mathrm{LQ}_{3}^{\mathrm{d}}$ signals with $m_{\mathrm{LQ}_{3}^{\mathrm{d}}}$=0.9, 1.1, and 1.3 TeV, assuming B=1. Events that satisfy the preselection requirements are considered. The yields correspond to an integrated luminosity of 139 fb$^{-1}$.

Cutflow of the signal region requirements in the 2$\ell$OS+$\geq 1\tau$ channel (see Table 4) for $\mathrm{LQ}_{3}^{\mathrm{d}}$ signals with $m_{\mathrm{LQ}_{3}^{\mathrm{d}}}$=0.9, 1.1, and 1.3 TeV, assuming B=1. Events that satisfy the preselection requirements are considered. The yields correspond to an integrated luminosity of 139 fb$^{-1}$.

Cutflow of the signal region requirements in the $2\ell$SS/$3\ell+\geq 1\tau$ channel (see Table 5) for $\mathrm{LQ}_{3}^{\mathrm{d}}$ signals with $m_{\mathrm{LQ}_{3}^{\mathrm{d}}}$=0.9, 1.1, and 1.3 TeV, assuming B=1. Events that satisfy the preselection requirements are considered. In this channel, two signal regions (SR-L and SR-H) are defined based on $p_{\mathrm{T}, 1}^{\tau}$, with SR-L and SR-H requiring $125< p_{\mathrm{T}, 1}^{\tau} < 225$ GeV and $p_{\mathrm{T}, 1}^{\tau}>225$ GeV, respectively. The yields correspond to an integrated luminosity of 139 fb$^{-1}$.

Cutflow of the signal region requirements in the $2\ell$SS/$3\ell+\geq 1\tau$ channel (see Table 5) for $\mathrm{LQ}_{3}^{\mathrm{d}}$ signals with $m_{\mathrm{LQ}_{3}^{\mathrm{d}}}$=0.9, 1.1, and 1.3 TeV, assuming B=1. Events that satisfy the preselection requirements are considered. In this channel, two signal regions (SR-L and SR-H) are defined based on $p_{\mathrm{T}, 1}^{\tau}$, with SR-L and SR-H requiring $125< p_{\mathrm{T}, 1}^{\tau} < 225$ GeV and $p_{\mathrm{T}, 1}^{\tau}>225$ GeV, respectively. The yields correspond to an integrated luminosity of 139 fb$^{-1}$.

Inclusive differential cross section for two electrons and two muons (Max = 1710~GeV).

Mass of the first dilepton pair in the $Z$ mass region (Max = 88.1~GeV).

Mass of the first dilepton pair in the $H$ mass region (Max = 107~GeV).

Mass of the first dilepton pair in the offshell region (Max = 146~GeV).

Mass of the first dilepton pair in the $ZZ$ mass region (Max = 232~GeV).

Mass of the second dilepton pair in the $Z$ mass region (Max = 44.2~GeV).

Mass of the second dilepton pair in the $H$ mass region (Max = 59.6~GeV).

Mass of the second dilepton pair in the offhell region (Max = 88.3~GeV).

Mass of the second dilepton pair in the $ZZ$ mass region (Max = 1000~GeV).

$p_T$ of the first dilepton pair in the $Z$ mass region (Max = 291~GeV).

$p_T$ of the first dilepton pair in the $H$ mass region (Max = 245~GeV).

$p_T$ of the first dilepton pair in the offshell region (Max = 216~GeV).

$p_T$ of the first dilepton pair in the $ZZ$ mass region (Max = 555~GeV).

$p_T$ of the second dilepton pair in the $Z$ mass region (Max = 117~GeV).

$p_T$ of the second dilepton pair in the $H$ mass region (Max = 124~GeV).

$p_T$ of the second dilepton pair in the offhell region (Max = 192~GeV).

$p_T$ of the second dilepton pair in the $ZZ$ mass region (Max = 555~GeV).

$\cos\theta^*$ between first dilepton pair in the $Z$ mass region.

$\cos\theta^*$ between first dilepton pair in the offshell region.

$\cos\theta^*$ between first dilepton pair in the $H$ mass region.

$\cos\theta^*$ between first dilepton pair in the $ZZ$ mass region.

$\cos\theta^*$ between second dilepton pair in the $Z$ mass region.

$\cos\theta^*$ between second dilepton pair in the $H$ mass region.

$\cos\theta^*$ between second dilepton pair in the offshell mass region.

$\cos\theta^*$ between second dilepton pair in the $ZZ$ mass region.

$\Delta y$ between the dilepton pairs in the $Z$ mass region.

$\Delta y$ between the dilepton pairs in the $H$ mass region.

$\Delta y$ between the dilepton pairs in the offshell region.

$\Delta y$ between the dilepton pairs in the $ZZ$ mass region.

$\Delta \phi$ between the dilepton pairs in the $Z$ mass region.

$\Delta \phi$ between the dilepton pairs in the $H$ mass region.

$\Delta \phi$ between the dilepton pairs in the offshell region.

$\Delta \phi$ between the dilepton pairs in the $ZZ$ mass region.

$\Delta \phi$ between the leading lepton pairs in the $Z$ mass region.

$\Delta \phi$ between the leading lepton pairs in the $H$ mass region.

$\Delta \phi$ between the leding lepton pairs in the offshell region.

$\Delta \phi$ between the leading lepton pairs in the $ZZ$ mass region.

Four lepton mass in first $p_T^{4l}$ slice (Max = 726~GeV)

Four lepton mass in second $p_T^{4l}$ slice (Max = 1079~GeV)

Four lepton mass in third $p_T^{4l}$ slice (Max = 1153~GeV)

Four lepton mass in fourth $p_T^{4l}$ slice (Max = 1320~GeV)

Four lepton mass in fifth $p_T^{4l}$ slice (Max = 1413~GeV)

Four lepton mass in first $y^{4l}$ slice (Max = 1710~GeV)

Four lepton mass in second $y^{4l}$ slice (Max = 1116~GeV)

Four lepton mass in third $y^{4l}$ slice (Max = 1413~GeV)

Four lepton mass in fourth $y^{4l}$ slice (Max = 1044~GeV)

Four lepton mass in fifth $y^{4l}$ slice (Max = 918~GeV)

Covariance matrix for M4l

Covariance matrix for M4lvChannel

Covariance matrix for m12_vs_M4l

Covariance matrix for m34_vs_M4l

Covariance matrix for pt12_vs_M4l

Covariance matrix for pt34_vs_M4l

Covariance matrix for CTS12_vs_M4l

Covariance matrix for CTS34_vs_M4l

Covariance matrix for dYpairs_vs_M4l

Covariance matrix for dPhiPairs_vs_M4l

Covariance matrix for dPhiLeadLep_vs_M4l

Covariance matrix for M4lvPt4lbin

Covariance matrix for M4lvRapiditybin

EFT results (linear)

EFT results (full)

The +1 $\sigma$ band of the observed exclusion contor for the $Z^{\prime}_{B}$ model in the $m_{\chi}$-$m_{Z^{\prime}_{B}}$ plane.

The -1 $\sigma$ band of the observed exclusion contor for the $Z^{\prime}_{B}$ model in the $m_{\chi}$-$m_{Z^{\prime}_{B}}$ plane.

A comparison of the inferred limits to the constraints from direct detection experiments on the spin-independent DM--nucleon cross section in the context of the $Z'_B$ simplified model with vector couplings. Limits are shown at 90% CL.

The observed exclusion contor for the $Z^{\prime}$-2HDM model in the $m_{A}$-$m_{Z^{\prime}}$ plane.

The expected exclusion contor for the $Z^{\prime}$-2HDM model in the $m_{A}$-$m_{Z^{\prime}}$ plane.

The +1 $\sigma$ band of the observed exclusion contor for the $m_{A}$-$Z^{\prime}$-2HDM model in the $m_{Z^{\prime}}$ plane.

The -1 $\sigma$ band of the observed exclusion contor for the $m_{A}$-$Z^{\prime}$-2HDM model in the $m_{Z^{\prime}}$ plane.

The observed exclusion contor for the 2HDM-a model in the $m_{A}$-$m_{a}$ plane.

The expected exclusion contor for the 2HDM-a model in the $m_{A}$-$m_{a}$ plane.

The +1 $\sigma$ band of the observed exclusion contor for the 2HDM-a model in the $m_{A}$-$m_{a}$ plane.

The -1 $\sigma$ band of the observed exclusion contor for the 2HDM-a model in the $m_{A}$-$m_{a}$ plane.

The observed exclusion contor for the 2HDM-a model in the $tan\beta$-$m_{a}$ plane.

The expected exclusion contor for the 2HDM-a model in the $tan\beta$-$m_{a}$ plane.

The +1 $\sigma$ band of the observed exclusion contor for the 2HDM-a model in the $tan\beta$-$m_{a}$ plane.

The -1 $\sigma$ band of the observed exclusion contor for the 2HDM-a model in the $tan\beta$-$m_{a}$ plane.

The exclusion limits at 95% CL for the 2HDM+a model as a function of $\sin \theta$ for $m_{A,H^{\pm},H}$= 600GeV, $m_a$ = 200GeV, $\tan \beta$ = 1.0$.

The exclusion limits at 95% CL for the 2HDM+a model as a function of $\sin \theta$ for $m_{A,H^{\pm},H}$= 1000GeV, $m_a$ = 350GeV, $\tan \beta$ = 1.0$.

Breakdown of the dominant systematic uncertainties.

Event yields in the range of 120 $<m_{\gamma\gamma}<$ 130 GeV for data, signal models, the SM Higgs boson background and non-resonant background in each analysis category, for an integrated luminosity of $139$fb$^{-1}$.

Detailed background contributions from the SM Higgs boson and continuum background for each cut

Detailed contributions from the signals for each cut.

Acceptance times efficiency for several signals in each category.

The \(m_{bb}\) distributions in the MediumBDT categories, overlaid with contributions from the \(H\gamma jj\) signal as well as the resonant \(Z\gamma jj\) and non-resonant \(b\bar{b} \gamma jj\) background fits. The combined \(\chi^2\) per degree of freedom is \(45.2/45\). The bottom panel in each plot presents the significance of the Higgs boson signal relative to the non-resonant \(b\bar{b} \gamma jj\) background in each bin.

The \(m_{bb}\) distributions in the MediumBDT categories, overlaid with contributions from the \(H\gamma jj\) signal as well as the resonant \(Z\gamma jj\) and non-resonant \(b\bar{b} \gamma jj\) background fits. The combined \(\chi^2\) per degree of freedom is \(45.2/45\). The bottom panel in each plot presents the significance of the Higgs boson signal relative to the non-resonant \(b\bar{b} \gamma jj\) background in each bin.

Comparisons of data and simulated event distributions of the BDT input variables \(\Delta R (b1, \gamma)\) in the two \(m_{bb}\) sidebands after kinematic reweighting of the non-resonant \(b\bar{b}\gamma jj\) background. The data are shown as black points, and the background contributions are stacked in coloured histograms. The Higgs boson signal contribution is scaled up and represented by the dashed red line. The bottom panel in each plot shows the ratio between the data and the SM prediction, where the uncertainty band corresponds to the statistical uncertainty only.

Comparisons of data and simulated event distributions of the BDT input variables \(m_{jj}\) in the two \(m_{bb}\) sidebands after kinematic reweighting of the non-resonant \(b\bar{b}\gamma jj\) background. The data are shown as black points, and the background contributions are stacked in coloured histograms. The Higgs boson signal contribution is scaled up and represented by the dashed red line. The bottom panel in each plot shows the ratio between the data and the SM prediction, where the uncertainty band corresponds to the statistical uncertainty only.

Comparisons of data and simulated event distributions of the BDT input variables \(\Delta R (b2, \gamma)\) in the two \(m_{bb}\) sidebands after kinematic reweighting of the non-resonant \(b\bar{b}\gamma jj\) background. The data are shown as black points, and the background contributions are stacked in coloured histograms. The Higgs boson signal contribution is scaled up and represented by the dashed red line. The bottom panel in each plot shows the ratio between the data and the SM prediction, where the uncertainty band corresponds to the statistical uncertainty only.

Comparisons of data and simulated event distributions of the BDT input variables \(\text{centrality}(\gamma, j1, j2)\) in the two \(m_{bb}\) sidebands after kinematic reweighting of the non-resonant \(b\bar{b}\gamma jj\) background. The data are shown as black points, and the background contributions are stacked in coloured histograms. The Higgs boson signal contribution is scaled up and represented by the dashed red line. The bottom panel in each plot shows the ratio between the data and the SM prediction, where the uncertainty band corresponds to the statistical uncertainty only.

Comparisons of data and simulated event distributions of the BDT input variables \(\Delta \phi(b\bar{b}, jj)\) in the two \(m_{bb}\) sidebands after kinematic reweighting of the non-resonant \(b\bar{b}\gamma jj\) background. The data are shown as black points, and the background contributions are stacked in coloured histograms. The Higgs boson signal contribution is scaled up and represented by the dashed red line. The bottom panel in each plot shows the ratio between the data and the SM prediction, where the uncertainty band corresponds to the statistical uncertainty only.

Comparisons of data and simulated event distributions of the BDT input variables \(p_{\mathrm{T}}^{jj}\) in the two \(m_{bb}\) sidebands after kinematic reweighting of the non-resonant \(b\bar{b}\gamma jj\) background. The data are shown as black points, and the background contributions are stacked in coloured histograms. The Higgs boson signal contribution is scaled up and represented by the dashed red line. The bottom panel in each plot shows the ratio between the data and the SM prediction, where the uncertainty band corresponds to the statistical uncertainty only.

Comparisons of data and simulated event distributions of the BDT input variables \(\cos \theta\) in the two \(m_{bb}\) sidebands after kinematic reweighting of the non-resonant \(b\bar{b}\gamma jj\) background. The data are shown as black points, and the background contributions are stacked in coloured histograms. The Higgs boson signal contribution is scaled up and represented by the dashed red line. The bottom panel in each plot shows the ratio between the data and the SM prediction, where the uncertainty band corresponds to the statistical uncertainty only.

Comparisons of data and simulated event distributions of the BDT input variables \(\Delta R(b1,j1)\) in the two \(m_{bb}\) sidebands after kinematic reweighting of the non-resonant \(b\bar{b}\gamma jj\) background. The data are shown as black points, and the background contributions are stacked in coloured histograms. The Higgs boson signal contribution is scaled up and represented by the dashed red line. The bottom panel in each plot shows the ratio between the data and the SM prediction, where the uncertainty band corresponds to the statistical uncertainty only.

Observed and expected limits on $\tan\beta$ as a function of $m_{H^+}$ in the $M_h^{125}(\tilde{\chi})$ scenario. Limits are shown for $\tan\beta$ values in the range of 0.5-60 due to the availability of the model prediction. The bands surrounding the expected limits show the 68% and 95% confidence intervals. Uncertainties in the predicted $H^+$ cross-sections or branching ratios are not considered.

Observed and expected limits on $\tan\beta$ as a function of $m_{H^+}$ in the $M_h^{125}(\tilde{\tau})$ scenario. Limits are shown for $\tan\beta$ values in the range of 0.5-60 due to the availability of the model prediction. The bands surrounding the expected limits show the 68% and 95% confidence intervals. Uncertainties in the predicted $H^+$ cross-sections or branching ratios are not considered.

Observed and expected limits on $\tan\beta$ as a function of $m_{H^+}$ in the $M_h^{125}$(alignment) scenario. Limits are shown for $\tan\beta$ values in the range of 0.5-60 due to the availability of the model prediction. The bands surrounding the expected limits show the 68% and 95% confidence intervals. Uncertainties in the predicted $H^+$ cross-sections or branching ratios are not considered.

Observed and expected limits on $\tan\beta$ as a function of $m_{H^+}$ in the $M_{h_1}^{125}$(CPV) scenario. Limits are shown for $\tan\beta$ values in the range of 0.5-60 due to the availability of the model prediction. The bands surrounding the expected limits show the 68% and 95% confidence intervals. Uncertainties in the predicted $H^+$ cross-sections or branching ratios are not considered.

Event yields of the SM background processes and the 800 GeV $H^{+}$ sample in the four analysis regions before the fit to the data. Uncertainties include both statistical and systematic uncertainties. The yields of the $H^{+}$ signal sample correspond to a cross-section times branching fraction of 10 pb.

Event acceptance for the different $H^+$ mass signal samples.

Efficiency in the Single-bin SR as a function of the M(Sbottom) vs. M(N2) with the $\Delta$ M(N2,N1) $= 130$ GeV. Keep in mind that the efficiency is given in units of $10^{-2}$.

Acceptance in the Multi-bin SR, $\min_{\Theta} < 0.5$ bin as a function of the M(Sbottom) vs. M(N2) with the $\Delta M$(N2,N1) = 130 GeV. Keep in mind that the acceptance is given in units of $10^{-4}$.

Efficiency in the Multi-bin SR, $\min_{\Theta} < 0.5$ bin as a function of the M(Sbottom) vs. M(N2) with the $\Delta M$(N2,N1) = 130 GeV. Keep in mind that the efficiency is given in units of $10^{-2}$.

Acceptance in the Multi-bin SR, $0.5 < \min_{\Theta} < 1.0$ bin as a function of the M(Sbottom) vs. M(N2) with the $\Delta M$(N2,N1) = 130 GeV. Keep in mind that the acceptance is given in units of $10^{-4}$.

Efficiency in the Multi-bin SR, $0.5 < \min_{\Theta} < 1.0$ bin as a function of the M(Sbottom) vs. M(N2) with the $\Delta M$(N2,N1) = 130 GeV. Keep in mind that the efficiency is given in units of $10^{-2}$.

Acceptance in the Multi-bin SR, $\min_{\Theta} > 1.0$ bin as a function of the M(Sbottom) vs. M(N2) with the $\Delta M$(N2,N1) = 130 GeV. Keep in mind that the acceptance is given in units of $10^{-4}$.

Efficiency in the Multi-bin SR, $\min_{\Theta} > 1.0$ bin as a function of the M(Sbottom) vs. M(N2) with the $\Delta M$(N2,N1) = 130 GeV. Keep in mind that the efficiency is given in units of $10^{-2}$.

Observed upper limits on the signal cross section as a function of the M(Sbottom) vs. M(N2) with the $\Delta M$(N2,N1) = 130 GeV.

Expected upper limits on the signal cross section as a function of the M(Sbottom) vs. M(N2) with the $\Delta M$(N2,N1) = 130 GeV.

Cutflows for the bechmarl signal point M(Sbottom) = 800 GeV, M(N2) = 180 GeV. Weighted event yields are reported starting with the "Preselection" line, normalized to an integrated luminosity of $139$ fb$^{−1}$.

Comparison of the expected and observed event yields in the signal regions. The top-quark and Z(mumu) background contributions are scaled with the normalization factors obtained from the background-only fit. The other contribution includes all the backgrounds not explicitly listed in the legend (V+jets except Z(mumu)+jets, di-/triboson, multijet). The hatched band indicates the total statistical and systematic uncertainties in the SM background. The contributions from three signal models to the signal regions are also displayed, where the masses M(Sbottom) and M(N2) are given in GeV in the legend. The lower panel shows the significance of the deviation of the observed yield from the expected background yield.

Dominant systematic uncertainties in the background prediction for the signal regions after the fit to the control regions. “Other” includes the uncertainties arising from muons, jet-vertex tagging, modeling of pile-up, the $E_{T}^{miss}$ computation, multijet background, and luminosity. The individual uncertainties can be correlated and do not necessarily add up quadratically to the total uncertainty.

Measured and predicted fiducial cross-sections in both the lllljj and ll$\nu\nu$jj channels for the inclusive ZZjj processes. Uncertainties due to different sources are presented.

Measured and predicted fiducial cross-sections in both the lllljj and ll$\nu\nu$jj channels for the inclusive ZZjj processes. Uncertainties due to different sources are presented.

Observed and expected multivariate discriminant distribution in the $\ell\ell\ell\ell jj$ QCD CR.

Observed and expected multivariate discriminant distribution in the $\ell\ell\ell\ell jj$ QCD CR.

Observed and expected multivariate discriminant distribution in the $\ell\ell\ell\ell jj$ SR.

Observed and expected multivariate discriminant distribution in the $\ell\ell\ell\ell jj$ SR.

Observed and expected multivariate discriminant distribution in the $\ell\ell\nu\nu jj$ SR.

Observed and expected multivariate discriminant distribution in the $\ell\ell\nu\nu jj$ SR.