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}$

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.

Parametrization of the $A_{X}$ factor, defined as the fraction of diphoton resonances satisfying the fiducial acceptance at the particle level, as function of $m_{X}$ obtained from the narrow width simulated signal samples produced in gluon fusion.

The correction factor, $C_{X}$, defined as the ratio between the number of reconstructed signal events passing the analysis cuts and the number of signal events at the particle level generated within the fiducial volume, and acceptance correction factor, $A_{X}$, defined as the fraction of diphoton resonances satisfying the fiducial acceptance at the particle level. Both are computed for NWA spin-0 models as a function of $m_{X}$.

Effect of event selections on a scalar MC signal sample generated for $m_{X}$ = 15 GeV and on the data. For the MC sample, the efficiencies are shown after applying event weights and a truth level filter that requires two photons with $p^{\gamma\gamma}_{T}>40$ GeV; for the data, the absolute yields are shown. The initial yields for data include a trigger preselection that is the OR of a list of single photon and diphoton triggers. The "2 $loose$ photons" step includes the kinematic acceptance cuts.

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

Observed and expected 95% CL upper limits on the product of the production cross section ($\sigma$) and the branching fraction, obtained after combining all categories with 138 $\mathrm{fb}^{−1}$ of data at $\sqrt{s}$ = 13 TeV for $\mathrm{V'}$ to $VV$ + $VH$ signal in HVT model C.

Total signal efficiency including all event categories as a function of the resonance mass $M_X$ for all DY/ggF signal models under consideration. The estimates from MC simulation (markers) are interpolated with a functional form (lines).

Total signal efficiency including all event categories as a function of the resonance mass $M_X$ for all VBF signal models under consideration. The estimates from MC simulation (markers) are interpolated with a functional form (lines).

The figure shows the signal acceptance (a) and efficiency (b) plots for the slepton pair production model, in the SR-0J $m_{\mathrm{T2}}^{100} \in[110,\infty)$ region. Acceptance is calculated by applying the signal region requirements to particle-level objects, which do not suffer from identification inefficiencies or mismeasurements. The efficiency is calculated with fully reconstructed objects with the acceptance divided out. Large acceptance and efficiency differences in neighbouring points are due to statistical fluctuations.

The figure shows the signal acceptance (a) and efficiency (b) plots for the slepton pair production model, in the SR-0J $m_{\mathrm{T2}}^{100} \in[110,\infty)$ region. Acceptance is calculated by applying the signal region requirements to particle-level objects, which do not suffer from identification inefficiencies or mismeasurements. The efficiency is calculated with fully reconstructed objects with the acceptance divided out. Large acceptance and efficiency differences in neighbouring points are due to statistical fluctuations.

The figure shows the signal acceptance (a) and efficiency (b) plots for the slepton pair production model, in the SR-0J $m_{\mathrm{T2}}^{100} \in[120,\infty)$ region. Acceptance is calculated by applying the signal region requirements to particle-level objects, which do not suffer from identification inefficiencies or mismeasurements. The efficiency is calculated with fully reconstructed objects with the acceptance divided out. Large acceptance and efficiency differences in neighbouring points are due to statistical fluctuations.

The figure shows the signal acceptance (a) and efficiency (b) plots for the slepton pair production model, in the SR-0J $m_{\mathrm{T2}}^{100} \in[120,\infty)$ region. Acceptance is calculated by applying the signal region requirements to particle-level objects, which do not suffer from identification inefficiencies or mismeasurements. The efficiency is calculated with fully reconstructed objects with the acceptance divided out. Large acceptance and efficiency differences in neighbouring points are due to statistical fluctuations.

The figure shows the signal acceptance (a) and efficiency (b) plots for the slepton pair production model, in the SR-0J $m_{\mathrm{T2}}^{100} \in[130,\infty)$ region. Acceptance is calculated by applying the signal region requirements to particle-level objects, which do not suffer from identification inefficiencies or mismeasurements. The efficiency is calculated with fully reconstructed objects with the acceptance divided out. Large acceptance and efficiency differences in neighbouring points are due to statistical fluctuations.

The figure shows the signal acceptance (a) and efficiency (b) plots for the slepton pair production model, in the SR-0J $m_{\mathrm{T2}}^{100} \in[130,\infty)$ region. Acceptance is calculated by applying the signal region requirements to particle-level objects, which do not suffer from identification inefficiencies or mismeasurements. The efficiency is calculated with fully reconstructed objects with the acceptance divided out. Large acceptance and efficiency differences in neighbouring points are due to statistical fluctuations.

The figure shows the signal acceptance (a) and efficiency (b) plots for the slepton pair production model, in the SR-0J $m_{\mathrm{T2}}^{100} \in[100,105)$ region. Acceptance is calculated by applying the signal region requirements to particle-level objects, which do not suffer from identification inefficiencies or mismeasurements. The efficiency is calculated with fully reconstructed objects with the acceptance divided out. Large acceptance and efficiency differences in neighbouring points are due to statistical fluctuations.

The figure shows the signal acceptance (a) and efficiency (b) plots for the slepton pair production model, in the SR-0J $m_{\mathrm{T2}}^{100} \in[100,105)$ region. Acceptance is calculated by applying the signal region requirements to particle-level objects, which do not suffer from identification inefficiencies or mismeasurements. The efficiency is calculated with fully reconstructed objects with the acceptance divided out. Large acceptance and efficiency differences in neighbouring points are due to statistical fluctuations.

The figure shows the signal acceptance (a) and efficiency (b) plots for the slepton pair production model, in the SR-0J $m_{\mathrm{T2}}^{100} \in[105,110)$ region. Acceptance is calculated by applying the signal region requirements to particle-level objects, which do not suffer from identification inefficiencies or mismeasurements. The efficiency is calculated with fully reconstructed objects with the acceptance divided out. Large acceptance and efficiency differences in neighbouring points are due to statistical fluctuations.

The figure shows the signal acceptance (a) and efficiency (b) plots for the slepton pair production model, in the SR-0J $m_{\mathrm{T2}}^{100} \in[105,110)$ region. Acceptance is calculated by applying the signal region requirements to particle-level objects, which do not suffer from identification inefficiencies or mismeasurements. The efficiency is calculated with fully reconstructed objects with the acceptance divided out. Large acceptance and efficiency differences in neighbouring points are due to statistical fluctuations.

The figure shows the signal acceptance (a) and efficiency (b) plots for the slepton pair production model, in the SR-0J $m_{\mathrm{T2}}^{100} \in[110,115)$ region. Acceptance is calculated by applying the signal region requirements to particle-level objects, which do not suffer from identification inefficiencies or mismeasurements. The efficiency is calculated with fully reconstructed objects with the acceptance divided out. Large acceptance and efficiency differences in neighbouring points are due to statistical fluctuations.

The figure shows the signal acceptance (a) and efficiency (b) plots for the slepton pair production model, in the SR-0J $m_{\mathrm{T2}}^{100} \in[110,115)$ region. Acceptance is calculated by applying the signal region requirements to particle-level objects, which do not suffer from identification inefficiencies or mismeasurements. The efficiency is calculated with fully reconstructed objects with the acceptance divided out. Large acceptance and efficiency differences in neighbouring points are due to statistical fluctuations.

The figure shows the signal acceptance (a) and efficiency (b) plots for the slepton pair production model, in the SR-0J $m_{\mathrm{T2}}^{100} \in[115,120)$ region. Acceptance is calculated by applying the signal region requirements to particle-level objects, which do not suffer from identification inefficiencies or mismeasurements. The efficiency is calculated with fully reconstructed objects with the acceptance divided out. Large acceptance and efficiency differences in neighbouring points are due to statistical fluctuations.

The figure shows the signal acceptance (a) and efficiency (b) plots for the slepton pair production model, in the SR-0J $m_{\mathrm{T2}}^{100} \in[115,120)$ region. Acceptance is calculated by applying the signal region requirements to particle-level objects, which do not suffer from identification inefficiencies or mismeasurements. The efficiency is calculated with fully reconstructed objects with the acceptance divided out. Large acceptance and efficiency differences in neighbouring points are due to statistical fluctuations.

The figure shows the signal acceptance (a) and efficiency (b) plots for the slepton pair production model, in the SR-0J $m_{\mathrm{T2}}^{100} \in[120,125)$ region. Acceptance is calculated by applying the signal region requirements to particle-level objects, which do not suffer from identification inefficiencies or mismeasurements. The efficiency is calculated with fully reconstructed objects with the acceptance divided out. Large acceptance and efficiency differences in neighbouring points are due to statistical fluctuations.

The figure shows the signal acceptance (a) and efficiency (b) plots for the slepton pair production model, in the SR-0J $m_{\mathrm{T2}}^{100} \in[120,125)$ region. Acceptance is calculated by applying the signal region requirements to particle-level objects, which do not suffer from identification inefficiencies or mismeasurements. The efficiency is calculated with fully reconstructed objects with the acceptance divided out. Large acceptance and efficiency differences in neighbouring points are due to statistical fluctuations.

The figure shows the signal acceptance (a) and efficiency (b) plots for the slepton pair production model, in the SR-0J $m_{\mathrm{T2}}^{100} \in[125,130)$ region. Acceptance is calculated by applying the signal region requirements to particle-level objects, which do not suffer from identification inefficiencies or mismeasurements. The efficiency is calculated with fully reconstructed objects with the acceptance divided out. Large acceptance and efficiency differences in neighbouring points are due to statistical fluctuations.

The figure shows the signal acceptance (a) and efficiency (b) plots for the slepton pair production model, in the SR-0J $m_{\mathrm{T2}}^{100} \in[125,130)$ region. Acceptance is calculated by applying the signal region requirements to particle-level objects, which do not suffer from identification inefficiencies or mismeasurements. The efficiency is calculated with fully reconstructed objects with the acceptance divided out. Large acceptance and efficiency differences in neighbouring points are due to statistical fluctuations.

The figure shows the signal acceptance (a) and efficiency (b) plots for the slepton pair production model, in the SR-0J $m_{\mathrm{T2}}^{100} \in[130,140)$ region. Acceptance is calculated by applying the signal region requirements to particle-level objects, which do not suffer from identification inefficiencies or mismeasurements. The efficiency is calculated with fully reconstructed objects with the acceptance divided out. Large acceptance and efficiency differences in neighbouring points are due to statistical fluctuations.

The figure shows the signal acceptance (a) and efficiency (b) plots for the slepton pair production model, in the SR-0J $m_{\mathrm{T2}}^{100} \in[130,140)$ region. Acceptance is calculated by applying the signal region requirements to particle-level objects, which do not suffer from identification inefficiencies or mismeasurements. The efficiency is calculated with fully reconstructed objects with the acceptance divided out. Large acceptance and efficiency differences in neighbouring points are due to statistical fluctuations.

The figure shows the signal acceptance (a) and efficiency (b) plots for the slepton pair production model, in the SR-0J $m_{\mathrm{T2}}^{100} \in[140,\infty)$ region. Acceptance is calculated by applying the signal region requirements to particle-level objects, which do not suffer from identification inefficiencies or mismeasurements. The efficiency is calculated with fully reconstructed objects with the acceptance divided out. Large acceptance and efficiency differences in neighbouring points are due to statistical fluctuations.

The figure shows the signal acceptance (a) and efficiency (b) plots for the slepton pair production model, in the SR-0J $m_{\mathrm{T2}}^{100} \in[140,\infty)$ region. Acceptance is calculated by applying the signal region requirements to particle-level objects, which do not suffer from identification inefficiencies or mismeasurements. The efficiency is calculated with fully reconstructed objects with the acceptance divided out. Large acceptance and efficiency differences in neighbouring points are due to statistical fluctuations.

The figure shows the signal acceptance (a) and efficiency (b) plots for the slepton pair production model, in the SR-1J $m_{\mathrm{T2}}^{100} \in[100,\infty)$ region. Acceptance is calculated by applying the signal region requirements to particle-level objects, which do not suffer from identification inefficiencies or mismeasurements. The efficiency is calculated with fully reconstructed objects with the acceptance divided out. Large acceptance and efficiency differences in neighbouring points are due to statistical fluctuations.

The figure shows the signal acceptance (a) and efficiency (b) plots for the slepton pair production model, in the SR-1J $m_{\mathrm{T2}}^{100} \in[100,\infty)$ region. Acceptance is calculated by applying the signal region requirements to particle-level objects, which do not suffer from identification inefficiencies or mismeasurements. The efficiency is calculated with fully reconstructed objects with the acceptance divided out. Large acceptance and efficiency differences in neighbouring points are due to statistical fluctuations.

The figure shows the signal acceptance (a) and efficiency (b) plots for the slepton pair production model, in the SR-1J $m_{\mathrm{T2}}^{100} \in[110,\infty)$ region. Acceptance is calculated by applying the signal region requirements to particle-level objects, which do not suffer from identification inefficiencies or mismeasurements. The efficiency is calculated with fully reconstructed objects with the acceptance divided out. Large acceptance and efficiency differences in neighbouring points are due to statistical fluctuations.

The figure shows the signal acceptance (a) and efficiency (b) plots for the slepton pair production model, in the SR-1J $m_{\mathrm{T2}}^{100} \in[110,\infty)$ region. Acceptance is calculated by applying the signal region requirements to particle-level objects, which do not suffer from identification inefficiencies or mismeasurements. The efficiency is calculated with fully reconstructed objects with the acceptance divided out. Large acceptance and efficiency differences in neighbouring points are due to statistical fluctuations.

The figure shows the signal acceptance (a) and efficiency (b) plots for the slepton pair production model, in the SR-1J $m_{\mathrm{T2}}^{100} \in[120,\infty)$ region. Acceptance is calculated by applying the signal region requirements to particle-level objects, which do not suffer from identification inefficiencies or mismeasurements. The efficiency is calculated with fully reconstructed objects with the acceptance divided out. Large acceptance and efficiency differences in neighbouring points are due to statistical fluctuations.

The figure shows the signal acceptance (a) and efficiency (b) plots for the slepton pair production model, in the SR-1J $m_{\mathrm{T2}}^{100} \in[120,\infty)$ region. Acceptance is calculated by applying the signal region requirements to particle-level objects, which do not suffer from identification inefficiencies or mismeasurements. The efficiency is calculated with fully reconstructed objects with the acceptance divided out. Large acceptance and efficiency differences in neighbouring points are due to statistical fluctuations.

The figure shows the signal acceptance (a) and efficiency (b) plots for the slepton pair production model, in the SR-1J $m_{\mathrm{T2}}^{100} \in[130,\infty)$ region. Acceptance is calculated by applying the signal region requirements to particle-level objects, which do not suffer from identification inefficiencies or mismeasurements. The efficiency is calculated with fully reconstructed objects with the acceptance divided out. Large acceptance and efficiency differences in neighbouring points are due to statistical fluctuations.

The figure shows the signal acceptance (a) and efficiency (b) plots for the slepton pair production model, in the SR-1J $m_{\mathrm{T2}}^{100} \in[130,\infty)$ region. Acceptance is calculated by applying the signal region requirements to particle-level objects, which do not suffer from identification inefficiencies or mismeasurements. The efficiency is calculated with fully reconstructed objects with the acceptance divided out. Large acceptance and efficiency differences in neighbouring points are due to statistical fluctuations.

The figure shows the signal acceptance (a) and efficiency (b) plots for the slepton pair production model, in the SR-1J $m_{\mathrm{T2}}^{100} \in[100,105)$ region. Acceptance is calculated by applying the signal region requirements to particle-level objects, which do not suffer from identification inefficiencies or mismeasurements. The efficiency is calculated with fully reconstructed objects with the acceptance divided out. Large acceptance and efficiency differences in neighbouring points are due to statistical fluctuations.

The figure shows the signal acceptance (a) and efficiency (b) plots for the slepton pair production model, in the SR-1J $m_{\mathrm{T2}}^{100} \in[100,105)$ region. Acceptance is calculated by applying the signal region requirements to particle-level objects, which do not suffer from identification inefficiencies or mismeasurements. The efficiency is calculated with fully reconstructed objects with the acceptance divided out. Large acceptance and efficiency differences in neighbouring points are due to statistical fluctuations.

The figure shows the signal acceptance (a) and efficiency (b) plots for the slepton pair production model, in the SR-1J $m_{\mathrm{T2}}^{100} \in[105,110)$ region. Acceptance is calculated by applying the signal region requirements to particle-level objects, which do not suffer from identification inefficiencies or mismeasurements. The efficiency is calculated with fully reconstructed objects with the acceptance divided out. Large acceptance and efficiency differences in neighbouring points are due to statistical fluctuations.

The figure shows the signal acceptance (a) and efficiency (b) plots for the slepton pair production model, in the SR-1J $m_{\mathrm{T2}}^{100} \in[105,110)$ region. Acceptance is calculated by applying the signal region requirements to particle-level objects, which do not suffer from identification inefficiencies or mismeasurements. The efficiency is calculated with fully reconstructed objects with the acceptance divided out. Large acceptance and efficiency differences in neighbouring points are due to statistical fluctuations.

The figure shows the signal acceptance (a) and efficiency (b) plots for the slepton pair production model, in the SR-1J $m_{\mathrm{T2}}^{100} \in[110,115)$ region. Acceptance is calculated by applying the signal region requirements to particle-level objects, which do not suffer from identification inefficiencies or mismeasurements. The efficiency is calculated with fully reconstructed objects with the acceptance divided out. Large acceptance and efficiency differences in neighbouring points are due to statistical fluctuations.

The figure shows the signal acceptance (a) and efficiency (b) plots for the slepton pair production model, in the SR-1J $m_{\mathrm{T2}}^{100} \in[110,115)$ region. Acceptance is calculated by applying the signal region requirements to particle-level objects, which do not suffer from identification inefficiencies or mismeasurements. The efficiency is calculated with fully reconstructed objects with the acceptance divided out. Large acceptance and efficiency differences in neighbouring points are due to statistical fluctuations.

The figure shows the signal acceptance (a) and efficiency (b) plots for the slepton pair production model, in the SR-1J $m_{\mathrm{T2}}^{100} \in[115,120)$ region. Acceptance is calculated by applying the signal region requirements to particle-level objects, which do not suffer from identification inefficiencies or mismeasurements. The efficiency is calculated with fully reconstructed objects with the acceptance divided out. Large acceptance and efficiency differences in neighbouring points are due to statistical fluctuations.

The figure shows the signal acceptance (a) and efficiency (b) plots for the slepton pair production model, in the SR-1J $m_{\mathrm{T2}}^{100} \in[115,120)$ region. Acceptance is calculated by applying the signal region requirements to particle-level objects, which do not suffer from identification inefficiencies or mismeasurements. The efficiency is calculated with fully reconstructed objects with the acceptance divided out. Large acceptance and efficiency differences in neighbouring points are due to statistical fluctuations.

The figure shows the signal acceptance (a) and efficiency (b) plots for the slepton pair production model, in the SR-1J $m_{\mathrm{T2}}^{100} \in[120,125)$ region. Acceptance is calculated by applying the signal region requirements to particle-level objects, which do not suffer from identification inefficiencies or mismeasurements. The efficiency is calculated with fully reconstructed objects with the acceptance divided out. Large acceptance and efficiency differences in neighbouring points are due to statistical fluctuations.

The figure shows the signal acceptance (a) and efficiency (b) plots for the slepton pair production model, in the SR-1J $m_{\mathrm{T2}}^{100} \in[120,125)$ region. Acceptance is calculated by applying the signal region requirements to particle-level objects, which do not suffer from identification inefficiencies or mismeasurements. The efficiency is calculated with fully reconstructed objects with the acceptance divided out. Large acceptance and efficiency differences in neighbouring points are due to statistical fluctuations.

The figure shows the signal acceptance (a) and efficiency (b) plots for the slepton pair production model, in the SR-1J $m_{\mathrm{T2}}^{100} \in[125,130)$ region. Acceptance is calculated by applying the signal region requirements to particle-level objects, which do not suffer from identification inefficiencies or mismeasurements. The efficiency is calculated with fully reconstructed objects with the acceptance divided out. Large acceptance and efficiency differences in neighbouring points are due to statistical fluctuations.

The figure shows the signal acceptance (a) and efficiency (b) plots for the slepton pair production model, in the SR-1J $m_{\mathrm{T2}}^{100} \in[125,130)$ region. Acceptance is calculated by applying the signal region requirements to particle-level objects, which do not suffer from identification inefficiencies or mismeasurements. The efficiency is calculated with fully reconstructed objects with the acceptance divided out. Large acceptance and efficiency differences in neighbouring points are due to statistical fluctuations.

The figure shows the signal acceptance (a) and efficiency (b) plots for the slepton pair production model, in the SR-1J $m_{\mathrm{T2}}^{100} \in[130,140)$ region. Acceptance is calculated by applying the signal region requirements to particle-level objects, which do not suffer from identification inefficiencies or mismeasurements. The efficiency is calculated with fully reconstructed objects with the acceptance divided out. Large acceptance and efficiency differences in neighbouring points are due to statistical fluctuations.

The figure shows the signal acceptance (a) and efficiency (b) plots for the slepton pair production model, in the SR-1J $m_{\mathrm{T2}}^{100} \in[130,140)$ region. Acceptance is calculated by applying the signal region requirements to particle-level objects, which do not suffer from identification inefficiencies or mismeasurements. The efficiency is calculated with fully reconstructed objects with the acceptance divided out. Large acceptance and efficiency differences in neighbouring points are due to statistical fluctuations.

The figure shows the signal acceptance (a) and efficiency (b) plots for the slepton pair production model, in the SR-1J $m_{\mathrm{T2}}^{100} \in[140,\infty)$ region. Acceptance is calculated by applying the signal region requirements to particle-level objects, which do not suffer from identification inefficiencies or mismeasurements. The efficiency is calculated with fully reconstructed objects with the acceptance divided out. Large acceptance and efficiency differences in neighbouring points are due to statistical fluctuations.

The figure shows the signal acceptance (a) and efficiency (b) plots for the slepton pair production model, in the SR-1J $m_{\mathrm{T2}}^{100} \in[140,\infty)$ region. Acceptance is calculated by applying the signal region requirements to particle-level objects, which do not suffer from identification inefficiencies or mismeasurements. The efficiency is calculated with fully reconstructed objects with the acceptance divided out. Large acceptance and efficiency differences in neighbouring points are due to statistical fluctuations.

Cutflow table for the slepton signal sample with $m(\tilde{\ell},\tilde{\chi}_1^0) = (100,70)$ GeV, in the SR-0J $m_{\mathrm{T2}}^{100} \in [100,\infty)$ region. The yields include the process cross section and are weighted to the 139 fb$^{-1}$ luminosity. 246000 events were generated for the sample.

Cutflow table for the slepton signal sample with $m(\tilde{\ell},\tilde{\chi}_1^0) = (100,70)$ GeV, in the SR-1J $m_{\mathrm{T2}}^{100} \in [100,\infty)$ region. The yields include the process cross section and are weighted to the 139 fb$^{-1}$ luminosity. 246000 events were generated for the sample.

Observed and expected exclusion limits on SUSY simplified models, with observed upper limits on signal cross-section (fb) overlaid, for slepton-pair production in the $m(\tilde{\ell})-m(\tilde{\chi}_1^0)$ plane. Only $\tilde{e}$ and $\tilde{\mu}$ are considered. The observed (solid thick line) and expected (thin dashed line) exclusion contours are indicated. The shaded band around the dashed line corresponds to the $\pm 1 \sigma$ variations in the expected limit, including all uncertainties except theoretical uncertainties in the signal cross-section. The dotted lines around the observed limit illustrate the change in the observed limit as the nominal signal cross-section is scaled up and down by the theoretical uncertainty. All limits are computed at 95% CL. The observed limits obtained at LEP for $\tilde{\mu}_{\textup{R}}$ and by the ATLAS experiment in previous searches are also shown.

Observed and expected exclusion limits on SUSY simplified models for slepton-pair production in the (a) $m(\tilde{\ell})-m(\tilde{\chi}_1^0)$ and (b) $m(\tilde{\ell})-\Delta m(\tilde{\ell},\tilde{\chi}_1^0)$ planes. Only $\tilde{e}$ and $\tilde{\mu}$ are considered. The observed (solid thick line) and expected (thin dashed line) exclusion contours are indicated. The shaded band around the dashed line corresponds to the $\pm 1 \sigma$ variations in the expected limit, including all uncertainties except theoretical uncertainties in the signal cross-section. The dotted lines around the observed limit illustrate the change in the observed limit as the nominal signal cross-section is scaled up and down by the theoretical uncertainty. All limits are computed at 95% CL. The observed limits obtained at LEP for $\tilde{\mu}_{\textup{R}}$ and by the ATLAS experiment in previous searches are also shown.

Observed and expected exclusion limits on SUSY simplified models for slepton-pair production in the (a) $m(\tilde{\ell})-m(\tilde{\chi}_1^0)$ and (b) $m(\tilde{\ell})-\Delta m(\tilde{\ell},\tilde{\chi}_1^0)$ planes. Only $\tilde{e}$ and $\tilde{\mu}$ are considered. The observed (solid thick line) and expected (thin dashed line) exclusion contours are indicated. The shaded band around the dashed line corresponds to the $\pm 1 \sigma$ variations in the expected limit, including all uncertainties except theoretical uncertainties in the signal cross-section. The dotted lines around the observed limit illustrate the change in the observed limit as the nominal signal cross-section is scaled up and down by the theoretical uncertainty. All limits are computed at 95% CL. The observed limits obtained at LEP for $\tilde{\mu}_{\textup{R}}$ and by the ATLAS experiment in previous searches are also shown.

Observed and expected exclusion limits on SUSY simplified models for slepton-pair production in the (a) $m(\tilde{\ell})-m(\tilde{\chi}_1^0)$ and (b) $m(\tilde{\ell})-\Delta m(\tilde{\ell},\tilde{\chi}_1^0)$ planes. Only $\tilde{e}$ and $\tilde{\mu}$ are considered. The observed (solid thick line) and expected (thin dashed line) exclusion contours are indicated. The shaded band around the dashed line corresponds to the $\pm 1 \sigma$ variations in the expected limit, including all uncertainties except theoretical uncertainties in the signal cross-section. The dotted lines around the observed limit illustrate the change in the observed limit as the nominal signal cross-section is scaled up and down by the theoretical uncertainty. All limits are computed at 95% CL. The observed limits obtained at LEP for $\tilde{\mu}_{\textup{R}}$ and by the ATLAS experiment in previous searches are also shown.

Observed and expected exclusion limits on SUSY simplified models for slepton-pair production in the (a) $m(\tilde{\ell})-m(\tilde{\chi}_1^0)$ and (b) $m(\tilde{\ell})-\Delta m(\tilde{\ell},\tilde{\chi}_1^0)$ planes. Only $\tilde{e}$ and $\tilde{\mu}$ are considered. The observed (solid thick line) and expected (thin dashed line) exclusion contours are indicated. The shaded band around the dashed line corresponds to the $\pm 1 \sigma$ variations in the expected limit, including all uncertainties except theoretical uncertainties in the signal cross-section. The dotted lines around the observed limit illustrate the change in the observed limit as the nominal signal cross-section is scaled up and down by the theoretical uncertainty. All limits are computed at 95% CL. The observed limits obtained at LEP for $\tilde{\mu}_{\textup{R}}$ and by the ATLAS experiment in previous searches are also shown.

Observed and expected exclusion limits on SUSY simplified models for slepton-pair production in the (a) $m(\tilde{\ell})-m(\tilde{\chi}_1^0)$ and (b) $m(\tilde{\ell})-\Delta m(\tilde{\ell},\tilde{\chi}_1^0)$ planes. Only $\tilde{e}$ and $\tilde{\mu}$ are considered. The observed (solid thick line) and expected (thin dashed line) exclusion contours are indicated. The shaded band around the dashed line corresponds to the $\pm 1 \sigma$ variations in the expected limit, including all uncertainties except theoretical uncertainties in the signal cross-section. The dotted lines around the observed limit illustrate the change in the observed limit as the nominal signal cross-section is scaled up and down by the theoretical uncertainty. All limits are computed at 95% CL. The observed limits obtained at LEP for $\tilde{\mu}_{\textup{R}}$ and by the ATLAS experiment in previous searches are also shown.

Observed and expected exclusion limits on SUSY simplified models for slepton-pair production in the (a) $m(\tilde{\ell})-m(\tilde{\chi}_1^0)$ and (b) $m(\tilde{\ell})-\Delta m(\tilde{\ell},\tilde{\chi}_1^0)$ planes. Only $\tilde{e}$ and $\tilde{\mu}$ are considered. The observed (solid thick line) and expected (thin dashed line) exclusion contours are indicated. The shaded band around the dashed line corresponds to the $\pm 1 \sigma$ variations in the expected limit, including all uncertainties except theoretical uncertainties in the signal cross-section. The dotted lines around the observed limit illustrate the change in the observed limit as the nominal signal cross-section is scaled up and down by the theoretical uncertainty. All limits are computed at 95% CL. The observed limits obtained at LEP for $\tilde{\mu}_{\textup{R}}$ and by the ATLAS experiment in previous searches are also shown.

Observed and expected exclusion limits on SUSY simplified models for slepton-pair production in the (a) $m(\tilde{\ell})-m(\tilde{\chi}_1^0)$ and (b) $m(\tilde{\ell})-\Delta m(\tilde{\ell},\tilde{\chi}_1^0)$ planes. Only $\tilde{e}$ and $\tilde{\mu}$ are considered. The observed (solid thick line) and expected (thin dashed line) exclusion contours are indicated. The shaded band around the dashed line corresponds to the $\pm 1 \sigma$ variations in the expected limit, including all uncertainties except theoretical uncertainties in the signal cross-section. The dotted lines around the observed limit illustrate the change in the observed limit as the nominal signal cross-section is scaled up and down by the theoretical uncertainty. All limits are computed at 95% CL. The observed limits obtained at LEP for $\tilde{\mu}_{\textup{R}}$ and by the ATLAS experiment in previous searches are also shown.

Observed and expected exclusion limits on SUSY simplified models for slepton-pair production in the (a) $m(\tilde{\ell})-m(\tilde{\chi}_1^0)$ and (b) $m(\tilde{\ell})-\Delta m(\tilde{\ell},\tilde{\chi}_1^0)$ planes. Only $\tilde{e}$ and $\tilde{\mu}$ are considered. The observed (solid thick line) and expected (thin dashed line) exclusion contours are indicated. The shaded band around the dashed line corresponds to the $\pm 1 \sigma$ variations in the expected limit, including all uncertainties except theoretical uncertainties in the signal cross-section. The dotted lines around the observed limit illustrate the change in the observed limit as the nominal signal cross-section is scaled up and down by the theoretical uncertainty. All limits are computed at 95% CL. The observed limits obtained at LEP for $\tilde{\mu}_{\textup{R}}$ and by the ATLAS experiment in previous searches are also shown.

Observed and expected exclusion limits on SUSY simplified models for slepton-pair production in the (a) $m(\tilde{\ell})-m(\tilde{\chi}_1^0)$ and (b) $m(\tilde{\ell})-\Delta m(\tilde{\ell},\tilde{\chi}_1^0)$ planes. Only $\tilde{e}$ and $\tilde{\mu}$ are considered. The observed (solid thick line) and expected (thin dashed line) exclusion contours are indicated. The shaded band around the dashed line corresponds to the $\pm 1 \sigma$ variations in the expected limit, including all uncertainties except theoretical uncertainties in the signal cross-section. The dotted lines around the observed limit illustrate the change in the observed limit as the nominal signal cross-section is scaled up and down by the theoretical uncertainty. All limits are computed at 95% CL. The observed limits obtained at LEP for $\tilde{\mu}_{\textup{R}}$ and by the ATLAS experiment in previous searches are also shown.

Observed and expected exclusion limits on SUSY simplified models for slepton-pair production in the (a) $m(\tilde{\ell})-m(\tilde{\chi}_1^0)$ and (b) $m(\tilde{\ell})-\Delta m(\tilde{\ell},\tilde{\chi}_1^0)$ planes. Only $\tilde{e}$ and $\tilde{\mu}$ are considered. The observed (solid thick line) and expected (thin dashed line) exclusion contours are indicated. The shaded band around the dashed line corresponds to the $\pm 1 \sigma$ variations in the expected limit, including all uncertainties except theoretical uncertainties in the signal cross-section. The dotted lines around the observed limit illustrate the change in the observed limit as the nominal signal cross-section is scaled up and down by the theoretical uncertainty. All limits are computed at 95% CL. The observed limits obtained at LEP for $\tilde{\mu}_{\textup{R}}$ and by the ATLAS experiment in previous searches are also shown.

Observed and expected exclusion limits on SUSY simplified models for slepton-pair production in the (a) $m(\tilde{\ell})-m(\tilde{\chi}_1^0)$ and (b) $m(\tilde{\ell})-\Delta m(\tilde{\ell},\tilde{\chi}_1^0)$ planes. Only $\tilde{e}$ and $\tilde{\mu}$ are considered. The observed (solid thick line) and expected (thin dashed line) exclusion contours are indicated. The shaded band around the dashed line corresponds to the $\pm 1 \sigma$ variations in the expected limit, including all uncertainties except theoretical uncertainties in the signal cross-section. The dotted lines around the observed limit illustrate the change in the observed limit as the nominal signal cross-section is scaled up and down by the theoretical uncertainty. All limits are computed at 95% CL. The observed limits obtained at LEP for $\tilde{\mu}_{\textup{R}}$ and by the ATLAS experiment in previous searches are also shown.

Observed and expected exclusion limits on SUSY simplified models for slepton-pair production in the (a) $m(\tilde{\ell})-m(\tilde{\chi}_1^0)$ and (b) $m(\tilde{\ell})-\Delta m(\tilde{\ell},\tilde{\chi}_1^0)$ planes. Only $\tilde{e}$ and $\tilde{\mu}$ are considered. The observed (solid thick line) and expected (thin dashed line) exclusion contours are indicated. The shaded band around the dashed line corresponds to the $\pm 1 \sigma$ variations in the expected limit, including all uncertainties except theoretical uncertainties in the signal cross-section. The dotted lines around the observed limit illustrate the change in the observed limit as the nominal signal cross-section is scaled up and down by the theoretical uncertainty. All limits are computed at 95% CL. The observed limits obtained at LEP for $\tilde{\mu}_{\textup{R}}$ and by the ATLAS experiment in previous searches are also shown.

Observed and expected exclusion limits on SUSY simplified models for direct selectron production in the (a) $m(\tilde{e})-m(\tilde{\chi}_1^0)$ and (c) $m(\tilde{e})-\Delta m(\tilde{e},\tilde{\chi}_1^0)$ planes, and for direct smuon production in the (b) $m(\tilde{\mu})-m(\tilde{\chi}_1^0)$ and (d) $m(\tilde{\mu})-\Delta m(\tilde{\mu},\tilde{\chi}_1^0)$ planes. In Figure (a) and (c) the observed (solid thick lines) and expected (dashed lines) exclusion contours are indicated for combined $\tilde{e}_{\textup{L,R}}$ and for $\tilde{e}_{\textup{L}}$ and $\tilde{e}_{\textup{R}}$. In Figure (b) and (d) the observed (solid thick lines) and expected (dashed lines) exclusion contours are indicated for combined $\tilde{\mu}_{\textup{L,R}}$ and for $\tilde{\mu}_{\textup{L}}$. No unique sensitivity to $\tilde{\mu}_{\textup{R}}$ is observed. All limits are computed at 95% CL. The observed limits obtained at LEP and by the ATLAS experiment in previous searches are also shown in the shaded areas.

Observed and expected exclusion limits on SUSY simplified models for direct selectron production in the (a) $m(\tilde{e})-m(\tilde{\chi}_1^0)$ and (c) $m(\tilde{e})-\Delta m(\tilde{e},\tilde{\chi}_1^0)$ planes, and for direct smuon production in the (b) $m(\tilde{\mu})-m(\tilde{\chi}_1^0)$ and (d) $m(\tilde{\mu})-\Delta m(\tilde{\mu},\tilde{\chi}_1^0)$ planes. In Figure (a) and (c) the observed (solid thick lines) and expected (dashed lines) exclusion contours are indicated for combined $\tilde{e}_{\textup{L,R}}$ and for $\tilde{e}_{\textup{L}}$ and $\tilde{e}_{\textup{R}}$. In Figure (b) and (d) the observed (solid thick lines) and expected (dashed lines) exclusion contours are indicated for combined $\tilde{\mu}_{\textup{L,R}}$ and for $\tilde{\mu}_{\textup{L}}$. No unique sensitivity to $\tilde{\mu}_{\textup{R}}$ is observed. All limits are computed at 95% CL. The observed limits obtained at LEP and by the ATLAS experiment in previous searches are also shown in the shaded areas.

Observed and expected exclusion limits on SUSY simplified models for direct selectron production in the (a) $m(\tilde{e})-m(\tilde{\chi}_1^0)$ and (c) $m(\tilde{e})-\Delta m(\tilde{e},\tilde{\chi}_1^0)$ planes, and for direct smuon production in the (b) $m(\tilde{\mu})-m(\tilde{\chi}_1^0)$ and (d) $m(\tilde{\mu})-\Delta m(\tilde{\mu},\tilde{\chi}_1^0)$ planes. In Figure (a) and (c) the observed (solid thick lines) and expected (dashed lines) exclusion contours are indicated for combined $\tilde{e}_{\textup{L,R}}$ and for $\tilde{e}_{\textup{L}}$ and $\tilde{e}_{\textup{R}}$. In Figure (b) and (d) the observed (solid thick lines) and expected (dashed lines) exclusion contours are indicated for combined $\tilde{\mu}_{\textup{L,R}}$ and for $\tilde{\mu}_{\textup{L}}$. No unique sensitivity to $\tilde{\mu}_{\textup{R}}$ is observed. All limits are computed at 95% CL. The observed limits obtained at LEP and by the ATLAS experiment in previous searches are also shown in the shaded areas.

Observed and expected exclusion limits on SUSY simplified models for direct selectron production in the (a) $m(\tilde{e})-m(\tilde{\chi}_1^0)$ and (c) $m(\tilde{e})-\Delta m(\tilde{e},\tilde{\chi}_1^0)$ planes, and for direct smuon production in the (b) $m(\tilde{\mu})-m(\tilde{\chi}_1^0)$ and (d) $m(\tilde{\mu})-\Delta m(\tilde{\mu},\tilde{\chi}_1^0)$ planes. In Figure (a) and (c) the observed (solid thick lines) and expected (dashed lines) exclusion contours are indicated for combined $\tilde{e}_{\textup{L,R}}$ and for $\tilde{e}_{\textup{L}}$ and $\tilde{e}_{\textup{R}}$. In Figure (b) and (d) the observed (solid thick lines) and expected (dashed lines) exclusion contours are indicated for combined $\tilde{\mu}_{\textup{L,R}}$ and for $\tilde{\mu}_{\textup{L}}$. No unique sensitivity to $\tilde{\mu}_{\textup{R}}$ is observed. All limits are computed at 95% CL. The observed limits obtained at LEP and by the ATLAS experiment in previous searches are also shown in the shaded areas.

Observed and expected exclusion limits on SUSY simplified models for direct selectron production in the (a) $m(\tilde{e})-m(\tilde{\chi}_1^0)$ and (c) $m(\tilde{e})-\Delta m(\tilde{e},\tilde{\chi}_1^0)$ planes, and for direct smuon production in the (b) $m(\tilde{\mu})-m(\tilde{\chi}_1^0)$ and (d) $m(\tilde{\mu})-\Delta m(\tilde{\mu},\tilde{\chi}_1^0)$ planes. In Figure (a) and (c) the observed (solid thick lines) and expected (dashed lines) exclusion contours are indicated for combined $\tilde{e}_{\textup{L,R}}$ and for $\tilde{e}_{\textup{L}}$ and $\tilde{e}_{\textup{R}}$. In Figure (b) and (d) the observed (solid thick lines) and expected (dashed lines) exclusion contours are indicated for combined $\tilde{\mu}_{\textup{L,R}}$ and for $\tilde{\mu}_{\textup{L}}$. No unique sensitivity to $\tilde{\mu}_{\textup{R}}$ is observed. All limits are computed at 95% CL. The observed limits obtained at LEP and by the ATLAS experiment in previous searches are also shown in the shaded areas.

Observed and expected exclusion limits on SUSY simplified models for direct selectron production in the (a) $m(\tilde{e})-m(\tilde{\chi}_1^0)$ and (c) $m(\tilde{e})-\Delta m(\tilde{e},\tilde{\chi}_1^0)$ planes, and for direct smuon production in the (b) $m(\tilde{\mu})-m(\tilde{\chi}_1^0)$ and (d) $m(\tilde{\mu})-\Delta m(\tilde{\mu},\tilde{\chi}_1^0)$ planes. In Figure (a) and (c) the observed (solid thick lines) and expected (dashed lines) exclusion contours are indicated for combined $\tilde{e}_{\textup{L,R}}$ and for $\tilde{e}_{\textup{L}}$ and $\tilde{e}_{\textup{R}}$. In Figure (b) and (d) the observed (solid thick lines) and expected (dashed lines) exclusion contours are indicated for combined $\tilde{\mu}_{\textup{L,R}}$ and for $\tilde{\mu}_{\textup{L}}$. No unique sensitivity to $\tilde{\mu}_{\textup{R}}$ is observed. All limits are computed at 95% CL. The observed limits obtained at LEP and by the ATLAS experiment in previous searches are also shown in the shaded areas.

Observed and expected exclusion limits on SUSY simplified models for direct selectron production in the (a) $m(\tilde{e})-m(\tilde{\chi}_1^0)$ and (c) $m(\tilde{e})-\Delta m(\tilde{e},\tilde{\chi}_1^0)$ planes, and for direct smuon production in the (b) $m(\tilde{\mu})-m(\tilde{\chi}_1^0)$ and (d) $m(\tilde{\mu})-\Delta m(\tilde{\mu},\tilde{\chi}_1^0)$ planes. In Figure (a) and (c) the observed (solid thick lines) and expected (dashed lines) exclusion contours are indicated for combined $\tilde{e}_{\textup{L,R}}$ and for $\tilde{e}_{\textup{L}}$ and $\tilde{e}_{\textup{R}}$. In Figure (b) and (d) the observed (solid thick lines) and expected (dashed lines) exclusion contours are indicated for combined $\tilde{\mu}_{\textup{L,R}}$ and for $\tilde{\mu}_{\textup{L}}$. No unique sensitivity to $\tilde{\mu}_{\textup{R}}$ is observed. All limits are computed at 95% CL. The observed limits obtained at LEP and by the ATLAS experiment in previous searches are also shown in the shaded areas.

Observed and expected exclusion limits on SUSY simplified models for direct selectron production in the (a) $m(\tilde{e})-m(\tilde{\chi}_1^0)$ and (c) $m(\tilde{e})-\Delta m(\tilde{e},\tilde{\chi}_1^0)$ planes, and for direct smuon production in the (b) $m(\tilde{\mu})-m(\tilde{\chi}_1^0)$ and (d) $m(\tilde{\mu})-\Delta m(\tilde{\mu},\tilde{\chi}_1^0)$ planes. In Figure (a) and (c) the observed (solid thick lines) and expected (dashed lines) exclusion contours are indicated for combined $\tilde{e}_{\textup{L,R}}$ and for $\tilde{e}_{\textup{L}}$ and $\tilde{e}_{\textup{R}}$. In Figure (b) and (d) the observed (solid thick lines) and expected (dashed lines) exclusion contours are indicated for combined $\tilde{\mu}_{\textup{L,R}}$ and for $\tilde{\mu}_{\textup{L}}$. No unique sensitivity to $\tilde{\mu}_{\textup{R}}$ is observed. All limits are computed at 95% CL. The observed limits obtained at LEP and by the ATLAS experiment in previous searches are also shown in the shaded areas.

Observed and expected exclusion limits on SUSY simplified models for direct selectron production in the (a) $m(\tilde{e})-m(\tilde{\chi}_1^0)$ and (c) $m(\tilde{e})-\Delta m(\tilde{e},\tilde{\chi}_1^0)$ planes, and for direct smuon production in the (b) $m(\tilde{\mu})-m(\tilde{\chi}_1^0)$ and (d) $m(\tilde{\mu})-\Delta m(\tilde{\mu},\tilde{\chi}_1^0)$ planes. In Figure (a) and (c) the observed (solid thick lines) and expected (dashed lines) exclusion contours are indicated for combined $\tilde{e}_{\textup{L,R}}$ and for $\tilde{e}_{\textup{L}}$ and $\tilde{e}_{\textup{R}}$. In Figure (b) and (d) the observed (solid thick lines) and expected (dashed lines) exclusion contours are indicated for combined $\tilde{\mu}_{\textup{L,R}}$ and for $\tilde{\mu}_{\textup{L}}$. No unique sensitivity to $\tilde{\mu}_{\textup{R}}$ is observed. All limits are computed at 95% CL. The observed limits obtained at LEP and by the ATLAS experiment in previous searches are also shown in the shaded areas.

Observed and expected exclusion limits on SUSY simplified models for direct selectron production in the (a) $m(\tilde{e})-m(\tilde{\chi}_1^0)$ and (c) $m(\tilde{e})-\Delta m(\tilde{e},\tilde{\chi}_1^0)$ planes, and for direct smuon production in the (b) $m(\tilde{\mu})-m(\tilde{\chi}_1^0)$ and (d) $m(\tilde{\mu})-\Delta m(\tilde{\mu},\tilde{\chi}_1^0)$ planes. In Figure (a) and (c) the observed (solid thick lines) and expected (dashed lines) exclusion contours are indicated for combined $\tilde{e}_{\textup{L,R}}$ and for $\tilde{e}_{\textup{L}}$ and $\tilde{e}_{\textup{R}}$. In Figure (b) and (d) the observed (solid thick lines) and expected (dashed lines) exclusion contours are indicated for combined $\tilde{\mu}_{\textup{L,R}}$ and for $\tilde{\mu}_{\textup{L}}$. No unique sensitivity to $\tilde{\mu}_{\textup{R}}$ is observed. All limits are computed at 95% CL. The observed limits obtained at LEP and by the ATLAS experiment in previous searches are also shown in the shaded areas.

Observed and expected exclusion limits on SUSY simplified models for direct selectron production in the (a) $m(\tilde{e})-m(\tilde{\chi}_1^0)$ and (c) $m(\tilde{e})-\Delta m(\tilde{e},\tilde{\chi}_1^0)$ planes, and for direct smuon production in the (b) $m(\tilde{\mu})-m(\tilde{\chi}_1^0)$ and (d) $m(\tilde{\mu})-\Delta m(\tilde{\mu},\tilde{\chi}_1^0)$ planes. In Figure (a) and (c) the observed (solid thick lines) and expected (dashed lines) exclusion contours are indicated for combined $\tilde{e}_{\textup{L,R}}$ and for $\tilde{e}_{\textup{L}}$ and $\tilde{e}_{\textup{R}}$. In Figure (b) and (d) the observed (solid thick lines) and expected (dashed lines) exclusion contours are indicated for combined $\tilde{\mu}_{\textup{L,R}}$ and for $\tilde{\mu}_{\textup{L}}$. No unique sensitivity to $\tilde{\mu}_{\textup{R}}$ is observed. All limits are computed at 95% CL. The observed limits obtained at LEP and by the ATLAS experiment in previous searches are also shown in the shaded areas.

Observed and expected exclusion limits on SUSY simplified models for direct selectron production in the (a) $m(\tilde{e})-m(\tilde{\chi}_1^0)$ and (c) $m(\tilde{e})-\Delta m(\tilde{e},\tilde{\chi}_1^0)$ planes, and for direct smuon production in the (b) $m(\tilde{\mu})-m(\tilde{\chi}_1^0)$ and (d) $m(\tilde{\mu})-\Delta m(\tilde{\mu},\tilde{\chi}_1^0)$ planes. In Figure (a) and (c) the observed (solid thick lines) and expected (dashed lines) exclusion contours are indicated for combined $\tilde{e}_{\textup{L,R}}$ and for $\tilde{e}_{\textup{L}}$ and $\tilde{e}_{\textup{R}}$. In Figure (b) and (d) the observed (solid thick lines) and expected (dashed lines) exclusion contours are indicated for combined $\tilde{\mu}_{\textup{L,R}}$ and for $\tilde{\mu}_{\textup{L}}$. No unique sensitivity to $\tilde{\mu}_{\textup{R}}$ is observed. All limits are computed at 95% CL. The observed limits obtained at LEP and by the ATLAS experiment in previous searches are also shown in the shaded areas.

Observed and expected exclusion limits on SUSY simplified models for direct selectron production in the (a) $m(\tilde{e})-m(\tilde{\chi}_1^0)$ and (c) $m(\tilde{e})-\Delta m(\tilde{e},\tilde{\chi}_1^0)$ planes, and for direct smuon production in the (b) $m(\tilde{\mu})-m(\tilde{\chi}_1^0)$ and (d) $m(\tilde{\mu})-\Delta m(\tilde{\mu},\tilde{\chi}_1^0)$ planes. In Figure (a) and (c) the observed (solid thick lines) and expected (dashed lines) exclusion contours are indicated for combined $\tilde{e}_{\textup{L,R}}$ and for $\tilde{e}_{\textup{L}}$ and $\tilde{e}_{\textup{R}}$. In Figure (b) and (d) the observed (solid thick lines) and expected (dashed lines) exclusion contours are indicated for combined $\tilde{\mu}_{\textup{L,R}}$ and for $\tilde{\mu}_{\textup{L}}$. No unique sensitivity to $\tilde{\mu}_{\textup{R}}$ is observed. All limits are computed at 95% CL. The observed limits obtained at LEP and by the ATLAS experiment in previous searches are also shown in the shaded areas.

Observed and expected exclusion limits on SUSY simplified models for direct selectron production in the (a) $m(\tilde{e})-m(\tilde{\chi}_1^0)$ and (c) $m(\tilde{e})-\Delta m(\tilde{e},\tilde{\chi}_1^0)$ planes, and for direct smuon production in the (b) $m(\tilde{\mu})-m(\tilde{\chi}_1^0)$ and (d) $m(\tilde{\mu})-\Delta m(\tilde{\mu},\tilde{\chi}_1^0)$ planes. In Figure (a) and (c) the observed (solid thick lines) and expected (dashed lines) exclusion contours are indicated for combined $\tilde{e}_{\textup{L,R}}$ and for $\tilde{e}_{\textup{L}}$ and $\tilde{e}_{\textup{R}}$. In Figure (b) and (d) the observed (solid thick lines) and expected (dashed lines) exclusion contours are indicated for combined $\tilde{\mu}_{\textup{L,R}}$ and for $\tilde{\mu}_{\textup{L}}$. No unique sensitivity to $\tilde{\mu}_{\textup{R}}$ is observed. All limits are computed at 95% CL. The observed limits obtained at LEP and by the ATLAS experiment in previous searches are also shown in the shaded areas.

Observed and expected exclusion limits on SUSY simplified models for direct selectron production in the (a) $m(\tilde{e})-m(\tilde{\chi}_1^0)$ and (c) $m(\tilde{e})-\Delta m(\tilde{e},\tilde{\chi}_1^0)$ planes, and for direct smuon production in the (b) $m(\tilde{\mu})-m(\tilde{\chi}_1^0)$ and (d) $m(\tilde{\mu})-\Delta m(\tilde{\mu},\tilde{\chi}_1^0)$ planes. In Figure (a) and (c) the observed (solid thick lines) and expected (dashed lines) exclusion contours are indicated for combined $\tilde{e}_{\textup{L,R}}$ and for $\tilde{e}_{\textup{L}}$ and $\tilde{e}_{\textup{R}}$. In Figure (b) and (d) the observed (solid thick lines) and expected (dashed lines) exclusion contours are indicated for combined $\tilde{\mu}_{\textup{L,R}}$ and for $\tilde{\mu}_{\textup{L}}$. No unique sensitivity to $\tilde{\mu}_{\textup{R}}$ is observed. All limits are computed at 95% CL. The observed limits obtained at LEP and by the ATLAS experiment in previous searches are also shown in the shaded areas.

Observed and expected exclusion limits on SUSY simplified models for direct selectron production in the (a) $m(\tilde{e})-m(\tilde{\chi}_1^0)$ and (c) $m(\tilde{e})-\Delta m(\tilde{e},\tilde{\chi}_1^0)$ planes, and for direct smuon production in the (b) $m(\tilde{\mu})-m(\tilde{\chi}_1^0)$ and (d) $m(\tilde{\mu})-\Delta m(\tilde{\mu},\tilde{\chi}_1^0)$ planes. In Figure (a) and (c) the observed (solid thick lines) and expected (dashed lines) exclusion contours are indicated for combined $\tilde{e}_{\textup{L,R}}$ and for $\tilde{e}_{\textup{L}}$ and $\tilde{e}_{\textup{R}}$. In Figure (b) and (d) the observed (solid thick lines) and expected (dashed lines) exclusion contours are indicated for combined $\tilde{\mu}_{\textup{L,R}}$ and for $\tilde{\mu}_{\textup{L}}$. No unique sensitivity to $\tilde{\mu}_{\textup{R}}$ is observed. All limits are computed at 95% CL. The observed limits obtained at LEP and by the ATLAS experiment in previous searches are also shown in the shaded areas.

Observed and expected exclusion limits on SUSY simplified models for direct selectron production in the (a) $m(\tilde{e})-m(\tilde{\chi}_1^0)$ and (c) $m(\tilde{e})-\Delta m(\tilde{e},\tilde{\chi}_1^0)$ planes, and for direct smuon production in the (b) $m(\tilde{\mu})-m(\tilde{\chi}_1^0)$ and (d) $m(\tilde{\mu})-\Delta m(\tilde{\mu},\tilde{\chi}_1^0)$ planes. In Figure (a) and (c) the observed (solid thick lines) and expected (dashed lines) exclusion contours are indicated for combined $\tilde{e}_{\textup{L,R}}$ and for $\tilde{e}_{\textup{L}}$ and $\tilde{e}_{\textup{R}}$. In Figure (b) and (d) the observed (solid thick lines) and expected (dashed lines) exclusion contours are indicated for combined $\tilde{\mu}_{\textup{L,R}}$ and for $\tilde{\mu}_{\textup{L}}$. No unique sensitivity to $\tilde{\mu}_{\textup{R}}$ is observed. All limits are computed at 95% CL. The observed limits obtained at LEP and by the ATLAS experiment in previous searches are also shown in the shaded areas.

Observed and expected exclusion limits on SUSY simplified models for direct selectron production in the (a) $m(\tilde{e})-m(\tilde{\chi}_1^0)$ and (c) $m(\tilde{e})-\Delta m(\tilde{e},\tilde{\chi}_1^0)$ planes, and for direct smuon production in the (b) $m(\tilde{\mu})-m(\tilde{\chi}_1^0)$ and (d) $m(\tilde{\mu})-\Delta m(\tilde{\mu},\tilde{\chi}_1^0)$ planes. In Figure (a) and (c) the observed (solid thick lines) and expected (dashed lines) exclusion contours are indicated for combined $\tilde{e}_{\textup{L,R}}$ and for $\tilde{e}_{\textup{L}}$ and $\tilde{e}_{\textup{R}}$. In Figure (b) and (d) the observed (solid thick lines) and expected (dashed lines) exclusion contours are indicated for combined $\tilde{\mu}_{\textup{L,R}}$ and for $\tilde{\mu}_{\textup{L}}$. No unique sensitivity to $\tilde{\mu}_{\textup{R}}$ is observed. All limits are computed at 95% CL. The observed limits obtained at LEP and by the ATLAS experiment in previous searches are also shown in the shaded areas.

Observed and expected exclusion limits on SUSY simplified models for direct selectron production in the (a) $m(\tilde{e})-m(\tilde{\chi}_1^0)$ and (c) $m(\tilde{e})-\Delta m(\tilde{e},\tilde{\chi}_1^0)$ planes, and for direct smuon production in the (b) $m(\tilde{\mu})-m(\tilde{\chi}_1^0)$ and (d) $m(\tilde{\mu})-\Delta m(\tilde{\mu},\tilde{\chi}_1^0)$ planes. In Figure (a) and (c) the observed (solid thick lines) and expected (dashed lines) exclusion contours are indicated for combined $\tilde{e}_{\textup{L,R}}$ and for $\tilde{e}_{\textup{L}}$ and $\tilde{e}_{\textup{R}}$. In Figure (b) and (d) the observed (solid thick lines) and expected (dashed lines) exclusion contours are indicated for combined $\tilde{\mu}_{\textup{L,R}}$ and for $\tilde{\mu}_{\textup{L}}$. No unique sensitivity to $\tilde{\mu}_{\textup{R}}$ is observed. All limits are computed at 95% CL. The observed limits obtained at LEP and by the ATLAS experiment in previous searches are also shown in the shaded areas.

Observed and expected exclusion limits on SUSY simplified models for direct selectron production in the (a) $m(\tilde{e})-m(\tilde{\chi}_1^0)$ and (c) $m(\tilde{e})-\Delta m(\tilde{e},\tilde{\chi}_1^0)$ planes, and for direct smuon production in the (b) $m(\tilde{\mu})-m(\tilde{\chi}_1^0)$ and (d) $m(\tilde{\mu})-\Delta m(\tilde{\mu},\tilde{\chi}_1^0)$ planes. In Figure (a) and (c) the observed (solid thick lines) and expected (dashed lines) exclusion contours are indicated for combined $\tilde{e}_{\textup{L,R}}$ and for $\tilde{e}_{\textup{L}}$ and $\tilde{e}_{\textup{R}}$. In Figure (b) and (d) the observed (solid thick lines) and expected (dashed lines) exclusion contours are indicated for combined $\tilde{\mu}_{\textup{L,R}}$ and for $\tilde{\mu}_{\textup{L}}$. No unique sensitivity to $\tilde{\mu}_{\textup{R}}$ is observed. All limits are computed at 95% CL. The observed limits obtained at LEP and by the ATLAS experiment in previous searches are also shown in the shaded areas.

Observed and expected exclusion limits on SUSY simplified models for direct selectron production in the (a) $m(\tilde{e})-m(\tilde{\chi}_1^0)$ and (c) $m(\tilde{e})-\Delta m(\tilde{e},\tilde{\chi}_1^0)$ planes, and for direct smuon production in the (b) $m(\tilde{\mu})-m(\tilde{\chi}_1^0)$ and (d) $m(\tilde{\mu})-\Delta m(\tilde{\mu},\tilde{\chi}_1^0)$ planes. In Figure (a) and (c) the observed (solid thick lines) and expected (dashed lines) exclusion contours are indicated for combined $\tilde{e}_{\textup{L,R}}$ and for $\tilde{e}_{\textup{L}}$ and $\tilde{e}_{\textup{R}}$. In Figure (b) and (d) the observed (solid thick lines) and expected (dashed lines) exclusion contours are indicated for combined $\tilde{\mu}_{\textup{L,R}}$ and for $\tilde{\mu}_{\textup{L}}$. No unique sensitivity to $\tilde{\mu}_{\textup{R}}$ is observed. All limits are computed at 95% CL. The observed limits obtained at LEP and by the ATLAS experiment in previous searches are also shown in the shaded areas.

Observed and expected exclusion limits on SUSY simplified models for direct selectron production in the (a) $m(\tilde{e})-m(\tilde{\chi}_1^0)$ and (c) $m(\tilde{e})-\Delta m(\tilde{e},\tilde{\chi}_1^0)$ planes, and for direct smuon production in the (b) $m(\tilde{\mu})-m(\tilde{\chi}_1^0)$ and (d) $m(\tilde{\mu})-\Delta m(\tilde{\mu},\tilde{\chi}_1^0)$ planes. In Figure (a) and (c) the observed (solid thick lines) and expected (dashed lines) exclusion contours are indicated for combined $\tilde{e}_{\textup{L,R}}$ and for $\tilde{e}_{\textup{L}}$ and $\tilde{e}_{\textup{R}}$. In Figure (b) and (d) the observed (solid thick lines) and expected (dashed lines) exclusion contours are indicated for combined $\tilde{\mu}_{\textup{L,R}}$ and for $\tilde{\mu}_{\textup{L}}$. No unique sensitivity to $\tilde{\mu}_{\textup{R}}$ is observed. All limits are computed at 95% CL. The observed limits obtained at LEP and by the ATLAS experiment in previous searches are also shown in the shaded areas.

Observed and expected exclusion limits on SUSY simplified models for direct selectron production in the (a) $m(\tilde{e})-m(\tilde{\chi}_1^0)$ and (c) $m(\tilde{e})-\Delta m(\tilde{e},\tilde{\chi}_1^0)$ planes, and for direct smuon production in the (b) $m(\tilde{\mu})-m(\tilde{\chi}_1^0)$ and (d) $m(\tilde{\mu})-\Delta m(\tilde{\mu},\tilde{\chi}_1^0)$ planes. In Figure (a) and (c) the observed (solid thick lines) and expected (dashed lines) exclusion contours are indicated for combined $\tilde{e}_{\textup{L,R}}$ and for $\tilde{e}_{\textup{L}}$ and $\tilde{e}_{\textup{R}}$. In Figure (b) and (d) the observed (solid thick lines) and expected (dashed lines) exclusion contours are indicated for combined $\tilde{\mu}_{\textup{L,R}}$ and for $\tilde{\mu}_{\textup{L}}$. No unique sensitivity to $\tilde{\mu}_{\textup{R}}$ is observed. All limits are computed at 95% CL. The observed limits obtained at LEP and by the ATLAS experiment in previous searches are also shown in the shaded areas.

Observed and expected exclusion limits on SUSY simplified models for direct selectron production in the (a) $m(\tilde{e})-m(\tilde{\chi}_1^0)$ and (c) $m(\tilde{e})-\Delta m(\tilde{e},\tilde{\chi}_1^0)$ planes, and for direct smuon production in the (b) $m(\tilde{\mu})-m(\tilde{\chi}_1^0)$ and (d) $m(\tilde{\mu})-\Delta m(\tilde{\mu},\tilde{\chi}_1^0)$ planes. In Figure (a) and (c) the observed (solid thick lines) and expected (dashed lines) exclusion contours are indicated for combined $\tilde{e}_{\textup{L,R}}$ and for $\tilde{e}_{\textup{L}}$ and $\tilde{e}_{\textup{R}}$. In Figure (b) and (d) the observed (solid thick lines) and expected (dashed lines) exclusion contours are indicated for combined $\tilde{\mu}_{\textup{L,R}}$ and for $\tilde{\mu}_{\textup{L}}$. No unique sensitivity to $\tilde{\mu}_{\textup{R}}$ is observed. All limits are computed at 95% CL. The observed limits obtained at LEP and by the ATLAS experiment in previous searches are also shown in the shaded areas.

Observed and expected exclusion limits on SUSY simplified models for slepton-pair production in the $m(\tilde{\ell})-m(\tilde{\chi}_1^0)$ plane. Only $\tilde{e}$ and $\tilde{\mu}$ are considered. The observed (solid thick line) and expected (thin dashed line) exclusion contours are indicated. The red contour shows the exclusion limits obtained using both the SR-0J and SR-1J region, as presented in Figure 6. The blue and green contours correspond to the result obtained considering only SR-0J and SR-1J region respectively. All limits are computed at 95% CL. The observed limits obtained by the ATLAS experiment in previous searches are also shown.

Observed and expected exclusion limits on SUSY simplified models for slepton-pair production in the $m(\tilde{\ell})-m(\tilde{\chi}_1^0)$ plane. Only $\tilde{e}$ and $\tilde{\mu}$ are considered. The observed (solid thick line) and expected (thin dashed line) exclusion contours are indicated. The red contour shows the exclusion limits obtained using both the SR-0J and SR-1J region, as presented in Figure 6. The blue and green contours correspond to the result obtained considering only SR-0J and SR-1J region respectively. All limits are computed at 95% CL. The observed limits obtained by the ATLAS experiment in previous searches are also shown.

Observed and expected exclusion limits on SUSY simplified models for slepton-pair production in the $m(\tilde{\ell})-m(\tilde{\chi}_1^0)$ plane. Only $\tilde{e}$ and $\tilde{\mu}$ are considered. The observed (solid thick line) and expected (thin dashed line) exclusion contours are indicated. The red contour shows the exclusion limits obtained using both the SR-0J and SR-1J region, as presented in Figure 6. The blue and green contours correspond to the result obtained considering only SR-0J and SR-1J region respectively. All limits are computed at 95% CL. The observed limits obtained by the ATLAS experiment in previous searches are also shown.

Observed and expected exclusion limits on SUSY simplified models for slepton-pair production in the $m(\tilde{\ell})-m(\tilde{\chi}_1^0)$ plane. Only $\tilde{e}$ and $\tilde{\mu}$ are considered. The observed (solid thick line) and expected (thin dashed line) exclusion contours are indicated. The red contour shows the exclusion limits obtained using both the SR-0J and SR-1J region, as presented in Figure 6. The blue and green contours correspond to the result obtained considering only SR-0J and SR-1J region respectively. All limits are computed at 95% CL. The observed limits obtained by the ATLAS experiment in previous searches are also shown.

The upper panel shows the observed number of events in each of the binned SRs defined in Table 3, together with the expected SM backgrounds obtained after applying the efficiency correction method to compute the number of expected FSB events. `Others' include the non-dominant background sources, e.g. $t \bar{t}$+$V$, Higgs boson and Drell--Yan events. The uncertainty band includes systematic and statistical errors from all sources. The distributions of two signal points with mass splittings $\Delta m(\tilde{\ell},\tilde{\chi}_1^0) = m(\tilde{\ell})-m(\tilde{\chi}_1^0) = 30$ GeV and $\Delta m(\tilde{\ell},\tilde{\chi}_1^0) = m(\tilde{\ell})-m(\tilde{\chi}_1^0) = 50$ GeV are overlaid. The lower panel shows the significance as defined in Ref. [115].

The figure shows the signal acceptance (a) and efficiency (b) plots for the $\tilde{\chi}_1^+\tilde{\chi}_1^-$ production with $W$-boson-mediated decay model, in the SR$^{\text{-DF BDT-signal}\in(0.81,1]}_{\text{-SF BDT-signal}\in(0.77,1]}$ inclusive region. Acceptance is calculated by applying the signal region requirements to particle-level objects, which do not suffer from identification inefficiencies or mismeasurements. The efficiency is calculated with fully reconstructed objects with the acceptance divided out. Large acceptance and efficiency differences in neighbouring points are due to statistical fluctuations.

The figure shows the signal acceptance (a) and efficiency (b) plots for the $\tilde{\chi}_1^+\tilde{\chi}_1^-$ production with $W$-boson-mediated decay model, in the SR$^{\text{-DF BDT-signal}\in(0.81,1]}_{\text{-SF BDT-signal}\in(0.77,1]}$ inclusive region. Acceptance is calculated by applying the signal region requirements to particle-level objects, which do not suffer from identification inefficiencies or mismeasurements. The efficiency is calculated with fully reconstructed objects with the acceptance divided out. Large acceptance and efficiency differences in neighbouring points are due to statistical fluctuations.

The figure shows the signal acceptance (a) and efficiency (b) plots for the $\tilde{\chi}_1^+\tilde{\chi}_1^-$ production with $W$-boson-mediated decay model, in the SR-DF BDT-signal$\in(0.81,1]$ inclusive region. Acceptance is calculated by applying the signal region requirements to particle-level objects, which do not suffer from identification inefficiencies or mismeasurements. The efficiency is calculated with fully reconstructed objects with the acceptance divided out. Large acceptance and efficiency differences in neighbouring points are due to statistical fluctuations.

The figure shows the signal acceptance (a) and efficiency (b) plots for the $\tilde{\chi}_1^+\tilde{\chi}_1^-$ production with $W$-boson-mediated decay model, in the SR-DF BDT-signal$\in(0.81,1]$ inclusive region. Acceptance is calculated by applying the signal region requirements to particle-level objects, which do not suffer from identification inefficiencies or mismeasurements. The efficiency is calculated with fully reconstructed objects with the acceptance divided out. Large acceptance and efficiency differences in neighbouring points are due to statistical fluctuations.

The figure shows the signal acceptance (a) and efficiency (b) plots for the $\tilde{\chi}_1^+\tilde{\chi}_1^-$ production with $W$-boson-mediated decay model, in the SR-DF BDT-signal$\in(0.82,1]$ inclusive region. Acceptance is calculated by applying the signal region requirements to particle-level objects, which do not suffer from identification inefficiencies or mismeasurements. The efficiency is calculated with fully reconstructed objects with the acceptance divided out. Large acceptance and efficiency differences in neighbouring points are due to statistical fluctuations.

The figure shows the signal acceptance (a) and efficiency (b) plots for the $\tilde{\chi}_1^+\tilde{\chi}_1^-$ production with $W$-boson-mediated decay model, in the SR-DF BDT-signal$\in(0.82,1]$ inclusive region. Acceptance is calculated by applying the signal region requirements to particle-level objects, which do not suffer from identification inefficiencies or mismeasurements. The efficiency is calculated with fully reconstructed objects with the acceptance divided out. Large acceptance and efficiency differences in neighbouring points are due to statistical fluctuations.

The figure shows the signal acceptance (a) and efficiency (b) plots for the $\tilde{\chi}_1^+\tilde{\chi}_1^-$ production with $W$-boson-mediated decay model, in the SR-DF BDT-signal$\in(0.83,1]$ inclusive region. Acceptance is calculated by applying the signal region requirements to particle-level objects, which do not suffer from identification inefficiencies or mismeasurements. The efficiency is calculated with fully reconstructed objects with the acceptance divided out. Large acceptance and efficiency differences in neighbouring points are due to statistical fluctuations.

The figure shows the signal acceptance (a) and efficiency (b) plots for the $\tilde{\chi}_1^+\tilde{\chi}_1^-$ production with $W$-boson-mediated decay model, in the SR-DF BDT-signal$\in(0.83,1]$ inclusive region. Acceptance is calculated by applying the signal region requirements to particle-level objects, which do not suffer from identification inefficiencies or mismeasurements. The efficiency is calculated with fully reconstructed objects with the acceptance divided out. Large acceptance and efficiency differences in neighbouring points are due to statistical fluctuations.

The figure shows the signal acceptance (a) and efficiency (b) plots for the $\tilde{\chi}_1^+\tilde{\chi}_1^-$ production with $W$-boson-mediated decay model, in the SR-DF BDT-signal$\in(0.84,1]$ inclusive region. Acceptance is calculated by applying the signal region requirements to particle-level objects, which do not suffer from identification inefficiencies or mismeasurements. The efficiency is calculated with fully reconstructed objects with the acceptance divided out. Large acceptance and efficiency differences in neighbouring points are due to statistical fluctuations.

The figure shows the signal acceptance (a) and efficiency (b) plots for the $\tilde{\chi}_1^+\tilde{\chi}_1^-$ production with $W$-boson-mediated decay model, in the SR-DF BDT-signal$\in(0.84,1]$ inclusive region. Acceptance is calculated by applying the signal region requirements to particle-level objects, which do not suffer from identification inefficiencies or mismeasurements. The efficiency is calculated with fully reconstructed objects with the acceptance divided out. Large acceptance and efficiency differences in neighbouring points are due to statistical fluctuations.

The figure shows the signal acceptance (a) and efficiency (b) plots for the $\tilde{\chi}_1^+\tilde{\chi}_1^-$ production with $W$-boson-mediated decay model, in the SR-DF BDT-signal$\in(0.85,1]$ inclusive region. Acceptance is calculated by applying the signal region requirements to particle-level objects, which do not suffer from identification inefficiencies or mismeasurements. The efficiency is calculated with fully reconstructed objects with the acceptance divided out. Large acceptance and efficiency differences in neighbouring points are due to statistical fluctuations.

The figure shows the signal acceptance (a) and efficiency (b) plots for the $\tilde{\chi}_1^+\tilde{\chi}_1^-$ production with $W$-boson-mediated decay model, in the SR-DF BDT-signal$\in(0.85,1]$ inclusive region. Acceptance is calculated by applying the signal region requirements to particle-level objects, which do not suffer from identification inefficiencies or mismeasurements. The efficiency is calculated with fully reconstructed objects with the acceptance divided out. Large acceptance and efficiency differences in neighbouring points are due to statistical fluctuations.

The figure shows the signal acceptance (a) and efficiency (b) plots for the $\tilde{\chi}_1^+\tilde{\chi}_1^-$ production with $W$-boson-mediated decay model, in the SR-DF BDT-signal$\in(0.81,0.8125]$ inclusive region. Acceptance is calculated by applying the signal region requirements to particle-level objects, which do not suffer from identification inefficiencies or mismeasurements. The efficiency is calculated with fully reconstructed objects with the acceptance divided out. Large acceptance and efficiency differences in neighbouring points are due to statistical fluctuations.

The figure shows the signal acceptance (a) and efficiency (b) plots for the $\tilde{\chi}_1^+\tilde{\chi}_1^-$ production with $W$-boson-mediated decay model, in the SR-DF BDT-signal$\in(0.81,0.8125]$ inclusive region. Acceptance is calculated by applying the signal region requirements to particle-level objects, which do not suffer from identification inefficiencies or mismeasurements. The efficiency is calculated with fully reconstructed objects with the acceptance divided out. Large acceptance and efficiency differences in neighbouring points are due to statistical fluctuations.

The figure shows the signal acceptance (a) and efficiency (b) plots for the $\tilde{\chi}_1^+\tilde{\chi}_1^-$ production with $W$-boson-mediated decay model, in the SR-DF BDT-signal$\in(0.8125,0.815]$ inclusive region. Acceptance is calculated by applying the signal region requirements to particle-level objects, which do not suffer from identification inefficiencies or mismeasurements. The efficiency is calculated with fully reconstructed objects with the acceptance divided out. Large acceptance and efficiency differences in neighbouring points are due to statistical fluctuations.

The figure shows the signal acceptance (a) and efficiency (b) plots for the $\tilde{\chi}_1^+\tilde{\chi}_1^-$ production with $W$-boson-mediated decay model, in the SR-DF BDT-signal$\in(0.8125,0.815]$ inclusive region. Acceptance is calculated by applying the signal region requirements to particle-level objects, which do not suffer from identification inefficiencies or mismeasurements. The efficiency is calculated with fully reconstructed objects with the acceptance divided out. Large acceptance and efficiency differences in neighbouring points are due to statistical fluctuations.

The figure shows the signal acceptance (a) and efficiency (b) plots for the $\tilde{\chi}_1^+\tilde{\chi}_1^-$ production with $W$-boson-mediated decay model, in the SR-DF BDT-signal$\in(0.815,0.8175]$ inclusive region. Acceptance is calculated by applying the signal region requirements to particle-level objects, which do not suffer from identification inefficiencies or mismeasurements. The efficiency is calculated with fully reconstructed objects with the acceptance divided out. Large acceptance and efficiency differences in neighbouring points are due to statistical fluctuations.

The figure shows the signal acceptance (a) and efficiency (b) plots for the $\tilde{\chi}_1^+\tilde{\chi}_1^-$ production with $W$-boson-mediated decay model, in the SR-DF BDT-signal$\in(0.815,0.8175]$ inclusive region. Acceptance is calculated by applying the signal region requirements to particle-level objects, which do not suffer from identification inefficiencies or mismeasurements. The efficiency is calculated with fully reconstructed objects with the acceptance divided out. Large acceptance and efficiency differences in neighbouring points are due to statistical fluctuations.

The figure shows the signal acceptance (a) and efficiency (b) plots for the $\tilde{\chi}_1^+\tilde{\chi}_1^-$ production with $W$-boson-mediated decay model, in the SR-DF BDT-signal$\in(0.8175,0.82]$ inclusive region. Acceptance is calculated by applying the signal region requirements to particle-level objects, which do not suffer from identification inefficiencies or mismeasurements. The efficiency is calculated with fully reconstructed objects with the acceptance divided out. Large acceptance and efficiency differences in neighbouring points are due to statistical fluctuations.

The figure shows the signal acceptance (a) and efficiency (b) plots for the $\tilde{\chi}_1^+\tilde{\chi}_1^-$ production with $W$-boson-mediated decay model, in the SR-DF BDT-signal$\in(0.8175,0.82]$ inclusive region. Acceptance is calculated by applying the signal region requirements to particle-level objects, which do not suffer from identification inefficiencies or mismeasurements. The efficiency is calculated with fully reconstructed objects with the acceptance divided out. Large acceptance and efficiency differences in neighbouring points are due to statistical fluctuations.

The figure shows the signal acceptance (a) and efficiency (b) plots for the $\tilde{\chi}_1^+\tilde{\chi}_1^-$ production with $W$-boson-mediated decay model, in the SR-DF BDT-signal$\in(0.82,0.8225]$ inclusive region. Acceptance is calculated by applying the signal region requirements to particle-level objects, which do not suffer from identification inefficiencies or mismeasurements. The efficiency is calculated with fully reconstructed objects with the acceptance divided out. Large acceptance and efficiency differences in neighbouring points are due to statistical fluctuations.

The figure shows the signal acceptance (a) and efficiency (b) plots for the $\tilde{\chi}_1^+\tilde{\chi}_1^-$ production with $W$-boson-mediated decay model, in the SR-DF BDT-signal$\in(0.82,0.8225]$ inclusive region. Acceptance is calculated by applying the signal region requirements to particle-level objects, which do not suffer from identification inefficiencies or mismeasurements. The efficiency is calculated with fully reconstructed objects with the acceptance divided out. Large acceptance and efficiency differences in neighbouring points are due to statistical fluctuations.

The figure shows the signal acceptance (a) and efficiency (b) plots for the $\tilde{\chi}_1^+\tilde{\chi}_1^-$ production with $W$-boson-mediated decay model, in the SR-DF BDT-signal$\in(0.8225,0.825]$ inclusive region. Acceptance is calculated by applying the signal region requirements to particle-level objects, which do not suffer from identification inefficiencies or mismeasurements. The efficiency is calculated with fully reconstructed objects with the acceptance divided out. Large acceptance and efficiency differences in neighbouring points are due to statistical fluctuations.

The figure shows the signal acceptance (a) and efficiency (b) plots for the $\tilde{\chi}_1^+\tilde{\chi}_1^-$ production with $W$-boson-mediated decay model, in the SR-DF BDT-signal$\in(0.8225,0.825]$ inclusive region. Acceptance is calculated by applying the signal region requirements to particle-level objects, which do not suffer from identification inefficiencies or mismeasurements. The efficiency is calculated with fully reconstructed objects with the acceptance divided out. Large acceptance and efficiency differences in neighbouring points are due to statistical fluctuations.

The figure shows the signal acceptance (a) and efficiency (b) plots for the $\tilde{\chi}_1^+\tilde{\chi}_1^-$ production with $W$-boson-mediated decay model, in the SR-DF BDT-signal$\in(0.825,0.8275]$ inclusive region. Acceptance is calculated by applying the signal region requirements to particle-level objects, which do not suffer from identification inefficiencies or mismeasurements. The efficiency is calculated with fully reconstructed objects with the acceptance divided out. Large acceptance and efficiency differences in neighbouring points are due to statistical fluctuations.

The figure shows the signal acceptance (a) and efficiency (b) plots for the $\tilde{\chi}_1^+\tilde{\chi}_1^-$ production with $W$-boson-mediated decay model, in the SR-DF BDT-signal$\in(0.825,0.8275]$ inclusive region. Acceptance is calculated by applying the signal region requirements to particle-level objects, which do not suffer from identification inefficiencies or mismeasurements. The efficiency is calculated with fully reconstructed objects with the acceptance divided out. Large acceptance and efficiency differences in neighbouring points are due to statistical fluctuations.

The figure shows the signal acceptance (a) and efficiency (b) plots for the $\tilde{\chi}_1^+\tilde{\chi}_1^-$ production with $W$-boson-mediated decay model, in the SR-DF BDT-signal$\in(0.8275,0.83]$ inclusive region. Acceptance is calculated by applying the signal region requirements to particle-level objects, which do not suffer from identification inefficiencies or mismeasurements. The efficiency is calculated with fully reconstructed objects with the acceptance divided out. Large acceptance and efficiency differences in neighbouring points are due to statistical fluctuations.

The figure shows the signal acceptance (a) and efficiency (b) plots for the $\tilde{\chi}_1^+\tilde{\chi}_1^-$ production with $W$-boson-mediated decay model, in the SR-DF BDT-signal$\in(0.8275,0.83]$ inclusive region. Acceptance is calculated by applying the signal region requirements to particle-level objects, which do not suffer from identification inefficiencies or mismeasurements. The efficiency is calculated with fully reconstructed objects with the acceptance divided out. Large acceptance and efficiency differences in neighbouring points are due to statistical fluctuations.

The figure shows the signal acceptance (a) and efficiency (b) plots for the $\tilde{\chi}_1^+\tilde{\chi}_1^-$ production with $W$-boson-mediated decay model, in the SR-DF BDT-signal$\in(0.83,0.8325]$ inclusive region. Acceptance is calculated by applying the signal region requirements to particle-level objects, which do not suffer from identification inefficiencies or mismeasurements. The efficiency is calculated with fully reconstructed objects with the acceptance divided out. Large acceptance and efficiency differences in neighbouring points are due to statistical fluctuations.

The figure shows the signal acceptance (a) and efficiency (b) plots for the $\tilde{\chi}_1^+\tilde{\chi}_1^-$ production with $W$-boson-mediated decay model, in the SR-DF BDT-signal$\in(0.83,0.8325]$ inclusive region. Acceptance is calculated by applying the signal region requirements to particle-level objects, which do not suffer from identification inefficiencies or mismeasurements. The efficiency is calculated with fully reconstructed objects with the acceptance divided out. Large acceptance and efficiency differences in neighbouring points are due to statistical fluctuations.

The figure shows the signal acceptance (a) and efficiency (b) plots for the $\tilde{\chi}_1^+\tilde{\chi}_1^-$ production with $W$-boson-mediated decay model, in the SR-DF BDT-signal$\in(0.8325,0.835]$ inclusive region. Acceptance is calculated by applying the signal region requirements to particle-level objects, which do not suffer from identification inefficiencies or mismeasurements. The efficiency is calculated with fully reconstructed objects with the acceptance divided out. Large acceptance and efficiency differences in neighbouring points are due to statistical fluctuations.

The figure shows the signal acceptance (a) and efficiency (b) plots for the $\tilde{\chi}_1^+\tilde{\chi}_1^-$ production with $W$-boson-mediated decay model, in the SR-DF BDT-signal$\in(0.8325,0.835]$ inclusive region. Acceptance is calculated by applying the signal region requirements to particle-level objects, which do not suffer from identification inefficiencies or mismeasurements. The efficiency is calculated with fully reconstructed objects with the acceptance divided out. Large acceptance and efficiency differences in neighbouring points are due to statistical fluctuations.

The figure shows the signal acceptance (a) and efficiency (b) plots for the $\tilde{\chi}_1^+\tilde{\chi}_1^-$ production with $W$-boson-mediated decay model, in the SR-DF BDT-signal$\in(0.835,0.8375]$ inclusive region. Acceptance is calculated by applying the signal region requirements to particle-level objects, which do not suffer from identification inefficiencies or mismeasurements. The efficiency is calculated with fully reconstructed objects with the acceptance divided out. Large acceptance and efficiency differences in neighbouring points are due to statistical fluctuations.

The figure shows the signal acceptance (a) and efficiency (b) plots for the $\tilde{\chi}_1^+\tilde{\chi}_1^-$ production with $W$-boson-mediated decay model, in the SR-DF BDT-signal$\in(0.835,0.8375]$ inclusive region. Acceptance is calculated by applying the signal region requirements to particle-level objects, which do not suffer from identification inefficiencies or mismeasurements. The efficiency is calculated with fully reconstructed objects with the acceptance divided out. Large acceptance and efficiency differences in neighbouring points are due to statistical fluctuations.

The figure shows the signal acceptance (a) and efficiency (b) plots for the $\tilde{\chi}_1^+\tilde{\chi}_1^-$ production with $W$-boson-mediated decay model, in the SR-DF BDT-signal$\in(0.8375,0.84]$ inclusive region. Acceptance is calculated by applying the signal region requirements to particle-level objects, which do not suffer from identification inefficiencies or mismeasurements. The efficiency is calculated with fully reconstructed objects with the acceptance divided out. Large acceptance and efficiency differences in neighbouring points are due to statistical fluctuations.

The figure shows the signal acceptance (a) and efficiency (b) plots for the $\tilde{\chi}_1^+\tilde{\chi}_1^-$ production with $W$-boson-mediated decay model, in the SR-DF BDT-signal$\in(0.8375,0.84]$ inclusive region. Acceptance is calculated by applying the signal region requirements to particle-level objects, which do not suffer from identification inefficiencies or mismeasurements. The efficiency is calculated with fully reconstructed objects with the acceptance divided out. Large acceptance and efficiency differences in neighbouring points are due to statistical fluctuations.

The figure shows the signal acceptance (a) and efficiency (b) plots for the $\tilde{\chi}_1^+\tilde{\chi}_1^-$ production with $W$-boson-mediated decay model, in the SR-DF BDT-signal$\in(0.84,0.845]$ inclusive region. Acceptance is calculated by applying the signal region requirements to particle-level objects, which do not suffer from identification inefficiencies or mismeasurements. The efficiency is calculated with fully reconstructed objects with the acceptance divided out. Large acceptance and efficiency differences in neighbouring points are due to statistical fluctuations.

The figure shows the signal acceptance (a) and efficiency (b) plots for the $\tilde{\chi}_1^+\tilde{\chi}_1^-$ production with $W$-boson-mediated decay model, in the SR-DF BDT-signal$\in(0.84,0.845]$ inclusive region. Acceptance is calculated by applying the signal region requirements to particle-level objects, which do not suffer from identification inefficiencies or mismeasurements. The efficiency is calculated with fully reconstructed objects with the acceptance divided out. Large acceptance and efficiency differences in neighbouring points are due to statistical fluctuations.

The figure shows the signal acceptance (a) and efficiency (b) plots for the $\tilde{\chi}_1^+\tilde{\chi}_1^-$ production with $W$-boson-mediated decay model, in the SR-DF BDT-signal$\in(0.845,0.85]$ inclusive region. Acceptance is calculated by applying the signal region requirements to particle-level objects, which do not suffer from identification inefficiencies or mismeasurements. The efficiency is calculated with fully reconstructed objects with the acceptance divided out. Large acceptance and efficiency differences in neighbouring points are due to statistical fluctuations.

The figure shows the signal acceptance (a) and efficiency (b) plots for the $\tilde{\chi}_1^+\tilde{\chi}_1^-$ production with $W$-boson-mediated decay model, in the SR-DF BDT-signal$\in(0.845,0.85]$ inclusive region. Acceptance is calculated by applying the signal region requirements to particle-level objects, which do not suffer from identification inefficiencies or mismeasurements. The efficiency is calculated with fully reconstructed objects with the acceptance divided out. Large acceptance and efficiency differences in neighbouring points are due to statistical fluctuations.

The figure shows the signal acceptance (a) and efficiency (b) plots for the $\tilde{\chi}_1^+\tilde{\chi}_1^-$ production with $W$-boson-mediated decay model, in the SR-DF BDT-signal$\in(0.85,0.86]$ inclusive region. Acceptance is calculated by applying the signal region requirements to particle-level objects, which do not suffer from identification inefficiencies or mismeasurements. The efficiency is calculated with fully reconstructed objects with the acceptance divided out. Large acceptance and efficiency differences in neighbouring points are due to statistical fluctuations.

The figure shows the signal acceptance (a) and efficiency (b) plots for the $\tilde{\chi}_1^+\tilde{\chi}_1^-$ production with $W$-boson-mediated decay model, in the SR-DF BDT-signal$\in(0.85,0.86]$ inclusive region. Acceptance is calculated by applying the signal region requirements to particle-level objects, which do not suffer from identification inefficiencies or mismeasurements. The efficiency is calculated with fully reconstructed objects with the acceptance divided out. Large acceptance and efficiency differences in neighbouring points are due to statistical fluctuations.

The figure shows the signal acceptance (a) and efficiency (b) plots for the $\tilde{\chi}_1^+\tilde{\chi}_1^-$ production with $W$-boson-mediated decay model, in the SR-DF BDT-signal$\in(0.86,1]$ inclusive region. Acceptance is calculated by applying the signal region requirements to particle-level objects, which do not suffer from identification inefficiencies or mismeasurements. The efficiency is calculated with fully reconstructed objects with the acceptance divided out. Large acceptance and efficiency differences in neighbouring points are due to statistical fluctuations.

The figure shows the signal acceptance (a) and efficiency (b) plots for the $\tilde{\chi}_1^+\tilde{\chi}_1^-$ production with $W$-boson-mediated decay model, in the SR-DF BDT-signal$\in(0.86,1]$ inclusive region. Acceptance is calculated by applying the signal region requirements to particle-level objects, which do not suffer from identification inefficiencies or mismeasurements. The efficiency is calculated with fully reconstructed objects with the acceptance divided out. Large acceptance and efficiency differences in neighbouring points are due to statistical fluctuations.

The figure shows the signal acceptance (a) and efficiency (b) plots for the $\tilde{\chi}_1^+\tilde{\chi}_1^-$ production with $W$-boson-mediated decay model, in the SR-SF BDT-signal$\in(0.77,1]$ inclusive region. Acceptance is calculated by applying the signal region requirements to particle-level objects, which do not suffer from identification inefficiencies or mismeasurements. The efficiency is calculated with fully reconstructed objects with the acceptance divided out. Large acceptance and efficiency differences in neighbouring points are due to statistical fluctuations.

The figure shows the signal acceptance (a) and efficiency (b) plots for the $\tilde{\chi}_1^+\tilde{\chi}_1^-$ production with $W$-boson-mediated decay model, in the SR-SF BDT-signal$\in(0.77,1]$ inclusive region. Acceptance is calculated by applying the signal region requirements to particle-level objects, which do not suffer from identification inefficiencies or mismeasurements. The efficiency is calculated with fully reconstructed objects with the acceptance divided out. Large acceptance and efficiency differences in neighbouring points are due to statistical fluctuations.

The figure shows the signal acceptance (a) and efficiency (b) plots for the $\tilde{\chi}_1^+\tilde{\chi}_1^-$ production with $W$-boson-mediated decay model, in the SR-SF BDT-signal$\in(0.78,1]$ inclusive region. Acceptance is calculated by applying the signal region requirements to particle-level objects, which do not suffer from identification inefficiencies or mismeasurements. The efficiency is calculated with fully reconstructed objects with the acceptance divided out. Large acceptance and efficiency differences in neighbouring points are due to statistical fluctuations.

The figure shows the signal acceptance (a) and efficiency (b) plots for the $\tilde{\chi}_1^+\tilde{\chi}_1^-$ production with $W$-boson-mediated decay model, in the SR-SF BDT-signal$\in(0.78,1]$ inclusive region. Acceptance is calculated by applying the signal region requirements to particle-level objects, which do not suffer from identification inefficiencies or mismeasurements. The efficiency is calculated with fully reconstructed objects with the acceptance divided out. Large acceptance and efficiency differences in neighbouring points are due to statistical fluctuations.

The figure shows the signal acceptance (a) and efficiency (b) plots for the $\tilde{\chi}_1^+\tilde{\chi}_1^-$ production with $W$-boson-mediated decay model, in the SR-SF BDT-signal$\in(0.79,1]$ inclusive region. Acceptance is calculated by applying the signal region requirements to particle-level objects, which do not suffer from identification inefficiencies or mismeasurements. The efficiency is calculated with fully reconstructed objects with the acceptance divided out. Large acceptance and efficiency differences in neighbouring points are due to statistical fluctuations.

The figure shows the signal acceptance (a) and efficiency (b) plots for the $\tilde{\chi}_1^+\tilde{\chi}_1^-$ production with $W$-boson-mediated decay model, in the SR-SF BDT-signal$\in(0.79,1]$ inclusive region. Acceptance is calculated by applying the signal region requirements to particle-level objects, which do not suffer from identification inefficiencies or mismeasurements. The efficiency is calculated with fully reconstructed objects with the acceptance divided out. Large acceptance and efficiency differences in neighbouring points are due to statistical fluctuations.

The figure shows the signal acceptance (a) and efficiency (b) plots for the $\tilde{\chi}_1^+\tilde{\chi}_1^-$ production with $W$-boson-mediated decay model, in the SR-SF BDT-signal$\in(0.80,1]$ inclusive region. Acceptance is calculated by applying the signal region requirements to particle-level objects, which do not suffer from identification inefficiencies or mismeasurements. The efficiency is calculated with fully reconstructed objects with the acceptance divided out. Large acceptance and efficiency differences in neighbouring points are due to statistical fluctuations.

The figure shows the signal acceptance (a) and efficiency (b) plots for the $\tilde{\chi}_1^+\tilde{\chi}_1^-$ production with $W$-boson-mediated decay model, in the SR-SF BDT-signal$\in(0.80,1]$ inclusive region. Acceptance is calculated by applying the signal region requirements to particle-level objects, which do not suffer from identification inefficiencies or mismeasurements. The efficiency is calculated with fully reconstructed objects with the acceptance divided out. Large acceptance and efficiency differences in neighbouring points are due to statistical fluctuations.

The figure shows the signal acceptance (a) and efficiency (b) plots for the $\tilde{\chi}_1^+\tilde{\chi}_1^-$ production with $W$-boson-mediated decay model, in the SR-SF BDT-signal$\in(0.77,0.775]$ inclusive region. Acceptance is calculated by applying the signal region requirements to particle-level objects, which do not suffer from identification inefficiencies or mismeasurements. The efficiency is calculated with fully reconstructed objects with the acceptance divided out. Large acceptance and efficiency differences in neighbouring points are due to statistical fluctuations.

The figure shows the signal acceptance (a) and efficiency (b) plots for the $\tilde{\chi}_1^+\tilde{\chi}_1^-$ production with $W$-boson-mediated decay model, in the SR-SF BDT-signal$\in(0.77,0.775]$ inclusive region. Acceptance is calculated by applying the signal region requirements to particle-level objects, which do not suffer from identification inefficiencies or mismeasurements. The efficiency is calculated with fully reconstructed objects with the acceptance divided out. Large acceptance and efficiency differences in neighbouring points are due to statistical fluctuations.

The figure shows the signal acceptance (a) and efficiency (b) plots for the $\tilde{\chi}_1^+\tilde{\chi}_1^-$ production with $W$-boson-mediated decay model, in the SR-SF BDT-signal$\in(0.775,0.78]$ inclusive region. Acceptance is calculated by applying the signal region requirements to particle-level objects, which do not suffer from identification inefficiencies or mismeasurements. The efficiency is calculated with fully reconstructed objects with the acceptance divided out. Large acceptance and efficiency differences in neighbouring points are due to statistical fluctuations.

The figure shows the signal acceptance (a) and efficiency (b) plots for the $\tilde{\chi}_1^+\tilde{\chi}_1^-$ production with $W$-boson-mediated decay model, in the SR-SF BDT-signal$\in(0.775,0.78]$ inclusive region. Acceptance is calculated by applying the signal region requirements to particle-level objects, which do not suffer from identification inefficiencies or mismeasurements. The efficiency is calculated with fully reconstructed objects with the acceptance divided out. Large acceptance and efficiency differences in neighbouring points are due to statistical fluctuations.

The figure shows the signal acceptance (a) and efficiency (b) plots for the $\tilde{\chi}_1^+\tilde{\chi}_1^-$ production with $W$-boson-mediated decay model, in the SR-SF BDT-signal$\in(0.78,0.785]$ inclusive region. Acceptance is calculated by applying the signal region requirements to particle-level objects, which do not suffer from identification inefficiencies or mismeasurements. The efficiency is calculated with fully reconstructed objects with the acceptance divided out. Large acceptance and efficiency differences in neighbouring points are due to statistical fluctuations.

The figure shows the signal acceptance (a) and efficiency (b) plots for the $\tilde{\chi}_1^+\tilde{\chi}_1^-$ production with $W$-boson-mediated decay model, in the SR-SF BDT-signal$\in(0.78,0.785]$ inclusive region. Acceptance is calculated by applying the signal region requirements to particle-level objects, which do not suffer from identification inefficiencies or mismeasurements. The efficiency is calculated with fully reconstructed objects with the acceptance divided out. Large acceptance and efficiency differences in neighbouring points are due to statistical fluctuations.

The figure shows the signal acceptance (a) and efficiency (b) plots for the $\tilde{\chi}_1^+\tilde{\chi}_1^-$ production with $W$-boson-mediated decay model, in the SR-SF BDT-signal$\in(0.785,0.79]$ inclusive region. Acceptance is calculated by applying the signal region requirements to particle-level objects, which do not suffer from identification inefficiencies or mismeasurements. The efficiency is calculated with fully reconstructed objects with the acceptance divided out. Large acceptance and efficiency differences in neighbouring points are due to statistical fluctuations.

The figure shows the signal acceptance (a) and efficiency (b) plots for the $\tilde{\chi}_1^+\tilde{\chi}_1^-$ production with $W$-boson-mediated decay model, in the SR-SF BDT-signal$\in(0.785,0.79]$ inclusive region. Acceptance is calculated by applying the signal region requirements to particle-level objects, which do not suffer from identification inefficiencies or mismeasurements. The efficiency is calculated with fully reconstructed objects with the acceptance divided out. Large acceptance and efficiency differences in neighbouring points are due to statistical fluctuations.

The figure shows the signal acceptance (a) and efficiency (b) plots for the $\tilde{\chi}_1^+\tilde{\chi}_1^-$ production with $W$-boson-mediated decay model, in the SR-SF BDT-signal$\in(0.79,0.795]$ inclusive region. Acceptance is calculated by applying the signal region requirements to particle-level objects, which do not suffer from identification inefficiencies or mismeasurements. The efficiency is calculated with fully reconstructed objects with the acceptance divided out. Large acceptance and efficiency differences in neighbouring points are due to statistical fluctuations.

The figure shows the signal acceptance (a) and efficiency (b) plots for the $\tilde{\chi}_1^+\tilde{\chi}_1^-$ production with $W$-boson-mediated decay model, in the SR-SF BDT-signal$\in(0.79,0.795]$ inclusive region. Acceptance is calculated by applying the signal region requirements to particle-level objects, which do not suffer from identification inefficiencies or mismeasurements. The efficiency is calculated with fully reconstructed objects with the acceptance divided out. Large acceptance and efficiency differences in neighbouring points are due to statistical fluctuations.

The figure shows the signal acceptance (a) and efficiency (b) plots for the $\tilde{\chi}_1^+\tilde{\chi}_1^-$ production with $W$-boson-mediated decay model, in the SR-SF BDT-signal$\in(0.795,0.80]$ inclusive region. Acceptance is calculated by applying the signal region requirements to particle-level objects, which do not suffer from identification inefficiencies or mismeasurements. The efficiency is calculated with fully reconstructed objects with the acceptance divided out. Large acceptance and efficiency differences in neighbouring points are due to statistical fluctuations.

The figure shows the signal acceptance (a) and efficiency (b) plots for the $\tilde{\chi}_1^+\tilde{\chi}_1^-$ production with $W$-boson-mediated decay model, in the SR-SF BDT-signal$\in(0.795,0.80]$ inclusive region. Acceptance is calculated by applying the signal region requirements to particle-level objects, which do not suffer from identification inefficiencies or mismeasurements. The efficiency is calculated with fully reconstructed objects with the acceptance divided out. Large acceptance and efficiency differences in neighbouring points are due to statistical fluctuations.

The figure shows the signal acceptance (a) and efficiency (b) plots for the $\tilde{\chi}_1^+\tilde{\chi}_1^-$ production with $W$-boson-mediated decay model, in the SR-SF BDT-signal$\in(0.80,0.81]$ inclusive region. Acceptance is calculated by applying the signal region requirements to particle-level objects, which do not suffer from identification inefficiencies or mismeasurements. The efficiency is calculated with fully reconstructed objects with the acceptance divided out. Large acceptance and efficiency differences in neighbouring points are due to statistical fluctuations.

The figure shows the signal acceptance (a) and efficiency (b) plots for the $\tilde{\chi}_1^+\tilde{\chi}_1^-$ production with $W$-boson-mediated decay model, in the SR-SF BDT-signal$\in(0.80,0.81]$ inclusive region. Acceptance is calculated by applying the signal region requirements to particle-level objects, which do not suffer from identification inefficiencies or mismeasurements. The efficiency is calculated with fully reconstructed objects with the acceptance divided out. Large acceptance and efficiency differences in neighbouring points are due to statistical fluctuations.

The figure shows the signal acceptance (a) and efficiency (b) plots for the $\tilde{\chi}_1^+\tilde{\chi}_1^-$ production with $W$-boson-mediated decay model, in the SR-SF BDT-signal$\in(0.81,1]$ inclusive region. Acceptance is calculated by applying the signal region requirements to particle-level objects, which do not suffer from identification inefficiencies or mismeasurements. The efficiency is calculated with fully reconstructed objects with the acceptance divided out. Large acceptance and efficiency differences in neighbouring points are due to statistical fluctuations.

The figure shows the signal acceptance (a) and efficiency (b) plots for the $\tilde{\chi}_1^+\tilde{\chi}_1^-$ production with $W$-boson-mediated decay model, in the SR-SF BDT-signal$\in(0.81,1]$ inclusive region. Acceptance is calculated by applying the signal region requirements to particle-level objects, which do not suffer from identification inefficiencies or mismeasurements. The efficiency is calculated with fully reconstructed objects with the acceptance divided out. Large acceptance and efficiency differences in neighbouring points are due to statistical fluctuations.

Cutflow table for the chargino signal sample with $m\tilde{\chi}_1^{\pm},\tilde{\chi}_1^0=(125,25)$ GeV, in the SR-SF BDT-signal$\in (0.77,1]$ and SR-DF BDT-signal$\in (0.81,1]$ regions. The yields include the process cross-section and are weighted to the 139 fb$^{-1}$ luminosity. 170000 events were generated for the sample.

Observed and expected exclusion limits on SUSY simplified models, with observed upper limits on signal cross-section (fb) overlaid, for chargino-pair production with $W$-boson-mediated decays in the $m(\tilde{\chi}_1^{\pm})-m(\tilde{\chi}_1^0)$ plane. The observed (solid thick line) and expected (thin dashed line) exclusion contours are indicated. The shaded band around the dashed line corresponds to the $\pm 1 \sigma$ variations in the expected limit, including all uncertainties except theoretical uncertainties in the signal cross-section. The dotted lines around the observed limit illustrate the change in the observed limit as the nominal signal cross-section is scaled up and down by the theoretical uncertainty. All limits are computed at 95% CL. The observed limits obtained at LEP and by the ATLAS experiment in previous searches are also shown. In case of the search performed on ATLAS Run 1 data at $\sqrt{s} = 8$ TeV no sensitivity was expected for the exclusion in the mass plane.

Observed and expected exclusion limits on SUSY simplified models for chargino-pair production with $W$-boson-mediated decays in the (a) $m(\tilde{\chi}_1^{\pm})-m(\tilde{\chi}_1^0)$ and (b) $m(\tilde{\chi}_1^{\pm})-\Delta m(\tilde{\chi}_1^{\pm},\tilde{\chi}_1^0)$ planes. The observed (solid thick line) and expected (thin dashed line) exclusion contours are indicated. The shaded band around the dashed line corresponds to the $\pm 1 \sigma$ variations in the expected limit, including all uncertainties except theoretical uncertainties in the signal cross-section. The dotted lines around the observed limit illustrate the change in the observed limit as the nominal signal cross-section is scaled up and down by the theoretical uncertainty. All limits are computed at 95% CL. The observed limits obtained at LEP and by the ATLAS experiment in previous searches are also shown. In case of the search performed on ATLAS Run 1 data at $\sqrt{s} = 8$ TeV no sensitivity was expected for the exclusion in the mass plane.

Observed and expected exclusion limits on SUSY simplified models for chargino-pair production with $W$-boson-mediated decays in the (a) $m(\tilde{\chi}_1^{\pm})-m(\tilde{\chi}_1^0)$ and (b) $m(\tilde{\chi}_1^{\pm})-\Delta m(\tilde{\chi}_1^{\pm},\tilde{\chi}_1^0)$ planes. The observed (solid thick line) and expected (thin dashed line) exclusion contours are indicated. The shaded band around the dashed line corresponds to the $\pm 1 \sigma$ variations in the expected limit, including all uncertainties except theoretical uncertainties in the signal cross-section. The dotted lines around the observed limit illustrate the change in the observed limit as the nominal signal cross-section is scaled up and down by the theoretical uncertainty. All limits are computed at 95% CL. The observed limits obtained at LEP and by the ATLAS experiment in previous searches are also shown. In case of the search performed on ATLAS Run 1 data at $\sqrt{s} = 8$ TeV no sensitivity was expected for the exclusion in the mass plane.

Observed and expected exclusion limits on SUSY simplified models for chargino-pair production with $W$-boson-mediated decays in the (a) $m(\tilde{\chi}_1^{\pm})-m(\tilde{\chi}_1^0)$ and (b) $m(\tilde{\chi}_1^{\pm})-\Delta m(\tilde{\chi}_1^{\pm},\tilde{\chi}_1^0)$ planes. The observed (solid thick line) and expected (thin dashed line) exclusion contours are indicated. The shaded band around the dashed line corresponds to the $\pm 1 \sigma$ variations in the expected limit, including all uncertainties except theoretical uncertainties in the signal cross-section. The dotted lines around the observed limit illustrate the change in the observed limit as the nominal signal cross-section is scaled up and down by the theoretical uncertainty. All limits are computed at 95% CL. The observed limits obtained at LEP and by the ATLAS experiment in previous searches are also shown. In case of the search performed on ATLAS Run 1 data at $\sqrt{s} = 8$ TeV no sensitivity was expected for the exclusion in the mass plane.

Observed and expected exclusion limits on SUSY simplified models for chargino-pair production with $W$-boson-mediated decays in the (a) $m(\tilde{\chi}_1^{\pm})-m(\tilde{\chi}_1^0)$ and (b) $m(\tilde{\chi}_1^{\pm})-\Delta m(\tilde{\chi}_1^{\pm},\tilde{\chi}_1^0)$ planes. The observed (solid thick line) and expected (thin dashed line) exclusion contours are indicated. The shaded band around the dashed line corresponds to the $\pm 1 \sigma$ variations in the expected limit, including all uncertainties except theoretical uncertainties in the signal cross-section. The dotted lines around the observed limit illustrate the change in the observed limit as the nominal signal cross-section is scaled up and down by the theoretical uncertainty. All limits are computed at 95% CL. The observed limits obtained at LEP and by the ATLAS experiment in previous searches are also shown. In case of the search performed on ATLAS Run 1 data at $\sqrt{s} = 8$ TeV no sensitivity was expected for the exclusion in the mass plane.

Observed and expected exclusion limits on SUSY simplified models for chargino-pair production with $W$-boson-mediated decays in the (a) $m(\tilde{\chi}_1^{\pm})-m(\tilde{\chi}_1^0)$ and (b) $m(\tilde{\chi}_1^{\pm})-\Delta m(\tilde{\chi}_1^{\pm},\tilde{\chi}_1^0)$ planes. The observed (solid thick line) and expected (thin dashed line) exclusion contours are indicated. The shaded band around the dashed line corresponds to the $\pm 1 \sigma$ variations in the expected limit, including all uncertainties except theoretical uncertainties in the signal cross-section. The dotted lines around the observed limit illustrate the change in the observed limit as the nominal signal cross-section is scaled up and down by the theoretical uncertainty. All limits are computed at 95% CL. The observed limits obtained at LEP and by the ATLAS experiment in previous searches are also shown. In case of the search performed on ATLAS Run 1 data at $\sqrt{s} = 8$ TeV no sensitivity was expected for the exclusion in the mass plane.

Observed and expected exclusion limits on SUSY simplified models for chargino-pair production with $W$-boson-mediated decays in the (a) $m(\tilde{\chi}_1^{\pm})-m(\tilde{\chi}_1^0)$ and (b) $m(\tilde{\chi}_1^{\pm})-\Delta m(\tilde{\chi}_1^{\pm},\tilde{\chi}_1^0)$ planes. The observed (solid thick line) and expected (thin dashed line) exclusion contours are indicated. The shaded band around the dashed line corresponds to the $\pm 1 \sigma$ variations in the expected limit, including all uncertainties except theoretical uncertainties in the signal cross-section. The dotted lines around the observed limit illustrate the change in the observed limit as the nominal signal cross-section is scaled up and down by the theoretical uncertainty. All limits are computed at 95% CL. The observed limits obtained at LEP and by the ATLAS experiment in previous searches are also shown. In case of the search performed on ATLAS Run 1 data at $\sqrt{s} = 8$ TeV no sensitivity was expected for the exclusion in the mass plane.

Observed and expected exclusion limits on SUSY simplified models for chargino-pair production with $W$-boson-mediated decays in the (a) $m(\tilde{\chi}_1^{\pm})-m(\tilde{\chi}_1^0)$ and (b) $m(\tilde{\chi}_1^{\pm})-\Delta m(\tilde{\chi}_1^{\pm},\tilde{\chi}_1^0)$ planes. The observed (solid thick line) and expected (thin dashed line) exclusion contours are indicated. The shaded band around the dashed line corresponds to the $\pm 1 \sigma$ variations in the expected limit, including all uncertainties except theoretical uncertainties in the signal cross-section. The dotted lines around the observed limit illustrate the change in the observed limit as the nominal signal cross-section is scaled up and down by the theoretical uncertainty. All limits are computed at 95% CL. The observed limits obtained at LEP and by the ATLAS experiment in previous searches are also shown. In case of the search performed on ATLAS Run 1 data at $\sqrt{s} = 8$ TeV no sensitivity was expected for the exclusion in the mass plane.

Observed and expected exclusion limits on SUSY simplified models for chargino-pair production with $W$-boson-mediated decays in the (a) $m(\tilde{\chi}_1^{\pm})-m(\tilde{\chi}_1^0)$ and (b) $m(\tilde{\chi}_1^{\pm})-\Delta m(\tilde{\chi}_1^{\pm},\tilde{\chi}_1^0)$ planes. The observed (solid thick line) and expected (thin dashed line) exclusion contours are indicated. The shaded band around the dashed line corresponds to the $\pm 1 \sigma$ variations in the expected limit, including all uncertainties except theoretical uncertainties in the signal cross-section. The dotted lines around the observed limit illustrate the change in the observed limit as the nominal signal cross-section is scaled up and down by the theoretical uncertainty. All limits are computed at 95% CL. The observed limits obtained at LEP and by the ATLAS experiment in previous searches are also shown. In case of the search performed on ATLAS Run 1 data at $\sqrt{s} = 8$ TeV no sensitivity was expected for the exclusion in the mass plane.

Observed and expected exclusion limits on SUSY simplified models for chargino-pair production with $W$-boson-mediated decays in the (a) $m(\tilde{\chi}_1^{\pm})-m(\tilde{\chi}_1^0)$ and (b) $m(\tilde{\chi}_1^{\pm})-\Delta m(\tilde{\chi}_1^{\pm},\tilde{\chi}_1^0)$ planes. The observed (solid thick line) and expected (thin dashed line) exclusion contours are indicated. The shaded band around the dashed line corresponds to the $\pm 1 \sigma$ variations in the expected limit, including all uncertainties except theoretical uncertainties in the signal cross-section. The dotted lines around the observed limit illustrate the change in the observed limit as the nominal signal cross-section is scaled up and down by the theoretical uncertainty. All limits are computed at 95% CL. The observed limits obtained at LEP and by the ATLAS experiment in previous searches are also shown. In case of the search performed on ATLAS Run 1 data at $\sqrt{s} = 8$ TeV no sensitivity was expected for the exclusion in the mass plane.

Observed and expected exclusion limits on SUSY simplified models for chargino-pair production with $W$-boson-mediated decays in the (a) $m(\tilde{\chi}_1^{\pm})-m(\tilde{\chi}_1^0)$ and (b) $m(\tilde{\chi}_1^{\pm})-\Delta m(\tilde{\chi}_1^{\pm},\tilde{\chi}_1^0)$ planes. The observed (solid thick line) and expected (thin dashed line) exclusion contours are indicated. The shaded band around the dashed line corresponds to the $\pm 1 \sigma$ variations in the expected limit, including all uncertainties except theoretical uncertainties in the signal cross-section. The dotted lines around the observed limit illustrate the change in the observed limit as the nominal signal cross-section is scaled up and down by the theoretical uncertainty. All limits are computed at 95% CL. The observed limits obtained at LEP and by the ATLAS experiment in previous searches are also shown. In case of the search performed on ATLAS Run 1 data at $\sqrt{s} = 8$ TeV no sensitivity was expected for the exclusion in the mass plane.

Observed and expected exclusion limits on SUSY simplified models for chargino-pair production with $W$-boson-mediated decays in the (a) $m(\tilde{\chi}_1^{\pm})-m(\tilde{\chi}_1^0)$ and (b) $m(\tilde{\chi}_1^{\pm})-\Delta m(\tilde{\chi}_1^{\pm},\tilde{\chi}_1^0)$ planes. The observed (solid thick line) and expected (thin dashed line) exclusion contours are indicated. The shaded band around the dashed line corresponds to the $\pm 1 \sigma$ variations in the expected limit, including all uncertainties except theoretical uncertainties in the signal cross-section. The dotted lines around the observed limit illustrate the change in the observed limit as the nominal signal cross-section is scaled up and down by the theoretical uncertainty. All limits are computed at 95% CL. The observed limits obtained at LEP and by the ATLAS experiment in previous searches are also shown. In case of the search performed on ATLAS Run 1 data at $\sqrt{s} = 8$ TeV no sensitivity was expected for the exclusion in the mass plane.

Observed and expected exclusion limits on SUSY simplified models for chargino-pair production with $W$-boson-mediated decays in the (a) $m(\tilde{\chi}_1^{\pm})-m(\tilde{\chi}_1^0)$ and (b) $m(\tilde{\chi}_1^{\pm})-\Delta m(\tilde{\chi}_1^{\pm},\tilde{\chi}_1^0)$ planes. The observed (solid thick line) and expected (thin dashed line) exclusion contours are indicated. The shaded band around the dashed line corresponds to the $\pm 1 \sigma$ variations in the expected limit, including all uncertainties except theoretical uncertainties in the signal cross-section. The dotted lines around the observed limit illustrate the change in the observed limit as the nominal signal cross-section is scaled up and down by the theoretical uncertainty. All limits are computed at 95% CL. The observed limits obtained at LEP and by the ATLAS experiment in previous searches are also shown. In case of the search performed on ATLAS Run 1 data at $\sqrt{s} = 8$ TeV no sensitivity was expected for the exclusion in the mass plane.

The upper panel shows the observed number of events in the SRs defined in Table 3, together with the expected SM backgrounds obtained after the background fit in the CRs. `Others' include the non-dominant background sources, e.g.$t \bar{t}$+$V$, Higgs boson and Drell--Yan events. The uncertainty band includes systematic and statistical errors from all sources. Distributions for three benchmark signal points are overlaid for comparison. The lower panel shows the significance as defined in Ref. [115].

Observed (points with error bars) and expected background (solid histograms) distributions for $E_{T}^{miss}$ in the signal region (a) SRL, (b) SRM and (c) SRH after the background-only fit applied to the CRs. The predicted signal distributions for the two models with a gluino mass of 2000 GeV and neutralino mass of 250 GeV (SRL), 1050 GeV (SRM) or 1950 GeV (SRH) are also shown for comparison. The uncertainties in the SM background are only statistical.

Observed and expected exclusion limit in the gluino-neutralino mass plane at 95% CL combined using the signal region with the best expected sensitivity at each point, for the full Run-2 dataset corresponding to an integrated luminosity of $139~\mathrm{fb}^{-1}$, for $\gamma/Z$ (a) and $\gamma/h$ (b) signal models. The black solid line corresponds to the expected limits at 95% CL, with the light (yellow) bands indicating the 1$\sigma$ exclusions due to experimental and background-theory uncertainties. The observed limits are indicated by medium (red) curves, the solid contour represents the nominal limit, and the dotted lines are obtained by varying the signal cross section by the theoretical scale and PDF uncertainties. For each point in the higgsino-bino parameter space, the labels indicate the best-expected signal region, where L, M and H mean SRL, SRM and SRH, respectively.

Observed and expected exclusion limit in the gluino-neutralino mass plane at 95% CL combined using the signal region with the best expected sensitivity at each point, for the full Run-2 dataset corresponding to an integrated luminosity of $139~\mathrm{fb}^{-1}$, for $\gamma/Z$ (a) and $\gamma/h$ (b) signal models. The black solid line corresponds to the expected limits at 95% CL, with the light (yellow) bands indicating the 1$\sigma$ exclusions due to experimental and background-theory uncertainties. The observed limits are indicated by medium (red) curves, the solid contour represents the nominal limit, and the dotted lines are obtained by varying the signal cross section by the theoretical scale and PDF uncertainties. For each point in the higgsino-bino parameter space, the labels indicate the best-expected signal region, where L, M and H mean SRL, SRM and SRH, respectively.

Acceptance (left) and efficiency (right) for the $\gamma/Z$ model signal grid for SRL (top), SRM (middle) and SRH (bottom).

Acceptance (left) and efficiency (right) for the $\gamma/Z$ model signal grid for SRL (top), SRM (middle) and SRH (bottom).

Acceptance (left) and efficiency (right) for the $\gamma/Z$ model signal grid for SRL (top), SRM (middle) and SRH (bottom).

Acceptance (left) and efficiency (right) for the $\gamma/Z$ model signal grid for SRL (top), SRM (middle) and SRH (bottom).

Acceptance (left) and efficiency (right) for the $\gamma/Z$ model signal grid for SRL (top), SRM (middle) and SRH (bottom).

Acceptance (left) and efficiency (right) for the $\gamma/Z$ model signal grid for SRL (top), SRM (middle) and SRH (bottom).

Acceptance (left) and efficiency (right) for the $\gamma/h$ model signal grid for SRL (top), SRM (middle) and SRH (bottom).

Acceptance (left) and efficiency (right) for the $\gamma/h$ model signal grid for SRL (top), SRM (middle) and SRH (bottom).

Acceptance (left) and efficiency (right) for the $\gamma/h$ model signal grid for SRL (top), SRM (middle) and SRH (bottom).

Acceptance (left) and efficiency (right) for the $\gamma/h$ model signal grid for SRL (top), SRM (middle) and SRH (bottom).

Acceptance (left) and efficiency (right) for the $\gamma/h$ model signal grid for SRL (top), SRM (middle) and SRH (bottom).

Acceptance (left) and efficiency (right) for the $\gamma/h$ model signal grid for SRL (top), SRM (middle) and SRH (bottom).

Cutflow for the SRL selection, for two relevant signal points for both $\gamma/Z$ and $\gamma/h$ models, where the gluinos have mass of 2000 GeV and the neutralinos have a mass of 250 GeV (10000 generated events). The numbers are normalized to a luminosity of 139 $fb^{-1}$.

Cutflow for the SRM selection, for two relevant signal points for both $\gamma/Z$ and $\gamma/h$ models, where the gluinos have mass of 2000 GeV and the neutralinos have a mass of 1050 GeV (10000 generated events). The numbers are normalized to a luminosity of 139 $fb^{-1}$.

Cutflow for the SRH selection, for two relevant signal points for both $\gamma/Z$ and $\gamma/h$ models, where the gluinos have mass of 2000 GeV and the neutralinos have a mass of 1950 GeV (10000 generated events). The numbers are normalized to a luminosity of 139 $fb^{-1}$.

Observed and expected exclusion limits in the gluino–neutralino mass plane at 95% CL for the full Run-2 dataset corresponding to an integrated luminosity of 139 fb−1 , for the (a) $\gamma/Z$ and (b) $\gamma/h$ signal models. They are obtained by combining limits from the signal region with the best expected sensitivity at each point. The dashed (black) line corresponds to the expected limits at 95% CL, with the light (yellow) band indicating the $\pm 1\sigma$ excursions due to experimental and background-theory uncertainties. The observed limits are indicated by medium (red) curves: the solid contour represents the nominal limit, and the dotted lines are obtained by varying the signal cross section by the theoretical scale and PDF uncertainties.

Observed and expected exclusion limits in the gluino–neutralino mass plane at 95% CL for the full Run-2 dataset corresponding to an integrated luminosity of 139 fb−1 , for the (a) $\gamma/Z$ and (b) $\gamma/h$ signal models. They are obtained by combining limits from the signal region with the best expected sensitivity at each point. The dashed (black) line corresponds to the expected limits at 95% CL, with the light (yellow) band indicating the $\pm 1\sigma$ excursions due to experimental and background-theory uncertainties. The observed limits are indicated by medium (red) curves: the solid contour represents the nominal limit, and the dotted lines are obtained by varying the signal cross section by the theoretical scale and PDF uncertainties.

Observed and expected exclusion limits in the gluino–neutralino mass plane at 95% CL for the full Run-2 dataset corresponding to an integrated luminosity of 139 fb−1 , for the (a) $\gamma/Z$ and (b) $\gamma/h$ signal models. They are obtained by combining limits from the signal region with the best expected sensitivity at each point. The dashed (black) line corresponds to the expected limits at 95% CL, with the light (yellow) band indicating the $\pm 1\sigma$ excursions due to experimental and background-theory uncertainties. The observed limits are indicated by medium (red) curves: the solid contour represents the nominal limit, and the dotted lines are obtained by varying the signal cross section by the theoretical scale and PDF uncertainties.

Observed and expected exclusion limits in the gluino–neutralino mass plane at 95% CL for the full Run-2 dataset corresponding to an integrated luminosity of 139 fb−1 , for the (a) $\gamma/Z$ and (b) $\gamma/h$ signal models. They are obtained by combining limits from the signal region with the best expected sensitivity at each point. The dashed (black) line corresponds to the expected limits at 95% CL, with the light (yellow) band indicating the $\pm 1\sigma$ excursions due to experimental and background-theory uncertainties. The observed limits are indicated by medium (red) curves: the solid contour represents the nominal limit, and the dotted lines are obtained by varying the signal cross section by the theoretical scale and PDF uncertainties.

Observed and expected exclusion limits in the gluino–neutralino mass plane at 95% CL for the full Run-2 dataset corresponding to an integrated luminosity of 139 fb−1 , for the (a) $\gamma/Z$ and (b) $\gamma/h$ signal models. They are obtained by combining limits from the signal region with the best expected sensitivity at each point. The dashed (black) line corresponds to the expected limits at 95% CL, with the light (yellow) band indicating the $\pm 1\sigma$ excursions due to experimental and background-theory uncertainties. The observed limits are indicated by medium (red) curves: the solid contour represents the nominal limit, and the dotted lines are obtained by varying the signal cross section by the theoretical scale and PDF uncertainties.

Observed and expected exclusion limits in the gluino–neutralino mass plane at 95% CL for the full Run-2 dataset corresponding to an integrated luminosity of 139 fb−1 , for the (a) $\gamma/Z$ and (b) $\gamma/h$ signal models. They are obtained by combining limits from the signal region with the best expected sensitivity at each point. The dashed (black) line corresponds to the expected limits at 95% CL, with the light (yellow) band indicating the $\pm 1\sigma$ excursions due to experimental and background-theory uncertainties. The observed limits are indicated by medium (red) curves: the solid contour represents the nominal limit, and the dotted lines are obtained by varying the signal cross section by the theoretical scale and PDF uncertainties.

Observed and expected exclusion limits in the gluino–neutralino mass plane at 95% CL for the full Run-2 dataset corresponding to an integrated luminosity of 139 fb−1 , for the (a) $\gamma/Z$ and (b) $\gamma/h$ signal models. They are obtained by combining limits from the signal region with the best expected sensitivity at each point. The dashed (black) line corresponds to the expected limits at 95% CL, with the light (yellow) band indicating the $\pm 1\sigma$ excursions due to experimental and background-theory uncertainties. The observed limits are indicated by medium (red) curves: the solid contour represents the nominal limit, and the dotted lines are obtained by varying the signal cross section by the theoretical scale and PDF uncertainties.

Observed and expected exclusion limits in the gluino–neutralino mass plane at 95% CL for the full Run-2 dataset corresponding to an integrated luminosity of 139 fb−1 , for the (a) $\gamma/Z$ and (b) $\gamma/h$ signal models. They are obtained by combining limits from the signal region with the best expected sensitivity at each point. The dashed (black) line corresponds to the expected limits at 95% CL, with the light (yellow) band indicating the $\pm 1\sigma$ excursions due to experimental and background-theory uncertainties. The observed limits are indicated by medium (red) curves: the solid contour represents the nominal limit, and the dotted lines are obtained by varying the signal cross section by the theoretical scale and PDF uncertainties.

Observed and expected exclusion limits in the gluino–neutralino mass plane at 95% CL for the full Run-2 dataset corresponding to an integrated luminosity of 139 fb−1 , for the (a) $\gamma/Z$ and (b) $\gamma/h$ signal models. They are obtained by combining limits from the signal region with the best expected sensitivity at each point. The dashed (black) line corresponds to the expected limits at 95% CL, with the light (yellow) band indicating the $\pm 1\sigma$ excursions due to experimental and background-theory uncertainties. The observed limits are indicated by medium (red) curves: the solid contour represents the nominal limit, and the dotted lines are obtained by varying the signal cross section by the theoretical scale and PDF uncertainties.

Observed and expected exclusion limits in the gluino–neutralino mass plane at 95% CL for the full Run-2 dataset corresponding to an integrated luminosity of 139 fb−1 , for the (a) $\gamma/Z$ and (b) $\gamma/h$ signal models. They are obtained by combining limits from the signal region with the best expected sensitivity at each point. The dashed (black) line corresponds to the expected limits at 95% CL, with the light (yellow) band indicating the $\pm 1\sigma$ excursions due to experimental and background-theory uncertainties. The observed limits are indicated by medium (red) curves: the solid contour represents the nominal limit, and the dotted lines are obtained by varying the signal cross section by the theoretical scale and PDF uncertainties.

Observed and expected exclusion limits in the gluino–neutralino mass plane at 95% CL for the full Run-2 dataset corresponding to an integrated luminosity of 139 fb−1 , for the (a) $\gamma/Z$ and (b) $\gamma/h$ signal models. They are obtained by combining limits from the signal region with the best expected sensitivity at each point. The dashed (black) line corresponds to the expected limits at 95% CL, with the light (yellow) band indicating the $\pm 1\sigma$ excursions due to experimental and background-theory uncertainties. The observed limits are indicated by medium (red) curves: the solid contour represents the nominal limit, and the dotted lines are obtained by varying the signal cross section by the theoretical scale and PDF uncertainties.

Observed and expected exclusion limits in the gluino–neutralino mass plane at 95% CL for the full Run-2 dataset corresponding to an integrated luminosity of 139 fb−1 , for the (a) $\gamma/Z$ and (b) $\gamma/h$ signal models. They are obtained by combining limits from the signal region with the best expected sensitivity at each point. The dashed (black) line corresponds to the expected limits at 95% CL, with the light (yellow) band indicating the $\pm 1\sigma$ excursions due to experimental and background-theory uncertainties. The observed limits are indicated by medium (red) curves: the solid contour represents the nominal limit, and the dotted lines are obtained by varying the signal cross section by the theoretical scale and PDF uncertainties.

Upper limits at 95% CL on $\sigma\times\mathcal{B}^2(\mathrm{H}_1\rightarrow b\bar{b})$ as a function of $m_{\mathrm{H_1}}$ for $m_{\mathrm{SUSY}}=1600\;\mathrm{GeV}$.

Upper limits at 95% CL on $\sigma\times\mathcal{B}^2(\mathrm{H}_1\rightarrow b\bar{b})$ as a function of $m_{\mathrm{H_1}}$ for $m_{\mathrm{SUSY}}=2000\;\mathrm{GeV}$.

Upper limits at 95% CL on $\sigma\times\mathcal{B}^2(\mathrm{H}_1\rightarrow b\bar{b})$ as a function of $m_{\mathrm{H_1}}$ for $m_{\mathrm{SUSY}}=2200\;\mathrm{GeV}$.

Upper limits at 95% CL on $\sigma\times\mathcal{B}^2(\mathrm{H}_1\rightarrow b\bar{b})$ as a function of $m_{\mathrm{H_1}}$ for $m_{\mathrm{SUSY}}=2400\;\mathrm{GeV}$.

Upper limits at 95% CL on $\sigma\times\mathcal{B}^2(\mathrm{H}_1\rightarrow b\bar{b})$ as a function of $m_{\mathrm{H_1}}$ for $m_{\mathrm{SUSY}}=2600\;\mathrm{GeV}$.

Upper limits at 95% CL on $\sigma\times\mathcal{B}^2(\mathrm{H}_1\rightarrow b\bar{b})$ as a function of $m_{\mathrm{H_1}}$ for $m_{\mathrm{SUSY}}=2800\;\mathrm{GeV}$.

Product of efficiency and acceptance of the event selection for T signal events as a function of the particle mass $m_\mathrm{T}$ and width $\Gamma$ for the different hypotheses considered.

Product of efficiency and acceptance of the event selection for T signal events as a function of the particle mass $m_\mathrm{T}$ and width $\Gamma$ for the different hypotheses considered.

Product of efficiency and acceptance of the event selection for T signal events as a function of the particle mass $m_\mathrm{T}$ and width $\Gamma$ for the different hypotheses considered.

Distribution of the transverse mass $M_T$ reconstructed from the top quark and the missing transverse energy in the merged validation region with 0 forward jets. Data are compared to backgrounds are after the data-driven extraction is performed.

Distribution of the transverse mass $M_T$ reconstructed from the top quark and the missing transverse energy in the merged validation region with at least 1 forward jet. Data are compared to backgrounds are after the data-driven extraction is performed.

Distribution of the transverse mass $M_T$ reconstructed from the top quark and the missing transverse energy in the merged validation region with 0 forward jets. Data are compared to backgrounds are after the data-driven extraction is performed.

Distribution of the transverse mass $M_T$ reconstructed from the top quark and the missing transverse energy in the merged validation region with at least 1 forward jet. Data are compared to backgrounds are after the data-driven extraction is performed.

Distribution of the transverse mass $M_T$ reconstructed from the top quark and the missing transverse energy in the merged validation region with 0 forward jets. Data are compared to backgrounds are after the data-driven extraction is performed.

Distribution of the transverse mass $M_T$ reconstructed from the top quark and the missing transverse energy in the merged validation region with at least 1 forward jet. Data are compared to backgrounds are after the data-driven extraction is performed.

Distribution of the transverse mass $M_T$ reconstructed from the top quark and the missing transverse energy in the resolved validation region with 0 forward jets. Data are compared to backgrounds are after the data-driven extraction is performed.

Distribution of the transverse mass $M_T$ reconstructed from the top quark and the missing transverse energy in the resolved validation region with at least 1 forward jet. Data are compared to backgrounds are after the data-driven extraction is performed.

Distribution of the transverse mass $M_T$ reconstructed from the top quark and the missing transverse energy in the resolved validation region with 0 forward jets. Data are compared to backgrounds are after the data-driven extraction is performed.

Distribution of the transverse mass $M_T$ reconstructed from the top quark and the missing transverse energy in the resolved validation region with at least 1 forward jet. Data are compared to backgrounds are after the data-driven extraction is performed.

Distribution of the transverse mass $M_T$ reconstructed from the top quark and the missing transverse energy in the resolved validation region with 0 forward jets. Data are compared to backgrounds are after the data-driven extraction is performed.

Distribution of the transverse mass $M_T$ reconstructed from the top quark and the missing transverse energy in the resolved validation region with at least 1 forward jet. Data are compared to backgrounds are after the data-driven extraction is performed.

Distribution of the transverse mass $M_T$ reconstructed from the top quark and the missing transverse energy in the merged region with 0 forward jets. Data are compared to backgrounds obtained post-fit.

Distribution of the transverse mass $M_T$ reconstructed from the top quark and the missing transverse energy in the merged region with at least 1 forward jet. Data are compared to backgrounds obtained post-fit.

Distribution of the transverse mass $M_T$ reconstructed from the top quark and the missing transverse energy in the merged region with 0 forward jets. Data are compared to backgrounds obtained post-fit.

Distribution of the transverse mass $M_T$ reconstructed from the top quark and the missing transverse energy in the merged region with at least 1 forward jet. Data are compared to backgrounds obtained post-fit.

Distribution of the transverse mass $M_T$ reconstructed from the top quark and the missing transverse energy in the merged region with 0 forward jets. Data are compared to backgrounds obtained post-fit.

Distribution of the transverse mass $M_T$ reconstructed from the top quark and the missing transverse energy in the merged region with at least 1 forward jet. Data are compared to backgrounds obtained post-fit.

Distribution of the transverse mass $M_T$ reconstructed from the top quark and the missing transverse energy in the partially merged region with 0 forward jets. Data are compared to backgrounds obtained post-fit.

Distribution of the transverse mass $M_T$ reconstructed from the top quark and the missing transverse energy in the partially merged region with at least 1 forward jet. Data are compared to backgrounds obtained post-fit.

Distribution of the transverse mass $M_T$ reconstructed from the top quark and the missing transverse energy in the partially merged region with 0 forward jets. Data are compared to backgrounds obtained post-fit.

Distribution of the transverse mass $M_T$ reconstructed from the top quark and the missing transverse energy in the partially merged region with at least 1 forward jet. Data are compared to backgrounds obtained post-fit.

Distribution of the transverse mass $M_T$ reconstructed from the top quark and the missing transverse energy in the partially merged region with 0 forward jets. Data are compared to backgrounds obtained post-fit.

Distribution of the transverse mass $M_T$ reconstructed from the top quark and the missing transverse energy in the partially merged region with at least 1 forward jet. Data are compared to backgrounds obtained post-fit.

Distribution of the transverse mass $M_T$ reconstructed from the top quark and the missing transverse energy in the resolved region with 0 forward jets. Data are compared to backgrounds obtained post-fit.

Distribution of the transverse mass $M_T$ reconstructed from the top quark and the missing transverse energy in the resolved region with at least 1 forward jet. Data are compared to backgrounds obtained post-fit.

Distribution of the transverse mass $M_T$ reconstructed from the top quark and the missing transverse energy in the resolved region with 0 forward jets. Data are compared to backgrounds obtained post-fit.

Distribution of the transverse mass $M_T$ reconstructed from the top quark and the missing transverse energy in the resolved region with at least 1 forward jet. Data are compared to backgrounds obtained post-fit.

Distribution of the transverse mass $M_T$ reconstructed from the top quark and the missing transverse energy in the resolved region with 0 forward jets. Data are compared to backgrounds obtained post-fit.

Distribution of the transverse mass $M_T$ reconstructed from the top quark and the missing transverse energy in the resolved region with at least 1 forward jet. Data are compared to backgrounds obtained post-fit.

Observed and expected 95% CL upper limits on the product of the cross section and branching fraction of a T decaying to a hadronic top quark and a Z boson to neutrinos as a function of the T mass, for a T with a width of 0.01 $ imes$ its mass.

Observed and expected 95% CL upper limits on the product of the cross section and branching fraction of a T decaying to a hadronic top quark and a Z boson to neutrinos as a function of the T mass, for a T with a width of 0.1 $ imes$ its mass.

Observed and expected 95% CL upper limits on the product of the cross section and branching fraction of a T decaying to a hadronic top quark and a Z boson to neutrinos as a function of the T mass, for a T with a width of 0.2 $ imes$ its mass.

Observed and expected 95% CL upper limits on the product of the cross section and branching fraction of a T decaying to a hadronic top quark and a Z boson to neutrinos as a function of the T mass, for a T with a width of 0.3 $ imes$ its mass.

Model independent observed 95% CL upper limits on the product of the cross section and branching fraction of a T decaying to a top quark and a Z boson as a function of the T mass and width.

Singlet model 95% CL excluded values of the product of the cross section and branching fraction for a T decaying to a top quark and a Z boson.

The normalized distributions of the BDT input variable $D_{p_{\text{T}}}$ for the two highest $p_{\text{T}}$ jets from the simulated SM backgrounds and several signal models. Each sample's jet $p_{\text{T}}$ distribution is weighted to match a reference distribution (see text). The last bin of each histogram includes the overflow events.

The normalized BDT discriminator distribution for the two highest $p_{\text{T}}$ jets from the simulated SM backgrounds and several signal models.

The BDT ROC curves for the two highest $p_{\text{T}}$ jets, comparing the simulated SM backgrounds with one signal model with $m_{\text{Z}^{\prime}} = 3.1\,\text{TeV}$, $m_{\text{dark}} = 20\,\text{GeV}$, $r_{\text{inv}} = 0.3$, $\alpha_{\text{dark}} = \alpha_{\text{dark}}^{\text{peak}}$.

The BDT ROC curves for the two highest $p_{\text{T}}$ jets, comparing the simulated SM backgrounds with one signal model with $m_{\text{Z}^{\prime}} = 3.1\,\text{TeV}$, $m_{\text{dark}} = 20\,\text{GeV}$, $r_{\text{inv}} = 0.3$, $\alpha_{\text{dark}} = \alpha_{\text{dark}}^{\text{peak}}$.

The BDT ROC curves for the two highest $p_{\text{T}}$ jets, comparing the simulated SM backgrounds with one signal model with $m_{\text{Z}^{\prime}} = 3.1\,\text{TeV}$, $m_{\text{dark}} = 20\,\text{GeV}$, $r_{\text{inv}} = 0.3$, $\alpha_{\text{dark}} = \alpha_{\text{dark}}^{\text{peak}}$.

The BDT ROC curves for the two highest $p_{\text{T}}$ jets, comparing the simulated SM backgrounds with one signal model with $m_{\text{Z}^{\prime}} = 3.1\,\text{TeV}$, $m_{\text{dark}} = 20\,\text{GeV}$, $r_{\text{inv}} = 0.3$, $\alpha_{\text{dark}} = \alpha_{\text{dark}}^{\text{peak}}$.

The $m_{\text{T}}$ distribution for the high-$R_{\text{T}}$ signal region, comparing the observed data to the background prediction from the analytic fit ($g_{3}(x) = \exp(p_{1}x)x^{p_{2}(1+p_{3}\ln(x))}$, $x = m_{\text{T}}/\sqrt{s}$). The distributions from several example signal models, with cross sections corresponding to the observed limits, are superimposed.

The $m_{\text{T}}$ distribution for the low-$R_{\text{T}}$ signal region, comparing the observed data to the background prediction from the analytic fit ($g_{2}(x) = \exp(p_{1}x)x^{p_{2}}$, $x = m_{\text{T}}/\sqrt{s}$). The distributions from several example signal models, with cross sections corresponding to the observed limits, are superimposed.

The $m_{\text{T}}$ distribution for the high-SVJ2 signal region, comparing the observed data to the background prediction from the analytic fit ($g_{2}(x) = \exp(p_{1}x)x^{p_{2}}$, $x = m_{\text{T}}/\sqrt{s}$). The distributions from several example signal models, with cross sections corresponding to the observed limits, are superimposed.

The $m_{\text{T}}$ distribution for the low-SVJ2 signal region, comparing the observed data to the background prediction from the analytic fit ($g_{2}(x) = \exp(p_{1}x)x^{p_{2}}$, $x = m_{\text{T}}/\sqrt{s}$). The distributions from several example signal models, with cross sections corresponding to the observed limits, are superimposed.

The 95% CL observed upper limits on the product of the cross section and branching fraction from the inclusive search for variations of the mediator mass and the dark hadron mass.

The expected exclusion for the nominal $\text{Z}^{\prime}$ cross section.

The expected exclusion for the nominal $\text{Z}^{\prime}$ cross section.

The observed exclusion for the nominal $\text{Z}^{\prime}$ cross section.

The observed exclusion for the nominal $\text{Z}^{\prime}$ cross section.

The lower 68% expected exclusion for the nominal $\text{Z}^{\prime}$ cross section.

The lower 68% expected exclusion for the nominal $\text{Z}^{\prime}$ cross section.

The upper 68% expected exclusion for the nominal $\text{Z}^{\prime}$ cross section.

The upper 68% expected exclusion for the nominal $\text{Z}^{\prime}$ cross section.

The lower 95% expected exclusion for the nominal $\text{Z}^{\prime}$ cross section.

The lower 95% expected exclusion for the nominal $\text{Z}^{\prime}$ cross section.

The upper 95% expected exclusion for the nominal $\text{Z}^{\prime}$ cross section.

The 95% CL observed upper limits on the product of the cross section and branching fraction from the inclusive search for variations of the mediator mass and the invisible fraction.

The expected exclusion for the nominal $\text{Z}^{\prime}$ cross section.

The observed exclusion for the nominal $\text{Z}^{\prime}$ cross section.

The lower 68% expected exclusion for the nominal $\text{Z}^{\prime}$ cross section.

The lower 68% expected exclusion for the nominal $\text{Z}^{\prime}$ cross section.

The upper 68% expected exclusion for the nominal $\text{Z}^{\prime}$ cross section.

The lower 95% expected exclusion for the nominal $\text{Z}^{\prime}$ cross section.

The upper 95% expected exclusion for the nominal $\text{Z}^{\prime}$ cross section.

The upper 95% expected exclusion for the nominal $\text{Z}^{\prime}$ cross section.

The 95% CL observed upper limits on the product of the cross section and branching fraction from the BDT-based search for variations of the mediator mass and the dark hadron mass.

The expected exclusion for the nominal $\text{Z}^{\prime}$ cross section.

The expected exclusion for the nominal $\text{Z}^{\prime}$ cross section.

The expected exclusion for the nominal $\text{Z}^{\prime}$ cross section.

The observed exclusion for the nominal $\text{Z}^{\prime}$ cross section.

The observed exclusion for the nominal $\text{Z}^{\prime}$ cross section.

The observed exclusion for the nominal $\text{Z}^{\prime}$ cross section.

The observed exclusion for the nominal $\text{Z}^{\prime}$ cross section.

The lower 68% expected exclusion for the nominal $\text{Z}^{\prime}$ cross section.

The lower 68% expected exclusion for the nominal $\text{Z}^{\prime}$ cross section.

The lower 68% expected exclusion for the nominal $\text{Z}^{\prime}$ cross section.

The upper 68% expected exclusion for the nominal $\text{Z}^{\prime}$ cross section.

The upper 68% expected exclusion for the nominal $\text{Z}^{\prime}$ cross section.

The upper 68% expected exclusion for the nominal $\text{Z}^{\prime}$ cross section.

The lower 95% expected exclusion for the nominal $\text{Z}^{\prime}$ cross section.

The lower 95% expected exclusion for the nominal $\text{Z}^{\prime}$ cross section.

The upper 95% expected exclusion for the nominal $\text{Z}^{\prime}$ cross section.

The upper 95% expected exclusion for the nominal $\text{Z}^{\prime}$ cross section.

The upper 95% expected exclusion for the nominal $\text{Z}^{\prime}$ cross section.

The 95% CL observed upper limits on the product of the cross section and branching fraction from the BDT-based search for variations of the mediator mass and the invisible fraction.

The expected exclusion for the nominal $\text{Z}^{\prime}$ cross section.

The expected exclusion for the nominal $\text{Z}^{\prime}$ cross section.

The expected exclusion for the nominal $\text{Z}^{\prime}$ cross section.

The observed exclusion for the nominal $\text{Z}^{\prime}$ cross section.

The observed exclusion for the nominal $\text{Z}^{\prime}$ cross section.

The lower 68% expected exclusion for the nominal $\text{Z}^{\prime}$ cross section.

The lower 68% expected exclusion for the nominal $\text{Z}^{\prime}$ cross section.

The lower 68% expected exclusion for the nominal $\text{Z}^{\prime}$ cross section.

The lower 68% expected exclusion for the nominal $\text{Z}^{\prime}$ cross section.

The upper 68% expected exclusion for the nominal $\text{Z}^{\prime}$ cross section.

The lower 95% expected exclusion for the nominal $\text{Z}^{\prime}$ cross section.

The lower 95% expected exclusion for the nominal $\text{Z}^{\prime}$ cross section.

The lower 95% expected exclusion for the nominal $\text{Z}^{\prime}$ cross section.

The lower 95% expected exclusion for the nominal $\text{Z}^{\prime}$ cross section.

The upper 95% expected exclusion for the nominal $\text{Z}^{\prime}$ cross section.

The 95% CL upper limits on the product of the cross section and branching fraction from the inclusive search for the $\alpha_{\text{dark}}$ variations.

The 95% CL upper limits on the product of the cross section and branching fraction from the BDT-based search for the $\alpha_{\text{dark}}$ variations.

The three two-dimensional signal model parameter scans.

Metrics representing the performance of the BDT for the benchmark signal model ($m_{\text{Z}^{\prime}} = 3.1\,\text{TeV}$, $m_{\text{dark}} = 20\,\text{GeV}$, $r_{\text{inv}} = 0.3$, $\alpha_{\text{dark}} = \alpha_{\text{dark}}^{\text{peak}}$), compared to each of the major SM background processes.

The range of effects on the signal yield for each systematic uncertainty and the total. Values less than 0.01% are rounded to 0.0%.

The normalized distribution of the variable $m_{\text{T}}$ for the simulated SM backgrounds and several signal models. The requirement on this variable is omitted, but all other preselection requirements are applied. The last bin of each histogram includes the overflow events.

The normalized distribution of the variable $\Delta\eta(\text{J}_{1},\text{J}_{2})$ for the simulated SM backgrounds and several signal models. The requirement on this variable is omitted, but all other preselection requirements are applied. The last bin of each histogram includes the overflow events.

The normalized distribution of the variable $p_{\text{T}}^{\text{miss}}$ for the simulated SM backgrounds and several signal models. The $R_{\text{T}}$ requirement is omitted, but all other preselection requirements are applied. The last bin of each histogram includes the overflow events.

The normalized distribution of the variable $N_{\text{e}}$ for the simulated SM backgrounds and several signal models. The requirement on this variable is omitted, but all other preselection requirements are applied. The last bin of each histogram includes the overflow events.

The normalized distribution of the variable $N_{\mu}$ for the simulated SM backgrounds and several signal models. The requirement on this variable is omitted, but all other preselection requirements are applied. The last bin of each histogram includes the overflow events.

The normalized distribution of $\Delta\eta(\text{J}_{1},\text{J}_{2})$ vs. $R_{\text{T}}$ for the simulated QCD background. The preselection requirements on both variables are omitted, but all other preselection requirements are applied.

The normalized distribution of $p_{\text{T}}^{\text{miss}}$ vs. $m_{\text{T}}$ for the simulated QCD background. All selection requirements are omitted, except for the requirement of two high-$p_{\text{T}}$ wide jets.

The normalized distribution of $R_{\text{T}}$ vs. $m_{\text{T}}$ for the simulated QCD background. All selection requirements are omitted, except for the requirement of two high-$p_{\text{T}}$ wide jets.

The normalized distributions of the BDT input variable $\tau_{21}$ for the two highest $p_{\text{T}}$ wide jets from the simulated SM backgrounds and several signal models with varying $m_{\text{dark}}$ values. Each sample's jet $p_{\text{T}}$ distribution is weighted to match a reference distribution (see text). The last bin of each histogram includes the overflow events.

The normalized distributions of the BDT input variable $\tau_{32}$ for the two highest $p_{\text{T}}$ wide jets from the simulated SM backgrounds and several signal models with varying $m_{\text{dark}}$ values. Each sample's jet $p_{\text{T}}$ distribution is weighted to match a reference distribution (see text). The last bin of each histogram includes the overflow events.

The normalized distributions of the BDT input variable $N_{2}^{(1)}$ for the two highest $p_{\text{T}}$ wide jets from the simulated SM backgrounds and several signal models with varying $m_{\text{dark}}$ values. Each sample's jet $p_{\text{T}}$ distribution is weighted to match a reference distribution (see text). The last bin of each histogram includes the overflow events.

The normalized distributions of the BDT input variable $N_{3}^{(1)}$ for the two highest $p_{\text{T}}$ wide jets from the simulated SM backgrounds and several signal models with varying $m_{\text{dark}}$ values. Each sample's jet $p_{\text{T}}$ distribution is weighted to match a reference distribution (see text). The last bin of each histogram includes the overflow events.

The normalized distributions of the BDT input variable $g_{\text{jet}}$ for the two highest $p_{\text{T}}$ wide jets from the simulated SM backgrounds and several signal models with varying $r_{\text{inv}}$ values. Each sample's jet $p_{\text{T}}$ distribution is weighted to match a reference distribution (see text). The last bin of each histogram includes the overflow events.

The normalized distributions of the BDT input variable $\sigma_{\text{major}}$ for the two highest $p_{\text{T}}$ wide jets from the simulated SM backgrounds and several signal models with varying $r_{\text{inv}}$ values. Each sample's jet $p_{\text{T}}$ distribution is weighted to match a reference distribution (see text). The last bin of each histogram includes the overflow events.

The normalized distributions of the BDT input variable $\sigma_{\text{minor}}$ for the two highest $p_{\text{T}}$ wide jets from the simulated SM backgrounds and several signal models with varying $r_{\text{inv}}$ values. Each sample's jet $p_{\text{T}}$ distribution is weighted to match a reference distribution (see text). The last bin of each histogram includes the overflow events.

The normalized distributions of the BDT input variable $\Delta\phi(\vec{J},\vec{p}_{\text{T}}^{\text{miss}})$ for the two highest $p_{\text{T}}$ wide jets from the simulated SM backgrounds and several signal models with varying $r_{\text{inv}}$ values. Each sample's jet $p_{\text{T}}$ distribution is weighted to match a reference distribution (see text). The last bin of each histogram includes the overflow events.

The normalized distributions of the BDT input variable $f_{\text{h}^{\pm}}$ for the two highest $p_{\text{T}}$ wide jets from the simulated SM backgrounds and several signal models with varying $m_{\text{Z}^{\prime}}$ values. Each sample's jet $p_{\text{T}}$ distribution is weighted to match a reference distribution (see text). The last bin of each histogram includes the overflow events.

The normalized distributions of the BDT input variable $f_{\text{e}}$ for the two highest $p_{\text{T}}$ wide jets from the simulated SM backgrounds and several signal models with varying $m_{\text{Z}^{\prime}}$ values. Each sample's jet $p_{\text{T}}$ distribution is weighted to match a reference distribution (see text). The last bin of each histogram includes the overflow events.

The normalized distributions of the BDT input variable $f_{\mu}$ for the two highest $p_{\text{T}}$ wide jets from the simulated SM backgrounds and several signal models with varying $m_{\text{Z}^{\prime}}$ values. Each sample's jet $p_{\text{T}}$ distribution is weighted to match a reference distribution (see text). The last bin of each histogram includes the overflow events.

The normalized distributions of the BDT input variable $f_{\text{h}^{0}}$ for the two highest $p_{\text{T}}$ wide jets from the simulated SM backgrounds and several signal models with varying $m_{\text{Z}^{\prime}}$ values. Each sample's jet $p_{\text{T}}$ distribution is weighted to match a reference distribution (see text). The last bin of each histogram includes the overflow events.

The normalized distributions of the BDT input variable $f_{\gamma}$ for the two highest $p_{\text{T}}$ wide jets from the simulated SM backgrounds and several signal models with varying $m_{\text{Z}^{\prime}}$ values. Each sample's jet $p_{\text{T}}$ distribution is weighted to match a reference distribution (see text). The last bin of each histogram includes the overflow events.

The product of signal acceptance and efficiency in the high-$R_{\text{T}}$ signal region, for variations of the mediator mass and the dark hadron mass.

The product of signal acceptance and efficiency in the high-$R_{\text{T}}$ signal region, for variations of the mediator mass and the invisible fraction.

The product of signal acceptance and efficiency in the high-$R_{\text{T}}$ signal region, for variations of the mediator mass and the dark coupling strength.

The product of signal acceptance and efficiency in the low-$R_{\text{T}}$ signal region, for variations of the mediator mass and the dark hadron mass.

The product of signal acceptance and efficiency in the low-$R_{\text{T}}$ signal region, for variations of the mediator mass and the invisible fraction.

The product of signal acceptance and efficiency in the low-$R_{\text{T}}$ signal region, for variations of the mediator mass and the dark coupling strength.

The product of signal acceptance and efficiency in the high-SVJ2 signal region, for variations of the mediator mass and the dark hadron mass.

The product of signal acceptance and efficiency in the high-SVJ2 signal region, for variations of the mediator mass and the invisible fraction.

The product of signal acceptance and efficiency in the high-SVJ2 signal region, for variations of the mediator mass and the dark coupling strength.

The product of signal acceptance and efficiency in the low-SVJ2 signal region, for variations of the mediator mass and the dark hadron mass.

The product of signal acceptance and efficiency in the low-SVJ2 signal region, for variations of the mediator mass and the invisible fraction.

The product of signal acceptance and efficiency in the low-SVJ2 signal region, for variations of the mediator mass and the dark coupling strength.

Comparison of different the dijet mass $m_{\text{J}\text{J}}$, the transverse mass $m_{\text{T}}$, and the Monte Carlo (MC) mass $m_{\text{MC}}$ for a signal model with $m_{\text{Z}^{\prime}} = 2.1\,\text{TeV}$, $m_{\text{dark}} = 20\,\text{GeV}$, $r_{\text{inv}} = 0.3$, $\alpha_{\text{dark}} = \alpha_{\text{dark}}^{\text{peak}}$. No selection is applied, except that there must be at least two jets. $m_{\text{MC}}$ is computed by adding the generator-level four-vectors for invisible particles to the dijet system, to represent the achievable resolution if the invisible component were fully measured. The last bin of each histogram includes the overflow events.

$m_{\text{T}}$ distributions for signal models with different $m_{\text{dark}}$ values for the high-$R_{\text{T}}$ inclusive signal region.

$m_{\text{T}}$ distributions for signal models with different $m_{\text{dark}}$ values for the low-$R_{\text{T}}$ inclusive signal region.

$m_{\text{T}}$ distributions for signal models with different $r_{\text{inv}}$ values for the high-$R_{\text{T}}$ inclusive signal region.

$m_{\text{T}}$ distributions for signal models with different $r_{\text{inv}}$ values for the low-$R_{\text{T}}$ inclusive signal region.

$m_{\text{T}}$ distributions for signal models with different $\alpha_{\text{dark}}$ values for the high-$R_{\text{T}}$ inclusive signal region.

$m_{\text{T}}$ distributions for signal models with different $\alpha_{\text{dark}}$ values for the low-$R_{\text{T}}$ inclusive signal region.

The proportions of each SM background process in the high-$R_{\text{T}}$ signal region.

The proportions of each SM background process in the low-$R_{\text{T}}$ signal region.

The proportions of each SM background process in the high-SVJ2 signal region.

The proportions of each SM background process in the low-SVJ2 signal region.

The 95% CL expected upper limits on the product of the cross section and branching fraction from the inclusive search for variations of the mediator mass and the dark hadron mass.

The 95% CL expected upper limits on the product of the cross section and branching fraction from the inclusive search for variations of the mediator mass and the invisible fraction.

The 95% CL expected upper limits on the product of the cross section and branching fraction from the BDT-based search for variations of the mediator mass and the dark hadron mass.

The 95% CL expected upper limits on the product of the cross section and branching fraction from the BDT-based search for variations of the mediator mass and the invisible fraction.

Relative efficiencies in % for each step of the event selection process for the major background processes. Statistical uncertainties, at most 1.8%, are omitted. The line "Efficiency [%]" is the absolute efficiency after the final selection. The subsequent lines show the efficiency for each signal region, relative to the final selection.

Relative efficiencies in % for each step of the event selection process for signals with $m_{\text{Z}^{\prime}} = 2.1\,\text{TeV}$, varying $m_{\text{dark}}$ values, $r_{\text{inv}} = 0.3$, and $\alpha_{\text{dark}} = \alpha_{\text{dark}}^{\text{peak}}$. Statistical uncertainties, at most 0.5%, are omitted. The line "Efficiency [%]" is the absolute efficiency after the final selection. The subsequent lines show the efficiency for each signal region, relative to the final selection.

Relative efficiencies in % for each step of the event selection process for signals with $m_{\text{Z}^{\prime}} = 3.1\,\text{TeV}$, varying $m_{\text{dark}}$ values, $r_{\text{inv}} = 0.3$, and $\alpha_{\text{dark}} = \alpha_{\text{dark}}^{\text{peak}}$. Statistical uncertainties, at most 0.4%, are omitted. The line "Efficiency [%]" is the absolute efficiency after the final selection. The subsequent lines show the efficiency for each signal region, relative to the final selection.

Relative efficiencies in % for each step of the event selection process for signals with $m_{\text{Z}^{\prime}} = 4.1\,\text{TeV}$, varying $m_{\text{dark}}$ values, $r_{\text{inv}} = 0.3$, and $\alpha_{\text{dark}} = \alpha_{\text{dark}}^{\text{peak}}$. Statistical uncertainties, at most 0.4%, are omitted. The line "Efficiency [%]" is the absolute efficiency after the final selection. The subsequent lines show the efficiency for each signal region, relative to the final selection.

Relative efficiencies in % for each step of the event selection process for signals with $m_{\text{Z}^{\prime}} = 2.1\,\text{TeV}$, $m_{\text{dark}} = 20\,\text{GeV}$, varying $r_{\text{inv}}$ values, and $\alpha_{\text{dark}} = \alpha_{\text{dark}}^{\text{peak}}$. Statistical uncertainties, at most 2.6%, are omitted. The line "Efficiency [%]" is the absolute efficiency after the final selection. The subsequent lines show the efficiency for each signal region, relative to the final selection.

Relative efficiencies in % for each step of the event selection process for signals with $m_{\text{Z}^{\prime}} = 3.1\,\text{TeV}$, $m_{\text{dark}} = 20\,\text{GeV}$, varying $r_{\text{inv}}$ values, and $\alpha_{\text{dark}} = \alpha_{\text{dark}}^{\text{peak}}$. Statistical uncertainties, at most 1.2%, are omitted. The line "Efficiency [%]" is the absolute efficiency after the final selection. The subsequent lines show the efficiency for each signal region, relative to the final selection.

Relative efficiencies in % for each step of the event selection process for signals with $m_{\text{Z}^{\prime}} = 4.1\,\text{TeV}$, $m_{\text{dark}} = 20\,\text{GeV}$, varying $r_{\text{inv}}$ values, and $\alpha_{\text{dark}} = \alpha_{\text{dark}}^{\text{peak}}$. Statistical uncertainties, at most 0.9%, are omitted. The line "Efficiency [%]" is the absolute efficiency after the final selection. The subsequent lines show the efficiency for each signal region, relative to the final selection.

Relative efficiencies in % for each step of the event selection process for signals with $m_{\text{Z}^{\prime}} = 2.1\,\text{TeV}$, $m_{\text{dark}} = 20\,\text{GeV}$, $r_{\text{inv}} = 0.3$, and varying $\alpha_{\text{dark}}$ values. Statistical uncertainties, at most 0.4%, are omitted. The line "Efficiency [%]" is the absolute efficiency after the final selection. The subsequent lines show the efficiency for each signal region, relative to the final selection.

Relative efficiencies in % for each step of the event selection process for signals with $m_{\text{Z}^{\prime}} = 3.1\,\text{TeV}$, $m_{\text{dark}} = 20\,\text{GeV}$, $r_{\text{inv}} = 0.3$, and varying $\alpha_{\text{dark}}$ values. Statistical uncertainties, at most 0.4%, are omitted. The line "Efficiency [%]" is the absolute efficiency after the final selection. The subsequent lines show the efficiency for each signal region, relative to the final selection.

Relative efficiencies in % for each step of the event selection process for signals with $m_{\text{Z}^{\prime}} = 4.1\,\text{TeV}$, $m_{\text{dark}} = 20\,\text{GeV}$, $r_{\text{inv}} = 0.3$, and varying $\alpha_{\text{dark}}$ values. Statistical uncertainties, at most 0.4%, are omitted. The line "Efficiency [%]" is the absolute efficiency after the final selection. The subsequent lines show the efficiency for each signal region, relative to the final selection.

Pre-fit distribution of the reconstructed Higgs boson candidate $p_T^H$ for the dilepton $SR^{\geq 4j}_{\geq 4b}$ signal region. The dashed line shows the ${t\bar {t}H}$ signal distribution normalised to the total background prediction. The uncertainty band includes all uncertainties and their correlations, except for the uncertainty in the $k({t\bar {t}+{\geq }1b})$ normalisation factor which is not defined pre-fit. The last bin includes the overflow.

Pre-fit distribution of the reconstructed Higgs boson candidate $p_T^H$ for the single-lepton resolved $SR^{\geq 6j}_{\geq 4b}$ signal region. The dashed line shows the ${t\bar {t}H}$ signal distribution normalised to the total background prediction. The uncertainty band includes all uncertainties and their correlations, except for the uncertainty in the $k({t\bar {t}+{\geq }1b})$ normalisation factor which is not defined pre-fit. The last bin includes the overflow.

Pre-fit distribution of the reconstructed Higgs boson candidate $p_T^H$ for the single-lepton boosted ${{SR}_{{boosted}}}$ signal region. The dashed line shows the ${t\bar {t}H}$ signal distribution normalised to the total background prediction. The uncertainty band includes all uncertainties and their correlations, except for the uncertainty in the $k({t\bar {t}+{\geq }1b})$ normalisation factor which is not defined pre-fit. The last bin includes the overflow.

Comparison of predicted and observed event yields in each of the control and signal regions in the dilepton channel after the fit to the data. The uncertainty band includes all uncertainties and their correlations.

Comparison of predicted and observed event yields in each of the control and signal regions in the single-lepton channels after the fit to the data. The uncertainty band includes all uncertainties and their correlations.

Comparison between data and prediction for the BDT discriminant in the dilepton SRs after the inclusive fit to the data for $0\le p_T^H<120$ GeV. The ${t\bar {t}H}$ signal yield (solid red) is normalised to the fitted $\mu $ value from the inclusive fit. The dashed line shows the ${t\bar {t}H}$ signal distribution normalised to the total background prediction. The uncertainty band includes all uncertainties and their correlations.

Comparison between data and prediction for the BDT discriminant in the dilepton SRs after the inclusive fit to the data for $120\le p_T^H<200$ GeV. The ${t\bar {t}H}$ signal yield (solid red) is normalised to the fitted $\mu $ value from the inclusive fit. The dashed line shows the ${t\bar {t}H}$ signal distribution normalised to the total background prediction. The uncertainty band includes all uncertainties and their correlations.

Comparison between data and prediction for the BDT discriminant in the dilepton SRs after the inclusive fit to the data for $200\le p_T^H<300$ GeV. The ${t\bar {t}H}$ signal yield (solid red) is normalised to the fitted $\mu $ value from the inclusive fit. The dashed line shows the ${t\bar {t}H}$ signal distribution normalised to the total background prediction. The uncertainty band includes all uncertainties and their correlations.

Comparison between data and prediction for the BDT discriminant in the dilepton SRs after the inclusive fit to the data for $p_{{T}}^{H}\ge 300$ GeV. The ${t\bar {t}H}$ signal yield (solid red) is normalised to the fitted $\mu $ value from the inclusive fit. The dashed line shows the ${t\bar {t}H}$ signal distribution normalised to the total background prediction. The uncertainty band includes all uncertainties and their correlations.

Comparison between data and prediction for the BDT discriminant in the single-lepton resolved SRs after the inclusive fit to the data for $0\le p_T^H<120$ GeV. The ${t\bar {t}H}$ signal yield (solid red) is normalised to the fitted $\mu $ value from the inclusive fit. The dashed line shows the ${t\bar {t}H}$ signal distribution normalised to the total background prediction. The uncertainty band includes all uncertainties and their correlations.

Comparison between data and prediction for the BDT discriminant in the single-lepton resolved SRs after the inclusive fit to the data for $120\le p_T^H<200$ GeV. The ${t\bar {t}H}$ signal yield (solid red) is normalised to the fitted $\mu $ value from the inclusive fit. The dashed line shows the ${t\bar {t}H}$ signal distribution normalised to the total background prediction. The uncertainty band includes all uncertainties and their correlations.

Comparison between data and prediction for the BDT discriminant in the single-lepton resolved SRs after the inclusive fit to the data for $200\le p_T^H<300$ GeV. The ${t\bar {t}H}$ signal yield (solid red) is normalised to the fitted $\mu $ value from the inclusive fit. The dashed line shows the ${t\bar {t}H}$ signal distribution normalised to the total background prediction. The uncertainty band includes all uncertainties and their correlations.

Comparison between data and prediction for the BDT discriminant in the single-lepton resolved SRs after the inclusive fit to the data for $300\le p_T^H<450$ GeV. The ${t\bar {t}H}$ signal yield (solid red) is normalised to the fitted $\mu $ value from the inclusive fit. The dashed line shows the ${t\bar {t}H}$ signal distribution normalised to the total background prediction. The uncertainty band includes all uncertainties and their correlations.

Comparison between data and prediction for the BDT discriminant in the single-lepton resolved SRs after the inclusive fit to the data for $p_{{T}}^{H}\ge 450$ GeV (yield only). The ${t\bar {t}H}$ signal yield (solid red) is normalised to the fitted $\mu $ value from the inclusive fit. The dashed line shows the ${t\bar {t}H}$ signal distribution normalised to the total background prediction. The uncertainty band includes all uncertainties and their correlations.

Comparison between data and prediction for the BDT discriminant in the single-lepton boosted SRs after the inclusive fit to the data for $300\le p_T^H<450$ GeV. The ${t\bar {t}H}$ signal yield (solid red) is normalised to the fitted $\mu $ value from the inclusive fit. The dashed line shows the ${t\bar {t}H}$ signal distribution normalised to the total background prediction. The uncertainty band includes all uncertainties and their correlations.

Comparison between data and prediction for the BDT discriminant in the single-lepton boosted SRs after the inclusive fit to the data for $p_{{T}}^{H}\ge 450$ GeV. The ${t\bar {t}H}$ signal yield (solid red) is normalised to the fitted $\mu $ value from the inclusive fit. The dashed line shows the ${t\bar {t}H}$ signal distribution normalised to the total background prediction. The uncertainty band includes all uncertainties and their correlations.

Comparison between data and prediction for ${\Delta R^{{avg}}_{bb}}$ after the inclusive fit to the data in the single-lepton $CR^{5j}_{{\geq}4b\ lo}$ control region. The ${t\bar {t}H}$ signal yield (solid red) is normalised to the fitted $\mu $ value from the inclusive fit. The dashed line shows the ${t\bar {t}H}$ signal distribution normalised to the total background prediction. The uncertainty band includes all uncertainties and their correlations. The first (last) bin includes the underflow (overflow).

Comparison between data and prediction for ${\Delta R^{{avg}}_{bb}}$ after the inclusive fit to the data in the single-lepton $CR^{5j}_{{\geq}4b\ hi}$ control region. The ${t\bar {t}H}$ signal yield (solid red) is normalised to the fitted $\mu $ value from the inclusive fit. The dashed line shows the ${t\bar {t}H}$ signal distribution normalised to the total background prediction. The uncertainty band includes all uncertainties and their correlations. The first (last) bin includes the underflow (overflow).

Post-fit yields of signal ($S$) and total background ($B$) as a function of $\log (S/B)$, compared with data. Final-discriminant bins in all dilepton and single-lepton analysis regions are combined into bins of $\log (S/B)$, with the signal normalised to the SM prediction used for the computation of $\log (S/B)$. The signal is then shown normalised to the best-fit value and the SM prediction. The lower frame reports the ratio of data to background, and this is compared with the expected ${t\bar {t}H}$-signal-plus-background yield divided by the background-only yield for the best-fit signal strength (solid red line) and the SM prediction (dashed orange line).

Comparison between data and prediction for the reconstruction BDT score for the Higgs boson candidate identified using Higgs boson information, after the inclusive fit to the data in the dilepton resolved channel for $0\le p_T^H<120$ GeV. The ${t\bar {t}H}$ signal yield (solid red) is normalised to the fitted $\mu $ value from the inclusive fit. The dashed line shows the ${t\bar {t}H}$ signal distribution normalised to the total background prediction. The uncertainty band includes all uncertainties and their correlations.

Comparison between data and prediction for the average $\Delta \eta $ between $b$-tagged jets, after the inclusive fit to the data in the dilepton resolved channel for $0\le p_T^H<120$ GeV. The ${t\bar {t}H}$ signal yield (solid red) is normalised to the fitted $\mu $ value from the inclusive fit. The dashed line shows the ${t\bar {t}H}$ signal distribution normalised to the total background prediction. The uncertainty band includes all uncertainties and their correlations.

Comparison between data and prediction for the likelihood discriminant, after the inclusive fit to the data in the single-lepton resolved channel for $0\le p_T^H<120$ GeV. The ${t\bar {t}H}$ signal yield (solid red) is normalised to the fitted $\mu $ value from the inclusive fit. The dashed line shows the ${t\bar {t}H}$ signal distribution normalised to the total background prediction. The uncertainty band includes all uncertainties and their correlations.

Comparison between data and prediction for the average $\Delta R$ for all possible combinations of $b$-tagged jet pairs, after the inclusive fit to the data in the single-lepton resolved channel for $0\le p_T^H<120$ GeV. The ${t\bar {t}H}$ signal yield (solid red) is normalised to the fitted $\mu $ value from the inclusive fit. The dashed line shows the ${t\bar {t}H}$ signal distribution normalised to the total background prediction. The uncertainty band includes all uncertainties and their correlations.

Comparison between data and prediction for the DNN $P(H)$ output for the Higgs boson candidate after the inclusive fit to the data in the single-lepton boosted channel for $300\le p_T^H<450$ GeV. The ${t\bar {t}H}$ signal yield (solid red) is normalised to the fitted $\mu $ value from the inclusive fit. The dashed line shows the ${t\bar {t}H}$ signal distribution normalised to the total background prediction. The uncertainty band includes all uncertainties and their correlations.

Comparison between data and prediction for the DNN $P(H)$ output for the Higgs boson candidate after the inclusive fit to the data in the single-lepton boosted channel for $p_{{T}}^{H}\ge 450$ GeV. The ${t\bar {t}H}$ signal yield (solid red) is normalised to the fitted $\mu $ value from the inclusive fit. The dashed line shows the ${t\bar {t}H}$ signal distribution normalised to the total background prediction. The uncertainty band includes all uncertainties and their correlations.

Post-fit distribution of the reconstructed Higgs boson candidate mass for the dilepton $SR^{\geq 4j}_{\geq 4b}$ signal region. The ${t\bar {t}H}$ signal yield (solid red) is normalised to the fitted $\mu $ value from the inclusive fit. The dashed line shows the ${t\bar {t}H}$ signal distribution normalised to the total background prediction. The uncertainty band includes all uncertainties and their correlations. The first (last) bin includes the underflow (overflow).

Post-fit distribution of the reconstructed Higgs boson candidate mass for the single-lepton resolved $SR^{\geq 6j}_{\geq 4b}$ signal region. The ${t\bar {t}H}$ signal yield (solid red) is normalised to the fitted $\mu $ value from the inclusive fit. The dashed line shows the ${t\bar {t}H}$ signal distribution normalised to the total background prediction. The uncertainty band includes all uncertainties and their correlations. The first (last) bin includes the underflow (overflow).

Post-fit distribution of the reconstructed Higgs boson candidate mass for the single-lepton boosted ${{SR}_{{boosted}}}$ signal region. The ${t\bar {t}H}$ signal yield (solid red) is normalised to the fitted $\mu $ value from the inclusive fit. The dashed line shows the ${t\bar {t}H}$ signal distribution normalised to the total background prediction. The uncertainty band includes all uncertainties and their correlations. The first (last) bin includes the underflow (overflow).

Fitted values of the ${t\bar {t}H}$ signal strength parameter in the individual channels and in the inclusive signal-strength measurement.

Ranking of the 20 nuisance parameters with the largest post-fit impact on $\mu $ in the fit. Nuisance parameters corresponding to statistical uncertainties in the simulated event samples are not included. The empty blue rectangles correspond to the pre-fit impact on $\mu $ and the filled blue ones to the post-fit impact on $\mu $, both referring to the upper scale. The impact of each nuisance parameter, $\Delta \mu $, is computed by comparing the nominal best-fit value of $\mu $ with the result of the fit when fixing the considered nuisance parameter to its best-fit value, $\hat{\theta }$, shifted by its pre-fit (post-fit) uncertainties $\pm \Delta \theta $ ($\pm \Delta \hat{\theta }$). The black points show the pulls of the nuisance parameters relative to their nominal values, $\theta _0$. These pulls and their relative post-fit errors, $\Delta \hat{\theta }/\Delta \theta $, refer to the lower scale. The `ljets' (`dilep') label refers to the single-lepton (dilepton) channel.

Pre-fit distribution of the number of jets in the dilepton $SR^{\geq 4j}_{\geq 4b}$ signal region. The ${t\bar {t}H}$ signal yield (solid red) is normalised to the Standard Model expectation. The dashed line shows the ${t\bar {t}H}$ signal distribution normalised to the total background prediction. The uncertainty band includes all uncertainties and their correlations, except the uncertainty in the $k({t\bar {t}+{\geq }1b})$ normalisation factor that is not defined pre-fit.

Pre-fit distribution of the number of jets in the single-lepton resolved $SR^{\geq 6j}_{\geq 4b}$ signal region. The ${t\bar {t}H}$ signal yield (solid red) is normalised to the Standard Model expectation. The dashed line shows the ${t\bar {t}H}$ signal distribution normalised to the total background prediction. The uncertainty band includes all uncertainties and their correlations, except the uncertainty in the $k({t\bar {t}+{\geq }1b})$ normalisation factor that is not defined pre-fit.

Pre-fit distribution of the number of jets in the single-lepton boosted ${{SR}_{{boosted}}}$ signal region. The ${t\bar {t}H}$ signal yield (solid red) is normalised to the Standard Model expectation. The dashed line shows the ${t\bar {t}H}$ signal distribution normalised to the total background prediction. The uncertainty band includes all uncertainties and their correlations, except the uncertainty in the $k({t\bar {t}+{\geq }1b})$ normalisation factor that is not defined pre-fit.

Post-fit distribution of the number of jets in the dilepton $SR^{\geq 4j}_{\geq 4b}$ signal region. The ${t\bar {t}H}$ signal yield (solid red) is normalised to the fitted $\mu $ value from the inclusive fit. The dashed line shows the ${t\bar {t}H}$ signal distribution normalised to the total background prediction. The uncertainty band includes all uncertainties and their correlations.

Post-fit distribution of the number of jets in the single-lepton resolved $SR^{\geq 6j}_{\geq 4b}$ signal region. The ${t\bar {t}H}$ signal yield (solid red) is normalised to the fitted $\mu $ value from the inclusive fit. The dashed line shows the ${t\bar {t}H}$ signal distribution normalised to the total background prediction. The uncertainty band includes all uncertainties and their correlations.

Post-fit distribution of the number of jets in the single-lepton boosted ${{SR}_{{boosted}}}$ signal region. The ${t\bar {t}H}$ signal yield (solid red) is normalised to the fitted $\mu $ value from the inclusive fit. The dashed line shows the ${t\bar {t}H}$ signal distribution normalised to the total background prediction. The uncertainty band includes all uncertainties and their correlations.

Post-fit distribution of the reconstructed Higgs boson candidate $p_T^H$ for the dilepton $SR^{\geq 4j}_{\geq 4b}$ signal region. The ${t\bar {t}H}$ signal yield (solid red) is normalised to the fitted $\mu $ value from the inclusive fit. The dashed line shows the ${t\bar {t}H}$ signal distribution normalised to the total background prediction. The uncertainty band includes all uncertainties and their correlations. The last bin includes the overflow.

Post-fit distribution of the reconstructed Higgs boson candidate $p_T^H$ for the single-lepton resolved $SR^{\geq 6j}_{\geq 4b}$ signal region. The ${t\bar {t}H}$ signal yield (solid red) is normalised to the fitted $\mu $ value from the inclusive fit. The dashed line shows the ${t\bar {t}H}$ signal distribution normalised to the total background prediction. The uncertainty band includes all uncertainties and their correlations. The last bin includes the overflow.

Post-fit distribution of the reconstructed Higgs boson candidate $p_T^H$ for the single-lepton boosted ${{SR}_{{boosted}}}$ signal region. The ${t\bar {t}H}$ signal yield (solid red) is normalised to the fitted $\mu $ value from the inclusive fit. The dashed line shows the ${t\bar {t}H}$ signal distribution normalised to the total background prediction. The uncertainty band includes all uncertainties and their correlations. The last bin includes the overflow.

Signal-strength measurements in the individual STXS ${\hat{p}_{{T}}^{H}}$ bins, as well as the inclusive signal strength.

95% CL simplified template cross-section upper limits in the individual STXS ${\hat{p}_{{T}}^{H}}$ bins, as well as the inclusive limit. The observed limits are shown (solid black lines), together with the expected limits both in the background-only hypothesis (dotted black lines) and in the SM hypothesis (dotted red lines). In the case of the expected limits in the background-only hypothesis, one- and two-standard-deviation uncertainty bands are also shown. The hatched uncertainty bands correspond to the theory uncertainty in the fiducial cross-section prediction in each bin.

The ratios $S/B$ (black solid line, referring to the vertical axis on the left) and $S/\sqrt{B}$ (red dashed line, referring to the vertical axis on the right) for each category in the inclusive analysis in the dilepton channel (left) and in the single-lepton channels (right), where $S$ ($B$) is the number of selected signal (background) events predicted by the simulation and normalised to a luminosity of 139 fb$^{-1}$ .

Comparison between data and prediction for the $\Delta R$ between the Higgs candidate and the ${t\bar {t}}$ candidate system, after the inclusive fit to the data in the dilepton resolved channel for $0\le p_T^H<120$ GeV. The ${t\bar {t}H}$ signal yield (solid red) is normalised to the fitted $\mu $ value from the inclusive fit. The dashed line shows the ${t\bar {t}H}$ signal distribution normalised to the total background prediction. The uncertainty band includes all uncertainties and their correlations.

Comparison between data and prediction for the number of $b$-tagged jet pairs with an invariant mass within 30 GeV of 125 GeV, after the inclusive fit to the data in the dilepton resolved channel for $0\le p_T^H<120$ GeV. The ${t\bar {t}H}$ signal yield (solid red) is normalised to the fitted $\mu $ value from the inclusive fit. The dashed line shows the ${t\bar {t}H}$ signal distribution normalised to the total background prediction. The uncertainty band includes all uncertainties and their correlations.

Comparison between data and prediction for the reconstruction BDT score for the Higgs boson candidate identified using Higgs boson information, after the inclusive fit to the data in the single-lepton resolved channel for $0\le p_T^H<120$ GeV. The ${t\bar {t}H}$ signal yield (solid red) is normalised to the fitted $\mu $ value from the inclusive fit. The dashed line shows the ${t\bar {t}H}$ signal distribution normalised to the total background prediction. The uncertainty band includes all uncertainties and their correlations.

Comparison between data and prediction for the $\Delta R$ between the two highest ${p_{{T}}}$ $b$-tagged jets, after the inclusive fit to the data in the single-lepton resolved channel for $0\le p_T^H<120$ GeV. The ${t\bar {t}H}$ signal yield (solid red) is normalised to the fitted $\mu $ value from the inclusive fit. The dashed line shows the ${t\bar {t}H}$ signal distribution normalised to the total background prediction. The uncertainty band includes all uncertainties and their correlations.

Comparison between data and prediction for the sum of $b$-tagging discriminants of jets from Higgs, hadronic top and leptonic top candidates, after the inclusive fit to the data in the single-lepton boosted channel for $300\le p_T^H<450$ GeV. The ${t\bar {t}H}$ signal yield (solid red) is normalised to the fitted $\mu $ value from the inclusive fit. The dashed line shows the ${t\bar {t}H}$ signal distribution normalised to the total background prediction. The uncertainty band includes all uncertainties and their correlations. The first (last) bin includes the underflow (overflow).

Comparison between data and prediction for the sum of $b$-tagging discriminants of jets from Higgs, hadronic top and leptonic top candidates, after the inclusive fit to the data in the single-lepton boosted channel for $p_{{T}}^{H}\ge 450$ GeV. The ${t\bar {t}H}$ signal yield (solid red) is normalised to the fitted $\mu $ value from the inclusive fit. The dashed line shows the ${t\bar {t}H}$ signal distribution normalised to the total background prediction. The uncertainty band includes all uncertainties and their correlations. The first (last) bin includes the underflow (overflow).

Comparison between data and prediction for the hadronic top candidate invariant mass, after the inclusive fit to the data in the single-lepton boosted channel for $300\le p_T^H<450$ GeV. The ${t\bar {t}H}$ signal yield (solid red) is normalised to the fitted $\mu $ value from the inclusive fit. The dashed line shows the ${t\bar {t}H}$ signal distribution normalised to the total background prediction. The uncertainty band includes all uncertainties and their correlations. The first (last) bin includes the underflow (overflow).

Comparison between data and prediction for the hadronic top candidate invariant mass, after the inclusive fit to the data in the single-lepton boosted channel for $p_{{T}}^{H}\ge 450$ GeV. The ${t\bar {t}H}$ signal yield (solid red) is normalised to the fitted $\mu $ value from the inclusive fit. The dashed line shows the ${t\bar {t}H}$ signal distribution normalised to the total background prediction. The uncertainty band includes all uncertainties and their correlations. The first (last) bin includes the underflow (overflow).

Comparison between data and prediction for the fraction of the sum of $b$-tagging discriminants of all jets not associated to the Higgs or hadronic top candidates, after the inclusive fit to the data in the single-lepton boosted channel for $300\le p_T^H<450$ GeV. The ${t\bar {t}H}$ signal yield (solid red) is normalised to the fitted $\mu $ value from the inclusive fit. The dashed line shows the ${t\bar {t}H}$ signal distribution normalised to the total background prediction. The uncertainty band includes all uncertainties and their correlations. The first (last) bin includes the underflow (overflow).

Comparison between data and prediction for the fraction of the sum of $b$-tagging discriminants of all jets not associated to the Higgs or hadronic top candidates, after the inclusive fit to the data in the single-lepton boosted channel for $p_{{T}}^{H}\ge 450$ GeV. The ${t\bar {t}H}$ signal yield (solid red) is normalised to the fitted $\mu $ value from the inclusive fit. The dashed line shows the ${t\bar {t}H}$ signal distribution normalised to the total background prediction. The uncertainty band includes all uncertainties and their correlations. The first (last) bin includes the underflow (overflow).

Ranking of the 20 nuisance parameters with the largest post-fit impact on $\mu $ in the STXS fit for $0\le {\hat{p}_{{T}}^{H}}<120$ GeV. Nuisance parameters corresponding to statistical uncertainties in the simulated event samples are not included. The empty blue rectangles correspond to the pre-fit impact on $\mu $ and the filled blue ones to the post-fit impact on $\mu $, both referring to the upper scale. The impact of each nuisance parameter, $\Delta \mu $, is computed by comparing the nominal best-fit value of $\mu $ with the result of the fit when fixing the considered nuisance parameter to its best-fit value, $\hat{\theta }$, shifted by its pre-fit (post-fit) uncertainties $\pm \Delta \theta $ ($\pm \Delta \hat{\theta }$). The black points show the pulls of the nuisance parameters relative to their nominal values, $\theta _0$. These pulls and their relative post-fit errors, $\Delta \hat{\theta }/\Delta \theta $, refer to the lower scale. For experimental uncertainties that are decomposed into several independent sources, NP X corresponds to the X$^{th}$ nuisance parameter, ordered by their impact on $\mu $. The `ljets' (`dilep') label refers to the single-lepton (dilepton) channel.

Ranking of the 20 nuisance parameters with the largest post-fit impact on $\mu $ in the STXS fit for $120\le {\hat{p}_{{T}}^{H}}<200$ GeV. Nuisance parameters corresponding to statistical uncertainties in the simulated event samples are not included. The empty blue rectangles correspond to the pre-fit impact on $\mu $ and the filled blue ones to the post-fit impact on $\mu $, both referring to the upper scale. The impact of each nuisance parameter, $\Delta \mu $, is computed by comparing the nominal best-fit value of $\mu $ with the result of the fit when fixing the considered nuisance parameter to its best-fit value, $\hat{\theta }$, shifted by its pre-fit (post-fit) uncertainties $\pm \Delta \theta $ ($\pm \Delta \hat{\theta }$). The black points show the pulls of the nuisance parameters relative to their nominal values, $\theta _0$. These pulls and their relative post-fit errors, $\Delta \hat{\theta }/\Delta \theta $, refer to the lower scale. For experimental uncertainties that are decomposed into several independent sources, NP X corresponds to the X$^{th}$ nuisance parameter, ordered by their impact on $\mu $. The `ljets' (`dilep') label refers to the single-lepton (dilepton) channel.

Ranking of the 20 nuisance parameters with the largest post-fit impact on $\mu $ in the STXS fit for $200\le {\hat{p}_{{T}}^{H}}<300$ GeV. Nuisance parameters corresponding to statistical uncertainties in the simulated event samples are not included. The empty blue rectangles correspond to the pre-fit impact on $\mu $ and the filled blue ones to the post-fit impact on $\mu $, both referring to the upper scale. The impact of each nuisance parameter, $\Delta \mu $, is computed by comparing the nominal best-fit value of $\mu $ with the result of the fit when fixing the considered nuisance parameter to its best-fit value, $\hat{\theta }$, shifted by its pre-fit (post-fit) uncertainties $\pm \Delta \theta $ ($\pm \Delta \hat{\theta }$). The black points show the pulls of the nuisance parameters relative to their nominal values, $\theta _0$. These pulls and their relative post-fit errors, $\Delta \hat{\theta }/\Delta \theta $, refer to the lower scale. For experimental uncertainties that are decomposed into several independent sources, NP X corresponds to the X$^{th}$ nuisance parameter, ordered by their impact on $\mu $. The `ljets' (`dilep') label refers to the single-lepton (dilepton) channel.

Ranking of the 20 nuisance parameters with the largest post-fit impact on $\mu $ in the STXS fit for $300\le {\hat{p}_{{T}}^{H}}<450$ GeV. Nuisance parameters corresponding to statistical uncertainties in the simulated event samples are not included. The empty blue rectangles correspond to the pre-fit impact on $\mu $ and the filled blue ones to the post-fit impact on $\mu $, both referring to the upper scale. The impact of each nuisance parameter, $\Delta \mu $, is computed by comparing the nominal best-fit value of $\mu $ with the result of the fit when fixing the considered nuisance parameter to its best-fit value, $\hat{\theta }$, shifted by its pre-fit (post-fit) uncertainties $\pm \Delta \theta $ ($\pm \Delta \hat{\theta }$). The black points show the pulls of the nuisance parameters relative to their nominal values, $\theta _0$. These pulls and their relative post-fit errors, $\Delta \hat{\theta }/\Delta \theta $, refer to the lower scale. For experimental uncertainties that are decomposed into several independent sources, NP X corresponds to the X$^{th}$ nuisance parameter, ordered by their impact on $\mu $. The `ljets' (`dilep') label refers to the single-lepton (dilepton) channel.

Ranking of the 20 nuisance parameters with the largest post-fit impact on $\mu $ in the STXS fit for ${\hat{p}_{{T}}^{H}}\ge 450$ GeV. Nuisance parameters corresponding to statistical uncertainties in the simulated event samples are not included. The empty blue rectangles correspond to the pre-fit impact on $\mu $ and the filled blue ones to the post-fit impact on $\mu $, both referring to the upper scale. The impact of each nuisance parameter, $\Delta \mu $, is computed by comparing the nominal best-fit value of $\mu $ with the result of the fit when fixing the considered nuisance parameter to its best-fit value, $\hat{\theta }$, shifted by its pre-fit (post-fit) uncertainties $\pm \Delta \theta $ ($\pm \Delta \hat{\theta }$). The black points show the pulls of the nuisance parameters relative to their nominal values, $\theta _0$. These pulls and their relative post-fit errors, $\Delta \hat{\theta }/\Delta \theta $, refer to the lower scale. For experimental uncertainties that are decomposed into several independent sources, NP X corresponds to the X$^{th}$ nuisance parameter, ordered by their impact on $\mu $. The `ljets' (`dilep') label refers to the single-lepton (dilepton) channel.

95% confidence level upper limits on signal-strength measurements in the individual STXS ${\hat{p}_{{T}}^{H}}$ bins, as well as the inclusive signal-strength limit, after the fit used to extract multiple signal-strength parameters. The observed limits are shown (solid black lines), together with the expected limits both in the background-only hypothesis (dotted black lines) and in the SM hypothesis (dotted red lines). In the case of the expected limits in the background-only hypothesis, one- and two-standard-deviation uncertainty bands are also shown.

Post-fit correlation matrix (in percentages) between the $\mu $ values obtained in the STXS bins.

Performance of the Higgs boson reconstruction algorithms. For each row of `truth' ${\hat{p}_{{T}}^{H}}$, the matrix shows (in percentages) the fraction of Higgs boson candidates which are truth-matched to ${b\bar {b}}$ decays, with reconstructed $p_T^H$ in the various bins of the dilepton (left), single lepton resolved (middle) and boosted (right) channels.

Pre-fit event yields in the dilepton signal regions and control regions. All uncertainties are included except the $k({t\bar {t}+{\geq }1b})$ uncertainty that is not defined pre-fit. For the ${t\bar {t}H}$ signal, the pre-fit yield values correspond to the theoretical prediction and corresponding uncertainties. `Other sources' refers to s-channel, t-channel, $tW$, $tWZ$, $tZq$, $Z+$ jets and diboson events.

Post-fit event yields in the dilepton signal regions and control regions, after the inclusive fit in all channels. All uncertainties are included, taking into account correlations. For the ${t\bar {t}H}$ signal, the post-fit yield and uncertainties correspond to those in the inclusive signal-strength measurement. `Other sources' refers to s-channel, t-channel, $tW$, $tWZ$, $tZq$, $Z+$ jets and diboson events.

Pre-fit event yields in the single-lepton resolved and boosted signal regions and control regions. All uncertainties are included except the $k({t\bar {t}+{\geq }1b})$ uncertainty that is not defined pre-fit. For the ${t\bar {t}H}$ signal, the pre-fit yield values correspond to the theoretical prediction and corresponding uncertainties. `Other top sources' refers to s-channel, t-channel, $tWZ$ and $tZq$ events.

Post-fit event yields in the single-lepton resolved and boosted signal regions and control regions, after the inclusive fit in all channels. All uncertainties are included, taking into account correlations. For the ${t\bar {t}H}$ signal, the post-fit yield and uncertainties correspond to those in the inclusive signal-strength measurement. `Other top sources' refers to s-channel, t-channel, $tWZ$ and $tZq$ events.

Breakdown of the contributions to the uncertainties in $\mu$. The contributions from the different sources of uncertainty are evaluated after the fit. The $\Delta \mu $ values are obtained by repeating the fit after having fixed a certain set of nuisance parameters corresponding to a group of systematic uncertainties, and then evaluating $(\Delta \mu)^2$ by subtracting the resulting squared uncertainty of $\mu $ from its squared uncertainty found in the full fit. The same procedure is followed when quoting the effect of the ${t\bar {t}+{\geq }1b}$ normalisation. The total uncertainty is different from the sum in quadrature of the different components due to correlations between nuisance parameters existing in the fit.

Fraction (in percentages) of signal events, after SR and CR selections, originating from $b\bar {b}$, $WW$ and other remaining Higgs boson decay modes in the dilepton channel.

Fraction (in percentages) of signal events, after SR and CR selections, originating from $b\bar {b}$, $WW$ and other remaining Higgs boson decay modes in the single-lepton channels.

Predicted SM ${t\bar {t}H}$ cross-section in each of the five STXS ${\hat{p}_{{T}}^{H}}$ bins and signal acceptance times efficiency (including all event selection criteria) in each STXS bin as well as for the inclusive ${\hat{p}_{{T}}^{H}}$ range.

Number of expected signal events before the fit, after each selection requirement applied to enter the dilepton channel $SR^{\geq 4j}_{\geq 4b}$ region. All ${t\bar {t}H}$ signal events are included, regardless of the $H$ or ${t\bar {t}H}$ decay mode. All object corrections are applied, except for the initial number of events which is calculated using the NLO QCD+EW theoretical prediction. All quoted numbers are rounded to unity. More details on the selection criteria can be found in the text.

Number of expected signal events before the fit, after each selection requirement applied to enter the single-lepton channel resolved $SR^{\geq 6j}_{\geq 4b}$ region. All ${t\bar {t}H}$ signal events are included, regardless of the $H$ or ${t\bar {t}H}$ decay mode. All object corrections are applied, except for the initial number of events which is calculated using the NLO QCD+EW theoretical prediction. All quoted numbers are rounded to unity. More details on the selection criteria can be found in the text.

Number of expected signal events before the fit, after each selection requirement applied to enter the single-lepton channel boosted $SR_{boosted}$ region. All ${t\bar {t}H}$ signal events are included, regardless of the $H$ or ${t\bar {t}H}$ decay mode. All object corrections are applied, except for the initial number of events which is calculated using the NLO QCD+EW theoretical prediction. All quoted numbers are rounded to unity. More details on the selection criteria can be found in the text.