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

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

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

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

Truth-level signal acceptances for each production process ($\tilde{\chi}_1^\pm \tilde{\chi}_1^0$, $\tilde{\chi}_1^\pm \tilde{\chi}_2^0$, $\tilde{\chi}_1^+ \tilde{\chi}_1^-$, and $\tilde{\chi}_2^0 \tilde{\chi}_1^0$) in a SR with the $S(d_0)$ requirement removed. The acceptance is defined as the fraction of accepted events divided by the total number of events in the generator-level signal Monte Carlo simulation, where the signal candidate track is identified as the charged particle with the largest distance between the interaction vertex and the secondary vertex of the higgsino decays.

Truth-level signal acceptances for each production process ($\tilde{\chi}_1^\pm \tilde{\chi}_1^0$, $\tilde{\chi}_1^\pm \tilde{\chi}_2^0$, $\tilde{\chi}_1^+ \tilde{\chi}_1^-$, and $\tilde{\chi}_2^0 \tilde{\chi}_1^0$) in a SR with the $S(d_0)$ requirement removed. The acceptance is defined as the fraction of accepted events divided by the total number of events in the generator-level signal Monte Carlo simulation, where the signal candidate track is identified as the charged particle with the largest distance between the interaction vertex and the secondary vertex of the higgsino decays.

Truth-level signal acceptances for each production process ($\tilde{\chi}_1^\pm \tilde{\chi}_1^0$, $\tilde{\chi}_1^\pm \tilde{\chi}_2^0$, $\tilde{\chi}_1^+ \tilde{\chi}_1^-$, and $\tilde{\chi}_2^0 \tilde{\chi}_1^0$) in a SR with the $S(d_0)$ requirement removed. The acceptance is defined as the fraction of accepted events divided by the total number of events in the generator-level signal Monte Carlo simulation, where the signal candidate track is identified as the charged particle with the largest distance between the interaction vertex and the secondary vertex of the higgsino decays.

Truth-level signal acceptances for each production process ($\tilde{\chi}_1^\pm \tilde{\chi}_1^0$, $\tilde{\chi}_1^\pm \tilde{\chi}_2^0$, $\tilde{\chi}_1^+ \tilde{\chi}_1^-$, and $\tilde{\chi}_2^0 \tilde{\chi}_1^0$) in a SR with the $S(d_0)$ requirement removed. The acceptance is defined as the fraction of accepted events divided by the total number of events in the generator-level signal Monte Carlo simulation, where the signal candidate track is identified as the charged particle with the largest distance between the interaction vertex and the secondary vertex of the higgsino decays.

Truth-level signal acceptances for each production process ($\tilde{\chi}_1^\pm \tilde{\chi}_1^0$, $\tilde{\chi}_1^\pm \tilde{\chi}_2^0$, $\tilde{\chi}_1^+ \tilde{\chi}_1^-$, and $\tilde{\chi}_2^0 \tilde{\chi}_1^0$) in a SR with the $S(d_0)$ requirement removed. The acceptance is defined as the fraction of accepted events divided by the total number of events in the generator-level signal Monte Carlo simulation, where the signal candidate track is identified as the charged particle with the largest distance between the interaction vertex and the secondary vertex of the higgsino decays.

Truth-level signal acceptances for each production process ($\tilde{\chi}_1^\pm \tilde{\chi}_1^0$, $\tilde{\chi}_1^\pm \tilde{\chi}_2^0$, $\tilde{\chi}_1^+ \tilde{\chi}_1^-$, and $\tilde{\chi}_2^0 \tilde{\chi}_1^0$) in a SR with the $S(d_0)$ requirement removed. The acceptance is defined as the fraction of accepted events divided by the total number of events in the generator-level signal Monte Carlo simulation, where the signal candidate track is identified as the charged particle with the largest distance between the interaction vertex and the secondary vertex of the higgsino decays.

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

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

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

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

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

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

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

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

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

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

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

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

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

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

Genuine correlation function for $D^{+}K^{+}$ in high-multiplicity pp collisions at $\sqrt{s}=13$ TeV.

Genuine correlation function for $D^{*+}\uppi^{-}$ in high-multiplicity pp collisions at $\sqrt{s}=13$ TeV.

Genuine correlation function for $D^{*+}\uppi^{+}$ in high-multiplicity pp collisions at $\sqrt{s}=13$ TeV.

Genuine correlation function for $D^{*+}K^{-}$ in high-multiplicity pp collisions at $\sqrt{s}=13$ TeV.

Genuine correlation function for $D^{*+}K^{+}$ in high-multiplicity pp collisions at $\sqrt{s}=13$ TeV.

The JEC, JER, and total systematic uncertainties in unfolded cross sections as functions of transverse momentum, for 1.5<|y|<2. The total systematic uncertainty includes also the luminosity, jet identification and trigger efficiency uncertainties.

The unfolded measured particle-level inclusive jet cross section as functions of jet pT (markers), for |y|<0.5, compared to the NLO perturbative QCD prediction (histogram), using the CT14NLO PDF set, with muR = muF = HT, and corrected for the NP effects. The experimental and theoretical systematic uncertainties are shown.

The unfolded measured particle-level inclusive jet cross section as functions of jet pT (markers), for 0.5<|y|<1, compared to the NLO perturbative QCD prediction (histogram), using the CT14NLO PDF set, with muR = muF = HT, and corrected for the NP effects. The experimental and theoretical systematic uncertainties are shown.

The unfolded measured particle-level inclusive jet cross section as functions of jet pT (markers), for 1<|y|<1.5, compared to the NLO perturbative QCD prediction (histogram), using the CT14NLO PDF set, with muR = muF = HT, and corrected for the NP effects. The experimental and theoretical systematic uncertainties are shown.

The unfolded measured particle-level inclusive jet cross section as functions of jet pT (markers), for 1.5<|y|<2, compared to the NLO perturbative QCD prediction (histogram), using the CT14NLO PDF set, with muR = muF = HT, and corrected for the NP effects. The experimental and theoretical systematic uncertainties are shown.

Ratios (points) of the unfolded measured cross sections to the NLO theoretical predictions, using the CT14NLO PDF set, with mu = pT, for |y|<0.5. The vertical error bars show the statistical experimental uncertainty. The systematic experimental uncertainty, the total theoretical uncertainty, and the individual sources of theoretical uncertainty are shown.

Ratios (points) of the unfolded measured cross sections to the NLO theoretical predictions, using the CT14NLO PDF set, with mu = pT, for 0.5<|y|<1. The vertical error bars show the statistical experimental uncertainty. The systematic experimental uncertainty, the total theoretical uncertainty, and the individual sources of theoretical uncertainty are shown.

Ratios (points) of the unfolded measured cross sections to the NLO theoretical predictions, using the CT14NLO PDF set, with mu = pT, for 1<|y|<1.5. The vertical error bars show the statistical experimental uncertainty. The systematic experimental uncertainty, the total theoretical uncertainty, and the individual sources of theoretical uncertainty are shown.

Ratios (points) of the unfolded measured cross sections to the NLO theoretical predictions, using the CT14NLO PDF set, with mu = pT, for 1.5<|y|<2. The vertical error bars show the statistical experimental uncertainty. The systematic experimental uncertainty, the total theoretical uncertainty, and the individual sources of theoretical uncertainty are shown.

Ratios (points) of the unfolded measured cross sections to the NLO theoretical predictions, using the CT14NLO PDF set, with mu = HT, for |y|<0.5. The vertical error bars show the statistical experimental uncertainty. The systematic experimental uncertainty, the total theoretical uncertainty, and the individual sources of theoretical uncertainty are shown.

Ratios (points) of the unfolded measured cross sections to the NLO theoretical predictions, using the CT14NLO PDF set, with mu = HT, for 0.5<|y|<1. The vertical error bars show the statistical experimental uncertainty. The systematic experimental uncertainty, the total theoretical uncertainty, and the individual sources of theoretical uncertainty are shown.

Ratios (points) of the unfolded measured cross sections to the NLO theoretical predictions, using the CT14NLO PDF set, with mu = HT, for 1<|y|<1.5. The vertical error bars show the statistical experimental uncertainty. The systematic experimental uncertainty, the total theoretical uncertainty, and the individual sources of theoretical uncertainty are shown.

Ratios (points) of the unfolded measured cross sections to the NLO theoretical predictions, using the CT14NLO PDF set, with mu = HT, for 1.5<|y|<2. The vertical error bars show the statistical experimental uncertainty. The systematic experimental uncertainty, the total theoretical uncertainty, and the individual sources of theoretical uncertainty are shown.

Ratios (points) of the unfolded measured cross sections to the NNLO theoretical predictions, using the CT14NNLO PDF set, with mu = HT, for |y|<0.5. The vertical error bars show the statistical experimental uncertainty. The systematic experimental uncertainty, the total theoretical uncertainty, and the individual sources of theoretical uncertainty are shown.

Ratios (points) of the unfolded measured cross sections to the NNLO theoretical predictions, using the CT14NNLO PDF set, with mu = HT, for 0.5<|y|<1. The vertical error bars show the statistical experimental uncertainty. The systematic experimental uncertainty, the total theoretical uncertainty, and the individual sources of theoretical uncertainty are shown.

Ratios (points) of the unfolded measured cross sections to the NNLO theoretical predictions, using the CT14NNLO PDF set, with mu = HT, for 1<|y|<1.5. The vertical error bars show the statistical experimental uncertainty. The systematic experimental uncertainty, the total theoretical uncertainty, and the individual sources of theoretical uncertainty are shown.

Ratios (points) of the unfolded measured cross sections to the NNLO theoretical predictions, using the CT14NNLO PDF set, with mu = HT, for 1.5<|y|<2. The vertical error bars show the statistical experimental uncertainty. The systematic experimental uncertainty, the total theoretical uncertainty, and the individual sources of theoretical uncertainty are shown.

Ratios (points) of the unfolded measured cross sections to the NNLO theoretical predictions, using the NNPDF31NNLO PDF set, with mu = HT, for |y|<0.5. The vertical error bars show the statistical experimental uncertainty. The systematic experimental uncertainty, the total theoretical uncertainty, and the individual sources of theoretical uncertainty are shown.

Ratios (points) of the unfolded measured cross sections to the NNLO theoretical predictions, using the NNPDF31NNLO PDF set, with mu = HT, for 0.5<|y|<1. The vertical error bars show the statistical experimental uncertainty. The systematic experimental uncertainty, the total theoretical uncertainty, and the individual sources of theoretical uncertainty are shown.

Ratios (points) of the unfolded measured cross sections to the NNLO theoretical predictions, using the NNPDF31NNLO PDF set, with mu = HT, for 1<|y|<1.5. The vertical error bars show the statistical experimental uncertainty. The systematic experimental uncertainty, the total theoretical uncertainty, and the individual sources of theoretical uncertainty are shown.

Ratios (points) of the unfolded measured cross sections to the NNLO theoretical predictions, using the NNPDF31NNLO PDF set, with mu = HT, for 1.5<|y|<2. The vertical error bars show the statistical experimental uncertainty. The systematic experimental uncertainty, the total theoretical uncertainty, and the individual sources of theoretical uncertainty are shown.

The effect of aS(MZ) variation, for |y|<0.5. The NNLO theoretical cross section predictions using the NNPDF31NNLO PDF with m = HT, calculated for different choices of aS (0.108, 0.110, 0.112, 0.114, 0.116, 0.117, 0.118, 0.119, 0.120, 0.122, and 0.124), are divided by the benchmark NNLO prediction for aS = 0.118 and the same choice of PDF set, muR, and muF. Also shown is the experimental unfolded measurement divided by the same benchmark prediction. The width of the unity line corresponds to the statistical uncertainty from the MC integration for the determination of the NNLO prediction. The error bars on the unfolded data correspond to the total experimental statistical and systematic uncertainty added in quadrature.

The effect of aS(MZ) variation, for 0.5<|y|<1. The NNLO theoretical cross section predictions using the NNPDF31NNLO PDF with m = HT, calculated for different choices of aS (0.108, 0.110, 0.112, 0.114, 0.116, 0.117, 0.118, 0.119, 0.120, 0.122, and 0.124), are divided by the benchmark NNLO prediction for aS = 0.118 and the same choice of PDF set, muR, and muF. Also shown is the experimental unfolded measurement divided by the same benchmark prediction. The width of the unity line corresponds to the statistical uncertainty from the MC integration for the determination of the NNLO prediction. The error bars on the unfolded data correspond to the total experimental statistical and systematic uncertainty added in quadrature.

The effect of aS(MZ) variation, for 1<|y|<1.5. The NNLO theoretical cross section predictions using the NNPDF31NNLO PDF with m = HT, calculated for different choices of aS (0.108, 0.110, 0.112, 0.114, 0.116, 0.117, 0.118, 0.119, 0.120, 0.122, and 0.124), are divided by the benchmark NNLO prediction for aS = 0.118 and the same choice of PDF set, muR, and muF. Also shown is the experimental unfolded measurement divided by the same benchmark prediction. The width of the unity line corresponds to the statistical uncertainty from the MC integration for the determination of the NNLO prediction. The error bars on the unfolded data correspond to the total experimental statistical and systematic uncertainty added in quadrature.

The effect of aS(MZ) variation, for 1.5<|y|<2. The NNLO theoretical cross section predictions using the NNPDF31NNLO PDF with m = HT, calculated for different choices of aS (0.108, 0.110, 0.112, 0.114, 0.116, 0.117, 0.118, 0.119, 0.120, 0.122, and 0.124), are divided by the benchmark NNLO prediction for aS = 0.118 and the same choice of PDF set, muR, and muF. Also shown is the experimental unfolded measurement divided by the same benchmark prediction. The width of the unity line corresponds to the statistical uncertainty from the MC integration for the determination of the NNLO prediction. The error bars on the unfolded data correspond to the total experimental statistical and systematic uncertainty added in quadrature.

Unfolded correlation matrix for |y|<0.5.

Unfolded correlation matrix for 0.5<|y|<1.

Unfolded correlation matrix for 1<|y|<1.5.

Unfolded correlation matrix for 1.5<|y|<2.

Data from Figure 2, panel b, U+U, 0-0.5% Centrality, 0.2<p_{T}<3 GeV/c

Data from Figure 3, panel a, 0.2<p_{T}<3 GeV/c

Data from Figure 3, panel b, 0.2<p_{T}<3 GeV/c

Data from Figure 3, panel c, 0.2<p_{T}<3 GeV/c

Data from Figure 3, panel d, 0-5% Centrality, 0.2<p_{T}<3 GeV/c

Data from Figure 3, panel e, 0-5% Centrality, 0.2<p_{T}<3 GeV/c

Data from Figure 3, panel f, 0-5% Centrality, 0.2<p_{T}<3 GeV/c

Data from Figure 4, panel a, Au+Au, 2subevent method, 0.2<p_{T}<3 GeV/c

Data from Figure 4, panel a, U+U, 2subevent method, 0.2<p_{T}<3 GeV/c

Data from Figure 4, panel b, Au+Au, 0.2<p_{T}<3 GeV/c

Data from Figure 4, panel b, U+U, 0.2<p_{T}<3 GeV/c

Data from Figure 4, panel c, Au+Au, 0.2<p_{T}<3 GeV/c

Data from Figure 4, panel c, U+U, 0.2<p_{T}<3 GeV/c

Data from Figure 4, panel d, Au+Au, 0.2<p_{T}<3 GeV/c

Data from Figure 4, panel d, U+U, 0.2<p_{T}<3 GeV/c

Data from Figure 5, panel a, Au+Au, 2subevent method, 0.2<p_{T}<3 GeV/c

Data from Figure 5, panel a, Au+Au, standard method, 0.2<p_{T}<3 GeV/c

Data from Figure 5, panel a, Au+Au, difference, 0.2<p_{T}<3 GeV/c

Data from Figure 5, panel b, Au+Au, 2subevent method, 0.2<p_{T}<3 GeV/c

Data from Figure 5, panel b, Au+Au, standard, 0.2<p_{T}<3 GeV/c

Data from Figure 5, panel b, Au+Au, difference, 0.2<p_{T}<3 GeV/c

Data from Figure 5, panel c, Au+Au, 2subevent method, 0.2<p_{T}<3 GeV/c

Data from Figure 5, panel c, Au+Au, standard, 0.2<p_{T}<3 GeV/c

Data from Figure 5, panel c, Au+Au, difference, 0.2<p_{T}<3 GeV/c