Observed 95% CL upper limits on the $\tilde{t}$ production cross section, as functions of $m_{\tilde{t}}$ and $\Delta m$, for $\mathcal{B}(\tilde{t} \to bf\overline{f}'\tilde{\chi}^{0}_{1})$ of 50%. The observed (solid black) and expected (dashed red) exclusion curves are overlaid on the plots. The search excludes the region to the left of the exclusion curves.

Observed 95% CL upper limits on the $\tilde{t}$ production cross section, as functions of $m_{\tilde{t}}$ and $\Delta m$, for $\mathcal{B}(\tilde{t} \to bf\overline{f}'\tilde{\chi}^{0}_{1})$ of 50%. The observed (solid black) and expected (dashed red) exclusion curves are overlaid on the plots. The search excludes the region to the left of the exclusion curves.

Observed 95% CL upper limits on the $\tilde{t}$ production cross section, as functions of $m_{\tilde{t}}$ and $\Delta m$, for $\mathcal{B}(\tilde{t} \to bf\overline{f}'\tilde{\chi}^{0}_{1})$ of 100%. The observed (solid black) and expected (dashed red) exclusion curves are overlaid on the plots. The search excludes the region to the left of the exclusion curves.

Observed 95% CL upper limits on the $\tilde{t}$ production cross section, as functions of $m_{\tilde{t}}$ and $\Delta m$, for $\mathcal{B}(\tilde{t} \to bf\overline{f}'\tilde{\chi}^{0}_{1})$ of 100%. The observed (solid black) and expected (dashed red) exclusion curves are overlaid on the plots. The search excludes the region to the left of the exclusion curves.

Observed 95% CL upper limits on the production cross section for the bino-wino NLSP model, as functions of $m_{ ext{LLP}}$ and $c\tau$, for $\Delta m$ of 12 GeV. The observed (solid black) and expected (dashed red) exclusion curves are overlaid on the plots. The search excludes the region to the left of the exclusion curves.

Observed 95% CL upper limits on the production cross section for the bino-wino NLSP model, as functions of $m_{ ext{LLP}}$ and $c\tau$, for $\Delta m$ of 12 GeV. The observed (solid black) and expected (dashed red) exclusion curves are overlaid on the plots. The search excludes the region to the left of the exclusion curves.

Observed 95% CL upper limits on the production cross section for the bino-wino NLSP model, as functions of the $m_{ ext{LLP}}$ and $c\tau$, for $\Delta m$ of 15 GeV. The observed (solid black) and expected (dashed red) exclusion curves are overlaid on the plots. The search excludes the region to the left of the exclusion curves.

Observed 95% CL upper limits on the production cross section for the bino-wino NLSP model, as functions of the $m_{ ext{LLP}}$ and $c\tau$, for $\Delta m$ of 15 GeV. The observed (solid black) and expected (dashed red) exclusion curves are overlaid on the plots. The search excludes the region to the left of the exclusion curves.

Observed 95% CL upper limits on the production cross section for the bino-wino NLSP model, as functions of the $m_{ ext{LLP}}$ and $c\tau$, for $\Delta m$ of 20 GeV. The observed (solid black) and expected (dashed red) exclusion curves are overlaid on the plots. The search excludes the region to the left of the exclusion curves.

Observed 95% CL upper limits on the production cross section for the bino-wino NLSP model, as functions of the $m_{ ext{LLP}}$ and $c\tau$, for $\Delta m$ of 20 GeV. The observed (solid black) and expected (dashed red) exclusion curves are overlaid on the plots. The search excludes the region to the left of the exclusion curves.

Observed 95% CL upper limits on the production cross section for the bino-wino NLSP model, as functions of the $m_{ ext{LLP}}$ and $c\tau$, for $\Delta m$ of 25 GeV. The observed (solid black) and expected (dashed red) exclusion curves are overlaid on the plots. The search excludes the region to the left of the exclusion curves.

Observed 95% CL upper limits on the production cross section for the bino-wino NLSP model, as functions of the $m_{ ext{LLP}}$ and $c\tau$, for $\Delta m$ of 25 GeV. The observed (solid black) and expected (dashed red) exclusion curves are overlaid on the plots. The search excludes the region to the left of the exclusion curves.

Example description

Example description

Normalized differential distributions of the $\Delta R (\gamma,\mathrm{tt})$. The nominal MC prediction used to compare the experimental results to is obtained with Madgraph5 at NLO in QCD for photons from the production part of the process and Madgraph5 at LO in QCD for photons from the decay part of the process. The alternative prediction is obtained with Madgraph5 at NLO in QCD for photons from the production part of the process and a POWHEG+Pythia $\mathrm{tt}$ simulation at NLO in QCD for photons from the decay part of the process.

Absolute differential distributions of the min. $\Delta R (\gamma,\mathrm{t})$. The nominal MC prediction used to compare the experimental results to is obtained with Madgraph5 at NLO in QCD for photons from the production part of the process and Madgraph5 at LO in QCD for photons from the decay part of the process. The alternative prediction is obtained with Madgraph5 at NLO in QCD for photons from the production part of the process and a POWHEG+Pythia $\mathrm{tt}$ simulation at NLO in QCD for photons from the decay part of the process.

Normalized differential distributions of the min. $\Delta R (\gamma,\mathrm{t})$. The nominal MC prediction used to compare the experimental results to is obtained with Madgraph5 at NLO in QCD for photons from the production part of the process and Madgraph5 at LO in QCD for photons from the decay part of the process. The alternative prediction is obtained with Madgraph5 at NLO in QCD for photons from the production part of the process and a POWHEG+Pythia $\mathrm{tt}$ simulation at NLO in QCD for photons from the decay part of the process.

Absolute differential distributions of the $m(\mathrm{tt})$. The nominal MC prediction used to compare the experimental results to is obtained with Madgraph5 at NLO in QCD for photons from the production part of the process and Madgraph5 at LO in QCD for photons from the decay part of the process. The alternative prediction is obtained with Madgraph5 at NLO in QCD for photons from the production part of the process and a POWHEG+Pythia $\mathrm{tt}$ simulation at NLO in QCD for photons from the decay part of the process.

Normalized differential distributions of the $m(\mathrm{tt})$. The nominal MC prediction used to compare the experimental results to is obtained with Madgraph5 at NLO in QCD for photons from the production part of the process and Madgraph5 at LO in QCD for photons from the decay part of the process. The alternative prediction is obtained with Madgraph5 at NLO in QCD for photons from the production part of the process and a POWHEG+Pythia $\mathrm{tt}$ simulation at NLO in QCD for photons from the decay part of the process.

Absolute differential distributions of the leading lepton $p_{\mathrm{T}}$. The nominal MC prediction used to compare the experimental results to is obtained with Madgraph5 at NLO in QCD for photons from the production part of the process and Madgraph5 at LO in QCD for photons from the decay part of the process. The alternative prediction is obtained with Madgraph5 at NLO in QCD for photons from the production part of the process and a POWHEG+Pythia $\mathrm{tt}$ simulation at NLO in QCD for photons from the decay part of the process. The theory predictions are computed at fixed order, at NLO in QCD, as described in Section 9.2 of the paper.

Normalized differential distributions of the leading lepton $p_{\mathrm{T}}$. The nominal MC prediction used to compare the experimental results to is obtained with Madgraph5 at NLO in QCD for photons from the production part of the process and Madgraph5 at LO in QCD for photons from the decay part of the process. The alternative prediction is obtained with Madgraph5 at NLO in QCD for photons from the production part of the process and a POWHEG+Pythia $\mathrm{tt}$ simulation at NLO in QCD for photons from the decay part of the process. The theory predictions are computed at fixed order, at NLO in QCD, as described in Section 9.2 of the paper.

Absolute differential distributions of the photon $p_{\mathrm{T}}$. The nominal MC prediction used to compare the experimental results to is obtained with Madgraph5 at NLO in QCD for photons from the production part of the process and Madgraph5 at LO in QCD for photons from the decay part of the process. The alternative prediction is obtained with Madgraph5 at NLO in QCD for photons from the production part of the process and a POWHEG+Pythia $\mathrm{tt}$ simulation at NLO in QCD for photons from the decay part of the process. The theory predictions are computed at fixed order, at NLO in QCD, as described in Section 9.2 of the paper.

Normalized differential distributions of the photon $p_{\mathrm{T}}$. The nominal MC prediction used to compare the experimental results to is obtained with Madgraph5 at NLO in QCD for photons from the production part of the process and Madgraph5 at LO in QCD for photons from the decay part of the process. The alternative prediction is obtained with Madgraph5 at NLO in QCD for photons from the production part of the process and a POWHEG+Pythia $\mathrm{tt}$ simulation at NLO in QCD for photons from the decay part of the process. The theory predictions are computed at fixed order, at NLO in QCD, as described in Section 9.2 of the paper.

Absolute differential distributions of the $\Delta \phi (\ell,\ell)$. The nominal MC prediction used to compare the experimental results to is obtained with Madgraph5 at NLO in QCD for photons from the production part of the process and Madgraph5 at LO in QCD for photons from the decay part of the process. The alternative prediction is obtained with Madgraph5 at NLO in QCD for photons from the production part of the process and a POWHEG+Pythia $\mathrm{tt}$ simulation at NLO in QCD for photons from the decay part of the process. The theory predictions are computed at fixed order, at NLO in QCD, as described in Section 9.2 of the paper.

Normalized differential distributions of the $\Delta \phi (\ell,\ell)$. The nominal MC prediction used to compare the experimental results to is obtained with Madgraph5 at NLO in QCD for photons from the production part of the process and Madgraph5 at LO in QCD for photons from the decay part of the process. The alternative prediction is obtained with Madgraph5 at NLO in QCD for photons from the production part of the process and a POWHEG+Pythia $\mathrm{tt}$ simulation at NLO in QCD for photons from the decay part of the process. The theory predictions are computed at fixed order, at NLO in QCD, as described in Section 9.2 of the paper.

Absolute differential distributions of the ratio between the cross sections of $\mathrm{tt\gamma}$ and $\mathrm{tt}$, as a function of the leading top quark $p_{\mathrm{T}}$. The nominal MC prediction used to compare the experimental results to is obtained with Madgraph5 at NLO in QCD for photons from the production part of the process and Madgraph5 at LO in QCD for photons from the decay part of the process. The alternative prediction is obtained with Madgraph5 at NLO in QCD for photons from the production part of the process and a POWHEG+Pythia $\mathrm{tt}$ simulation at NLO in QCD for photons from the decay part of the process.

Absolute differential distributions of the ratio between the cross sections of $\mathrm{tt\gamma}$ and $\mathrm{tt}$, as a function of the leading lepton $p_{\mathrm{T}}$. The nominal MC prediction used to compare the experimental results to is obtained with Madgraph5 at NLO in QCD for photons from the production part of the process and Madgraph5 at LO in QCD for photons from the decay part of the process. The alternative prediction is obtained with Madgraph5 at NLO in QCD for photons from the production part of the process and a POWHEG+Pythia $\mathrm{tt}$ simulation at NLO in QCD for photons from the decay part of the process. The theory predictions are computed at fixed order, at NLO in QCD, as described in Section 9.2 of the paper.

Correlation matrix of the total uncertainty in the absolute differential cross section as a function of the leading top quark $p_{\mathrm{T}}$.

Correlation matrix of the statistical uncertainty in the absolute differential cross section as a function of the leading top quark $p_{\mathrm{T}}$.

Correlation matrix of the total uncertainty in the absolute differential cross section as a function of the $\Delta R (\gamma,\mathrm{tt})$.

Correlation matrix of the statistical uncertainty in the absolute differential cross section as a function of the $\Delta R (\gamma,\mathrm{tt})$.

Correlation matrix of the total uncertainty in the absolute differential cross section as a function of the min. $\Delta R (\gamma,\mathrm{t})$.

Correlation matrix of the statistical uncertainty in the absolute differential cross section as a function of the min. $\Delta R (\gamma,\mathrm{t})$.

Correlation matrix of the total uncertainty in the absolute differential cross section as a function of the $m(\mathrm{tt})$.

Correlation matrix of the statistical uncertainty in the absolute differential cross section as a function of the $m(\mathrm{tt})$.

Correlation matrix of the total uncertainty in the absolute differential cross section as a function of the leading lepton $p_{\mathrm{T}}$.

Correlation matrix of the statistical uncertainty in the absolute differential cross section as a function of the leading lepton $p_{\mathrm{T}}$.

Correlation matrix of the total uncertainty in the absolute differential cross section as a function of the photon $p_{\mathrm{T}}$.

Correlation matrix of the statistical uncertainty in the absolute differential cross section as a function of the photon $p_{\mathrm{T}}$.

Correlation matrix of the total uncertainty in the absolute differential cross section as a function of the $\Delta \phi (\ell,\ell)$.

Correlation matrix of the statistical uncertainty in the absolute differential cross section as a function of the $\Delta \phi (\ell,\ell)$.

Correlation matrix of the total uncertainty in the absolute differential $\mathrm{tt\gamma}$ / $\mathrm{tt}$ cross section ratio as a function of the leading top $p_{\mathrm{T}}$.

Correlation matrix of the statistical uncertainty in the absolute differential $\mathrm{tt\gamma}$ / $\mathrm{tt}$ cross section ratio as a function of the leading top $p_{\mathrm{T}}$.

Correlation matrix of the total uncertainty in the absolute differential $\mathrm{tt\gamma}$ / $\mathrm{tt}$ cross section ratio as a function of the leading lepton $p_{\mathrm{T}}$.

Correlation matrix of the statistical uncertainty in the absolute differential $\mathrm{tt\gamma}$ / $\mathrm{tt}$ cross section ratio as a function of the leading lepton $p_{\mathrm{T}}$.

Distributions in the trailing lepton pt for the DY CR after a fit to the data.

Distributions in the trailing lepton pt for the tt CR after a fit to the data.

Likelihood scan of the signal strengths for tWZ and ttZ production. The black cross shows the best fit value, while the black diamond indicates the SM expected value. The con- tours in different shades of blue correspond to the 68, 95, and 99% confidence level.

The 95% CL upper limits on the inclusive signal strength $ r = \sigma_{HH} /\sigma_{HH} ^\text{SM} $ for each channel and their combination. The inner (green) band and the outer (yellow) band indicate the 68 and 95% CL intervals, respectively, under the background-only hypothesis. The $ b\overline{b}b\overline{b} $ and $ b\overline{b}WW $ contributions have been combined in order to simplify the presentation of results.

The 95$\%$ CL upper limits on the VBF signal strength $r_{VBF HH} = \sigma_{VBF HH}/\sigma_{VBF HH}^{SM}$ for each channel and their combination. The inner (green) band and the outer (yellow) band indicate the 68 and 95$\%$ CL intervals, respectively, under the background-only hypothesis. Not all of the eight channels are used to extract the combined VBF limit. The contributing channels are indicated in the figure. The b\overline{b}b\overline{b} $ and $ b\overline{b}WW $ contributions have been combined in order to simplify the presentation of results.

The 95$\%$ CL upper limits on the inclusive HH cross section as a function of $\kappa_\lambda$. All other couplings are set to the values predicted by the SM. The theoretical uncertainties in the HH GGF and VBF signal cross sections are not considered because here we directly constrain the measured cross section.The inner (green) band and the outer (yellow) band indicate the 68 and 95$\%$ CL intervals, respectively, under the background-only hypothesis. The star shows the limit at the SM value for $\kappa_\lambda$.

The 95$\%$ CL upper limits on the inclusive HH cross section as a function of $\kappa_{2V}$. All other couplings are set to the values predicted by the SM. The theoretical uncertainties in the HH VBF signal cross sections are not considered because here we directly constrain the measured cross section.The inner (green) band and the outer (yellow) band indicate the 68 and 95$\%$ CL intervals, respectively, under the background-only hypothesis. The star shows the limit at the SM value for $\kappa_{2V}$.

The profile likelihood ratio test statistic $-2\Delta\log(L)$ as a function of coupling modifiers $\kappa_\lambda$ for the combination of all channels, when all the other parameters are fixed to their SM values.

The profile likelihood ratio test statistic $-2\Delta\log(L)$ as a function of coupling modifiers $\kappa_{2V}$ for the combination of all channels, when all the other parameters are fixed to their SM values.

The 68 $\%$ CL expected contours of $-2\Delta\log(L)$ in the ($\kappa_\lambda$, $\kappa_{2V}$) plane for the combination of all channels when all the other parameters are fixed to their SM values.

The 95$\%$ CL expected contours of $-2\Delta\log(L)$ in the ($\kappa_\lambda$, $\kappa_{2V}$) plane for the combination of all channels when all the other parameters are fixed to their SM values.

The 68 $\%$ CL observed contours of $-2\Delta\log(L)$ in the ($\kappa_\lambda$, $\kappa_{2V}$) plane for the combination of all channels when all the other parameters are fixed to their SM values.

The 95$\%$ observed CL contours of $-2\Delta\log(L)$ in the ($\kappa_\lambda$, $\kappa_{2V}$) plane for the combination of all channels when all the other parameters are fixed to their SM values.

The 5 standard deviation expected contours of $-2\Delta\log(L)$ in the ($\kappa_\lambda$, $\kappa_{2V}$) plane for the combination of all channels when all the other parameters are fixed to their SM values.

The 5 standard deviation observed contours of $-2\Delta\log(L)$ in the ($\kappa_\lambda$, $\kappa_{2V}$) plane for the combination of all channels when all the other parameters are fixed to their SM values.

The observed best-fit values in the ($\kappa_\lambda$, $\kappa_{2V}$) plane for the combination of all channels when all the other parameters are fixed to their SM values.

The 68 $\%$ expected CL contours of $-2\Delta\log(L)$ in the ($\kappa_{V}$, $\kappa_{2V}$) plane for the combination of all channels when all the other parameters are fixed to their SM values.

The 68$\%$ expected CL contours of $-2\Delta\log(L)$ in the ($\kappa_{V}$, $\kappa_{2V}$) plane for the combination of all channels when all the other parameters are fixed to their SM values.

The 95$\%$ CL expected contours of $-2\Delta\log(L)$ in the ($\kappa_{V}$, $\kappa_{2V}$) plane for the combination of all channels when all the other parameters are fixed to their SM values.

The 95$\%$ CL observed contours of $-2\Delta\log(L)$ in the ($\kappa_{V}$, $\kappa_{2V}$) plane for the combination of all channels when all the other parameters are fixed to their SM values.

The 5 standard deviation expected contours of $-2\Delta\log(L)$ in the ($\kappa_{V}$, $\kappa_{2V}$) plane for the combination of all channels when all the other parameters are fixed to their SM values.

The 5 standard deviation observed contours of $-2\Delta\log(L)$ in the ($\kappa_{V}$, $\kappa_{2V}$) plane for the combination of all channels when all the other parameters are fixed to their SM values.

The observed best-fit values in the ($\kappa_{V}$, $\kappa_{2V}$) plane for the combination of all channels when all the other parameters are fixed to their SM values.

The 68 $\%$ CL expected contours of $-2\Delta\log(L)$ in the ($\kappa_\lambda$, $\kappa_t$) plane for the combination of all channels when all the other parameters are fixed to their SM values.

The 68 $\%$ CL observed contours of $-2\Delta\log(L)$ in the ($\kappa_\lambda$, $\kappa_t$) plane for the combination of all channels when all the other parameters are fixed to their SM values.

The 95$\%$ CL expected contours of $-2\Delta\log(L)$ in the ($\kappa_\lambda$, $\kappa_t$) plane for the combination of all channels when all the other parameters are fixed to their SM values.

The 95$\%$ CL observed contours of $-2\Delta\log(L)$ in the ($\kappa_\lambda$, $\kappa_t$) plane for the combination of all channels when all the other parameters are fixed to their SM values.

The observed best-fit values in the ($\kappa_\lambda$, $\kappa_t$) plane for the combination of all channels when all the other parameters are fixed to their SM values.

Upper limits on the HH production cross section at 95$\%$ CL for the two sets of HEFT benchmarks. The theoretical uncertainties in the HH GGF signal cross section are not considered because we directly constrain the measured cross section.

Upper limits on the HH cross section as a function of the $C_2$ coupling modifier. The theoretical uncertainties in the HH GGF signal cross section are not considered because we directly constrain the measured cross section.

The profile likelihood ratio test statistic $-2\Delta\log(L)$ as a function of the $c_2$ coupling modifier.

The 68$\%$CL expected contours of $-2\Delta\log(L)$ in the ($c_2$,$\kappa_t$) for the combination of all channels when all the other parameters are fixed to their SM value.

The 68$\%$CL observed contours of $-2\Delta\log(L)$ in the ($c_2$,$\kappa_t$) for the combination of all channels when all the other parameters are fixed to their SM value.

The 95$\%$CL expected contours of $-2\Delta\log(L)$ in the ($c_2$,$\kappa_t$) for the combination of all channels when all the other parameters are fixed to their SM value.

The 95$\%$CL observed contours of $-2\Delta\log(L)$ in the ($c_2$,$\kappa_t$) for the combination of all channels when all the other parameters are fixed to their SM value.

The 68 $\%$ CL expected contours of $-2\Delta\log(L)$ in the ($c_2$, $\kappa_\lambda$) for the combination of all channels when all the other parameters are fixed to their SM value.

The 68 $\%$ CL observed contours of $-2\Delta\log(L)$ in the ($c_2$, $\kappa_\lambda$) for the combination of all channels when all the other parameters are fixed to their SM value.

The 95 $\%$ CL expected contours of $-2\Delta\log(L)$ in the ($c_2$, $\kappa_\lambda$) for the combination of all channels when all the other parameters are fixed to their SM value.

The 95 $\%$ CL observed contours of $-2\Delta\log(L)$ in the ($c_2$, $\kappa_\lambda$) for the combination of all channels when all the other parameters are fixed to their SM value.

The $-2\Delta\log(L)$ scan for the combination of all channels as a function of $|\cos{\alpha}|$ for singlet extensions of the SM with spontaneous $Z_2$ symmetry breaking. The $|\cos{\alpha}|$ is constrained at 95\% CL between 0.79 and 1. The ranges of $|\cos{\alpha}|$ and $\lambda_\text{eff}$ are chosen to guarantee the validity of the models. The excluded regions are below the curves shown.

The 95$\%$ CL expected contour of $-2\Delta\log(L)$ in the ($|\cos{\alpha}|$, $\lambda_{eff}$) plane, where $\lambda_{eff}$=$\lambda_\alpha$-$\tan{(\alpha)\frac{m_{2}}{\nu}}$. The ranges of $|\cos{\alpha}|$ and $\lambda_{eff}$ are chosen to guarantee the validity of the models. The excluded regions are to the left of the curves shown.

The 95$\%$ CL observed contour of $-2\Delta\log(L)$ in the ($|\cos{\alpha}|$, $\lambda_{eff}$) plane, where $\lambda_{eff}$=$\lambda_\alpha$-$\tan{(\alpha)\frac{m_{2}}{\nu}}$. The ranges of $|\cos{\alpha}|$ and $\lambda_{eff}$ are chosen to guarantee the validity of the models. The excluded regions are to the left of the curves shown.

The 95% CL expected contour of $-2\Delta\log(L)$ for the combination of all channels in the ($m_{H}$, $\tan\beta$) plane for a fixed value of $|\cos(\beta-\alpha)|=0.2$ in the 2HDM-I model. In these and all other cases considered in this paper, $m_H$ = $m_A$. The ranges of $m_{H}$ and $\tan\beta$ are chosen to guarantee the validity of the models. The excluded regions are below the curves shown. The value $\tan\beta$ = 0.5 is excluded for $m_{H}$ > 800 GeV for all models considered.

The 95$\%$ CL observed contour of $-2\Delta\log(L)$ for the combination of all channels in the ($m_{H}$, $\tan\beta$) plane for a fixed value of $|\cos(\beta-\alpha)|=0.2$ in the 2HDM-I model. In these and all other cases considered in this paper, $m_H$ = $m_A$. The ranges of $m_{H}$ and $\tan\beta$ are chosen to guarantee the validity of the models. The excluded regions are below the curves shown. The value $\tan\beta$ = 0.5 is excluded for $m_{H}$ > 800 GeV for all models considered.

The 95$\%$ CL contour of $-2\Delta\log(L)$ for the combination of all channels in the ($m_{H}$, $\tan\beta$) plane for a fixed value of $|\cos(\beta-\alpha)|=0.2$ for the 2HDM I models. In these and all other cases considered in this paper, $m_H$ = $m_A$. The ranges of $m_{H}$ and $\tan\beta$ are chosen to guarantee the validity of the models. The excluded regions are below the curves shown. The value $\tan\beta$ = 0.5 is excluded for $m_{H}$ > 800 GeV for all models considered.

The 95$\%$ CL contour of $-2\Delta\log(L)$ for the combination of all channels in the ($m_{H}$, $\tan\beta$) plane for a fixed value of $|\cos(\beta-\alpha)|=0.2$ for the 2HDM II models. In these and all other cases considered in this paper, $m_H$ = $m_A$. The ranges of $m_{H}$ and $\tan\beta$ are chosen to guarantee the validity of the models. The excluded regions are below the curves shown. The value $\tan\beta$ = 0.5 is excluded for $m_{H}$ > 800 GeV for all models considered.

The 95$\%$ CL contour of $-2\Delta\log(L)$ for the combination of all channels in the ($m_{H}$, $\tan\beta$) plane for a fixed value of $|\cos(\beta-\alpha)|=0.2$ for the 2HDM X models. In these and all other cases considered in this paper, $m_H$ = $m_A$. The ranges of $m_{H}$ and $\tan\beta$ are chosen to guarantee the validity of the models. The excluded regions are below the curves shown. The value $\tan\beta$ = 0.5 is excluded for $m_{H}$ > 800 GeV for all models considered.

The 95$\%$ CL contour of $-2\Delta\log(L)$ for the combination of all channels in the ($m_{H}$, $\tan\beta$) plane for a fixed value of $|\cos(\beta-\alpha)|=0.2$ for the 2HDM Y models. In these and all other cases considered in this paper, $m_H$ = $m_A$. The ranges of $m_{H}$ and $\tan\beta$ are chosen to guarantee the validity of the models. The excluded regions are below the curves shown. The value $\tan\beta$ = 0.5 is excluded for $m_{H}$ > 800 GeV for all models considered.

The 95$\%$ CL expected contour of $-2\Delta\log(L)$ for the combination of all channels in the ($|\cos(\beta-\alpha)|$, $\tan\beta$) plane for a fixed value of $m_{H}=1000GeV$ in the 2HDM-I model. In these and all other cases considered in this paper, $m_H$ = $m_A$. The ranges of $|\cos(\beta-\alpha)|$ and $\tan\beta$ are chosen to guarantee the validity of the models. The excluded regions are below the curves shown. The value $\tan\beta$ = 0.5 is excluded for $|\cos(\beta-\alpha)|$ > 0.16 and $|\cos(\beta-\alpha)|$ < -0.13 for all models considered.

The 95$\%$ CL observed contour of $-2\Delta\log(L)$ for the combination of all channels in the ($|\cos(\beta-\alpha)|$, $\tan\beta$) plane for a fixed value of $m_{H}=1000GeV$ in the 2HDM-I model. In these and all other cases considered in this paper, $m_H$ = $m_A$. The ranges of $|\cos(\beta-\alpha)|$ and $\tan\beta$ are chosen to guarantee the validity of the models. The excluded regions are below the curves shown. The value $\tan\beta$ = 0.5 is excluded for $|\cos(\beta-\alpha)|$ > 0.16 and $|\cos(\beta-\alpha)|$ < -0.13 for all models considered.

The 95% CL contour of $-2\Delta\log(L)$ for the combination of all channels in the ($|\cos(\beta-\alpha)|$, $\tan\beta$) plane for a fixed value of $m_{H}=1000GeV$ for the 2HDM I models. In these and all other cases considered in this paper, $m_H$ = $m_A$. The ranges of $|\cos(\beta-\alpha)|$ and $\tan\beta$ are chosen to guarantee the validity of the models. The excluded regions are below the curves shown. The value $\tan\beta$ = 0.5 is excluded for $|\cos(\beta-\alpha)|$ > 0.16 and $|\cos(\beta-\alpha)|$ < -0.13 for all models considered.

The 95% CL contour of $-2\Delta\log(L)$ for the combination of all channels in the ($|\cos(\beta-\alpha)|$, $\tan\beta$) plane for a fixed value of $m_{H}=1000GeV$ for the 2HDM II models. In these and all other cases considered in this paper, $m_H$ = $m_A$. The ranges of $|\cos(\beta-\alpha)|$ and $\tan\beta$ are chosen to guarantee the validity of the models. The excluded regions are below the curves shown. The value $\tan\beta$ = 0.5 is excluded for $|\cos(\beta-\alpha)|$ > 0.16 and $|\cos(\beta-\alpha)|$ < -0.13 for all models considered.

The 95% CL contour of $-2\Delta\log(L)$ for the combination of all channels in the ($|\cos(\beta-\alpha)|$, $\tan\beta$) plane for a fixed value of $m_{H}=1000GeV$ for the 2HDM X models. In these and all other cases considered in this paper, $m_H$ = $m_A$. The ranges of $|\cos(\beta-\alpha)|$ and $\tan\beta$ are chosen to guarantee the validity of the models. The excluded regions are below the curves shown. The value $\tan\beta$ = 0.5 is excluded for $|\cos(\beta-\alpha)|$ > 0.16 and $|\cos(\beta-\alpha)|$ < -0.13 for all models considered.

The 95% CL contour of $-2\Delta\log(L)$ for the combination of all channels in the ($|\cos(\beta-\alpha)|$, $\tan\beta$) plane for a fixed value of $m_{H}=1000GeV$ for the 2HDM Y models. In these and all other cases considered in this paper, $m_H$ = $m_A$. The ranges of $|\cos(\beta-\alpha)|$ and $\tan\beta$ are chosen to guarantee the validity of the models. The excluded regions are below the curves shown. The value $\tan\beta$ = 0.5 is excluded for $|\cos(\beta-\alpha)|$ > 0.16 and $|\cos(\beta-\alpha)|$ < -0.13 for all models considered.

The $-2\Delta\log(L)$ scan for the combination of all channels as a function of $\xi$ for MCHM$_4$. The range of $\xi$ is chosen to guarantee the validity of the model. At 95% CL, $\xi$ is constrained to be between 0 and 0.45 in MCHM$_4$.

The $-2\Delta\log(L)$ scan for the combination of all channels as a function of $\xi$ for MCHM$_5$. The range of $\xi$ is chosen to guarantee the validity of the model. At 95% CL, $\xi$ is constrained to be between 0 and 0.26 in MCHM$_5$.

Treatment of most important common systematic uncertainties in the S2 scenario.

The expected upper limits at 95% CL on the HH signal strength from the combination of all the considered channels at different integrated luminosities for the S2 scenario.

The expected upper limits at 95% CL on the HH signal strength from the combination of all the considered channels under different assumptions on the systematic uncertainties for an integrated luminosity of 3000fb$^{-1}$ (right).

The expected $-2\Delta\log(L)$ scan as a function of coupling modifier $\kappa_\lambda$ for the combination of all contributing channels, projected to different integrated luminosities. All other parameters are fixed to their SM values.

The expected $-2\Delta\log(L)$ scan as a function of coupling modifier $\kappa_\lambda$ for the combination of all contributing channels, projected under different assumptions on the systematic uncertainties for an integrated luminosity of 3000fb$^{-1}$. All other parameters are fixed to their SM values.

Expected significance for the HH signal projected to 2000 or 3000 fb$^{-1}$ under different assumptions of systematic uncertainties.

The expected signal significance as a function of integrated luminosity for the nominal scenario of systematic uncertainties S2 and for the scenario with statistical uncertainties only.

The expected signal significance as a function of $\kappa_\lambda$ under different assumptions on the systematic uncertainties for an integrated luminosity of 3000fb$^{-1}$.

The 95$\%$ CL upper limits on the inclusive HH signal strenngth as a function of $\kappa_{\lambda}$ coupling modifier. All other couplings are set to the values predicted by the SM. The theoretical uncertainties in the HH VBF signal cross sections are not considered because here we directly constrain the measured cross section.The inner (green) band and the outer (yellow) band indicate the 68 and 95$\%$ CL intervals, respectively, under the background-only hypothesis.

The 95$\%$ CL upper limits on the inclusive HH signal strenngth as a function of $\kappa_{2V}$ coupling modifier. All other couplings are set to the values predicted by the SM. The theoretical uncertainties in the HH VBF signal cross sections are not considered because here we directly constrain the measured cross section.The inner (green) band and the outer (yellow) band indicate the 68 and 95$\%$ CL intervals, respectively, under the background-only hypothesis.

The $-2\Delta\log(L)$ scan as a function of $\kappa_\lambda$ for all channels, when all the other parameters are fixed to their SM value.

The $-2\Delta\log(L)$ scan as a function of $\kappa_{2V}$ for all channels, when all the other parameters are fixed to their SM value.

The $-2\Delta\log(L)$ scan as a function of $\kappa_\lambda$ for all channels, when all the other parameters are fixed to their SM value.

The 95$\%$ CL upper limits on the inclusive HH cross section as a function of $\kappa_\lambda$ coupling modifier for all channels. All other couplings are set to the values predicted by the SM. The theoretical uncertainties in the HH GGF and VBF signal cross sections are not considered because here we directly constrain the measured cross section.

The 95$\%$ CL upper limits on the inclusive VBF HH cross section as a function of $\kappa_{2V}$ coupling modifier for all channels. All other couplings are set to the values predicted by the SM. The theoretical uncertainties in the HH GGF and VBF signal cross sections are not considered because here we directly constrain the measured cross section.

The 95$\%$ CL upper limits on the inclusive HH signal strength as a function of $\kappa_\lambda$ coupling modifier for all channels. All other couplings are set to the values predicted by the SM. The theoretical uncertainties in the HH GGF and VBF signal cross sections are not considered because here we directly constrain the measured cross section.

The 95$\%$ CL upper limits on the inclusive VBF HH signal strength as a function of $\kappa_{2V}$ coupling modifier for all channels. All other couplings are set to the values predicted by the SM. The theoretical uncertainties in the HH GGF and VBF signal cross sections are not considered because here we directly constrain the measured cross section.

The best fit value for $\kappa_\lambda$ compared to the SM expectation for all channels and their combination, when all the other parameters are fixed to their SM values.

The best fit value for $\kappa_2V$ compared to the SM expectation for all channels and their combination, when all the other parameters are fixed to their SM values.

The 68$\%$ CL expected contours of $-2\Delta\log(L)$ in the ($\kappa_{V}$, $\kappa_{2V}$) plane for the combination of all channels when all the other parameters are fixed to their SM values.

The 68$\%$ CL observed contours of $-2\Delta\log(L)$ in the ($\kappa_{V}$, $\kappa_{2V}$) plane for the combination of all channels when all the other parameters are fixed to their SM values.

The 95$\%$ CL expected contours of $-2\Delta\log(L)$ in the ($\kappa_{V}$, $\kappa_{2V}$) plane for the combination of all channels when all the other parameters are fixed to their SM values.

The 95$\%$ CL observed contours of $-2\Delta\log(L)$ in the ($\kappa_{V}$, $\kappa_{2V}$) plane for the combination of all channels when all the other parameters are fixed to their SM values.

The 5 $\sigma$ expected contours of $-2\Delta\log(L)$ in the ($\kappa_{V}$, $\kappa_{2V}$) plane for the combination of all channels when all the other parameters are fixed to their SM values.

The 5 $\sigma$ observed contours of $-2\Delta\log(L)$ in the ($\kappa_{V}$, $\kappa_{2V}$) plane for the combination of all channels when all the other parameters are fixed to their SM values.

The best-fit values of ($\kappa_{V}$, $\kappa_{2V}$) for the combination of all channels when all the other parameters are fixed to their SM values.

The 95$\%$ CL upper limits on the inclusive VBF HH signal strength as a function of $\kappa_{2V}$ coupling modifier for all channels. All other couplings are set to the values predicted by the SM. The theoretical uncertainties in the HH GGF and VBF signal cross sections are not considered because here we directly constrain the measured cross section. The inner (green) band and the outer (yellow) band indicate the 68 and 95% CL intervals, respectively, under the background-only hypothesis.

The $-2\Delta\log(L)$ scan as a function of the $c_2$ coupling modifier for all channels, when all the other parameters are fixed to their SM values.

The 95% CL upper limits on the inclusive HH cross section as a function of the $c_2$ coupling modifier for all channels. All other couplings are set to the values predicted by the SM. The theoretical uncertainties in the HH ggF signal cross sections are not considered because we directly constrain the measured cross section.

The 95% CL upper limits on the inclusive HH signal strength as a function of the $c_2$ coupling modifier for all channels. All other couplings are set to the values predicted by the SM. The theoretical uncertainties in the HH ggF signal cross sections are considered in this case

The best fit value for the $c_2$ coupling modifier compared to the SM expectation for all channels and their combination, when all the other parameters are fixed to their SM values.

The observed 68 $\%$ CL observed contours ($\kappa_\lambda$ =1.0) of $-2\Delta\log(L)$ in the ($c_2$, $\kappa_t$) plane for the combination of all channels when $\kappa_t$ is allowed to vary. The range of $\kappa_t$ is set between 0.5 and 1.5 to avoid unphysical areas of the phase space. All the other parameters are fixed to their SM values.

The observed 68 $\%$ CL observed contours (floating $\kappa_\lambda$) of $-2\Delta\log(L)$ in the ($c_2$, $\kappa_t$) plane for the combination of all channels when $\kappa_t$ is allowed to vary. The range of $\kappa_t$ is set between 0.5 and 1.5 to avoid unphysical areas of the phase space. All the other parameters are fixed to their SM values.

The observed 95 $\%$ CL observed contours ($\kappa_\lambda$ =1.0) of $-2\Delta\log(L)$ in the ($c_2$, $\kappa_t$) plane for the combination of all channels when $\kappa_t$ is allowed to vary. The range of $\kappa_t$ is set between 0.5 and 1.5 to avoid unphysical areas of the phase space. All the other parameters are fixed to their SM values.

The observed 95 $\%$ CL observed contours (floating $\kappa_\lambda$) of $-2\Delta\log(L)$ in the ($c_2$, $\kappa_t$) plane for the combination of all channels when $\kappa_t$ is allowed to vary. The range of $\kappa_t$ is set between 0.5 and 1.5 to avoid unphysical areas of the phase space. All the other parameters are fixed to their SM values.

The best fit values ($\kappa_\lambda$ =1.0) of ($c_2$, $\kappa_t$) for the combination of all channels when $\kappa_t$ is allowed to vary. The range of $\kappa_t$ is set between 0.5 and 1.5 to avoid unphysical areas of the phase space. All the other parameters are fixed to their SM values.

The best fit values (floating $\kappa_\lambda$) of ($c_2$, $\kappa_t$) for the combination of all channels when $\kappa_t$ is allowed to vary. The range of $\kappa_t$ is set between 0.5 and 1.5 to avoid unphysical areas of the phase space. All the other parameters are fixed to their SM values.

The observed 68$\%$ CL cobsevred contours ($\kappa_{t}$ =1.0 (SM)) of $-2\Delta\log(L)$ in the ($c_2$, $\kappa_\lambda$) plane for the combination of all channels when $\kappa_\lambda$ is allowed to vary. The range of $\kappa_\lambda$ is set between -15 and 15 to avoid unphysical areas of the phase space. All the other parameters are fixed to their SM values.

The observed 68$\%$ CL cobsevred contours (floating $\kappa_{t}$) of $-2\Delta\log(L)$ in the ($c_2$, $\kappa_\lambda$) plane for the combination of all channels when $\kappa_\lambda$ is allowed to vary. The range of $\kappa_\lambda$ is set between -15 and 15 to avoid unphysical areas of the phase space. All the other parameters are fixed to their SM values.

The observed 95$\%$ CL cobsevred contours ($\kappa_{t}$ =1.0 (SM)) of $-2\Delta\log(L)$ in the ($c_2$, $\kappa_\lambda$) plane for the combination of all channels when $\kappa_\lambda$ is allowed to vary. The range of $\kappa_\lambda$ is set between -15 and 15 to avoid unphysical areas of the phase space. All the other parameters are fixed to their SM values.

The observed 95$\%$ CL cobsevred contours (floating $\kappa_{t}$) of $-2\Delta\log(L)$ in the ($c_2$, $\kappa_\lambda$) plane for the combination of all channels when $\kappa_\lambda$ is allowed to vary. The range of $\kappa_\lambda$ is set between -15 and 15 to avoid unphysical areas of the phase space. All the other parameters are fixed to their SM values.

The best fit values ($\kappa_{t}$ =1.0 (SM)) of ($c_2$, $\kappa_\lambda$) for the combination of all channels when $\kappa_\lambda$ is allowed to vary. The range of $\kappa_\lambda$ is set between -15 and 15 to avoid unphysical areas of the phase space. All the other parameters are fixed to their SM values.

The best fit values (floating $\kappa_{t}$) for ($c_2$, $\kappa_\lambda$) for the combination of all channels when $\kappa_\lambda$ is allowed to vary. The range of $\kappa_\lambda$ is set between -15 and 15 to avoid unphysical areas of the phase space. All the other parameters are fixed to their SM values.

The observed 68% CL contours of $-2\Delta\log(L)$ in the ($c_2$, $\kappa_\lambda$) plane for the combination of all channels when all the other parameters are fixed to their SM value. The best fit value and 68% CL contour for the $b\bar{b}\tau\tau$ channel are not within the range of the figure.

The observed best fit values in the ($c_2$, $\kappa_\lambda$) plane for the combination of all channels when all the other parameters are fixed to their SM value. The best fit value and 68% CL contour for the $b\bar{b}\tau\tau$ channel are not within the range of the figure.

The observed 68% CL contours of $-2\Delta\log(L)$ in the ($c_2$, $\kappa_\lambda$) plane for the $b\bar{b}b\bar{b}$ resolved channel when all the other parameters are fixed to their SM value. The best fit value and 68% CL contour for the $b\bar{b}\tau\tau$ channel are not within the range of the figure.

The observed best fit values in the ($c_2$, $\kappa_\lambda$) plane for the $b\bar{b}b\bar{b}$ resolved channel when all the other parameters are fixed to their SM value. The best fit value and 68% CL contour for the $b\bar{b}\tau\tau$ channel are not within the range of the figure.

The observed 68% CL contours of $-2\Delta\log(L)$ in the ($c_2$, $\kappa_\lambda$) plane for the $b\bar{b}b\bar{b}$ merged channel when all the other parameters are fixed to their SM value. The best fit value and 68% CL contour for the $b\bar{b}\tau\tau$ channel are not within the range of the figure.

The observed best fit values in the ($c_2$, $\kappa_\lambda$) plane for the $b\bar{b}b\bar{b}$ merged channel when all the other parameters are fixed to their SM value. The best fit value and 68% CL contour for the $b\bar{b}\tau\tau$ channel are not within the range of the figure.

The observed 68% CL contours of $-2\Delta\log(L)$ in the ($c_2$, $\kappa_\lambda$) plane for the $b\bar{b}\gamma\gamma$ channel when all the other parameters are fixed to their SM value. The best fit value and 68% CL contour for the $b\bar{b}\tau\tau$ channel are not within the range of the figure.

The observed best fit values in the ($c_2$, $\kappa_\lambda$) plane for the $b\bar{b}\gamma\gamma$ channel when all the other parameters are fixed to their SM value. The best fit value and 68% CL contour for the $b\bar{b}\tau\tau$ channel are not within the range of the figure.

The observed 68% CL contours of $-2\Delta\log(L)$ in in the ($c_2$, $\kappa_\lambda$) plane for the $b\bar{b}W^{+}W^{-}$ channel when all the other parameters are fixed to their SM value. The best fit value and 68% CL contour for the $b\bar{b}\tau\tau$ channel are not within the range of the figure.

The best fit values in in the ($c_2$, $\kappa_\lambda$) plane for the $b\bar{b}W^{+}W^{-}$ channel when all the other parameters are fixed to their SM value. The best fit value and 68% CL contour for the $b\bar{b}\tau\tau$ channel are not within the range of the figure.

The observed 68% CL contours of $-2\Delta\log(L)$ in in the ($c_2$, $\kappa_\lambda$) plane for the multilepton channel when all the other parameters are fixed to their SM value. The best fit value and 68% CL contour for the $b\bar{b}\tau\tau$ channel are not within the range of the figure.

The best fit values in in the ($c_2$, $\kappa_\lambda$) plane for the multilepton channel when all the other parameters are fixed to their SM value. The best fit value and 68% CL contour for the $b\bar{b}\tau\tau$ channel are not within the range of the figure.

The observed 68% CL contours of $-2\Delta\log(L)$ in in the ($c_2$, $\kappa_\lambda$) plane for the $W^{+}W^{-}\gamma\gamma$ channel when all the other parameters are fixed to their SM value. The best fit value and 68% CL contour for the $b\bar{b}\tau\tau$ channel are not within the range of the figure.

The best fit values in in the ($c_2$, $\kappa_\lambda$) plane for the $W^{+}W^{-}\gamma\gamma$ channel when all the other parameters are fixed to their SM value. The best fit value and 68% CL contour for the $b\bar{b}\tau\tau$ channel are not within the range of the figure.

The $v_{3}\{2,\lvert\Delta\eta\rvert>2\}$ ratios for charged particles as functions of centrality in NeNe to OO collisions at 5.36 TeV.

Signal strength and significance for $\text{t}\overline{\text{t}}\text{Z}(\text{Z}\to\text{b}\overline{\text{b}})$ decays with respect to the standard model expectation of unity.

Signal strength and significance for $\text{t}\overline{\text{t}}\text{Z}(\text{Z}\to\text{c}\overline{\text{c}})$ decays with respect to the standard model expectation of unity.

2D likelihood scans for $\kappa_{b}$ and $\kappa_{c}$ for $\text{t}\overline{\text{t}}\text{H}(\text{H}\to\text{c}\overline{\text{c}})$ and $\text{t}\overline{\text{t}}\text{H}(\text{H}\to\text{b}\overline{\text{b}})$.

The correlation matrix corresponding to the statistical uncertainties on the differential cross-section ($d\sigma/dp_T$) fit results for $W^-$. To combine with $W^+$ results (Fig8a), use the rows and columns ordered as $W^+$ and then $W^-$. Assume no correlation in the statistical uncertainties between $W^+$ and $W^-$ (zero entries in the off-diagonal blocks).

The correlation matrix corresponding to all of the systematic uncertainties on the differential cross-section ($d\sigma/dp_T$) fit results for $W^+$ and $W^-$ combined together, with rows and columns ordered as $W^+$ and then $W^-$.

The $m_{jj}$ and $m_{J}$ projections for the number of observed events (black markers) compared with the backgrounds estimated in the fit to the data (filled histograms) in the CR. Pass and Fail categories are shown. The high level of agreement between the model and the data in the Fail region is due to the nature of the background estimate. The lower panels show the ``Pull'' defined as $(\text{observed events}{-}\text{expected events})/\sqrt{\smash[b]{\sigma_\text{obs}^{2} + \sigma_\text{exp}^{2}}}$, where $\sigma_\text{obs}$ and $\sigma_\text{exp}$ are the total uncertainties in the observation and the background estimation, respectively.

The $m_{jj}$ and $m_{J}$ projections for the number of observed events (black markers) compared with the backgrounds estimated in the fit to the data (filled histograms) in the MR. Pass and Fail categories are shown. The expected contribution from the signal benchmark with $M_{X} = 2200\,\mathrm{GeV}$ and $M_{Y} = 250\,\mathrm{GeV}$ is overlaid in the MR Pass, assuming a production cross section of 5 fb. The high level of agreement between the model and the data in the Fail region is due to the nature of the background estimate.

The $m_{jj}$ and $m_{J}$ projections for the number of observed events (black markers) compared with the backgrounds estimated in the fit to the data (filled histograms) in the MR. Pass and Fail categories are shown. The expected contribution from the signal benchmark with $M_{X} = 2200\,\mathrm{GeV}$ and $M_{Y} = 250\,\mathrm{GeV}$ is overlaid in the MR Pass, assuming a production cross section of 5 fb. The high level of agreement between the model and the data in the Fail region is due to the nature of the background estimate.

The $m_{jj}$ and $m_{J}$ projections for the number of observed events (black markers) compared with the backgrounds estimated in the fit to the data (filled histograms) in the MR. Pass and Fail categories are shown.

The $m_{jj}$ and $m_{J}$ projections for the number of observed events (black markers) compared with the backgrounds estimated in the fit to the data (filled histograms) in the MR. Pass and Fail categories are shown.

The expected and observed 95% confidence level upper limits on the cross section $\sigma(\mathrm{pp} \to \mathrm{X} \to (\mathrm{H} \to b\bar{b})(\mathrm{Y} \to WW \to 4q))$ are shown for different values of $M_X$ and $M_Y$. The limits have been evaluated in discrete steps corresponding to the centers of the boxes. The numbers in the boxes are given in fb. In addition to the observed and median expected limits, the expected limits at $\pm1\sigma$ and $\pm2\sigma$ quantiles are also provided. These represent the regions containing 68% and 95%, respectively, of the distribution of limits expected under the background-only hypothesis.

The median expected (dashed line) and observed (solid line) 95% confidence level upper limits on the main and three alternative signal scenarios as a function of $M_{X}$. The inner (green) band and outer (yellow) band represent the regions containing 68% and 95%, respectively, of the distribution of limits expected under the background-only hypothesis.

The median expected (dashed line) and observed (solid line) 95% confidence level upper limits on the main and three alternative signal scenarios as a function of $M_{X}$. The inner (green) band and outer (yellow) band represent the regions containing 68% and 95%, respectively, of the distribution of limits expected under the background-only hypothesis.

The median expected (dashed line) and observed (solid line) 95% confidence level upper limits on the main and three alternative signal scenarios as a function of $M_{X}$. The inner (green) band and outer (yellow) band represent the regions containing 68% and 95%, respectively, of the distribution of limits expected under the background-only hypothesis.

The median expected (dashed line) and observed (solid line) 95% confidence level upper limits on the main and three alternative signal scenarios as a function of $M_{X}$. The inner (green) band and outer (yellow) band represent the regions containing 68% and 95%, respectively, of the distribution of limits expected under the background-only hypothesis.

Selection efficiency as a function of MX and MY

Selection efficiency as a function of MX for fixed MY

Selection efficiency as a function of MX for fixed MY

Selection efficiency as a function of MX for fixed MY

Response-corrected resolutions of $u_{\parallel}$ vs $N_\mathrm{vtx}$ of different $p_{T}^{miss}$ estimators in data and MC simulations after the $Z\to\mu\mu$ selections.

Response-corrected resolutions of $u_{\perp}$ vs $N_\mathrm{vtx}$ of different $p_{T}^{miss}$ estimators in data and MC simulations after the $Z\to\mu\mu$ selections.

Response-corrected resolutions of $u_{\parallel}$ vs $N_\mathrm{vtx}$ of different $p_{T}^{miss}$ estimators in data after the $Z\to\mu\mu$ selections, in the region with $q_T < 50$ GeV.

Response-corrected resolutions of $u_{\perp}$ vs $N_\mathrm{vtx}$ of different $p_{T}^{miss}$ estimators in data after the $Z\to\mu\mu$ selections, in the region with $q_T < 50$ GeV.

Response-corrected resolutions of $u_{\parallel}$ vs $N_\mathrm{vtx}$ of different $p_{T}^{miss}$ estimators in data after the $Z\to\mu\mu$ selections, in the region with $q_T > 50$ GeV.

Response-corrected resolutions of $u_{\perp}$ vs $N_\mathrm{vtx}$ of different $p_{T}^{miss}$ estimators in data after the $Z\to\mu\mu$ selections, in the region with $q_T > 50$ GeV.

Distribution of $p_{T}^{miss}$ using different $\vec{p}_T^{miss}$ estimators in data after the $W\to\mu\nu$ selections.

Distribution of $m_{T}$ using different $\vec{p}_T^{miss}$ estimators in data after the $W\to\mu\nu$ selections.