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.

Expected (dashed black line) and observed (solid red line) 95% CL exclusion limits on the higgsino simplified model being considered. These are shown with $\pm 1\sigma_{\mathrm{exp}}$ (yellow band) from experimental systematic and statistical uncertainties, and with $\pm 1\sigma_{\mathrm{theory}}^{\mathrm{SUSY}}$ (red dotted lines) from signal cross-section uncertainties, respectively. The limits set by the latest ATLAS searches using the soft lepton and disappearing track signatures are illustrated by the blue and green regions, respectively, while the limit imposed by the LEP experiments is shown in gray. The dot-dashed gray line indicates the predicted mass-splitting for the pure higgsino scenario.

Expected (dashed black line) and observed (solid red line) 95% CL exclusion limits on the higgsino simplified model being considered. These are shown with $\pm 1\sigma_{\mathrm{exp}}$ (yellow band) from experimental systematic and statistical uncertainties, and with $\pm 1\sigma_{\mathrm{theory}}^{\mathrm{SUSY}}$ (red dotted lines) from signal cross-section uncertainties, respectively. The limits set by the latest ATLAS searches using the soft lepton and disappearing track signatures are illustrated by the blue and green regions, respectively, while the limit imposed by the LEP experiments is shown in gray. The dot-dashed gray line indicates the predicted mass-splitting for the pure higgsino scenario.

Expected (dashed black line) and observed (solid red line) 95% CL exclusion limits on the higgsino simplified model being considered. These are shown with $\pm 1\sigma_{\mathrm{exp}}$ (yellow band) from experimental systematic and statistical uncertainties, and with $\pm 1\sigma_{\mathrm{theory}}^{\mathrm{SUSY}}$ (red dotted lines) from signal cross-section uncertainties, respectively. The limits set by the latest ATLAS searches using the soft lepton and disappearing track signatures are illustrated by the blue and green regions, respectively, while the limit imposed by the LEP experiments is shown in gray. The dot-dashed gray line indicates the predicted mass-splitting for the pure higgsino scenario.

Expected (dashed black line) and observed (solid red line) 95% CL exclusion limits on the higgsino simplified model being considered. These are shown with $\pm 1\sigma_{\mathrm{exp}}$ (yellow band) from experimental systematic and statistical uncertainties, and with $\pm 1\sigma_{\mathrm{theory}}^{\mathrm{SUSY}}$ (red dotted lines) from signal cross-section uncertainties, respectively. The limits set by the latest ATLAS searches using the soft lepton and disappearing track signatures are illustrated by the blue and green regions, respectively, while the limit imposed by the LEP experiments is shown in gray. The dot-dashed gray line indicates the predicted mass-splitting for the pure higgsino scenario.

Expected and observed CLs values per signal point represented by the grey numbers. The expected (dashed) and observed (solid) 95% CL exclusion limits are overlaid along with $\pm 1\sigma_{\mathrm{exp}}$ (yellow band) from experimental systematic and statistical uncertainties, and with $\pm 1\sigma_{\mathrm{theory}}^{\mathrm{SUSY}}$ (red dotted lines) from signal cross-section uncertainties, respectively.

Expected and observed CLs values per signal point represented by the grey numbers. The expected (dashed) and observed (solid) 95% CL exclusion limits are overlaid along with $\pm 1\sigma_{\mathrm{exp}}$ (yellow band) from experimental systematic and statistical uncertainties, and with $\pm 1\sigma_{\mathrm{theory}}^{\mathrm{SUSY}}$ (red dotted lines) from signal cross-section uncertainties, respectively.

Expected and observed cross-section upper-limit per signal point represented by the grey numbers. The expected (dashed) and observed (solid) 95% CL exclusion limits are overlaid along with $\pm 1\sigma_{\mathrm{exp}}$ (yellow band) from experimental systematic and statistical uncertainties, and with $\pm 1\sigma_{\mathrm{theory}}^{\mathrm{SUSY}}$ (red dotted lines) from signal cross-section uncertainties, respectively.

Expected and observed cross-section upper-limit per signal point represented by the grey numbers. The expected (dashed) and observed (solid) 95% CL exclusion limits are overlaid along with $\pm 1\sigma_{\mathrm{exp}}$ (yellow band) from experimental systematic and statistical uncertainties, and with $\pm 1\sigma_{\mathrm{theory}}^{\mathrm{SUSY}}$ (red dotted lines) from signal cross-section uncertainties, respectively.

Truth-level signal acceptances for each production process ($\tilde{\chi}_1^\pm \tilde{\chi}_1^0$, $\tilde{\chi}_1^\pm \tilde{\chi}_2^0$, $\tilde{\chi}_1^+ \tilde{\chi}_1^-$, and $\tilde{\chi}_2^0 \tilde{\chi}_1^0$) in a SR with the $S(d_0)$ requirement removed. 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, where the signal candidate track is identified as the charged particle with the largest distance between the interaction vertex and the secondary vertex of the higgsino decays.

Truth-level signal acceptances for each production process ($\tilde{\chi}_1^\pm \tilde{\chi}_1^0$, $\tilde{\chi}_1^\pm \tilde{\chi}_2^0$, $\tilde{\chi}_1^+ \tilde{\chi}_1^-$, and $\tilde{\chi}_2^0 \tilde{\chi}_1^0$) in a SR with the $S(d_0)$ requirement removed. 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, where the signal candidate track is identified as the charged particle with the largest distance between the interaction vertex and the secondary vertex of the higgsino decays.

Truth-level signal acceptances for each production process ($\tilde{\chi}_1^\pm \tilde{\chi}_1^0$, $\tilde{\chi}_1^\pm \tilde{\chi}_2^0$, $\tilde{\chi}_1^+ \tilde{\chi}_1^-$, and $\tilde{\chi}_2^0 \tilde{\chi}_1^0$) in a SR with the $S(d_0)$ requirement removed. 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, where the signal candidate track is identified as the charged particle with the largest distance between the interaction vertex and the secondary vertex of the higgsino decays.

Truth-level signal acceptances for each production process ($\tilde{\chi}_1^\pm \tilde{\chi}_1^0$, $\tilde{\chi}_1^\pm \tilde{\chi}_2^0$, $\tilde{\chi}_1^+ \tilde{\chi}_1^-$, and $\tilde{\chi}_2^0 \tilde{\chi}_1^0$) in a SR with the $S(d_0)$ requirement removed. 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, where the signal candidate track is identified as the charged particle with the largest distance between the interaction vertex and the secondary vertex of the higgsino decays.

Truth-level signal acceptances for each production process ($\tilde{\chi}_1^\pm \tilde{\chi}_1^0$, $\tilde{\chi}_1^\pm \tilde{\chi}_2^0$, $\tilde{\chi}_1^+ \tilde{\chi}_1^-$, and $\tilde{\chi}_2^0 \tilde{\chi}_1^0$) in a SR with the $S(d_0)$ requirement removed. 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, where the signal candidate track is identified as the charged particle with the largest distance between the interaction vertex and the secondary vertex of the higgsino decays.

Truth-level signal acceptances for each production process ($\tilde{\chi}_1^\pm \tilde{\chi}_1^0$, $\tilde{\chi}_1^\pm \tilde{\chi}_2^0$, $\tilde{\chi}_1^+ \tilde{\chi}_1^-$, and $\tilde{\chi}_2^0 \tilde{\chi}_1^0$) in a SR with the $S(d_0)$ requirement removed. 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, where the signal candidate track is identified as the charged particle with the largest distance between the interaction vertex and the secondary vertex of the higgsino decays.

Signal efficiencies in SR-Low for each production process ($\tilde{\chi}_1^\pm \tilde{\chi}_1^0$, $\tilde{\chi}_1^\pm \tilde{\chi}_2^0$, $\tilde{\chi}_1^+ \tilde{\chi}_1^-$, and $\tilde{\chi}_2^0 \tilde{\chi}_1^0$), defined by the number of events of reconstructed-level signal simulation divided by the number of events obtained at generator level, where the $S(d_0)$ selecton efficiency has the largest impact. The higgsino decay products from $\Delta \mathrm{m}(\tilde{\chi}_1^\pm,\tilde{\chi}_1^0) < 0.4$ GeV signal have $p_{\mathrm{T}}$ too low to be reconstructed as the signal candidate tracks, and therefore the identified signal candidate tracks are typically from pile-up collisions or underlying events similar to the QCD track background, causing a low $S(d_0)$ selection efficiency in these plots.

Signal efficiencies in SR-Low for each production process ($\tilde{\chi}_1^\pm \tilde{\chi}_1^0$, $\tilde{\chi}_1^\pm \tilde{\chi}_2^0$, $\tilde{\chi}_1^+ \tilde{\chi}_1^-$, and $\tilde{\chi}_2^0 \tilde{\chi}_1^0$), defined by the number of events of reconstructed-level signal simulation divided by the number of events obtained at generator level, where the $S(d_0)$ selecton efficiency has the largest impact. The higgsino decay products from $\Delta \mathrm{m}(\tilde{\chi}_1^\pm,\tilde{\chi}_1^0) < 0.4$ GeV signal have $p_{\mathrm{T}}$ too low to be reconstructed as the signal candidate tracks, and therefore the identified signal candidate tracks are typically from pile-up collisions or underlying events similar to the QCD track background, causing a low $S(d_0)$ selection efficiency in these plots.

Signal efficiencies in SR-Low for each production process ($\tilde{\chi}_1^\pm \tilde{\chi}_1^0$, $\tilde{\chi}_1^\pm \tilde{\chi}_2^0$, $\tilde{\chi}_1^+ \tilde{\chi}_1^-$, and $\tilde{\chi}_2^0 \tilde{\chi}_1^0$), defined by the number of events of reconstructed-level signal simulation divided by the number of events obtained at generator level, where the $S(d_0)$ selecton efficiency has the largest impact. The higgsino decay products from $\Delta \mathrm{m}(\tilde{\chi}_1^\pm,\tilde{\chi}_1^0) < 0.4$ GeV signal have $p_{\mathrm{T}}$ too low to be reconstructed as the signal candidate tracks, and therefore the identified signal candidate tracks are typically from pile-up collisions or underlying events similar to the QCD track background, causing a low $S(d_0)$ selection efficiency in these plots.

Signal efficiencies in SR-Low for each production process ($\tilde{\chi}_1^\pm \tilde{\chi}_1^0$, $\tilde{\chi}_1^\pm \tilde{\chi}_2^0$, $\tilde{\chi}_1^+ \tilde{\chi}_1^-$, and $\tilde{\chi}_2^0 \tilde{\chi}_1^0$), defined by the number of events of reconstructed-level signal simulation divided by the number of events obtained at generator level, where the $S(d_0)$ selecton efficiency has the largest impact. The higgsino decay products from $\Delta \mathrm{m}(\tilde{\chi}_1^\pm,\tilde{\chi}_1^0) < 0.4$ GeV signal have $p_{\mathrm{T}}$ too low to be reconstructed as the signal candidate tracks, and therefore the identified signal candidate tracks are typically from pile-up collisions or underlying events similar to the QCD track background, causing a low $S(d_0)$ selection efficiency in these plots.

Signal efficiencies in SR-Low for each production process ($\tilde{\chi}_1^\pm \tilde{\chi}_1^0$, $\tilde{\chi}_1^\pm \tilde{\chi}_2^0$, $\tilde{\chi}_1^+ \tilde{\chi}_1^-$, and $\tilde{\chi}_2^0 \tilde{\chi}_1^0$), defined by the number of events of reconstructed-level signal simulation divided by the number of events obtained at generator level, where the $S(d_0)$ selecton efficiency has the largest impact. The higgsino decay products from $\Delta \mathrm{m}(\tilde{\chi}_1^\pm,\tilde{\chi}_1^0) < 0.4$ GeV signal have $p_{\mathrm{T}}$ too low to be reconstructed as the signal candidate tracks, and therefore the identified signal candidate tracks are typically from pile-up collisions or underlying events similar to the QCD track background, causing a low $S(d_0)$ selection efficiency in these plots.

Signal efficiencies in SR-Low for each production process ($\tilde{\chi}_1^\pm \tilde{\chi}_1^0$, $\tilde{\chi}_1^\pm \tilde{\chi}_2^0$, $\tilde{\chi}_1^+ \tilde{\chi}_1^-$, and $\tilde{\chi}_2^0 \tilde{\chi}_1^0$), defined by the number of events of reconstructed-level signal simulation divided by the number of events obtained at generator level, where the $S(d_0)$ selecton efficiency has the largest impact. The higgsino decay products from $\Delta \mathrm{m}(\tilde{\chi}_1^\pm,\tilde{\chi}_1^0) < 0.4$ GeV signal have $p_{\mathrm{T}}$ too low to be reconstructed as the signal candidate tracks, and therefore the identified signal candidate tracks are typically from pile-up collisions or underlying events similar to the QCD track background, causing a low $S(d_0)$ selection efficiency in these plots.

Signal efficiencies in SR-High for each production process ($\tilde{\chi}_1^\pm \tilde{\chi}_1^0$, $\tilde{\chi}_1^\pm \tilde{\chi}_2^0$, $\tilde{\chi}_1^+ \tilde{\chi}_1^-$, and $\tilde{\chi}_2^0 \tilde{\chi}_1^0$), defined by the number of events of reconstructed-level signal simulation divided by the number of events obtained at generator level, where the $S(d_0)$ selecton efficiency has the largest impact. The higgsino decay products from $\Delta \mathrm{m}(\tilde{\chi}_1^\pm,\tilde{\chi}_1^0) < 0.4$ GeV signal have $p_{\mathrm{T}}$ too low to be reconstructed as the signal candidate tracks, and therefore the identified signal candidate tracks are typically from pile-up collisions or underlying events similar to the QCD track background, causing a low $S(d_0)$ selection efficiency in these plots.

Signal efficiencies in SR-High for each production process ($\tilde{\chi}_1^\pm \tilde{\chi}_1^0$, $\tilde{\chi}_1^\pm \tilde{\chi}_2^0$, $\tilde{\chi}_1^+ \tilde{\chi}_1^-$, and $\tilde{\chi}_2^0 \tilde{\chi}_1^0$), defined by the number of events of reconstructed-level signal simulation divided by the number of events obtained at generator level, where the $S(d_0)$ selecton efficiency has the largest impact. The higgsino decay products from $\Delta \mathrm{m}(\tilde{\chi}_1^\pm,\tilde{\chi}_1^0) < 0.4$ GeV signal have $p_{\mathrm{T}}$ too low to be reconstructed as the signal candidate tracks, and therefore the identified signal candidate tracks are typically from pile-up collisions or underlying events similar to the QCD track background, causing a low $S(d_0)$ selection efficiency in these plots.

Signal efficiencies in SR-High for each production process ($\tilde{\chi}_1^\pm \tilde{\chi}_1^0$, $\tilde{\chi}_1^\pm \tilde{\chi}_2^0$, $\tilde{\chi}_1^+ \tilde{\chi}_1^-$, and $\tilde{\chi}_2^0 \tilde{\chi}_1^0$), defined by the number of events of reconstructed-level signal simulation divided by the number of events obtained at generator level, where the $S(d_0)$ selecton efficiency has the largest impact. The higgsino decay products from $\Delta \mathrm{m}(\tilde{\chi}_1^\pm,\tilde{\chi}_1^0) < 0.4$ GeV signal have $p_{\mathrm{T}}$ too low to be reconstructed as the signal candidate tracks, and therefore the identified signal candidate tracks are typically from pile-up collisions or underlying events similar to the QCD track background, causing a low $S(d_0)$ selection efficiency in these plots.

Signal efficiencies in SR-High for each production process ($\tilde{\chi}_1^\pm \tilde{\chi}_1^0$, $\tilde{\chi}_1^\pm \tilde{\chi}_2^0$, $\tilde{\chi}_1^+ \tilde{\chi}_1^-$, and $\tilde{\chi}_2^0 \tilde{\chi}_1^0$), defined by the number of events of reconstructed-level signal simulation divided by the number of events obtained at generator level, where the $S(d_0)$ selecton efficiency has the largest impact. The higgsino decay products from $\Delta \mathrm{m}(\tilde{\chi}_1^\pm,\tilde{\chi}_1^0) < 0.4$ GeV signal have $p_{\mathrm{T}}$ too low to be reconstructed as the signal candidate tracks, and therefore the identified signal candidate tracks are typically from pile-up collisions or underlying events similar to the QCD track background, causing a low $S(d_0)$ selection efficiency in these plots.

Signal efficiencies in SR-High for each production process ($\tilde{\chi}_1^\pm \tilde{\chi}_1^0$, $\tilde{\chi}_1^\pm \tilde{\chi}_2^0$, $\tilde{\chi}_1^+ \tilde{\chi}_1^-$, and $\tilde{\chi}_2^0 \tilde{\chi}_1^0$), defined by the number of events of reconstructed-level signal simulation divided by the number of events obtained at generator level, where the $S(d_0)$ selecton efficiency has the largest impact. The higgsino decay products from $\Delta \mathrm{m}(\tilde{\chi}_1^\pm,\tilde{\chi}_1^0) < 0.4$ GeV signal have $p_{\mathrm{T}}$ too low to be reconstructed as the signal candidate tracks, and therefore the identified signal candidate tracks are typically from pile-up collisions or underlying events similar to the QCD track background, causing a low $S(d_0)$ selection efficiency in these plots.

Signal efficiencies in SR-High for each production process ($\tilde{\chi}_1^\pm \tilde{\chi}_1^0$, $\tilde{\chi}_1^\pm \tilde{\chi}_2^0$, $\tilde{\chi}_1^+ \tilde{\chi}_1^-$, and $\tilde{\chi}_2^0 \tilde{\chi}_1^0$), defined by the number of events of reconstructed-level signal simulation divided by the number of events obtained at generator level, where the $S(d_0)$ selecton efficiency has the largest impact. The higgsino decay products from $\Delta \mathrm{m}(\tilde{\chi}_1^\pm,\tilde{\chi}_1^0) < 0.4$ GeV signal have $p_{\mathrm{T}}$ too low to be reconstructed as the signal candidate tracks, and therefore the identified signal candidate tracks are typically from pile-up collisions or underlying events similar to the QCD track background, causing a low $S(d_0)$ selection efficiency in these plots.

Event selection cutflows for signal samples with $m(\tilde{\chi}_{1}^0)$ = 150 GeV and $\Delta m(\tilde{\chi}_{1}^\pm, \tilde{\chi}_{1}^0)$ = 1.5, 1.0, and 0.75 GeV, including all six production processes ($\tilde{\chi}_1^\pm \tilde{\chi}_1^0$, $\tilde{\chi}_1^\pm \tilde{\chi}_2^0$, $\tilde{\chi}_1^+ \tilde{\chi}_1^-$, and $\tilde{\chi}_2^0 \tilde{\chi}_1^0$). The cross-section used to obtain the initial number of events ($\sigma(\mathrm{n}_{\mathrm{jets}}) \geq 1$) refers to an emission of at least one gluon or quark with $p_{\mathrm{T}} > 50$ GeV at the parton level.

Event selection cutflows for signal samples with $m(\tilde{\chi}_{1}^0)$ = 150 GeV and $\Delta m(\tilde{\chi}_{1}^\pm, \tilde{\chi}_{1}^0)$ = 0.5, 0.35, and 0.25 GeV, including all six production processes ($\tilde{\chi}_1^\pm \tilde{\chi}_1^0$, $\tilde{\chi}_1^\pm \tilde{\chi}_2^0$, $\tilde{\chi}_1^+ \tilde{\chi}_1^-$, and $\tilde{\chi}_2^0 \tilde{\chi}_1^0$). The cross-section used to obtain the initial number of events ($\sigma(\mathrm{n}_{\mathrm{jets}}) \geq 1$) refers to an emission of at least one gluon or quark with $p_{\mathrm{T}} > 50$ GeV at the parton level.

Numbers of observed and background events for SR3L low for each bin of the distribution in Meff. The table corresponds to Fig. 4(d). The statistical and systematic uncertainties are combined for the predicted numbers.

Numbers of observed and background events for SR3L high for each bin of the distribution in Meff. The table corresponds to Fig. 4(e). The statistical and systematic uncertainties are combined for the predicted numbers.

The efficiencies are calculated for all simplified models. For each model, the values are given for the five signal regions and their combination. The simplified model is for direct production of squarks that decay into two steps into q q W Z W Z chi1^0 chi1^0 (see Fig. 6c in the paper).

The efficiencies are calculated for all simplified extra dimension model (see Fig. 8d in the paper). For each model, the values are given for the five signal regions and their combination.

The efficiencies are calculated for all simplified models. For each model, the values are given for the five signal regions and their combination. The simplified model is for direct production of gluinos that decay via sleptons into q q q q l l (l l) chi1^0 chi1^0 + neutrinos (see Fig. 6d in the paper).

The efficiencies are calculated for all simplified models. For each model, the values are given for the five signal regions and their combination. The simplified model is for direct production of gluinos that decay into q q q q W W chi1^0 chi1^0 (see Fig. 6a in the paper).

The efficiencies are calculated for all simplified models. For each model, the values are given for the five signal regions and their combination. The simplified model is for direct production of gluinos that decay into t tbar t tbar chi1^0 chi1^0 (see Fig. 5a in the paper). This particular model assumes that top quark is much heavier than gluino.

The efficiencies are calculated for all mSUGRA models (see Fig. 8a in the paper). For each model, the values are given for the five signal regions and their combination. The model assumes tan(beta)=30, A0=2m0, and mu>0.

The efficiencies are calculated for all simplified models. For each model, the values are given for the five signal regions and their combination. The simplified model is for direct production of gluinos. A gluino decays into t c chi1^0 (see Fig. 5c in the paper). This particular model assumes that m(chi1^0) = m(stop) - 20 GeV.

The efficiencies are calculated for all GMSB models (see Fig. 8c in the paper). For each model, the values are given for the five signal regions and their combination. The model assumes mmess=250 TeV, m5=3, mu>0, and Cgrav=1.

The efficiencies are calculated for all simplified models. For each model, the values are given for the five signal regions and their combination. The simplified model is for direct production of bottom squarks. A bottom squark decays into t chi1^(+-) and chi1^(+-) --> W^(+-) chi1^0 (see Fig. 7a in the paper). This particular model assumes that m(chi1^0)=60 GeV.

The efficiencies are calculated for all simplified models. For each model, the values are given for the five signal regions and their combination. The simplified model is for direct production of gluinos and top squarks. Top squarks undergo R-parity violating decays into b s and gluinos decay into t stop (see Fig. 5d in the paper).

The efficiencies are calculated for all simplified models. For each model, the values are given for the five signal regions and their combination. The simplified model is for direct production of bottom squarks. A bottom squark decays into t chi1^(+-) and chi1^(+-) --> W^(+-) chi1^0 (see Fig. 7b in the paper). This particular model assumes that m(chi1^0)=2(chi1^0).

The efficiencies are calculated for all mSUGRA/CMSSM models with bRPV (see Fig. 8b in the paper). For each model, the values are given for the five signal regions and their combination. The model assumes tan(beta)=30, A0=2m0, mu>0, and bRPV.

The efficiencies are calculated for all simplified models. For each model, the values are given for the five signal regions and their combination. The simplified model is for direct production of squarks. Squarks decay into q q l l (l l) chi1^0 chi1^0 + neutrinos (see Fig. 6e in the paper).

The efficiencies are calculated for all simplified models. For each model, the values are given for the five signal regions and their combination. The simplified model is for direct pair-production of gluinos that decay via a two-step process into q q q q W Z W Z chi1^0 chi1^0 (see Fig. 6b in the paper).

The efficiencies are calculated for all simplified models. For each model, the values are given for the five signal regions and their combination. The simplified model is for direct pair production of gluinos. A gluino decays into t stop. Consequently, a top squark squark decays into b chi1^(+-) and chi1^(+-) --> W^(+-) chi1^0 (see Fig. 5b in the paper). This particular model assumes that m(stop) < m(gluino), m(chi1^0)=6 GeV, and m(chi1^(+-))=118 GeV.

The acceptances (in percent, %) are calculated for all simplified models. For each model, the values are given for the five signal regions and their combination. The simplified model is for direct production of squarks that decay into two steps into q q W Z W Z chi1^0 chi1^0 (see Fig. 6c in the paper).

The acceptances (in percent, %) are calculated for all simplified extra dimension model (see Fig. 8d in the paper). For each model, the values are given for the five signal regions and their combination.

The acceptances (in percent, %) are calculated for all simplified models. For each model, the values are given for the five signal regions and their combination. The simplified model is for direct production of gluinos that decay via sleptons into q q q q l l (l l) chi1^0 chi1^0 + neutrinos (see Fig. 6d in the paper).

The acceptances (in percent, %) are calculated for all simplified models. For each model, the values are given for the five signal regions and their combination. The simplified model is for direct production of gluinos that decay into q q q q W W chi1^0 chi1^0 (see Fig. 6a in the paper).

The acceptances (in percent, %) are calculated for all simplified models. For each model, the values are given for the five signal regions and their combination. The simplified model is for direct production of gluinos that decay into t tbar t tbar chi1^0 chi1^0 (see Fig. 5a in the paper). This particular model assumes that top quark is much heavier than gluino.

The acceptances (in percent, %) are calculated for all mSUGRA models (see Fig. 8a in the paper). For each model, the values are given for the five signal regions and their combination. The model assumes tan(beta)=30, A0=2m0, and mu>0.

The acceptances (in percent, %) are calculated for all simplified models. For each model, the values are given for the five signal regions and their combination. The simplified model is for direct production of gluinos. A gluino decays into t c chi1^0 (see Fig. 5c in the paper). This particular model assumes that m(chi1^0) = m(stop) - 20 GeV.

The acceptances (in percent, %) are calculated for all GMSB models (see Fig. 8c in the paper). For each model, the values are given for the five signal regions and their combination. The model assumes mmess=250 TeV, m5=3, mu>0, and Cgrav=1.

The acceptances (in percent, %) are calculated for all simplified models. For each model, the values are given for the five signal regions and their combination. The simplified model is for direct production of bottom squarks. A bottom squark decays into t chi1^(+-) and chi1^(+-) --> W^(+-) chi1^0 (see Fig. 7a in the paper). This particular model assumes that m(chi1^0)=60 GeV.

The acceptances (in percent, %) are calculated for all simplified models. For each model, the values are given for the five signal regions and their combination. The simplified model is for direct production of gluinos and top squarks. Top squarks undergo R-parity violating decays into bs and gluinos decay into t stop (see Fig. 5d in the paper).

The acceptances (in percent, %) are calculated for all simplified models. For each model, the values are given for the five signal regions and their combination. The simplified model is for direct production of bottom squarks. A bottom squark decays into t chi1^(+-) and chi1^(+-) --> W chi1^0 (see Fig. 7b in the paper). This particular model assumes that m(chi1^0)=2(chi1^0).

The acceptances (in percent, %) are calculated for all mSUGRA/CMSSM models with bRPV (see Fig. 8b in the paper). For each model, the values are given for the five signal regions and their combination. The model assumes tan(beta)=30, A0=2m0, mu>0, and bRPV.

The acceptances (in percent, %) are calculated for all simplified models. For each model, the values are given for the five signal regions and their combination. The simplified model is for direct production of squarks. Squarks decay into q q l l (l l) chi1^0 chi1^0 + neutrinos (see Fig. 6e in the paper).

The acceptances (in percent, %) are calculated for all simplified models. For each model, the values are given for the five signal regions and their combination. The simplified model is for direct pair-production of gluinos that decay via a two-step process into q q q q W Z W Z chi1^0 chi1^0 (see Fig. 6b in the paper).

The acceptances (in percent, %) are calculated for all simplified models. For each model, the values are given for the five signal regions and their combination. The simplified model is for direct pair production of gluinos. A gluino decays into t stop. Consequently, a top squark squark decays into b chi1^(+-) and chi1^(+-) --> W^(+-) chi1^0 (see Fig. 5b in the paper). This particular model assumes that m(stop) < m(gluino), m(chi1^0)=6 GeV, and m(chi1^(+-))=118 GeV.

The limits on observed cross section are calculated for all simplified models. The simplified models are for direct pair production of squarks that decay into two steps into q q W Z W Z chi1^0 chi1^0 (see Fig. 6c in the paper).

The limits on observed cross sections are calculated for all simplified models. The simplified models are for direct pair-production of gluinos that decay via sleptons into q q q q l l (l l) chi1^0 chi1^0 + neutrinos (see Fig. 6d in the paper).

The limits on observed cross sections are calculated for all simplified models. The simplified models are for direct production of gluinos that decay into q q q q W W chi1^0 chi1^0 (see Fig. 6a in the paper).

The limits on observed cross sections are calculated for all simplified models. The simplified models are for direct production of gluinos that decay into t tbar t tbar chi1^0 chi1^0 (see Fig. 5a in the paper). This particular model assumes that top quark is much heavier than gluino.

The limits on observed cross sections are calculated for all simplified models. The simplified models are for direct pair production of gluinos. A gluino decays into t c chi1^0 (see Fig. 5c in the paper). This particular model assumes that m(chi1^0) = m(stop) - 20 GeV.

The limits on observed cross sections are calculated for all simplified models. The simplified models are for direct production of bottom squarks. A bottom squark decays into t chi1^(+-) and chi1^(+-) --> W^(+-) chi1^0 (see Fig. 7a in the paper). This particular model assumes that m(chi1^0)=60 GeV.

The limits on observed cross sections are calculated for all simplified models. The simplified models are for direct production of gluinos and top squarks. Top squarks undergo R-parity violating decays into bs and gluinos decay into t stop (see Fig. 5d in the paper).

The limits on observed cross sections are calculated for all simplified models. The simplified models are for direct production of bottom squarks. A bottom squark decays into t chi1^(+-) and chi1^(+-) --> W^(+-) chi1^0 (see Fig. 7b in the paper). This particular model assumes that m(chi1^0)=2(chi1^0).

The limits on observed cross sections are calculated for all simplified models. The simplified models are for direct production of squarks. Squarks decay into q q l l (l l) chi1^0 chi1^0 + neutrinos (see Fig. 6e in the paper).

The limits on observed cross sections are calculated for all simplified models. The simplified models are for direct pair-production of gluinos that decay via a two-step process into q q q q W Z W Z chi1^0 chi1^0 (see Fig. 6b in the paper).

The limits on observed cross sections are calculated for all simplified models. The simplified models are for direct pair production of gluinos. A gluino decays into t stop. Consequently, a top squark squark decays into b chi1^(+-) and chi1^(+-) --> W^(+-) chi1^0 (see Fig. 5b in the paper). This particular model assumes that m(stop) < m(gluino), m(chi1^0)=6 GeV, and m(chi1^(+-))=118 GeV.

The signal event yields are calculated for all simplified models. For each model, the values are given for the five signal regions and their combination. The simplified model is for direct production of squarks that decay into two steps into q q W Z W Z chi1^0 chi1^0 (see Fig. 6c in the paper).

The signal event yields are calculated for all simplified extra dimension model (see Fig. 8d in the paper). For each model, the values are given for the five signal regions and their combination.

The signal event yields are calculated for all simplified models. For each model, the values are given for the five signal regions and their combination. The simplified model is for direct production of gluinos that decay via sleptons into q q q q l l (l l) chi1^0 chi1^0 + neutrinos (see Fig. 6d in the paper).

The signal event yields are calculated for all simplified models. For each model, the values are given for the five signal regions and their combination. The simplified model is for direct production of gluinos that decay into q q q q W W chi1^0 chi1^0 (see Fig. 6a in the paper).

The signal event yields are calculated for all simplified models. For each model, the values are given for the five signal regions and their combination. The simplified model is for direct production of gluinos that decay into t tbar t tbar chi1^0 chi1^0 (see Fig. 5a in the paper). This particular model assumes that top quark is much heavier than gluino.

The signal event yields are calculated for all mSUGRA models (see Fig. 8a in the paper). For each model, the values are given for the five signal regions and their combination. The model assumes tan(beta)=30, A0=2m0, and mu>0.

The signal event yields are calculated for all simplified models. For each model, the values are given for the five signal regions and their combination. The simplified model is for direct production of gluinos. A gluino decays into t c chi1^0 (see Fig. 5c in the paper). This particular model assumes that m(chi1^0) = m(stop)-20 GeV.

The signal event yields are calculated for all GMSB models (see Fig. 8c in the paper). For each model, the values are given for the five signal regions and their combination. The model assumes mmess=250 TeV, m5=3, mu>0, and Cgrav=1.

The signal event yields are calculated for all simplified models. For each model, the values are given for the five signal regions and their combination. The simplified model is for direct production of bottom squarks. A bottom squark decays into t chi1^(+-) and chi1^(+-) --> W^(+-) chi1^0 (see Fig. 7a in the paper). This particular model assumes that m(chi1^0)=60 GeV.

The signal event yields are calculated for all simplified models. For each model, the values are given for the five signal regions and their combination. The simplified model is for direct production of gluinos and top squarks. Top squarks undergo R-parity violating decays into bs and gluinos decay into t stop (see Fig. 5d in the paper).

The signal event yields are calculated for all simplified models. For each model, the values are given for the five signal regions and their combination. The simplified model is for direct production of bottom squarks. A bottom squark decays into t chi1^(+-) and chi1^(+-) --> W^(+-) chi1^0 (see Fig. 7b in the paper). This particular model assumes that m(chi1^0)=2(chi1^0).

The signal event yields are calculated for all mSUGRA/CMSSM models with bRPV (see Fig. 8b in the paper). For each model, the values are given for the five signal regions and their combination. The model assumes tan(beta)=30, A0=2m0, mu>0, and bRPV.

The signal event yields are calculated for all simplified models. For each model, the values are given for the five signal regions and their combination. The simplified model is for direct production of squarks. Squarks decay into q q l l (l l) chi1^0 chi1^0 + neutrinos (see Fig. 6e in the paper).

The signal event yields are calculated for all simplified models. For each model, the values are given for the five signal regions and their combination. The simplified model is for direct pair-production of gluinos that decay via a two-step process into q q q q W Z W Z chi1^0 chi1^0 (see Fig. 6b in the paper).

The signal event yields are calculated for all simplified models. For each model, the values are given for the five signal regions and their combination. The simplified model is for direct pair-production of gluinos. A gluino decays into t stop. Consequently, a top squark squark decays into b chi1^(+-) and chi1^(+-) --> W^(+-) chi1^0 (see Fig. 5b in the paper). This particular model assumes that m(stop) < m(gluino), m(chi1^0)=6 GeV, and m(chi1^(+-))=118 GeV.

Experimental uncertainties on the signal event yields are calculated for all simplified models. For each model, the values are given for the five signal regions and their combination. The simplified model is for direct production of squarks that decay into two steps into q q W Z W Z chi1^0 chi1^0 (see Fig. 6c in the paper).

Experimental uncertainties on the signal event yields are calculated for all simplified extra dimension model (see Fig. 8d in the paper). For each model, the values are given for the five signal regions and their combination.

Experimental uncertainties on the signal event yields are calculated for all simplified models. For each model, the values are given for the five signal regions and their combination. The simplified model is for direct production of gluinos that decay via sleptons into q q q q l l (l l) chi1^0 chi1^0 + neutrinos (see Fig. 6d in the paper).

Experimental uncertainties on the signal event yields are calculated for all simplified models. For each model, the values are given for the five signal regions and their combination. The simplified model is for direct production of gluinos that decay into q q q q W W chi1^0 chi1^0 (see Fig. 6a in the paper).

Experimental uncertainties on the signal event yields are calculated for all simplified models. For each model, the values are given for the five signal regions and their combination. The simplified model is for direct production of gluinos that decay into t tbar t tbar chi1^0 chi1^0 (see Fig. 5a in the paper). This particular model assumes that top quark is much heavier than gluino.

Experimental uncertainties on the signal event yields are calculated for all mSUGRA models (see Fig. 8a in the paper). For each model, the values are given for the five signal regions and their combination. The model assumes tan(beta)=30, A0=2m0, and mu>0.

Experimental uncertainties on the signal event yields are calculated for all simplified models. For each model, the values are given for the five signal regions and their combination. The simplified model is for direct production of gluinos. A gluino decays into t c chi1^0 (see Fig. 5c in the paper). This particular model assumes that m(chi1^0) = m(stop) - 20 GeV.

Experimental uncertainties on the signal event yields are calculated for all GMSB models (see Fig. 8c in the paper). For each model, the values are given for the five signal regions and their combination. The model assumes mmess=250 TeV, m5=3, mu>0, and Cgrav=1.

Experimental uncertainties on the signal event yields are calculated for all simplified models. For each model, the values are given for the five signal regions and their combination. The simplified model is for direct production of bottom squarks. A bottom squark decays into t chi1^(+-) and chi1^(+-) --> W^(+-) chi1^0 (see Fig. 7a in the paper). This particular model assumes that m(chi1^0)=60 GeV.

Experimental uncertainties on the signal event yields are calculated for all simplified models. For each model, the values are given for the five signal regions and their combination. The simplified model is for direct production of gluinos and top squarks. Top squarks undergo R-parity violating decays into bs and gluinos decay into t stop (see Fig. 5d in the paper).

Experimental uncertainties on the signal event yields are calculated for all simplified models. For each model, the values are given for the five signal regions and their combination. The simplified model is for direct production of bottom squarks. A bottom squark decays into t chi1^(+-) and chi1^(+-) --> W^(+-) chi1^0 (see Fig. 7b in the paper). This particular model assumes that m(chi1^0)=2(chi1^0).

Experimental uncertainties on the signal event yields are calculated for all mSUGRA/CMSSM models with bRPV (see Fig. 8b in the paper). For each model, the values are given for the five signal regions and their combination. The model assumes tan(beta)=30, A0=2m0, mu>0, and bRPV.

Experimental uncertainties on the signal event yields are calculated for all simplified models. For each model, the values are given for the five signal regions and their combination. The simplified model is for direct production of squarks. Squarks decay into q q l l (l l) chi1^0 chi1^0 + neutrinos (see Fig. 6e in the paper).

Experimental uncertainties on the signal event yields are calculated for all simplified models. For each model, the values are given for the five signal regions and their combination. The simplified model is for direct pair-production of gluinos that decay via a two-step process into q q q q W Z W Z chi1^0 chi1^0 (see Fig. 6b in the paper).

Experimental uncertainties on the signal event yields are calculated for all simplified models. For each model, the values are given for the five signal regions and their combination. The simplified model is for direct pair-production of gluinos. A gluino decays into t stop. Consequently, a top squark squark decays into b chi1^(+-) and chi1^(+-) --> W^(+-) chi1^0 (see Fig. 5b in the paper). This particular model assumes that m(stop) < m(gluino), m(chi1^0)=6 GeV, and m(chi1^(+-))=118 GeV.

Statistical uncertainties on the signal event yields are calculated for all simplified models. For each model, the values are given for the five signal regions and their combination. The simplified model is for direct production of squarks that decay into two steps into q q W Z W Z chi1^0 chi1^0 (see Fig. 6c in the paper).

Statistical uncertainties on the signal event yields are calculated for all simplified extra dimension model (see Fig. 8d in the paper). For each model, the values are given for the five signal regions and their combination.

Statistical uncertainties on the signal event yields are calculated for all simplified models. For each model, the values are given for the five signal regions and their combination. The simplified model is for direct production of gluinos that decay via sleptons into q q q q l l (l l) chi1^0 chi1^0 + neutrinos (see Fig. 6d in the paper).

Statistical uncertainties on the signal event yields are calculated for all simplified models. For each model, the values are given for the five signal regions and their combination. The simplified model is for direct production of gluinos that decay into q q q q W W chi1^0 chi1^0 (see Fig. 6a in the paper).

Statistical uncertainties on the signal event yields are calculated for all simplified models. For each model, the values are given for the five signal regions and their combination. The simplified model is for direct production of gluinos that decay into t tbar t tbar chi1^0 chi1^0 (see Fig. 5a in the paper). This particular model assumes that top quark is much heavier than gluino.

Statistical uncertainties on the signal event yields are calculated for all mSUGRA models (see Fig. 8a in the paper). For each model, the values are given for the five signal regions and their combination. The model assumes tan(beta)=30, A0=2m0, and mu>0.

Statistical uncertainties on the signal event yields are calculated for all simplified models. For each model, the values are given for the five signal regions and their combination. The simplified model is for direct production of gluinos. A gluino decays into t c chi1^0 (see Fig. 5c in the paper). This particular model assumes that m(chi1^0) = m(stop) - 20 GeV.

Statistical uncertainties on the signal event yields are calculated for all GMSB models (see Fig. 8c in the paper). For each model, the values are given for the five signal regions and their combination. The model assumes mmess=250 TeV, m5=3, mu>0, and Cgrav=1.

Statistical uncertainties on the signal event yields are calculated for all simplified models. For each model, the values are given for the five signal regions and their combination. The simplified model is for direct production of bottom squarks. A bottom squark decays into t chi1^(+-) and chi1^(+-) --> W^(+-) chi1^0 (see Fig. 7a in the paper). This particular model assumes that m(chi1^0)=60 GeV.

Statistical uncertainties on the signal event yields are calculated for all simplified models. For each model, the values are given for the five signal regions and their combination. The simplified model is for direct production of gluinos and top squarks. Top squarks undergo R-parity violating decays into bs and gluinos decay into t stop (see Fig. 5d in the paper).

Statistical uncertainties on the signal event yields are calculated for all simplified models. For each model, the values are given for the five signal regions and their combination. The simplified model is for direct production of bottom squarks. A bottom squark decays into t chi1^(+-) and chi1^(+-) --> W^(+-) chi1^0 (see Fig. 7b in the paper). This particular model assumes that m(chi1^0)=2(chi1^0).

Statistical uncertainties on the signal event yields are calculated for all mSUGRA/CMSSM models with bRPV (see Fig. 8b in the paper). For each model, the values are given for the five signal regions and their combination. The model assumes tan(beta)=30, A0=2m0, mu>0, and bRPV.

Statistical uncertainties on the signal event yields are calculated for all simplified models. For each model, the values are given for the five signal regions and their combination. The simplified model is for direct production of squarks. Squarks decay into q q l l (l l) chi1^0 chi1^0 + neutrinos (see Fig. 6e in the paper).

Statistical uncertainties on the signal event yields are calculated for all simplified models. For each model, the values are given for the five signal regions and their combination. The simplified model is for direct pair-production of gluinos that decay via a two-step process into q q q q W Z W Z chi1^0 chi1^0 (see Fig. 6b in the paper).

Statistical uncertainties on the signal event yields are calculated for all simplified models. For each model, the values are given for the five signal regions and their combination. The simplified model is for direct pair-production of gluinos. A gluino decays into t stop. Consequently, a top squark squark decays into b chi1^(+-) and chi1^(+-) --> W ^(+-) chi1^0 (see Fig. 5b in the paper). This particular model assumes that m(stop) < m(gluino), m(chi1^0)=6 GeV, and m(chi1^(+-))=118 GeV.

The confidence levels are calculated for all simplified models. For each model, the observed and expected values are given. The simplified model is for direct production of gluinos that decay into t tbar t tbar chi1^0 chi1^0 (see Fig. 5a in the paper). This particular model assumes that top quark is much heavier than gluino.

The confidence levels are calculated for all simplified models. For each model, the observed and expected values are given. The simplified model is for direct production of squarks that decay into two steps into q q W Z W Z chi1^0 chi1^0 (see Fig. 6c in the paper).

The confidence levels are calculated for all simplified models. For each model, the values are given for the five signal regions and their combination. The simplified model is for direct pair-production of gluinos that decay via a two-step process into q q q q W Z W Z chi1^0 chi1^0 (see Fig. 6b in the paper).

The confidence levels are calculated for all simplified models. For each model, the expected and observed values are given. The simplified model is for direct production of gluinos that decay via sleptons into q q q q l l (l l) chi1^0 chi1^0 + neutrinos (see Fig. 6d in the paper).

The confidence levels are calculated for all simplified models. For each model, the expected and observed values are given. The simplified model is for direct pair-production of gluinos. A gluino decays into t stop. Consequently, a top squark squark decays into b chi1^(+-) and chi1^(+-) --> W^(+-) chi1^0 (see Fig. 5b in the paper). This particular model assumes that m(stop) < m(gluion), m(chi1^0)=6 GeV, and m(chi1^(+-))=118 GeV.

The confidence levels are calculated for all simplified models. For each model, the expected and observed values are given. The simplified model is for direct production of gluinos. A gluino decays into t c chi1^0 (see Fig. 5c in the paper). This particular model assumes that m(chi1^0) = m(stop) - 20 GeV.

The confidence levels are calculated for all simplified models. For each model, the expected and observed values are given. The simplified model is for direct production of bottom squarks. A bottom squark decays into t chi1^(+-) and chi1^(+-) --> W^(+-) chi1^0 (see Fig. 7b in the paper). This particular model assumes that m(chi1^0)=2(chi1^0).

The confidence levels are calculated for all simplified models. For each model, the expected and observed values are given. The simplified model is for direct production of bottom squarks. A bottom squark decays into t chi1^(+-) and chi1^(+-) --> W^(+-) chi1^0 (see Fig. 7a in the paper). This particular model assumes that m(chi1^0)=60 GeV.

The confidence levels are calculated for all simplified models. For each model, the expected and observed values are given. The simplified model is for direct production of squarks. Squarks decay into q q l l (l l) chi1^0 chi1^0 + neutrinos (see Fig. 6e in the paper).

The confidence levels are calculated for all GMSB models (see Fig. 8c in the paper). For each model, the expected and observed values are given. The model assumes mmess=250 TeV, m5=3, mu>0, and Cgrav=1.

The confidence levels are calculated for all simplified models. For each model, the expected and observed values are given. The simplified model is for direct production of gluinos and top squarks. Top squarks undergo R-parity violating decays into bs and gluinos decay into t stop (see Fig. 5d in the paper).

The confidence levels are calculated for all mSUGRA/CMSSM models with bRPV (see Fig. 8b in the paper). For each model, the expected and observed values are given. The model assumes tan(beta)=30, A0=2m0, mu>0, and bRPV.

The confidence levels are calculated for all simplified extra dimension model (see Fig. 8d in the paper). For each model, the expected and observed values are given.

The confidence levels are calculated for all simplified models. For each model, the expected and observed values are given. The simplified model is for direct production of gluinos that decay into q q q q W W chi1^0 chi1^0 (see Fig. 6a in the paper).

The confidence levels are calculated for all mSUGRA models (see Fig. 8a in the paper). For each model, the expected and observed values are given. The model assumes tan(beta)=30, A0=2m0, and mu>0.

The effective mass distribution in VR2Z.

The ETmiss distribution in SR0noZa.

The effective mass distribution in SR0noZa.

The ETmiss distribution in SR1noZa.

The effective mass distribution in SR1noZa.

The ETmiss distribution in SR2noZa.

The effective mass distribution in SR2noZa.

The ETmiss distribution in SR0noZb.

The effective mass distribution in SR0noZb.

The ETmiss distribution in SR1noZb.

The effective mass distribution in SR1noZb.

The ETmiss distribution in SR2noZb.

The effective mass distribution in SR2noZb.

The ETmiss distribution in SR0Z.

The effective mass distribution in SR0Z.

The ETmiss distribution in SR1Z.

The effective mass distribution in SR1Z.

The ETmiss distribution in SR2Z.

The effective mass distribution in SR2Z.

Observed 95% CL exclusion contour for the RPV chargino NLSP model with lambda_121 != 0.

Expected 95% CL exclusion contour for the RPV chargino NLSP model with lambda_121 != 0.

Observed 95% CL exclusion contour for the RPV chargino NLSP model with lambda_122 != 0.

Expected 95% CL exclusion contour for the RPV chargino NLSP model with lambda_122 != 0.

Observed 95% CL exclusion contour for the RPV chargino NLSP model with lambda_133 != 0.

Expected 95% CL exclusion contour for the RPV chargino NLSP model with lambda_133 != 0.

Observed 95% CL exclusion contour for the RPV chargino NLSP model with lambda_233 != 0.

Expected 95% CL exclusion contour for the RPV chargino NLSP model with lambda_233 != 0.

Observed 95% CL exclusion contour for the RPV gluino NLSP model with lambda_121 != 0.

Expected 95% CL exclusion contour for the RPV gluino NLSP model with lambda_121 != 0.

Observed 95% CL exclusion contour for the RPV gluino NLSP model with lambda_122 != 0.

Expected 95% CL exclusion contour for the RPV gluino NLSP model with lambda_122 != 0.

Observed 95% CL exclusion contour for the RPV gluino NLSP model with lambda_133 != 0.

Expected 95% CL exclusion contour for the RPV gluino NLSP model with lambda_133 != 0.

Observed 95% CL exclusion contour for the RPV gluino NLSP model with lambda_233 != 0.

Expected 95% CL exclusion contour for the RPV gluino NLSP model with lambda_233 != 0.

Observed 95% CL exclusion contour for the RPV Lslepton NLSP model with lambda_121 != 0.

Expected 95% CL exclusion contour for the RPV Lslepton NLSP model with lambda_121 != 0.

Observed 95% CL exclusion contour for the RPV Lslepton NLSP model with lambda_122 != 0.

Expected 95% CL exclusion contour for the RPV Lslepton NLSP model with lambda_122 != 0.

Observed 95% CL exclusion contour for the RPV Lslepton NLSP model with lambda_133 != 0.

Expected 95% CL exclusion contour for the RPV Lslepton NLSP model with lambda_133 != 0.

Observed 95% CL exclusion contour for the RPV Lslepton NLSP model with lambda_233 != 0.

Expected 95% CL exclusion contour for the RPV Lslepton NLSP model with lambda_233 != 0.

Observed 95% CL exclusion contour for the RPV Rslepton NLSP model with lambda_121 != 0.

Expected 95% CL exclusion contour for the RPV Rslepton NLSP model with lambda_121 != 0.

Observed 95% CL exclusion contour for the RPV Rslepton NLSP model with lambda_122 != 0.

Expected 95% CL exclusion contour for the RPV Rslepton NLSP model with lambda_122 != 0.

Observed 95% CL exclusion contour for the RPV Rslepton NLSP model with lambda_133 != 0.

Expected 95% CL exclusion contour for the RPV Rslepton NLSP model with lambda_133 != 0.

Observed 95% CL exclusion contour for the RPV Rslepton NLSP model with lambda_233 != 0.

Expected 95% CL exclusion contour for the RPV Rslepton NLSP model with lambda_233 != 0.

Observed 95% CL exclusion contour for the RPV sneutrino NLSP model with lambda_121 != 0.

Expected 95% CL exclusion contour for the RPV sneutrino NLSP model with lambda_121 != 0.

Observed 95% CL exclusion contour for the RPV sneutrino NLSP model with lambda_122 != 0.

Expected 95% CL exclusion contour for the RPV sneutrino NLSP model with lambda_122 != 0.

Observed 95% CL exclusion contour for the RPV sneutrino NLSP model with lambda_133 != 0.

Expected 95% CL exclusion contour for the RPV sneutrino NLSP model with lambda_133 != 0.

Observed 95% CL exclusion contour for the RPV sneutrino NLSP model with lambda_233 != 0.

Expected 95% CL exclusion contour for the RPV sneutrino NLSP model with lambda_233 != 0.

Observed 95% CL exclusion contour for the R-slepton RPC model.

Expected 95% CL exclusion contour for the R-slepton RPC model.

Observed and expected 95% CL cross-section upper limits for the Stau RPC model, together with the theoretically predicted cross-section.

Observed and expected 95% CL cross-section upper limits for the Z RPC model, together with the theoretically predicted cross-section.

Observed 95% CL exclusion contour for the GGM tan beta = 1.5 model.

Expected 95% CL exclusion contour for the GGM tan beta = 1.5 model.

Observed 95% CL exclusion contour for the GGM tan beta = 30 model.

Expected 95% CL exclusion contour for the GGM tan beta = 30 model.

Cut flows for a representative selection of SUSY signal points in the Z-veto signal regions. In each case, m2 and m1 refer to the axes of the plots in Sec. XI, where m2 is the larger of the two masses. The number of events expected for a luminosity of 20.3 fb-1 is quoted at each step of the selection. The preselection requires four baseline leptons, at least two of which are light leptons; the signal lepton selection is made at the ``Lepton Multiplicity' stage. ``Event Cleaning' refers to the selection criteria applied to remove non-collision backgrounds and detector noise.

Cut flows for a representative selection of SUSY signal points in the Z-rich signal regions. In each case, m2 and m1 refer to the axes of the plots in Sec. XI, where m2 is the larger of the two masses (or the value of mu in the case of GGM models). The number of events expected for a luminosity of 20.3 fb-1 is quoted at each step of the selection. The preselection requires four baseline leptons, at least two of which are light leptons; the signal lepton selection is made at the ``Lepton Multiplicity' stage. ``Event Cleaning' refers to the selection criteria applied to remove non-collision backgrounds and detector noise.

Cut flows by lepton channel for a representative selection of SUSY signal points in the SR0noZa signal region. In each case, m2 and m1 refer to the axes of the plots in Sec. XI, where m2 is the larger of the two masses. The number of events expected for a luminosity of 20.3 fb-1 is quoted at each step of the selection. The preselection requires four baseline leptons, at least two of which are light leptons; the signal lepton selection is made at the ``Lepton Multiplicity' stage. ``Event Cleaning' refers to the selection criteria applied to remove non-collision backgrounds and detector noise. The RPC R-slepton model is used, with (m2,m1) = (450,300) GeV.

Cut flows by lepton channel for a representative selection of SUSY signal points in the SR1noZb signal region. In each case, m2 and m1 refer to the axes of the plots in Sec. XI, where m2 is the larger of the two masses. The number of events expected for a luminosity of 20.3 fb-1 is quoted at each step of the selection. The preselection requires four baseline leptons, at least two of which are light leptons; the signal lepton selection is made at the ``Lepton Multiplicity' stage. ``Event Cleaning' refers to the selection criteria applied to remove non-collision backgrounds and detector noise. The RPV gluino NLSP model is used, with lambda_133 != 0 and (m2,m1) = (800,400) GeV.

Cut flows by lepton channel for a representative selection of SUSY signal points in the SR0Z signal region. In each case, m2 and m1 refer to the axes of the plots in Sec. XI, where m2 is the value of mu. The number of events expected for a luminosity of 20.3 fb-1 is quoted at each step of the selection. The preselection requires four baseline leptons, at least two of which are light leptons; the signal lepton selection is made at the ``Lepton Multiplicity' stage. ``Event Cleaning' refers to the selection criteria applied to remove non-collision backgrounds and detector noise. The GGM tan beta = 30 model is used, with (m2,m1) = (200,1000) GeV.

The observed and expected $E_T^{miss}$ distribution in the dimuon SR-Z. The negigible estimated contribution from Z+jets is omitted in these distributions. The last bin contains the overflow.

The observed and expected dielectron invariant mass distribution in SR-loose. The last bin contains the overflow.

The observed and expected dimuon invariant mass distribution in SR-loose. The last bin contains the overflow.

The observed and expected dielectron invariant mass distribution in the two-jet $b$-veto SR. The last bin contains the overflow.

The observed and expected dimuon invariant mass distribution in the two-jet $b$-veto SR. The last bin contains the overflow.

The observed and expected dielectron invariant mass distribution in the four jet b-veto SR. The last bin contains the overflow.

The observed and expected dimuon invariant mass distribution in the four-jet $b$-veto SR. The last bin contains the overflow.

The observed and expected dielectron invariant mass distribution in the two-jet $b$-tag SR. The last bin contains the overflow.

The observed and expected dimuon invariant mass distribution in two-jet $b$-tag SR. The last bin contains the overflow.

The observed and expected dielectron invariant mass distribution in the four-jet $b$-tag SR. The last bin contains the overflow.

The observed and expected dimuon invariant mass distribution in the four-jet $b$-tag SR. The last bin contains the overflow.

Expected 95% exclusion contour for the GGM model with $\tan(\beta)=1.5$ in SR-Z.

Observed 95% exclusion contour for the GGM model with $\tan(\beta)=1.5$ in SR-Z.

Expected 95% exclusion contour for the GGM model with $\tan(\beta)=30$ in SR-Z.

Observed 95% exclusion contour for the GGM model with $\tan(\beta)=30$ in SR-Z.

Expected 95% exclusion contour for the two-step first- and second-generation squark simplified model with sleptons in the two-jet $b$-veto SR.

Observed 95% exclusion contour for the two-step first- and second-generation squark simplified model with sleptons in the two-jet $b$-veto SR.

Expected 95% exclusion contour for the two-step gluino simplified model with sleptons in the four-jet $b$-veto SR.

Observed 95% exclusion contour for the two-step gluino simplified model with sleptons in the four-jet $b$-veto SR.

Number of generated events in the two-step gluino simplified model with sleptons.

Production cross-section in the two-step gluino simplified model with sleptons.

Number of generated events in the two-step first- and second-generation squark simplified model with sleptons.

Production cross-section in the two-step first- and second-generation squark simplified model with sleptons.

Number of generated events in the GGM model with $\tan(\beta)=1.5$.

Production cross-section in the GGM model with $\tan(\beta)=1.5$.

Number of generated events in the GGM model with $\tan(\beta)=30$.

Production cross-section in the GGM model with $\tan(\beta)=30$.

Total experimental uncertainty [%] for the two-step gluino simplified model with sleptons.

Total experimental uncertainty [%] for the GGM model with $\tan(\beta)=1.5$.

Total experimental uncertainty for the GGM model with $\tan(\beta)=30$.

Signal acceptance for the GGM model with $\tan(\beta)=1.5$ in the combined electron and muon SR-Z.

Signal acceptance for the GGM model with $\tan(\beta)=30$ in the combined electron and muon SR-Z.

Signal efficiency for the GGM model with $\tan(\beta)=1.5$ in the dielectron SR-Z.

Signal efficiency for the GGM model with $\tan(\beta)=30$ in the dielectron SR-Z.

Signal efficiency for the GGM model with $\tan(\beta)=1.5$ in the dimuon SR-Z.

Signal efficiency for the GGM model with $\tan(\beta)=30$ in the dimuon SR-Z.

Signal efficiency for the GGM model with $\tan(\beta)=1.5$ in the electron and muon combined SR-Z.

Signal efficiency for the GGM model with $\tan(\beta)=30$ in the the electron and muon combined SR-Z.

Signal acceptance for the two-step first- and second-generation squarks simplified model with sleptons in the two-jet $b$-veto SR.

Signal acceptance for the two-step gluino simplified model with sleptons in the four-jet $b$-veto SR.

Signal efficiency for the two-step first- and second-generation squarks simplified model with sleptons in the two-jet $b$-veto SR.

Signal efficiency for the two-step gluino simplified model with sleptons in the four-jet $b$-veto SR.

Upper limits on the signal cross-section at 95% CL for the GGM model with $\tan(\beta)=1.5$.

Observed CL$_{\text{S}}$ for the GGM model with $\tan(\beta)=1.5$.

Expected CL$_{\text{S}}$ for the GGM model with $\tan(\beta)=1.5$.

Upper limits on the signal cross-section at 95% CL for the GGM model with $\tan(\beta)=30$.

Observed CL$_{\text{S}}$ for the GGM model with $\tan(\beta)=30$.

Expected CL$_{\text{S}}$ for the GGM model with $\tan(\beta)=30$.

Upper limits on the signal stength at 95% CL for the two-step first- and second-generation squark simplified model with sleptons. The excluded signal strength is defined as the ratio of the observed excluded production cross section to the expected production cross section calculated at NLO+NLL.

Upper limits on the signal stength at 95% CL for the two-step gluino simplified model with sleptons. The excluded signal strength is defined as the ratio of the observed excluded production cross section to the expected production cross section calculated at NLO+NLL.

Observed CL$_{\text{S}}$ for the two-step first- and second-generation squark simplified model with sleptons.

Expected CL$_{\text{S}}$ for the two-step first- and second-generation squark simplified model with sleptons.

Observed CL$_{\text{S}}$ for the two-step gluino simplified model with sleptons.

Expected CL$_{\text{S}}$ for the two-step gluino simplified model with sleptons.

Number of data events (20.3 fb$^{-1}$), number of predicted events from the fit, statistical uncertainty on the fit, systematic uncertainty on the choice of control region, and on the choice of fit function versus inclusive $H_{\textrm{T}}^{\textrm{min}}$ lower bin edge for inclusive jet multiplicity $N_{\textrm{Jet}} \geq 6$. The total uncertainty is obtained by adding the three uncertainties linearly.

Number of data events (20.3 fb$^{-1}$), number of predicted events from the fit, statistical uncertainty on the fit, systematic uncertainty on the choice of control region, and on the choice of fit function versus inclusive $H_{\textrm{T}}^{\textrm{min}}$ lower bin edge for inclusive jet multiplicity $N_{\textrm{Jet}} \geq 7$. The total uncertainty is obtained by adding the three uncertainties linearly.

Number of data events (20.3 fb$^{-1}$), number of predicted events from the fit, statistical uncertainty on the fit, systematic uncertainty on the choice of control region, and on the choice of fit function versus inclusive $H_{\textrm{T}}^{\textrm{min}}$ lower bin edge for inclusive jet multiplicity $N_{\textrm{Jet}} \geq 8$. The total uncertainty is obtained by adding the three uncertainties linearly.

Efficiency for the calorimetric MET>80 GeV trigger as a function of the metastable R-hadron mass. The R-hadron decays to g/qq plus neutralino of mass = m(gluino) - 100 GeV with a lifetime of 1 ns.

Efficiency for the calorimetric MET>80 GeV trigger as a function of the metastable R-hadron mass. The R-hadron decays to g/qq plus neutralino of mass 100 GeV with a lifetime of 1 ns.

Efficiency for the calorimetric MET>80 GeV trigger as a function of the stable chargino mass.

Total selection efficiency as a function of the stable R-hadron mass.

Total selection efficiency as a function of the metastable R-hadron mass. The R-hadron decays to g/qq plus neutralino of mass 100 GeV with a lifetime of 10 ns.

Total selection efficiency as a function of the metastable R-hadron mass. The R-hadron decays to g/qq plus neutralino of mass = m(gluino) - 100 GeV with a lifetime of 10 ns.

Total selection efficiency as a function of the metastable R-hadron mass. The R-hadron decays to g/qq plus neutralino of mass 100 GeV with a lifetime of 1 ns.

Total selection efficiency as a function of the stable chargino mass.

Ionisation distribution of all the CR2 tracks, and those not matched to a reconstructed muon. The two distributions are normalised to their total number of entries.

Distribution of the mass of selected candidates, derived from the specific ionisation loss, for an example of gluino R-hadron signal, for searches for stable particles. The signal distributions are stacked on the expected background, and a narrower binning is used for them to allow the signal shape to be seen more clearly. The number of signal events is that expected according to the theoretical cross sections.

Distribution of the mass of selected candidates, derived from the specific ionisation loss, for one example of chargino signal, for searches for stable particles. The signal distributions are stacked on the expected background, and a narrower binning is used for them to allow the signal shape to be seen more clearly. The number of signal events is that expected according to the theoretical cross sections.

Distribution of the mass of selected candidates, derived from the specific ionisation loss, for background and data, for searches for stable particles. The expected background is shown with its total uncertainty (sum in quadrature of statistical, normalisation and systematic errors).

Distribution of the mass of selected candidates, derived from the specific ionisation loss, for an example of gluino R-hadron signal, for searches for metastable particles. The signal distributions are stacked on the expected background, and a narrower binning is used for them to allow the signal shape to be seen more clearly. The number of signal events is that expected according to the theoretical cross sections.

Distribution of the mass of selected candidates, derived from the specific ionisation loss, for an example of chargino signal, for searches for metastable particles. The signal distributions are stacked on the expected background, and a narrower binning is used for them to allow the signal shape to be seen more clearly. The number of signal events is that expected according to the theoretical cross sections.

Distribution of the mass of selected candidates, derived from the specific ionisation loss, for background and data. The expected background is shown with its total uncertainty (sum in quadrature of statistical, normalisation and systematic errors).

Theoretical values for the cross section of gluino pairs production with their uncertainty.

Expected upper limits on the production cross section as a function of mass for metastable gluino R-hadrons, with lifetime tau =10 ns, decaying into g/qq plus a light neutralino of mass 100 GeV, in the background-only case, with its 1 sigma band.

Observed 95 PCT upper limits on the production cross section as a function of mass for metastable gluino R-hadrons, with lifetime tau =10 ns, decaying into g/qq plus a light neutralino of mass 100 GeV.

Expected upper limits on the production cross section as a function of mass for metastable gluino R-hadrons, with lifetime tau =10 ns, decaying into g/qq plus a heavy neutralino of mass(gluino) - 100 GeV, in the background-only case, with its 1 sigma band.

Observed 95 PCT upper limits on the production cross section as a function of mass for metastable gluino R-hadrons, with lifetime tau =10 ns, decaying into g/qq plus a heavy neutralino of mass(gluino) - 100 GeV.

The expected excluded range of lifetimes as a function of gluino mass for gluino R-hadrons decaying into g/qq plus a light neutralino of mass 100 GeV, with respect to the nominal theoretical cross section.

The expected excluded range of lifetimes as a function of gluino mass for gluino R-hadrons decaying into g/qq plus a light neutralino of mass 100 GeV, with respect to the nominal theoretical cross section, plus 1 experimental sigma.

The expected excluded range of lifetimes as a function of gluino mass for gluino R-hadrons decaying into g/qq plus a light neutralino of mass 100 GeV, with respect to the nominal theoretical cross section, minus 1 experimental sigma.

The observed excluded range of lifetimes as a function of gluino mass for gluino R-hadrons decaying into g/qq plus a light neutralino of mass 100 GeV, with respect to the nominal theoretical cross section.

The observed excluded range of lifetimes as a function of gluino mass for gluino R-hadrons decaying into g/qq plus a light neutralino of mass 100 GeV, with respect to the nominal theoretical cross section minus its uncertainty.

The observed excluded range of lifetimes as a function of gluino mass for gluino R-hadrons decaying into g/qq plus a light neutralino of mass 100 GeV, with respect to the nominal theoretical cross section plus its uncertainty.

The expected excluded range of lifetimes as a function of gluino mass for gluino R-hadrons decaying into g/qq plus a heavy neutralino of mass = m(gluino) - 100 GeV, with respect to the nominal theoretical cross section.

The expected excluded range of lifetimes as a function of gluino mass for gluino R-hadrons decaying into g/qq plus a heavy neutralino of mass = m(gluino) - 100 GeV, with respect to the nominal theoretical cross section, plus 1 experimental sigma.

The expected excluded range of lifetimes as a function of gluino mass for gluino R-hadrons decaying into g/qq plus a heavy neutralino of mass = m(gluino) - 100 GeV, with respect to the nominal theoretical cross section, minus 1 experimental sigma.

The observed excluded range of lifetimes as a function of gluino mass for gluino R-hadrons decaying into g/qq plus a heavy neutralino of mass = m(gluino) - 100 GeV, with respect to the nominal theoretical cross section.

The observed excluded range of lifetimes as a function of gluino mass for gluino R-hadrons decaying into g/qq plus a heavy neutralino of mass = m(gluino) - 100 GeV, with respect to the nominal theoretical cross section minus its uncertainty.

The observed excluded range of lifetimes as a function of gluino mass for gluino R-hadrons decaying into g/qq plus a heavy neutralino of mass = m(gluino) - 100 GeV, with respect to the nominal theoretical cross section plus its uncertainty.

Expected upper limits on the production cross section as a function of mass for metastable gluino R-hadrons, with lifetime tau =10 ns, decaying into tt plus a light neutralino of mass 100 GeV, in the background-only case, with its 1 sigma band.

Observed 95 PCT upper limits on the production cross section as a function of mass for metastable gluino R-hadrons, with lifetime tau =10 ns, decaying into tt plus a light neutralino of mass 100 GeV.

Expected upper limits on the production cross section as a function of mass for metastable gluino R-hadrons, with lifetime tau =10 ns, decaying into tt plus a heavy neutralino of mass(gluino) - 100 GeV, in the background-only case, with its 1 sigma band.

Observed 95 PCT upper limits on the production cross section as a function of mass for metastable gluino R-hadrons, with lifetime tau =10 ns, decaying into tt plus a heavy neutralino of mass(gluino) - 100 GeV.

The expected excluded range of lifetimes as a function of gluino mass for gluino R-hadrons decaying into tt plus a light neutralino of mass 100 GeV, with respect to the nominal theoretical cross section.

The expected excluded range of lifetimes as a function of gluino mass for gluino R-hadrons decaying into tt plus a light neutralino of mass 100 GeV, with respect to the nominal theoretical cross section, plus 1 experimental sigma.

The expected excluded range of lifetimes as a function of gluino mass for gluino R-hadrons decaying into tt plus a light neutralino of mass 100 GeV, with respect to the nominal theoretical cross section, minus 1 experimental sigma.

The observed excluded range of lifetimes as a function of gluino mass for gluino R-hadrons decaying into tt plus a light neutralino of mass 100 GeV, with respect to the nominal theoretical cross section.

The observed excluded range of lifetimes as a function of gluino mass for gluino R-hadrons decaying into tt plus a light neutralino of mass 100 GeV, with respect to the nominal theoretical cross section minus its uncertainty.

The observed excluded range of lifetimes as a function of gluino mass for gluino R-hadrons decaying into tt plus a light neutralino of mass 100 GeV, with respect to the nominal theoretical cross section plus its uncertainty.

The expected excluded range of lifetimes as a function of gluino mass for gluino R-hadrons decaying into tt plus a heavy neutralino of mass = m(gluino) - 100 GeV, with respect to the nominal theoretical cross section.

The expected excluded range of lifetimes as a function of gluino mass for gluino R-hadrons decaying into tt plus a heavy neutralino of mass = m(gluino) - 100 GeV, with respect to the nominal theoretical cross section, plus 1 experimental sigma.

The expected excluded range of lifetimes as a function of gluino mass for gluino R-hadrons decaying into tt plus a heavy neutralino of mass = m(gluino) - 100 GeV, with respect to the nominal theoretical cross section, minus 1 experimental sigma.

The observed excluded range of lifetimes as a function of gluino mass for gluino R-hadrons decaying into tt plus a heavy neutralino of mass = m(gluino) - 100 GeV, with respect to the nominal theoretical cross section.

The observed excluded range of lifetimes as a function of gluino mass for gluino R-hadrons decaying into tt plus a heavy neutralino of mass = m(gluino) - 100 GeV, with respect to the nominal theoretical cross section minus its uncertainty.

The observed excluded range of lifetimes as a function of gluino mass for gluino R-hadrons decaying into tt plus a heavy neutralino of mass = m(gluino) - 100 GeV, with respect to the nominal theoretical cross section plus its uncertainty.

Theoretical values for the production cross section of charginos or chargino/neutralino pairs, with their uncertainty.

Expected upper limits on the production cross section as a function of mass for metastable charginos, with lifetime tau =1.0 ns, decaying into neutralino + pion, in the background-only case, with its 1 sigma band.

Observed 95 PCT upper limits on the production cross section as a function of mass for metastable charginos, with lifetime tau =1.0 ns, decaying into neutralino + pion.

The expected excluded range of lifetimes as a function of chargino mass for charginos decaying into neutralino plus pion, with respect to the nominal theoretical cross section.

The expected excluded range of lifetimes as a function of chargino mass for charginos decaying into neutralino plus pion, with respect to the nominal theoretical cross section, plus 1 experimental sigma.

The expected excluded range of lifetimes as a function of chargino mass for charginos decaying into neutralino plus pion, with respect to the nominal theoretical cross section, minus 1 experimental sigma.

The observed excluded range of lifetimes as a function of chargino mass for charginos decaying into neutralino plus pion, with respect to the nominal theoretical cross section.

The observed excluded range of lifetimes as a function of chargino mass for charginos decaying into neutralino plus pion, with respect to the nominal theoretical cross section minus its uncertainty.

The observed excluded range of lifetimes as a function of gluino mass for chargino mass for charginos decaying into neutralino plus pion, with respect to the nominal theoretical cross section plus its uncertainty.

dEdx ionization for data, 1 TeV gluino R-hadrons stable and decaying in 100 GeV neutralinos with a 10 ns lifetime and for charginos of 350 GeV. Tracks that fulfil all the requirements up to including the High-m_T (see Tab.1 in the paper) are considered at this stage and normalised to an integrated luminosity of 18.4 fb^-1.

Expected upper limits on the production cross section as a function of mass for stable gluino R-hadrons, in case of background only, with its 1 sigma band.

Observed 95 PCT upper limits on the production cross section as a function of mass for stable gluino R-hadrons.

Theoretical values for the cross section of squark pairs production with their uncertainty.

Expected upper limits on the production cross section as a function of mass for stable sbottom R-hadrons, in case of background only, with its 1 sigma band.

Observed 95 PCT upper limits on the production cross section as a function of mass for stable sbottom $R$-hadrons. Cross section IN PB.

Expected upper limits on the production cross section as a function of mass for stop R-hadrons, in case of background only, with its 1 sigma band.

Observed 95 PCT upper limits on the production cross section as a function of mass for stop R-hadrons.

Expected upper limits on the production cross section as a function of mass for stable charginos, in case of background only, with its 1 sigma band.

Observed 95 PCT upper limits on the production cross section as a function of mass for stable charginos.

Expected upper limits on the production cross section as a function of mass for metastable gluino R-hadrons, with lifetime tau=1.0 ns, decaying to g/qq plus a light neutralino of mass 100 GeV.

Observed 95 PCT upper limits on the production cross section as a function of mass for metastable gluino R-hadrons, with lifetime tau=1.0 ns, decaying to g/qq plus a light neutralino of mass 100 GeV.

Expected upper limits on the production cross section as a function of mass for metastable gluino R-hadrons, with lifetime tau=1.0 ns, decaying to g/qq plus a heavy neutralino of mass = m(gluino) - 100 GeV.

Observed 95 PCT upper limits on the production cross section as a function of mass for metastable gluino R-hadrons, with lifetime tau=1.0 ns, decaying to g/qq plus a heavy neutralino of mass = m(gluino) - 100 GeV.

Expected upper limits on the production cross section as a function of mass for metastable gluino R-hadrons, with lifetime tau=1.0 ns, decaying to tt plus a light neutralino of mass 100 GeV.

Observed 95 PCT upper limits on the production cross section as a function of mass for metastable gluino R-hadrons, with lifetime tau=1.0 ns, decaying to tt plus a light neutralino of mass 100 GeV.

Expected upper limits on the production cross section as a function of mass for metastable gluino R-hadrons, with lifetime tau=1.0 ns, decaying to tt plus a heavy neutralino of mass = m(gluino) - 100 GeV.

Observed 95 PCT upper limits on the production cross section as a function of mass for metastable gluino R-hadrons, with lifetime tau=1.0 ns, decaying to tt plus a heavy neutralino of mass = m(gluino) - 100 GeV.

Expected upper limits on the production cross section as a function of mass for metastable charginos, with lifetime tau =15 ns, decaying to neutralino and pion, in case of background only, with its 1 sigma band.

Observed 95 PCT upper limits on the production cross section as a function of mass for metastable charginos, with lifetime tau =15 ns, decaying to neutralino and pion, in case of background only, with its 1 sigma band.

$E_{\mathrm{T}}^{\mathrm{miss}} / \sqrt{H_{\mathrm{T}}}$ distribution in signal region 10j50 2b. Two benchmark signal models are overlaid on the plot for comparison. Labelled `pMSSM' and `2-step', they show signal distributions from the example SUSY models (as described in the paper): a pMSSM slice model with ($m \tilde{g}$, $m \tilde{\chi_{1}^{\pm}}$) = (1300, 200) GeV and a cascade decay model with ($m \tilde{g}$, $m \tilde{\chi_{1}^{0}}$) = (1300, 200) GeV.

$E_{\mathrm{T}}^{\mathrm{miss}} / \sqrt{H_{\mathrm{T}}}$ distribution in signal region 8j80 0b. Two benchmark signal models are overlaid on the plot for comparison. Labelled `pMSSM' and `2-step', they show signal distributions from the example SUSY models (as described in the paper): a pMSSM slice model with ($m \tilde{g}$, $m \tilde{\chi_{1}^{\pm}}$) = (1300, 200) GeV and a cascade decay model with ($m \tilde{g}$, $m \tilde{\chi_{1}^{0}}$) = (1300, 200) GeV.

$E_{\mathrm{T}}^{\mathrm{miss}} / \sqrt{H_{\mathrm{T}}}$ distribution in signal region 8j80 2b. Two benchmark signal models are overlaid on the plot for comparison. Labelled `pMSSM' and `2-step', they show signal distributions from the example SUSY models (as described in the paper): a pMSSM slice model with ($m \tilde{g}$, $m \tilde{\chi_{1}^{\pm}}$) = (1300, 200) GeV and a cascade decay model with ($m \tilde{g}$, $m \tilde{\chi_{1}^{0}}$) = (1300, 200) GeV.

Observed 95% CL limit for the pMSSM grid.

Observed 95% CL limit for the pMSSM grid when the signal cross section is increased by one standard deviation.

Observed 95% CL limit for the pMSSM grid when the signal cross section is decreased by one standard deviation.

Expected 95% CL limit for the pMSSM grid.

+1 sigma excursion of the expected 95% CL limit for the pMSSM grid.

-1 sigma excursion of the expected 95% CL limit for the pMSSM grid.

Observed 95% CL limit for the 2Step grid.

Observed 95% CL limit for the 2Step grid when the signal cross section is increased by one standard deviation.

Observed 95% CL limit for the 2Step grid when the signal cross section is decreased by one standard deviation.

Expected 95% CL limit for the 2Step grid.

+1 sigma excursion of the expected 95% CL limit for the 2Step grid.

-1 sigma excursion of the expected 95% CL limit for the 2Step grid.

$E_{\mathrm{T}}^{\mathrm{miss}} / \sqrt{H_{\mathrm{T}}}$ distribution in signal region 8j50 0b. Two benchmark signal models are overlaid on the plot for comparison. Labelled `pMSSM' and `2-step', they show signal distributions from the example SUSY models (as described in the paper): a pMSSM slice model with ($m \tilde{g}$, $m \tilde{\chi_{1}^{\pm}}$) = (1300, 200) GeV and a cascade decay model with ($m \tilde{g}$, $m \tilde{\chi_{1}^{0}}$) = (1300, 200) GeV.

$E_{\mathrm{T}}^{\mathrm{miss}} / \sqrt{H_{\mathrm{T}}}$ distribution in signal region 8j50 1b. Two benchmark signal models are overlaid on the plot for comparison. Labelled `pMSSM' and `2-step', they show signal distributions from the example SUSY models (as described in the paper): a pMSSM slice model with ($m \tilde{g}$, $m \tilde{\chi_{1}^{\pm}}$) = (1300, 200) GeV and a cascade decay model with ($m \tilde{g}$, $m \tilde{\chi_{1}^{0}}$) = (1300, 200) GeV.

$E_{\mathrm{T}}^{\mathrm{miss}} / \sqrt{H_{\mathrm{T}}}$ distribution in signal region 8j50 2b. Two benchmark signal models are overlaid on the plot for comparison. Labelled `pMSSM' and `2-step', they show signal distributions from the example SUSY models (as described in the paper): a pMSSM slice model with ($m \tilde{g}$, $m \tilde{\chi_{1}^{\pm}}$) = (1300, 200) GeV and a cascade decay model with ($m \tilde{g}$, $m \tilde{\chi_{1}^{0}}$) = (1300, 200) GeV.

$E_{\mathrm{T}}^{\mathrm{miss}} / \sqrt{H_{\mathrm{T}}}$ distribution in signal region 9j50 0b. Two benchmark signal models are overlaid on the plot for comparison. Labelled `pMSSM' and `2-step', they show signal distributions from the example SUSY models (as described in the paper): a pMSSM slice model with ($m \tilde{g}$, $m \tilde{\chi_{1}^{\pm}}$) = (1300, 200) GeV and a cascade decay model with ($m \tilde{g}$, $m \tilde{\chi_{1}^{0}}$) = (1300, 200) GeV.

$E_{\mathrm{T}}^{\mathrm{miss}} / \sqrt{H_{\mathrm{T}}}$ distribution in signal region 9j50 1b. Two benchmark signal models are overlaid on the plot for comparison. Labelled `pMSSM' and `2-step', they show signal distributions from the example SUSY models (as described in the paper): a pMSSM slice model with ($m \tilde{g}$, $m \tilde{\chi_{1}^{\pm}}$) = (1300, 200) GeV and a cascade decay model with ($m \tilde{g}$, $m \tilde{\chi_{1}^{0}}$) = (1300, 200) GeV.

$E_{\mathrm{T}}^{\mathrm{miss}} / \sqrt{H_{\mathrm{T}}}$ distribution in signal region 9j50 2b. Two benchmark signal models are overlaid on the plot for comparison. Labelled `pMSSM' and `2-step', they show signal distributions from the example SUSY models (as described in the paper): a pMSSM slice model with ($m \tilde{g}$, $m \tilde{\chi_{1}^{\pm}}$) = (1300, 200) GeV and a cascade decay model with ($m \tilde{g}$, $m \tilde{\chi_{1}^{0}}$) = (1300, 200) GeV.

$E_{\mathrm{T}}^{\mathrm{miss}} / \sqrt{H_{\mathrm{T}}}$ distribution in signal region 10j50 0b. Two benchmark signal models are overlaid on the plot for comparison. Labelled `pMSSM' and `2-step', they show signal distributions from the example SUSY models (as described in the paper): a pMSSM slice model with ($m \tilde{g}$, $m \tilde{\chi_{1}^{\pm}}$) = (1300, 200) GeV and a cascade decay model with ($m \tilde{g}$, $m \tilde{\chi_{1}^{0}}$) = (1300, 200) GeV.

$E_{\mathrm{T}}^{\mathrm{miss}} / \sqrt{H_{\mathrm{T}}}$ distribution in signal region 10j50 1b. Two benchmark signal models are overlaid on the plot for comparison. Labelled `pMSSM' and `2-step', they show signal distributions from the example SUSY models (as described in the paper): a pMSSM slice model with ($m \tilde{g}$, $m \tilde{\chi_{1}^{\pm}}$) = (1300, 200) GeV and a cascade decay model with ($m \tilde{g}$, $m \tilde{\chi_{1}^{0}}$) = (1300, 200) GeV.

$E_{\mathrm{T}}^{\mathrm{miss}} / \sqrt{H_{\mathrm{T}}}$ distribution in signal region 10j50 2b. Two benchmark signal models are overlaid on the plot for comparison. Labelled `pMSSM' and `2-step', they show signal distributions from the example SUSY models (as described in the paper): a pMSSM slice model with ($m \tilde{g}$, $m \tilde{\chi_{1}^{\pm}}$) = (1300, 200) GeV and a cascade decay model with ($m \tilde{g}$, $m \tilde{\chi_{1}^{0}}$) = (1300, 200) GeV.

$E_{\mathrm{T}}^{\mathrm{miss}} / \sqrt{H_{\mathrm{T}}}$ distribution in signal region 7j80 0b. Two benchmark signal models are overlaid on the plot for comparison. Labelled `pMSSM' and `2-step', they show signal distributions from the example SUSY models (as described in the paper): a pMSSM slice model with ($m \tilde{g}$, $m \tilde{\chi_{1}^{\pm}}$) = (1300, 200) GeV and a cascade decay model with ($m \tilde{g}$, $m \tilde{\chi_{1}^{0}}$) = (1300, 200) GeV.

$E_{\mathrm{T}}^{\mathrm{miss}} / \sqrt{H_{\mathrm{T}}}$ distribution in signal region 7j80 1b. Two benchmark signal models are overlaid on the plot for comparison. Labelled `pMSSM' and `2-step', they show signal distributions from the example SUSY models (as described in the paper): a pMSSM slice model with ($m \tilde{g}$, $m \tilde{\chi_{1}^{\pm}}$) = (1300, 200) GeV and a cascade decay model with ($m \tilde{g}$, $m \tilde{\chi_{1}^{0}}$) = (1300, 200) GeV.

$E_{\mathrm{T}}^{\mathrm{miss}} / \sqrt{H_{\mathrm{T}}}$ distribution in signal region 7j80 2b. Two benchmark signal models are overlaid on the plot for comparison. Labelled `pMSSM' and `2-step', they show signal distributions from the example SUSY models (as described in the paper): a pMSSM slice model with ($m \tilde{g}$, $m \tilde{\chi_{1}^{\pm}}$) = (1300, 200) GeV and a cascade decay model with ($m \tilde{g}$, $m \tilde{\chi_{1}^{0}}$) = (1300, 200) GeV.

$E_{\mathrm{T}}^{\mathrm{miss}} / \sqrt{H_{\mathrm{T}}}$ distribution in signal region 8j80 0b. Two benchmark signal models are overlaid on the plot for comparison. Labelled `pMSSM' and `2-step', they show signal distributions from the example SUSY models (as described in the paper): a pMSSM slice model with ($m \tilde{g}$, $m \tilde{\chi_{1}^{\pm}}$) = (1300, 200) GeV and a cascade decay model with ($m \tilde{g}$, $m \tilde{\chi_{1}^{0}}$) = (1300, 200) GeV.

$E_{\mathrm{T}}^{\mathrm{miss}} / \sqrt{H_{\mathrm{T}}}$ distribution in signal region 8j80 1b. Two benchmark signal models are overlaid on the plot for comparison. Labelled `pMSSM' and `2-step', they show signal distributions from the example SUSY models (as described in the paper): a pMSSM slice model with ($m \tilde{g}$, $m \tilde{\chi_{1}^{\pm}}$) = (1300, 200) GeV and a cascade decay model with ($m \tilde{g}$, $m \tilde{\chi_{1}^{0}}$) = (1300, 200) GeV.

$E_{\mathrm{T}}^{\mathrm{miss}} / \sqrt{H_{\mathrm{T}}}$ distribution in signal region 8j80 2b. Two benchmark signal models are overlaid on the plot for comparison. Labelled `pMSSM' and `2-step', they show signal distributions from the example SUSY models (as described in the paper): a pMSSM slice model with ($m \tilde{g}$, $m \tilde{\chi_{1}^{\pm}}$) = (1300, 200) GeV and a cascade decay model with ($m \tilde{g}$, $m \tilde{\chi_{1}^{0}}$) = (1300, 200) GeV.

Degree of multijet closure for signal and vaidation regions with at no b-jet requirement. The solid lines are the pre-fit predicted numbers of events and the points are the observed numbers. The blue hatched band shows only the statistical (MC and data) uncertainty on the background estimate. The bins labelled in bold are signal regions, while the others are validation regions. The template closure uncertainty for each SR bin is given by the maximal deviation of data from prediction in any non-SR bin to its left on this plot (although those for 80 GeV regions are independent of deviations in 50 GeV regions).

Degree of multijet closure for signal and vaidation regions with at least 1 b-jet. The solid lines are the pre-fit predicted numbers of events and the points are the observed numbers. The blue hatched band shows only the statistical (MC and data) uncertainty on the background estimate. The bins labelled in bold are signal regions, while the others are validation regions. The template closure uncertainty for each SR bin is given by the maximal deviation of data from prediction in any non-SR bin to its left on this plot (although those for 80 GeV regions are independent of deviations in 50 GeV regions).

Degree of multijet closure for signal and vaidation regions with at least 2 b-jets. The solid lines are the pre-fit predicted numbers of events and the points are the observed numbers. The blue hatched band shows only the statistical (MC and data) uncertainty on the background estimate. The bins labelled in bold are signal regions, while the others are validation regions. The template closure uncertainty for each SR bin is given by the maximal deviation of data from prediction in any non-SR bin to its left on this plot (although those for 80 GeV regions are independent of deviations in 50 GeV regions).

Summary of all 15 signal regions (post-fit).

95% CLs observed upper limit on model cross-section for 2-step signal points for the best-expected signal region.

Data (bold boxes) and background estimates (colour fill) for m_beta vs. m_betagamma for the gluino R-hadron search (1000 GeV). The blue thin-line boxes illustrate the expected signal (on top of background) for the given R-hadron mass hypothesis. The black dashed vertical/horizontal lines at 500 GeV show the mass selection (signal region in the top-right). Two events pass this selection.

Expected (dashed black line) and observed (solid red line) 95% CL upper limits on the cross section as a function of mass for the production of long-lived gluino R-hadrons. The theory prediction along with its +-1sigma uncertainty is show as a black line and a blue band, respectively. The observed 8 TeV Run-1 limit and theory prediction [arXiv:1411.6795] are shown in dash-dotted and dotted lines, respectively.

Expected (dashed black line) and observed (solid red line) 95% CL upper limits on the cross section as a function of mass for the production of bottom-squark R-hadrons. The theory prediction along with its +-1sigma uncertainty is show as a black line and a blue band, respectively. The observed 8 TeV Run-1 limit and theory prediction [arXiv:1411.6795] are shown in dash-dotted and dotted lines, respectively.

Expected (dashed black line) and observed (solid red line) 95% CL upper limits on the cross section as a function of mass for the production of top-squark R-hadrons. The theory prediction along with its +-1sigma uncertainty is show as a black line and a blue band, respectively. The observed 8 TeV Run-1 limit and theory prediction [arXiv:1411.6795] are shown in dash-dotted and dotted lines, respectively.

Final selection requirements as a function of the simulated R-hadron mass.

Summary of all studied systematic uncertainties. Ranges indicate a dependency on the R-hadron mass hypothesis (from low to high masses).

Expected signal yield (Nsig) and efficiency (eff.), estimated background (Nbkg) and observed number of events in data (Nobs) for the full mass range after the final selection using 3.2/fb of data. The stated uncertainties include both the statistical and systematic contribution.

Distribution of the truth-level beta for gluino R-hadrons in exemplary signal MC samples and muons in a Zmumu MC sample. All distributions have been normalised to one. The last bin contains the overflow of the histograms. The distributions illustrate the good discriminating power of the variables.

Distribution of the truth-level betagamma for gluino R-hadrons in exemplary signal MC samples and muons in a Zmumu MC sample. All distributions have been normalised to one. The last bin contains the overflow of the histograms. The distributions illustrate the good discriminating power of the variables.

Expected (dashed black line) and observed (solid red line) 95% confidence level upper limits on the cross section as a function of mass for the production of long-lived gluino R-hadrons. The theory prediction along with its +-1sigma uncertainty is show as a black line and a blue band, respectively. For meta-stable gluinos with a lifetime of 50 ns. (mass exclusion: about 1660 GeV expected, 1520 GeV observed).

Expected (dashed black line) and observed (solid red line) 95% confidence level upper limits on the cross section as a function of mass for the production of long-lived gluino R-hadrons. The theory prediction along with its +-1sigma uncertainty is show as a black line and a blue band, respectively. For meta-stable gluinos with a lifetime of 30 ns. (mass exclusion: about 1660 GeV expected, 1520 GeV observed).

Expected (dashed black line) and observed (solid red line) 95% confidence level upper limits on the cross section as a function of mass for the production of long-lived gluino R-hadrons. The theory prediction along with its +-1sigma uncertainty is show as a black line and a blue band, respectively. For meta-stable gluinos with a lifetime of 10 ns. (mass exclusion: about 1660 GeV expected, 1520 GeV observed).

Object-quality selection cut-flow with observed data and exemplary expected events (scaled to 3.2/fb for MC) in the gluino R-hadron search.

Object-quality selection cut-flow with observed data and exemplary expected events (scaled to 3.2/fb for MC) in the squark R-hadron search.

Expected signal yield (Nsig) and efficiency (eff.), estimated background (Nbkg) and observed number of events in data (Nobs) for the full mass range in the meta-stable gluino R-hadron search using 3.2/fb of data. The stated uncertainties include both the statistical and systematic contribution.

Comparison of data with estimated backgrounds in the $\tilde{E}_\text{T}^\text{miss}$ distribution with the TZCR1 event selection except for the requirement on $\tilde{E}_\text{T}^\text{miss}$. The variables $\tilde{E}_\text{T}^\text{miss}$ and $\tilde{m}_\text{T}$ are constructed in the same way as $E_\text{T}^\text{miss}$ and $m_\text{T}$ but treating the leading photon transverse momentum as invisible. The predicted backgrounds are scaled with normalization factors. The uncertainty band includes statistical and all experimental systematic uncertainties. The last bin includes overflow.

Comparison of data with estimated backgrounds in the $\tilde{m}_\text{T}$ distribution with the TZCR1 event selection except for the requirement on $\tilde{m}_\text{T}$. The variables $\tilde{E}_\text{T}^\text{miss}$ and $\tilde{m}_\text{T}$ are constructed in the same way as $E_\text{T}^\text{miss}$ and $m_\text{T}$ but treating the leading photon transverse momentum as invisible. The predicted backgrounds are scaled with normalization factors. The uncertainty band includes statistical and all experimental systematic uncertainties. The last bin includes overflow.

Comparison of the observed data ($n_\text{obs}$) with the predicted background ($n_\text{exp}$) in the validation and signal regions. The background predictions are obtained using the background-only fit configuration. The bottom panel shows the significance of the difference between data and predicted background, where the significance is based on the total uncertainty ($\sigma_\text{tot}$).

Jet multiplicity distributions for events where exactly two signal leptons are selected. No correction factors are included in the background normalizations. The uncertainty band includes statistical and all experimental systematic uncertainties. The last bin includes overflow.

Jet multiplicity distributions for events where exactly one lepton plus one $\tau$ candidate are selected. No correction factors are included in the background normalizations. The uncertainty band includes statistical and all experimental systematic uncertainties. The last bin includes overflow.

The $E_\text{T}^\text{miss}$ distribution in SR1. In the plot, the full event selection in the corresponding signal region is applied, except for the requirement on $E_\text{T}^\text{miss}$. The predicted backgrounds are scaled with normalization factors. The uncertainty band includes statistical and all experimental systematic uncertainties. The last bin contains the overflow. Benchmark signal models are overlaid for comparison. The benchmark models are specified by the gluino and stop masses, given in TeV in the table.

The $m_\text{T}$ distribution in SR1. In the plot, the full event selection in the corresponding signal region is applied, except for the requirement on $m_\text{T}$. The predicted backgrounds are scaled with normalization factors. The uncertainty band includes statistical and all experimental systematic uncertainties. The last bin contains the overflow. Benchmark signal models are overlaid for comparison. The benchmark models are specified by the gluino and stop masses, given in TeV in the table.

Expected (black dashed) 95% excluded regions in the plane of $m_{\tilde{g}}$ versus $m_{\tilde{t}_1}$ for gluino-mediated stop production.

Observed (red solid) 95% excluded regions in the plane of $m_{\tilde{g}}$ versus $m_{\tilde{t}_1}$ for gluino-mediated stop production.

Expected (black dashed) 95% excluded regions in the plane of $m_{\tilde{t}_1}$ versus $m_{\tilde{\chi}_1^0}$ for direct stop production.

Observed (red solid) 95% excluded regions in the plane of $m_{\tilde{t}_1}$ versus $m_{\tilde{\chi}_1^0}$ for direct stop production.

The expected upper limits on $T$ quark pair production times the squared branching ratio for $T \rightarrow tZ$ as a function of the $T$ quark mass.

The observed upper limits on $T$ quark pair production times the squared branching ratio for $T \rightarrow tZ$ as a function of the $T$ quark mass.

The expected limits on $T$ quarks as a function of the branching ratios $B\left(T \rightarrow bW\right)$ and $B\left(T \rightarrow tH\right)$ for a $T$ quark with a mass of 800 GeV. The $T$ is assumed to decay in three possible ways: $T \to tZ$, $T \to tH$, and $T \to bW$.

The observed limits on $T$ quarks as a function of the branching ratios $B\left(T \rightarrow bW\right)$ and $B\left(T \rightarrow tH\right)$ for a $T$ quark with a mass of 800 GeV. The $T$ is assumed to decay in three possible ways: $T \to tZ$, $T \to tH$, and $T \to bW$.

The $m_\text{T}$ distribution in the WVR2-tail validation region which has the same preselection and jet $p_\text{T}$ requirements as SR2.

The $am_\text{T2}$ distribution in the WVR2-tail validation region which has the same preselection and jet $p_\text{T}$ requirements as SR2.

Large-radius jet mass ($R=1.2$), decomposed into the number of small-radius jet constituents. The lower panel shows the ratio of the total data to the total prediction (summed over all jet multiplicities). Events are required to have one lepton, four jets with $p_\text{T}>80,50,40,40$ GeV, at least one $b$-tagged jet, $E_\text{T}^\text{miss}>200$ GeV, and $m_\text{T}>30$ GeV.

Distribution of $m_\text{T2}^\tau$ in data for a selection enriched in $t\bar{t}$ events with one hadronically decaying $\tau$. Events that have no hadronic $\tau$ candidate (that passes the Loose identification criteria, as well as other requirements) are not shown in the plot.

Upper limits on the model cross-section in units of pb for the gluino-mediated stop models.

Upper limits on the model cross-section in units of pb for the models with direct stop pair production.

Illustration of the best expected signal region per signal grid point for the gluino-mediated stop models. This mapping is used for the final combined exclusion limits.

Illustration of the best expected signal region per signal grid point for models with direct stop pair production. This mapping is used for the final combined exclusion limits.

Expected $CL_s$ values for the gluino-mediated stop models.

Observed $CL_s$ values for the gluino-mediated stop models.

Expected $CL_s$ values for the direct stop pair production models.

Observed $CL_s$ values for the direct stop pair production models.

Expected limit using SR1 for models with direct stop pair production and an unpolarized stop (and bino LSP).

Expected limit using SR1 for models with direct stop pair production with $\tilde{t}_1=\tilde{t}_L$ (and bino LSP).

Expected limit using SR1 for models with direct stop pair production with $\tilde{t}_1\sim\tilde{t}_R$ (and bino LSP).

Observed limit using SR1 for models with direct stop pair production and an unpolarized stop (and bino LSP).

Observed limit using SR1 for models with direct stop pair production with $\tilde{t}_1=\tilde{t}_L$ (and bino LSP).

Observed limit using SR1 for models with direct stop pair production with $\tilde{t}_1\sim\tilde{t}_R$ (and bino LSP).

Expected limit using SR2 for models with direct stop pair production and an unpolarized stop (and bino LSP).

Expected limit using SR2 for models with direct stop pair production with $\tilde{t}_1=\tilde{t}_L$ (and bino LSP).

Expected limit using SR2 for models with direct stop pair production with $\tilde{t}_1\sim\tilde{t}_R$ (and bino LSP).

Observed limit using SR2 for models with direct stop pair production and an unpolarized stop (and bino LSP).

Observed limit using SR2 for models with direct stop pair production with $\tilde{t}_1=\tilde{t}_L$ (and bino LSP).

Observed limit using SR2 for models with direct stop pair production with $\tilde{t}_1\sim\tilde{t}_R$ (and bino LSP).

Expected limit using SR1+SR2 (best expected) for models with direct stop pair production and an unpolarized stop (and bino LSP).

Expected limit using SR1+SR2 (best expected) for models with direct stop pair production with $\tilde{t}_1=\tilde{t}_L$ (and bino LSP).

Expected limit using SR1+SR2 (best expected) for models with direct stop pair production with $\tilde{t}_1\sim\tilde{t}_R$ (and bino LSP).

Observed limit using SR1+SR2 (best expected) for models with direct stop pair production and an unpolarized stop (and bino LSP).

Observed limit using SR1+SR2 (best expected) for models with direct stop pair production with $\tilde{t}_1=\tilde{t}_L$ (and bino LSP).

Observed limit using SR1+SR2 (best expected) for models with direct stop pair production with $\tilde{t}_1\sim\tilde{t}_R$ (and bino LSP).

Acceptance for SR1 in the gluino-mediated stop models. The acceptance is defined as the fraction of signal events that pass the analysis selection performed on generator-level objects, therefore emulating an ideal detector with perfect particle identification and no measurement resolution effects.

Acceptance for SR1 in the direct stop pair production. The acceptance is defined as the fraction of signal events that pass the analysis selection performed on generator-level objects, therefore emulating an ideal detector with perfect particle identification and no measurement resolution effects.

Acceptance for SR2 in the gluino-mediated stop models. The acceptance is defined as the fraction of signal events that pass the analysis selection performed on generator-level objects, therefore emulating an ideal detector with perfect particle identification and no measurement resolution effects.

Acceptance for SR2 in the direct stop pair production. The acceptance is defined as the fraction of signal events that pass the analysis selection performed on generator-level objects, therefore emulating an ideal detector with perfect particle identification and no measurement resolution effects.

Acceptance for SR3 in the gluino-mediated stop models. The acceptance is defined as the fraction of signal events that pass the analysis selection performed on generator-level objects, therefore emulating an ideal detector with perfect particle identification and no measurement resolution effects.

Acceptance for SR3 in the direct stop pair production. The acceptance is defined as the fraction of signal events that pass the analysis selection performed on generator-level objects, therefore emulating an ideal detector with perfect particle identification and no measurement resolution effects.

Efficiency for SR1 in the gluino-mediated stop models. The efficiency is the ratio between the expected signal rate calculated with simulated data passing all the reconstruction level cuts applied to reconstructed objects, and the signal rate for an ideal detector (with perfect particle identification and no measurement resolution effects).

Efficiency for SR1 in the direct stop pair production. The efficiency is the ratio between the expected signal rate calculated with simulated data passing all the reconstruction level cuts applied to reconstructed objects, and the signal rate for an ideal detector (with perfect particle identification and no measurement resolution effects).

Efficiency for SR2 in the gluino-mediated stop models. The efficiency is the ratio between the expected signal rate calculated with simulated data passing all the reconstruction level cuts applied to reconstructed objects, and the signal rate for an ideal detector (with perfect particle identification and no measurement resolution effects).

Efficiency for SR2 in the direct stop pair production. The efficiency is the ratio between the expected signal rate calculated with simulated data passing all the reconstruction level cuts applied to reconstructed objects, and the signal rate for an ideal detector (with perfect particle identification and no measurement resolution effects).

Efficiency for SR3 in the gluino-mediated stop models. The efficiency is the ratio between the expected signal rate calculated with simulated data passing all the reconstruction level cuts applied to reconstructed objects, and the signal rate for an ideal detector (with perfect particle identification and no measurement resolution effects).

Efficiency for SR3 in the direct stop pair production. The efficiency is the ratio between the expected signal rate calculated with simulated data passing all the reconstruction level cuts applied to reconstructed objects, and the signal rate for an ideal detector (with perfect particle identification and no measurement resolution effects).