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.

Definition of the tetralepton signal regions.

Definition of the fiducial volumes at particle- and parton-level. Leptons refer exclusively to electrons and muons - they are dressed with additional radiation at particle-level, but not at parton-level.

Definition of the dilepton $t\bar{t}$ validation regions.

Pre-fit distribution of the number of $b$-jets in 2L-$e\mu$-6j2b, this distribution is not used in the fit.

Pre-fit distribution of the DNN output 2L-$e\mu$-6j1b, this distribution is not used in the fit.

Pre-fit distribution of the DNN output 2L-$e\mu$-5j2b, this distribution is not used in the fit.

Pre-fit distribution of the DNN output 2L-$e\mu$-6j2b, this distribution is not used in the fit.

Definition of the tetralepton control region.

Definition of the trilepton fakes control regions.

Pre-fit distribution of jet multiplicity in CR-$t\bar{t}$-e region.

Pre-fit distribution of loose lepton transverse momentum in CR-$t\bar{t}$-$\mu$ region.

Pre-fit distribution of the transverse mass of the trailing lepton and the missing transverse momentum in CR-Z-e region.

Post-fit distribution of jet multiplicity in CR-$t\bar{t}$-e region

Post-fit distribution of loose lepton transverse momentum in CR-$t\bar{t}$-$\mu$ region

Post-fit distribution of the transverse mass of the trailing lepton and the missing transverse momentum in CR-Z-e region

Post-fit distribution of NN output in SR-2L-5j2b region.

Post-fit distribution of NN output in SR-2L-6j1b region.

Post-fit distribution of NN output in SR-2L-6j2b region.

Post-fit distribution of DNN-$t\bar{t}Z$ output in 3L-SR-ttZ region.

Post-fit distribution of DNN-$t\bar{t}Z$ outputt in 3L-SR-tZq region.

Post fit events yields in 3L-SR-WZ region.

Post-fit distribution of NN output in 4L-SR-SF region.

Post-fit distribution of NN output in 4L-SR-DF region.

Post-fit distribution of b-tagger output for leading b-jet in 4L-CR-ZZ region.

Measured values of the background normalizations obtained from the combined fit. The uncertainties include statistical and systematic sources.

Measured $\sigma_{t\bar{t}\text{Z}}$ cross sections obtained from the fits in the different lepton channels. The uncertainties include statistical and systematic sources.

Grouped impact of systematic uncertainties in the combined inclusive fit to data.

Unfolded absolute cross section as a function of $p^{Z}_{T}$ in the combination of $3\ell$ and $4\ell$ channels at particle-level (Figure 8 top-left).

Unfolded absolute cross section as a function of $p^{Z}_{T}$ in the combination of $3\ell$ and $4\ell$ channels at parton-level (Figure 8 top-right).

Unfolded normalized cross section as a function of $p^{Z}_{T}$ in the combination of $3\ell$ and $4\ell$ channels at particle-level (Figure 8 bottom-left).

Unfolded normalized cross section as a function of $p^{Z}_{T}$ in the combination of $3\ell$ and $4\ell$ channels at parton-level (Figure 8 bottom-right).

Unfolded absolute cross section as a function of $|y^{Z}$| in the combination of $3\ell$ and $4\ell$ channels at particle-level (Figure 17 top-left and Figure 11 top-left).

Unfolded absolute cross section as a function of $|y^{Z}$| in the combination of $3\ell$ and $4\ell$ channels at parton-level (Figure 17 top-right).

Unfolded normalized cross section as a function of $|y^{Z}$| in the combination of $3\ell$ and $4\ell$ channels at particle-level (Figure 17 bottom-left).

Unfolded normalized cross section as a function of $|y^{Z}$| in the combination of $3\ell$ and $4\ell$ channels at parton-level (Figure 17 bottom-right).

Unfolded absolute cross section as a function of cos $\theta_{Z}^{*}$ in the combination of $3\ell$ and $4\ell$ channels at particle-level (Figure 18 top-left and Figure 11 top-right).

Unfolded absolute cross section as a function of cos $\theta_{Z}^{*}$ in the combination of $3\ell$ and $4\ell$ channels at parton-level (Figure 18 top-right).

Unfolded normalized cross section as a function of cos $\theta_{Z}^{*}$ in the combination of $3\ell$ and $4\ell$ channels at particle-level (Figure 18 bottom-left).

Unfolded normalized cross section as a function of cos $\theta_{Z}^{*}$ in the combination of $3\ell$ and $4\ell$ channels at parton-level (Figure 18 bottom-right).

Unfolded absolute cross section as a function of $p_{T}^{\mathrm{top}}$ in the combination of $3\ell$ and $4\ell$ channels at particle-level (Figure 19 top-left and Figure 11 bottom-left).

Unfolded absolute cross section as a function of $p_{T}^{\mathrm{top}}$ in the combination of $3\ell$ and $4\ell$ channels at parton-level (Figure 19, top-right).

Unfolded normalized cross section as a function of $p_{T}^{\mathrm{top}}$ in the combination of $3\ell$ and $4\ell$ channels at particle-level (Figure 19, bottom-left).

Unfolded normalized cross section as a function of $p_{T}^{\mathrm{top}}$ in the combination of $3\ell$ and $4\ell$ channels at parton-level (Figure 19, bottom-right).

Unfolded absolute cross section as a function of $p_{T}^{t\bar{t}}$ in the combination of $3\ell$ and $4\ell$ channels at particle-level (Figure 20 top-left and Figure 11 bottom-right).

Unfolded absolute cross section as a function of $p_{T}^{t\bar{t}}$ in the combination of $3\ell$ and $4\ell$ channels at parton-level (Figure 20, top-right).

Unfolded normalized cross section as a function of $p_{T}^{t\bar{t}}$ in the combination of $3\ell$ and $4\ell$ channels at particle-level (Figure 20, bottom-left)

Unfolded normalized cross section as a function of $p_{T}^{t\bar{t}}$ in the combination of $3\ell$ and $4\ell$ channels at parton-level (Figure 20, bottom-right)

Unfolded absolute cross section as a function of $|\Delta\Phi(t\bar{t}, Z)|/\pi$ in the combination of $3\ell$ and $4\ell$ channels at particle-level (Figure 21 top-left and Figure 12 top-left).

Unfolded absolute cross section as a function of $|\Delta\Phi(t\bar{t}, Z)|/\pi$ in the combination of $3\ell$ and $4\ell$ channels at parton-level (Figure 21, top-right).

Unfolded normalized cross section as a function of $|\Delta\Phi(t\bar{t}, Z)|/\pi$ in the combination of $3\ell$ and $4\ell$ channels at particle-level (Figure 21, bottom-left).

Unfolded normalized cross section as a function of $|\Delta\Phi(t\bar{t}, Z)|/\pi$ in the combination of $3\ell$ and $4\ell$ channels at parton-level (Figure 21, top-right).

Unfolded absolute cross section as a function of $m^{t\bar{t}Z}$ in the combination of $3\ell$ and $4\ell$ channels at particle-level (Figure 22 top-left and Figure 12 bottom-left).

Unfolded absolute cross section as a function of $m^{t\bar{t}Z}$ in the combination of $3\ell$ and $4\ell$ channels at parton-level (Figure 22, top-right).

Unfolded normalized cross section as a function of $m^{t\bar{t}Z}$ in the combination of $3\ell$ and $4\ell$ channels at particle-level (Figure 22, bottom-left).

Unfolded normalized cross section as a function of $m^{t\bar{t}Z}$ in the combination of $3\ell$ and $4\ell$ channels at parton-level (Figure 22, bottom-right).

Unfolded absolute cross section as a function of $m^{t\bar{t}}$ in the combination of $3\ell$ and $4\ell$ channels at particle-level (Figure 23 top-left and Figure 12 bottom-right).

Unfolded absolute cross section as a function of $m^{t\bar{t}}$ in the combination of $3\ell$ and $4\ell$ channels at parton-level (Figure 23, top-right).

Unfolded normalized cross section as a function of $m^{t\bar{t}}$ in the combination of $3\ell$ and $4\ell$ channels at particle-level (Figure 23, bottom-left).

Unfolded normalized cross section as a function of $m^{t\bar{t}}$ in the combination of $3\ell$ and $4\ell$ channels at parton-level (Figure 23, bottom-right).

Unfolded absolute cross section as a function of $|y^{t\bar{t}Z}|$ in the combination of $3\ell$ and $4\ell$ channels at particle-level (Figure 24 top-left and Figure 12 top-right).

Unfolded absolute cross section as a function of $|y^{t\bar{t}Z}|$ in the combination of $3\ell$ and $4\ell$ channels at parton-level (Figure 24, top-right).

Unfolded normalized cross section as a function of $|y^{t\bar{t}Z}|$ in the combination of $3\ell$ and $4\ell$ channels at particle-level (Figure 24, bottom-left).

Unfolded normalized cross section as a function of $|y^{t\bar{t}Z}|$ in the combination of $3\ell$ and $4\ell$ channels at parton-level (Figure 24, bottom-right).

Unfolded absolute cross section as a function of $H_{\text{T}}^{\text{l}}$ in the trilepton channel at particle-level (Figure 25 top-left and Figure 9 top-left).

Unfolded absolute cross section as a function of $H_{\text{T}}^{\text{l}}$ in the trilepton channel at parton-level (Figure 25 top-right).

Unfolded normalized cross section as a function of $H_{\text{T}}^{\text{l}}$ in the trilepton channel at particle-level (Figure 25 bottom-left).

Unfolded normalized cross section as a function of $H_{\text{T}}^{\text{l}}$ in the trilepton channel at parton-level (Figure 25 bottom-right).

Unfolded absolute cross section as a function of $|\Delta\Phi(Z, t_{lep})|/\pi$ in the trilepton channel at particle-level (Figure 26 top-left and Figure 10 bottom-left).

Unfolded absolute cross section as a function of $|\Delta\Phi(Z, t_{lep})|/\pi$ in the trilepton channel at parton-level (Figure 26 top-right).

Unfolded normalized cross section as a function of $|\Delta\Phi(Z, t_{lep})|/\pi$ in the trilepton channel at particle-level (Figure 26 bottom-left).

Unfolded normalized cross section as a function of $|\Delta\Phi(Z, t_{lep})|/\pi$ in the trilepton channel at parton-level (Figure 26 bottom-right).

Unfolded absolute cross section as a function of $|\Delta y(Z, t_{lep})|$ in the trilepton channel at particle-level (Figure 27 top-left and Figure 10 bottom-right).

Unfolded absolute cross section as a function of $|\Delta y(Z, t_{lep})|$ in the trilepton channel at parton-level (Figure 27 top-right).

Unfolded normalized cross section as a function of $|\Delta y(Z, t_{lep})|$ in the trilepton channel at particle-level (Figure 27 bottom-left).

Unfolded normalized cross section as a function of $|\Delta y(Z, t_{lep})|$ in the trilepton channel at parton-level (Figure 27 bottom-right).

Unfolded absolute cross section as a function of $p_{\text{T}}^{\ell, non-Z}$ in the trilepton channel at particle-level (Figure 28 top-left and Figure 10 top-left).

Unfolded absolute cross section as a function of $p_{\text{T}}^{\ell, non-Z}$ in the trilepton channel at parton-level (Figure 28 top-right).

Unfolded normalized cross section as a function of $p_{\text{T}}^{\ell, non-Z}$ in the trilepton channel at particle-level (Figure 28 bottom-left).

Unfolded normalized cross section as a function of $p_{\text{T}}^{\ell, non-Z}$ in the trilepton channel at parton-level (Figure 28 bottom-right).

Unfolded absolute cross section as a function of $N_{\text{jets}}$ in the trilepton channel at particle-level (Figure 29 left and Figure 9 bottom-left).

Unfolded normalized cross section as a function of $N_{\text{jets}}$ in the trilepton channel at particle-level (Figure 29 right).

Unfolded absolute cross section as a function of $H_{\text{T}}^{\text{l}}$ in the tetralepton channel at particle-level (Figure 30 top-left and Figure 9 top-right).

Unfolded absolute cross section as a function of $H_{\text{T}}^{\text{l}}$ in the tetralepton channel at parton-level (Figure 30 top-right).

Unfolded normalized cross section as a function of $H_{\text{T}}^{\text{l}}$ in the tetralepton channel at particle-level (Figure 30 bottom-left).

Unfolded normalized cross section as a function of $H_{\text{T}}^{\text{l}}$ in the tetralepton channel at parton-level (Figure 30 bottom-right).

Unfolded absolute cross section as a function of $|\Delta\Phi(l_{t}^{+}, l_{\bar{t}}^{-})|/\pi$ in the tetralepton channel at particle-level (Figure 31 top-left and Figure 10 top-right).

Unfolded absolute cross section as a function of $|\Delta\Phi(l_{t}^{+}, l_{\bar{t}}^{-})|/\pi$ in the tetralepton channel at parton-level (Figure 31 top-right).

Unfolded normalized cross section as a function of $|\Delta\Phi(l_{t}^{+}, l_{\bar{t}}^{-})|/\pi$ in the tetralepton channel at particle-level (Figure 31 bottom-left).

Unfolded normalized cross section as a function of $|\Delta\Phi(l_{t}^{+}, l_{\bar{t}}^{-})|/\pi$ in the tetralepton channel at parton-level (Figure 31 bottom-right).

Unfolded absolute cross section as a function of $N_{\text{jets}}$ in the tetralepton channel at particle-level (Figure 32 left and Figure 9 bottom-right).

Unfolded normalized cross section as a function of $N_{\text{jets}}$ in the tetralepton channel at particle-level (Figure 32 right).

Bootstrap replicas (0-499) for data in all regions used in inclusive cross section measurement. The used bootstrap method is described in ATL-PHYS-PUB-2021-011 (https://atlas.web.cern.ch/Atlas/GROUPS/PHYSICS/PUBNOTES/ATL-PHYS-PUB-2021-011/).

Bootstrap replicas (500-999) for data in all regions used in inclusive cross section measurement. The used bootstrap method is described in ATL-PHYS-PUB-2021-011 (https://atlas.web.cern.ch/Atlas/GROUPS/PHYSICS/PUBNOTES/ATL-PHYS-PUB-2021-011/).

Bootstrap replicas (0-499) for data, variable $|\Delta\Phi(t\bar{t}, Z)|/\pi$. The used bootstrap method is described in ATL-PHYS-PUB-2021-011 (https://atlas.web.cern.ch/Atlas/GROUPS/PHYSICS/PUBNOTES/ATL-PHYS-PUB-2021-011/).

Bootstrap replicas (500-999) for data, variable $|\Delta\Phi(t\bar{t}, Z)|/\pi$. The used bootstrap method is described in ATL-PHYS-PUB-2021-011 (https://atlas.web.cern.ch/Atlas/GROUPS/PHYSICS/PUBNOTES/ATL-PHYS-PUB-2021-011/).

Bootstrap replicas (0-499) for data, variable $|\Delta\Phi(Z, t_{lep})|/\pi$. The used bootstrap method is described in ATL-PHYS-PUB-2021-011 (https://atlas.web.cern.ch/Atlas/GROUPS/PHYSICS/PUBNOTES/ATL-PHYS-PUB-2021-011/).

Bootstrap replicas (500-999) for data, variable $|\Delta\Phi(Z, t_{lep})|/\pi$. The used bootstrap method is described in ATL-PHYS-PUB-2021-011 (https://atlas.web.cern.ch/Atlas/GROUPS/PHYSICS/PUBNOTES/ATL-PHYS-PUB-2021-011/).

Bootstrap replicas (0-499) for data, variable $m^{t\bar{t}}$. The used bootstrap method is described in ATL-PHYS-PUB-2021-011 (https://atlas.web.cern.ch/Atlas/GROUPS/PHYSICS/PUBNOTES/ATL-PHYS-PUB-2021-011/).

Bootstrap replicas (500-999) for data, variable $m^{t\bar{t}}$. The used bootstrap method is described in ATL-PHYS-PUB-2021-011 (https://atlas.web.cern.ch/Atlas/GROUPS/PHYSICS/PUBNOTES/ATL-PHYS-PUB-2021-011/).

Bootstrap replicas (0-499) for data, variable $N_{\text{jets}}$ in $3\ell$ channel. The used bootstrap method is described in ATL-PHYS-PUB-2021-011 (https://atlas.web.cern.ch/Atlas/GROUPS/PHYSICS/PUBNOTES/ATL-PHYS-PUB-2021-011/).

Bootstrap replicas (500-999) for data, variable $N_{\text{jets}}$ in $3\ell$ channel. The used bootstrap method is described in ATL-PHYS-PUB-2021-011 (https://atlas.web.cern.ch/Atlas/GROUPS/PHYSICS/PUBNOTES/ATL-PHYS-PUB-2021-011/).

Bootstrap replicas (0-499) for data, variable $|y^{t\bar{t}Z}|$. The used bootstrap method is described in ATL-PHYS-PUB-2021-011 (https://atlas.web.cern.ch/Atlas/GROUPS/PHYSICS/PUBNOTES/ATL-PHYS-PUB-2021-011/).

Bootstrap replicas (500-999) for data, variable $|y^{t\bar{t}Z}|$. The used bootstrap method is described in ATL-PHYS-PUB-2021-011 (https://atlas.web.cern.ch/Atlas/GROUPS/PHYSICS/PUBNOTES/ATL-PHYS-PUB-2021-011/).

Bootstrap replicas (0-499) for data, variable $H_{\text{T}}^{\text{l}}$ in $3\ell$ channel. The used bootstrap method is described in ATL-PHYS-PUB-2021-011 (https://atlas.web.cern.ch/Atlas/GROUPS/PHYSICS/PUBNOTES/ATL-PHYS-PUB-2021-011/).

Bootstrap replicas (500-999) for data, variable $H_{\text{T}}^{\text{l}}$ in $3\ell$ channel. The used bootstrap method is described in ATL-PHYS-PUB-2021-011 (https://atlas.web.cern.ch/Atlas/GROUPS/PHYSICS/PUBNOTES/ATL-PHYS-PUB-2021-011/).

Bootstrap replicas (0-499) for data, variable $y^{Z}$. The used bootstrap method is described in ATL-PHYS-PUB-2021-011 (https://atlas.web.cern.ch/Atlas/GROUPS/PHYSICS/PUBNOTES/ATL-PHYS-PUB-2021-011/).

Bootstrap replicas (500-999) for data, variable $y^{Z}$. The used bootstrap method is described in ATL-PHYS-PUB-2021-011 (https://atlas.web.cern.ch/Atlas/GROUPS/PHYSICS/PUBNOTES/ATL-PHYS-PUB-2021-011/).

Bootstrap replicas (0-499) for data, variable $p_{T}^{\mathrm{top}}$. The used bootstrap method is described in ATL-PHYS-PUB-2021-011 (https://atlas.web.cern.ch/Atlas/GROUPS/PHYSICS/PUBNOTES/ATL-PHYS-PUB-2021-011/).

Bootstrap replicas (500-999) for data, variable $p_{T}^{\mathrm{top}}$. The used bootstrap method is described in ATL-PHYS-PUB-2021-011 (https://atlas.web.cern.ch/Atlas/GROUPS/PHYSICS/PUBNOTES/ATL-PHYS-PUB-2021-011/).

Bootstrap replicas (0-499) for data, variable cos $\theta^{*}_{Z}$. The used bootstrap method is described in ATL-PHYS-PUB-2021-011 (https://atlas.web.cern.ch/Atlas/GROUPS/PHYSICS/PUBNOTES/ATL-PHYS-PUB-2021-011/).

Bootstrap replicas (500-999) for data, variable cos $\theta^{*}_{Z}$. The used bootstrap method is described in ATL-PHYS-PUB-2021-011 (https://atlas.web.cern.ch/Atlas/GROUPS/PHYSICS/PUBNOTES/ATL-PHYS-PUB-2021-011/).

Bootstrap replicas (0-499) for data, variable $p_{\text{T}}^{\ell, non-Z}$. The used bootstrap method is described in ATL-PHYS-PUB-2021-011 (https://atlas.web.cern.ch/Atlas/GROUPS/PHYSICS/PUBNOTES/ATL-PHYS-PUB-2021-011/).

Bootstrap replicas (500-999) for data, variable $p_{\text{T}}^{\ell, non-Z}$. The used bootstrap method is described in ATL-PHYS-PUB-2021-011 (https://atlas.web.cern.ch/Atlas/GROUPS/PHYSICS/PUBNOTES/ATL-PHYS-PUB-2021-011/).

Bootstrap replicas (0-499) for data, variable $|\Delta\Phi(l_{t}^{+}, l_{\bar{t}}^{-})|/\pi$. The used bootstrap method is described in ATL-PHYS-PUB-2021-011 (https://atlas.web.cern.ch/Atlas/GROUPS/PHYSICS/PUBNOTES/ATL-PHYS-PUB-2021-011/).

Bootstrap replicas (500-999) for data, variable $|\Delta\Phi(l_{t}^{+}, l_{\bar{t}}^{-})|/\pi$. The used bootstrap method is described in ATL-PHYS-PUB-2021-011 (https://atlas.web.cern.ch/Atlas/GROUPS/PHYSICS/PUBNOTES/ATL-PHYS-PUB-2021-011/).

Bootstrap replicas (0-499) for data, variable $H_{\text{T}}^{\text{l}}$ in $4\ell$ channel. The used bootstrap method is described in ATL-PHYS-PUB-2021-011 (https://atlas.web.cern.ch/Atlas/GROUPS/PHYSICS/PUBNOTES/ATL-PHYS-PUB-2021-011/).

Bootstrap replicas (500-999) for data, variable $H_{\text{T}}^{\text{l}}$ in $4\ell$ channel. The used bootstrap method is described in ATL-PHYS-PUB-2021-011 (https://atlas.web.cern.ch/Atlas/GROUPS/PHYSICS/PUBNOTES/ATL-PHYS-PUB-2021-011/).

Bootstrap replicas (0-499) for data, variable $m^{t\bar{t}Z}$. The used bootstrap method is described in ATL-PHYS-PUB-2021-011 (https://atlas.web.cern.ch/Atlas/GROUPS/PHYSICS/PUBNOTES/ATL-PHYS-PUB-2021-011/).

Bootstrap replicas (500-999) for data, variable $m^{t\bar{t}Z}$. The used bootstrap method is described in ATL-PHYS-PUB-2021-011 (https://atlas.web.cern.ch/Atlas/GROUPS/PHYSICS/PUBNOTES/ATL-PHYS-PUB-2021-011/).

Bootstrap replicas (0-499) for data, variable $N_{\text{jets}}$ in $4\ell$ channel. The used bootstrap method is described in ATL-PHYS-PUB-2021-011 (https://atlas.web.cern.ch/Atlas/GROUPS/PHYSICS/PUBNOTES/ATL-PHYS-PUB-2021-011/).

Bootstrap replicas (500-999) for data, variable $N_{\text{jets}}$ in $4\ell$ channel. The used bootstrap method is described in ATL-PHYS-PUB-2021-011 (https://atlas.web.cern.ch/Atlas/GROUPS/PHYSICS/PUBNOTES/ATL-PHYS-PUB-2021-011/).

Bootstrap replicas (0-499) for data, variable $|\Delta y(Z, t_{lep})|$. The used bootstrap method is described in ATL-PHYS-PUB-2021-011 (https://atlas.web.cern.ch/Atlas/GROUPS/PHYSICS/PUBNOTES/ATL-PHYS-PUB-2021-011/).

Bootstrap replicas (500-999) for data, variable $|\Delta y(Z, t_{lep})|$. The used bootstrap method is described in ATL-PHYS-PUB-2021-011 (https://atlas.web.cern.ch/Atlas/GROUPS/PHYSICS/PUBNOTES/ATL-PHYS-PUB-2021-011/).

Bootstrap replicas (0-499) for data, variable $p^{Z}_{T}$. The used bootstrap method is described in ATL-PHYS-PUB-2021-011 (https://atlas.web.cern.ch/Atlas/GROUPS/PHYSICS/PUBNOTES/ATL-PHYS-PUB-2021-011/).

Bootstrap replicas (500-999) for data, variable $p^{Z}_{T}$. The used bootstrap method is described in ATL-PHYS-PUB-2021-011 (https://atlas.web.cern.ch/Atlas/GROUPS/PHYSICS/PUBNOTES/ATL-PHYS-PUB-2021-011/).

Bootstrap replicas (0-499) for data, variable $p_{T}^{t\bar{t}}$. The used bootstrap method is described in ATL-PHYS-PUB-2021-011 (https://atlas.web.cern.ch/Atlas/GROUPS/PHYSICS/PUBNOTES/ATL-PHYS-PUB-2021-011/).

Bootstrap replicas (500-999) for data, variable $p_{T}^{t\bar{t}}$. The used bootstrap method is described in ATL-PHYS-PUB-2021-011 (https://atlas.web.cern.ch/Atlas/GROUPS/PHYSICS/PUBNOTES/ATL-PHYS-PUB-2021-011/).

Parton-level acceptance and selection efficiency histograms for $|\Delta\Phi(Z, t_{lep})|/\pi$ variable.

Parton-level acceptance and selection efficiency histograms for $|\Delta y(Z, t_{lep})|$ variable.

Parton-level acceptance and selection efficiency histograms for $H_{\text{T}}^{\text{ l}}$ variable.

Parton-level acceptance and selection efficiency histograms for $p_{\text{T}}^{\ell, non-Z}$ variable.

Parton-level acceptance and selection efficiency histograms for $|\Delta\Phi(l_{t}^{+}, l_{\bar{t}}^{-})|/\pi$ variable.

Parton-level acceptance and selection efficiency histograms for $H_{\text{T}}^{\text{ l}}$ variable.

Parton-level acceptance and selection efficiency histograms for cos $\theta_{Z}^{*}$ variable.

Parton-level acceptance and selection efficiency histograms for $p^{Z}_{T}$ variable.

Parton-level acceptance and selection efficiency histograms for $|y^{Z}$| variable.

Parton-level acceptance and selection efficiency histograms for $|\Delta\Phi(t\bar{t}, Z)|/\pi$ variable.

Parton-level acceptance and selection efficiency histograms for $m^{t\bar{t}}$ variable.

Parton-level acceptance and selection efficiency histograms for $m^{t\bar{t}Z}$ variable.

Parton-level acceptance and selection efficiency histograms for $p_{T}^{\mathrm{top}}$ variable.

Parton-level acceptance and selection efficiency histograms for $p_{T}^{t\bar{t}}$ variable.

Parton-level acceptance and selection efficiency histograms for $|y^{t\bar{t}Z}|$ variable.

Particle-level acceptance and selection efficiency histograms for $|\Delta\Phi(Z, t_{lep})|/\pi$ variable.

Particle-level acceptance and selection efficiency histograms for $|\Delta y(Z, t_{lep})|$ variable.

Particle-level acceptance and selection efficiency histograms for $H_{\text{T}}^{\text{ l}}$ variable.

Particle-level acceptance and selection efficiency histograms for $N_{\text{jets}}$ variable.

Particle-level acceptance and selection efficiency histograms for $p_{\text{T}}^{\ell, non-Z}$ variable.

Particle-level acceptance and selection efficiency histograms for $|\Delta\Phi(l_{t}^{+}, l_{\bar{t}}^{-})|/\pi$ variable.

Particle-level acceptance and selection efficiency histograms for $H_{\text{T}}^{\text{ l}}$ variable.

Particle-level acceptance and selection efficiency histograms for $N_{\text{jets}}$ variable.

Particle-level acceptance and selection efficiency histograms for cos $\theta_{Z}^{*}$ variable.

Particle-level acceptance and selection efficiency histograms for $p^{Z}_{T}$ variable.

Particle-level acceptance and selection efficiency histograms for $|y^{Z}$| variable.

Particle-level acceptance and selection efficiency histograms for $|\Delta\Phi(t\bar{t}, Z)|/\pi$ variable.

Particle-level acceptance and selection efficiency histograms for $m^{t\bar{t}}$ variable.

Particle-level acceptance and selection efficiency histograms for $m^{t\bar{t}Z}$ variable.

Particle-level acceptance and selection efficiency histograms for $p_{T}^{\mathrm{top}}$ variable.

Particle-level acceptance and selection efficiency histograms for $p_{T}^{t\bar{t}}$ variable.

Particle-level acceptance and selection efficiency histograms for $|y^{t\bar{t}Z}|$ variable.

Migration matrix for cos $\theta_{Z}^{*}$ variable at particle-level in region SR-3L-ttZ.

Migration matrix for cos $\theta_{Z}^{*}$ variable at particle-level in region SR-3L-tZq.

Migration matrix for cos $\theta_{Z}^{*}$ variable at particle-level in region SR-3L-WZ.

Migration matrix for cos $\theta_{Z}^{*}$ variable at particle-level in region SR-4L-DF.

Migration matrix for cos $\theta_{Z}^{*}$ variable at particle-level in region SR-4L-SF.

Migration matrix for cos $\theta_{Z}^{*}$ variable at particle-level in region CR-4L-ZZ.

Migration matrix for cos $\theta_{Z}^{*}$ variable at parton-level in region SR-3L-ttZ.

Migration matrix for cos $\theta_{Z}^{*}$ variable at parton-level in region SR-3L-tZq.

Migration matrix for cos $\theta_{Z}^{*}$ variable at parton-level in region SR-3L-WZ.

Migration matrix for cos $\theta_{Z}^{*}$ variable at parton-level in region SR-4L-DF.

Migration matrix for cos $\theta_{Z}^{*}$ variable at parton-level in region SR-4L-SF.

Migration matrix for cos $\theta_{Z}^{*}$ variable at parton-level in region CR-4L-ZZ.

Migration matrix for $|\Delta\Phi(t\bar{t}, Z)|/\pi$ variable at particle-level in region SR-3L-ttZ.

Migration matrix for $|\Delta\Phi(t\bar{t}, Z)|/\pi$ variable at particle-level in region SR-3L-tZq.

Migration matrix for $|\Delta\Phi(t\bar{t}, Z)|/\pi$ variable at particle-level in region SR-3L-WZ.

Migration matrix for $|\Delta\Phi(t\bar{t}, Z)|/\pi$ variable at particle-level in region SR-4L-DF.

Migration matrix for $|\Delta\Phi(t\bar{t}, Z)|/\pi$ variable at particle-level in region SR-4L-SF.

Migration matrix for $|\Delta\Phi(t\bar{t}, Z)|/\pi$ variable at particle-level in region CR-4L-ZZ.

Migration matrix for $|\Delta\Phi(t\bar{t}, Z)|/\pi$ variable at parton-level in region SR-3L-ttZ.

Migration matrix for $|\Delta\Phi(t\bar{t}, Z)|/\pi$ variable at parton-level in region SR-3L-tZq.

Migration matrix for $|\Delta\Phi(t\bar{t}, Z)|/\pi$ variable at parton-level in region SR-3L-WZ.

Migration matrix for $|\Delta\Phi(t\bar{t}, Z)|/\pi$ variable at parton-level in region SR-4L-DF.

Migration matrix for $|\Delta\Phi(t\bar{t}, Z)|/\pi$ variable at parton-level in region SR-4L-SF.

Migration matrix for $|\Delta\Phi(t\bar{t}, Z)|/\pi$ variable at parton-level in region CR-4L-ZZ.

Migration matrix for $|\Delta\Phi(l_{t}^{+}, l_{\bar{t}}^{-})|/\pi$ variable at particle-level in region SR-4L-DF.

Migration matrix for $|\Delta\Phi(l_{t}^{+}, l_{\bar{t}}^{-})|/\pi$ variable at particle-level in region SR-4L-SF.

Migration matrix for $|\Delta\Phi(l_{t}^{+}, l_{\bar{t}}^{-})|/\pi$ variable at particle-level in region CR-4L-ZZ.

Migration matrix for $|\Delta\Phi(l_{t}^{+}, l_{\bar{t}}^{-})|/\pi$ variable at parton-level in region SR-4L-DF.

Migration matrix for $|\Delta\Phi(l_{t}^{+}, l_{\bar{t}}^{-})|/\pi$ variable at parton-level in region SR-4L-SF.

Migration matrix for $|\Delta\Phi(l_{t}^{+}, l_{\bar{t}}^{-})|/\pi$ variable at parton-level in region CR-4L-ZZ.

Migration matrix for $|\Delta\Phi(Z, t_{lep})|/\pi$ variable at particle-level in region SR-3L-ttZ.

Migration matrix for $|\Delta\Phi(Z, t_{lep})|/\pi$ variable at particle-level in region SR-3L-tZq.

Migration matrix for $|\Delta\Phi(Z, t_{lep})|/\pi$ variable at particle-level in region SR-3L-WZ.

Migration matrix for $|\Delta\Phi(Z, t_{lep})|/\pi$ variable at parton-level in region SR-3L-ttZ.

Migration matrix for $|\Delta\Phi(Z, t_{lep})|/\pi$ variable at parton-level in region SR-3L-tZq.

Migration matrix for $|\Delta\Phi(Z, t_{lep})|/\pi$ variable at parton-level in region SR-3L-WZ.

Migration matrix for $|\Delta y(Z, t_{lep})|$ variable at particle-level in region SR-3L-ttZ.

Migration matrix for $|\Delta y(Z, t_{lep})|$ variable at particle-level in region SR-3L-tZq.

Migration matrix for $|\Delta y(Z, t_{lep})|$ variable at particle-level in region SR-3L-WZ.

Migration matrix for $|\Delta y(Z, t_{lep})|$ variable at parton-level in region SR-3L-ttZ.

Migration matrix for $|\Delta y(Z, t_{lep})|$ variable at parton-level in region SR-3L-tZq.

Migration matrix for $|\Delta y(Z, t_{lep})|$ variable at parton-level in region SR-3L-WZ.

Migration matrix for $H_{\text{T}}^{\text{ l}}$ variable at particle-level in region SR-4L-DF.

Migration matrix for $H_{\text{T}}^{\text{ l}}$ variable at particle-level in region SR-4L-SF.

Migration matrix for $H_{\text{T}}^{\text{ l}}$ variable at particle-level in region CR-4L-ZZ.

Migration matrix for $H_{\text{T}}^{\text{ l}}$ variable at parton-level in region SR-4L-DF.

Migration matrix for $H_{\text{T}}^{\text{ l}}$ variable at parton-level in region SR-4L-SF.

Migration matrix for $H_{\text{T}}^{\text{ l}}$ variable at parton-level in region CR-4L-ZZ.

Migration matrix for $H_{\text{T}}^{\text{ l}}$ variable at particle-level in region SR-3L-ttZ.

Migration matrix for $H_{\text{T}}^{\text{ l}}$ variable at particle-level in region SR-3L-tZq.

Migration matrix for $H_{\text{T}}^{\text{ l}}$ variable at particle-level in region SR-3L-WZ.

Migration matrix for $H_{\text{T}}^{\text{ l}}$ variable at parton-level in region SR-3L-ttZ.

Migration matrix for $H_{\text{T}}^{\text{ l}}$ variable at parton-level in region SR-3L-tZq.

Migration matrix for $H_{\text{T}}^{\text{ l}}$ variable at parton-level in region SR-3L-WZ.

Migration matrix for $m^{t\bar{t}Z}$ variable at particle-level in region SR-3L-ttZ.

Migration matrix for $m^{t\bar{t}Z}$ variable at particle-level in region SR-3L-tZq.

Migration matrix for $m^{t\bar{t}Z}$ variable at particle-level in region SR-3L-WZ.

Migration matrix for $m^{t\bar{t}Z}$ variable at particle-level in region SR-4L-DF.

Migration matrix for $m^{t\bar{t}Z}$ variable at particle-level in region SR-4L-SF.

Migration matrix for $m^{t\bar{t}Z}$ variable at particle-level in region CR-4L-ZZ.

Migration matrix for $m^{t\bar{t}Z}$ variable at parton-level in region SR-3L-ttZ.

Migration matrix for $m^{t\bar{t}Z}$ variable at parton-level in region SR-3L-tZq.

Migration matrix for $m^{t\bar{t}Z}$ variable at parton-level in region SR-3L-WZ.

Migration matrix for $m^{t\bar{t}Z}$ variable at parton-level in region SR-4L-DF.

Migration matrix for $m^{t\bar{t}Z}$ variable at parton-level in region SR-4L-SF.

Migration matrix for $m^{t\bar{t}Z}$ variable at parton-level in region CR-4L-ZZ.

Migration matrix for $m^{t\bar{t}}$ variable at particle-level in region SR-3L-ttZ.

Migration matrix for $m^{t\bar{t}}$ variable at particle-level in region SR-3L-tZq.

Migration matrix for $m^{t\bar{t}}$ variable at particle-level in region SR-3L-WZ.

Migration matrix for $m^{t\bar{t}}$ variable at particle-level in region SR-4L-DF.

Migration matrix for $m^{t\bar{t}}$ variable at particle-level in region SR-4L-SF.

Migration matrix for $m^{t\bar{t}}$ variable at particle-level in region CR-4L-ZZ.

Migration matrix for $m^{t\bar{t}}$ variable at parton-level in region SR-3L-ttZ.

Migration matrix for $m^{t\bar{t}}$ variable at parton-level in region SR-3L-tZq.

Migration matrix for $m^{t\bar{t}}$ variable at parton-level in region SR-3L-WZ.

Migration matrix for $m^{t\bar{t}}$ variable at parton-level in region SR-4L-DF.

Migration matrix for $m^{t\bar{t}}$ variable at parton-level in region SR-4L-SF.

Migration matrix for $m^{t\bar{t}}$ variable at parton-level in region CR-4L-ZZ.

Migration matrix for $N_{\text{jets}}$ variable at particle-level in region SR-4L-DF.

Migration matrix for $N_{\text{jets}}$ variable at particle-level in region SR-4L-SF.

Migration matrix for $N_{\text{jets}}$ variable at particle-level in region CR-4L-ZZ.

Migration matrix for $N_{\text{jets}}$ variable at particle-level in region SR-3L-ttZ.

Migration matrix for $N_{\text{jets}}$ variable at particle-level in region SR-3L-tZq.

Migration matrix for $N_{\text{jets}}$ variable at particle-level in region SR-3L-WZ.

Migration matrix for $p^{Z}_{T}$ variable at particle-level in region SR-3L-ttZ.

Migration matrix for $p^{Z}_{T}$ variable at particle-level in region SR-3L-tZq.

Migration matrix for $p^{Z}_{T}$ variable at particle-level in region SR-3L-WZ.

Migration matrix for $p^{Z}_{T}$ variable at particle-level in region SR-4L-DF.

Migration matrix for $p^{Z}_{T}$ variable at particle-level in region SR-4L-SF.

Migration matrix for $p^{Z}_{T}$ variable at particle-level in region CR-4L-ZZ.

Migration matrix for $p^{Z}_{T}$ variable at parton-level in region SR-3L-ttZ.

Migration matrix for $p^{Z}_{T}$ variable at parton-level in region SR-3L-tZq.

Migration matrix for $p^{Z}_{T}$ variable at parton-level in region SR-3L-WZ.

Migration matrix for $p^{Z}_{T}$ variable at parton-level in region SR-4L-DF.

Migration matrix for $p^{Z}_{T}$ variable at parton-level in region SR-4L-SF.

Migration matrix for $p^{Z}_{T}$ variable at parton-level in region CR-4L-ZZ.

Migration matrix for $p_{T}^{\mathrm{top}}$ variable at particle-level in region SR-3L-ttZ.

Migration matrix for $p_{T}^{\mathrm{top}}$ variable at particle-level in region SR-3L-tZq.

Migration matrix for $p_{T}^{\mathrm{top}}$ variable at particle-level in region SR-3L-WZ.

Migration matrix for $p_{T}^{\mathrm{top}}$ variable at particle-level in region SR-4L-DF.

Migration matrix for $p_{T}^{\mathrm{top}}$ variable at particle-level in region SR-4L-SF.

Migration matrix for $p_{T}^{\mathrm{top}}$ variable at particle-level in region CR-4L-ZZ.

Migration matrix for $p_{T}^{\mathrm{top}}$ variable at parton-level in region SR-3L-ttZ.

Migration matrix for $p_{T}^{\mathrm{top}}$ variable at parton-level in region SR-3L-tZq.

Migration matrix for $p_{T}^{\mathrm{top}}$ variable at parton-level in region SR-3L-WZ.

Migration matrix for $p_{T}^{\mathrm{top}}$ variable at parton-level in region SR-4L-DF.

Migration matrix for $p_{T}^{\mathrm{top}}$ variable at parton-level in region SR-4L-SF.

Migration matrix for $p_{T}^{\mathrm{top}}$ variable at parton-level in region CR-4L-ZZ.

Migration matrix for $p_{T}^{t\bar{t}}$ variable at particle-level in region SR-3L-ttZ.

Migration matrix for $p_{T}^{t\bar{t}}$ variable at particle-level in region SR-3L-tZq.

Migration matrix for $p_{T}^{t\bar{t}}$ variable at particle-level in region SR-3L-WZ.

Migration matrix for $p_{T}^{t\bar{t}}$ variable at particle-level in region SR-4L-DF.

Migration matrix for $p_{T}^{t\bar{t}}$ variable at particle-level in region SR-4L-SF.

Migration matrix for $p_{T}^{t\bar{t}}$ variable at particle-level in region CR-4L-ZZ.

Migration matrix for $p_{T}^{t\bar{t}}$ variable at parton-level in region SR-3L-ttZ.

Migration matrix for $p_{T}^{t\bar{t}}$ variable at parton-level in region SR-3L-tZq.

Migration matrix for $p_{T}^{t\bar{t}}$ variable at parton-level in region SR-3L-WZ.

Migration matrix for $p_{T}^{t\bar{t}}$ variable at parton-level in region SR-4L-DF.

Migration matrix for $p_{T}^{t\bar{t}}$ variable at parton-level in region SR-4L-SF.

Migration matrix for $p_{T}^{t\bar{t}}$ variable at parton-level in region CR-4L-ZZ.

Migration matrix for $p_{\text{T}}^{\ell, non-Z}$ variable at particle-level in region SR-3L-ttZ.

Migration matrix for $p_{\text{T}}^{\ell, non-Z}$ variable at particle-level in region SR-3L-tZq.

Migration matrix for $p_{\text{T}}^{\ell, non-Z}$ variable at particle-level in region SR-3L-WZ.

Migration matrix for $p_{\text{T}}^{\ell, non-Z}$ variable at parton-level in region SR-3L-ttZ.

Migration matrix for $p_{\text{T}}^{\ell, non-Z}$ variable at parton-level in region SR-3L-tZq.

Migration matrix for $p_{\text{T}}^{\ell, non-Z}$ variable at parton-level in region SR-3L-WZ.

Migration matrix for $|y^{Z}$| variable at particle-level in region SR-3L-ttZ.

Migration matrix for $|y^{Z}$| variable at particle-level in region SR-3L-tZq.

Migration matrix for $|y^{Z}$| variable at particle-level in region SR-3L-WZ.

Migration matrix for $|y^{Z}$| variable at particle-level in region SR-4L-DF.

Migration matrix for $|y^{Z}$| variable at particle-level in region SR-4L-SF.

Migration matrix for $|y^{Z}$| variable at particle-level in region CR-4L-ZZ.

Migration matrix for $|y^{Z}$| variable at parton-level in region SR-3L-ttZ.

Migration matrix for $|y^{Z}$| variable at parton-level in region SR-3L-tZq.

Migration matrix for $|y^{Z}$| variable at parton-level in region SR-3L-WZ.

Migration matrix for $|y^{Z}$| variable at parton-level in region SR-4L-DF.

Migration matrix for $|y^{Z}$| variable at parton-level in region SR-4L-SF.

Migration matrix for $|y^{Z}$| variable at parton-level in region CR-4L-ZZ.

Migration matrix for $|y^{t\bar{t}Z}|$ variable at particle-level in region SR-3L-ttZ.

Migration matrix for $|y^{t\bar{t}Z}|$ variable at particle-level in region SR-3L-tZq.

Migration matrix for $|y^{t\bar{t}Z}|$ variable at particle-level in region SR-3L-WZ.

Migration matrix for $|y^{t\bar{t}Z}|$ variable at particle-level in region SR-4L-DF.

Migration matrix for $|y^{t\bar{t}Z}|$ variable at particle-level in region SR-4L-SF.

Migration matrix for $|y^{t\bar{t}Z}|$ variable at particle-level in region CR-4L-ZZ.

Migration matrix for $|y^{t\bar{t}Z}|$ variable at parton-level in region SR-3L-ttZ.

Migration matrix for $|y^{t\bar{t}Z}|$ variable at parton-level in region SR-3L-tZq.

Migration matrix for $|y^{t\bar{t}Z}|$ variable at parton-level in region SR-3L-WZ.

Migration matrix for $|y^{t\bar{t}Z}|$ variable at parton-level in region SR-4L-DF.

Migration matrix for $|y^{t\bar{t}Z}|$ variable at parton-level in region SR-4L-SF.

Migration matrix for $|y^{t\bar{t}Z}|$ variable at parton-level in region CR-4L-ZZ.

Covariance matrix for absolute cross section as a function of $p_{T}^{\mathrm{top}}$ at particle-level.

Covariance matrix for normalized cross section as a function of $p_{T}^{\mathrm{top}}$ at particle-level.

Covariance matrix for absolute cross section as a function of $p_{T}^{\mathrm{top}}$ at parton-level.

Covariance matrix for normalized cross section as a function of $p_{T}^{\mathrm{top}}$ at parton-level.

Covariance matrix for absolute cross section as a function of $p_{T}^{t\bar{t}}$ at particle-level.

Covariance matrix for normalized cross section as a function of $p_{T}^{t\bar{t}}$ at particle-level.

Covariance matrix for absolute cross section as a function of $p_{T}^{t\bar{t}}$ at parton-level.

Covariance matrix for normalized cross section as a function of $p_{T}^{t\bar{t}}$ at parton-level.

Covariance matrix for absolute cross section as a function of $|\Delta\Phi(t\bar{t}, Z)|/\pi$ at particle-level.

Covariance matrix for normalized cross section as a function of $|\Delta\Phi(t\bar{t}, Z)|/\pi$ at particle-level.

Covariance matrix for absolute cross section as a function of $|\Delta\Phi(t\bar{t}, Z)|/\pi$ at parton-level.

Covariance matrix for normalized cross section as a function of $|\Delta\Phi(t\bar{t}, Z)|/\pi$ at parton-level.

Covariance matrix for absolute cross section as a function of $m^{t\bar{t}Z}$ at particle-level.

Covariance matrix for normalized cross section as a function of $m^{t\bar{t}Z}$ at particle-level.

Covariance matrix for absolute cross section as a function of $m^{t\bar{t}Z}$ at parton-level.

Covariance matrix for normalized cross section as a function of $m^{t\bar{t}Z}$ at parton-level.

Covariance matrix for absolute cross section as a function of $m^{t\bar{t}}$ at particle-level.

Covariance matrix for normalized cross section as a function of $m^{t\bar{t}}$ at particle-level.

Covariance matrix for absolute cross section as a function of $m^{t\bar{t}}$ at parton-level.

Covariance matrix for normalized cross section as a function of $m^{t\bar{t}}$ at parton-level.

Covariance matrix for absolute cross section as a function of $|y^{t\bar{t}Z}|$ at particle-level.

Covariance matrix for normalized cross section as a function of $|y^{t\bar{t}Z}|$ at particle-level.

Covariance matrix for absolute cross section as a function of $|y^{t\bar{t}Z}|$ at parton-level.

Covariance matrix for normalized cross section as a function of $|y^{t\bar{t}Z}|$ at parton-level.

Covariance matrix for absolute cross section as a function of cos $\theta_{Z}^{*}$ at particle-level.

Covariance matrix for normalized cross section as a function of cos $\theta_{Z}^{*}$ at particle-level.

Covariance matrix for absolute cross section as a function of cos $\theta_{Z}^{*}$ at parton-level.

Covariance matrix for normalized cross section as a function of cos $\theta_{Z}^{*}$ at parton-level.

Covariance matrix for absolute cross section as a function of $|\Delta\Phi(l_{t}^{+}, l_{\bar{t}}^{-})|/\pi$ at particle-level.

Covariance matrix for normalized cross section as a function of $|\Delta\Phi(l_{t}^{+}, l_{\bar{t}}^{-})|/\pi$ at particle-level.

Covariance matrix for absolute cross section as a function of $|\Delta\Phi(l_{t}^{+}, l_{\bar{t}}^{-})|/\pi$ at parton-level.

Covariance matrix for normalized cross section as a function of $|\Delta\Phi(l_{t}^{+}, l_{\bar{t}}^{-})|/\pi$ at parton-level.

Covariance matrix for absolute cross section as a function of $|\Delta\Phi(Z, t_{lep})|/\pi$ at particle-level.

Covariance matrix for normalized cross section as a function of $|\Delta\Phi(Z, t_{lep})|/\pi$ at particle-level.

Covariance matrix for absolute cross section as a function of $|\Delta\Phi(Z, t_{lep})|/\pi$ at parton-level.

Covariance matrix for normalized cross section as a function of $|\Delta\Phi(Z, t_{lep})|/\pi$ at parton-level.

Covariance matrix for absolute cross section as a function of $|\Delta y(Z, t_{lep})|$ at particle-level.

Covariance matrix for normalized cross section as a function of $|\Delta y(Z, t_{lep})|$ at particle-level.

Covariance matrix for absolute cross section as a function of $|\Delta y(Z, t_{lep})|$ at parton-level.

Covariance matrix for normalized cross section as a function of $|\Delta y(Z, t_{lep})|$ at parton-level.

Covariance matrix for absolute cross section as a function of $H_{\text{T}}^{\text{l}}$ at in the tetralepton channel particle-level.

Covariance matrix for normalized cross section as a function of $H_{\text{T}}^{\text{l}}$ at in the tetralepton channel particle-level.

Covariance matrix for absolute cross section as a function of $H_{\text{T}}^{\text{l}}$ at in the tetralepton channel parton-level.

Covariance matrix for normalized cross section as a function of $H_{\text{T}}^{\text{l}}$ in the tetralepton channel at parton-level.

Covariance matrix for absolute cross section as a function of $H_{\text{T}}^{\text{l}}$ in the trilepton channel at particle-level.

Covariance matrix for normalized cross section as a function of $H_{\text{T}}^{\text{l}}$ in the trilepton channel at particle-level.

Covariance matrix for absolute cross section as a function of $H_{\text{T}}^{\text{l}}$ in the trilepton channel at parton-level.

Covariance matrix for normalized cross section as a function of $H_{\text{T}}^{\text{l}}$ in the trilepton channel at parton-level.

Covariance matrix for absolute cross section as a function of $N_{\text{jets}}$ in the tetralepton channel at particle-level.

Covariance matrix for normalized cross section as a function of $N_{\text{jets}}$ in the tetralepton channel at particle-level.

Covariance matrix for absolute cross section as a function of $N_{\text{jets}}$ in the trilepton channel at particle-level.

Covariance matrix for normalized cross section as a function of $N_{\text{jets}}$ in the trilepton channel at particle-level.

Covariance matrix for absolute cross section as a function of $p^{Z}_{T}$ at particle-level.

Covariance matrix for normalized cross section as a function of $p^{Z}_{T}$ at particle-level.

Covariance matrix for absolute cross section as a function of $p^{Z}_{T}$ at parton-level.

Covariance matrix for normalized cross section as a function of $p^{Z}_{T}$ at parton-level.

Covariance matrix for absolute cross section as a function of $p_{\text{T}}^{\ell, non-Z}$ at particle-level.

Covariance matrix for normalized cross section as a function of $p_{\text{T}}^{\ell, non-Z}$ at particle-level.

Covariance matrix for absolute cross section as a function of $p_{\text{T}}^{\ell, non-Z}$ at parton-level.

Covariance matrix for normalized cross section as a function of $p_{\text{T}}^{\ell, non-Z}$ at parton-level.

Covariance matrix for absolute cross section as a function of $|y^{Z}$| at particle-level.

Covariance matrix for normalized cross section as a function of $|y^{Z}$| at particle-level.

Covariance matrix for absolute cross section as a function of $|y^{Z}$| at parton-level.

Covariance matrix for normalized cross section as a function of $|y^{Z}$| at parton-level.

Ranking of nuisance parameters and background normalizations on signal strength for inclusive cross section measurement in combination of all channels

Observed and expected 68% and 95% credible intervals for the top-boson operators, in the marginalised linear fit.

Observed and expected 68% and 95% credible intervals for the top-boson operators, in the marginalised quadratic fit.

Observed and expected 68% and 95% credible intervals for the top-boson operators, in the independent quadratic fits (allowing only one Wilson Coefficient to be non-zero).

Observed and expected 68% and 95% credible intervals for the four-quark operators, in the marginalised linear fit.

Observed and expected 68% and 95% credible intervals for the four-quark operators, in the marginalised quadratic fit.

Observed and expected 68% and 95% credible intervals for the four-quark operators, in the independent quadratic fits (allowing only one Wilson Coefficient to be non-zero).

Observed and expected 68% and 95% credible intervals for Fisher-rotated directions of EFT sensitivity, in the marginalised linear fit.

Correlation matrix of the input particle-level observables used in the EFT fit.

Expected and observed 95% CL upper limits on the branching ratio of the Higgs boson decay to $Za$ times the branching ratio $a $->$ \gamma \gamma$ as a function of the $a$ particle mass in the merged ($m_a \le 2$ GeV) and the resolved ($m_a > 2$ GeV) categories. The figure uses $pp$ collision data at $\sqrt{s}=13$ TeV corresponding to 139 fb$^{-1}$.

ATLAS observed 95% CL exclusion contours limits in terms of the ALP mass and its effective coupling to photons, $|C_{\gamma\gamma}|/\Lambda$, for different values of the Higgs coupling to $Za$, $|C_{Zh}|/\Lambda$. The figure uses $pp$ collision data at $\sqrt{s}=13$ TeV corresponding to 139 fb$^{-1}$. The value 100 in the tables is presentational, denoting the upper edge of the plot.

The distribution of the fit discriminant mTVH in the 3+ b-tag signal region of the vvbb channel for the mH=300 GeV hypothesis. <br><br><a href="?table=overview">return to overview</a>

The m(tt) distribution in the Lhi3_Zin region of the lltt channel. <br><br><a href="?table=overview">return to overview</a>

The m(bb) distribution in the 2 b-tag 0L region of the vvbb channel. <br><br><a href="?table=overview">return to overview</a>

The m(bb) distribution in the 3+ b-tag 0L region of the vvbb channel. <br><br><a href="?table=overview">return to overview</a>

Expected and observed upper limits at 95% CL on sigma(gg->A)*B(A->ZH)*B(H->tt) in the (mA,mH) plane. The limits are shown for tanbeta=1. The tanbeta value is relevant only for the choice of the A boson width. <br><a href="?table=overview">return to overview</a>

Expected and observed upper limits at 95% CL on sigma(bbA)*B(A->ZH)*B(H->tt)$ in the (mA,mH) plane. The limits are shown for tanbeta=10. The tanbeta value is relevant only for the choice of the A boson width. <br><br><a href="?table=overview">return to overview</a>

Expected and observed upper limits at 95% CL on sigma(gg->A)*B(A->ZH)*B(H->bb) in the (mA,mH) plane. The limits are shown for tanbeta=1. The tanbeta value is relevant only for the choice of the A boson width. <br><br><a href="?table=overview">return to overview</a>

Expected and observed upper limits at 95% CL on sigma(bbA)*B(A->ZH)*B(H->bb) in the (mA,mH) plane. The limits are shown for tanbeta=10. The tanbeta value is relevant only for the choice of the A boson width. <br><br><a href="?table=overview">return to overview</a>

Expected and observed upper limits at 95% CL on sigma(gg->A)*B(A->ZH)*B(H->tt) in he (mA,mH) plane. The limits are shown for tanbeta=0.5. The tanbeta value is relevant only for the choice of the A boson width. <br><a href="?table=overview">return to overview</a>

Expected and observed upper limits at 95% CL on sigma(gg->A)*B(A->ZH)*B(H->tt) in the (mA,mH) plane. The limits are shown for tanbeta=5. The tanbeta value is relevant only for the choice of the A boson width. <br><a href="?table=overview">return to overview</a>

Expected and observed upper limits at 95% CL on sigma(bbA)*B(A->ZH)*B(H->tt) in the (mA,mH) plane. The limits are shown for tanbeta=1. The tanbeta value is relevant only for the choice of the A boson width. <br><br><a href="?table=overview">return to overview</a>

Expected and observed upper limits at 95% CL on sigma(bbA)*B(A->ZH)*B(H->tt)$ in the (mA,mH) plane. The limits are shown for tanbeta=5. The tanbeta value is relevant only for the choice of the A boson width. <br><br><a href="?table=overview">return to overview</a>

Expected and observed upper limits at 95% CL on sigma(gg->A)*B(A->ZH)*B(H->bb) in the (mA,mH) plane. The limits are shown for tanbeta=0.5. The tanbeta value is relevant only for the choice of the A boson width. <br><br><a href="?table=overview">return to overview</a>

Expected and observed upper limits at 95% CL on sigma(gg->A)*B(A->ZH)*B(H->bb) in the (mA,mH) plane. The limits are shown for tanbeta=5. The tanbeta value is relevant only for the choice of the A boson width. <br><br><a href="?table=overview">return to overview</a>

Expected and observed upper limits at 95% CL on sigma(bbA)*B(A->ZH)*B(H->bb) in the (mA,mH) plane. The limits are shown for tanbeta=1. The tanbeta value is relevant only for the choice of the A boson width. <br><br><a href="?table=overview">return to overview</a>

Expected and observed upper limits at 95% CL on sigma(bbA)*B(A->ZH)*B(H->bb) in the (mA,mH) plane. The limits are shown for tanbeta=5. The tanbeta value is relevant only for the choice of the A boson width. <br><br><a href="?table=overview">return to overview</a>

Expected and observed upper limits at 95% CL on sigma(bbA)*B(A->ZH)*B(H->bb) in the (mA,mH) plane. The limits are shown for tanbeta=20. The tanbeta value is relevant only for the choice of the A boson width. <br><br><a href="?table=overview">return to overview</a>

The distribution of the fit discriminant m(lltt)-m(tt) in the signal region of the lltt channel for the mH=450 GeV hypothesis with the bbA signal shown for comparison. <br><br><a href="?table=overview">return to overview</a>

The m(tt) distribution in the L3hi_Zin region of the lltt channel for the mH=450 GeV hypothesis with the bbA signal shown for comparison. <br><br><a href="?table=overview">return to overview</a>

The distribution of the fit discriminant m(lltt)-m(tt) in the signal region of the lltt channel for the mH=350 GeV hypothesis. <br><br><a href="?table=overview">return to overview</a>

The distribution of the fit discriminant m(lltt)-m(tt) in the signal region of the lltt channel for the mH=400 GeV hypothesis. <br><br><a href="?table=overview">return to overview</a>

The distribution of the fit discriminant m(lltt)-m(tt) in the signal region of the lltt channel for the mH=500 GeV hypothesis. <br><br><a href="?table=overview">return to overview</a>

The distribution of the fit discriminant m(lltt)-m(tt) in the signal region of the lltt channel for the mH=550 GeV hypothesis. <br><br><a href="?table=overview">return to overview</a>

The distribution of the fit discriminant m(lltt)-m(tt) in the signal region of the lltt channel for the mH=600 GeV hypothesis. <br><br><a href="?table=overview">return to overview</a>

The distribution of the fit discriminant m(lltt)-m(tt) in the signal region of the lltt channel for the mH=700 GeV hypothesis. <br><br><a href="?table=overview">return to overview</a>

The distribution of the fit discriminant m(lltt)-m(tt) in the signal region of the lltt channel for the mH=800 GeV hypothesis. <br><br><a href="?table=overview">return to overview</a>

The distribution of the fit discriminant mTVH in the 2 b-tag signal region of the vvbb channel for the mH=130 GeV hypothesis. <br><br><a href="?table=overview">return to overview</a>

The distribution of the fit discriminant mTVH in the 2 b-tag signal region of the vvbb channel for the mH=150 GeV hypothesis. <br><br><a href="?table=overview">return to overview</a>

The distribution of the fit discriminant mTVH in the 2 b-tag signal region of the vvbb channel for the mH=200 GeV hypothesis. <br><br><a href="?table=overview">return to overview</a>

The distribution of the fit discriminant mTVH in the 2 b-tag signal region of the vvbb channel for the mH=250 GeV hypothesis. <br><br><a href="?table=overview">return to overview</a>

The distribution of the fit discriminant mTVH in the 2 b-tag signal region of the vvbb channel for the mH=350 GeV hypothesis. <br><br><a href="?table=overview">return to overview</a>

The distribution of the fit discriminant mTVH in the 2 b-tag signal region of the vvbb channel for the mH=400 GeV hypothesis. <br><br><a href="?table=overview">return to overview</a>

The distribution of the fit discriminant mTVH in the 2 b-tag signal region of the vvbb channel for the mH=450 GeV hypothesis. <br><br><a href="?table=overview">return to overview</a>

The distribution of the fit discriminant mTVH in the 2 b-tag signal region of the vvbb channel for the mH=500 GeV hypothesis. <br><br><a href="?table=overview">return to overview</a>

The distribution of the fit discriminant mTVH in the 2 b-tag signal region of the vvbb channel for the mH=600 GeV hypothesis. <br><br><a href="?table=overview">return to overview</a>

The distribution of the fit discriminant mTVH in the 2 b-tag signal region of the vvbb channel for the mH=700 GeV hypothesis. <br><br><a href="?table=overview">return to overview</a>

The distribution of the fit discriminant mTVH in the 2 b-tag signal region of the vvbb channel for the mH=800 GeV hypothesis. <br><br><a href="?table=overview">return to overview</a>

The distribution of the fit discriminant mTVH i+ the 3p b-tag signal region of the vvbb channel for the mH=130 GeV hypothesis. <br><br><a href="?table=overview">return to overview</a>

The distribution of the fit discriminant mTVH in the 3+ b-tag signal region of the vvbb channel for the mH=150 GeV hypothesis. <br><br><a href="?table=overview">return to overview</a>

The distribution of the fit discriminant mTVH in the 3+ b-tag signal region of the vvbb channel for the mH=200 GeV hypothesis. <br><br><a href="?table=overview">return to overview</a>

The distribution of the fit discriminant mTVH in the 3+ b-tag signal region of the vvbb channel for the mH=250 GeV hypothesis. <br><br><a href="?table=overview">return to overview</a>

The distribution of the fit discriminant mTVH in the 3+ b-tag signal region of the vvbb channel for the mH=350 GeV hypothesis. <br><br><a href="?table=overview">return to overview</a>

The distribution of the fit discriminant mTVH in the 3+ b-tag signal region of the vvbb channel for the mH=400 GeV hypothesis. <br><br><a href="?table=overview">return to overview</a>

The distribution of the fit discriminant mTVH in the 3+ b-tag signal region of the vvbb channel for the mH=450 GeV hypothesis. <br><br><a href="?table=overview">return to overview</a>

The distribution of the fit discriminant mTVH in the 3+ b-tag signal region of the vvbb channel for the mH=500 GeV hypothesis. <br><br><a href="?table=overview">return to overview</a>

The distribution of the fit discriminant mTVH in the 3+ b-tag signal region of the vvbb channel for the mH=600 GeV hypothesis. <br><br><a href="?table=overview">return to overview</a>

The distribution of the fit discriminant mTVH in the 3+ b-tag signal region of the vvbb channel for the mH=700 GeV hypothesis. <br><br><a href="?table=overview">return to overview</a>

The distribution of the fit discriminant mTVH in the 3+ b-tag signal region of the vvbb channel for the mH=800 GeV hypothesis. <br><br><a href="?table=overview">return to overview</a>

The distribution of the fit discriminant mTVH in the 2 b-tag 2L region of the vvbb channel. <br><br><a href="?table=overview">return to overview</a>

The distribution of the fit discriminant mTVH in the 2 b-tag em region of the vvbb channel. <br><br><a href="?table=overview">return to overview</a>

The distribution of the fit discriminant mTVH in the 3+ b-tag 2L region of the vvbb channel. <br><br><a href="?table=overview">return to overview</a>

The distribution of the fit discriminant mTVH in the 3+ b-tag em region of the vvbb channel. <br><br><a href="?table=overview">return to overview</a>

Local significance for lltt ggF signals for tanbeta=1. <br><br><a href="?table=overview">return to overview</a>

Local significance for lltt bbA signals for tanbeta=20. <br><br><a href="?table=overview">return to overview</a>

Local significance for vvbb ggF signals for tanbeta=1. <br><br><a href="?table=overview">return to overview</a>

Local significance for vvbb bbA signals for tanbeta=20. <br><br><a href="?table=overview">return to overview</a>

Signal acceptance times efficiency for lltt ggF signals. <br><br><a href="?table=overview">return to overview</a>

Signal acceptance times efficiency for lltt bbA signals. <br><br><a href="?table=overview">return to overview</a>

Signal acceptance times efficiency for vvbb ggF signals. <br><br><a href="?table=overview">return to overview</a>

Signal acceptance times efficiency for vvbb bbA signals. <br><br><a href="?table=overview">return to overview</a>

The pT distribution of the third lepton in the L3hi_Zin region of the lltt channel. <br><br><a href="?table=overview">return to overview</a>

The eta distribution of the H boson candidate in the rest frame of the A boson candidate in the L3hi_Zin region of the lltt channel. <br><br><a href="?table=overview">return to overview</a>

The ETmiss distribution in the 2 b-tag 0L region of the vvbb channel. <br><br><a href="?table=overview">return to overview</a>

The m_top_near distribution in the 2 b-tag 0L region of the vvbb channel. <br><br><a href="?table=overview">return to overview</a>

The ETmiss distribution in the 3+ b-tag 0L region of the vvbb channel. <br><br><a href="?table=overview">return to overview</a>

The m_top_near distribution in the 3+ b-tag 0L region of the vvbb channel. <br><br><a href="?table=overview">return to overview</a>

Predicted and observed yields as a function of $m_{4\ell}$ in the VBS-Suppressed region. Overflow events are included in the last bin of the distribution.

Relative systematic uncertainties in the unfolded differential cross section for 4ljj production in the VBS-Enhanced region as a function of m_{jj}

Relative systematic uncertainties in the unfolded differential cross section for 4ljj production in the VBS-Enhanced region as a function of m_{4\ell}

Differential cross-sections for inclusive 4ljj production in the VBS-Enhanced region as a function of $m_{jj}$. Overflow events are included in the last bin of the distribution.

Differential cross-sections for inclusive 4ljj production in the VBS-Suppressed region as a function of $m_{jj}$. Overflow events are included in the last bin of the distribution.

Differential cross-sections for inclusive 4ljj production in the VBS-Enhanced region as a function of $m_{4\ell}$. Overflow events are included in the last bin of the distribution.

Differential cross-sections for inclusive 4ljj production in the VBS-Suppressed region as a function of $m_{4\ell}$. Overflow events are included in the last bin of the distribution.

Differential cross-sections for inclusive 4ljj production in the VBS-Enhanced region as a function of $ p_{T,4l}$. Overflow events are included in the last bin of the distribution.

Differential cross-sections for inclusive 4ljj production in the VBS-Suppressed region as a function of $ p_{T,4l}$. Overflow events are included in the last bin of the distribution.

Differential cross-sections for inclusive 4ljj production in the VBS-Enhanced region as a function of $\Delta\phi_{jj}$. Overflow events are included in the last bin of the distribution.

Differential cross-sections for inclusive 4ljj production in the VBS-Suppressed region as a function of $\Delta\phi_{jj}$. Overflow events are included in the last bin of the distribution.

Differential cross-sections for inclusive 4ljj production in the VBS-Enhanced region as a function of $|\Delta y_{jj}|$. Overflow events are included in the last bin of the distribution.

Differential cross-sections for inclusive 4ljj production in the VBS-Suppressed region as a function of $|\Delta y_{jj}|$. Overflow events are included in the last bin of the distribution.

Differential cross-sections for inclusive 4ljj production in the VBS-Enhanced region as a function of $cos(\theta^{*}_{12})$. Overflow events are included in the last bin of the distribution.

Differential cross-sections for inclusive 4ljj production in the VBS-Suppressed region as a function of $cos(\theta^{*}_{12})$. Overflow events are included in the last bin of the distribution.

Differential cross-sections for inclusive 4ljj production in the VBS-Enhanced region as a function of $cos(\theta^{*}_{34})$. Overflow events are included in the last bin of the distribution.

Differential cross-sections for inclusive 4ljj production in the VBS-Suppressed region as a function of $cos(\theta^{*}_{34})$. Overflow events are included in the last bin of the distribution.

Differential cross-sections for inclusive 4ljj production in the VBS-Enhanced region as a function of $p_{T,jj}$. Overflow events are included in the last bin of the distribution.

Differential cross-sections for inclusive 4ljj production in the VBS-Suppressed region as a function of $p_{T,jj}$. Overflow events are included in the last bin of the distribution.

Differential cross-sections for inclusive 4ljj production in the VBS-Enhanced region as a function of $p_{T,4ljj}$. Overflow events are included in the last bin of the distribution.

Differential cross-sections for inclusive 4ljj production in the VBS-Suppressed region as a function of $p_{T,4ljj}$. Overflow events are included in the last bin of the distribution.

Differential cross-sections for inclusive 4ljj production in the VBS-Enhanced region as a function of $S_{T,4ljj}$. Overflow events are included in the last bin of the distribution.

Differential cross-sections for inclusive 4ljj production in the VBS-Suppressed region as a function of $S_{T,4ljj}$. Overflow events are included in the last bin of the distribution.

Statistical correlation between different bins of the unfolded cross-section for all the measured observables in the VBS-Enhanced region.

Statistical correlation between different bins of the unfolded cross-section for all the measured observables in the VBS-Suppressed region.

Table with post-fit yields in the four regions used in the measurement. The correlations of the nuisance parameters, as obtained by the fit, are considered for the calculation of the uncertainties.

Ratio of the $t\bar{t}$ to the $Z$-boson cross section compared to the predictions for several sets of parton distribution functions. The total cross section for $t\bar{t}$ production and the fiducial cross section for $Z$-boson production are used. The uncertainty on the prediction represents the PDF+$\alpha_{s}$ variations.

Nuisance parameter ranking plot for the $t\bar{t}$ signal-strength. Only the 10 most important nuisance parameters are shown. The impact of each nuisance parameter on the measured cross section, is defined as the shift induced in measured cross-section as the nuisance parameter is shifted by its pre-fit and post-fit uncertainties from its nominal best-fit value and then fixed while all other nuisance parameters profiled.

Nuisance parameter ranking plot for the ratio of the $t\bar{t}$ and $Z$-boson signal-strengths. Only the 10 most important nuisance parameters are shown. The impact of each nuisance parameter on the measured ratio, is defined as the shift induced in measured ratio as the nuisance parameter is shifted by its pre-fit and post-fit uncertainties from its nominal best-fit value and then fixed while all other nuisance parameters profiled.

Nuisance parameter ranking plot for the $Z$-boson signal-strength. Only the 10 most important nuisance parameters are shown. The impact of each nuisance parameter on the measured cross section, is defined as the shift induced in measured cross-section as the nuisance parameter is shifted by its pre-fit and post-fit uncertainties from its nominal best-fit value and then fixed while all other nuisance parameters profiled.