Pre-fit background composition of the SR$3\ell$ Prod. The table shows the event yields as opposed to just the percentages of the relevant background processes.

Post-fit plot of $H_\text{T}(\text{jets})$ in the SR$2\ell$ Dec from a signal-blinded background-only fit.

Post-fit plot of $m(t_\text{SM}, b\text{-jet}_0)$ in the SR$2\ell$ Prod from a signal-blinded background-only fit.

Post-fit plot of $m(\ell_\text{OS},\ell_\text{SS,1})$ in the SR$3\ell$ Dec from a signal-blinded background-only fit.

Post-fit plot of $m(\ell_\text{OS},\ell_\text{SS,1})$ in the SR$3\ell$ Prod from a signal-blinded background-only fit.

Post-fit plot of $D_\text{NN}(tHc)$ in the SR$2\ell$ Dec from the full fit to data.

Post-fit plot of $D_\text{NN}(tHc)$ in the SR$2\ell$ Prod from the full fit to data.

Post-fit plot of $D_\text{NN}(tHc)$ in the SR$3\ell$ Dec from the full fit to data.

Post-fit plot of $D_\text{NN}(tHc)$ in the SR$3\ell$ Prod from the full fit to data.

Post-fit plot of $p_\text{T}(\ell_1)$ in the CR$2\ell$ HF$e$ from the full fit to data.

Post-fit plot of $p_\text{T}(\ell_1)$ in the CR$2\ell$ HF$\mu$ from the full fit to data.

Post-fit plot of $p_\text{T}(\ell_1)$ in the CR$2\ell$ $t\bar{t}V$ from the full fit to data.

Post-fit plot of $p_\text{T}(\ell_2)$ in the CR$3\ell$ HF$e$ from the full fit to data.

Post-fit plot of $p_\text{T}(\ell_2)$ in the CR$3\ell$ HF$\mu$ from the full fit to data.

Post-fit plot of $p_\text{T}(b\text{-jet}_0)$ in the CR$3\ell$ $t\bar{t}W$ from the full fit to data.

Post-fit plot of $p_\text{T}(b\text{-jet}_0)$ in the CR$3\ell$ $t\bar{t}Z$ from the full fit to data.

Observed and expected upper exclusion limits on the branching ratio $\mathcal{B}(t\to Hu)$ for different analyses and their statistical combination.

Observed and expected upper exclusion limits on the branching ratio $\mathcal{B}(t\to Hc)$ for different analyses and their statistical combination.

Post-fit normalisation factors of free-floating background processes and the signal normalisation.

Post-fit predicted and observed yields in all $2\ell$SS signal and control regions. Pre-fit signal contributions for a signal cross section equivalent to $\mathcal{B}(t\to Hq)=0.1\,\%$ are given as well.

Post-fit predicted and observed yields in all $3\ell$ signal and control regions. Pre-fit signal contributions for a signal cross section equivalent to $\mathcal{B}(t\to Hq)=0.1\,\%$ are given as well.

Expected upper limits on $\mathcal{B}(t\to Hq)$ for the nominal fit and alternative fit configurations. One contains the full phase space but only considers statistical uncertainties. Two other configurations consider the full set of systematic uncertainties, but only encompass one final state.

Expected and observed upper limits on $\mathcal{B}(t\to Hq)$ and $|C_{u\phi}^{qt,tq}|$ for the full fit containing all systematic uncertainties.

Pre-fit plot of $H_\text{T}(\text{jets})$ in the SR$2\ell$ Dec from a signal-blinded background-only fit.

Pre-fit plot of $m(t_\text{SM}, b\text{-jet}_0)$ in the SR$2\ell$ Prod from a signal-blinded background-only fit.

Pre-fit plot of $m(\ell_\text{OS},\ell_\text{SS,1})$ in the SR$3\ell$ Dec from a signal-blinded background-only fit.

Pre-fit plot of $m(\ell_\text{OS},\ell_\text{SS,1})$ in the SR$3\ell$ Prod from a signal-blinded background-only fit.

Post-fit plot of $D_\text{NN}(tHu)$ in the SR$2\ell$ Dec from the full fit to data.

Post-fit plot of $D_\text{NN}(tHu)$ in the SR$2\ell$ Prod from the full fit to data.

Post-fit plot of $D_\text{NN}(tHu)$ in the SR$3\ell$ Dec from the full fit to data.

Post-fit plot of $D_\text{NN}(tHu)$ in the SR$3\ell$ Prod from the full fit to data.

Pre-fit plot of $D_\text{NN}(tHc)$ in the SR$2\ell$ Dec from the full fit to data.

Pre-fit plot of $D_\text{NN}(tHc)$ in the SR$2\ell$ Prod from the full fit to data.

Pre-fit plot of $D_\text{NN}(tHc)$ in the SR$3\ell$ Dec from the full fit to data.

Pre-fit plot of $D_\text{NN}(tHc)$ in the SR$3\ell$ Prod from the full fit to data.

Pre-fit plot of $D_\text{NN}(tHu)$ in the SR$2\ell$ Dec from the full fit to data.

Pre-fit plot of $D_\text{NN}(tHu)$ in the SR$2\ell$ Prod from the full fit to data.

Pre-fit plot of $D_\text{NN}(tHu)$ in the SR$3\ell$ Dec from the full fit to data.

Pre-fit plot of $D_\text{NN}(tHu)$ in the SR$3\ell$ Prod from the full fit to data.

Pre-fit plot of $p_\text{T}(\ell_1)$ in the CR$2\ell$ HF$e$ from the full fit to data.

Pre-fit plot of $p_\text{T}(\ell_1)$ in the CR$2\ell$ HF$\mu$ from the full fit to data.

Pre-fit plot of $p_\text{T}(\ell_1)$ in the CR$2\ell$ $t\bar{t}V$ from the full fit to data.

Pre-fit plot of $p_\text{T}(\ell_2)$ in the CR$3\ell$ HF$e$ from the full fit to data.

Pre-fit plot of $p_\text{T}(\ell_2)$ in the CR$3\ell$ HF$\mu$ from the full fit to data.

Pre-fit plot of $p_\text{T}(b\text{-jet}_0)$ in the CR$3\ell$ $t\bar{t}W$ from the full fit to data.

Pre-fit plot of $p_\text{T}(b\text{-jet}_0)$ in the CR$3\ell$ $t\bar{t}Z$ from the full fit to data.

Ranking of fit nuisance parameters according to their impact on the post-fit $tHu$ signal normalisation when fixed to $\pm1\sigma$

Ranking of fit nuisance parameters according to their impact on the post-fit $tHc$ signal normalisation when fixed to $\pm1\sigma$

Expected upper exclusion limits on the branching ratio $\mathcal{B}(t\to Hu)$ for each individual final state and the full analysis.

Expected upper exclusion limits on the branching ratio $\mathcal{B}(t\to Hc)$ for each individual final state and the full analysis.

The impact of the eight most important systematic uncertainties on \(R_t\), in %.

The post-fit yields in the two SRs. All uncertainties applied in the analysis are included.

Limits on EFT parameters and Vtb extracted from the analysis. The upper and lower limits correspond to 95% CL.

Results of the main analysis.

Pre-fit distribution of the NN discriminant \(D_{\text{nn}}\) in SR plus.

Post-fit distribution of the NN discriminant \(D_{\text{nn}}\) in SR plus.

Post-fit distribution of the \(\Delta\phi(\vec{p}_{T}^{miss},\mu)\) variable in the muon-plus CR.

Pre-fit distribution of the NN discriminant \(D_{\text{nn}}\) in SR minus.

Post-fit distribution of the NN discriminant \(D_{\text{nn}}\) in SR minus.

Post-fit distribution of the \(\Delta\phi(\vec{p}_{T}^{miss},\mu)\) variable in the muon-minus CR.

Pre-fit distribution of the invariant mass of the untagged jet and the b-tagged jet, \(m(jb)\), in SR plus.

Pre-fit distribution of the absolute value of the pseudorapidity of the untagged jet, \(|\eta(j)|\), in SR plus.

Pre-fit distribution of the absolute value of the difference in \(p_{\text{T}}\) between the reconstructed W boson and the jet pair, \(|\Delta p_{\text{T}}(W, jb)|\), in SR plus.

Pre-fit distribution of the difference in azimuth angle between the reconstructed W boson and the jet pair, \(|\Delta\phi(W, jb)|\), in SR plus.

Pre-fit distribution of the invariant mass of the reconstructed top quark, \(m(t)\), in SR plus.

Pre-fit distribution of the absolute value of the difference in pseudorapidity between the charged lepton and the untagged jet, \(|\Delta\eta(\ell ,j)|\), SR plus.

Pre-fit distribution of the angular distance of the charged lepton and the untagged jet, \(\Delta R(\ell ,j)\), in SR plus.

Pre-fit distribution of the absolute value of the difference in pseudorapidity between the b-tagged jet and the charged lepton, \(|\Delta\eta(b,\ell )|\), in SR plus.

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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.

Bayesian upper exclusion limits on the ratio $g_{W'}/g_{W}$ for an SSM-like W$\prime$ boson are shown. The unity coupling ratio (blue dotted curve) corresponds to the SSM common benchmark. The 68 and 95% quantiles of the limits are represented by the green and yellow bands, respectively.

Bayesian lower exclusion limits on the NUGIM G(221) mixing angle $\cot(\theta_{E})$ are shown as a function of the W$\prime$ boson mass. The theoretically excluded region is shaded in grey. The 68 and 95% quantiles of the limits are represented by the green and yellow bands, respectively.

Bayesian upper exclusion limits at 95% CL on the product of the production cross section and branching fraction of a QBH in an associated $ au$ lepton and neutrino final state. Minimum threshold masses $m_{th}$ of up to 6.6 TeV are excluded at 95% CL. The observed limit (solid line) is obtained from the intersection with the LO QBH cross section (blue dotted curve). The 68 and 95% quantiles of the limits are represented by the green and yellow bands, respectively.

Bayesian upper limits at 95% CL on the cross section of the process $pp\rightarrow\tau\nu$ mediated via LQ exchange in the t-channel. The inner (green) band and the outer (yellow) band indicate the regions containing 68 and 95%, respectively, of the distribution of limits expected under the background-only hypothesis. The predicted LQ cross section at LO in the three coupling benchmark scenarios is depicted in different colors for $g_{U}=1$. The uncertainty bandy correspond to the sum in quadrature of PDF and scale variations. The first benchmark scenario considers only couplings to left-handed SM fermions (i.e. $\beta_{\text{R}}^{ij}=0$) and is referred to as "best fit LH". The second benchmark, referred to as "best fit LH+RH", considers $|\beta_{\text{R}}^{\text{b}\tau}|=1$ and all other $\beta_{\text{R}}^{ij}=0$. In the third "democratic" benchmark, equal couplings only to LH fermions are assumed, i.e. $\beta_{\text{L}}^{ij}=1$ and $\beta_{\text{R}}^{ij}=0$ for all $i$ and $j$.

Expected and observed lower limits of the LQ mass as a function of the coupling $g_{U}$ in the LH scenario. The blue band shows the 68% and 95% regions of $g_{U}$ preferred by the fit to the b anomalies data.

Expected and observed lower limits of the LQ mass as a function of the coupling $g_{U}$ in the LH+RH scenario. The blue band shows the 68% and 95% regions of $g_{U}$ preferred by the fit to the b anomalies data.

Expected and observed lower limits of the LQ mass as a function of the coupling $g_{U}$ in the democratic scenario. The blue band shows the 68% and 95% regions of $g_{U}$ preferred by the fit to the b anomalies data.

Bayesian upper exclusion limits at 95% CL on each of the Wilson coefficients described by the EFT model based on 2016-2018 data. The three different coupling types represent a left-handed vector coupling ($\epsilon^{cb}_{L}$), tensor-like coupling ($\epsilon^{cb}_{T}$), and scalar-tensor-like coupling ($\epsilon^{cb}_{S_{L}}$). The 68 and 95% quantiles of the limits are represented by the green and yellow bands, respectively.

Summary of exclusion limits (expected and observed) calculated at 95% CL for full Run-2 CMS data.

Background prediction and observed data yields in the signal region bins. The background yields are obtained from the background-only fit and serve as input to the simplified likelihood reinterpretation scheme. The naming of the bins is "year_binnumber", following the binning from Figure 4.

Matrix of covariance coefficients between signal region bins. The coefficients are obtained from the background-only fit and serve as input to the simplified likelihood reinterpretation scheme. The naming of the bins is "year_binnumber", following the binning used in Figure 4.

Predicted signal yields for the 2017 data-taking period, corresponding to 41.3 fb$^{-1}$, after the application of each search requirement (cumulative) for various signal hypothesis. The requirements listed are presented as total efficiencies w.r.t. the previous selection step.

Observed 95% CL exclusion limits on the simplified models of chargino-pair production decaying via W bosons.

$+1\sigma$ observed 95% CL exclusion limits on the simplified models of chargino-pair production decaying via W bosons.

$-1\sigma$ observed 95% CL exclusion limits on the simplified models of chargino-pair production decaying via W bosons.

Expected 95% CL exclusion limits on the simplified models of chargino--neutralino production decaying via W and Z bosons.

$+1\sigma$ expected 95% CL exclusion limits on the simplified models of chargino--neutralino production decaying via W and Z bosons.

$-1\sigma$ expected 95% CL exclusion limits on the simplified models of chargino--neutralino production decaying via W and Z bosons.

Observed 95% CL exclusion limits on the simplified models of chargino--neutralino production decaying via W and Z bosons.

$+1\sigma$ observed 95% CL exclusion limits on the simplified models of chargino--neutralino production decaying via W and Z bosons.

$-1\sigma$ observed 95% CL exclusion limits on the simplified models of chargino--neutralino production decaying via W and Z bosons.

Expected 95% CL exclusion limits on the simplified models of chargino--neutralino production decaying via W and h bosons.

$+1\sigma$ expected 95% CL exclusion limits on the simplified models of chargino--neutralino production decaying via W and h bosons.

$-1\sigma$ expected 95% CL exclusion limits on the simplified models of chargino--neutralino production decaying via W and h bosons.

Observed 95% CL exclusion limits on the simplified models of chargino--neutralino production decaying via W and h bosons.

$+1\sigma$ observed 95% CL exclusion limits on the simplified models of chargino--neutralino production decaying via W and h bosons.

$-1\sigma$ observed 95% CL exclusion limits on the simplified models of chargino--neutralino production decaying via W and h bosons.

Expected 95% CL exclusion limits on the simplified models of higgsino GGM scenarios.

$+1\sigma$ expected 95% CL exclusion limits on the simplified models of higgsino GGM scenarios.

$-1\sigma$ expected 95% CL exclusion limits on the simplified models of higgsino GGM scenarios.

Observed 95% CL exclusion limits on the simplified models of higgsino GGM scenarios.

$+1\sigma$ observed 95% CL exclusion limits on the simplified models of higgsino GGM scenarios.

$-1\sigma$ observed 95% CL exclusion limits on the simplified models of higgsino GGM scenarios.

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

The analyses used in combination for each scenario to set limits in models of the production of $\tilde{\chi}_1^{+}\tilde{\chi}_{1}^{-}$.

Observed upper limit on the signal cross section in fb for chargino--neutralino production decaying via W and Z bosons.

The analyses used in combination for each scenario to set limits in models of chargino--neutralino production decaying via W and Z bosons.

Expected 95% CL exclusion limits on the simplified models of chargino--neutralino production decaying via W and Z bosons.

$+1\sigma$ expected 95% CL exclusion limits on the simplified models of chargino--neutralino production decaying via W and Z bosons.

$-1\sigma$ expected 95% CL exclusion limits on the simplified models of chargino--neutralino production decaying via W and Z bosons.

Observed 95% CL exclusion limits on the simplified models of chargino--neutralino production decaying via W and Z bosons.

$+1\sigma$ observed 95% CL exclusion limits on the simplified models of chargino--neutralino production decaying via W and Z bosons.

$-1\sigma$ observed 95% CL exclusion limits on the simplified models of chargino--neutralino production decaying via W and Z bosons.

Observed upper limit on the signal cross section in fb for chargino--neutralino production decaying via W and h bosons.

The analyses used in combination for each scenario to set limits in models of chargino--neutralino production decaying via W and h bosons.

Observed upper limit on the signal cross section in fb for higgsino GGM scenarios.

The analyses used in combination for each scenario to set limits in higgsino GGM scenarios.

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.

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.

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.

Parametrization of the $A_{X}$ factor, defined as the fraction of diphoton resonances satisfying the fiducial acceptance at the particle level, as function of $m_{X}$ obtained from the narrow width simulated signal samples produced in gluon fusion.

The correction factor, $C_{X}$, defined as the ratio between the number of reconstructed signal events passing the analysis cuts and the number of signal events at the particle level generated within the fiducial volume, and acceptance correction factor, $A_{X}$, defined as the fraction of diphoton resonances satisfying the fiducial acceptance at the particle level. Both are computed for NWA spin-0 models as a function of $m_{X}$.

Effect of event selections on a scalar MC signal sample generated for $m_{X}$ = 15 GeV and on the data. For the MC sample, the efficiencies are shown after applying event weights and a truth level filter that requires two photons with $p^{\gamma\gamma}_{T}>40$ GeV; for the data, the absolute yields are shown. The initial yields for data include a trigger preselection that is the OR of a list of single photon and diphoton triggers. The "2 $loose$ photons" step includes the kinematic acceptance cuts.

Parameterization of the Double Sided Crystal Ball function parameters describing the scalar mass resolution model as a function of $m_{X}$ [GeV].