Post-fit SR yields of the high-mass channel. The upper panel shows the observed number of events, as well the post-fit background predictions in each region. The bottom panel shows the ratio of the observed data and the total background prediction. The shaded areas correspond to the total statistical and systematic uncertainties obtained after the fit and described in Section 6.

Post-fit SR yields of the high-mass channel. The upper panel shows the observed number of events, as well the post-fit background predictions in each region. The bottom panel shows the ratio of the observed data and the total background prediction. The shaded areas correspond to the total statistical and systematic uncertainties obtained after the fit and described in Section 6.

Post-fit SR yields of the high-mass channel. The upper panel shows the observed number of events, as well the post-fit background predictions in each region. The bottom panel shows the ratio of the observed data and the total background prediction. The shaded areas correspond to the total statistical and systematic uncertainties obtained after the fit and described in Section 6.

Post-fit SR yields of the high-mass channel. The upper panel shows the observed number of events, as well the post-fit background predictions in each region. The bottom panel shows the ratio of the observed data and the total background prediction. The shaded areas correspond to the total statistical and systematic uncertainties obtained after the fit and described in Section 6.

Post-fit SR yields of the high-mass channel. The upper panel shows the observed number of events, as well the post-fit background predictions in each region. The bottom panel shows the ratio of the observed data and the total background prediction. The shaded areas correspond to the total statistical and systematic uncertainties obtained after the fit and described in Section 6.

Post-fit SR yields of the high-mass channel. The upper panel shows the observed number of events, as well the post-fit background predictions in each region. The bottom panel shows the ratio of the observed data and the total background prediction. The shaded areas correspond to the total statistical and systematic uncertainties obtained after the fit and described in Section 6.

Post-fit SR yields of the high-mass channel. The upper panel shows the observed number of events, as well the post-fit background predictions in each region. The bottom panel shows the ratio of the observed data and the total background prediction. The shaded areas correspond to the total statistical and systematic uncertainties obtained after the fit and described in Section 6.

Pre-fit data and background (reweighted $2b$) predictions for each $4b$ SR $E_\text{T}^\text{miss}$ and $m_\text{eff}$ bin of the low-mass channel for the 2016 data-taking period. The bottom panel shows the significance of any differences between the observed $4b$ data and the background prediction. The $1\sigma$ and $2\sigma$ bands are shown in green and yellow, respectively. All systematics are included except the background normalization, which is 2.3%.

Pre-fit data and background (reweighted $2b$) predictions for each $4b$ SR $E_\text{T}^\text{miss}$ and $m_\text{eff}$ bin of the low-mass channel for the 2017 data-taking period. The bottom panel shows the significance of any differences between the observed $4b$ data and the background prediction. The $1\sigma$ and $2\sigma$ bands are shown in green and yellow, respectively. All systematics are included except the background normalization, which is 3.7%.

Pre-fit data and background (reweighted $2b$) predictions for each $4b$ SR $E_\text{T}^\text{miss}$ and $m_\text{eff}$ bin of the low-mass channel for the 2018 data-taking period. The bottom panel shows the significance of any differences between the observed $4b$ data and the background prediction. The $1\sigma$ and $2\sigma$ bands are shown in green and yellow, respectively. All systematics are included except the background normalization, which is 1.8%.

Exclusion limits of the low-mass and high-mass channels. The low-mass channel is used for $m_{\tilde{H}}<250$ GeV while the high-mass channel is used for $m_{\tilde{H}}\ge250$ GeV. The plot shows the observed (solid) and expected (dashed) 95% CL upper limits on the cross section of higgsino pair production, assuming a higgsino decay branching ratio of $\mathcal{B}(\tilde{H}\rightarrow h + \tilde{G})=100\%$. The theory cross section and its uncertainty are shown by the solid red line and red shading. Results from a previous ATLAS search using 24.3-36.1 fb$^{-1}$ [13] are shown by the solid (observed) and dashed (expected) blue lines. The bottom panel shows the ratio of the limits to the theory cross section. The phase space above the lines is excluded.

Exclusion limits of the low-mass and high-mass channels. The low-mass channel is used for $m_{\tilde{H}}<250$ GeV while the high-mass channel is used for $m_{\tilde{H}}\ge250$ GeV. The plot shows the observed (solid) and expected (dashed) 95% CL upper limits on the cross section of higgsino pair production, assuming a higgsino decay branching ratio of $\mathcal{B}(\tilde{H}\rightarrow h + \tilde{G})=100\%$. The theory cross section and its uncertainty are shown by the solid red line and red shading. Results from a previous ATLAS search using 24.3-36.1 fb$^{-1}$ [13] are shown by the solid (observed) and dashed (expected) blue lines. The bottom panel shows the ratio of the limits to the theory cross section. The phase space above the lines is excluded.

Exclusion limits of the low-mass and high-mass channels. The low-mass channel is used for $m_{\tilde{H}}<250$ GeV while the high-mass channel is used for $m_{\tilde{H}}\ge250$ GeV. The plot shows the observed (solid) and expected (dashed) 95% CL upper limits on the cross section of higgsino pair production, assuming a higgsino decay branching ratio of $\mathcal{B}(\tilde{H}\rightarrow h + \tilde{G})=100\%$. The theory cross section and its uncertainty are shown by the solid red line and red shading. Results from a previous ATLAS search using 24.3-36.1 fb$^{-1}$ [13] are shown by the solid (observed) and dashed (expected) blue lines. The bottom panel shows the ratio of the limits to the theory cross section. The phase space above the lines is excluded.

Exclusion limits of the low-mass and high-mass channels. The low-mass channel is used for $m_{\tilde{H}}<250$ GeV while the high-mass channel is used for $m_{\tilde{H}}\ge250$ GeV. The plot shows the observed (solid) and expected (dashed) 95% CL upper limits on the cross section of higgsino pair production, assuming a higgsino decay branching ratio of $\mathcal{B}(\tilde{H}\rightarrow h + \tilde{G})=100\%$. The theory cross section and its uncertainty are shown by the solid red line and red shading. Results from a previous ATLAS search using 24.3-36.1 fb$^{-1}$ [13] are shown by the solid (observed) and dashed (expected) blue lines. The bottom panel shows the ratio of the limits to the theory cross section. The phase space above the lines is excluded.

Exclusion limits of the low-mass and high-mass channels. The low-mass channel is used for $m_{\tilde{H}}<250$ GeV while the high-mass channel is used for $m_{\tilde{H}}\ge250$ GeV. The plot shows the observed (solid) and expected (dashed) 95% CL upper limits on the cross section of higgsino pair production, assuming a higgsino decay branching ratio of $\mathcal{B}(\tilde{H}\rightarrow h + \tilde{G})=100\%$. The theory cross section and its uncertainty are shown by the solid red line and red shading. Results from a previous ATLAS search using 24.3-36.1 fb$^{-1}$ [13] are shown by the solid (observed) and dashed (expected) blue lines. The bottom panel shows the ratio of the limits to the theory cross section. The phase space above the lines is excluded.

Exclusion limits of the low-mass and high-mass channels. The low-mass channel is used for $m_{\tilde{H}}<250$ GeV while the high-mass channel is used for $m_{\tilde{H}}\ge250$ GeV. The plot shows the observed (solid) and expected (dashed) 95% CL upper limits on the cross section of higgsino pair production, assuming a higgsino decay branching ratio of $\mathcal{B}(\tilde{H}\rightarrow h + \tilde{G})=100\%$. The theory cross section and its uncertainty are shown by the solid red line and red shading. Results from a previous ATLAS search using 24.3-36.1 fb$^{-1}$ [13] are shown by the solid (observed) and dashed (expected) blue lines. The bottom panel shows the ratio of the limits to the theory cross section. The phase space above the lines is excluded.

Exclusion limits of the low-mass and high-mass channels. The low-mass channel is used for $m_{\tilde{H}}<250$ GeV while the high-mass channel is used for $m_{\tilde{H}}\ge250$ GeV. The plot shows the observed (solid) and expected (dashed) 95% CL upper limits on the cross section of higgsino pair production, assuming a higgsino decay branching ratio of $\mathcal{B}(\tilde{H}\rightarrow h + \tilde{G})=100\%$. The theory cross section and its uncertainty are shown by the solid red line and red shading. Results from a previous ATLAS search using 24.3-36.1 fb$^{-1}$ [13] are shown by the solid (observed) and dashed (expected) blue lines. The bottom panel shows the ratio of the limits to the theory cross section. The phase space above the lines is excluded.

Exclusion limits of the low-mass and high-mass channels. The low-mass channel is used for $m_{\tilde{H}}<250$ GeV while the high-mass channel is used for $m_{\tilde{H}}\ge250$ GeV. The plot shows the observed (solid) and expected (dashed) 95% CL upper limits on the cross section of higgsino pair production, assuming a higgsino decay branching ratio of $\mathcal{B}(\tilde{H}\rightarrow h + \tilde{G})=100\%$. The theory cross section and its uncertainty are shown by the solid red line and red shading. Results from a previous ATLAS search using 24.3-36.1 fb$^{-1}$ [13] are shown by the solid (observed) and dashed (expected) blue lines. The bottom panel shows the ratio of the limits to the theory cross section. The phase space above the lines is excluded.

Exclusion limits of the low-mass and high-mass channels. The low-mass channel is used for $m_{\tilde{H}}<250$ GeV while the high-mass channel is used for $m_{\tilde{H}}\ge250$ GeV. The plot shows the observed (solid) and expected (dashed) 95% CL upper limits on the cross section of higgsino pair production, assuming a higgsino decay branching ratio of $\mathcal{B}(\tilde{H}\rightarrow h + \tilde{G})=100\%$. The theory cross section and its uncertainty are shown by the solid red line and red shading. Results from a previous ATLAS search using 24.3-36.1 fb$^{-1}$ [13] are shown by the solid (observed) and dashed (expected) blue lines. The bottom panel shows the ratio of the limits to the theory cross section. The phase space above the lines is excluded.

Exclusion limits of the low-mass and high-mass channels. The low-mass channel is used for $m_{\tilde{H}}<250$ GeV while the high-mass channel is used for $m_{\tilde{H}}\ge250$ GeV. The plot shows the observed (solid) and expected (dashed) 95% CL upper limits on the cross section of higgsino pair production, assuming a higgsino decay branching ratio of $\mathcal{B}(\tilde{H}\rightarrow h + \tilde{G})=100\%$. The theory cross section and its uncertainty are shown by the solid red line and red shading. Results from a previous ATLAS search using 24.3-36.1 fb$^{-1}$ [13] are shown by the solid (observed) and dashed (expected) blue lines. The bottom panel shows the ratio of the limits to the theory cross section. The phase space above the lines is excluded.

Exclusion limits of the low-mass and high-mass channels. The low-mass channel is used for $m_{\tilde{H}}<250$ GeV while the high-mass channel is used for $m_{\tilde{H}}\ge250$ GeV. The plot shows the 95% CL observed (solid) and expected (dashed) upper limits on $\mathcal{B}(\tilde{H}\rightarrow h + \tilde{G})$, assuming the theory cross section for higgsino pair production. The higgsinos are assumed to decay as $\tilde{H}\rightarrow h + \tilde{G}$ or $\tilde{H}\rightarrow Z + \tilde{G}$. The phase space above the lines is excluded.

Exclusion limits of the low-mass and high-mass channels. The low-mass channel is used for $m_{\tilde{H}}<250$ GeV while the high-mass channel is used for $m_{\tilde{H}}\ge250$ GeV. The plot shows the 95% CL observed (solid) and expected (dashed) upper limits on $\mathcal{B}(\tilde{H}\rightarrow h + \tilde{G})$, assuming the theory cross section for higgsino pair production. The higgsinos are assumed to decay as $\tilde{H}\rightarrow h + \tilde{G}$ or $\tilde{H}\rightarrow Z + \tilde{G}$. The phase space above the lines is excluded.

Exclusion limits of the low-mass and high-mass channels. The low-mass channel is used for $m_{\tilde{H}}<250$ GeV while the high-mass channel is used for $m_{\tilde{H}}\ge250$ GeV. The plot shows the 95% CL observed (solid) and expected (dashed) upper limits on $\mathcal{B}(\tilde{H}\rightarrow h + \tilde{G})$, assuming the theory cross section for higgsino pair production. The higgsinos are assumed to decay as $\tilde{H}\rightarrow h + \tilde{G}$ or $\tilde{H}\rightarrow Z + \tilde{G}$. The phase space above the lines is excluded.

Exclusion limits of the low-mass and high-mass channels. The low-mass channel is used for $m_{\tilde{H}}<250$ GeV while the high-mass channel is used for $m_{\tilde{H}}\ge250$ GeV. The plot shows the 95% CL observed (solid) and expected (dashed) upper limits on $\mathcal{B}(\tilde{H}\rightarrow h + \tilde{G})$, assuming the theory cross section for higgsino pair production. The higgsinos are assumed to decay as $\tilde{H}\rightarrow h + \tilde{G}$ or $\tilde{H}\rightarrow Z + \tilde{G}$. The phase space above the lines is excluded.

Exclusion limits of the low-mass and high-mass channels. The low-mass channel is used for $m_{\tilde{H}}<250$ GeV while the high-mass channel is used for $m_{\tilde{H}}\ge250$ GeV. The plot shows the 95% CL observed (solid) and expected (dashed) upper limits on $\mathcal{B}(\tilde{H}\rightarrow h + \tilde{G})$, assuming the theory cross section for higgsino pair production. The higgsinos are assumed to decay as $\tilde{H}\rightarrow h + \tilde{G}$ or $\tilde{H}\rightarrow Z + \tilde{G}$. The phase space above the lines is excluded.

Exclusion limits of the low-mass channel. The plot shows the observed (solid) and expected (dashed) 95% CL upper limits on the cross section of higgsino pair production, assuming a higgsino decay branching ratio of $\mathcal{B}(\tilde{H}\rightarrow h + \tilde{G})=100\%$. The theory cross section and its uncertainty are shown by the solid red line and red shading. The bottom panel shows the ratio of the limits to the theory cross section. The phase space above the lines is excluded.

Exclusion limits of the low-mass channel. The plot shows the observed (solid) and expected (dashed) 95% CL upper limits on the cross section of higgsino pair production, assuming a higgsino decay branching ratio of $\mathcal{B}(\tilde{H}\rightarrow h + \tilde{G})=100\%$. The theory cross section and its uncertainty are shown by the solid red line and red shading. The bottom panel shows the ratio of the limits to the theory cross section. The phase space above the lines is excluded.

Exclusion limits of the high-mass channel. The plot shows the observed (solid) and expected (dashed) 95% CL upper limits on the cross section of higgsino pair production, assuming a higgsino decay branching ratio of $\mathcal{B}(\tilde{H}\rightarrow h + \tilde{G})=100\%$. The theory cross section and its uncertainty are shown by the solid red line and red shading. The bottom panel shows the ratio of the limits to the theory cross section. The phase space above the lines is excluded.

Exclusion limits of the high-mass channel. The plot shows the observed (solid) and expected (dashed) 95% CL upper limits on the cross section of higgsino pair production, assuming a higgsino decay branching ratio of $\mathcal{B}(\tilde{H}\rightarrow h + \tilde{G})=100\%$. The theory cross section and its uncertainty are shown by the solid red line and red shading. The bottom panel shows the ratio of the limits to the theory cross section. The phase space above the lines is excluded.

Results of the background-only fit in the low-mass channel discovery region SR_LM_150. Both pre-fit and post-fit values are shown.

Results of the background-only fit in the low-mass channel discovery region SR_LM_300. Both pre-fit and post-fit values are shown.

The experimental efficiency of the low-mass channel for the exclusion and discovery signal regions as a function of higgsino mass. The experimental efficiency is defined as the number of events passing the detector-level event selections divided by the number of events passing the event selections for a perfect detector. The denominator is obtained by implementing particle-level event selections that emulate the detector-level selections. This treats the lack of availability of $b$-jet triggers as an inefficiency.

The particle-level acceptance for the low-mass exclusion and discovery signal regions, shown as a function of higgsino mass. The acceptance is defined as the fraction of signal events passing the particle-level event selection that emulates the detector-level selection. The acceptance calculation considers only those signal events where both higgsinos decay to Higgs bosons.

The experimental efficiency of the high-mass channel discovery regions as a function of higgsino mass. For each higgsino mass, the efficiency is shown for the SR-1 region corresponding to the mass. For masses above 1100 GeV, SR-1-1100 is used. The experimental efficiency is defined as the number of events passing the detector-level event selections divided by the number of events passing the event selections for a perfect detector. The denominator is obtained by implementing particle-level event selections that emulate the detector-level selections. The efficiency calculation considers only those signal events where both higgsinos decay to Higgs bosons.

The particle-level acceptance for the high-mass signal regions, shown as a function of higgsino mass. For each higgsino mass, the acceptance is shown for the SR-1 region corresponding to the mass. For masses above 1100 GeV, SR-1-1100 is used. The acceptance is defined as the fraction of signal events passing the particle-level event selection that emulates the detector-level selection. The acceptance calculation considers only those signal events where both higgsinos decay to Higgs bosons.

Cutflow for the low-mass channel for a representative 130 GeV signal. The preselection requires 4 or more jets with $p_\text{T}>25$ GeV and 2 or more $b$-jets with $p_\text{T}>25$ GeV. The $b$-jet cut requires 4 or more $b$-jets with $p_\text{T}>40$ GeV. As the samples are generated with $\mathcal{B}(\tilde{H}\rightarrow h\tilde{G})$=50%, $\mathcal{B}(\tilde{H}\rightarrow Z\tilde{G})$=50% to allow for both decays to be studied, the $hh$ events selection is used to select the events where each of the higgsinos decays to a Higgs boson. Expected yields are normalized to a luminosity of 139 fb$^{-1}$, with the availability of $b$-jet triggers lowering the luminosity to 126 fb$^{-1}$. All selections are cumulative, with the exception of the SR cuts, which are each applied separately.

Cutflow for the low-mass channel for a representative 150 GeV signal. The preselection requires 4 or more jets with $p_\text{T}>25$ GeV and 2 or more $b$-jets with $p_\text{T}>25$ GeV. The $b$-jet cut requires 4 or more $b$-jets with $p_\text{T}>40$ GeV. As the samples are generated with $\mathcal{B}(\tilde{H}\rightarrow h\tilde{G})$=50%, $\mathcal{B}(\tilde{H}\rightarrow Z\tilde{G})$=50% to allow for both decays to be studied, the $hh$ events selection is used to select the events where each of the higgsinos decays to a Higgs boson. Expected yields are normalized to a luminosity of 139 fb$^{-1}$, with the availability of $b$-jet triggers lowering the luminosity to 126 fb$^{-1}$. All selections are cumulative, with the exception of the SR cuts, which are each applied separately.

Cutflow for the low-mass channel for a representative 200 GeV signal. The preselection requires 4 or more jets with $p_\text{T}>25$ GeV and 2 or more $b$-jets with $p_\text{T}>25$ GeV. The $b$-jet cut requires 4 or more $b$-jets with $p_\text{T}>40$ GeV. As the samples are generated with $\mathcal{B}(\tilde{H}\rightarrow h\tilde{G})$=50%, $\mathcal{B}(\tilde{H}\rightarrow Z\tilde{G})$=50% to allow for both decays to be studied, the $hh$ events selection is used to select the events where each of the higgsinos decays to a Higgs boson. Expected yields are normalized to a luminosity of 139 fb$^{-1}$, with the availability of $b$-jet triggers lowering the luminosity to 126 fb$^{-1}$. All selections are cumulative, with the exception of the SR cuts, which are each applied separately.

Cutflow for the low-mass channel for a representative 250 GeV signal. The preselection requires 4 or more jets with $p_\text{T}>25$ GeV and 2 or more $b$-jets with $p_\text{T}>25$ GeV. The $b$-jet cut requires 4 or more $b$-jets with $p_\text{T}>40$ GeV. As the samples are generated with $\mathcal{B}(\tilde{H}\rightarrow h\tilde{G})$=50%, $\mathcal{B}(\tilde{H}\rightarrow Z\tilde{G})$=50% to allow for both decays to be studied, the $hh$ events selection is used to select the events where each of the higgsinos decays to a Higgs boson. Expected yields are normalized to a luminosity of 139 fb$^{-1}$, with the availability of $b$-jet triggers lowering the luminosity to 126 fb$^{-1}$. All selections are cumulative, with the exception of the SR cuts, which are each applied separately.

Cutflow for the low-mass channel for a representative 300 GeV signal. The preselection requires 4 or more jets with $p_\text{T}>25$ GeV and 2 or more $b$-jets with $p_\text{T}>25$ GeV. The $b$-jet cut requires 4 or more $b$-jets with $p_\text{T}>40$ GeV. As the samples are generated with $\mathcal{B}(\tilde{H}\rightarrow h\tilde{G})$=50%, $\mathcal{B}(\tilde{H}\rightarrow Z\tilde{G})$=50% to allow for both decays to be studied, the $hh$ events selection is used to select the events where each of the higgsinos decays to a Higgs boson. Expected yields are normalized to a luminosity of 139 fb$^{-1}$, with the availability of $b$-jet triggers lowering the luminosity to 126 fb$^{-1}$. All selections are cumulative, with the exception of the SR cuts, which are each applied separately.

Cutflow for the low-mass channel for a representative 400 GeV signal. The preselection requires 4 or more jets with $p_\text{T}>25$ GeV and 2 or more $b$-jets with $p_\text{T}>25$ GeV. The $b$-jet cut requires 4 or more $b$-jets with $p_\text{T}>40$ GeV. As the samples are generated with $\mathcal{B}(\tilde{H}\rightarrow h\tilde{G})$=50%, $\mathcal{B}(\tilde{H}\rightarrow Z\tilde{G})$=50% to allow for both decays to be studied, the $hh$ events selection is used to select the events where each of the higgsinos decays to a Higgs boson. Expected yields are normalized to a luminosity of 139 fb$^{-1}$, with the availability of $b$-jet triggers lowering the luminosity to 126 fb$^{-1}$. All selections are cumulative, with the exception of the SR cuts, which are each applied separately.

Cutflow for the low-mass channel for a representative 500 GeV signal. The preselection requires 4 or more jets with $p_\text{T}>25$ GeV and 2 or more $b$-jets with $p_\text{T}>25$ GeV. The $b$-jet cut requires 4 or more $b$-jets with $p_\text{T}>40$ GeV. As the samples are generated with $\mathcal{B}(\tilde{H}\rightarrow h\tilde{G})$=50%, $\mathcal{B}(\tilde{H}\rightarrow Z\tilde{G})$=50% to allow for both decays to be studied, the $hh$ events selection is used to select the events where each of the higgsinos decays to a Higgs boson. Expected yields are normalized to a luminosity of 139 fb$^{-1}$, with the availability of $b$-jet triggers lowering the luminosity to 126 fb$^{-1}$. All selections are cumulative, with the exception of the SR cuts, which are each applied separately.

Cutflow for the low-mass channel for a representative 600 GeV signal. The preselection requires 4 or more jets with $p_\text{T}>25$ GeV and 2 or more $b$-jets with $p_\text{T}>25$ GeV. The $b$-jet cut requires 4 or more $b$-jets with $p_\text{T}>40$ GeV. As the samples are generated with $\mathcal{B}(\tilde{H}\rightarrow h\tilde{G})$=50%, $\mathcal{B}(\tilde{H}\rightarrow Z\tilde{G})$=50% to allow for both decays to be studied, the $hh$ events selection is used to select the events where each of the higgsinos decays to a Higgs boson. Expected yields are normalized to a luminosity of 139 fb$^{-1}$, with the availability of $b$-jet triggers lowering the luminosity to 126 fb$^{-1}$. All selections are cumulative, with the exception of the SR cuts, which are each applied separately.

Cutflow for the low-mass channel for a representative 700 GeV signal. The preselection requires 4 or more jets with $p_\text{T}>25$ GeV and 2 or more $b$-jets with $p_\text{T}>25$ GeV. The $b$-jet cut requires 4 or more $b$-jets with $p_\text{T}>40$ GeV. As the samples are generated with $\mathcal{B}(\tilde{H}\rightarrow h\tilde{G})$=50%, $\mathcal{B}(\tilde{H}\rightarrow Z\tilde{G})$=50% to allow for both decays to be studied, the $hh$ events selection is used to select the events where each of the higgsinos decays to a Higgs boson. Expected yields are normalized to a luminosity of 139 fb$^{-1}$, with the availability of $b$-jet triggers lowering the luminosity to 126 fb$^{-1}$. All selections are cumulative, with the exception of the SR cuts, which are each applied separately.

Cutflow for the low-mass channel for a representative 800 GeV signal. The preselection requires 4 or more jets with $p_\text{T}>25$ GeV and 2 or more $b$-jets with $p_\text{T}>25$ GeV. The $b$-jet cut requires 4 or more $b$-jets with $p_\text{T}>40$ GeV. As the samples are generated with $\mathcal{B}(\tilde{H}\rightarrow h\tilde{G})$=50%, $\mathcal{B}(\tilde{H}\rightarrow Z\tilde{G})$=50% to allow for both decays to be studied, the $hh$ events selection is used to select the events where each of the higgsinos decays to a Higgs boson. Expected yields are normalized to a luminosity of 139 fb$^{-1}$, with the availability of $b$-jet triggers lowering the luminosity to 126 fb$^{-1}$. All selections are cumulative, with the exception of the SR cuts, which are each applied separately.

Cutflow for the low-mass channel for a representative 900 GeV signal. The preselection requires 4 or more jets with $p_\text{T}>25$ GeV and 2 or more $b$-jets with $p_\text{T}>25$ GeV. The $b$-jet cut requires 4 or more $b$-jets with $p_\text{T}>40$ GeV. As the samples are generated with $\mathcal{B}(\tilde{H}\rightarrow h\tilde{G})$=50%, $\mathcal{B}(\tilde{H}\rightarrow Z\tilde{G})$=50% to allow for both decays to be studied, the $hh$ events selection is used to select the events where each of the higgsinos decays to a Higgs boson. Expected yields are normalized to a luminosity of 139 fb$^{-1}$, with the availability of $b$-jet triggers lowering the luminosity to 126 fb$^{-1}$. All selections are cumulative, with the exception of the SR cuts, which are each applied separately.

Cutflow for the low-mass channel for a representative 1000 GeV signal. The preselection requires 4 or more jets with $p_\text{T}>25$ GeV and 2 or more $b$-jets with $p_\text{T}>25$ GeV. The $b$-jet cut requires 4 or more $b$-jets with $p_\text{T}>40$ GeV. As the samples are generated with $\mathcal{B}(\tilde{H}\rightarrow h\tilde{G})$=50%, $\mathcal{B}(\tilde{H}\rightarrow Z\tilde{G})$=50% to allow for both decays to be studied, the $hh$ events selection is used to select the events where each of the higgsinos decays to a Higgs boson. Expected yields are normalized to a luminosity of 139 fb$^{-1}$, with the availability of $b$-jet triggers lowering the luminosity to 126 fb$^{-1}$. All selections are cumulative, with the exception of the SR cuts, which are each applied separately.

Cutflow for the low-mass channel for a representative 1100 GeV signal. The preselection requires 4 or more jets with $p_\text{T}>25$ GeV and 2 or more $b$-jets with $p_\text{T}>25$ GeV. The $b$-jet cut requires 4 or more $b$-jets with $p_\text{T}>40$ GeV. As the samples are generated with $\mathcal{B}(\tilde{H}\rightarrow h\tilde{G})$=50%, $\mathcal{B}(\tilde{H}\rightarrow Z\tilde{G})$=50% to allow for both decays to be studied, the $hh$ events selection is used to select the events where each of the higgsinos decays to a Higgs boson. Expected yields are normalized to a luminosity of 139 fb$^{-1}$, with the availability of $b$-jet triggers lowering the luminosity to 126 fb$^{-1}$. All selections are cumulative, with the exception of the SR cuts, which are each applied separately.

Cutflow for the high-mass channel for a representative 200 GeV signal. The preselection requires 4 or more jets with $p_\text{T}>25$ GeV and 2 or more $b$-jets with $p_\text{T}>25$ GeV. As the samples are generated with $\mathcal{B}(\tilde{H}\rightarrow h\tilde{G})$=50%, $\mathcal{B}(\tilde{H}\rightarrow Z\tilde{G})$=50% to allow for both decays to be studied, the $hh$ events selection is used to select the events where each of the higgsinos decays to a Higgs boson. Expected yields are normalized to a luminosity of 139 fb$^{-1}$. All selections are cumulative, with the exception of the SR cuts, which are each applied separately.

Cutflow for the high-mass channel for a representative 250 GeV signal. The preselection requires 4 or more jets with $p_\text{T}>25$ GeV and 2 or more $b$-jets with $p_\text{T}>25$ GeV. As the samples are generated with $\mathcal{B}(\tilde{H}\rightarrow h\tilde{G})$=50%, $\mathcal{B}(\tilde{H}\rightarrow Z\tilde{G})$=50% to allow for both decays to be studied, the $hh$ events selection is used to select the events where each of the higgsinos decays to a Higgs boson. Expected yields are normalized to a luminosity of 139 fb$^{-1}$. All selections are cumulative, with the exception of the SR cuts, which are each applied separately.

Cutflow for the high-mass channel for a representative 300 GeV signal. The preselection requires 4 or more jets with $p_\text{T}>25$ GeV and 2 or more $b$-jets with $p_\text{T}>25$ GeV. As the samples are generated with $\mathcal{B}(\tilde{H}\rightarrow h\tilde{G})$=50%, $\mathcal{B}(\tilde{H}\rightarrow Z\tilde{G})$=50% to allow for both decays to be studied, the $hh$ events selection is used to select the events where each of the higgsinos decays to a Higgs boson. Expected yields are normalized to a luminosity of 139 fb$^{-1}$. All selections are cumulative, with the exception of the SR cuts, which are each applied separately.

Cutflow for the high-mass channel for a representative 400 GeV signal. The preselection requires 4 or more jets with $p_\text{T}>25$ GeV and 2 or more $b$-jets with $p_\text{T}>25$ GeV. As the samples are generated with $\mathcal{B}(\tilde{H}\rightarrow h\tilde{G})$=50%, $\mathcal{B}(\tilde{H}\rightarrow Z\tilde{G})$=50% to allow for both decays to be studied, the $hh$ events selection is used to select the events where each of the higgsinos decays to a Higgs boson. Expected yields are normalized to a luminosity of 139 fb$^{-1}$. All selections are cumulative, with the exception of the SR cuts, which are each applied separately.

Cutflow for the high-mass channel for a representative 500 GeV signal. The preselection requires 4 or more jets with $p_\text{T}>25$ GeV and 2 or more $b$-jets with $p_\text{T}>25$ GeV. As the samples are generated with $\mathcal{B}(\tilde{H}\rightarrow h\tilde{G})$=50%, $\mathcal{B}(\tilde{H}\rightarrow Z\tilde{G})$=50% to allow for both decays to be studied, the $hh$ events selection is used to select the events where each of the higgsinos decays to a Higgs boson. Expected yields are normalized to a luminosity of 139 fb$^{-1}$. All selections are cumulative, with the exception of the SR cuts, which are each applied separately.

Cutflow for the high-mass channel for a representative 600 GeV signal. The preselection requires 4 or more jets with $p_\text{T}>25$ GeV and 2 or more $b$-jets with $p_\text{T}>25$ GeV. As the samples are generated with $\mathcal{B}(\tilde{H}\rightarrow h\tilde{G})$=50%, $\mathcal{B}(\tilde{H}\rightarrow Z\tilde{G})$=50% to allow for both decays to be studied, the $hh$ events selection is used to select the events where each of the higgsinos decays to a Higgs boson. Expected yields are normalized to a luminosity of 139 fb$^{-1}$. All selections are cumulative, with the exception of the SR cuts, which are each applied separately.

Cutflow for the high-mass channel for a representative 700 GeV signal. The preselection requires 4 or more jets with $p_\text{T}>25$ GeV and 2 or more $b$-jets with $p_\text{T}>25$ GeV. As the samples are generated with $\mathcal{B}(\tilde{H}\rightarrow h\tilde{G})$=50%, $\mathcal{B}(\tilde{H}\rightarrow Z\tilde{G})$=50% to allow for both decays to be studied, the $hh$ events selection is used to select the events where each of the higgsinos decays to a Higgs boson. Expected yields are normalized to a luminosity of 139 fb$^{-1}$. All selections are cumulative, with the exception of the SR cuts, which are each applied separately.

Cutflow for the high-mass channel for a representative 800 GeV signal. The preselection requires 4 or more jets with $p_\text{T}>25$ GeV and 2 or more $b$-jets with $p_\text{T}>25$ GeV. As the samples are generated with $\mathcal{B}(\tilde{H}\rightarrow h\tilde{G})$=50%, $\mathcal{B}(\tilde{H}\rightarrow Z\tilde{G})$=50% to allow for both decays to be studied, the $hh$ events selection is used to select the events where each of the higgsinos decays to a Higgs boson. Expected yields are normalized to a luminosity of 139 fb$^{-1}$. All selections are cumulative, with the exception of the SR cuts, which are each applied separately.

Cutflow for the high-mass channel for a representative 900 GeV signal. The preselection requires 4 or more jets with $p_\text{T}>25$ GeV and 2 or more $b$-jets with $p_\text{T}>25$ GeV. As the samples are generated with $\mathcal{B}(\tilde{H}\rightarrow h\tilde{G})$=50%, $\mathcal{B}(\tilde{H}\rightarrow Z\tilde{G})$=50% to allow for both decays to be studied, the $hh$ events selection is used to select the events where each of the higgsinos decays to a Higgs boson. Expected yields are normalized to a luminosity of 139 fb$^{-1}$. All selections are cumulative, with the exception of the SR cuts, which are each applied separately.

Cutflow for the high-mass channel for a representative 1000 GeV signal. The preselection requires 4 or more jets with $p_\text{T}>25$ GeV and 2 or more $b$-jets with $p_\text{T}>25$ GeV. As the samples are generated with $\mathcal{B}(\tilde{H}\rightarrow h\tilde{G})$=50%, $\mathcal{B}(\tilde{H}\rightarrow Z\tilde{G})$=50% to allow for both decays to be studied, the $hh$ events selection is used to select the events where each of the higgsinos decays to a Higgs boson. Expected yields are normalized to a luminosity of 139 fb$^{-1}$. All selections are cumulative, with the exception of the SR cuts, which are each applied separately.

Cutflow for the high-mass channel for a representative 1100 GeV signal. The preselection requires 4 or more jets with $p_\text{T}>25$ GeV and 2 or more $b$-jets with $p_\text{T}>25$ GeV. As the samples are generated with $\mathcal{B}(\tilde{H}\rightarrow h\tilde{G})$=50%, $\mathcal{B}(\tilde{H}\rightarrow Z\tilde{G})$=50% to allow for both decays to be studied, the $hh$ events selection is used to select the events where each of the higgsinos decays to a Higgs boson. Expected yields are normalized to a luminosity of 139 fb$^{-1}$. All selections are cumulative, with the exception of the SR cuts, which are each applied separately.

Cutflow for the high-mass channel for a representative 1200 GeV signal. The preselection requires 4 or more jets with $p_\text{T}>25$ GeV and 2 or more $b$-jets with $p_\text{T}>25$ GeV. As the samples are generated with $\mathcal{B}(\tilde{H}\rightarrow h\tilde{G})$=50%, $\mathcal{B}(\tilde{H}\rightarrow Z\tilde{G})$=50% to allow for both decays to be studied, the $hh$ events selection is used to select the events where each of the higgsinos decays to a Higgs boson. Expected yields are normalized to a luminosity of 139 fb$^{-1}$. All selections are cumulative, with the exception of the SR cuts, which are each applied separately.

Cutflow for the high-mass channel for a representative 1300 GeV signal. The preselection requires 4 or more jets with $p_\text{T}>25$ GeV and 2 or more $b$-jets with $p_\text{T}>25$ GeV. As the samples are generated with $\mathcal{B}(\tilde{H}\rightarrow h\tilde{G})$=50%, $\mathcal{B}(\tilde{H}\rightarrow Z\tilde{G})$=50% to allow for both decays to be studied, the $hh$ events selection is used to select the events where each of the higgsinos decays to a Higgs boson. Expected yields are normalized to a luminosity of 139 fb$^{-1}$. All selections are cumulative, with the exception of the SR cuts, which are each applied separately.

Cutflow for the high-mass channel for a representative 1400 GeV signal. The preselection requires 4 or more jets with $p_\text{T}>25$ GeV and 2 or more $b$-jets with $p_\text{T}>25$ GeV. As the samples are generated with $\mathcal{B}(\tilde{H}\rightarrow h\tilde{G})$=50%, $\mathcal{B}(\tilde{H}\rightarrow Z\tilde{G})$=50% to allow for both decays to be studied, the $hh$ events selection is used to select the events where each of the higgsinos decays to a Higgs boson. Expected yields are normalized to a luminosity of 139 fb$^{-1}$. All selections are cumulative, with the exception of the SR cuts, which are each applied separately.

Cutflow for the high-mass channel for a representative 1500 GeV signal. The preselection requires 4 or more jets with $p_\text{T}>25$ GeV and 2 or more $b$-jets with $p_\text{T}>25$ GeV. As the samples are generated with $\mathcal{B}(\tilde{H}\rightarrow h\tilde{G})$=50%, $\mathcal{B}(\tilde{H}\rightarrow Z\tilde{G})$=50% to allow for both decays to be studied, the $hh$ events selection is used to select the events where each of the higgsinos decays to a Higgs boson. Expected yields are normalized to a luminosity of 139 fb$^{-1}$. All selections are cumulative, with the exception of the SR cuts, which are each applied separately.

Distributions of MaxMin alternative algorithm $m(h_{\mathrm{cand1}},h_{\mathrm{cand2}})_{\mathrm{avg}}$ after the background-only fit. The backgrounds which contribute only a small amount (diboson, W+jets and ttbar+W/Z/h) are grouped and labelled as `Other'.

Distributions of Leading jet $p_T$ after the background-only fit. The backgrounds which contribute only a small amount (diboson, W+jets and ttbar+W/Z/h) are grouped and labelled as `Other'.

Distributions of MaxMin algorithm $m_{hcand}$ after the background-only fit. The backgrounds which contribute only a small amount (diboson, W+jets and ttbar+W/Z/h) are grouped and labelled as `Other'.

Signal efficiency in SRA_M_m60 for simplified models with '$\widetilde{b}\widetilde{b}$ $\rightarrow$ $b\bar{b} \widetilde{\chi}_2^0 \widetilde{\chi}_2^0$ $\rightarrow$ hh$\widetilde{\chi}_1^0 \widetilde{\chi}_1^0$, h $\rightarrow$ $b\bar{b}$' production

Signal acceptance in SRC_28 for simplified models with $\widetilde{b}\widetilde{b}$ $\rightarrow$ $b\bar{b} \widetilde{\chi}_2^0 \widetilde{\chi}_2^0$ $\rightarrow$ hh$\widetilde{\chi}_1^0 \widetilde{\chi}_1^0$, h $\rightarrow$ $b\bar{b}$ production

Signal acceptance in SRC_26 for simplified models with $\widetilde{b}\widetilde{b}$ $\rightarrow$ $b\bar{b} \widetilde{\chi}_2^0 \widetilde{\chi}_2^0$ $\rightarrow$ hh$\widetilde{\chi}_1^0 \widetilde{\chi}_1^0$, h $\rightarrow$ $b\bar{b}$ production

Signal acceptance in SRC_24 for simplified models with $\widetilde{b}\widetilde{b}$ $\rightarrow$ $b\bar{b} \widetilde{\chi}_2^0 \widetilde{\chi}_2^0$ $\rightarrow$ hh$\widetilde{\chi}_1^0 \widetilde{\chi}_1^0$, h $\rightarrow$ $b\bar{b}$ production

Signal acceptance in SRA_M_dm130 for simplified models with $\widetilde{b}\widetilde{b}$ $\rightarrow$ $b\bar{b} \widetilde{\chi}_2^0 \widetilde{\chi}_2^0$ $\rightarrow$ hh$\widetilde{\chi}_1^0 \widetilde{\chi}_1^0$, h $\rightarrow$ $b\bar{b}$ production

Signal acceptance in SRB for simplified models with $\widetilde{b}\widetilde{b}$ $\rightarrow$ $b\bar{b} \widetilde{\chi}_2^0 \widetilde{\chi}_2^0$ $\rightarrow$ hh$\widetilde{\chi}_1^0 \widetilde{\chi}_1^0$, h $\rightarrow$ $b\bar{b}$ production

Signal acceptance in SRA_L_dm130 for simplified models with $\widetilde{b}\widetilde{b}$ $\rightarrow$ $b\bar{b} \widetilde{\chi}_2^0 \widetilde{\chi}_2^0$ $\rightarrow$ hh$\widetilde{\chi}_1^0 \widetilde{\chi}_1^0$, h $\rightarrow$ $b\bar{b}$ production

Signal acceptance in SRC_incl for simplified models with $\widetilde{b}\widetilde{b}$ $\rightarrow$ $b\bar{b} \widetilde{\chi}_2^0 \widetilde{\chi}_2^0$ $\rightarrow$ hh$\widetilde{\chi}_1^0 \widetilde{\chi}_1^0$, h $\rightarrow$ $b\bar{b}$ production

Signal acceptance in SRA_L_m60 for simplified models with $\widetilde{b}\widetilde{b}$ $\rightarrow$ $b\bar{b} \widetilde{\chi}_2^0 \widetilde{\chi}_2^0$ $\rightarrow$ hh$\widetilde{\chi}_1^0 \widetilde{\chi}_1^0$, h $\rightarrow$ $b\bar{b}$ production

Signal acceptance in SRA_incl_dm130 for simplified models with $\widetilde{b}\widetilde{b}$ $\rightarrow$ $b\bar{b} \widetilde{\chi}_2^0 \widetilde{\chi}_2^0$ $\rightarrow$ hh$\widetilde{\chi}_1^0 \widetilde{\chi}_1^0$, h $\rightarrow$ $b\bar{b}$ production

Signal acceptance in SRA_incl_m60 for simplified models with $\widetilde{b}\widetilde{b}$ $\rightarrow$ $b\bar{b} \widetilde{\chi}_2^0 \widetilde{\chi}_2^0$ $\rightarrow$ hh$\widetilde{\chi}_1^0 \widetilde{\chi}_1^0$, h $\rightarrow$ $b\bar{b}$ production

Signal efficiency in SRA_H_m60 for simplified models with '$\widetilde{b}\widetilde{b}$ $\rightarrow$ $b\bar{b} \widetilde{\chi}_2^0 \widetilde{\chi}_2^0$ $\rightarrow$ hh$\widetilde{\chi}_1^0 \widetilde{\chi}_1^0$, h $\rightarrow$ $b\bar{b}$' production

Signal efficiency in SRA_L_dm130 for simplified models with '$\widetilde{b}\widetilde{b}$ $\rightarrow$ $b\bar{b} \widetilde{\chi}_2^0 \widetilde{\chi}_2^0$ $\rightarrow$ hh$\widetilde{\chi}_1^0 \widetilde{\chi}_1^0$, h $\rightarrow$ $b\bar{b}$' production

Signal efficiency in SRB for simplified models with '$\widetilde{b}\widetilde{b}$ $\rightarrow$ $b\bar{b} \widetilde{\chi}_2^0 \widetilde{\chi}_2^0$ $\rightarrow$ hh$\widetilde{\chi}_1^0 \widetilde{\chi}_1^0$, h $\rightarrow$ $b\bar{b}$' production

Signal acceptance in SRC_22 for simplified models with $\widetilde{b}\widetilde{b}$ $\rightarrow$ $b\bar{b} \widetilde{\chi}_2^0 \widetilde{\chi}_2^0$ $\rightarrow$ hh$\widetilde{\chi}_1^0 \widetilde{\chi}_1^0$, h $\rightarrow$ $b\bar{b}$ production

Signal efficiency in SRA_H_dm130 for simplified models with '$\widetilde{b}\widetilde{b}$ $\rightarrow$ $b\bar{b} \widetilde{\chi}_2^0 \widetilde{\chi}_2^0$ $\rightarrow$ hh$\widetilde{\chi}_1^0 \widetilde{\chi}_1^0$, h $\rightarrow$ $b\bar{b}$' production

Signal efficiency in SRC_24 for simplified models with '$\widetilde{b}\widetilde{b}$ $\rightarrow$ $b\bar{b} \widetilde{\chi}_2^0 \widetilde{\chi}_2^0$ $\rightarrow$ hh$\widetilde{\chi}_1^0 \widetilde{\chi}_1^0$, h $\rightarrow$ $b\bar{b}$' production

Signal efficiency in SRC_26 for simplified models with '$\widetilde{b}\widetilde{b}$ $\rightarrow$ $b\bar{b} \widetilde{\chi}_2^0 \widetilde{\chi}_2^0$ $\rightarrow$ hh$\widetilde{\chi}_1^0 \widetilde{\chi}_1^0$, h $\rightarrow$ $b\bar{b}$' production

Signal acceptance in SRA_H_m60 for simplified models with $\widetilde{b}\widetilde{b}$ $\rightarrow$ $b\bar{b} \widetilde{\chi}_2^0 \widetilde{\chi}_2^0$ $\rightarrow$ hh$\widetilde{\chi}_1^0 \widetilde{\chi}_1^0$, h $\rightarrow$ $b\bar{b}$ production

Signal efficiency in SRA_incl_m60 for simplified models with '$\widetilde{b}\widetilde{b}$ $\rightarrow$ $b\bar{b} \widetilde{\chi}_2^0 \widetilde{\chi}_2^0$ $\rightarrow$ hh$\widetilde{\chi}_1^0 \widetilde{\chi}_1^0$, h $\rightarrow$ $b\bar{b}$' production

Signal efficiency in SRC_22 for simplified models with '$\widetilde{b}\widetilde{b}$ $\rightarrow$ $b\bar{b} \widetilde{\chi}_2^0 \widetilde{\chi}_2^0$ $\rightarrow$ hh$\widetilde{\chi}_1^0 \widetilde{\chi}_1^0$, h $\rightarrow$ $b\bar{b}$' production

Signal acceptance in SRA_M_m60 for simplified models with $\widetilde{b}\widetilde{b}$ $\rightarrow$ $b\bar{b} \widetilde{\chi}_2^0 \widetilde{\chi}_2^0$ $\rightarrow$ hh$\widetilde{\chi}_1^0 \widetilde{\chi}_1^0$, h $\rightarrow$ $b\bar{b}$ production

Signal efficiency in SRC_28 for simplified models with '$\widetilde{b}\widetilde{b}$ $\rightarrow$ $b\bar{b} \widetilde{\chi}_2^0 \widetilde{\chi}_2^0$ $\rightarrow$ hh$\widetilde{\chi}_1^0 \widetilde{\chi}_1^0$, h $\rightarrow$ $b\bar{b}$' production

Signal acceptance in SRA_H_dm130 for simplified models with $\widetilde{b}\widetilde{b}$ $\rightarrow$ $b\bar{b} \widetilde{\chi}_2^0 \widetilde{\chi}_2^0$ $\rightarrow$ hh$\widetilde{\chi}_1^0 \widetilde{\chi}_1^0$, h $\rightarrow$ $b\bar{b}$ production

Signal efficiency in SRA_incl_dm130 for simplified models with '$\widetilde{b}\widetilde{b}$ $\rightarrow$ $b\bar{b} \widetilde{\chi}_2^0 \widetilde{\chi}_2^0$ $\rightarrow$ hh$\widetilde{\chi}_1^0 \widetilde{\chi}_1^0$, h $\rightarrow$ $b\bar{b}$' production

Signal efficiency in SRA_L_m60 for simplified models with '$\widetilde{b}\widetilde{b}$ $\rightarrow$ $b\bar{b} \widetilde{\chi}_2^0 \widetilde{\chi}_2^0$ $\rightarrow$ hh$\widetilde{\chi}_1^0 \widetilde{\chi}_1^0$, h $\rightarrow$ $b\bar{b}$' production

Signal efficiency in SRA_M_dm130 for simplified models with '$\widetilde{b}\widetilde{b}$ $\rightarrow$ $b\bar{b} \widetilde{\chi}_2^0 \widetilde{\chi}_2^0$ $\rightarrow$ hh$\widetilde{\chi}_1^0 \widetilde{\chi}_1^0$, h $\rightarrow$ $b\bar{b}$' production

Signal efficiency in SRC_incl for simplified models with '$\widetilde{b}\widetilde{b}$ $\rightarrow$ $b\bar{b} \widetilde{\chi}_2^0 \widetilde{\chi}_2^0$ $\rightarrow$ hh$\widetilde{\chi}_1^0 \widetilde{\chi}_1^0$, h $\rightarrow$ $b\bar{b}$' production

Observed 95% CLs exclusion limit for the $\Delta M(\widetilde{\chi}_{2}^{0},\widetilde{\chi}_{1}^{0})$=130GeV signal grid for the best combined signal regions.

Expected 95% CLs exclusion limit for the $\Delta M(\widetilde{\chi}_{2}^{0},\widetilde{\chi}_{1}^{0})$=130GeV signal grid for the best combined signal regions.

Observed 95% CLs exclusion limit for the $M(\widetilde{\chi}_{1}^{0})$=60GeV signal grid for the best combined signal regions.

Expected 95% CLs exclusion limit for the $M(\widetilde{\chi}_{1}^{0})$=60GeV signal grid for the best combined signal regions.

Model dependent upper limit on the best combined signal regions considered in the $\Delta M(\widetilde{\chi}_{2}^{0},\widetilde{\chi}_{1}^{0})$=130GeV signal grid

Model dependet upper limits on the best combined signal regions considered in the $\Delta M(\widetilde{\chi}_{2}^{0},\widetilde{\chi}_{1}^{0})$=130GeV signal grid

Result of background only fit applied to signal regions. Event yields from the signal regions compared with SM MC predictions for the 3 highest contributing backgrounds separately and combined minor backgrounds.

Expected background event yields and dominant systematic uncertainties on background estimates in the A-type (inclusive), B-type and C-type (inclusive) regions.

Background-only fit results for the A- and B-type regions performed using 139$fb^{-1}$ of data. The quoted uncertainties on the fitted SM background include both the statistical and systematic uncertainties.

Background-only fit results for the C-type region performed using 139$fb^{-1}$ of data. The quoted uncertainties on the fitted SM background include both the statistical and systematic uncertainties.

Observed 95% CL upper limits on the visible cross sections σvis, the observed (S95obs) and expected (S95exp) 95% CL upper limits on the number of signal events with ± 1 σ excursions of the expectation, the CL of the background-only hypothesis, CLB, the discovery p-value (p0), truncated at 0.5, and the associated significance.

Cutflow of the MC events scaled to 139 $fb^{-1}$ for the SRA selections, with a scalar bottom signal of m$(\widetilde{b}_{1},\widetilde{\chi}_2^0,\widetilde{\chi}_1^0) = (1100, 330, 200)$ GeV, considered.

Cutflow of the MC events scaled to 139 $fb^{-1}$ for the SRB selections, with a scalar bottom signal of m$(\widetilde{b}_{1},\widetilde{\chi}_2^0,\widetilde{\chi}_1^0) = (700, 680, 550)$ GeV, considered.

Cutflow of the MC events scaled to 139 $fb^{-1}$ for the SRC selections, with a scalar bottom signal of m$(\widetilde{b}_{1},\widetilde{\chi}_2^0,\widetilde{\chi}_1^0) = (1200, 1150, 60)$ GeV, considered.

Distribution of the missing transverse momentum $E_{\mathrm{T}}^{\mathrm{miss}}$ for the sample satisfying all requirements of the $\mathrm{SR}^{\gamma\gamma}_{W-L}$ selection except the $E_{\mathrm{T}}^{\mathrm{miss}}$ requirement itself. Overlaid are the expected SM backgrounds, separated into the various contributing sources. Also shown are the signal expectations for the ($m_{\tilde{W}}$, $m_{\tilde{\chi}^{0}_{1}}$) = (1000,100) GeV and ($m_{\tilde{W}}$, $m_{\tilde{\chi}^{0}_{1}}$) = (1000,800) GeV models. The vertical dashed lines and right-pointing arrows show the region of the $E_{\mathrm{T}}^{\mathrm{miss}}$ observable selected for inclusion in $\mathrm{SR}^{\gamma\gamma}_{W-L}$ and $\mathrm{SR}^{\gamma\gamma}_{W-H}$. The lower panels show the ratio of observed data to the combined SM expectation. For these plots, the band represents the range of combined statistical and systematic uncertainty in the SM expectation. Events outside the range of the displayed region are included in the highest-value bin.

Distribution of the missing transverse momentum $E_{\mathrm{T}}^{\mathrm{miss}}$ for the sample satisfying all requirements of the $\mathrm{SR}^{\gamma\gamma}_{W-H}$ selection except the $E_{\mathrm{T}}^{\mathrm{miss}}$ requirement itself. Overlaid are the expected SM backgrounds, separated into the various contributing sources. Also shown are the signal expectations for the ($m_{\tilde{W}}$, $m_{\tilde{\chi}^{0}_{1}}$) = (1000,100) GeV and ($m_{\tilde{W}}$, $m_{\tilde{\chi}^{0}_{1}}$) = (1000,800) GeV models. The vertical dashed lines and right-pointing arrows show the region of the $E_{\mathrm{T}}^{\mathrm{miss}}$ observable selected for inclusion in $\mathrm{SR}^{\gamma\gamma}_{W-L}$ and $\mathrm{SR}^{\gamma\gamma}_{W-H}$. The lower panels show the ratio of observed data to the combined SM expectation. For these plots, the band represents the range of combined statistical and systematic uncertainty in the SM expectation. Events outside the range of the displayed region are included in the highest-value bin.

Comparisons between expected and observed content of the validation and signal regions for the photon+jets analysis. The uncertainties in the expected numbers of events are the combined statistical and systematic uncertainties. The lower panel shows the pull (difference between observed and expected event counts normalized by the uncertainty) for each region.

Distribution of the missing transverse momentum $E_{\mathrm{T}}^{\mathrm{miss}}$ for the sample satisfying all requirements of the $\mathrm{SR}^{\gamma j}_{H}$ selection except the $E_{\mathrm{T}}^{\mathrm{miss}}$ requirement itself. Overlaid are the expected SM backgrounds, separated into the various contributing sources. Also shown are the signal expectations for points in the $m_{\tilde{g}}-m_{\tilde{\chi}^{0}_{1}}$ parameter space of the GGM model relevant to the photon+jets analysis (mass values in GeV). The value of the gluino mass arises from the choice $M_{3}$ = 1900 GeV. The $\tilde{\chi}^{0}_{1}$ mass values of 1868, 1920, 442 and 652 GeV arise from the choices $\mu$ = 1810, 1868, 400 and 600 GeV, respectively, combined with the constraint that the branching fraction of $\tilde{\chi}^{0}_{1}$ $\to \gamma \tilde{G}$ be 50%. The vertical dashed lines and right-pointing arrows show the region of the $E_{\mathrm{T}}^{\mathrm{miss}}$ observable selected for inclusion in $\mathrm{SR}^{\gamma j}_{H}$ and $\mathrm{SR}^{\gamma j}_{L}$ for $\mathrm{SR}^{\gamma j}_{L200}$, the $E_{\mathrm{T}}^{\mathrm{miss}}$ requirement is 200 GeV rather than 300 GeV. The lower panels show the ratio of observed data to the combined SM expectation. For these plots, the band represents the range of statistical uncertainty in the SM expectation. Events outside the range of the displayed region are included in the highest-value bin.

Distribution of the missing transverse momentum $E_{\mathrm{T}}^{\mathrm{miss}}$ for the sample satisfying all requirements of the $\mathrm{SR}^{\gamma j}_{L}$ or $\mathrm{SR}^{\gamma j}_{L200}$ selection except the $E_{\mathrm{T}}^{\mathrm{miss}}$ requirement itself. Overlaid are the expected SM backgrounds, separated into the various contributing sources. Also shown are the signal expectations for points in the $m_{\tilde{g}}-m_{\tilde{\chi}^{0}_{1}}$ parameter space of the GGM model relevant to the photon+jets analysis (mass values in GeV). The value of the gluino mass arises from the choice $M_{3}$ = 1900 GeV. The $\tilde{\chi}^{0}_{1}$ mass values of 1868, 1920, 442 and 652 GeV arise from the choices $\mu$ = 1810, 1868, 400 and 600 GeV, respectively, combined with the constraint that the branching fraction of $\tilde{\chi}^{0}_{1} \to \gamma \tilde{G}$ be 50%. The vertical dashed lines and right-pointing arrows show the region of the $E_{\mathrm{T}}^{\mathrm{miss}}$ observable selected for inclusion in $\mathrm{SR}^{\gamma j}_{H}$ and $\mathrm{SR}^{\gamma j}_{L}$ for $\mathrm{SR}^{\gamma j}_{L200}$, the $E_{\mathrm{T}}^{\mathrm{miss}}$ requirement is 200 GeV rather than 300 GeV. The lower panels show the ratio of observed data to the combined SM expectation. For these plots, the band represents the range of statistical uncertainty in the SM expectation. Events outside the range of the displayed region are included in the highest-value bin.

Expected exclusion limits in the gluino-bino mass plane, using the $\mathrm{SR}^{\gamma\gamma}_{S-H}$ analysis for $m_{\tilde{\chi}^{0}_{1}} > 1600$ GeV and the $\mathrm{SR}^{\gamma\gamma}_{S-L}$ analysis for $m_{\tilde{\chi}^{0}_{1}} < 1600$ GeV.

Observed exclusion limits in the gluino--bino mass plane, using the $\mathrm{SR}^{\gamma\gamma}_{S-H}$ analysis for $m_{\tilde{\chi}^{0}_{1}} > 1600$ GeV and the $\mathrm{SR}^{\gamma\gamma}_{S-L}$ analysis for $m_{\tilde{\chi}^{0}_{1}} < 1600$ GeV.

Expected exclusion limit in the squark-bino mass plane, using the $\mathrm{SR}^{\gamma\gamma}_{S-H}$ analysis for $m_{\tilde{\chi}^{0}_{1}} > 900$ GeV and the $\mathrm{SR}^{\gamma\gamma}_{S-L}$ analysis for $m_{\tilde{\chi}^{0}_{1}} < 900$ GeV.

Observed exclusion limit in the squark--bino mass plane, using the $\mathrm{SR}^{\gamma\gamma}_{S-H}$ analysis for $m_{\tilde{\chi}^{0}_{1}} > 900$ GeV and the $\mathrm{SR}^{\gamma\gamma}_{S-L}$ analysis for $m_{\tilde{\chi}^{0}_{1}} < 900$ GeV.

Expected exclusion limit in the wino-bino mass plane, using the $\mathrm{SR}^{\gamma\gamma}_{W-H}$ analysis for $m_{\tilde{\chi}^{0}_{1}} > 400$ GeV and the $\mathrm{SR}^{\gamma\gamma}_{W-L}$ analysis for $m_{\tilde{\chi}^{0}_{1}}$ < 400$ GeV.

Observed exclusion limit in the wino-bino mass plane, using the $\mathrm{SR}^{\gamma\gamma}_{W-H}$ analysis for $m_{\tilde{\chi}^{0}_{1}} > 400$ GeV and the $\mathrm{SR}^{\gamma\gamma}_{W-L}$ analysis for $m_{\tilde{\chi}^{0}_{1}}$ < 400$ GeV.

Expected exclusion limits for the $\mu > 0$ higgsino-bino GGM model explored by the photon+jets analysis.

Observed exclusion limits for the $\mu > 0$ higgsino-bino GGM model explored by the photon+jets analysis.

Distribution of the transverse momentum $p_{\mathrm{T}} (\ell\gamma\gamma)$ of events in the $\ell\gamma\gamma$ control region (except without a cut on $p_{\mathrm{T}} (\ell\gamma\gamma)$). Also shown is the expected contribution from various SM sources, including $W(\to\ell\nu) + \gamma\gamma$ production itself. The displayed uncertainties are a combination of those from all SM sources except $W(\to\ell\nu) + \gamma\gamma$ production, and include statistical and systematic uncertainties.

Distribution of $E_{\mathrm{T}}^{\mathrm{miss}}$ for diphoton events in a validation region defined by a requirement of $H_{\mathrm{T}} > 1750$ GeV. Also shown is the expected contribution from various SM sources, as well as their combined statistical and systematic uncertainties.

Distribution of $H_{\mathrm{T}}$ for diphoton events in a validation region defined by requirement of $E_{\mathrm{T}}^{\mathrm{miss}} > 100$ GeV. Also shown is the expected contribution from various SM sources, as well as their combined statistical and systematic uncertainties.

Distribution of $m_{\mathrm{eff}}$ for events satisfying all requirements $\mathrm{SR}^{\gamma j}_{H}$ save the $m_{\mathrm{eff}}$ requirement itself. Also shown is the expected contribution from various SM sources, and their combined statistical uncertainties.

Distribution of $m_{\mathrm{eff}}$ for events satisfying all requirements $\mathrm{SR}^{\gamma j}_{L}$ save the $m_{\mathrm{eff}}$ requirement itself. Also shown is the expected contribution from various SM sources, and their combined statistical uncertainties.

Derived exclusion limits for the gluino-bino GGM model explored by the diphoton analysis. For each point in the gluino-bino parameter space, the SR ($\mathrm{SR}^{\gamma\gamma}_{S-L}$ or $\mathrm{SR}^{\gamma\gamma}_{S-H}$) that provides the best expected sensitivity is used to estimate the exclusion likelihood. The model dependent upper limits on cross-section (fb) are shown by grey numbers for each signal point.

Derived exclusion limits for the squark-bino GGM model explored by the diphoton analysis. For each point in the squark-bino parameter space, the SR ($\mathrm{SR}^{\gamma\gamma}_{S-L}$ or $\mathrm{SR}^{\gamma\gamma}_{S-H}$) that provides the best expected sensitivity is used to estimate the exclusion likelihood. The model dependent upper limits on cross-section (fb) are shown by grey numbers for each signal point.

Derived exclusion limits for the wino-bino GGM model explored by the diphoton analysis. For each point in the wino-bino parameter space, the SR ($\mathrm{SR}^{\gamma\gamma}_{W-L}$ or $\mathrm{SR}^{\gamma\gamma}_{W-H}$) that provides the best expected sensitivity is used to estimate the exclusion likelihood. The model dependent upper limits on cross-section (fb) are shown by grey numbers for each signal point.

Derived exclusion limits for the $\mu > 0$ higgsino-bino GGM model explored by the photon+jets analysis. For each point in the higgsino-bino parameter space, the SR ($\mathrm{SR}^{\gamma j}_{L}$ or $\mathrm{SR}^{\gamma j}_{H}$) that provides the best expected sensitivity is used to estimate the exclusion likelihood. The model dependent upper limits on cross-section (fb) are shown by grey numbers for each signal point.

Derived exclusion limits for the gluino-bino GGM model explored by the diphoton analysis. For each point in the gluino-bino parameter space, the SR ($\mathrm{SR}^{\gamma\gamma}_{S-L}$ or $\mathrm{SR}^{\gamma\gamma}_{S-H}$) that provides the best expected sensitivity is used to estimate the exclusion likelihood. The labels indicate the best-expected signal region for each point, where SL and SH mean $\mathrm{SR}^{\gamma\gamma}_{S-L}$ and $\mathrm{SR}^{\gamma\gamma}_{S-H}$, respectively.

Derived exclusion limits for the squark-bino GGM model explored by the diphoton analysis. For each point in the squark-bino parameter space, the SR ($\mathrm{SR}^{\gamma\gamma}_{S-L}$ or $\mathrm{SR}^{\gamma\gamma}_{S-H}$) that provides the best expected sensitivity is used to estimate the exclusion likelihood. The labels indicate the best-expected signal region for each point, where SL and SH mean $\mathrm{SR}^{\gamma\gamma}_{S-L}$ and $\mathrm{SR}^{\gamma\gamma}_{S-H}$, respectively.

Derived exclusion limits for the wino--bino GGM model explored by the diphoton analysis. For each point in the wino-bino parameter space, the SR ($\mathrm{SR}^{\gamma\gamma}_{W-L}$ or $\mathrm{SR}^{\gamma\gamma}_{W-H}$) that provides the best expected sensitivity is used to estimate the exclusion likelihood. The labels indicate the best-expected signal region for each point, where WL and WH mean $\mathrm{SR}^{\gamma\gamma}_{W-L}$ and $\mathrm{SR}^{\gamma\gamma}_{W-H}$, respectively.

Derived exclusion limits for the $\mu > 0$ higgsino-bino GGM model explored by the photon+jets analysis. For each point in the higgsino-bino parameter space, the SR ($\mathrm{SR}^{\gamma j}_{L}$ or $\mathrm{SR}^{\gamma j}_{H}$) that provides the best expected sensitivity is used to estimate the exclusion likelihood. The labels indicate the best-expected signal region for each point, where L and H mean $\mathrm{SR}^{\gamma j}_{L}$ and $\mathrm{SR}^{\gamma j}_{H}$, respectively.

Acceptance and efficiency for $\mathrm{SR}^{\gamma\gamma}_{S-L}$ for the signal models of the gluino-bino GGM grid.

Acceptance and efficiency for $\mathrm{SR}^{\gamma\gamma}_{S-H}$ for the signal models of the gluino-bino GGM grid.

Acceptance and efficiency for $\mathrm{SR}^{\gamma\gamma}_{S-L}$ for the signal models of the squark-bino GGM grid.

Acceptance and efficiency for $\mathrm{SR}^{\gamma\gamma}_{S-H}$ for the signal models of the squark-bino GGM grid.

Acceptance and efficiency for $\mathrm{SR}^{\gamma\gamma}_{W-L}$ for the signal models of the wino-bino GGM grid.

Acceptance and efficiency for $\mathrm{SR}^{\gamma\gamma}_{W-H}$ for the signal models of the wino-bino GGM grid.

Acceptance and efficiency for $\mathrm{SR}^{\gamma j}_{L}$ for the signal models of the photon+jets GGM grid.

Acceptance and efficiency for $\mathrm{SR}^{\gamma j}_{H}$ for the signal models of the photon+jets GGM grid.

Cutflow for the $\mathrm{SR}^{\gamma\gamma}_{S-L}$ selection for one relevant signal point in the gluino-bino model, where the gluinos have mass of 1900 GeV and the $\tilde{\chi}^{0}_{1}$ has a mass of 300 GeV (10000 generated events). The numbers are normalized to a luminosity of 36.1 fb$^{-1}$.

Cutflow for the $\mathrm{SR}^{\gamma\gamma}_{S-H}$ selection for one relevant signal point in the gluino-bino model, where the gluinos have mass of 1900 GeV and the $\tilde{\chi}^{0}_{1}$ has a mass of 1700 GeV (10000 generated events). The numbers are normalized to a luminosity of 36.1 fb$^{-1}$.

Cutflow for the $\mathrm{SR}^{\gamma\gamma}_{S-L}$ selection for one relevant signal point in the squark-bino model, where the squarks have mass of 1700 GeV and the $\tilde{\chi}^{0}_{1}$ has a mass of 200 GeV (10000 generated events). The numbers are normalized to a luminosity of 36.1 fb$^{-1}$.

Cutflow for the $\mathrm{SR}^{\gamma\gamma}_{S-H}$ selection for one relevant signal point in the squark-bino model, where the squarks have mass of 1700 GeV and the $\tilde{\chi}^{0}_{1}$ has a mass of 1600 GeV (10000 generated events). The numbers are normalized to a luminosity of 36.1 fb$^{-1}$.

Cutflow for the $\mathrm{SR}^{\gamma\gamma}_{W-L}$ selection for one relevant signal point in the wino-bino model, where the winos have mass of 1000 GeV and the $\tilde{\chi}^{0}_{1}$ has a mass of 200 GeV (10000 generated events). The numbers are normalized to a luminosity of 36.1 fb$^{-1}$.

Cutflow for the $\mathrm{SR}^{\gamma\gamma}_{W-H}$ selection for one relevant signal point in the wino-bino model, where the winos have mass of 1000 GeV and the $\tilde{\chi}^{0}_{1}$ has a mass of 800 GeV (10000 generated events). The numbers are normalized to a luminosity of 36.1 fb$^{-1}$.

Cutflow for the $\mathrm{SR}^{\gamma j}_{L}$ selection, for two relevant signal points in the higgsino-bino model, where the gluinos have mass of 1974 GeV and the $\tilde{\chi}^{0}_{1}$ has a mass of 442 GeV (10000 generated events), and 652 GeV (10000 generated events). The numbers are normalized to a luminosity of 36.1 fb$^{-1}$.

Cutflow for the $\mathrm{SR}^{\gamma j}_{H}$ selection, for two relevant signal points in the higgsino-bino model, where the gluinos have mass of 1974 GeV and the $\tilde{\chi}^{0}_{1}$ has a mass of 1868 GeV (10000 generated events), and 1920 GeV (10000 generated events). The numbers are normalized to a luminosity of 36.1 fb$^{-1}$.

Post-fit yields for data and prediction in each of the single-bin signal regions of the analysis.

Observed 95% confidence level limit for the two-step signal grid.

Observed 95% confidence level limit for the two-step signal grid with the signal cross section increased by one sigma.

Observed 95% confidence level limit for the two-step signal grid with the signal cross section decreased by one sigma.

Expected 95% confidence level limit for the two-step signal grid.

Expected 95% confidence level limit for the two-step signal grid plus one sigma from experimental systematics.

Expected 95% confidence level limit for the two-step signal grid minus one sigma from experimental systematics.

Observed 95% confidence level limit for the Gtt signal grid.

Observed 95% confidence level limit for the Gtt signal grid with the signal cross section increased by one sigma.

Observed 95% confidence level limit for the Gtt signal grid with the signal cross section decreased by one sigma.

Expected 95% confidence level limit for the Gtt signal grid.

Expected 95% confidence level limit for the Gtt signal grid plus one sigma from experimental systematics.

Expected 95% confidence level limit for the Gtt signal grid minus one sigma from experimental systematics.

Observed 95% confidence level limit for the RPV signal grid.

Observed 95% confidence level limit for the RPV signal grid with the signal cross section increased by one sigma.

Observed 95% confidence level limit for the RPV signal grid with the signal cross section decreased by one sigma.

Expected 95% confidence level limit for the RPV signal grid.

Expected 95% confidence level limit for the RPV signal grid plus one sigma from experimental systematics.

Expected 95% confidence level limit for the RPV signal grid minus one sigma from experimental systematics.

Observed 95% confidence level limit for the two-step signal grid.

Expected 95% confidence level limit for the two-step signal grid.

Observed 95% confidence level limit for the Gtt signal grid.

Expected 95% confidence level limit for the Gtt signal grid.

Observed 95% confidence level limit for the RPV signal grid.

Expected 95% confidence level limit for the RPV signal grid.

$\mathcal{S}(E_{\mathrm{T}}^{\mathrm{miss}})$ distribution in the signal region SR-10ij50-0ib-MJ340. Two benchmark signal models are shown along with the background yields. These models, each representing a single mass point, are labelled 'RPV' with $(m_{\tilde{g}}, m_{\tilde{t}}) = (1600, 600) \, \mathrm{GeV}$ and 'two-step' with $(m_{\tilde{g}}, m_{\tilde{\chi^{0}_{1}}}) = (1600, 100) \, \mathrm{GeV}$.

$\mathcal{S}(E_{\mathrm{T}}^{\mathrm{miss}})$ distribution in the signal region SR-12ij50-2ib. Two benchmark signal models are shown along with the background yields. These models, each representing a single mass point, are labelled 'RPV' with $(m_{\tilde{g}}, m_{\tilde{t}}) = (1600, 600) \, \mathrm{GeV}$ and 'two-step' with $(m_{\tilde{g}}, m_{\tilde{\chi^{0}_{1}}}) = (1600, 100) \, \mathrm{GeV}$.

$\mathcal{S}(E_{\mathrm{T}}^{\mathrm{miss}})$ distribution in the signal region SR-9ij80-0ib. Two benchmark signal models are shown along with the background yields. These models, each representing a single mass point, are labelled 'RPV' with $(m_{\tilde{g}}, m_{\tilde{t}}) = (1600, 600) \, \mathrm{GeV}$ and 'two-step' with $(m_{\tilde{g}}, m_{\tilde{\chi^{0}_{1}}}) = (1600, 100) \, \mathrm{GeV}$.

Number of signal events expected for $139 \, \mathrm{fb}^{-1} $ after different analysis selections in the signal region SR-8ij50-0ib-MJ500. This 'two-step' model requires that a strongly produced gluino decays into quarks, the W and Z bosons, and the lightest stable neutralino where $(m_{\tilde{g}}, m_{\tilde{\chi^{0}_{1}}}) = (1600, 100) \, \mathrm{GeV}$.

Number of signal events expected for $139 \, \mathrm{fb}^{-1} $ after different analysis selections in the signal region SR-9ij50-0ib-MJ340. This 'two-step' model requires that a strongly produced gluino decays into quarks, the W and Z bosons, and the lightest stable neutralino where $(m_{\tilde{g}}, m_{\tilde{\chi^{0}_{1}}}) = (1600, 100) \, \mathrm{GeV}$.

Number of signal events expected for $139 \, \mathrm{fb}^{-1} $ after different analysis selections in the signal region SR-10ij50-0ib-MJ340. This 'two-step' model requires that a strongly produced gluino decays into quarks, the W and Z bosons, and the lightest stable neutralino where $(m_{\tilde{g}}, m_{\tilde{\chi^{0}_{1}}}) = (1600, 100) \, \mathrm{GeV}$.

Number of signal events expected for $139 \, \mathrm{fb}^{-1} $ after different analysis selections in the signal region SR-10ij50-0ib-MJ500. This 'two-step' model requires that a strongly produced gluino decays into quarks, the W and Z bosons, and the lightest stable neutralino where $(m_{\tilde{g}}, m_{\tilde{\chi^{0}_{1}}}) = (1600, 100) \, \mathrm{GeV}$.

Number of signal events expected for $139 \, \mathrm{fb}^{-1} $ after different analysis selections in the signal region SR-10ij50-1ib-MJ500. This 'two-step' model requires that a strongly produced gluino decays into quarks, the W and Z bosons, and the lightest stable neutralino where $(m_{\tilde{g}}, m_{\tilde{\chi^{0}_{1}}}) = (1600, 100) \, \mathrm{GeV}$.

Number of signal events expected for $139 \, \mathrm{fb}^{-1} $ after different analysis selections in the signal region SR-11ij50-0ib. This 'two-step' model requires that a strongly produced gluino decays into quarks, the W and Z bosons, and the lightest stable neutralino where $(m_{\tilde{g}}, m_{\tilde{\chi^{0}_{1}}}) = (1600, 100) \, \mathrm{GeV}$.

Number of signal events expected for $139 \, \mathrm{fb}^{-1} $ after different analysis selections in the signal region SR-12ij50-2ib. This 'two-step' model requires that a strongly produced gluino decays into quarks, the W and Z bosons, and the lightest stable neutralino where $(m_{\tilde{g}}, m_{\tilde{\chi^{0}_{1}}}) = (1600, 100) \, \mathrm{GeV}$.

Number of signal events expected for $139 \, \mathrm{fb}^{-1} $ after different analysis selections in the signal region SR-9ij80-0ib. This 'two-step' model requires that a strongly produced gluino decays into quarks, the W and Z bosons, and the lightest stable neutralino where $(m_{\tilde{g}}, m_{\tilde{\chi^{0}_{1}}}) = (1600, 100) \, \mathrm{GeV}$.

Acceptance for the signal region SR-8ij50-0ib-MJ500 showing the acceptance for the complete two-step signal grid.

Efficiency for the signal region SR-8ij50-0ib-MJ500 showing the efficiency for the complete two-step signal grid.

Acceptance for the signal region SR-9ij50-0ib-MJ340 showing the acceptance for the complete two-step signal grid.

Efficiency for the signal region SR-9ij50-0ib-MJ340 showing the efficiency for the complete two-step signal grid.

Acceptance for the signal region SR-10ij50-0ib-MJ340 showing the acceptance for the complete two-step signal grid.

Efficiency for the signal region SR-10ij50-0ib-MJ340 showing the efficiency for the complete two-step signal grid.

Acceptance for the signal region SR-10ij50-0ib-MJ500 showing the acceptance for the complete two-step signal grid.

Efficiency for the signal region SR-10ij50-0ib-MJ500 showing the efficiency for the complete two-step signal grid.

Acceptance for the signal region SR-10ij50-1ib-MJ500 showing the acceptance for the complete two-step signal grid.

Efficiency for the signal region SR-10ij50-1ib-MJ500 showing the efficiency for the complete two-step signal grid.

Acceptance for the signal region SR-11ij50-0ib showing the acceptance for the complete two-step signal grid.

Efficiency for the signal region SR-11ij50-0ib showing the efficiency for the complete two-step signal grid.

Acceptance for the signal region SR-12ij50-2ib showing the acceptance for the complete two-step signal grid.

Efficiency for the signal region SR-12ij50-2ib showing the efficiency for the complete two-step signal grid.

Acceptance for the signal region SR-9ij80-0ib showing the acceptance for the complete two-step signal grid.

Efficiency for the signal region SR-9ij80-0ib showing the efficiency for the complete two-step signal grid.

The normalisation factors for the dominant backgrounds of the analysis in each of the multi-bin and single-bin regions.

Post-fit yields for data and prediction in each of the single-bin validation regions to test the $N_{\mathrm{jet}}$ extraction.

Post-fit yields for data and prediction in each of the single-bin validation regions to test the $\mathcal{S}(E_{\mathrm{T}}^{\mathrm{miss}})$ extrapolation.

Post-fit yields for data and prediction in each of the multi-bin validation regions to test the $N_{\mathrm{jet}}$ extraction.

Post-fit yields for data and prediction in each of the multi-bin validation regions to test the $\mathcal{S}(E_{\mathrm{T}}^{\mathrm{miss}})$ extrapolation.

The observed Cls from the best expected signal regions for the two-step decay.

The observed Cls from the best expected signal regions for the Gtt decay.

The observed Cls from the best expected signal regions for the RPV decay.

Number of events in each signal region broken down by background type and the number of observed data events.

From left to right; the $95\%$ CL upper limits on the visible cross section (${\langle \epsilon\sigma \rangle}^{95}_{obs}$) and on the number of signal events. Next is the $95\%$ CL upper limit on the number of signal events, given the expected number of background events. The last two columns show the confidence level for the background only hypothesis ($CL_{b}$) and the dicovery $p$-value along with the Gaussian significance (Z).

Visualisation of the highest jet multiplicity event selected in signal regions targeting long cascade decays of pair-produced gluinos. This event was recorded by ATLAS on 23 October 2016, and contains 16 jets, illustrated by cones. Yellow blocks represent the calorimeter energy measured in noise-suppressed clusters. Of the reconstructed jets, 13 (11) have transverse momenta above 50 GeV (80 GeV), with 3 (2) being b-tagged. The leading jet has a transverse momentum of 507 GeV, and the sum of jet transverse momenta $H_T=2.9$ TeV. A value of 343 GeV is observed for the $E_{T}^{miss}$, whose direction is shown by the dashed red line, producing a significance $S(E_{T}^{miss})=6.4$. The sum of the masses of large-radius jets is evaluated as $M_{J}^{\Sigma}=1070$ GeV.

Visualisation of the highest jet multiplicity event selected in a control region used to make predictions of the background from multijet production. This event was recorded by ATLAS on 18 July 2018, and contains 19 jets, illustrated by cones. Yellow blocks represent the calorimeter energy measured in in noise-suppressed clusters. Of the reconstructed jets, 16 (10) have transverse momenta above 50 GeV (80 GeV). No jets were b-tagged. The leading et has a transverse momentum of 371 GeV, and the sum of jet transverse momenta $H_T=2.2$ TeV. A value of 8 GeV is observed for the $E_{T}^{miss}$, whose direction is shown by the dashed red line, producing a significance $S(E_{T}^{miss})=0.2$. The sum of the masses of large-radius jets is evaluated as $M_{J}^{\Sigma}=767$ GeV.

Data bootstraps post unfolding for the normalised differential cross-section in the phase space with at least three $b$-jets as a function of the number of $b$-jets compared with predictions. The replicas are obtained by reweighting each observed data event by a random integer generated according to Poisson statistics, using the BootstrapGenerator software package (https://gitlab.cern.ch/atlas-physics/sm/StandardModelTools_BootstrapGenerator/BootstrapGenerator), which implements a technique described in ATL-PHYS-PUB-2021-011 (https://cds.cern.ch/record/2759945). The ATLAS event number and run number of each event are used as seed to uniquely but reproducibly initialise the random number generator for each event. The last bin contains the overflow.

Data bootstraps post unfolding for the normalised differential cross-section in the phase space with at least three $b$-jets as a function of the number of $l/c$-jets compared with predictions. The replicas are obtained by reweighting each observed data event by a random integer generated according to Poisson statistics, using the BootstrapGenerator software package (https://gitlab.cern.ch/atlas-physics/sm/StandardModelTools_BootstrapGenerator/BootstrapGenerator), which implements a technique described in ATL-PHYS-PUB-2021-011 (https://cds.cern.ch/record/2759945). The ATLAS event number and run number of each event are used as seed to uniquely but reproducibly initialise the random number generator for each event. The last bin contains the overflow.

Data bootstraps post unfolding for the normalised differential cross-section in the phase space with at least three $b$-jets as a function of $H_{\text{T}}^{\text{had}}$ compared with predictions. The replicas are obtained by reweighting each observed data event by a random integer generated according to Poisson statistics, using the BootstrapGenerator software package (https://gitlab.cern.ch/atlas-physics/sm/StandardModelTools_BootstrapGenerator/BootstrapGenerator), which implements a technique described in ATL-PHYS-PUB-2021-011 (https://cds.cern.ch/record/2759945). The ATLAS event number and run number of each event are used as seed to uniquely but reproducibly initialise the random number generator for each event. The last bin contains the overflow.

Data bootstraps post unfolding for the normalised differential cross-section in the phase space with at least three $b$-jets as a function of $\Delta R_{\text{avg}}^{bb}$ compared with predictions. The replicas are obtained by reweighting each observed data event by a random integer generated according to Poisson statistics, using the BootstrapGenerator software package (https://gitlab.cern.ch/atlas-physics/sm/StandardModelTools_BootstrapGenerator/BootstrapGenerator), which implements a technique described in ATL-PHYS-PUB-2021-011 (https://cds.cern.ch/record/2759945). The ATLAS event number and run number of each event are used as seed to uniquely but reproducibly initialise the random number generator for each event. The last bin contains the overflow.

Data bootstraps post unfolding for the normalised differential cross-section in the phase space with at least three $b$-jets as a function of $p_{\text{T}}(b_{1})$ compared with predictions. The replicas are obtained by reweighting each observed data event by a random integer generated according to Poisson statistics, using the BootstrapGenerator software package (https://gitlab.cern.ch/atlas-physics/sm/StandardModelTools_BootstrapGenerator/BootstrapGenerator), which implements a technique described in ATL-PHYS-PUB-2021-011 (https://cds.cern.ch/record/2759945). The ATLAS event number and run number of each event are used as seed to uniquely but reproducibly initialise the random number generator for each event. The last bin contains the overflow.

Data bootstraps post unfolding for the normalised differential cross-section in the phase space with at least three $b$-jets as a function of $p_{\text{T}}(b_{2})$ compared with predictions. The replicas are obtained by reweighting each observed data event by a random integer generated according to Poisson statistics, using the BootstrapGenerator software package (https://gitlab.cern.ch/atlas-physics/sm/StandardModelTools_BootstrapGenerator/BootstrapGenerator), which implements a technique described in ATL-PHYS-PUB-2021-011 (https://cds.cern.ch/record/2759945). The ATLAS event number and run number of each event are used as seed to uniquely but reproducibly initialise the random number generator for each event. The last bin contains the overflow.

Data bootstraps post unfolding for the normalised differential cross-section in the phase space with at least three $b$-jets as a function of $p_{\text{T}}(b_{1}^{\text{top}})$ compared with predictions. The replicas are obtained by reweighting each observed data event by a random integer generated according to Poisson statistics, using the BootstrapGenerator software package (https://gitlab.cern.ch/atlas-physics/sm/StandardModelTools_BootstrapGenerator/BootstrapGenerator), which implements a technique described in ATL-PHYS-PUB-2021-011 (https://cds.cern.ch/record/2759945). The ATLAS event number and run number of each event are used as seed to uniquely but reproducibly initialise the random number generator for each event. The last bin contains the overflow.

Data bootstraps post unfolding for the normalised differential cross-section in the phase space with at least three $b$-jets as a function of $p_{\text{T}}(b_{2}^{\text{top}})$ compared with predictions. The replicas are obtained by reweighting each observed data event by a random integer generated according to Poisson statistics, using the BootstrapGenerator software package (https://gitlab.cern.ch/atlas-physics/sm/StandardModelTools_BootstrapGenerator/BootstrapGenerator), which implements a technique described in ATL-PHYS-PUB-2021-011 (https://cds.cern.ch/record/2759945). The ATLAS event number and run number of each event are used as seed to uniquely but reproducibly initialise the random number generator for each event. The last bin contains the overflow.

Data bootstraps post unfolding for the normalised differential cross-section in the phase space with at least three $b$-jets as a function of $p_{\text{T}}(b_{3})$ compared with predictions. The replicas are obtained by reweighting each observed data event by a random integer generated according to Poisson statistics, using the BootstrapGenerator software package (https://gitlab.cern.ch/atlas-physics/sm/StandardModelTools_BootstrapGenerator/BootstrapGenerator), which implements a technique described in ATL-PHYS-PUB-2021-011 (https://cds.cern.ch/record/2759945). The ATLAS event number and run number of each event are used as seed to uniquely but reproducibly initialise the random number generator for each event. The last bin contains the overflow.

Data bootstraps post unfolding for the normalised differential cross-section in the phase space with at least three $b$-jets as a function of $p_{\text{T}}(b_{1}^{\text{add}})$ compared with predictions. The replicas are obtained by reweighting each observed data event by a random integer generated according to Poisson statistics, using the BootstrapGenerator software package (https://gitlab.cern.ch/atlas-physics/sm/StandardModelTools_BootstrapGenerator/BootstrapGenerator), which implements a technique described in ATL-PHYS-PUB-2021-011 (https://cds.cern.ch/record/2759945). The ATLAS event number and run number of each event are used as seed to uniquely but reproducibly initialise the random number generator for each event. The last bin contains the overflow.

Data bootstraps post unfolding for the normalised differential cross-section in the phase space with at least three $b$-jets as a function of $m(b_{1}b_{2})$ compared with predictions. The replicas are obtained by reweighting each observed data event by a random integer generated according to Poisson statistics, using the BootstrapGenerator software package (https://gitlab.cern.ch/atlas-physics/sm/StandardModelTools_BootstrapGenerator/BootstrapGenerator), which implements a technique described in ATL-PHYS-PUB-2021-011 (https://cds.cern.ch/record/2759945). The ATLAS event number and run number of each event are used as seed to uniquely but reproducibly initialise the random number generator for each event. The last bin contains the overflow.

Data bootstraps post unfolding for the normalised differential cross-section in the phase space with at least three $b$-jets as a function of $p_{\text{T}}(b_{1}b_{2})$ compared with predictions. The replicas are obtained by reweighting each observed data event by a random integer generated according to Poisson statistics, using the BootstrapGenerator software package (https://gitlab.cern.ch/atlas-physics/sm/StandardModelTools_BootstrapGenerator/BootstrapGenerator), which implements a technique described in ATL-PHYS-PUB-2021-011 (https://cds.cern.ch/record/2759945). The ATLAS event number and run number of each event are used as seed to uniquely but reproducibly initialise the random number generator for each event. The last bin contains the overflow.

Data bootstraps post unfolding for the normalised differential cross-section in the phase space with at least three $b$-jets as a function of $m(bb^{\text{top}})$ compared with predictions. The replicas are obtained by reweighting each observed data event by a random integer generated according to Poisson statistics, using the BootstrapGenerator software package (https://gitlab.cern.ch/atlas-physics/sm/StandardModelTools_BootstrapGenerator/BootstrapGenerator), which implements a technique described in ATL-PHYS-PUB-2021-011 (https://cds.cern.ch/record/2759945). The ATLAS event number and run number of each event are used as seed to uniquely but reproducibly initialise the random number generator for each event. The last bin contains the overflow.

Data bootstraps post unfolding for the normalised differential cross-section in the phase space with at least three $b$-jets as a function of $p_{\text{T}}(bb^{\text{top}})$ compared with predictions. The replicas are obtained by reweighting each observed data event by a random integer generated according to Poisson statistics, using the BootstrapGenerator software package (https://gitlab.cern.ch/atlas-physics/sm/StandardModelTools_BootstrapGenerator/BootstrapGenerator), which implements a technique described in ATL-PHYS-PUB-2021-011 (https://cds.cern.ch/record/2759945). The ATLAS event number and run number of each event are used as seed to uniquely but reproducibly initialise the random number generator for each event. The last bin contains the overflow.

Data bootstraps post unfolding for the normalised differential cross-section in the phase space with at least three $b$-jets as a function of $\Delta R(e\mu bb^{\text{top}}, b_{1}^{\text{add}})$ compared with predictions. The replicas are obtained by reweighting each observed data event by a random integer generated according to Poisson statistics, using the BootstrapGenerator software package (https://gitlab.cern.ch/atlas-physics/sm/StandardModelTools_BootstrapGenerator/BootstrapGenerator), which implements a technique described in ATL-PHYS-PUB-2021-011 (https://cds.cern.ch/record/2759945). The ATLAS event number and run number of each event are used as seed to uniquely but reproducibly initialise the random number generator for each event. The last bin contains the overflow.

Data bootstraps post unfolding for the normalised differential cross-section in the phase space with at least three $b$-jets and at least one $l/c$-jet as a function of $\Delta R(e\mu bb^{\text{top}}, l/c\text{-jet}_{1})$ compared with predictions. The replicas are obtained by reweighting each observed data event by a random integer generated according to Poisson statistics, using the BootstrapGenerator software package (https://gitlab.cern.ch/atlas-physics/sm/StandardModelTools_BootstrapGenerator/BootstrapGenerator), which implements a technique described in ATL-PHYS-PUB-2021-011 (https://cds.cern.ch/record/2759945). The ATLAS event number and run number of each event are used as seed to uniquely but reproducibly initialise the random number generator for each event. The last bin contains the overflow.

Data bootstraps post unfolding for the normalised differential cross-section in the phase space with at least three $b$-jets and at least one $l/c$-jet as a function of $p_{\text{T}}(l/c\text{-jet}_{1})$ compared with predictions. The replicas are obtained by reweighting each observed data event by a random integer generated according to Poisson statistics, using the BootstrapGenerator software package (https://gitlab.cern.ch/atlas-physics/sm/StandardModelTools_BootstrapGenerator/BootstrapGenerator), which implements a technique described in ATL-PHYS-PUB-2021-011 (https://cds.cern.ch/record/2759945). The ATLAS event number and run number of each event are used as seed to uniquely but reproducibly initialise the random number generator for each event. The last bin contains the overflow.

Data bootstraps post unfolding for the normalised differential cross-section in the phase space with at least three $b$-jets and at least one $l/c$-jet as a function of $p_{\text{T}}(l/c\text{-jet}_{1}) - p_{\text{T}}(b_{1}^{\text{add}})$ compared with predictions. The replicas are obtained by reweighting each observed data event by a random integer generated according to Poisson statistics, using the BootstrapGenerator software package (https://gitlab.cern.ch/atlas-physics/sm/StandardModelTools_BootstrapGenerator/BootstrapGenerator), which implements a technique described in ATL-PHYS-PUB-2021-011 (https://cds.cern.ch/record/2759945). The ATLAS event number and run number of each event are used as seed to uniquely but reproducibly initialise the random number generator for each event. The last bin contains the overflow.

Data bootstraps post unfolding for the normalised differential cross-section in the phase space with at least four $b$-jets as a function of $m(bb^{\text{min}\Delta R})$ compared with predictions. The replicas are obtained by reweighting each observed data event by a random integer generated according to Poisson statistics, using the BootstrapGenerator software package (https://gitlab.cern.ch/atlas-physics/sm/StandardModelTools_BootstrapGenerator/BootstrapGenerator), which implements a technique described in ATL-PHYS-PUB-2021-011 (https://cds.cern.ch/record/2759945). The ATLAS event number and run number of each event are used as seed to uniquely but reproducibly initialise the random number generator for each event. The last bin contains the overflow.

Data bootstraps post unfolding for the normalised differential cross-section in the phase space with at least four $b$-jets as a function of $p_{\text{T}}(bb^{\text{min}\Delta R})$ compared with predictions. The replicas are obtained by reweighting each observed data event by a random integer generated according to Poisson statistics, using the BootstrapGenerator software package (https://gitlab.cern.ch/atlas-physics/sm/StandardModelTools_BootstrapGenerator/BootstrapGenerator), which implements a technique described in ATL-PHYS-PUB-2021-011 (https://cds.cern.ch/record/2759945). The ATLAS event number and run number of each event are used as seed to uniquely but reproducibly initialise the random number generator for each event. The last bin contains the overflow.

Data bootstraps post unfolding for the normalised differential cross-section in the phase space with at least four $b$-jets as a function of $m(bb^{\text{add}})$ compared with predictions. The replicas are obtained by reweighting each observed data event by a random integer generated according to Poisson statistics, using the BootstrapGenerator software package (https://gitlab.cern.ch/atlas-physics/sm/StandardModelTools_BootstrapGenerator/BootstrapGenerator), which implements a technique described in ATL-PHYS-PUB-2021-011 (https://cds.cern.ch/record/2759945). The ATLAS event number and run number of each event are used as seed to uniquely but reproducibly initialise the random number generator for each event. The last bin contains the overflow.

Data bootstraps post unfolding for the normalised differential cross-section in the phase space with at least four $b$-jets as a function of $p_{\text{T}}(bb^{\text{add}})$ compared with predictions. The replicas are obtained by reweighting each observed data event by a random integer generated according to Poisson statistics, using the BootstrapGenerator software package (https://gitlab.cern.ch/atlas-physics/sm/StandardModelTools_BootstrapGenerator/BootstrapGenerator), which implements a technique described in ATL-PHYS-PUB-2021-011 (https://cds.cern.ch/record/2759945). The ATLAS event number and run number of each event are used as seed to uniquely but reproducibly initialise the random number generator for each event. The last bin contains the overflow.

Data bootstraps post unfolding for the normalised differential cross-section in the phase space with at least three $b$-jets as a function of $|\eta(b_{3})|$ compared with predictions. The replicas are obtained by reweighting each observed data event by a random integer generated according to Poisson statistics, using the BootstrapGenerator software package (https://gitlab.cern.ch/atlas-physics/sm/StandardModelTools_BootstrapGenerator/BootstrapGenerator), which implements a technique described in ATL-PHYS-PUB-2021-011 (https://cds.cern.ch/record/2759945). The ATLAS event number and run number of each event are used as seed to uniquely but reproducibly initialise the random number generator for each event. The last bin contains the overflow.

Data bootstraps post unfolding for the normalised differential cross-section in the phase space with at least three $b$-jets as a function of $|\eta(b_{1}^{\text{add}})|$ compared with predictions. The replicas are obtained by reweighting each observed data event by a random integer generated according to Poisson statistics, using the BootstrapGenerator software package (https://gitlab.cern.ch/atlas-physics/sm/StandardModelTools_BootstrapGenerator/BootstrapGenerator), which implements a technique described in ATL-PHYS-PUB-2021-011 (https://cds.cern.ch/record/2759945). The ATLAS event number and run number of each event are used as seed to uniquely but reproducibly initialise the random number generator for each event. The last bin contains the overflow.

Data bootstraps post unfolding for the normalised differential cross-section in the phase space with at least three $b$-jets as a function of $\Delta R(b_{1}b_{2})$ compared with predictions. The replicas are obtained by reweighting each observed data event by a random integer generated according to Poisson statistics, using the BootstrapGenerator software package (https://gitlab.cern.ch/atlas-physics/sm/StandardModelTools_BootstrapGenerator/BootstrapGenerator), which implements a technique described in ATL-PHYS-PUB-2021-011 (https://cds.cern.ch/record/2759945). The ATLAS event number and run number of each event are used as seed to uniquely but reproducibly initialise the random number generator for each event. The last bin contains the overflow.

Data bootstraps post unfolding for the normalised differential cross-section in the phase space with at least three $b$-jets as a function of $m(e\mu bb^{\text{top}})$ compared with predictions. The replicas are obtained by reweighting each observed data event by a random integer generated according to Poisson statistics, using the BootstrapGenerator software package (https://gitlab.cern.ch/atlas-physics/sm/StandardModelTools_BootstrapGenerator/BootstrapGenerator), which implements a technique described in ATL-PHYS-PUB-2021-011 (https://cds.cern.ch/record/2759945). The ATLAS event number and run number of each event are used as seed to uniquely but reproducibly initialise the random number generator for each event. The last bin contains the overflow.

Data bootstraps post unfolding for the normalised differential cross-section in the phase space with at least three $b$-jets as a function of $|\eta(l/c\text{-jet}_{1})|$ compared with predictions. The replicas are obtained by reweighting each observed data event by a random integer generated according to Poisson statistics, using the BootstrapGenerator software package (https://gitlab.cern.ch/atlas-physics/sm/StandardModelTools_BootstrapGenerator/BootstrapGenerator), which implements a technique described in ATL-PHYS-PUB-2021-011 (https://cds.cern.ch/record/2759945). The ATLAS event number and run number of each event are used as seed to uniquely but reproducibly initialise the random number generator for each event. The last bin contains the overflow.

Data bootstraps post unfolding for the normalised differential cross-section in the phase space with at least three $b$-jets as a function of $\Delta\eta_{\text{max}}^{jj}$ compared with predictions. The replicas are obtained by reweighting each observed data event by a random integer generated according to Poisson statistics, using the BootstrapGenerator software package (https://gitlab.cern.ch/atlas-physics/sm/StandardModelTools_BootstrapGenerator/BootstrapGenerator), which implements a technique described in ATL-PHYS-PUB-2021-011 (https://cds.cern.ch/record/2759945). The ATLAS event number and run number of each event are used as seed to uniquely but reproducibly initialise the random number generator for each event. The last bin contains the overflow.

Data bootstraps post unfolding for the normalised differential cross-section in the phase space with at least three $b$-jets as a function of $H_{\text{T}}^{\text{all}}$ compared with predictions. The replicas are obtained by reweighting each observed data event by a random integer generated according to Poisson statistics, using the BootstrapGenerator software package (https://gitlab.cern.ch/atlas-physics/sm/StandardModelTools_BootstrapGenerator/BootstrapGenerator), which implements a technique described in ATL-PHYS-PUB-2021-011 (https://cds.cern.ch/record/2759945). The ATLAS event number and run number of each event are used as seed to uniquely but reproducibly initialise the random number generator for each event. The last bin contains the overflow.

Data bootstraps post unfolding for the normalised differential cross-section in the phase space with at least three $b$-jets as a function of $m(e\mu b_{1}b_{2})$ compared with predictions. The replicas are obtained by reweighting each observed data event by a random integer generated according to Poisson statistics, using the BootstrapGenerator software package (https://gitlab.cern.ch/atlas-physics/sm/StandardModelTools_BootstrapGenerator/BootstrapGenerator), which implements a technique described in ATL-PHYS-PUB-2021-011 (https://cds.cern.ch/record/2759945). The ATLAS event number and run number of each event are used as seed to uniquely but reproducibly initialise the random number generator for each event. The last bin contains the overflow.

Data bootstraps post unfolding for the normalised differential cross-section in the phase space with at least three $b$-jets as a function of $|\eta(b_{1})|$ compared with predictions. The replicas are obtained by reweighting each observed data event by a random integer generated according to Poisson statistics, using the BootstrapGenerator software package (https://gitlab.cern.ch/atlas-physics/sm/StandardModelTools_BootstrapGenerator/BootstrapGenerator), which implements a technique described in ATL-PHYS-PUB-2021-011 (https://cds.cern.ch/record/2759945). The ATLAS event number and run number of each event are used as seed to uniquely but reproducibly initialise the random number generator for each event. The last bin contains the overflow.

Data bootstraps post unfolding for the normalised differential cross-section in the phase space with at least three $b$-jets as a function of $|\eta(b_{2})|$ compared with predictions. The replicas are obtained by reweighting each observed data event by a random integer generated according to Poisson statistics, using the BootstrapGenerator software package (https://gitlab.cern.ch/atlas-physics/sm/StandardModelTools_BootstrapGenerator/BootstrapGenerator), which implements a technique described in ATL-PHYS-PUB-2021-011 (https://cds.cern.ch/record/2759945). The ATLAS event number and run number of each event are used as seed to uniquely but reproducibly initialise the random number generator for each event. The last bin contains the overflow.

Data bootstraps post unfolding for the normalised differential cross-section in the phase space with at least three $b$-jets as a function of $|\eta(b_{1}^{\text{top}})|$ compared with predictions. The replicas are obtained by reweighting each observed data event by a random integer generated according to Poisson statistics, using the BootstrapGenerator software package (https://gitlab.cern.ch/atlas-physics/sm/StandardModelTools_BootstrapGenerator/BootstrapGenerator), which implements a technique described in ATL-PHYS-PUB-2021-011 (https://cds.cern.ch/record/2759945). The ATLAS event number and run number of each event are used as seed to uniquely but reproducibly initialise the random number generator for each event. The last bin contains the overflow.

Data bootstraps post unfolding for the normalised differential cross-section in the phase space with at least three $b$-jets as a function of $|\eta(b_{2}^{\text{top}})|$ compared with predictions. The replicas are obtained by reweighting each observed data event by a random integer generated according to Poisson statistics, using the BootstrapGenerator software package (https://gitlab.cern.ch/atlas-physics/sm/StandardModelTools_BootstrapGenerator/BootstrapGenerator), which implements a technique described in ATL-PHYS-PUB-2021-011 (https://cds.cern.ch/record/2759945). The ATLAS event number and run number of each event are used as seed to uniquely but reproducibly initialise the random number generator for each event. The last bin contains the overflow.

Data bootstraps post unfolding for the normalised differential cross-section in the phase space with at least four $b$-jets as a function of $p_{\text{T}}(b_{1})$ compared with predictions. The replicas are obtained by reweighting each observed data event by a random integer generated according to Poisson statistics, using the BootstrapGenerator software package (https://gitlab.cern.ch/atlas-physics/sm/StandardModelTools_BootstrapGenerator/BootstrapGenerator), which implements a technique described in ATL-PHYS-PUB-2021-011 (https://cds.cern.ch/record/2759945). The ATLAS event number and run number of each event are used as seed to uniquely but reproducibly initialise the random number generator for each event. The last bin contains the overflow.

Data bootstraps post unfolding for the normalised differential cross-section in the phase space with at least four $b$-jets as a function of $p_{\text{T}}(b_{2})$ compared with predictions. The replicas are obtained by reweighting each observed data event by a random integer generated according to Poisson statistics, using the BootstrapGenerator software package (https://gitlab.cern.ch/atlas-physics/sm/StandardModelTools_BootstrapGenerator/BootstrapGenerator), which implements a technique described in ATL-PHYS-PUB-2021-011 (https://cds.cern.ch/record/2759945). The ATLAS event number and run number of each event are used as seed to uniquely but reproducibly initialise the random number generator for each event. The last bin contains the overflow.

Data bootstraps post unfolding for the normalised differential cross-section in the phase space with at least four $b$-jets as a function of $p_{\text{T}}(b_{1}^{\text{top}})$ compared with predictions. The replicas are obtained by reweighting each observed data event by a random integer generated according to Poisson statistics, using the BootstrapGenerator software package (https://gitlab.cern.ch/atlas-physics/sm/StandardModelTools_BootstrapGenerator/BootstrapGenerator), which implements a technique described in ATL-PHYS-PUB-2021-011 (https://cds.cern.ch/record/2759945). The ATLAS event number and run number of each event are used as seed to uniquely but reproducibly initialise the random number generator for each event. The last bin contains the overflow.

Data bootstraps post unfolding for the normalised differential cross-section in the phase space with at least four $b$-jets as a function of $p_{\text{T}}(b_{2}^{\text{top}})$ compared with predictions. The replicas are obtained by reweighting each observed data event by a random integer generated according to Poisson statistics, using the BootstrapGenerator software package (https://gitlab.cern.ch/atlas-physics/sm/StandardModelTools_BootstrapGenerator/BootstrapGenerator), which implements a technique described in ATL-PHYS-PUB-2021-011 (https://cds.cern.ch/record/2759945). The ATLAS event number and run number of each event are used as seed to uniquely but reproducibly initialise the random number generator for each event. The last bin contains the overflow.

Data bootstraps post unfolding for the normalised differential cross-section in the phase space with at least four $b$-jets as a function of $p_{\text{T}}(b_{3})$ compared with predictions. The replicas are obtained by reweighting each observed data event by a random integer generated according to Poisson statistics, using the BootstrapGenerator software package (https://gitlab.cern.ch/atlas-physics/sm/StandardModelTools_BootstrapGenerator/BootstrapGenerator), which implements a technique described in ATL-PHYS-PUB-2021-011 (https://cds.cern.ch/record/2759945). The ATLAS event number and run number of each event are used as seed to uniquely but reproducibly initialise the random number generator for each event. The last bin contains the overflow.

Data bootstraps post unfolding for the normalised differential cross-section in the phase space with at least four $b$-jets as a function of $p_{\text{T}}(b_{4})$ compared with predictions. The replicas are obtained by reweighting each observed data event by a random integer generated according to Poisson statistics, using the BootstrapGenerator software package (https://gitlab.cern.ch/atlas-physics/sm/StandardModelTools_BootstrapGenerator/BootstrapGenerator), which implements a technique described in ATL-PHYS-PUB-2021-011 (https://cds.cern.ch/record/2759945). The ATLAS event number and run number of each event are used as seed to uniquely but reproducibly initialise the random number generator for each event. The last bin contains the overflow.

Data bootstraps post unfolding for the normalised differential cross-section in the phase space with at least four $b$-jets as a function of $p_{\text{T}}(b_{1}^{\text{add}})$ compared with predictions. The replicas are obtained by reweighting each observed data event by a random integer generated according to Poisson statistics, using the BootstrapGenerator software package (https://gitlab.cern.ch/atlas-physics/sm/StandardModelTools_BootstrapGenerator/BootstrapGenerator), which implements a technique described in ATL-PHYS-PUB-2021-011 (https://cds.cern.ch/record/2759945). The ATLAS event number and run number of each event are used as seed to uniquely but reproducibly initialise the random number generator for each event. The last bin contains the overflow.

Data bootstraps post unfolding for the normalised differential cross-section in the phase space with at least four $b$-jets as a function of $p_{\text{T}}(b_{2}^{\text{add}})$ compared with predictions. The replicas are obtained by reweighting each observed data event by a random integer generated according to Poisson statistics, using the BootstrapGenerator software package (https://gitlab.cern.ch/atlas-physics/sm/StandardModelTools_BootstrapGenerator/BootstrapGenerator), which implements a technique described in ATL-PHYS-PUB-2021-011 (https://cds.cern.ch/record/2759945). The ATLAS event number and run number of each event are used as seed to uniquely but reproducibly initialise the random number generator for each event. The last bin contains the overflow.

Data bootstraps post unfolding for the normalised differential cross-section in the phase space with at least four $b$-jets as a function of $m(b_{1}b_{2})$ compared with predictions. The replicas are obtained by reweighting each observed data event by a random integer generated according to Poisson statistics, using the BootstrapGenerator software package (https://gitlab.cern.ch/atlas-physics/sm/StandardModelTools_BootstrapGenerator/BootstrapGenerator), which implements a technique described in ATL-PHYS-PUB-2021-011 (https://cds.cern.ch/record/2759945). The ATLAS event number and run number of each event are used as seed to uniquely but reproducibly initialise the random number generator for each event. The last bin contains the overflow.

Data bootstraps post unfolding for the normalised differential cross-section in the phase space with at least four $b$-jets as a function of $p_{\text{T}}(b_{1}b_{2})$ compared with predictions. The replicas are obtained by reweighting each observed data event by a random integer generated according to Poisson statistics, using the BootstrapGenerator software package (https://gitlab.cern.ch/atlas-physics/sm/StandardModelTools_BootstrapGenerator/BootstrapGenerator), which implements a technique described in ATL-PHYS-PUB-2021-011 (https://cds.cern.ch/record/2759945). The ATLAS event number and run number of each event are used as seed to uniquely but reproducibly initialise the random number generator for each event. The last bin contains the overflow.

Data bootstraps post unfolding for the normalised differential cross-section in the phase space with at least four $b$-jets as a function of $m(bb^{\text{top}})$ compared with predictions. The replicas are obtained by reweighting each observed data event by a random integer generated according to Poisson statistics, using the BootstrapGenerator software package (https://gitlab.cern.ch/atlas-physics/sm/StandardModelTools_BootstrapGenerator/BootstrapGenerator), which implements a technique described in ATL-PHYS-PUB-2021-011 (https://cds.cern.ch/record/2759945). The ATLAS event number and run number of each event are used as seed to uniquely but reproducibly initialise the random number generator for each event. The last bin contains the overflow.

Data bootstraps post unfolding for the normalised differential cross-section in the phase space with at least four $b$-jets as a function of $p_{\text{T}}(bb^{\text{top}})$ compared with predictions. The replicas are obtained by reweighting each observed data event by a random integer generated according to Poisson statistics, using the BootstrapGenerator software package (https://gitlab.cern.ch/atlas-physics/sm/StandardModelTools_BootstrapGenerator/BootstrapGenerator), which implements a technique described in ATL-PHYS-PUB-2021-011 (https://cds.cern.ch/record/2759945). The ATLAS event number and run number of each event are used as seed to uniquely but reproducibly initialise the random number generator for each event. The last bin contains the overflow.

Data bootstraps post unfolding for the normalised differential cross-section in the phase space with at least four $b$-jets as a function of $H_{\text{T}}^{\text{all}}$ compared with predictions. The replicas are obtained by reweighting each observed data event by a random integer generated according to Poisson statistics, using the BootstrapGenerator software package (https://gitlab.cern.ch/atlas-physics/sm/StandardModelTools_BootstrapGenerator/BootstrapGenerator), which implements a technique described in ATL-PHYS-PUB-2021-011 (https://cds.cern.ch/record/2759945). The ATLAS event number and run number of each event are used as seed to uniquely but reproducibly initialise the random number generator for each event. The last bin contains the overflow.

Data bootstraps post unfolding for the normalised differential cross-section in the phase space with at least four $b$-jets as a function of $m(e\mu b_{1}b_{2})$ compared with predictions. The replicas are obtained by reweighting each observed data event by a random integer generated according to Poisson statistics, using the BootstrapGenerator software package (https://gitlab.cern.ch/atlas-physics/sm/StandardModelTools_BootstrapGenerator/BootstrapGenerator), which implements a technique described in ATL-PHYS-PUB-2021-011 (https://cds.cern.ch/record/2759945). The ATLAS event number and run number of each event are used as seed to uniquely but reproducibly initialise the random number generator for each event. The last bin contains the overflow.

Data bootstraps post unfolding for the normalised differential cross-section in the phase space with at least four $b$-jets as a function of $m(e\mu bb^{\text{top}})$ compared with predictions. The replicas are obtained by reweighting each observed data event by a random integer generated according to Poisson statistics, using the BootstrapGenerator software package (https://gitlab.cern.ch/atlas-physics/sm/StandardModelTools_BootstrapGenerator/BootstrapGenerator), which implements a technique described in ATL-PHYS-PUB-2021-011 (https://cds.cern.ch/record/2759945). The ATLAS event number and run number of each event are used as seed to uniquely but reproducibly initialise the random number generator for each event. The last bin contains the overflow.

Data bootstraps post unfolding for the normalised differential cross-section in the phase space with at least four $b$-jets as a function of $H_{\text{T}}^{\text{had}}$ compared with predictions. The replicas are obtained by reweighting each observed data event by a random integer generated according to Poisson statistics, using the BootstrapGenerator software package (https://gitlab.cern.ch/atlas-physics/sm/StandardModelTools_BootstrapGenerator/BootstrapGenerator), which implements a technique described in ATL-PHYS-PUB-2021-011 (https://cds.cern.ch/record/2759945). The ATLAS event number and run number of each event are used as seed to uniquely but reproducibly initialise the random number generator for each event. The last bin contains the overflow.

Data bootstraps post unfolding for the normalised differential cross-section in the phase space with at least four $b$-jets as a function of $\text{min}\Delta R(bb)$ compared with predictions. The replicas are obtained by reweighting each observed data event by a random integer generated according to Poisson statistics, using the BootstrapGenerator software package (https://gitlab.cern.ch/atlas-physics/sm/StandardModelTools_BootstrapGenerator/BootstrapGenerator), which implements a technique described in ATL-PHYS-PUB-2021-011 (https://cds.cern.ch/record/2759945). The ATLAS event number and run number of each event are used as seed to uniquely but reproducibly initialise the random number generator for each event. The last bin contains the overflow.

Data bootstraps post unfolding for the normalised differential cross-section in the phase space with at least four $b$-jets as a function of $\Delta R(e\mu bb^{\text{top}}, b_{1}^{\text{add}})$ compared with predictions. The replicas are obtained by reweighting each observed data event by a random integer generated according to Poisson statistics, using the BootstrapGenerator software package (https://gitlab.cern.ch/atlas-physics/sm/StandardModelTools_BootstrapGenerator/BootstrapGenerator), which implements a technique described in ATL-PHYS-PUB-2021-011 (https://cds.cern.ch/record/2759945). The ATLAS event number and run number of each event are used as seed to uniquely but reproducibly initialise the random number generator for each event. The last bin contains the overflow.

Data bootstraps post unfolding for the normalised differential cross-section in the phase space with at least four $b$-jets as a function of $\Delta R_{\text{avg}}^{bb}$ compared with predictions. The replicas are obtained by reweighting each observed data event by a random integer generated according to Poisson statistics, using the BootstrapGenerator software package (https://gitlab.cern.ch/atlas-physics/sm/StandardModelTools_BootstrapGenerator/BootstrapGenerator), which implements a technique described in ATL-PHYS-PUB-2021-011 (https://cds.cern.ch/record/2759945). The ATLAS event number and run number of each event are used as seed to uniquely but reproducibly initialise the random number generator for each event. The last bin contains the overflow.

Data bootstraps post unfolding for the normalised differential cross-section in the phase space with at least four $b$-jets as a function of $\Delta\eta_{\text{max}}^{jj}$ compared with predictions. The replicas are obtained by reweighting each observed data event by a random integer generated according to Poisson statistics, using the BootstrapGenerator software package (https://gitlab.cern.ch/atlas-physics/sm/StandardModelTools_BootstrapGenerator/BootstrapGenerator), which implements a technique described in ATL-PHYS-PUB-2021-011 (https://cds.cern.ch/record/2759945). The ATLAS event number and run number of each event are used as seed to uniquely but reproducibly initialise the random number generator for each event. The last bin contains the overflow.

Data bootstraps post unfolding for the normalised differential cross-section in the phase space with at least four $b$-jets as a function of the number of $l/c$-jets compared with predictions. The replicas are obtained by reweighting each observed data event by a random integer generated according to Poisson statistics, using the BootstrapGenerator software package (https://gitlab.cern.ch/atlas-physics/sm/StandardModelTools_BootstrapGenerator/BootstrapGenerator), which implements a technique described in ATL-PHYS-PUB-2021-011 (https://cds.cern.ch/record/2759945). The ATLAS event number and run number of each event are used as seed to uniquely but reproducibly initialise the random number generator for each event. The last bin contains the overflow.

Data bootstraps post unfolding for the normalised differential cross-section in the phase space with at least four $b$-jets and at least one $l/c$-jet as a function of $p_{\text{T}}(l/c\text{-jet}_{1})$ compared with predictions. The replicas are obtained by reweighting each observed data event by a random integer generated according to Poisson statistics, using the BootstrapGenerator software package (https://gitlab.cern.ch/atlas-physics/sm/StandardModelTools_BootstrapGenerator/BootstrapGenerator), which implements a technique described in ATL-PHYS-PUB-2021-011 (https://cds.cern.ch/record/2759945). The ATLAS event number and run number of each event are used as seed to uniquely but reproducibly initialise the random number generator for each event. The last bin contains the overflow.

Data bootstraps post unfolding for the normalised differential cross-section in the phase space with at least four $b$-jets and at least one $l/c$-jet as a function of $|\eta(l/c\text{-jet}_{1})|$ compared with predictions. The replicas are obtained by reweighting each observed data event by a random integer generated according to Poisson statistics, using the BootstrapGenerator software package (https://gitlab.cern.ch/atlas-physics/sm/StandardModelTools_BootstrapGenerator/BootstrapGenerator), which implements a technique described in ATL-PHYS-PUB-2021-011 (https://cds.cern.ch/record/2759945). The ATLAS event number and run number of each event are used as seed to uniquely but reproducibly initialise the random number generator for each event. The last bin contains the overflow.

Data bootstraps post unfolding for the normalised differential cross-section in the phase space with at least four $b$-jets and at least one $l/c$-jet as a function of $\Delta R(e\mu bb^{\text{top}}, l/c\text{-jet}_{1})$ compared with predictions. The replicas are obtained by reweighting each observed data event by a random integer generated according to Poisson statistics, using the BootstrapGenerator software package (https://gitlab.cern.ch/atlas-physics/sm/StandardModelTools_BootstrapGenerator/BootstrapGenerator), which implements a technique described in ATL-PHYS-PUB-2021-011 (https://cds.cern.ch/record/2759945). The ATLAS event number and run number of each event are used as seed to uniquely but reproducibly initialise the random number generator for each event. The last bin contains the overflow.

Data bootstraps post unfolding for the normalised differential cross-section in the phase space with at least four $b$-jets and at least one $l/c$-jet as a function of $p_{\text{T}}(l/c\text{-jet}_{1}) - p_{\text{T}}(b_{1}^{\text{add}})$ compared with predictions. The replicas are obtained by reweighting each observed data event by a random integer generated according to Poisson statistics, using the BootstrapGenerator software package (https://gitlab.cern.ch/atlas-physics/sm/StandardModelTools_BootstrapGenerator/BootstrapGenerator), which implements a technique described in ATL-PHYS-PUB-2021-011 (https://cds.cern.ch/record/2759945). The ATLAS event number and run number of each event are used as seed to uniquely but reproducibly initialise the random number generator for each event. The last bin contains the overflow.

Data bootstraps post unfolding for the normalised differential cross-section in the phase space with at least four $b$-jets as a function of $|\eta(b_{1})|$ compared with predictions. The replicas are obtained by reweighting each observed data event by a random integer generated according to Poisson statistics, using the BootstrapGenerator software package (https://gitlab.cern.ch/atlas-physics/sm/StandardModelTools_BootstrapGenerator/BootstrapGenerator), which implements a technique described in ATL-PHYS-PUB-2021-011 (https://cds.cern.ch/record/2759945). The ATLAS event number and run number of each event are used as seed to uniquely but reproducibly initialise the random number generator for each event. The last bin contains the overflow.

Data bootstraps post unfolding for the normalised differential cross-section in the phase space with at least four $b$-jets as a function of $|\eta(b_{2})|$ compared with predictions. The replicas are obtained by reweighting each observed data event by a random integer generated according to Poisson statistics, using the BootstrapGenerator software package (https://gitlab.cern.ch/atlas-physics/sm/StandardModelTools_BootstrapGenerator/BootstrapGenerator), which implements a technique described in ATL-PHYS-PUB-2021-011 (https://cds.cern.ch/record/2759945). The ATLAS event number and run number of each event are used as seed to uniquely but reproducibly initialise the random number generator for each event. The last bin contains the overflow.

Data bootstraps post unfolding for the normalised differential cross-section in the phase space with at least four $b$-jets as a function of $|\eta(b_{1}^{\text{top}})|$ compared with predictions. The replicas are obtained by reweighting each observed data event by a random integer generated according to Poisson statistics, using the BootstrapGenerator software package (https://gitlab.cern.ch/atlas-physics/sm/StandardModelTools_BootstrapGenerator/BootstrapGenerator), which implements a technique described in ATL-PHYS-PUB-2021-011 (https://cds.cern.ch/record/2759945). The ATLAS event number and run number of each event are used as seed to uniquely but reproducibly initialise the random number generator for each event. The last bin contains the overflow.

Data bootstraps post unfolding for the normalised differential cross-section in the phase space with at leastfour $b$-jets as a function of $|\eta(b_{2}^{\text{top}})|$ compared with predictions. The replicas are obtained by reweighting each observed data event by a random integer generated according to Poisson statistics, using the BootstrapGenerator software package (https://gitlab.cern.ch/atlas-physics/sm/StandardModelTools_BootstrapGenerator/BootstrapGenerator), which implements a technique described in ATL-PHYS-PUB-2021-011 (https://cds.cern.ch/record/2759945). The ATLAS event number and run number of each event are used as seed to uniquely but reproducibly initialise the random number generator for each event. The last bin contains the overflow.

Data bootstraps post unfolding for the normalised differential cross-section in the phase space with at least four $b$-jets as a function of $|\eta(b_{3})|$ compared with predictions. The replicas are obtained by reweighting each observed data event by a random integer generated according to Poisson statistics, using the BootstrapGenerator software package (https://gitlab.cern.ch/atlas-physics/sm/StandardModelTools_BootstrapGenerator/BootstrapGenerator), which implements a technique described in ATL-PHYS-PUB-2021-011 (https://cds.cern.ch/record/2759945). The ATLAS event number and run number of each event are used as seed to uniquely but reproducibly initialise the random number generator for each event. The last bin contains the overflow.

Data bootstraps post unfolding for the normalised differential cross-section in the phase space with at least four $b$-jets as a function of $|\eta(b_{4})|$ compared with predictions. The replicas are obtained by reweighting each observed data event by a random integer generated according to Poisson statistics, using the BootstrapGenerator software package (https://gitlab.cern.ch/atlas-physics/sm/StandardModelTools_BootstrapGenerator/BootstrapGenerator), which implements a technique described in ATL-PHYS-PUB-2021-011 (https://cds.cern.ch/record/2759945). The ATLAS event number and run number of each event are used as seed to uniquely but reproducibly initialise the random number generator for each event. The last bin contains the overflow.

Data bootstraps post unfolding for the normalised differential cross-section in the phase space with at least four $b$-jets as a function of $|\eta(b_{1}^{\text{add}})|$ compared with predictions. The replicas are obtained by reweighting each observed data event by a random integer generated according to Poisson statistics, using the BootstrapGenerator software package (https://gitlab.cern.ch/atlas-physics/sm/StandardModelTools_BootstrapGenerator/BootstrapGenerator), which implements a technique described in ATL-PHYS-PUB-2021-011 (https://cds.cern.ch/record/2759945). The ATLAS event number and run number of each event are used as seed to uniquely but reproducibly initialise the random number generator for each event. The last bin contains the overflow.

Data bootstraps post unfolding for the normalised differential cross-section in the phase space with at least four $b$-jets as a function of $|\eta(b_{2}^{\text{add}})|$ compared with predictions. The replicas are obtained by reweighting each observed data event by a random integer generated according to Poisson statistics, using the BootstrapGenerator software package (https://gitlab.cern.ch/atlas-physics/sm/StandardModelTools_BootstrapGenerator/BootstrapGenerator), which implements a technique described in ATL-PHYS-PUB-2021-011 (https://cds.cern.ch/record/2759945). The ATLAS event number and run number of each event are used as seed to uniquely but reproducibly initialise the random number generator for each event. The last bin contains the overflow.

The measured normalised differential cross-section as a function of $N_{b-\text{jets}}$ in the $e\mu+\geq2b$ phase space. The overflow is included in the last bin.

The measured normalised differential cross-section as a function of $H_{\text{T}}^{\text{had}}$ in the $e\mu+\geq3b$ phase space. The overflow is included in the last bin.

The measured normalised differential cross-section as a function of $H_{\text{T}}^{\text{all}}$ in the $e\mu+\geq3b$ phase space. The overflow is included in the last bin.

The measured normalised differential cross-section as a function of $\Delta R_{\text{avg}}^{bb}$ in the $e\mu+\geq3b$ phase space. The overflow is included in the last bin.

The measured normalised differential cross-section as a function of $\Delta\eta_{\text{max}}^{jj}$ in the $e\mu+\geq3b$ phase space. The overflow is included in the last bin.

The measured normalised differential cross-section as a function of $p_{\text{T}}(b_{1})$ in the $e\mu+\geq3b$ phase space. The overflow is included in the last bin.

The measured normalised differential cross-section as a function of $p_{\text{T}}(b_{1}^{\text{top}})$ in the $e\mu+\geq3b$ phase space. The overflow is included in the last bin.

The measured normalised differential cross-section as a function of $p_{\text{T}}(b_{2})$ in the $e\mu+\geq3b$ phase space. The overflow is included in the last bin.

The measured normalised differential cross-section as a function of $p_{\text{T}}(b_{2}^{\text{top}})$ in the $e\mu+\geq3b$ phase space. The overflow is included in the last bin.

The measured normalised differential cross-section as a function of $p_{\text{T}}(b_{3})$ in the $e\mu+\geq3b$ phase space. The overflow is included in the last bin.

The measured normalised differential cross-section as a function of $p_{\text{T}}(b_{1}^{\text{add}})$ in the $e\mu+\geq3b$ phase space. The overflow is included in the last bin.

The measured normalised differential cross-section as a function of $|\eta(b_{1})|$ in the $e\mu+\geq3b$ phase space. The overflow is included in the last bin.

The measured normalised differential cross-section as a function of $|\eta(b_{1}^{\text{top}})|$ in the $e\mu+\geq3b$ phase space. The overflow is included in the last bin.

The measured normalised differential cross-section as a function of $|\eta(b_{2})|$ in the $e\mu+\geq3b$ phase space. The overflow is included in the last bin.

The measured normalised differential cross-section as a function of $|\eta(b_{2}^{\text{top}})|$ in the $e\mu+\geq3b$ phase space. The overflow is included in the last bin.

The measured normalised differential cross-section as a function of $|\eta(b_{3})|$ in the $e\mu+\geq3b$ phase space. The overflow is included in the last bin.

The measured normalised differential cross-section as a function of $|\eta(b_{1}^{\text{add}})|$ in the $e\mu+\geq3b$ phase space. The overflow is included in the last bin.

The measured normalised differential cross-section as a function of $m(b_{1}b_{2})$ in the $e\mu+\geq3b$ phase space. The overflow is included in the last bin.

The measured normalised differential cross-section as a function of $p_{\text{T}}(b_{1}b_{2})$ in the $e\mu+\geq3b$ phase space. The overflow is included in the last bin.

The measured normalised differential cross-section as a function of $m(bb^{\text{top}})$ in the $e\mu+\geq3b$ phase space. The overflow is included in the last bin.

The measured normalised differential cross-section as a function of $p_{\text{T}}(bb^{\text{top}})$ in the $e\mu+\geq3b$ phase space. The overflow is included in the last bin.

The measured normalised differential cross-section as a function of $m(e\mu b_{1}b_{2})$ in the $e\mu+\geq3b$ phase space. The overflow is included in the last bin.

The measured normalised differential cross-section as a function of $m(e\mu bb^{\text{top}})$ in the $e\mu+\geq3b$ phase space. The overflow is included in the last bin.

The measured normalised differential cross-section as a function of $\Delta R(b_{1}b_{2})$ in the $e\mu+\geq3b$ phase space. The overflow is included in the last bin.

The measured normalised differential cross-section as a function of $N_{l/c-\text{jets}}$ in the $e\mu+\geq3b$ phase space. The overflow is included in the last bin.

The measured normalised differential cross-section as a function of $\Delta R(e\mu b_{1}b_{2},b_{3})$ in the $e\mu+\geq3b$ phase space. The overflow is included in the last bin.

The measured normalised differential cross-section as a function of $\Delta R(e\mu bb^{\text{top}}, b_{1}^{\text{add}})$ in the $e\mu+\geq3b$ phase space. The overflow is included in the last bin.

The measured normalised differential cross-section as a function of $\Delta R(e\mu bb^{\text{top}},l/c-\text{jet})$ in the $e\mu+\geq3b+\geq1l/c-\text{jet}$ phase space. The overflow is included in the last bin.

The measured normalised differential cross-section as a function of $p_{\text{T}}(l/c\text{-jet}_{1}) - p_{\text{T}}(b_{1}^{\text{add}})$ in the $e\mu+\geq3b+\geq1l/c-\text{jet}$ phase space. The overflow is included in the last bin.

The measured normalised differential cross-section as a function of $|\eta(l/c\text{-jet}_{1})|$ in the $e\mu+\geq3b+\geq1l/c-\text{jet}$ phase space. The overflow is included in the last bin.

The measured normalised differential cross-section as a function of $p_{\text{T}}(l/c\text{-jet}_{1})$ in the $e\mu+\geq3b+\geq1l/c-\text{jet}$ phase space. The overflow is included in the last bin.

The measured normalised differential cross-section as a function of $H_{\text{T}}^{\text{had}}$ in the $e\mu+\geq4b$ phase space. The overflow is included in the last bin.

The measured normalised differential cross-section as a function of $H_{\text{T}}^{\text{all}}$ in the $e\mu+\geq4b$ phase space. The overflow is included in the last bin.

The measured normalised differential cross-section as a function of $\Delta R_{\text{avg}}^{bb}$ in the $e\mu+\geq4b$ phase space. The overflow is included in the last bin.

The measured normalised differential cross-section as a function of $\Delta\eta_{\text{max}}^{jj}$ in the $e\mu+\geq4b$ phase space. The overflow is included in the last bin.

The measured normalised differential cross-section as a function of $p_{\text{T}}(b_{1})$ in the $e\mu+\geq4b$ phase space. The overflow is included in the last bin.

The measured normalised differential cross-section as a function of $p_{\text{T}}(b_{1}^{\text{top}})$ in the $e\mu+\geq4b$ phase space. The overflow is included in the last bin.

The measured normalised differential cross-section as a function of $p_{\text{T}}(b_{2})$ in the $e\mu+\geq4b$ phase space. The overflow is included in the last bin.

The measured normalised differential cross-section as a function of $p_{\text{T}}(b_{2}^{\text{top}})$ in the $e\mu+\geq4b$ phase space. The overflow is included in the last bin.

The measured normalised differential cross-section as a function of $p_{\text{T}}(b_{3})$ in the $e\mu+\geq4b$ phase space. The overflow is included in the last bin.

The measured normalised differential cross-section as a function of $p_{\text{T}}(b_{1}^{\text{add}})$ in the $e\mu+\geq4b$ phase space. The overflow is included in the last bin.

The measured normalised differential cross-section as a function of $p_{\text{T}}(b_{4})$ in the $e\mu+\geq4b$ phase space. The overflow is included in the last bin.

The measured normalised differential cross-section as a function of $p_{\text{T}}(b_{2}^{\text{add}})$ in the $e\mu+\geq4b$ phase space. The overflow is included in the last bin.

The measured normalised differential cross-section as a function of $|\eta(b_{1})|$ in the $e\mu+\geq4b$ phase space. The overflow is included in the last bin.

The measured normalised differential cross-section as a function of $|\eta(b_{1}^{\text{top}})|$ in the $e\mu+\geq4b$ phase space. The overflow is included in the last bin.

The measured normalised differential cross-section as a function of $|\eta(b_{2})|$ in the $e\mu+\geq4b$ phase space. The overflow is included in the last bin.

The measured normalised differential cross-section as a function of $|\eta(b_{2}^{\text{top}})|$ in the $e\mu+\geq4b$ phase space. The overflow is included in the last bin.

The measured normalised differential cross-section as a function of $|\eta(b_{3})|$ in the $e\mu+\geq4b$ phase space. The overflow is included in the last bin.

The measured normalised differential cross-section as a function of $|\eta(b_{1}^{\text{add}})|$ in the $e\mu+\geq4b$ phase space. The overflow is included in the last bin.

The measured normalised differential cross-section as a function of $|\eta(b_{4})|$ in the $e\mu+\geq4b$ phase space. The overflow is included in the last bin.

The measured normalised differential cross-section as a function of $|\eta(b_{2}^{\text{add}})|$ in the $e\mu+\geq4b$ phase space. The overflow is included in the last bin.

The measured normalised differential cross-section as a function of $m(b_{1}b_{2})$ in the $e\mu+\geq4b$ phase space. The overflow is included in the last bin.

The measured normalised differential cross-section as a function of $p_{\text{T}}(b_{1}b_{2})$ in the $e\mu+\geq4b$ phase space. The overflow is included in the last bin.

The measured normalised differential cross-section as a function of $m(bb^{\text{top}})$ in the $e\mu+\geq4b$ phase space. The overflow is included in the last bin.

The measured normalised differential cross-section as a function of $p_{\text{T}}(bb^{\text{top}})$ in the $e\mu+\geq4b$ phase space. The overflow is included in the last bin.

The measured normalised differential cross-section as a function of $m(e\mu b_{1}b_{2})$ in the $e\mu+\geq4b$ phase space. The overflow is included in the last bin.

The measured normalised differential cross-section as a function of $m(e\mu bb^{\text{top}})$ in the $e\mu+\geq4b$ phase space. The overflow is included in the last bin.

The measured normalised differential cross-section as a function of $m(bb^{\text{min}\Delta R})$ in the $e\mu+\geq4b$ phase space. The overflow is included in the last bin.

The measured normalised differential cross-section as a function of $p_{\text{T}}(bb^{\text{min}\Delta R})$ in the $e\mu+\geq4b$ phase space. The overflow is included in the last bin.

The measured normalised differential cross-section as a function of $m(bb^{\text{add}})$ in the $e\mu+\geq4b$ phase space. The overflow is included in the last bin.

The measured normalised differential cross-section as a function of $p_{\text{T}}(bb^{\text{add}})$ in the $e\mu+\geq4b$ phase space. The overflow is included in the last bin.

The measured normalised differential cross-section as a function of $\text{min}\Delta R(bb)$ in the $e\mu+\geq4b$ phase space. The overflow is included in the last bin.

The measured normalised differential cross-section as a function of $\Delta R(b_{1}b_{2})$ in the $e\mu+\geq4b$ phase space. The overflow is included in the last bin.

The measured normalised differential cross-section as a function of $N_{l/c-\text{jets}}$ in the $e\mu+\geq4b$ phase space. The overflow is included in the last bin.

The measured normalised differential cross-section as a function of $\Delta R(e\mu b_{1}b_{2},b_{3})$ in the $e\mu+\geq4b$ phase space. The overflow is included in the last bin.

The measured normalised differential cross-section as a function of $\Delta R(e\mu bb^{\text{top}}, b_{1}^{\text{add}})$ in the $e\mu+\geq4b$ phase space. The overflow is included in the last bin.

The measured normalised differential cross-section as a function of $\Delta R(e\mu bb^{\text{top}}, l/c\text{-jet}_{1})$ in the $e\mu+\geq4b+\geq1l/c-\text{jet}$ phase space. The overflow is included in the last bin.

The measured normalised differential cross-section as a function of $p_{\text{T}}(l/c\text{-jet}_{1}) - p_{\text{T}}(b_{1}^{\text{add}})$ in the $e\mu+\geq4b+\geq1l/c-\text{jet}$ phase space. The overflow is included in the last bin.

The measured normalised differential cross-section as a function of $|\eta(l/c\text{-jet}_{1})|$ in the $e\mu+\geq4b+\geq1l/c-\text{jet}$ phase space. The overflow is included in the last bin.

The measured normalised differential cross-section as a function of $p_{\text{T}}(l/c\text{-jet}_{1})$ in the $e\mu+\geq4b+\geq1l/c-\text{jet}$ phase space. The overflow is included in the last bin.

The correlation matrix for the measured normalised differential cross-section in terms of $N_{b-\text{jets}}$ in the phase space with at least two b-jets. The overflow is included in the last bin.

The correlation matrix for the measured normalised differential cross-section in terms of $N_{b-\text{jets}}$ in the phase space with at least three b-jets. The overflow is included in the last bin.

The correlation matrix for the measured normalised differential cross-section in terms of $H_{\text{T}}^{\text{had}}$ in the phase space with at least three b-jets. The overflow is included in the last bin.

The correlation matrix for the measured normalised differential cross-section in terms of $H_{\text{T}}^{\text{all}}$ in the phase space with at least three b-jets. The overflow is included in the last bin.

The correlation matrix for the measured normalised differential cross-section in terms of $\Delta R_{\text{avg}}^{bb}$ in the phase space with at least three b-jets. The overflow is included in the last bin.

The correlation matrix for the measured normalised differential cross-section in terms of $\Delta\eta_{\text{max}}^{jj}$ in the phase space with at least three b-jets. The overflow is included in the last bin.

The correlation matrix for the measured normalised differential cross-section in terms of $p_{\text{T}}(b_{1})$ in the phase space with at least three b-jets. The overflow is included in the last bin.

The correlation matrix for the measured normalised differential cross-section in terms of $p_{\text{T}}(b_{1}^{\text{top}})$ in the phase space with at least three b-jets. The overflow is included in the last bin.

The correlation matrix for the measured normalised differential cross-section in terms of $p_{\text{T}}(b_{2})$ in the phase space with at least three b-jets. The overflow is included in the last bin.

The correlation matrix for the measured normalised differential cross-section in terms of $p_{\text{T}}(b_{2}^{\text{top}})$ in the phase space with at least three b-jets. The overflow is included in the last bin.

The correlation matrix for the measured normalised differential cross-section in terms of $p_{\text{T}}(b_{3})$ in the phase space with at least three b-jets. The overflow is included in the last bin.

The correlation matrix for the measured normalised differential cross-section in terms of $p_{\text{T}}(b_{1}^{\text{add}})$ in the phase space with at least three b-jets. The overflow is included in the last bin.

The correlation matrix for the measured normalised differential cross-section in terms of $|\eta(b_{1})|$ in the phase space with at least three b-jets. The overflow is included in the last bin.

The correlation matrix for the measured normalised differential cross-section in terms of $|\eta(b_{1}^{\text{top}})|$ in the phase space with at least three b-jets. The overflow is included in the last bin.

The correlation matrix for the measured normalised differential cross-section in terms of $|\eta(b_{2}^{\text{top}})|$ in the phase space with at least three b-jets. The overflow is included in the last bin.

The correlation matrix for the measured normalised differential cross-section in terms of $|\eta(b_{2}^{\text{top}})|$ in the phase space with at least three b-jets. The overflow is included in the last bin.

The correlation matrix for the measured normalised differential cross-section in terms of $|\eta(b_{3})|$ in the phase space with at least three b-jets. The overflow is included in the last bin.

The correlation matrix for the measured normalised differential cross-section in terms of $|\eta(b_{1}^{\text{add}})|$ in the phase space with at least three b-jets. The overflow is included in the last bin.

The correlation matrix for the measured normalised differential cross-section in terms of $m(b_{1}b_{2})$ in the phase space with at least three b-jets. The overflow is included in the last bin.

The correlation matrix for the measured normalised differential cross-section in terms of $p_{\text{T}}(b_{1}b_{2})$ in the phase space with at least three b-jets. The overflow is included in the last bin.

The correlation matrix for the measured normalised differential cross-section in terms of $m(bb^{\text{top}})$ in the phase space with at least three b-jets. The overflow is included in the last bin.

The correlation matrix for the measured normalised differential cross-section in terms of $p_{\text{T}}(bb^{\text{top}})$ in the phase space with at least three b-jets. The overflow is included in the last bin.

The correlation matrix for the measured normalised differential cross-section in terms of $m(e\mu b_{1}b_{2})$ in the phase space with at least three b-jets. The overflow is included in the last bin.

The correlation matrix for the measured normalised differential cross-section in terms of $m(e\mu bb^{\text{top}})$ in the phase space with at least three b-jets. The overflow is included in the last bin.

The correlation matrix for the measured normalised differential cross-section in terms of $\Delta R(b_{1}b_{2})$ in the phase space with at least three b-jets. The overflow is included in the last bin.

The correlation matrix for the measured normalised differential cross-section in terms of $N_{l/c-\text{jets}}$ in the phase space with at least three b-jets. The overflow is included in the last bin.

The correlation matrix for the measured normalised differential cross-section in terms of $\Delta R(e\mu b_{1}b_{2},b_{3})$ in the phase space with at least three b-jets. The overflow is included in the last bin.

The correlation matrix for the measured normalised differential cross-section in terms of $\Delta R(e\mu bb^{\text{top}}, b_{1}^{\text{add}})$ in the phase space with at least three b-jets. The overflow is included in the last bin.

The correlation matrix for the measured normalised differential cross-section in terms of $\Delta R(e\mu bb^{\text{top}},l/c-\text{jet})$ in the phase space with at least three b-jets and at least one $l/c$-jet. The overflow is included in the last bin.

The correlation matrix for the measured normalised differential cross-section in terms of $p_{\text{T}}(l/c\text{-jet}_{1}) - p_{\text{T}}(b_{1}^{\text{add}})$ in the phase space with at least three b-jets and at least one $l/c$-jet. The overflow is included in the last bin.

The correlation matrix for the measured normalised differential cross-section in terms of $|\eta(l/c\text{-jet}_{1})|$ in the phase space with at least three b-jets and at least one $l/c$-jet. The overflow is included in the last bin.

The correlation matrix for the measured normalised differential cross-section in terms of $p_{\text{T}}(l/c\text{-jet}_{1})$ in the phase space with at least three b-jets and at least one $l/c$-jet. The overflow is included in the last bin.

The correlation matrix for the measured normalised differential cross-section in terms of $H_{\text{T}}^{\text{had}}$ in the phase space with at least four b-jets. The overflow is included in the last bin.

The correlation matrix for the measured normalised differential cross-section in terms of $H_{\text{T}}^{\text{all}}$ in the phase space with at least four b-jets. The overflow is included in the last bin.

The correlation matrix for the measured normalised differential cross-section in terms of $\Delta R_{\text{avg}}^{bb}$ in the phase space with at least four b-jets. The overflow is included in the last bin.

The correlation matrix for the measured normalised differential cross-section in terms of $\Delta\eta_{\text{max}}^{jj}$ in the phase space with at least four b-jets. The overflow is included in the last bin.

The correlation matrix for the measured normalised differential cross-section in terms of $p_{\text{T}}(b_{1})$ in the phase space with at least four b-jets. The overflow is included in the last bin.

The correlation matrix for the measured normalised differential cross-section in terms of $p_{\text{T}}(b_{1}^{\text{top}})$ in the phase space with at least four b-jets. The overflow is included in the last bin.

The correlation matrix for the measured normalised differential cross-section in terms of $p_{\text{T}}(b_{2})$ in the phase space with at least four b-jets. The overflow is included in the last bin.

The correlation matrix for the measured normalised differential cross-section in terms of $p_{\text{T}}(b_{2}^{\text{top}})$ in the phase space with at least four b-jets. The overflow is included in the last bin.

The correlation matrix for the measured normalised differential cross-section in terms of $p_{\text{T}}(b_{3})$ in the phase space with at least four b-jets. The overflow is included in the last bin.

The correlation matrix for the measured normalised differential cross-section in terms of $p_{\text{T}}(b_{1}^{\text{add}})$ in the phase space with at least four b-jets. The overflow is included in the last bin.

The correlation matrix for the measured normalised differential cross-section in terms of $p_{\text{T}}(b_{4})$ in the phase space with at least four b-jets. The overflow is included in the last bin.

The correlation matrix for the measured normalised differential cross-section in terms of $p_{\text{T}}(b_{2}^{\text{add}})$ in the phase space with at least four b-jets. The overflow is included in the last bin.

The correlation matrix for the measured normalised differential cross-section in terms of $|\eta(b_{1})|$ in the phase space with at least four b-jets. The overflow is included in the last bin.

The correlation matrix for the measured normalised differential cross-section in terms of $|\eta(b_{1}^{\text{top}})|$ in the phase space with at least four b-jets. The overflow is included in the last bin.

The correlation matrix for the measured normalised differential cross-section in terms of $|\eta(b_{2})|$ in the phase space with at least four b-jets. The overflow is included in the last bin.

The correlation matrix for the measured normalised differential cross-section in terms of $|\eta(b_{2}^{\text{top}})|$ in the phase space with at least four b-jets. The overflow is included in the last bin.

The correlation matrix for the measured normalised differential cross-section in terms of $|\eta(b_{3})|$ in the phase space with at least four b-jets. The overflow is included in the last bin.

The correlation matrix for the measured normalised differential cross-section in terms of $|\eta(b_{1}^{\text{add}})|$ in the phase space with at least four b-jets. The overflow is included in the last bin.

The correlation matrix for the measured normalised differential cross-section in terms of $|\eta(b_{4})|$ in the phase space with at least four b-jets. The overflow is included in the last bin.

The correlation matrix for the measured normalised differential cross-section in terms of $|\eta(b_{2}^{\text{add}})|$ in the phase space with at least four b-jets. The overflow is included in the last bin.

The correlation matrix for the measured normalised differential cross-section in terms of $m(b_{1}b_{2})$ in the phase space with at least four b-jets. The overflow is included in the last bin.

The correlation matrix for the measured normalised differential cross-section in terms of $p_{\text{T}}(b_{1}b_{2})$ in the phase space with at least four b-jets. The overflow is included in the last bin.

The correlation matrix for the measured normalised differential cross-section in terms of $m(bb^{\text{top}})$ in the phase space with at least four b-jets. The overflow is included in the last bin.

The correlation matrix for the measured normalised differential cross-section in terms of $p_{\text{T}}(bb^{\text{top}})$ in the phase space with at least four b-jets. The overflow is included in the last bin.

The correlation matrix for the measured normalised differential cross-section in terms of $m(e\mu b_{1}b_{2})$ in the phase space with at least four b-jets. The overflow is included in the last bin.

The correlation matrix for the measured normalised differential cross-section in terms of $m(e\mu bb^{\text{top}})$ in the phase space with at least four b-jets. The overflow is included in the last bin.

The correlation matrix for the measured normalised differential cross-section in terms of $m(bb^{\text{min}\Delta R})$ in the phase space with at least four b-jets. The overflow is included in the last bin.

The correlation matrix for the measured normalised differential cross-section in terms of $p_{\text{T}}(bb^{\text{min}\Delta R})$ in the phase space with at least four b-jets. The overflow is included in the last bin.

The correlation matrix for the measured normalised differential cross-section in terms of $m(bb^{\text{add}})$ in the phase space with at least four b-jets. The overflow is included in the last bin.

The correlation matrix for the measured normalised differential cross-section in terms of $p_{\text{T}}(bb^{\text{add}})$ in the phase space with at least four b-jets. The overflow is included in the last bin.

The correlation matrix for the measured normalised differential cross-section in terms of $\text{min}\Delta R(bb)$ in the phase space with at least four b-jets. The overflow is included in the last bin.

The correlation matrix for the measured normalised differential cross-section in terms of $\Delta R(b_{1}b_{2})$ in the phase space with at least four b-jets. The overflow is included in the last bin.

The correlation matrix for the measured normalised differential cross-section in terms of $N_{l/c-\text{jets}}$ in the phase space with at least four b-jets. The overflow is included in the last bin.

The correlation matrix for the measured normalised differential cross-section in terms of $\Delta R(e\mu b_{1}b_{2},b_{3})$ in the phase space with at least four b-jets. The overflow is included in the last bin.

The correlation matrix for the measured normalised differential cross-section in terms of $\Delta R(e\mu bb^{\text{top}}, b_{1}^{\text{add}})$ in the phase space with at least four b-jets. The overflow is included in the last bin.

The correlation matrix for the measured normalised differential cross-section in terms of $\Delta R(e\mu bb^{\text{top}}, l/c\text{-jet}_{1})$ in the phase space with at least four b-jets and at least one $l/c$-jet. The overflow is included in the last bin.

The correlation matrix for the measured normalised differential cross-section in terms of $p_{\text{T}}(l/c\text{-jet}_{1}) - p_{\text{T}}(b_{1}^{\text{add}})$ in the phase space with at least four b-jets and at least one $l/c$-jet. The overflow is included in the last bin.

The correlation matrix for the measured normalised differential cross-section in terms of $|\eta(l/c\text{-jet}_{1})|$ in the phase space with at least four b-jets and at least one $l/c$-jet. The overflow is included in the last bin.

The correlation matrix for the measured normalised differential cross-section in terms of $p_{\text{T}}(l/c\text{-jet}_{1})$ in the phase space with at least four b-jets and at least one $l/c$-jet. The overflow is included in the last bin.

Expected and observed combined limits at 95% CL on the production cross-section of vector-like T in an $(X,T)$ doublet (50% $T\rightarrow Zt$, 50% $T\rightarrow Ht$).

Expected and observed combined limits at 95% CL on the production cross-section of vector-like T in a $T$ singlet (50% $T\rightarrow Wb$, 25% $T\rightarrow Zt$, 25% $T\rightarrow Ht$).

Expected and observed combined limits at 95% CL on the production cross-section of vector-like B for 100% $B\rightarrow Zb$.

Expected and observed combined limits at 95% CL on the production cross-section of vector-like B in a $(B,Y)$ doublet (50% $B\rightarrow Zb$, 50% $B\rightarrow Hb$).

Expected and observed combined limits at 95% CL on the production cross-section of vector-like B in a B singlet (50% $B\rightarrow Wt$, 25% $B\rightarrow Zb$, 25% $B\rightarrow Hb$).

Distribution of the final discriminant after the combined background-only fit in the 2$\ell$ 2$b$ double-tag 1 signal region. The signal contribution is the expected distribution for a singlet TT with a mass of 1.2 TeV. The last bin contains overflow.

Distribution of the final discriminant after the combined background-only fit in the 2$\ell$ 2$b$ top-tag signal region. The signal contribution is the expected distribution for a singlet TT with a mass of 1.2 TeV. The last bin contains overflow.

Distribution of the final discriminant after the combined background-only fit in the 2$\ell$ 1$b$ $V$-tag signal region. The signal contribution is the expected distribution for a singlet TT with a mass of 1.2 TeV. The last bin contains overflow.

Distribution of the final discriminant after the combined background-only fit in the 2$\ell$ 1$b$ no-tag signal region. The signal contribution is the expected distribution for a singlet TT with a mass of 1.2 TeV. The last bin contains overflow.

Distribution of the final discriminant after the combined background-only fit in the 3$\ell$ $VV$ control region. The signal contribution is the expected distribution for a singlet TT with a mass of 1.2 TeV. The last bin contains overflow.

Distribution of the final discriminant after the combined background-only fit in the 3$\ell$ $H$-tag signal region. The signal contribution is the expected distribution for a singlet TT with a mass of 1.2 TeV. The last bin contains overflow.

Distribution of the final discriminant after the combined background-only fit in the 3$\ell$ top-tag signal region. The signal contribution is the expected distribution for a singlet TT with a mass of 1.2 TeV. The last bin contains overflow.

Distribution of the final discriminant after the combined background-only fit in the 2$\ell$ 1$b$ control region. The signal contribution is the expected distribution for a singlet TT with a mass of 1.2 TeV. The last bin contains overflow.

Distribution of the final discriminant after the combined background-only fit in the 2$\ell$ 2$b$ control region. The signal contribution is the expected distribution for a singlet TT with a mass of 1.2 TeV. The last bin contains overflow.

Summary of observed and predicted yields in the control regions and signal region categories. The background prediction is shown after the combined likelihood fit to data under the background-only hypothesis across all control region and signal region categories. The uncertainties include statistical and systematic uncertainties. Due to correlations, the total background uncertainty is not necessarily equal to the sum in quadrature of the uncertainties of the individual background processes. The expected yields for benchmark signals with mVLQ=1.2 TeV obtained by using their theoretical cross-sections are also shown with their pre-fit uncertainties. Dashes indicate negligible yields.

Summary of observed and predicted yields in the control regions and signal region categories. The background prediction is shown before the combined likelihood fit to data under the background-only hypothesis across all control region and signal region categories. The uncertainties include statistical and systematic uncertainties. The expected yields for benchmark signals with mVLQ=1.2 TeV obtained by using their theoretical cross sections are also shown with their pre-fit uncertainties. Dashes indicate negligible yields.

Summary of data and background yields in all analysis categories after the background-only fit. The expected yields for singlet TT̄ signal with mT=1.2 TeV are also shown for comparison. The bottom panel shows the ratio of the data over the total background prediction. The shaded band includes statistical and post-fit systematic uncertainties.

Unfolded $R(\rho_{s})$ observable in the 2-to-3 measurement (normalized and divided by bin width). The parton level is defined with two stable top-quarks and a jet with $p_{T}>50$ GeV and $|\eta|<2.5$.

Covariance matrix for statistical effects of the measured $R(\rho_{s})$ observable after unfolding, for the 2-to-3 measurement (normalized)

Covariance matrix for statistical and systematic effects of the measured $R(\rho_{s})$ observable after unfolding, for the 2-to-3 measurement (normalized)

Impact of systematic uncertainties on the 2-to-3 unfolded observable. Values are given in percentage of bin content.

Central value and breakdown of the uncertainties affecting the top-quark pole mass extraction from the 2-to-3 unfolded observable.

Unfolded number of events in the 2-to-7measurement (not normalized). The parton level is defined with two neutrinos, one electron and one muon of opposite electric charges, two $b$-jets and an additional jet (extrajet). The four-momentum of the sum of neutrinos has transverse component larger than 30 GeV. The $p_{T}$-leading lepton has $p_{T}>28$ GeV, while the sub-leading has $p_{T}>20$ GeV. The two $b$-jets with have $p_{T}>30$ GeV and the extrajet has $p^\text{extrajet}_{T}>60$ GeV. All the leptons and jets are separated by $\Delta R >0.4$ and have $|\eta|<2.5$.

Covariance matrix for statistical effects of the measured number of events after unfolding, for the 2-to-7 measurement (not normalized)

Covariance matrix for statistical and systematic effects of the measured number of events after unfolding, for the 2-to-7 measurement (not normalized)

Unfolded $R(\rho_{s})$ observable in the 2-to-7 measurement (normalized and divided by bin width). The parton level is defined with two neutrinos, one electron and one muon of opposite electric charges, two $b$-jets and an additional jet (extrajet). The four-momentum of the sum of neutrinos has transverse component larger than 30 GeV. The $p_{T}$-leading lepton has $p_{T}>28$ GeV, while the sub-leading has $p_{T}>20$ GeV. The two $b$-jets with have $p_{T}>30$ GeV and the extrajet has $p^\text{extrajet}_{T}>60$ GeV. All the leptons and jets are separated by $\Delta R >0.4$ and have $|\eta|<2.5$.

Covariance matrix for statistical effects of the measured $R(\rho_{s})$ observable after unfolding, for the 2-to-7 measurement (normalized)

Covariance matrix for statistical and systematic effects of the measured $R(\rho_{s})$ observable after unfolding, for the 2-to-7 measurement (normalized)

Impact of systematic uncertainties on the 2-to-7 unfolded observable. Values are given in percentage of bin content.

Central value and breakdown of the uncertainties affecting the top-quark pole mass extraction from the 2-to-7 unfolded observable.

Three-body selection. Distributions of $M_{\Delta}^R$ in (a,b) $SR_{W}^{3-body}$ and (c,d) $SR_{T}^{3-body}$ for (left) same-flavour and (right) different-flavour events satisfying the selection criteria of the given SR, except the one for the presented variable, after the background fit. The contributions from all SM backgrounds are shown as a histogram stack. ''Others'' includes contributions from $VVV$, $t\bar{t} t$, $t\bar{t}t\bar{t}$, $t\bar{t} W$, $t\bar{t} WW$, $t\bar{t} WZ$, $t\bar{t} H$, and $tZ$ processes. The hatched bands represent the total statistical and systematic uncertainty. The rightmost bin of each plot includes overflow events. Reference top squark pair production signal models are overlayed for comparison. Red arrows in the upper panels indicate the signal region selection criteria. The bottom panels show the ratio of the observed data to the total SM background prediction, with hatched bands representing the total uncertainty in the background prediction; red arrows show data outside the vertical-axis range.

Three-body selection. Distributions of $M_{\Delta}^R$ in (a,b) $SR_{W}^{3-body}$ and (c,d) $SR_{T}^{3-body}$ for (left) same-flavour and (right) different-flavour events satisfying the selection criteria of the given SR, except the one for the presented variable, after the background fit. The contributions from all SM backgrounds are shown as a histogram stack. ''Others'' includes contributions from $VVV$, $t\bar{t} t$, $t\bar{t}t\bar{t}$, $t\bar{t} W$, $t\bar{t} WW$, $t\bar{t} WZ$, $t\bar{t} H$, and $tZ$ processes. The hatched bands represent the total statistical and systematic uncertainty. The rightmost bin of each plot includes overflow events. Reference top squark pair production signal models are overlayed for comparison. Red arrows in the upper panels indicate the signal region selection criteria. The bottom panels show the ratio of the observed data to the total SM background prediction, with hatched bands representing the total uncertainty in the background prediction; red arrows show data outside the vertical-axis range.

Three-body selection. Distributions of $M_{\Delta}^R$ in (a,b) $SR_{W}^{3-body}$ and (c,d) $SR_{T}^{3-body}$ for (left) same-flavour and (right) different-flavour events satisfying the selection criteria of the given SR, except the one for the presented variable, after the background fit. The contributions from all SM backgrounds are shown as a histogram stack. ''Others'' includes contributions from $VVV$, $t\bar{t} t$, $t\bar{t}t\bar{t}$, $t\bar{t} W$, $t\bar{t} WW$, $t\bar{t} WZ$, $t\bar{t} H$, and $tZ$ processes. The hatched bands represent the total statistical and systematic uncertainty. The rightmost bin of each plot includes overflow events. Reference top squark pair production signal models are overlayed for comparison. Red arrows in the upper panels indicate the signal region selection criteria. The bottom panels show the ratio of the observed data to the total SM background prediction, with hatched bands representing the total uncertainty in the background prediction; red arrows show data outside the vertical-axis range.

Four-body selection. (a) distributions of $E_{T}^{miss}$ in $SR^{4-body}_{Small\,\Delta m}$ and (b) distribution of $R_{2\ell 4j}$ in $SR^{4-body}_{Large\,\Delta m}$ for events satisfying the selection criteria of the given SR, except the one for the presented variable, after the background fit. The contributions from all SM backgrounds are shown as a histogram stack. ''Others'' includes contributions from $VVV$, $t\bar{t} t$, $t\bar{t}$, $t\bar{t} W$, $t\bar{t} WW$, $t\bar{t} WZ$, $t\bar{t} H$, and $tZ$ processes. The hatched bands represent the total statistical and systematic uncertainty. The rightmost bin of each plot includes overflow events. Reference top squark pair production signal models are overlayed for comparison. Red arrows in the upper panel indicate the signal region selection criteria. The bottom panels show the ratio of the observed data to the total SM background prediction, with hatched bands representing the total uncertainty in the background prediction; red arrows show data outside the vertical-axis range.

Four-body selection. (a) distributions of $E_{T}^{miss}$ in $SR^{4-body}_{Small\,\Delta m}$ and (b) distribution of $R_{2\ell 4j}$ in $SR^{4-body}_{Large\,\Delta m}$ for events satisfying the selection criteria of the given SR, except the one for the presented variable, after the background fit. The contributions from all SM backgrounds are shown as a histogram stack. ''Others'' includes contributions from $VVV$, $t\bar{t} t$, $t\bar{t}$, $t\bar{t} W$, $t\bar{t} WW$, $t\bar{t} WZ$, $t\bar{t} H$, and $tZ$ processes. The hatched bands represent the total statistical and systematic uncertainty. The rightmost bin of each plot includes overflow events. Reference top squark pair production signal models are overlayed for comparison. Red arrows in the upper panel indicate the signal region selection criteria. The bottom panels show the ratio of the observed data to the total SM background prediction, with hatched bands representing the total uncertainty in the background prediction; red arrows show data outside the vertical-axis range.

Exclusion limit contour (95% CL) for a simplified model assuming $\tilde{t}_1$ pair production, decaying via $\tilde{t}_1 \rightarrow t^{(*)}\tilde{\chi}_1^0$ with 100% branching ratio, in the (a) $m(\tilde{t}_1)$--$m(\tilde{\chi}_1^0)$ and (b) $m(\tilde{t}_1)$--$\Delta m(\tilde{t}_1,\tilde{\chi}_1^0)$ planes. The dashed lines and the shaded bands are the expected limits and their $\pm1\sigma$ uncertainties. The thick solid lines are the observed limits for the central value of the signal cross-section. The expected and observed limits do not include the effect of the theoretical uncertainties in the signal cross-section. The dotted lines show the effect on the observed limit when varying the signal cross-section by $\pm1\sigma$ of the theoretical uncertainty.

Exclusion limit contour (95% CL) for a simplified model assuming $\tilde{t}_1$ pair production, decaying via $\tilde{t}_1 \rightarrow t^{(*)}\tilde{\chi}_1^0$ with 100% branching ratio, in the (a) $m(\tilde{t}_1)$--$m(\tilde{\chi}_1^0)$ and (b) $m(\tilde{t}_1)$--$\Delta m(\tilde{t}_1,\tilde{\chi}_1^0)$ planes. The dashed lines and the shaded bands are the expected limits and their $\pm1\sigma$ uncertainties. The thick solid lines are the observed limits for the central value of the signal cross-section. The expected and observed limits do not include the effect of the theoretical uncertainties in the signal cross-section. The dotted lines show the effect on the observed limit when varying the signal cross-section by $\pm1\sigma$ of the theoretical uncertainty.

Exclusion limit contour (95% CL) for a simplified model assuming $\tilde{t}_1$ pair production, decaying via $\tilde{t}_1 \rightarrow t^{(*)}\tilde{\chi}_1^0$ with 100% branching ratio, in the (a) $m(\tilde{t}_1)$--$m(\tilde{\chi}_1^0)$ and (b) $m(\tilde{t}_1)$--$\Delta m(\tilde{t}_1,\tilde{\chi}_1^0)$ planes. The dashed lines and the shaded bands are the expected limits and their $\pm1\sigma$ uncertainties. The thick solid lines are the observed limits for the central value of the signal cross-section. The expected and observed limits do not include the effect of the theoretical uncertainties in the signal cross-section. The dotted lines show the effect on the observed limit when varying the signal cross-section by $\pm1\sigma$ of the theoretical uncertainty.

Exclusion limit contour (95% CL) for a simplified model assuming $\tilde{t}_1$ pair production, decaying via $\tilde{t}_1 \rightarrow t^{(*)}\tilde{\chi}_1^0$ with 100% branching ratio, in the (a) $m(\tilde{t}_1)$--$m(\tilde{\chi}_1^0)$ and (b) $m(\tilde{t}_1)$--$\Delta m(\tilde{t}_1,\tilde{\chi}_1^0)$ planes. The dashed lines and the shaded bands are the Observed limits and their $\pm1\sigma$ uncertainties. The thick solid lines are the observed limits for the central value of the signal cross-section. The expected and observed limits do not include the effect of the theoretical uncertainties in the signal cross-section. The dotted lines show the effect on the observed limit when varying the signal cross-section by $\pm1\sigma$ of the theoretical uncertainty.

Exclusion limit contour (95% CL) for a simplified model assuming $\tilde{t}_1$ pair production, decaying via $\tilde{t}_1 \rightarrow t^{(*)}\tilde{\chi}_1^0$ with 100% branching ratio, in the (a) $m(\tilde{t}_1)$--$m(\tilde{\chi}_1^0)$ and (b) $m(\tilde{t}_1)$--$\Delta m(\tilde{t}_1,\tilde{\chi}_1^0)$ planes. The dashed lines and the shaded bands are the expected limits and their $\pm1\sigma$ uncertainties. The thick solid lines are the observed limits for the central value of the signal cross-section. The expected and observed limits do not include the effect of the theoretical uncertainties in the signal cross-section. The dotted lines show the effect on the observed limit when varying the signal cross-section by $\pm1\sigma$ of the theoretical uncertainty.

Exclusion limit contour (95% CL) for a simplified model assuming $\tilde{t}_1$ pair production, decaying via $\tilde{t}_1 \rightarrow t^{(*)}\tilde{\chi}_1^0$ with 100% branching ratio, in the (a) $m(\tilde{t}_1)$--$m(\tilde{\chi}_1^0)$ and (b) $m(\tilde{t}_1)$--$\Delta m(\tilde{t}_1,\tilde{\chi}_1^0)$ planes. The dashed lines and the shaded bands are the expected limits and their $\pm1\sigma$ uncertainties. The thick solid lines are the observed limits for the central value of the signal cross-section. The expected and observed limits do not include the effect of the theoretical uncertainties in the signal cross-section. The dotted lines show the effect on the observed limit when varying the signal cross-section by $\pm1\sigma$ of the theoretical uncertainty.

Exclusion limit contour (95% CL) for a simplified model assuming $\tilde{t}_1$ pair production, decaying via $\tilde{t}_1 \rightarrow t^{(*)}\tilde{\chi}_1^0$ with 100\% branching ratio, in the (a) $m(\tilde{t}_1)$--$m(\tilde{\chi}_1^0)$ and (b) $m(\tilde{t}_1)$--$\Delta m(\tilde{t}_1,\tilde{\chi}_1^0)$ planes. The dashed lines and the shaded bands are the expected limits and their $\pm1\sigma$ uncertainties. The thick solid lines are the observed limits for the central value of the signal cross-section. The expected and observed limits do not include the effect of the theoretical uncertainties in the signal cross-section. The dotted lines show the effect on the observed limit when varying the signal cross-section by $\pm1\sigma$ of the theoretical uncertainty.

Exclusion limit contour (95% CL) for a simplified model assuming $\tilde{t}_1$ pair production, decaying via $\tilde{t}_1 \rightarrow t^{(*)}\tilde{\chi}_1^0$ with 100\% branching ratio, in the (a) $m(\tilde{t}_1)$--$m(\tilde{\chi}_1^0)$ and (b) $m(\tilde{t}_1)$--$\Delta m(\tilde{t}_1,\tilde{\chi}_1^0)$ planes. The dashed lines and the shaded bands are the expected limits and their $\pm1\sigma$ uncertainties. The thick solid lines are the observed limits for the central value of the signal cross-section. The expected and observed limits do not include the effect of the theoretical uncertainties in the signal cross-section. The dotted lines show the effect on the observed limit when varying the signal cross-section by $\pm1\sigma$ of the theoretical uncertainty.

Exclusion limit contour (95% CL) for a simplified model assuming $\tilde{t}_1$ pair production, decaying via $\tilde{t}_1 \rightarrow t^{(*)}\tilde{\chi}_1^0$ with 100\% branching ratio, in the (a) $m(\tilde{t}_1)$--$m(\tilde{\chi}_1^0)$ and (b) $m(\tilde{t}_1)$--$\Delta m(\tilde{t}_1,\tilde{\chi}_1^0)$ planes. The dashed lines and the shaded bands are the expected limits and their $\pm1\sigma$ uncertainties. The thick solid lines are the observed limits for the central value of the signal cross-section. The expected and observed limits do not include the effect of the theoretical uncertainties in the signal cross-section. The dotted lines show the effect on the observed limit when varying the signal cross-section by $\pm1\sigma$ of the theoretical uncertainty.

Exclusion limit contour (95% CL) for a simplified model assuming $\tilde{t}_1$ pair production, decaying via $\tilde{t}_1 \rightarrow t^{(*)}\tilde{\chi}_1^0$ with 100\% branching ratio, in the (a) $m(\tilde{t}_1)$--$m(\tilde{\chi}_1^0)$ and (b) $m(\tilde{t}_1)$--$\Delta m(\tilde{t}_1,\tilde{\chi}_1^0)$ planes. The dashed lines and the shaded bands are the expected limits and their $\pm1\sigma$ uncertainties. The thick solid lines are the observed limits for the central value of the signal cross-section. The expected and observed limits do not include the effect of the theoretical uncertainties in the signal cross-section. The dotted lines show the effect on the observed limit when varying the signal cross-section by $\pm1\sigma$ of the theoretical uncertainty.

Exclusion limit contour (95% CL) for a simplified model assuming $\tilde{t}_1$ pair production, decaying via $\tilde{t}_1 \rightarrow t^{(*)}\tilde{\chi}_1^0$ with 100\% branching ratio, in the (a) $m(\tilde{t}_1)$--$m(\tilde{\chi}_1^0)$ and (b) $m(\tilde{t}_1)$--$\Delta m(\tilde{t}_1,\tilde{\chi}_1^0)$ planes. The dashed lines and the shaded bands are the expected limits and their $\pm1\sigma$ uncertainties. The thick solid lines are the observed limits for the central value of the signal cross-section. The expected and observed limits do not include the effect of the theoretical uncertainties in the signal cross-section. The dotted lines show the effect on the observed limit when varying the signal cross-section by $\pm1\sigma$ of the theoretical uncertainty.

Exclusion limit contour (95% CL) for a simplified model assuming $\tilde{t}_1$ pair production, decaying via $\tilde{t}_1 \rightarrow t^{(*)}\tilde{\chi}_1^0$ with 100\% branching ratio, in the (a) $m(\tilde{t}_1)$--$m(\tilde{\chi}_1^0)$ and (b) $m(\tilde{t}_1)$--$\Delta m(\tilde{t}_1,\tilde{\chi}_1^0)$ planes. The dashed lines and the shaded bands are the expected limits and their $\pm1\sigma$ uncertainties. The thick solid lines are the observed limits for the central value of the signal cross-section. The expected and observed limits do not include the effect of the theoretical uncertainties in the signal cross-section. The dotted lines show the effect on the observed limit when varying the signal cross-section by $\pm1\sigma$ of the theoretical uncertainty.

Exclusion limits for (a) $t\bar{t} + \phi $ scalar and (b) $t\bar{t} + a $ pseudoscalar models as a function of the mediator mass for a DM particle mass of $m(\chi)=1$ GeV. 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 shows the observed (expected) exclusion limits.

Exclusion limits for (a) $t\bar{t} + \phi $ scalar and (b) $t\bar{t} + a $ pseudoscalar models as a function of the mediator mass for a DM particle mass of $m(\chi)=1$ GeV. 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 shows the observed (expected) exclusion limits.

Exclusion limits for (a) $t\bar{t} + \phi $ scalar and (b) $t\bar{t} + a $ pseudoscalar models as a function of the mediator mass for a DM particle mass of $m(\chi)=1$ GeV. 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 shows the observed (expected) exclusion limits.

Exclusion limits for (a) $t\bar{t} + \phi $ scalar and (b) $t\bar{t} + a $ pseudoscalar models as a function of the mediator mass for a DM particle mass of $m(\chi)=1$ GeV. 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 shows the observed (expected) exclusion limits.

Two-body selection. Background fit results for $\mathrm{CR}^{\mathrm{2-body}}_{t\bar{t}}$, $\mathrm{CR}^{\mathrm{2-body}}_{t\bar{t}Z}$, $\mathrm{VR}^{\mathrm{2-body}}_{t\bar{t}, DF}$, $\mathrm{VR}^{\mathrm{2-body}}_{t\bar{t}, SF}$ and $\mathrm{VR}^{\mathrm{2-body}}_{t\bar{t} Z}$. ''Others'' includes contributions from $VVV$, $t\bar{t} t$, $t\bar{t}t\bar{t}$, $t\bar{t} W$, $t\bar{t} WW$, $t\bar{t} WZ$, $t\bar{t} H$, and $tZ$ processes. Combined statistical and systematic uncertainties are given. Entries marked `--' indicate a negligible background contribution (less than 0.001 events). The individual uncertainties can be correlated, and do not necessarily add up in quadrature to the total background uncertainty.

Three-body selection. Background fit results for $\mathrm{CR}^{\mathrm{3-body}}_{t\bar{t}}$, $\mathrm{CR}^{\mathrm{3-body}}_{VV}$, $\mathrm{CR}^{\mathrm{2-body}}_{t\bar{t}Z}$, $\mathrm{VR}^{\mathrm{3-body}}_{VV}$, $\mathrm{VR(1)}^{\mathrm{3-body}}_{t\bar{t}}$ and $\mathrm{VR(2)}^{\mathrm{3-body}}_{t\bar{t}}$. ''Others'' includes contributions from $VVV$, $t\bar{t} t$, $t\bar{t}t\bar{t}$, $t\bar{t} W$, $t\bar{t} WW$, $t\bar{t} WZ$, $t\bar{t} H$, and $tZ$ processes. Combined statistical and systematic uncertainties are given. Entries marked `--' indicate a negligible background contribution (less than 0.001 events). The individual uncertainties can be correlated, and do not necessarily add up in quadrature to the total background uncertainty.

Four-body selection. Background fit results for $\mathrm{CR}^{\mathrm{4-body}}_{t\bar{t}}$,$\mathrm{CR}^{\mathrm{4-body}}_{VV}$, $\mathrm{VR}^{\mathrm{4-body}}_{t\bar{t}}$, $VR^{4-body}_{VV}$ and $\mathrm{VR}^{\mathrm{4-body}}_{VV,lll}$. The ''Others'' category contains the contributions from $VVV$, $t\bar{t} t$, $t\bar{t}t\bar{t}$, $t\bar{t} W$, $t\bar{t} WW$, $t\bar{t} WZ$, $t\bar{t} H$, and $tZ$. Combined statistical and systematic uncertainties are given. Entries marked `--' indicate a negligible background contribution (less than 0.001 events). The individual uncertainties can be correlated, and do not necessarily add up in quadrature to the total background uncertainty.

Two-body selection. Background fit results for the different-flavour leptons binned SRs. The ''Others'' category contains the contributions from $VVV$, $t\bar{t} t$, $t\bar{t}t\bar{t}$, $t\bar{t} W$, $t\bar{t} WW$, $t\bar{t} WZ$, $t\bar{t} H$, and $tZ$. Combined statistical and systematic uncertainties are given. Entries marked `--' indicate a negligible background contribution (less than 0.001 events). The individual uncertainties can be correlated, and do not necessarily add up in quadrature to the total background uncertainty.

Two-body selection. Background fit results for the same-flavour leptons binned SRs. The ''Others'' category contains the contributions from $VVV$, $t\bar{t} t$, $t\bar{t}t\bar{t}$, $t\bar{t} W$, $t\bar{t} WW$, $t\bar{t} WZ$, $t\bar{t} H$, and $tZ$. Combined statistical and systematic uncertainties are given. The individual uncertainties can be correlated, and do not necessarily add up in quadrature to the total background uncertainty.

Three-body selection. Observed event yields and background fit results for the three-body selection SRs. The ''Others'' category contains contributions from $VVV$, $t\bar{t} t$, $t\bar{t}t\bar{t}$, $t\bar{t} W$, $t\bar{t} WW$, $t\bar{t} WZ$, $t\bar{t} H$, and $tZ$. Combined statistical and systematic uncertainties are given. Entries marked `--' indicate a negligible background contribution (less than 0.001 events). The individual uncertainties can be correlated, and do not necessarily add up in quadrature to the total background uncertainty.

Four-body selection. Observed event yields and background fit results for SR$^{\mathrm{4-body}}_{\mathrm{Small}\,\Delta m}$ and SR$^{\mathrm{4-body}}_{\mathrm{Large}\,\Delta m}$. The ''Others'' category contains the contributions from $VVV$, $t\bar{t} t$, $t\bar{t}t\bar{t}$, $t\bar{t} W$, $t\bar{t} WW$, $t\bar{t} WZ$, $t\bar{t} H$, and $tZ$. Combined statistical and systematic uncertainties are given. The individual uncertainties can be correlated, and do not necessarily add up in quadrature to the total background uncertainty.

Exclusion limits contours (95% CL) for a simplified model assuming $\tilde{t}_1$ pair production, decaying via $\tilde{t}_1 \rightarrow t^{(*)}\tilde{\chi}^0_1$ with 100% branching ratio in $\tilde{t}_1--\tilde{\chi}^0_1$ masses planes. The dashed lines and the shaded bands are the expected limit and its $\pm 1\sigma$ uncertainty. The thick solid lines are the observed limits for the central value of the signal cross-section. The expected and observed limits do not include the effect of the theoretical uncertainties in the signal cross-section. The dotted lines show the effect on the observed limit when varying the signal cross-section by $\pm1\sigma$ of the theoretical uncertainty. The exclusion limits contours for the two-body, three-body and four-body selections are respectively shown in blue, green and red.

Exclusion limits contours (95% CL) for a simplified model assuming $\tilde{t}_1$ pair production, decaying via $\tilde{t}_1 \rightarrow t^{(*)}\tilde{\chi}^0_1$ with 100% branching ratio in $\tilde{t}_1--\tilde{\chi}^0_1$ masses planes. The dashed lines and the shaded bands are the expected limit and its $\pm 1\sigma$ uncertainty. The thick solid lines are the observed limits for the central value of the signal cross-section. The expected and observed limits do not include the effect of the theoretical uncertainties in the signal cross-section. The dotted lines show the effect on the observed limit when varying the signal cross-section by $\pm1\sigma$ of the theoretical uncertainty. The exclusion limits contours for the two-body, three-body and four-body selections are respectively shown in blue, green and red.

Exclusion limit contour (95% CL) for a simplified model assuming $\tilde{t}_1$ pair production, decaying via $\tilde{t}_1 \rightarrow t \tilde{\chi}_1^0$ with 100% branching ratio, in $\tilde{t}_1$--$\tilde{\chi}_1^0$ masses plane. The dashed lines and the shaded bands are the expected limit and its $\pm1\sigma$ uncertainty. The thick solid lines are the observed limits for the central value of the signal cross-section. The expected and observed limits do not include the effect of the theoretical uncertainties in the signal cross-section. The dotted lines show the effect on the observed limit when varying the signal cross-section by $\pm1\sigma$ of the theoretical uncertainty. The observed (a) and expected (b) CLs values are respectively shown.

Exclusion limit contour (95% CL) for a simplified model assuming $\tilde{t}_1$ pair production, decaying via $\tilde{t}_1 \rightarrow t \tilde{\chi}_1^0$ with 100% branching ratio, in $\tilde{t}_1$--$\tilde{\chi}_1^0$ masses plane. The dashed lines and the shaded bands are the expected limit and its $\pm1\sigma$ uncertainty. The thick solid lines are the observed limits for the central value of the signal cross-section. The expected and observed limits do not include the effect of the theoretical uncertainties in the signal cross-section. The dotted lines show the effect on the observed limit when varying the signal cross-section by $\pm1\sigma$ of the theoretical uncertainty. The observed (a) and expected (b) CLs values are respectively shown.

Exclusion limit contour (95% CL) for a simplified model assuming $\tilde{t}_1$ pair production, decaying via $\tilde{t}_1 \rightarrow t \tilde{\chi}_1^0$ with 100% branching ratio, in $\tilde{t}_1$--$\tilde{\chi}_1^0$ masses plane. The dashed lines and the shaded bands are the expected limit and its $\pm1\sigma$ uncertainty. The thick solid lines are the observed limits for the central value of the signal cross-section. The expected and observed limits do not include the effect of the theoretical uncertainties in the signal cross-section. The dotted lines show the effect on the observed limit when varying the signal cross-section by $\pm1\sigma$ of the theoretical uncertainty. The observed (a) and expected (b) CLs values are respectively shown.

Exclusion limit contour (95% CL) for a simplified model assuming $\tilde{t}_1$ pair production, decaying via $\tilde{t}_1 \rightarrow t \tilde{\chi}_1^0$ with 100% branching ratio, in $\tilde{t}_1$--$\tilde{\chi}_1^0$ masses plane. The dashed lines and the shaded bands are the expected limit and its $\pm1\sigma$ uncertainty. The thick solid lines are the observed limits for the central value of the signal cross-section. The expected and observed limits do not include the effect of the theoretical uncertainties in the signal cross-section. The dotted lines show the effect on the observed limit when varying the signal cross-section by $\pm1\sigma$ of the theoretical uncertainty. The observed (a) and expected (b) CLs values are respectively shown.

Exclusion limit contour (95% CL) for a simplified model assuming $\tilde{t}_1$ pair production, decaying via $\tilde{t}_1 \rightarrow b W \tilde{\chi}_1^0$ with 100% branching ratio, in $\tilde{t}_1$--$\tilde{\chi}_1^0$ masses plane. The dashed lines and the shaded bands are the expected limit and its $\pm1\sigma$ uncertainty. The thick solid lines are the observed limits for the central value of the signal cross-section. The expected and observed limits do not include the effect of the theoretical uncertainties in the signal cross-section. The dotted lines show the effect on the observed limit when varying the signal cross-section by $\pm1\sigma$ of the theoretical uncertainty. The observed (a) and expected (b) CLs values are respectively shown.

Exclusion limit contour (95% CL) for a simplified model assuming $\tilde{t}_1$ pair production, decaying via $\tilde{t}_1 \rightarrow b W \tilde{\chi}_1^0$ with 100% branching ratio, in $\tilde{t}_1$--$\tilde{\chi}_1^0$ masses plane. The dashed lines and the shaded bands are the expected limit and its $\pm1\sigma$ uncertainty. The thick solid lines are the observed limits for the central value of the signal cross-section. The expected and observed limits do not include the effect of the theoretical uncertainties in the signal cross-section. The dotted lines show the effect on the observed limit when varying the signal cross-section by $\pm1\sigma$ of the theoretical uncertainty. The observed (a) and expected (b) CLs values are respectively shown.

Exclusion limit contour (95% CL) for a simplified model assuming $\tilde{t}_1$ pair production, decaying via $\tilde{t}_1 \rightarrow b W \tilde{\chi}_1^0$ with 100% branching ratio, in $\tilde{t}_1$--$\tilde{\chi}_1^0$ masses plane. The dashed lines and the shaded bands are the expected limit and its $\pm1\sigma$ uncertainty. The thick solid lines are the observed limits for the central value of the signal cross-section. The expected and observed limits do not include the effect of the theoretical uncertainties in the signal cross-section. The dotted lines show the effect on the observed limit when varying the signal cross-section by $\pm 1\sigma$ of the theoretical uncertainty. The observed (a) and expected (b) CLs values are respectively shown.

Exclusion limit contour (95% CL) for a simplified model assuming $\tilde{t}_1$ pair production, decaying via $\tilde{t}_1 \rightarrow b W \tilde{\chi}_1^0$ with 100% branching ratio, in $\tilde{t}_1$--$\tilde{\chi}_1^0$ masses plane. The dashed lines and the shaded bands are the expected limit and its $\pm1\sigma$ uncertainty. The thick solid lines are the observed limits for the central value of the signal cross-section. The expected and observed limits do not include the effect of the theoretical uncertainties in the signal cross-section. The dotted lines show the effect on the observed limit when varying the signal cross-section by $\pm1\sigma$ of the theoretical uncertainty. The observed (a) and expected (b) CLs values are respectively shown.

Exclusion limit contour (95% CL) for a simplified model assuming $\tilde{t}_1$ pair production, decaying via $\tilde{t}_1 \rightarrow b l \nu \tilde{\chi}_1^0$ with 100% branching ratio, in $\tilde{t}_1$--$\tilde{\chi}_1^0$ masses plane. The dashed lines and the shaded bands are the expected limit and its $\pm1\sigma$ uncertainty. The thick solid lines are the observed limits for the central value of the signal cross-section. The expected and observed limits do not include the effect of the theoretical uncertainties in the signal cross-section. The dotted lines show the effect on the observed limit when varying the signal cross-section by $\pm1\sigma$ of the theoretical uncertainty. The observed (a) and expected (b) CLs values are respectively shown.

Exclusion limit contour (95% CL) for a simplified model assuming $\tilde{t}_1$ pair production, decaying via $\tilde{t}_1 \rightarrow b l \nu \tilde{\chi}_1^0$ with 100% branching ratio, in $\tilde{t}_1$--$\tilde{\chi}_1^0$ masses plane. The dashed lines and the shaded bands are the expected limit and its $\pm1\sigma$ uncertainty. The thick solid lines are the observed limits for the central value of the signal cross-section. The expected and observed limits do not include the effect of the theoretical uncertainties in the signal cross-section. The dotted lines show the effect on the observed limit when varying the signal cross-section by $\pm1\sigma$ of the theoretical uncertainty. The observed (a) and expected (b) CLs values are respectively shown.

Exclusion limit contour (95% CL) for a simplified model assuming $\tilde{t}_1$ pair production, decaying via $\tilde{t}_1 \rightarrow b l \nu \tilde{\chi}_1^0$ with 100% branching ratio, in $\tilde{t}_1$--$\tilde{\chi}_1^0$ masses plane. The dashed lines and the shaded bands are the expected limit and its $\pm1\sigma$ uncertainty.The thick solid lines are the observed limits for the central value of the signal cross-section. The expected and observed limits do not include the effect of the theoretical uncertainties in the signal cross-section. The dotted lines show the effect on the observed limit when varying the signal cross-section by $\pm1\sigma$ of the theoretical uncertainty. The observed (a) and expected (b) CLs values are respectively shown.

Exclusion limit contour (95% CL) for a simplified model assuming $\tilde{t}_1$ pair production, decaying via $\tilde{t}_1 \rightarrow b l \nu \tilde{\chi}_1^0$ with 100% branching ratio, in $\tilde{t}_1$--$\tilde{\chi}_1^0$ masses plane. The dashed lines and the shaded bands are the expected limit and its $\pm1\sigma$ uncertainty.The thick solid lines are the observed limits for the central value of the signal cross-section. The expected and observed limits do not include the effect of the theoretical uncertainties in the signal cross-section. The dotted lines show the effect on the observed limit when varying the signal cross-section by $\pm1\sigma$ of the theoretical uncertainty. The observed (a) and expected (b) CLs values are respectively shown.

Exclusion limits for (a) $t\bar{t} + \phi $ scalar and (b) $t\bar{t} + a $ pseudoscalar models as a function of the DM particle mass for a mediator mass of 10 GeV. 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 shows the observed (expected) exclusion limits.

Exclusion limits for (a) $t\bar{t} + \phi $ scalar and (b) $t\bar{t} + a $ pseudoscalar models as a function of the DM particle mass for a mediator mass of 10 GeV. 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 shows the observed (expected) exclusion limits.

Exclusion limits for (a) $t\bar{t} + \phi $ scalar and (b) $t\bar{t} + a $ pseudoscalar models as a function of the DM particle mass for a mediator mass of 10 GeV. 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 shows the observed (expected) exclusion limits.

Exclusion limits for (a) $t\bar{t} + \phi $ scalar and (b) $t\bar{t} + a $ pseudoscalar models as a function of the DM particle mass for a mediator mass of 10 GeV. 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 shows the observed (expected) exclusion limits.

Two-body selection efficiency (a) SR-DF$^{2-body}_{[110,120)}$, (b) SR-DF1$^{2-body}_{[120,140)}$, (c) SR-DF2$^{2-body}_{[140,160)}$, (d) SR-DF3$^{2-body}_{[160,180)}$, (e) SR-DF4$^{2-body}_{[180,220)}$, (f) SR-DF5$^{2-body}_{[220,\infty)}$ for a simplified model assuming $ \tilde{t}_1$ pair production.

Two-body selection efficiency (a) SR-DF$^{2-body}_{[110,120)}$, (b) SR-DF1$^{2-body}_{[120,140)}$, (c) SR-DF2$^{2-body}_{[140,160)}$, (d) SR-DF3$^{2-body}_{[160,180)}$, (e) SR-DF4$^{2-body}_{[180,220)}$, (f) SR-DF5$^{2-body}_{[220,\infty)}$ for a simplified model assuming $ \tilde{t}_1$ pair production.

Two-body selection efficiency (a) SR-DF$^{2-body}_{[110,120)}$, (b) SR-DF1$^{2-body}_{[120,140)}$, (c) SR-DF2$^{2-body}_{[140,160)}$, (d) SR-DF3$^{2-body}_{[160,180)}$, (e) SR-DF4$^{2-body}_{[180,220)}$, (f) SR-DF5$^{2-body}_{[220,\infty)}$ for a simplified model assuming $ \tilde{t}_1$ pair production.

Two-body selection efficiency (a) SR-DF$^{2-body}_{[110,120)}$, (b) SR-DF1$^{2-body}_{[120,140)}$, (c) SR-DF2$^{2-body}_{[140,160)}$, (d) SR-DF3$^{2-body}_{[160,180)}$, (e) SR-DF4$^{2-body}_{[180,220)}$, (f) SR-DF5$^{2-body}_{[220,\infty)}$ for a simplified model assuming $ \tilde{t}_1$ pair production.

Two-body selection efficiency (a) SR-DF$^{2-body}_{[110,120)}$, (b) SR-DF1$^{2-body}_{[120,140)}$, (c) SR-DF2$^{2-body}_{[140,160)}$, (d) SR-DF3$^{2-body}_{[160,180)}$, (e) SR-DF4$^{2-body}_{[180,220)}$, (f) SR-DF5$^{2-body}_{[220,\infty)}$ for a simplified model assuming $ \tilde{t}_1$ pair production.

Two-body selection efficiency (a) SR-DF$^{2-body}_{[110,120)}$, (b) SR-DF1$^{2-body}_{[120,140)}$, (c) SR-DF2$^{2-body}_{[140,160)}$, (d) SR-DF3$^{2-body}_{[160,180)}$, (e) SR-DF4$^{2-body}_{[180,220)}$, (f) SR-DF5$^{2-body}_{[220,\infty)}$ for a simplified model assuming $ \tilde{t}_1$ pair production.

Two-body selection efficiency (a) SR-SF$^{2-body}_{[110,120)}$, (b) SR-SF1$^{2-body}_{[120,140)}$, (c) SR-SF2$^{2-body}_{[140,160)}$, (d) SR-SF3$^{2-body}_{[160,180)}$, (e) SR-SF4$^{2-body}_{[180,220)}$, (f) SR-SF5$^{2-body}_{[220,\infty)}$ for a simplified model assuming $ \tilde{t}_1$ pair production.

Two-body selection efficiency (a) SR-SF$^{2-body}_{[110,120)}$, (b) SR-SF1$^{2-body}_{[120,140)}$, (c) SR-SF2$^{2-body}_{[140,160)}$, (d) SR-SF3$^{2-body}_{[160,180)}$, (e) SR-SF4$^{2-body}_{[180,220)}$, (f) SR-SF5$^{2-body}_{[220,\infty)}$ for a simplified model assuming $ \tilde{t}_1$ pair production.

Two-body selection efficiency (a) SR-SF$^{2-body}_{[110,120)}$, (b) SR-SF1$^{2-body}_{[120,140)}$, (c) SR-SF2$^{2-body}_{[140,160)}$, (d) SR-SF3$^{2-body}_{[160,180)}$, (e) SR-SF4$^{2-body}_{[180,220)}$, (f) SR-SF5$^{2-body}_{[220,\infty)}$ for a simplified model assuming $ \tilde{t}_1$ pair production.

Two-body selection efficiency (a) SR-SF$^{2-body}_{[110,120)}$, (b) SR-SF1$^{2-body}_{[120,140)}$, (c) SR-SF2$^{2-body}_{[140,160)}$, (d) SR-SF3$^{2-body}_{[160,180)}$, (e) SR-SF4$^{2-body}_{[180,220)}$, (f) SR-SF5$^{2-body}_{[220,\infty)}$ for a simplified model assuming $ \tilde{t}_1$ pair production.

Two-body selection efficiency (a) SR-SF$^{2-body}_{[110,120)}$, (b) SR-SF1$^{2-body}_{[120,140)}$, (c) SR-SF2$^{2-body}_{[140,160)}$, (d) SR-SF3$^{2-body}_{[160,180)}$, (e) SR-SF4$^{2-body}_{[180,220)}$, (f) SR-SF5$^{2-body}_{[220,\infty)}$ for a simplified model assuming $ \tilde{t}_1$ pair production.

Two-body selection efficiency (a) SR-SF$^{2-body}_{[110,120)}$, (b) SR-SF1$^{2-body}_{[120,140)}$, (c) SR-SF2$^{2-body}_{[140,160)}$, (d) SR-SF3$^{2-body}_{[160,180)}$, (e) SR-SF4$^{2-body}_{[180,220)}$, (f) SR-SF5$^{2-body}_{[220,\infty)}$ for a simplified model assuming $ \tilde{t}_1$ pair production.

Two-body selection efficiency (a) $SR^{2-body}_{[110,\infty)}$ , (b) $SR^{2-body}_{[120,\infty)}$ , (c) $SR^{2-body}_{[140,\infty)}$ , (d) $SR^{2-body}_{[160,\infty)}$ , (e) $SR^{2-body}_{[180,\infty)}$ , (f) $SR^{2-body}_{[200,\infty)}$ , (g) $SR^{2-body}_{[220,\infty)}$ for a simplified model assuming $ \tilde{t}_1$ pair production.

Two-body selection efficiency (a) $SR^{2-body}_{[110,\infty)}$ , (b) $SR^{2-body}_{[120,\infty)}$ , (c) $SR^{2-body}_{[140,\infty)}$ , (d) $SR^{2-body}_{[160,\infty)}$ , (e) $SR^{2-body}_{[180,\infty)}$ , (f) $SR^{2-body}_{[200,\infty)}$ , (g) $SR^{2-body}_{[220,\infty)}$ for a simplified model assuming $ \tilde{t}_1$ pair production.

Two-body selection efficiency (a) $SR^{2-body}_{[110,\infty)}$ , (b) $SR^{2-body}_{[120,\infty)}$ , (c) $SR^{2-body}_{[140,\infty)}$ , (d) $SR^{2-body}_{[160,\infty)}$ , (e) $SR^{2-body}_{[180,\infty)}$ , (f) $SR^{2-body}_{[200,\infty)}$ , (g) $SR^{2-body}_{[220,\infty)}$ for a simplified model assuming $ \tilde{t}_1$ pair production.

Two-body selection efficiency (a) $SR^{2-body}_{[110,\infty)}$ , (b) $SR^{2-body}_{[120,\infty)}$ , (c) $SR^{2-body}_{[140,\infty)}$ , (d) $SR^{2-body}_{[160,\infty)}$ , (e) $SR^{2-body}_{[180,\infty)}$ , (f) $SR^{2-body}_{[200,\infty)}$ , (g) $SR^{2-body}_{[220,\infty)}$ for a simplified model assuming $ \tilde{t}_1$ pair production.

Two-body selection efficiency (a) $SR^{2-body}_{[110,\infty)}$ , (b) $SR^{2-body}_{[120,\infty)}$ , (c) $SR^{2-body}_{[140,\infty)}$ , (d) $SR^{2-body}_{[160,\infty)}$ , (e) $SR^{2-body}_{[180,\infty)}$ , (f) $SR^{2-body}_{[200,\infty)}$ , (g) $SR^{2-body}_{[220,\infty)}$ for a simplified model assuming $ \tilde{t}_1$ pair production.

Two-body selection efficiency (a) $SR^{2-body}_{[110,\infty)}$ , (b) $SR^{2-body}_{[120,\infty)}$ , (c) $SR^{2-body}_{[140,\infty)}$ , (d) $SR^{2-body}_{[160,\infty)}$ , (e) $SR^{2-body}_{[180,\infty)}$ , (f) $SR^{2-body}_{[200,\infty)}$ , (g) $SR^{2-body}_{[220,\infty)}$ for a simplified model assuming $ \tilde{t}_1$ pair production.

Two-body selection efficiency (a) $SR^{2-body}_{[110,\infty)}$ , (b) $SR^{2-body}_{[120,\infty)}$ , (c) $SR^{2-body}_{[140,\infty)}$ , (d) $SR^{2-body}_{[160,\infty)}$ , (e) $SR^{2-body}_{[180,\infty)}$ , (f) $SR^{2-body}_{[200,\infty)}$ , (g) $SR^{2-body}_{[220,\infty)}$ for a simplified model assuming $ \tilde{t}_1$ pair production.

Two-body selection efficiency (a) SR-DF$^{2-body}_{[110,120)}$, (b) SR-DF1$^{2-body}_{[120,140)}$, (c) SR-DF2$^{2-body}_{[140,160)}$, (d) SR-DF3$^{2-body}_{[160,180)}$, (e) SR-DF4$^{2-body}_{[180,220)}$, (f) SR-DF5$^{2-body}_{[220,\infty)}$ for a simplified model assuming $t \tilde{t} + \phi$ pair production.

Two-body selection efficiency (a) SR-DF$^{2-body}_{[110,120)}$, (b) SR-DF1$^{2-body}_{[120,140)}$, (c) SR-DF2$^{2-body}_{[140,160)}$, (d) SR-DF3$^{2-body}_{[160,180)}$, (e) SR-DF4$^{2-body}_{[180,220)}$, (f) SR-DF5$^{2-body}_{[220,\infty)}$ for a simplified model assuming $t \tilde{t} + \phi$ pair production.

Two-body selection efficiency (a) SR-DF$^{2-body}_{[110,120)}$, (b) SR-DF1$^{2-body}_{[120,140)}$, (c) SR-DF2$^{2-body}_{[140,160)}$, (d) SR-DF3$^{2-body}_{[160,180)}$, (e) SR-DF4$^{2-body}_{[180,220)}$, (f) SR-DF5$^{2-body}_{[220,\infty)}$ for a simplified model assuming $t \tilde{t} + \phi$ pair production.

Two-body selection efficiency (a) SR-DF$^{2-body}_{[110,120)}$, (b) SR-DF1$^{2-body}_{[120,140)}$, (c) SR-DF2$^{2-body}_{[140,160)}$, (d) SR-DF3$^{2-body}_{[160,180)}$, (e) SR-DF4$^{2-body}_{[180,220)}$, (f) SR-DF5$^{2-body}_{[220,\infty)}$ for a simplified model assuming $t \tilde{t} + \phi$ pair production.

Two-body selection efficiency (a) SR-DF$^{2-body}_{[110,120)}$, (b) SR-DF1$^{2-body}_{[120,140)}$, (c) SR-DF2$^{2-body}_{[140,160)}$, (d) SR-DF3$^{2-body}_{[160,180)}$, (e) SR-DF4$^{2-body}_{[180,220)}$, (f) SR-DF5$^{2-body}_{[220,\infty)}$ for a simplified model assuming $t \tilde{t} + \phi$ pair production.

Two-body selection efficiency (a) SR-DF$^{2-body}_{[110,120)}$, (b) SR-DF1$^{2-body}_{[120,140)}$, (c) SR-DF2$^{2-body}_{[140,160)}$, (d) SR-DF3$^{2-body}_{[160,180)}$, (e) SR-DF4$^{2-body}_{[180,220)}$, (f) SR-DF5$^{2-body}_{[220,\infty)}$ for a simplified model assuming $t \tilde{t} + \phi$ pair production.

Two-body selection efficiency (a) SR-SF$^{2-body}_{[110,120)}$, (b) SR-SF1$^{2-body}_{[120,140)}$, (c) SR-SF2$^{2-body}_{[140,160)}$, (d) SR-SF3$^{2-body}_{[160,180)}$, (e) SR-SF4$^{2-body}_{[180,220)}$, (f) SR-SF5$^{2-body}_{[220,\infty)}$ for a simplified model assuming $t \tilde{t} + \phi$ pair production.

Two-body selection efficiency (a) SR-SF$^{2-body}_{[110,120)}$, (b) SR-SF1$^{2-body}_{[120,140)}$, (c) SR-SF2$^{2-body}_{[140,160)}$, (d) SR-SF3$^{2-body}_{[160,180)}$, (e) SR-SF4$^{2-body}_{[180,220)}$, (f) SR-SF5$^{2-body}_{[220,\infty)}$ for a simplified model assuming $t \tilde{t} + \phi$ pair production.

Two-body selection efficiency (a) SR-SF$^{2-body}_{[110,120)}$, (b) SR-SF1$^{2-body}_{[120,140)}$, (c) SR-SF2$^{2-body}_{[140,160)}$, (d) SR-SF3$^{2-body}_{[160,180)}$, (e) SR-SF4$^{2-body}_{[180,220)}$, (f) SR-SF5$^{2-body}_{[220,\infty)}$ for a simplified model assuming $t \tilde{t} + \phi$ pair production.

Two-body selection efficiency (a) SR-SF$^{2-body}_{[110,120)}$, (b) SR-SF1$^{2-body}_{[120,140)}$, (c) SR-SF2$^{2-body}_{[140,160)}$, (d) SR-SF3$^{2-body}_{[160,180)}$, (e) SR-SF4$^{2-body}_{[180,220)}$, (f) SR-SF5$^{2-body}_{[220,\infty)}$ for a simplified model assuming $t \tilde{t} + \phi$ pair production.

Two-body selection efficiency (a) SR-SF$^{2-body}_{[110,120)}$, (b) SR-SF1$^{2-body}_{[120,140)}$, (c) SR-SF2$^{2-body}_{[140,160)}$, (d) SR-SF3$^{2-body}_{[160,180)}$, (e) SR-SF4$^{2-body}_{[180,220)}$, (f) SR-SF5$^{2-body}_{[220,\infty)}$ for a simplified model assuming $t \tilde{t} + \phi$ pair production.

Two-body selection efficiency (a) SR-SF$^{2-body}_{[110,120)}$, (b) SR-SF1$^{2-body}_{[120,140)}$, (c) SR-SF2$^{2-body}_{[140,160)}$, (d) SR-SF3$^{2-body}_{[160,180)}$, (e) SR-SF4$^{2-body}_{[180,220)}$, (f) SR-SF5$^{2-body}_{[220,\infty)}$ for a simplified model assuming $t \tilde{t} + \phi$ pair production.

Two-body selection efficiency (a) $SR^{2-body}_{[110,\infty)}$ , (b) $SR^{2-body}_{[120,\infty)}$ , (c) $SR^{2-body}_{[140,\infty)}$ , (d) $SR^{2-body}_{[160,\infty)}$ , (e) $SR^{2-body}_{[180,\infty)}$ , (f) $SR^{2-body}_{[200,\infty)}$ , (g) $SR^{2-body}_{[220,\infty)}$ for a simplified model assuming $t \tilde{t} + \phi$ pair production.

Two-body selection efficiency (a) $SR^{2-body}_{[110,\infty)}$ , (b) $SR^{2-body}_{[120,\infty)}$ , (c) $SR^{2-body}_{[140,\infty)}$ , (d) $SR^{2-body}_{[160,\infty)}$ , (e) $SR^{2-body}_{[180,\infty)}$ , (f) $SR^{2-body}_{[200,\infty)}$ , (g) $SR^{2-body}_{[220,\infty)}$ for a simplified model assuming $t \tilde{t} + \phi$ pair production.

Two-body selection efficiency (a) $SR^{2-body}_{[110,\infty)}$ , (b) $SR^{2-body}_{[120,\infty)}$ , (c) $SR^{2-body}_{[140,\infty)}$ , (d) $SR^{2-body}_{[160,\infty)}$ , (e) $SR^{2-body}_{[180,\infty)}$ , (f) $SR^{2-body}_{[200,\infty)}$ , (g) $SR^{2-body}_{[220,\infty)}$ for a simplified model assuming $t \tilde{t} + \phi$ pair production.

Two-body selection efficiency (a) $SR^{2-body}_{[110,\infty)}$ , (b) $SR^{2-body}_{[120,\infty)}$ , (c) $SR^{2-body}_{[140,\infty)}$ , (d) $SR^{2-body}_{[160,\infty)}$ , (e) $SR^{2-body}_{[180,\infty)}$ , (f) $SR^{2-body}_{[200,\infty)}$ , (g) $SR^{2-body}_{[220,\infty)}$ for a simplified model assuming $t \tilde{t} + \phi$ pair production.

Two-body selection efficiency (a) $SR^{2-body}_{[110,\infty)}$ , (b) $SR^{2-body}_{[120,\infty)}$ , (c) $SR^{2-body}_{[140,\infty)}$ , (d) $SR^{2-body}_{[160,\infty)}$ , (e) $SR^{2-body}_{[180,\infty)}$ , (f) $SR^{2-body}_{[200,\infty)}$ , (g) $SR^{2-body}_{[220,\infty)}$ for a simplified model assuming $t \tilde{t} + \phi$ pair production.

Two-body selection efficiency (a) $SR^{2-body}_{[110,\infty)}$ , (b) $SR^{2-body}_{[120,\infty)}$ , (c) $SR^{2-body}_{[140,\infty)}$ , (d) $SR^{2-body}_{[160,\infty)}$ , (e) $SR^{2-body}_{[180,\infty)}$ , (f) $SR^{2-body}_{[200,\infty)}$ , (g) $SR^{2-body}_{[220,\infty)}$ for a simplified model assuming $t \tilde{t} + \phi$ pair production.

Two-body selection efficiency (a) $SR^{2-body}_{[110,\infty)}$ , (b) $SR^{2-body}_{[120,\infty)}$ , (c) $SR^{2-body}_{[140,\infty)}$ , (d) $SR^{2-body}_{[160,\infty)}$ , (e) $SR^{2-body}_{[180,\infty)}$ , (f) $SR^{2-body}_{[200,\infty)}$ , (g) $SR^{2-body}_{[220,\infty)}$ for a simplified model assuming $t \tilde{t} + \phi$ pair production.

Two-body selection efficiency (a) SR-DF$^{2-body}_{[110,120)}$, (b) SR-DF1$^{2-body}_{[120,140)}$, (c) SR-DF2$^{2-body}_{[140,160)}$, (d) SR-DF3$^{2-body}_{[160,180)}$, (e) SR-DF4$^{2-body}_{[180,220)}$, (f) SR-DF5$^{2-body}_{[220,\infty)}$ for a simplified model assuming $t \tilde{t} + a$ pair production.

Two-body selection efficiency (a) SR-DF$^{2-body}_{[110,120)}$, (b) SR-DF1$^{2-body}_{[120,140)}$, (c) SR-DF2$^{2-body}_{[140,160)}$, (d) SR-DF3$^{2-body}_{[160,180)}$, (e) SR-DF4$^{2-body}_{[180,220)}$, (f) SR-DF5$^{2-body}_{[220,\infty)}$ for a simplified model assuming $t \tilde{t} + a$ pair production.

Two-body selection efficiency (a) SR-DF$^{2-body}_{[110,120)}$, (b) SR-DF1$^{2-body}_{[120,140)}$, (c) SR-DF2$^{2-body}_{[140,160)}$, (d) SR-DF3$^{2-body}_{[160,180)}$, (e) SR-DF4$^{2-body}_{[180,220)}$, (f) SR-DF5$^{2-body}_{[220,\infty)}$ for a simplified model assuming $t \tilde{t} + a$ pair production.

Two-body selection efficiency (a) SR-DF$^{2-body}_{[110,120)}$, (b) SR-DF1$^{2-body}_{[120,140)}$, (c) SR-DF2$^{2-body}_{[140,160)}$, (d) SR-DF3$^{2-body}_{[160,180)}$, (e) SR-DF4$^{2-body}_{[180,220)}$, (f) SR-DF5$^{2-body}_{[220,\infty)}$ for a simplified model assuming $t \tilde{t} + a$ pair production.

Two-body selection efficiency (a) SR-DF$^{2-body}_{[110,120)}$, (b) SR-DF1$^{2-body}_{[120,140)}$, (c) SR-DF2$^{2-body}_{[140,160)}$, (d) SR-DF3$^{2-body}_{[160,180)}$, (e) SR-DF4$^{2-body}_{[180,220)}$, (f) SR-DF5$^{2-body}_{[220,\infty)}$ for a simplified model assuming $t \tilde{t} + a$ pair production.

Two-body selection efficiency (a) SR-DF$^{2-body}_{[110,120)}$, (b) SR-DF1$^{2-body}_{[120,140)}$, (c) SR-DF2$^{2-body}_{[140,160)}$, (d) SR-DF3$^{2-body}_{[160,180)}$, (e) SR-DF4$^{2-body}_{[180,220)}$, (f) SR-DF5$^{2-body}_{[220,\infty)}$ for a simplified model assuming $t \tilde{t} + a$ pair production.

Two-body selection efficiency (a) SR-SF$^{2-body}_{[110,120)}$, (b) SR-SF1$^{2-body}_{[120,140)}$, (c) SR-SF2$^{2-body}_{[140,160)}$, (d) SR-SF3$^{2-body}_{[160,180)}$, (e) SR-SF4$^{2-body}_{[180,220)}$, (f) SR-SF5$^{2-body}_{[220,\infty)}$ for a simplified model assuming $t \tilde{t} + a$ pair production.

Two-body selection efficiency (a) SR-SF$^{2-body}_{[110,120)}$, (b) SR-SF1$^{2-body}_{[120,140)}$, (c) SR-SF2$^{2-body}_{[140,160)}$, (d) SR-SF3$^{2-body}_{[160,180)}$, (e) SR-SF4$^{2-body}_{[180,220)}$, (f) SR-SF5$^{2-body}_{[220,\infty)}$ for a simplified model assuming $t \tilde{t} + a$ pair production.

Two-body selection efficiency (a) SR-SF$^{2-body}_{[110,120)}$, (b) SR-SF1$^{2-body}_{[120,140)}$, (c) SR-SF2$^{2-body}_{[140,160)}$, (d) SR-SF3$^{2-body}_{[160,180)}$, (e) SR-SF4$^{2-body}_{[180,220)}$, (f) SR-SF5$^{2-body}_{[220,\infty)}$ for a simplified model assuming $t \tilde{t} + a$ pair production.

Two-body selection efficiency (a) SR-SF$^{2-body}_{[110,120)}$, (b) SR-SF1$^{2-body}_{[120,140)}$, (c) SR-SF2$^{2-body}_{[140,160)}$, (d) SR-SF3$^{2-body}_{[160,180)}$, (e) SR-SF4$^{2-body}_{[180,220)}$, (f) SR-SF5$^{2-body}_{[220,\infty)}$ for a simplified model assuming $t \tilde{t} + a$ pair production.

Two-body selection efficiency (a) SR-SF$^{2-body}_{[110,120)}$, (b) SR-SF1$^{2-body}_{[120,140)}$, (c) SR-SF2$^{2-body}_{[140,160)}$, (d) SR-SF3$^{2-body}_{[160,180)}$, (e) SR-SF4$^{2-body}_{[180,220)}$, (f) SR-SF5$^{2-body}_{[220,\infty)}$ for a simplified model assuming $t \tilde{t} + a$ pair production.

Two-body selection efficiency (a) SR-SF$^{2-body}_{[110,120)}$, (b) SR-SF1$^{2-body}_{[120,140)}$, (c) SR-SF2$^{2-body}_{[140,160)}$, (d) SR-SF3$^{2-body}_{[160,180)}$, (e) SR-SF4$^{2-body}_{[180,220)}$, (f) SR-SF5$^{2-body}_{[220,\infty)}$ for a simplified model assuming $t \tilde{t} + a$ pair production.

Two-body selection efficiency (a) $SR^{2-body}_{[110,\infty)}$ , (b) $SR^{2-body}_{[120,\infty)}$ , (c) $SR^{2-body}_{[140,\infty)}$ , (d) $SR^{2-body}_{[160,\infty)}$ , (e) $SR^{2-body}_{[180,\infty)}$ , (f) $SR^{2-body}_{[200,\infty)}$ , (g) $SR^{2-body}_{[220,\infty)}$ for a simplified model assuming $t \tilde{t} + a$ pair production.

Two-body selection efficiency (a) $SR^{2-body}_{[110,\infty)}$ , (b) $SR^{2-body}_{[120,\infty)}$ , (c) $SR^{2-body}_{[140,\infty)}$ , (d) $SR^{2-body}_{[160,\infty)}$ , (e) $SR^{2-body}_{[180,\infty)}$ , (f) $SR^{2-body}_{[200,\infty)}$ , (g) $SR^{2-body}_{[220,\infty)}$ for a simplified model assuming $t \tilde{t} + a$ pair production.

Two-body selection efficiency (a) $SR^{2-body}_{[110,\infty)}$ , (b) $SR^{2-body}_{[120,\infty)}$ , (c) $SR^{2-body}_{[140,\infty)}$ , (d) $SR^{2-body}_{[160,\infty)}$ , (e) $SR^{2-body}_{[180,\infty)}$ , (f) $SR^{2-body}_{[200,\infty)}$ , (g) $SR^{2-body}_{[220,\infty)}$ for a simplified model assuming $t \tilde{t} + a$ pair production.

Two-body selection efficiency (a) $SR^{2-body}_{[110,\infty)}$ , (b) $SR^{2-body}_{[120,\infty)}$ , (c) $SR^{2-body}_{[140,\infty)}$ , (d) $SR^{2-body}_{[160,\infty)}$ , (e) $SR^{2-body}_{[180,\infty)}$ , (f) $SR^{2-body}_{[200,\infty)}$ , (g) $SR^{2-body}_{[220,\infty)}$ for a simplified model assuming $t \tilde{t} + a$ pair production.

Two-body selection efficiency (a) $SR^{2-body}_{[110,\infty)}$ , (b) $SR^{2-body}_{[120,\infty)}$ , (c) $SR^{2-body}_{[140,\infty)}$ , (d) $SR^{2-body}_{[160,\infty)}$ , (e) $SR^{2-body}_{[180,\infty)}$ , (f) $SR^{2-body}_{[200,\infty)}$ , (g) $SR^{2-body}_{[220,\infty)}$ for a simplified model assuming $t \tilde{t} + a$ pair production.

Two-body selection efficiency (a) $SR^{2-body}_{[110,\infty)}$ , (b) $SR^{2-body}_{[120,\infty)}$ , (c) $SR^{2-body}_{[140,\infty)}$ , (d) $SR^{2-body}_{[160,\infty)}$ , (e) $SR^{2-body}_{[180,\infty)}$ , (f) $SR^{2-body}_{[200,\infty)}$ , (g) $SR^{2-body}_{[220,\infty)}$ for a simplified model assuming $t \tilde{t} + a$ pair production.

Two-body selection efficiency (a) $SR^{2-body}_{[110,\infty)}$ , (b) $SR^{2-body}_{[120,\infty)}$ , (c) $SR^{2-body}_{[140,\infty)}$ , (d) $SR^{2-body}_{[160,\infty)}$ , (e) $SR^{2-body}_{[180,\infty)}$ , (f) $SR^{2-body}_{[200,\infty)}$ , (g) $SR^{2-body}_{[220,\infty)}$ for a simplified model assuming $t \tilde{t} + a$ pair production.

Three-body selection efficiency (a) SR-DF$^{3-body}_{t}$, (b) SR-SF$^{3-body}_{t}$, (c) SR-DF$^{3-body}_{W}$, (d) SR-SF$^{3-body}_{W}$ for a simplified model assuming $ \tilde{t}_1$ pair production.

Three-body selection efficiency (a) SR-DF$^{3-body}_{t}$, (b) SR-SF$^{3-body}_{t}$, (c) SR-DF$^{3-body}_{W}$, (d) SR-SF$^{3-body}_{W}$ for a simplified model assuming $ \tilde{t}_1$ pair production.

Three-body selection efficiency (a) SR-DF$^{3-body}_{t}$, (b) SR-SF$^{3-body}_{t}$, (c) SR-DF$^{3-body}_{W}$, (d) SR-SF$^{3-body}_{W}$ for a simplified model assuming $ \tilde{t}_1$ pair production.

Three-body selection efficiency (a) SR-DF$^{3-body}_{t}$, (b) SR-SF$^{3-body}_{t}$, (c) SR-DF$^{3-body}_{W}$, (d) SR-SF$^{3-body}_{W}$ for a simplified model assuming $ \tilde{t}_1$ pair production.

Four-body selection Efficiency (a) SR$^{4-body}_{Small \Delta m}$ , (b) $SR^{4-body}_{Large \Delta m}$ for a simplified model assuming $\tilde{t}_1$ pair production.

Four-body selection Efficiency (a) SR$^{4-body}_{Small \Delta m}$ , (b) $SR^{4-body}_{Large \Delta\ m}$ for a simplified model assuming $\tilde{t}_1$ pair production.

Two-body selection acceptance (a) SR-DF$^{2-body}_{[110,120)}$, (b) SR-DF1$^{2-body}_{[120,140)}$, (c) SR-DF2$^{2-body}_{[140,160)}$, (d) SR-DF3$^{2-body}_{[160,180)}$, (e) SR-DF4$^{2-body}_{[180,220)}$, (f) SR-DF5$^{2-body}_{[220,\infty)}$ for a simplified model assuming $ \tilde{t}_1$ pair production.

Two-body selection acceptance (a) SR-DF$^{2-body}_{[110,120)}$, (b) SR-DF1$^{2-body}_{[120,140)}$, (c) SR-DF2$^{2-body}_{[140,160)}$, (d) SR-DF3$^{2-body}_{[160,180)}$, (e) SR-DF4$^{2-body}_{[180,220)}$, (f) SR-DF5$^{2-body}_{[220,\infty)}$ for a simplified model assuming $ \tilde{t}_1$ pair production.

Two-body selection acceptance (a) SR-DF$^{2-body}_{[110,120)}$, (b) SR-DF1$^{2-body}_{[120,140)}$, (c) SR-DF2$^{2-body}_{[140,160)}$, (d) SR-DF3$^{2-body}_{[160,180)}$, (e) SR-DF4$^{2-body}_{[180,220)}$, (f) SR-DF5$^{2-body}_{[220,\infty)}$ for a simplified model assuming $ \tilde{t}_1$ pair production.

Two-body selection acceptance (a) SR-DF$^{2-body}_{[110,120)}$, (b) SR-DF1$^{2-body}_{[120,140)}$, (c) SR-DF2$^{2-body}_{[140,160)}$, (d) SR-DF3$^{2-body}_{[160,180)}$, (e) SR-DF4$^{2-body}_{[180,220)}$, (f) SR-DF5$^{2-body}_{[220,\infty)}$ for a simplified model assuming $ \tilde{t}_1$ pair production.

Two-body selection acceptance (a) SR-DF$^{2-body}_{[110,120)}$, (b) SR-DF1$^{2-body}_{[120,140)}$, (c) SR-DF2$^{2-body}_{[140,160)}$, (d) SR-DF3$^{2-body}_{[160,180)}$, (e) SR-DF4$^{2-body}_{[180,220)}$, (f) SR-DF5$^{2-body}_{[220,\infty)}$ for a simplified model assuming $ \tilde{t}_1$ pair production.

Two-body selection acceptance (a) SR-DF$^{2-body}_{[110,120)}$, (b) SR-DF1$^{2-body}_{[120,140)}$, (c) SR-DF2$^{2-body}_{[140,160)}$, (d) SR-DF3$^{2-body}_{[160,180)}$, (e) SR-DF4$^{2-body}_{[180,220)}$, (f) SR-DF5$^{2-body}_{[220,\infty)}$ for a simplified model assuming $ \tilde{t}_1$ pair production.

Two-body selection acceptance (a) SR-SF$^{2-body}_{[110,120)}$, (b) SR-SF1$^{2-body}_{[120,140)}$, (c) SR-SF2$^{2-body}_{[140,160)}$, (d) SR-SF3$^{2-body}_{[160,180)}$, (e) SR-SF4$^{2-body}_{[180,220)}$, (f) SR-SF5$^{2-body}_{[220,\infty)}$ for a simplified model assuming $ \tilde{t}_1$ pair production.

Two-body selection acceptance (a) SR-SF$^{2-body}_{[110,120)}$, (b) SR-SF1$^{2-body}_{[120,140)}$, (c) SR-SF2$^{2-body}_{[140,160)}$, (d) SR-SF3$^{2-body}_{[160,180)}$, (e) SR-SF4$^{2-body}_{[180,220)}$, (f) SR-SF5$^{2-body}_{[220,\infty)}$ for a simplified model assuming $ \tilde{t}_1$ pair production.

Two-body selection acceptance (a) SR-SF$^{2-body}_{[110,120)}$, (b) SR-SF1$^{2-body}_{[120,140)}$, (c) SR-SF2$^{2-body}_{[140,160)}$, (d) SR-SF3$^{2-body}_{[160,180)}$, (e) SR-SF4$^{2-body}_{[180,220)}$, (f) SR-SF5$^{2-body}_{[220,\infty)}$ for a simplified model assuming $ \tilde{t}_1$ pair production.

Two-body selection acceptance (a) SR-SF$^{2-body}_{[110,120)}$, (b) SR-SF1$^{2-body}_{[120,140)}$, (c) SR-SF2$^{2-body}_{[140,160)}$, (d) SR-SF3$^{2-body}_{[160,180)}$, (e) SR-SF4$^{2-body}_{[180,220)}$, (f) SR-SF5$^{2-body}_{[220,\infty)}$ for a simplified model assuming $ \tilde{t}_1$ pair production.

Two-body selection acceptance (a) SR-SF$^{2-body}_{[110,120)}$, (b) SR-SF1$^{2-body}_{[120,140)}$, (c) SR-SF2$^{2-body}_{[140,160)}$, (d) SR-SF3$^{2-body}_{[160,180)}$, (e) SR-SF4$^{2-body}_{[180,220)}$, (f) SR-SF5$^{2-body}_{[220,\infty)}$ for a simplified model assuming $ \tilde{t}_1$ pair production.

Two-body selection acceptance (a) SR-SF$^{2-body}_{[110,120)}$, (b) SR-SF1$^{2-body}_{[120,140)}$, (c) SR-SF2$^{2-body}_{[140,160)}$, (d) SR-SF3$^{2-body}_{[160,180)}$, (e) SR-SF4$^{2-body}_{[180,220)}$, (f) SR-SF5$^{2-body}_{[220,\infty)}$ for a simplified model assuming $ \tilde{t}_1$ pair production.

Two-body selection acceptance (a) $SR^{2-body}_{[110,\infty)}$ , (b) $SR^{2-body}_{[120,\infty)}$ , (c) $SR^{2-body}_{[140,\infty)}$ , (d) $SR^{2-body}_{[160,\infty)}$ , (e) $SR^{2-body}_{[180,\infty)}$ , (f) $SR^{2-body}_{[200,\infty)}$ , (g) $SR^{2-body}_{[220,\infty)}$ for a simplified model assuming $ \tilde{t}_1$ pair production.

Two-body selection acceptance (a) $SR^{2-body}_{[110,\infty)}$ , (b) $SR^{2-body}_{[120,\infty)}$ , (c) $SR^{2-body}_{[140,\infty)}$ , (d) $SR^{2-body}_{[160,\infty)}$ , (e) $SR^{2-body}_{[180,\infty)}$ , (f) $SR^{2-body}_{[200,\infty)}$ , (g) $SR^{2-body}_{[220,\infty)}$ for a simplified model assuming $ \tilde{t}_1$ pair production.

Two-body selection acceptance (a) $SR^{2-body}_{[110,\infty)}$ , (b) $SR^{2-body}_{[120,\infty)}$ , (c) $SR^{2-body}_{[140,\infty)}$ , (d) $SR^{2-body}_{[160,\infty)}$ , (e) $SR^{2-body}_{[180,\infty)}$ , (f) $SR^{2-body}_{[200,\infty)}$ , (g) $SR^{2-body}_{[220,\infty)}$ for a simplified model assuming $ \tilde{t}_1$ pair production.

Two-body selection acceptance (a) $SR^{2-body}_{[110,\infty)}$ , (b) $SR^{2-body}_{[120,\infty)}$ , (c) $SR^{2-body}_{[140,\infty)}$ , (d) $SR^{2-body}_{[160,\infty)}$ , (e) $SR^{2-body}_{[180,\infty)}$ , (f) $SR^{2-body}_{[200,\infty)}$ , (g) $SR^{2-body}_{[220,\infty)}$ for a simplified model assuming $ \tilde{t}_1$ pair production.

Two-body selection acceptance (a) $SR^{2-body}_{[110,\infty)}$ , (b) $SR^{2-body}_{[120,\infty)}$ , (c) $SR^{2-body}_{[140,\infty)}$ , (d) $SR^{2-body}_{[160,\infty)}$ , (e) $SR^{2-body}_{[180,\infty)}$ , (f) $SR^{2-body}_{[200,\infty)}$ , (g) $SR^{2-body}_{[220,\infty)}$ for a simplified model assuming $ \tilde{t}_1$ pair production.

Two-body selection acceptance (a) $SR^{2-body}_{[110,\infty)}$ , (b) $SR^{2-body}_{[120,\infty)}$ , (c) $SR^{2-body}_{[140,\infty)}$ , (d) $SR^{2-body}_{[160,\infty)}$ , (e) $SR^{2-body}_{[180,\infty)}$ , (f) $SR^{2-body}_{[200,\infty)}$ , (g) $SR^{2-body}_{[220,\infty)}$ for a simplified model assuming $ \tilde{t}_1$ pair production.

Two-body selection acceptance (a) $SR^{2-body}_{[110,\infty)}$ , (b) $SR^{2-body}_{[120,\infty)}$ , (c) $SR^{2-body}_{[140,\infty)}$ , (d) $SR^{2-body}_{[160,\infty)}$ , (e) $SR^{2-body}_{[180,\infty)}$ , (f) $SR^{2-body}_{[200,\infty)}$ , (g) $SR^{2-body}_{[220,\infty)}$ for a simplified model assuming $ \tilde{t}_1$ pair production.

Two-body selection acceptance (a) SR-DF$^{2-body}_{[110,120)}$, (b) SR-DF1$^{2-body}_{[120,140)}$, (c) SR-DF2$^{2-body}_{[140,160)}$, (d) SR-DF3$^{2-body}_{[160,180)}$, (e) SR-DF4$^{2-body}_{[180,220)}$, (f) SR-DF5$^{2-body}_{[220,\infty)}$ for a simplified model assuming $t \tilde{t} + \phi$ pair production.

Two-body selection acceptance (a) SR-DF$^{2-body}_{[110,120)}$, (b) SR-DF1$^{2-body}_{[120,140)}$, (c) SR-DF2$^{2-body}_{[140,160)}$, (d) SR-DF3$^{2-body}_{[160,180)}$, (e) SR-DF4$^{2-body}_{[180,220)}$, (f) SR-DF5$^{2-body}_{[220,\infty)}$ for a simplified model assuming $t \tilde{t} + \phi$ pair production.

Two-body selection acceptance (a) SR-DF$^{2-body}_{[110,120)}$, (b) SR-DF1$^{2-body}_{[120,140)}$, (c) SR-DF2$^{2-body}_{[140,160)}$, (d) SR-DF3$^{2-body}_{[160,180)}$, (e) SR-DF4$^{2-body}_{[180,220)}$, (f) SR-DF5$^{2-body}_{[220,\infty)}$ for a simplified model assuming $t \tilde{t} + \phi$ pair production.

Two-body selection acceptance (a) SR-DF$^{2-body}_{[110,120)}$, (b) SR-DF1$^{2-body}_{[120,140)}$, (c) SR-DF2$^{2-body}_{[140,160)}$, (d) SR-DF3$^{2-body}_{[160,180)}$, (e) SR-DF4$^{2-body}_{[180,220)}$, (f) SR-DF5$^{2-body}_{[220,\infty)}$ for a simplified model assuming $t \tilde{t} + \phi$ pair production.

Two-body selection acceptance (a) SR-DF$^{2-body}_{[110,120)}$, (b) SR-DF1$^{2-body}_{[120,140)}$, (c) SR-DF2$^{2-body}_{[140,160)}$, (d) SR-DF3$^{2-body}_{[160,180)}$, (e) SR-DF4$^{2-body}_{[180,220)}$, (f) SR-DF5$^{2-body}_{[220,\infty)}$ for a simplified model assuming $t \tilde{t} + \phi$ pair production.

Two-body selection acceptance (a) SR-DF$^{2-body}_{[110,120)}$, (b) SR-DF1$^{2-body}_{[120,140)}$, (c) SR-DF2$^{2-body}_{[140,160)}$, (d) SR-DF3$^{2-body}_{[160,180)}$, (e) SR-DF4$^{2-body}_{[180,220)}$, (f) SR-DF5$^{2-body}_{[220,\infty)}$ for a simplified model assuming $t \tilde{t} + \phi$ pair production.

Two-body selection acceptance (a) SR-SF$^{2-body}_{[110,120)}$, (b) SR-SF1$^{2-body}_{[120,140)}$, (c) SR-SF2$^{2-body}_{[140,160)}$, (d) SR-SF3$^{2-body}_{[160,180)}$, (e) SR-SF4$^{2-body}_{[180,220)}$, (f) SR-SF5$^{2-body}_{[220,\infty)}$ for a simplified model assuming $t \tilde{t} +\phi$ pair production.

Two-body selection acceptance (a) SR-SF$^{2-body}_{[110,120)}$, (b) SR-SF1$^{2-body}_{[120,140)}$, (c) SR-SF2$^{2-body}_{[140,160)}$, (d) SR-SF3$^{2-body}_{[160,180)}$, (e) SR-SF4$^{2-body}_{[180,220)}$, (f) SR-SF5$^{2-body}_{[220,\infty)}$ for a simplified model assuming $ t \tilde{t} +\phi$ pair production.

Two-body selection acceptance (a) SR-SF$^{2-body}_{[110,120)}$, (b) SR-SF1$^{2-body}_{[120,140)}$, (c) SR-SF2$^{2-body}_{[140,160)}$, (d) SR-SF3$^{2-body}_{[160,180)}$, (e) SR-SF4$^{2-body}_{[180,220)}$, (f) SR-SF5$^{2-body}_{[220,\infty)}$ for a simplified model assuming $ t \tilde{t} +\phi$ pair production.

Two-body selection acceptance (a) SR-SF$^{2-body}_{[110,120)}$, (b) SR-SF1$^{2-body}_{[120,140)}$, (c) SR-SF2$^{2-body}_{[140,160)}$, (d) SR-SF3$^{2-body}_{[160,180)}$, (e) SR-SF4$^{2-body}_{[180,220)}$, (f) SR-SF5$^{2-body}_{[220,\infty)}$ for a simplified model assuming $ t \tilde{t} +\phi$ pair production.

Two-body selection acceptance (a) SR-SF$^{2-body}_{[110,120)}$, (b) SR-SF1$^{2-body}_{[120,140)}$, (c) SR-SF2$^{2-body}_{[140,160)}$, (d) SR-SF3$^{2-body}_{[160,180)}$, (e) SR-SF4$^{2-body}_{[180,220)}$, (f) SR-SF5$^{2-body}_{[220,\infty)}$ for a simplified model assuming $ t \tilde{t} +\phi$ pair production.

Two-body selection acceptance (a) SR-SF$^{2-body}_{[110,120)}$, (b) SR-SF1$^{2-body}_{[120,140)}$, (c) SR-SF2$^{2-body}_{[140,160)}$, (d) SR-SF3$^{2-body}_{[160,180)}$, (e) SR-SF4$^{2-body}_{[180,220)}$, (f) SR-SF5$^{2-body}_{[220,\infty)}$ for a simplified model assuming $ t \tilde{t} +\phi$ pair production.

Two-body selection acceptance (a) $SR^{2-body}_{[110,\infty)}$ , (b) $SR^{2-body}_{[120,\infty)}$ , (c) $SR^{2-body}_{[140,\infty)}$ , (d) $SR^{2-body}_{[160,\infty)}$ , (e) $SR^{2-body}_{[180,\infty)}$ , (f) $SR^{2-body}_{[200,\infty)}$ , (g) $SR^{2-body}_{[220,\infty)}$ for a simplified model assuming $t \tilde{t} + \phi$ pair production.

Two-body selection acceptance (a) $SR^{2-body}_{[110,\infty)}$ , (b) $SR^{2-body}_{[120,\infty)}$ , (c) $SR^{2-body}_{[140,\infty)}$ , (d) $SR^{2-body}_{[160,\infty)}$ , (e) $SR^{2-body}_{[180,\infty)}$ , (f) $SR^{2-body}_{[200,\infty)}$ , (g) $SR^{2-body}_{[220,\infty)}$ for a simplified model assuming $t \tilde{t} + \phi$ pair production.

Two-body selection acceptance (a) $SR^{2-body}_{[110,\infty)}$ , (b) $SR^{2-body}_{[120,\infty)}$ , (c) $SR^{2-body}_{[140,\infty)}$ , (d) $SR^{2-body}_{[160,\infty)}$ , (e) $SR^{2-body}_{[180,\infty)}$ , (f) $SR^{2-body}_{[200,\infty)}$ , (g) $SR^{2-body}_{[220,\infty)}$ for a simplified model assuming $t \tilde{t} + \phi$ pair production.

Two-body selection acceptance (a) $SR^{2-body}_{[110,\infty)}$ , (b) $SR^{2-body}_{[120,\infty)}$ , (c) $SR^{2-body}_{[140,\infty)}$ , (d) $SR^{2-body}_{[160,\infty)}$ , (e) $SR^{2-body}_{[180,\infty)}$ , (f) $SR^{2-body}_{[200,\infty)}$ , (g) $SR^{2-body}_{[220,\infty)}$ for a simplified model assuming $t \tilde{t} + \phi$ pair production.

Two-body selection acceptance (a) $SR^{2-body}_{[110,\infty)}$ , (b) $SR^{2-body}_{[120,\infty)}$ , (c) $SR^{2-body}_{[140,\infty)}$ , (d) $SR^{2-body}_{[160,\infty)}$ , (e) $SR^{2-body}_{[180,\infty)}$ , (f) $SR^{2-body}_{[200,\infty)}$ , (g) $SR^{2-body}_{[220,\infty)}$ for a simplified model assuming $t \tilde{t} + \phi$ pair production.

Two-body selection acceptance (a) $SR^{2-body}_{[110,\infty)}$ , (b) $SR^{2-body}_{[120,\infty)}$ , (c) $SR^{2-body}_{[140,\infty)}$ , (d) $SR^{2-body}_{[160,\infty)}$ , (e) $SR^{2-body}_{[180,\infty)}$ , (f) $SR^{2-body}_{[200,\infty)}$ , (g) $SR^{2-body}_{[220,\infty)}$ for a simplified model assuming $t \tilde{t} + \phi$ pair production.

Two-body selection acceptance (a) $SR^{2-body}_{[110,\infty)}$ , (b) $SR^{2-body}_{[120,\infty)}$ , (c) $SR^{2-body}_{[140,\infty)}$ , (d) $SR^{2-body}_{[160,\infty)}$ , (e) $SR^{2-body}_{[180,\infty)}$ , (f) $SR^{2-body}_{[200,\infty)}$ , (g) $SR^{2-body}_{[220,\infty)}$ for a simplified model assuming $t \tilde{t} + \phi$ pair production.

Two-body selection acceptance (a) SR-DF$^{2-body}_{[110,120)}$, (b) SR-DF1$^{2-body}_{[120,140)}$, (c) SR-DF2$^{2-body}_{[140,160)}$, (d) SR-DF3$^{2-body}_{[160,180)}$, (e) SR-DF4$^{2-body}_{[180,220)}$, (f) SR-DF5$^{2-body}_{[220,\infty)}$ for a simplified model assuming $t \tilde{t} + a$ pair production.

Two-body selection acceptance (a) SR-DF$^{2-body}_{[110,120)}$, (b) SR-DF1$^{2-body}_{[120,140)}$, (c) SR-DF2$^{2-body}_{[140,160)}$, (d) SR-DF3$^{2-body}_{[160,180)}$, (e) SR-DF4$^{2-body}_{[180,220)}$, (f) SR-DF5$^{2-body}_{[220,\infty)}$ for a simplified model assuming $t \tilde{t} + a$ pair production.

Two-body selection acceptance (a) SR-DF$^{2-body}_{[110,120)}$, (b) SR-DF1$^{2-body}_{[120,140)}$, (c) SR-DF2$^{2-body}_{[140,160)}$, (d) SR-DF3$^{2-body}_{[160,180)}$, (e) SR-DF4$^{2-body}_{[180,220)}$, (f) SR-DF5$^{2-body}_{[220,\infty)}$ for a simplified model assuming $t \tilde{t} + a$ pair production.

Two-body selection acceptance (a) SR-DF$^{2-body}_{[110,120)}$, (b) SR-DF1$^{2-body}_{[120,140)}$, (c) SR-DF2$^{2-body}_{[140,160)}$, (d) SR-DF3$^{2-body}_{[160,180)}$, (e) SR-DF4$^{2-body}_{[180,220)}$, (f) SR-DF5$^{2-body}_{[220,\infty)}$ for a simplified model assuming $t \tilde{t} + a$ pair production.

Two-body selection acceptance (a) SR-DF$^{2-body}_{[110,120)}$, (b) SR-DF1$^{2-body}_{[120,140)}$, (c) SR-DF2$^{2-body}_{[140,160)}$, (d) SR-DF3$^{2-body}_{[160,180)}$, (e) SR-DF4$^{2-body}_{[180,220)}$, (f) SR-DF5$^{2-body}_{[220,\infty)}$ for a simplified model assuming $t \tilde{t} + a$ pair production.

Two-body selection acceptance (a) SR-DF$^{2-body}_{[110,120)}$, (b) SR-DF1$^{2-body}_{[120,140)}$, (c) SR-DF2$^{2-body}_{[140,160)}$, (d) SR-DF3$^{2-body}_{[160,180)}$, (e) SR-DF4$^{2-body}_{[180,220)}$, (f) SR-DF5$^{2-body}_{[220,\infty)}$ for a simplified model assuming $t \tilde{t} + a$ pair production.

Two-body selection acceptance (a) SR-SF$^{2-body}_{[110,120)}$, (b) SR-SF1$^{2-body}_{[120,140)}$, (c) SR-SF2$^{2-body}_{[140,160)}$, (d) SR-SF3$^{2-body}_{[160,180)}$, (e) SR-SF4$^{2-body}_{[180,220)}$, (f) SR-SF5$^{2-body}_{[220,\infty)}$ for a simplified model assuming t \tilde{t} + a$ pair production.

Two-body selection acceptance (a) SR-SF$^{2-body}_{[110,120)}$, (b) SR-SF1$^{2-body}_{[120,140)}$, (c) SR-SF2$^{2-body}_{[140,160)}$, (d) SR-SF3$^{2-body}_{[160,180)}$, (e) SR-SF4$^{2-body}_{[180,220)}$, (f) SR-SF5$^{2-body}_{[220,\infty)}$ for a simplified model assuming $t \tilde{t} + a$ pair production.

Two-body selection acceptance (a) SR-SF$^{2-body}_{[110,120)}$, (b) SR-SF1$^{2-body}_{[120,140)}$, (c) SR-SF2$^{2-body}_{[140,160)}$, (d) SR-SF3$^{2-body}_{[160,180)}$, (e) SR-SF4$^{2-body}_{[180,220)}$, (f) SR-SF5$^{2-body}_{[220,\infty)}$ for a simplified model assuming $t \tilde{t} + a$ pair production.

Two-body selection acceptance (a) SR-SF$^{2-body}_{[110,120)}$, (b) SR-SF1$^{2-body}_{[120,140)}$, (c) SR-SF2$^{2-body}_{[140,160)}$, (d) SR-SF3$^{2-body}_{[160,180)}$, (e) SR-SF4$^{2-body}_{[180,220)}$, (f) SR-SF5$^{2-body}_{[220,\infty)}$ for a simplified model assuming $t \tilde{t} + a$ pair production.

Two-body selection acceptance (a) SR-SF$^{2-body}_{[110,120)}$, (b) SR-SF1$^{2-body}_{[120,140)}$, (c) SR-SF2$^{2-body}_{[140,160)}$, (d) SR-SF3$^{2-body}_{[160,180)}$, (e) SR-SF4$^{2-body}_{[180,220)}$, (f) SR-SF5$^{2-body}_{[220,\infty)}$ for a simplified model assuming $t \tilde{t} + a$ pair production.

Two-body selection acceptance (a) SR-SF$^{2-body}_{[110,120)}$, (b) SR-SF1$^{2-body}_{[120,140)}$, (c) SR-SF2$^{2-body}_{[140,160)}$, (d) SR-SF3$^{2-body}_{[160,180)}$, (e) SR-SF4$^{2-body}_{[180,220)}$, (f) SR-SF5$^{2-body}_{[220,\infty)}$ for a simplified model assuming $t \tilde{t} + a$ pair production.

Two-body selection acceptance (a) $SR^{2-body}_{[110,\infty)}$ , (b) $SR^{2-body}_{[120,\infty)}$ , (c) $SR^{2-body}_{[140,\infty)}$ , (d) $SR^{2-body}_{[160,\infty)}$ , (e) $SR^{2-body}_{[180,\infty)}$ , (f) $SR^{2-body}_{[200,\infty)}$ , (g) $SR^{2-body}_{[220,\infty)}$ for a simplified model assuming $ \tilde{t}_1$ pair production.

Two-body selection acceptance (a) $SR^{2-body}_{[110,\infty)}$ , (b) $SR^{2-body}_{[120,\infty)}$ , (c) $SR^{2-body}_{[140,\infty)}$ , (d) $SR^{2-body}_{[160,\infty)}$ , (e) $SR^{2-body}_{[180,\infty)}$ , (f) $SR^{2-body}_{[200,\infty)}$ , (g) $SR^{2-body}_{[220,\infty)}$ for a simplified model assuming $ \tilde{t}_1$ pair production.

Two-body selection acceptance (a) $SR^{2-body}_{[110,\infty)}$ , (b) $SR^{2-body}_{[120,\infty)}$ , (c) $SR^{2-body}_{[140,\infty)}$ , (d) $SR^{2-body}_{[160,\infty)}$ , (e) $SR^{2-body}_{[180,\infty)}$ , (f) $SR^{2-body}_{[200,\infty)}$ , (g) $SR^{2-body}_{[220,\infty)}$ for a simplified model assuming $ \tilde{t}_1$ pair production.

Two-body selection acceptance (a) $SR^{2-body}_{[110,\infty)}$ , (b) $SR^{2-body}_{[120,\infty)}$ , (c) $SR^{2-body}_{[140,\infty)}$ , (d) $SR^{2-body}_{[160,\infty)}$ , (e) $SR^{2-body}_{[180,\infty)}$ , (f) $SR^{2-body}_{[200,\infty)}$ , (g) $SR^{2-body}_{[220,\infty)}$ for a simplified model assuming $ \tilde{t}_1$ pair production.

Two-body selection acceptance (a) $SR^{2-body}_{[110,\infty)}$ , (b) $SR^{2-body}_{[120,\infty)}$ , (c) $SR^{2-body}_{[140,\infty)}$ , (d) $SR^{2-body}_{[160,\infty)}$ , (e) $SR^{2-body}_{[180,\infty)}$ , (f) $SR^{2-body}_{[200,\infty)}$ , (g) $SR^{2-body}_{[220,\infty)}$ for a simplified model assuming $ \tilde{t}_1$ pair production.

Two-body selection acceptance (a) $SR^{2-body}_{[110,\infty)}$ , (b) $SR^{2-body}_{[120,\infty)}$ , (c) $SR^{2-body}_{[140,\infty)}$ , (d) $SR^{2-body}_{[160,\infty)}$ , (e) $SR^{2-body}_{[180,\infty)}$ , (f) $SR^{2-body}_{[200,\infty)}$ , (g) $SR^{2-body}_{[220,\infty)}$ for a simplified model assuming $ \tilde{t}_1$ pair production.

Two-body selection acceptance (a) $SR^{2-body}_{[110,\infty)}$ , (b) $SR^{2-body}_{[120,\infty)}$ , (c) $SR^{2-body}_{[140,\infty)}$ , (d) $SR^{2-body}_{[160,\infty)}$ , (e) $SR^{2-body}_{[180,\infty)}$ , (f) $SR^{2-body}_{[200,\infty)}$ , (g) $SR^{2-body}_{[220,\infty)}$ for a simplified model assuming $ \tilde{t}_1$ pair production.

Three-body selection acceptance (a) SR-DF$^{3-body}_{t}$, (b) SR-SF$^{3-body}_{t}$, (c) SR-DF$^{3-body}_{W}$, (d) SR-SF$^{3-body}_{W}$ for a simplified model assuming $ \tilde{t}_1$ pair production.

Three-body selection acceptance (a) SR-DF$^{3-body}_{t}$, (b) SR-SF$^{3-body}_{t}$, (c) SR-DF$^{3-body}_{W}$, (d) SR-SF$^{3-body}_{W}$ for a simplified model assuming $ \tilde{t}_1$ pair production.

Three-body selection acceptance (a) SR-DF$^{3-body}_{t}$, (b) SR-SF$^{3-body}_{t}$, (c) SR-DF$^{3-body}_{W}$, (d) SR-SF$^{3-body}_{W}$ for a simplified model assuming $ \tilde{t}_1$ pair production.

Three-body selection acceptance (a) SR-DF$^{3-body}_{t}$, (b) SR-SF$^{3-body}_{t}$, (c) SR-DF$^{3-body}_{W}$, (d) SR-SF$^{3-body}_{W}$ for a simplified model assuming $ \tilde{t}_1$ pair production.

Four-body selection acceptance (a) SR$^{4-body}_{Small \Delta m}$ , (b) $SR^{4-body}_{Large \Delta m}$ for a simplified model assuming $\tilde{t}_1$ pair production.

Four-body selection acceptance (a) SR$^{4-body}_{Small \Delta m}$ , (b) $SR^{4-body}_{Large \Delta m}$ for a simplified model assuming $\tilde{t}_1$ pair production.

Two-body selection The numbers indicate the observed upper limits on the signal strenght for (a) a simplified model assuming $\tilde{t}_1$ pair production, (b) for $t\tilde{t} + a $ pseudoscalar models, (c) for $t\tilde{t} + \phi $ scalar models. In Figure (a), the red line corresponds to the observed limit.

Two-body selection The numbers indicate the observed upper limits on the signal strenght for (a) a simplified model assuming $\tilde{t}_1$ pair production, (b) for $t\tilde{t} + a $ pseudoscalar models, (c) for $t\tilde{t} + \phi $ scalar models. In Figure (a), the red line corresponds to the observed limit.

Two-body selection The numbers indicate the observed upper limits on the signal strenght for (a) a simplified model assuming $\tilde{t}_1$ pair production, (b) for $t\tilde{t} + a $ pseudoscalar models, (c) for $t\tilde{t} + \phi $ scalar models. In Figure (a), the red line corresponds to the observed limit.

Three-body selection The numbers indicate the upper limits on the signal strenght for a simplified model assuming $\tilde{t}_1$ pair production. For comparison, the red line corresponds to the observed limit.

Four-body selection The numbers indicate the upper limits on the signal strenght for a simplified model assuming $\tilde{t}_1$ pair production. For comparison, the red line corresponds to the observed limit.

Two-body selection The numbers indicate the upper limits on the signal cross-section for (a) a simplified model assuming $\tilde{t}_1$ pair production, (b) for $t\tilde{t} + a $ pseudoscalar models, (c) for $t\tilde{t} + \phi $ scalar models. In Figure (a), the red line corresponds to the observed limit.

Two-body selection The numbers indicate the upper limits on the signal cross-section for (a) a simplified model assuming $\tilde{t}_1$ pair production, (b) for $t\tilde{t} + a $ pseudoscalar models, (c) for $t\tilde{t} + \phi $ scalar models. In Figure (a), the red line corresponds to the observed limit.

Two-body selection The numbers indicate the upper limits on the signal cross-section for (a) a simplified model assuming $\tilde{t}_1$ pair production, (b) for $t\tilde{t} + a $ pseudoscalar models, (c) for $t\tilde{t} + \phi $ scalar models. In Figure (a), the red line corresponds to the observed limit.

Three-body selection The numbers indicate the upper limits on the signal cross-section for a simplified model assuming $\tilde{t}_1$ pair production. For comparison, the red line corresponds to the observed limit.

Four-body selection The numbers indicate the upper limits on the signal cross-section for a simplified model assuming $\tilde{t}_1$ pair production. For comparison, the red line corresponds to the observed limit.

Two-body selection. Background fit results for the $inclusive$ SRs. The Others category contains the contributions from $VVV$, $t\bar{t} t$, $t\bar{t}t\bar{t}$, $t\bar{t} W$, $t\bar{t} WW$, $t\bar{t} WZ$, $t\bar{t} H$, and $tZ$. Combined statistical and systematic uncertainties are given. Note that the individual uncertainties can be correlated, and do not necessarily add up quadratically to the total background uncertainty.

Cut flow for the simplified signal model $\tilde{t}_1 \rightarrow t^{(*)}\tilde{\chi}^0_1$ with $m(\tilde{t}_1)=600~ GeV$ and $m(\tilde{\chi}^0_1)=400~ GeV$ in the SRs for the two-body selection. The number of events is normalized to the cross-section and to an integrated luminosity of $139~fb^{-1}$.

Cut flow for the scalar signal model $t\bar{t} + \phi $ with $m(\phi)=150~ GeV$ and $m(\chi)=1~ GeV$ in the SRs for the two-body selection. The number of events is normalized to the cross-section and to an integrated luminosity of $139~fb^{-1}$.

Cut flow for the pseudoscalar signal model $t\bar{t} + a $ with $m(a)=150~ GeV$ and $m(\chi)=1~ GeV$ in the SRs for the two-body selection. The number of events is normalized to the cross-section and to an integrated luminosity of $139~fb^{-1}$.

Cut flow for the simplified signal model $\tilde{t}_1 \rightarrow bW\tilde{\chi}^0_1$ with $m(\tilde{t}_1)=550~ GeV$ and $m(\tilde{\chi}^0_1)=385~ GeV$ in the SRs for the three-body selection. The number of events is normalized to the cross-section and to an integrated luminosity of $139~fb^{-1}$.

Cut flow for the simplified signal model $\tilde{t}_1 \rightarrow bW\tilde{\chi}^0_1$ with $m(\tilde{t}_1)=550~ GeV$ and $m(\tilde{\chi}^0_1)=400~ GeV$ in the SRs for the three-body selection. The number of events is normalized to the cross-section and to an integrated luminosity of $139~fb^{-1}$.

Cut flow for the simplified signal model $\tilde{t}_1 \rightarrow bW\tilde{\chi}^0_1$ with $m(\tilde{t}_1)=550~ GeV$ and $m(\tilde{\chi}^0_1)=430~ GeV$ in the SRs for the three-body selection. The number of events is normalized to the cross-section and to an integrated luminosity of $139~fb^{-1}$.

Cut flow for the simplified signal model $\tilde{t}_1 \rightarrow bW\tilde{\chi}^0_1$ with $m(\tilde{t}_1)=550~ GeV$ and $m(\tilde{\chi}^0_1)=460~ GeV$ in the SRs for the three-body selection. The number of events is normalized to the cross-section and to an integrated luminosity of $139~fb^{-1}$.

Cut flow for the simplified signal model $\tilde{t}_1 \rightarrow b l \nu \tilde{\chi}^0_1$ with $m(\tilde{t}_1)=400~ GeV$ and $m(\tilde{\chi}^0_1)=380~ GeV$ in the SRs for the four-body selection. The number of events is normalized to the cross-section and to an integrated luminosity of $139~fb^{-1}$.

Cut flow for the simplified signal model $\tilde{t}_1 \rightarrow b l \nu \tilde{\chi}^0_1$ with $m(\tilde{t}_1)=460~ GeV$ and $m(\tilde{\chi}^0_1)=415~ GeV$ in the SRs for the four-body selection. The number of events is normalized to the cross-section and to an integrated luminosity of $139~fb^{-1}$.

Cut flow for the simplified signal model $\tilde{t}_1 \rightarrow b l \nu \tilde{\chi}^0_1$ with $m(\tilde{t}_1)=400~ GeV$ and $m(\tilde{\chi}^0_1)=320~ GeV$ in the SRs for the four-body selection. The number of events is normalized to the cross-section and to an integrated luminosity of $139~fb^{-1}$.

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

Expected (dashed black line) and observed (solid red line) 95% CL exclusion limits on the compressed SUSY simplified model with a bino-like LSP and wino-like NLSPs being considered. These are shown with $\pm1\sigma_\text{exp}$ (yellow band) from experimental systematic and statistical uncertainties, and with $\pm1\sigma^{\text{SUSY}}_{\text{theory}}$ (red dotted lines) from signal cross-section uncertainties, respectively. The limits set by the ATLAS searches using the soft lepton signature is illustrated by the blue region while the limit imposed by the LEP experiments is shown in grey.

Expected (dashed black line) and observed (solid red line) 95% CL exclusion limits on the compressed SUSY simplified model with a bino-like LSP and wino-like NLSPs being considered. These are shown with $\pm1\sigma_ ext{exp}$ (yellow band) from experimental systematic and statistical uncertainties, and with $\pm1\sigma^{ ext{SUSY}}_{ ext{theory}}$ (red dotted lines) from signal cross-section uncertainties, respectively. The limits set by the ATLAS searches using the soft lepton signature is illustrated by the blue region while the limit imposed by the LEP experiments is shown in grey.

The expected upper limits on the cross-sections at 95% CL for the simplified SUSY model featuring the production of mass-degenerate pairs of wino-like $\widetilde{\chi}_{2}^{0}$ and $\widetilde{\chi}_{1}^{\pm}$ particles in association with at least two jets. The production modes considered are $\widetilde{\chi}_{2}^{0}\widetilde{\chi}_{1}^{\pm}$, $\widetilde{\chi}_{1}^{\pm}\widetilde{\chi}_{1}^{\mp}$, $\widetilde{\chi}_{2}^{0}\widetilde{\chi}_{2}^{0}$, and $\widetilde{\chi}_{1}^{\pm}\widetilde{\chi}_{1}^{\pm}$.

The observed upper limits on the cross-sections at 95% CL for the simplified SUSY model featuring the production of mass-degenerate pairs of wino-like $\widetilde{\chi}_{2}^{0}$ and $\widetilde{\chi}_{1}^{\pm}$ particles in association with at least two jets. The production modes considered are $\widetilde{\chi}_{2}^{0}\widetilde{\chi}_{1}^{\pm}$, $\widetilde{\chi}_{1}^{\pm}\widetilde{\chi}_{1}^{\mp}$, $\widetilde{\chi}_{2}^{0}\widetilde{\chi}_{2}^{0}$, and $\widetilde{\chi}_{1}^{\pm}\widetilde{\chi}_{1}^{\pm}$.

Truth-level signal acceptances in $\text{SR}_\text{2j}$ with BDT score $\in [0.88, 1.0]$. All of the considered production modes ($\widetilde{\chi}_{2}^{0}\widetilde{\chi}_{1}^{\pm}$, $\widetilde{\chi}_{1}^{\pm}\widetilde{\chi}_{1}^{\mp}$, $\widetilde{\chi}_{2}^{0}\widetilde{\chi}_{2}^{0}$, and $\widetilde{\chi}_{1}^{\pm}\widetilde{\chi}_{1}^{\pm}$) are included. The acceptance is defined as the fraction of accepted events divided by the total number of events in the generator-level signal Monte Carlo simulation. The generator-level selections are the same as the ones at reconstruction level.

Truth-level signal acceptances in $\text{SR}_{\geq\text{3j}}$ with BDT score $\in [0.88, 1.0]$. All of the considered production modes ($\widetilde{\chi}_{2}^{0}\widetilde{\chi}_{1}^{\pm}$, $\widetilde{\chi}_{1}^{\pm}\widetilde{\chi}_{1}^{\mp}$, $\widetilde{\chi}_{2}^{0}\widetilde{\chi}_{2}^{0}$, and $\widetilde{\chi}_{1}^{\pm}\widetilde{\chi}_{1}^{\pm}$) are included. The acceptance is defined as the fraction of accepted events divided by the total number of events in the generator-level signal Monte Carlo simulation. The generator-level selections are the same as the ones at reconstruction level.

Signal efficiencies in $\text{SR}_\text{2j}$ with BDT score $\in [0.88, 1.0]$. All of the considered production modes ($\widetilde{\chi}_{2}^{0}\widetilde{\chi}_{1}^{\pm}$, $\widetilde{\chi}_{1}^{\pm}\widetilde{\chi}_{1}^{\mp}$, $\widetilde{\chi}_{2}^{0}\widetilde{\chi}_{2}^{0}$, and $\widetilde{\chi}_{1}^{\pm}\widetilde{\chi}_{1}^{\pm}$) are included. The efficiency is defined by the number of events of reconstructed-level signal simulation divided by the number of events obtained at generator level. The generator-level selections are the same as the ones at reconstruction level.

Signal efficiencies in $\text{SR}_{\geq\text{3j}}$ with BDT score $\in [0.88, 1.0]$. All of the considered production modes ($\widetilde{\chi}_{2}^{0}\widetilde{\chi}_{1}^{\pm}$, $\widetilde{\chi}_{1}^{\pm}\widetilde{\chi}_{1}^{\mp}$, $\widetilde{\chi}_{2}^{0}\widetilde{\chi}_{2}^{0}$, and $\widetilde{\chi}_{1}^{\pm}\widetilde{\chi}_{1}^{\pm}$) are included. The efficiency is defined by the number of events of reconstructed-level signal simulation divided by the number of events obtained at generator level. The generator-level selections are the same as the ones at reconstruction level.

Event selection cutflows for signal samples with $m(\tilde\chi^0_2) =100~\text{GeV}$ and $\Delta m(\tilde\chi^0_2 / \tilde\chi^\pm_1,\tilde\chi^0_1) = 0.2~\text{GeV}, 1.0~\text{GeV},$ and $5.0~\text{GeV}$. The first number in the column for each mass point corresponds to the signal event yield after applying the associated cut while the second number in parentheses represents the efficiency with respect to the previous cut. All of the production modes are included. The total cross-section used to obtain the initial number of events ($\sigma$) includes the cuts applied at parton level as described in the main text.