The post-fit BDT score distribution for the direct stau channel, showing the scores for BDT4, before the selections on the BDT score is made. The black arrow depicts the BDT score selection for the SR-BDT. A few example SUSY scenarios targeted by each BDT are overlaid for illustration.

The post-fit kinematic distribution for the Intermediate stau channel, showing the $m_{\mathrm{T}2}$ distribution in SR-C1C1-LM, before the selection on $m_{\mathrm{T}2}$ is made. The black arrow depicts the selection for the signal region. An example SUSY scenario is overlaid for illustration.

The post-fit kinematic distribution for the Intermediate stau channel, showing the $m_{\mathrm{T}2}$ distribution in SR-C1N2OS-LM, before the selection on $m_{\mathrm{T}2}$ is made. The black arrow depicts the selection for the signal region. An example SUSY scenario is overlaid for illustration.

The post-fit kinematic distribution for the Intermediate stau channel, showing the $m_{\mathrm{T}2}$ distribution in SR-C1N2SS-LM, before the selection on $m_{\mathrm{T}2}$ is made. The black arrow depicts the selection for the signal region. An example SUSY scenario is overlaid for illustration.

The post-fit kinematic distribution for the Intermediate $Wh$ channel, showing the $m_{\mathrm{T}2}$ distribution in SR-Wh-LM, before the selection on $m_{\mathrm{T}2}$ is made. The black arrow depicts the selection for the signal region. An example SUSY scenario is overlaid for illustration.

The post-fit kinematic distribution for the Intermediate $Wh$ channel, showing the $m_{\mathrm{Tsum}}$ distribution in SR-Wh-HM, before the selection on $m_{\mathrm{Tsum}}$ is made. The black arrow depicts the selection for the signal region. An example SUSY scenario is overlaid for illustration.

Expected 95\% CL exclusion limits on the simplified models of $\tilde{\tau}_{L,R}\tilde{\tau}_{L,R}$ production.

$+1\sigma$ expected 95\% CL exclusion limits on the simplified models of $\tilde{\tau}_{L,R}\tilde{\tau}_{L,R}$ production.

$-1\sigma$ expected 95\% CL exclusion limits on the simplified models of $\tilde{\tau}_{L,R}\tilde{\tau}_{L,R}$ production.

Observed 95\% CL exclusion limits on the simplified models of $\tilde{\tau}_{L,R}\tilde{\tau}_{L,R}$ production.

$+1\sigma$ observed 95\% CL exclusion limits on the simplified models of $\tilde{\tau}_{L,R}\tilde{\tau}_{L,R}$ production.

$-1\sigma$ observed 95\% CL exclusion limits on the simplified models of $\tilde{\tau}_{L,R}\tilde{\tau}_{L,R}$ production.

Expected 95\% CL exclusion limits on the simplified models of $\tilde{\tau}_{L}\tilde{\tau}_{L}$ production.

$+1\sigma$ expected 95\% CL exclusion limits on the simplified models of $\tilde{\tau}_{L}\tilde{\tau}_{L}$ production.

$-1\sigma$ expected 95\% CL exclusion limits on the simplified models of $\tilde{\tau}_{L}\tilde{\tau}_{L}$ production.

Observed 95\% CL exclusion limits on the simplified models of $\tilde{\tau}_{L}\tilde{\tau}_{L}$ production.

$+1\sigma$ observed 95\% CL exclusion limits on the simplified models of $\tilde{\tau}_{L}\tilde{\tau}_{L}$ production.

$-1\sigma$ observed 95\% CL exclusion limits on the simplified models of $\tilde{\tau}_{L}\tilde{\tau}_{L}$ production.

Expected 95\% CL exclusion limits on the simplified models of $\tilde{\tau}_{R}\tilde{\tau}_{R}$ production.

$+1\sigma$ expected 95\% CL exclusion limits on the simplified models of $\tilde{\tau}_{R}\tilde{\tau}_{R}$ production.

$-1\sigma$ expected 95\% CL exclusion limits on the simplified models of $\tilde{\tau}_{R}\tilde{\tau}_{R}$ production.

Observed 95\% CL exclusion limits on the simplified models of $\tilde{\tau}_{R}\tilde{\tau}_{R}$ production.

$+1\sigma$ observed 95\% CL exclusion limits on the simplified models of $\tilde{\tau}_{R}\tilde{\tau}_{R}$ production.

$-1\sigma$ observed 95\% CL exclusion limits on the simplified models of $\tilde{\tau}_{R}\tilde{\tau}_{R}$ production.

Expected 95\% CL exclusion limits on the simplified models of $\tilde{\chi}_1^{+}\tilde{\chi}_1^{-}$ production decaying via staus.

$+1\sigma$ expected 95\% CL exclusion limits on the simplified models of $\tilde{\chi}_1^{+}\tilde{\chi}_1^{-}$ production decaying via staus.

$-1\sigma$ expected 95\% CL exclusion limits on the simplified models of $\tilde{\chi}_1^{+}\tilde{\chi}_1^{-}$ production decaying via staus.

Observed 95\% CL exclusion limits on the simplified models of $\tilde{\chi}_1^{+}\tilde{\chi}_1^{-}$ production decaying via staus.

Expected 95\% CL exclusion limits on the simplified models of $\tilde{\chi}_1^{\pm}\tilde{\chi}_2^0$ production decaying via staus (parameterised as $m(\tilde{\chi}_1^{\pm}/\tilde{\chi}_2^0)$ and $m(\tilde{\chi}_1^0)$).

$+1\sigma$ expected 95\% CL exclusion limits on the simplified models of $\tilde{\chi}_1^{\pm}\tilde{\chi}_2^0$ production decaying via staus (parameterised as $m(\tilde{\chi}_1^{\pm}/\tilde{\chi}_2^0)$ and $m(\tilde{\chi}_1^0)$).

$-1\sigma$ expected 95\% CL exclusion limits on the simplified models of $\tilde{\chi}_1^{\pm}\tilde{\chi}_2^0$ production decaying via staus (parameterised as $m(\tilde{\chi}_1^{\pm}/\tilde{\chi}_2^0)$ and $m(\tilde{\chi}_1^0)$).

Observed 95\% CL exclusion limits on the simplified models of $\tilde{\chi}_1^{\pm}\tilde{\chi}_2^0$ production decaying via staus (parameterised as $m(\tilde{\chi}_1^{\pm}/\tilde{\chi}_2^0)$ and $m(\tilde{\chi}_1^0)$).

Expected 95\% CL exclusion limits on the simplified models of $\tilde{\chi}_1^{\pm}\tilde{\chi}_2^0$ production decaying via staus (parameterised as $m(\tilde{\chi}_1^{\pm}/\tilde{\chi}_2^0)$ and $m(\tilde{\chi}_1^0)$) using the OS signal regions.

Observed 95\% CL exclusion limits on the simplified models of $\tilde{\chi}_1^{\pm}\tilde{\chi}_2^0$ production decaying via staus (parameterised as $m(\tilde{\chi}_1^{\pm}/\tilde{\chi}_2^0)$ and $m(\tilde{\chi}_1^0)$) using the OS signal regions.

Expected 95\% CL exclusion limits on the simplified models of $\tilde{\chi}_1^{\pm}\tilde{\chi}_2^0$ production decaying via staus (parameterised as $m(\tilde{\chi}_1^{\pm}/\tilde{\chi}_2^0)$ and $m(\tilde{\chi}_1^0)$) using the SS signal regions.

Observed 95\% CL exclusion limits on the simplified models of $\tilde{\chi}_1^{\pm}\tilde{\chi}_2^0$ production decaying via staus (parameterised as $m(\tilde{\chi}_1^{\pm}/\tilde{\chi}_2^0)$ and $m(\tilde{\chi}_1^0)$) using the SS signal regions.

Expected 95\% CL exclusion limits on the simplified models of $\tilde{\chi}_1^{\pm}\tilde{\chi}_2^0$ production decaying via $Wh$ (parameterised as $m(\tilde{\chi}_1^{\pm}/\tilde{\chi}_2^0)$ and $m(\tilde{\chi}_2^0)-m(\tilde{\chi}_1^0)$).

$+1\sigma$ expected 95\% CL exclusion limits on the simplified models of $\tilde{\chi}_1^{\pm}\tilde{\chi}_2^0$ production decaying via $Wh$ (parameterised as $m(\tilde{\chi}_1^{\pm}/\tilde{\chi}_2^0)$ and $m(\tilde{\chi}_2^0)-m(\tilde{\chi}_1^0)$).

$-1\sigma$ expected 95\% CL exclusion limits on the simplified models of $\tilde{\chi}_1^{\pm}\tilde{\chi}_2^0$ production decaying via $Wh$ (parameterised as $m(\tilde{\chi}_1^{\pm}/\tilde{\chi}_2^0)$ and $m(\tilde{\chi}_2^0)-m(\tilde{\chi}_1^0)$).

Observed 95\% CL exclusion limits on the simplified models of $\tilde{\chi}_1^{\pm}\tilde{\chi}_2^0$ production decaying via $Wh$ (parameterised as $m(\tilde{\chi}_1^{\pm}/\tilde{\chi}_2^0)$ and $m(\tilde{\chi}_2^0)-m(\tilde{\chi}_1^0)$).

Expected 95\% CL exclusion limits on the simplified models of $\tilde{\chi}_1^{\pm}\tilde{\chi}_2^0$ production decaying via $Wh$ (parameterised as $m(\tilde{\chi}_1^{\pm}/\tilde{\chi}_2^0)$ and $m(\tilde{\chi}_2^0)-m(\tilde{\chi}_1^0)$) using the OS signal regions.

Observed 95\% CL exclusion limits on the simplified models of $\tilde{\chi}_1^{\pm}\tilde{\chi}_2^0$ production decaying via $Wh$ (parameterised as $m(\tilde{\chi}_1^{\pm}/\tilde{\chi}_2^0)$ and $m(\tilde{\chi}_2^0)-m(\tilde{\chi}_1^0)$) using the OS signal regions.

Expected 95\% CL exclusion limits on the simplified models of $\tilde{\chi}_1^{\pm}\tilde{\chi}_2^0$ production decaying via $Wh$ (parameterised as $m(\tilde{\chi}_1^{\pm}/\tilde{\chi}_2^0)$ and $m(\tilde{\chi}_2^0)-m(\tilde{\chi}_1^0)$) using the SS signal regions.

Observed 95\% CL exclusion limits on the simplified models of $\tilde{\chi}_1^{\pm}\tilde{\chi}_2^0$ production decaying via $Wh$ (parameterised as $m(\tilde{\chi}_1^{\pm}/\tilde{\chi}_2^0)$ and $m(\tilde{\chi}_2^0)-m(\tilde{\chi}_1^0)$) using the SS signal regions.

Expected 95\% CL exclusion limits on the simplified models of $\tilde{\chi}_1^{\pm}\tilde{\chi}_2^0$ production decaying via $Wh$.

$+1\sigma$ expected 95\% CL exclusion limits on the simplified models of $\tilde{\chi}_1^{\pm}\tilde{\chi}_2^0$ production decaying via $Wh$.

$-1\sigma$ expected 95\% CL exclusion limits on the simplified models of $\tilde{\chi}_1^{\pm}\tilde{\chi}_2^0$ production decaying via $Wh$.

Observed 95\% CL exclusion limits on the simplified models of $\tilde{\chi}_1^{\pm}\tilde{\chi}_2^0$ production decaying via $Wh$.

Observed upper limit on the signal cross section in fb for the production of $\tilde{\tau}_{L,R}\tilde{\tau}_{L,R}$.

Observed upper limit on the signal cross section in fb for the production of $\tilde{\tau}_{L}\tilde{\tau}_{L}$.

Observed upper limit on the signal cross section in fb for the production of $\tilde{\tau}_{R}\tilde{\tau}_{R}$.

The best expected signal region used to set limits in models of the production of $\tilde{\tau}_{L,R}\tilde{\tau}_{L,R}$.

The best expected signal region used to set limits in models of the production of $\tilde{\tau}_{L}\tilde{\tau}_{L}$.

The best expected signal region used to set limits in models of the production of $\tilde{\tau}_{R}\tilde{\tau}_{R}$.

The acceptance $A$ of the direct stau production signal region SR-BDT1 in models of the production of $\tilde{\tau}_{L,R}\tilde{\tau}_{L,R}$.

The efficiency $\epsilon$ of the direct stau production signal region SR-BDT1 in models of the production of $\tilde{\tau}_{L,R}\tilde{\tau}_{L,R}$.

The acceptance $A$ of the direct stau production signal region SR-BDT2 in models of the production of $\tilde{\tau}_{L,R}\tilde{\tau}_{L,R}$.

The efficiency $\epsilon$ of the direct stau production signal region SR-BDT2 in models of the production of $\tilde{\tau}_{L,R}\tilde{\tau}_{L,R}$.

The acceptance $A$ of the direct stau production signal region SR-BDT3 in models of the production of $\tilde{\tau}_{L,R}\tilde{\tau}_{L,R}$.

The efficiency $\epsilon$ of the direct stau production signal region SR-BDT3 in models of the production of $\tilde{\tau}_{L,R}\tilde{\tau}_{L,R}$.

The acceptance $A$ of the direct stau production signal region SR-BDT4 in models of the production of $\tilde{\tau}_{L,R}\tilde{\tau}_{L,R}$.

The efficiency $\epsilon$ of the direct stau production signal region SR-BDT4 in models of the production of $\tilde{\tau}_{L,R}\tilde{\tau}_{L,R}$.

Cutflow event yields of the direct stau production SRs in models of the production of $\tilde{\tau}_{L,R}\tilde{\tau}_{L,R}$. All yields correspond to weighted events, so that effects from reconstruction efficiencies, trigger corrections, pileup reweighting, etc., are included. They are normalised to the integrated luminosity of the data sample, $\int L dt = 139$ fb$^{-1}$. The preliminary event reduction is a centralised stage requiring the event to: have at least two hadronically-decaying taus; be selected and matched by the asymmetric ditau trigger (with plateau cuts); have at least two medium taus of opposite sign and have $m_{\mathrm{T2}} > 30$ GeV. Cuts except the BDT score cuts are applied consecutively, while the BDT score cuts are applied one at a time.

Observed upper limit on the signal cross section in fb for the production of $\tilde{\chi}_1^{+}\tilde{\chi}_1^{-}$.

Observed upper limit on the signal cross section in fb for the production of $\tilde{\chi}_1^{\pm}\tilde{\chi}_2^{0}$.

The acceptance $A$ of the Intermediate stau signal region SR-C1C1-LM in models of the production of $\tilde{\chi}_1^{+}\tilde{\chi}_1^{-}$ decaying via staus.

The efficiency $\epsilon$ of the Intermediate stau signal region SR-C1C1-LM in models of the production of $\tilde{\chi}_1^{+}\tilde{\chi}_1^{-}$ decaying via staus.

The acceptance $A$ of the Intermediate stau signal region SR-C1C1-HM in models of the production of $\tilde{\chi}_1^{+}\tilde{\chi}_1^{-}$ decaying via staus.

The efficiency $\epsilon$ of the Intermediate stau signal region SR-C1C1-HM in models of the production of $\tilde{\chi}_1^{+}\tilde{\chi}_1^{-}$ decaying via staus.

The acceptance $A$ of the Intermediate stau signal region SR-C1N2OS-LM in models of the production of $\tilde{\chi}_1^{\pm}\tilde{\chi}_2^{0}$ decaying via staus.

The efficiency $\epsilon$ of the Intermediate stau signal region SR-C1N2OS-LM in models of the production of $\tilde{\chi}_1^{\pm}\tilde{\chi}_2^{0}$ decaying via staus.

The acceptance $A$ of the Intermediate stau signal region SR-C1N2SS-LM in models of the production of $\tilde{\chi}_1^{\pm}\tilde{\chi}_2^{0}$ decaying via staus.

The efficiency $\epsilon$ of the Intermediate stau signal region SR-C1N2SS-LM in models of the production of $\tilde{\chi}_1^{\pm}\tilde{\chi}_2^{0}$ decaying via staus.

The acceptance $A$ of the Intermediate stau signal region SR-C1N2OS-HM in models of the production of $\tilde{\chi}_1^{\pm}\tilde{\chi}_2^{0}$ decaying via staus.

The efficiency $\epsilon$ of the Intermediate stau signal region SR-C1N2OS-HM in models of the production of $\tilde{\chi}_1^{\pm}\tilde{\chi}_2^{0}$ decaying via staus.

The acceptance $A$ of the Intermediate stau signal region SR-C1N2SS-HM in models of the production of $\tilde{\chi}_1^{\pm}\tilde{\chi}_2^{0}$ decaying via staus.

The efficiency $\epsilon$ of the Intermediate stau signal region SR-C1N2SS-HM in models of the production of $\tilde{\chi}_1^{\pm}\tilde{\chi}_2^{0}$ decaying via staus.

Cutflow event yields of the Intermediate stau SRs in models of the production of $\tilde{\chi}_1^{+}\tilde{\chi}_1^{-}$. All yields correspond to weighted events, so that effects from reconstruction efficiencies, trigger corrections, pileup reweighting, etc., are included. They are normalised to the integrated luminosity of the data sample, $\int L dt = 139$ fb$^{-1}$. The preliminary event reduction is a centralised stage requiring the event to: have at least two hadronically-decaying medium taus; be selected and matched by the asymmetric ditau trigger or ditau+MET trigger with offline cuts using HLT threshold values. The different yields of preliminary event reduction are caused by different trigger scale factor.

Cutflow event yields of the Intermediate stau SRs in models of the production of $\tilde{\chi}_1^{\pm}\tilde{\chi}_2^{0}$. All yields correspond to weighted events, so that effects from reconstruction efficiencies, trigger corrections, pileup reweighting, etc., are included. They are normalised to the integrated luminosity of the data sample, $\int L dt = 139$ fb$^{-1}$. The preliminary event reduction is a centralised stage requiring the event to: have at least two hadronically-decaying medium taus; be selected and matched by the asymmetric ditau trigger or ditau+MET trigger with offline cuts using HLT threshold values. The different yields of preliminary event reduction are caused by different trigger scale factor.

Observed upper limit on the signal cross section in fb for the production of $\tilde{\chi}_1^{\pm}\tilde{\chi}_2^0$ decaying via $Wh$.

The best expected signal region used to set limits in models of the production of $\tilde{\chi}_1^{\pm}\tilde{\chi}_2^0$ decaying via $Wh$.

The acceptance $A$ of the Intermediate $Wh$ signal region SR-Wh-LM in models of the production of $\tilde{\chi}_1^{\pm}\tilde{\chi}_2^0$ decaying via $Wh$.

The efficiency $\epsilon$ of the Intermediate $Wh$ signal region SR-Wh-LM in models of the production of $\tilde{\chi}_1^{\pm}\tilde{\chi}_2^0$ decaying via $Wh$.

The acceptance $A$ of the Intermediate $Wh$ signal region SR-Wh-HM in models of the production of $\tilde{\chi}_1^{\pm}\tilde{\chi}_2^0$ decaying via $Wh$.

The efficiency $\epsilon$ of the Intermediate $Wh$ signal region SR-Wh-HM in models of the production of $\tilde{\chi}_1^{\pm}\tilde{\chi}_2^0$ decaying via $Wh$.

Cutflow event yields of the Intermediate $Wh$ SRs in models of the production of $\tilde{\chi}_1^{\pm}\tilde{\chi}_2^0$ decaying via $Wh$. All yields correspond to weighted events, so that effects from reconstruction efficiencies, trigger corrections, pileup reweighting, etc., are included. They are normalised to the integrated luminosity of the data sample, $\int L dt = 139$ fb$^{-1}$. The preliminary event reduction is a centralised stage requiring the event to: have one single lepton and two hadronically decaying medium taus; be selected and matched by the single lepton trigger with plateau cuts.

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.

Covariance matrix of the two parameters $A$ and $m$ of the linear 2-to-5-jets-bin adjustment in the background Monte Carlo. Please see the text of the paper for an explanation of $A$ and $m$. Best fit values of $A$ and m are available as a separate table in this HEPData record.

Best fit values of the two parameters $A$ and $m$ of the linear 2-to-6-jets-bin adjustment in the background Monte Carlo. Please see the text of the paper for an explanation of $A$ and $m$. Uncertainties on $A$ and $m$ can be obtained from the covariance matrix of the fit, available as a separate table in this HEPData Record.

Covariance matrix of the two parameters $A$ and $m$ of the linear 2-to-6-jets-bin adjustment in the background Monte Carlo. Please see the text of the paper for an explanation of $A$ and $m$. Best fit values of $A$ and m are available as a separate table in this HEPData record.

This table contains the pre-fit background model (i.e. setting all nuisance parameters to 0), the post-fit background model (i.e. setting all nuisance parameters to their best fit values given the observed data), and the observed data in the 4 Jets bin. Uncertainties are provided, but per-bin uncertainties do not capture the correlations between bins with nuisance parameter variations. For a full analysis, the full covariance matrix between the search regions, also present in this HEPData record, should be used.

This table contains the pre-fit background model (i.e. setting all nuisance parameters to 0), the post-fit background model (i.e. setting all nuisance parameters to their best fit values given the observed data), and the observed data in the 5 Jets bin. Uncertainties are provided, but per-bin uncertainties do not capture the correlations between bins with nuisance parameter variations. For a full analysis, the full covariance matrix between the search regions, also present in this HEPData record, should be used.

This table contains the pre-fit background model (i.e. setting all nuisance parameters to 0), the post-fit background model (i.e. setting all nuisance parameters to their best fit values given the observed data), and the observed data in the $\geq 6$ Jets bin. Uncertainties are provided, but per-bin uncertainties do not capture the correlations between bins with nuisance parameter variations. For a full analysis, the full covariance matrix between the search regions, also present in this HEPData record, should be used.

Covariance matrix between search regions, estimated by toy MC runs, varying the nuisance parameters according to their covariance matrix before performing any fit. Note that the two columns ($S_T$ Bin (axis 1)) and ($N_{jets}$ Bin (axis 1)) together define the x-axis of the covariance, and the two columns ($S_T$ Bin (axis 2)) and ($N_{jets}$ Bin (axis 2)) together define the y-axis of the covariance.

Covariance matrix between search regions, estimated by toy MC runs, varying the nuisance parameters according to their covariance matrix in the background-only fit. Note that the two columns ($S_T$ Bin (axis 1)) and (NJets Bin (axis 1)) together define the x-axis of the covariance, and the two columns ($S_T$ Bin (axis 2)) and (NJets Bin (axis 2)) together define the y-axis of the covariance.

Scan through (gluino mass, neutralino mass) parameter-space. This table has the gluino and neutralino mass as the independent parameters, and the following as dependent parameters: 1. Expected upper limits (95% $CL_s$) on signal strength, with $\pm 1\sigma$ and $\pm 2\sigma$ uncertainty terms from limited statistics. 2. Observed upper limits (95% $CL_s$) on signal strength, with a $\pm 1\sigma$ uncertainty from the theoretical prediction of the gluino production cross section. 3. Observed upper limits (95% $CL_s$) on gluino production cross section.

Scan through (squark mass, neutralino mass) parameter-space. This table has the squark and neutralino mass as the independent parameters, and the following as dependent parameters: 1. Expected upper limits (95% $CL_s$) on signal strength, with $\pm 1\sigma$ and $\pm 2\sigma$ uncertainty terms from limited statistics. 2. Observed upper limits (95% $CL_s$) on signal strength, with a $\pm 1\sigma$ uncertainty from the theoretical prediction of the squark production cross section. 3. Observed upper limits (95% $CL_s$) on squark production cross section.

The 95% CL expected upper limits on the production cross sections for $m_{\tilde{g}}$ in T5bbbbZg.

Gluino - neutralino mass points lying on the expected exclusion contour of T5bbbbZg model.

The 95% CL observed upper limits on the production cross sections for $m_{\tilde{g}}$ in T5qqqqHg.

Gluino - neutralino mass points lying on the observed exclusion contour of T5qqqqHg model.

The 95% CL expected upper limits on the production cross sections for $m_{\tilde{g}}$ in T5qqqqHg.

Gluino - neutralino mass points lying on the expected exclusion contour of T5qqqqHg model.

The 95% CL observed upper limits on the production cross sections for $m_{\tilde{g}}$ in T5ttttZg.

Gluino - neutralino mass points lying on the observed exclusion contour of T5ttttZg model.

The 95% CL expected upper limits on the production cross sections for $m_{\tilde{g}}$ in T5ttttZg.

Gluino - neutralino mass points lying on the expected exclusion contour of T5ttttZg model.

The 95% CL observed upper limits on the production cross sections for $m_{\tilde{t}}$ in T6ttZg.

Stop - neutralino mass points lying on the observed exclusion contour of T6ttZg model.

The 95% CL expected upper limits on the production cross sections for $m_{\tilde{t}}$ in T6ttZg.

Stop - neutralino mass points lying on the expected exclusion contour of T6ttZg model.

observed exclusion limit on the neutralino mass points in the TChiNGnn model.

expected exclusion limit on the neutralino mass points in the TChiNGnn model.

observed exclusion limit on the neutralino mass points in the TChiNG model.

expected exclusion limit on the neutralino mass points in the TChiNG model.

observed exclusion limit on the neutralino mass points in the TChiWG model.

expected exclusion limit on the neutralino mass points in the TChiWG model.

Acceptance times efficiency values with statistical uncertainties for T5bbbbZg in the SR region.

Acceptance times efficiency values with statistical uncertainties for T5qqqqHg in the SR region.

Acceptance times efficiency values with statistical uncertainties for T5ttttZg in the SR region.

Acceptance times efficiency values with statistical uncertainties for T6ttZg in the SR region.

Acceptance times efficiency values with statistical uncertainties for T5bbbbZg inclusive in the 45 bins.

Acceptance times efficiency values with statistical uncertainties for T5qqqqHg inclusive in the 45 bins.

Acceptance times efficiency values with statistical uncertainties for T5ttttZg inclusive in the 45 bins.

Acceptance times efficiency values with statistical uncertainties for T6ttZg inclusive in the 45 bins.

Acceptance times efficiency values with statistical uncertainties for TChiWG in the 45 optimised bins region.

Acceptance times efficiency values with statistical uncertainties for TChiNG in the 45 optimised bins region.

Event yield for T5bbbbZg and T5qqqqHg for each bin after SR baseline selections.

Event yield for T5ttttZg and T6ttZg for each bin after SR baseline selections.

Event yield for TChiWg and TChiNg for each bin after SR baseline selections.

Expected exclusion limits at 95% CL from Fig 7(a) for $\tilde{g}$ decays into SM gauge bosons and $\tilde{\chi}^{0}_{1}$

One $\sigma$ band of expected exclusion limits at 95% CL from Fig 7(a) for $\tilde{g}$ decays into SM gauge bosons and $\tilde{\chi}^{0}_{1}$

Observed exclusion limits at 95% CL from Fig 7(c) for $\tilde{g}$ decays into sleptons and subsequently to SM leptons and $\tilde{\chi}^{0}_{1}$

Positive one $\sigma$ observed exclusion limits at 95% CL from Fig 7(c) for $\tilde{g}$ decays into sleptons and subsequently to SM leptons and $\tilde{\chi}^{0}_{1}$

Negative one $\sigma$ observed exclusion limits at 95% CL from Fig 7(c) for $\tilde{g}$ decays into sleptons and subsequently to SM leptons and $\tilde{\chi}^{0}_{1}$

Expected exclusion limits at 95% CL from Fig 7(c) for $\tilde{g}$ decays into sleptons and subsequently to SM leptons and $\tilde{\chi}^{0}_{1}$

One $\sigma$ band of expected exclusion limits at 95% CL from Fig 7(c) for $\tilde{g}$ decays into sleptons and subsequently to SM leptons and $\tilde{\chi}^{0}_{1}$

Observed exclusion limits at 95% CL from Fig 7(f) for $\tilde{g}$ decays into anti-top and $\tilde{t}$ and $\tilde{t}$ decays via a non-zero RPV coupling $\lambda''$

Positive one $\sigma$ observed exclusion limits at 95% CL from Fig 7(f) for $\tilde{g}$ decays into anti-top and $\tilde{t}$ and $\tilde{t}$ decays via a non-zero RPV coupling $\lambda''$

Negative one $\sigma$ observed exclusion limits at 95% CL from Fig 7(f) for $\tilde{g}$ decays into anti-top and $\tilde{t}$ and $\tilde{t}$ decays via a non-zero RPV coupling $\lambda''$

Expected exclusion limits at 95% CL from Fig 7(f) for $\tilde{g}$ decays into anti-top and $\tilde{t}$ and $\tilde{t}$ decays via a non-zero RPV coupling $\lambda''$

One $\sigma$ band of expected exclusion limits at 95% CL from Fig 7(f) for $\tilde{g}$ decays into anti-top and $\tilde{t}$ and $\tilde{t}$ decays via a non-zero RPV coupling $\lambda''$

Observed exclusion limits at 95% CL from Fig 7(e) for direct $\tilde{\chi_{1}^{0}}$ decay into SM leptons and quarks via a non-zero RPV coupling $\lambda'$

Positive one $\sigma$ observed exclusion limits at 95% CL from Fig 7(e) for direct $\tilde{\chi_{1}^{0}}$ decay into SM leptons and quarks via a non-zero RPV coupling $\lambda'$

Negative one $\sigma$ observed exclusion limits at 95% CL from Fig 7(e) for direct $\tilde{\chi_{1}^{0}}$ decay into SM leptons and quarks via a non-zero RPV coupling $\lambda'$

Expected exclusion limits at 95% CL from Fig 7(e) for direct $\tilde{\chi_{1}^{0}}$ decay into SM leptons and quarks via a non-zero RPV coupling $\lambda'$

One $\sigma$ band of expected exclusion limits at 95% CL from Fig 7(e) for direct $\tilde{\chi_{1}^{0}}$ decay into SM leptons and quarks via a non-zero RPV coupling $\lambda'$

Observed exclusion limits at 95% CL from Fig 7(b) for $\tilde{q}$ decays into SM gauge bosons and $\tilde{\chi}^{0}_{1}$

Positive one $\sigma$ observed exclusion limits at 95% CL from Fig 7(b) for $\tilde{q}$ decays into SM gauge bosons and $\tilde{\chi}^{0}_{1}$

Negative one $\sigma$ observed exclusion limits at 95% CL from Fig 7(b) for $\tilde{q}$ decays into SM gauge bosons and $\tilde{\chi}^{0}_{1}$

Expected exclusion limits at 95% CL from Fig 7(b) for $\tilde{q}$ decays into SM gauge bosons and $\tilde{\chi}^{0}_{1}$

One $\sigma$ band of expected exclusion limits at 95% CL from Fig 7(b) for $\tilde{q}$ decays into SM gauge bosons and $\tilde{\chi}^{0}_{1}$

Observed exclusion limits at 95% CL from Fig 7(d) for $\tilde{q}$ decays into sleptons and subsequently to SM leptons and $\tilde{\chi}^{0}_{1}$

Positive one $\sigma$ observed exclusion limits at 95% CL from Fig 7(d) for $\tilde{q}$ decays into sleptons and subsequently to SM leptons and $\tilde{\chi}^{0}_{1}$

Negative one $\sigma$ observed exclusion limits at 95% CL from Fig 7(d) for $\tilde{q}$ decays into sleptons and subsequently to SM leptons and $\tilde{\chi}^{0}_{1}$

Expected exclusion limits at 95% CL from Fig 7(d) for $\tilde{q}$ decays into sleptons and subsequently to SM leptons and $\tilde{\chi}^{0}_{1}$

One $\sigma$ band of expected exclusion limits at 95% CL from Fig 7(d) for $\tilde{q}$ decays into sleptons and subsequently to SM leptons and $\tilde{\chi}^{0}_{1}$

N-1 distribution for $m_{\mathrm{eff}}$of observed data and expected background in SRGGWZ-H.

N-1 distribution for $E_{\mathrm{T}}^{\mathrm{miss}}$of observed data and expected background in SRGGSlep-M.

N-1 distribution for $\sum{p_{\mathrm{T}}^\mathrm{jet}}$of observed data and expected background in SRUDD-ge2b.

N-1 distribution for $m_{\mathrm{eff}}$of observed data and expected background in SRLQD.

N-1 distribution for $m_{\mathrm{eff}}$of observed data and expected background in SRSSWZ-H.

N-1 distribution for $m_{\mathrm{eff}}$of observed data and expected background in SRSSSlep-H(loose).

Signal acceptance for SRGGWZ-H signal region from Fig 10(c) in a SUSY scenario where $\tilde{g}$ decays into SM gauge bosons and $\tilde{\chi}^{0}_{1}$

Signal efficiency for SRGGWZ-H signal region from Fig 15(c) in a SUSY scenario where $\tilde{g}$ decays into SM gauge bosons and $\tilde{\chi}^{0}_{1}$

Signal acceptance for SRGGWZ-M signal region from Fig 10(b) in a SUSY scenario where $\tilde{g}$ decays into SM gauge bosons and $\tilde{\chi}^{0}_{1}$

Signal efficiency for SRGGWZ-M signal region from Fig 15(b) in a SUSY scenario where $\tilde{g}$ decays into SM gauge bosons and $\tilde{\chi}^{0}_{1}$

Signal acceptance for SRGGWZ-L signal region from Fig 10(a) in a SUSY scenario where $\tilde{g}$ decays into SM gauge bosons and $\tilde{\chi}^{0}_{1}$

Signal efficiency for SRGGWZ-L signal region from Fig 15(a) in a SUSY scenario where $\tilde{g}$ decays into SM gauge bosons and $\tilde{\chi}^{0}_{1}$

Signal acceptance for SRGGSlep-L signal region from Fig 12(a) in a SUSY scenario where $\tilde{g}$ decays into sleptons and subsequently to SM leptons and $\tilde{\chi}^{0}_{1}$

Signal efficiency for SRGGSlep-L signal region from Fig 17(a) in a SUSY scenario where $\tilde{g}$ decays into sleptons and subsequently to SM leptons and $\tilde{\chi}^{0}_{1}$

Signal acceptance for SRGGSlep-M signal region from Fig 12(b) in a SUSY scenario where $\tilde{g}$ decays into sleptons and subsequently to SM leptons and $\tilde{\chi}^{0}_{1}$

Signal efficiency for SRGGSlep-M signal region from Fig 17(b) in a SUSY scenario where $\tilde{g}$ decays into sleptons and subsequently to SM leptons and $\tilde{\chi}^{0}_{1}$

Signal acceptance for SRGGSlep-H signal region from Fig 12(c) in a SUSY scenario where $\tilde{g}$ decays into sleptons and subsequently to SM leptons and $\tilde{\chi}^{0}_{1}$

Signal efficiency for SRGGSlep-H signal region from Fig 17(c) in a SUSY scenario where $\tilde{g}$ decays into sleptons and subsequently to SM leptons and $\tilde{\chi}^{0}_{1}$

Signal acceptance for SRUDD-1b signal region from Fig 14(b) in a SUSY scenario where $\tilde{g}$ decays into anti-top and $\tilde{t}$ and $\tilde{t}$ decays via a non-zero RPV coupling $\lambda''$

Signal efficiency for SRUDD-1b signal region from Fig 19(b) in a SUSY scenario where $\tilde{g}$ decays into anti-top and $\tilde{t}$ and $\tilde{t}$ decays via a non-zero RPV coupling $\lambda''$

Signal acceptance for SRUDD-2b signal region from Fig 14(c) in a SUSY scenario where $\tilde{g}$ decays into anti-top and $\tilde{t}$ and $\tilde{t}$ decays via a non-zero RPV coupling $\lambda''$

Signal efficiency for SRUDD-2b signal region from Fig 19(c) in a SUSY scenario where $\tilde{g}$ decays into anti-top and $\tilde{t}$ and $\tilde{t}$ decays via a non-zero RPV coupling $\lambda''$

Signal acceptance for SRUDD-ge2b signal region from Fig 14(d) in a SUSY scenario where $\tilde{g}$ decays into anti-top and $\tilde{t}$ and $\tilde{t}$ decays via a non-zero RPV coupling $\lambda''$

Signal efficiency for SRUDD-ge2b signal region from Fig 19(d) in a SUSY scenario where $\tilde{g}$ decays into anti-top and $\tilde{t}$ and $\tilde{t}$ decays via a non-zero RPV coupling $\lambda''$

Signal acceptance for SRUDD-ge3b signal region from Fig 14(e) in a SUSY scenario where $\tilde{g}$ decays into anti-top and $\tilde{t}$ and $\tilde{t}$ decays via a non-zero RPV coupling $\lambda''$

Signal efficiency for SRUDD-ge3b signal region from Fig 19(e) in a SUSY scenario where $\tilde{g}$ decays into anti-top and $\tilde{t}$ and $\tilde{t}$ decays via a non-zero RPV coupling $\lambda''$

Signal acceptance for SRLQD signal region from Fig 14(a) in a SUSY scenario where direct $\tilde{\chi_{1}^{0}}$ decay into SM leptons and quarks via a non-zero RPV coupling $\lambda'$

Signal efficiency for SRLQD signal region from Fig 19(a) in a SUSY scenario where direct $\tilde{\chi_{1}^{0}}$ decay into SM leptons and quarks via a non-zero RPV coupling $\lambda'$

Signal acceptance for SRSSWZ-L signal region from Fig 11(a) in a SUSY scenario where $\tilde{q}$ decays into SM gauge bosons and $\tilde{\chi}^{0}_{1}$

Signal efficiency for SRSSWZ-L signal region from Fig 16(a) in a SUSY scenario where $\tilde{q}$ decays into SM gauge bosons and $\tilde{\chi}^{0}_{1}$

Signal acceptance for SRSSWZ-ML signal region from Fig 11(b) in a SUSY scenario where $\tilde{q}$ decays into SM gauge bosons and $\tilde{\chi}^{0}_{1}$

Signal efficiency for SRSSWZ-ML signal region from Fig 16(b) in a SUSY scenario where $\tilde{q}$ decays into SM gauge bosons and $\tilde{\chi}^{0}_{1}$

Signal acceptance for SRSSWZ-MH signal region from Fig 11(c) in a SUSY scenario where $\tilde{q}$ decays into SM gauge bosons and $\tilde{\chi}^{0}_{1}$

Signal efficiency for SRSSWZ-MH signal region from Fig 16(c) in a SUSY scenario where $\tilde{q}$ decays into SM gauge bosons and $\tilde{\chi}^{0}_{1}$

Signal acceptance for SRSSWZ-H signal region from Fig 11(d) in a SUSY scenario where $\tilde{q}$ decays into SM gauge bosons and $\tilde{\chi}^{0}_{1}$

Signal efficiency for SRSSWZ-H signal region from Fig 16(d) in a SUSY scenario where $\tilde{q}$ decays into SM gauge bosons and $\tilde{\chi}^{0}_{1}$

Signal acceptance for SRSSSlep-H signal region from Fig 13(d) in a SUSY scenario where $\tilde{q}$ decays into sleptons and subsequently to SM leptons and $\tilde{\chi}^{0}_{1}$

Signal efficiency for SRSSSlep-H signal region from Fig 18(d) in a SUSY scenario where $\tilde{q}$ decays into sleptons and subsequently to SM leptons and $\tilde{\chi}^{0}_{1}$

Signal acceptance for SRSSSlep-MH signal region from Fig 13(c) in a SUSY scenario where $\tilde{q}$ decays into sleptons and subsequently to SM leptons and $\tilde{\chi}^{0}_{1}$

Signal efficiency for SRSSSlep-MH signal region from Fig 18(c) in a SUSY scenario where $\tilde{q}$ decays into sleptons and subsequently to SM leptons and $\tilde{\chi}^{0}_{1}$

Signal acceptance for SRSSSlep-L signal region from Fig 13(a) in a SUSY scenario where $\tilde{q}$ decays into sleptons and subsequently to SM leptons and $\tilde{\chi}^{0}_{1}$

Signal efficiency for SRSSSlep-L signal region from Fig 18(a) in a SUSY scenario where $\tilde{q}$ decays into sleptons and subsequently to SM leptons and $\tilde{\chi}^{0}_{1}$

Signal acceptance for SRSSSlep-ML signal region from Fig 13(b) in a SUSY scenario where $\tilde{q}$ decays into sleptons and subsequently to SM leptons and $\tilde{\chi}^{0}_{1}$

Signal efficiency for SRSSSlep-ML signal region from Fig 18(b) in a SUSY scenario where $\tilde{q}$ decays into sleptons and subsequently to SM leptons and $\tilde{\chi}^{0}_{1}$

Signal acceptance for SRSSSlep-H(loose) signal region from Fig 13(e) in a SUSY scenario where $\tilde{q}$ decays into sleptons and subsequently to SM leptons and $\tilde{\chi}^{0}_{1}$

Signal efficiency for SRSSSlep-H(loose) signal region from Fig 18(e) in a SUSY scenario where $\tilde{q}$ decays into sleptons and subsequently to SM leptons and $\tilde{\chi}^{0}_{1}$

Number of signal events expected for 139 fb$^{-1}$ at different stages of the event selection for the signal region SRGGWZ-H in a susy scenario where $\tilde{g}$ decays into SM gauge bosons and $\tilde{\chi}^{0}_{1}$. The masses of the superpartners involved in the process are set to $m(\tilde{g})$ = 1400 GeV, $m(\tilde{\chi_{1}^{0}})$ = 1000 GeV. Only statistical uncertainties are shown.

Number of signal events expected for 139 fb$^{-1}$ at different stages of the event selection for the signal region SRGGWZ-M in a susy scenario where $\tilde{g}$ decays into SM gauge bosons and $\tilde{\chi}^{0}_{1}$. The masses of the superpartners involved in the process are set to $m(\tilde{g})$ = 1400 GeV, $m(\tilde{\chi_{1}^{0}})$ = 1000 GeV. Only statistical uncertainties are shown.

Number of signal events expected for 139 fb$^{-1}$ at different stages of the event selection for the signal region SRGGWZ-L in a susy scenario where $\tilde{g}$ decays into SM gauge bosons and $\tilde{\chi}^{0}_{1}$. The masses of the superpartners involved in the process are set to $m(\tilde{g})$ = 1400 GeV, $m(\tilde{\chi_{1}^{0}})$ = 1000 GeV. Only statistical uncertainties are shown.

Number of signal events expected for 139 fb$^{-1}$ at different stages of the event selection for the signal region SRGGSlep-L in a susy scenario where $\tilde{g}$ decays into sleptons and subsequently to SM leptons and $\tilde{\chi}^{0}_{1}$. The masses of the superpartners involved in the process are set to $m(\tilde{g})$ = 2000 GeV, $m(\tilde{\chi_{1}^{0}})$ = 500 GeV. Only statistical uncertainties are shown.

Number of signal events expected for 139 fb$^{-1}$ at different stages of the event selection for the signal region SRGGSlep-M in a susy scenario where $\tilde{g}$ decays into sleptons and subsequently to SM leptons and $\tilde{\chi}^{0}_{1}$. The masses of the superpartners involved in the process are set to $m(\tilde{g})$ = 2000 GeV, $m(\tilde{\chi_{1}^{0}})$ = 500 GeV. Only statistical uncertainties are shown.

Number of signal events expected for 139 fb$^{-1}$ at different stages of the event selection for the signal region SRGGSlep-H in a susy scenario where $\tilde{g}$ decays into sleptons and subsequently to SM leptons and $\tilde{\chi}^{0}_{1}$. The masses of the superpartners involved in the process are set to $m(\tilde{g})$ = 2000 GeV, $m(\tilde{\chi_{1}^{0}})$ = 500 GeV. Only statistical uncertainties are shown.

Number of signal events expected for 139 fb$^{-1}$ at different stages of the event selection for the signal region SRUDD-1b in a susy scenario where $\tilde{g}$ decays into anti-top and $\tilde{t}$ and $\tilde{t}$ decays via a non-zero RPV coupling $\lambda''$. The masses of the superpartners involved in the process are set to $m(\tilde{g})$ = 1600 GeV, $m(\tilde{t})$ = 600 GeV. Only statistical uncertainties are shown.

Number of signal events expected for 139 fb$^{-1}$ at different stages of the event selection for the signal region SRUDD-2b in a susy scenario where $\tilde{g}$ decays into anti-top and $\tilde{t}$ and $\tilde{t}$ decays via a non-zero RPV coupling $\lambda''$. The masses of the superpartners involved in the process are set to $m(\tilde{g})$ = 1600 GeV, $m(\tilde{t})$ = 600 GeV. Only statistical uncertainties are shown.

Number of signal events expected for 139 fb$^{-1}$ at different stages of the event selection for the signal region SRUDD-ge2b in a susy scenario where $\tilde{g}$ decays into anti-top and $\tilde{t}$ and $\tilde{t}$ decays via a non-zero RPV coupling $\lambda''$. The masses of the superpartners involved in the process are set to $m(\tilde{g})$ = 1600 GeV, $m(\tilde{t})$ = 600 GeV. Only statistical uncertainties are shown.

Number of signal events expected for 139 fb$^{-1}$ at different stages of the event selection for the signal region SRUDD-ge3b in a susy scenario where $\tilde{g}$ decays into anti-top and $\tilde{t}$ and $\tilde{t}$ decays via a non-zero RPV coupling $\lambda''$. The masses of the superpartners involved in the process are set to $m(\tilde{g})$ = 1600 GeV, $m(\tilde{t})$ = 600 GeV. Only statistical uncertainties are shown.

Number of signal events expected for 139 fb$^{-1}$ at different stages of the event selection for the signal region SRLQD in a susy scenario where direct $\tilde{\chi_{1}^{0}}$ decay into SM leptons and quarks via a non-zero RPV coupling $\lambda'$. The masses of the superpartners involved in the process are set to $m(\tilde{g})$ = 2200 GeV, $m(\tilde{\chi_{1}^{0}})$ = 1870 GeV. Only statistical uncertainties are shown.

Number of signal events expected for 139 fb$^{-1}$ at different stages of the event selection for the signal region SRSSWZ-L in a susy scenario where $\tilde{q}$ decays into SM gauge bosons and $\tilde{\chi}^{0}_{1}$. The masses of the superpartners involved in the process are set to $m(\tilde{q})$ = 800 GeV, $m(\tilde{\chi_{1}^{0}})$ = 600 GeV. Only statistical uncertainties are shown.

Number of signal events expected for 139 fb$^{-1}$ at different stages of the event selection for the signal region SRSSWZ-ML in a susy scenario where $\tilde{q}$ decays into SM gauge bosons and $\tilde{\chi}^{0}_{1}$. The masses of the superpartners involved in the process are set to $m(\tilde{q})$ = 800 GeV, $m(\tilde{\chi_{1}^{0}})$ = 600 GeV. Only statistical uncertainties are shown.

Number of signal events expected for 139 fb$^{-1}$ at different stages of the event selection for the signal region SRSSWZ-MH in a susy scenario where $\tilde{q}$ decays into SM gauge bosons and $\tilde{\chi}^{0}_{1}$. The masses of the superpartners involved in the process are set to $m(\tilde{q})$ = 800 GeV, $m(\tilde{\chi_{1}^{0}})$ = 600 GeV. Only statistical uncertainties are shown.

Number of signal events expected for 139 fb$^{-1}$ at different stages of the event selection for the signal region SRSSWZ-H in a susy scenario where $\tilde{q}$ decays into SM gauge bosons and $\tilde{\chi}^{0}_{1}$. The masses of the superpartners involved in the process are set to $m(\tilde{q})$ = 800 GeV, $m(\tilde{\chi_{1}^{0}})$ = 600 GeV. Only statistical uncertainties are shown.

Number of signal events expected for 139 fb$^{-1}$ at different stages of the event selection for the signal region SRSSSlep-H in a susy scenario where $\tilde{q}$ decays into sleptons and subsequently to SM leptons and $\tilde{\chi}^{0}_{1}$. The masses of the superpartners involved in the process are set to $m(\tilde{q})$ = 1000 GeV, $m(\tilde{\chi_{1}^{0}})$ = 800 GeV. Only statistical uncertainties are shown.

Number of signal events expected for 139 fb$^{-1}$ at different stages of the event selection for the signal region SRSSSlep-MH in a susy scenario where $\tilde{q}$ decays into sleptons and subsequently to SM leptons and $\tilde{\chi}^{0}_{1}$. The masses of the superpartners involved in the process are set to $m(\tilde{q})$ = 1000 GeV, $m(\tilde{\chi_{1}^{0}})$ = 800 GeV. Only statistical uncertainties are shown.

Number of signal events expected for 139 fb$^{-1}$ at different stages of the event selection for the signal region SRSSSlep-L in a susy scenario where $\tilde{q}$ decays into sleptons and subsequently to SM leptons and $\tilde{\chi}^{0}_{1}$. The masses of the superpartners involved in the process are set to $m(\tilde{q})$ = 1000 GeV, $m(\tilde{\chi_{1}^{0}})$ = 800 GeV. Only statistical uncertainties are shown.

Number of signal events expected for 139 fb$^{-1}$ at different stages of the event selection for the signal region SRSSSlep-ML in a susy scenario where $\tilde{q}$ decays into sleptons and subsequently to SM leptons and $\tilde{\chi}^{0}_{1}$. The masses of the superpartners involved in the process are set to $m(\tilde{q})$ = 1000 GeV, $m(\tilde{\chi_{1}^{0}})$ = 800 GeV. Only statistical uncertainties are shown.

Number of signal events expected for 139 fb$^{-1}$ at different stages of the event selection for the signal region SRSSSlep-H(loose) in a susy scenario where $\tilde{q}$ decays into sleptons and subsequently to SM leptons and $\tilde{\chi}^{0}_{1}$. The masses of the superpartners involved in the process are set to $m(\tilde{q})$ = 1000 GeV, $m(\tilde{\chi_{1}^{0}})$ = 800 GeV. Only statistical uncertainties are shown.

Cross-section upper limits at 95% CL from Fig1(a) for $\tilde{g}$ decays into SM gauge bosons and $\tilde{\chi}^{0}_{1}$

Cross-section upper limits at 95% CL from Fig1(c) for $\tilde{g}$ decays into sleptons and subsequently to SM leptons and $\tilde{\chi}^{0}_{1}$

Cross-section upper limits at 95% CL from Fig1(f) for $\tilde{g}$ decays into anti-top and $\tilde{t}$ and $\tilde{t}$ decays via a non-zero RPV coupling $\lambda''$

Cross-section upper limits at 95% CL from Fig1(e) for direct $\tilde{\chi_{1}^{0}}$ decay into SM leptons and quarks via a non-zero RPV coupling $\lambda'$

Cross-section upper limits at 95% CL from Fig1(b) for $\tilde{q}$ decays into SM gauge bosons and $\tilde{\chi}^{0}_{1}$

Cross-section upper limits at 95% CL from Fig1(d) for $\tilde{q}$ decays into sleptons and subsequently to SM leptons and $\tilde{\chi}^{0}_{1}$

Observed excluded cross-section at 95% CL for the WZ-mediated simplified model of wino $\tilde{\chi}^{\pm}_{1}/\tilde{\chi}^{0}_{2}$ production from Fig 8(aux).

Expected exclusion limits at 95% CL for the WZ-mediated simplified model of wino $\tilde{\chi}^{\pm}_{1}/\tilde{\chi}^{0}_{2}$ production.

Observed exclusion limits at 95% CL for the Wh-mediated simplified model of wino $\tilde{\chi}^{\pm}_{1}/\tilde{\chi}^{0}_{2}$ production from from Fig 13(a) and from Fig 7 and Fig 10(aux).

Observed excluded cross-section at 95% CL for the Wh-mediated simplified model of wino $\tilde{\chi}^{\pm}_{1}/\tilde{\chi}^{0}_{2}$ production from Fig 7(aux) and Fig 10(aux).

positive one $\sigma$ observed exclusion limits at 95% CL for the Wh-mediated simplified model of wino $\tilde{\chi}^{\pm}_{1}/\tilde{\chi}^{0}_{2}$ production from from Fig 13(a) and from Fig 7 and Fig 10(aux).

negative one $\sigma$ observed exclusion limits at 95% CL for the Wh-mediated simplified model of wino $\tilde{\chi}^{\pm}_{1}/\tilde{\chi}^{0}_{2}$ production from from Fig 13(a) and from Fig 7 and Fig 10(aux).

Expected exclusion limits at 95% CL for the Wh-mediated simplified model of wino $\tilde{\chi}^{\pm}_{1}/\tilde{\chi}^{0}_{2}$ production.

Expected exclusion limits at 95% CL for the Wh-mediated simplified model of wino $\tilde{\chi}^{\pm}_{1}/\tilde{\chi}^{0}_{2}$ production.

Expected exclusion limits at 95% CL for the Wh-mediated simplified model of wino $\tilde{\chi}^{\pm}_{1}/\tilde{\chi}^{0}_{2}$ production.

Number of signal events expected for 139 fb$^{-1}$ at different stages of the event selection for the signal region $SR^{bRPV}_{2l-SS}$. in a susy scenario where $\tilde{\chi}^{0} _{1} and \tilde{\chi}^{0} _{2}$ are produced in pairs and decay to all possible allowed bRPV decays. The masses of the superpartners involved in the process are set to $m(\tilde{\chi}^{0} _{1}/\tilde{\chi}^{0} _{2})$ = 200 GeV, tan$\beta$=5. Only statistical uncertainties are shown.

Number of signal events expected for 139 fb$^{-1}$ at different stages of the event selection for the signal region $SR^{bRPV}_{3l}$. in a susy scenario where $\tilde{\chi}^{0} _{1} and \tilde{\chi}^{0} _{2}$ are produced in pairs and decay to all possible allowed bRPV decays. The masses of the superpartners involved in the process are set to $m(\tilde{\chi}^{0} _{1}/\tilde{\chi}^{0} _{2})$ = 200 GeV, tan$\beta$=5. Only statistical uncertainties are shown.

Number of signal events expected for 139 fb$^{-1}$ at different stages of the event selection for the signal region $SR^{WZ}_{high-m_{T2}}$. The wino-like doublet pair ($\tilde{\chi}^{\pm} _{1} and \tilde{\chi}^{0} _{2}$) were produced and then decays into $bino-like \tilde{\chi}^{0} _{1}$ which is the lightest SUSY particle (LSP) accompanied by mass on-shell or mass off-shell W and Z bosons. The masses of the superpartners involved in the process are set to $m(\tilde{\chi}^{\pm} _{1}/\tilde{\chi}^{0} _{2})$ = 150 GeV, $m(\tilde{\chi}^{0} _{1})$ = 50 GeV. Only statistical uncertainties are shown.

Number of signal events expected for 139 fb$^{-1}$ at different stages of the event selection for the signal region $SR^{WZ}_{low-m_{T2}}$. The wino-like doublet pair ($\tilde{\chi}^{\pm} _{1} and \tilde{\chi}^{0} _{2}$) were produced and then decays into $bino-like \tilde{\chi}^{0} _{1}$ which is the lightest SUSY particle (LSP) accompanied by mass on-shell or mass off-shell W and Z bosons. The masses of the superpartners involved in the process are set to $m(\tilde{\chi}^{\pm} _{1}/\tilde{\chi}^{0} _{2})$ = 150 GeV, $m(\tilde{\chi}^{0} _{1})$ = 50 GeV. Only statistical uncertainties are shown.

Number of signal events expected for 139 fb$^{-1}$ at different stages of the event selection for the low mass $SR^{RPV}_{2l1b}$, where the $\tilde{\chi}^{0} _{1} and \tilde{\chi}^{0} _{2}$ are directly produced and undergoes prompt RPV decays. The masses of the superpartners involved in the process are set to $m(\tilde{\chi}^{0} _{1}/\tilde{\chi}^{0} _{2})$ = 200 GeV. Only statistical uncertainties are shown.

Number of signal events expected for 139 fb$^{-1}$ at different stages of the event selection for the medium mass $SR^{RPV}_{2l1b}$, where the $\tilde{\chi}^{0} _{1} and \tilde{\chi}^{0} _{2}$ are directly produced and undergoes prompt RPV decays. The masses of the superpartners involved in the process are set to $m(\tilde{\chi}^{0} _{1}/\tilde{\chi}^{0} _{2})$ = 200 GeV. Only statistical uncertainties are shown.

Number of signal events expected for 139 fb$^{-1}$ at different stages of the event selection for the low mass $SR^{RPV}_{2l2b}$, where the $\tilde{\chi}^{0} _{1} and \tilde{\chi}^{0} _{2}$ are directly produced and undergoes prompt RPV decays. The masses of the superpartners involved in the process are set to $m(\tilde{\chi}^{0} _{1}/\tilde{\chi}^{0} _{2})$ = 200 GeV. Only statistical uncertainties are shown.

Number of signal events expected for 139 fb$^{-1}$ at different stages of the event selection for the medium mass $SR^{RPV}_{2l2b}$, where the $\tilde{\chi}^{0} _{1} and \tilde{\chi}^{0} _{2}$ are directly produced and undergoes prompt RPV decays. The masses of the superpartners involved in the process are set to $m(\tilde{\chi}^{0} _{1}/\tilde{\chi}^{0} _{2})$ = 200 GeV. Only statistical uncertainties are shown.

Number of signal events expected for 139 fb$^{-1}$ at different stages of the event selection for the high mass $SR^{RPV}_{2l2b}$, where the $\tilde{\chi}^{0} _{1} and \tilde{\chi}^{0} _{2}$ are directly produced and undergoes prompt RPV decays. The masses of the superpartners involved in the process are set to $m(\tilde{\chi}^{0} _{1}/\tilde{\chi}^{0} _{2})$ = 200 GeV. Only statistical uncertainties are shown.

Number of signal events expected for 139 fb$^{-1}$ at different stages of the event selection for the low mass $SR^{RPV}_{2l3b}$, where the $\tilde{\chi}^{0} _{1} and \tilde{\chi}^{0} _{2}$ are directly produced and undergoes prompt RPV decays. The masses of the superpartners involved in the process are set to $m(\tilde{\chi}^{0} _{1}/\tilde{\chi}^{0} _{2})$ = 200 GeV. Only statistical uncertainties are shown.

Number of signal events expected for 139 fb$^{-1}$ at different stages of the event selection for the medium mass $SR^{RPV}_{2l3b}$, where the $\tilde{\chi}^{0} _{1} and \tilde{\chi}^{0} _{2}$ are directly produced and undergoes prompt RPV decays. The masses of the superpartners involved in the process are set to $m(\tilde{\chi}^{0} _{1}/\tilde{\chi}^{0} _{2})$ = 200 GeV. Only statistical uncertainties are shown.

Number of signal events expected for 139 fb$^{-1}$ at different stages of the event selection for the high mass $SR^{RPV}_{2l3b}$, where the $\tilde{\chi}^{0} _{1} and \tilde{\chi}^{0} _{2}$ are directly produced and undergoes prompt RPV decays. The masses of the superpartners involved in the process are set to $m(\tilde{\chi}^{0} _{1}/\tilde{\chi}^{0} _{2})$ = 200 GeV. Only statistical uncertainties are shown.

Number of signal events expected for 139 fb$^{-1}$ at different stages of the event selection for the $SR^{Wh}_{low-m_{T2} }$. The wino-like doublet pair ($\tilde{\chi}^{\pm} _{1} and \tilde{\chi}^{0} _{2}$) were produced and then decays into $bino-like \tilde{\chi}^{0} _{1}$ which is the lightest SUSY particle (LSP) accompanied by mass on-shell or mass off-shell W and Higgs bosons. The masses of the superpartners involved in the process are set to $m(\tilde{\chi}^{\pm} _{1}/\tilde{\chi}^{0} _{2})$ = 300 GeV, $m(\tilde{\chi}^{0} _{1})$ = 100 GeV. Only statistical uncertainties are shown.

Number of signal events expected for 139 fb$^{-1}$ at different stages of the event selection for the $SR^{Wh}_{high-m_{T2} }$. The wino-like doublet pair ($\tilde{\chi}^{\pm} _{1} and \tilde{\chi}^{0} _{2}$) were produced and then decays into $bino-like \tilde{\chi}^{0} _{1}$ which is the lightest SUSY particle (LSP) accompanied by mass on-shell or mass off-shell W and Higgs bosons. The masses of the superpartners involved in the process are set to $m(\tilde{\chi}^{\pm} _{1}/\tilde{\chi}^{0} _{2})$ = 300 GeV, $m(\tilde{\chi}^{0} _{1})$ = 100 GeV. Only statistical uncertainties are shown.

Signal Hepdataeptance for $SR^{bRPV}_{2l-SS}$ signal region from Fig 13(a)(aux) in a SUSY scenario where $\tilde{\chi}^{0} _{1} and \tilde{\chi}^{0} _{2}$ are produced in pairs and decay to all possible allowed bRPV decays.

Signal Hepdataeptance for $SR^{bRPV}_{3l}$ signal region from Fig 13(b)(aux) in a SUSY scenario where $\tilde{\chi}^{0} _{1} and \tilde{\chi}^{0} _{2}$ are produced in pairs and decay to all possible allowed bRPV decays.

Signal acceptance for $SR^{WZ}_{high-m_{T2}}$ in a SUSY scenario where the wino-like doublet pair ($\tilde{\chi}^{\pm} _{1} and \tilde{\chi}^{0} _{2}$) were produced and then decays into $bino-like \tilde{\chi}^{0} _{1}$ which is the lightest SUSY particle (LSP) accompanied by mass on-shell or mass off-shell W and Z bosons.

Signal acceptance for $SR^{WZ}_{low-m_{T2}}$ in a SUSY scenario where the wino-like doublet pair ($\tilde{\chi}^{\pm} _{1} and \tilde{\chi}^{0} _{2}$) were produced and then decays into $bino-like \tilde{\chi}^{0} _{1}$ which is the lightest SUSY particle (LSP) accompanied by mass on-shell or mass off-shell W and Z bosons.

Signal acceptance for $SR^{RPV}_{2l1b}-L$ signal region in a SUSY scenario where the $\tilde{\chi}^{0} _{1} and \tilde{\chi}^{0} _{2}$ are directly produced and undergoes prompt RPV decays.

Signal acceptance for $SR^{RPV}_{2l1b}-M$ signal region in a SUSY scenario where the $\tilde{\chi}^{0} _{1} and \tilde{\chi}^{0} _{2}$ are directly produced and undergoes prompt RPV decays.

Signal acceptance for $SR^{RPV}_{2l2b}-L$ signal region in a SUSY scenario where the $\tilde{\chi}^{0} _{1} and \tilde{\chi}^{0} _{2}$ are directly produced and undergoes prompt RPV decays.

Signal acceptance for $SR^{RPV}_{2l2b}-M$ signal region in a SUSY scenario where the $\tilde{\chi}^{0} _{1} and \tilde{\chi}^{0} _{2}$ are directly produced and undergoes prompt RPV decays.

Signal acceptance for $SR^{RPV}_{2l2b}-H$ signal region in a SUSY scenario where the $\tilde{\chi}^{0} _{1} and \tilde{\chi}^{0} _{2}$ are directly produced and undergoes prompt RPV decays.

Signal acceptance for $SR^{RPV}_{2l3b}-L$ signal region in a SUSY scenario where the $\tilde{\chi}^{0} _{1} and \tilde{\chi}^{0} _{2}$ are directly produced and undergoes prompt RPV decays.

Signal acceptance for $SR^{RPV}_{2l3b}-M$ signal region in a SUSY scenario where the $\tilde{\chi}^{0} _{1} and \tilde{\chi}^{0} _{2}$ are directly produced and undergoes prompt RPV decays.

Signal acceptance for $SR^{RPV}_{2l3b}-H$ signal region in a SUSY scenario where the $\tilde{\chi}^{0} _{1} and \tilde{\chi}^{0} _{2}$ are directly produced and undergoes prompt RPV decays.

Signal efficiency for $SR^{bRPV}_{2l-SS}$ signal region in a SUSY scenario where $\tilde{\chi}^{0} _{1} and \tilde{\chi}^{0} _{2}$ are produced in pairs and decay to all possible allowed bRPV decays.

Signal efficiency for $SR^{bRPV}_{3l}$ signal region in a SUSY scenario where $\tilde{\chi}^{0} _{1} and \tilde{\chi}^{0} _{2}$ are produced in pairs and decay to all possible allowed bRPV decays.

Signal efficiency for $SR^{WZ}_{high-m_{T2}}$ in a SUSY scenario where the wino-like doublet pair ($\tilde{\chi}^{\pm} _{1} and \tilde{\chi}^{0} _{2}$) were produced and then decays into $bino-like \tilde{\chi}^{0} _{1}$ which is the lightest SUSY particle (LSP) accompanied by mass on-shell or mass off-shell W and Z bosons.

Signal efficiency for $SR^{WZ}_{low-m_{T2}}$ in a SUSY scenario where the wino-like doublet pair ($\tilde{\chi}^{\pm} _{1} and \tilde{\chi}^{0} _{2}$) were produced and then decays into $bino-like \tilde{\chi}^{0} _{1}$ which is the lightest SUSY particle (LSP) accompanied by mass on-shell or mass off-shell W and Z bosons.

Signal efficiency for $SR^{RPV}_{2l1b}-L$ signal region in a SUSY scenario where the $\tilde{\chi}^{0} _{1} and \tilde{\chi}^{0} _{2}$ are directly produced and undergoes prompt RPV decays.

Signal efficiency for $SR^{RPV}_{2l1b}-M$ signal region in a SUSY scenario where the $\tilde{\chi}^{0} _{1} and \tilde{\chi}^{0} _{2}$ are directly produced and undergoes prompt RPV decays.

Signal efficiency for $SR^{RPV}_{2l2b}-L$ signal region in a SUSY scenario where the $\tilde{\chi}^{0} _{1} and \tilde{\chi}^{0} _{2}$ are directly produced and undergoes prompt RPV decays.

Signal efficiency for $SR^{RPV}_{2l2b}-M$ signal region in a SUSY scenario where the $\tilde{\chi}^{0} _{1} and \tilde{\chi}^{0} _{2}$ are directly produced and undergoes prompt RPV decays.

Signal efficiency for $SR^{RPV}_{2l2b}-H$ signal region in a SUSY scenario where the $\tilde{\chi}^{0} _{1} and \tilde{\chi}^{0} _{2}$ are directly produced and undergoes prompt RPV decays.

Signal efficiency for $SR^{RPV}_{2l3b}-L$ signal region in a SUSY scenario where the $\tilde{\chi}^{0} _{1} and \tilde{\chi}^{0} _{2}$ are directly produced and undergoes prompt RPV decays.

Signal efficiency for $SR^{RPV}_{2l3b}-M$ signal region in a SUSY scenario where the $\tilde{\chi}^{0} _{1} and \tilde{\chi}^{0} _{2}$ are directly produced and undergoes prompt RPV decays.

Signal efficiency for $SR^{RPV}_{2l3b}-H$ signal region in a SUSY scenario where the $\tilde{\chi}^{0} _{1} and \tilde{\chi}^{0} _{2}$ are directly produced and undergoes prompt RPV decays.

Signal acceptance for $SR^{Wh}_{high-m_{T2} }$ signal region from Fig 11(a)(aux) in a SUSY scenario where direct production of a lightest $\tilde{\chi}^{\pm} _{1} and \tilde{\chi}^{0} _{2}$ , decay with 100% branching ratio to a final state with a same sign light lepton (e or $\mu$) pair and two lightest neutralino1, via the on-shell emission of SM W and Higgs bosons,

Signal acceptance for $SR^{Wh}_{low-m_{T2} }$ signal region from Fig 11(b)(aux) in a SUSY scenario where direct production of a lightest $\tilde{\chi}^{\pm} _{1} and \tilde{\chi}^{0} _{2}$ , decay with 100% branching ratio to a final state with a same sign light lepton (e or $\mu$) pair and two lightest neutralino1, via the on-shell emission of SM W and Higgs bosons,

Signal efficiency for $SR^{Wh}_{high-m_{T2} }$ signal region from Fig 15(a)(aux) in a SUSY scenario where direct production of a lightest $\tilde{\chi}^{\pm} _{1} and \tilde{\chi}^{0} _{2}$ , decay with 100% branching ratio to a final state with a same sign light lepton (e or $\mu$) pair and two lightest neutralino1, via the on-shell emission of SM W and Higgs bosons,

Signal efficiency for $SR^{Wh}_{low-m_{T2} }$ signal region from Fig 15(b)(aux) in a SUSY scenario where direct production of a lightest $\tilde{\chi}^{\pm} _{1} and \tilde{\chi}^{0} _{2}$ , decay with 100% branching ratio to a final state with a same sign light lepton (e or $\mu$) pair and two lightest neutralino1, via the on-shell emission of SM W and Higgs bosons,

Observed 95% X-section upper limits as a function of higgsino $\tilde{\chi}^{\pm}_{1}/\tilde{\chi}^{0}_{1}/\tilde{\chi}^{0}_{2}$ mass in the bilinear RPV model from Fig 14.

Observed 95% X-section upper limits as a function of higgsino $\tilde{\chi}^{0}_{1}/\tilde{\chi}^{0}_{2}$ mass in the UDD RPV model from Fig 18.

Observed 95% X-section upper limits as a function of wino $\tilde{\chi}^{\pm}_{1}/\tilde{\chi}^{0}_{2}$ mass in the WZ-mediated simplified model of wino $\tilde{\chi}^{\pm}_{1}/\tilde{\chi}^{0}_{2}$ production from Fig 9(aux).

N-1 distributions for $m_{T2}$ of observed data and expected background towards $SR^{WZ}_{high-m_{T2}}$ from publication's Figure 11(a) . The last bin is inclusive.

N-1 distributions for $m_{T2}$ of observed data and expected background towards $SR^{WZ}_{low-m_{T2}}$ from publication's Figure 11(b) . The last bin is inclusive.

N-1 distributions for $m_{T2}$ of observed data and expected background towards $SR^{bRPV}_{2l-SS}$ from publication's Figure 11(c) . The last bin is inclusive.

N-1 distributions for $m_{T2}$ of observed data and expected background towards $SR^{bRPV}_{3l}$ from publication's Figure 11(d) . The last bin is inclusive.

N-1 distributions for $\sum p^{b-jet}_{T}/\sum p^{jet}_{T}$ of observed data and expected background towards $SR^{RPV}_{2l1b}-L$ from publication's Figure 16(a) . The last bin is inclusive.

N-1 distributions for $\sum p^{b-jet}_{T}/\sum p^{jet}_{T}$ of observed data and expected background towards $SR^{RPV}_{2l2b}-M$ from publication's Figure 16(b) . The last bin is inclusive.

N-1 distributions for $\sum p^{b-jet}_{T}/\sum p^{jet}_{T}$ of observed data and expected background towards $SR^{RPV}_{2l3b}-H$ from publication's Figure 16(c) . The last bin is inclusive.

N-1 distribution for $E_{T}^{miss}$ in $SR^{Wh}_{high-m_{T2} }$ in ee channel

N-1 distribution for $E_{T}^{miss}$ in $SR^{Wh}_{high-m_{T2} }$ in e$\mu$ channel

N-1 distribution for $E_{T}^{miss}$ in $SR^{Wh}_{high-m_{T2} }$ in $\mu\mu$ channel

N-1 distribution for $\mathcal{S}(E_{T}^{miss})$ in $SR^{Wh}_{low-m_{T2} }$ in ee channel

N-1 distribution for $\mathcal{S}(E_{T}^{miss})$ in $SR^{Wh}_{low-m_{T2} }$ in e$\mu$ channel

N-1 distribution for $\mathcal{S}(E_{T}^{miss})$ in $SR^{Wh}_{low-m_{T2} }$ in $\mu\mu$ channel

Average timing distributions for SR data and the estimated background as determined by the background-only fit, in each of the five exclusive $\rho$ categories. For comparison, the expected timing shapes for a few different signal models are superimposed, with each model labeled by the values of the $\tilde\chi^0_1$ mass and lifetime, as well as decay mode. To provide some indication of the variations in signal yield and shape, three signal models are shown for each of the $\tilde\chi^0_1$ decay modes, namely $\tilde\chi^0_1$ $\rightarrow$ $H \tilde G$ and $\tilde\chi^0_1$ $\rightarrow$ $Z \tilde G$. The models shown include a rather low $\tilde\chi^0_1$ mass value of 135 GeV for lifetimes of either 2 ns or 10 ns, and a higher $\tilde\chi^0_1$ mass value which is near the 95% CL exclusion limit for each decay mode for a lifetime of 2 ns. Each signal model is shown with the signal normalization corresponding to a BR value of unity for the decay mode in question.

Average timing distributions for SR data and the estimated background as determined by the background-only fit, in each of the five exclusive $\rho$ categories. For comparison, the expected timing shapes for a few different signal models are superimposed, with each model labeled by the values of the $\tilde\chi^0_1$ mass and lifetime, as well as decay mode. To provide some indication of the variations in signal yield and shape, three signal models are shown for each of the $\tilde\chi^0_1$ decay modes, namely $\tilde\chi^0_1$ $\rightarrow$ $H \tilde G$ and $\tilde\chi^0_1$ $\rightarrow$ $Z \tilde G$. The models shown include a rather low $\tilde\chi^0_1$ mass value of 135 GeV for lifetimes of either 2 ns or 10 ns, and a higher $\tilde\chi^0_1$ mass value which is near the 95% CL exclusion limit for each decay mode for a lifetime of 2 ns. Each signal model is shown with the signal normalization corresponding to a BR value of unity for the decay mode in question.

The 95% CL limits on $\sigma(pp \rightarrow \tilde\chi^0_1 \tilde\chi^0_1$) in fb as a function of $\tilde\chi^0_1$ mass (left) and $\tilde\chi^0_1$ lifetime (right), for the different decay modes of $\mathcal{B}$($\tilde\chi^0_1$ $\rightarrow H + \tilde{G}$) = 1 (top) and $\mathcal{B}$($\tilde\chi^0_1$ $\rightarrow Z +\tilde{G}$) = 1 (bottom). For the limits as a function of mass (lifetime), several signal models with varying lifetime (mass) are overlaid for comparison. Included are the theoretical expectations from higgsino production for each mass hypothesis, calculated from a GMSB SUSY model that assumes nearly degenerate $\tilde\chi^0_1$, $\tilde\chi^\pm_1$, and $\tilde\chi^0_2$.

The 95% CL limits on $\sigma(pp \rightarrow \tilde\chi^0_1 \tilde\chi^0_1$) in fb as a function of $\tilde\chi^0_1$ mass (left) and $\tilde\chi^0_1$ lifetime (right), for the different decay modes of $\mathcal{B}$($\tilde\chi^0_1$ $\rightarrow H + \tilde{G}$) = 1 (top) and $\mathcal{B}$($\tilde\chi^0_1$ $\rightarrow Z +\tilde{G}$) = 1 (bottom). For the limits as a function of mass (lifetime), several signal models with varying lifetime (mass) are overlaid for comparison. Included are the theoretical expectations from higgsino production for each mass hypothesis, calculated from a GMSB SUSY model that assumes nearly degenerate $\tilde\chi^0_1$, $\tilde\chi^\pm_1$, and $\tilde\chi^0_2$.

The 95% CL limits on $\sigma(pp \rightarrow \tilde\chi^0_1 \tilde\chi^0_1$) in fb as a function of $\tilde\chi^0_1$ mass (left) and $\tilde\chi^0_1$ lifetime (right), for the different decay modes of $\mathcal{B}$($\tilde\chi^0_1$ $\rightarrow H + \tilde{G}$) = 1 (top) and $\mathcal{B}$($\tilde\chi^0_1$ $\rightarrow Z +\tilde{G}$) = 1 (bottom). For the limits as a function of mass (lifetime), several signal models with varying lifetime (mass) are overlaid for comparison. Included are the theoretical expectations from higgsino production for each mass hypothesis, calculated from a GMSB SUSY model that assumes nearly degenerate $\tilde\chi^0_1$, $\tilde\chi^\pm_1$, and $\tilde\chi^0_2$.

The 95% CL limits on $\sigma(pp \rightarrow \tilde\chi^0_1 \tilde\chi^0_1$) in fb as a function of $\tilde\chi^0_1$ mass (left) and $\tilde\chi^0_1$ lifetime (right), for the different decay modes of $\mathcal{B}$($\tilde\chi^0_1$ $\rightarrow H + \tilde{G}$) = 1 (top) and $\mathcal{B}$($\tilde\chi^0_1$ $\rightarrow Z +\tilde{G}$) = 1 (bottom). For the limits as a function of mass (lifetime), several signal models with varying lifetime (mass) are overlaid for comparison. Included are the theoretical expectations from higgsino production for each mass hypothesis, calculated from a GMSB SUSY model that assumes nearly degenerate $\tilde\chi^0_1$, $\tilde\chi^\pm_1$, and $\tilde\chi^0_2$.

The 95% CL limits on $\sigma(pp \rightarrow \tilde\chi^0_1 \tilde\chi^0_1$) in fb as a function of $\tilde\chi^0_1$ mass (left) and $\tilde\chi^0_1$ lifetime (right), for the different decay modes of $\mathcal{B}$($\tilde\chi^0_1$ $\rightarrow H + \tilde{G}$) = 1 (top) and $\mathcal{B}$($\tilde\chi^0_1$ $\rightarrow Z +\tilde{G}$) = 1 (bottom). For the limits as a function of mass (lifetime), several signal models with varying lifetime (mass) are overlaid for comparison. Included are the theoretical expectations from higgsino production for each mass hypothesis, calculated from a GMSB SUSY model that assumes nearly degenerate $\tilde\chi^0_1$, $\tilde\chi^\pm_1$, and $\tilde\chi^0_2$.

The 95% CL limits on $\sigma(pp \rightarrow \tilde\chi^0_1 \tilde\chi^0_1$) in fb as a function of $\tilde\chi^0_1$ mass (left) and $\tilde\chi^0_1$ lifetime (right), for the different decay modes of $\mathcal{B}$($\tilde\chi^0_1$ $\rightarrow H + \tilde{G}$) = 1 (top) and $\mathcal{B}$($\tilde\chi^0_1$ $\rightarrow Z +\tilde{G}$) = 1 (bottom). For the limits as a function of mass (lifetime), several signal models with varying lifetime (mass) are overlaid for comparison. Included are the theoretical expectations from higgsino production for each mass hypothesis, calculated from a GMSB SUSY model that assumes nearly degenerate $\tilde\chi^0_1$, $\tilde\chi^\pm_1$, and $\tilde\chi^0_2$.

The 95% CL limits on $\sigma(pp \rightarrow \tilde\chi^0_1 \tilde\chi^0_1$) in fb as a function of $\tilde\chi^0_1$ mass (left) and $\tilde\chi^0_1$ lifetime (right), for the different decay modes of $\mathcal{B}$($\tilde\chi^0_1$ $\rightarrow H + \tilde{G}$) = 1 (top) and $\mathcal{B}$($\tilde\chi^0_1$ $\rightarrow Z +\tilde{G}$) = 1 (bottom). For the limits as a function of mass (lifetime), several signal models with varying lifetime (mass) are overlaid for comparison. Included are the theoretical expectations from higgsino production for each mass hypothesis, calculated from a GMSB SUSY model that assumes nearly degenerate $\tilde\chi^0_1$, $\tilde\chi^\pm_1$, and $\tilde\chi^0_2$.

The 95% CL limits on $\sigma(pp \rightarrow \tilde\chi^0_1 \tilde\chi^0_1$) in fb as a function of $\tilde\chi^0_1$ mass (left) and $\tilde\chi^0_1$ lifetime (right), for the different decay modes of $\mathcal{B}$($\tilde\chi^0_1$ $\rightarrow H + \tilde{G}$) = 1 (top) and $\mathcal{B}$($\tilde\chi^0_1$ $\rightarrow Z +\tilde{G}$) = 1 (bottom). For the limits as a function of mass (lifetime), several signal models with varying lifetime (mass) are overlaid for comparison. Included are the theoretical expectations from higgsino production for each mass hypothesis, calculated from a GMSB SUSY model that assumes nearly degenerate $\tilde\chi^0_1$, $\tilde\chi^\pm_1$, and $\tilde\chi^0_2$.

The 95% CL limits on $\sigma(pp \rightarrow \tilde\chi^0_1 \tilde\chi^0_1$) in fb as a function of $\tilde\chi^0_1$ mass (left) and $\tilde\chi^0_1$ lifetime (right), for the different decay modes of $\mathcal{B}$($\tilde\chi^0_1$ $\rightarrow H + \tilde{G}$) = 1 (top) and $\mathcal{B}$($\tilde\chi^0_1$ $\rightarrow Z +\tilde{G}$) = 1 (bottom). For the limits as a function of mass (lifetime), several signal models with varying lifetime (mass) are overlaid for comparison. Included are the theoretical expectations from higgsino production for each mass hypothesis, calculated from a GMSB SUSY model that assumes nearly degenerate $\tilde\chi^0_1$, $\tilde\chi^\pm_1$, and $\tilde\chi^0_2$.

The 95% CL limits on $\sigma(pp \rightarrow \tilde\chi^0_1 \tilde\chi^0_1$) in fb as a function of $\tilde\chi^0_1$ mass (left) and $\tilde\chi^0_1$ lifetime (right), for the different decay modes of $\mathcal{B}$($\tilde\chi^0_1$ $\rightarrow H + \tilde{G}$) = 1 (top) and $\mathcal{B}$($\tilde\chi^0_1$ $\rightarrow Z +\tilde{G}$) = 1 (bottom). For the limits as a function of mass (lifetime), several signal models with varying lifetime (mass) are overlaid for comparison. Included are the theoretical expectations from higgsino production for each mass hypothesis, calculated from a GMSB SUSY model that assumes nearly degenerate $\tilde\chi^0_1$, $\tilde\chi^\pm_1$, and $\tilde\chi^0_2$.

The 95% CL limits on $\sigma(pp \rightarrow \tilde\chi^0_1 \tilde\chi^0_1$) in fb as a function of $\tilde\chi^0_1$ mass (left) and $\tilde\chi^0_1$ lifetime (right), for the different decay modes of $\mathcal{B}$($\tilde\chi^0_1$ $\rightarrow H + \tilde{G}$) = 1 (top) and $\mathcal{B}$($\tilde\chi^0_1$ $\rightarrow Z +\tilde{G}$) = 1 (bottom). For the limits as a function of mass (lifetime), several signal models with varying lifetime (mass) are overlaid for comparison. Included are the theoretical expectations from higgsino production for each mass hypothesis, calculated from a GMSB SUSY model that assumes nearly degenerate $\tilde\chi^0_1$, $\tilde\chi^\pm_1$, and $\tilde\chi^0_2$.

The 95% CL limits on $\sigma(pp \rightarrow \tilde\chi^0_1 \tilde\chi^0_1$) in fb as a function of $\tilde\chi^0_1$ mass (left) and $\tilde\chi^0_1$ lifetime (right), for the different decay modes of $\mathcal{B}$($\tilde\chi^0_1$ $\rightarrow H + \tilde{G}$) = 1 (top) and $\mathcal{B}$($\tilde\chi^0_1$ $\rightarrow Z +\tilde{G}$) = 1 (bottom). For the limits as a function of mass (lifetime), several signal models with varying lifetime (mass) are overlaid for comparison. Included are the theoretical expectations from higgsino production for each mass hypothesis, calculated from a GMSB SUSY model that assumes nearly degenerate $\tilde\chi^0_1$, $\tilde\chi^\pm_1$, and $\tilde\chi^0_2$.

The 95% CL limits on $\sigma(pp \rightarrow \tilde\chi^0_1 \tilde\chi^0_1$) in fb as a function of $\tilde\chi^0_1$ mass (left) and $\tilde\chi^0_1$ lifetime (right), for the different decay modes of $\mathcal{B}$($\tilde\chi^0_1$ $\rightarrow H + \tilde{G}$) = 1 (top) and $\mathcal{B}$($\tilde\chi^0_1$ $\rightarrow Z +\tilde{G}$) = 1 (bottom). For the limits as a function of mass (lifetime), several signal models with varying lifetime (mass) are overlaid for comparison. Included are the theoretical expectations from higgsino production for each mass hypothesis, calculated from a GMSB SUSY model that assumes nearly degenerate $\tilde\chi^0_1$, $\tilde\chi^\pm_1$, and $\tilde\chi^0_2$.

The 95% CL limits on $\sigma(pp \rightarrow \tilde\chi^0_1 \tilde\chi^0_1$) in fb as a function of $\tilde\chi^0_1$ mass (left) and $\tilde\chi^0_1$ lifetime (right), for the different decay modes of $\mathcal{B}$($\tilde\chi^0_1$ $\rightarrow H + \tilde{G}$) = 1 (top) and $\mathcal{B}$($\tilde\chi^0_1$ $\rightarrow Z +\tilde{G}$) = 1 (bottom). For the limits as a function of mass (lifetime), several signal models with varying lifetime (mass) are overlaid for comparison. Included are the theoretical expectations from higgsino production for each mass hypothesis, calculated from a GMSB SUSY model that assumes nearly degenerate $\tilde\chi^0_1$, $\tilde\chi^\pm_1$, and $\tilde\chi^0_2$.

The 95% CL limits on $\sigma(pp \rightarrow \tilde\chi^0_1 \tilde\chi^0_1$) in fb as a function of $\tilde\chi^0_1$ mass (left) and $\tilde\chi^0_1$ lifetime (right), for the different decay modes of $\mathcal{B}$($\tilde\chi^0_1$ $\rightarrow H + \tilde{G}$) = 1 (top) and $\mathcal{B}$($\tilde\chi^0_1$ $\rightarrow Z +\tilde{G}$) = 1 (bottom). For the limits as a function of mass (lifetime), several signal models with varying lifetime (mass) are overlaid for comparison. Included are the theoretical expectations from higgsino production for each mass hypothesis, calculated from a GMSB SUSY model that assumes nearly degenerate $\tilde\chi^0_1$, $\tilde\chi^\pm_1$, and $\tilde\chi^0_2$.

The 95% CL limits on $\sigma(pp \rightarrow \tilde\chi^0_1 \tilde\chi^0_1$) in fb as a function of $\tilde\chi^0_1$ mass (left) and $\tilde\chi^0_1$ lifetime (right), for the different decay modes of $\mathcal{B}$($\tilde\chi^0_1$ $\rightarrow H + \tilde{G}$) = 1 (top) and $\mathcal{B}$($\tilde\chi^0_1$ $\rightarrow Z +\tilde{G}$) = 1 (bottom). For the limits as a function of mass (lifetime), several signal models with varying lifetime (mass) are overlaid for comparison. Included are the theoretical expectations from higgsino production for each mass hypothesis, calculated from a GMSB SUSY model that assumes nearly degenerate $\tilde\chi^0_1$, $\tilde\chi^\pm_1$, and $\tilde\chi^0_2$.

The 95% CL limits on $\sigma(pp \rightarrow \tilde\chi^0_1 \tilde\chi^0_1$) in fb as a function of $\tilde\chi^0_1$ branching ratio to the SM Higgs boson, where the assumed cross-section is for higgsino production, and $\mathcal{B}$($\tilde\chi^0_1$ $\rightarrow Z +\tilde{G}$) = 1 - $\mathcal{B}$($\tilde\chi^0_1$ $\rightarrow H + \tilde{G}$). Several signal hypotheses are overlaid that are labelled by the $\tilde\chi^0_1$ mass, all with a fixed $\tilde\chi^0_1$ lifetime of 2 ns.

The 95% CL limits on $\sigma(pp \rightarrow \tilde\chi^0_1 \tilde\chi^0_1$) in fb as a function of $\tilde\chi^0_1$ branching ratio to the SM Higgs boson, where the assumed cross-section is for higgsino production, and $\mathcal{B}$($\tilde\chi^0_1$ $\rightarrow Z +\tilde{G}$) = 1 - $\mathcal{B}$($\tilde\chi^0_1$ $\rightarrow H + \tilde{G}$). Several signal hypotheses are overlaid that are labelled by the $\tilde\chi^0_1$ mass, all with a fixed $\tilde\chi^0_1$ lifetime of 2 ns.

The 95% CL limits on $\sigma(pp \rightarrow \tilde\chi^0_1 \tilde\chi^0_1$) in fb as a function of $\tilde\chi^0_1$ branching ratio to the SM Higgs boson, where the assumed cross-section is for higgsino production, and $\mathcal{B}$($\tilde\chi^0_1$ $\rightarrow Z +\tilde{G}$) = 1 - $\mathcal{B}$($\tilde\chi^0_1$ $\rightarrow H + \tilde{G}$). Several signal hypotheses are overlaid that are labelled by the $\tilde\chi^0_1$ mass, all with a fixed $\tilde\chi^0_1$ lifetime of 2 ns.

The 95% CL limits on $\sigma(pp \rightarrow \tilde\chi^0_1 \tilde\chi^0_1$) in fb as a function of $\tilde\chi^0_1$ branching ratio to the SM Higgs boson, where the assumed cross-section is for higgsino production, and $\mathcal{B}$($\tilde\chi^0_1$ $\rightarrow Z +\tilde{G}$) = 1 - $\mathcal{B}$($\tilde\chi^0_1$ $\rightarrow H + \tilde{G}$). Several signal hypotheses are overlaid that are labelled by the $\tilde\chi^0_1$ mass, all with a fixed $\tilde\chi^0_1$ lifetime of 2 ns.

The 95% CL exclusion limits on the target signal hypothesis, for $\tilde\chi^0_1$ lifetime in ns as a function of $\tilde\chi^0_1$ mass in GeV. The overlaid curves correspond to different decay hypotheses, where the assumed cross-section is for higgsino production, and the $\tilde\chi^0_1$ decays to $H + \tilde{G}$ or $Z + \tilde{G}$ such that $\mathcal{B}(H + \tilde{G}) + \mathcal{B}(Z + \tilde{G})$ = 100%. The curve shown in red represents the decay hypothesis where the $\tilde\chi^0_1$ decays to $Z + \tilde{G}$ with 100% branching ratio. The curve shown in blue represents the decay hypothesis where the $\tilde\chi^0_1$ decays to $H + \tilde{G}$ with 100% branching ratio.

The 95% CL exclusion limits on the target signal hypothesis, for $\tilde\chi^0_1$ lifetime in ns as a function of $\tilde\chi^0_1$ mass in GeV. The overlaid curves correspond to different decay hypotheses, where the assumed cross-section is for higgsino production, and the $\tilde\chi^0_1$ decays to $H + \tilde{G}$ or $Z + \tilde{G}$ such that $\mathcal{B}(H + \tilde{G}) + \mathcal{B}(Z + \tilde{G})$ = 100%. The curve shown in red represents the decay hypothesis where the $\tilde\chi^0_1$ decays to $Z + \tilde{G}$ with 100% branching ratio. The curve shown in blue represents the decay hypothesis where the $\tilde\chi^0_1$ decays to $H + \tilde{G}$ with 100% branching ratio.

The 95% CL exclusion limits on the target signal hypothesis, for $\tilde\chi^0_1$ lifetime in ns as a function of $\tilde\chi^0_1$ mass in GeV. The overlaid curves correspond to different decay hypotheses, where the assumed cross-section is for higgsino production, and the $\tilde\chi^0_1$ decays to $H + \tilde{G}$ or $Z + \tilde{G}$ such that $\mathcal{B}(H + \tilde{G}) + \mathcal{B}(Z + \tilde{G})$ = 100%. The curve shown in red represents the decay hypothesis where the $\tilde\chi^0_1$ decays to $Z + \tilde{G}$ with 100% branching ratio. The curve shown in blue represents the decay hypothesis where the $\tilde\chi^0_1$ decays to $H + \tilde{G}$ with 100% branching ratio.

The 95% CL exclusion limits on the target signal hypothesis, for $\tilde\chi^0_1$ lifetime in ns as a function of $\tilde\chi^0_1$ mass in GeV. The overlaid curves correspond to different decay hypotheses, where the assumed cross-section is for higgsino production, and the $\tilde\chi^0_1$ decays to $H + \tilde{G}$ or $Z + \tilde{G}$ such that $\mathcal{B}(H + \tilde{G}) + \mathcal{B}(Z + \tilde{G})$ = 100%. The curve shown in red represents the decay hypothesis where the $\tilde\chi^0_1$ decays to $Z + \tilde{G}$ with 100% branching ratio. The curve shown in blue represents the decay hypothesis where the $\tilde\chi^0_1$ decays to $H + \tilde{G}$ with 100% branching ratio.

The 95% CL exclusion limits on the target signal hypothesis, for $\tilde\chi^0_1$ lifetime in ns as a function of $\tilde\chi^0_1$ mass in GeV. The overlaid curves correspond to different decay hypotheses, where the assumed cross-section is for higgsino production, and the $\tilde\chi^0_1$ decays to $H + \tilde{G}$ or $Z + \tilde{G}$ such that $\mathcal{B}(H + \tilde{G}) + \mathcal{B}(Z + \tilde{G})$ = 100%. The curve shown in red represents the decay hypothesis where the $\tilde\chi^0_1$ decays to $Z + \tilde{G}$ with 100% branching ratio. The curve shown in blue represents the decay hypothesis where the $\tilde\chi^0_1$ decays to $H + \tilde{G}$ with 100% branching ratio.

The 95% CL exclusion limits on the target signal hypothesis, for $\tilde\chi^0_1$ lifetime in ns as a function of $\tilde\chi^0_1$ mass in GeV. The overlaid curves correspond to different decay hypotheses, where the assumed cross-section is for higgsino production, and the $\tilde\chi^0_1$ decays to $H + \tilde{G}$ or $Z + \tilde{G}$ such that $\mathcal{B}(H + \tilde{G}) + \mathcal{B}(Z + \tilde{G})$ = 100%. The curve shown in red represents the decay hypothesis where the $\tilde\chi^0_1$ decays to $Z + \tilde{G}$ with 100% branching ratio. The curve shown in blue represents the decay hypothesis where the $\tilde\chi^0_1$ decays to $H + \tilde{G}$ with 100% branching ratio.

The 95% CL exclusion limits on the target signal hypothesis, for $\tilde\chi^0_1$ lifetime in ns as a function of $\tilde\chi^0_1$ mass in GeV. The overlaid curves correspond to different decay hypotheses, where the assumed cross-section is for higgsino production, and the $\tilde\chi^0_1$ decays to $H + \tilde{G}$ or $Z + \tilde{G}$ such that $\mathcal{B}(H + \tilde{G}) + \mathcal{B}(Z + \tilde{G})$ = 100%. The curve shown in red represents the decay hypothesis where the $\tilde\chi^0_1$ decays to $Z + \tilde{G}$ with 100% branching ratio. The curve shown in blue represents the decay hypothesis where the $\tilde\chi^0_1$ decays to $H + \tilde{G}$ with 100% branching ratio.

The 95% CL exclusion limits on the target signal hypothesis, for $\tilde\chi^0_1$ lifetime in ns as a function of $\tilde\chi^0_1$ mass in GeV. The overlaid curves correspond to different decay hypotheses, where the assumed cross-section is for higgsino production, and the $\tilde\chi^0_1$ decays to $H + \tilde{G}$ or $Z + \tilde{G}$ such that $\mathcal{B}(H + \tilde{G}) + \mathcal{B}(Z + \tilde{G})$ = 100%. The curve shown in red represents the decay hypothesis where the $\tilde\chi^0_1$ decays to $Z + \tilde{G}$ with 100% branching ratio. The curve shown in blue represents the decay hypothesis where the $\tilde\chi^0_1$ decays to $H + \tilde{G}$ with 100% branching ratio.

The 95% CL exclusion limits on the target signal hypothesis, for $\tilde\chi^0_1$ lifetime in ns as a function of $\tilde\chi^0_1$ mass in GeV. The overlaid curves correspond to different decay hypotheses, where the assumed cross-section is for higgsino production, and the $\tilde\chi^0_1$ decays to $H + \tilde{G}$ or $Z + \tilde{G}$ such that $\mathcal{B}(H + \tilde{G}) + \mathcal{B}(Z + \tilde{G})$ = 100%. The curve shown in red represents the decay hypothesis where the $\tilde\chi^0_1$ decays to $Z + \tilde{G}$ with 100% branching ratio. The curve shown in blue represents the decay hypothesis where the $\tilde\chi^0_1$ decays to $H + \tilde{G}$ with 100% branching ratio.

The 95% CL exclusion limits on the target signal hypothesis, for $\tilde\chi^0_1$ lifetime in ns as a function of $\tilde\chi^0_1$ mass in GeV. The overlaid curves correspond to different decay hypotheses, where the assumed cross-section is for higgsino production, and the $\tilde\chi^0_1$ decays to $H + \tilde{G}$ or $Z + \tilde{G}$ such that $\mathcal{B}(H + \tilde{G}) + \mathcal{B}(Z + \tilde{G})$ = 100%. The curve shown in red represents the decay hypothesis where the $\tilde\chi^0_1$ decays to $Z + \tilde{G}$ with 100% branching ratio. The curve shown in blue represents the decay hypothesis where the $\tilde\chi^0_1$ decays to $H + \tilde{G}$ with 100% branching ratio.

The 95% CL exclusion limits on the target signal hypothesis, for $\tilde\chi^0_1$ lifetime in ns as a function of $\tilde\chi^0_1$ mass in GeV. The overlaid curves correspond to different decay hypotheses, where the assumed cross-section is for higgsino production, and the $\tilde\chi^0_1$ decays to $H + \tilde{G}$ or $Z + \tilde{G}$ such that $\mathcal{B}(H + \tilde{G}) + \mathcal{B}(Z + \tilde{G})$ = 100%. The curve shown in red represents the decay hypothesis where the $\tilde\chi^0_1$ decays to $Z + \tilde{G}$ with 100% branching ratio. The curve shown in blue represents the decay hypothesis where the $\tilde\chi^0_1$ decays to $H + \tilde{G}$ with 100% branching ratio.

The 95% CL exclusion limits on the target signal hypothesis, for $\tilde\chi^0_1$ lifetime in ns as a function of $\tilde\chi^0_1$ mass in GeV. The overlaid curves correspond to different decay hypotheses, where the assumed cross-section is for higgsino production, and the $\tilde\chi^0_1$ decays to $H + \tilde{G}$ or $Z + \tilde{G}$ such that $\mathcal{B}(H + \tilde{G}) + \mathcal{B}(Z + \tilde{G})$ = 100%. The curve shown in red represents the decay hypothesis where the $\tilde\chi^0_1$ decays to $Z + \tilde{G}$ with 100% branching ratio. The curve shown in blue represents the decay hypothesis where the $\tilde\chi^0_1$ decays to $H + \tilde{G}$ with 100% branching ratio.

Cutflow for an example higgsino signal with mass 225 GeV and lifetime 10 ns, in the H decay mode. Acceptance is defined at truth level, and efficiency compares the events passing at reconstruction level with respect to truth.

Cutflow for an example higgsino signal with mass 225 GeV and lifetime 10 ns, in the Z decay mode. Acceptance is defined at truth level, and efficiency compares the events passing at reconstruction level with respect to truth.

Cutflow for an example higgsino signal with mass 225 GeV and lifetime 2 ns, in the H decay mode. Acceptance is defined at truth level, and efficiency compares the events passing at reconstruction level with respect to truth.

Cutflow for an example higgsino signal with mass 225 GeV and lifetime 2 ns, in the Z decay mode. Acceptance is defined at truth level, and efficiency compares the events passing at reconstruction level with respect to truth.

Acceptance across the H decay mode signal grid, calculated using truth information. The selection applied corresponds to the model-independent signal region (i.e. the standard SR with $t_{\text{avg}$ > 0.9 ns).

Acceptance across the Z decay mode signal grid, calculated using truth information. The selection applied corresponds to the model-independent signal region (i.e. the standard SR with $t_{\text{avg}$ > 0.9 ns).

Efficiency across the H decay mode signal grid, calculated using reco information. The selection applied corresponds to the model-independent signal region (i.e. the standard SR with $t_{\text{avg}$ > 0.9 ns). Here, the numerator is the signal yield passing the reco selection and the denominator is the signal yield passing the truth selection.

Efficiency across the Z decay mode signal grid, calculated using reco information. The selection applied corresponds to the model-independent signal region (i.e. the standard SR with $t_{\text{avg}$ > 0.9 ns). Here, the numerator is the signal yield passing the reco selection and the denominator is the signal yield passing the truth selection.

Two-dimensional distribution of the invariant mass $m_{DV}$ and the track multiplicity in the High-pT jet SR for observed data events

Two-dimensional distribution of the invariant mass $m_{DV}$ and the track multiplicity in the High-pT jet SR for expected signal events in the strong gluino pair pair production model with m(gluino)=1.8 TeV, m(chi0)=0.2 TeV, tau(chi0)=0.1 ns

Two-dimensional distribution of the invariant mass $m_{DV}$ and the track multiplicity in the Trackless jet SR for observed data events

Two-dimensional distribution of the invariant mass $m_{DV}$ and the track multiplicity in the Trackless jet SR for expected signal events in the electroweak pair production model

Expected exclusion limits at 95% CL on the lifetime and mass of the neutralino in electroweakino pair production models

Expected (+1 sigma) exclusion limits at 95% CL on the lifetime and mass of the neutralino in electroweakino pair production models

Expected (-1 sigma) exclusion limits at 95% CL on the lifetime and mass of the neutralino in electroweakino pair production models

Observed exclusion limits at 95% CL on the lifetime and mass of the neutralino in electroweakino pair production models

Observed (+1 sigma) exclusion limits at 95% CL on the lifetime and mass of the neutralino in electroweakino pair production models

Observed (-1 sigma) exclusion limits at 95% CL on the lifetime and mass of the neutralino in electroweakino pair production models

Expected exclusion limits at 95% CL on the lifetime and mass of the neutralino in strong gluino pair production models and m(gluino)=2.4 TeV

Expected (+1 sigma) exclusion limits at 95% CL on the lifetime and mass of the neutralino in strong gluino pair production models and m(gluino)=2.4 TeV

Expected (-1 sigma) exclusion limits at 95% CL on the lifetime and mass of the neutralino in strong gluino pair production models and m(gluino)=2.4 TeV

Observed exclusion limits at 95% CL on the lifetime and mass of the neutralino in strong gluino pair production models and m(gluino)=2.4 TeV

Observed (+1 sigma) exclusion limits at 95% CL on the lifetime and mass of the neutralino in strong gluino pair production models and m(gluino)=2.4 TeV

Observed (-1 sigma) exclusion limits at 95% CL on the lifetime and mass of the neutralino in strong gluino pair production models and m(gluino)=2.4 TeV

Exclusion limits at 95% CL on the production cross section in the electroweak pair production model.

Exclusion limits at 95% CL on the production cross section in the strong gluino pair production models and m(gluino)=2.4 TeV

Expected exclusion limits at 95% CL on the lifetime and mass of the neutralino in strong gluino pair production models and m(gluino)=2.0 TeV

Expected (+1 sigma) exclusion limits at 95% CL on the lifetime and mass of the neutralino in strong gluino pair production models and m(gluino)=2.0 TeV

Expected (-1 sigma) exclusion limits at 95% CL on the lifetime and mass of the neutralino in strong gluino pair production models and m(gluino)=2.0 TeV

Observed exclusion limits at 95% CL on the lifetime and mass of the neutralino in strong gluino pair production models and m(gluino)=2.0 TeV

Observed (+1 sigma) exclusion limits at 95% CL on the lifetime and mass of the neutralino in strong gluino pair production models and m(gluino)=2.0 TeV

Observed (-1 sigma) exclusion limits at 95% CL on the lifetime and mass of the neutralino in strong gluino pair production models and m(gluino)=2.0 TeV

Expected exclusion limits at 95% CL on the lifetime and mass of the neutralino in strong gluino pair production models and m(gluino)=2.2 TeV

Expected (+1 sigma) exclusion limits at 95% CL on the lifetime and mass of the neutralino in strong gluino pair production models and m(gluino)=2.2 TeV

Expected (-1 sigma) exclusion limits at 95% CL on the lifetime and mass of the neutralino in strong gluino pair production models and m(gluino)=2.2 TeV

Observed exclusion limits at 95% CL on the lifetime and mass of the neutralino in strong gluino pair production models and m(gluino)=2.2 TeV

Observed (+1 sigma) exclusion limits at 95% CL on the lifetime and mass of the neutralino in strong gluino pair production models and m(gluino)=2.2 TeV

Observed (-1 sigma) exclusion limits at 95% CL on the lifetime and mass of the neutralino in strong gluino pair production models and m(gluino)=2.2 TeV

Expected exclusion limits at 95% CL on the lifetime and mass of the gluino in strong gluino pair production models and m(chi0)=50 GeV

Expected (+1 sigma) exclusion limits at 95% CL on the lifetime and mass of the gluino in strong gluino pair production models and m(chi0)=50 GeV

Expected (-1 sigma) exclusion limits at 95% CL on the lifetime and mass of the gluino in strong gluino pair production models and m(chi0)=50 GeV

Observed exclusion limits at 95% CL on the lifetime and mass of the gluino in strong gluino pair production models and m(chi0)=50 GeV

Observed (+1 sigma) exclusion limits at 95% CL on the lifetime and mass of the gluino in strong gluino pair production models and m(chi0)=50 GeV

Observed (-1 sigma) exclusion limits at 95% CL on the lifetime and mass of the gluino in strong gluino pair production models and m(chi0)=50 GeV

Expected exclusion limits at 95% CL on the lifetime and mass of the gluino in strong gluino pair production models and m(chi0)=450 GeV

Expected (+1 sigma) exclusion limits at 95% CL on the lifetime and mass of the gluino in strong gluino pair production models and m(chi0)=450 GeV

Expected (-1 sigma) exclusion limits at 95% CL on the lifetime and mass of the gluino in strong gluino pair production models and m(chi0)=450 GeV

Observed exclusion limits at 95% CL on the lifetime and mass of the gluino in strong gluino pair production models and m(chi0)=450 GeV

Observed (+1 sigma) exclusion limits at 95% CL on the lifetime and mass of the gluino in strong gluino pair production models and m(chi0)=450 GeV

Observed (-1 sigma) exclusion limits at 95% CL on the lifetime and mass of the gluino in strong gluino pair production models and m(chi0)=450 GeV

Expected exclusion limits at 95% CL on the mass of the gluino and neutralino in strong gluino pair production models and tau(chi0)=0.01 ns

Expected (+1 sigma) exclusion limits at 95% CL on the mass of the gluino and neutralino in strong gluino pair production models and tau(chi0)=0.01 ns

Expected (-1 sigma) exclusion limits at 95% CL on the mass of the gluino and neutralino in strong gluino pair production models and tau(chi0)=0.01 ns

Observed exclusion limits at 95% CL on the mass of the gluino and neutralino in strong gluino pair production models and tau(chi0)=0.01 ns

Observed (+1 sigma) exclusion limits at 95% CL on the mass of the gluino and neutralino in strong gluino pair production models and tau(chi0)=0.01 ns

Observed (-1 sigma) exclusion limits at 95% CL on the mass of the gluino and neutralino in strong gluino pair production models and tau(chi0)=0.01 ns

Expected exclusion limits at 95% CL on the mass of the gluino and neutralino in strong gluino pair production models and tau(chi0)=0.1 ns

Expected (+1 sigma) exclusion limits at 95% CL on the mass of the gluino and neutralino in strong gluino pair production models and tau(chi0)=0.1 ns

Expected (-1 sigma) exclusion limits at 95% CL on the mass of the gluino and neutralino in strong gluino pair production models and tau(chi0)=0.1 ns

Observed exclusion limits at 95% CL on the mass of the gluino and neutralino in strong gluino pair production models and tau(chi0)=0.1 ns

Observed (+1 sigma) exclusion limits at 95% CL on the mass of the gluino and neutralino in strong gluino pair production models and tau(chi0)=0.1 ns

Observed (-1 sigma) exclusion limits at 95% CL on the mass of the gluino and neutralino in strong gluino pair production models and tau(chi0)=0.1 ns

Expected exclusion limits at 95% CL on the mass of the gluino and neutralino in strong gluino pair production models and tau(chi0)=1 ns

Expected (+1 sigma) exclusion limits at 95% CL on the mass of the gluino and neutralino in strong gluino pair production models and tau(chi0)=1 ns

Expected (-1 sigma) exclusion limits at 95% CL on the mass of the gluino and neutralino in strong gluino pair production models and tau(chi0)=1 ns

Observed exclusion limits at 95% CL on the mass of the gluino and neutralino in strong gluino pair production models and tau(chi0)=1 ns

Observed (+1 sigma) exclusion limits at 95% CL on the mass of the gluino and neutralino in strong gluino pair production models and tau(chi0)=1 ns

Observed (-1 sigma) exclusion limits at 95% CL on the mass of the gluino and neutralino in strong gluino pair production models and tau(chi0)=1 ns

Expected exclusion limits at 95% CL on the mass of the gluino and neutralino in strong gluino pair production models and tau(chi0)=10 ns

Expected (+1 sigma) exclusion limits at 95% CL on the mass of the gluino and neutralino in strong gluino pair production models and tau(chi0)=10 ns

Expected (-1 sigma) exclusion limits at 95% CL on the mass of the gluino and neutralino in strong gluino pair production models and tau(chi0)=10 ns

Observed exclusion limits at 95% CL on the mass of the gluino and neutralino in strong gluino pair production models and tau(chi0)=10 ns

Observed (+1 sigma) exclusion limits at 95% CL on the mass of the gluino and neutralino in strong gluino pair production models and tau(chi0)=10 ns

Observed (-1 sigma) exclusion limits at 95% CL on the mass of the gluino and neutralino in strong gluino pair production models and tau(chi0)=10 ns

Exclusion limits at 95% CL on the production cross section in the strong gluino pair production models and m($ ilde{\chi}^0_1$)=1.25 TeV

Acceptance cutflow for the High-pT SR for representative points in the strong gluino pair production model. See additional resources for more information.

Acceptance cutflow for the Trackless SR for representative points in the electroweak pair production model. See additional resources for more information.

Acceptance cutflow for the Trackless SR for representative points in the electroweak pair production model with heavy-flavor quarks final state. See additional resources for more information.

Acceptance cutflow for the High-pT SR for representative points in the electroweak pair production model with heavy-flavor quarks final state. See additional resources for more information.

Reinterpretation Material: Event-level Efficiency for HighPt SR selections, R < 1150 mm

Reinterpretation Material: Event-level Efficiency for HighPt SR selections, R [1150, 3870] mm

Reinterpretation Material: Event-level Efficiency for HighPt SR selections, R > 3870 mm

Reinterpretation Material: Event-level Efficiency for Trackless SR selections, R < 1150 mm

Reinterpretation Material: Event-level Efficiency for Trackless SR selections, R [1150, 3870] mm

Reinterpretation Material: Event-level Efficiency for Trackless SR selections, R > 3870 mm

Reinterpretation Material: Vertex-level Efficiency for R < 22 mm

Reinterpretation Material: Vertex-level Efficiency for R [22, 25] mm

Reinterpretation Material: Vertex-level Efficiency for R [25, 29] mm

Reinterpretation Material: Vertex-level Efficiency for R [29, 38] mm

Reinterpretation Material: Vertex-level Efficiency for R [38, 46] mm

Reinterpretation Material: Vertex-level Efficiency for R [46, 73] mm

Reinterpretation Material: Vertex-level Efficiency for R [73, 84] mm

Reinterpretation Material: Vertex-level Efficiency for R [84, 111] mm

Reinterpretation Material: Vertex-level Efficiency for R [111, 120] mm

Reinterpretation Material: Vertex-level Efficiency for R [120, 145] mm

Reinterpretation Material: Vertex-level Efficiency for R [145, 180] mm

Reinterpretation Material: Vertex-level Efficiency for R [180, 300] mm

Cutflow (acceptance x efficiency) for the High-pT SR for representative points in the strong gluino pair production model. See additional resources for more information.

Cutflow (acceptance x efficiency) for the Trackless SR for representative points in the electroweak pair production model. See additional resources for more information.

Cutflow (acceptance x efficiency) for the Trackless SR for representative points in the electroweak pair production model with heavy-flavor quarks. See additional resources for more information.

Cutflow (acceptance x efficiency) for the High-pT SR for representative points in the electroweak pair production model with heavy-flavor quarks. See additional resources for more information.

Feynman diagrams of the effective neutralino pair production in the GMSB simplified model in which the two neutralinos decay into two gravitinos ($\tilde{G}$) and two $Z$ bosons (left), a $Z$ and a Higgs boson ($H$) (center), or two Higgs bosons (right).

The distributions of the most impactful input variables to the TD jet tagger for signal (red, lighter) and collision background (blue, darker). They include the charged (upper left) and neutral (upper right) hadron energy fractions, the number of track constituents in the jet (middle left), the $\Delta R$ between the jet axis and the closest track associated with the PV (middle right), and the jet time (lower).

The distributions of the most impactful input variables to the TD jet tagger for signal (red, lighter) and collision background (blue, darker). They include the charged (upper left) and neutral (upper right) hadron energy fractions, the number of track constituents in the jet (middle left), the $\Delta R$ between the jet axis and the closest track associated with the PV (middle right), and the jet time (lower).

The distributions of the most impactful input variables to the TD jet tagger for signal (red, lighter) and collision background (blue, darker). They include the charged (upper left) and neutral (upper right) hadron energy fractions, the number of track constituents in the jet (middle left), the $\Delta R$ between the jet axis and the closest track associated with the PV (middle right), and the jet time (lower).

The distributions of the most impactful input variables to the TD jet tagger for signal (red, lighter) and collision background (blue, darker). They include the charged (upper left) and neutral (upper right) hadron energy fractions, the number of track constituents in the jet (middle left), the $\Delta R$ between the jet axis and the closest track associated with the PV (middle right), and the jet time (lower).

The distributions of the most impactful input variables to the TD jet tagger for signal (red, lighter) and collision background (blue, darker). They include the charged (upper left) and neutral (upper right) hadron energy fractions, the number of track constituents in the jet (middle left), the $\Delta R$ between the jet axis and the closest track associated with the PV (middle right), and the jet time (lower).

TD jet tagger score distributions (left) for signal (red, lighter) and collision background (blue, darker). Identification probability for the signal versus the misidentification probability for the background (right) with the tagger working point (w.~p.) used in the analysis shown as a blue marker. Both are evaluated using an independent sample of testing events.

TD jet tagger score distributions (left) for signal (red, lighter) and collision background (blue, darker). Identification probability for the signal versus the misidentification probability for the background (right) with the tagger working point (w.~p.) used in the analysis shown as a blue marker. Both are evaluated using an independent sample of testing events.

The efficiency of the TD jet tagger working point used in the analysis is shown as a function of the lab frame transverse decay length for simulated signals with $\chi$ mass of 400 GeV. The uncertainties shown account for lifetime dependence and statistical uncertainty.

The TD jet tagger score distributions for simulation (shaded histogram) and data (black markers) when using electrons from $W\to e\nu_e$ events as proxy objects for signal jets. The histograms and data points have been normalized to unit area. The last bin contains jets with tagger scores greater than 0.996, the threshold used to tag signal jets. Similar levels of agreement are observed for photon proxy objects from the $Z\to\ell^+\ell^-\gamma$ sample.

The TD jet tagger misidentification probability measured using the nominal $W$+jets MR (black round markers) is shown along with the systematic uncertainty (gray band), quantifying the degree of process dependence measured from alternative MRs. The measurements in the alternative MRs are displayed as well ($Z$+jets MR as green round markers, $t\bar{t}$ MR as red squared markers, QCD MR as blue triangular markers) along with their respective statistical uncertainty. On the left, this probability is shown for the first 19.9 fb$^{-1}$ of data collected in 2016, while on the right it is shown for the last 16.4 fb$^{-1}$ of data collected in 2016combined with data collected in 2017-2018.

The TD jet tagger misidentification probability measured using the nominal $W$+jets MR (black round markers) is shown along with the systematic uncertainty (gray band), quantifying the degree of process dependence measured from alternative MRs. The measurements in the alternative MRs are displayed as well ($Z$+jets MR as green round markers, $t\bar{t}$ MR as red squared markers, QCD MR as blue triangular markers) along with their respective statistical uncertainty. On the left, this probability is shown for the first 19.9 fb$^{-1}$ of data collected in 2016, while on the right it is shown for the last 16.4 fb$^{-1}$ of data collected in 2016combined with data collected in 2017-2018.

Distribution of the number of TD tagged jets for the $m_{\chi} = 400$ GeVsimulated signal samples with $c\tau_{\chi} = 0.5$ m (solid red line) and $c\tau_{\chi} = 3.0$ m (dotted green line), estimated background (blue square markers), and data (black round markers). The signal distributions are normalized to the expected cross section limit. The blue shaded region indicates the systematic uncertainty in the background prediction. No background prediction is shown for the bin with zero TD tagged jets as it is the main control region used to predict the background for the other two bins. There are zero observed events in the bin with two or more TD tagged jets.

Expected and observed 95% CL upper limits on $\sigma_{\chi\chi}$ as functions of $m_\chi$ in a scenario with $\mathcal{B}(\chi\to HG) = 0.5$ and $c\tau = 0.5$ m (left) or 3 m (right).

Expected and observed 95% CL upper limits on $\sigma_{\chi\chi}$ as functions of $m_\chi$ in a scenario with $\mathcal{B}(\chi\to HG) = 0.5$ and $c\tau = 0.5$ m (left) or 3 m (right).

Expected and observed 95% CL upper limits on $\sigma_{\chi\chi}$ as functions of $c\tau_{\chi}$ in a scenario with $\mathcal{B}(\chi\to H\tilde{G}) = 0.5$ and $m_{\chi} = 400$ GeV (left) or 1000 GeV (right).

Expected and observed 95% CL upper limits on $\sigma_{\chi\chi}$ as functions of $c\tau_{\chi}$ in a scenario with $\mathcal{B}(\chi\to H\tilde{G}) = 0.5$ and $m_{\chi} = 400$ GeV (left) or 1000 GeV (right).

The observed 95% CL upper limit on $\sigma_{\chi\chi}$ as a function of $m_{\chi}$ and $c\tau_{\chi}$ in a scenario with $\mathcal{B}(\chi\to H\tilde{G}) = 0.5$. The area enclosed by the dotted black line corresponds to the observed excluded region.

The distribution of the jet charged hadron energy fraction, a variable used as input to the TD jet tagger score, for simulation (shaded histogram) and data (black markers) when using electrons from $W\to e\nu_e$ events as proxy objects for signal jets. The histograms and data points have been normalized to unit area. Similar levels of agreement are observed for photon proxy objects from the $Z\to\ell^+\ell^-\gamma$ sample.

The distribution of the jet neutral hadron energy fraction, a variable used as input to the TD jet tagger score, for simulation (shaded histogram) and data (black markers) when using electrons from $W\to e\nu_e$ events as proxy objects for signal jets. The histograms and data points have been normalized to unit area. Similar levels of agreement are observed for photon proxy objects from the $Z\to\ell^+\ell^-\gamma$ sample.

The distribution of the number of track constituents in the jet, a variable used as input to the TD jet tagger score, for simulation (shaded histogram) and data (black markers) when using electrons from $W\to e\nu_e$ events as proxy objects for signal jets. The histograms and data points have been normalized to unit area. Similar levels of agreement are observed for photon proxy objects from the $Z\to\ell^+\ell^-\gamma$ sample.

The $\eta$ distribution of TD-tagged jets in a background-enriched control region that comprises events with exactly one TD-tagged jet. Observed data (black round markers) and the corresponding prediction based on control samples in data (empty squared markers), measured using the nominal $W$+jets MR, are compared. The prediction uncertainty (gray band) includes the systematic uncertainty quantifying the degree of process dependence measured from alternative MRs. The predictions for the shape and the normalization of the $\eta$ distribution are consistent with the data.

Jet time distribution in a sample of b-tagged jets from dilepton $t \bar{t}$ events in 2017 data-taking period (black round markers) and simulation (filled histogram). A Gaussian smearing procedure is applied to the jet time in the $t \bar{t}$ sample (green line) to correct for effects that are difficult to simulate (timing drift caused by crystals transparency loss due to detector aging, electronics jitter).

The observed 95% CL upper limit on $\sigma_{\chi\chi}$ as a function of $m_{\chi}$ and $c\tau_{\chi}$ in a scenario with $\mathcal{B}(\chi\to H\tilde{G}) = 1$. The area enclosed by the dotted black line corresponds to the observed excluded region.

The observed 95% CL upper limit on $\sigma_{\chi\chi}$ as a function of $m_{\chi}$ and $c\tau_{\chi}$ in a scenario with $\mathcal{B}(\chi\to H\tilde{G}) = 0.75$, $\mathcal{B}(\chi\to Z\tilde{G}) = 0.25$. The area enclosed by the dotted black line corresponds to the observed excluded region.

The observed 95% CL upper limit on $\sigma_{\chi\chi}$ as a function of $m_{\chi}$ and $c\tau_{\chi}$ in a scenario with $\mathcal{B}(\chi\to H\tilde{G}) = 0.25$, $\mathcal{B}(\chi\to Z\tilde{G}) = 0.75$. The area enclosed by the dotted black line corresponds to the observed excluded region.

The observed 95% CL upper limit on $\sigma_{\chi\chi}$ as a function of $m_{\chi}$ and $c\tau_{\chi}$ in a scenario with $\mathcal{B}(\chi\to Z\tilde{G}) = 1$. The area enclosed by the dotted black line corresponds to the observed excluded region.

The efficiency of identifying a LLP decay as a TD-tagged jet in bins of the LLP transverse and longitudinal decay position. The sample used to compute the efficiency contains events with pair production of $\chi$ with a mass of 127 GeV and a lifetime of 0.5 and 3 m, and considering the combinations of the $\chi$ decay modes considered in this search ($H \tilde{G} \rightarrow b\bar{b} \tilde{G}$ or $Z\tilde{G} \rightarrow q\bar{q} \tilde{G}$). The efficiency is calculated on top of the acceptance definition for the merged topology, namely, the H (or Z) decay products are produced with an angular separation $\Delta R < 0.8$, and the H (or Z) has $p_T > 30$ GeV and $|\eta|<1$. The full simulation signal yield prediction can be reproduced within 3\%. This nonclosure uncertainty is added in quadrature to the statistical uncertainty of each bin.

The efficiency of identifying a LLP decay as a TD-tagged jet in bins of the LLP transverse and longitudinal decay position. The sample used to compute the efficiency contains events with pair production of $\chi$ with a mass of 127 GeV and a lifetime of 0.5 and 3 m, and considering the combinations of the $\chi$ decay modes considered in this search ($H \tilde{G} \rightarrow b\bar{b} \tilde{G}$ or $Z\tilde{G} \rightarrow q\bar{q} \tilde{G}$). The efficiency is calculated on top of the acceptance definition for the resolved topology with exactly one quark in acceptance, namely, the H (or Z) decay products are produced with an angular separation $\Delta R \geq 0.8$, and only one b-quark (or light quark) has $p_T > 30$ GeV and $|\eta|<1$. The full simulation signal yield prediction can be reproduced within 5\%. This nonclosure uncertainty is added in quadrature to the statistical uncertainty of each bin.

The efficiency of identifying a LLP decay as a TD-tagged jet in bins of the LLP transverse and longitudinal decay position. The sample used to compute the efficiency contains events with pair production of $\chi$ with a mass of 127 GeV and a lifetime of 0.5 and 3 m, and considering the combinations of the $\chi$ decay modes considered in this search ($H \tilde{G} \rightarrow b\bar{b} \tilde{G}$ or $Z\tilde{G} \rightarrow q\bar{q} \tilde{G}$). The efficiency is calculated on top of the acceptance definition for the resolved topology with two quarks in acceptance, namely, the H (or Z) decay products are produced with an angular separation $\Delta R \geq 0.8$, and both b-quarks (or light quarks) have $p_T > 30$ GeV and $|\eta|<1$. The full simulation signal yield prediction can be reproduced within 7\%. This nonclosure uncertainty is added in quadrature to the statistical uncertainty of each bin.

The efficiency of identifying a LLP decay as a TD-tagged jet in bins of the LLP transverse and longitudinal decay position. The sample used to compute the efficiency contains events with pair production of $\chi$ with a mass of 150 GeV and a lifetime of 0.5 and 3 m, and considering the combinations of the $\chi$ decay modes considered in this search ($H \tilde{G} \rightarrow b\bar{b} \tilde{G}$ or $Z\tilde{G} \rightarrow q\bar{q} \tilde{G}$). The efficiency is calculated on top of the acceptance definition for the merged topology, namely, the H (or Z) decay products are produced with an angular separation $\Delta R < 0.8$, and the H (or Z) has $p_T > 30$ GeV and $|\eta|<1$. The full simulation signal yield prediction can be reproduced within 3\%. This nonclosure uncertainty is added in quadrature to the statistical uncertainty of each bin.

The efficiency of identifying a LLP decay as a TD-tagged jet in bins of the LLP transverse and longitudinal decay position. The sample used to compute the efficiency contains events with pair production of $\chi$ with a mass of 150 GeV and a lifetime of 0.5 and 3 m, and considering the combinations of the $\chi$ decay modes considered in this search ($H \tilde{G} \rightarrow b\bar{b} \tilde{G}$ or $Z\tilde{G} \rightarrow q\bar{q} \tilde{G}$). The efficiency is calculated on top of the acceptance definition for the resolved topology with exactly one quark in acceptance, namely, the H (or Z) decay products are produced with an angular separation $\Delta R \geq 0.8$, and only one b-quark (or light quark) has $p_T > 30$ GeV and $|\eta|<1$. The full simulation signal yield prediction can be reproduced within 5\%. This nonclosure uncertainty is added in quadrature to the statistical uncertainty of each bin.

The efficiency of identifying a LLP decay as a TD-tagged jet in bins of the LLP transverse and longitudinal decay position. The sample used to compute the efficiency contains events with pair production of $\chi$ with a mass of 150 GeV and a lifetime of 0.5 and 3 m, and considering the combinations of the $\chi$ decay modes considered in this search ($H \tilde{G} \rightarrow b\bar{b} \tilde{G}$ or $Z\tilde{G} \rightarrow q\bar{q} \tilde{G}$). The efficiency is calculated on top of the acceptance definition for the resolved topology with two quarks in acceptance, namely, the H (or Z) decay products are produced with an angular separation $\Delta R \geq 0.8$, and both b-quarks (or light quarks) have $p_T > 30$ GeV and $|\eta|<1$. The full simulation signal yield prediction can be reproduced within 7\%. This nonclosure uncertainty is added in quadrature to the statistical uncertainty of each bin.

The efficiency of identifying a LLP decay as a TD-tagged jet in bins of the LLP transverse and longitudinal decay position. The sample used to compute the efficiency contains events with pair production of $\chi$ with a mass of 175 GeV and a lifetime of 0.5 and 3 m, and considering the combinations of the $\chi$ decay modes considered in this search ($H \tilde{G} \rightarrow b\bar{b} \tilde{G}$ or $Z\tilde{G} \rightarrow q\bar{q} \tilde{G}$). The efficiency is calculated on top of the acceptance definition for the merged topology, namely, the H (or Z) decay products are produced with an angular separation $\Delta R < 0.8$, and the H (or Z) has $p_T > 30$ GeV and $|\eta|<1$. The full simulation signal yield prediction can be reproduced within 3\%. This nonclosure uncertainty is added in quadrature to the statistical uncertainty of each bin.

The efficiency of identifying a LLP decay as a TD-tagged jet in bins of the LLP transverse and longitudinal decay position. The sample used to compute the efficiency contains events with pair production of $\chi$ with a mass of 175 GeV and a lifetime of 0.5 and 3 m, and considering the combinations of the $\chi$ decay modes considered in this search ($H \tilde{G} \rightarrow b\bar{b} \tilde{G}$ or $Z\tilde{G} \rightarrow q\bar{q} \tilde{G}$). The efficiency is calculated on top of the acceptance definition for the resolved topology with exactly one quark in acceptance, namely, the H (or Z) decay products are produced with an angular separation $\Delta R \geq 0.8$, and only one b-quark (or light quark) has $p_T > 30$ GeV and $|\eta|<1$. The full simulation signal yield prediction can be reproduced within 5\%. This nonclosure uncertainty is added in quadrature to the statistical uncertainty of each bin.

The efficiency of identifying a LLP decay as a TD-tagged jet in bins of the LLP transverse and longitudinal decay position. The sample used to compute the efficiency contains events with pair production of $\chi$ with a mass of 175 GeV and a lifetime of 0.5 and 3 m, and considering the combinations of the $\chi$ decay modes considered in this search ($H \tilde{G} \rightarrow b\bar{b} \tilde{G}$ or $Z\tilde{G} \rightarrow q\bar{q} \tilde{G}$). The efficiency is calculated on top of the acceptance definition for the resolved topology with two quarks in acceptance, namely, the H (or Z) decay products are produced with an angular separation $\Delta R \geq 0.8$, and both b-quarks (or light quarks) have $p_T > 30$ GeV and $|\eta|<1$. The full simulation signal yield prediction can be reproduced within 7\%. This nonclosure uncertainty is added in quadrature to the statistical uncertainty of each bin.

The efficiency of identifying a LLP decay as a TD-tagged jet in bins of the LLP transverse and longitudinal decay position. The sample used to compute the efficiency contains events with pair production of $\chi$ with a mass of 200 GeV and a lifetime of 0.5 and 3 m, and considering the combinations of the $\chi$ decay modes considered in this search ($H \tilde{G} \rightarrow b\bar{b} \tilde{G}$ or $Z\tilde{G} \rightarrow q\bar{q} \tilde{G}$). The efficiency is calculated on top of the acceptance definition for the merged topology, namely, the H (or Z) decay products are produced with an angular separation $\Delta R < 0.8$, and the H (or Z) has $p_T > 30$ GeV and $|\eta|<1$. The full simulation signal yield prediction can be reproduced within 3\%. This nonclosure uncertainty is added in quadrature to the statistical uncertainty of each bin.

The efficiency of identifying a LLP decay as a TD-tagged jet in bins of the LLP transverse and longitudinal decay position. The sample used to compute the efficiency contains events with pair production of $\chi$ with a mass of 200 GeV and a lifetime of 0.5 and 3 m, and considering the combinations of the $\chi$ decay modes considered in this search ($H \tilde{G} \rightarrow b\bar{b} \tilde{G}$ or $Z\tilde{G} \rightarrow q\bar{q} \tilde{G}$). The efficiency is calculated on top of the acceptance definition for the resolved topology with exactly one quark in acceptance, namely, the H (or Z) decay products are produced with an angular separation $\Delta R \geq 0.8$, and only one b-quark (or light quark) has $p_T > 30$ GeV and $|\eta|<1$. The full simulation signal yield prediction can be reproduced within 5\%. This nonclosure uncertainty is added in quadrature to the statistical uncertainty of each bin.

The efficiency of identifying a LLP decay as a TD-tagged jet in bins of the LLP transverse and longitudinal decay position. The sample used to compute the efficiency contains events with pair production of $\chi$ with a mass of 200 GeV and a lifetime of 0.5 and 3 m, and considering the combinations of the $\chi$ decay modes considered in this search ($H \tilde{G} \rightarrow b\bar{b} \tilde{G}$ or $Z\tilde{G} \rightarrow q\bar{q} \tilde{G}$). The efficiency is calculated on top of the acceptance definition for the resolved topology with two quarks in acceptance, namely, the H (or Z) decay products are produced with an angular separation $\Delta R \geq 0.8$, and both b-quarks (or light quarks) have $p_T > 30$ GeV and $|\eta|<1$. The full simulation signal yield prediction can be reproduced within 7\%. This nonclosure uncertainty is added in quadrature to the statistical uncertainty of each bin.

The efficiency of identifying a LLP decay as a TD-tagged jet in bins of the LLP transverse and longitudinal decay position. The sample used to compute the efficiency contains events with pair production of $\chi$ with a mass of 250 GeV and a lifetime of 0.5 and 3 m, and considering the combinations of the $\chi$ decay modes considered in this search ($H \tilde{G} \rightarrow b\bar{b} \tilde{G}$ or $Z\tilde{G} \rightarrow q\bar{q} \tilde{G}$). The efficiency is calculated on top of the acceptance definition for the merged topology, namely, the H (or Z) decay products are produced with an angular separation $\Delta R < 0.8$, and the H (or Z) has $p_T > 30$ GeV and $|\eta|<1$. The full simulation signal yield prediction can be reproduced within 3\%. This nonclosure uncertainty is added in quadrature to the statistical uncertainty of each bin.

The efficiency of identifying a LLP decay as a TD-tagged jet in bins of the LLP transverse and longitudinal decay position. The sample used to compute the efficiency contains events with pair production of $\chi$ with a mass of 250 GeV and a lifetime of 0.5 and 3 m, and considering the combinations of the $\chi$ decay modes considered in this search ($H \tilde{G} \rightarrow b\bar{b} \tilde{G}$ or $Z\tilde{G} \rightarrow q\bar{q} \tilde{G}$). The efficiency is calculated on top of the acceptance definition for the resolved topology with exactly one quark in acceptance, namely, the H (or Z) decay products are produced with an angular separation $\Delta R \geq 0.8$, and only one b-quark (or light quark) has $p_T > 30$ GeV and $|\eta|<1$. The full simulation signal yield prediction can be reproduced within 5\%. This nonclosure uncertainty is added in quadrature to the statistical uncertainty of each bin.

The efficiency of identifying a LLP decay as a TD-tagged jet in bins of the LLP transverse and longitudinal decay position. The sample used to compute the efficiency contains events with pair production of $\chi$ with a mass of 250 GeV and a lifetime of 0.5 and 3 m, and considering the combinations of the $\chi$ decay modes considered in this search ($H \tilde{G} \rightarrow b\bar{b} \tilde{G}$ or $Z\tilde{G} \rightarrow q\bar{q} \tilde{G}$). The efficiency is calculated on top of the acceptance definition for the resolved topology with two quarks in acceptance, namely, the H (or Z) decay products are produced with an angular separation $\Delta R \geq 0.8$, and both b-quarks (or light quarks) have $p_T > 30$ GeV and $|\eta|<1$. The full simulation signal yield prediction can be reproduced within 7\%. This nonclosure uncertainty is added in quadrature to the statistical uncertainty of each bin.

The efficiency of identifying a LLP decay as a TD-tagged jet in bins of the LLP transverse and longitudinal decay position. The sample used to compute the efficiency contains events with pair production of $\chi$ with a mass of 300 GeV and a lifetime of 0.5 and 3 m, and considering the combinations of the $\chi$ decay modes considered in this search ($H \tilde{G} \rightarrow b\bar{b} \tilde{G}$ or $Z\tilde{G} \rightarrow q\bar{q} \tilde{G}$). The efficiency is calculated on top of the acceptance definition for the merged topology, namely, the H (or Z) decay products are produced with an angular separation $\Delta R < 0.8$, and the H (or Z) has $p_T > 30$ GeV and $|\eta|<1$. The full simulation signal yield prediction can be reproduced within 3\%. This nonclosure uncertainty is added in quadrature to the statistical uncertainty of each bin.

The efficiency of identifying a LLP decay as a TD-tagged jet in bins of the LLP transverse and longitudinal decay position. The sample used to compute the efficiency contains events with pair production of $\chi$ with a mass of 300 GeV and a lifetime of 0.5 and 3 m, and considering the combinations of the $\chi$ decay modes considered in this search ($H \tilde{G} \rightarrow b\bar{b} \tilde{G}$ or $Z\tilde{G} \rightarrow q\bar{q} \tilde{G}$). The efficiency is calculated on top of the acceptance definition for the resolved topology with exactly one quark in acceptance, namely, the H (or Z) decay products are produced with an angular separation $\Delta R \geq 0.8$, and only one b-quark (or light quark) has $p_T > 30$ GeV and $|\eta|<1$. The full simulation signal yield prediction can be reproduced within 5\%. This nonclosure uncertainty is added in quadrature to the statistical uncertainty of each bin.

The efficiency of identifying a LLP decay as a TD-tagged jet in bins of the LLP transverse and longitudinal decay position. The sample used to compute the efficiency contains events with pair production of $\chi$ with a mass of 300 GeV and a lifetime of 0.5 and 3 m, and considering the combinations of the $\chi$ decay modes considered in this search ($H \tilde{G} \rightarrow b\bar{b} \tilde{G}$ or $Z\tilde{G} \rightarrow q\bar{q} \tilde{G}$). The efficiency is calculated on top of the acceptance definition for the resolved topology with two quarks in acceptance, namely, the H (or Z) decay products are produced with an angular separation $\Delta R \geq 0.8$, and both b-quarks (or light quarks) have $p_T > 30$ GeV and $|\eta|<1$. The full simulation signal yield prediction can be reproduced within 7\%. This nonclosure uncertainty is added in quadrature to the statistical uncertainty of each bin.

The efficiency of identifying a LLP decay as a TD-tagged jet in bins of the LLP transverse and longitudinal decay position. The sample used to compute the efficiency contains events with pair production of $\chi$ with a mass of 400 GeV and a lifetime of 0.5 and 3 m, and considering the combinations of the $\chi$ decay modes considered in this search ($H \tilde{G} \rightarrow b\bar{b} \tilde{G}$ or $Z\tilde{G} \rightarrow q\bar{q} \tilde{G}$). The efficiency is calculated on top of the acceptance definition for the merged topology, namely, the H (or Z) decay products are produced with an angular separation $\Delta R < 0.8$, and the H (or Z) has $p_T > 30$ GeV and $|\eta|<1$. The full simulation signal yield prediction can be reproduced within 3\%. This nonclosure uncertainty is added in quadrature to the statistical uncertainty of each bin.

The efficiency of identifying a LLP decay as a TD-tagged jet in bins of the LLP transverse and longitudinal decay position. The sample used to compute the efficiency contains events with pair production of $\chi$ with a mass of 400 GeV and a lifetime of 0.5 and 3 m, and considering the combinations of the $\chi$ decay modes considered in this search ($H \tilde{G} \rightarrow b\bar{b} \tilde{G}$ or $Z\tilde{G} \rightarrow q\bar{q} \tilde{G}$). The efficiency is calculated on top of the acceptance definition for the resolved topology with exactly one quark in acceptance, namely, the H (or Z) decay products are produced with an angular separation $\Delta R \geq 0.8$, and only one b-quark (or light quark) has $p_T > 30$ GeV and $|\eta|<1$. The full simulation signal yield prediction can be reproduced within 5\%. This nonclosure uncertainty is added in quadrature to the statistical uncertainty of each bin.

The efficiency of identifying a LLP decay as a TD-tagged jet in bins of the LLP transverse and longitudinal decay position. The sample used to compute the efficiency contains events with pair production of $\chi$ with a mass of 400 GeV and a lifetime of 0.5 and 3 m, and considering the combinations of the $\chi$ decay modes considered in this search ($H \tilde{G} \rightarrow b\bar{b} \tilde{G}$ or $Z\tilde{G} \rightarrow q\bar{q} \tilde{G}$). The efficiency is calculated on top of the acceptance definition for the resolved topology with two quarks in acceptance, namely, the H (or Z) decay products are produced with an angular separation $\Delta R \geq 0.8$, and both b-quarks (or light quarks) have $p_T > 30$ GeV and $|\eta|<1$. The full simulation signal yield prediction can be reproduced within 7\%. This nonclosure uncertainty is added in quadrature to the statistical uncertainty of each bin.

The efficiency of identifying a LLP decay as a TD-tagged jet in bins of the LLP transverse and longitudinal decay position. The sample used to compute the efficiency contains events with pair production of $\chi$ with a mass of 600 GeV and a lifetime of 0.5 and 3 m, and considering the combinations of the $\chi$ decay modes considered in this search ($H \tilde{G} \rightarrow b\bar{b} \tilde{G}$ or $Z\tilde{G} \rightarrow q\bar{q} \tilde{G}$). The efficiency is calculated on top of the acceptance definition for the merged topology, namely, the H (or Z) decay products are produced with an angular separation $\Delta R < 0.8$, and the H (or Z) has $p_T > 30$ GeV and $|\eta|<1$. The full simulation signal yield prediction can be reproduced within 3\%. This nonclosure uncertainty is added in quadrature to the statistical uncertainty of each bin.

The efficiency of identifying a LLP decay as a TD-tagged jet in bins of the LLP transverse and longitudinal decay position. The sample used to compute the efficiency contains events with pair production of $\chi$ with a mass of 600 GeV and a lifetime of 0.5 and 3 m, and considering the combinations of the $\chi$ decay modes considered in this search ($H \tilde{G} \rightarrow b\bar{b} \tilde{G}$ or $Z\tilde{G} \rightarrow q\bar{q} \tilde{G}$). The efficiency is calculated on top of the acceptance definition for the resolved topology with exactly one quark in acceptance, namely, the H (or Z) decay products are produced with an angular separation $\Delta R \geq 0.8$, and only one b-quark (or light quark) has $p_T > 30$ GeV and $|\eta|<1$. The full simulation signal yield prediction can be reproduced within 5\%. This nonclosure uncertainty is added in quadrature to the statistical uncertainty of each bin.

The efficiency of identifying a LLP decay as a TD-tagged jet in bins of the LLP transverse and longitudinal decay position. The sample used to compute the efficiency contains events with pair production of $\chi$ with a mass of 600 GeV and a lifetime of 0.5 and 3 m, and considering the combinations of the $\chi$ decay modes considered in this search ($H \tilde{G} \rightarrow b\bar{b} \tilde{G}$ or $Z\tilde{G} \rightarrow q\bar{q} \tilde{G}$). The efficiency is calculated on top of the acceptance definition for the resolved topology with two quarks in acceptance, namely, the H (or Z) decay products are produced with an angular separation $\Delta R \geq 0.8$, and both b-quarks (or light quarks) have $p_T > 30$ GeV and $|\eta|<1$. The full simulation signal yield prediction can be reproduced within 7\%. This nonclosure uncertainty is added in quadrature to the statistical uncertainty of each bin.

The efficiency of identifying a LLP decay as a TD-tagged jet in bins of the LLP transverse and longitudinal decay position. The sample used to compute the efficiency contains events with pair production of $\chi$ with a mass of 800 GeV and a lifetime of 0.5 and 3 m, and considering the combinations of the $\chi$ decay modes considered in this search ($H \tilde{G} \rightarrow b\bar{b} \tilde{G}$ or $Z\tilde{G} \rightarrow q\bar{q} \tilde{G}$). The efficiency is calculated on top of the acceptance definition for the merged topology, namely, the H (or Z) decay products are produced with an angular separation $\Delta R < 0.8$, and the H (or Z) has $p_T > 30$ GeV and $|\eta|<1$. The full simulation signal yield prediction can be reproduced within 3\%. This nonclosure uncertainty is added in quadrature to the statistical uncertainty of each bin.

The efficiency of identifying a LLP decay as a TD-tagged jet in bins of the LLP transverse and longitudinal decay position. The sample used to compute the efficiency contains events with pair production of $\chi$ with a mass of 800 GeV and a lifetime of 0.5 and 3 m, and considering the combinations of the $\chi$ decay modes considered in this search ($H \tilde{G} \rightarrow b\bar{b} \tilde{G}$ or $Z\tilde{G} \rightarrow q\bar{q} \tilde{G}$). The efficiency is calculated on top of the acceptance definition for the resolved topology with exactly one quark in acceptance, namely, the H (or Z) decay products are produced with an angular separation $\Delta R \geq 0.8$, and only one b-quark (or light quark) has $p_T > 30$ GeV and $|\eta|<1$. The full simulation signal yield prediction can be reproduced within 5\%. This nonclosure uncertainty is added in quadrature to the statistical uncertainty of each bin.

The efficiency of identifying a LLP decay as a TD-tagged jet in bins of the LLP transverse and longitudinal decay position. The sample used to compute the efficiency contains events with pair production of $\chi$ with a mass of 800 GeV and a lifetime of 0.5 and 3 m, and considering the combinations of the $\chi$ decay modes considered in this search ($H \tilde{G} \rightarrow b\bar{b} \tilde{G}$ or $Z\tilde{G} \rightarrow q\bar{q} \tilde{G}$). The efficiency is calculated on top of the acceptance definition for the resolved topology with two quarks in acceptance, namely, the H (or Z) decay products are produced with an angular separation $\Delta R \geq 0.8$, and both b-quarks (or light quarks) have $p_T > 30$ GeV and $|\eta|<1$. The full simulation signal yield prediction can be reproduced within 7\%. This nonclosure uncertainty is added in quadrature to the statistical uncertainty of each bin.

The efficiency of identifying a LLP decay as a TD-tagged jet in bins of the LLP transverse and longitudinal decay position. The sample used to compute the efficiency contains events with pair production of $\chi$ with a mass of 1000 GeV and a lifetime of 0.5 and 3 m, and considering the combinations of the $\chi$ decay modes considered in this search ($H \tilde{G} \rightarrow b\bar{b} \tilde{G}$ or $Z\tilde{G} \rightarrow q\bar{q} \tilde{G}$). The efficiency is calculated on top of the acceptance definition for the merged topology, namely, the H (or Z) decay products are produced with an angular separation $\Delta R < 0.8$, and the H (or Z) has $p_T > 30$ GeV and $|\eta|<1$. The full simulation signal yield prediction can be reproduced within 3\%. This nonclosure uncertainty is added in quadrature to the statistical uncertainty of each bin.

The efficiency of identifying a LLP decay as a TD-tagged jet in bins of the LLP transverse and longitudinal decay position. The sample used to compute the efficiency contains events with pair production of $\chi$ with a mass of 1000 GeV and a lifetime of 0.5 and 3 m, and considering the combinations of the $\chi$ decay modes considered in this search ($H \tilde{G} \rightarrow b\bar{b} \tilde{G}$ or $Z\tilde{G} \rightarrow q\bar{q} \tilde{G}$). The efficiency is calculated on top of the acceptance definition for the resolved topology with exactly one quark in acceptance, namely, the H (or Z) decay products are produced with an angular separation $\Delta R \geq 0.8$, and only one b-quark (or light quark) has $p_T > 30$ GeV and $|\eta|<1$. The full simulation signal yield prediction can be reproduced within 5\%. This nonclosure uncertainty is added in quadrature to the statistical uncertainty of each bin.

The efficiency of identifying a LLP decay as a TD-tagged jet in bins of the LLP transverse and longitudinal decay position. The sample used to compute the efficiency contains events with pair production of $\chi$ with a mass of 1000 GeV and a lifetime of 0.5 and 3 m, and considering the combinations of the $\chi$ decay modes considered in this search ($H \tilde{G} \rightarrow b\bar{b} \tilde{G}$ or $Z\tilde{G} \rightarrow q\bar{q} \tilde{G}$). The efficiency is calculated on top of the acceptance definition for the resolved topology with two quarks in acceptance, namely, the H (or Z) decay products are produced with an angular separation $\Delta R \geq 0.8$, and both b-quarks (or light quarks) have $p_T > 30$ GeV and $|\eta|<1$. The full simulation signal yield prediction can be reproduced within 7\%. This nonclosure uncertainty is added in quadrature to the statistical uncertainty of each bin.

The efficiency of identifying a LLP decay as a TD-tagged jet in bins of the LLP transverse and longitudinal decay position. The sample used to compute the efficiency contains events with pair production of $\chi$ with a mass of 1250 GeV and a lifetime of 0.5 and 3 m, and considering the combinations of the $\chi$ decay modes considered in this search ($H \tilde{G} \rightarrow b\bar{b} \tilde{G}$ or $Z\tilde{G} \rightarrow q\bar{q} \tilde{G}$). The efficiency is calculated on top of the acceptance definition for the merged topology, namely, the H (or Z) decay products are produced with an angular separation $\Delta R < 0.8$, and the H (or Z) has $p_T > 30$ GeV and $|\eta|<1$. The full simulation signal yield prediction can be reproduced within 3\%. This nonclosure uncertainty is added in quadrature to the statistical uncertainty of each bin.

The efficiency of identifying a LLP decay as a TD-tagged jet in bins of the LLP transverse and longitudinal decay position. The sample used to compute the efficiency contains events with pair production of $\chi$ with a mass of 1250 GeV and a lifetime of 0.5 and 3 m, and considering the combinations of the $\chi$ decay modes considered in this search ($H \tilde{G} \rightarrow b\bar{b} \tilde{G}$ or $Z\tilde{G} \rightarrow q\bar{q} \tilde{G}$). The efficiency is calculated on top of the acceptance definition for the resolved topology with exactly one quark in acceptance, namely, the H (or Z) decay products are produced with an angular separation $\Delta R \geq 0.8$, and only one b-quark (or light quark) has $p_T > 30$ GeV and $|\eta|<1$. The full simulation signal yield prediction can be reproduced within 5\%. This nonclosure uncertainty is added in quadrature to the statistical uncertainty of each bin.

The efficiency of identifying a LLP decay as a TD-tagged jet in bins of the LLP transverse and longitudinal decay position. The sample used to compute the efficiency contains events with pair production of $\chi$ with a mass of 1250 GeV and a lifetime of 0.5 and 3 m, and considering the combinations of the $\chi$ decay modes considered in this search ($H \tilde{G} \rightarrow b\bar{b} \tilde{G}$ or $Z\tilde{G} \rightarrow q\bar{q} \tilde{G}$). The efficiency is calculated on top of the acceptance definition for the resolved topology with two quarks in acceptance, namely, the H (or Z) decay products are produced with an angular separation $\Delta R \geq 0.8$, and both b-quarks (or light quarks) have $p_T > 30$ GeV and $|\eta|<1$. The full simulation signal yield prediction can be reproduced within 7\%. This nonclosure uncertainty is added in quadrature to the statistical uncertainty of each bin.

The efficiency of identifying a LLP decay as a TD-tagged jet in bins of the LLP transverse and longitudinal decay position. The sample used to compute the efficiency contains events with pair production of $\chi$ with a mass of 1500 GeV and a lifetime of 0.5 and 3 m, and considering the combinations of the $\chi$ decay modes considered in this search ($H \tilde{G} \rightarrow b\bar{b} \tilde{G}$ or $Z\tilde{G} \rightarrow q\bar{q} \tilde{G}$). The efficiency is calculated on top of the acceptance definition for the merged topology, namely, the H (or Z) decay products are produced with an angular separation $\Delta R < 0.8$, and the H (or Z) has $p_T > 30$ GeV and $|\eta|<1$. The full simulation signal yield prediction can be reproduced within 3\%. This nonclosure uncertainty is added in quadrature to the statistical uncertainty of each bin.

The efficiency of identifying a LLP decay as a TD-tagged jet in bins of the LLP transverse and longitudinal decay position. The sample used to compute the efficiency contains events with pair production of $\chi$ with a mass of 1500 GeV and a lifetime of 0.5 and 3 m, and considering the combinations of the $\chi$ decay modes considered in this search ($H \tilde{G} \rightarrow b\bar{b} \tilde{G}$ or $Z\tilde{G} \rightarrow q\bar{q} \tilde{G}$). The efficiency is calculated on top of the acceptance definition for the resolved topology with exactly one quark in acceptance, namely, the H (or Z) decay products are produced with an angular separation $\Delta R \geq 0.8$, and only one b-quark (or light quark) has $p_T > 30$ GeV and $|\eta|<1$. The full simulation signal yield prediction can be reproduced within 5\%. This nonclosure uncertainty is added in quadrature to the statistical uncertainty of each bin.

The efficiency of identifying a LLP decay as a TD-tagged jet in bins of the LLP transverse and longitudinal decay position. The sample used to compute the efficiency contains events with pair production of $\chi$ with a mass of 1500 GeV and a lifetime of 0.5 and 3 m, and considering the combinations of the $\chi$ decay modes considered in this search ($H \tilde{G} \rightarrow b\bar{b} \tilde{G}$ or $Z\tilde{G} \rightarrow q\bar{q} \tilde{G}$). The efficiency is calculated on top of the acceptance definition for the resolved topology with two quarks in acceptance, namely, the H (or Z) decay products are produced with an angular separation $\Delta R \geq 0.8$, and both b-quarks (or light quarks) have $p_T > 30$ GeV and $|\eta|<1$. The full simulation signal yield prediction can be reproduced within 7\%. This nonclosure uncertainty is added in quadrature to the statistical uncertainty of each bin.

The efficiency of identifying a LLP decay as a TD-tagged jet in bins of the LLP transverse and longitudinal decay position. The sample used to compute the efficiency contains events with pair production of $\chi$ with a mass of 1800 GeV and a lifetime of 0.5 and 3 m, and considering the combinations of the $\chi$ decay modes considered in this search ($H \tilde{G} \rightarrow b\bar{b} \tilde{G}$ or $Z\tilde{G} \rightarrow q\bar{q} \tilde{G}$). The efficiency is calculated on top of the acceptance definition for the merged topology, namely, the H (or Z) decay products are produced with an angular separation $\Delta R < 0.8$, and the H (or Z) has $p_T > 30$ GeV and $|\eta|<1$. The full simulation signal yield prediction can be reproduced within 3\%. This nonclosure uncertainty is added in quadrature to the statistical uncertainty of each bin.

The efficiency of identifying a LLP decay as a TD-tagged jet in bins of the LLP transverse and longitudinal decay position. The sample used to compute the efficiency contains events with pair production of $\chi$ with a mass of 1800 GeV and a lifetime of 0.5 and 3 m, and considering the combinations of the $\chi$ decay modes considered in this search ($H \tilde{G} \rightarrow b\bar{b} \tilde{G}$ or $Z\tilde{G} \rightarrow q\bar{q} \tilde{G}$). The efficiency is calculated on top of the acceptance definition for the resolved topology with exactly one quark in acceptance, namely, the H (or Z) decay products are produced with an angular separation $\Delta R \geq 0.8$, and only one b-quark (or light quark) has $p_T > 30$ GeV and $|\eta|<1$. The full simulation signal yield prediction can be reproduced within 5\%. This nonclosure uncertainty is added in quadrature to the statistical uncertainty of each bin.

The efficiency of identifying a LLP decay as a TD-tagged jet in bins of the LLP transverse and longitudinal decay position. The sample used to compute the efficiency contains events with pair production of $\chi$ with a mass of 1800 GeV and a lifetime of 0.5 and 3 m, and considering the combinations of the $\chi$ decay modes considered in this search ($H \tilde{G} \rightarrow b\bar{b} \tilde{G}$ or $Z\tilde{G} \rightarrow q\bar{q} \tilde{G}$). The efficiency is calculated on top of the acceptance definition for the resolved topology with two quarks in acceptance, namely, the H (or Z) decay products are produced with an angular separation $\Delta R \geq 0.8$, and both b-quarks (or light quarks) have $p_T > 30$ GeV and $|\eta|<1$. The full simulation signal yield prediction can be reproduced within 7\%. This nonclosure uncertainty is added in quadrature to the statistical uncertainty of each bin.

Cutflow table for a $\tilde{\chi}_{1}^{0}$ signal sample with a mass of 127 GeV.

Cutflow table for a $\tilde{\chi}_{1}^{0}$ signal sample with a mass of 127 GeV.

Cutflow table for a $\tilde{\chi}_{1}^{0}$ signal sample with a mass of 150 GeV.

Cutflow table for a $\tilde{\chi}_{1}^{0}$ signal sample with a mass of 150 GeV.

Cutflow table for a $\tilde{\chi}_{1}^{0}$ signal sample with a mass of 175 GeV.

Cutflow table for a $\tilde{\chi}_{1}^{0}$ signal sample with a mass of 175 GeV.

Cutflow table for a $\tilde{\chi}_{1}^{0}$ signal sample with a mass of 200 GeV.

Cutflow table for a $\tilde{\chi}_{1}^{0}$ signal sample with a mass of 200 GeV.

Cutflow table for a $\tilde{\chi}_{1}^{0}$ signal sample with a mass of 250 GeV.

Cutflow table for a $\tilde{\chi}_{1}^{0}$ signal sample with a mass of 250 GeV.

Cutflow table for a $\tilde{\chi}_{1}^{0}$ signal sample with a mass of 300 GeV.

Cutflow table for a $\tilde{\chi}_{1}^{0}$ signal sample with a mass of 300 GeV.

Cutflow table for a $\tilde{\chi}_{1}^{0}$ signal sample with a mass of 400 GeV.