Fraction of non-prompt J/$\psi$ in p--Pb collisions at $\sqrt{s_\mathrm{NN}}$ = 8.16 TeV as a function of $p_\mathrm{T}$.

d$^2\sigma$/d$y$d$p_{\rm T}$ in bins of $p_{\mathrm{T}}^{J/\psi}$ for prompt J/$\psi$ in p--Pb collisions at $\sqrt{s_{NN}}$ = 8.16 TeV.

d$^2\sigma$/d$y$d$p_{\rm T}$ in bins of $p_{\mathrm{T}}^{J/\psi}$ for non-prompt J/$\psi$ in p--Pb collisions at $\sqrt{s_{NN}}$ = 8.16 TeV.

Nuclear modification factor ($R_{pPb}$) of prompt J/$\psi$ in p--Pb collisions at $\sqrt{s_{NN}}$ = 8.16 TeV at midrapidity.

Nuclear modification factor ($R_{pPb}$) of non-prompt J/$\psi$ in p--Pb collisions at $\sqrt{s_{NN}}$ = 8.16 TeV at midrapidity.

Measured unfolded differential cross-section as a function of the dilepton transverse momentum $p^{ll}_{\mathrm{T}}$. NLO predictions from Sherpa 2.2.10 and MadGraph5_aMC@NLO 2.7.3 are also shown. The uncertainty in the predictions is divided into statistical and theoretical uncertainties (scale and PDF+$\alpha_{s}$).

Measured unfolded differential cross-section as a function of the the four-body transverse momentum $p^{ll\gamma\gamma}_{\mathrm{T}}$. NLO predictions from Sherpa 2.2.10 and MadGraph5_aMC@NLO 2.7.3 are also shown. The uncertainty in the predictions is divided into statistical and theoretical uncertainties (scale and PDF+$\alpha_{s}$).

Measured unfolded differential cross-section as a function of the diphoton invariant mass $m_{\gamma\gamma}$. NLO predictions from Sherpa 2.2.10 and MadGraph5_aMC@NLO 2.7.3 are also shown. The uncertainty in the predictions is divided into statistical and theoretical uncertainties (scale and PDF+$\alpha_{s}$).

Measured unfolded differential cross-section as a function of the four-body invariant mass $m_{ll\gamma\gamma}$. NLO predictions from Sherpa 2.2.10 and MadGraph5_aMC@NLO 2.7.3 are also shown. The uncertainty in the predictions is divided into statistical and theoretical uncertainties (scale and PDF+$\alpha_{s}$).

Expected and observed $95\%$ confidence intervals for the coupling parameters $f_{T,j}/\Lambda^{4}$ of transverse dimension-8 operators. All parameter values outside of the stated range are excluded at the chosen confidence level. No unitarity constraints are applied.

Expected and observed unitarised $95\%$ confidence intervals for the coupling parameter $f_{T,8}/\Lambda^{4}$ in the clipping energy range between 1.1 and 5 TeV. The non-unitarised limits ($E_c = \infty$) are also shown. All parameter values outside of the stated range are excluded at the chosen confidence level.

Expected and observed unitarised $95\%$ confidence intervals for the coupling parameter $f_{T,0}/\Lambda^{4}$ in the clipping energy range between 1.1 and 5 TeV. The non-unitarised limits ($E_c = \infty$) are also shown. All parameter values outside of the stated range are excluded at the chosen confidence level.

Expected and observed unitarised $95\%$ confidence intervals for the coupling parameter $f_{T,1}/\Lambda^{4}$ in the clipping energy range between 1.1 and 5 TeV. The non-unitarised limits ($E_c = \infty$) are also shown. All parameter values outside of the stated range are excluded at the chosen confidence level.

Expected and observed unitarised $95\%$ confidence intervals for the coupling parameter $f_{T,2}/\Lambda^{4}$ in the clipping energy range between 1.1 and 5 TeV. The non-unitarised limits ($E_c = \infty$) are also shown. All parameter values outside of the stated range are excluded at the chosen confidence level.

Expected and observed unitarised $95\%$ confidence intervals for the coupling parameter $f_{T,5}/\Lambda^{4}$ in the clipping energy range between 1.1 and 5 TeV. The non-unitarised limits ($E_c = \infty$) are also shown. All parameter values outside of the stated range are excluded at the chosen confidence level.

Expected and observed unitarised $95\%$ confidence intervals for the coupling parameter $f_{T,6}/\Lambda^{4}$ in the clipping energy range between 1.1 and 5 TeV. The non-unitarised limits ($E_c = \infty$) are also shown. All parameter values outside of the stated range are excluded at the chosen confidence level.

Expected and observed unitarised $95\%$ confidence intervals for the coupling parameter $f_{T,7}/\Lambda^{4}$ in the clipping energy range between 1.1 and 5 TeV. The non-unitarised limits ($E_c = \infty$) are also shown. All parameter values outside of the stated range are excluded at the chosen confidence level.

Expected and observed unitarised $95\%$ confidence intervals for the coupling parameter $f_{T,9}/\Lambda^{4}$ in the clipping energy range between 1.1 and 5 TeV. The non-unitarised limits ($E_c = \infty$) are also shown. All parameter values outside of the stated range are excluded at the chosen confidence level.

The observed exclusion contour at 95% CL as a function of the $m_{H^{\pm}}$ vs. tan$\beta$ and assuming $m_a$ = 150 $\mathrm{GeV}$, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.7$. Masses that are within the contours are excluded. Only signals simulating the tW+DM final states are considered in this contour.

The expected exclusion contour at 95% CL as a function of the $m_{H^{\pm}}$ vs. tan$\beta$ and assuming $m_a$ = 150 $\mathrm{GeV}$, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.7$. Masses that are within the contours are excluded. Only signals simulating the tW+DM final states are considered in this contour.

The observed exclusion contour at 95% CL as a function of the $m_{H^{\pm}}$ vs. tan$\beta$ and assuming $m_a$ = 250 $\mathrm{GeV}$, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.7$. Masses that are within the contours are excluded. Only signals simulating the tW+DM final states are considered in this contour.

The expected exclusion contour at 95% CL as a function of the $m_{H^{\pm}}$ vs. tan$\beta$ and assuming $m_a$ = 250 $\mathrm{GeV}$, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.7$. Masses that are within the contours are excluded. Only signals simulating the tW+DM final states are considered in this contour.

The observed exclusion contour at 95% CL as a function of the $m_a$ vs. $m_{H^{\pm}}$ and assuming tan$\beta$ = 1, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.7$. Masses that are within the contours are excluded. Signals simulating the tW+DM + tt+DM final states are considered in this contour.

The expected exclusion contour at 95% CL as a function of the $m_a$ vs. $m_{H^{\pm}}$ and assuming tan$\beta$ = 1, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.7$. Masses that are within the contours are excluded. Signals simulating the tW+DM + tt+DM final states are considered in this contour.

The observed exclusion contour at 95% CL as a function of the $m_{H^{\pm}}$ vs. tan$\beta$ and assuming $m_a$ = 150 $\mathrm{GeV}$, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.7$. Masses that are within the contours are excluded. Signals simulating the tW+DM + tt+DM final states are considered in this contour.

The expected exclusion contour at 95% CL as a function of the $m_{H^{\pm}}$ vs. tan$\beta$ and assuming $m_a$ = 150 $\mathrm{GeV}$, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.7$. Masses that are within the contours are excluded. Signals simulating the tW+DM + tt+DM final states are considered in this contour.

The observed exclusion contour at 95% CL as a function of the $m_{H^{\pm}}$ vs. tan$\beta$ and assuming $m_a$ = 250 $\mathrm{GeV}$, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.7$. Masses that are within the contours are excluded. Signals simulating the tW+DM + tt+DM final states are considered in this contour.

The expected exclusion contour at 95% CL as a function of the $m_{H^{\pm}}$ vs. tan$\beta$ and assuming $m_a$ = 250 $\mathrm{GeV}$, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.7$. Masses that are within the contours are excluded. Signals simulating the tW+DM + tt+DM final states are considered in this contour.

The observed exclusion contour at 95% CL as a function of the $m_a$ vs. $m_{H^{\pm}}$ and assuming tan$\beta$ = 1, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.35$. Masses that are within the contours are excluded. Only signals simulating the tW+DM final states are considered in this contour.

The expected exclusion contour at 95% CL as a function of the $m_a$ vs. $m_{H^{\pm}}$ and assuming tan$\beta$ = 1, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.35$. Masses that are within the contours are excluded. Only signals simulating the tW+DM final states are considered in this contour.

The observed exclusion contour at 95% CL as a function of the $m_{H^{\pm}}$ vs. tan$\beta$ and assuming $m_a$ = 150 $\mathrm{GeV}$, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.35$. Masses that are within the contours are excluded. Only signals simulating the tW+DM final states are considered in this contour.

The expected exclusion contour at 95% CL as a function of the $m_{H^{\pm}}$ vs. tan$\beta$ and assuming $m_a$ = 150 $\mathrm{GeV}$, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.35$. Masses that are within the contours are excluded. Only signals simulating the tW+DM final states are considered in this contour.

The observed exclusion contour at 95% CL as a function of the $m_{H^{\pm}}$ vs. tan$\beta$ and assuming $m_a$ = 250 $\mathrm{GeV}$, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.35$. Masses that are within the contours are excluded. Only signals simulating the tW+DM final states are considered in this contour.

The expected exclusion contour at 95% CL as a function of the $m_{H^{\pm}}$ vs. tan$\beta$ and assuming $m_a$ = 250 $\mathrm{GeV}$, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.35$. Masses that are within the contours are excluded. Only signals simulating the tW+DM final states are considered in this contour.

The observed exclusion contour at 95% CL as a function of the $m_a$ vs. $m_{H^{\pm}}$ and assuming tan$\beta$ = 1, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.35$. Masses that are within the contours are excluded. Signals simulating the tW+DM + tt+DM final states are considered in this contour.

The expected exclusion contour at 95% CL as a function of the $m_a$ vs. $m_{H^{\pm}}$ and assuming tan$\beta$ = 1, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.35$. Masses that are within the contours are excluded. Signals simulating the tW+DM + tt+DM final states are considered in this contour.

The observed exclusion contour at 95% CL as a function of the $m_{H^{\pm}}$ vs. tan$\beta$ and assuming $m_a$ = 150 $\mathrm{GeV}$, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.35$. Masses that are within the contours are excluded. Signals simulating the tW+DM + tt+DM final states are considered in this contour.

The expected exclusion contour at 95% CL as a function of the $m_{H^{\pm}}$ vs. tan$\beta$ and assuming $m_a$ = 150 $\mathrm{GeV}$, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.35$. Masses that are within the contours are excluded. Signals simulating the tW+DM + tt+DM final states are considered in this contour.

The observed exclusion contour at 95% CL as a function of the $m_{H^{\pm}}$ vs. tan$\beta$ and assuming $m_a$ = 250 $\mathrm{GeV}$, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.35$. Masses that are within the contours are excluded. Signals simulating the tW+DM + tt+DM final states are considered in this contour.

The expected exclusion contour at 95% CL as a function of the $m_{H^{\pm}}$ vs. tan$\beta$ and assuming $m_a$ = 250 $\mathrm{GeV}$, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.35$. Masses that are within the contours are excluded. Signals simulating the tW+DM + tt+DM final states are considered in this contour.

The observed exclusion contour at 95% CL as a function of the $m_a$ vs. $m_{H^{\pm}}$ and assuming tan$\beta$ = 1, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.7$. Masses that are within the contours are excluded. Only signals simulating the tW+DM final states are considered in this contour. These exclusion contours are derived using the 0L channel only.

The expected exclusion contour at 95% CL as a function of the $m_a$ vs. $m_{H^{\pm}}$ and assuming tan$\beta$ = 1, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.7$. Masses that are within the contours are excluded. Only signals simulating the tW+DM final states are considered in this contour. These exclusion contours are derived using the 0L channel only.

The observed exclusion contour at 95% CL as a function of the $m_{H^{\pm}}$ vs. tan$\beta$ and assuming $m_a$ = 150 $\mathrm{GeV}$, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.7$. Masses that are within the contours are excluded. Only signals simulating the tW+DM final states are considered in this contour. These exclusion contours are derived using the 0L channel only.

The expected exclusion contour at 95% CL as a function of the $m_{H^{\pm}}$ vs. tan$\beta$ and assuming $m_a$ = 150 $\mathrm{GeV}$, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.7$. Masses that are within the contours are excluded. Only signals simulating the tW+DM final states are considered in this contour. These exclusion contours are derived using the 0L channel only.

The observed exclusion contour at 95% CL as a function of the $m_{H^{\pm}}$ vs. tan$\beta$ and assuming $m_a$ = 250 $\mathrm{GeV}$, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.7$. Masses that are within the contours are excluded. Only signals simulating the tW+DM final states are considered in this contour. These exclusion contours are derived using the 0L channel only.

The expected exclusion contour at 95% CL as a function of the $m_{H^{\pm}}$ vs. tan$\beta$ and assuming $m_a$ = 250 $\mathrm{GeV}$, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.7$. Masses that are within the contours are excluded. Only signals simulating the tW+DM final states are considered in this contour. These exclusion contours are derived using the 0L channel only.

The observed exclusion contour at 95% CL as a function of the $m_a$ vs. $m_{H^{\pm}}$ and assuming tan$\beta$ = 1, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.7$. Masses that are within the contours are excluded. Signals simulating the tW+DM + tt+DM final states are considered in this contour. These exclusion contours are derived using the 0L channel only.

The expected exclusion contour at 95% CL as a function of the $m_a$ vs. $m_{H^{\pm}}$ and assuming tan$\beta$ = 1, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.7$. Masses that are within the contours are excluded. Signals simulating the tW+DM + tt+DM final states are considered in this contour. These exclusion contours are derived using the 0L channel only.

The observed exclusion contour at 95% CL as a function of the $m_{H^{\pm}}$ vs. tan$\beta$ and assuming $m_a$ = 150 $\mathrm{GeV}$, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.7$. Masses that are within the contours are excluded. Signals simulating the tW+DM + tt+DM final states are considered in this contour. These exclusion contours are derived using the 0L channel only.

The expected exclusion contour at 95% CL as a function of the $m_{H^{\pm}}$ vs. tan$\beta$ and assuming $m_a$ = 150 $\mathrm{GeV}$, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.7$. Masses that are within the contours are excluded. Signals simulating the tW+DM + tt+DM final states are considered in this contour. These exclusion contours are derived using the 0L channel only.

The observed exclusion contour at 95% CL as a function of the $m_{H^{\pm}}$ vs. tan$\beta$ and assuming $m_a$ = 250 $\mathrm{GeV}$, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.7$. Masses that are within the contours are excluded. Signals simulating the tW+DM + tt+DM final states are considered in this contour. These exclusion contours are derived using the 0L channel only.

The expected exclusion contour at 95% CL as a function of the $m_{H^{\pm}}$ vs. tan$\beta$ and assuming $m_a$ = 250 $\mathrm{GeV}$, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.7$. Masses that are within the contours are excluded. Signals simulating the tW+DM + tt+DM final states are considered in this contour. These exclusion contours are derived using the 0L channel only.

The observed exclusion contour at 95% CL as a function of the $m_a$ vs. $m_{H^{\pm}}$ and assuming tan$\beta$ = 1, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.35$. Masses that are within the contours are excluded. Only signals simulating the tW+DM final states are considered in this contour. These exclusion contours are derived using the 0L channel only.

The expected exclusion contour at 95% CL as a function of the $m_a$ vs. $m_{H^{\pm}}$ and assuming tan$\beta$ = 1, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.35$. Masses that are within the contours are excluded. Only signals simulating the tW+DM final states are considered in this contour. These exclusion contours are derived using the 0L channel only.

The observed exclusion contour at 95% CL as a function of the $m_{H^{\pm}}$ vs. tan$\beta$ and assuming $m_a$ = 150 $\mathrm{GeV}$, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.35$. Masses that are within the contours are excluded. Only signals simulating the tW+DM final states are considered in this contour. These exclusion contours are derived using the 0L channel only.

The expected exclusion contour at 95% CL as a function of the $m_{H^{\pm}}$ vs. tan$\beta$ and assuming $m_a$ = 150 $\mathrm{GeV}$, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.35$. Masses that are within the contours are excluded. Only signals simulating the tW+DM final states are considered in this contour. These exclusion contours are derived using the 0L channel only.

The observed exclusion contour at 95% CL as a function of the $m_{H^{\pm}}$ vs. tan$\beta$ and assuming $m_a$ = 250 $\mathrm{GeV}$, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.35$. Masses that are within the contours are excluded. Only signals simulating the tW+DM final states are considered in this contour. These exclusion contours are derived using the 0L channel only.

The expected exclusion contour at 95% CL as a function of the $m_{H^{\pm}}$ vs. tan$\beta$ and assuming $m_a$ = 250 $\mathrm{GeV}$, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.35$. Masses that are within the contours are excluded. Only signals simulating the tW+DM final states are considered in this contour. These exclusion contours are derived using the 0L channel only.

The observed exclusion contour at 95% CL as a function of the $m_a$ vs. $m_{H^{\pm}}$ and assuming tan$\beta$ = 1, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.35$. Masses that are within the contours are excluded. Signals simulating the tW+DM + tt+DM final states are considered in this contour. These exclusion contours are derived using the 0L channel only.

The expected exclusion contour at 95% CL as a function of the $m_a$ vs. $m_{H^{\pm}}$ and assuming tan$\beta$ = 1, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.35$. Masses that are within the contours are excluded. Signals simulating the tW+DM + tt+DM final states are considered in this contour. These exclusion contours are derived using the 0L channel only.

The observed exclusion contour at 95% CL as a function of the $m_{H^{\pm}}$ vs. tan$\beta$ and assuming $m_a$ = 150 $\mathrm{GeV}$, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.35$. Masses that are within the contours are excluded. Signals simulating the tW+DM + tt+DM final states are considered in this contour. These exclusion contours are derived using the 0L channel only.

The expected exclusion contour at 95% CL as a function of the $m_{H^{\pm}}$ vs. tan$\beta$ and assuming $m_a$ = 150 $\mathrm{GeV}$, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.35$. Masses that are within the contours are excluded. Signals simulating the tW+DM + tt+DM final states are considered in this contour. These exclusion contours are derived using the 0L channel only.

The observed exclusion contour at 95% CL as a function of the $m_{H^{\pm}}$ vs. tan$\beta$ and assuming $m_a$ = 250 $\mathrm{GeV}$, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.35$. Masses that are within the contours are excluded. Signals simulating the tW+DM + tt+DM final states are considered in this contour. These exclusion contours are derived using the 0L channel only.

The expected exclusion contour at 95% CL as a function of the $m_{H^{\pm}}$ vs. tan$\beta$ and assuming $m_a$ = 250 $\mathrm{GeV}$, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.35$. Masses that are within the contours are excluded. Signals simulating the tW+DM + tt+DM final states are considered in this contour. These exclusion contours are derived using the 0L channel only.

The observed exclusion contour at 95% CL as a function of the $m_a$ vs. $m_{H^{\pm}}$ and assuming tan$\beta$ = 1, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.7$. Masses that are within the contours are excluded. Only signals simulating the tW+DM final states are considered in this contour. These exclusion contours are derived using the 1L channel only.

The expected exclusion contour at 95% CL as a function of the $m_a$ vs. $m_{H^{\pm}}$ and assuming tan$\beta$ = 1, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.7$. Masses that are within the contours are excluded. Only signals simulating the tW+DM final states are considered in this contour. These exclusion contours are derived using the 1L channel only.

The observed exclusion contour at 95% CL as a function of the $m_{H^{\pm}}$ vs. tan$\beta$ and assuming $m_a$ = 150 $\mathrm{GeV}$, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.7$. Masses that are within the contours are excluded. Only signals simulating the tW+DM final states are considered in this contour. These exclusion contours are derived using the 1L channel only.

The expected exclusion contour at 95% CL as a function of the $m_{H^{\pm}}$ vs. tan$\beta$ and assuming $m_a$ = 150 $\mathrm{GeV}$, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.7$. Masses that are within the contours are excluded. Only signals simulating the tW+DM final states are considered in this contour. These exclusion contours are derived using the 1L channel only.

The observed exclusion contour at 95% CL as a function of the $m_{H^{\pm}}$ vs. tan$\beta$ and assuming $m_a$ = 250 $\mathrm{GeV}$, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.7$. Masses that are within the contours are excluded. Only signals simulating the tW+DM final states are considered in this contour. These exclusion contours are derived using the 1L channel only.

The expected exclusion contour at 95% CL as a function of the $m_{H^{\pm}}$ vs. tan$\beta$ and assuming $m_a$ = 250 $\mathrm{GeV}$, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.7$. Masses that are within the contours are excluded. Only signals simulating the tW+DM final states are considered in this contour. These exclusion contours are derived using the 1L channel only.

The observed exclusion contour at 95% CL as a function of the $m_a$ vs. $m_{H^{\pm}}$ and assuming tan$\beta$ = 1, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.7$. Masses that are within the contours are excluded. Signals simulating the tW+DM + tt+DM final states are considered in this contour. These exclusion contours are derived using the 1L channel only.

The expected exclusion contour at 95% CL as a function of the $m_a$ vs. $m_{H^{\pm}}$ and assuming tan$\beta$ = 1, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.7$. Masses that are within the contours are excluded. Signals simulating the tW+DM + tt+DM final states are considered in this contour. These exclusion contours are derived using the 1L channel only.

The observed exclusion contour at 95% CL as a function of the $m_{H^{\pm}}$ vs. tan$\beta$ and assuming $m_a$ = 150 $\mathrm{GeV}$, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.7$. Masses that are within the contours are excluded. Signals simulating the tW+DM + tt+DM final states are considered in this contour. These exclusion contours are derived using the 1L channel only.

The expected exclusion contour at 95% CL as a function of the $m_{H^{\pm}}$ vs. tan$\beta$ and assuming $m_a$ = 150 $\mathrm{GeV}$, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.7$. Masses that are within the contours are excluded. Signals simulating the tW+DM + tt+DM final states are considered in this contour. These exclusion contours are derived using the 1L channel only.

The observed exclusion contour at 95% CL as a function of the $m_{H^{\pm}}$ vs. tan$\beta$ and assuming $m_a$ = 250 $\mathrm{GeV}$, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.7$. Masses that are within the contours are excluded. Signals simulating the tW+DM + tt+DM final states are considered in this contour. These exclusion contours are derived using the 1L channel only.

The expected exclusion contour at 95% CL as a function of the $m_{H^{\pm}}$ vs. tan$\beta$ and assuming $m_a$ = 250 $\mathrm{GeV}$, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.7$. Masses that are within the contours are excluded. Signals simulating the tW+DM + tt+DM final states are considered in this contour. These exclusion contours are derived using the 1L channel only.

The observed exclusion contour at 95% CL as a function of the $m_a$ vs. $m_{H^{\pm}}$ and assuming tan$\beta$ = 1, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.35$. Masses that are within the contours are excluded. Only signals simulating the tW+DM final states are considered in this contour. These exclusion contours are derived using the 1L channel only.

The expected exclusion contour at 95% CL as a function of the $m_a$ vs. $m_{H^{\pm}}$ and assuming tan$\beta$ = 1, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.35$. Masses that are within the contours are excluded. Only signals simulating the tW+DM final states are considered in this contour. These exclusion contours are derived using the 1L channel only.

The observed exclusion contour at 95% CL as a function of the $m_{H^{\pm}}$ vs. tan$\beta$ and assuming $m_a$ = 150 $\mathrm{GeV}$, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.35$. Masses that are within the contours are excluded. Only signals simulating the tW+DM final states are considered in this contour. These exclusion contours are derived using the 1L channel only.

The expected exclusion contour at 95% CL as a function of the $m_{H^{\pm}}$ vs. tan$\beta$ and assuming $m_a$ = 150 $\mathrm{GeV}$, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.35$. Masses that are within the contours are excluded. Only signals simulating the tW+DM final states are considered in this contour. These exclusion contours are derived using the 1L channel only.

The expected exclusion contour at 95% CL as a function of the $m_{H^{\pm}}$ vs. tan$\beta$ and assuming $m_a$ = 250 $\mathrm{GeV}$, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.35$. Masses that are within the contours are excluded. Only signals simulating the tW+DM final states are considered in this contour. These exclusion contours are derived using the 1L channel only.

The observed exclusion contour at 95% CL as a function of the $m_a$ vs. $m_{H^{\pm}}$ and assuming tan$\beta$ = 1, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.35$. Masses that are within the contours are excluded. Signals simulating the tW+DM + tt+DM final states are considered in this contour. These exclusion contours are derived using the 1L channel only.

The expected exclusion contour at 95% CL as a function of the $m_a$ vs. $m_{H^{\pm}}$ and assuming tan$\beta$ = 1, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.35$. Masses that are within the contours are excluded. Signals simulating the tW+DM + tt+DM final states are considered in this contour. These exclusion contours are derived using the 1L channel only.

The observed exclusion contour at 95% CL as a function of the $m_{H^{\pm}}$ vs. tan$\beta$ and assuming $m_a$ = 150 $\mathrm{GeV}$, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.35$. Masses that are within the contours are excluded. Signals simulating the tW+DM + tt+DM final states are considered in this contour. These exclusion contours are derived using the 1L channel only.

The expected exclusion contour at 95% CL as a function of the $m_{H^{\pm}}$ vs. tan$\beta$ and assuming $m_a$ = 150 $\mathrm{GeV}$, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.35$. Masses that are within the contours are excluded. Signals simulating the tW+DM + tt+DM final states are considered in this contour. These exclusion contours are derived using the 1L channel only.

The observed exclusion contour at 95% CL as a function of the $m_{H^{\pm}}$ vs. tan$\beta$ and assuming $m_a$ = 250 $\mathrm{GeV}$, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.35$. Masses that are within the contours are excluded. Signals simulating the tW+DM + tt+DM final states are considered in this contour. These exclusion contours are derived using the 1L channel only.

The expected exclusion contour at 95% CL as a function of the $m_{H^{\pm}}$ vs. tan$\beta$ and assuming $m_a$ = 250 $\mathrm{GeV}$, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.35$. Masses that are within the contours are excluded. Signals simulating the tW+DM + tt+DM final states are considered in this contour. These exclusion contours are derived using the 1L channel only.

The observed exclusion contour at 95% CL as a function of the $m_a$ vs. $m_{H^{\pm}}$ and assuming tan$\beta$ = 1, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.7$. Masses that are within the contours are excluded. Only signals simulating the tW+DM final states are considered in this contour. These exclusion contours are derived using the 2L channel only.

The expected exclusion contour at 95% CL as a function of the $m_a$ vs. $m_{H^{\pm}}$ and assuming tan$\beta$ = 1, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.7$. Masses that are within the contours are excluded. Only signals simulating the tW+DM final states are considered in this contour. These exclusion contours are derived using the 2L channel only.

The observed exclusion contour at 95% CL as a function of the $m_{H^{\pm}}$ vs. tan$\beta$ and assuming $m_a$ = 150 $\mathrm{GeV}$, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.7$. Masses that are within the contours are excluded. Only signals simulating the tW+DM final states are considered in this contour. These exclusion contours are derived using the 2L channel only.

The expected exclusion contour at 95% CL as a function of the $m_{H^{\pm}}$ vs. tan$\beta$ and assuming $m_a$ = 150 $\mathrm{GeV}$, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.7$. Masses that are within the contours are excluded. Only signals simulating the tW+DM final states are considered in this contour. These exclusion contours are derived using the 2L channel only.

The observed exclusion contour at 95% CL as a function of the $m_{H^{\pm}}$ vs. tan$\beta$ and assuming $m_a$ = 250 $\mathrm{GeV}$, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.7$. Masses that are within the contours are excluded. Only signals simulating the tW+DM final states are considered in this contour. These exclusion contours are derived using the 2L channel only.

The expected exclusion contour at 95% CL as a function of the $m_{H^{\pm}}$ vs. tan$\beta$ and assuming $m_a$ = 250 $\mathrm{GeV}$, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.7$. Masses that are within the contours are excluded. Only signals simulating the tW+DM final states are considered in this contour. These exclusion contours are derived using the 2L channel only.

The observed exclusion contour at 95% CL as a function of the $m_a$ vs. $m_{H^{\pm}}$ and assuming tan$\beta$ = 1, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.7$. Masses that are within the contours are excluded. Signals simulating the tW+DM + tt+DM final states are considered in this contour. These exclusion contours are derived using the 2L channel only.

The expected exclusion contour at 95% CL as a function of the $m_a$ vs. $m_{H^{\pm}}$ and assuming tan$\beta$ = 1, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.7$. Masses that are within the contours are excluded. Signals simulating the tW+DM + tt+DM final states are considered in this contour. These exclusion contours are derived using the 2L channel only.

The observed exclusion contour at 95% CL as a function of the $m_{H^{\pm}}$ vs. tan$\beta$ and assuming $m_a$ = 150 $\mathrm{GeV}$, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.7$. Masses that are within the contours are excluded. Signals simulating the tW+DM + tt+DM final states are considered in this contour. These exclusion contours are derived using the 2L channel only.

The expected exclusion contour at 95% CL as a function of the $m_{H^{\pm}}$ vs. tan$\beta$ and assuming $m_a$ = 150 $\mathrm{GeV}$, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.7$. Masses that are within the contours are excluded. Signals simulating the tW+DM + tt+DM final states are considered in this contour. These exclusion contours are derived using the 2L channel only.

The observed exclusion contour at 95% CL as a function of the $m_{H^{\pm}}$ vs. tan$\beta$ and assuming $m_a$ = 250 $\mathrm{GeV}$, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.7$. Masses that are within the contours are excluded. Signals simulating the tW+DM + tt+DM final states are considered in this contour. These exclusion contours are derived using the 2L channel only.

The expected exclusion contour at 95% CL as a function of the $m_{H^{\pm}}$ vs. tan$\beta$ and assuming $m_a$ = 250 $\mathrm{GeV}$, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.7$. Masses that are within the contours are excluded. Signals simulating the tW+DM + tt+DM final states are considered in this contour. These exclusion contours are derived using the 2L channel only.

The expected exclusion contour at 95% CL as a function of the $m_a$ vs. $m_{H^{\pm}}$ and assuming tan$\beta$ = 1, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.35$. Masses that are within the contours are excluded. Only signals simulating the tW+DM final states are considered in this contour. These exclusion contours are derived using the 2L channel only.

The expected exclusion contour at 95% CL as a function of the $m_{H^{\pm}}$ vs. tan$\beta$ and assuming $m_a$ = 150 $\mathrm{GeV}$, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.35$. Masses that are within the contours are excluded. Only signals simulating the tW+DM final states are considered in this contour. These exclusion contours are derived using the 2L channel only.

The expected exclusion contour at 95% CL as a function of the $m_{H^{\pm}}$ vs. tan$\beta$ and assuming $m_a$ = 250 $\mathrm{GeV}$, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.35$. Masses that are within the contours are excluded. Only signals simulating the tW+DM final states are considered in this contour. These exclusion contours are derived using the 2L channel only.

The expected exclusion contour at 95% CL as a function of the $m_a$ vs. $m_{H^{\pm}}$ and assuming tan$\beta$ = 1, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.35$. Masses that are within the contours are excluded. Signals simulating the tW+DM + tt+DM final states are considered in this contour. These exclusion contours are derived using the 2L channel only.

The observed exclusion contour at 95% CL as a function of the $m_{H^{\pm}}$ vs. tan$\beta$ and assuming $m_a$ = 150 $\mathrm{GeV}$, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.35$. Masses that are within the contours are excluded. Signals simulating the tW+DM + tt+DM final states are considered in this contour. These exclusion contours are derived using the 2L channel only.

The expected exclusion contour at 95% CL as a function of the $m_{H^{\pm}}$ vs. tan$\beta$ and assuming $m_a$ = 150 $\mathrm{GeV}$, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.35$. Masses that are within the contours are excluded. Signals simulating the tW+DM + tt+DM final states are considered in this contour. These exclusion contours are derived using the 2L channel only.

The observed exclusion contour at 95% CL as a function of the $m_{H^{\pm}}$ vs. tan$\beta$ and assuming $m_a$ = 250 $\mathrm{GeV}$, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.35$. Masses that are within the contours are excluded. Signals simulating the tW+DM + tt+DM final states are considered in this contour. These exclusion contours are derived using the 2L channel only.

The expected exclusion contour at 95% CL as a function of the $m_{H^{\pm}}$ vs. tan$\beta$ and assuming $m_a$ = 250 $\mathrm{GeV}$, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.35$. Masses that are within the contours are excluded. Signals simulating the tW+DM + tt+DM final states are considered in this contour. These exclusion contours are derived using the 2L channel only.

Model dependent upper limit on the cross section for the $m_a$ vs. $ m_{H^{\pm}}$ signal grid assuming tan$\beta$ = 1, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.7$. Only signals simulating the tW+DM final states are considered. Upper limits with large $\mu_{\mathrm{sig}}$ for the observed limit are capped at 500.

Model dependent upper limit on the cross section for the $m_{H^{\pm}}$ vs. tan$\beta$ signal grid assuming $m_a$ = 150 GeV, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.7$. Only signals simulating the tW+DM final states are considered. Upper limits with large $\mu_{\mathrm{sig}}$ for the observed limit are capped at 500.

Model dependent upper limit on the cross section for the $m_{H^{\pm}}$ vs. tan$\beta$ signal grid assuming $m_a$ = 250 GeV, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.7$. Only signals simulating the tW+DM final states are considered. Upper limits with large $\mu_{\mathrm{sig}}$ for the observed limit are capped at 500.

Model dependent upper limit on the cross section for the $m_a$ vs. $ m_{H^{\pm}}$ signal grid assuming tan$\beta$ = 1, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.7$. Signals simulating the tW+DM + tt+DM final states are considered. Upper limits with large $\mu_{\mathrm{sig}}$ for the observed limit are capped at 500.

Model dependent upper limit on the cross section for the $m_{H^{\pm}}$ vs. tan$\beta$ signal grid assuming $m_a$ = 150 GeV, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.7$. Signals simulating the tW+DM + tt+DM final states are considered. Upper limits with large $\mu_{\mathrm{sig}}$ for the observed limit are capped at 500.

Model dependent upper limit on the cross section for the $m_{H^{\pm}}$ vs. tan$\beta$ signal grid assuming $m_a$ = 250 GeV, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.7$. Signals simulating the tW+DM + tt+DM final states are considered. Upper limits with large $\mu_{\mathrm{sig}}$ for the observed limit are capped at 500.

Model dependent upper limit on the cross section for the $m_a$ vs. $ m_{H^{\pm}}$ signal grid assuming tan$\beta$ = 1, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.35$. Only signals simulating the tW+DM final states are considered. Upper limits with large $\mu_{\mathrm{sig}}$ for the observed limit are capped at 500.

Model dependent upper limit on the cross section for the $m_{H^{\pm}}$ vs. tan$\beta$ signal grid assuming $m_a$ = 150 GeV, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.35$. Only signals simulating the tW+DM final states are considered. Upper limits with large $\mu_{\mathrm{sig}}$ for the observed limit are capped at 500.

Model dependent upper limit on the cross section for the $m_{H^{\pm}}$ vs. tan$\beta$ signal grid assuming $m_a$ = 250 GeV, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.35$. Only signals simulating the tW+DM final states are considered. Upper limits with large $\mu_{\mathrm{sig}}$ for the observed limit are capped at 500.

Model dependent upper limit on the cross section for the $m_a$ vs. $ m_{H^{\pm}}$ signal grid assuming tan$\beta$ = 1, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.35$. Signals simulating the tW+DM + tt+DM final states are considered. Upper limits with large $\mu_{\mathrm{sig}}$ for the observed limit are capped at 500.

Model dependent upper limit on the cross section for the $m_{H^{\pm}}$ vs. tan$\beta$ signal grid assuming $m_a$ = 150 GeV, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.35$. Signals simulating the tW+DM + tt+DM final states are considered. Upper limits with large $\mu_{\mathrm{sig}}$ for the observed limit are capped at 500.

Model dependent upper limit on the cross section for the $m_{H^{\pm}}$ vs. tan$\beta$ signal grid assuming $m_a$ = 250 GeV, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.35$. Signals simulating the tW+DM + tt+DM final states are considered. Upper limits with large $\mu_{\mathrm{sig}}$ for the observed limit are capped at 500.

Model dependent upper limit on the cross section for the $m_a$ vs. $ m_{H^{\pm}}$ signal grid assuming tan$\beta$ = 1, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.7$. Only signals simulating the tW+DM final states are considered. Upper limits with large $\mu_{\mathrm{sig}}$ for the observed limit are capped at 500.

Model dependent upper limit on the cross section for the $m_{H^{\pm}}$ vs. tan$\beta$ signal grid assuming $m_a$ = 150 GeV, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.7$. Only signals simulating the tW+DM final states are considered. Upper limits with large $\mu_{\mathrm{sig}}$ for the observed limit are capped at 500.

Model dependent upper limit on the cross section for the $m_{H^{\pm}}$ vs. tan$\beta$ signal grid assuming $m_a$ = 250 GeV, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.7$. Only signals simulating the tW+DM final states are considered. Upper limits with large $\mu_{\mathrm{sig}}$ for the observed limit are capped at 500.

Model dependent upper limit on the cross section for the $m_a$ vs. $ m_{H^{\pm}}$ signal grid assuming tan$\beta$ = 1, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.7$. Signals simulating the tW+DM + tt+DM final states are considered. Upper limits with large $\mu_{\mathrm{sig}}$ for the observed limit are capped at 500.

Model dependent upper limit on the cross section for the $m_{H^{\pm}}$ vs. tan$\beta$ signal grid assuming $m_a$ = 150 GeV, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.7$. Signals simulating the tW+DM + tt+DM final states are considered. Upper limits with large $\mu_{\mathrm{sig}}$ for the observed limit are capped at 500.

Model dependent upper limit on the cross section for the $m_{H^{\pm}}$ vs. tan$\beta$ signal grid assuming $m_a$ = 250 GeV, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.7$. Signals simulating the tW+DM + tt+DM final states are considered. Upper limits with large $\mu_{\mathrm{sig}}$ for the observed limit are capped at 500.

Model dependent upper limit on the cross section for the $m_a$ vs. $ m_{H^{\pm}}$ signal grid assuming tan$\beta$ = 1, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.35$. Only signals simulating the tW+DM final states are considered. Upper limits with large $\mu_{\mathrm{sig}}$ for the observed limit are capped at 500.

Model dependent upper limit on the cross section for the $m_{H^{\pm}}$ vs. tan$\beta$ signal grid assuming $m_a$ = 150 GeV, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.35$. Only signals simulating the tW+DM final states are considered. Upper limits with large $\mu_{\mathrm{sig}}$ for the observed limit are capped at 500.

Model dependent upper limit on the cross section for the $m_{H^{\pm}}$ vs. tan$\beta$ signal grid assuming $m_a$ = 250 GeV, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.35$. Only signals simulating the tW+DM final states are considered. Upper limits with large $\mu_{\mathrm{sig}}$ for the observed limit are capped at 500.

Model dependent upper limit on the cross section for the $m_a$ vs. $ m_{H^{\pm}}$ signal grid assuming tan$\beta$ = 1, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.35$. Signals simulating the tW+DM + tt+DM final states are considered. Upper limits with large $\mu_{\mathrm{sig}}$ for the observed limit are capped at 500.

Model dependent upper limit on the cross section for the $m_{H^{\pm}}$ vs. tan$\beta$ signal grid assuming $m_a$ = 150 GeV, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.35$. Signals simulating the tW+DM + tt+DM final states are considered. Upper limits with large $\mu_{\mathrm{sig}}$ for the observed limit are capped at 500.

Model dependent upper limit on the cross section for the $m_{H^{\pm}}$ vs. tan$\beta$ signal grid assuming $m_a$ = 250 GeV, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.35$. Signals simulating the tW+DM + tt+DM final states are considered. Upper limits with large $\mu_{\mathrm{sig}}$ for the observed limit are capped at 500.

Model dependent upper limit on the cross section for the $m_a$ vs. $ m_{H^{\pm}}$ signal grid assuming tan$\beta$ = 1, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.7$. Only signals simulating the tW+DM final states are considered. Upper limits with large $\mu_{\mathrm{sig}}$ for the observed limit are capped at 500.

Model dependent upper limit on the cross section for the $m_{H^{\pm}}$ vs. tan$\beta$ signal grid assuming $m_a$ = 150 GeV, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.7$. Only signals simulating the tW+DM final states are considered. Upper limits with large $\mu_{\mathrm{sig}}$ for the observed limit are capped at 500.

Model dependent upper limit on the cross section for the $m_{H^{\pm}}$ vs. tan$\beta$ signal grid assuming $m_a$ = 250 GeV, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.7$. Only signals simulating the tW+DM final states are considered. Upper limits with large $\mu_{\mathrm{sig}}$ for the observed limit are capped at 500.

Model dependent upper limit on the cross section for the $m_a$ vs. $ m_{H^{\pm}}$ signal grid assuming tan$\beta$ = 1, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.7$. Signals simulating the tW+DM + tt+DM final states are considered. Upper limits with large $\mu_{\mathrm{sig}}$ for the observed limit are capped at 500.

Model dependent upper limit on the cross section for the $m_{H^{\pm}}$ vs. tan$\beta$ signal grid assuming $m_a$ = 150 GeV, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.7$. Signals simulating the tW+DM + tt+DM final states are considered. Upper limits with large $\mu_{\mathrm{sig}}$ for the observed limit are capped at 500.

Model dependent upper limit on the cross section for the $m_{H^{\pm}}$ vs. tan$\beta$ signal grid assuming $m_a$ = 250 GeV, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.7$. Signals simulating the tW+DM + tt+DM final states are considered. Upper limits with large $\mu_{\mathrm{sig}}$ for the observed limit are capped at 500.

Model dependent upper limit on the cross section for the $m_a$ vs. $ m_{H^{\pm}}$ signal grid assuming tan$\beta$ = 1, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.35$. Only signals simulating the tW+DM final states are considered. Upper limits with large $\mu_{\mathrm{sig}}$ for the observed limit are capped at 500.

Model dependent upper limit on the cross section for the $m_{H^{\pm}}$ vs. tan$\beta$ signal grid assuming $m_a$ = 150 GeV, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.35$. Only signals simulating the tW+DM final states are considered. Upper limits with large $\mu_{\mathrm{sig}}$ for the observed limit are capped at 500.

Model dependent upper limit on the cross section for the $m_{H^{\pm}}$ vs. tan$\beta$ signal grid assuming $m_a$ = 250 GeV, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.35$. Only signals simulating the tW+DM final states are considered. Upper limits with large $\mu_{\mathrm{sig}}$ for the observed limit are capped at 500.

Model dependent upper limit on the cross section for the $m_a$ vs. $ m_{H^{\pm}}$ signal grid assuming tan$\beta$ = 1, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.35$. Signals simulating the tW+DM + tt+DM final states are considered. Upper limits with large $\mu_{\mathrm{sig}}$ for the observed limit are capped at 500.

Model dependent upper limit on the cross section for the $m_{H^{\pm}}$ vs. tan$\beta$ signal grid assuming $m_a$ = 150 GeV, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.35$. Signals simulating the tW+DM + tt+DM final states are considered. Upper limits with large $\mu_{\mathrm{sig}}$ for the observed limit are capped at 500.

Model dependent upper limit on the cross section for the $m_{H^{\pm}}$ vs. tan$\beta$ signal grid assuming $m_a$ = 250 GeV, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.35$. Signals simulating the tW+DM + tt+DM final states are considered. Upper limits with large $\mu_{\mathrm{sig}}$ for the observed limit are capped at 500.

The distributions of $m_{\mathrm{b1},\mathrm{W-tagged}}$ in the 0L inclusive signal region. For each bin yields for the data and total SM prediction are provided. The SM prediction is provided with the total uncertainty, including the MC statistical uncertainty, detector-related systematic uncertainties and theoretical uncertainties. The rightmost bin includes overflow events.

The distributions of $m_{\mathrm{T}}^{\mathrm{b,E_{\mathrm{T}^{\mathrm{miss}}}}}$ in the 0L inclusive signal region. For each bin yields for the data and total SM prediction are provided. The SM prediction is provided with the total uncertainty, including the MC statistical uncertainty, detector-related systematic uncertainties and theoretical uncertainties. The rightmost bin includes overflow events.

The distributions of $N_{\mathrm{W-tagged}}$ in the 0L inclusive signal region. For each bin yields for the data and total SM prediction are provided. The SM prediction is provided with the total uncertainty, including the MC statistical uncertainty, detector-related systematic uncertainties and theoretical uncertainties. The rightmost bin includes overflow events.

The distributions of $m_{\mathrm{b1},\mathrm{\cancel{b1}}}$ in the hadronic top inclusive signal region. For each bin yields for the data and total SM prediction are provided. The SM prediction is provided with the total uncertainty, including the MC statistical uncertainty, detector-related systematic uncertainties and theoretical uncertainties. The rightmost bin includes overflow events.

The distributions of $m_{\mathrm{b1},\mathrm{\cancel{b1}}}$ in the leptonic top inclusive signal region. For each bin yields for the data and total SM prediction are provided. The SM prediction is provided with the total uncertainty, including the MC statistical uncertainty, detector-related systematic uncertainties and theoretical uncertainties. The rightmost bin includes overflow events.

The distributions of $m_{\mathrm{b1},\mathrm{\cancel{b1}}}$ in the leptonic top inclusive signal region. For each bin yields for the data and total SM prediction are provided. The SM prediction is provided with the total uncertainty, including the MC statistical uncertainty, detector-related systematic uncertainties and theoretical uncertainties. The rightmost bin includes overflow events.

Cutflow for the reference point $(\it{m}_{\mathrm{H^{\pm}}}, \it{m}_{a}, tan\beta, sin\theta )=$ (500,100,1,0.7) , (800,150,20,0.7), (600,250,30,0.7), (1000,400,1,0.7) in 0L regions. Results are shown including all correction factors applied to simulation, and is normalised to 139 fb$^{-1}$.

Cutflow for the reference point $(\it{m}_{\mathrm{H^{\pm}}}, \it{m}_{a}, tan\beta, sin\theta )=$ (500,100,1,0.7) , (800,150,20,0.7), (600,250,30,0.7), (1000,400,1,0.7) in 1L leptonic top regions. Results are shown including all correction factors applied to simulation, and is normalised to 139 fb$^{-1}$.

Cutflow for the reference point $(\it{m}_{\mathrm{H^{\pm}}}, \it{m}_{a}, tan\beta, sin\theta )=$ (500,100,1,0.7) , (800,150,20,0.7), (600,250,30,0.7), (1000,400,1,0.7) in 1L hadronic top regions. Results are shown including all correction factors applied to simulation, and is normalised to 139 fb$^{-1}$.

Signal acceptance in the 0L region for 2HDM+a model DM signals on the plane defined by m$_a$--m$_{H^{\pm}}$ assuming tan$\beta$ = 1, m$_{\chi}$= 10 GeV and sin$\theta$ = 0.7. Please mind that the acceptance given in the table is multiplied by factor of $10^{3}$

Signal efficiency in the 0L region for 2HDM+a model DM signals on the plane defined by m$_a$--m$_{H^{\pm}}$ assuming tan$\beta$ = 1, m$_{\chi}$= 10 GeV and sin$\theta$ = 0.7. Please mind that the efficiency given in the table is multiplied by factor of $10^{2}$

Signal acceptance in the 0L region for 2HDM+a model DM signals on the plane defined by m$_{H^{\pm}}$--tan$\beta$ assuming m$_a$ = 150 GeV, m$_{\chi}$= 10 GeV and sin$\theta$ = 0.7. Please mind that the acceptance given in the table is multiplied by factor of $10^{3}$

Signal efficiency in the 0L region for 2HDM+a model DM signals on the plane defined by m$_{H^{\pm}}$--tan$\beta$ assuming m$_a$ = 150 GeV, m$_{\chi}$= 10 GeV and sin$\theta$ = 0.7. Please mind that the efficiency given in the table is multiplied by factor of $10^{2}$

Signal acceptance in the 0L region for 2HDM+a model DM signals on the plane defined by m$_{H^{\pm}}$--tan$\beta$ assuming m$_a$ = 250 GeV, m$_{\chi}$= 10 GeV and sin$\theta$ = 0.7. Please mind that the acceptance given in the table is multiplied by factor of $10^{3}$

Signal efficiency in the 0L region for 2HDM+a model DM signals on the plane defined by m$_{H^{\pm}}$--tan$\beta$ assuming m$_a$ = 250 GeV, m$_{\chi}$= 10 GeV and sin$\theta$ = 0.7. Please mind that the efficiency given in the table is multiplied by factor of $10^{2}$

Signal acceptance in the 1L region for 2HDM+a model DM signals on the plane defined by m$_a$--m$_{H^{\pm}}$ assuming tan$\beta$ = 1, m$_{\chi}$= 10 GeV and sin$\theta$ = 0.7. Please mind that the acceptance given in the table is multiplied by factor of $10^{3}$

Signal efficiency in the 1L region for 2HDM+a model DM signals on the plane defined by m$_a$--m$_{H^{\pm}}$ assuming tan$\beta$ = 1, m$_{\chi}$= 10 GeV and sin$\theta$ = 0.7. Please mind that the efficiency given in the table is multiplied by factor of $10^{2}$

Signal acceptance in the 1L region for 2HDM+a model DM signals on the plane defined by m$_{H^{\pm}}$--tan$\beta$ assuming m$_a$ = 150 GeV, m$_{\chi}$= 10 GeV and sin$\theta$ = 0.7. Please mind that the acceptance given in the table is multiplied by factor of $10^{3}$

Signal efficiency in the 1L region for 2HDM+a model DM signals on the plane defined by m$_{H^{\pm}}$--tan$\beta$ assuming m$_a$ = 150 GeV, m$_{\chi}$= 10 GeV and sin$\theta$ = 0.7. Please mind that the efficiency given in the table is multiplied by factor of $10^{2}$

Signal acceptance in the 1L region for 2HDM+a model DM signals on the plane defined by m$_{H^{\pm}}$--tan$\beta$ assuming m$_a$ = 250 GeV, m$_{\chi}$= 10 GeV and sin$\theta$ = 0.7. Please mind that the acceptance given in the table is multiplied by factor of $10^{3}$

Signal efficiency in the 1L region for 2HDM+a model DM signals on the plane defined by m$_{H^{\pm}}$--tan$\beta$ assuming m$_a$ = 250 GeV, m$_{\chi}$= 10 GeV and sin$\theta$ = 0.7. Please mind that the efficiency given in the table is multiplied by factor of $10^{2}$

Signal acceptance in the 0L region for 2HDM+a model DM signals on the plane defined by m$_a$--m$_{H^{\pm}}$ assuming tan$\beta$ = 1, m$_{\chi}$= 10 GeV and sin$\theta$ = 0.35. Please mind that the acceptance given in the table is multiplied by factor of $10^{3}$

Signal efficiency in the 0L region for 2HDM+a model DM signals on the plane defined by m$_a$--m$_{H^{\pm}}$ assuming tan$\beta$ = 1, m$_{\chi}$= 10 GeV and sin$\theta$ = 0.35. Please mind that the efficiency given in the table is multiplied by factor of $10^{2}$

Signal acceptance in the 0L region for 2HDM+a model DM signals on the plane defined by m$_{H^{\pm}}$--tan$\beta$ assuming m$_a$ = 150 GeV, m$_{\chi}$= 10 GeV and sin$\theta$ = 0.35. Please mind that the acceptance given in the table is multiplied by factor of $10^{3}$

Signal efficiency in the 0L region for 2HDM+a model DM signals on the plane defined by m$_{H^{\pm}}$--tan$\beta$ assuming m$_a$ = 150 GeV, m$_{\chi}$= 10 GeV and sin$\theta$ = 0.35. Please mind that the efficiency given in the table is multiplied by factor of $10^{2}$

Signal acceptance in the 0L region for 2HDM+a model DM signals on the plane defined by m$_{H^{\pm}}$--tan$\beta$ assuming m$_a$ = 250 GeV, m$_{\chi}$= 10 GeV and sin$\theta$ = 0.35. Please mind that the acceptance given in the table is multiplied by factor of $10^{3}$

Signal efficiency in the 0L region for 2HDM+a model DM signals on the plane defined by m$_{H^{\pm}}$--tan$\beta$ assuming m$_a$ = 250 GeV, m$_{\chi}$= 10 GeV and sin$\theta$ = 0.35. Please mind that the efficiency given in the table is multiplied by factor of $10^{2}$

Signal acceptance in the 1L region for 2HDM+a model DM signals on the plane defined by m$_a$--m$_{H^{\pm}}$ assuming tan$\beta$ = 1, m$_{\chi}$= 10 GeV and sin$\theta$ = 0.35. Please mind that the acceptance given in the table is multiplied by factor of $10^{3}$

Signal efficiency in the 1L region for 2HDM+a model DM signals on the plane defined by m$_a$--m$_{H^{\pm}}$ assuming tan$\beta$ = 1, m$_{\chi}$= 10 GeV and sin$\theta$ = 0.35. Please mind that the efficiency given in the table is multiplied by factor of $10^{2}$

Signal acceptance in the 1L region for 2HDM+a model DM signals on the plane defined by m$_{H^{\pm}}$--tan$\beta$ assuming m$_a$ = 150 GeV, m$_{\chi}$= 10 GeV and sin$\theta$ = 0.35. Please mind that the acceptance given in the table is multiplied by factor of $10^{3}$

Signal efficiency in the 1L region for 2HDM+a model DM signals on the plane defined by m$_{H^{\pm}}$--tan$\beta$ assuming m$_a$ = 150 GeV, m$_{\chi}$= 10 GeV and sin$\theta$ = 0.35. Please mind that the efficiency given in the table is multiplied by factor of $10^{2}$

Signal acceptance in the 1L region for 2HDM+a model DM signals on the plane defined by m$_{H^{\pm}}$--tan$\beta$ assuming m$_a$ = 250 GeV, m$_{\chi}$= 10 GeV and sin$\theta$ = 0.35. Please mind that the acceptance given in the table is multiplied by factor of $10^{3}$

Signal efficiency in the 1L region for 2HDM+a model DM signals on the plane defined by m$_{H^{\pm}}$--tan$\beta$ assuming m$_a$ = 250 GeV, m$_{\chi}$= 10 GeV and sin$\theta$ = 0.35. Please mind that the efficiency given in the table is multiplied by factor of $10^{2}$

$\Delta R_{\rm axis}$ distribution for WTA$\textendash$SD with grooming setting ($z_{\rm cut}=0.1,\beta=2$) for jets of $R=0.2$, in the interval $20<p_{\rm T}^{\rm ch \ jet}<40 \ {\rm GeV}/c$.

$\Delta R_{\rm axis}$ distribution for WTA$\textendash$SD with grooming setting ($z_{\rm cut}=0.1,\beta=3$) for jets of $R=0.2$, in the interval $20<p_{\rm T}^{\rm ch \ jet}<40 \ {\rm GeV}/c$.

$\Delta R_{\rm axis}$ distribution for WTA$\textendash$SD with grooming setting ($z_{\rm cut}=0.2,\beta=1$) for jets of $R=0.2$, in the interval $20<p_{\rm T}^{\rm ch \ jet}<40 \ {\rm GeV}/c$.

$\Delta R_{\rm axis}$ distribution for WTA$\textendash$SD with grooming setting ($z_{\rm cut}=0.3,\beta=1$) for jets of $R=0.2$, in the interval $20<p_{\rm T}^{\rm ch \ jet}<40 \ {\rm GeV}/c$.

$\Delta R_{\rm axis}$ distribution for Standard$\textendash$SD with grooming setting ($z_{\rm cut}=0.3,\beta=1$) for jets of $R=0.2$, in the interval $20<p_{\rm T}^{\rm ch \ jet}<40 \ {\rm GeV}/c$.

$\Delta R_{\rm axis}$ distribution for Standard$\textendash$SD with grooming setting ($z_{\rm cut}=0.2,\beta=1$) for jets of $R=0.2$, in the interval $20<p_{\rm T}^{\rm ch \ jet}<40 \ {\rm GeV}/c$.

$\Delta R_{\rm axis}$ distribution for Standard$\textendash$SD with grooming setting ($z_{\rm cut}=0.1,\beta=0$) for jets of $R=0.2$, in the interval $20<p_{\rm T}^{\rm ch \ jet}<40 \ {\rm GeV}/c$.

$\Delta R_{\rm axis}$ distribution for Standard$\textendash$SD with grooming setting ($z_{\rm cut}=0.1,\beta=1$) for jets of $R=0.2$, in the interval $20<p_{\rm T}^{\rm ch \ jet}<40 \ {\rm GeV}/c$.

$\Delta R_{\rm axis}$ distribution for Standard$\textendash$SD with grooming setting ($z_{\rm cut}=0.1,\beta=2$) for jets of $R=0.2$, in the interval $20<p_{\rm T}^{\rm ch \ jet}<40 \ {\rm GeV}/c$.

$\Delta R_{\rm axis}$ distribution for Standard$\textendash$SD with grooming setting ($z_{\rm cut}=0.1,\beta=3$) for jets of $R=0.2$, in the interval $20<p_{\rm T}^{\rm ch \ jet}<40 \ {\rm GeV}/c$.

$\Delta R_{\rm axis}$ distribution for WTA$\textendash$Standard for jets of $R=0.2$, in the interval $40<p_{\rm T}^{\rm ch \ jet}<60 \ {\rm GeV}/c$.

$\Delta R_{\rm axis}$ distribution for WTA$\textendash$SD with grooming setting ($z_{\rm cut}=0.1,\beta=0$) for jets of $R=0.2$, in the interval $40<p_{\rm T}^{\rm ch \ jet}<60 \ {\rm GeV}/c$.

$\Delta R_{\rm axis}$ distribution for WTA$\textendash$SD with grooming setting ($z_{\rm cut}=0.1,\beta=1$) for jets of $R=0.2$, in the interval $40<p_{\rm T}^{\rm ch \ jet}<60 \ {\rm GeV}/c$.

$\Delta R_{\rm axis}$ distribution for WTA$\textendash$SD with grooming setting ($z_{\rm cut}=0.1,\beta=2$) for jets of $R=0.2$, in the interval $40<p_{\rm T}^{\rm ch \ jet}<60 \ {\rm GeV}/c$.

$\Delta R_{\rm axis}$ distribution for WTA$\textendash$SD with grooming setting ($z_{\rm cut}=0.1,\beta=3$) for jets of $R=0.2$, in the interval $40<p_{\rm T}^{\rm ch \ jet}<60 \ {\rm GeV}/c$.

$\Delta R_{\rm axis}$ distribution for WTA$\textendash$SD with grooming setting ($z_{\rm cut}=0.2,\beta=1$) for jets of $R=0.2$, in the interval $40<p_{\rm T}^{\rm ch \ jet}<60 \ {\rm GeV}/c$.

$\Delta R_{\rm axis}$ distribution for WTA$\textendash$SD with grooming setting ($z_{\rm cut}=0.3,\beta=1$) for jets of $R=0.2$, in the interval $40<p_{\rm T}^{\rm ch \ jet}<60 \ {\rm GeV}/c$.

$\Delta R_{\rm axis}$ distribution for Standard$\textendash$SD with grooming setting ($z_{\rm cut}=0.3,\beta=1$) for jets of $R=0.2$, in the interval $40<p_{\rm T}^{\rm ch \ jet}<60 \ {\rm GeV}/c$.

$\Delta R_{\rm axis}$ distribution for Standard$\textendash$SD with grooming setting ($z_{\rm cut}=0.2,\beta=1$) for jets of $R=0.2$, in the interval $40<p_{\rm T}^{\rm ch \ jet}<60 \ {\rm GeV}/c$.

$\Delta R_{\rm axis}$ distribution for Standard$\textendash$SD with grooming setting ($z_{\rm cut}=0.1,\beta=0$) for jets of $R=0.2$, in the interval $40<p_{\rm T}^{\rm ch \ jet}<60 \ {\rm GeV}/c$.

$\Delta R_{\rm axis}$ distribution for Standard$\textendash$SD with grooming setting ($z_{\rm cut}=0.1,\beta=1$) for jets of $R=0.2$, in the interval $40<p_{\rm T}^{\rm ch \ jet}<60 \ {\rm GeV}/c$.

$\Delta R_{\rm axis}$ distribution for Standard$\textendash$SD with grooming setting ($z_{\rm cut}=0.1,\beta=2$) for jets of $R=0.2$, in the interval $40<p_{\rm T}^{\rm ch \ jet}<60 \ {\rm GeV}/c$.

$\Delta R_{\rm axis}$ distribution for Standard$\textendash$SD with grooming setting ($z_{\rm cut}=0.1,\beta=3$) for jets of $R=0.2$, in the interval $40<p_{\rm T}^{\rm ch \ jet}<60 \ {\rm GeV}/c$.

$\Delta R_{\rm axis}$ distribution for WTA$\textendash$Standard for jets of $R=0.2$, in the interval $60<p_{\rm T}^{\rm ch \ jet}<80 \ {\rm GeV}/c$.

$\Delta R_{\rm axis}$ distribution for WTA$\textendash$SD with grooming setting ($z_{\rm cut}=0.1,\beta=0$) for jets of $R=0.2$, in the interval $60<p_{\rm T}^{\rm ch \ jet}<80 \ {\rm GeV}/c$.

$\Delta R_{\rm axis}$ distribution for WTA$\textendash$SD with grooming setting ($z_{\rm cut}=0.1,\beta=1$) for jets of $R=0.2$, in the interval $60<p_{\rm T}^{\rm ch \ jet}<80 \ {\rm GeV}/c$.

$\Delta R_{\rm axis}$ distribution for WTA$\textendash$SD with grooming setting ($z_{\rm cut}=0.1,\beta=2$) for jets of $R=0.2$, in the interval $60<p_{\rm T}^{\rm ch \ jet}<80 \ {\rm GeV}/c$.

$\Delta R_{\rm axis}$ distribution for WTA$\textendash$SD with grooming setting ($z_{\rm cut}=0.1,\beta=3$) for jets of $R=0.2$, in the interval $60<p_{\rm T}^{\rm ch \ jet}<80 \ {\rm GeV}/c$.

$\Delta R_{\rm axis}$ distribution for WTA$\textendash$SD with grooming setting ($z_{\rm cut}=0.2,\beta=1$) for jets of $R=0.2$, in the interval $60<p_{\rm T}^{\rm ch \ jet}<80 \ {\rm GeV}/c$.

$\Delta R_{\rm axis}$ distribution for WTA$\textendash$SD with grooming setting ($z_{\rm cut}=0.3,\beta=1$) for jets of $R=0.2$, in the interval $60<p_{\rm T}^{\rm ch \ jet}<80 \ {\rm GeV}/c$.

$\Delta R_{\rm axis}$ distribution for Standard$\textendash$SD with grooming setting ($z_{\rm cut}=0.3,\beta=1$) for jets of $R=0.2$, in the interval $60<p_{\rm T}^{\rm ch \ jet}<80 \ {\rm GeV}/c$.

$\Delta R_{\rm axis}$ distribution for Standard$\textendash$SD with grooming setting ($z_{\rm cut}=0.2,\beta=1$) for jets of $R=0.2$, in the interval $60<p_{\rm T}^{\rm ch \ jet}<80 \ {\rm GeV}/c$.

$\Delta R_{\rm axis}$ distribution for Standard$\textendash$SD with grooming setting ($z_{\rm cut}=0.1,\beta=0$) for jets of $R=0.2$, in the interval $60<p_{\rm T}^{\rm ch \ jet}<80 \ {\rm GeV}/c$.

$\Delta R_{\rm axis}$ distribution for Standard$\textendash$SD with grooming setting ($z_{\rm cut}=0.1,\beta=1$) for jets of $R=0.2$, in the interval $60<p_{\rm T}^{\rm ch \ jet}<80 \ {\rm GeV}/c$.

$\Delta R_{\rm axis}$ distribution for Standard$\textendash$SD with grooming setting ($z_{\rm cut}=0.1,\beta=2$) for jets of $R=0.2$, in the interval $60<p_{\rm T}^{\rm ch \ jet}<80 \ {\rm GeV}/c$.

$\Delta R_{\rm axis}$ distribution for Standard$\textendash$SD with grooming setting ($z_{\rm cut}=0.1,\beta=3$) for jets of $R=0.2$, in the interval $60<p_{\rm T}^{\rm ch \ jet}<80 \ {\rm GeV}/c$.

$\Delta R_{\rm axis}$ distribution for WTA$\textendash$Standard for jets of $R=0.2$, in the interval $80<p_{\rm T}^{\rm ch \ jet}<100 \ {\rm GeV}/c$.

$\Delta R_{\rm axis}$ distribution for WTA$\textendash$SD with grooming setting ($z_{\rm cut}=0.1,\beta=0$) for jets of $R=0.2$, in the interval $80<p_{\rm T}^{\rm ch \ jet}<100 \ {\rm GeV}/c$.

$\Delta R_{\rm axis}$ distribution for WTA$\textendash$SD with grooming setting ($z_{\rm cut}=0.1,\beta=1$) for jets of $R=0.2$, in the interval $80<p_{\rm T}^{\rm ch \ jet}<100 \ {\rm GeV}/c$.

$\Delta R_{\rm axis}$ distribution for WTA$\textendash$SD with grooming setting ($z_{\rm cut}=0.1,\beta=2$) for jets of $R=0.2$, in the interval $80<p_{\rm T}^{\rm ch \ jet}<100 \ {\rm GeV}/c$.

$\Delta R_{\rm axis}$ distribution for WTA$\textendash$SD with grooming setting ($z_{\rm cut}=0.1,\beta=3$) for jets of $R=0.2$, in the interval $80<p_{\rm T}^{\rm ch \ jet}<100 \ {\rm GeV}/c$.

$\Delta R_{\rm axis}$ distribution for WTA$\textendash$SD with grooming setting ($z_{\rm cut}=0.2,\beta=1$) for jets of $R=0.2$, in the interval $80<p_{\rm T}^{\rm ch \ jet}<100 \ {\rm GeV}/c$.

$\Delta R_{\rm axis}$ distribution for WTA$\textendash$SD with grooming setting ($z_{\rm cut}=0.3,\beta=1$) for jets of $R=0.2$, in the interval $80<p_{\rm T}^{\rm ch \ jet}<100 \ {\rm GeV}/c$.

$\Delta R_{\rm axis}$ distribution for Standard$\textendash$SD with grooming setting ($z_{\rm cut}=0.3,\beta=1$) for jets of $R=0.2$, in the interval $80<p_{\rm T}^{\rm ch \ jet}<100 \ {\rm GeV}/c$.

$\Delta R_{\rm axis}$ distribution for Standard$\textendash$SD with grooming setting ($z_{\rm cut}=0.2,\beta=1$) for jets of $R=0.2$, in the interval $80<p_{\rm T}^{\rm ch \ jet}<100 \ {\rm GeV}/c$.

$\Delta R_{\rm axis}$ distribution for Standard$\textendash$SD with grooming setting ($z_{\rm cut}=0.1,\beta=0$) for jets of $R=0.2$, in the interval $80<p_{\rm T}^{\rm ch \ jet}<100 \ {\rm GeV}/c$.

$\Delta R_{\rm axis}$ distribution for Standard$\textendash$SD with grooming setting ($z_{\rm cut}=0.1,\beta=1$) for jets of $R=0.2$, in the interval $80<p_{\rm T}^{\rm ch \ jet}<100 \ {\rm GeV}/c$.

$\Delta R_{\rm axis}$ distribution for Standard$\textendash$SD with grooming setting ($z_{\rm cut}=0.1,\beta=2$) for jets of $R=0.2$, in the interval $80<p_{\rm T}^{\rm ch \ jet}<100 \ {\rm GeV}/c$.

$\Delta R_{\rm axis}$ distribution for Standard$\textendash$SD with grooming setting ($z_{\rm cut}=0.1,\beta=3$) for jets of $R=0.2$, in the interval $80<p_{\rm T}^{\rm ch \ jet}<100 \ {\rm GeV}/c$.

$\Delta R_{\rm axis}$ distribution for WTA$\textendash$Standard for jets of $R=0.4$, in the interval $20<p_{\rm T}^{\rm ch \ jet}<40 \ {\rm GeV}/c$.

$\Delta R_{\rm axis}$ distribution for WTA$\textendash$SD with grooming setting ($z_{\rm cut}=0.1,\beta=0$) for jets of $R=0.4$, in the interval $20<p_{\rm T}^{\rm ch \ jet}<40 \ {\rm GeV}/c$.

$\Delta R_{\rm axis}$ distribution for WTA$\textendash$SD with grooming setting ($z_{\rm cut}=0.1,\beta=1$) for jets of $R=0.4$, in the interval $20<p_{\rm T}^{\rm ch \ jet}<40 \ {\rm GeV}/c$.

$\Delta R_{\rm axis}$ distribution for WTA$\textendash$SD with grooming setting ($z_{\rm cut}=0.1,\beta=2$) for jets of $R=0.4$, in the interval $20<p_{\rm T}^{\rm ch \ jet}<40 \ {\rm GeV}/c$.

$\Delta R_{\rm axis}$ distribution for WTA$\textendash$SD with grooming setting ($z_{\rm cut}=0.1,\beta=3$) for jets of $R=0.4$, in the interval $20<p_{\rm T}^{\rm ch \ jet}<40 \ {\rm GeV}/c$.

$\Delta R_{\rm axis}$ distribution for WTA$\textendash$SD with grooming setting ($z_{\rm cut}=0.2,\beta=1$) for jets of $R=0.4$, in the interval $20<p_{\rm T}^{\rm ch \ jet}<40 \ {\rm GeV}/c$.

$\Delta R_{\rm axis}$ distribution for WTA$\textendash$SD with grooming setting ($z_{\rm cut}=0.3,\beta=1$) for jets of $R=0.4$, in the interval $20<p_{\rm T}^{\rm ch \ jet}<40 \ {\rm GeV}/c$.

$\Delta R_{\rm axis}$ distribution for Standard$\textendash$SD with grooming setting ($z_{\rm cut}=0.3,\beta=1$) for jets of $R=0.4$, in the interval $20<p_{\rm T}^{\rm ch \ jet}<40 \ {\rm GeV}/c$.

$\Delta R_{\rm axis}$ distribution for Standard$\textendash$SD with grooming setting ($z_{\rm cut}=0.2,\beta=1$) for jets of $R=0.4$, in the interval $20<p_{\rm T}^{\rm ch \ jet}<40 \ {\rm GeV}/c$.

$\Delta R_{\rm axis}$ distribution for Standard$\textendash$SD with grooming setting ($z_{\rm cut}=0.1,\beta=0$) for jets of $R=0.4$, in the interval $20<p_{\rm T}^{\rm ch \ jet}<40 \ {\rm GeV}/c$.

$\Delta R_{\rm axis}$ distribution for Standard$\textendash$SD with grooming setting ($z_{\rm cut}=0.1,\beta=1$) for jets of $R=0.4$, in the interval $20<p_{\rm T}^{\rm ch \ jet}<40 \ {\rm GeV}/c$.

$\Delta R_{\rm axis}$ distribution for Standard$\textendash$SD with grooming setting ($z_{\rm cut}=0.1,\beta=2$) for jets of $R=0.4$, in the interval $20<p_{\rm T}^{\rm ch \ jet}<40 \ {\rm GeV}/c$.

$\Delta R_{\rm axis}$ distribution for Standard$\textendash$SD with grooming setting ($z_{\rm cut}=0.1,\beta=3$) for jets of $R=0.4$, in the interval $20<p_{\rm T}^{\rm ch \ jet}<40 \ {\rm GeV}/c$.

$\Delta R_{\rm axis}$ distribution for WTA$\textendash$Standard for jets of $R=0.4$, in the interval $40<p_{\rm T}^{\rm ch \ jet}<60 \ {\rm GeV}/c$.

$\Delta R_{\rm axis}$ distribution for WTA$\textendash$SD with grooming setting ($z_{\rm cut}=0.1,\beta=0$) for jets of $R=0.4$, in the interval $40<p_{\rm T}^{\rm ch \ jet}<60 \ {\rm GeV}/c$.

$\Delta R_{\rm axis}$ distribution for WTA$\textendash$SD with grooming setting ($z_{\rm cut}=0.1,\beta=1$) for jets of $R=0.4$, in the interval $40<p_{\rm T}^{\rm ch \ jet}<60 \ {\rm GeV}/c$.

$\Delta R_{\rm axis}$ distribution for WTA$\textendash$SD with grooming setting ($z_{\rm cut}=0.1,\beta=2$) for jets of $R=0.4$, in the interval $40<p_{\rm T}^{\rm ch \ jet}<60 \ {\rm GeV}/c$.

$\Delta R_{\rm axis}$ distribution for WTA$\textendash$SD with grooming setting ($z_{\rm cut}=0.1,\beta=3$) for jets of $R=0.4$, in the interval $40<p_{\rm T}^{\rm ch \ jet}<60 \ {\rm GeV}/c$.

$\Delta R_{\rm axis}$ distribution for WTA$\textendash$SD with grooming setting ($z_{\rm cut}=0.2,\beta=1$) for jets of $R=0.4$, in the interval $40<p_{\rm T}^{\rm ch \ jet}<60 \ {\rm GeV}/c$.

$\Delta R_{\rm axis}$ distribution for WTA$\textendash$SD with grooming setting ($z_{\rm cut}=0.3,\beta=1$) for jets of $R=0.4$, in the interval $40<p_{\rm T}^{\rm ch \ jet}<60 \ {\rm GeV}/c$.

$\Delta R_{\rm axis}$ distribution for Standard$\textendash$SD with grooming setting ($z_{\rm cut}=0.3,\beta=1$) for jets of $R=0.4$, in the interval $40<p_{\rm T}^{\rm ch \ jet}<60 \ {\rm GeV}/c$.

$\Delta R_{\rm axis}$ distribution for Standard$\textendash$SD with grooming setting ($z_{\rm cut}=0.2,\beta=1$) for jets of $R=0.4$, in the interval $40<p_{\rm T}^{\rm ch \ jet}<60 \ {\rm GeV}/c$.

$\Delta R_{\rm axis}$ distribution for Standard$\textendash$SD with grooming setting ($z_{\rm cut}=0.1,\beta=0$) for jets of $R=0.4$, in the interval $40<p_{\rm T}^{\rm ch \ jet}<60 \ {\rm GeV}/c$.

$\Delta R_{\rm axis}$ distribution for Standard$\textendash$SD with grooming setting ($z_{\rm cut}=0.1,\beta=1$) for jets of $R=0.4$, in the interval $40<p_{\rm T}^{\rm ch \ jet}<60 \ {\rm GeV}/c$.

$\Delta R_{\rm axis}$ distribution for Standard$\textendash$SD with grooming setting ($z_{\rm cut}=0.1,\beta=2$) for jets of $R=0.4$, in the interval $40<p_{\rm T}^{\rm ch \ jet}<60 \ {\rm GeV}/c$.

$\Delta R_{\rm axis}$ distribution for Standard$\textendash$SD with grooming setting ($z_{\rm cut}=0.1,\beta=3$) for jets of $R=0.4$, in the interval $40<p_{\rm T}^{\rm ch \ jet}<60 \ {\rm GeV}/c$.

$\Delta R_{\rm axis}$ distribution for WTA$\textendash$Standard for jets of $R=0.4$, in the interval $60<p_{\rm T}^{\rm ch \ jet}<80 \ {\rm GeV}/c$.

$\Delta R_{\rm axis}$ distribution for WTA$\textendash$SD with grooming setting ($z_{\rm cut}=0.1,\beta=0$) for jets of $R=0.4$, in the interval $60<p_{\rm T}^{\rm ch \ jet}<80 \ {\rm GeV}/c$.

$\Delta R_{\rm axis}$ distribution for WTA$\textendash$SD with grooming setting ($z_{\rm cut}=0.1,\beta=1$) for jets of $R=0.4$, in the interval $60<p_{\rm T}^{\rm ch \ jet}<80 \ {\rm GeV}/c$.

$\Delta R_{\rm axis}$ distribution for WTA$\textendash$SD with grooming setting ($z_{\rm cut}=0.1,\beta=2$) for jets of $R=0.4$, in the interval $60<p_{\rm T}^{\rm ch \ jet}<80 \ {\rm GeV}/c$.

$\Delta R_{\rm axis}$ distribution for WTA$\textendash$SD with grooming setting ($z_{\rm cut}=0.1,\beta=3$) for jets of $R=0.4$, in the interval $60<p_{\rm T}^{\rm ch \ jet}<80 \ {\rm GeV}/c$.

$\Delta R_{\rm axis}$ distribution for WTA$\textendash$SD with grooming setting ($z_{\rm cut}=0.2,\beta=1$) for jets of $R=0.4$, in the interval $60<p_{\rm T}^{\rm ch \ jet}<80 \ {\rm GeV}/c$.

$\Delta R_{\rm axis}$ distribution for WTA$\textendash$SD with grooming setting ($z_{\rm cut}=0.3,\beta=1$) for jets of $R=0.4$, in the interval $60<p_{\rm T}^{\rm ch \ jet}<80 \ {\rm GeV}/c$.

$\Delta R_{\rm axis}$ distribution for Standard$\textendash$SD with grooming setting ($z_{\rm cut}=0.3,\beta=1$) for jets of $R=0.4$, in the interval $60<p_{\rm T}^{\rm ch \ jet}<80 \ {\rm GeV}/c$.

$\Delta R_{\rm axis}$ distribution for Standard$\textendash$SD with grooming setting ($z_{\rm cut}=0.2,\beta=1$) for jets of $R=0.4$, in the interval $60<p_{\rm T}^{\rm ch \ jet}<80 \ {\rm GeV}/c$.

$\Delta R_{\rm axis}$ distribution for Standard$\textendash$SD with grooming setting ($z_{\rm cut}=0.1,\beta=0$) for jets of $R=0.4$, in the interval $60<p_{\rm T}^{\rm ch \ jet}<80 \ {\rm GeV}/c$.

$\Delta R_{\rm axis}$ distribution for Standard$\textendash$SD with grooming setting ($z_{\rm cut}=0.1,\beta=1$) for jets of $R=0.4$, in the interval $60<p_{\rm T}^{\rm ch \ jet}<80 \ {\rm GeV}/c$.

$\Delta R_{\rm axis}$ distribution for Standard$\textendash$SD with grooming setting ($z_{\rm cut}=0.1,\beta=2$) for jets of $R=0.4$, in the interval $60<p_{\rm T}^{\rm ch \ jet}<80 \ {\rm GeV}/c$.

$\Delta R_{\rm axis}$ distribution for Standard$\textendash$SD with grooming setting ($z_{\rm cut}=0.1,\beta=3$) for jets of $R=0.4$, in the interval $60<p_{\rm T}^{\rm ch \ jet}<80 \ {\rm GeV}/c$.

$\Delta R_{\rm axis}$ distribution for WTA$\textendash$Standard for jets of $R=0.4$, in the interval $80<p_{\rm T}^{\rm ch \ jet}<100 \ {\rm GeV}/c$.

$\Delta R_{\rm axis}$ distribution for WTA$\textendash$SD with grooming setting ($z_{\rm cut}=0.1,\beta=0$) for jets of $R=0.4$, in the interval $80<p_{\rm T}^{\rm ch \ jet}<100 \ {\rm GeV}/c$.

$\Delta R_{\rm axis}$ distribution for WTA$\textendash$SD with grooming setting ($z_{\rm cut}=0.1,\beta=1$) for jets of $R=0.4$, in the interval $80<p_{\rm T}^{\rm ch \ jet}<100 \ {\rm GeV}/c$.

$\Delta R_{\rm axis}$ distribution for WTA$\textendash$SD with grooming setting ($z_{\rm cut}=0.1,\beta=2$) for jets of $R=0.4$, in the interval $80<p_{\rm T}^{\rm ch \ jet}<100 \ {\rm GeV}/c$.

$\Delta R_{\rm axis}$ distribution for WTA$\textendash$SD with grooming setting ($z_{\rm cut}=0.1,\beta=3$) for jets of $R=0.4$, in the interval $80<p_{\rm T}^{\rm ch \ jet}<100 \ {\rm GeV}/c$.

$\Delta R_{\rm axis}$ distribution for WTA$\textendash$SD with grooming setting ($z_{\rm cut}=0.2,\beta=1$) for jets of $R=0.4$, in the interval $80<p_{\rm T}^{\rm ch \ jet}<100 \ {\rm GeV}/c$.

$\Delta R_{\rm axis}$ distribution for WTA$\textendash$SD with grooming setting ($z_{\rm cut}=0.3,\beta=1$) for jets of $R=0.4$, in the interval $80<p_{\rm T}^{\rm ch \ jet}<100 \ {\rm GeV}/c$.

$\Delta R_{\rm axis}$ distribution for Standard$\textendash$SD with grooming setting ($z_{\rm cut}=0.3,\beta=1$) for jets of $R=0.4$, in the interval $80<p_{\rm T}^{\rm ch \ jet}<100 \ {\rm GeV}/c$.

$\Delta R_{\rm axis}$ distribution for Standard$\textendash$SD with grooming setting ($z_{\rm cut}=0.2,\beta=1$) for jets of $R=0.4$, in the interval $80<p_{\rm T}^{\rm ch \ jet}<100 \ {\rm GeV}/c$.

$\Delta R_{\rm axis}$ distribution for Standard$\textendash$SD with grooming setting ($z_{\rm cut}=0.1,\beta=0$) for jets of $R=0.4$, in the interval $80<p_{\rm T}^{\rm ch \ jet}<100 \ {\rm GeV}/c$.

$\Delta R_{\rm axis}$ distribution for Standard$\textendash$SD with grooming setting ($z_{\rm cut}=0.1,\beta=1$) for jets of $R=0.4$, in the interval $80<p_{\rm T}^{\rm ch \ jet}<100 \ {\rm GeV}/c$.

$\Delta R_{\rm axis}$ distribution for Standard$\textendash$SD with grooming setting ($z_{\rm cut}=0.1,\beta=2$) for jets of $R=0.4$, in the interval $80<p_{\rm T}^{\rm ch \ jet}<100 \ {\rm GeV}/c$.

$\Delta R_{\rm axis}$ distribution for Standard$\textendash$SD with grooming setting ($z_{\rm cut}=0.1,\beta=3$) for jets of $R=0.4$, in the interval $80<p_{\rm T}^{\rm ch \ jet}<100 \ {\rm GeV}/c$.

Distribution of the most-probable value of the WTA$\textendash$Standard spectra as a function of $p_{\rm T}^{\rm ch \ jet}$ for $R=0.2$.

Distribution of the most-probable value of the WTA$\textendash$Standard spectra as a function of $p_{\rm T}^{\rm ch \ jet}$ for $R=0.4$.

$p_{\rm T}$-differential density of $\Xi^{\pm}$ in jets and the underlying event in pp collisions at $\sqrt{s} = 13$ TeV.

$p_{\rm T}$-differential density of inclusive $\Omega^{\pm}$ in pp collisions at $\sqrt{s} = 13$ TeV.

$p_{\rm T}$-differential density of $\Omega^{\pm}$ in jets and the underlying event in pp collisions at $\sqrt{s} = 13$ TeV.

$p_{\rm T}$-differential $(\Lambda + \overline{\Lambda})/2{\rm K}_{\rm S}^{0}$ ratio of inclusive particles in pp collisions at $\sqrt{s} = 13$ TeV.

$p_{\rm T}$-differential $(\Lambda + \overline{\Lambda})/2{\rm K}_{\rm S}^{0}$ ratio in jets and the underlying event in pp collisions at $\sqrt{s} = 13$ TeV.

$p_{\rm T}$-differential $(\Xi^{-} + \Xi^{+}) / 2{\rm K}_{\rm S}^{0}$ and $(\Xi^{-} + \Xi^{+}) / (\Lambda + \overline{\Lambda})$ ratios of inclusive particles in pp collisions at $\sqrt{s} = 13$ TeV.

$p_{\rm T}$-differential $(\Xi^{-} + \Xi^{+}) / 2{\rm K}_{\rm S}^{0}$ and $(\Xi^{-} + \Xi^{+}) / (\Lambda + \overline{\Lambda})$ ratios in jets and the underlying event in pp collisions at $\sqrt{s} = 13$ TeV.

$p_{\rm T}$-differential $(\Omega^{-} + \Omega^{+}) / 2{\rm K}_{\rm S}^{0}$, $(\Omega^{-} + \Omega^{+}) / (\Lambda + \overline{\Lambda})$ and $(\Omega^{-} + \Omega^{+}) / (\Xi^{-} + \Xi^{+})$ ratios of inclusive particles in pp collisions at $\sqrt{s} = 13$ TeV.

$p_{\rm T}$-differential $(\Omega^{-} + \Omega^{+}) / 2{\rm K}_{\rm S}^{0}$, $(\Omega^{-} + \Omega^{+}) / (\Lambda + \overline{\Lambda})$ and $(\Omega^{-} + \Omega^{+}) / (\Xi^{-} + \Xi^{+})$ ratios in jets and the underlying event in pp collisions at $\sqrt{s} = 13$ TeV.

$p_{\rm T}$-differential density of inclusive ${\rm K}_{\rm S}^{0}$ and $\Lambda$ ($\overline{\Lambda}$) in p-Pb collisions at $\sqrt{s_{\rm NN}} = 5.02$ TeV.

$p_{\rm T}$-differential densities of ${\rm K}_{\rm S}^{0}$ and $\Lambda$ ($\overline{\Lambda}$) in jets and the underlying event in p-Pb collisions at $\sqrt{s_{\rm NN}} = 5.02$ TeV.

$p_{\rm T}$-differential density of inclusive $\Xi^{\pm}$ in p-Pb collisions at $\sqrt{s_{\rm NN}} = 5.02$ TeV.

$p_{\rm T}$-differential density of $\Xi^{\pm}$ in jets and the underlying event in p-Pb collisions at $\sqrt{s_{\rm NN}} = 5.02$ TeV.

$p_{\rm T}$-differential density of inclusive $\Omega^{\pm}$ in p-Pb collisions at $\sqrt{s} = 5.02$ TeV.

$p_{\rm T}$-differential density of $\Omega^{\pm}$ in jets and the underlying event in p-Pb collisions at $\sqrt{s_{\rm NN}} = 5.02$ TeV.

$p_{\rm T}$-differential $(\Lambda + \overline{\Lambda})/2{\rm K}_{\rm S}^{0}$ ratio in jets for various event activity classes in p-Pb collisions at $\sqrt{s_{\rm NN}} = 5.02$ TeV.

$p_{\rm T}$-differential $(\Xi^{-} + \Xi^{+}) / 2{\rm K}_{\rm S}^{0}$ and $(\Xi^{-} + \Xi^{+}) / (\Lambda + \overline{\Lambda})$ ratios in jets for various event activity classes in p-Pb collisions at $\sqrt{s_{\rm NN}} = 5.02$ TeV.

$p_{\rm T}$-differential $(\Omega^{-} + \Omega^{+}) / 2{\rm K}_{\rm S}^{0}$, $(\Omega^{-} + \Omega^{+}) / (\Lambda + \overline{\Lambda})$ and $(\Omega^{-} + \Omega^{+}) / (\Xi^{-} + \Xi^{+})$ ratios in jets in MB p-Pb collisions at $\sqrt{s_{\rm NN}} = 5.02$ TeV.

$p_{\rm T}$-differential density of inclusive ${\rm K}_{\rm S}^{0}$ and $\Lambda$ ($\overline{\Lambda}$) for $0$-$10\%$ event activity class in p-Pb collisions at $\sqrt{s_{\rm NN}} = 5.02$ TeV.

$p_{\rm T}$-differential densities of ${\rm K}_{\rm S}^{0}$ and $\Lambda$ ($\overline{\Lambda}$) in jets and the underlying event for $0$-$10\%$ event activity class in p-Pb collisions at $\sqrt{s_{\rm NN}} = 5.02$ TeV.

$p_{\rm T}$-differential density of inclusive $\Xi^{\pm}$ for $0$-$10\%$ event activity class in p-Pb collisions at $\sqrt{s_{\rm NN}} = 5.02$ TeV.

$p_{\rm T}$-differential density of $\Xi^{\pm}$ in jets and the underlying event for $0$-$10\%$ event activity class in p-Pb collisions at $\sqrt{s_{\rm NN}} = 5.02$ TeV.

$p_{\rm T}$-differential density of inclusive ${\rm K}_{\rm S}^{0}$ and $\Lambda$ ($\overline{\Lambda}$) for $10$-$40\%$ event activity class in p-Pb collisions at $\sqrt{s_{\rm NN}} = 5.02$ TeV.

$p_{\rm T}$-differential densities of ${\rm K}_{\rm S}^{0}$ and $\Lambda$ ($\overline{\Lambda}$) in jets and the underlying event for $10$-$40\%$ event activity class in p-Pb collisions at $\sqrt{s_{\rm NN}} = 5.02$ TeV.

$p_{\rm T}$-differential density of inclusive $\Xi^{\pm}$ for $10$-$40\%$ event activity class in p-Pb collisions at $\sqrt{s_{\rm NN}} = 5.02$ TeV.

$p_{\rm T}$-differential density of $\Xi^{\pm}$ in jets and the underlying event for $10$-$40\%$ event activity class in p-Pb collisions at $\sqrt{s_{\rm NN}} = 5.02$ TeV.

$p_{\rm T}$-differential density of inclusive ${\rm K}_{\rm S}^{0}$ and $\Lambda$ ($\overline{\Lambda}$) for $40$-$100\%$ event activity class in p-Pb collisions at $\sqrt{s_{\rm NN}} = 5.02$ TeV.

$p_{\rm T}$-differential densities of ${\rm K}_{\rm S}^{0}$ and $\Lambda$ ($\overline{\Lambda}$) in jets and the underlying event for $40$-$100\%$ event activity class in p-Pb collisions at $\sqrt{s_{\rm NN}} = 5.02$ TeV.

$p_{\rm T}$-differential density of inclusive $\Xi^{\pm}$ for $40$-$100\%$ event activity class in p-Pb collisions at $\sqrt{s_{\rm NN}} = 5.02$ TeV.

$p_{\rm T}$-differential density of $\Xi^{\pm}$ in jets and the underlying event for $40$-$100\%$ event activity class in p-Pb collisions at $\sqrt{s_{\rm NN}} = 5.02$ TeV.

$p_{\rm T}$-differential $(\Lambda + \overline{\Lambda})/2{\rm K}_{\rm S}^{0}$ ratio of inclusive particles for various event activity classes in p-Pb collisions at $\sqrt{s_{\rm NN}} = 5.02$ TeV.

$p_{\rm T}$-differential $(\Lambda + \overline{\Lambda})/2{\rm K}_{\rm S}^{0}$ ratio in the underlying event for various event activity classes in p-Pb collisions at $\sqrt{s_{\rm NN}} = 5.02$ TeV.

$p_{\rm T}$-differential $(\Lambda + \overline{\Lambda})/2{\rm K}_{\rm S}^{0}$ ratio in jets for various event activity classes in p-Pb collisions at $\sqrt{s_{\rm NN}} = 5.02$ TeV.

$p_{\rm T}$-differential $(\Xi^{-} + \Xi^{+}) / 2{\rm K}_{\rm S}^{0}$ ratio of inclusive particles for various event activity classes in p-Pb collisions at $\sqrt{s_{\rm NN}} = 5.02$ TeV.

$p_{\rm T}$-differential $(\Xi^{-} + \Xi^{+}) / 2{\rm K}_{\rm S}^{0}$ ratio in the underlying event for $0$-$10\%$ event activity class in p-Pb collisions at $\sqrt{s_{\rm NN}} = 5.02$ TeV.

$p_{\rm T}$-differential $(\Xi^{-} + \Xi^{+}) / 2{\rm K}_{\rm S}^{0}$ ratio in the underlying event for $10$-$40\%$ and $40$-$100\%$ event activity classes in p-Pb collisions at $\sqrt{s_{\rm NN}} = 5.02$ TeV.

$p_{\rm T}$-differential $(\Xi^{-} + \Xi^{+}) / 2{\rm K}_{\rm S}^{0}$ ratio in jets for various event activity classes in p-Pb collisions at $\sqrt{s_{\rm NN}} = 5.02$ TeV.

$p_{\rm T}$-differential $(\Xi^{-} + \Xi^{+}) / (\Lambda + \overline{\Lambda})$ ratio of inclusive particles for various event activity classes in p-Pb collisions at $\sqrt{s_{\rm NN}} = 5.02$ TeV.

$p_{\rm T}$-differential $(\Xi^{-} + \Xi^{+}) / (\Lambda + \overline{\Lambda})$ ratio in the underlying event for $0$-$10\%$ event activity class in p-Pb collisions at $\sqrt{s_{\rm NN}} = 5.02$ TeV.

$p_{\rm T}$-differential $(\Xi^{-} + \Xi^{+}) / (\Lambda + \overline{\Lambda})$ ratio in the underlying event for $10$-$40\%$ and $40$-$100\%$ event activity classes in p-Pb collisions at $\sqrt{s_{\rm NN}} = 5.02$ TeV.

$p_{\rm T}$-differential $(\Xi^{-} + \Xi^{+}) / (\Lambda + \overline{\Lambda})$ ratio in jets for various event activity classes in p-Pb collisions at $\sqrt{s_{\rm NN}} = 5.02$ TeV.

Distributions of $\Delta\phi$ as obtained from simulation, requiring various $\textrm{t}$ tag multiplicities for the signal in two representative combinations of (gluino, neutralino) masses with large (2.2, 0.1)$\mathrm{TeV}$ and small (1.8, 1.3)$\mathrm{TeV}$ mass difference.

Results of fits to the $n_{\textrm{b}}$ multiplicity for control regions for the muon channel and with the requirements $3\leq n_{\textrm{jet}}\leq4$, $250<L_T<350~\mathrm{GeV}$, $500<H_T<750~\mathrm{GeV}$, $n_{\textrm{W}}\geq1$, $\Delta \phi<1$. The shaded area shows the fit uncertainty of the total background.

Results of fits to the $n_{\textrm{b}}$ multiplicity for control regions for the muon channel and with the requirements $3\leq n_{\textrm{jet}}\leq4$, $350<L_T<450~\mathrm{GeV}$, $H_T>1000~\mathrm{GeV}$, $n_{\textrm{W}}\geq0$, $\Delta \phi<1$. The shaded area shows the fit uncertainty of the total background.

Jet multiplicity distribution after the single-lepton baseline selection excluding the SRs for the multi-b analysis. The simulation is normalized to data with the SF mentioned in the plot.

Jet multiplicity distribution after the single-lepton baseline selection excluding the SRs for the the zero-b analysis (right). The simulation is normalized to data with the SF mentioned in the plot.

Jet multiplicity distribution in the dilepton CRs for the multi-b analysis. The simulation is normalized to data with the SF mentioned in the plot.

Jet multiplicity distribution in the dilepton CRs for the zero-b analysis. The simulation is normalized to data with the SF mentioned in the plot.

The double ratio of the single-lepton and dilepton ratio between data and simulation together with fit results and their uncertainties is shown for the multi-b analysis. The fits are performed for each data-taking year; 2018 is shown as an example.

The double ratio of the single-lepton and dilepton ratio between data and simulation together with fit results and their uncertainties is shown for the zero-b analysis. The fits are performed for each data-taking year; 2018 is shown as an example.

The prefit $L_{\mathrm{P}}$ distribution for selected electron candidates in the baseline QCD selection, with modified requirements of $n_{\text{jet}}\in[3,4]$ and $n_{\mathrm{b}}=0$.

The prefit $L_{\mathrm{P}}$ distribution for anti-selected electron candidates in the baseline QCD selection, with modified requirements of $n_{\text{jet}}\in[3,4]$ and $n_{\mathrm{b}}=0$.

Observed event yields in the MB SRs of the multi-b analysis compared to signal and background predictions. The relative fraction of the different SM EW background contributions determined in simulation is shown by the stacked, colored histograms, normalized so that their sum is equal to the background estimated using data control regions. The QCD background is predicted using the $L_p$ method. The signal is shown for two representative combinations of (gluino, neutralino) masses with large (2.2, 0.1) $\mathrm{TeV}$ and small (1.8, 1.3) $\mathrm{TeV}$ mass differences.

Observed event yields in the MB SRs of the zero-b analysis compared to signal and background predictions. The $\textrm{W}$+jets, $\textrm{t}\bar{\textrm{t}}$+jets, and QCD predictions are extracted from data control samples, while the other background contributions are estimated from simulation. The signal is shown for two representative combinations of (gluino, neutralino) masses with large (2.2, 0.1) $\mathrm{TeV}$ and small (1.8, 1.3) $\mathrm{TeV}$ mass differences.

Cross section limits at 95% CL for the T1tttt (left) model, as functions of the gluino and LSP masses, assuming a branching fraction of 100%. The mass of the intermediate chargino is taken to be halfway between the gluino and the neutralino masses. The solid black (dashed red) lines correspond to the observed (expected) mass limits, with the thicker lines representing the central values and the thinner lines representing the $\pm1\sigma$ uncertainty bands related to the theoretical (experimental) uncertainties.

Cross section limits at 95% CL for the T1tttt model, as functions of the gluino and LSP masses, assuming a branching fraction of 100%. The mass of the intermediate chargino is taken to be halfway between the gluino and the neutralino masses. The solid black (dashed red) lines correspond to the observed (expected) mass limits, with the thicker lines representing the central values and the thinner lines representing the $\pm1\sigma$ uncertainty bands related to the theoretical (experimental) uncertainties.

Cross section limits at 95% CL for the T1tttt model, as functions of the gluino and LSP masses, assuming a branching fraction of 100%. The mass of the intermediate chargino is taken to be halfway between the gluino and the neutralino masses. The solid black (dashed red) lines correspond to the observed (expected) mass limits, with the thicker lines representing the central values and the thinner lines representing the $\pm1\sigma$ uncertainty bands related to the theoretical (experimental) uncertainties.

Cross section limits at 95% CL for the T1tttt model, as functions of the gluino and LSP masses, assuming a branching fraction of 100%. The mass of the intermediate chargino is taken to be halfway between the gluino and the neutralino masses. The solid black (dashed red) lines correspond to the observed (expected) mass limits, with the thicker lines representing the central values and the thinner lines representing the $\pm1\sigma$ uncertainty bands related to the theoretical (experimental) uncertainties.

Cross section limits at 95% CL for the T1tttt model, as functions of the gluino and LSP masses, assuming a branching fraction of 100%. The mass of the intermediate chargino is taken to be halfway between the gluino and the neutralino masses. The solid black (dashed red) lines correspond to the experved (expected) mass limits, with the thicker lines representing the central values and the thinner lines representing the $\pm1\sigma$ uncertainty bands related to the theoretical (experimental) uncertainties.

Cross section limits at 95% CL for the T1tttt model, as functions of the gluino and LSP masses, assuming a branching fraction of 100%. The mass of the intermediate chargino is taken to be halfway between the gluino and the neutralino masses. The solid black (dashed red) lines correspond to the experved (expected) mass limits, with the thicker lines representing the central values and the thinner lines representing the $\pm1\sigma$ uncertainty bands related to the theoretical (experimental) uncertainties.

Cross section limits at 95% CL for the T1tttt model, as functions of the gluino and LSP masses, assuming a branching fraction of 100%. The mass of the intermediate chargino is taken to be halfway between the gluino and the neutralino masses. The solid black (dashed red) lines correspond to the experved (expected) mass limits, with the thicker lines representing the central values and the thinner lines representing the $\pm1\sigma$ uncertainty bands related to the theoretical (experimental) uncertainties.

Cross section limits at 95% CL for the T5qqqqWW model, as functions of the gluino and LSP masses, assuming a branching fraction of 100%. The mass of the intermediate chargino is taken to be halfway between the gluino and the neutralino masses. The solid black (dashed red) lines correspond to the observed (expected) mass limits, with the thicker lines representing the central values and the thinner lines representing the $\pm1\sigma$ uncertainty bands related to the theoretical (experimental) uncertainties.

Cross section limits at 95% CL for the T5qqqqWW model, as functions of the gluino and LSP masses, assuming a branching fraction of 100%. The mass of the intermediate chargino is taken to be halfway between the gluino and the neutralino masses. The solid black (dashed red) lines correspond to the observed (expected) mass limits, with the thicker lines representing the central values and the thinner lines representing the $\pm1\sigma$ uncertainty bands related to the theoretical (experimental) uncertainties.

Cross section limits at 95% CL for the T5qqqqWW model, as functions of the gluino and LSP masses, assuming a branching fraction of 100%. The mass of the intermediate chargino is taken to be halfway between the gluino and the neutralino masses. The solid black (dashed red) lines correspond to the observed (expected) mass limits, with the thicker lines representing the central values and the thinner lines representing the $\pm1\sigma$ uncertainty bands related to the theoretical (experimental) uncertainties.

Cross section limits at 95% CL for the T5qqqqWW model, as functions of the gluino and LSP masses, assuming a branching fraction of 100%. The mass of the intermediate chargino is taken to be halfway between the gluino and the neutralino masses. The solid black (dashed red) lines correspond to the observed (expected) mass limits, with the thicker lines representing the central values and the thinner lines representing the $\pm1\sigma$ uncertainty bands related to the theoretical (experimental) uncertainties.

Cross section limits at 95% CL for the T5qqqqWW model, as functions of the gluino and LSP masses, assuming a branching fraction of 100%. The mass of the intermediate chargino is taken to be halfway between the gluino and the neutralino masses. The solid black (dashed red) lines correspond to the experved (expected) mass limits, with the thicker lines representing the central values and the thinner lines representing the $\pm1\sigma$ uncertainty bands related to the theoretical (experimental) uncertainties.

Cross section limits at 95% CL for the T5qqqqWW model, as functions of the gluino and LSP masses, assuming a branching fraction of 100%. The mass of the intermediate chargino is taken to be halfway between the gluino and the neutralino masses. The solid black (dashed red) lines correspond to the experved (expected) mass limits, with the thicker lines representing the central values and the thinner lines representing the $\pm1\sigma$ uncertainty bands related to the theoretical (experimental) uncertainties.

Cross section limits at 95% CL for the T5qqqqWW model, as functions of the gluino and LSP masses, assuming a branching fraction of 100%. The mass of the intermediate chargino is taken to be halfway between the gluino and the neutralino masses. The solid black (dashed red) lines correspond to the experved (expected) mass limits, with the thicker lines representing the central values and the thinner lines representing the $\pm1\sigma$ uncertainty bands related to the theoretical (experimental) uncertainties.

Observed number of events in the MB SR bins of the multi-b analysis. All bins are defined with $\Delta\phi > 0.75$.

Observed number of events in the MB SR bins of the zero-b analysis.

The simplified likelihood can be constructed using this covariance matrix of nuisance parameters in the signal region (MB SR). The bin contents have been aggregated accross the three years to simplify the likelihood.

Selection efficiency for the multi-b signal for the SRs.

The simplified likelihood can be constructed using this covariance matrix of nuisance parameters in the signal region (MB SR). The bins were aggregated across HT to stabilize the simplified likelihood.

Selection efficiency for the zero-b signal for the SRs. The bins correspond to the aggregated bins that were used for the covariance matrix.

A summary of the uncertainties in the background estimates for SR-Gtt-0L-C. The individual experimental and theoretical uncertainties are assumed to be uncorrelated and are combined by adding in quadrature.

A summary of the uncertainties in the background estimates for SR-Gtt-1L-B. The individual experimental and theoretical uncertainties are assumed to be uncorrelated and are combined by adding in quadrature.

A summary of the uncertainties in the background estimates for SR-Gtt-1L-M1. The individual experimental and theoretical uncertainties are assumed to be uncorrelated and are combined by adding in quadrature.

A summary of the uncertainties in the background estimates for SR-Gtt-1L-M2. The individual experimental and theoretical uncertainties are assumed to be uncorrelated and are combined by adding in quadrature.

A summary of the uncertainties in the background estimates for SR-Gtt-1L-C. The individual experimental and theoretical uncertainties are assumed to be uncorrelated and are combined by adding in quadrature.

A summary of the uncertainties in the background estimates for SR-Gbb-B. The individual experimental and theoretical uncertainties are assumed to be uncorrelated and are combined by adding in quadrature.

A summary of the uncertainties in the background estimates for SR-Gbb-M. The individual experimental and theoretical uncertainties are assumed to be uncorrelated and are combined by adding in quadrature.

A summary of the uncertainties in the background estimates for SR-Gbb-C. The individual experimental and theoretical uncertainties are assumed to be uncorrelated and are combined by adding in quadrature.

A summary of the uncertainties in the background estimates for SR-Gtb-B. The individual experimental and theoretical uncertainties are assumed to be uncorrelated and are combined by adding in quadrature.

A summary of the uncertainties in the background estimates for SR-Gtb-M. The individual experimental and theoretical uncertainties are assumed to be uncorrelated and are combined by adding in quadrature.

A summary of the uncertainties in the background estimates for SR-Gtb-C. The individual experimental and theoretical uncertainties are assumed to be uncorrelated and are combined by adding in quadrature.

A summary of the uncertainties in the background estimates for SR-Gtt-2100-1. The individual experimental and theoretical uncertainties are assumed to be uncorrelated and are combined by adding in quadrature.

A summary of the uncertainties in the background estimates for SR-Gtt-1800-1. The individual experimental and theoretical uncertainties are assumed to be uncorrelated and are combined by adding in quadrature.

A summary of the uncertainties in the background estimates for SR-Gtt-2300-1200. The individual experimental and theoretical uncertainties are assumed to be uncorrelated and are combined by adding in quadrature.

A summary of the uncertainties in the background estimates for SR-Gtt-1900-1400. The individual experimental and theoretical uncertainties are assumed to be uncorrelated and are combined by adding in quadrature.

A summary of the uncertainties in the background estimates for SR-Gbb-2800-1400. The individual experimental and theoretical uncertainties are assumed to be uncorrelated and are combined by adding in quadrature.

A summary of the uncertainties in the background estimates for SR-Gbb-2300-1000. The individual experimental and theoretical uncertainties are assumed to be uncorrelated and are combined by adding in quadrature.

A summary of the uncertainties in the background estimates for SR-Gbb-2100-1600. The individual experimental and theoretical uncertainties are assumed to be uncorrelated and are combined by adding in quadrature.

A summary of the uncertainties in the background estimates for SR-Gbb-2000-1800. The individual experimental and theoretical uncertainties are assumed to be uncorrelated and are combined by adding in quadrature.

Results of the background-only fit extrapolated to SR_Gtt_0L_B in the CC analysis, for both the total expected background yields and the main contributing background processes. The quoted uncertainties include both experimental and theoretical systematics. The data in the SRs are not included in the fit. The background category $t\bar{t}+X$ includes $t\bar{t} W/Z$, $t\bar{t} H$ and $t\bar{t} t\bar{t}$ events. The row ``Pre-fit background'' provides the total background prediction when the $t\bar{t}$ and $Z+$jets normalisations are obtained from theoretical calculation, taking into account the kinematic weights described in Section 5.

Results of the background-only fit extrapolated to SR_Gtt_0L_M1 in the CC analysis, for both the total expected background yields and the main contributing background processes. The quoted uncertainties include both experimental and theoretical systematics. The data in the SRs are not included in the fit. The background category $t\bar{t}+X$ includes $t\bar{t} W/Z$, $t\bar{t} H$ and $t\bar{t} t\bar{t}$ events. The row ``Pre-fit background'' provides the total background prediction when the $t\bar{t}$ and $Z+$jets normalisations are obtained from theoretical calculation, taking into account the kinematic weights described in Section 5.

Results of the background-only fit extrapolated to SR_Gtt_0L_M2 in the CC analysis, for both the total expected background yields and the main contributing background processes. The quoted uncertainties include both experimental and theoretical systematics. The data in the SRs are not included in the fit. The background category $t\bar{t}+X$ includes $t\bar{t} W/Z$, $t\bar{t} H$ and $t\bar{t} t\bar{t}$ events. The row ``Pre-fit background'' provides the total background prediction when the $t\bar{t}$ and $Z+$jets normalisations are obtained from theoretical calculation, taking into account the kinematic weights described in Section 5.

Results of the background-only fit extrapolated to SR_Gtt_0L_C in the CC analysis, for both the total expected background yields and the main contributing background processes. The quoted uncertainties include both experimental and theoretical systematics. The data in the SRs are not included in the fit. The background category $t\bar{t}+X$ includes $t\bar{t} W/Z$, $t\bar{t} H$ and $t\bar{t} t\bar{t}$ events. The row ``Pre-fit background'' provides the total background prediction when the $t\bar{t}$ and $Z+$jets normalisations are obtained from theoretical calculation, taking into account the kinematic weights described in Section 5.

Results of the background-only fit extrapolated to SR_Gtt_1L_B in the CC analysis, for both the total expected background yields and the main contributing background processes. The quoted uncertainties include both experimental and theoretical systematics. The data in the SRs are not included in the fit. The background category $t\bar{t}+X$ includes $t\bar{t} W/Z$, $t\bar{t} H$ and $t\bar{t} t\bar{t}$ events. The row ``Pre-fit background'' provides the total background prediction when the $t\bar{t}$ and $Z+$jets normalisations are obtained from theoretical calculation, taking into account the kinematic weights described in Section 5.

Results of the background-only fit extrapolated to SR_Gtt_1L_M1 in the CC analysis, for both the total expected background yields and the main contributing background processes. The quoted uncertainties include both experimental and theoretical systematics. The data in the SRs are not included in the fit. The background category $t\bar{t}+X$ includes $t\bar{t} W/Z$, $t\bar{t} H$ and $t\bar{t} t\bar{t}$ events. The row ``Pre-fit background'' provides the total background prediction when the $t\bar{t}$ and $Z+$jets normalisations are obtained from theoretical calculation, taking into account the kinematic weights described in Section 5.

Results of the background-only fit extrapolated to SR_Gtt_1L_M2 in the CC analysis, for both the total expected background yields and the main contributing background processes. The quoted uncertainties include both experimental and theoretical systematics. The data in the SRs are not included in the fit. The background category $t\bar{t}+X$ includes $t\bar{t} W/Z$, $t\bar{t} H$ and $t\bar{t} t\bar{t}$ events. The row ``Pre-fit background'' provides the total background prediction when the $t\bar{t}$ and $Z+$jets normalisations are obtained from theoretical calculation, taking into account the kinematic weights described in Section 5.

Results of the background-only fit extrapolated to SR_Gtt_1L_C in the CC analysis, for both the total expected background yields and the main contributing background processes. The quoted uncertainties include both experimental and theoretical systematics. The data in the SRs are not included in the fit. The background category $t\bar{t}+X$ includes $t\bar{t} W/Z$, $t\bar{t} H$ and $t\bar{t} t\bar{t}$ events. The row ``Pre-fit background'' provides the total background prediction when the $t\bar{t}$ and $Z+$jets normalisations are obtained from theoretical calculation, taking into account the kinematic weights described in Section 5.

Results of the background-only fit extrapolated to SR_Gbb_B in the CC analysis, for both the total expected background yields and the main contributing background processes. The quoted uncertainties include both experimental and theoretical systematics. The data in the SRs are not included in the fit. The background category $t\bar{t}+X$ includes $t\bar{t} W/Z$, $t\bar{t} H$ and $t\bar{t} t\bar{t}$ events. The row ``Pre-fit background'' provides the total background prediction when the $t\bar{t}$ and $Z+$jets normalisations are obtained from theoretical calculation, taking into account the kinematic weights described in Section 5.

Results of the background-only fit extrapolated to SR_Gbb_M in the CC analysis, for both the total expected background yields and the main contributing background processes. The quoted uncertainties include both experimental and theoretical systematics. The data in the SRs are not included in the fit. The background category $t\bar{t}+X$ includes $t\bar{t} W/Z$, $t\bar{t} H$ and $t\bar{t} t\bar{t}$ events. The row ``Pre-fit background'' provides the total background prediction when the $t\bar{t}$ and $Z+$jets normalisations are obtained from theoretical calculation, taking into account the kinematic weights described in Section 5.

Results of the background-only fit extrapolated to SR_Gbb_C in the CC analysis, for both the total expected background yields and the main contributing background processes. The quoted uncertainties include both experimental and theoretical systematics. The data in the SRs are not included in the fit. The background category $t\bar{t}+X$ includes $t\bar{t} W/Z$, $t\bar{t} H$ and $t\bar{t} t\bar{t}$ events. The row ``Pre-fit background'' provides the total background prediction when the $t\bar{t}$ and $Z+$jets normalisations are obtained from theoretical calculation, taking into account the kinematic weights described in Section 5.

Results of the background-only fit extrapolated to SR_Gtb_B in the CC analysis, for both the total expected background yields and the main contributing background processes. The quoted uncertainties include both experimental and theoretical systematics. The data in the SRs are not included in the fit. The background category $t\bar{t}+X$ includes $t\bar{t} W/Z$, $t\bar{t} H$ and $t\bar{t} t\bar{t}$ events. The row ``Pre-fit background'' provides the total background prediction when the $t\bar{t}$ and $Z+$jets normalisations are obtained from theoretical calculation, taking into account the kinematic weights described in Section 5.

Results of the background-only fit extrapolated to SR_Gtb_M in the CC analysis, for both the total expected background yields and the main contributing background processes. The quoted uncertainties include both experimental and theoretical systematics. The data in the SRs are not included in the fit. The background category $t\bar{t}+X$ includes $t\bar{t} W/Z$, $t\bar{t} H$ and $t\bar{t} t\bar{t}$ events. The row ``Pre-fit background'' provides the total background prediction when the $t\bar{t}$ and $Z+$jets normalisations are obtained from theoretical calculation, taking into account the kinematic weights described in Section 5.

Results of the background-only fit extrapolated to SR_Gtb_C in the CC analysis, for both the total expected background yields and the main contributing background processes. The quoted uncertainties include both experimental and theoretical systematics. The data in the SRs are not included in the fit. The background category $t\bar{t}+X$ includes $t\bar{t} W/Z$, $t\bar{t} H$ and $t\bar{t} t\bar{t}$ events. The row ``Pre-fit background'' provides the total background prediction when the $t\bar{t}$ and $Z+$jets normalisations are obtained from theoretical calculation, taking into account the kinematic weights described in Section 5.

Results of the background-only fit extrapolated to SR_Gtt_2100_1 in the NN analysis, for both the total expected background yields and the main contributing background processes. The quoted uncertainties include both experimental and theoretical systematics. The data in the SRs are not included in the fit. The background category $t\bar{t}+X$ includes $t\bar{t} W/Z$, $t\bar{t} H$ and $t\bar{t} t\bar{t}$ events. The row ``Pre-fit background'' provides the total background prediction when the $t\bar{t}$ and $Z+$jets normalisations are obtained from theoretical calculation, taking into account the kinematic weights described in Section 5.

Results of the background-only fit extrapolated to SR_Gtt_1800_1 in the NN analysis, for both the total expected background yields and the main contributing background processes. The quoted uncertainties include both experimental and theoretical systematics. The data in the SRs are not included in the fit. The background category $t\bar{t}+X$ includes $t\bar{t} W/Z$, $t\bar{t} H$ and $t\bar{t} t\bar{t}$ events. The row ``Pre-fit background'' provides the total background prediction when the $t\bar{t}$ and $Z+$jets normalisations are obtained from theoretical calculation, taking into account the kinematic weights described in Section 5.

Results of the background-only fit extrapolated to SR_Gtt_2300_1200 in the NN analysis, for both the total expected background yields and the main contributing background processes. The quoted uncertainties include both experimental and theoretical systematics. The data in the SRs are not included in the fit. The background category $t\bar{t}+X$ includes $t\bar{t} W/Z$, $t\bar{t} H$ and $t\bar{t} t\bar{t}$ events. The row ``Pre-fit background'' provides the total background prediction when the $t\bar{t}$ and $Z+$jets normalisations are obtained from theoretical calculation, taking into account the kinematic weights described in Section 5.

Results of the background-only fit extrapolated to SR_Gtt_1900_1400 in the NN analysis, for both the total expected background yields and the main contributing background processes. The quoted uncertainties include both experimental and theoretical systematics. The data in the SRs are not included in the fit. The background category $t\bar{t}+X$ includes $t\bar{t} W/Z$, $t\bar{t} H$ and $t\bar{t} t\bar{t}$ events. The row ``Pre-fit background'' provides the total background prediction when the $t\bar{t}$ and $Z+$jets normalisations are obtained from theoretical calculation, taking into account the kinematic weights described in Section 5.

Results of the background-only fit extrapolated to SR_Gbb_2800_1400 in the NN analysis, for both the total expected background yields and the main contributing background processes. The quoted uncertainties include both experimental and theoretical systematics. The data in the SRs are not included in the fit. The background category $t\bar{t}+X$ includes $t\bar{t} W/Z$, $t\bar{t} H$ and $t\bar{t} t\bar{t}$ events. The row ``Pre-fit background'' provides the total background prediction when the $t\bar{t}$ and $Z+$jets normalisations are obtained from theoretical calculation, taking into account the kinematic weights described in Section 5.

Results of the background-only fit extrapolated to SR_Gbb_2300_1000 in the NN analysis, for both the total expected background yields and the main contributing background processes. The quoted uncertainties include both experimental and theoretical systematics. The data in the SRs are not included in the fit. The background category $t\bar{t}+X$ includes $t\bar{t} W/Z$, $t\bar{t} H$ and $t\bar{t} t\bar{t}$ events. The row ``Pre-fit background'' provides the total background prediction when the $t\bar{t}$ and $Z+$jets normalisations are obtained from theoretical calculation, taking into account the kinematic weights described in Section 5.

Results of the background-only fit extrapolated to SR_Gbb_2100_1600 in the NN analysis, for both the total expected background yields and the main contributing background processes. The quoted uncertainties include both experimental and theoretical systematics. The data in the SRs are not included in the fit. The background category $t\bar{t}+X$ includes $t\bar{t} W/Z$, $t\bar{t} H$ and $t\bar{t} t\bar{t}$ events. The row ``Pre-fit background'' provides the total background prediction when the $t\bar{t}$ and $Z+$jets normalisations are obtained from theoretical calculation, taking into account the kinematic weights described in Section 5.

Results of the background-only fit extrapolated to SR_Gbb_2000_1800 in the NN analysis, for both the total expected background yields and the main contributing background processes. The quoted uncertainties include both experimental and theoretical systematics. The data in the SRs are not included in the fit. The background category $t\bar{t}+X$ includes $t\bar{t} W/Z$, $t\bar{t} H$ and $t\bar{t} t\bar{t}$ events. The row ``Pre-fit background'' provides the total background prediction when the $t\bar{t}$ and $Z+$jets normalisations are obtained from theoretical calculation, taking into account the kinematic weights described in Section 5.

Observed exclusion limit in the $\tilde{g}$--$\tilde\chi^0_1$ mass plane for the Gtt models obtained from the NN analysis. The dashed and solid bold lines show the 95\% CL expected and observed limits, respectively. The shaded bands around the expected limits show the impact of the experimental and background uncertainties. The dotted lines show the impact on the observed limit of the variation of the nominal signal cross-section by $\pm 1 \sigma$ of its theoretical uncertainty.

Observed exclusion limit in the $\tilde{g}$--$\tilde\chi^0_1$ mass plane for the Gtt models obtained from the NN analysis. The dashed and solid bold lines show the 95\% CL expected and observed limits, respectively. The shaded bands around the expected limits show the impact of the experimental and background uncertainties. The dotted lines show the impact on the observed limit of the variation of the nominal signal cross-section by $\pm 1 \sigma$ of its theoretical uncertainty.

Observed exclusion limit in the $\tilde{g}$--$\tilde\chi^0_1$ mass plane for the Gtt models obtained from the NN analysis. The dashed and solid bold lines show the 95\% CL expected and observed limits, respectively. The shaded bands around the expected limits show the impact of the experimental and background uncertainties. The dotted lines show the impact on the observed limit of the variation of the nominal signal cross-section by $\pm 1 \sigma$ of its theoretical uncertainty.

Expected exclusion limit in the $\tilde{g}$--$\tilde\chi^0_1$ mass plane for the Gtt models obtained from the NN analysis. The dashed and solid bold lines show the 95\% CL expected and observed limits, respectively. The shaded bands around the expected limits show the impact of the experimental and background uncertainties. The dotted lines show the impact on the observed limit of the variation of the nominal signal cross-section by $\pm 1 \sigma$ of its theoretical uncertainty.

Expected exclusion limit in the $\tilde{g}$--$\tilde\chi^0_1$ mass plane for the Gtt models obtained from the NN analysis. The dashed and solid bold lines show the 95\% CL expected and observed limits, respectively. The shaded bands around the expected limits show the impact of the experimental and background uncertainties. The dotted lines show the impact on the observed limit of the variation of the nominal signal cross-section by $\pm 1 \sigma$ of its theoretical uncertainty.

Expected exclusion limit in the $\tilde{g}$--$\tilde\chi^0_1$ mass plane for the Gtt models obtained from the NN analysis. The dashed and solid bold lines show the 95\% CL expected and observed limits, respectively. The shaded bands around the expected limits show the impact of the experimental and background uncertainties. The dotted lines show the impact on the observed limit of the variation of the nominal signal cross-section by $\pm 1 \sigma$ of its theoretical uncertainty.

Observed exclusion limit in the $\tilde{g}$--$\tilde\chi^0_1$ mass plane for the Gbb models obtained from the NN analysis. The dashed and solid bold lines show the 95\% CL expected and observed limits, respectively. The shaded bands around the expected limits show the impact of the experimental and background uncertainties. The dotted lines show the impact on the observed limit of the variation of the nominal signal cross-section by $\pm 1 \sigma$ of its theoretical uncertainty.

Observed exclusion limit in the $\tilde{g}$--$\tilde\chi^0_1$ mass plane for the Gbb models obtained from the NN analysis. The dashed and solid bold lines show the 95\% CL expected and observed limits, respectively. The shaded bands around the expected limits show the impact of the experimental and background uncertainties. The dotted lines show the impact on the observed limit of the variation of the nominal signal cross-section by $\pm 1 \sigma$ of its theoretical uncertainty.

Observed exclusion limit in the $\tilde{g}$--$\tilde\chi^0_1$ mass plane for the Gbb models obtained from the NN analysis. The dashed and solid bold lines show the 95\% CL expected and observed limits, respectively. The shaded bands around the expected limits show the impact of the experimental and background uncertainties. The dotted lines show the impact on the observed limit of the variation of the nominal signal cross-section by $\pm 1 \sigma$ of its theoretical uncertainty.

Expected exclusion limit in the $\tilde{g}$--$\tilde\chi^0_1$ mass plane for the Gbb models obtained from the NN analysis. The dashed and solid bold lines show the 95\% CL expected and observed limits, respectively. The shaded bands around the expected limits show the impact of the experimental and background uncertainties. The dotted lines show the impact on the observed limit of the variation of the nominal signal cross-section by $\pm 1 \sigma$ of its theoretical uncertainty.

Expected exclusion limit in the $\tilde{g}$--$\tilde\chi^0_1$ mass plane for the Gbb models obtained from the NN analysis. The dashed and solid bold lines show the 95\% CL expected and observed limits, respectively. The shaded bands around the expected limits show the impact of the experimental and background uncertainties. The dotted lines show the impact on the observed limit of the variation of the nominal signal cross-section by $\pm 1 \sigma$ of its theoretical uncertainty.

Expected exclusion limit in the $\tilde{g}$--$\tilde\chi^0_1$ mass plane for the Gbb models obtained from the NN analysis. The dashed and solid bold lines show the 95\% CL expected and observed limits, respectively. The shaded bands around the expected limits show the impact of the experimental and background uncertainties. The dotted lines show the impact on the observed limit of the variation of the nominal signal cross-section by $\pm 1 \sigma$ of its theoretical uncertainty.

Observed (left) 95\%~CL exclusion limits on the gluino mass as a function of BR$(\tilde{g} \to b\bar{b}\tilde\chi^{0}_{1}$) (vertical) and BR$(\tilde{g} \to t\bar{t}\tilde\chi^{0}_{1}$) (horizontal) for Gtb models with $m(\tilde\chi^{0}_{1}) = 1$~GeV, obtained from the CC analysis.

Expected (right) 95\%~CL exclusion limits on the gluino mass as a function of BR$(\tilde{g} \to b\bar{b}\tilde\chi^{0}_{1}$) (vertical) and BR$(\tilde{g} \to t\bar{t}\tilde\chi^{0}_{1}$) (horizontal) for Gtb models with $m(\tilde\chi^{0}_{1}) = 1$~GeV, obtained from the CC analysis.

Observed (left) 95\%~CL exclusion limits on the gluino mass as a function of BR$(\tilde{g} \to b\bar{b}\tilde\chi^{0}_{1}$) (vertical) and BR$(\tilde{g} \to t\bar{t}\tilde\chi^{0}_{1}$) (horizontal) for Gtb models with $m(\tilde\chi^{0}_{1}) = 600$~GeV, obtained from the CC analysis.

Expected (right) 95\%~CL exclusion limits on the gluino mass as a function of BR$(\tilde{g} \to b\bar{b}\tilde\chi^{0}_{1}$) (vertical) and BR$(\tilde{g} \to t\bar{t}\tilde\chi^{0}_{1}$) (horizontal) for Gtb models with $m(\tilde\chi^{0}_{1}) = 600$~GeV, obtained from the CC analysis.

Observed (left) 95\%~CL exclusion limits on the gluino mass as a function of BR$(\tilde{g} \to b\bar{b}\tilde\chi^{0}_{1}$) (vertical) and BR$(\tilde{g} \to t\bar{t}\tilde\chi^{0}_{1}$) (horizontal) for Gtb models with $m(\tilde\chi^{0}_{1}) = 1$~TeV, obtained from the CC analysis.

Expected (right) 95\%~CL exclusion limits on the gluino mass as a function of BR$(\tilde{g} \to b\bar{b}\tilde\chi^{0}_{1}$) (vertical) and BR$(\tilde{g} \to t\bar{t}\tilde\chi^{0}_{1}$) (horizontal) for Gtb models with $m(\tilde\chi^{0}_{1}) = 1$~TeV, obtained from the CC analysis.

Observed exclusion limit in the $\tilde{g}$--$\tilde\chi^0_1$ mass plane for the Gtt models obtained from the CC analysis. The dashed and solid bold lines show the 95\% CL expected and observed limits, respectively. The shaded bands around the expected limits show the impact of the experimental and background uncertainties. The dotted lines show the impact on the observed limit of the variation of the nominal signal cross-section by $\pm 1 \sigma$ of its theoretical uncertainty.

Observed exclusion limit in the $\tilde{g}$--$\tilde\chi^0_1$ mass plane for the Gtt models obtained from the CC analysis. The dashed and solid bold lines show the 95\% CL expected and observed limits, respectively. The shaded bands around the expected limits show the impact of the experimental and background uncertainties. The dotted lines show the impact on the observed limit of the variation of the nominal signal cross-section by $\pm 1 \sigma$ of its theoretical uncertainty.

Observed exclusion limit in the $\tilde{g}$--$\tilde\chi^0_1$ mass plane for the Gtt models obtained from the CC analysis. The dashed and solid bold lines show the 95\% CL expected and observed limits, respectively. The shaded bands around the expected limits show the impact of the experimental and background uncertainties. The dotted lines show the impact on the observed limit of the variation of the nominal signal cross-section by $\pm 1 \sigma$ of its theoretical uncertainty.

Expected exclusion limit in the $\tilde{g}$--$\tilde\chi^0_1$ mass plane for the Gtt models obtained from the CC analysis. The dashed and solid bold lines show the 95\% CL expected and observed limits, respectively. The shaded bands around the expected limits show the impact of the experimental and background uncertainties. The dotted lines show the impact on the observed limit of the variation of the nominal signal cross-section by $\pm 1 \sigma$ of its theoretical uncertainty.

Expected exclusion limit in the $\tilde{g}$--$\tilde\chi^0_1$ mass plane for the Gtt models obtained from the CC analysis. The dashed and solid bold lines show the 95\% CL expected and observed limits, respectively. The shaded bands around the expected limits show the impact of the experimental and background uncertainties. The dotted lines show the impact on the observed limit of the variation of the nominal signal cross-section by $\pm 1 \sigma$ of its theoretical uncertainty.

Expected exclusion limit in the $\tilde{g}$--$\tilde\chi^0_1$ mass plane for the Gtt models obtained from the CC analysis. The dashed and solid bold lines show the 95\% CL expected and observed limits, respectively. The shaded bands around the expected limits show the impact of the experimental and background uncertainties. The dotted lines show the impact on the observed limit of the variation of the nominal signal cross-section by $\pm 1 \sigma$ of its theoretical uncertainty.

Observed exclusion limit in the $\tilde{g}$--$\tilde\chi^0_1$ mass plane for the Gbb models obtained from the CC analysis. The dashed and solid bold lines show the 95\% CL expected and observed limits, respectively. The shaded bands around the expected limits show the impact of the experimental and background uncertainties. The dotted lines show the impact on the observed limit of the variation of the nominal signal cross-section by $\pm 1 \sigma$ of its theoretical uncertainty.

Observed exclusion limit in the $\tilde{g}$--$\tilde\chi^0_1$ mass plane for the Gbb models obtained from the CC analysis. The dashed and solid bold lines show the 95\% CL expected and observed limits, respectively. The shaded bands around the expected limits show the impact of the experimental and background uncertainties. The dotted lines show the impact on the observed limit of the variation of the nominal signal cross-section by $\pm 1 \sigma$ of its theoretical uncertainty.

Observed exclusion limit in the $\tilde{g}$--$\tilde\chi^0_1$ mass plane for the Gbb models obtained from the CC analysis. The dashed and solid bold lines show the 95\% CL expected and observed limits, respectively. The shaded bands around the expected limits show the impact of the experimental and background uncertainties. The dotted lines show the impact on the observed limit of the variation of the nominal signal cross-section by $\pm 1 \sigma$ of its theoretical uncertainty.

Expected exclusion limit in the $\tilde{g}$--$\tilde\chi^0_1$ mass plane for the Gbb models obtained from the CC analysis. The dashed and solid bold lines show the 95\% CL expected and observed limits, respectively. The shaded bands around the expected limits show the impact of the experimental and background uncertainties. The dotted lines show the impact on the observed limit of the variation of the nominal signal cross-section by $\pm 1 \sigma$ of its theoretical uncertainty.

Expected exclusion limit in the $\tilde{g}$--$\tilde\chi^0_1$ mass plane for the Gbb models obtained from the CC analysis. The dashed and solid bold lines show the 95\% CL expected and observed limits, respectively. The shaded bands around the expected limits show the impact of the experimental and background uncertainties. The dotted lines show the impact on the observed limit of the variation of the nominal signal cross-section by $\pm 1 \sigma$ of its theoretical uncertainty.

Expected exclusion limit in the $\tilde{g}$--$\tilde\chi^0_1$ mass plane for the Gbb models obtained from the CC analysis. The dashed and solid bold lines show the 95\% CL expected and observed limits, respectively. The shaded bands around the expected limits show the impact of the experimental and background uncertainties. The dotted lines show the impact on the observed limit of the variation of the nominal signal cross-section by $\pm 1 \sigma$ of its theoretical uncertainty.

Upper limit at 95\% CL on the cross-section times branching ratio (fb) in the $\tilde{g}$--$\tilde\chi^0_1$ mass plane for the Gbb (right) models obtained from the CC analysis. The numbers give the observed 95\% CL upper limit on the cross section in fb, with the label colour matching the associated best-expected region. Only a lower limit on the excluded cross section (>0.7 fb) is given at some points due to the very small number events expected and observed in the chosen SR. The dashed and solid bold lines show the 95\% CL expected and observed limits, respectively. The shaded bands around the expected limits show the impact of the experimental and background theoretical uncertainties. The dotted lines show the impact on the observed limit of the variation of the nominal signal cross-section by $\pm1\sigma$ of its theoretical uncertainty.

Upper limit at 95\% CL on the cross-section times branching ratio (fb) in the $\tilde{g}$--$\tilde\chi^0_1$ mass plane for the Gbb (right) models obtained from the NN analysis. The numbers give the observed 95\% CL upper limit on the cross section in fb, with the label colour matching the associated best-expected region. Only a lower limit on the excluded cross section (>0.7 fb) is given at some points due to the very small number events expected and observed in the chosen SR. The dashed and solid bold lines show the 95\% CL expected and observed limits, respectively. The shaded bands around the expected limits show the impact of the experimental and background theoretical uncertainties. The dotted lines show the impact on the observed limit of the variation of the nominal signal cross-section by $\pm1\sigma$ of its theoretical uncertainty.

Upper limit at 95\% CL on the cross-section times branching ratio (fb) in the $\tilde{g}$--$\tilde\chi^0_1$ mass plane for the Gtt (left) models obtained from the CC analysis. The numbers give the observed 95\% CL upper limit on the cross section in fb, with the label colour matching the associated best-expected region. Only a lower limit on the excluded cross section (>0.7 fb) is given at some points due to the very small number events expected and observed in the chosen SR. The dashed and solid bold lines show the 95\% CL expected and observed limits, respectively. The shaded bands around the expected limits show the impact of the experimental and background theoretical uncertainties. The dotted lines show the impact on the observed limit of the variation of the nominal signal cross-section by $\pm1\sigma$ of its theoretical uncertainty.

Upper limit at 95\% CL on the cross-section times branching ratio (fb) in the $\tilde{g}$--$\tilde\chi^0_1$ mass plane for the Gtt (left) models obtained from the NN analysis. The numbers give the observed 95\% CL upper limit on the cross section in fb, with the label colour matching the associated best-expected region. Only a lower limit on the excluded cross section (>0.7 fb) is given at some points due to the very small number events expected and observed in the chosen SR. The dashed and solid bold lines show the 95\% CL expected and observed limits, respectively. The shaded bands around the expected limits show the impact of the experimental and background theoretical uncertainties. The dotted lines show the impact on the observed limit of the variation of the nominal signal cross-section by $\pm1\sigma$ of its theoretical uncertainty.

Acceptance for SR-Gtt-0L-B and the $\tilde{g}\rightarrow t\bar{t}\tilde\chi^0_1$ signal process.

Efficiency for SR-Gtt-0L-B and the $\tilde{g}\rightarrow t\bar{t}\tilde\chi^0_1$ signal process.

Acceptance for SR-Gtt-0L-M1 and the $\tilde{g}\rightarrow t\bar{t}\tilde\chi^0_1$ signal process.

Efficiency for SR-Gtt-0L-M1 and the $\tilde{g}\rightarrow t\bar{t}\tilde\chi^0_1$ signal process.

Acceptance for SR-Gtt-0L-M2 and the $\tilde{g}\rightarrow t\bar{t}\tilde\chi^0_1$ signal process.

Efficiency for SR-Gtt-0L-M2 and the $\tilde{g}\rightarrow t\bar{t}\tilde\chi^0_1$ signal process.

Acceptance for SR-Gtt-0L-C and the $\tilde{g}\rightarrow t\bar{t}\tilde\chi^0_1$ signal process.

Efficiency for SR-Gtt-0L-C and the $\tilde{g}\rightarrow t\bar{t}\tilde\chi^0_1$ signal process.

Acceptance for SR-Gtt-1L-B and the $\tilde{g}\rightarrow t\bar{t}\tilde\chi^0_1$ signal process.

Efficiency for SR-Gtt-1L-B and the $\tilde{g}\rightarrow t\bar{t}\tilde\chi^0_1$ signal process.

Acceptance for SR-Gtt-1L-M1 and the $\tilde{g}\rightarrow t\bar{t}\tilde\chi^0_1$ signal process.

Efficiency for SR-Gtt-1L-M1 and the $\tilde{g}\rightarrow t\bar{t}\tilde\chi^0_1$ signal process.

Acceptance for SR-Gtt-1L-M2 and the $\tilde{g}\rightarrow t\bar{t}\tilde\chi^0_1$ signal process.

Efficiency for SR-Gtt-1L-M2 and the $\tilde{g}\rightarrow t\bar{t}\tilde\chi^0_1$ signal process.

Acceptance for SR-Gtt-1L-C and the $\tilde{g}\rightarrow t\bar{t}\tilde\chi^0_1$ signal process.

Efficiency for SR-Gtt-1L-C and the $\tilde{g}\rightarrow t\bar{t}\tilde\chi^0_1$ signal process.

Acceptance for SR-Gbb-B and the $\tilde{g}\rightarrow b\bar{b}\tilde\chi^0_1$ signal process.

Efficiency for SR-Gbb-B and the $\tilde{g}\rightarrow b\bar{b}\tilde\chi^0_1$ signal process.

Acceptance for SR-Gbb-M and the $\tilde{g}\rightarrow b\bar{b}\tilde\chi^0_1$ signal process.

Efficiency for SR-Gbb-M and the $\tilde{g}\rightarrow b\bar{b}\tilde\chi^0_1$ signal process.

Acceptance for SR-Gbb-C and the $\tilde{g}\rightarrow b\bar{b}\tilde\chi^0_1$ signal process.

Efficiency for SR-Gbb-C and the $\tilde{g}\rightarrow b\bar{b}\tilde\chi^0_1$ signal process.

Acceptance for SR-Gtt-2100-1 and the $\tilde{g}\rightarrow t\bar{t}\tilde\chi^0_1$ signal process.

Efficiency for SR-Gtt-2100-1 and the $\tilde{g}\rightarrow t\bar{t}\tilde\chi^0_1$ signal process.

Acceptance for SR-Gtt-1800-1 and the $\tilde{g}\rightarrow t\bar{t}\tilde\chi^0_1$ signal process.

Efficiency for SR-Gtt-1800-1 and the $\tilde{g}\rightarrow t\bar{t}\tilde\chi^0_1$ signal process.

Acceptance for SR-Gtt-2300-1200 and the $\tilde{g}\rightarrow t\bar{t}\tilde\chi^0_1$ signal process.

Efficiency for SR-Gtt-2300-1200 and the $\tilde{g}\rightarrow t\bar{t}\tilde\chi^0_1$ signal process.

Acceptance for SR-Gtt-1900-1400 and the $\tilde{g}\rightarrow t\bar{t}\tilde\chi^0_1$ signal process.

Efficiency for SR-Gtt-1900-1400 and the $\tilde{g}\rightarrow t\bar{t}\tilde\chi^0_1$ signal process.

Acceptance for SR-Gbb-2800-1400 and the $\tilde{g}\rightarrow b\bar{b}\tilde\chi^0_1$ signal process.

Efficiency for SR-Gbb-2800-1400 and the $\tilde{g}\rightarrow b\bar{b}\tilde\chi^0_1$ signal process.

Acceptance for SR-Gbb-2300-1000 and the $\tilde{g}\rightarrow b\bar{b}\tilde\chi^0_1$ signal process.

Efficiency for SR-Gbb-2300-1000 and the $\tilde{g}\rightarrow b\bar{b}\tilde\chi^0_1$ signal process.

Acceptance for SR-Gbb-2100-1600 and the $\tilde{g}\rightarrow b\bar{b}\tilde\chi^0_1$ signal process.

Efficiency for SR-Gbb-2100-1600 and the $\tilde{g}\rightarrow b\bar{b}\tilde\chi^0_1$ signal process.

Acceptance for SR-Gbb-2000-1800 and the $\tilde{g}\rightarrow b\bar{b}\tilde\chi^0_1$ signal process.

Efficiency for SR-Gbb-2000-1800 and the $\tilde{g}\rightarrow b\bar{b}\tilde\chi^0_1$ signal process.

Cutflow for the SR-Gtt-0L-B for a representative Gtt signal. Signal was generated with 30000 events. Expected yields are normalised to a luminosity of 139~fb$^{-1}$.

Cutflow for the SR-Gtt-0L-M1 for a representative Gtt signal. Signal was generated with 30000 events. Expected yields are normalised to a luminosity of 139~fb$^{-1}$.

Cutflow for the SR-Gtt-0L-M2 for a representative Gtt signal. Signal was generated with 30000 events. Expected yields are normalised to a luminosity of 139~fb$^{-1}$.

Cutflow for the SR-Gtt-0L-C for a representative Gtt signal. Signal was generated with 30000 events. Expected yields are normalised to a luminosity of 139~fb$^{-1}$.

Cutflow for the SR-Gtt-1L-B for a representative Gtt signal. Signal was generated with 30000 events. Expected yields are normalised to a luminosity of 139~fb$^{-1}$.

Cutflow for the SR-Gtt-1L-M1 for a representative Gtt signal. Signal was generated with 30000 events. Expected yields are normalised to a luminosity of 139~fb$^{-1}$.

Cutflow for the SR-Gtt-1L-M2 for a representative Gtt signal. Signal was generated with 30000 events. Expected yields are normalised to a luminosity of 139~fb$^{-1}$.

Cutflow for the SR-Gtt-1L-C for a representative Gtt signal. Signal was generated with 30000 events. Expected yields are normalised to a luminosity of 139~fb$^{-1}$.

Cutflow for the SR-Gbb-B for a representative Gbb signal. Signal was generated with 30000 events. Expected yields are normalised to a luminosity of 139~fb$^{-1}$.

Cutflow for the SR-Gbb-M for a representative Gbb signal. Signal was generated with 30000 events. Expected yields are normalised to a luminosity of 139~fb$^{-1}$.

Cutflow for the SR-Gbb-C for a representative Gbb signal. Signal was generated with 30000 events. Expected yields are normalised to a luminosity of 139~fb$^{-1}$.

Cutflow for the SR-Gtb-B for a representative Gtb signal. Signal was generated with 30000 events. Expected yields are normalised to a luminosity of 139~fb$^{-1}$.

Cutflow for the SR-Gtb-M for a representative Gtb signal. Signal was generated with 30000 events. Expected yields are normalised to a luminosity of 139~fb$^{-1}$.

Cutflow for the SR-Gtb-C for a representative Gtb signal. Signal was generated with 30000 events. Expected yields are normalised to a luminosity of 139~fb$^{-1}$.

Cutflow for the SR-Gtt-2100-1 for a representative Gtt signal. Signal was generated with 30000 events. Expected yields are normalised to a luminosity of 139~fb$^{-1}$.

Cutflow for the SR-Gtt-1800-1 for a representative Gtt signal. Signal was generated with 30000 events. Expected yields are normalised to a luminosity of 139~fb$^{-1}$.

Cutflow for the SR-Gtt-2300-1200 for a representative Gtt signal. Signal was generated with 30000 events. Expected yields are normalised to a luminosity of 139~fb$^{-1}$.

Cutflow for the SR-Gtt-1900-1400 for a representative Gtt signal. Signal was generated with 30000 events. Expected yields are normalised to a luminosity of 139~fb$^{-1}$.

Cutflow for the SR-Gbb-2800-1400 for a representative Gbb signal. Signal was generated with 30000 events. Expected yields are normalised to a luminosity of 139~fb$^{-1}$.

Cutflow for the SR-Gbb-2300-1000 for a representative Gbb signal. Signal was generated with 30000 events. Expected yields are normalised to a luminosity of 139~fb$^{-1}$.

Cutflow for the SR-Gbb-2100-1600 for a representative Gbb signal. Signal was generated with 30000 events. Expected yields are normalised to a luminosity of 139~fb$^{-1}$.

Cutflow for the SR-Gbb-2000-1800 for a representative Gbb signal. Signal was generated with 30000 events. Expected yields are normalised to a luminosity of 139~fb$^{-1}$.

Distribution of $m(e^{\pm},\mu^{\pm})_{\mathrm{lead}}$ in the electron-muon signal region after the background-only fit.

Distribution of $m(\mu^{\pm},\mu^{\pm})_{\mathrm{lead}}$ in the muon-muon signal region after the background-only fit.

Distribution of $m(\ell^{\pm},\ell^{\prime \pm})_{\mathrm{lead}}$ in the three-lepton signal region after the background-only fit.

Event yield in the four-lepton signal region after the background-only fit.

Two-lepton signal region weighted cutflows for electron-electron channel. Data and a few representative signal samples are presented alongside the main background contributions. Yield stands for the observed yield in data and the expected yields based on $\sigma \cdot \mathcal{L}$ for the signal and background samples. Top quark and multiboson process backgrounds are merged under the 'other' category.

Two-lepton signal region weighted cutflows for electron-muon channel. Data and a few representative signal samples are presented alongside the main background contributions. Yield stands for the observed yield in data and the expected yields based on $\sigma \cdot \mathcal{L}$ for the signal and background samples. Top quark and multiboson process backgrounds are merged under the 'other' category.

Two-lepton signal region weighted cutflows for muon-muon channel. Data and a few representative signal samples are presented alongside the main background contributions. Yield stands for the observed yield in data and the expected yields based on $\sigma \cdot \mathcal{L}$ for the signal and background samples. Top quark and multiboson process backgrounds are merged under the 'other' category.

Three-lepton signal region weighted cutflow for $\ell \ell \ell$ channel. Data and a few representative signal samples are presented alongside the main background contributions. Yield stands for the observed yield in data and the expected yields based on $\sigma \cdot \mathcal{L}$ for the signal and background samples. Top quark and multiboson process backgrounds are merged under the 'other' category.

Four-lepton signal region weighted cutflow for $\ell \ell \ell \ell$ channel. Data and a few representative signal samples are presented alongside the main background contributions. Yield stands for the observed yield in data and the expected yields based on $\sigma \cdot \mathcal{L}$ for the signal and background samples. Top quark, Drell-Yan, and multiboson process backgrounds are merged under the 'other' category.

Probability of finding at least one TAR jet, where the p_T-leading TAR jet passes the m_Wcand and D₂^&beta;=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&#37; C.L. for the dark Higgs model in the (m_Z', m_s) plane for g_q=0.25, g_&chi;=1, m_&chi;=200 GeV, and sin&theta;=0.01.

Expected exclusion contour at 95&#37; C.L. for the dark Higgs model in the (m_Z', m_s) plane for g_q=0.25, g_&chi;=1, m_&chi;=200 GeV, and sin&theta;=0.01.

Expected+1&sigma; exclusion contour at 95&#37; C.L. for the dark Higgs model in the (m_Z', m_s) plane for g_q=0.25, g_&chi;=1, m_&chi;=200 GeV, and sin&theta;=0.01.

Expected-1&sigma; exclusion contour at 95&#37; C.L. for the dark Higgs model in the (m_Z', m_s) plane for g_q=0.25, g_&chi;=1, m_&chi;=200 GeV, and sin&theta;=0.01.

Expected+2&sigma; exclusion contour at 95&#37; C.L. for the dark Higgs model in the (m_Z', m_s) plane for g_q=0.25, g_&chi;=1, m_&chi;=200 GeV, and sin&theta;=0.01.

Expected-2&sigma; exclusion contour at 95&#37; C.L. for the dark Higgs model in the (m_Z', m_s) plane for g_q=0.25, g_&chi;=1, m_&chi;=200 GeV, and sin&theta;=0.01.

Observed upper limits at 95&#37; C.L. on &sigma;(pp &rarr; s &chi;&chi;) &times; B(s &rarr; W^&pm; W^&mp;) 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 &sigma;(pp &rarr; s &chi;&chi;) &times; B(s &rarr; W^&pm; W^&mp;) process for m_Z'=0.5 TeV, shown in dashed blue.

Observed upper limits at 95&#37; C.L. on &sigma;(pp &rarr; s &chi;&chi;) &times; B(s &rarr; W^&pm; W^&mp;) 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 &sigma;(pp &rarr; s &chi;&chi;) &times; B(s &rarr; W^&pm; W^&mp;) process for m_Z'=1 TeV, shown in dashed blue.

Observed upper limits at 95&#37; C.L. on &sigma;(pp &rarr; s &chi;&chi;) &times; B(s &rarr; W^&pm; W^&mp;) 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 &sigma;(pp &rarr; s &chi;&chi;) &times; B(s &rarr; W^&pm; W^&mp;) process for m_Z'=1.7 TeV, shown in dashed blue.

Observed upper limits at 95&#37; C.L. on &sigma;(pp &rarr; s &chi;&chi;) &times; B(s &rarr; W^&pm; W^&mp;) 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 &sigma;(pp &rarr; s &chi;&chi;) &times; B(s &rarr; W^&pm; W^&mp;) 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_&chi;=1.0, sin&theta;=0.01 and m_&chi;=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_&chi;=1.0, sin&theta;=0.01, m_Z' = 1 TeV, and m_&chi;=200 GeV considering s&rarr;W(&ell;&nu;)W(qq) decays only.

Cumulative efficiencies in the resolved category for three representative dark Higgs signal scenarios with g_q=0.25, g_&chi;=1.0, sin&theta;=0.01, m_Z' = 1 TeV, and m_&chi;=200 GeV considering s&rarr;W(&ell;&nu;)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;).

Summary of the total uncertainty in the background prediction for each SR of the tt0L-low, tt0L-high, tt1L and tt2L analysis channels in the statistical combination. Their dominant contributions are indicated by individual lines. Individual uncertainties can be correlated, and do not necessarily add up in quadrature to the total background uncertainty.

Exclusion limits for colour-neutral scalar mediator dark matter models as a function of the mediator mass $m(\phi)$ for a DM mass $m_{\chi} = 1$ GeV. Associated production of DM with both single top quarks ($tW$ and $tj$ channels) and top quark pairs is considered. The limits are calculated at 95% CL and are expressed in terms of the ratio of the excluded cross section to the cross section for a coupling assumption of $g = g_q = g_{\chi} = 1$. The solid (dashed) lines show the observed (expected) exclusion limits for each individual channel and their statistical combination.

Exclusion limits for colour-neutral pseudoscalar mediator dark matter models as a function of the mediator mass $m(a)$ for a DM mass $m_{\chi} = 1$ GeV. Associated production of DM with both single top quarks ($tW$ and $tj$ channels) and top quark pairs is considered. The limits are calculated at 95% CL and are expressed in terms of the ratio of the excluded cross section to the cross section for a coupling assumption of $g = g_q = g_{\chi} = 1$. The solid (dashed) lines show the observed (expected) exclusion limits for each individual channel and their statistical combination.

$E_{\text{T}}^{\text{miss}}$ distribution in SR0X for the tt0L-low analysis. The contributions from all SM backgrounds are shown after the profile likelihood simultaneous fit to all tt0L-low CRs, with the hatched bands representing the total uncertainty. The category '$t\bar{t}$ (other)' represents $t\bar{t}$ events without extra jets or events with extra light-flavour jets. 'Other' includes contributions from $t\bar{t}W$, $tZ$ and $tWZ$ processes. The expected distributions for selected signal models are shown as dashed lines. The overflow events are included in the last bin. The bottom panels show the ratio of the observed data to the total SM background prediction, with the hatched area representing the total uncertainty in the background prediction and the red arrows marking data outside the vertical-axis range.

$E_{\text{T}}^{\text{miss}}$ distribution in SRWX for the tt0L-low analysis. The contributions from all SM backgrounds are shown after the profile likelihood simultaneous fit to all tt0L-low CRs, with the hatched bands representing the total uncertainty. The category '$t\bar{t}$ (other)' represents $t\bar{t}$ events without extra jets or events with extra light-flavour jets. 'Other' includes contributions from $t\bar{t}W$, $tZ$ and $tWZ$ processes. The expected distributions for selected signal models are shown as dashed lines. The overflow events are included in the last bin. The bottom panels show the ratio of the observed data to the total SM background prediction, with the hatched area representing the total uncertainty in the background prediction and the red arrows marking data outside the vertical-axis range.

$E_{\text{T}}^{\text{miss}}$ distribution in SRTX for the tt0L-low analysis. The contributions from all SM backgrounds are shown after the profile likelihood simultaneous fit to all tt0L-low CRs, with the hatched bands representing the total uncertainty. The category '$t\bar{t}$ (other)' represents $t\bar{t}$ events without extra jets or events with extra light-flavour jets. 'Other' includes contributions from $t\bar{t}W$, $tZ$ and $tWZ$ processes. The expected distributions for selected signal models are shown as dashed lines. The overflow events are included in the last bin. The bottom panels show the ratio of the observed data to the total SM background prediction, with the hatched area representing the total uncertainty in the background prediction and the red arrows marking data outside the vertical-axis range.

Exclusion limits for colour-neutral scalar mediator dark matter models as a function of the mediator mass $m(\phi)$ for a DM mass $m_{\chi} = 1$ GeV. Associated production of DM with both single top quarks ($tW$ and $tj$ channels) and top quark pairs is considered. The limits are calculated at 95% CL and are expressed in terms of the ratio of the excluded cross section to the nominal cross section for a coupling assumption of $g = g_q = g_{\chi} = 1$. The solid (dashed) lines show the observed (expected) exclusion limits for the tt0L-high and tt0L-low analyses and their statistical combination.

Exclusion limits for colour-neutral pseudoscalar mediator dark matter models as a function of the mediator mass $m(a)$ for a DM mass $m_{\chi} = 1$ GeV. Associated production of DM with both single top quarks ($tW$ and $tj$ channels) and top quark pairs is considered. The limits are calculated at 95% CL and are expressed in terms of the ratio of the excluded cross section to the nominal cross section for a coupling assumption of $g = g_q = g_{\chi} = 1$. The solid (dashed) lines show the observed (expected) exclusion limits for the tt0L-high and tt0L-low analyses and their statistical combination.

Exclusion limits for colour-neutral scalar mediator dark matter models as a function of the mediator mass $m(\phi)$ for a DM mass $m_{\chi} = 1$ GeV. Only associated production of DM with top quark pairs is considered for this interpretation. The limits are calculated at 95% CL and are expressed in terms of the ratio of the excluded cross section to the cross section for a coupling assumption of $g = g_q = g_{\chi} = 1$. The solid (dashed) lines show the observed (expected) exclusion limits for each individual channel and their statistical combination.

Exclusion limits for colour-neutral pseudoscalar mediator dark matter models as a function of the mediator mass $m(a)$ for a DM mass $m_{\chi} = 1$ GeV. Only associated production of DM with top quark pairs is considered for this interpretation. The limits are calculated at 95% CL and are expressed in terms of the ratio of the excluded cross section to the cross section for a coupling assumption of $g = g_q = g_{\chi} = 1$. The solid (dashed) lines show the observed (expected) exclusion limits for each individual channel and their statistical combination.

Exclusion limits for colour-neutral scalar mediator dark matter models as a function of the mediator mass $m(\phi)$ for a DM mass $m_{\chi} = 1$ GeV. Only associated production of DM with top quark pairs is considered for this interpretation. The limits are calculated at 95% CL and are expressed in terms of the ratio of the excluded cross section to the nominal cross section for a coupling assumption of $g = g_q = g_{\chi} = 1$. The solid (dashed) lines show the observed (expected) exclusion limits for the tt0L-high and tt0L-low analyses and their statistical combination.

Exclusion limits for colour-neutral pseudoscalar mediator dark matter models as a function of the mediator mass $m(a)$ for a DM mass $m_{\chi} = 1$ GeV. Only associated production of DM with top quark pairs is considered for this interpretation. The limits are calculated at 95% CL and are expressed in terms of the ratio of the excluded cross section to the nominal cross section for a coupling assumption of $g = g_q = g_{\chi} = 1$. The solid (dashed) lines show the observed (expected) exclusion limits for the tt0L-high and tt0L-low analyses and their statistical combination.

Representative fit distribution in the different flavour leptons signal region for the tt2L analysis: each bin of such distribution, starting from the red arrow, corresponds to a single SR included in the fit. 'FNP' includes the contribution from fake/non-prompt lepton background arising from jets (mainly $\pi/K$, heavy-flavour hadron decays and photon conversion) misidentified as leptons, estimated in a purely data-driven way. 'Other' includes contributions from $t\bar{t}W$, $tZ$ and $tWZ$ processes. The total uncertainty in the SM expectation is represented with hatched bands and the expected distributions for selected signal models are shown as dashed lines.

Signal acceptance in SR0X, SRWX and SRTX for simplified DM+$t\bar{t}$ model, defined as the number of accepted events at generator level in signal Monte Carlo simulation divided by the total number of events in the sample.

Signal acceptance in SR0X, SRWX and SRTX for simplified DM+$tW$ model, defined as the number of accepted events at generator level in signal Monte Carlo simulation divided by the total number of events in the sample.

Signal acceptance in SR0X, SRWX and SRTX for simplified DM+$tj$ model, defined as the number of accepted events at generator level in signal Monte Carlo simulation divided by the total number of events in the sample.

Signal efficiency in SR0X, SRWX and SRTX for simplified DM+$t\bar{t}$ model, defined as the number of selected reconstructed events divided by the acceptance.

Signal efficiency in SR0X, SRWX and SRTX for simplified DM+$tW$ model, defined as the number of selected reconstructed events divided by the acceptance.

Signal efficiency in SR0X, SRWX and SRTX for simplified DM+$tj$ model, defined as the number of selected reconstructed events divided by the acceptance.

Cutflow for the reference point DM+$t\bar{t}$ $m(\phi, \chi) = (10, 1)$ GeV in signal region SR0X. The column labelled 'weighted' shows the event yield including all correction factors applied to simulation, and is normalised to 139 fb$^{-1}$. A notable exception concerns the 'weighted' numbers in the first and the second row, labelled 'Total' and 'Filtered', which correspond to $\mathcal{L}\cdot\sigma$ and $\mathcal{L}\cdot\sigma\cdot\epsilon$ expected, respectively. The 'Skim' selection requires the $p_{\text{T}}$ of the leading four jets to be above (80, 60, 40, 40) GeV, the missing transverse momentum $E_{\text{T}}^{\text{miss}} > 140$ GeV, the missing momentum significance $\mathcal{S} > 8$, $\Delta\phi_{\min}(\vec{p}_{\text{T,1-4}},\vec{p}_{\text{T}}^{\text{miss}}) > 0.4$ and a lepton veto. The 'Orthogonalisation' selection is defined in the main body. In total 2045000 raw MC events were generated prior to the specified cuts, with the column 'Unweighted yield' collecting the numbers after each cut.

Cutflow for the reference point DM+$t\bar{t}$ $m(\phi, \chi) = (10, 1)$ GeV in signal region SRWX. The column labelled 'weighted' shows the event yield including all correction factors applied to simulation, and is normalised to 139 fb$^{-1}$. A notable exception concerns the 'weighted' numbers in the first and the second row, labelled 'Total' and 'Filtered', which correspond to $\mathcal{L}\cdot\sigma$ and $\mathcal{L}\cdot\sigma\cdot\epsilon$ expected, respectively. The 'Skim' selection requires the $p_{\text{T}}$ of the leading four jets to be above (80, 60, 40, 40) GeV, the missing transverse momentum $E_{\text{T}}^{\text{miss}} > 140$ GeV, the missing momentum significance $\mathcal{S} > 8$, $\Delta\phi_{\min}(\vec{p}_{\text{T,1-4}},\vec{p}_{\text{T}}^{\text{miss}}) > 0.4$ and a lepton veto. The 'Orthogonalisation' selection is defined in the main body. In total 2045000 raw MC events were generated prior to the specified cuts, with the column 'Unweighted yield' collecting the numbers after each cut.

Cutflow for the reference point DM+$t\bar{t}$ $m(\phi, \chi) = (10, 1)$ GeV in signal region SRTX. The column labelled 'weighted' shows the event yield including all correction factors applied to simulation, and is normalised to 139 fb$^{-1}$. A notable exception concerns the 'weighted' numbers in the first and the second row, labelled 'Total' and 'Filtered', which correspond to $\mathcal{L}\cdot\sigma$ and $\mathcal{L}\cdot\sigma\cdot\epsilon$ expected, respectively. The 'Skim' selection requires the $p_{\text{T}}$ of the leading four jets to be above (80, 60, 40, 40) GeV, the missing transverse momentum $E_{\text{T}}^{\text{miss}} > 140$ GeV, the missing momentum significance $\mathcal{S} > 8$, $\Delta\phi_{\min}(\vec{p}_{\text{T,1-4}},\vec{p}_{\text{T}}^{\text{miss}}) > 0.4$ and a lepton veto. The 'Orthogonalisation' selection is defined in the main body. In total 2045000 raw MC events were generated prior to the specified cuts, with the column 'Unweighted yield' collecting the numbers after each cut.

Cutflow for the reference point DM+$t\bar{t}$ $m(a, \chi) = (10, 1)$ GeV in signal region SR0X. The column labelled 'weighted' shows the event yield including all correction factors applied to simulation, and is normalised to 139 fb$^{-1}$. A notable exception concerns the 'weighted' numbers in the first and the second row, labelled 'Total' and 'Filtered', which correspond to $\mathcal{L}\cdot\sigma$ and $\mathcal{L}\cdot\sigma\cdot\epsilon$ expected, respectively. The 'Skim' selection requires the $p_{\text{T}}$ of the leading four jets to be above (80, 60, 40, 40) GeV, the missing transverse momentum $E_{\text{T}}^{\text{miss}} > 140$ GeV, the missing momentum significance $\mathcal{S} > 8$, $\Delta\phi_{\min}(\vec{p}_{\text{T,1-4}},\vec{p}_{\text{T}}^{\text{miss}}) > 0.4$ and a lepton veto. The 'Orthogonalisation' selection is defined in the main body. In total 400000 raw MC events were generated prior to the specified cuts, with the column 'Unweighted yield' collecting the numbers after each cut.

Cutflow for the reference point DM+$t\bar{t}$ $m(a, \chi) = (10, 1)$ GeV in signal region SRWX. The column labelled 'weighted' shows the event yield including all correction factors applied to simulation, and is normalised to 139 fb$^{-1}$. A notable exception concerns the 'weighted' numbers in the first and the second row, labelled 'Total' and 'Filtered', which correspond to $\mathcal{L}\cdot\sigma$ and $\mathcal{L}\cdot\sigma\cdot\epsilon$ expected, respectively. The 'Skim' selection requires the $p_{\text{T}}$ of the leading four jets to be above (80, 60, 40, 40) GeV, the missing transverse momentum $E_{\text{T}}^{\text{miss}} > 140$ GeV, the missing momentum significance $\mathcal{S} > 8$, $\Delta\phi_{\min}(\vec{p}_{\text{T,1-4}},\vec{p}_{\text{T}}^{\text{miss}}) > 0.4$ and a lepton veto. The 'Orthogonalisation' selection is defined in the main body. In total 400000 raw MC events were generated prior to the specified cuts, with the column 'Unweighted yield' collecting the numbers after each cut.

Cutflow for the reference point DM+$t\bar{t}$ $m(a, \chi) = (10, 1)$ GeV in signal region SRTX. The column labelled 'weighted' shows the event yield including all correction factors applied to simulation, and is normalised to 139 fb$^{-1}$. A notable exception concerns the 'weighted' numbers in the first and the second row, labelled 'Total' and 'Filtered', which correspond to $\mathcal{L}\cdot\sigma$ and $\mathcal{L}\cdot\sigma\cdot\epsilon$ expected, respectively. The 'Skim' selection requires the $p_{\text{T}}$ of the leading four jets to be above (80, 60, 40, 40) GeV, the missing transverse momentum $E_{\text{T}}^{\text{miss}} > 140$ GeV, the missing momentum significance $\mathcal{S} > 8$, $\Delta\phi_{\min}(\vec{p}_{\text{T,1-4}},\vec{p}_{\text{T}}^{\text{miss}}) > 0.4$ and a lepton veto. The 'Orthogonalisation' selection is defined in the main body. In total 400000 raw MC events were generated prior to the specified cuts, with the column 'Unweighted yield' collecting the numbers after each cut.

Cutflow for the reference point DM+$tW$ $m(\phi, \chi) = (10, 1)$ GeV in signal region SR0X. The column labelled 'weighted' shows the event yield including all correction factors applied to simulation, and is normalised to 139 fb$^{-1}$. A notable exception concerns the 'weighted' numbers in the first and the second row, labelled 'Total' and 'Filtered', which correspond to $\mathcal{L}\cdot\sigma$ and $\mathcal{L}\cdot\sigma\cdot\epsilon$ expected, respectively. The 'Skim' selection requires the $p_{\text{T}}$ of the leading four jets to be above (80, 60, 40, 40) GeV, the missing transverse momentum $E_{\text{T}}^{\text{miss}} > 140$ GeV, the missing momentum significance $\mathcal{S} > 8$, $\Delta\phi_{\min}(\vec{p}_{\text{T,1-4}},\vec{p}_{\text{T}}^{\text{miss}}) > 0.4$ and a lepton veto. The 'Orthogonalisation' selection is defined in the main body. In total 120000 raw MC events were generated prior to the specified cuts, with the column 'Unweighted yield' collecting the numbers after each cut.

Cutflow for the reference point DM+$tW$ $m(\phi, \chi) = (10, 1)$ GeV in signal region SRWX. The column labelled 'weighted' shows the event yield including all correction factors applied to simulation, and is normalised to 139 fb$^{-1}$. A notable exception concerns the 'weighted' numbers in the first and the second row, labelled 'Total' and 'Filtered', which correspond to $\mathcal{L}\cdot\sigma$ and $\mathcal{L}\cdot\sigma\cdot\epsilon$ expected, respectively. The 'Skim' selection requires the $p_{\text{T}}$ of the leading four jets to be above (80, 60, 40, 40) GeV, the missing transverse momentum $E_{\text{T}}^{\text{miss}} > 140$ GeV, the missing momentum significance $\mathcal{S} > 8$, $\Delta\phi_{\min}(\vec{p}_{\text{T,1-4}},\vec{p}_{\text{T}}^{\text{miss}}) > 0.4$ and a lepton veto. The 'Orthogonalisation' selection is defined in the main body. In total 120000 raw MC events were generated prior to the specified cuts, with the column 'Unweighted yield' collecting the numbers after each cut.

Cutflow for the reference point DM+$tW$ $m(\phi, \chi) = (10, 1)$ GeV in signal region SRTX. The column labelled 'weighted' shows the event yield including all correction factors applied to simulation, and is normalised to 139 fb$^{-1}$. A notable exception concerns the 'weighted' numbers in the first and the second row, labelled 'Total' and 'Filtered', which correspond to $\mathcal{L}\cdot\sigma$ and $\mathcal{L}\cdot\sigma\cdot\epsilon$ expected, respectively. The 'Skim' selection requires the $p_{\text{T}}$ of the leading four jets to be above (80, 60, 40, 40) GeV, the missing transverse momentum $E_{\text{T}}^{\text{miss}} > 140$ GeV, the missing momentum significance $\mathcal{S} > 8$, $\Delta\phi_{\min}(\vec{p}_{\text{T,1-4}},\vec{p}_{\text{T}}^{\text{miss}}) > 0.4$ and a lepton veto. The 'Orthogonalisation' selection is defined in the main body. In total 120000 raw MC events were generated prior to the specified cuts, with the column 'Unweighted yield' collecting the numbers after each cut.

Cutflow for the reference point DM+$tW$ $m(a, \chi) = (10, 1)$ GeV in signal region SR0X. The column labelled 'weighted' shows the event yield including all correction factors applied to simulation, and is normalised to 139 fb$^{-1}$. A notable exception concerns the 'weighted' numbers in the first and the second row, labelled 'Total' and 'Filtered', which correspond to $\mathcal{L}\cdot\sigma$ and $\mathcal{L}\cdot\sigma\cdot\epsilon$ expected, respectively. The 'Skim' selection requires the $p_{\text{T}}$ of the leading four jets to be above (80, 60, 40, 40) GeV, the missing transverse momentum $E_{\text{T}}^{\text{miss}} > 140$ GeV, the missing momentum significance $\mathcal{S} > 8$, $\Delta\phi_{\min}(\vec{p}_{\text{T,1-4}},\vec{p}_{\text{T}}^{\text{miss}}) > 0.4$ and a lepton veto. The 'Orthogonalisation' selection is defined in the main body. In total 100000 raw MC events were generated prior to the specified cuts, with the column 'Unweighted yield' collecting the numbers after each cut.

Cutflow for the reference point DM+$tW$ $m(a, \chi) = (10, 1)$ GeV in signal region SRWX. The column labelled 'weighted' shows the event yield including all correction factors applied to simulation, and is normalised to 139 fb$^{-1}$. A notable exception concerns the 'weighted' numbers in the first and the second row, labelled 'Total' and 'Filtered', which correspond to $\mathcal{L}\cdot\sigma$ and $\mathcal{L}\cdot\sigma\cdot\epsilon$ expected, respectively. The 'Skim' selection requires the $p_{\text{T}}$ of the leading four jets to be above (80, 60, 40, 40) GeV, the missing transverse momentum $E_{\text{T}}^{\text{miss}} > 140$ GeV, the missing momentum significance $\mathcal{S} > 8$, $\Delta\phi_{\min}(\vec{p}_{\text{T,1-4}},\vec{p}_{\text{T}}^{\text{miss}}) > 0.4$ and a lepton veto. The 'Orthogonalisation' selection is defined in the main body. In total 100000 raw MC events were generated prior to the specified cuts, with the column 'Unweighted yield' collecting the numbers after each cut.

Cutflow for the reference point DM+$tW$ $m(a, \chi) = (10, 1)$ GeV in signal region SRTX. The column labelled 'weighted' shows the event yield including all correction factors applied to simulation, and is normalised to 139 fb$^{-1}$. A notable exception concerns the 'weighted' numbers in the first and the second row, labelled 'Total' and 'Filtered', which correspond to $\mathcal{L}\cdot\sigma$ and $\mathcal{L}\cdot\sigma\cdot\epsilon$ expected, respectively. The 'Skim' selection requires the $p_{\text{T}}$ of the leading four jets to be above (80, 60, 40, 40) GeV, the missing transverse momentum $E_{\text{T}}^{\text{miss}} > 140$ GeV, the missing momentum significance $\mathcal{S} > 8$, $\Delta\phi_{\min}(\vec{p}_{\text{T,1-4}},\vec{p}_{\text{T}}^{\text{miss}}) > 0.4$ and a lepton veto. The 'Orthogonalisation' selection is defined in the main body. In total 100000 raw MC events were generated prior to the specified cuts, with the column 'Unweighted yield' collecting the numbers after each cut.

Cutflow for the reference point DM+$tj$ $m(\phi, \chi) = (10, 1)$ GeV in signal region SR0X. The column labelled 'weighted' shows the event yield including all correction factors applied to simulation, and is normalised to 139 fb$^{-1}$. A notable exception concerns the 'weighted' numbers in the first and the second row, labelled 'Total' and 'Filtered', which correspond to $\mathcal{L}\cdot\sigma$ and $\mathcal{L}\cdot\sigma\cdot\epsilon$ expected, respectively. The 'Skim' selection requires the $p_{\text{T}}$ of the leading four jets to be above (80, 60, 40, 40) GeV, the missing transverse momentum $E_{\text{T}}^{\text{miss}} > 140$ GeV, the missing momentum significance $\mathcal{S} > 8$, $\Delta\phi_{\min}(\vec{p}_{\text{T,1-4}},\vec{p}_{\text{T}}^{\text{miss}}) > 0.4$ and a lepton veto. The 'Orthogonalisation' selection is defined in the main body. In total 169000 raw MC events were generated prior to the specified cuts, with the column 'Unweighted yield' collecting the numbers after each cut.

Cutflow for the reference point DM+$tj$ $m(\phi, \chi) = (10, 1)$ GeV in signal region SRWX. The column labelled 'weighted' shows the event yield including all correction factors applied to simulation, and is normalised to 139 fb$^{-1}$. A notable exception concerns the 'weighted' numbers in the first and the second row, labelled 'Total' and 'Filtered', which correspond to $\mathcal{L}\cdot\sigma$ and $\mathcal{L}\cdot\sigma\cdot\epsilon$ expected, respectively. The 'Skim' selection requires the $p_{\text{T}}$ of the leading four jets to be above (80, 60, 40, 40) GeV, the missing transverse momentum $E_{\text{T}}^{\text{miss}} > 140$ GeV, the missing momentum significance $\mathcal{S} > 8$, $\Delta\phi_{\min}(\vec{p}_{\text{T,1-4}},\vec{p}_{\text{T}}^{\text{miss}}) > 0.4$ and a lepton veto. The 'Orthogonalisation' selection is defined in the main body. In total 169000 raw MC events were generated prior to the specified cuts, with the column 'Unweighted yield' collecting the numbers after each cut.

Cutflow for the reference point DM+$tj$ $m(\phi, \chi) = (10, 1)$ GeV in signal region SRTX. The column labelled 'weighted' shows the event yield including all correction factors applied to simulation, and is normalised to 139 fb$^{-1}$. A notable exception concerns the 'weighted' numbers in the first and the second row, labelled 'Total' and 'Filtered', which correspond to $\mathcal{L}\cdot\sigma$ and $\mathcal{L}\cdot\sigma\cdot\epsilon$ expected, respectively. The 'Skim' selection requires the $p_{\text{T}}$ of the leading four jets to be above (80, 60, 40, 40) GeV, the missing transverse momentum $E_{\text{T}}^{\text{miss}} > 140$ GeV, the missing momentum significance $\mathcal{S} > 8$, $\Delta\phi_{\min}(\vec{p}_{\text{T,1-4}},\vec{p}_{\text{T}}^{\text{miss}}) > 0.4$ and a lepton veto. The 'Orthogonalisation' selection is defined in the main body. In total 169000 raw MC events were generated prior to the specified cuts, with the column 'Unweighted yield' collecting the numbers after each cut.

Cutflow for the reference point DM+$tj$ $m(a, \chi) = (10, 1)$ GeV in signal region SR0X. The column labelled 'weighted' shows the event yield including all correction factors applied to simulation, and is normalised to 139 fb$^{-1}$. A notable exception concerns the 'weighted' numbers in the first and the second row, labelled 'Total' and 'Filtered', which correspond to $\mathcal{L}\cdot\sigma$ and $\mathcal{L}\cdot\sigma\cdot\epsilon$ expected, respectively. The 'Skim' selection requires the $p_{\text{T}}$ of the leading four jets to be above (80, 60, 40, 40) GeV, the missing transverse momentum $E_{\text{T}}^{\text{miss}} > 140$ GeV, the missing momentum significance $\mathcal{S} > 8$, $\Delta\phi_{\min}(\vec{p}_{\text{T,1-4}},\vec{p}_{\text{T}}^{\text{miss}}) > 0.4$ and a lepton veto. The 'Orthogonalisation' selection is defined in the main body. In total 140000 raw MC events were generated prior to the specified cuts, with the column 'Unweighted yield' collecting the numbers after each cut.

Cutflow for the reference point DM+$tj$ $m(a, \chi) = (10, 1)$ GeV in signal region SRWX. The column labelled 'weighted' shows the event yield including all correction factors applied to simulation, and is normalised to 139 fb$^{-1}$. A notable exception concerns the 'weighted' numbers in the first and the second row, labelled 'Total' and 'Filtered', which correspond to $\mathcal{L}\cdot\sigma$ and $\mathcal{L}\cdot\sigma\cdot\epsilon$ expected, respectively. The 'Skim' selection requires the $p_{\text{T}}$ of the leading four jets to be above (80, 60, 40, 40) GeV, the missing transverse momentum $E_{\text{T}}^{\text{miss}} > 140$ GeV, the missing momentum significance $\mathcal{S} > 8$, $\Delta\phi_{\min}(\vec{p}_{\text{T,1-4}},\vec{p}_{\text{T}}^{\text{miss}}) > 0.4$ and a lepton veto. The 'Orthogonalisation' selection is defined in the main body. In total 140000 raw MC events were generated prior to the specified cuts, with the column 'Unweighted yield' collecting the numbers after each cut.

Cutflow for the reference point DM+$tj$ $m(a, \chi) = (10, 1)$ GeV in signal region SRTX. The column labelled 'weighted' shows the event yield including all correction factors applied to simulation, and is normalised to 139 fb$^{-1}$. A notable exception concerns the 'weighted' numbers in the first and the second row, labelled 'Total' and 'Filtered', which correspond to $\mathcal{L}\cdot\sigma$ and $\mathcal{L}\cdot\sigma\cdot\epsilon$ expected, respectively. The 'Skim' selection requires the $p_{\text{T}}$ of the leading four jets to be above (80, 60, 40, 40) GeV, the missing transverse momentum $E_{\text{T}}^{\text{miss}} > 140$ GeV, the missing momentum significance $\mathcal{S} > 8$, $\Delta\phi_{\min}(\vec{p}_{\text{T,1-4}},\vec{p}_{\text{T}}^{\text{miss}}) > 0.4$ and a lepton veto. The 'Orthogonalisation' selection is defined in the main body. In total 140000 raw MC events were generated prior to the specified cuts, with the column 'Unweighted yield' collecting the numbers after each cut.