Data compared to the SM prediction in each region of the positive \lam{} analysis, after the fit to data. The shaded bands represent the total post-fit uncertainty on the prediction.

Fit results in the ($\kappa_{Z}$, $\kappa_{W}$) plane, using the negative-$\lambda_{WZ}$ analysis. The results are overlaid with the confidence regions (shown in blue) from a separate fit using a combination of the Higgs boson measurements collected in Nature 607 (2022) 52. This fit assumes that all Higgs boson couplings besides the ones represented are positive, and that only SM particles contribute to loop processes. Confidence regions are constructed from the log-likelihood ratio $\Lambda_{LR}=-2\ln(L/L_{max})$, where $L_{\max}$ is the likelihood for the best-fit point, which is shown as a black dot for the VBF $WH$ analysis or a blue dot for the Higgs combination. The $1\sigma$, $2\sigma$, and $5\sigma$ regions are defined by $\Lambda_{LR}$ values smaller than 2.30, 6.18, and 28.7, respectively. The SM value is marked with a star, while green triangles mark the points with $\kappa_{Z} = \pm1$, $\kappa_{W} = \mp1$.

Post-fit SR yields of the high-mass channel. The upper panel shows the observed number of events, as well the post-fit background predictions in each region. The bottom panel shows the ratio of the observed data and the total background prediction. The shaded areas correspond to the total statistical and systematic uncertainties obtained after the fit and described in Section 6.

Post-fit SR yields of the high-mass channel. The upper panel shows the observed number of events, as well the post-fit background predictions in each region. The bottom panel shows the ratio of the observed data and the total background prediction. The shaded areas correspond to the total statistical and systematic uncertainties obtained after the fit and described in Section 6.

Post-fit SR yields of the high-mass channel. The upper panel shows the observed number of events, as well the post-fit background predictions in each region. The bottom panel shows the ratio of the observed data and the total background prediction. The shaded areas correspond to the total statistical and systematic uncertainties obtained after the fit and described in Section 6.

Post-fit SR yields of the high-mass channel. The upper panel shows the observed number of events, as well the post-fit background predictions in each region. The bottom panel shows the ratio of the observed data and the total background prediction. The shaded areas correspond to the total statistical and systematic uncertainties obtained after the fit and described in Section 6.

Post-fit SR yields of the high-mass channel. The upper panel shows the observed number of events, as well the post-fit background predictions in each region. The bottom panel shows the ratio of the observed data and the total background prediction. The shaded areas correspond to the total statistical and systematic uncertainties obtained after the fit and described in Section 6.

Post-fit SR yields of the high-mass channel. The upper panel shows the observed number of events, as well the post-fit background predictions in each region. The bottom panel shows the ratio of the observed data and the total background prediction. The shaded areas correspond to the total statistical and systematic uncertainties obtained after the fit and described in Section 6.

Post-fit SR yields of the high-mass channel. The upper panel shows the observed number of events, as well the post-fit background predictions in each region. The bottom panel shows the ratio of the observed data and the total background prediction. The shaded areas correspond to the total statistical and systematic uncertainties obtained after the fit and described in Section 6.

Pre-fit data and background (reweighted $2b$) predictions for each $4b$ SR $E_\text{T}^\text{miss}$ and $m_\text{eff}$ bin of the low-mass channel for the 2016 data-taking period. The bottom panel shows the significance of any differences between the observed $4b$ data and the background prediction. The $1\sigma$ and $2\sigma$ bands are shown in green and yellow, respectively. All systematics are included except the background normalization, which is 2.3%.

Pre-fit data and background (reweighted $2b$) predictions for each $4b$ SR $E_\text{T}^\text{miss}$ and $m_\text{eff}$ bin of the low-mass channel for the 2017 data-taking period. The bottom panel shows the significance of any differences between the observed $4b$ data and the background prediction. The $1\sigma$ and $2\sigma$ bands are shown in green and yellow, respectively. All systematics are included except the background normalization, which is 3.7%.

Pre-fit data and background (reweighted $2b$) predictions for each $4b$ SR $E_\text{T}^\text{miss}$ and $m_\text{eff}$ bin of the low-mass channel for the 2018 data-taking period. The bottom panel shows the significance of any differences between the observed $4b$ data and the background prediction. The $1\sigma$ and $2\sigma$ bands are shown in green and yellow, respectively. All systematics are included except the background normalization, which is 1.8%.

Exclusion limits of the low-mass and high-mass channels. The low-mass channel is used for $m_{\tilde{H}}<250$ GeV while the high-mass channel is used for $m_{\tilde{H}}\ge250$ GeV. The plot shows the observed (solid) and expected (dashed) 95% CL upper limits on the cross section of higgsino pair production, assuming a higgsino decay branching ratio of $\mathcal{B}(\tilde{H}\rightarrow h + \tilde{G})=100\%$. The theory cross section and its uncertainty are shown by the solid red line and red shading. Results from a previous ATLAS search using 24.3-36.1 fb$^{-1}$ [13] are shown by the solid (observed) and dashed (expected) blue lines. The bottom panel shows the ratio of the limits to the theory cross section. The phase space above the lines is excluded.

Exclusion limits of the low-mass and high-mass channels. The low-mass channel is used for $m_{\tilde{H}}<250$ GeV while the high-mass channel is used for $m_{\tilde{H}}\ge250$ GeV. The plot shows the observed (solid) and expected (dashed) 95% CL upper limits on the cross section of higgsino pair production, assuming a higgsino decay branching ratio of $\mathcal{B}(\tilde{H}\rightarrow h + \tilde{G})=100\%$. The theory cross section and its uncertainty are shown by the solid red line and red shading. Results from a previous ATLAS search using 24.3-36.1 fb$^{-1}$ [13] are shown by the solid (observed) and dashed (expected) blue lines. The bottom panel shows the ratio of the limits to the theory cross section. The phase space above the lines is excluded.

Exclusion limits of the low-mass and high-mass channels. The low-mass channel is used for $m_{\tilde{H}}<250$ GeV while the high-mass channel is used for $m_{\tilde{H}}\ge250$ GeV. The plot shows the observed (solid) and expected (dashed) 95% CL upper limits on the cross section of higgsino pair production, assuming a higgsino decay branching ratio of $\mathcal{B}(\tilde{H}\rightarrow h + \tilde{G})=100\%$. The theory cross section and its uncertainty are shown by the solid red line and red shading. Results from a previous ATLAS search using 24.3-36.1 fb$^{-1}$ [13] are shown by the solid (observed) and dashed (expected) blue lines. The bottom panel shows the ratio of the limits to the theory cross section. The phase space above the lines is excluded.

Exclusion limits of the low-mass and high-mass channels. The low-mass channel is used for $m_{\tilde{H}}<250$ GeV while the high-mass channel is used for $m_{\tilde{H}}\ge250$ GeV. The plot shows the observed (solid) and expected (dashed) 95% CL upper limits on the cross section of higgsino pair production, assuming a higgsino decay branching ratio of $\mathcal{B}(\tilde{H}\rightarrow h + \tilde{G})=100\%$. The theory cross section and its uncertainty are shown by the solid red line and red shading. Results from a previous ATLAS search using 24.3-36.1 fb$^{-1}$ [13] are shown by the solid (observed) and dashed (expected) blue lines. The bottom panel shows the ratio of the limits to the theory cross section. The phase space above the lines is excluded.

Exclusion limits of the low-mass and high-mass channels. The low-mass channel is used for $m_{\tilde{H}}<250$ GeV while the high-mass channel is used for $m_{\tilde{H}}\ge250$ GeV. The plot shows the observed (solid) and expected (dashed) 95% CL upper limits on the cross section of higgsino pair production, assuming a higgsino decay branching ratio of $\mathcal{B}(\tilde{H}\rightarrow h + \tilde{G})=100\%$. The theory cross section and its uncertainty are shown by the solid red line and red shading. Results from a previous ATLAS search using 24.3-36.1 fb$^{-1}$ [13] are shown by the solid (observed) and dashed (expected) blue lines. The bottom panel shows the ratio of the limits to the theory cross section. The phase space above the lines is excluded.

Exclusion limits of the low-mass and high-mass channels. The low-mass channel is used for $m_{\tilde{H}}<250$ GeV while the high-mass channel is used for $m_{\tilde{H}}\ge250$ GeV. The plot shows the observed (solid) and expected (dashed) 95% CL upper limits on the cross section of higgsino pair production, assuming a higgsino decay branching ratio of $\mathcal{B}(\tilde{H}\rightarrow h + \tilde{G})=100\%$. The theory cross section and its uncertainty are shown by the solid red line and red shading. Results from a previous ATLAS search using 24.3-36.1 fb$^{-1}$ [13] are shown by the solid (observed) and dashed (expected) blue lines. The bottom panel shows the ratio of the limits to the theory cross section. The phase space above the lines is excluded.

Exclusion limits of the low-mass and high-mass channels. The low-mass channel is used for $m_{\tilde{H}}<250$ GeV while the high-mass channel is used for $m_{\tilde{H}}\ge250$ GeV. The plot shows the observed (solid) and expected (dashed) 95% CL upper limits on the cross section of higgsino pair production, assuming a higgsino decay branching ratio of $\mathcal{B}(\tilde{H}\rightarrow h + \tilde{G})=100\%$. The theory cross section and its uncertainty are shown by the solid red line and red shading. Results from a previous ATLAS search using 24.3-36.1 fb$^{-1}$ [13] are shown by the solid (observed) and dashed (expected) blue lines. The bottom panel shows the ratio of the limits to the theory cross section. The phase space above the lines is excluded.

Exclusion limits of the low-mass and high-mass channels. The low-mass channel is used for $m_{\tilde{H}}<250$ GeV while the high-mass channel is used for $m_{\tilde{H}}\ge250$ GeV. The plot shows the observed (solid) and expected (dashed) 95% CL upper limits on the cross section of higgsino pair production, assuming a higgsino decay branching ratio of $\mathcal{B}(\tilde{H}\rightarrow h + \tilde{G})=100\%$. The theory cross section and its uncertainty are shown by the solid red line and red shading. Results from a previous ATLAS search using 24.3-36.1 fb$^{-1}$ [13] are shown by the solid (observed) and dashed (expected) blue lines. The bottom panel shows the ratio of the limits to the theory cross section. The phase space above the lines is excluded.

Exclusion limits of the low-mass and high-mass channels. The low-mass channel is used for $m_{\tilde{H}}<250$ GeV while the high-mass channel is used for $m_{\tilde{H}}\ge250$ GeV. The plot shows the observed (solid) and expected (dashed) 95% CL upper limits on the cross section of higgsino pair production, assuming a higgsino decay branching ratio of $\mathcal{B}(\tilde{H}\rightarrow h + \tilde{G})=100\%$. The theory cross section and its uncertainty are shown by the solid red line and red shading. Results from a previous ATLAS search using 24.3-36.1 fb$^{-1}$ [13] are shown by the solid (observed) and dashed (expected) blue lines. The bottom panel shows the ratio of the limits to the theory cross section. The phase space above the lines is excluded.

Exclusion limits of the low-mass and high-mass channels. The low-mass channel is used for $m_{\tilde{H}}<250$ GeV while the high-mass channel is used for $m_{\tilde{H}}\ge250$ GeV. The plot shows the observed (solid) and expected (dashed) 95% CL upper limits on the cross section of higgsino pair production, assuming a higgsino decay branching ratio of $\mathcal{B}(\tilde{H}\rightarrow h + \tilde{G})=100\%$. The theory cross section and its uncertainty are shown by the solid red line and red shading. Results from a previous ATLAS search using 24.3-36.1 fb$^{-1}$ [13] are shown by the solid (observed) and dashed (expected) blue lines. The bottom panel shows the ratio of the limits to the theory cross section. The phase space above the lines is excluded.

Exclusion limits of the low-mass and high-mass channels. The low-mass channel is used for $m_{\tilde{H}}<250$ GeV while the high-mass channel is used for $m_{\tilde{H}}\ge250$ GeV. The plot shows the 95% CL observed (solid) and expected (dashed) upper limits on $\mathcal{B}(\tilde{H}\rightarrow h + \tilde{G})$, assuming the theory cross section for higgsino pair production. The higgsinos are assumed to decay as $\tilde{H}\rightarrow h + \tilde{G}$ or $\tilde{H}\rightarrow Z + \tilde{G}$. The phase space above the lines is excluded.

Exclusion limits of the low-mass and high-mass channels. The low-mass channel is used for $m_{\tilde{H}}<250$ GeV while the high-mass channel is used for $m_{\tilde{H}}\ge250$ GeV. The plot shows the 95% CL observed (solid) and expected (dashed) upper limits on $\mathcal{B}(\tilde{H}\rightarrow h + \tilde{G})$, assuming the theory cross section for higgsino pair production. The higgsinos are assumed to decay as $\tilde{H}\rightarrow h + \tilde{G}$ or $\tilde{H}\rightarrow Z + \tilde{G}$. The phase space above the lines is excluded.

Exclusion limits of the low-mass and high-mass channels. The low-mass channel is used for $m_{\tilde{H}}<250$ GeV while the high-mass channel is used for $m_{\tilde{H}}\ge250$ GeV. The plot shows the 95% CL observed (solid) and expected (dashed) upper limits on $\mathcal{B}(\tilde{H}\rightarrow h + \tilde{G})$, assuming the theory cross section for higgsino pair production. The higgsinos are assumed to decay as $\tilde{H}\rightarrow h + \tilde{G}$ or $\tilde{H}\rightarrow Z + \tilde{G}$. The phase space above the lines is excluded.

Exclusion limits of the low-mass and high-mass channels. The low-mass channel is used for $m_{\tilde{H}}<250$ GeV while the high-mass channel is used for $m_{\tilde{H}}\ge250$ GeV. The plot shows the 95% CL observed (solid) and expected (dashed) upper limits on $\mathcal{B}(\tilde{H}\rightarrow h + \tilde{G})$, assuming the theory cross section for higgsino pair production. The higgsinos are assumed to decay as $\tilde{H}\rightarrow h + \tilde{G}$ or $\tilde{H}\rightarrow Z + \tilde{G}$. The phase space above the lines is excluded.

Exclusion limits of the low-mass and high-mass channels. The low-mass channel is used for $m_{\tilde{H}}<250$ GeV while the high-mass channel is used for $m_{\tilde{H}}\ge250$ GeV. The plot shows the 95% CL observed (solid) and expected (dashed) upper limits on $\mathcal{B}(\tilde{H}\rightarrow h + \tilde{G})$, assuming the theory cross section for higgsino pair production. The higgsinos are assumed to decay as $\tilde{H}\rightarrow h + \tilde{G}$ or $\tilde{H}\rightarrow Z + \tilde{G}$. The phase space above the lines is excluded.

Exclusion limits of the low-mass channel. The plot shows the observed (solid) and expected (dashed) 95% CL upper limits on the cross section of higgsino pair production, assuming a higgsino decay branching ratio of $\mathcal{B}(\tilde{H}\rightarrow h + \tilde{G})=100\%$. The theory cross section and its uncertainty are shown by the solid red line and red shading. The bottom panel shows the ratio of the limits to the theory cross section. The phase space above the lines is excluded.

Exclusion limits of the low-mass channel. The plot shows the observed (solid) and expected (dashed) 95% CL upper limits on the cross section of higgsino pair production, assuming a higgsino decay branching ratio of $\mathcal{B}(\tilde{H}\rightarrow h + \tilde{G})=100\%$. The theory cross section and its uncertainty are shown by the solid red line and red shading. The bottom panel shows the ratio of the limits to the theory cross section. The phase space above the lines is excluded.

Exclusion limits of the high-mass channel. The plot shows the observed (solid) and expected (dashed) 95% CL upper limits on the cross section of higgsino pair production, assuming a higgsino decay branching ratio of $\mathcal{B}(\tilde{H}\rightarrow h + \tilde{G})=100\%$. The theory cross section and its uncertainty are shown by the solid red line and red shading. The bottom panel shows the ratio of the limits to the theory cross section. The phase space above the lines is excluded.

Exclusion limits of the high-mass channel. The plot shows the observed (solid) and expected (dashed) 95% CL upper limits on the cross section of higgsino pair production, assuming a higgsino decay branching ratio of $\mathcal{B}(\tilde{H}\rightarrow h + \tilde{G})=100\%$. The theory cross section and its uncertainty are shown by the solid red line and red shading. The bottom panel shows the ratio of the limits to the theory cross section. The phase space above the lines is excluded.

Results of the background-only fit in the low-mass channel discovery region SR_LM_150. Both pre-fit and post-fit values are shown.

Results of the background-only fit in the low-mass channel discovery region SR_LM_300. Both pre-fit and post-fit values are shown.

The experimental efficiency of the low-mass channel for the exclusion and discovery signal regions as a function of higgsino mass. The experimental efficiency is defined as the number of events passing the detector-level event selections divided by the number of events passing the event selections for a perfect detector. The denominator is obtained by implementing particle-level event selections that emulate the detector-level selections. This treats the lack of availability of $b$-jet triggers as an inefficiency.

The particle-level acceptance for the low-mass exclusion and discovery signal regions, shown as a function of higgsino mass. The acceptance is defined as the fraction of signal events passing the particle-level event selection that emulates the detector-level selection. The acceptance calculation considers only those signal events where both higgsinos decay to Higgs bosons.

The experimental efficiency of the high-mass channel discovery regions as a function of higgsino mass. For each higgsino mass, the efficiency is shown for the SR-1 region corresponding to the mass. For masses above 1100 GeV, SR-1-1100 is used. The experimental efficiency is defined as the number of events passing the detector-level event selections divided by the number of events passing the event selections for a perfect detector. The denominator is obtained by implementing particle-level event selections that emulate the detector-level selections. The efficiency calculation considers only those signal events where both higgsinos decay to Higgs bosons.

The particle-level acceptance for the high-mass signal regions, shown as a function of higgsino mass. For each higgsino mass, the acceptance is shown for the SR-1 region corresponding to the mass. For masses above 1100 GeV, SR-1-1100 is used. The acceptance is defined as the fraction of signal events passing the particle-level event selection that emulates the detector-level selection. The acceptance calculation considers only those signal events where both higgsinos decay to Higgs bosons.

Cutflow for the low-mass channel for a representative 130 GeV signal. The preselection requires 4 or more jets with $p_\text{T}>25$ GeV and 2 or more $b$-jets with $p_\text{T}>25$ GeV. The $b$-jet cut requires 4 or more $b$-jets with $p_\text{T}>40$ GeV. As the samples are generated with $\mathcal{B}(\tilde{H}\rightarrow h\tilde{G})$=50%, $\mathcal{B}(\tilde{H}\rightarrow Z\tilde{G})$=50% to allow for both decays to be studied, the $hh$ events selection is used to select the events where each of the higgsinos decays to a Higgs boson. Expected yields are normalized to a luminosity of 139 fb$^{-1}$, with the availability of $b$-jet triggers lowering the luminosity to 126 fb$^{-1}$. All selections are cumulative, with the exception of the SR cuts, which are each applied separately.

Cutflow for the low-mass channel for a representative 150 GeV signal. The preselection requires 4 or more jets with $p_\text{T}>25$ GeV and 2 or more $b$-jets with $p_\text{T}>25$ GeV. The $b$-jet cut requires 4 or more $b$-jets with $p_\text{T}>40$ GeV. As the samples are generated with $\mathcal{B}(\tilde{H}\rightarrow h\tilde{G})$=50%, $\mathcal{B}(\tilde{H}\rightarrow Z\tilde{G})$=50% to allow for both decays to be studied, the $hh$ events selection is used to select the events where each of the higgsinos decays to a Higgs boson. Expected yields are normalized to a luminosity of 139 fb$^{-1}$, with the availability of $b$-jet triggers lowering the luminosity to 126 fb$^{-1}$. All selections are cumulative, with the exception of the SR cuts, which are each applied separately.

Cutflow for the low-mass channel for a representative 200 GeV signal. The preselection requires 4 or more jets with $p_\text{T}>25$ GeV and 2 or more $b$-jets with $p_\text{T}>25$ GeV. The $b$-jet cut requires 4 or more $b$-jets with $p_\text{T}>40$ GeV. As the samples are generated with $\mathcal{B}(\tilde{H}\rightarrow h\tilde{G})$=50%, $\mathcal{B}(\tilde{H}\rightarrow Z\tilde{G})$=50% to allow for both decays to be studied, the $hh$ events selection is used to select the events where each of the higgsinos decays to a Higgs boson. Expected yields are normalized to a luminosity of 139 fb$^{-1}$, with the availability of $b$-jet triggers lowering the luminosity to 126 fb$^{-1}$. All selections are cumulative, with the exception of the SR cuts, which are each applied separately.

Cutflow for the low-mass channel for a representative 250 GeV signal. The preselection requires 4 or more jets with $p_\text{T}>25$ GeV and 2 or more $b$-jets with $p_\text{T}>25$ GeV. The $b$-jet cut requires 4 or more $b$-jets with $p_\text{T}>40$ GeV. As the samples are generated with $\mathcal{B}(\tilde{H}\rightarrow h\tilde{G})$=50%, $\mathcal{B}(\tilde{H}\rightarrow Z\tilde{G})$=50% to allow for both decays to be studied, the $hh$ events selection is used to select the events where each of the higgsinos decays to a Higgs boson. Expected yields are normalized to a luminosity of 139 fb$^{-1}$, with the availability of $b$-jet triggers lowering the luminosity to 126 fb$^{-1}$. All selections are cumulative, with the exception of the SR cuts, which are each applied separately.

Cutflow for the low-mass channel for a representative 300 GeV signal. The preselection requires 4 or more jets with $p_\text{T}>25$ GeV and 2 or more $b$-jets with $p_\text{T}>25$ GeV. The $b$-jet cut requires 4 or more $b$-jets with $p_\text{T}>40$ GeV. As the samples are generated with $\mathcal{B}(\tilde{H}\rightarrow h\tilde{G})$=50%, $\mathcal{B}(\tilde{H}\rightarrow Z\tilde{G})$=50% to allow for both decays to be studied, the $hh$ events selection is used to select the events where each of the higgsinos decays to a Higgs boson. Expected yields are normalized to a luminosity of 139 fb$^{-1}$, with the availability of $b$-jet triggers lowering the luminosity to 126 fb$^{-1}$. All selections are cumulative, with the exception of the SR cuts, which are each applied separately.

Cutflow for the low-mass channel for a representative 400 GeV signal. The preselection requires 4 or more jets with $p_\text{T}>25$ GeV and 2 or more $b$-jets with $p_\text{T}>25$ GeV. The $b$-jet cut requires 4 or more $b$-jets with $p_\text{T}>40$ GeV. As the samples are generated with $\mathcal{B}(\tilde{H}\rightarrow h\tilde{G})$=50%, $\mathcal{B}(\tilde{H}\rightarrow Z\tilde{G})$=50% to allow for both decays to be studied, the $hh$ events selection is used to select the events where each of the higgsinos decays to a Higgs boson. Expected yields are normalized to a luminosity of 139 fb$^{-1}$, with the availability of $b$-jet triggers lowering the luminosity to 126 fb$^{-1}$. All selections are cumulative, with the exception of the SR cuts, which are each applied separately.

Cutflow for the low-mass channel for a representative 500 GeV signal. The preselection requires 4 or more jets with $p_\text{T}>25$ GeV and 2 or more $b$-jets with $p_\text{T}>25$ GeV. The $b$-jet cut requires 4 or more $b$-jets with $p_\text{T}>40$ GeV. As the samples are generated with $\mathcal{B}(\tilde{H}\rightarrow h\tilde{G})$=50%, $\mathcal{B}(\tilde{H}\rightarrow Z\tilde{G})$=50% to allow for both decays to be studied, the $hh$ events selection is used to select the events where each of the higgsinos decays to a Higgs boson. Expected yields are normalized to a luminosity of 139 fb$^{-1}$, with the availability of $b$-jet triggers lowering the luminosity to 126 fb$^{-1}$. All selections are cumulative, with the exception of the SR cuts, which are each applied separately.

Cutflow for the low-mass channel for a representative 600 GeV signal. The preselection requires 4 or more jets with $p_\text{T}>25$ GeV and 2 or more $b$-jets with $p_\text{T}>25$ GeV. The $b$-jet cut requires 4 or more $b$-jets with $p_\text{T}>40$ GeV. As the samples are generated with $\mathcal{B}(\tilde{H}\rightarrow h\tilde{G})$=50%, $\mathcal{B}(\tilde{H}\rightarrow Z\tilde{G})$=50% to allow for both decays to be studied, the $hh$ events selection is used to select the events where each of the higgsinos decays to a Higgs boson. Expected yields are normalized to a luminosity of 139 fb$^{-1}$, with the availability of $b$-jet triggers lowering the luminosity to 126 fb$^{-1}$. All selections are cumulative, with the exception of the SR cuts, which are each applied separately.

Cutflow for the low-mass channel for a representative 700 GeV signal. The preselection requires 4 or more jets with $p_\text{T}>25$ GeV and 2 or more $b$-jets with $p_\text{T}>25$ GeV. The $b$-jet cut requires 4 or more $b$-jets with $p_\text{T}>40$ GeV. As the samples are generated with $\mathcal{B}(\tilde{H}\rightarrow h\tilde{G})$=50%, $\mathcal{B}(\tilde{H}\rightarrow Z\tilde{G})$=50% to allow for both decays to be studied, the $hh$ events selection is used to select the events where each of the higgsinos decays to a Higgs boson. Expected yields are normalized to a luminosity of 139 fb$^{-1}$, with the availability of $b$-jet triggers lowering the luminosity to 126 fb$^{-1}$. All selections are cumulative, with the exception of the SR cuts, which are each applied separately.

Cutflow for the low-mass channel for a representative 800 GeV signal. The preselection requires 4 or more jets with $p_\text{T}>25$ GeV and 2 or more $b$-jets with $p_\text{T}>25$ GeV. The $b$-jet cut requires 4 or more $b$-jets with $p_\text{T}>40$ GeV. As the samples are generated with $\mathcal{B}(\tilde{H}\rightarrow h\tilde{G})$=50%, $\mathcal{B}(\tilde{H}\rightarrow Z\tilde{G})$=50% to allow for both decays to be studied, the $hh$ events selection is used to select the events where each of the higgsinos decays to a Higgs boson. Expected yields are normalized to a luminosity of 139 fb$^{-1}$, with the availability of $b$-jet triggers lowering the luminosity to 126 fb$^{-1}$. All selections are cumulative, with the exception of the SR cuts, which are each applied separately.

Cutflow for the low-mass channel for a representative 900 GeV signal. The preselection requires 4 or more jets with $p_\text{T}>25$ GeV and 2 or more $b$-jets with $p_\text{T}>25$ GeV. The $b$-jet cut requires 4 or more $b$-jets with $p_\text{T}>40$ GeV. As the samples are generated with $\mathcal{B}(\tilde{H}\rightarrow h\tilde{G})$=50%, $\mathcal{B}(\tilde{H}\rightarrow Z\tilde{G})$=50% to allow for both decays to be studied, the $hh$ events selection is used to select the events where each of the higgsinos decays to a Higgs boson. Expected yields are normalized to a luminosity of 139 fb$^{-1}$, with the availability of $b$-jet triggers lowering the luminosity to 126 fb$^{-1}$. All selections are cumulative, with the exception of the SR cuts, which are each applied separately.

Cutflow for the low-mass channel for a representative 1000 GeV signal. The preselection requires 4 or more jets with $p_\text{T}>25$ GeV and 2 or more $b$-jets with $p_\text{T}>25$ GeV. The $b$-jet cut requires 4 or more $b$-jets with $p_\text{T}>40$ GeV. As the samples are generated with $\mathcal{B}(\tilde{H}\rightarrow h\tilde{G})$=50%, $\mathcal{B}(\tilde{H}\rightarrow Z\tilde{G})$=50% to allow for both decays to be studied, the $hh$ events selection is used to select the events where each of the higgsinos decays to a Higgs boson. Expected yields are normalized to a luminosity of 139 fb$^{-1}$, with the availability of $b$-jet triggers lowering the luminosity to 126 fb$^{-1}$. All selections are cumulative, with the exception of the SR cuts, which are each applied separately.

Cutflow for the low-mass channel for a representative 1100 GeV signal. The preselection requires 4 or more jets with $p_\text{T}>25$ GeV and 2 or more $b$-jets with $p_\text{T}>25$ GeV. The $b$-jet cut requires 4 or more $b$-jets with $p_\text{T}>40$ GeV. As the samples are generated with $\mathcal{B}(\tilde{H}\rightarrow h\tilde{G})$=50%, $\mathcal{B}(\tilde{H}\rightarrow Z\tilde{G})$=50% to allow for both decays to be studied, the $hh$ events selection is used to select the events where each of the higgsinos decays to a Higgs boson. Expected yields are normalized to a luminosity of 139 fb$^{-1}$, with the availability of $b$-jet triggers lowering the luminosity to 126 fb$^{-1}$. All selections are cumulative, with the exception of the SR cuts, which are each applied separately.

Cutflow for the high-mass channel for a representative 200 GeV signal. The preselection requires 4 or more jets with $p_\text{T}>25$ GeV and 2 or more $b$-jets with $p_\text{T}>25$ GeV. As the samples are generated with $\mathcal{B}(\tilde{H}\rightarrow h\tilde{G})$=50%, $\mathcal{B}(\tilde{H}\rightarrow Z\tilde{G})$=50% to allow for both decays to be studied, the $hh$ events selection is used to select the events where each of the higgsinos decays to a Higgs boson. Expected yields are normalized to a luminosity of 139 fb$^{-1}$. All selections are cumulative, with the exception of the SR cuts, which are each applied separately.

Cutflow for the high-mass channel for a representative 250 GeV signal. The preselection requires 4 or more jets with $p_\text{T}>25$ GeV and 2 or more $b$-jets with $p_\text{T}>25$ GeV. As the samples are generated with $\mathcal{B}(\tilde{H}\rightarrow h\tilde{G})$=50%, $\mathcal{B}(\tilde{H}\rightarrow Z\tilde{G})$=50% to allow for both decays to be studied, the $hh$ events selection is used to select the events where each of the higgsinos decays to a Higgs boson. Expected yields are normalized to a luminosity of 139 fb$^{-1}$. All selections are cumulative, with the exception of the SR cuts, which are each applied separately.

Cutflow for the high-mass channel for a representative 300 GeV signal. The preselection requires 4 or more jets with $p_\text{T}>25$ GeV and 2 or more $b$-jets with $p_\text{T}>25$ GeV. As the samples are generated with $\mathcal{B}(\tilde{H}\rightarrow h\tilde{G})$=50%, $\mathcal{B}(\tilde{H}\rightarrow Z\tilde{G})$=50% to allow for both decays to be studied, the $hh$ events selection is used to select the events where each of the higgsinos decays to a Higgs boson. Expected yields are normalized to a luminosity of 139 fb$^{-1}$. All selections are cumulative, with the exception of the SR cuts, which are each applied separately.

Cutflow for the high-mass channel for a representative 400 GeV signal. The preselection requires 4 or more jets with $p_\text{T}>25$ GeV and 2 or more $b$-jets with $p_\text{T}>25$ GeV. As the samples are generated with $\mathcal{B}(\tilde{H}\rightarrow h\tilde{G})$=50%, $\mathcal{B}(\tilde{H}\rightarrow Z\tilde{G})$=50% to allow for both decays to be studied, the $hh$ events selection is used to select the events where each of the higgsinos decays to a Higgs boson. Expected yields are normalized to a luminosity of 139 fb$^{-1}$. All selections are cumulative, with the exception of the SR cuts, which are each applied separately.

Cutflow for the high-mass channel for a representative 500 GeV signal. The preselection requires 4 or more jets with $p_\text{T}>25$ GeV and 2 or more $b$-jets with $p_\text{T}>25$ GeV. As the samples are generated with $\mathcal{B}(\tilde{H}\rightarrow h\tilde{G})$=50%, $\mathcal{B}(\tilde{H}\rightarrow Z\tilde{G})$=50% to allow for both decays to be studied, the $hh$ events selection is used to select the events where each of the higgsinos decays to a Higgs boson. Expected yields are normalized to a luminosity of 139 fb$^{-1}$. All selections are cumulative, with the exception of the SR cuts, which are each applied separately.

Cutflow for the high-mass channel for a representative 600 GeV signal. The preselection requires 4 or more jets with $p_\text{T}>25$ GeV and 2 or more $b$-jets with $p_\text{T}>25$ GeV. As the samples are generated with $\mathcal{B}(\tilde{H}\rightarrow h\tilde{G})$=50%, $\mathcal{B}(\tilde{H}\rightarrow Z\tilde{G})$=50% to allow for both decays to be studied, the $hh$ events selection is used to select the events where each of the higgsinos decays to a Higgs boson. Expected yields are normalized to a luminosity of 139 fb$^{-1}$. All selections are cumulative, with the exception of the SR cuts, which are each applied separately.

Cutflow for the high-mass channel for a representative 700 GeV signal. The preselection requires 4 or more jets with $p_\text{T}>25$ GeV and 2 or more $b$-jets with $p_\text{T}>25$ GeV. As the samples are generated with $\mathcal{B}(\tilde{H}\rightarrow h\tilde{G})$=50%, $\mathcal{B}(\tilde{H}\rightarrow Z\tilde{G})$=50% to allow for both decays to be studied, the $hh$ events selection is used to select the events where each of the higgsinos decays to a Higgs boson. Expected yields are normalized to a luminosity of 139 fb$^{-1}$. All selections are cumulative, with the exception of the SR cuts, which are each applied separately.

Cutflow for the high-mass channel for a representative 800 GeV signal. The preselection requires 4 or more jets with $p_\text{T}>25$ GeV and 2 or more $b$-jets with $p_\text{T}>25$ GeV. As the samples are generated with $\mathcal{B}(\tilde{H}\rightarrow h\tilde{G})$=50%, $\mathcal{B}(\tilde{H}\rightarrow Z\tilde{G})$=50% to allow for both decays to be studied, the $hh$ events selection is used to select the events where each of the higgsinos decays to a Higgs boson. Expected yields are normalized to a luminosity of 139 fb$^{-1}$. All selections are cumulative, with the exception of the SR cuts, which are each applied separately.

Cutflow for the high-mass channel for a representative 900 GeV signal. The preselection requires 4 or more jets with $p_\text{T}>25$ GeV and 2 or more $b$-jets with $p_\text{T}>25$ GeV. As the samples are generated with $\mathcal{B}(\tilde{H}\rightarrow h\tilde{G})$=50%, $\mathcal{B}(\tilde{H}\rightarrow Z\tilde{G})$=50% to allow for both decays to be studied, the $hh$ events selection is used to select the events where each of the higgsinos decays to a Higgs boson. Expected yields are normalized to a luminosity of 139 fb$^{-1}$. All selections are cumulative, with the exception of the SR cuts, which are each applied separately.

Cutflow for the high-mass channel for a representative 1000 GeV signal. The preselection requires 4 or more jets with $p_\text{T}>25$ GeV and 2 or more $b$-jets with $p_\text{T}>25$ GeV. As the samples are generated with $\mathcal{B}(\tilde{H}\rightarrow h\tilde{G})$=50%, $\mathcal{B}(\tilde{H}\rightarrow Z\tilde{G})$=50% to allow for both decays to be studied, the $hh$ events selection is used to select the events where each of the higgsinos decays to a Higgs boson. Expected yields are normalized to a luminosity of 139 fb$^{-1}$. All selections are cumulative, with the exception of the SR cuts, which are each applied separately.

Cutflow for the high-mass channel for a representative 1100 GeV signal. The preselection requires 4 or more jets with $p_\text{T}>25$ GeV and 2 or more $b$-jets with $p_\text{T}>25$ GeV. As the samples are generated with $\mathcal{B}(\tilde{H}\rightarrow h\tilde{G})$=50%, $\mathcal{B}(\tilde{H}\rightarrow Z\tilde{G})$=50% to allow for both decays to be studied, the $hh$ events selection is used to select the events where each of the higgsinos decays to a Higgs boson. Expected yields are normalized to a luminosity of 139 fb$^{-1}$. All selections are cumulative, with the exception of the SR cuts, which are each applied separately.

Cutflow for the high-mass channel for a representative 1200 GeV signal. The preselection requires 4 or more jets with $p_\text{T}>25$ GeV and 2 or more $b$-jets with $p_\text{T}>25$ GeV. As the samples are generated with $\mathcal{B}(\tilde{H}\rightarrow h\tilde{G})$=50%, $\mathcal{B}(\tilde{H}\rightarrow Z\tilde{G})$=50% to allow for both decays to be studied, the $hh$ events selection is used to select the events where each of the higgsinos decays to a Higgs boson. Expected yields are normalized to a luminosity of 139 fb$^{-1}$. All selections are cumulative, with the exception of the SR cuts, which are each applied separately.

Cutflow for the high-mass channel for a representative 1300 GeV signal. The preselection requires 4 or more jets with $p_\text{T}>25$ GeV and 2 or more $b$-jets with $p_\text{T}>25$ GeV. As the samples are generated with $\mathcal{B}(\tilde{H}\rightarrow h\tilde{G})$=50%, $\mathcal{B}(\tilde{H}\rightarrow Z\tilde{G})$=50% to allow for both decays to be studied, the $hh$ events selection is used to select the events where each of the higgsinos decays to a Higgs boson. Expected yields are normalized to a luminosity of 139 fb$^{-1}$. All selections are cumulative, with the exception of the SR cuts, which are each applied separately.

Cutflow for the high-mass channel for a representative 1400 GeV signal. The preselection requires 4 or more jets with $p_\text{T}>25$ GeV and 2 or more $b$-jets with $p_\text{T}>25$ GeV. As the samples are generated with $\mathcal{B}(\tilde{H}\rightarrow h\tilde{G})$=50%, $\mathcal{B}(\tilde{H}\rightarrow Z\tilde{G})$=50% to allow for both decays to be studied, the $hh$ events selection is used to select the events where each of the higgsinos decays to a Higgs boson. Expected yields are normalized to a luminosity of 139 fb$^{-1}$. All selections are cumulative, with the exception of the SR cuts, which are each applied separately.

Cutflow for the high-mass channel for a representative 1500 GeV signal. The preselection requires 4 or more jets with $p_\text{T}>25$ GeV and 2 or more $b$-jets with $p_\text{T}>25$ GeV. As the samples are generated with $\mathcal{B}(\tilde{H}\rightarrow h\tilde{G})$=50%, $\mathcal{B}(\tilde{H}\rightarrow Z\tilde{G})$=50% to allow for both decays to be studied, the $hh$ events selection is used to select the events where each of the higgsinos decays to a Higgs boson. Expected yields are normalized to a luminosity of 139 fb$^{-1}$. All selections are cumulative, with the exception of the SR cuts, which are each applied separately.

Expected (dashed black line) and observed (solid red line) 95% CL exclusion limits on the higgsino simplified model being considered. These are shown with $\pm 1\sigma_{\mathrm{exp}}$ (yellow band) from experimental systematic and statistical uncertainties, and with $\pm 1\sigma_{\mathrm{theory}}^{\mathrm{SUSY}}$ (red dotted lines) from signal cross-section uncertainties, respectively. The limits set by the latest ATLAS searches using the soft lepton and disappearing track signatures are illustrated by the blue and green regions, respectively, while the limit imposed by the LEP experiments is shown in gray. The dot-dashed gray line indicates the predicted mass-splitting for the pure higgsino scenario.

Expected (dashed black line) and observed (solid red line) 95% CL exclusion limits on the higgsino simplified model being considered. These are shown with $\pm 1\sigma_{\mathrm{exp}}$ (yellow band) from experimental systematic and statistical uncertainties, and with $\pm 1\sigma_{\mathrm{theory}}^{\mathrm{SUSY}}$ (red dotted lines) from signal cross-section uncertainties, respectively. The limits set by the latest ATLAS searches using the soft lepton and disappearing track signatures are illustrated by the blue and green regions, respectively, while the limit imposed by the LEP experiments is shown in gray. The dot-dashed gray line indicates the predicted mass-splitting for the pure higgsino scenario.

Expected (dashed black line) and observed (solid red line) 95% CL exclusion limits on the higgsino simplified model being considered. These are shown with $\pm 1\sigma_{\mathrm{exp}}$ (yellow band) from experimental systematic and statistical uncertainties, and with $\pm 1\sigma_{\mathrm{theory}}^{\mathrm{SUSY}}$ (red dotted lines) from signal cross-section uncertainties, respectively. The limits set by the latest ATLAS searches using the soft lepton and disappearing track signatures are illustrated by the blue and green regions, respectively, while the limit imposed by the LEP experiments is shown in gray. The dot-dashed gray line indicates the predicted mass-splitting for the pure higgsino scenario.

Expected (dashed black line) and observed (solid red line) 95% CL exclusion limits on the higgsino simplified model being considered. These are shown with $\pm 1\sigma_{\mathrm{exp}}$ (yellow band) from experimental systematic and statistical uncertainties, and with $\pm 1\sigma_{\mathrm{theory}}^{\mathrm{SUSY}}$ (red dotted lines) from signal cross-section uncertainties, respectively. The limits set by the latest ATLAS searches using the soft lepton and disappearing track signatures are illustrated by the blue and green regions, respectively, while the limit imposed by the LEP experiments is shown in gray. The dot-dashed gray line indicates the predicted mass-splitting for the pure higgsino scenario.

Expected (dashed black line) and observed (solid red line) 95% CL exclusion limits on the higgsino simplified model being considered. These are shown with $\pm 1\sigma_{\mathrm{exp}}$ (yellow band) from experimental systematic and statistical uncertainties, and with $\pm 1\sigma_{\mathrm{theory}}^{\mathrm{SUSY}}$ (red dotted lines) from signal cross-section uncertainties, respectively. The limits set by the latest ATLAS searches using the soft lepton and disappearing track signatures are illustrated by the blue and green regions, respectively, while the limit imposed by the LEP experiments is shown in gray. The dot-dashed gray line indicates the predicted mass-splitting for the pure higgsino scenario.

Expected (dashed black line) and observed (solid red line) 95% CL exclusion limits on the higgsino simplified model being considered. These are shown with $\pm 1\sigma_{\mathrm{exp}}$ (yellow band) from experimental systematic and statistical uncertainties, and with $\pm 1\sigma_{\mathrm{theory}}^{\mathrm{SUSY}}$ (red dotted lines) from signal cross-section uncertainties, respectively. The limits set by the latest ATLAS searches using the soft lepton and disappearing track signatures are illustrated by the blue and green regions, respectively, while the limit imposed by the LEP experiments is shown in gray. The dot-dashed gray line indicates the predicted mass-splitting for the pure higgsino scenario.

Expected (dashed black line) and observed (solid red line) 95% CL exclusion limits on the higgsino simplified model being considered. These are shown with $\pm 1\sigma_{\mathrm{exp}}$ (yellow band) from experimental systematic and statistical uncertainties, and with $\pm 1\sigma_{\mathrm{theory}}^{\mathrm{SUSY}}$ (red dotted lines) from signal cross-section uncertainties, respectively. The limits set by the latest ATLAS searches using the soft lepton and disappearing track signatures are illustrated by the blue and green regions, respectively, while the limit imposed by the LEP experiments is shown in gray. The dot-dashed gray line indicates the predicted mass-splitting for the pure higgsino scenario.

Expected (dashed black line) and observed (solid red line) 95% CL exclusion limits on the higgsino simplified model being considered. These are shown with $\pm 1\sigma_{\mathrm{exp}}$ (yellow band) from experimental systematic and statistical uncertainties, and with $\pm 1\sigma_{\mathrm{theory}}^{\mathrm{SUSY}}$ (red dotted lines) from signal cross-section uncertainties, respectively. The limits set by the latest ATLAS searches using the soft lepton and disappearing track signatures are illustrated by the blue and green regions, respectively, while the limit imposed by the LEP experiments is shown in gray. The dot-dashed gray line indicates the predicted mass-splitting for the pure higgsino scenario.

Expected (dashed black line) and observed (solid red line) 95% CL exclusion limits on the higgsino simplified model being considered. These are shown with $\pm 1\sigma_{\mathrm{exp}}$ (yellow band) from experimental systematic and statistical uncertainties, and with $\pm 1\sigma_{\mathrm{theory}}^{\mathrm{SUSY}}$ (red dotted lines) from signal cross-section uncertainties, respectively. The limits set by the latest ATLAS searches using the soft lepton and disappearing track signatures are illustrated by the blue and green regions, respectively, while the limit imposed by the LEP experiments is shown in gray. The dot-dashed gray line indicates the predicted mass-splitting for the pure higgsino scenario.

Expected (dashed black line) and observed (solid red line) 95% CL exclusion limits on the higgsino simplified model being considered. These are shown with $\pm 1\sigma_{\mathrm{exp}}$ (yellow band) from experimental systematic and statistical uncertainties, and with $\pm 1\sigma_{\mathrm{theory}}^{\mathrm{SUSY}}$ (red dotted lines) from signal cross-section uncertainties, respectively. The limits set by the latest ATLAS searches using the soft lepton and disappearing track signatures are illustrated by the blue and green regions, respectively, while the limit imposed by the LEP experiments is shown in gray. The dot-dashed gray line indicates the predicted mass-splitting for the pure higgsino scenario.

Expected (dashed black line) and observed (solid red line) 95% CL exclusion limits on the higgsino simplified model being considered. These are shown with $\pm 1\sigma_{\mathrm{exp}}$ (yellow band) from experimental systematic and statistical uncertainties, and with $\pm 1\sigma_{\mathrm{theory}}^{\mathrm{SUSY}}$ (red dotted lines) from signal cross-section uncertainties, respectively. The limits set by the latest ATLAS searches using the soft lepton and disappearing track signatures are illustrated by the blue and green regions, respectively, while the limit imposed by the LEP experiments is shown in gray. The dot-dashed gray line indicates the predicted mass-splitting for the pure higgsino scenario.

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

Absolute differential cross-section as a function of RC large-R jet $\tau_{21}$ at particle level in the $\ell$+jets channel. The measured differential cross-section is compared with the prediction obtained with the Powheg+Pythia8 Monte Carlo generator.

Covariance matrix of the Absolute differential cross-section as function of RC large-R jet $\tau_{21}$ at particle level in the $\ell$+jets channel, accounting for the statistical uncertainty.

Relative differential cross-section as a function of RC large-R jet $\tau_{21}$ at particle level in the $\ell$+jets channel. The measured differential cross-section is compared with the prediction obtained with the Powheg+Pythia8 Monte Carlo generator.

Absolute differential cross-section as a function of RC large-R jet $\tau_{3}$ at particle level in the $\ell$+jets channel. The measured differential cross-section is compared with the prediction obtained with the Powheg+Pythia8 Monte Carlo generator.

Covariance matrix of the Absolute differential cross-section as function of RC large-R jet $\tau_{3}$ at particle level in the $\ell$+jets channel, accounting for the statistical uncertainty.

Relative differential cross-section as a function of RC large-R jet $\tau_{3}$ at particle level in the $\ell$+jets channel. The measured differential cross-section is compared with the prediction obtained with the Powheg+Pythia8 Monte Carlo generator.

Absolute differential cross-section as a function of RC large-R jet ECF2 at particle level in the $\ell$+jets channel. The measured differential cross-section is compared with the prediction obtained with the Powheg+Pythia8 Monte Carlo generator.

Covariance matrix of the Absolute differential cross-section as function of RC large-R jet ECF2 at particle level in the $\ell$+jets channel, accounting for the statistical uncertainty.

Relative differential cross-section as a function of RC large-R jet ECF2 at particle level in the $\ell$+jets channel. The measured differential cross-section is compared with the prediction obtained with the Powheg+Pythia8 Monte Carlo generator.

Absolute differential cross-section as a function of RC large-R jet $D_{2}$ at particle level in the $\ell$+jets channel. The measured differential cross-section is compared with the prediction obtained with the Powheg+Pythia8 Monte Carlo generator.

Covariance matrix of the Absolute differential cross-section as function of RC large-R jet $D_{2}$ at particle level in the $\ell$+jets channel, accounting for the statistical uncertainty.

Relative differential cross-section as a function of RC large-R jet $D_{2}$ at particle level in the $\ell$+jets channel. The measured differential cross-section is compared with the prediction obtained with the Powheg+Pythia8 Monte Carlo generator.

Absolute differential cross-section as a function of RC large-R jet $C_{3}$ at particle level in the $\ell$+jets channel. The measured differential cross-section is compared with the prediction obtained with the Powheg+Pythia8 Monte Carlo generator.

Covariance matrix of the Absolute differential cross-section as function of RC large-R jet $C_{3}$ at particle level in the $\ell$+jets channel, accounting for the statistical uncertainty.

Relative differential cross-section as a function of RC large-R jet $C_{3}$ at particle level in the $\ell$+jets channel. The measured differential cross-section is compared with the prediction obtained with the Powheg+Pythia8 Monte Carlo generator.

Absolute differential cross-section as a function of RC large-R jet $p_\mathrm{T}^{d,*}$ (jet $p_T$ dispersion) at particle level in the $\ell$+jets channel. The measured differential cross-section is compared with the prediction obtained with the Powheg+Pythia8 Monte Carlo generator.

Covariance matrix of the Absolute differential cross-section as function of RC large-R jet $p_\mathrm{T}^{d,*}$ (jet $p_T$ dispersion) at particle level in the $\ell$+jets channel, accounting for the statistical uncertainty.

Relative differential cross-section as a function of RC large-R jet $p_\mathrm{T}^{d,*}$ (jet $p_T$ dispersion) at particle level in the $\ell$+jets channel. The measured differential cross-section is compared with the prediction obtained with the Powheg+Pythia8 Monte Carlo generator.

Absolute differential cross-section as a function of RC large-R jet LHA at particle level in the $\ell$+jets channel. The measured differential cross-section is compared with the prediction obtained with the Powheg+Pythia8 Monte Carlo generator.

Covariance matrix of the Absolute differential cross-section as function of RC large-R jet LHA at particle level in the $\ell$+jets channel, accounting for the statistical uncertainty.

Relative differential cross-section as a function of RC large-R jet LHA at particle level in the $\ell$+jets channel. The measured differential cross-section is compared with the prediction obtained with the Powheg+Pythia8 Monte Carlo generator.

Absolute double-differential cross-section as a function of RC large-R jet $D_{2}$ vs $m^{top}$ at particle level in the $\ell$+jets channel in 122.5 GeV < $m^{top}$ < 157.5 GeV. The measured differential cross-section is compared with the prediction obtained with the Powheg+Pythia8 Monte Carlo generator.

Absolute double-differential cross-section as a function of RC large-R jet $D_{2}$ vs $m^{top}$ at particle level in the $\ell$+jets channel in 157.5 GeV < $m^{top}$ < 187.5 GeV. The measured differential cross-section is compared with the prediction obtained with the Powheg+Pythia8 Monte Carlo generator.

Absolute double-differential cross-section as a function of RC large-R jet $D_{2}$ vs $m^{top}$ at particle level in the $\ell$+jets channel in 187.5 GeV < $m^{top}$ < 222.5 GeV. The measured differential cross-section is compared with the prediction obtained with the Powheg+Pythia8 Monte Carlo generator.

Covariance matrix between the Absolute double-differential cross-section as function of RC large-R jet $D_{2}$ vs $m^{top}$ in 122.5 GeV < $m^{top}$ < 157.5 GeV and the Absolute double-differential cross-section as function of RC large-R jet $D_{2}$ vs $m^{top}$ in 122.5 GeV < $m^{top}$ < 157.5 GeV at particle level in the $\ell$+jets channel, accounting for the statistical uncertainty.

Covariance matrix between the Absolute double-differential cross-section as function of RC large-R jet $D_{2}$ vs $m^{top}$ in 157.5 GeV < $m^{top}$ < 187.5 GeV and the Absolute double-differential cross-section as function of RC large-R jet $D_{2}$ vs $m^{top}$ in 122.5 GeV < $m^{top}$ < 157.5 GeV at particle level in the $\ell$+jets channel, accounting for the statistical uncertainty.

Covariance matrix between the Absolute double-differential cross-section as function of RC large-R jet $D_{2}$ vs $m^{top}$ in 157.5 GeV < $m^{top}$ < 187.5 GeV and the Absolute double-differential cross-section as function of RC large-R jet $D_{2}$ vs $m^{top}$ in 157.5 GeV < $m^{top}$ < 187.5 GeV at particle level in the $\ell$+jets channel, accounting for the statistical uncertainty.

Covariance matrix between the Absolute double-differential cross-section as function of RC large-R jet $D_{2}$ vs $m^{top}$ in 187.5 GeV < $m^{top}$ < 222.5 GeV and the Absolute double-differential cross-section as function of RC large-R jet $D_{2}$ vs $m^{top}$ in 122.5 GeV < $m^{top}$ < 157.5 GeV at particle level in the $\ell$+jets channel, accounting for the statistical uncertainty.

Covariance matrix between the Absolute double-differential cross-section as function of RC large-R jet $D_{2}$ vs $m^{top}$ in 187.5 GeV < $m^{top}$ < 222.5 GeV and the Absolute double-differential cross-section as function of RC large-R jet $D_{2}$ vs $m^{top}$ in 157.5 GeV < $m^{top}$ < 187.5 GeV at particle level in the $\ell$+jets channel, accounting for the statistical uncertainty.

Covariance matrix between the Absolute double-differential cross-section as function of RC large-R jet $D_{2}$ vs $m^{top}$ in 187.5 GeV < $m^{top}$ < 222.5 GeV and the Absolute double-differential cross-section as function of RC large-R jet $D_{2}$ vs $m^{top}$ in 187.5 GeV < $m^{top}$ < 222.5 GeV at particle level in the $\ell$+jets channel, accounting for the statistical uncertainty.

Relative double-differential cross-section as a function of RC large-R jet $D_{2}$ vs $m^{top}$ at particle level in the $\ell$+jets channel in 122.5 GeV < $m^{top}$ < 157.5 GeV. The measured differential cross-section is compared with the prediction obtained with the Powheg+Pythia8 Monte Carlo generator.

Relative double-differential cross-section as a function of RC large-R jet $D_{2}$ vs $m^{top}$ at particle level in the $\ell$+jets channel in 157.5 GeV < $m^{top}$ < 187.5 GeV. The measured differential cross-section is compared with the prediction obtained with the Powheg+Pythia8 Monte Carlo generator.

Relative double-differential cross-section as a function of RC large-R jet $D_{2}$ vs $m^{top}$ at particle level in the $\ell$+jets channel in 187.5 GeV < $m^{top}$ < 222.5 GeV. The measured differential cross-section is compared with the prediction obtained with the Powheg+Pythia8 Monte Carlo generator.

Absolute double-differential cross-section as a function of RC large-R jet $\tau_{32}$ vs $m^{top}$ at particle level in the $\ell$+jets channel in 122.5 GeV < $m^{top}$ < 157.5 GeV. The measured differential cross-section is compared with the prediction obtained with the Powheg+Pythia8 Monte Carlo generator.

Absolute double-differential cross-section as a function of RC large-R jet $\tau_{32}$ vs $m^{top}$ at particle level in the $\ell$+jets channel in 157.5 GeV < $m^{top}$ < 187.5 GeV. The measured differential cross-section is compared with the prediction obtained with the Powheg+Pythia8 Monte Carlo generator.

Absolute double-differential cross-section as a function of RC large-R jet $\tau_{32}$ vs $m^{top}$ at particle level in the $\ell$+jets channel in 187.5 GeV < $m^{top}$ < 222.5 GeV. The measured differential cross-section is compared with the prediction obtained with the Powheg+Pythia8 Monte Carlo generator.

Covariance matrix between the Absolute double-differential cross-section as function of RC large-R jet $\tau_{32}$ vs $m^{top}$ in 122.5 GeV < $m^{top}$ < 157.5 GeV and the Absolute double-differential cross-section as function of RC large-R jet $\tau_{32}$ vs $m^{top}$ in 122.5 GeV < $m^{top}$ < 157.5 GeV at particle level in the $\ell$+jets channel, accounting for the statistical uncertainty.

Covariance matrix between the Absolute double-differential cross-section as function of RC large-R jet $\tau_{32}$ vs $m^{top}$ in 157.5 GeV < $m^{top}$ < 187.5 GeV and the Absolute double-differential cross-section as function of RC large-R jet $\tau_{32}$ vs $m^{top}$ in 122.5 GeV < $m^{top}$ < 157.5 GeV at particle level in the $\ell$+jets channel, accounting for the statistical uncertainty.

Covariance matrix between the Absolute double-differential cross-section as function of RC large-R jet $\tau_{32}$ vs $m^{top}$ in 157.5 GeV < $m^{top}$ < 187.5 GeV and the Absolute double-differential cross-section as function of RC large-R jet $\tau_{32}$ vs $m^{top}$ in 157.5 GeV < $m^{top}$ < 187.5 GeV at particle level in the $\ell$+jets channel, accounting for the statistical uncertainty.

Covariance matrix between the Absolute double-differential cross-section as function of RC large-R jet $\tau_{32}$ vs $m^{top}$ in 187.5 GeV < $m^{top}$ < 222.5 GeV and the Absolute double-differential cross-section as function of RC large-R jet $\tau_{32}$ vs $m^{top}$ in 122.5 GeV < $m^{top}$ < 157.5 GeV at particle level in the $\ell$+jets channel, accounting for the statistical uncertainty.

Covariance matrix between the Absolute double-differential cross-section as function of RC large-R jet $\tau_{32}$ vs $m^{top}$ in 187.5 GeV < $m^{top}$ < 222.5 GeV and the Absolute double-differential cross-section as function of RC large-R jet $\tau_{32}$ vs $m^{top}$ in 157.5 GeV < $m^{top}$ < 187.5 GeV at particle level in the $\ell$+jets channel, accounting for the statistical uncertainty.

Covariance matrix between the Absolute double-differential cross-section as function of RC large-R jet $\tau_{32}$ vs $m^{top}$ in 187.5 GeV < $m^{top}$ < 222.5 GeV and the Absolute double-differential cross-section as function of RC large-R jet $\tau_{32}$ vs $m^{top}$ in 187.5 GeV < $m^{top}$ < 222.5 GeV at particle level in the $\ell$+jets channel, accounting for the statistical uncertainty.

Relative double-differential cross-section as a function of RC large-R jet $\tau_{32}$ vs $m^{top}$ at particle level in the $\ell$+jets channel in 122.5 GeV < $m^{top}$ < 157.5 GeV. The measured differential cross-section is compared with the prediction obtained with the Powheg+Pythia8 Monte Carlo generator.

Relative double-differential cross-section as a function of RC large-R jet $\tau_{32}$ vs $m^{top}$ at particle level in the $\ell$+jets channel in 157.5 GeV < $m^{top}$ < 187.5 GeV. The measured differential cross-section is compared with the prediction obtained with the Powheg+Pythia8 Monte Carlo generator.

Relative double-differential cross-section as a function of RC large-R jet $\tau_{32}$ vs $m^{top}$ at particle level in the $\ell$+jets channel in 187.5 GeV < $m^{top}$ < 222.5 GeV. The measured differential cross-section is compared with the prediction obtained with the Powheg+Pythia8 Monte Carlo generator.

Absolute double-differential cross-section as a function of RC large-R jet $\tau_{32}$ vs $p_{T}$ at particle level in the $\ell$+jets channel in 350.0 GeV < $p_{T}$ < 500.0 GeV. The measured differential cross-section is compared with the prediction obtained with the Powheg+Pythia8 Monte Carlo generator.

Absolute double-differential cross-section as a function of RC large-R jet $\tau_{32}$ vs $p_{T}$ at particle level in the $\ell$+jets channel in 500.0 GeV < $p_{T}$ < 550.0 GeV. The measured differential cross-section is compared with the prediction obtained with the Powheg+Pythia8 Monte Carlo generator.

Absolute double-differential cross-section as a function of RC large-R jet $\tau_{32}$ vs $p_{T}$ at particle level in the $\ell$+jets channel in 550.0 GeV < $p_{T}$ < 650.0 GeV. The measured differential cross-section is compared with the prediction obtained with the Powheg+Pythia8 Monte Carlo generator.

Absolute double-differential cross-section as a function of RC large-R jet $\tau_{32}$ vs $p_{T}$ at particle level in the $\ell$+jets channel in 650.0 GeV < $p_{T}$ < 750.0 GeV. The measured differential cross-section is compared with the prediction obtained with the Powheg+Pythia8 Monte Carlo generator.

Absolute double-differential cross-section as a function of RC large-R jet $\tau_{32}$ vs $p_{T}$ at particle level in the $\ell$+jets channel in 750.0 GeV < $p_{T}$ < 2000.0 GeV. The measured differential cross-section is compared with the prediction obtained with the Powheg+Pythia8 Monte Carlo generator.

Covariance matrix between the Absolute double-differential cross-section as function of RC large-R jet $\tau_{32}$ vs $p_{T}$ in 350.0 GeV < $p_{T}$ < 500.0 GeV and the Absolute double-differential cross-section as function of RC large-R jet $\tau_{32}$ vs $p_{T}$ in 350.0 GeV < $p_{T}$ < 500.0 GeV at particle level in the $\ell$+jets channel, accounting for the statistical uncertainty.

Covariance matrix between the Absolute double-differential cross-section as function of RC large-R jet $\tau_{32}$ vs $p_{T}$ in 500.0 GeV < $p_{T}$ < 550.0 GeV and the Absolute double-differential cross-section as function of RC large-R jet $\tau_{32}$ vs $p_{T}$ in 350.0 GeV < $p_{T}$ < 500.0 GeV at particle level in the $\ell$+jets channel, accounting for the statistical uncertainty.

Covariance matrix between the Absolute double-differential cross-section as function of RC large-R jet $\tau_{32}$ vs $p_{T}$ in 500.0 GeV < $p_{T}$ < 550.0 GeV and the Absolute double-differential cross-section as function of RC large-R jet $\tau_{32}$ vs $p_{T}$ in 500.0 GeV < $p_{T}$ < 550.0 GeV at particle level in the $\ell$+jets channel, accounting for the statistical uncertainty.

Covariance matrix between the Absolute double-differential cross-section as function of RC large-R jet $\tau_{32}$ vs $p_{T}$ in 550.0 GeV < $p_{T}$ < 650.0 GeV and the Absolute double-differential cross-section as function of RC large-R jet $\tau_{32}$ vs $p_{T}$ in 350.0 GeV < $p_{T}$ < 500.0 GeV at particle level in the $\ell$+jets channel, accounting for the statistical uncertainty.

Covariance matrix between the Absolute double-differential cross-section as function of RC large-R jet $\tau_{32}$ vs $p_{T}$ in 550.0 GeV < $p_{T}$ < 650.0 GeV and the Absolute double-differential cross-section as function of RC large-R jet $\tau_{32}$ vs $p_{T}$ in 500.0 GeV < $p_{T}$ < 550.0 GeV at particle level in the $\ell$+jets channel, accounting for the statistical uncertainty.

Covariance matrix between the Absolute double-differential cross-section as function of RC large-R jet $\tau_{32}$ vs $p_{T}$ in 550.0 GeV < $p_{T}$ < 650.0 GeV and the Absolute double-differential cross-section as function of RC large-R jet $\tau_{32}$ vs $p_{T}$ in 550.0 GeV < $p_{T}$ < 650.0 GeV at particle level in the $\ell$+jets channel, accounting for the statistical uncertainty.

Covariance matrix between the Absolute double-differential cross-section as function of RC large-R jet $\tau_{32}$ vs $p_{T}$ in 650.0 GeV < $p_{T}$ < 750.0 GeV and the Absolute double-differential cross-section as function of RC large-R jet $\tau_{32}$ vs $p_{T}$ in 350.0 GeV < $p_{T}$ < 500.0 GeV at particle level in the $\ell$+jets channel, accounting for the statistical uncertainty.

Covariance matrix between the Absolute double-differential cross-section as function of RC large-R jet $\tau_{32}$ vs $p_{T}$ in 650.0 GeV < $p_{T}$ < 750.0 GeV and the Absolute double-differential cross-section as function of RC large-R jet $\tau_{32}$ vs $p_{T}$ in 500.0 GeV < $p_{T}$ < 550.0 GeV at particle level in the $\ell$+jets channel, accounting for the statistical uncertainty.

Covariance matrix between the Absolute double-differential cross-section as function of RC large-R jet $\tau_{32}$ vs $p_{T}$ in 650.0 GeV < $p_{T}$ < 750.0 GeV and the Absolute double-differential cross-section as function of RC large-R jet $\tau_{32}$ vs $p_{T}$ in 550.0 GeV < $p_{T}$ < 650.0 GeV at particle level in the $\ell$+jets channel, accounting for the statistical uncertainty.

Covariance matrix between the Absolute double-differential cross-section as function of RC large-R jet $\tau_{32}$ vs $p_{T}$ in 650.0 GeV < $p_{T}$ < 750.0 GeV and the Absolute double-differential cross-section as function of RC large-R jet $\tau_{32}$ vs $p_{T}$ in 650.0 GeV < $p_{T}$ < 750.0 GeV at particle level in the $\ell$+jets channel, accounting for the statistical uncertainty.

Covariance matrix between the Absolute double-differential cross-section as function of RC large-R jet $\tau_{32}$ vs $p_{T}$ in 750.0 GeV < $p_{T}$ < 2000.0 GeV and the Absolute double-differential cross-section as function of RC large-R jet $\tau_{32}$ vs $p_{T}$ in 350.0 GeV < $p_{T}$ < 500.0 GeV at particle level in the $\ell$+jets channel, accounting for the statistical uncertainty.

Covariance matrix between the Absolute double-differential cross-section as function of RC large-R jet $\tau_{32}$ vs $p_{T}$ in 750.0 GeV < $p_{T}$ < 2000.0 GeV and the Absolute double-differential cross-section as function of RC large-R jet $\tau_{32}$ vs $p_{T}$ in 500.0 GeV < $p_{T}$ < 550.0 GeV at particle level in the $\ell$+jets channel, accounting for the statistical uncertainty.

Covariance matrix between the Absolute double-differential cross-section as function of RC large-R jet $\tau_{32}$ vs $p_{T}$ in 750.0 GeV < $p_{T}$ < 2000.0 GeV and the Absolute double-differential cross-section as function of RC large-R jet $\tau_{32}$ vs $p_{T}$ in 550.0 GeV < $p_{T}$ < 650.0 GeV at particle level in the $\ell$+jets channel, accounting for the statistical uncertainty.

Covariance matrix between the Absolute double-differential cross-section as function of RC large-R jet $\tau_{32}$ vs $p_{T}$ in 750.0 GeV < $p_{T}$ < 2000.0 GeV and the Absolute double-differential cross-section as function of RC large-R jet $\tau_{32}$ vs $p_{T}$ in 650.0 GeV < $p_{T}$ < 750.0 GeV at particle level in the $\ell$+jets channel, accounting for the statistical uncertainty.

Covariance matrix between the Absolute double-differential cross-section as function of RC large-R jet $\tau_{32}$ vs $p_{T}$ in 750.0 GeV < $p_{T}$ < 2000.0 GeV and the Absolute double-differential cross-section as function of RC large-R jet $\tau_{32}$ vs $p_{T}$ in 750.0 GeV < $p_{T}$ < 2000.0 GeV at particle level in the $\ell$+jets channel, accounting for the statistical uncertainty.

Relative double-differential cross-section as a function of RC large-R jet $\tau_{32}$ vs $p_{T}$ at particle level in the $\ell$+jets channel in 350.0 GeV < $p_{T}$ < 500.0 GeV. The measured differential cross-section is compared with the prediction obtained with the Powheg+Pythia8 Monte Carlo generator.

Relative double-differential cross-section as a function of RC large-R jet $\tau_{32}$ vs $p_{T}$ at particle level in the $\ell$+jets channel in 500.0 GeV < $p_{T}$ < 550.0 GeV. The measured differential cross-section is compared with the prediction obtained with the Powheg+Pythia8 Monte Carlo generator.

Relative double-differential cross-section as a function of RC large-R jet $\tau_{32}$ vs $p_{T}$ at particle level in the $\ell$+jets channel in 550.0 GeV < $p_{T}$ < 650.0 GeV. The measured differential cross-section is compared with the prediction obtained with the Powheg+Pythia8 Monte Carlo generator.

Relative double-differential cross-section as a function of RC large-R jet $\tau_{32}$ vs $p_{T}$ at particle level in the $\ell$+jets channel in 650.0 GeV < $p_{T}$ < 750.0 GeV. The measured differential cross-section is compared with the prediction obtained with the Powheg+Pythia8 Monte Carlo generator.

Relative double-differential cross-section as a function of RC large-R jet $\tau_{32}$ vs $p_{T}$ at particle level in the $\ell$+jets channel in 750.0 GeV < $p_{T}$ < 2000.0 GeV. The measured differential cross-section is compared with the prediction obtained with the Powheg+Pythia8 Monte Carlo generator.

Absolute double-differential cross-section as a function of RC large-R jet $D_{2}$ vs $p_{T}$ at particle level in the $\ell$+jets channel in 350.0 GeV < $p_{T}$ < 500.0 GeV. The measured differential cross-section is compared with the prediction obtained with the Powheg+Pythia8 Monte Carlo generator.

Absolute double-differential cross-section as a function of RC large-R jet $D_{2}$ vs $p_{T}$ at particle level in the $\ell$+jets channel in 500.0 GeV < $p_{T}$ < 550.0 GeV. The measured differential cross-section is compared with the prediction obtained with the Powheg+Pythia8 Monte Carlo generator.

Absolute double-differential cross-section as a function of RC large-R jet $D_{2}$ vs $p_{T}$ at particle level in the $\ell$+jets channel in 550.0 GeV < $p_{T}$ < 650.0 GeV. The measured differential cross-section is compared with the prediction obtained with the Powheg+Pythia8 Monte Carlo generator.

Absolute double-differential cross-section as a function of RC large-R jet $D_{2}$ vs $p_{T}$ at particle level in the $\ell$+jets channel in 650.0 GeV < $p_{T}$ < 750.0 GeV. The measured differential cross-section is compared with the prediction obtained with the Powheg+Pythia8 Monte Carlo generator.

Absolute double-differential cross-section as a function of RC large-R jet $D_{2}$ vs $p_{T}$ at particle level in the $\ell$+jets channel in 750.0 GeV < $p_{T}$ < 2000.0 GeV. The measured differential cross-section is compared with the prediction obtained with the Powheg+Pythia8 Monte Carlo generator.

Covariance matrix between the Absolute double-differential cross-section as function of RC large-R jet $D_{2}$ vs $p_{T}$ in 350.0 GeV < $p_{T}$ < 500.0 GeV and the Absolute double-differential cross-section as function of RC large-R jet $D_{2}$ vs $p_{T}$ in 350.0 GeV < $p_{T}$ < 500.0 GeV at particle level in the $\ell$+jets channel, accounting for the statistical uncertainty.

Covariance matrix between the Absolute double-differential cross-section as function of RC large-R jet $D_{2}$ vs $p_{T}$ in 500.0 GeV < $p_{T}$ < 550.0 GeV and the Absolute double-differential cross-section as function of RC large-R jet $D_{2}$ vs $p_{T}$ in 350.0 GeV < $p_{T}$ < 500.0 GeV at particle level in the $\ell$+jets channel, accounting for the statistical uncertainty.

Covariance matrix between the Absolute double-differential cross-section as function of RC large-R jet $D_{2}$ vs $p_{T}$ in 500.0 GeV < $p_{T}$ < 550.0 GeV and the Absolute double-differential cross-section as function of RC large-R jet $D_{2}$ vs $p_{T}$ in 500.0 GeV < $p_{T}$ < 550.0 GeV at particle level in the $\ell$+jets channel, accounting for the statistical uncertainty.

Covariance matrix between the Absolute double-differential cross-section as function of RC large-R jet $D_{2}$ vs $p_{T}$ in 550.0 GeV < $p_{T}$ < 650.0 GeV and the Absolute double-differential cross-section as function of RC large-R jet $D_{2}$ vs $p_{T}$ in 350.0 GeV < $p_{T}$ < 500.0 GeV at particle level in the $\ell$+jets channel, accounting for the statistical uncertainty.

Covariance matrix between the Absolute double-differential cross-section as function of RC large-R jet $D_{2}$ vs $p_{T}$ in 550.0 GeV < $p_{T}$ < 650.0 GeV and the Absolute double-differential cross-section as function of RC large-R jet $D_{2}$ vs $p_{T}$ in 500.0 GeV < $p_{T}$ < 550.0 GeV at particle level in the $\ell$+jets channel, accounting for the statistical uncertainty.

Covariance matrix between the Absolute double-differential cross-section as function of RC large-R jet $D_{2}$ vs $p_{T}$ in 550.0 GeV < $p_{T}$ < 650.0 GeV and the Absolute double-differential cross-section as function of RC large-R jet $D_{2}$ vs $p_{T}$ in 550.0 GeV < $p_{T}$ < 650.0 GeV at particle level in the $\ell$+jets channel, accounting for the statistical uncertainty.

Covariance matrix between the Absolute double-differential cross-section as function of RC large-R jet $D_{2}$ vs $p_{T}$ in 650.0 GeV < $p_{T}$ < 750.0 GeV and the Absolute double-differential cross-section as function of RC large-R jet $D_{2}$ vs $p_{T}$ in 350.0 GeV < $p_{T}$ < 500.0 GeV at particle level in the $\ell$+jets channel, accounting for the statistical uncertainty.

Covariance matrix between the Absolute double-differential cross-section as function of RC large-R jet $D_{2}$ vs $p_{T}$ in 650.0 GeV < $p_{T}$ < 750.0 GeV and the Absolute double-differential cross-section as function of RC large-R jet $D_{2}$ vs $p_{T}$ in 500.0 GeV < $p_{T}$ < 550.0 GeV at particle level in the $\ell$+jets channel, accounting for the statistical uncertainty.

Covariance matrix between the Absolute double-differential cross-section as function of RC large-R jet $D_{2}$ vs $p_{T}$ in 650.0 GeV < $p_{T}$ < 750.0 GeV and the Absolute double-differential cross-section as function of RC large-R jet $D_{2}$ vs $p_{T}$ in 550.0 GeV < $p_{T}$ < 650.0 GeV at particle level in the $\ell$+jets channel, accounting for the statistical uncertainty.

Covariance matrix between the Absolute double-differential cross-section as function of RC large-R jet $D_{2}$ vs $p_{T}$ in 650.0 GeV < $p_{T}$ < 750.0 GeV and the Absolute double-differential cross-section as function of RC large-R jet $D_{2}$ vs $p_{T}$ in 650.0 GeV < $p_{T}$ < 750.0 GeV at particle level in the $\ell$+jets channel, accounting for the statistical uncertainty.

Covariance matrix between the Absolute double-differential cross-section as function of RC large-R jet $D_{2}$ vs $p_{T}$ in 750.0 GeV < $p_{T}$ < 2000.0 GeV and the Absolute double-differential cross-section as function of RC large-R jet $D_{2}$ vs $p_{T}$ in 350.0 GeV < $p_{T}$ < 500.0 GeV at particle level in the $\ell$+jets channel, accounting for the statistical uncertainty.

Covariance matrix between the Absolute double-differential cross-section as function of RC large-R jet $D_{2}$ vs $p_{T}$ in 750.0 GeV < $p_{T}$ < 2000.0 GeV and the Absolute double-differential cross-section as function of RC large-R jet $D_{2}$ vs $p_{T}$ in 500.0 GeV < $p_{T}$ < 550.0 GeV at particle level in the $\ell$+jets channel, accounting for the statistical uncertainty.

Covariance matrix between the Absolute double-differential cross-section as function of RC large-R jet $D_{2}$ vs $p_{T}$ in 750.0 GeV < $p_{T}$ < 2000.0 GeV and the Absolute double-differential cross-section as function of RC large-R jet $D_{2}$ vs $p_{T}$ in 550.0 GeV < $p_{T}$ < 650.0 GeV at particle level in the $\ell$+jets channel, accounting for the statistical uncertainty.

Covariance matrix between the Absolute double-differential cross-section as function of RC large-R jet $D_{2}$ vs $p_{T}$ in 750.0 GeV < $p_{T}$ < 2000.0 GeV and the Absolute double-differential cross-section as function of RC large-R jet $D_{2}$ vs $p_{T}$ in 650.0 GeV < $p_{T}$ < 750.0 GeV at particle level in the $\ell$+jets channel, accounting for the statistical uncertainty.

Covariance matrix between the Absolute double-differential cross-section as function of RC large-R jet $D_{2}$ vs $p_{T}$ in 750.0 GeV < $p_{T}$ < 2000.0 GeV and the Absolute double-differential cross-section as function of RC large-R jet $D_{2}$ vs $p_{T}$ in 750.0 GeV < $p_{T}$ < 2000.0 GeV at particle level in the $\ell$+jets channel, accounting for the statistical uncertainty.

Relative double-differential cross-section as a function of RC large-R jet $D_{2}$ vs $p_{T}$ at particle level in the $\ell$+jets channel in 350.0 GeV < $p_{T}$ < 500.0 GeV. The measured differential cross-section is compared with the prediction obtained with the Powheg+Pythia8 Monte Carlo generator.

Relative double-differential cross-section as a function of RC large-R jet $D_{2}$ vs $p_{T}$ at particle level in the $\ell$+jets channel in 500.0 GeV < $p_{T}$ < 550.0 GeV. The measured differential cross-section is compared with the prediction obtained with the Powheg+Pythia8 Monte Carlo generator.

Relative double-differential cross-section as a function of RC large-R jet $D_{2}$ vs $p_{T}$ at particle level in the $\ell$+jets channel in 550.0 GeV < $p_{T}$ < 650.0 GeV. The measured differential cross-section is compared with the prediction obtained with the Powheg+Pythia8 Monte Carlo generator.

Relative double-differential cross-section as a function of RC large-R jet $D_{2}$ vs $p_{T}$ at particle level in the $\ell$+jets channel in 650.0 GeV < $p_{T}$ < 750.0 GeV. The measured differential cross-section is compared with the prediction obtained with the Powheg+Pythia8 Monte Carlo generator.

Relative double-differential cross-section as a function of RC large-R jet $D_{2}$ vs $p_{T}$ at particle level in the $\ell$+jets channel in 750.0 GeV < $p_{T}$ < 2000.0 GeV. The measured differential cross-section is compared with the prediction obtained with the Powheg+Pythia8 Monte Carlo generator.

Upper limits on $\mathcal{B}(H\rightarrow aa\rightarrow 4\gamma)$ at 95% CL as a function of the axion mass and for ALP-photon coupling $C_{a\gamma\gamma}=0.01$.

Upper limits on $\mathcal{B}(H\rightarrow aa\rightarrow 4\gamma)$ at 95% CL as a function of the axion mass and for ALP-photon coupling $C_{a\gamma\gamma}=5\times10^{-4}$

Upper limits on $\mathcal{B}(H\rightarrow aa\rightarrow 4\gamma)$ at 95% CL as a function of the axion mass and for ALP-photon coupling $C_{a\gamma\gamma}=10^{-5}$

Upper limits on $B(H\rightarrow 4\gamma)$ at 95% CL as a function of the signal mass hypothesis and for the assumption of promptly decaying ALPs.

Identification efficiency for isolated photons that pass the tight photon ID requirements in dependence of the opening angle in $\Delta R$ of both decay photons of the axion and the decay radius of the originating axion from the primary vertex for axions with masses between 100 MeV and 300 MeV. All four axion decay photons on MC truth level are required to fulfil $|\eta|<2.5$ and $p_T>5$ GeV. Only a mild dependency on the decay radius can be seen, since the distance between the two photons is small and hence the combined decay signature in the calorimeter still points to the primary vertex.

Efficiency for classified merged isolated photons in dependence of the opening angle in $\Delta R$ of both decay photons of the axion and the decay radius of the originating axion from the primary vertex for axions with masses between 100 MeV and 300 MeV. All four axion decay photons on MC truth level are required to fulfil $|\eta|<2.5$ and $p_T>5$ GeV. Only a mild dependency on the decay radius can be seen, since the distance between the two photons is small and hence the combined decay signature in the calorimeter still points to the primary vertex.

Limits on the ALP mass and coupling to photons at 95\%~CL, assuming $\mathcal{B}(a\rightarrow\gamma\gamma) = 1$, $\Lambda = 1\,TeV$ with $|C_{aH}^{\text{eff}}|$ = $1$ (solid line) and $|C_{aH}^{\text{eff}}|$ = $0.1$.

The number of data and estimated background events in the signal region of the most sensitive categories..

Observed 95% CL upper limits on B(H&rarr;&nbsp;2&gamma;_d+X) for different &gamma;_d masses and a 125&nbsp;GeV Higgs boson, as a function of the dark-photon mean proper decay length c&tau;. The limits are shown for all search channels combined, assuming an FRVZ signal model. The hatched band denotes the region in which the branching ratio is larger than unity.

Observed 95% CL upper limits on B(H&rarr;&nbsp;2&gamma;_d+X) for a &gamma;_d mass of 100&nbsp;MeV, as a function of the dark-photon mean proper decay length c&tau;. The solid black curve shows the observed exclusion limit from their statistical combination. The hatched band denotes the region in which the branching ratio is larger than unity.

Observed 95% CL upper limits on B(H&rarr;&nbsp;2&gamma;_d+X) for a &gamma;_d mass of 10&nbsp;GeV, as a function of the dark-photon mean proper decay length c&tau;. The solid black curve shows the observed exclusion limit from their statistical combination. The hatched band denotes the region in which the branching ratio is larger than unity.

Acceptance times efficiency extrapolated as a function of c&tau;_&gamma;d for the &mu;DPJ channel (SR_&mu;), obtained for the different FRVZ MC signal mass points.

Acceptance times efficiency extrapolated as a function of c&tau;_&gamma;d for the low E_T^miss caloDPJ channel (SR_c^L), obtained for the different FRVZ MC signal mass points.

Acceptance times efficiency extrapolated as a function of c&tau;_&gamma;d for the high E_T^miss caloDPJ channel (SR_c^H), obtained for the different FRVZ MC signal mass points.

Observed (solid line) and expected (dashed line) 95% CL upper limits on B(H&rarr;&nbsp;2&gamma;_d+X) for FRVZ processes with a 125&nbsp;GeV Higgs boson as a function of dark photon mean proper decay length for the ggF, WH and VBF analyses. The full combination is shown as solid or dashed black lines for a dark photon mass of 17&nbsp;MeV.

Observed (solid line) and expected (dashed line) 95% CL upper limits on B(H&rarr;&nbsp;2&gamma;_d+X) for FRVZ processes with a 125&nbsp;GeV Higgs boson as a function of dark photon mean proper decay length for the ggF, WH and VBF analyses. The full combination is shown as solid or dashed black lines for a dark photon mass of 50&nbsp;MeV.

Observed (solid line) and expected (dashed line) 95% CL upper limits on B(H&rarr;&nbsp;2&gamma;_d+X) for FRVZ processes with a 125&nbsp;GeV Higgs boson as a function of dark photon mean proper decay length for the ggF, WH and VBF analyses. The full combination is shown as solid or dashed black lines for a dark photon mass of 400&nbsp;MeV.

Observed (solid line) and expected (dashed line) 95% CL upper limits on B(H&rarr;&nbsp;2&gamma;_d+X) for FRVZ processes with a 125&nbsp;GeV Higgs boson as a function of dark photon mean proper decay length for the ggF, WH and VBF analyses. The full combination is shown as solid or dashed black lines for a dark photon mass of 900&nbsp;MeV.

Observed (solid line) and expected (dashed line) 95% CL upper limits on B(H&rarr;&nbsp;2&gamma;_d+X) for FRVZ processes with a 125&nbsp;GeV Higgs boson as a function of dark photon mean proper decay length for the ggF, WH and VBF analyses. The full combination is shown as solid or dashed black lines for a dark photon mass of 2&nbsp;GeV.

Observed (solid line) and expected (dashed line) 95% CL upper limits on B(H&rarr;&nbsp;2&gamma;_d+X) for FRVZ processes with a 125&nbsp;GeV Higgs boson as a function of dark photon mean proper decay length for the ggF, WH and VBF analyses. The full combination is shown as solid or dashed black lines for a dark photon mass of 6&nbsp;GeV.

Observed (solid line) and expected (dashed line) 95% CL upper limits on B(H&rarr;&nbsp;2&gamma;_d+X) for FRVZ processes with a 125&nbsp;GeV Higgs boson as a function of dark photon mean proper decay length for the ggF, WH and VBF analyses. The full combination is shown as solid or dashed black lines for a dark photon mass of 15&nbsp;GeV.

Observed 95% CL upper limits on B(H&rarr;&nbsp;2&gamma;_d+X) for different &gamma;_d masses and a 125&nbsp;GeV Higgs boson, as a function of the dark-photon mean proper decay length c&tau;. The limits are shown for the SR_&mu; search channel, assuming an FRVZ signal model. The hatched band denotes the region in which the branching ratio is larger than unity.

Observed 95% CL upper limits on B(H&rarr;&nbsp;2&gamma;_d+X) for different &gamma;_d masses and a 125&nbsp;GeV Higgs boson, as a function of the dark-photon mean proper decay length c&tau;. The limits are shown for the SR_c^L search channel, assuming an FRVZ signal model. The hatched band denotes the region in which the branching ratio is larger than unity.

Observed 95% CL upper limits on B(H&rarr;&nbsp;2&gamma;_d+X) for different &gamma;_d masses and a 125&nbsp;GeV Higgs boson, as a function of the dark-photon mean proper decay length c&tau;. The limits are shown for the SR_c^H search channel, assuming an FRVZ signal model. The hatched band denotes the region in which the branching ratio is larger than unity.

Observed 95% CL upper limits on B(H&rarr;&nbsp;2&gamma;_d+X) for different &gamma;_d masses and a 125&nbsp;GeV Higgs boson, as a function of the dark-photon mean proper decay length c&tau;. The limits are shown for all search channels combined, assuming an FRVZ signal model. The hatched band denotes the region in which the branching ratio is larger than unity.

Unfolded $M(l^{-}\gamma)$ distribution for $Z \to \mu\mu\gamma$ process with dressed leptons and bkg subtraction. $M_{ll}>20$ GeV. Nexp.un f. = 64822.9 $\pm$ 254.6, NPowHeg truth =634276.

Unfolded dR distribution for $Z \to ee\gamma$ process with dressed leptons and bkg subtraction. $M_{ll}>20$ GeV. Nexp.un f. = 64273.2 $\pm$ 253.5, NPowHeg truth =338671.

Unfolded dR distribution for $Z \to \mu\mu\gamma$ process with dressed leptons and bkg subtraction. $M_{ll}>20$ GeV. Nexp.un f. = 65362.4 $\pm$ 255.7 , NPowHeg truth =634214.

Unfolded $P_{T}^{\gamma}$ distribution for $Z \to ee\gamma$ process with dressed leptons and bkg subtraction. $M_{ll}>20$ GeV. Nexp.un f. = 63486.8 $\pm$ 252.0, NPowHeg truth =333348.

Unfolded $P_{T}^{\gamma}$ distribution for $Z \to \mu\mu\gamma$ process with dressed leptons and bkg subtraction. $M_{ll}>20$ GeV. Nexp.un f. = 64432.6 $\pm$ 253.8, NPowHeg truth =624059.

Unfolded $M(l^{+}\gamma)$ distribution for $Z \to ee\gamma$ process with dressed leptons and bkg subtraction. $M_{ll}>45$ GeV. Nexp.un f. = 57095.9 $\pm$ 254 , NPowHeg truth =301666, N Sherpa truth =76327.

Unfolded $M(l^{-}\gamma)$ distribution for $Z \to ee\gamma$ process with dressed leptons and bkg subtraction. $M_{ll}>45$ GeV. Nexp.un f. = 57130 $\pm$ 254 , NPowHeg truth =301666, N Sherpa truth =763273.

Unfolded $M(l^{+}\gamma)$ distribution for $Z \to \mu\mu\gamma$ process with dressed leptons and bkg subtraction. $M_{ll}>45$ GeV. Nexp.un f. = 58310.9 $\pm$ 245, NPowHeg truth =564672, N Sherpa truth =814016.

Unfolded $M(l^{-}\gamma)$ distribution for $Z \to \mu\mu\gamma$ process with dressed leptons and bkg subtraction. $M_{ll}>45$ GeV. Nexp.un f. = 58320.9 $\pm$ 245, NPowHeg truth =564672, N Sherpa truth =814016.

Unfolded dR distribution for $Z \to ee\gamma$ process with dressed leptons and bkg subtraction. $M_{ll}>45$ GeV. Nexp.un f. = 57705.4 $\pm$ 254 , NPowHeg truth =301655, N Sherpa truth =763259.

Unfolded dR distribution for $Z \to \mu\mu\gamma$ process with dressed leptons and bkg subtraction. $M_{ll}>45$ GeV. Nexp.un f. = 58869.2 $\pm$ 256, NPowHeg truth =564660, NS herpa truth =813995.

Unfolded $P_{T}^{\gamma}$ distribution for $Z \to ee\gamma$ process with dressed leptons and bkg subtraction. $M_{ll}>45$ GeV. Nexp.un f. = 57358.6 $\pm$ 254 , NPowHeg truth =298662, N Sherpa truth =756609.

Unfolded $P_{T}^{\gamma}$ distribution for $Z \to \mu\mu\gamma$ process with dressed leptons and bkg subtraction. $M_{ll}>45$ GeV. Nexp.un f. = 58226.4 $\pm$ 254, NPowHeg truth =558949, N Sherpa truth =806988.

Unfolded $M(l^{+}\gamma)$ distribution for $Z \to ee\gamma$ process with bare leptons and bkg subtraction. $M_{ll}>20$ GeV. Nexp.un f. = 63717.4 $\pm$ 252.4, NPowHeg truth =338714.

Unfolded $M(l^{-}\gamma)$ distribution for $Z \to ee\gamma$ process with bare leptons and bkg subtraction. $M_{ll}>20$ GeV. Nexp.un f. = 63855.8 $\pm$ 252.7 , NPowHeg truth =338708.

Unfolded $M(l^{+}\gamma)$ distribution for $Z \to \mu\mu\gamma$ process with bare leptons and bkg subtraction. $M_{ll}>20$ GeV. Nexp.un f. = 64809.8 $\pm$ 254.6, NPowHeg truth =634285.

Unfolded $M(l^{-}\gamma)$ distribution for $Z \to \mu\mu\gamma$ process with bare leptons and bkg subtraction. $M_{ll}>20$ GeV. Nexp.un f. = 64822.9 $\pm$ 254.6, NPowHeg truth =634276.

Unfolded dR distribution for $Z \to ee\gamma$ process with bare leptons and bkg subtraction. $M_{ll}>20$ GeV. Nexp.un f. = 64273.2 $\pm$ 253.5, NPowHeg truth =338671.

Unfolded dR distribution for $Z \to \mu\mu\gamma$ process with bare leptons and bkg subtraction. $M_{ll}>20$ GeV. Nexp.un f. = 65362.4 $\pm$ 255.7 , NPowHeg truth =634214.

Unfolded $P_{T}^{\gamma}$ distribution for $Z \to ee\gamma$ process with bare leptons and bkg subtraction. $M_{ll}>20$ GeV. Nexp.un f. = 63486.8 $\pm$ 252.0, NPowHeg truth =333348.

Unfolded $P_{T}^{\gamma}$ distribution for $Z \to \mu\mu\gamma$ process with bare leptons and bkg subtraction. $M_{ll}>20$ GeV. Nexp.un f. = 64432.6 $\pm$ 253.8, NPowHeg truth =624059.

Unfolded $M(l^{+}\gamma)$ distribution for $Z \to ee\gamma$ process with bare leptons and bkg subtraction. $M_{ll}>45$ GeV. Nexp.un f. = 57095.9 $\pm$ 254 , NPowHeg truth =301666, N Sherpa truth =76327.

Unfolded $M(l^{-}\gamma)$ distribution for $Z \to ee\gamma$ process with bare leptons and bkg subtraction. $M_{ll}>45$ GeV. Nexp.un f. = 57130 $\pm$ 254 , NPowHeg truth =301666, N Sherpa truth =763273.

Unfolded $M(l^{+}\gamma)$ distribution for $Z \to \mu\mu\gamma$ process with bare leptons and bkg subtraction. $M_{ll}>45$ GeV. Nexp.un f. = 58310.9 $\pm$ 245, NPowHeg truth =564672, N Sherpa truth =814016.

Unfolded $M(l^{-}\gamma)$ distribution for $Z \to \mu\mu\gamma$ process with bare leptons and bkg subtraction. $M_{ll}>45$ GeV. Nexp.un f. = 58320.9 $\pm$ 245, NPowHeg truth =564672, N Sherpa truth =814016.

Unfolded dR distribution for $Z \to ee\gamma$ process with bare leptons and bkg subtraction. $M_{ll}>45$ GeV. Nexp.un f. = 57705.4 $\pm$ 254 , NPowHeg truth =301655, N Sherpa truth =763259.

Unfolded dR distribution for $Z \to \mu\mu\gamma$ process with bare leptons and bkg subtraction. $M_{ll}>45$ GeV. Nexp.un f. = 58869.2 $\pm$ 256, NPowHeg truth =564660, NS herpa truth =813995.

Unfolded $P_{T}^{\gamma}$ distribution for $Z \to ee\gamma$ process with bare leptons and bkg subtraction. $M_{ll}>45$ GeV. Nexp.un f. = 57358.6 $\pm$ 254 , NPowHeg truth =298662, N Sherpa truth =756609.

Unfolded $P_{T}^{\gamma}$ distribution for $Z \to \mu\mu\gamma$ process with bare leptons and bkg subtraction. $M_{ll}>45$ GeV. Nexp.un f. = 58226.4 $\pm$ 254, NPowHeg truth =558949, N Sherpa truth =806988.

Truth $P_{T}^{Z}$ distribution after reweghting to data.

Combined Covariance Matrix for Unfolded $M(l^{+}\gamma)$ distribution for $Z \to ee\gamma$ process with dressed leptons and bkg subtraction. $M_{ll}>20$ GeV. Nexp.un f. = 63717.4 $\pm$ 252.4, NPowHeg truth =338714.

Combined Covariance Matrix for Unfolded $M(l^{-}\gamma)$ distribution for $Z \to ee\gamma$ process with dressed leptons and bkg subtraction. $M_{ll}>20$ GeV. Nexp.un f. = 63855.8 $\pm$ 252.7 , NPowHeg truth =338708.

Combined Covariance Matrix for Unfolded $M(l^{+}\gamma)$ distribution for $Z \to \mu\mu\gamma$ process with dressed leptons and bkg subtraction. $M_{ll}>20$ GeV. Nexp.un f. = 64809.8 $\pm$ 254.6, NPowHeg truth =634285.

Combined Covariance Matrix for Unfolded $M(l^{-}\gamma)$ distribution for $Z \to \mu\mu\gamma$ process with dressed leptons and bkg subtraction. $M_{ll}>20$ GeV. Nexp.un f. = 64822.9 $\pm$ 254.6, NPowHeg truth =634276.

Combined Covariance Matrix for Unfolded dR distribution for $Z \to ee\gamma$ process with dressed leptons and bkg subtraction. $M_{ll}>20$ GeV. Nexp.un f. = 64273.2 $\pm$ 253.5, NPowHeg truth =338671.

Combined Covariance Matrix for Unfolded dR distribution for $Z \to \mu\mu\gamma$ process with dressed leptons and bkg subtraction. $M_{ll}>20$ GeV. Nexp.un f. = 65362.4 $\pm$ 255.7 , NPowHeg truth =634214.

Combined Covariance Matrix for Unfolded $P_{T}^{\gamma}$ distribution for $Z \to ee\gamma$ process with dressed leptons and bkg subtraction. $M_{ll}>20$ GeV. Nexp.un f. = 63486.8 $\pm$ 252.0, NPowHeg truth =333348.

Combined Covariance Matrix for Unfolded $P_{T}^{\gamma}$ distribution for $Z \to \mu\mu\gamma$ process with dressed leptons and bkg subtraction. $M_{ll}>20$ GeV. Nexp.un f. = 64432.6 $\pm$ 253.8, NPowHeg truth =624059.

Combined Covariance Matrix for Unfolded $M(l^{+}\gamma)$ distribution for $Z \to ee\gamma$ process with dressed leptons and bkg subtraction. $M_{ll}>45$ GeV. Nexp.un f. = 57095.9 $\pm$ 254 , NPowHeg truth =301666, N Sherpa truth =76327.

Combined Covariance Matrix for Unfolded $M(l^{-}\gamma)$ distribution for $Z \to ee\gamma$ process with dressed leptons and bkg subtraction. $M_{ll}>45$ GeV. Nexp.un f. = 57130 $\pm$ 254 , NPowHeg truth =301666, N Sherpa truth =763273.

Combined Covariance Matrix for Unfolded $M(l^{+}\gamma)$ distribution for $Z \to \mu\mu\gamma$ process with dressed leptons and bkg subtraction. $M_{ll}>45$ GeV. Nexp.un f. = 58310.9 $\pm$ 245, NPowHeg truth =564672, N Sherpa truth =814016.

Combined Covariance Matrix for Unfolded $M(l^{-}\gamma)$ distribution for $Z \to \mu\mu\gamma$ process with dressed leptons and bkg subtraction. $M_{ll}>45$ GeV. Nexp.un f. = 58320.9 $\pm$ 245, NPowHeg truth =564672, N Sherpa truth =814016.

Combined Covariance Matrix for Unfolded dR distribution for $Z \to ee\gamma$ process with dressed leptons and bkg subtraction. $M_{ll}>45$ GeV. Nexp.un f. = 57705.4 $\pm$ 254 , NPowHeg truth =301655, N Sherpa truth =763259.

Combined Covariance Matrix for Unfolded dR distribution for $Z \to \mu\mu\gamma$ process with dressed leptons and bkg subtraction. $M_{ll}>45$ GeV. Nexp.un f. = 58869.2 $\pm$ 256, NPowHeg truth =564660, NS herpa truth =813995.

Combined Covariance Matrix for Unfolded $P_{T}^{\gamma}$ distribution for $Z \to ee\gamma$ process with dressed leptons and bkg subtraction. $M_{ll}>45$ GeV. Nexp.un f. = 57358.6 $\pm$ 254 , NPowHeg truth =298662, N Sherpa truth =756609.

Combined Covariance Matrix for Unfolded $P_{T}^{\gamma}$ distribution for $Z \to \mu\mu\gamma$ process with dressed leptons and bkg subtraction. $M_{ll}>45$ GeV. Nexp.un f. = 58226.4 $\pm$ 254, NPowHeg truth =558949, N Sherpa truth =806988.

Combined Covariance Matrix for Unfolded $M(l^{+}\gamma)$ distribution for $Z \to ee\gamma$ process with bare leptons and bkg subtraction. $M_{ll}>20$ GeV. Nexp.un f. = 63717.4 $\pm$ 252.4, NPowHeg truth =338714.

Combined Covariance Matrix for Unfolded $M(l^{-}\gamma)$ distribution for $Z \to ee\gamma$ process with bare leptons and bkg subtraction. $M_{ll}>20$ GeV. Nexp.un f. = 63855.8 $\pm$ 252.7 , NPowHeg truth =338708.

Combined Covariance Matrix for Unfolded $M(l^{+}\gamma)$ distribution for $Z \to \mu\mu\gamma$ process with bare leptons and bkg subtraction. $M_{ll}>20$ GeV. Nexp.un f. = 64809.8 $\pm$ 254.6, NPowHeg truth =634285.

Combined Covariance Matrix for Unfolded $M(l^{-}\gamma)$ distribution for $Z \to \mu\mu\gamma$ process with bare leptons and bkg subtraction. $M_{ll}>20$ GeV. Nexp.un f. = 64822.9 $\pm$ 254.6, NPowHeg truth =634276.

Combined Covariance Matrix for Unfolded dR distribution for $Z \to ee\gamma$ process with bare leptons and bkg subtraction. $M_{ll}>20$ GeV. Nexp.un f. = 64273.2 $\pm$ 253.5, NPowHeg truth =338671.

Combined Covariance Matrix for Unfolded dR distribution for $Z \to \mu\mu\gamma$ process with bare leptons and bkg subtraction. $M_{ll}>20$ GeV. Nexp.un f. = 65362.4 $\pm$ 255.7 , NPowHeg truth =634214.

Combined Covariance Matrix for Unfolded $P_{T}^{\gamma}$ distribution for $Z \to ee\gamma$ process with bare leptons and bkg subtraction. $M_{ll}>20$ GeV. Nexp.un f. = 63486.8 $\pm$ 252.0, NPowHeg truth =333348.

Combined Covariance Matrix for Unfolded $P_{T}^{\gamma}$ distribution for $Z \to \mu\mu\gamma$ process with bare leptons and bkg subtraction. $M_{ll}>20$ GeV. Nexp.un f. = 64432.6 $\pm$ 253.8, NPowHeg truth =624059.

Combined Covariance Matrix for Unfolded $M(l^{+}\gamma)$ distribution for $Z \to ee\gamma$ process with bare leptons and bkg subtraction. $M_{ll}>45$ GeV. Nexp.un f. = 57095.9 $\pm$ 254 , NPowHeg truth =301666, N Sherpa truth =76327.

Combined Covariance Matrix for Unfolded $M(l^{-}\gamma)$ distribution for $Z \to ee\gamma$ process with bare leptons and bkg subtraction. $M_{ll}>45$ GeV. Nexp.un f. = 57130 $\pm$ 254 , NPowHeg truth =301666, N Sherpa truth =763273.

Combined Covariance Matrix for Unfolded $M(l^{+}\gamma)$ distribution for $Z \to \mu\mu\gamma$ process with bare leptons and bkg subtraction. $M_{ll}>45$ GeV. Nexp.un f. = 58310.9 $\pm$ 245, NPowHeg truth =564672, N Sherpa truth =814016.

Combined Covariance Matrix for Unfolded $M(l^{-}\gamma)$ distribution for $Z \to \mu\mu\gamma$ process with bare leptons and bkg subtraction. $M_{ll}>45$ GeV. Nexp.un f. = 58320.9 $\pm$ 245, NPowHeg truth =564672, N Sherpa truth =814016.

Combined Covariance Matrix for Unfolded dR distribution for $Z \to ee\gamma$ process with bare leptons and bkg subtraction. $M_{ll}>45$ GeV. Nexp.un f. = 57705.4 $\pm$ 254 , NPowHeg truth =301655, N Sherpa truth =763259.

Combined Covariance Matrix for Unfolded dR distribution for $Z \to \mu\mu\gamma$ process with bare leptons and bkg subtraction. $M_{ll}>45$ GeV. Nexp.un f. = 58869.2 $\pm$ 256, NPowHeg truth =564660, NS herpa truth =813995.

Combined Covariance Matrix for Unfolded $P_{T}^{\gamma}$ distribution for $Z \to ee\gamma$ process with bare leptons and bkg subtraction. $M_{ll}>45$ GeV. Nexp.un f. = 57358.6 $\pm$ 254 , NPowHeg truth =298662, N Sherpa truth =756609.

Combined Covariance Matrix for Unfolded $P_{T}^{\gamma}$ distribution for $Z \to \mu\mu\gamma$ process with bare leptons and bkg subtraction. $M_{ll}>45$ GeV. Nexp.un f. = 58226.4 $\pm$ 254, NPowHeg truth =558949, N Sherpa truth =806988.

Unfolded $M(l^{+}\gamma\gamma)$ distribution for $Z \to ll\gamma\gamma$ process with dressed leptons.

Unfolded $M(l^{+}\gamma\gamma)$ distribution for $Z \to ll\gamma\gamma$ process with bare leptons.

Unfolded dR distribution for $Z \to ll\gamma\gamma$ process with dressed leptons for leading photon.

Unfolded dR distribution for $Z \to ll\gamma\gamma$ process with bare leptons for leading photon.

Unfolded dR distribution for $Z \to ll\gamma\gamma$ process with dressed leptons for second photon.

Unfolded dR distribution for $Z \to ll\gamma\gamma$ process with bare leptons for second photon.

Unfolded dRgg distribution for $Z \to ll\gamma\gamma$ process with dressed leptons.

Unfolded dRgg distribution for $Z \to ll\gamma\gamma$ process with bare leptons.

Unfolded $P_{T}^{\gamma 1}$ distribution for $Z \to ll\gamma\gamma$ process with dressed leptons for leading photon.

Unfolded $P_{T}^{\gamma 1}$ distribution for $Z \to ll\gamma\gamma$ process with bare leptons for leading photon.

Unfolded $P_{T}^{\gamma 2}$ distribution for $Z \to ll\gamma\gamma$ process with dressed leptons for second photon.

Unfolded $P_{T}^{\gamma 2}$ distribution for $Z \to ll\gamma\gamma$ process with bare leptons for second photon.

Summary of results for cross-section of non-prompt $\psi(2S)$ decaying to a muon pair for 13 TeV data in fb/GeV. Uncertainties are statistical and systematic, respectively.

Summary of results for non-prompt fraction of $J/\psi$ decaying to a muon pair as a percentage for 13 TeV data. Uncertainties are statistical and systematic, respectively.

Summary of results for non-prompt fraction of $\psi(2S)$ decaying to a muon pair as a percentage for 13 TeV data. Uncertainties are statistical and systematic, respectively.

Summary of results for prompt $\psi(2S)$ over $J/\psi$ ratio as a percentage for 13 TeV data. Uncertainties are statistical and systematic, respectively.

Summary of results for non-prompt $\psi(2S)$ over $J/\psi$ ratio as a percentage for 13 TeV data. Uncertainties are statistical and systematic, respectively.

Correction factors for an assumed angular dependence $\propto 1+\lambda_{\theta} \cos^2\theta $ in the helicity frame, for several values of the parameter $\lambda_{\theta}$. The correction factors were found to be the same (within 1% - 2%) for $J/\psi$ and $\psi(2S)$ mesons, for prompt and non-prompt production mechanisms, and for the three rapidity intervals considered in this paper.

Measured $p_T$ cross sections in full-lepton phase space for 1.2 < |y| < 1.6.

Measured $p_T$ cross sections in full-lepton phase space for 1.6 < |y| < 2.0.

Measured $p_T$ cross sections in full-lepton phase space for 2.0 < |y| < 2.4.

Measured $p_T$ cross sections in full-lepton phase space for 2.4 < |y| < 2.8.

Measured $p_T$ cross sections in full-lepton phase space for 2.8 < |y| < 3.6.

Measured cross sections in full-lepton phase space as a function of |y|.

Normalised measured $p_T$ cross sections in full-lepton phase space for |y| < 1.6.