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 (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 (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 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.

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.

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-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.

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)$ = 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.

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.

The observed limit contour in the SQUARK-NEUTRALINO plane.

The expected limit contour in the SQUARK-NEUTRALINO plane.

Transverse mass after requiring one muon with pT>20 GeV, at least three jets with pT>60,25,25 GeV and dphi(jets,Etmiss)>0.2.

Effective mass after requiring one electron with pT>25 GeV, at least three jets with pT>60,25,25 GeV and dphi(jets,Etmiss)>0.2.

Effective mass after requiring one muon with pT>20 GeV, at least three jets with pT>60,25,25 GeV and dphi(jets,Etmiss)>0.2.

Missing transverse energy after requiring one electron with pT>25 GeV, at least four jets with pT>60,25,25,25 GeV and dphi(jets,Etmiss)>0.2.

Missing transverse energy after requiring one muon with pT>20 GeV, at least four jets with pT>60,25,25,25 GeV and dphi(jets,Etmiss)>0.2.

Transverse mass after requiring one electron with pT>25 GeV, at least four jets with pT>60,25,25,25 GeV and dphi(jets,Etmiss)>0.2.

Transverse mass after requiring one muon with pT>20 GeV, at least four jets with pT>60,25,25,25 GeV and dphi(jets,Etmiss)>0.2.

Effective mass after requiring one electron with pT>25 GeV, at least four jets with pT>60,25,25,25 GeV and dphi(jets,Etmiss)>0.2.

Effective mass after requiring one muon with pT>20 GeV, at least four jets with pT>60,25,25,25 GeV and dphi(jets,Etmiss)>0.2.

Effective mass in the electron plus three jets W+jets control region.

Effective mass in the muon plus three jets W+jets control region.

Effective mass in the electron plus three jets top control region.

Effective mass in the muon plus three jets top control region.

Number of b-tagged jets in the combined electron plus three jets W+jets and top control region.

Number of b-tagged jets in the combined muon plus three jets W+jets and top control region.

Effective mass in the electron plus four jets W+jets control region.

Effective mass in the muon plus four jets W+jets control region.

Effective mass in the electron plus four jets top control region.

Effective mass in the muon plus four jets top control region.

Effective mass in the 3-jet loose signal region (electron channel).

Effective mass in the 3-jet loose signal region (muon channel).

Effective mass in the 3-jet tight signal region (electron channel).

Effective mass in the 3-jet tight signal region (muon channel).

Effective mass in the 4-jet loose signal region (electron channel).

Effective mass in the 4-jet loose signal region (muon channel).

Effective mass in the 4-jet tight signal region (electron channel).

Effective mass in the 4-jet tight signal region (muon channel).

Expected 95% CL exclusion limits in the MSUGRA/CMSSM framework with the combined electron and muon channels.

Observed 95% CL exclusion limits in the MSUGRA/CMSSM framework with the combined electron and muon channels.

Expected 95% CL exclusion contour for the Gbb signal.

Observed 95% CL exclusion contour for the Gbb signal.

Expected 95% CL exclusion contour for the Gtt combination.

Observed 95% CL exclusion contour for the Gtt combination.

Acceptances for the Gbb model in SR-Gbb-A. Acceptance is evaluated at truth level, with only leptons from heavy bosons and taus considered, and no further quality or isolation criteria applied in their selection.

Acceptances for the Gbb model in SR-Gbb-B. Acceptance is evaluated at truth level, with only leptons from heavy bosons and taus considered, and no further quality or isolation criteria applied in their selection.

Acceptances for the Gbb model in SR-Gbb-C. Acceptance is evaluated at truth level, with only leptons from heavy bosons and taus considered, and no further quality or isolation criteria applied in their selection.

Acceptances for the Gtt model in SR-Gtt-0L-A. Acceptance is evaluated at truth level, with only leptons from heavy bosons and taus considered, and no further quality or isolation criteria applied in their selection.

Acceptances for the Gtt model in SR-Gtt-0L-B. Acceptance is evaluated at truth level, with only leptons from heavy bosons and taus considered, and no further quality or isolation criteria applied in their selection.

Acceptances for the Gtt model in SR-Gtt-0L-C. Acceptance is evaluated at truth level, with only leptons from heavy bosons and taus considered, and no further quality or isolation criteria applied in their selection.

Acceptances for the Gtt model in SR-Gtt-1L-A. Acceptance is evaluated at truth level, with only leptons from heavy bosons and taus considered, and no further quality or isolation criteria applied in their selection.

Acceptances for the Gtt model in SR-Gtt-1L-B. Acceptance is evaluated at truth level, with only leptons from heavy bosons and taus considered, and no further quality or isolation criteria applied in their selection.

Acceptance times efficiency for the Gbb model in SR-Gbb-A.

Acceptance times efficiency for the Gbb model in SR-Gbb-B.

Acceptance times efficiency for the Gbb model in SR-Gbb-C.

Acceptance times efficiency for the Gtt model in SR-Gtt-0L-A.

Acceptance times efficiency for the Gtt model in SR-Gtt-0L-B.

Acceptance times efficiency for the Gtt model in SR-Gtt-0L-C.

Acceptance times efficiency for the Gtt model in SR-Gtt-1L-A.

Acceptance times efficiency for the Gtt model in SR-Gtt-1L-B.

95% CL upper limit on the cross-section times branching ratio (in fb) for the Gbb model in SR-Gbb-A.

95% CL upper limit on the cross-section times branching ratio (in fb) for the Gbb model in SR-Gbb-B.

95% CL upper limit on the cross-section times branching ratio (in fb) for the Gbb model in SR-Gbb-C.

95% CL upper limit on the cross-section times branching ratio (in fb) for the Gtt model in SR-Gtt-0L-A.

95% CL upper limit on the cross-section times branching ratio (in fb) for the Gtt model in SR-Gtt-0L-B.

95% CL upper limit on the cross-section times branching ratio (in fb) for the Gtt model in SR-Gtt-0L-C.

95% CL upper limit on the cross-section times branching ratio (in fb) for the Gtt model in SR-Gtt-1L-A.

95% CL upper limit on the cross-section times branching ratio (in fb) for the Gtt model in SR-Gtt-1L-B.

Signal region yielding the best expected sensitivity for each point of the parameter space in the Gbb model.

Signal region yielding the best expected sensitivity for each point of the parameter space in the Gtt model for the 0-lepton channel.

Signal region yielding the best expected sensitivity for each point of the parameter space in the Gtt model for the 1-lepton channel.

Combination of two 0-lepton and 1-lepton signal regions yielding the best expected sensitivity for each point of the parameter space in the Gtt model.

The transverse momentum distribution of the 4th-leading lepton for events with at least 4 leptons and no Z-boson candidate.

Distribution of missing transverse momentum for events with at least 4 leptons and no Z-boson candidate.

Distribution of effective mass for events with at least 4 leptons and no Z-boson candidate.

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.

Distribution of HT after the MTtau requirement for the 1-tau ''Tight'' SR. The SM prediction includes the data-driven corrections discussed in the paper. MC events are normalized to data in the CRs corresponding to MTtau below 130 GeV. Also shown is the expected signal from typical mSUGRA, GMSB and bRPV samples. The last bin in the expected background distribution is an overflow bin.

Distribution of MTtau1 + MTtau2 in the 2tau channel after all analysis requirements but the final SR requirements on MTtau1 + MTtau2 and HT2j. To reduce the contributions from events with Z bosons decaying into tau leptons, the requirement MTtau1 + MTtau2 > 150 GeV is applied to all distributions. The SM prediction includes the data-driven corrections discussed in the paper. MC events are normalized to data in the CRs corresponding to HT2j below 550 GeV.

Distribution of HT2j in the 2tau channel after all analysis requirements but the final SR requirements on MTtau1 + MTtau2 and HT2j. To reduce the contributions from events with Z bosons decaying into tau leptons, the requirement MTtau1 + MTtau2 > 150 GeV is applied to all distributions. The SM prediction includes the data-driven corrections discussed in the paper. MC events are normalized to data in the CRs corresponding to HT2j below 550 GeV.

Distribution of the jet multiplicity in the 2tau channel after all analysis requirements but the final SR requirements on MTtau1 + MTtau2 and HT2j. To reduce the contributions from events with Z bosons decaying into tau leptons, the requirement MTtau1 + MTtau2 > 150 GeV is applied to all distributions. The SM prediction includes the data-driven corrections discussed in the paper. MC events are normalized to data in the CRs corresponding to HT2j below 550 GeV.

Distribution of the final kinematic variables in the tau+e channel after all analysis requirements but the final SR selections on Meff for the bRPV model. The SM prediction includes the data-driven corrections discussed in the paper. MC events are normalized to data in specific CRs. The last bin in the expected background distribution is an overflow bin. There are no data events in the overflow bin after all analysis requirements are applied.

Distribution of the final kinematic variables in the tau+e channel after all analysis requirements but the final SR selections on MEFF for the GMSB model. The SM prediction includes the data-driven corrections discussed in the paper. MC events are normalized to data in specific CRs. The last bin in the expected background distribution is an overflow bin. There are no data events in the overflow bin after all analysis requirements are applied.

Distribution of the final kinematic variables in the tau+e channel after all analysis requirements but the final SR selections on MET for the mSUGRA model. The SM prediction includes the data-driven corrections discussed in the paper. MC events are normalized to data in specific CRs. The last bin in the expected background distribution is an overflow bin. There are no data events in the overflow bin after all analysis requirements are applied.

Distribution of the final kinematic variables in the tau+e channel after all analysis requirements but the final SR selections on MET for the nGM model. The SM prediction includes the data-driven corrections discussed in the paper. MC events are normalized to data in specific CRs. The last bin in the expected background distribution is an overflow bin. There are no data events in the overflow bin after all analysis requirements are applied.

Distribution of the final kinematic variables in the tau+mu channel after all analysis requirements but the final SR selections on MEFF for the bRPV model. The SM prediction includes the data-driven corrections discussed in the paper. MC events are normalized to data in specific CRs. The last bin in the expected background distribution is an overflow bin. There are no data events in the overflow bin after all analysis requirements are applied.

Distribution of the final kinematic variables in the tau+mu channel after all analysis requirements but the final SR selections on MEFF for the GMSB model. The SM prediction includes the data-driven corrections discussed in the paper. MC events are normalized to data in specific CRs. The last bin in the expected background distribution is an overflow bin. There are no data events in the overflow bin after all analysis requirements are applied.

Distribution of the final kinematic variables in the tau+mu channel after all analysis requirements but the final SR selections on MET for the mSUGRA model. The SM prediction includes the data-driven corrections discussed in the paper. MC events are normalized to data in specific CRs. The last bin in the expected background distribution is an overflow bin. There are no data events in the overflow bin after all analysis requirements are applied.

Distribution of the final kinematic variables in the tau+mu channel after all analysis requirements but the final SR selections on MET for the nGM model. The SM prediction includes the data-driven corrections discussed in the paper. MC events are normalized to data in specific CRs. The last bin in the expected background distribution is an overflow bin. There are no data events in the overflow bin after all analysis requirements are applied.

Observed 95% CL lower limits on the minimal GMSB model parameters Lambda and tan(beta) using a combination of all channels. The result is obtained using 20.3 fb-1 of sqrt(s) = 8 TeV ATLAS data. Additional model parameters are M(mess) = 250 TeV, N5 = 3, mu>0 and Cgrav =1.

Expected 95% CL lower limits on the minimal GMSB model parameters Lambda and tan(beta) using a combination of all channels. The result is obtained using 20.3 fb-1 of sqrt(s) = 8 TeV ATLAS data. Additional model parameters are M(mess) = 250 TeV, N5 = 3, mu>0 and Cgrav =1.

Observed 95% CL lower limits on the mSUGRA/CMSSM model parameters m0 and m1/2 for the combination of the 1tau, tau+e and tau+mu channels. Additional model parameters are A0 = -2m0, tan(beta) = 30 and sign(mu) = +1.

Expected 95% CL lower limits on the mSUGRA/CMSSM model parameters m0 and m1/2 for the combination of the 1tau, tau+e and tau+mu channels. Additional model parameters are A0 = -2m0, tan(beta) = 30 and sign(mu) = +1.

Observed 95% CL lower limits on the simplified nGM model parameters m(stau) and m(gluino) for the combination of the 2tau, tau+e and tau+mu channels. Additional squark and slepton mass parameters are set to 2.5 TeV, M1 = M2 = 2.5 TeV, and all trilinear coupling terms are set to zero. Also, the parameter mu is fixed to mu = 400 GeV.

Expected 95% CL lower limits on the simplified nGM model parameters m(stau) and m(gluino) for the combination of the 2tau, tau+e and tau+mu channels. Additional squark and slepton mass parameters are set to 2.5 TeV, M1 = M2 = 2.5 TeV, and all trilinear coupling terms are set to zero. Also, the parameter mu is fixed to mu = 400 GeV.

Observed 95% CL lower limits on the bRPV model parameters m0 and m1/2 for the combination of all channels. Additional model parameters are A0 = -2m0 , tan(beta) = 30 and sign(mu) = +1.

Expected 95% CL lower limits on the bRPV model parameters m0 and m1/2 for the combination of all channels. Additional model parameters are A0 = -2m0 , tan(beta) = 30 and sign(mu) = +1.

Cross section predictions for the nGM grid. For each signal point 25000 MC events have been generated.

Observed upper cross section limits for the nGM grid. For each signal point 25000 MC events have been generated. The limit is derived for the combination of the 2tau and the tau+lepton channels.

Systematic uncertainty for the GMSB grid in the 1-tau analysis.

Acceptance for the GMSB grid in the 1-tau analysis.

Efficiency for the GMSB grid in the 1-tau analysis.

The product of acceptance and efficiency for the GMSB grid in the 1-tau analysis.

Expected CLs values for the GMSB grid in the 1tau analysis.

Observed CLs values for the GMSB grid in the 1tau analysis.

Systematic uncertainty for the mSUGRA grid in the 1-tau analysis.

Acceptance for the mSUGRA grid in the 1-tau analysis.

Efficiency for the mSUGRA grid in the 1-tau analysis.

Product of acceptance and efficiency for the mSUGRA grid in the 1-tau analysis.

Expected CLs values for the mSUGRA grid in the 1tau analysis.

Observed CLs values for the mSUGRA grid in the 1tau analysis.

Systematic uncertainty for the bRPV grid in the 1-tau analysis.

Acceptance for the bRPV grid in the 1-tau analysis.

Efficiency for the bRPV grid in the 1-tau analysis.

Product of acceptance and efficiency for the bRPV grid in the 1-tau analysis.

Expected CLs values for the bRPV grid in the 1tau analysis.

Observed CLs values for the bRPV grid in the 1tau analysis.

Systematic uncertainty for the GMSB grid in the 2tau analysis.

Acceptance for the GMSB grid in the 2tau analysis.

Efficiency for the GMSB grid in the 2tau analysis.

Product of acceptance and efficiency for the GMSB grid in the 2tau analysis.

Expected CLs values for the GMSB grid in the 2tau analysis.

Observed CLs values for the GMSB grid in the 2tau analysis.

Systematic uncertainty for the nGM grid in the 2tau analysis.

Acceptance for the nGM grid in the 2tau analysis.

Efficiency for the nGM grid in the 2tau analysis.

Product of acceptance and efficiency for the nGM grid in the 2tau analysis.

Expected CLs values for the nGM grid in the 2tau analysis.

Observed CLs values for the nGM grid in the 2tau analysis.

Systematic uncertainty for the bRPV grid in the 2tau analysis.

Acceptance for the bRPV grid in the 2tau analysis.

Efficiency for the bRPV grid in the 2tau analysis.

Product of acceptance and efficiency for the bRPV grid in the 2tau analysis.

Expected CLs values for the bRPV grid in the 2tau analysis.

Observed CLs values for the bRPV grid in the 2tau analysis.

Systematic Uncertainty for the GMSB grid in the tau+e analysis.

Acceptance for the GMSB grid in the tau+e analysis.

Efficiency for the GMSB grid in the tau+e analysis.

Product of acceptance and efficiency for the GMSB grid in the tau+e analysis.

Expected CLs values for the GMSB grid in the tau+e analysis.

Observed CLs values for the GMSB grid in the tau+e analysis.

Systematic Uncertainty for the nGM grid in the tau+e analysis.

Acceptance for the nGM grid in the tau+e analysis.

Efficiency for the nGM grid in the tau+e analysis.

Product of acceptance and efficiency for the nGM grid in the tau+e analysis.

Expected CLs values for the nGM grid in the tau+e analysis.

Observed CLs values for the nGM grid in the tau+e analysis.

Systematic Uncertainty for the bRPV grid in the tau+e analysis.

Acceptance for the bRPV grid in the tau+e analysis.

Efficiency for the bRPV grid in the tau+e analysis.

Product of acceptance and efficiency for the bRPV grid in the tau+e analysis.

Expected CLs values for the bRPV grid in the tau+e analysis.

Observed CLs values for the bRPV grid in the tau+e analysis.

Systematic Uncertainty for the mSUGRA grid in the tau+e analysis.

Acceptance for the mSUGRA grid in the tau+e analysis.

Efficiency for the mSUGRA grid in the tau+e analysis.

Product of acceptance and efficiency for the mSUGRA grid in the tau+e analysis.

Expected CLs values for the mSUGRA grid in the tau+e analysis.

Observed CLs values for the mSUGRA grid in the tau+e analysis.

Systematic Uncertainty for the GMSB grid in the tau+mu analysis.

Acceptance for the GMSB grid in the tau+mu analysis.

Efficiency for the GMSB grid in the tau+mu analysis.

Product of acceptance and efficiency for the GMSB grid in the tau+mu analysis.

Expected CLs values for the GMSB grid in the tau+mu analysis.

Observed CLs values for the GMSB grid in the tau+mu analysis.

Systematic Uncertainty for the nGM grid in the tau+mu analysis.

Acceptance for the nGM grid in the tau+mu analysis.

Efficiency for the nGM grid in the tau+mu analysis.

Product of acceptance and efficiency for the nGM grid in the tau+mu analysis.

Expected CLs values for the nGM grid in the tau+mu analysis.

Observed CLs values for the nGM grid in the tau+mu analysis.

Systematic Uncertainty for the bRPV grid in the tau+mu analysis.

Acceptance for the bRPV grid in the tau+mu analysis.

Efficiency for the bRPV grid in the tau+mu analysis.

Product of acceptance and efficiency for the bRPV grid in the tau+mu analysis.

Expected CLs values for the bRPV grid in the tau+mu analysis.

Observed CLs values for the bRPV grid in the tau+mu analysis.

Systematic Uncertainty for the mSUGRA grid in the tau+mu analysis.

Acceptance for the mSUGRA grid in the tau+mu analysis.

Efficiency for the mSUGRA grid in the tau+mu analysis.

Product of acceptance and efficiency for the mSUGRA grid in the tau+mu analysis.

Expected CLs values for the mSUGRA grid in the tau+mu analysis.

Observed CLs values for the mSUGRA grid in the tau+mu analysis.

Example cutflow for three benchmark signal points in the 1-tau channel. For the GMSB point 50000 events have been generated, 20000 for nGM and 20000 for mSUGRA respectively. These event numbers are then normalised to 21 fb^{-1} luminosity.

Example cutflow for three benchmark signal points in the 2tau analysis. Specific SRs are indicated at the respective cut. Event numbers are normalised to 21 fb^{-1} luminosity. For the GMSB point 50000 events have been generated, 20000 for nGM and 25000 for bRPV respectively.

Example cutflow for four benchmark signal points in the tau+e analysis. Event numbers are normalised to 21 fb^{-1} luminosity. For the GMSB point 50000 events have been generated, 20000 for nGM, 20000 for mSUGRA and 25000 for bRPV respectively.

Example cutflow for four benchmark signal points in the tau+mu analysis. Event numbers are normalised to 21 fb^{-1} luminosity. For the GMSB point 50000 events have been generated, 20000 for nGM, 20000 for mSUGRA and 25000 for bRPV respectively.

Missing transverse momentum distribution after SR3b selection, beside the $E_\mathrm{T}^\mathrm{miss}$ requirement. The results in the signal region correspond to the last inclusive bin. The systematic uncertainties include theory uncertainties for the backgrounds with prompt SS/3L and the full systematic uncertainties for data-driven backgrounds. For illustration the distribution for a benchmark SUSY scenario ($pp\to \tilde g\tilde g$, $\tilde g\to t\bar t\tilde\chi_1^0$, $m_{\tilde g}=1.2$ TeV, $m_{\tilde\chi_1^0}=0.7$ TeV) is also shown.

Observed exclusion limits on the $\tilde g$ and $\tilde\chi_1^0$ masses in the context of SUSY scenarios with simplified mass spectra featuring $\tilde g\tilde g$ pair production with exclusive $\tilde g\to qq(\tilde\ell\ell/\tilde\nu\nu)$ decays. All limits are computed at 95% CL.

Expected exclusion limits on the $\tilde g$ and $\tilde\chi_1^0$ masses in the context of SUSY scenarios with simplified mass spectra featuring $\tilde g\tilde g$ pair production with exclusive $\tilde g\to qq(\tilde\ell\ell/\tilde\nu\nu)$ decays. All limits are computed at 95% CL.

Upper limits on signal cross-sections as function of the $\tilde g$ and $\tilde\chi_1^0$ masses in the context of SUSY scenarios with simplified mass spectra featuring $\tilde g\tilde g$ pair production with exclusive $\tilde g\to qq(\tilde\ell\ell/\tilde\nu\nu)$ decays, obtained using the signal efficiency and acceptance specific to each model. All limits are computed at 95% CL.

Observed exclusion limits on the $\tilde g$ and $\tilde\chi_1^0$ masses in the context of SUSY scenarios with simplified mass spectra featuring $\tilde g\tilde g$ pair production with exclusive $\tilde g\to qqWZ\tilde\chi_1^0$ decays. All limits are computed at 95% CL.

Expected exclusion limits on the $\tilde g$ and $\tilde\chi_1^0$ masses in the context of SUSY scenarios with simplified mass spectra featuring $\tilde g\tilde g$ pair production with exclusive $\tilde g\to qqWZ\tilde\chi_1^0$ decays. All limits are computed at 95% CL.

Upper limits on signal cross-sections as function of the $\tilde g$ and $\tilde\chi_1^0$ masses in the context of SUSY scenarios with simplified mass spectra featuring $\tilde g\tilde g$ pair production with exclusive $\tilde g\to qqWZ\tilde\chi_1^0$ decays, obtained using the signal efficiency and acceptance specific to each model. All limits are computed at 95% CL.

Observed exclusion limits on the $\tilde b_1$ and $\tilde\chi_1^0$ masses in the context of SUSY scenarios with simplified mass spectra featuring $\tilde b_1\tilde b_1^*$ pair production with exclusive $\tilde b_1\to t\tilde\chi_1^-$ decays. All limits are computed at 95% CL.

Expected exclusion limits on the $\tilde b_1$ and $\tilde\chi_1^0$ masses in the context of SUSY scenarios with simplified mass spectra featuring $\tilde b_1\tilde b_1^*$ pair production with exclusive $\tilde b_1\to t\tilde\chi_1^-$ decays. All limits are computed at 95% CL.

Upper limits on signal cross-sections as function of the $\tilde b_1$ and $\tilde\chi_1^0$ masses in the context of SUSY scenarios with simplified mass spectra featuring $\tilde b_1\tilde b_1^*$ pair production with exclusive $\tilde b_1\to t\tilde\chi_1^-$ decays, obtained using the signal efficiency and acceptance specific to each model. All limits are computed at 95% CL.

Observed exclusion limits on the $\tilde g$ and $\tilde\chi_1^0$ masses in the context of SUSY scenarios with simplified mass spectra featuring $\tilde g\tilde g$ pair production with exclusive $\tilde g\to t\bar t\tilde\chi_1^0$ decays. All limits are computed at 95% CL.

Expected exclusion limits on the $\tilde g$ and $\tilde\chi_1^0$ masses in the context of SUSY scenarios with simplified mass spectra featuring $\tilde g\tilde g$ pair production with exclusive $\tilde g\to t\bar t\tilde\chi_1^0$ decays. All limits are computed at 95% CL.

Upper limits on signal cross-sections as function of the $\tilde g$ and $\tilde\chi_1^0$ masses in the context of SUSY scenarios with simplified mass spectra featuring $\tilde g\tilde g$ pair production with exclusive $\tilde g\to t\bar t\tilde\chi_1^0$ decays, obtained using the signal efficiency and acceptance specific to each model. All limits are computed at 95% CL.

SUSY scenario with $\tilde g\tilde g$ production and $\tilde g\to q\bar q(\tilde\ell\ell/\tilde\nu\nu)$ decay: signal acceptance (in %) in the signal region SR0b3j. The benchmark scenarios used to set exclusion limits are materialized by black dot markers. Acceptance and efficiency are defined as in appendix A of [JHEP 06 (2014) 124, arXiv: 1403.4853v1 [hep-ex]].

SUSY scenario with $\tilde g\tilde g$ production and $\tilde g\to q\bar q(\tilde\ell\ell/\tilde\nu\nu)$ decay: reconstruction efficiency (in %) in the signal region SR0b3j. The benchmark scenarios used to set exclusion limits are materialized by black dot markers. Acceptance and efficiency are defined as in appendix A of [JHEP 06 (2014) 124, arXiv: 1403.4853v1 [hep-ex]].

SUSY scenario with $\tilde g\tilde g$ production and $\tilde g\to q\bar qWZ\tilde\chi_1^0$ decay: signal acceptance (in %) in the signal region SR0b5j. The benchmark scenarios used to set exclusion limits are materialized by black dot markers. Acceptance and efficiency are defined as in appendix A of [JHEP 06 (2014) 124, arXiv: 1403.4853v1 [hep-ex]].

SUSY scenario with $\tilde g\tilde g$ production and $\tilde g\to q\bar qWZ\tilde\chi_1^0$ decay: reconstruction efficiency (in %) in the signal region SR0b5j. The benchmark scenarios used to set exclusion limits are materialized by black dot markers. Acceptance and efficiency are defined as in appendix A of [JHEP 06 (2014) 124, arXiv: 1403.4853v1 [hep-ex]].

SUSY scenario with $\tilde b_1\tilde b_1^*$ production and $\tilde b_1\to tW\tilde\chi_1^0$ decay: signal acceptance (in %) in the signal region SR1b. The benchmark scenarios used to set exclusion limits are materialized by black dot markers. Acceptance and efficiency are defined as in appendix A of [JHEP 06 (2014) 124, arXiv: 1403.4853v1 [hep-ex]].

SUSY scenario with $\tilde b_1\tilde b_1^*$ production and $\tilde b_1\to tW\tilde\chi_1^0$ decay: reconstruction efficiency (in %) in the signal region SR1b. The benchmark scenarios used to set exclusion limits are materialized by black dot markers. Acceptance and efficiency are defined as in appendix A of [JHEP 06 (2014) 124, arXiv: 1403.4853v1 [hep-ex]].

SUSY scenario with $\tilde g\tilde g$ production and $\tilde g\to t\bar t\tilde\chi_1^0$ decay: signal acceptance (in %) in the signal region SR3b. The benchmark scenarios used to set exclusion limits are materialized by black dot markers. Acceptance and efficiency are defined as in appendix A of [JHEP 06 (2014) 124, arXiv: 1403.4853v1 [hep-ex]].

SUSY scenario with $\tilde g\tilde g$ production and $\tilde g\to t\bar t\tilde\chi_1^0$ decay: reconstruction efficiency (in %) in the signal region SR3b. The benchmark scenarios used to set exclusion limits are materialized by black dot markers. Acceptance and efficiency are defined as in appendix A of [JHEP 06 (2014) 124, arXiv: 1403.4853v1 [hep-ex]].