CHI distribution for mass bin > 1200 GeV.

Centrality Ratio.

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

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.

Signal region detector-level distribution for the observable $p_{\mathrm{T}}^{\mathrm{lead}\, \ell}$.

Signal region detector-level distribution for the observable $m_{e\mu}$.

Signal region detector-level distribution for the observable $p_{\mathrm{T}}^{e\mu}$.

Measured fiducial differential cross-section of $WW \rightarrow e^{\pm}\nu\mu^{\mp}\nu$ production for the observable $|y_{e\mu}|$

Relative systematic uncertainties for the fiducial differential cross-section of $WW \rightarrow e^{\pm}\nu\mu^{\mp}\nu$ production for the observable $|y_{e\mu}|$

Measured fiducial differential cross-section of $WW \rightarrow e^{\pm}\nu\mu^{\mp}\nu$ production for the observable $|\Delta \phi(e \mu)|$

Relative systematic uncertainties for the fiducial differential cross-section of $WW \rightarrow e^{\pm}\nu\mu^{\mp}\nu$ production for the observable $|\Delta \phi(e \mu)|$

Measured fiducial differential cross-section of $WW \rightarrow e^{\pm}\nu\mu^{\mp}\nu$ production for the observable $ \cos\theta^{\ast}$

Relative systematic uncertainties for the fiducial differential cross-section of $WW \rightarrow e^{\pm}\nu\mu^{\mp}\nu$ production for the observable $ \cos\theta^{\ast}$

Measured fiducial differential cross-section of $WW \rightarrow e^{\pm}\nu\mu^{\mp}\nu$ production for the observable $p_{\mathrm{T}}^{\mathrm{lead}\, \ell}$

Relative systematic uncertainties for the fiducial differential cross-section of $WW \rightarrow e^{\pm}\nu\mu^{\mp}\nu$ production for the observable $p_{\mathrm{T}}^{\mathrm{lead}\, \ell}$

Measured fiducial differential cross-section of $WW \rightarrow e^{\pm}\nu\mu^{\mp}\nu$ production for the observable $m_{e\mu}$

Relative systematic uncertainties for the fiducial differential cross-section of $WW \rightarrow e^{\pm}\nu\mu^{\mp}\nu$ production for the observable $m_{e\mu}$

Measured fiducial differential cross-section of $WW \rightarrow e^{\pm}\nu\mu^{\mp}\nu$ production for the observable $p_{\mathrm{T}}^{e\mu}$

Relative systematic uncertainties for the fiducial differential cross-section of $WW \rightarrow e^{\pm}\nu\mu^{\mp}\nu$ production for the observable $p_{\mathrm{T}}^{e\mu}$

The statistical correlation coefficients (in percentage) between bins for the measured fiducial differential cross-section of $WW \rightarrow e^{\pm}\nu\mu^{\mp}\nu$ production for the observable $|y_{e\mu}|$

The total correlation coefficients (in percentage) between bins for the measured fiducial differential cross-section of $WW \rightarrow e^{\pm}\nu\mu^{\mp}\nu$ production for the observable $|y_{e\mu}|$

The statistical correlation coefficients (in percentage) between bins for the measured fiducial differential cross-section of $WW \rightarrow e^{\pm}\nu\mu^{\mp}\nu$ production for the observable $|\Delta \phi(e \mu)|$

The total correlation coefficients (in percentage) between bins for the measured fiducial differential cross-section of $WW \rightarrow e^{\pm}\nu\mu^{\mp}\nu$ production for the observable $|\Delta \phi(e \mu)|$

The statistical correlation coefficients (in percentage) between bins for the measured fiducial differential cross-section of $WW \rightarrow e^{\pm}\nu\mu^{\mp}\nu$ production for the observable $ \cos\theta^{\ast}$

The total correlation coefficients (in percentage) between bins for the measured fiducial differential cross-section of $WW \rightarrow e^{\pm}\nu\mu^{\mp}\nu$ production for the observable $ \cos\theta^{\ast}$

The statistical correlation coefficients (in percentage) between bins for the measured fiducial differential cross-section of $WW \rightarrow e^{\pm}\nu\mu^{\mp}\nu$ production for the observable $p_{\mathrm{T}}^{\mathrm{lead}\, \ell}$

The total correlation coefficients (in percentage) between bins for the measured fiducial differential cross-section of $WW \rightarrow e^{\pm}\nu\mu^{\mp}\nu$ production for the observable $p_{\mathrm{T}}^{\mathrm{lead}\, \ell}$

The statistical correlation coefficients (in percentage) between bins for the measured fiducial differential cross-section of $WW \rightarrow e^{\pm}\nu\mu^{\mp}\nu$ production for the observable $m_{e\mu}$

The total correlation coefficients (in percentage) between bins for the measured fiducial differential cross-section of $WW \rightarrow e^{\pm}\nu\mu^{\mp}\nu$ production for the observable $m_{e\mu}$

The statistical correlation coefficients (in percentage) between bins for the measured fiducial differential cross-section of $WW \rightarrow e^{\pm}\nu\mu^{\mp}\nu$ production for the observable $p_{\mathrm{T}}^{e\mu}$

The total correlation coefficients (in percentage) between bins for the measured fiducial differential cross-section of $WW \rightarrow e^{\pm}\nu\mu^{\mp}\nu$ production for the observable $p_{\mathrm{T}}^{e\mu}$

The BDT distribution of the di-Higgs ggF signal for two different values of $\kappa_{\lambda}$ and the main backgrounds in the low mass region ($m^{*}_{b\bar{b}\gamma\gamma} < 350$ GeV). Distributions are normalized to unit area. The dotted lines denote the category boundaries. Events with a BDT score below 0.881 in the low mass region are discarded.

The BDT distribution (with x-axis zoomed in) of the di-Higgs ggF signal for two different values of $\kappa_{\lambda}$ and the main backgrounds in the low mass region ($m^{*}_{b\bar{b}\gamma\gamma} < 350$ GeV). Distributions are normalized to unit area. The dotted lines denote the category boundaries. Events with a BDT score below 0.881 in the low mass region are discarded. The range of BDT scores is from 0.8 to 1.

The BDT distribution (with x-axis zoomed in) of the di-Higgs ggF signal for two different values of $\kappa_{\lambda}$ and the main backgrounds in the low mass region ($m^{*}_{b\bar{b}\gamma\gamma} < 350$ GeV). Distributions are normalized to unit area. The dotted lines denote the category boundaries. Events with a BDT score below 0.881 in the low mass region are discarded. The range of BDT scores is from 0.8 to 1.

The BDT distribution (with x-axis zoomed in) of the di-Higgs ggF signal for two different values of $\kappa_{\lambda}$ and the main backgrounds in the low mass region ($m^{*}_{b\bar{b}\gamma\gamma} < 350$ GeV). Distributions are normalized to unit area. The dotted lines denote the category boundaries. Events with a BDT score below 0.881 in the low mass region are discarded. The range of BDT scores is from 0.8 to 1.

The BDT distribution (with x-axis zoomed in) of the di-Higgs ggF signal for two different values of $\kappa_{\lambda}$ and the main backgrounds in the low mass region ($m^{*}_{b\bar{b}\gamma\gamma} < 350$ GeV). Distributions are normalized to unit area. The dotted lines denote the category boundaries. Events with a BDT score below 0.881 in the low mass region are discarded. The range of BDT scores is from 0.8 to 1.

The BDT distribution of the di-Higgs ggF signal for two different values of $\kappa_{\lambda}$ and the main backgrounds in the high mass region ($m^{*}_{b\bar{b}\gamma\gamma} > 350$ GeV). Distributions are normalized to unit area. The dotted lines denote the category boundaries. Events with a BDT score below 0.857 in the high mass region are discarded.

The BDT distribution of the di-Higgs ggF signal for two different values of $\kappa_{\lambda}$ and the main backgrounds in the high mass region ($m^{*}_{b\bar{b}\gamma\gamma} > 350$ GeV). Distributions are normalized to unit area. The dotted lines denote the category boundaries. Events with a BDT score below 0.857 in the high mass region are discarded.

The BDT distribution of the di-Higgs ggF signal for two different values of $\kappa_{\lambda}$ and the main backgrounds in the high mass region ($m^{*}_{b\bar{b}\gamma\gamma} > 350$ GeV). Distributions are normalized to unit area. The dotted lines denote the category boundaries. Events with a BDT score below 0.857 in the high mass region are discarded.

The BDT distribution of the di-Higgs ggF signal for two different values of $\kappa_{\lambda}$ and the main backgrounds in the high mass region ($m^{*}_{b\bar{b}\gamma\gamma} > 350$ GeV). Distributions are normalized to unit area. The dotted lines denote the category boundaries. Events with a BDT score below 0.857 in the high mass region are discarded.

The BDT distribution (with x-axis zoomed in) of the di-Higgs ggF signal for two different values of $\kappa_{\lambda}$ and the main backgrounds in the high mass region ($m^{*}_{b\bar{b}\gamma\gamma} > 350$ GeV). Distributions are normalized to unit area. The dotted lines denote the category boundaries. Events with a BDT score below 0.857 in the high mass region are discarded. The range of BDT scores is from 0.8 to 1.

The BDT distribution (with x-axis zoomed in) of the di-Higgs ggF signal for two different values of $\kappa_{\lambda}$ and the main backgrounds in the high mass region ($m^{*}_{b\bar{b}\gamma\gamma} > 350$ GeV). Distributions are normalized to unit area. The dotted lines denote the category boundaries. Events with a BDT score below 0.857 in the high mass region are discarded. The range of BDT scores is from 0.8 to 1.

The BDT distribution (with x-axis zoomed in) of the di-Higgs ggF signal for two different values of $\kappa_{\lambda}$ and the main backgrounds in the high mass region ($m^{*}_{b\bar{b}\gamma\gamma} > 350$ GeV). Distributions are normalized to unit area. The dotted lines denote the category boundaries. Events with a BDT score below 0.857 in the high mass region are discarded. The range of BDT scores is from 0.8 to 1.

The BDT distribution (with x-axis zoomed in) of the di-Higgs ggF signal for two different values of $\kappa_{\lambda}$ and the main backgrounds in the high mass region ($m^{*}_{b\bar{b}\gamma\gamma} > 350$ GeV). Distributions are normalized to unit area. The dotted lines denote the category boundaries. Events with a BDT score below 0.857 in the high mass region are discarded. The range of BDT scores is from 0.8 to 1.

The $BDT_{tot}$ score for the benchmark signal $m_{X}$ = 300 GeV and for the main backgrounds. Distributions are normalized to unit area. The dotted line denotes the event selection threshold. Events with a $BDT_{tot}$ score below 0.85 for $m_{X}$ = 300 GeV are discarded.

The $BDT_{tot}$ score for the benchmark signal $m_{X}$ = 300 GeV and for the main backgrounds. Distributions are normalized to unit area. The dotted line denotes the event selection threshold. Events with a $BDT_{tot}$ score below 0.85 for $m_{X}$ = 300 GeV are discarded.

The $BDT_{tot}$ score for the benchmark signal $m_{X}$ = 300 GeV and for the main backgrounds. Distributions are normalized to unit area. The dotted line denotes the event selection threshold. Events with a $BDT_{tot}$ score below 0.85 for $m_{X}$ = 300 GeV are discarded.

The $BDT_{tot}$ score for the benchmark signal $m_{X}$ = 300 GeV and for the main backgrounds. Distributions are normalized to unit area. The dotted line denotes the event selection threshold. Events with a $BDT_{tot}$ score below 0.85 for $m_{X}$ = 300 GeV are discarded.

The $BDT_{tot}$ score for the benchmark signal $m_{X}$ = 500 GeV and for the main backgrounds. Distributions are normalized to unit area. The dotted line denotes the event selection threshold. Events with a $BDT_{tot}$ score below 0.75 for $m_{X}$ = 500 GeV are discarded.

The $BDT_{tot}$ score for the benchmark signal $m_{X}$ = 500 GeV and for the main backgrounds. Distributions are normalized to unit area. The dotted line denotes the event selection threshold. Events with a $BDT_{tot}$ score below 0.75 for $m_{X}$ = 500 GeV are discarded.

The $BDT_{tot}$ score for the benchmark signal $m_{X}$ = 500 GeV and for the main backgrounds. Distributions are normalized to unit area. The dotted line denotes the event selection threshold. Events with a $BDT_{tot}$ score below 0.75 for $m_{X}$ = 500 GeV are discarded.

The $BDT_{tot}$ score for the benchmark signal $m_{X}$ = 500 GeV and for the main backgrounds. Distributions are normalized to unit area. The dotted line denotes the event selection threshold. Events with a $BDT_{tot}$ score below 0.75 for $m_{X}$ = 500 GeV are discarded.

Distributions of $m_{\gamma\gamma}$ in high mass BDT tight category for the nonresonant $HH$ search. The data-derived fractions of nonresonant $\gamma\gamma$, $\gamma$-jet or jet-$\gamma$, and dijet background are applied and the total background is normalized to the data sideband.

Distributions of $m_{\gamma\gamma}$ in high mass BDT tight category for the nonresonant $HH$ search. The data-derived fractions of nonresonant $\gamma\gamma$, $\gamma$-jet or jet-$\gamma$, and dijet background are applied and the total background is normalized to the data sideband.

Distributions of $m_{\gamma\gamma}$ in high mass BDT tight category for the nonresonant $HH$ search. The data-derived fractions of nonresonant $\gamma\gamma$, $\gamma$-jet or jet-$\gamma$, and dijet background are applied and the total background is normalized to the data sideband.

Distributions of $m_{\gamma\gamma}$ in high mass BDT tight category for the nonresonant $HH$ search. The data-derived fractions of nonresonant $\gamma\gamma$, $\gamma$-jet or jet-$\gamma$, and dijet background are applied and the total background is normalized to the data sideband.

Distributions of $m_{\gamma\gamma}$ in high mass BDT loose category for the nonresonant $HH$ search. The data-derived fractions of nonresonant $\gamma\gamma$, $\gamma$-jet or jet-$\gamma$, and dijet background are applied and the total background is normalized to the data sideband.

Distributions of $m_{\gamma\gamma}$ in high mass BDT loose category for the nonresonant $HH$ search. The data-derived fractions of nonresonant $\gamma\gamma$, $\gamma$-jet or jet-$\gamma$, and dijet background are applied and the total background is normalized to the data sideband.

Distributions of $m_{\gamma\gamma}$ in high mass BDT loose category for the nonresonant $HH$ search. The data-derived fractions of nonresonant $\gamma\gamma$, $\gamma$-jet or jet-$\gamma$, and dijet background are applied and the total background is normalized to the data sideband.

Distributions of $m_{\gamma\gamma}$ in high mass BDT loose category for the nonresonant $HH$ search. The data-derived fractions of nonresonant $\gamma\gamma$, $\gamma$-jet or jet-$\gamma$, and dijet background are applied and the total background is normalized to the data sideband.

Distributions of $m_{\gamma\gamma}$ in low mass BDT tight category for the nonresonant $HH$ search. The data-derived fractions of nonresonant $\gamma\gamma$, $\gamma$-jet or jet-$\gamma$, and dijet background are applied and the total background is normalized to the data sideband.

Distributions of $m_{\gamma\gamma}$ in low mass BDT tight category for the nonresonant $HH$ search. The data-derived fractions of nonresonant $\gamma\gamma$, $\gamma$-jet or jet-$\gamma$, and dijet background are applied and the total background is normalized to the data sideband.

Distributions of $m_{\gamma\gamma}$ in low mass BDT tight category for the nonresonant $HH$ search. The data-derived fractions of nonresonant $\gamma\gamma$, $\gamma$-jet or jet-$\gamma$, and dijet background are applied and the total background is normalized to the data sideband.

Distributions of $m_{\gamma\gamma}$ in low mass BDT tight category for the nonresonant $HH$ search. The data-derived fractions of nonresonant $\gamma\gamma$, $\gamma$-jet or jet-$\gamma$, and dijet background are applied and the total background is normalized to the data sideband.

Distributions of $m_{\gamma\gamma}$ in low mass BDT loose category for the nonresonant $HH$ search. The data-derived fractions of nonresonant $\gamma\gamma$, $\gamma$-jet or jet-$\gamma$, and dijet background are applied and the total background is normalized to the data sideband.

Distributions of $m_{\gamma\gamma}$ in low mass BDT loose category for the nonresonant $HH$ search. The data-derived fractions of nonresonant $\gamma\gamma$, $\gamma$-jet or jet-$\gamma$, and dijet background are applied and the total background is normalized to the data sideband.

Distributions of $m_{\gamma\gamma}$ in low mass BDT loose category for the nonresonant $HH$ search. The data-derived fractions of nonresonant $\gamma\gamma$, $\gamma$-jet or jet-$\gamma$, and dijet background are applied and the total background is normalized to the data sideband.

Distributions of $m_{\gamma\gamma}$ in low mass BDT loose category for the nonresonant $HH$ search. The data-derived fractions of nonresonant $\gamma\gamma$, $\gamma$-jet or jet-$\gamma$, and dijet background are applied and the total background is normalized to the data sideband.

Distributions of $m_{\gamma\gamma}$ for the selections used for the resonance mass point $m_{X}$ = 300 GeV for the resonant search. The data-derived fractions of nonresonant $\gamma\gamma$, $\gamma$-jet or jet-$\gamma$, and dijet background are applied and the total background is normalized to the data sideband. The scalar resonance signal is scaled to a total production cross section $\sigma(pp \rightarrow X \rightarrow HH)$ = 370 fb.

Distributions of $m_{\gamma\gamma}$ for the selections used for the resonance mass point $m_{X}$ = 300 GeV for the resonant search. The data-derived fractions of nonresonant $\gamma\gamma$, $\gamma$-jet or jet-$\gamma$, and dijet background are applied and the total background is normalized to the data sideband. The scalar resonance signal is scaled to a total production cross section $\sigma(pp \rightarrow X \rightarrow HH)$ = 370 fb.

Distributions of $m_{\gamma\gamma}$ for the selections used for the resonance mass point $m_{X}$ = 300 GeV for the resonant search. The data-derived fractions of nonresonant $\gamma\gamma$, $\gamma$-jet or jet-$\gamma$, and dijet background are applied and the total background is normalized to the data sideband. The scalar resonance signal is scaled to a total production cross section $\sigma(pp \rightarrow X \rightarrow HH)$ = 370 fb.

Distributions of $m_{\gamma\gamma}$ for the selections used for the resonance mass point $m_{X}$ = 300 GeV for the resonant search. The data-derived fractions of nonresonant $\gamma\gamma$, $\gamma$-jet or jet-$\gamma$, and dijet background are applied and the total background is normalized to the data sideband. The scalar resonance signal is scaled to a total production cross section $\sigma(pp \rightarrow X \rightarrow HH)$ = 370 fb.

Distributions of $m_{\gamma\gamma}$ for the selections used for the resonance mass point $m_{X}$ = 500 GeV for the resonant search. The data-derived fractions of nonresonant $\gamma\gamma$, $\gamma$-jet or jet-$\gamma$, and dijet background are applied and the total background is normalized to the data sideband. The scalar resonance signal is scaled to a total production cross section $\sigma(pp \rightarrow X \rightarrow HH)$ = 67 fb.

Distributions of $m_{\gamma\gamma}$ for the selections used for the resonance mass point $m_{X}$ = 500 GeV for the resonant search. The data-derived fractions of nonresonant $\gamma\gamma$, $\gamma$-jet or jet-$\gamma$, and dijet background are applied and the total background is normalized to the data sideband. The scalar resonance signal is scaled to a total production cross section $\sigma(pp \rightarrow X \rightarrow HH)$ = 67 fb.

Distributions of $m_{\gamma\gamma}$ for the selections used for the resonance mass point $m_{X}$ = 500 GeV for the resonant search. The data-derived fractions of nonresonant $\gamma\gamma$, $\gamma$-jet or jet-$\gamma$, and dijet background are applied and the total background is normalized to the data sideband. The scalar resonance signal is scaled to a total production cross section $\sigma(pp \rightarrow X \rightarrow HH)$ = 67 fb.

Distributions of $m_{\gamma\gamma}$ for the selections used for the resonance mass point $m_{X}$ = 500 GeV for the resonant search. The data-derived fractions of nonresonant $\gamma\gamma$, $\gamma$-jet or jet-$\gamma$, and dijet background are applied and the total background is normalized to the data sideband. The scalar resonance signal is scaled to a total production cross section $\sigma(pp \rightarrow X \rightarrow HH)$ = 67 fb.

The number of data events observed in the 120 GeV < $m_{\gamma\gamma}$ < 130 GeV window, the number of $HH$ signal events expected for $\kappa_{\lambda}$ = 1 and for $\kappa_{\lambda}$ = 10, and events expected for single Higgs boson production (estimated using MC simulation), as well as for continuum background. For the single Higgs boson, Rest includes VBF, $WH$, $tHqb$, and $tHW$. The values are obtained from a fit of the Asimov data set generated under the SM signal-plus-background hypothesis, $\kappa_{\lambda}$ = 1. The continuum background component of the Asimov data set is obtained from the fit of the data sideband. The uncertainties in $HH$ signals and single Higgs boson background include the systematic uncertainties discussed in Section 6. The uncertainty in the continuum background is given by the sum in quadrature of the statistical uncertainty from the fit to the data and the spurious-signal uncertainty.

The number of data events observed in the 120 GeV < $m_{\gamma\gamma}$ < 130 GeV window, the number of $HH$ signal events expected for $\kappa_{\lambda}$ = 1 and for $\kappa_{\lambda}$ = 10, and events expected for single Higgs boson production (estimated using MC simulation), as well as for continuum background. For the single Higgs boson, Rest includes VBF, $WH$, $tHqb$, and $tHW$. The values are obtained from a fit of the Asimov data set generated under the SM signal-plus-background hypothesis, $\kappa_{\lambda}$ = 1. The continuum background component of the Asimov data set is obtained from the fit of the data sideband. The uncertainties in $HH$ signals and single Higgs boson background include the systematic uncertainties discussed in Section 6. The uncertainty in the continuum background is given by the sum in quadrature of the statistical uncertainty from the fit to the data and the spurious-signal uncertainty.

The number of data events observed in the 120 GeV < $m_{\gamma\gamma}$ < 130 GeV window, the number of $HH$ signal events expected for $\kappa_{\lambda}$ = 1 and for $\kappa_{\lambda}$ = 10, and events expected for single Higgs boson production (estimated using MC simulation), as well as for continuum background. For the single Higgs boson, Rest includes VBF, $WH$, $tHqb$, and $tHW$. The values are obtained from a fit of the Asimov data set generated under the SM signal-plus-background hypothesis, $\kappa_{\lambda}$ = 1. The continuum background component of the Asimov data set is obtained from the fit of the data sideband. The uncertainties in $HH$ signals and single Higgs boson background include the systematic uncertainties discussed in Section 6. The uncertainty in the continuum background is given by the sum in quadrature of the statistical uncertainty from the fit to the data and the spurious-signal uncertainty.

The number of data events observed in the 120 GeV < $m_{\gamma\gamma}$ < 130 GeV window, the number of $HH$ signal events expected for $\kappa_{\lambda}$ = 1 and for $\kappa_{\lambda}$ = 10, and events expected for single Higgs boson production (estimated using MC simulation), as well as for continuum background. For the single Higgs boson, Rest includes VBF, $WH$, $tHqb$, and $tHW$. The values are obtained from a fit of the Asimov data set generated under the SM signal-plus-background hypothesis, $\kappa_{\lambda}$ = 1. The continuum background component of the Asimov data set is obtained from the fit of the data sideband. The uncertainties in $HH$ signals and single Higgs boson background include the systematic uncertainties discussed in Section 6. The uncertainty in the continuum background is given by the sum in quadrature of the statistical uncertainty from the fit to the data and the spurious-signal uncertainty.

Observed and expected limits at 95% CL on the cross section of nonresonant Higgs boson pair production as a function of the Higgs boson self-coupling modifier $\kappa_{\lambda}= \lambda_{HHH}/\lambda^{\textrm{SM}}_{HHH}$. The expected constraints on $\kappa_{\lambda}$ are obtained with a background hypothesis excluding $pp \rightarrow HH$ production. The $\pm 1\sigma$ and $\pm 2\sigma$ variations about the expected limit due to statistical and systematic uncertainties are also shown. The theory prediction curve represents the scenario where all parameters and couplings are set to their SM values except for $\kappa_{\lambda}$. The uncertainty band of the theory prediction curve shows the cross-section uncertainty.

Observed and expected limits at 95% CL on the cross section of nonresonant Higgs boson pair production as a function of the Higgs boson self-coupling modifier $\kappa_{\lambda}= \lambda_{HHH}/\lambda^{\textrm{SM}}_{HHH}$. The expected constraints on $\kappa_{\lambda}$ are obtained with a background hypothesis excluding $pp \rightarrow HH$ production. The $\pm 1\sigma$ and $\pm 2\sigma$ variations about the expected limit due to statistical and systematic uncertainties are also shown. The theory prediction curve represents the scenario where all parameters and couplings are set to their SM values except for $\kappa_{\lambda}$. The uncertainty band of the theory prediction curve shows the cross-section uncertainty.

Observed and expected limits at 95% CL on the cross section of nonresonant Higgs boson pair production as a function of the Higgs boson self-coupling modifier $\kappa_{\lambda}= \lambda_{HHH}/\lambda^{\textrm{SM}}_{HHH}$. The expected constraints on $\kappa_{\lambda}$ are obtained with a background hypothesis excluding $pp \rightarrow HH$ production. The $\pm 1\sigma$ and $\pm 2\sigma$ variations about the expected limit due to statistical and systematic uncertainties are also shown. The theory prediction curve represents the scenario where all parameters and couplings are set to their SM values except for $\kappa_{\lambda}$. The uncertainty band of the theory prediction curve shows the cross-section uncertainty.

Observed and expected limits at 95% CL on the cross section of nonresonant Higgs boson pair production as a function of the Higgs boson self-coupling modifier $\kappa_{\lambda}= \lambda_{HHH}/\lambda^{\textrm{SM}}_{HHH}$. The expected constraints on $\kappa_{\lambda}$ are obtained with a background hypothesis excluding $pp \rightarrow HH$ production. The $\pm 1\sigma$ and $\pm 2\sigma$ variations about the expected limit due to statistical and systematic uncertainties are also shown. The theory prediction curve represents the scenario where all parameters and couplings are set to their SM values except for $\kappa_{\lambda}$. The uncertainty band of the theory prediction curve shows the cross-section uncertainty.

Values of the negative log-profile-likelihood ratio ($-2ln\Lambda$) as a function of $\kappa_{\lambda}$ evaluated for the combination of all the categories of the nonresonant search. The coupling of the Higgs boson to fermions and gauge bosons is set to SM values in the profile likelihood calculation. The Asimov data set is generated under the SM signal-plus-background hypothesis, $\kappa_{\lambda}$= 1. All systematic uncertainties, including the theoretical uncertainties in the di-Higgs boson production cross section, are included. The intersections of the solid curves and the horizontal dashed lines indicate the 1$\sigma$ and 2$\sigma$ confidence-level intervals.

Values of the negative log-profile-likelihood ratio ($-2ln\Lambda$) as a function of $\kappa_{\lambda}$ evaluated for the combination of all the categories of the nonresonant search. The coupling of the Higgs boson to fermions and gauge bosons is set to SM values in the profile likelihood calculation. The Asimov data set is generated under the SM signal-plus-background hypothesis, $\kappa_{\lambda}$= 1. All systematic uncertainties, including the theoretical uncertainties in the di-Higgs boson production cross section, are included. The intersections of the solid curves and the horizontal dashed lines indicate the 1$\sigma$ and 2$\sigma$ confidence-level intervals.

Values of the negative log-profile-likelihood ratio ($-2ln\Lambda$) as a function of $\kappa_{\lambda}$ evaluated for the combination of all the categories of the nonresonant search. The coupling of the Higgs boson to fermions and gauge bosons is set to SM values in the profile likelihood calculation. The Asimov data set is generated under the SM signal-plus-background hypothesis, $\kappa_{\lambda}$= 1. All systematic uncertainties, including the theoretical uncertainties in the di-Higgs boson production cross section, are included. The intersections of the solid curves and the horizontal dashed lines indicate the 1$\sigma$ and 2$\sigma$ confidence-level intervals.

Values of the negative log-profile-likelihood ratio ($-2ln\Lambda$) as a function of $\kappa_{\lambda}$ evaluated for the combination of all the categories of the nonresonant search. The coupling of the Higgs boson to fermions and gauge bosons is set to SM values in the profile likelihood calculation. The Asimov data set is generated under the SM signal-plus-background hypothesis, $\kappa_{\lambda}$= 1. All systematic uncertainties, including the theoretical uncertainties in the di-Higgs boson production cross section, are included. The intersections of the solid curves and the horizontal dashed lines indicate the 1$\sigma$ and 2$\sigma$ confidence-level intervals.

The number of events observed in the 120 < $m_{\gamma\gamma}$ < 130 GeV window in data, the number of events expected for scalar resonance signals of masses $m_{X}$ = 300 GeV and $m_{X}$ = 500 GeV assuming a total production cross section $\sigma(pp \rightarrow X \rightarrow HH)$ equal to the observed exclusion limits of Figure 15, and events expected for SM $HH$ and single Higgs boson production (estimated using MC simulation), as well as for continuum background. The values are obtained from a fit of the Asimov data set generated under the signal-plus-background hypothesis. The continuum background component of the Asimov data set is obtained from the fit of the data sideband. The uncertainties in the resonant signals and the SM $HH$ and single-Higgs-boson backgrounds include the systematic uncertainties discussed in Section 6. The uncertainty in the continuum background is given by the sum in quadrature of the statistical uncertainty from the fit to the data and the spurious-signal uncertainty.

The number of events observed in the 120 < $m_{\gamma\gamma}$ < 130 GeV window in data, the number of events expected for scalar resonance signals of masses $m_{X}$ = 300 GeV and $m_{X}$ = 500 GeV assuming a total production cross section $\sigma(pp \rightarrow X \rightarrow HH)$ equal to the observed exclusion limits of Figure 15, and events expected for SM $HH$ and single Higgs boson production (estimated using MC simulation), as well as for continuum background. The values are obtained from a fit of the Asimov data set generated under the signal-plus-background hypothesis. The continuum background component of the Asimov data set is obtained from the fit of the data sideband. The uncertainties in the resonant signals and the SM $HH$ and single-Higgs-boson backgrounds include the systematic uncertainties discussed in Section 6. The uncertainty in the continuum background is given by the sum in quadrature of the statistical uncertainty from the fit to the data and the spurious-signal uncertainty.

The number of events observed in the 120 < $m_{\gamma\gamma}$ < 130 GeV window in data, the number of events expected for scalar resonance signals of masses $m_{X}$ = 300 GeV and $m_{X}$ = 500 GeV assuming a total production cross section $\sigma(pp \rightarrow X \rightarrow HH)$ equal to the observed exclusion limits of Figure 15, and events expected for SM $HH$ and single Higgs boson production (estimated using MC simulation), as well as for continuum background. The values are obtained from a fit of the Asimov data set generated under the signal-plus-background hypothesis. The continuum background component of the Asimov data set is obtained from the fit of the data sideband. The uncertainties in the resonant signals and the SM $HH$ and single-Higgs-boson backgrounds include the systematic uncertainties discussed in Section 6. The uncertainty in the continuum background is given by the sum in quadrature of the statistical uncertainty from the fit to the data and the spurious-signal uncertainty.

The number of events observed in the 120 < $m_{\gamma\gamma}$ < 130 GeV window in data, the number of events expected for scalar resonance signals of masses $m_{X}$ = 300 GeV and $m_{X}$ = 500 GeV assuming a total production cross section $\sigma(pp \rightarrow X \rightarrow HH)$ equal to the observed exclusion limits of Figure 15, and events expected for SM $HH$ and single Higgs boson production (estimated using MC simulation), as well as for continuum background. The values are obtained from a fit of the Asimov data set generated under the signal-plus-background hypothesis. The continuum background component of the Asimov data set is obtained from the fit of the data sideband. The uncertainties in the resonant signals and the SM $HH$ and single-Higgs-boson backgrounds include the systematic uncertainties discussed in Section 6. The uncertainty in the continuum background is given by the sum in quadrature of the statistical uncertainty from the fit to the data and the spurious-signal uncertainty.

Observed and expected limits at 95% CL on the production cross section of a narrow-width scalar resonance $X$ as a function of the mass $m_{X}$ of the hypothetical scalar particle. The black solid line represents the observed upper limits. The dashed line represents the expected upper limits. The $\pm 1\sigma$ and $\pm 2\sigma$ variations about the expected limit due to statistical and systematic uncertainties are also shown.

Observed and expected limits at 95% CL on the production cross section of a narrow-width scalar resonance $X$ as a function of the mass $m_{X}$ of the hypothetical scalar particle. The black solid line represents the observed upper limits. The dashed line represents the expected upper limits. The $\pm 1\sigma$ and $\pm 2\sigma$ variations about the expected limit due to statistical and systematic uncertainties are also shown.

Observed and expected limits at 95% CL on the production cross section of a narrow-width scalar resonance $X$ as a function of the mass $m_{X}$ of the hypothetical scalar particle. The black solid line represents the observed upper limits. The dashed line represents the expected upper limits. The $\pm 1\sigma$ and $\pm 2\sigma$ variations about the expected limit due to statistical and systematic uncertainties are also shown.

Observed and expected limits at 95% CL on the production cross section of a narrow-width scalar resonance $X$ as a function of the mass $m_{X}$ of the hypothetical scalar particle. The black solid line represents the observed upper limits. The dashed line represents the expected upper limits. The $\pm 1\sigma$ and $\pm 2\sigma$ variations about the expected limit due to statistical and systematic uncertainties are also shown.

Breakdown of the dominant systematic uncertainties. The impact of the uncertainties is defined according to the statistical analysis described in Section 7. It corresponds to the relative variation of the expected upper limit on the cross section when re-evaluating the profile likelihood ratio after fixing the nuisance parameter in question to its best-fit value, while all remaining nuisance parameters remain free to float. The impact is shown in %. Only systematic uncertainties with an impact of at least 0.2% are shown. Uncertainties of the "Norm. + Shape" type affect both the normalization and the parameters of the functional form. The rest of the uncertainties affect only the yields.

Breakdown of the dominant systematic uncertainties. The impact of the uncertainties is defined according to the statistical analysis described in Section 7. It corresponds to the relative variation of the expected upper limit on the cross section when re-evaluating the profile likelihood ratio after fixing the nuisance parameter in question to its best-fit value, while all remaining nuisance parameters remain free to float. The impact is shown in %. Only systematic uncertainties with an impact of at least 0.2% are shown. Uncertainties of the "Norm. + Shape" type affect both the normalization and the parameters of the functional form. The rest of the uncertainties affect only the yields.

Breakdown of the dominant systematic uncertainties. The impact of the uncertainties is defined according to the statistical analysis described in Section 7. It corresponds to the relative variation of the expected upper limit on the cross section when re-evaluating the profile likelihood ratio after fixing the nuisance parameter in question to its best-fit value, while all remaining nuisance parameters remain free to float. The impact is shown in %. Only systematic uncertainties with an impact of at least 0.2% are shown. Uncertainties of the "Norm. + Shape" type affect both the normalization and the parameters of the functional form. The rest of the uncertainties affect only the yields.

Breakdown of the dominant systematic uncertainties. The impact of the uncertainties is defined according to the statistical analysis described in Section 7. It corresponds to the relative variation of the expected upper limit on the cross section when re-evaluating the profile likelihood ratio after fixing the nuisance parameter in question to its best-fit value, while all remaining nuisance parameters remain free to float. The impact is shown in %. Only systematic uncertainties with an impact of at least 0.2% are shown. Uncertainties of the "Norm. + Shape" type affect both the normalization and the parameters of the functional form. The rest of the uncertainties affect only the yields.

Cutflow for nonresonant di-Higgs ggF signal sample, yields are normalized to 139 $fb^{-1}$.

Cutflow for nonresonant di-Higgs ggF signal sample, yields are normalized to 139 $fb^{-1}$.

Cutflow for nonresonant di-Higgs ggF signal sample, yields are normalized to 139 $fb^{-1}$.

Cutflow for nonresonant di-Higgs ggF signal sample, yields are normalized to 139 $fb^{-1}$.

Cutflow for resonant signal sample, with $m_{X}$ = 300 GeV, yields are normalized to 139 $fb^{-1}$.

Cutflow for resonant signal sample, with $m_{X}$ = 300 GeV, yields are normalized to 139 $fb^{-1}$.

Cutflow for resonant signal sample, with $m_{X}$ = 300 GeV, yields are normalized to 139 $fb^{-1}$.

Cutflow for resonant signal sample, with $m_{X}$ = 300 GeV, yields are normalized to 139 $fb^{-1}$.

Cutflow for resonant signal sample, with $m_{X}$ = 500 GeV, yields are normalized to 139 $fb^{-1}$.

Cutflow for resonant signal sample, with $m_{X}$ = 500 GeV, yields are normalized to 139 $fb^{-1}$.

Cutflow for resonant signal sample, with $m_{X}$ = 500 GeV, yields are normalized to 139 $fb^{-1}$.

Cutflow for resonant signal sample, with $m_{X}$ = 500 GeV, yields are normalized to 139 $fb^{-1}$.

Comparison of $m_{b\bar{b}}$ distributions when applying the specific b-jet energy calibration and the nominal jet energy calibration. The distributions are fitted using a Bukin function, and the values of the peak position, resolution and the relative improvement are reported in the legend.

Comparison of $m_{b\bar{b}}$ distributions when applying the specific b-jet energy calibration and the nominal jet energy calibration. The distributions are fitted using a Bukin function, and the values of the peak position, resolution and the relative improvement are reported in the legend.

Comparison of $m_{b\bar{b}}$ distributions when applying the specific b-jet energy calibration and the nominal jet energy calibration. The distributions are fitted using a Bukin function, and the values of the peak position, resolution and the relative improvement are reported in the legend.

Comparison of $m_{b\bar{b}}$ distributions when applying the specific b-jet energy calibration and the nominal jet energy calibration. The distributions are fitted using a Bukin function, and the values of the peak position, resolution and the relative improvement are reported in the legend.

Fit results of $m_{b\bar{b}}$ distributions when applying the specific b-jet energy calibration and the nominal jet energy calibration. The distributions are fitted using a Bukin function, and the values of the peak position, resolution and the relative improvement are reported in the legend.

Fit results of $m_{b\bar{b}}$ distributions when applying the specific b-jet energy calibration and the nominal jet energy calibration. The distributions are fitted using a Bukin function, and the values of the peak position, resolution and the relative improvement are reported in the legend.

Fit results of $m_{b\bar{b}}$ distributions when applying the specific b-jet energy calibration and the nominal jet energy calibration. The distributions are fitted using a Bukin function, and the values of the peak position, resolution and the relative improvement are reported in the legend.

Fit results of $m_{b\bar{b}}$ distributions when applying the specific b-jet energy calibration and the nominal jet energy calibration. The distributions are fitted using a Bukin function, and the values of the peak position, resolution and the relative improvement are reported in the legend.

The relative amount (purity) of expected events from SM $HH$ and single Higgs boson production processes for each of the four categories of the nonresonant search. The Higgs boson pair production with $k_{\lambda} = 1$ is considered as signal in (a), while the case with $k_{\lambda} = 10$ is considered as signal in (b).

The relative amount (purity) of expected events from SM $HH$ and single Higgs boson production processes for each of the four categories of the nonresonant search. The Higgs boson pair production with $k_{\lambda} = 1$ is considered as signal in (a), while the case with $k_{\lambda} = 10$ is considered as signal in (b).

The relative amount (purity) of expected events from SM $HH$ and single Higgs boson production processes for each of the four categories of the nonresonant search. The Higgs boson pair production with $k_{\lambda} = 1$ is considered as signal in (a), while the case with $k_{\lambda} = 10$ is considered as signal in (b).

The relative amount (purity) of expected events from SM $HH$ and single Higgs boson production processes for each of the four categories of the nonresonant search. The Higgs boson pair production with $k_{\lambda} = 1$ is considered as signal in (a), while the case with $k_{\lambda} = 10$ is considered as signal in (b).

The expected significance in each of the four categories of the nonresonant search. The Higgs boson pair production with $k_{\lambda} = 1$ is considered as signal in (a), while the case with $k_{\lambda} = 10$ is considered as signal in (b). The single Higgs boson processes and the di-photon continuum spectrum are considered as background.

The expected significance in each of the four categories of the nonresonant search. The Higgs boson pair production with $k_{\lambda} = 1$ is considered as signal in (a), while the case with $k_{\lambda} = 10$ is considered as signal in (b). The single Higgs boson processes and the di-photon continuum spectrum are considered as background.

The expected significance in each of the four categories of the nonresonant search. The Higgs boson pair production with $k_{\lambda} = 1$ is considered as signal in (a), while the case with $k_{\lambda} = 10$ is considered as signal in (b). The single Higgs boson processes and the di-photon continuum spectrum are considered as background.

The expected significance in each of the four categories of the nonresonant search. The Higgs boson pair production with $k_{\lambda} = 1$ is considered as signal in (a), while the case with $k_{\lambda} = 10$ is considered as signal in (b). The single Higgs boson processes and the di-photon continuum spectrum are considered as background.

Spurious signal result for the exponential pdf for the various ggF nonresonant di-Higgs categories. In each category, the spurious signal value ($n_{sp}$) and its ratio to the expected statistical error ($Z_{spur}$) from data are shown.

Spurious signal result for the exponential pdf for the various ggF nonresonant di-Higgs categories. In each category, the spurious signal value ($n_{sp}$) and its ratio to the expected statistical error ($Z_{spur}$) from data are shown.

Spurious signal result for the exponential pdf for the various ggF nonresonant di-Higgs categories. In each category, the spurious signal value ($n_{sp}$) and its ratio to the expected statistical error ($Z_{spur}$) from data are shown.

Spurious signal result for the exponential pdf for the various ggF nonresonant di-Higgs categories. In each category, the spurious signal value ($n_{sp}$) and its ratio to the expected statistical error ($Z_{spur}$) from data are shown.

Spurious signal result for the exponential pdf as function of the resonant di-Higgs signal mass. The spurious signal value and its ratio to the expected statistical error from data are shown.

Spurious signal result for the exponential pdf as function of the resonant di-Higgs signal mass. The spurious signal value and its ratio to the expected statistical error from data are shown.

Spurious signal result for the exponential pdf as function of the resonant di-Higgs signal mass. The spurious signal value and its ratio to the expected statistical error from data are shown.

Spurious signal result for the exponential pdf as function of the resonant di-Higgs signal mass. The spurious signal value and its ratio to the expected statistical error from data are shown.

Distributions of the signal efficiency as a function of $\kappa_{\lambda}$, for the di-Higgs boson ggF nonresonant production mode. The range of $\kappa_{\lambda}$ in the table is from -10 to 10.

Distributions of the signal efficiency as a function of $\kappa_{\lambda}$, for the di-Higgs boson ggF nonresonant production mode. The range of $\kappa_{\lambda}$ in the table is from -10 to 10.

Distributions of the signal efficiency as a function of $\kappa_{\lambda}$, for the di-Higgs boson ggF nonresonant production mode. The range of $\kappa_{\lambda}$ in the table is from -10 to 10.

Distributions of the signal efficiency as a function of $\kappa_{\lambda}$, for the di-Higgs boson ggF nonresonant production mode. The range of $\kappa_{\lambda}$ in the table is from -10 to 10.

Distributions of the signal efficiency as a function of $\kappa_{\lambda}$, for the di-Higgs boson VBF nonresonant production mode. The range of $\kappa_{\lambda}$ in the table is from -10 to 10.

Distributions of the signal efficiency as a function of $\kappa_{\lambda}$, for the di-Higgs boson VBF nonresonant production mode. The range of $\kappa_{\lambda}$ in the table is from -10 to 10.

Distributions of the signal efficiency as a function of $\kappa_{\lambda}$, for the di-Higgs boson VBF nonresonant production mode. The range of $\kappa_{\lambda}$ in the table is from -10 to 10.

Distributions of the signal efficiency as a function of $\kappa_{\lambda}$, for the di-Higgs boson VBF nonresonant production mode. The range of $\kappa_{\lambda}$ in the table is from -10 to 10.

Values of the negative log-profile-likelihood ($-2ln\Lambda$) as a function of $\kappa_{\lambda}$ evaluated for the combination of all the categories of the nonresonant search. The coupling of the Higgs boson to fermions and gauge bosons is set to SM values in the profile likelihood calculation. The Asimov data set is generated under the SM signal-plus-background hypothesis, $\kappa_{\lambda}$= 1. All systematic uncertainties, including the theoretical uncertainties on the di-Higgs boson production cross section, are included. The intersections of the solid curves and the horizontal dashed lines indicate the 1$\sigma$ and 2$\sigma$ confidence level intervals. The best fit value corresponds to $\kappa_{\lambda}$ = $2.8^{+2.0}_{-2.2}(^{+3.8}_{-4.3})$ for the 1$\sigma$(2$\sigma$) confidence interval. The expected value corresponds to $\kappa_{\lambda}$ = $1.0^{+5.5}_{-2.4}(^{+7.3}_{-4.2})$ for the 1$\sigma$(2$\sigma$) confidence interval. The dashed curves represent values of the negative log-profile-likelihood where the Higgs boson branching fractions and the cross section of the production modes are varied as a function of $\kappa_{\lambda}$. In this case,the best fit value corresponds to $\kappa_{\lambda}$ = $2.7^{+2.0}_{-2.2}(^{+3.8}_{-4.3})$ and the expected value corresponds to $\kappa_{\lambda}$ = $1.0^{+5.4}_{-2.5}(^{+7.3}_{-4.3})$ for the 1$\sigma$(2$\sigma$) confidence interval.

Values of the negative log-profile-likelihood ($-2ln\Lambda$) as a function of $\kappa_{\lambda}$ evaluated for the combination of all the categories of the nonresonant search. The coupling of the Higgs boson to fermions and gauge bosons is set to SM values in the profile likelihood calculation. The Asimov data set is generated under the SM signal-plus-background hypothesis, $\kappa_{\lambda}$= 1. All systematic uncertainties, including the theoretical uncertainties on the di-Higgs boson production cross section, are included. The intersections of the solid curves and the horizontal dashed lines indicate the 1$\sigma$ and 2$\sigma$ confidence level intervals. The best fit value corresponds to $\kappa_{\lambda}$ = $2.8^{+2.0}_{-2.2}(^{+3.8}_{-4.3})$ for the 1$\sigma$(2$\sigma$) confidence interval. The expected value corresponds to $\kappa_{\lambda}$ = $1.0^{+5.5}_{-2.4}(^{+7.3}_{-4.2})$ for the 1$\sigma$(2$\sigma$) confidence interval. The dashed curves represent values of the negative log-profile-likelihood where the Higgs boson branching fractions and the cross section of the production modes are varied as a function of $\kappa_{\lambda}$. In this case,the best fit value corresponds to $\kappa_{\lambda}$ = $2.7^{+2.0}_{-2.2}(^{+3.8}_{-4.3})$ and the expected value corresponds to $\kappa_{\lambda}$ = $1.0^{+5.4}_{-2.5}(^{+7.3}_{-4.3})$ for the 1$\sigma$(2$\sigma$) confidence interval.

Values of the negative log-profile-likelihood ($-2ln\Lambda$) as a function of $\kappa_{\lambda}$ evaluated for the combination of all the categories of the nonresonant search. The coupling of the Higgs boson to fermions and gauge bosons is set to SM values in the profile likelihood calculation. The Asimov data set is generated under the SM signal-plus-background hypothesis, $\kappa_{\lambda}$= 1. All systematic uncertainties, including the theoretical uncertainties on the di-Higgs boson production cross section, are included. The intersections of the solid curves and the horizontal dashed lines indicate the 1$\sigma$ and 2$\sigma$ confidence level intervals. The best fit value corresponds to $\kappa_{\lambda}$ = $2.8^{+2.0}_{-2.2}(^{+3.8}_{-4.3})$ for the 1$\sigma$(2$\sigma$) confidence interval. The expected value corresponds to $\kappa_{\lambda}$ = $1.0^{+5.5}_{-2.4}(^{+7.3}_{-4.2})$ for the 1$\sigma$(2$\sigma$) confidence interval. The dashed curves represent values of the negative log-profile-likelihood where the Higgs boson branching fractions and the cross section of the production modes are varied as a function of $\kappa_{\lambda}$. In this case,the best fit value corresponds to $\kappa_{\lambda}$ = $2.7^{+2.0}_{-2.2}(^{+3.8}_{-4.3})$ and the expected value corresponds to $\kappa_{\lambda}$ = $1.0^{+5.4}_{-2.5}(^{+7.3}_{-4.3})$ for the 1$\sigma$(2$\sigma$) confidence interval.

Values of the negative log-profile-likelihood ($-2ln\Lambda$) as a function of $\kappa_{\lambda}$ evaluated for the combination of all the categories of the nonresonant search. The coupling of the Higgs boson to fermions and gauge bosons is set to SM values in the profile likelihood calculation. The Asimov data set is generated under the SM signal-plus-background hypothesis, $\kappa_{\lambda}$= 1. All systematic uncertainties, including the theoretical uncertainties on the di-Higgs boson production cross section, are included. The intersections of the solid curves and the horizontal dashed lines indicate the 1$\sigma$ and 2$\sigma$ confidence level intervals. The best fit value corresponds to $\kappa_{\lambda}$ = $2.8^{+2.0}_{-2.2}(^{+3.8}_{-4.3})$ for the 1$\sigma$(2$\sigma$) confidence interval. The expected value corresponds to $\kappa_{\lambda}$ = $1.0^{+5.5}_{-2.4}(^{+7.3}_{-4.2})$ for the 1$\sigma$(2$\sigma$) confidence interval. The dashed curves represent values of the negative log-profile-likelihood where the Higgs boson branching fractions and the cross section of the production modes are varied as a function of $\kappa_{\lambda}$. In this case,the best fit value corresponds to $\kappa_{\lambda}$ = $2.7^{+2.0}_{-2.2}(^{+3.8}_{-4.3})$ and the expected value corresponds to $\kappa_{\lambda}$ = $1.0^{+5.4}_{-2.5}(^{+7.3}_{-4.3})$ for the 1$\sigma$(2$\sigma$) confidence interval.

Minimum BDT value of the events passing the selection criteria of the resonant search. The combined BDT score is formed using as coefficients $C_{1}$ = 0.65 and $C_{2}$ = 1 − $C_{1}$. The selection efficiency for the resonant $X \rightarrow HH$ signal is also shown.

Minimum BDT value of the events passing the selection criteria of the resonant search. The combined BDT score is formed using as coefficients $C_{1}$ = 0.65 and $C_{2}$ = 1 − $C_{1}$. The selection efficiency for the resonant $X \rightarrow HH$ signal is also shown.

Minimum BDT value of the events passing the selection criteria of the resonant search. The combined BDT score is formed using as coefficients $C_{1}$ = 0.65 and $C_{2}$ = 1 − $C_{1}$. The selection efficiency for the resonant $X \rightarrow HH$ signal is also shown.

Minimum BDT value of the events passing the selection criteria of the resonant search. The combined BDT score is formed using as coefficients $C_{1}$ = 0.65 and $C_{2}$ = 1 − $C_{1}$. The selection efficiency for the resonant $X \rightarrow HH$ signal is also shown.

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.

Breakdown of expected and observed yields in the four edge signal regions, integrated over the $m_{\ell\ell}$ distribution after a separate simultaneous fit to each signal region and control region pair. The uncertainties include both the statistical and systematic sources.

Breakdown of expected and observed yields in the three on-$Z$ signal regions after a separate simultaneous fit to each signal region and control region pair. The uncertainties include both the statistical and systematic sources.

Distributions of $\mathcal{S}(E_{\mathrm{T}}^{\mathrm{miss}})$ in VR-High-Sideband-EWK (top-left), VR-High-R-EWK (top-right), VR-1J-High-EWK (bottom-left), and VR-$\ell\ell bb$-EWK (bottom-right) from the EWK search after a simultaneous fit of the control regions. The hatched band includes both the systematic and statistical uncertainties. The last bin includes the overflow.

Distributions of $\mathcal{S}(E_{\mathrm{T}}^{\mathrm{miss}})$ in VR-High-Sideband-EWK (top-left), VR-High-R-EWK (top-right), VR-1J-High-EWK (bottom-left), and VR-$\ell\ell bb$-EWK (bottom-right) from the EWK search after a simultaneous fit of the control regions. The hatched band includes both the systematic and statistical uncertainties. The last bin includes the overflow.

Distributions of $\mathcal{S}(E_{\mathrm{T}}^{\mathrm{miss}})$ in VR-High-Sideband-EWK (top-left), VR-High-R-EWK (top-right), VR-1J-High-EWK (bottom-left), and VR-$\ell\ell bb$-EWK (bottom-right) from the EWK search after a simultaneous fit of the control regions. The hatched band includes both the systematic and statistical uncertainties. The last bin includes the overflow.

Distributions of $\mathcal{S}(E_{\mathrm{T}}^{\mathrm{miss}})$ in VR-High-Sideband-EWK (top-left), VR-High-R-EWK (top-right), VR-1J-High-EWK (bottom-left), and VR-$\ell\ell bb$-EWK (bottom-right) from the EWK search after a simultaneous fit of the control regions. The hatched band includes both the systematic and statistical uncertainties. The last bin includes the overflow.

Distributions of $\mathcal{S}(E_{\mathrm{T}}^{\mathrm{miss}})$ in VR-Int-EWK (top-left), VR-Low-EWK (top-right), VR-Low-2-EWK (bottom-left), and VR-OffShell-EWK (bottom-right) from the EWK search after a simultaneous fit of the control regions. The hatched band includes both the systematic and statistical uncertainties. The last bin includes the overflow.

Distributions of $\mathcal{S}(E_{\mathrm{T}}^{\mathrm{miss}})$ in VR-Int-EWK (top-left), VR-Low-EWK (top-right), VR-Low-2-EWK (bottom-left), and VR-OffShell-EWK (bottom-right) from the EWK search after a simultaneous fit of the control regions. The hatched band includes both the systematic and statistical uncertainties. The last bin includes the overflow.

Distributions of $\mathcal{S}(E_{\mathrm{T}}^{\mathrm{miss}})$ in VR-Int-EWK (top-left), VR-Low-EWK (top-right), VR-Low-2-EWK (bottom-left), and VR-OffShell-EWK (bottom-right) from the EWK search after a simultaneous fit of the control regions. The hatched band includes both the systematic and statistical uncertainties. The last bin includes the overflow.

Distributions of $\mathcal{S}(E_{\mathrm{T}}^{\mathrm{miss}})$ in VR-Int-EWK (top-left), VR-Low-EWK (top-right), VR-Low-2-EWK (bottom-left), and VR-OffShell-EWK (bottom-right) from the EWK search after a simultaneous fit of the control regions. The hatched band includes both the systematic and statistical uncertainties. The last bin includes the overflow.

Observed and expected dilepton mass distributions in VRC-STR (top-left), VRLow-STR (top-right), VRMed-STR (bottom-left), and VRHigh-STR (bottom-right). Each validation region is fit separately with the corresponding control region. All statistical and systematic uncertainties are included in the hatched band. The entries are normalized to the bin width, and the last bin is the overflow.

Observed and expected dilepton mass distributions in VRC-STR (top-left), VRLow-STR (top-right), VRMed-STR (bottom-left), and VRHigh-STR (bottom-right). Each validation region is fit separately with the corresponding control region. All statistical and systematic uncertainties are included in the hatched band. The entries are normalized to the bin width, and the last bin is the overflow.

Observed and expected dilepton mass distributions in VRC-STR (top-left), VRLow-STR (top-right), VRMed-STR (bottom-left), and VRHigh-STR (bottom-right). Each validation region is fit separately with the corresponding control region. All statistical and systematic uncertainties are included in the hatched band. The entries are normalized to the bin width, and the last bin is the overflow.

Observed and expected dilepton mass distributions in VRC-STR (top-left), VRLow-STR (top-right), VRMed-STR (bottom-left), and VRHigh-STR (bottom-right). Each validation region is fit separately with the corresponding control region. All statistical and systematic uncertainties are included in the hatched band. The entries are normalized to the bin width, and the last bin is the overflow.

Observed and expected jet multiplicity in VRLow-STR (top-left), VRMed-STR (top-right), and VRHigh-STR (bottom) after a fit performed on the $m_{\ell\ell}$ distribution and corresponding control region. All statistical and systematic uncertainties are included in the hatched band. The last bin contains the overflow.

Observed and expected jet multiplicity in VRLow-STR (top-left), VRMed-STR (top-right), and VRHigh-STR (bottom) after a fit performed on the $m_{\ell\ell}$ distribution and corresponding control region. All statistical and systematic uncertainties are included in the hatched band. The last bin contains the overflow.

Observed and expected jet multiplicity in VRLow-STR (top-left), VRMed-STR (top-right), and VRHigh-STR (bottom) after a fit performed on the $m_{\ell\ell}$ distribution and corresponding control region. All statistical and systematic uncertainties are included in the hatched band. The last bin contains the overflow.

Observed and expected dilepton mass distributions in VR3L-STR without a fit to the data. The 'Other' category includes the negligible contributions from $t\bar{t}$ and $Z/\gamma^*$+jets processes. The hatched band contains the statistical uncertainty and the theoretical systematic uncertainties of the $WZ$/$ZZ$ prediction, which are the dominant sources of uncertainty. No fit is performed. The last bin contains the overflow.

Observed and expected distributions in five EWK search regions after a simultaneous fit to the SR and CR. In the top row, left-to-right, are $\mathcal{S}(E_{\mathrm{T}}^{\mathrm{miss}})$ in SR-High_8-EWK and $m_{bb}$ in SR-$\ell\ell bb$-EWK. In the middle row, left-to-right, are $\mathcal{S}(E_{\mathrm{T}}^{\mathrm{miss}})$ in SR-Int-EWK and $\mathcal{S}(E_{\mathrm{T}}^{\mathrm{miss}})$ in SR-Low-EWK. In the bottom row is $m_{\ell\ell}$ in SR-OffShell-EWK. Overlaid are example C1N2 and GMSB signal models, where the numbers in the brackets indicate the masses, in $\mathrm{GeV}$, of the $\tilde{\chi}_1^\pm$ and $\tilde{\chi}_2^0$ or the mass of the $\tilde{\chi}_1^0$ and branching ratio to the Higgs boson respectively. All statistical and systematic uncertainties are included in the hatched bands. The last bin includes the overflow.

Observed and expected distributions in five EWK search regions after a simultaneous fit to the SR and CR. In the top row, left-to-right, are $\mathcal{S}(E_{\mathrm{T}}^{\mathrm{miss}})$ in SR-High_8-EWK and $m_{bb}$ in SR-$\ell\ell bb$-EWK. In the middle row, left-to-right, are $\mathcal{S}(E_{\mathrm{T}}^{\mathrm{miss}})$ in SR-Int-EWK and $\mathcal{S}(E_{\mathrm{T}}^{\mathrm{miss}})$ in SR-Low-EWK. In the bottom row is $m_{\ell\ell}$ in SR-OffShell-EWK. Overlaid are example C1N2 and GMSB signal models, where the numbers in the brackets indicate the masses, in $\mathrm{GeV}$, of the $\tilde{\chi}_1^\pm$ and $\tilde{\chi}_2^0$ or the mass of the $\tilde{\chi}_1^0$ and branching ratio to the Higgs boson respectively. All statistical and systematic uncertainties are included in the hatched bands. The last bin includes the overflow.

Observed and expected distributions in five EWK search regions after a simultaneous fit to the SR and CR. In the top row, left-to-right, are $\mathcal{S}(E_{\mathrm{T}}^{\mathrm{miss}})$ in SR-High_8-EWK and $m_{bb}$ in SR-$\ell\ell bb$-EWK. In the middle row, left-to-right, are $\mathcal{S}(E_{\mathrm{T}}^{\mathrm{miss}})$ in SR-Int-EWK and $\mathcal{S}(E_{\mathrm{T}}^{\mathrm{miss}})$ in SR-Low-EWK. In the bottom row is $m_{\ell\ell}$ in SR-OffShell-EWK. Overlaid are example C1N2 and GMSB signal models, where the numbers in the brackets indicate the masses, in $\mathrm{GeV}$, of the $\tilde{\chi}_1^\pm$ and $\tilde{\chi}_2^0$ or the mass of the $\tilde{\chi}_1^0$ and branching ratio to the Higgs boson respectively. All statistical and systematic uncertainties are included in the hatched bands. The last bin includes the overflow.

Observed and expected distributions in five EWK search regions after a simultaneous fit to the SR and CR. In the top row, left-to-right, are $\mathcal{S}(E_{\mathrm{T}}^{\mathrm{miss}})$ in SR-High_8-EWK and $m_{bb}$ in SR-$\ell\ell bb$-EWK. In the middle row, left-to-right, are $\mathcal{S}(E_{\mathrm{T}}^{\mathrm{miss}})$ in SR-Int-EWK and $\mathcal{S}(E_{\mathrm{T}}^{\mathrm{miss}})$ in SR-Low-EWK. In the bottom row is $m_{\ell\ell}$ in SR-OffShell-EWK. Overlaid are example C1N2 and GMSB signal models, where the numbers in the brackets indicate the masses, in $\mathrm{GeV}$, of the $\tilde{\chi}_1^\pm$ and $\tilde{\chi}_2^0$ or the mass of the $\tilde{\chi}_1^0$ and branching ratio to the Higgs boson respectively. All statistical and systematic uncertainties are included in the hatched bands. The last bin includes the overflow.

Observed and expected distributions in five EWK search regions after a simultaneous fit to the SR and CR. In the top row, left-to-right, are $\mathcal{S}(E_{\mathrm{T}}^{\mathrm{miss}})$ in SR-High_8-EWK and $m_{bb}$ in SR-$\ell\ell bb$-EWK. In the middle row, left-to-right, are $\mathcal{S}(E_{\mathrm{T}}^{\mathrm{miss}})$ in SR-Int-EWK and $\mathcal{S}(E_{\mathrm{T}}^{\mathrm{miss}})$ in SR-Low-EWK. In the bottom row is $m_{\ell\ell}$ in SR-OffShell-EWK. Overlaid are example C1N2 and GMSB signal models, where the numbers in the brackets indicate the masses, in $\mathrm{GeV}$, of the $\tilde{\chi}_1^\pm$ and $\tilde{\chi}_2^0$ or the mass of the $\tilde{\chi}_1^0$ and branching ratio to the Higgs boson respectively. All statistical and systematic uncertainties are included in the hatched bands. The last bin includes the overflow.

Observed and expected dilepton mass distributions in SRC-STR (top-left), SRLow-STR (top-right), SRMed-STR (bottom-left), and SRHigh-STR (bottom-right), with the binning used for interpretations after a separate simultaneous fit to each signal region and control region pair. The red dashed lines are example signal models overlaid on the figure. All statistical and systematic uncertainties are included in the hatched bands. The last bins are the overflow.

Observed and expected dilepton mass distributions in SRC-STR (top-left), SRLow-STR (top-right), SRMed-STR (bottom-left), and SRHigh-STR (bottom-right), with the binning used for interpretations after a separate simultaneous fit to each signal region and control region pair. The red dashed lines are example signal models overlaid on the figure. All statistical and systematic uncertainties are included in the hatched bands. The last bins are the overflow.

Observed and expected dilepton mass distributions in SRC-STR (top-left), SRLow-STR (top-right), SRMed-STR (bottom-left), and SRHigh-STR (bottom-right), with the binning used for interpretations after a separate simultaneous fit to each signal region and control region pair. The red dashed lines are example signal models overlaid on the figure. All statistical and systematic uncertainties are included in the hatched bands. The last bins are the overflow.

Observed and expected dilepton mass distributions in SRC-STR (top-left), SRLow-STR (top-right), SRMed-STR (bottom-left), and SRHigh-STR (bottom-right), with the binning used for interpretations after a separate simultaneous fit to each signal region and control region pair. The red dashed lines are example signal models overlaid on the figure. All statistical and systematic uncertainties are included in the hatched bands. The last bins are the overflow.

Expected and observed exclusion contours from the EWK analysis for the C1N2 model (left) and GMSB model (right). The dashed line indicates the expected limits at 95$\%$ CL and the surrounding band shows the $1\sigma$ variation of the expected limit as a consequence of the uncertainties on the background prediction and experimental uncertainties on the signal ($\pm1\sigma_\mathrm{exp}$). The red dotted lines surrounding the observed limit contours indicate the variation resulting from changing the signal cross-section within its uncertainty ($\pm1\sigma_\mathrm{theory}^\mathrm{SUSY}$). The gray shaded areas indicate observed limits on these models from the two lepton channels of Ref.~[arXiv: 1803.02762] and Ref.~[arXiv: 1403.5294].

Expected and observed exclusion contours from the EWK analysis for the C1N2 model (left) and GMSB model (right). The dashed line indicates the expected limits at 95$\%$ CL and the surrounding band shows the $1\sigma$ variation of the expected limit as a consequence of the uncertainties on the background prediction and experimental uncertainties on the signal ($\pm1\sigma_\mathrm{exp}$). The red dotted lines surrounding the observed limit contours indicate the variation resulting from changing the signal cross-section within its uncertainty ($\pm1\sigma_\mathrm{theory}^\mathrm{SUSY}$). The gray shaded areas indicate observed limits on these models from the two lepton channels of Ref.~[arXiv: 1803.02762] and Ref.~[arXiv: 1403.5294].

Expected and observed exclusion contours from the EWK analysis for the C1N2 model (left) and GMSB model (right). The dashed line indicates the expected limits at 95$\%$ CL and the surrounding band shows the $1\sigma$ variation of the expected limit as a consequence of the uncertainties on the background prediction and experimental uncertainties on the signal ($\pm1\sigma_\mathrm{exp}$). The red dotted lines surrounding the observed limit contours indicate the variation resulting from changing the signal cross-section within its uncertainty ($\pm1\sigma_\mathrm{theory}^\mathrm{SUSY}$). The gray shaded areas indicate observed limits on these models from the two lepton channels of Ref.~[arXiv: 1803.02762] and Ref.~[arXiv: 1403.5294].

Expected and observed exclusion contours from the EWK analysis for the C1N2 model (left) and GMSB model (right). The dashed line indicates the expected limits at 95$\%$ CL and the surrounding band shows the $1\sigma$ variation of the expected limit as a consequence of the uncertainties on the background prediction and experimental uncertainties on the signal ($\pm1\sigma_\mathrm{exp}$). The red dotted lines surrounding the observed limit contours indicate the variation resulting from changing the signal cross-section within its uncertainty ($\pm1\sigma_\mathrm{theory}^\mathrm{SUSY}$). The gray shaded areas indicate observed limits on these models from the two lepton channels of Ref.~[arXiv: 1803.02762] and Ref.~[arXiv: 1403.5294].

Expected and observed exclusion contours from the EWK analysis for the C1N2 model (left) and GMSB model (right). The dashed line indicates the expected limits at 95$\%$ CL and the surrounding band shows the $1\sigma$ variation of the expected limit as a consequence of the uncertainties on the background prediction and experimental uncertainties on the signal ($\pm1\sigma_\mathrm{exp}$). The red dotted lines surrounding the observed limit contours indicate the variation resulting from changing the signal cross-section within its uncertainty ($\pm1\sigma_\mathrm{theory}^\mathrm{SUSY}$). The gray shaded areas indicate observed limits on these models from the two lepton channels of Ref.~[arXiv: 1803.02762] and Ref.~[arXiv: 1403.5294].

Expected and observed exclusion contours from the EWK analysis for the C1N2 model (left) and GMSB model (right). The dashed line indicates the expected limits at 95$\%$ CL and the surrounding band shows the $1\sigma$ variation of the expected limit as a consequence of the uncertainties on the background prediction and experimental uncertainties on the signal ($\pm1\sigma_\mathrm{exp}$). The red dotted lines surrounding the observed limit contours indicate the variation resulting from changing the signal cross-section within its uncertainty ($\pm1\sigma_\mathrm{theory}^\mathrm{SUSY}$). The gray shaded areas indicate observed limits on these models from the two lepton channels of Ref.~[arXiv: 1803.02762] and Ref.~[arXiv: 1403.5294].

Expected and observed exclusion contours from the EWK analysis for the C1N2 model (left) and GMSB model (right). The dashed line indicates the expected limits at 95$\%$ CL and the surrounding band shows the $1\sigma$ variation of the expected limit as a consequence of the uncertainties on the background prediction and experimental uncertainties on the signal ($\pm1\sigma_\mathrm{exp}$). The red dotted lines surrounding the observed limit contours indicate the variation resulting from changing the signal cross-section within its uncertainty ($\pm1\sigma_\mathrm{theory}^\mathrm{SUSY}$). The gray shaded areas indicate observed limits on these models from the two lepton channels of Ref.~[arXiv: 1803.02762] and Ref.~[arXiv: 1403.5294].

Expected and observed exclusion contours from the EWK analysis for the C1N2 model (left) and GMSB model (right). The dashed line indicates the expected limits at 95$\%$ CL and the surrounding band shows the $1\sigma$ variation of the expected limit as a consequence of the uncertainties on the background prediction and experimental uncertainties on the signal ($\pm1\sigma_\mathrm{exp}$). The red dotted lines surrounding the observed limit contours indicate the variation resulting from changing the signal cross-section within its uncertainty ($\pm1\sigma_\mathrm{theory}^\mathrm{SUSY}$). The gray shaded areas indicate observed limits on these models from the two lepton channels of Ref.~[arXiv: 1803.02762] and Ref.~[arXiv: 1403.5294].

Expected and observed exclusion contours from the EWK analysis for the C1N2 model (left) and GMSB model (right). The dashed line indicates the expected limits at 95$\%$ CL and the surrounding band shows the $1\sigma$ variation of the expected limit as a consequence of the uncertainties on the background prediction and experimental uncertainties on the signal ($\pm1\sigma_\mathrm{exp}$). The red dotted lines surrounding the observed limit contours indicate the variation resulting from changing the signal cross-section within its uncertainty ($\pm1\sigma_\mathrm{theory}^\mathrm{SUSY}$). The gray shaded areas indicate observed limits on these models from the two lepton channels of Ref.~[arXiv: 1803.02762] and Ref.~[arXiv: 1403.5294].

Expected and observed exclusion contours from the EWK analysis for the C1N2 model (left) and GMSB model (right). The dashed line indicates the expected limits at 95$\%$ CL and the surrounding band shows the $1\sigma$ variation of the expected limit as a consequence of the uncertainties on the background prediction and experimental uncertainties on the signal ($\pm1\sigma_\mathrm{exp}$). The red dotted lines surrounding the observed limit contours indicate the variation resulting from changing the signal cross-section within its uncertainty ($\pm1\sigma_\mathrm{theory}^\mathrm{SUSY}$). The gray shaded areas indicate observed limits on these models from the two lepton channels of Ref.~[arXiv: 1803.02762] and Ref.~[arXiv: 1403.5294].

Expected and observed exclusion contours from the EWK analysis for the C1N2 model (left) and GMSB model (right). The dashed line indicates the expected limits at 95$\%$ CL and the surrounding band shows the $1\sigma$ variation of the expected limit as a consequence of the uncertainties on the background prediction and experimental uncertainties on the signal ($\pm1\sigma_\mathrm{exp}$). The red dotted lines surrounding the observed limit contours indicate the variation resulting from changing the signal cross-section within its uncertainty ($\pm1\sigma_\mathrm{theory}^\mathrm{SUSY}$). The gray shaded areas indicate observed limits on these models from the two lepton channels of Ref.~[arXiv: 1803.02762] and Ref.~[arXiv: 1403.5294].

Expected and observed exclusion contours from the EWK analysis for the C1N2 model (left) and GMSB model (right). The dashed line indicates the expected limits at 95$\%$ CL and the surrounding band shows the $1\sigma$ variation of the expected limit as a consequence of the uncertainties on the background prediction and experimental uncertainties on the signal ($\pm1\sigma_\mathrm{exp}$). The red dotted lines surrounding the observed limit contours indicate the variation resulting from changing the signal cross-section within its uncertainty ($\pm1\sigma_\mathrm{theory}^\mathrm{SUSY}$). The gray shaded areas indicate observed limits on these models from the two lepton channels of Ref.~[arXiv: 1803.02762] and Ref.~[arXiv: 1403.5294].

Expected and observed exclusion contours from the EWK analysis for the C1N2 model (left) and GMSB model (right). The dashed line indicates the expected limits at 95$\%$ CL and the surrounding band shows the $1\sigma$ variation of the expected limit as a consequence of the uncertainties on the background prediction and experimental uncertainties on the signal ($\pm1\sigma_\mathrm{exp}$). The red dotted lines surrounding the observed limit contours indicate the variation resulting from changing the signal cross-section within its uncertainty ($\pm1\sigma_\mathrm{theory}^\mathrm{SUSY}$). The gray shaded areas indicate observed limits on these models from the two lepton channels of Ref.~[arXiv: 1803.02762] and Ref.~[arXiv: 1403.5294]. The grey numbers indicate the observed 95\% CLs upper limit on the cross section.

Expected and observed exclusion contours from the EWK analysis for the C1N2 model (left) and GMSB model (right). The dashed line indicates the expected limits at 95$\%$ CL and the surrounding band shows the $1\sigma$ variation of the expected limit as a consequence of the uncertainties on the background prediction and experimental uncertainties on the signal ($\pm1\sigma_\mathrm{exp}$). The red dotted lines surrounding the observed limit contours indicate the variation resulting from changing the signal cross-section within its uncertainty ($\pm1\sigma_\mathrm{theory}^\mathrm{SUSY}$). The gray shaded areas indicate observed limits on these models from the two lepton channels of Ref.~[arXiv: 1803.02762] and Ref.~[arXiv: 1403.5294]. The grey numbers indicate the observed 95$\%$ CLs upper limit on the cross section.

Expected and observed exclusion contours derived from the combination of all of the Strong search SRs for the $\tilde{g}$--$ ilde{\ell}$ (top-left), $\tilde{g}$--$Z$ (top-right), and $\tilde{s}--Z$ (bottom) models. The dashed line indicates the expected limits at 95\% CL and the surrounding band shows the $1\sigma$ variation of the expected limit as a consequence of the uncertainties in the background prediction and experimental uncertainties of the signal ($\pm1\sigma_\mathrm{exp}$). The red dotted lines surrounding the observed limit contours indicate the variation resulting from changing the signal cross-section within its uncertainty ($\pm1\sigma_\mathrm{theory}^\mathrm{SUSY}$). The grey-shaded area indicates the observed limits on these models from Ref. [23].

Expected and observed exclusion contours derived from the combination of all of the Strong search SRs for the $\tilde{g}$--$ ilde{\ell}$ (top-left), $\tilde{g}$--$Z$ (top-right), and $\tilde{s}--Z$ (bottom) models. The dashed line indicates the expected limits at 95\% CL and the surrounding band shows the $1\sigma$ variation of the expected limit as a consequence of the uncertainties in the background prediction and experimental uncertainties of the signal ($\pm1\sigma_\mathrm{exp}$). The red dotted lines surrounding the observed limit contours indicate the variation resulting from changing the signal cross-section within its uncertainty ($\pm1\sigma_\mathrm{theory}^\mathrm{SUSY}$). The grey-shaded area indicates the observed limits on these models from Ref. [23].

Expected and observed exclusion contours derived from the combination of all of the Strong search SRs for the $\tilde{g}$--$ ilde{\ell}$ (top-left), $\tilde{g}$--$Z$ (top-right), and $\tilde{s}--Z$ (bottom) models. The dashed line indicates the expected limits at 95\% CL and the surrounding band shows the $1\sigma$ variation of the expected limit as a consequence of the uncertainties in the background prediction and experimental uncertainties of the signal ($\pm1\sigma_\mathrm{exp}$). The red dotted lines surrounding the observed limit contours indicate the variation resulting from changing the signal cross-section within its uncertainty ($\pm1\sigma_\mathrm{theory}^\mathrm{SUSY}$). The grey-shaded area indicates the observed limits on these models from Ref. [23].

Expected and observed exclusion contours derived from the combination of all of the Strong search SRs for the $\tilde{g}$--$ ilde{\ell}$ (top-left), $\tilde{g}$--$Z$ (top-right), and $\tilde{s}--Z$ (bottom) models. The dashed line indicates the expected limits at 95\% CL and the surrounding band shows the $1\sigma$ variation of the expected limit as a consequence of the uncertainties in the background prediction and experimental uncertainties of the signal ($\pm1\sigma_\mathrm{exp}$). The red dotted lines surrounding the observed limit contours indicate the variation resulting from changing the signal cross-section within its uncertainty ($\pm1\sigma_\mathrm{theory}^\mathrm{SUSY}$). The grey-shaded area indicates the observed limits on these models from Ref. [23].

Expected and observed exclusion contours derived from the combination of all of the Strong search SRs for the $\tilde{g}$--$\tilde{\ell}$ (top-left), $\tilde{g}$--$Z$ (top-right), and $\tilde{s}--Z$ (bottom) models. The dashed line indicates the expected limits at 95\% CL and the surrounding band shows the $1\sigma$ variation of the expected limit as a consequence of the uncertainties in the background prediction and experimental uncertainties of the signal ($\pm1\sigma_\mathrm{exp}$). The red dotted lines surrounding the observed limit contours indicate the variation resulting from changing the signal cross-section within its uncertainty ($\pm1\sigma_\mathrm{theory}^\mathrm{SUSY}$). The grey-shaded area indicates the observed limits on these models from Ref. [23].

Expected and observed exclusion contours derived from the combination of all of the Strong search SRs for the $\tilde{g}$--$\tilde{\ell}$ (top-left), $\tilde{g}$--$Z$ (top-right), and $\tilde{s}--Z$ (bottom) models. The dashed line indicates the expected limits at 95\% CL and the surrounding band shows the $1\sigma$ variation of the expected limit as a consequence of the uncertainties in the background prediction and experimental uncertainties of the signal ($\pm1\sigma_\mathrm{exp}$). The red dotted lines surrounding the observed limit contours indicate the variation resulting from changing the signal cross-section within its uncertainty ($\pm1\sigma_\mathrm{theory}^\mathrm{SUSY}$). The grey-shaded area indicates the observed limits on these models from Ref. [23].

Expected and observed exclusion contours derived from the combination of all of the Strong search SRs for the $\tilde{g}$--$\tilde{\ell}$ (top-left), $\tilde{g}$--$Z$ (top-right), and $\tilde{s}--Z$ (bottom) models. The dashed line indicates the expected limits at 95\% CL and the surrounding band shows the $1\sigma$ variation of the expected limit as a consequence of the uncertainties in the background prediction and experimental uncertainties of the signal ($\pm1\sigma_\mathrm{exp}$). The red dotted lines surrounding the observed limit contours indicate the variation resulting from changing the signal cross-section within its uncertainty ($\pm1\sigma_\mathrm{theory}^\mathrm{SUSY}$). The grey-shaded area indicates the observed limits on these models from Ref. [23].The grey numbers indicated the observed 95\% CL upper limit on the cross section.

Expected and observed exclusion contours derived from the combination of all of the Strong search SRs for the $\tilde{g}$--$\tilde{\ell}$ (top-left), $\tilde{g}$--$Z$ (top-right), and $\tilde{s}--Z$ (bottom) models. The dashed line indicates the expected limits at 95\% CL and the surrounding band shows the $1\sigma$ variation of the expected limit as a consequence of the uncertainties in the background prediction and experimental uncertainties of the signal ($\pm1\sigma_\mathrm{exp}$). The red dotted lines surrounding the observed limit contours indicate the variation resulting from changing the signal cross-section within its uncertainty ($\pm1\sigma_\mathrm{theory}^\mathrm{SUSY}$). The grey-shaded area indicates the observed limits on these models from Ref. [23].

Expected and observed exclusion contours derived from the combination of all of the Strong search SRs for the $\tilde{g}$--$\tilde{\ell}$ (top-left), $\tilde{g}$--$Z$ (top-right), and $\tilde{s}--Z$ (bottom) models. The dashed line indicates the expected limits at 95\% CL and the surrounding band shows the $1\sigma$ variation of the expected limit as a consequence of the uncertainties in the background prediction and experimental uncertainties of the signal ($\pm1\sigma_\mathrm{exp}$). The red dotted lines surrounding the observed limit contours indicate the variation resulting from changing the signal cross-section within its uncertainty ($\pm1\sigma_\mathrm{theory}^\mathrm{SUSY}$). The grey-shaded area indicates the observed limits on these models from Ref. [23].

Expected and observed exclusion contours derived from the combination of all of the Strong search SRs for the $\tilde{g}$--$\tilde{\ell}$ (top-left), $\tilde{g}$--$Z$ (top-right), and $\tilde{s}--Z$ (bottom) models. The dashed line indicates the expected limits at 95\% CL and the surrounding band shows the $1\sigma$ variation of the expected limit as a consequence of the uncertainties in the background prediction and experimental uncertainties of the signal ($\pm1\sigma_\mathrm{exp}$). The red dotted lines surrounding the observed limit contours indicate the variation resulting from changing the signal cross-section within its uncertainty ($\pm1\sigma_\mathrm{theory}^\mathrm{SUSY}$). The grey-shaded area indicates the observed limits on these models from Ref. [23].

Expected and observed exclusion contours derived from the combination of all of the Strong search SRs for the $\tilde{g}$--$\tilde{\ell}$ (top-left), $\tilde{g}$--$Z$ (top-right), and $\tilde{s}--Z$ (bottom) models. The dashed line indicates the expected limits at 95\% CL and the surrounding band shows the $1\sigma$ variation of the expected limit as a consequence of the uncertainties in the background prediction and experimental uncertainties of the signal ($\pm1\sigma_\mathrm{exp}$). The red dotted lines surrounding the observed limit contours indicate the variation resulting from changing the signal cross-section within its uncertainty ($\pm1\sigma_\mathrm{theory}^\mathrm{SUSY}$). The grey-shaded area indicates the observed limits on these models from Ref. [23].

Expected and observed exclusion contours derived from the combination of all of the Strong search SRs for the $\tilde{g}$--$\tilde{\ell}$ (top-left), $\tilde{g}$--$Z$ (top-right), and $\tilde{s}--Z$ (bottom) models. The dashed line indicates the expected limits at 95\% CL and the surrounding band shows the $1\sigma$ variation of the expected limit as a consequence of the uncertainties in the background prediction and experimental uncertainties of the signal ($\pm1\sigma_\mathrm{exp}$). The red dotted lines surrounding the observed limit contours indicate the variation resulting from changing the signal cross-section within its uncertainty ($\pm1\sigma_\mathrm{theory}^\mathrm{SUSY}$). The grey-shaded area indicates the observed limits on these models from Ref. [23].

Expected and observed exclusion contours derived from the combination of all of the Strong search SRs for the $\tilde{g}$--$\tilde{\ell}$ (top-left), $\tilde{g}$--$Z$ (top-right), and $\tilde{s}--Z$ (bottom) models. The dashed line indicates the expected limits at 95\% CL and the surrounding band shows the $1\sigma$ variation of the expected limit as a consequence of the uncertainties in the background prediction and experimental uncertainties of the signal ($\pm1\sigma_\mathrm{exp}$). The red dotted lines surrounding the observed limit contours indicate the variation resulting from changing the signal cross-section within its uncertainty ($\pm1\sigma_\mathrm{theory}^\mathrm{SUSY}$). The grey-shaded area indicates the observed limits on these models from Ref. [23].

Expected and observed exclusion contours derived from the combination of all of the Strong search SRs for the $\tilde{g}$--$\tilde{\ell}$ (top-left), $\tilde{g}$--$Z$ (top-right), and $\tilde{s}--Z$ (bottom) models. The dashed line indicates the expected limits at 95\% CL and the surrounding band shows the $1\sigma$ variation of the expected limit as a consequence of the uncertainties in the background prediction and experimental uncertainties of the signal ($\pm1\sigma_\mathrm{exp}$). The red dotted lines surrounding the observed limit contours indicate the variation resulting from changing the signal cross-section within its uncertainty ($\pm1\sigma_\mathrm{theory}^\mathrm{SUSY}$). The grey-shaded area indicates the observed limits on these models from Ref. [23].The grey numbers indicated the observed 95\% CL upper limit on the cross section.

Expected and observed exclusion contours derived from the combination of all of the Strong search SRs for the $\tilde{g}$--$\tilde{\ell}$ (top-left), $\tilde{g}$--$Z$ (top-right), and $\tilde{s}--Z$ (bottom) models. The dashed line indicates the expected limits at 95\% CL and the surrounding band shows the $1\sigma$ variation of the expected limit as a consequence of the uncertainties in the background prediction and experimental uncertainties of the signal ($\pm1\sigma_\mathrm{exp}$). The red dotted lines surrounding the observed limit contours indicate the variation resulting from changing the signal cross-section within its uncertainty ($\pm1\sigma_\mathrm{theory}^\mathrm{SUSY}$). The grey-shaded area indicates the observed limits on these models from Ref. [23].

Expected and observed exclusion contours derived from the combination of all of the Strong search SRs for the $\tilde{g}$--$\tilde{\ell}$ (top-left), $\tilde{g}$--$Z$ (top-right), and $\tilde{s}--Z$ (bottom) models. The dashed line indicates the expected limits at 95\% CL and the surrounding band shows the $1\sigma$ variation of the expected limit as a consequence of the uncertainties in the background prediction and experimental uncertainties of the signal ($\pm1\sigma_\mathrm{exp}$). The red dotted lines surrounding the observed limit contours indicate the variation resulting from changing the signal cross-section within its uncertainty ($\pm1\sigma_\mathrm{theory}^\mathrm{SUSY}$). The grey-shaded area indicates the observed limits on these models from Ref. [23].

Expected and observed exclusion contours derived from the combination of all of the Strong search SRs for the $\tilde{g}$--$\tilde{\ell}$ (top-left), $\tilde{g}$--$Z$ (top-right), and $\tilde{s}--Z$ (bottom) models. The dashed line indicates the expected limits at 95\% CL and the surrounding band shows the $1\sigma$ variation of the expected limit as a consequence of the uncertainties in the background prediction and experimental uncertainties of the signal ($\pm1\sigma_\mathrm{exp}$). The red dotted lines surrounding the observed limit contours indicate the variation resulting from changing the signal cross-section within its uncertainty ($\pm1\sigma_\mathrm{theory}^\mathrm{SUSY}$). The grey-shaded area indicates the observed limits on these models from Ref. [23].

Expected and observed exclusion contours derived from the combination of all of the Strong search SRs for the $\tilde{g}$--$\tilde{\ell}$ (top-left), $\tilde{g}$--$Z$ (top-right), and $\tilde{s}--Z$ (bottom) models. The dashed line indicates the expected limits at 95\% CL and the surrounding band shows the $1\sigma$ variation of the expected limit as a consequence of the uncertainties in the background prediction and experimental uncertainties of the signal ($\pm1\sigma_\mathrm{exp}$). The red dotted lines surrounding the observed limit contours indicate the variation resulting from changing the signal cross-section within its uncertainty ($\pm1\sigma_\mathrm{theory}^\mathrm{SUSY}$). The grey-shaded area indicates the observed limits on these models from Ref. [23].

Expected and observed exclusion contours derived from the combination of all of the Strong search SRs for the $\tilde{g}$--$\tilde{\ell}$ (top-left), $\tilde{g}$--$Z$ (top-right), and $\tilde{s}--Z$ (bottom) models. The dashed line indicates the expected limits at 95\% CL and the surrounding band shows the $1\sigma$ variation of the expected limit as a consequence of the uncertainties in the background prediction and experimental uncertainties of the signal ($\pm1\sigma_\mathrm{exp}$). The red dotted lines surrounding the observed limit contours indicate the variation resulting from changing the signal cross-section within its uncertainty ($\pm1\sigma_\mathrm{theory}^\mathrm{SUSY}$). The grey-shaded area indicates the observed limits on these models from Ref. [23].The grey numbers indicated the observed 95\% CL upper limit on the cross section.

Expected and observed exclusion contours derived from the combination of all of the Strong search SRs for the $\tilde{g}$--$\tilde{\ell}$ (top-left), $\tilde{g}$--$Z$ (top-right), and $\tilde{s}--Z$ (bottom) models. The dashed line indicates the expected limits at 95\% CL and the surrounding band shows the $1\sigma$ variation of the expected limit as a consequence of the uncertainties in the background prediction and experimental uncertainties of the signal ($\pm1\sigma_\mathrm{exp}$). The red dotted lines surrounding the observed limit contours indicate the variation resulting from changing the signal cross-section within its uncertainty ($\pm1\sigma_\mathrm{theory}^\mathrm{SUSY}$). The grey-shaded area indicates the observed limits on these models from Ref. [23].

Expected and observed exclusion contours derived from the combination of all of the Strong search SRs for the $\tilde{g}$--$\tilde{\ell}$ (top-left), $\tilde{g}$--$Z$ (top-right), and $\tilde{s}--Z$ (bottom) models. The dashed line indicates the expected limits at 95\% CL and the surrounding band shows the $1\sigma$ variation of the expected limit as a consequence of the uncertainties in the background prediction and experimental uncertainties of the signal ($\pm1\sigma_\mathrm{exp}$). The red dotted lines surrounding the observed limit contours indicate the variation resulting from changing the signal cross-section within its uncertainty ($\pm1\sigma_\mathrm{theory}^\mathrm{SUSY}$). The grey-shaded area indicates the observed limits on these models from Ref. [23].

The combined $E_{\mathrm{T}}^{\mathrm{miss}}$ distribution of VRLow-STR and SRLow-STR (left), and the same region with the $\Delta\phi(\boldsymbol{j}_{1,2},\boldsymbol{\mathit{p}}_{ ext{T}}^{ ext{miss}})<0.4$ requirement, used as a control region to normalize the $Z/\gamma^*+\mathrm{jets}$ process (right). Separate fits for the SR and VR are performed, as for the results in the paper, and the resulting distributions are merged. All statistical and systematic uncertainties are included in the hatched bands. The last bins contain the overflow.

Cutflow of expected events in the four Strong search edge signal regions. `Leptons' refers to electrons and muons only. The gluino-$Z^{(*)}$ model with $m_{ ilde{g}}=800~GeV$ and $m_{ ilde{\chi}_1^0}=700~GeV$ is used for SRC-STR with 60000 Monte Carlo (MC) events generated. The slepton-$Z^{(*)}$ model with $m_{ ilde{\ell}}=1200~GeV$ and $m_{ ilde{\chi}_1^0}=700~GeV$ is used for SRMed-STR with 30000 (MC) events generated. The gluino-slepton model with $m_{ ilde{g}}=2~TeV$ and $m_{ ilde{\ell}}=1.3~TeV$ is used for SRLow-STR and SRHigh-STR with 30000 MC events generated. The Generator Filter requires two 5~GeV leptons and 100~GeV of \met. The SUSY2 kernel requires at least two leptons with $p_{\mathrm{T}}>9~GeV$ or at least one lepton with $p_{\mathrm{T}}>25~GeV$ and a photon with $p_{\mathrm{T}}>40~GeV$, with all objects within $|\eta|=2.6$.

Cutflow of expected events in the four Strong search edge signal regions. `Leptons' refers to electrons and muons only. The gluino-$Z^{(*)}$ model with $m_{ ilde{g}}=800~GeV$ and $m_{ ilde{\chi}_1^0}=700~GeV$ is used for SRC-STR with 60000 Monte Carlo (MC) events generated. The slepton-$Z^{(*)}$ model with $m_{ ilde{\ell}}=1200~GeV$ and $m_{ ilde{\chi}_1^0}=700~GeV$ is used for SRMed-STR with 30000 (MC) events generated. The gluino-slepton model with $m_{ ilde{g}}=2~TeV$ and $m_{ ilde{\ell}}=1.3~TeV$ is used for SRLow-STR and SRHigh-STR with 30000 MC events generated. The Generator Filter requires two 5~GeV leptons and 100~GeV of \met. The SUSY2 kernel requires at least two leptons with $p_{\mathrm{T}}>9~GeV$ or at least one lepton with $p_{\mathrm{T}}>25~GeV$ and a photon with $p_{\mathrm{T}}>40~GeV$, with all objects within $|\eta|=2.6$.

Cutflow of expected events in the four Strong search edge signal regions. `Leptons' refers to electrons and muons only. The gluino-$Z^{(*)}$ model with $m_{ ilde{g}}=800~GeV$ and $m_{ ilde{\chi}_1^0}=700~GeV$ is used for SRC-STR with 60000 Monte Carlo (MC) events generated. The slepton-$Z^{(*)}$ model with $m_{ ilde{\ell}}=1200~GeV$ and $m_{ ilde{\chi}_1^0}=700~GeV$ is used for SRMed-STR with 30000 (MC) events generated. The gluino-slepton model with $m_{ ilde{g}}=2~TeV$ and $m_{ ilde{\ell}}=1.3~TeV$ is used for SRLow-STR and SRHigh-STR with 30000 MC events generated. The Generator Filter requires two 5~GeV leptons and 100~GeV of \met. The SUSY2 kernel requires at least two leptons with $p_{\mathrm{T}}>9~GeV$ or at least one lepton with $p_{\mathrm{T}}>25~GeV$ and a photon with $p_{\mathrm{T}}>40~GeV$, with all objects within $|\eta|=2.6$.

Cutflow of expected events in the four Strong search edge signal regions. `Leptons' refers to electrons and muons only. The gluino-$Z^{(*)}$ model with $m_{ ilde{g}}=800~GeV$ and $m_{ ilde{\chi}_1^0}=700~GeV$ is used for SRC-STR with 60000 Monte Carlo (MC) events generated. The slepton-$Z^{(*)}$ model with $m_{ ilde{\ell}}=1200~GeV$ and $m_{ ilde{\chi}_1^0}=700~GeV$ is used for SRMed-STR with 30000 (MC) events generated. The gluino-slepton model with $m_{ ilde{g}}=2~TeV$ and $m_{ ilde{\ell}}=1.3~TeV$ is used for SRLow-STR and SRHigh-STR with 30000 MC events generated. The Generator Filter requires two 5~GeV leptons and 100~GeV of \met. The SUSY2 kernel requires at least two leptons with $p_{\mathrm{T}}>9~GeV$ or at least one lepton with $p_{\mathrm{T}}>25~GeV$ and a photon with $p_{\mathrm{T}}>40~GeV$, with all objects within $|\eta|=2.6$.

Cutflow of expected events in the three Strong search on-$Z$ signal regions. The cutflow up to the signal region specific requirements is the same as in the Strong search edge cutflow. The slepton-$Z^{(*)}$ model with $m_{ ilde{\ell}}=1200~GeV$ and $m_{ ilde{\chi}_1^0}=700~GeV$ is used for all of the on-$Z$ signal regions with 30000 (MC) events generated.

Cutflow of expected events in the three Strong search on-$Z$ signal regions. The cutflow up to the signal region specific requirements is the same as in the Strong search edge cutflow. The slepton-$Z^{(*)}$ model with $m_{ ilde{\ell}}=1200~GeV$ and $m_{ ilde{\chi}_1^0}=700~GeV$ is used for all of the on-$Z$ signal regions with 30000 (MC) events generated.

Cutflow of expected events in the three Strong search on-$Z$ signal regions. The cutflow up to the signal region specific requirements is the same as in the Strong search edge cutflow. The slepton-$Z^{(*)}$ model with $m_{ ilde{\ell}}=1200~GeV$ and $m_{ ilde{\chi}_1^0}=700~GeV$ is used for all of the on-$Z$ signal regions with 30000 (MC) events generated.

Table 36: Cutflow of expected events in the region SR-OffShell_a-EWK. Requirements below the line are specific to this region. On the Generator Filter line, the total number of unweighted events simulated is given in brackets. `Leptons' refers to electrons and muons only. For C1N2 models, the Generator Filter requires at least two $7~\mathrm{GeV}$ leptons and for C1N2 models with mass splittings below the Z boson mass it also requires $75~\mathrm{GeV}$ of $E_{\mathrm{T}}^{\mathrm{miss}}$. For GMSB models, the Generator Filter requires at least two $3~\mathrm{GeV}$ leptons. For on-shell C1N2 models, the `Forced Decays' require each Z boson to decay to a charged lepton pair (electron, muon, or tau) and each W boson to decay hadronically. For off-shell C1N2 models, each neutralino is forced to produce a charged lepton pair in its decay, and each chargino can produce any fermion pair. The SUSY2 kernel requires at least two leptons with $p_{\mathrm{T}}>9~\mathrm{GeV}$ or at least one lepton with $p_{\mathrm{T}}>25~\mathrm{GeV}$ and a photon with $p_{\mathrm{T}}>40~\mathrm{GeV}$, with all objects within $|\eta|=2.6$.

Table 36: Cutflow of expected events in the region SR-OffShell_a-EWK. Requirements below the line are specific to this region. On the Generator Filter line, the total number of unweighted events simulated is given in brackets. `Leptons' refers to electrons and muons only. For C1N2 models, the Generator Filter requires at least two $7~\mathrm{GeV}$ leptons and for C1N2 models with mass splittings below the Z boson mass it also requires $75~\mathrm{GeV}$ of $E_{\mathrm{T}}^{\mathrm{miss}}$. For GMSB models, the Generator Filter requires at least two $3~\mathrm{GeV}$ leptons. For on-shell C1N2 models, the `Forced Decays' require each Z boson to decay to a charged lepton pair (electron, muon, or tau) and each W boson to decay hadronically. For off-shell C1N2 models, each neutralino is forced to produce a charged lepton pair in its decay, and each chargino can produce any fermion pair. The SUSY2 kernel requires at least two leptons with $p_{\mathrm{T}}>9~\mathrm{GeV}$ or at least one lepton with $p_{\mathrm{T}}>25~\mathrm{GeV}$ and a photon with $p_{\mathrm{T}}>40~\mathrm{GeV}$, with all objects within $|\eta|=2.6$.

Table 37: Cutflow of expected events in the region SR-OffShell_b-EWK. Requirements below the line are specific to this region. On the Generator Filter line, the total number of unweighted events simulated is given in brackets. `Leptons' refers to electrons and muons only. For C1N2 models, the Generator Filter requires at least two $7~\mathrm{GeV}$ leptons and for C1N2 models with mass splittings below the Z boson mass it also requires $75~\mathrm{GeV}$ of $E_{\mathrm{T}}^{\mathrm{miss}}$. For GMSB models, the Generator Filter requires at least two $3~\mathrm{GeV}$ leptons. For on-shell C1N2 models, the `Forced Decays' require each Z boson to decay to a charged lepton pair (electron, muon, or tau) and each W boson to decay hadronically. For off-shell C1N2 models, each neutralino is forced to produce a charged lepton pair in its decay, and each chargino can produce any fermion pair. The SUSY2 kernel requires at least two leptons with $p_{\mathrm{T}}>9~\mathrm{GeV}$ or at least one lepton with $p_{\mathrm{T}}>25~\mathrm{GeV}$ and a photon with $p_{\mathrm{T}}>40~\mathrm{GeV}$, with all objects within $|\eta|=2.6$.

Table 37: Cutflow of expected events in the region SR-OffShell_b-EWK. Requirements below the line are specific to this region. On the Generator Filter line, the total number of unweighted events simulated is given in brackets. `Leptons' refers to electrons and muons only. For C1N2 models, the Generator Filter requires at least two $7~\mathrm{GeV}$ leptons and for C1N2 models with mass splittings below the Z boson mass it also requires $75~\mathrm{GeV}$ of $E_{\mathrm{T}}^{\mathrm{miss}}$. For GMSB models, the Generator Filter requires at least two $3~\mathrm{GeV}$ leptons. For on-shell C1N2 models, the `Forced Decays' require each Z boson to decay to a charged lepton pair (electron, muon, or tau) and each W boson to decay hadronically. For off-shell C1N2 models, each neutralino is forced to produce a charged lepton pair in its decay, and each chargino can produce any fermion pair. The SUSY2 kernel requires at least two leptons with $p_{\mathrm{T}}>9~\mathrm{GeV}$ or at least one lepton with $p_{\mathrm{T}}>25~\mathrm{GeV}$ and a photon with $p_{\mathrm{T}}>40~\mathrm{GeV}$, with all objects within $|\eta|=2.6$.

Table 38: Cutflow of expected events in the region SR-Low_a-EWK. Requirements below the line are specific to this region. On the Generator Filter line, the total number of unweighted events simulated is given in brackets. `Leptons' refers to electrons and muons only. For C1N2 models, the Generator Filter requires at least two $7~\mathrm{GeV}$ leptons and for C1N2 models with mass splittings below the Z boson mass it also requires $75~\mathrm{GeV}$ of $E_{\mathrm{T}}^{\mathrm{miss}}$. For GMSB models, the Generator Filter requires at least two $3~\mathrm{GeV}$ leptons. For on-shell C1N2 models, the `Forced Decays' require each Z boson to decay to a charged lepton pair (electron, muon, or tau) and each W boson to decay hadronically. For off-shell C1N2 models, each neutralino is forced to produce a charged lepton pair in its decay, and each chargino can produce any fermion pair. The SUSY2 kernel requires at least two leptons with $p_{\mathrm{T}}>9~\mathrm{GeV}$ or at least one lepton with $p_{\mathrm{T}}>25~\mathrm{GeV}$ and a photon with $p_{\mathrm{T}}>40~\mathrm{GeV}$, with all objects within $|\eta|=2.6$.

Table 38: Cutflow of expected events in the region SR-Low_a-EWK. Requirements below the line are specific to this region. On the Generator Filter line, the total number of unweighted events simulated is given in brackets. `Leptons' refers to electrons and muons only. For C1N2 models, the Generator Filter requires at least two $7~\mathrm{GeV}$ leptons and for C1N2 models with mass splittings below the Z boson mass it also requires $75~\mathrm{GeV}$ of $E_{\mathrm{T}}^{\mathrm{miss}}$. For GMSB models, the Generator Filter requires at least two $3~\mathrm{GeV}$ leptons. For on-shell C1N2 models, the `Forced Decays' require each Z boson to decay to a charged lepton pair (electron, muon, or tau) and each W boson to decay hadronically. For off-shell C1N2 models, each neutralino is forced to produce a charged lepton pair in its decay, and each chargino can produce any fermion pair. The SUSY2 kernel requires at least two leptons with $p_{\mathrm{T}}>9~\mathrm{GeV}$ or at least one lepton with $p_{\mathrm{T}}>25~\mathrm{GeV}$ and a photon with $p_{\mathrm{T}}>40~\mathrm{GeV}$, with all objects within $|\eta|=2.6$.

Table 39: Cutflow of expected events in the region SR-Low_b-EWK. Requirements below the line are specific to this region. On the Generator Filter line, the total number of unweighted events simulated is given in brackets. `Leptons' refers to electrons and muons only. For C1N2 models, the Generator Filter requires at least two $7~\mathrm{GeV}$ leptons and for C1N2 models with mass splittings below the Z boson mass it also requires $75~\mathrm{GeV}$ of $E_{\mathrm{T}}^{\mathrm{miss}}$. For GMSB models, the Generator Filter requires at least two $3~\mathrm{GeV}$ leptons. For on-shell C1N2 models, the `Forced Decays' require each Z boson to decay to a charged lepton pair (electron, muon, or tau) and each W boson to decay hadronically. For off-shell C1N2 models, each neutralino is forced to produce a charged lepton pair in its decay, and each chargino can produce any fermion pair. The SUSY2 kernel requires at least two leptons with $p_{\mathrm{T}}>9~\mathrm{GeV}$ or at least one lepton with $p_{\mathrm{T}}>25~\mathrm{GeV}$ and a photon with $p_{\mathrm{T}}>40~\mathrm{GeV}$, with all objects within $|\eta|=2.6$.

Table 39: Cutflow of expected events in the region SR-Low_b-EWK. Requirements below the line are specific to this region. On the Generator Filter line, the total number of unweighted events simulated is given in brackets. `Leptons' refers to electrons and muons only. For C1N2 models, the Generator Filter requires at least two $7~\mathrm{GeV}$ leptons and for C1N2 models with mass splittings below the Z boson mass it also requires $75~\mathrm{GeV}$ of $E_{\mathrm{T}}^{\mathrm{miss}}$. For GMSB models, the Generator Filter requires at least two $3~\mathrm{GeV}$ leptons. For on-shell C1N2 models, the `Forced Decays' require each Z boson to decay to a charged lepton pair (electron, muon, or tau) and each W boson to decay hadronically. For off-shell C1N2 models, each neutralino is forced to produce a charged lepton pair in its decay, and each chargino can produce any fermion pair. The SUSY2 kernel requires at least two leptons with $p_{\mathrm{T}}>9~\mathrm{GeV}$ or at least one lepton with $p_{\mathrm{T}}>25~\mathrm{GeV}$ and a photon with $p_{\mathrm{T}}>40~\mathrm{GeV}$, with all objects within $|\eta|=2.6$.

Table 40: Cutflow of expected events in the region SR-Low-2-EWK. Requirements below the line are specific to this region. On the Generator Filter line, the total number of unweighted events simulated is given in brackets. `Leptons' refers to electrons and muons only. For C1N2 models, the Generator Filter requires at least two $7~\mathrm{GeV}$ leptons and for C1N2 models with mass splittings below the Z boson mass it also requires $75~\mathrm{GeV}$ of $E_{\mathrm{T}}^{\mathrm{miss}}$. For GMSB models, the Generator Filter requires at least two $3~\mathrm{GeV}$ leptons. For on-shell C1N2 models, the `Forced Decays' require each Z boson to decay to a charged lepton pair (electron, muon, or tau) and each W boson to decay hadronically. For off-shell C1N2 models, each neutralino is forced to produce a charged lepton pair in its decay, and each chargino can produce any fermion pair. The SUSY2 kernel requires at least two leptons with $p_{\mathrm{T}}>9~\mathrm{GeV}$ or at least one lepton with $p_{\mathrm{T}}>25~\mathrm{GeV}$ and a photon with $p_{\mathrm{T}}>40~\mathrm{GeV}$, with all objects within $|\eta|=2.6$.

Table 40: Cutflow of expected events in the region SR-Low-2-EWK. Requirements below the line are specific to this region. On the Generator Filter line, the total number of unweighted events simulated is given in brackets. `Leptons' refers to electrons and muons only. For C1N2 models, the Generator Filter requires at least two $7~\mathrm{GeV}$ leptons and for C1N2 models with mass splittings below the Z boson mass it also requires $75~\mathrm{GeV}$ of $E_{\mathrm{T}}^{\mathrm{miss}}$. For GMSB models, the Generator Filter requires at least two $3~\mathrm{GeV}$ leptons. For on-shell C1N2 models, the `Forced Decays' require each Z boson to decay to a charged lepton pair (electron, muon, or tau) and each W boson to decay hadronically. For off-shell C1N2 models, each neutralino is forced to produce a charged lepton pair in its decay, and each chargino can produce any fermion pair. The SUSY2 kernel requires at least two leptons with $p_{\mathrm{T}}>9~\mathrm{GeV}$ or at least one lepton with $p_{\mathrm{T}}>25~\mathrm{GeV}$ and a photon with $p_{\mathrm{T}}>40~\mathrm{GeV}$, with all objects within $|\eta|=2.6$.

Table 41: Cutflow of expected events in the region SR-Int_a-EWK. Requirements below the line are specific to this region. On the Generator Filter line, the total number of unweighted events simulated is given in brackets. `Leptons' refers to electrons and muons only. For C1N2 models, the Generator Filter requires at least two $7~\mathrm{GeV}$ leptons and for C1N2 models with mass splittings below the Z boson mass it also requires $75~\mathrm{GeV}$ of $E_{\mathrm{T}}^{\mathrm{miss}}$. For GMSB models, the Generator Filter requires at least two $3~\mathrm{GeV}$ leptons. For on-shell C1N2 models, the `Forced Decays' require each Z boson to decay to a charged lepton pair (electron, muon, or tau) and each W boson to decay hadronically. For off-shell C1N2 models, each neutralino is forced to produce a charged lepton pair in its decay, and each chargino can produce any fermion pair. The SUSY2 kernel requires at least two leptons with $p_{\mathrm{T}}>9~\mathrm{GeV}$ or at least one lepton with $p_{\mathrm{T}}>25~\mathrm{GeV}$ and a photon with $p_{\mathrm{T}}>40~\mathrm{GeV}$, with all objects within $|\eta|=2.6$.

Table 41: Cutflow of expected events in the region SR-Int_a-EWK. Requirements below the line are specific to this region. On the Generator Filter line, the total number of unweighted events simulated is given in brackets. `Leptons' refers to electrons and muons only. For C1N2 models, the Generator Filter requires at least two $7~\mathrm{GeV}$ leptons and for C1N2 models with mass splittings below the Z boson mass it also requires $75~\mathrm{GeV}$ of $E_{\mathrm{T}}^{\mathrm{miss}}$. For GMSB models, the Generator Filter requires at least two $3~\mathrm{GeV}$ leptons. For on-shell C1N2 models, the `Forced Decays' require each Z boson to decay to a charged lepton pair (electron, muon, or tau) and each W boson to decay hadronically. For off-shell C1N2 models, each neutralino is forced to produce a charged lepton pair in its decay, and each chargino can produce any fermion pair. The SUSY2 kernel requires at least two leptons with $p_{\mathrm{T}}>9~\mathrm{GeV}$ or at least one lepton with $p_{\mathrm{T}}>25~\mathrm{GeV}$ and a photon with $p_{\mathrm{T}}>40~\mathrm{GeV}$, with all objects within $|\eta|=2.6$.

Table 42: Cutflow of expected events in the region SR-Int_b-EWK. Requirements below the line are specific to this region. On the Generator Filter line, the total number of unweighted events simulated is given in brackets. `Leptons' refers to electrons and muons only. For C1N2 models, the Generator Filter requires at least two $7~\mathrm{GeV}$ leptons and for C1N2 models with mass splittings below the Z boson mass it also requires $75~\mathrm{GeV}$ of $E_{\mathrm{T}}^{\mathrm{miss}}$. For GMSB models, the Generator Filter requires at least two $3~\mathrm{GeV}$ leptons. For on-shell C1N2 models, the `Forced Decays' require each Z boson to decay to a charged lepton pair (electron, muon, or tau) and each W boson to decay hadronically. For off-shell C1N2 models, each neutralino is forced to produce a charged lepton pair in its decay, and each chargino can produce any fermion pair. The SUSY2 kernel requires at least two leptons with $p_{\mathrm{T}}>9~\mathrm{GeV}$ or at least one lepton with $p_{\mathrm{T}}>25~\mathrm{GeV}$ and a photon with $p_{\mathrm{T}}>40~\mathrm{GeV}$, with all objects within $|\eta|=2.6$.

Table 42: Cutflow of expected events in the region SR-Int_b-EWK. Requirements below the line are specific to this region. On the Generator Filter line, the total number of unweighted events simulated is given in brackets. `Leptons' refers to electrons and muons only. For C1N2 models, the Generator Filter requires at least two $7~\mathrm{GeV}$ leptons and for C1N2 models with mass splittings below the Z boson mass it also requires $75~\mathrm{GeV}$ of $E_{\mathrm{T}}^{\mathrm{miss}}$. For GMSB models, the Generator Filter requires at least two $3~\mathrm{GeV}$ leptons. For on-shell C1N2 models, the `Forced Decays' require each Z boson to decay to a charged lepton pair (electron, muon, or tau) and each W boson to decay hadronically. For off-shell C1N2 models, each neutralino is forced to produce a charged lepton pair in its decay, and each chargino can produce any fermion pair. The SUSY2 kernel requires at least two leptons with $p_{\mathrm{T}}>9~\mathrm{GeV}$ or at least one lepton with $p_{\mathrm{T}}>25~\mathrm{GeV}$ and a photon with $p_{\mathrm{T}}>40~\mathrm{GeV}$, with all objects within $|\eta|=2.6$.

Table 43: Cutflow of expected events in the region SR-High_16a-EWK. Requirements below the line are specific to this region. On the Generator Filter line, the total number of unweighted events simulated is given in brackets. `Leptons' refers to electrons and muons only. For C1N2 models, the Generator Filter requires at least two $7~\mathrm{GeV}$ leptons and for C1N2 models with mass splittings below the Z boson mass it also requires $75~\mathrm{GeV}$ of $E_{\mathrm{T}}^{\mathrm{miss}}$. For GMSB models, the Generator Filter requires at least two $3~\mathrm{GeV}$ leptons. For on-shell C1N2 models, the `Forced Decays' require each Z boson to decay to a charged lepton pair (electron, muon, or tau) and each W boson to decay hadronically. For off-shell C1N2 models, each neutralino is forced to produce a charged lepton pair in its decay, and each chargino can produce any fermion pair. The SUSY2 kernel requires at least two leptons with $p_{\mathrm{T}}>9~\mathrm{GeV}$ or at least one lepton with $p_{\mathrm{T}}>25~\mathrm{GeV}$ and a photon with $p_{\mathrm{T}}>40~\mathrm{GeV}$, with all objects within $|\eta|=2.6$.

Table 43: Cutflow of expected events in the region SR-High_16a-EWK. Requirements below the line are specific to this region. On the Generator Filter line, the total number of unweighted events simulated is given in brackets. `Leptons' refers to electrons and muons only. For C1N2 models, the Generator Filter requires at least two $7~\mathrm{GeV}$ leptons and for C1N2 models with mass splittings below the Z boson mass it also requires $75~\mathrm{GeV}$ of $E_{\mathrm{T}}^{\mathrm{miss}}$. For GMSB models, the Generator Filter requires at least two $3~\mathrm{GeV}$ leptons. For on-shell C1N2 models, the `Forced Decays' require each Z boson to decay to a charged lepton pair (electron, muon, or tau) and each W boson to decay hadronically. For off-shell C1N2 models, each neutralino is forced to produce a charged lepton pair in its decay, and each chargino can produce any fermion pair. The SUSY2 kernel requires at least two leptons with $p_{\mathrm{T}}>9~\mathrm{GeV}$ or at least one lepton with $p_{\mathrm{T}}>25~\mathrm{GeV}$ and a photon with $p_{\mathrm{T}}>40~\mathrm{GeV}$, with all objects within $|\eta|=2.6$.

Table 44: Cutflow of expected events in the region SR-High_16b-EWK. Requirements below the line are specific to this region. On the Generator Filter line, the total number of unweighted events simulated is given in brackets. `Leptons' refers to electrons and muons only. For C1N2 models, the Generator Filter requires at least two $7~\mathrm{GeV}$ leptons and for C1N2 models with mass splittings below the Z boson mass it also requires $75~\mathrm{GeV}$ of $E_{\mathrm{T}}^{\mathrm{miss}}$. For GMSB models, the Generator Filter requires at least two $3~\mathrm{GeV}$ leptons. For on-shell C1N2 models, the `Forced Decays' require each Z boson to decay to a charged lepton pair (electron, muon, or tau) and each W boson to decay hadronically. For off-shell C1N2 models, each neutralino is forced to produce a charged lepton pair in its decay, and each chargino can produce any fermion pair. The SUSY2 kernel requires at least two leptons with $p_{\mathrm{T}}>9~\mathrm{GeV}$ or at least one lepton with $p_{\mathrm{T}}>25~\mathrm{GeV}$ and a photon with $p_{\mathrm{T}}>40~\mathrm{GeV}$, with all objects within $|\eta|=2.6$.

Table 44: Cutflow of expected events in the region SR-High_16b-EWK. Requirements below the line are specific to this region. On the Generator Filter line, the total number of unweighted events simulated is given in brackets. `Leptons' refers to electrons and muons only. For C1N2 models, the Generator Filter requires at least two $7~\mathrm{GeV}$ leptons and for C1N2 models with mass splittings below the Z boson mass it also requires $75~\mathrm{GeV}$ of $E_{\mathrm{T}}^{\mathrm{miss}}$. For GMSB models, the Generator Filter requires at least two $3~\mathrm{GeV}$ leptons. For on-shell C1N2 models, the `Forced Decays' require each Z boson to decay to a charged lepton pair (electron, muon, or tau) and each W boson to decay hadronically. For off-shell C1N2 models, each neutralino is forced to produce a charged lepton pair in its decay, and each chargino can produce any fermion pair. The SUSY2 kernel requires at least two leptons with $p_{\mathrm{T}}>9~\mathrm{GeV}$ or at least one lepton with $p_{\mathrm{T}}>25~\mathrm{GeV}$ and a photon with $p_{\mathrm{T}}>40~\mathrm{GeV}$, with all objects within $|\eta|=2.6$.

Table 45: Cutflow of expected events in the region SR-High_8a-EWK. Requirements below the line are specific to this region. On the Generator Filter line, the total number of unweighted events simulated is given in brackets. `Leptons' refers to electrons and muons only. For C1N2 models, the Generator Filter requires at least two $7~\mathrm{GeV}$ leptons and for C1N2 models with mass splittings below the Z boson mass it also requires $75~\mathrm{GeV}$ of $E_{\mathrm{T}}^{\mathrm{miss}}$. For GMSB models, the Generator Filter requires at least two $3~\mathrm{GeV}$ leptons. For on-shell C1N2 models, the `Forced Decays' require each Z boson to decay to a charged lepton pair (electron, muon, or tau) and each W boson to decay hadronically. For off-shell C1N2 models, each neutralino is forced to produce a charged lepton pair in its decay, and each chargino can produce any fermion pair. The SUSY2 kernel requires at least two leptons with $p_{\mathrm{T}}>9~\mathrm{GeV}$ or at least one lepton with $p_{\mathrm{T}}>25~\mathrm{GeV}$ and a photon with $p_{\mathrm{T}}>40~\mathrm{GeV}$, with all objects within $|\eta|=2.6$.

Table 45: Cutflow of expected events in the region SR-High_8a-EWK. Requirements below the line are specific to this region. On the Generator Filter line, the total number of unweighted events simulated is given in brackets. `Leptons' refers to electrons and muons only. For C1N2 models, the Generator Filter requires at least two $7~\mathrm{GeV}$ leptons and for C1N2 models with mass splittings below the Z boson mass it also requires $75~\mathrm{GeV}$ of $E_{\mathrm{T}}^{\mathrm{miss}}$. For GMSB models, the Generator Filter requires at least two $3~\mathrm{GeV}$ leptons. For on-shell C1N2 models, the `Forced Decays' require each Z boson to decay to a charged lepton pair (electron, muon, or tau) and each W boson to decay hadronically. For off-shell C1N2 models, each neutralino is forced to produce a charged lepton pair in its decay, and each chargino can produce any fermion pair. The SUSY2 kernel requires at least two leptons with $p_{\mathrm{T}}>9~\mathrm{GeV}$ or at least one lepton with $p_{\mathrm{T}}>25~\mathrm{GeV}$ and a photon with $p_{\mathrm{T}}>40~\mathrm{GeV}$, with all objects within $|\eta|=2.6$.

Table 46: Cutflow of expected events in the region SR-High_8b-EWK. Requirements below the line are specific to this region. On the Generator Filter line, the total number of unweighted events simulated is given in brackets. `Leptons' refers to electrons and muons only. For C1N2 models, the Generator Filter requires at least two $7~\mathrm{GeV}$ leptons and for C1N2 models with mass splittings below the Z boson mass it also requires $75~\mathrm{GeV}$ of $E_{\mathrm{T}}^{\mathrm{miss}}$. For GMSB models, the Generator Filter requires at least two $3~\mathrm{GeV}$ leptons. For on-shell C1N2 models, the `Forced Decays' require each Z boson to decay to a charged lepton pair (electron, muon, or tau) and each W boson to decay hadronically. For off-shell C1N2 models, each neutralino is forced to produce a charged lepton pair in its decay, and each chargino can produce any fermion pair. The SUSY2 kernel requires at least two leptons with $p_{\mathrm{T}}>9~\mathrm{GeV}$ or at least one lepton with $p_{\mathrm{T}}>25~\mathrm{GeV}$ and a photon with $p_{\mathrm{T}}>40~\mathrm{GeV}$, with all objects within $|\eta|=2.6$.

Table 46: Cutflow of expected events in the region SR-High_8b-EWK. Requirements below the line are specific to this region. On the Generator Filter line, the total number of unweighted events simulated is given in brackets. `Leptons' refers to electrons and muons only. For C1N2 models, the Generator Filter requires at least two $7~\mathrm{GeV}$ leptons and for C1N2 models with mass splittings below the Z boson mass it also requires $75~\mathrm{GeV}$ of $E_{\mathrm{T}}^{\mathrm{miss}}$. For GMSB models, the Generator Filter requires at least two $3~\mathrm{GeV}$ leptons. For on-shell C1N2 models, the `Forced Decays' require each Z boson to decay to a charged lepton pair (electron, muon, or tau) and each W boson to decay hadronically. For off-shell C1N2 models, each neutralino is forced to produce a charged lepton pair in its decay, and each chargino can produce any fermion pair. The SUSY2 kernel requires at least two leptons with $p_{\mathrm{T}}>9~\mathrm{GeV}$ or at least one lepton with $p_{\mathrm{T}}>25~\mathrm{GeV}$ and a photon with $p_{\mathrm{T}}>40~\mathrm{GeV}$, with all objects within $|\eta|=2.6$.

Table 47: Cutflow of expected events in the region SR-1J-High-EWK. Requirements below the line are specific to this region. On the Generator Filter line, the total number of unweighted events simulated is given in brackets. `Leptons' refers to electrons and muons only. For C1N2 models, the Generator Filter requires at least two $7~\mathrm{GeV}$ leptons and for C1N2 models with mass splittings below the Z boson mass it also requires $75~\mathrm{GeV}$ of $E_{\mathrm{T}}^{\mathrm{miss}}$. For GMSB models, the Generator Filter requires at least two $3~\mathrm{GeV}$ leptons. For on-shell C1N2 models, the `Forced Decays' require each Z boson to decay to a charged lepton pair (electron, muon, or tau) and each W boson to decay hadronically. For off-shell C1N2 models, each neutralino is forced to produce a charged lepton pair in its decay, and each chargino can produce any fermion pair. The SUSY2 kernel requires at least two leptons with $p_{\mathrm{T}}>9~\mathrm{GeV}$ or at least one lepton with $p_{\mathrm{T}}>25~\mathrm{GeV}$ and a photon with $p_{\mathrm{T}}>40~\mathrm{GeV}$, with all objects within $|\eta|=2.6$.

Table 47: Cutflow of expected events in the region SR-1J-High-EWK. Requirements below the line are specific to this region. On the Generator Filter line, the total number of unweighted events simulated is given in brackets. `Leptons' refers to electrons and muons only. For C1N2 models, the Generator Filter requires at least two $7~\mathrm{GeV}$ leptons and for C1N2 models with mass splittings below the Z boson mass it also requires $75~\mathrm{GeV}$ of $E_{\mathrm{T}}^{\mathrm{miss}}$. For GMSB models, the Generator Filter requires at least two $3~\mathrm{GeV}$ leptons. For on-shell C1N2 models, the `Forced Decays' require each Z boson to decay to a charged lepton pair (electron, muon, or tau) and each W boson to decay hadronically. For off-shell C1N2 models, each neutralino is forced to produce a charged lepton pair in its decay, and each chargino can produce any fermion pair. The SUSY2 kernel requires at least two leptons with $p_{\mathrm{T}}>9~\mathrm{GeV}$ or at least one lepton with $p_{\mathrm{T}}>25~\mathrm{GeV}$ and a photon with $p_{\mathrm{T}}>40~\mathrm{GeV}$, with all objects within $|\eta|=2.6$.

Table 48: Cutflow of expected events in the region SR-llbb-EWK. Requirements below the line are specific to this region. On the Generator Filter line, the total number of unweighted events simulated is given in brackets. `Leptons' refers to electrons and muons only. For C1N2 models, the Generator Filter requires at least two $7~\mathrm{GeV}$ leptons and for C1N2 models with mass splittings below the Z boson mass it also requires $75~\mathrm{GeV}$ of $E_{\mathrm{T}}^{\mathrm{miss}}$. For GMSB models, the Generator Filter requires at least two $3~\mathrm{GeV}$ leptons. For on-shell C1N2 models, the `Forced Decays' require each Z boson to decay to a charged lepton pair (electron, muon, or tau) and each W boson to decay hadronically. For off-shell C1N2 models, each neutralino is forced to produce a charged lepton pair in its decay, and each chargino can produce any fermion pair. The SUSY2 kernel requires at least two leptons with $p_{\mathrm{T}}>9~\mathrm{GeV}$ or at least one lepton with $p_{\mathrm{T}}>25~\mathrm{GeV}$ and a photon with $p_{\mathrm{T}}>40~\mathrm{GeV}$, with all objects within $|\eta|=2.6$.

Table 48: Cutflow of expected events in the region SR-llbb-EWK. Requirements below the line are specific to this region. On the Generator Filter line, the total number of unweighted events simulated is given in brackets. `Leptons' refers to electrons and muons only. For C1N2 models, the Generator Filter requires at least two $7~\mathrm{GeV}$ leptons and for C1N2 models with mass splittings below the Z boson mass it also requires $75~\mathrm{GeV}$ of $E_{\mathrm{T}}^{\mathrm{miss}}$. For GMSB models, the Generator Filter requires at least two $3~\mathrm{GeV}$ leptons. For on-shell C1N2 models, the `Forced Decays' require each Z boson to decay to a charged lepton pair (electron, muon, or tau) and each W boson to decay hadronically. For off-shell C1N2 models, each neutralino is forced to produce a charged lepton pair in its decay, and each chargino can produce any fermion pair. The SUSY2 kernel requires at least two leptons with $p_{\mathrm{T}}>9~\mathrm{GeV}$ or at least one lepton with $p_{\mathrm{T}}>25~\mathrm{GeV}$ and a photon with $p_{\mathrm{T}}>40~\mathrm{GeV}$, with all objects within $|\eta|=2.6$.

The combined $m_{jj}$ distribution of CR-Z-EWK and SR-Low-EWK (left), and the $\mathcal{S}(E_{\mathrm{T}}^{\mathrm{miss}})$ distribution in CR-Z-met-EWK (right), which removes the upper limit of $\mathcal{S}(E_{\mathrm{T}}^{\mathrm{miss}}) < 9$ from the definition of CR-Z-EWK. This $\mathcal{S}(E_{\mathrm{T}}^{\mathrm{miss}})$ tail overlaps other control and validation regions, but not signal regions. The arrows indicate the signal region SR-Low-EWK (left), and the $\mathcal{S}(E_{\mathrm{T}}^{\mathrm{miss}})$ phase space which is not included in CR-Z-EWK (right). All EWK search control and signal regions are included in the fit. All statistical and systematic uncertainties are included in the hatched bands. The theoretical uncertainties from CR-Z-EWK are applied to CR-Z-met-EWK. The last bins contain the overflow.

The combined $m_{jj}$ distribution of CR-Z-EWK and SR-Low-EWK (left), and the $\mathcal{S}(E_{\mathrm{T}}^{\mathrm{miss}})$ distribution in CR-Z-met-EWK (right), which removes the upper limit of $\mathcal{S}(E_{\mathrm{T}}^{\mathrm{miss}}) < 9$ from the definition of CR-Z-EWK. This $\mathcal{S}(E_{\mathrm{T}}^{\mathrm{miss}})$ tail overlaps other control and validation regions, but not signal regions. The arrows indicate the signal region SR-Low-EWK (left), and the $\mathcal{S}(E_{\mathrm{T}}^{\mathrm{miss}})$ phase space which is not included in CR-Z-EWK (right). All EWK search control and signal regions are included in the fit. All statistical and systematic uncertainties are included in the hatched bands. The theoretical uncertainties from CR-Z-EWK are applied to CR-Z-met-EWK. The last bins contain the overflow.

Signal region acceptance (left) and efficiency (right) for the GMSB model in the regions SR-OffShell-EWK and SR-Low-EWK. Acceptance is calculated by applying the signal-region requirements to particle-level objects, which do not suffer from identification inefficiencies or mismeasurements. For models with mass splittings below the Z boson mass, this filter also requires $E_{\mathrm{T}}^{\mathrm{miss}} > 75~\mathrm{GeV}$. The efficiency is calculated with fully reconstructed objects with the acceptance divided out.

Signal region acceptance (left) and efficiency (right) for the GMSB model in the regions SR-OffShell-EWK and SR-Low-EWK. Acceptance is calculated by applying the signal-region requirements to particle-level objects, which do not suffer from identification inefficiencies or mismeasurements. For models with mass splittings below the Z boson mass, this filter also requires $E_{\mathrm{T}}^{\mathrm{miss}} > 75~\mathrm{GeV}$. The efficiency is calculated with fully reconstructed objects with the acceptance divided out.

Signal region acceptance (left) and efficiency (right) for the GMSB model in the regions SR-OffShell-EWK and SR-Low-EWK. Acceptance is calculated by applying the signal-region requirements to particle-level objects, which do not suffer from identification inefficiencies or mismeasurements. For models with mass splittings below the Z boson mass, this filter also requires $E_{\mathrm{T}}^{\mathrm{miss}} > 75~\mathrm{GeV}$. The efficiency is calculated with fully reconstructed objects with the acceptance divided out.

Signal region acceptance (left) and efficiency (right) for the GMSB model in the regions SR-OffShell-EWK and SR-Low-EWK. Acceptance is calculated by applying the signal-region requirements to particle-level objects, which do not suffer from identification inefficiencies or mismeasurements. For models with mass splittings below the Z boson mass, this filter also requires $E_{\mathrm{T}}^{\mathrm{miss}} > 75~\mathrm{GeV}$. The efficiency is calculated with fully reconstructed objects with the acceptance divided out.

Signal region acceptance (left) and efficiency (right) for the GMSB model in the regions SR-OffShell-EWK and SR-Low-EWK. Acceptance is calculated by applying the signal-region requirements to particle-level objects, which do not suffer from identification inefficiencies or mismeasurements. For models with mass splittings below the Z boson mass, this filter also requires $E_{\mathrm{T}}^{\mathrm{miss}} > 75~\mathrm{GeV}$. The efficiency is calculated with fully reconstructed objects with the acceptance divided out.

Signal region acceptance (left) and efficiency (right) for the GMSB model in the regions SR-OffShell-EWK and SR-Low-EWK. Acceptance is calculated by applying the signal-region requirements to particle-level objects, which do not suffer from identification inefficiencies or mismeasurements. For models with mass splittings below the Z boson mass, this filter also requires $E_{\mathrm{T}}^{\mathrm{miss}} > 75~\mathrm{GeV}$. The efficiency is calculated with fully reconstructed objects with the acceptance divided out.

Signal region acceptance (left) and efficiency (right) for the GMSB model in the regions SR-OffShell-EWK and SR-Low-EWK. Acceptance is calculated by applying the signal-region requirements to particle-level objects, which do not suffer from identification inefficiencies or mismeasurements. For models with mass splittings below the Z boson mass, this filter also requires $E_{\mathrm{T}}^{\mathrm{miss}} > 75~\mathrm{GeV}$. The efficiency is calculated with fully reconstructed objects with the acceptance divided out.

Signal region acceptance (left) and efficiency (right) for the GMSB model in the regions SR-OffShell-EWK and SR-Low-EWK. Acceptance is calculated by applying the signal-region requirements to particle-level objects, which do not suffer from identification inefficiencies or mismeasurements. For models with mass splittings below the Z boson mass, this filter also requires $E_{\mathrm{T}}^{\mathrm{miss}} > 75~\mathrm{GeV}$. The efficiency is calculated with fully reconstructed objects with the acceptance divided out.

Signal region acceptance (left) and efficiency (right) for the GMSB model in the regions SR-Low-2-EWK and SR-Int-EWK. Acceptance is calculated by applying the signal-region requirements to particle-level objects, which do not suffer from identification inefficiencies or mismeasurements. For models with mass splittings below the Z boson mass, this filter also requires $E_{\mathrm{T}}^{\mathrm{miss}} > 75~\mathrm{GeV}$. The efficiency is calculated with fully reconstructed objects with the acceptance divided out.

Signal region acceptance (left) and efficiency (right) for the GMSB model in the regions SR-Low-2-EWK and SR-Int-EWK. Acceptance is calculated by applying the signal-region requirements to particle-level objects, which do not suffer from identification inefficiencies or mismeasurements. For models with mass splittings below the Z boson mass, this filter also requires $E_{\mathrm{T}}^{\mathrm{miss}} > 75~\mathrm{GeV}$. The efficiency is calculated with fully reconstructed objects with the acceptance divided out.

Signal region acceptance (left) and efficiency (right) for the GMSB model in the regions SR-Low-2-EWK and SR-Int-EWK. Acceptance is calculated by applying the signal-region requirements to particle-level objects, which do not suffer from identification inefficiencies or mismeasurements. For models with mass splittings below the Z boson mass, this filter also requires $E_{\mathrm{T}}^{\mathrm{miss}} > 75~\mathrm{GeV}$. The efficiency is calculated with fully reconstructed objects with the acceptance divided out.

Signal region acceptance (left) and efficiency (right) for the GMSB model in the regions SR-Low-2-EWK and SR-Int-EWK. Acceptance is calculated by applying the signal-region requirements to particle-level objects, which do not suffer from identification inefficiencies or mismeasurements. For models with mass splittings below the Z boson mass, this filter also requires $E_{\mathrm{T}}^{\mathrm{miss}} > 75~\mathrm{GeV}$. The efficiency is calculated with fully reconstructed objects with the acceptance divided out.

Signal region acceptance (left) and efficiency (right) for the GMSB model in the regions SR-Low-2-EWK and SR-Int-EWK. Acceptance is calculated by applying the signal-region requirements to particle-level objects, which do not suffer from identification inefficiencies or mismeasurements. For models with mass splittings below the Z boson mass, this filter also requires $E_{\mathrm{T}}^{\mathrm{miss}} > 75~\mathrm{GeV}$. The efficiency is calculated with fully reconstructed objects with the acceptance divided out.

Signal region acceptance (left) and efficiency (right) for the GMSB model in the regions SR-Low-2-EWK and SR-Int-EWK. Acceptance is calculated by applying the signal-region requirements to particle-level objects, which do not suffer from identification inefficiencies or mismeasurements. For models with mass splittings below the Z boson mass, this filter also requires $E_{\mathrm{T}}^{\mathrm{miss}} > 75~\mathrm{GeV}$. The efficiency is calculated with fully reconstructed objects with the acceptance divided out.

Signal region acceptance (left) and efficiency (right) for the GMSB model in the regions SR-High-EWK. Acceptance is calculated by applying the signal-region requirements to particle-level objects, which do not suffer from identification inefficiencies or mismeasurements. For models with mass splittings below the Z boson mass, this filter also requires $E_{\mathrm{T}}^{\mathrm{miss}} > 75~\mathrm{GeV}$. The efficiency is calculated with fully reconstructed objects with the acceptance divided out.

Signal region acceptance (left) and efficiency (right) for the GMSB model in the regions SR-High-EWK. Acceptance is calculated by applying the signal-region requirements to particle-level objects, which do not suffer from identification inefficiencies or mismeasurements. For models with mass splittings below the Z boson mass, this filter also requires $E_{\mathrm{T}}^{\mathrm{miss}} > 75~\mathrm{GeV}$. The efficiency is calculated with fully reconstructed objects with the acceptance divided out.

Signal region acceptance (left) and efficiency (right) for the GMSB model in the regions SR-High-EWK. Acceptance is calculated by applying the signal-region requirements to particle-level objects, which do not suffer from identification inefficiencies or mismeasurements. For models with mass splittings below the Z boson mass, this filter also requires $E_{\mathrm{T}}^{\mathrm{miss}} > 75~\mathrm{GeV}$. The efficiency is calculated with fully reconstructed objects with the acceptance divided out.

Signal region acceptance (left) and efficiency (right) for the GMSB model in the regions SR-High-EWK. Acceptance is calculated by applying the signal-region requirements to particle-level objects, which do not suffer from identification inefficiencies or mismeasurements. For models with mass splittings below the Z boson mass, this filter also requires $E_{\mathrm{T}}^{\mathrm{miss}} > 75~\mathrm{GeV}$. The efficiency is calculated with fully reconstructed objects with the acceptance divided out.

Signal region acceptance (left) and efficiency (right) for the GMSB model in the regions SR-High-EWK. Acceptance is calculated by applying the signal-region requirements to particle-level objects, which do not suffer from identification inefficiencies or mismeasurements. For models with mass splittings below the Z boson mass, this filter also requires $E_{\mathrm{T}}^{\mathrm{miss}} > 75~\mathrm{GeV}$. The efficiency is calculated with fully reconstructed objects with the acceptance divided out.

Signal region acceptance (left) and efficiency (right) for the GMSB model in the regions SR-High-EWK. Acceptance is calculated by applying the signal-region requirements to particle-level objects, which do not suffer from identification inefficiencies or mismeasurements. For models with mass splittings below the Z boson mass, this filter also requires $E_{\mathrm{T}}^{\mathrm{miss}} > 75~\mathrm{GeV}$. The efficiency is calculated with fully reconstructed objects with the acceptance divided out.

Signal region acceptance (left) and efficiency (right) for the GMSB model in the regions SR-High-EWK. Acceptance is calculated by applying the signal-region requirements to particle-level objects, which do not suffer from identification inefficiencies or mismeasurements. For models with mass splittings below the Z boson mass, this filter also requires $E_{\mathrm{T}}^{\mathrm{miss}} > 75~\mathrm{GeV}$. The efficiency is calculated with fully reconstructed objects with the acceptance divided out.

Signal region acceptance (left) and efficiency (right) for the GMSB model in the regions SR-High-EWK. Acceptance is calculated by applying the signal-region requirements to particle-level objects, which do not suffer from identification inefficiencies or mismeasurements. For models with mass splittings below the Z boson mass, this filter also requires $E_{\mathrm{T}}^{\mathrm{miss}} > 75~\mathrm{GeV}$. The efficiency is calculated with fully reconstructed objects with the acceptance divided out.

Signal region acceptance (left) and efficiency (right) for the GMSB model in the regions SR-1J-High-EWK and SR-$\ell\ell bb$-EWK. Acceptance is calculated by applying the signal-region requirements to particle-level objects, which do not suffer from identification inefficiencies or mismeasurements. For models with mass splittings below the Z boson mass, this filter also requires $E_{\mathrm{T}}^{\mathrm{miss}} > 75~\mathrm{GeV}$. The efficiency is calculated with fully reconstructed objects with the acceptance divided out.

Signal region acceptance (left) and efficiency (right) for the GMSB model in the regions SR-1J-High-EWK and SR-$\ell\ell bb$-EWK. Acceptance is calculated by applying the signal-region requirements to particle-level objects, which do not suffer from identification inefficiencies or mismeasurements. For models with mass splittings below the Z boson mass, this filter also requires $E_{\mathrm{T}}^{\mathrm{miss}} > 75~\mathrm{GeV}$. The efficiency is calculated with fully reconstructed objects with the acceptance divided out.

Signal region acceptance (left) and efficiency (right) for the GMSB model in the regions SR-1J-High-EWK and SR-$\ell\ell bb$-EWK. Acceptance is calculated by applying the signal-region requirements to particle-level objects, which do not suffer from identification inefficiencies or mismeasurements. For models with mass splittings below the Z boson mass, this filter also requires $E_{\mathrm{T}}^{\mathrm{miss}} > 75~\mathrm{GeV}$. The efficiency is calculated with fully reconstructed objects with the acceptance divided out.

Signal region acceptance (left) and efficiency (right) for the GMSB model in the regions SR-1J-High-EWK and SR-$\ell\ell bb$-EWK. Acceptance is calculated by applying the signal-region requirements to particle-level objects, which do not suffer from identification inefficiencies or mismeasurements. For models with mass splittings below the Z boson mass, this filter also requires $E_{\mathrm{T}}^{\mathrm{miss}} > 75~\mathrm{GeV}$. The efficiency is calculated with fully reconstructed objects with the acceptance divided out.

Signal region acceptance (left) and efficiency (right) for the C1N2 model in the regions SR-OffShell-EWK and SR-Low-EWK. Acceptance is calculated by applying the signal-region requirements to particle-level objects, which do not suffer from identification inefficiencies or mismeasurements. For models with mass splittings below the Z boson mass, this filter also requires $E_{\mathrm{T}}^{\mathrm{miss}} > 75~\mathrm{GeV}$. The efficiency is calculated with fully reconstructed objects with the acceptance divided out.

Signal region acceptance (left) and efficiency (right) for the C1N2 model in the regions SR-OffShell-EWK and SR-Low-EWK. Acceptance is calculated by applying the signal-region requirements to particle-level objects, which do not suffer from identification inefficiencies or mismeasurements. For models with mass splittings below the Z boson mass, this filter also requires $E_{\mathrm{T}}^{\mathrm{miss}} > 75~\mathrm{GeV}$. The efficiency is calculated with fully reconstructed objects with the acceptance divided out.

Signal region acceptance (left) and efficiency (right) for the C1N2 model in the regions SR-OffShell-EWK and SR-Low-EWK. Acceptance is calculated by applying the signal-region requirements to particle-level objects, which do not suffer from identification inefficiencies or mismeasurements. For models with mass splittings below the Z boson mass, this filter also requires $E_{\mathrm{T}}^{\mathrm{miss}} > 75~\mathrm{GeV}$. The efficiency is calculated with fully reconstructed objects with the acceptance divided out.

Signal region acceptance (left) and efficiency (right) for the C1N2 model in the regions SR-OffShell-EWK and SR-Low-EWK. Acceptance is calculated by applying the signal-region requirements to particle-level objects, which do not suffer from identification inefficiencies or mismeasurements. For models with mass splittings below the Z boson mass, this filter also requires $E_{\mathrm{T}}^{\mathrm{miss}} > 75~\mathrm{GeV}$. The efficiency is calculated with fully reconstructed objects with the acceptance divided out.

Signal region acceptance (left) and efficiency (right) for the C1N2 model in the regions SR-OffShell-EWK and SR-Low-EWK. Acceptance is calculated by applying the signal-region requirements to particle-level objects, which do not suffer from identification inefficiencies or mismeasurements. For models with mass splittings below the Z boson mass, this filter also requires $E_{\mathrm{T}}^{\mathrm{miss}} > 75~\mathrm{GeV}$. The efficiency is calculated with fully reconstructed objects with the acceptance divided out.

Signal region acceptance (left) and efficiency (right) for the C1N2 model in the regions SR-OffShell-EWK and SR-Low-EWK. Acceptance is calculated by applying the signal-region requirements to particle-level objects, which do not suffer from identification inefficiencies or mismeasurements. For models with mass splittings below the Z boson mass, this filter also requires $E_{\mathrm{T}}^{\mathrm{miss}} > 75~\mathrm{GeV}$. The efficiency is calculated with fully reconstructed objects with the acceptance divided out.

Signal region acceptance (left) and efficiency (right) for the C1N2 model in the regions SR-OffShell-EWK and SR-Low-EWK. Acceptance is calculated by applying the signal-region requirements to particle-level objects, which do not suffer from identification inefficiencies or mismeasurements. For models with mass splittings below the Z boson mass, this filter also requires $E_{\mathrm{T}}^{\mathrm{miss}} > 75~\mathrm{GeV}$. The efficiency is calculated with fully reconstructed objects with the acceptance divided out.

Signal region acceptance (left) and efficiency (right) for the C1N2 model in the regions SR-OffShell-EWK and SR-Low-EWK. Acceptance is calculated by applying the signal-region requirements to particle-level objects, which do not suffer from identification inefficiencies or mismeasurements. For models with mass splittings below the Z boson mass, this filter also requires $E_{\mathrm{T}}^{\mathrm{miss}} > 75~\mathrm{GeV}$. The efficiency is calculated with fully reconstructed objects with the acceptance divided out.

Signal region acceptance (left) and efficiency (right) for the C1N2 model in the regions SR-Low-2-EWK and SR-Int-EWK. Acceptance is calculated by applying the signal-region requirements to particle-level objects, which do not suffer from identification inefficiencies or mismeasurements. For models with mass splittings below the Z boson mass, this filter also requires $E_{\mathrm{T}}^{\mathrm{miss}} > 75~\mathrm{GeV}$. The efficiency is calculated with fully reconstructed objects with the acceptance divided out.

Signal region acceptance (left) and efficiency (right) for the C1N2 model in the regions SR-Low-2-EWK and SR-Int-EWK. Acceptance is calculated by applying the signal-region requirements to particle-level objects, which do not suffer from identification inefficiencies or mismeasurements. For models with mass splittings below the Z boson mass, this filter also requires $E_{\mathrm{T}}^{\mathrm{miss}} > 75~\mathrm{GeV}$. The efficiency is calculated with fully reconstructed objects with the acceptance divided out.

Signal region acceptance (left) and efficiency (right) for the C1N2 model in the regions SR-Low-2-EWK and SR-Int-EWK. Acceptance is calculated by applying the signal-region requirements to particle-level objects, which do not suffer from identification inefficiencies or mismeasurements. For models with mass splittings below the Z boson mass, this filter also requires $E_{\mathrm{T}}^{\mathrm{miss}} > 75~\mathrm{GeV}$. The efficiency is calculated with fully reconstructed objects with the acceptance divided out.

Signal region acceptance (left) and efficiency (right) for the C1N2 model in the regions SR-Low-2-EWK and SR-Int-EWK. Acceptance is calculated by applying the signal-region requirements to particle-level objects, which do not suffer from identification inefficiencies or mismeasurements. For models with mass splittings below the Z boson mass, this filter also requires $E_{\mathrm{T}}^{\mathrm{miss}} > 75~\mathrm{GeV}$. The efficiency is calculated with fully reconstructed objects with the acceptance divided out.

Signal region acceptance (left) and efficiency (right) for the C1N2 model in the regions SR-Low-2-EWK and SR-Int-EWK. Acceptance is calculated by applying the signal-region requirements to particle-level objects, which do not suffer from identification inefficiencies or mismeasurements. For models with mass splittings below the Z boson mass, this filter also requires $E_{\mathrm{T}}^{\mathrm{miss}} > 75~\mathrm{GeV}$. The efficiency is calculated with fully reconstructed objects with the acceptance divided out.

Signal region acceptance (left) and efficiency (right) for the C1N2 model in the regions SR-Low-2-EWK and SR-Int-EWK. Acceptance is calculated by applying the signal-region requirements to particle-level objects, which do not suffer from identification inefficiencies or mismeasurements. For models with mass splittings below the Z boson mass, this filter also requires $E_{\mathrm{T}}^{\mathrm{miss}} > 75~\mathrm{GeV}$. The efficiency is calculated with fully reconstructed objects with the acceptance divided out.

Signal region acceptance (left) and efficiency (right) for the C1N2 model in the regions SR-High-EWK. Acceptance is calculated by applying the signal-region requirements to particle-level objects, which do not suffer from identification inefficiencies or mismeasurements. For models with mass splittings below the Z boson mass, this filter also requires $E_{\mathrm{T}}^{\mathrm{miss}} > 75~\mathrm{GeV}$. The efficiency is calculated with fully reconstructed objects with the acceptance divided out.

Signal region acceptance (left) and efficiency (right) for the C1N2 model in the regions SR-High-EWK. Acceptance is calculated by applying the signal-region requirements to particle-level objects, which do not suffer from identification inefficiencies or mismeasurements. For models with mass splittings below the Z boson mass, this filter also requires $E_{\mathrm{T}}^{\mathrm{miss}} > 75~\mathrm{GeV}$. The efficiency is calculated with fully reconstructed objects with the acceptance divided out.

Signal region acceptance (left) and efficiency (right) for the C1N2 model in the regions SR-High-EWK. Acceptance is calculated by applying the signal-region requirements to particle-level objects, which do not suffer from identification inefficiencies or mismeasurements. For models with mass splittings below the Z boson mass, this filter also requires $E_{\mathrm{T}}^{\mathrm{miss}} > 75~\mathrm{GeV}$. The efficiency is calculated with fully reconstructed objects with the acceptance divided out.

Signal region acceptance (left) and efficiency (right) for the C1N2 model in the regions SR-High-EWK. Acceptance is calculated by applying the signal-region requirements to particle-level objects, which do not suffer from identification inefficiencies or mismeasurements. For models with mass splittings below the Z boson mass, this filter also requires $E_{\mathrm{T}}^{\mathrm{miss}} > 75~\mathrm{GeV}$. The efficiency is calculated with fully reconstructed objects with the acceptance divided out.

Signal region acceptance (left) and efficiency (right) for the C1N2 model in the regions SR-High-EWK. Acceptance is calculated by applying the signal-region requirements to particle-level objects, which do not suffer from identification inefficiencies or mismeasurements. For models with mass splittings below the Z boson mass, this filter also requires $E_{\mathrm{T}}^{\mathrm{miss}} > 75~\mathrm{GeV}$. The efficiency is calculated with fully reconstructed objects with the acceptance divided out.

Signal region acceptance (left) and efficiency (right) for the C1N2 model in the regions SR-High-EWK. Acceptance is calculated by applying the signal-region requirements to particle-level objects, which do not suffer from identification inefficiencies or mismeasurements. For models with mass splittings below the Z boson mass, this filter also requires $E_{\mathrm{T}}^{\mathrm{miss}} > 75~\mathrm{GeV}$. The efficiency is calculated with fully reconstructed objects with the acceptance divided out.

Signal region acceptance (left) and efficiency (right) for the C1N2 model in the regions SR-High-EWK. Acceptance is calculated by applying the signal-region requirements to particle-level objects, which do not suffer from identification inefficiencies or mismeasurements. For models with mass splittings below the Z boson mass, this filter also requires $E_{\mathrm{T}}^{\mathrm{miss}} > 75~\mathrm{GeV}$. The efficiency is calculated with fully reconstructed objects with the acceptance divided out.

Signal region acceptance (left) and efficiency (right) for the C1N2 model in the regions SR-High-EWK. Acceptance is calculated by applying the signal-region requirements to particle-level objects, which do not suffer from identification inefficiencies or mismeasurements. For models with mass splittings below the Z boson mass, this filter also requires $E_{\mathrm{T}}^{\mathrm{miss}} > 75~\mathrm{GeV}$. The efficiency is calculated with fully reconstructed objects with the acceptance divided out.

Signal region acceptance (left) and efficiency (right) for the C1N2 model in the regions SR-1J-High-EWK and SR-$\ell\ell bb$-EWK. Acceptance is calculated by applying the signal-region requirements to particle-level objects, which do not suffer from identification inefficiencies or mismeasurements. For models with mass splittings below the Z boson mass, this filter also requires $E_{\mathrm{T}}^{\mathrm{miss}} > 75~\mathrm{GeV}$. The efficiency is calculated with fully reconstructed objects with the acceptance divided out.

Signal region acceptance (left) and efficiency (right) for the C1N2 model in the regions SR-1J-High-EWK and SR-$\ell\ell bb$-EWK. Acceptance is calculated by applying the signal-region requirements to particle-level objects, which do not suffer from identification inefficiencies or mismeasurements. For models with mass splittings below the Z boson mass, this filter also requires $E_{\mathrm{T}}^{\mathrm{miss}} > 75~\mathrm{GeV}$. The efficiency is calculated with fully reconstructed objects with the acceptance divided out.

Signal region acceptance (left) and efficiency (right) for the C1N2 model in the regions SR-1J-High-EWK and SR-$\ell\ell bb$-EWK. Acceptance is calculated by applying the signal-region requirements to particle-level objects, which do not suffer from identification inefficiencies or mismeasurements. For models with mass splittings below the Z boson mass, this filter also requires $E_{\mathrm{T}}^{\mathrm{miss}} > 75~\mathrm{GeV}$. The efficiency is calculated with fully reconstructed objects with the acceptance divided out.

Signal region acceptance (left) and efficiency (right) for the C1N2 model in the regions SR-1J-High-EWK and SR-$\ell\ell bb$-EWK. Acceptance is calculated by applying the signal-region requirements to particle-level objects, which do not suffer from identification inefficiencies or mismeasurements. For models with mass splittings below the Z boson mass, this filter also requires $E_{\mathrm{T}}^{\mathrm{miss}} > 75~\mathrm{GeV}$. The efficiency is calculated with fully reconstructed objects with the acceptance divided out.

Signal region acceptance (left) and efficiency (right) over the full \mll\ range for the GG_N2_SLN1 model in Strong search regions. 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.

Signal region acceptance (left) and efficiency (right) over the full \mll\ range for the GG_N2_SLN1 model in Strong search regions. 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.

Signal region acceptance (left) and efficiency (right) over the full \mll\ range for the GG_N2_ZN1 model in Strong search regions. 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.

Signal region acceptance (left) and efficiency (right) over the full \mll\ range for the GG_N2_ZN1 model in Strong search regions. 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.

Signal region acceptance (left) and efficiency (right) over the full \mll\ range for the SS_N2_ZN1 model in Strong search regions. 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.

Signal region acceptance (left) and efficiency (right) over the full \mll\ range for the SS_N2_ZN1 model in Strong search regions. 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.

Signal region acceptance (left) and efficiency (right) over the full \mll\ range for the GG_N2_SLN1 model in Strong search regions. 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.

Signal region acceptance (left) and efficiency (right) over the full \mll\ range for the GG_N2_SLN1 model in Strong search regions. 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.

Signal region acceptance (left) and efficiency (right) over the full \mll\ range for the GG_N2_ZN1 model in Strong search regions. 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.

Signal region acceptance (left) and efficiency (right) over the full \mll\ range for the GG_N2_ZN1 model in Strong search regions. 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.

Signal region acceptance (left) and efficiency (right) over the full \mll\ range for the SS_N2_ZN1 model in Strong search regions. 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.

Signal region acceptance (left) and efficiency (right) over the full \mll\ range for the SS_N2_ZN1 model in Strong search regions. 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.

Signal region acceptance (left) and efficiency (right) over the full \mll\ range for the GG_N2_SLN1 model in Strong search regions. 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.

Signal region acceptance (left) and efficiency (right) over the full \mll\ range for the GG_N2_SLN1 model in Strong search regions. 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.

Signal region acceptance (left) and efficiency (right) over the full \mll\ range for the GG_N2_ZN1 model in Strong search regions. 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.

Signal region acceptance (left) and efficiency (right) over the full \mll\ range for the GG_N2_ZN1 model in Strong search regions. 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.

Signal region acceptance (left) and efficiency (right) over the full \mll\ range for the SS_N2_ZN1 model in Strong search regions. 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.

Signal region acceptance (left) and efficiency (right) over the full \mll\ range for the SS_N2_ZN1 model in Strong search regions. 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.

Signal region acceptance (left) and efficiency (right) over the full \mll\ range for the GG_N2_SLN1 model in Strong search regions. 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.

Signal region acceptance (left) and efficiency (right) over the full \mll\ range for the GG_N2_SLN1 model in Strong search regions. 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.

Signal region acceptance (left) and efficiency (right) over the full \mll\ range for the GG_N2_ZN1 model in Strong search regions. 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.

Signal region acceptance (left) and efficiency (right) over the full \mll\ range for the GG_N2_ZN1 model in Strong search regions. 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.

Signal region acceptance (left) and efficiency (right) over the full \mll\ range for the SS_N2_ZN1 model in Strong search regions. 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.

Signal region acceptance (left) and efficiency (right) over the full \mll\ range for the SS_N2_ZN1 model in Strong search regions. 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.

Signal region acceptance (left) and efficiency (right) over the full \mll\ range for the GG_N2_ZN1 model in Strong search regions. 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.

Signal region acceptance (left) and efficiency (right) over the full \mll\ range for the GG_N2_ZN1 model in Strong search regions. 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.

Signal region acceptance (left) and efficiency (right) over the full \mll\ range for the SS_N2_ZN1 model in Strong search regions. 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.

Signal region acceptance (left) and efficiency (right) over the full \mll\ range for the SS_N2_ZN1 model in Strong search regions. 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.

Signal region acceptance (left) and efficiency (right) over the full \mll\ range for the GG_N2_ZN1 model in Strong search regions. 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.

Signal region acceptance (left) and efficiency (right) over the full \mll\ range for the GG_N2_ZN1 model in Strong search regions. 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.

Signal region acceptance (left) and efficiency (right) over the full \mll\ range for the SS_N2_ZN1 model in Strong search regions. 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.

Signal region acceptance (left) and efficiency (right) over the full \mll\ range for the SS_N2_ZN1 model in Strong search regions. 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.

Signal region acceptance (left) and efficiency (right) over the full \mll\ range for the GG_N2_ZN1 model in Strong search regions. 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.

Signal region acceptance (left) and efficiency (right) over the full \mll\ range for the GG_N2_ZN1 model in Strong search regions. 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.

Signal region acceptance (left) and efficiency (right) over the full \mll\ range for the SS_N2_ZN1 model in Strong search regions. 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.

Signal region acceptance (left) and efficiency (right) over the full \mll\ range for the SS_N2_ZN1 model in Strong search regions. 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.

The observed mass distribution in the SR-Inclusive_High signal-region bin. The band on the expected background indicates the total uncertainty of the estimation. Several representative signal models are overlaid. Events outside the shown range are accumulated in the rightmost bin indicated as 'Overflow'. Downward triangle markers at the bottom of the panels indicate that no events are observed in the corresponding mass bin, while upward triangle markers in the lower panels indicate that the observed data is beyond the range.

List of expected and observed events, $p_{0}$-value and the corresponding $Z$ local significance, as well as the 95% CLs upper limit of the expected and observed signal events ($S^{95}_ ext{exp} and $S^{95}_ ext{obs}$ ) in each mass window for SR-Inclusive bins of the short lifetime regime.

List of expected and observed events, $p_{0}$-value and the corresponding $Z$ local significance, as well as the 95% CLs upper limit of the expected and observed signal events ($S^{95}_ ext{exp} and $S^{95}_ ext{obs}$ ) in each mass window for SR-Inclusive bins of the long lifetime regime.

The observed $p_{\rm T$ distribution in the SR-Inclusive_Low signal-region bin. The band on the expected background indicates the total uncertainty of the estimation. Several representative signal models are overlaid. Events outside the shown range are accumulated in the rightmost bin indicated as 'Overflow'. Downward triangle markers at the bottom of the panels indicate that no events are observed in the corresponding mass bin, while upward triangle markers in the lower panels indicate that the observed data is beyond the range.

The observed $p_{\rm T$ distribution in the SR-Inclusive_High signal-region bin. The band on the expected background indicates the total uncertainty of the estimation. Several representative signal models are overlaid. Events outside the shown range are accumulated in the rightmost bin indicated as 'Overflow'. Downward triangle markers at the bottom of the panels indicate that no events are observed in the corresponding mass bin, while upward triangle markers in the lower panels indicate that the observed data is beyond the range.

The observed $|\eta|$ distribution in the SR-Inclusive_Low signal-region bin. The band on the expected background indicates the total uncertainty of the estimation. Several representative signal models are overlaid. Events outside the shown range are accumulated in the rightmost bin indicated as 'Overflow'. Downward triangle markers at the bottom of the panels indicate that no events are observed in the corresponding mass bin, while upward triangle markers in the lower panels indicate that the observed data is beyond the range.

The observed $|\eta|$ distribution in the SR-Inclusive_High signal-region bin. The band on the expected background indicates the total uncertainty of the estimation. Several representative signal models are overlaid. Events outside the shown range are accumulated in the rightmost bin indicated as 'Overflow'. Downward triangle markers at the bottom of the panels indicate that no events are observed in the corresponding mass bin, while upward triangle markers in the lower panels indicate that the observed data is beyond the range.

The observed dE/dx distribution in the SR-Inclusive_Low signal-region bin. The band on the expected background indicates the total uncertainty of the estimation. Several representative signal models are overlaid. Events outside the shown range are accumulated in the rightmost bin indicated as 'Overflow'. Downward triangle markers at the bottom of the panels indicate that no events are observed in the corresponding mass bin, while upward triangle markers in the lower panels indicate that the observed data is beyond the range.

The observed dE/dx distribution in the SR-Inclusive_High signal-region bin. The band on the expected background indicates the total uncertainty of the estimation. Several representative signal models are overlaid. Events outside the shown range are accumulated in the rightmost bin indicated as 'Overflow'. Downward triangle markers at the bottom of the panels indicate that no events are observed in the corresponding mass bin, while upward triangle markers in the lower panels indicate that the observed data is beyond the range.

The observed mass distribution in the SR-Trk-IBL0_Low signal-region bin. The band on the expected background indicates the total uncertainty of the estimation. Several representative signal models are overlaid. Events outside the shown range are accumulated in the rightmost bin indicated as 'Overflow'. Downward triangle markers at the bottom of the panels indicate that no events are observed in the corresponding mass bin, while upward triangle markers in the lower panels indicate that the observed data is beyond the range.

The observed mass distribution in the SR-Mu-IBL0_Low signal-region bin. The band on the expected background indicates the total uncertainty of the estimation. Several representative signal models are overlaid. Events outside the shown range are accumulated in the rightmost bin indicated as 'Overflow'. Downward triangle markers at the bottom of the panels indicate that no events are observed in the corresponding mass bin, while upward triangle markers in the lower panels indicate that the observed data is beyond the range.

The observed mass distribution in the SR-Trk-IBL0_High signal-region bin. The band on the expected background indicates the total uncertainty of the estimation. Several representative signal models are overlaid. Events outside the shown range are accumulated in the rightmost bin indicated as 'Overflow'. Downward triangle markers at the bottom of the panels indicate that no events are observed in the corresponding mass bin, while upward triangle markers in the lower panels indicate that the observed data is beyond the range.

The observed mass distribution in the SR-Mu-IBL0_High signal-region bin. The band on the expected background indicates the total uncertainty of the estimation. Several representative signal models are overlaid. Events outside the shown range are accumulated in the rightmost bin indicated as 'Overflow'. Downward triangle markers at the bottom of the panels indicate that no events are observed in the corresponding mass bin, while upward triangle markers in the lower panels indicate that the observed data is beyond the range.

The observed mass distribution in the SR-Trk-IBL1 signal-region bin. The band on the expected background indicates the total uncertainty of the estimation. Several representative signal models are overlaid. Events outside the shown range are accumulated in the rightmost bin indicated as 'Overflow'. Downward triangle markers at the bottom of the panels indicate that no events are observed in the corresponding mass bin, while upward triangle markers in the lower panels indicate that the observed data is beyond the range.

The observed mass distribution in the SR-Mu-IBL1 signal-region bin. The band on the expected background indicates the total uncertainty of the estimation. Several representative signal models are overlaid. Events outside the shown range are accumulated in the rightmost bin indicated as 'Overflow'. Downward triangle markers at the bottom of the panels indicate that no events are observed in the corresponding mass bin, while upward triangle markers in the lower panels indicate that the observed data is beyond the range.

Lower limits on the gluino mass, from gluino $R$-hadron pair production, as a function of gluino lifetime for two neutralino mass assumptions of (a) $m(\tilde{\chi}_{1}^{0}) = 100 \text{GeV}$ and (b) $\Delta m(\tilde{g}, \tilde{\chi}_{1}^{0}) = 30 \text{GeV}$. The upper $1 \sigma_\text{exp}$ expected bound is very close to the expected limit for some lifetime values due to the expected background getting very close to 0 events.

Lower limits on the gluino mass, from gluino $R$-hadron pair production, as a function of gluino lifetime for two neutralino mass assumptions of (a) $m(\tilde{\chi}_{1}^{0}) = 100 \text{GeV}$ and (b) $\Delta m(\tilde{g}, \tilde{\chi}_{1}^{0}) = 30 \text{GeV}$. The upper $1 \sigma_\text{exp}$ expected bound is very close to the expected limit for some lifetime values due to the expected background getting very close to 0 events.

Lower limits on the gluino mass, from gluino $R$-hadron pair production, as a function of gluino lifetime for two neutralino mass assumptions of (a) $m(\tilde{\chi}_{1}^{0}) = 100 \text{GeV}$ and (b) $\Delta m(\tilde{g}, \tilde{\chi}_{1}^{0}) = 30 \text{GeV}$. The upper $1 \sigma_\text{exp}$ expected bound is very close to the expected limit for some lifetime values due to the expected background getting very close to 0 events.

Lower limits on the gluino mass, from gluino $R$-hadron pair production, as a function of gluino lifetime for two neutralino mass assumptions of (a) $m(\tilde{\chi}_{1}^{0}) = 100 \text{GeV}$ and (b) $\Delta m(\tilde{g}, \tilde{\chi}_{1}^{0}) = 30 \text{GeV}$. The upper $1 \sigma_\text{exp}$ expected bound is very close to the expected limit for some lifetime values due to the expected background getting very close to 0 events.

(a) Lower limits on the chargino mass as a function of lifetime, and (b) the contours around the excluded mass-lifetime region for stau pair production.

(a) Lower limits on the chargino mass as a function of lifetime, and (b) the contours around the excluded mass-lifetime region for stau pair production.

(a) Lower limits on the chargino mass as a function of lifetime, and (b) the contours around the excluded mass-lifetime region for stau pair production.

(a) Lower limits on the chargino mass as a function of lifetime, and (b) the contours around the excluded mass-lifetime region for stau pair production.

Comparison of the observed and expected VAR distributionsin VR-LowPt-Trk-IBL0_Low. The band on the expected background estimation indicates the total uncertainty of the estimation. Downward triangle markers at the bottom of the panels indicate there is no events observed in the corresponding bin, while upward triangle markers at the bottom panel indicate the observed data is beyond the range.

Comparison of the observed and expected VAR distributionsin VR-LowPt-Mu-IBL0_Low. The band on the expected background estimation indicates the total uncertainty of the estimation. Downward triangle markers at the bottom of the panels indicate there is no events observed in the corresponding bin, while upward triangle markers at the bottom panel indicate the observed data is beyond the range.

Comparison of the observed and expected VAR distributionsin VR-LowPt-Trk-IBL0_High. The band on the expected background estimation indicates the total uncertainty of the estimation. Downward triangle markers at the bottom of the panels indicate there is no events observed in the corresponding bin, while upward triangle markers at the bottom panel indicate the observed data is beyond the range.

Comparison of the observed and expected VAR distributionsin VR-LowPt-Mu-IBL0_High. The band on the expected background estimation indicates the total uncertainty of the estimation. Downward triangle markers at the bottom of the panels indicate there is no events observed in the corresponding bin, while upward triangle markers at the bottom panel indicate the observed data is beyond the range.

Comparison of the observed and expected VAR distributionsin VR-LowPt-Trk-IBL1. The band on the expected background estimation indicates the total uncertainty of the estimation. Downward triangle markers at the bottom of the panels indicate there is no events observed in the corresponding bin, while upward triangle markers at the bottom panel indicate the observed data is beyond the range.

Comparison of the observed and expected VAR distributionsin VR-LowPt-Mu-IBL1. The band on the expected background estimation indicates the total uncertainty of the estimation. Downward triangle markers at the bottom of the panels indicate there is no events observed in the corresponding bin, while upward triangle markers at the bottom panel indicate the observed data is beyond the range.

Comparison of the observed and expected VAR distributionsin VR-HiEta-Trk-IBL0_Low. The band on the expected background estimation indicates the total uncertainty of the estimation. Downward triangle markers at the bottom of the panels indicate there is no events observed in the corresponding bin, while upward triangle markers at the bottom panel indicate the observed data is beyond the range.

Comparison of the observed and expected VAR distributionsin VR-HiEta-Mu-IBL0_Low. The band on the expected background estimation indicates the total uncertainty of the estimation. Downward triangle markers at the bottom of the panels indicate there is no events observed in the corresponding bin, while upward triangle markers at the bottom panel indicate the observed data is beyond the range.

Comparison of the observed and expected VAR distributionsin VR-HiEta-Trk-IBL0_High. The band on the expected background estimation indicates the total uncertainty of the estimation. Downward triangle markers at the bottom of the panels indicate there is no events observed in the corresponding bin, while upward triangle markers at the bottom panel indicate the observed data is beyond the range.

Comparison of the observed and expected VAR distributionsin VR-HiEta-Mu-IBL0_High. The band on the expected background estimation indicates the total uncertainty of the estimation. Downward triangle markers at the bottom of the panels indicate there is no events observed in the corresponding bin, while upward triangle markers at the bottom panel indicate the observed data is beyond the range.

Comparison of the observed and expected VAR distributionsin VR-HiEta-Trk-IBL1. The band on the expected background estimation indicates the total uncertainty of the estimation. Downward triangle markers at the bottom of the panels indicate there is no events observed in the corresponding bin, while upward triangle markers at the bottom panel indicate the observed data is beyond the range.

Comparison of the observed and expected VAR distributionsin VR-HiEta-Mu-IBL1. The band on the expected background estimation indicates the total uncertainty of the estimation. Downward triangle markers at the bottom of the panels indicate there is no events observed in the corresponding bin, while upward triangle markers at the bottom panel indicate the observed data is beyond the range.

The observed $p_{\rm T$ distribution in the SR-Trk-IBL0_Low signal-region bin. The band on the expected background indicates the total uncertainty of the estimation. Several representative signal models are overlaid. Events outside the shown range are accumulated in the rightmost bin indicated as 'Overflow'. Downward triangle markers at the bottom of the panels indicate that no events are observed in the corresponding mass bin, while upward triangle markers in the lower panels indicate that the observed data is beyond the range.

The observed $p_{\rm T$ distribution in the SR-Mu-IBL0_Low signal-region bin. The band on the expected background indicates the total uncertainty of the estimation. Several representative signal models are overlaid. Events outside the shown range are accumulated in the rightmost bin indicated as 'Overflow'. Downward triangle markers at the bottom of the panels indicate that no events are observed in the corresponding mass bin, while upward triangle markers in the lower panels indicate that the observed data is beyond the range.

The observed $p_{\rm T$ distribution in the SR-Trk-IBL0_High signal-region bin. The band on the expected background indicates the total uncertainty of the estimation. Several representative signal models are overlaid. Events outside the shown range are accumulated in the rightmost bin indicated as 'Overflow'. Downward triangle markers at the bottom of the panels indicate that no events are observed in the corresponding mass bin, while upward triangle markers in the lower panels indicate that the observed data is beyond the range.

The observed $p_{\rm T$ distribution in the SR-Mu-IBL0_High signal-region bin. The band on the expected background indicates the total uncertainty of the estimation. Several representative signal models are overlaid. Events outside the shown range are accumulated in the rightmost bin indicated as 'Overflow'. Downward triangle markers at the bottom of the panels indicate that no events are observed in the corresponding mass bin, while upward triangle markers in the lower panels indicate that the observed data is beyond the range.

The observed $p_{\rm T$ distribution in the SR-Trk-IBL1 signal-region bin. The band on the expected background indicates the total uncertainty of the estimation. Several representative signal models are overlaid. Events outside the shown range are accumulated in the rightmost bin indicated as 'Overflow'. Downward triangle markers at the bottom of the panels indicate that no events are observed in the corresponding mass bin, while upward triangle markers in the lower panels indicate that the observed data is beyond the range.

The observed $p_{\rm T$ distribution in the SR-Mu-IBL1 signal-region bin. The band on the expected background indicates the total uncertainty of the estimation. Several representative signal models are overlaid. Events outside the shown range are accumulated in the rightmost bin indicated as 'Overflow'. Downward triangle markers at the bottom of the panels indicate that no events are observed in the corresponding mass bin, while upward triangle markers in the lower panels indicate that the observed data is beyond the range.

The observed dE/dx distribution in the SR-Trk-IBL0_Low signal-region bin. The band on the expected background indicates the total uncertainty of the estimation. Several representative signal models are overlaid. Events outside the shown range are accumulated in the rightmost bin indicated as 'Overflow'. Downward triangle markers at the bottom of the panels indicate that no events are observed in the corresponding mass bin, while upward triangle markers in the lower panels indicate that the observed data is beyond the range.

The observed dE/dx distribution in the SR-Mu-IBL0_Low signal-region bin. The band on the expected background indicates the total uncertainty of the estimation. Several representative signal models are overlaid. Events outside the shown range are accumulated in the rightmost bin indicated as 'Overflow'. Downward triangle markers at the bottom of the panels indicate that no events are observed in the corresponding mass bin, while upward triangle markers in the lower panels indicate that the observed data is beyond the range.

The observed dE/dx distribution in the SR-Trk-IBL0_High signal-region bin. The band on the expected background indicates the total uncertainty of the estimation. Several representative signal models are overlaid. Events outside the shown range are accumulated in the rightmost bin indicated as 'Overflow'. Downward triangle markers at the bottom of the panels indicate that no events are observed in the corresponding mass bin, while upward triangle markers in the lower panels indicate that the observed data is beyond the range.

The observed dE/dx distribution in the SR-Mu-IBL0_High signal-region bin. The band on the expected background indicates the total uncertainty of the estimation. Several representative signal models are overlaid. Events outside the shown range are accumulated in the rightmost bin indicated as 'Overflow'. Downward triangle markers at the bottom of the panels indicate that no events are observed in the corresponding mass bin, while upward triangle markers in the lower panels indicate that the observed data is beyond the range.

The observed dE/dx distribution in the SR-Trk-IBL1 signal-region bin. The band on the expected background indicates the total uncertainty of the estimation. Several representative signal models are overlaid. Events outside the shown range are accumulated in the rightmost bin indicated as 'Overflow'. Downward triangle markers at the bottom of the panels indicate that no events are observed in the corresponding mass bin, while upward triangle markers in the lower panels indicate that the observed data is beyond the range.

The observed dE/dx distribution in the SR-Mu-IBL1 signal-region bin. The band on the expected background indicates the total uncertainty of the estimation. Several representative signal models are overlaid. Events outside the shown range are accumulated in the rightmost bin indicated as 'Overflow'. Downward triangle markers at the bottom of the panels indicate that no events are observed in the corresponding mass bin, while upward triangle markers in the lower panels indicate that the observed data is beyond the range.

Expected and observed distributions in SR-Inclusive_Low of missing transverse momentum. The expected background distribution is calculated for each |eta| slice using CR-kin control region as the template and applying the scale factor using the dE/dx distribution in CR-dEdx of the corresponding |eta| slice. The last bins of the plots include overflow events above the range.

Expected and observed distributions in SR-Inclusive_High of missing transverse momentum. The expected background distribution is calculated for each |eta| slice using CR-kin control region as the template and applying the scale factor using the dE/dx distribution in CR-dEdx of the corresponding |eta| slice. The last bins of the plots include overflow events above the range.

Expected and observed distributions in SR-Inclusive_Low of relative phi-angle between pTmiss and the signal candidate track. The expected background distribution is calculated for each |eta| slice using CR-kin control region as the template and applying the scale factor using the dE/dx distribution in CR-dEdx of the corresponding |eta| slice. The last bins of the plots include overflow events above the range.

Expected and observed distributions in SR-Inclusive_High of relative phi-angle between pTmiss and the signal candidate track. The expected background distribution is calculated for each |eta| slice using CR-kin control region as the template and applying the scale factor using the dE/dx distribution in CR-dEdx of the corresponding |eta| slice. The last bins of the plots include overflow events above the range.

Expected and observed distributions in SR-Inclusive_Low of the transverse mass of pTmiss and the signal candidate track. The expected background distribution is calculated for each |eta| slice using CR-kin control region as the template and applying the scale factor using the dE/dx distribution in CR-dEdx of the corresponding |eta| slice. The last bins of the plots include overflow events above the range.

Expected and observed distributions in SR-Inclusive_High of the transverse mass of pTmiss and the signal candidate track. The expected background distribution is calculated for each |eta| slice using CR-kin control region as the template and applying the scale factor using the dE/dx distribution in CR-dEdx of the corresponding |eta| slice. The last bins of the plots include overflow events above the range.

Expected and observed distributions in SR-Inclusive_Low of the leading jet pT, required to be separated by at least deltaR > 0.4 with respect to the signal candidate track. The expected background distribution is calculated for each |eta| slice using CR-kin control region as the template and applying the scale factor using the dE/dx distribution in CR-dEdx of the corresponding |eta| slice. The last bins of the plots include overflow events above the range.

Expected and observed distributions in SR-Inclusive_High of the leading jet pT, required to be separated by at least deltaR > 0.4 with respect to the signal candidate track. The expected background distribution is calculated for each |eta| slice using CR-kin control region as the template and applying the scale factor using the dE/dx distribution in CR-dEdx of the corresponding |eta| slice. The last bins of the plots include overflow events above the range.

Expected and observed distributions in SR-Inclusive_Low of the relative phi-angle between the leading jet pT, required to be separated by at least deltaR > 0.4 with respect to the signal candidate track, and the signal candidate track. The expected background distribution is calculated for each |eta| slice using CR-kin control region as the template and applying the scale factor using the dE/dx distribution in CR-dEdx of the corresponding |eta| slice. The last bins of the plots include overflow events above the range.

Expected and observed distributions in SR-Inclusive_High of the relative phi-angle between the leading jet pT, required to be separated by at least deltaR > 0.4 with respect to the signal candidate track, and the signal candidate track. The expected background distribution is calculated for each |eta| slice using CR-kin control region as the template and applying the scale factor using the dE/dx distribution in CR-dEdx of the corresponding |eta| slice. The last bins of the plots include overflow events above the range.

Expected and observed distributions in SR-Inclusive_Low of the relative phi-angle between pTmiss and the leading jet pT, required to be separated by at least deltaR > 0.4 with respect to the signal candidate track. The expected background distribution is calculated for each |eta| slice using CR-kin control region as the template and applying the scale factor using the dE/dx distribution in CR-dEdx of the corresponding |eta| slice. The last bins of the plots include overflow events above the range.

Expected and observed distributions in SR-Inclusive_High of the relative phi-angle between pTmiss and the leading jet pT, required to be separated by at least deltaR > 0.4 with respect to the signal candidate track. The expected background distribution is calculated for each |eta| slice using CR-kin control region as the template and applying the scale factor using the dE/dx distribution in CR-dEdx of the corresponding |eta| slice. The last bins of the plots include overflow events above the range.

Expected and observed distributions in SR-Inclusive_Low of the transverse mass of pTmiss and the leading jet pT, required to be separated by at least deltaR > 0.4 with respect to the signal candidate track. The expected background distribution is calculated for each |eta| slice using CR-kin control region as the template and applying the scale factor using the dE/dx distribution in CR-dEdx of the corresponding |eta| slice. The last bins of the plots include overflow events above the range.

Expected and observed distributions in SR-Inclusive_High of the transverse mass of pTmiss and the leading jet pT, required to be separated by at least deltaR > 0.4 with respect to the signal candidate track. The expected background distribution is calculated for each |eta| slice using CR-kin control region as the template and applying the scale factor using the dE/dx distribution in CR-dEdx of the corresponding |eta| slice. The last bins of the plots include overflow events above the range.

Expected and observed distributions in SR-Inclusive_Low of the effective mass, defined as the scalar sum pT of the signal candidate track, jets satisfying pT > 30 GeV, excluding ones within deltaR < 0.4 with respect to the signal candidate track, and pTmiss. The expected background distribution is calculated for each |eta| slice using CR-kin control region as the template and applying the scale factor using the dE/dx distribution in CR-dEdx of the corresponding |eta| slice. The last bins of the plots include overflow events above the range.

Expected and observed distributions in SR-Inclusive_High of the effective mass, defined as the scalar sum pT of the signal candidate track, jets satisfying pT > 30 GeV, excluding ones within deltaR < 0.4 with respect to the signal candidate track, and pTmiss. The expected background distribution is calculated for each |eta| slice using CR-kin control region as the template and applying the scale factor using the dE/dx distribution in CR-dEdx of the corresponding |eta| slice. The last bins of the plots include overflow events above the range.

The expected upper limits on cross-section for gluinos with $m(\tilde{\chi}_{1}^{0}) = 100 \text{GeV}$, with lifetime with lifetime (a) 1 ns, (b) 3 ns, (c) 10 ns, (d) 30 ns, and (e) stable.

The expected upper limits on cross-section for gluinos with $m(\tilde{\chi}_{1}^{0}) = 100 \text{GeV}$, with lifetime with lifetime (a) 1 ns, (b) 3 ns, (c) 10 ns, (d) 30 ns, and (e) stable.

The expected upper limits on cross-section for gluinos with $m(\tilde{\chi}_{1}^{0}) = 100 \text{GeV}$, with lifetime with lifetime (a) 1 ns, (b) 3 ns, (c) 10 ns, (d) 30 ns, and (e) stable.

The expected upper limits on cross-section for gluinos with $m(\tilde{\chi}_{1}^{0}) = 100 \text{GeV}$, with lifetime with lifetime (a) 1 ns, (b) 3 ns, (c) 10 ns, (d) 30 ns, and (e) stable.

The expected upper limits on cross-section for gluinos with $m(\tilde{\chi}_{1}^{0}) = 100 \text{GeV}$, with lifetime with lifetime (a) 1 ns, (b) 3 ns, (c) 10 ns, (d) 30 ns, and (e) stable.

The expected upper limits on cross-section for gluinos with $\Delta m(\tilde{g}, \tilde{\chi}_{1}^{0}) = 30 \text{GeV}$, with lifetime (a) 1 ns, (b) 3 ns, (c) 10 ns, and (d) 30 ns.

The expected upper limits on cross-section for gluinos with $\Delta m(\tilde{g}, \tilde{\chi}_{1}^{0}) = 30 \text{GeV}$, with lifetime (a) 1 ns, (b) 3 ns, (c) 10 ns, and (d) 30 ns.

The expected upper limits on cross-section for gluinos with $\Delta m(\tilde{g}, \tilde{\chi}_{1}^{0}) = 30 \text{GeV}$, with lifetime (a) 1 ns, (b) 3 ns, (c) 10 ns, and (d) 30 ns.

The expected upper limits on cross-section for gluinos with $\Delta m(\tilde{g}, \tilde{\chi}_{1}^{0}) = 30 \text{GeV}$, with lifetime (a) 1 ns, (b) 3 ns, (c) 10 ns, and (d) 30 ns.

The expected upper limits on cross-section for charginos with lifetime (c) 10 ns, (d) 30 ns, and (e) stable.

The expected upper limits on cross-section for charginos with lifetime (c) 10 ns, (d) 30 ns, and (e) stable.

The expected upper limits on cross-section for charginos with lifetime (c) 10 ns, (d) 30 ns, and (e) stable.

The expected upper limits on cross-section for charginos with lifetime (c) 10 ns, (d) 30 ns, and (e) stable.

The expected upper limits on cross-section for charginos with lifetime (c) 10 ns, (d) 30 ns, and (e) stable.

The expected upper limits on cross-section for sleptons with lifetime (a) 1 ns, (b) 3 ns, (c) 10 ns, (d) 30 ns, and (e) stable.

The expected upper limits on cross-section for sleptons with lifetime (a) 1 ns, (b) 3 ns, (c) 10 ns, (d) 30 ns, and (e) stable.

The expected upper limits on cross-section for sleptons with lifetime (a) 1 ns, (b) 3 ns, (c) 10 ns, (d) 30 ns, and (e) stable.

The expected upper limits on cross-section for sleptons with lifetime (a) 1 ns, (b) 3 ns, (c) 10 ns, (d) 30 ns, and (e) stable.

The expected upper limits on cross-section for sleptons with lifetime (a) 1 ns, (b) 3 ns, (c) 10 ns, (d) 30 ns, and (e) stable.

Muon reconstruction efficiency as a function of β and |η| for (a) stable charginos and (b) stable charged R-hadrons. For weakly interacting LLPs with calorimeter materials the efficiency for the chargino is recommended to refer to. The muon reconstruction efficiency for R-hadrons is significantly lower due to having QCD interactions with materials.

Muon reconstruction efficiency as a function of β and |η| for (a) stable charginos and (b) stable charged R-hadrons. For weakly interacting LLPs with calorimeter materials the efficiency for the chargino is recommended to refer to. The muon reconstruction efficiency for R-hadrons is significantly lower due to having QCD interactions with materials.

Trigger and event selection efficiencies. The band on the marker indicates a typical size of fluctuation by the LLP mass and lifetime observed by the samples used in efficiency derivation, but it does not indicate the full envelope of model dependence.

Trigger and event selection efficiencies. The band on the marker indicates a typical size of fluctuation by the LLP mass and lifetime observed by the samples used in efficiency derivation, but it does not indicate the full envelope of model dependence.

Signal track selection efficiency as a function of CLLP $\beta\gamma$ for SR-Inclusive_Low and SR-Inclusive_High bins. The band on the marker indicates a typical size of fluctuation by the LLP mass and lifetime observed by the samples used in efficiency derivation, but it does not indicate the full envelope of model dependence.

Signal selection efficiency by the mass window for SR-Inclusive_Low and SR-Inclusive_High bins.

Acceptance for the R-hadron pair-production model with m(N1) = 100 GeV for various masses and lifetimes. The acceptance is defined as the fraction of events having at least one charged LLP satisfying pT > 120 GeV, |\eta| < 1.8 and r_decay > 500 mm.

Acceptance for the R-hadron pair-production model with DeltaM(gluino, N1) = 30 GeV for various masses and lifetimes. The acceptance is defined as the fraction of events having at least one charged LLP satisfying pT > 120 GeV, |eta| < 1.8 and r_decay > 500 mm.

Acceptance for the chargino pair-production model for various masses and lifetimes. The acceptance is defined as the fraction of events having at least one charged LLP satisfying pT > 120 GeV, |\eta| < 1.8 and r_decay > 500 mm.

Acceptance for the stau pair-production model for various masses and lifetimes. The acceptance is defined as the fraction of events having at least one charged LLP satisfying pT > 120 GeV, |\eta| < 1.8 and r_decay > 500 mm.

Event-level efficiency for the R-hadron pair-production model with m(N1) = 100 GeV for various masses and lifetimes. The efficiency is defined as the fraction of events satisfying the selection of trigger, event and jet cleaning, ETmiss and primary vertex requirements per events satisfying the acceptance criteria.

Event-level efficiency for the R-hadron pair-production model with DeltaM(gluino, N1) = 30 GeV for various masses and lifetimes. The efficiency is defined as the fraction of events satisfying the selection of trigger, event and jet cleaning, ETmiss and primary vertex requirements per events satisfying the acceptance criteria.

Event-level efficiency for the chargino pair-production model for various masses and lifetimes. The efficiency is defined as the fraction of events satisfying the selection of trigger, event and jet cleaning, ETmiss and primary vertex requirements per events satisfying the acceptance criteria.

Event-level efficiency for the stau pair-production model for various masses and lifetimes. The efficiency is defined as the fraction of events satisfying the selection of trigger, event and jet cleaning, ETmiss and primary vertex requirements per events satisfying the acceptance criteria.

Efficiency of SR-Inclusive_Highfor the R-hadron pair-production model with m(N1) = 100 GeV for various masses and lifetimes. The efficiency is defined as the ratio of events satisfying the signal region selection to those satisfying the acceptance criteria. The mass window is not applied for the presented numbers.

Efficiency of SR-Inclusive_Highfor the R-hadron pair-production model with DeltaM(gluino, N1) = 30 GeV for various masses and lifetimes. The efficiency is defined as the ratio of events satisfying the signal region selection to those satisfying the acceptance criteria. The mass window is not applied for the presented numbers.

Efficiency of SR-Inclusive_Highfor the chargino pair-production model for various masses and lifetimes. The efficiency is defined as the ratio of events satisfying the signal region selection to those satisfying the acceptance criteria. The mass window is not applied for the presented numbers.

Efficiency of SR-Inclusive_Highfor the stau pair-production model for various masses and lifetimes. The efficiency is defined as the ratio of events satisfying the signal region selection to those satisfying the acceptance criteria. The mass window is not applied for the presented numbers.

Efficiency of SR-Inclusive_Low for the R-hadron pair-production model with m(N1) = 100 GeV for various masses and lifetimes. The efficiency is defined as the ratio of events satisfying the signal region selection to those satisfying the acceptance criteria. The mass window is not applied for the presented numbers.

Efficiency of SR-Inclusive_Low for the R-hadron pair-production model with DeltaM(gluino, N1) = 30 GeV for various masses and lifetimes. The efficiency is defined as the ratio of events satisfying the signal region selection to those satisfying the acceptance criteria. The mass window is not applied for the presented numbers.

Efficiency of SR-Inclusive_Low for the chargino pair-production model for various masses and lifetimes. The efficiency is defined as the ratio of events satisfying the signal region selection to those satisfying the acceptance criteria. The mass window is not applied for the presented numbers.

Efficiency of SR-Inclusive_Low for the stau pair-production model for various masses and lifetimes. The efficiency is defined as the ratio of events satisfying the signal region selection to those satisfying the acceptance criteria. The mass window is not applied for the presented numbers.

Passing events in event selection steps for the R-hadron pair-production model with m(N1) = 100 GeV for various masses and lifetimes.

Passing events in event selection steps for the R-hadron pair-production model with DeltaM(gluino, N1) = 30 GeV for various masses and lifetimes.

Passing events in event selection steps for the chargino pair-production model for various masses and lifetimes.

Passing events in event selection steps for the stau pair-production model for various masses and lifetimes.