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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

The observed upper limits at $95\%$ CL for the gluino pair production cross section in the GMSB model, shown in the plane of $m_{\tilde{g}}$ and $c\tau_{0}$. A branching fraction of $100\%$ for the gluino decay to a gluon and a gravitino is assumed. The area below the thick black curve represents the observed exclusion region, while the dashed red lines indicate the expected limits and their $\pm 1$ standard deviation ranges. The thin black lines show the effect of the theoretical uncertainties on the signal cross section.

The observed and expected upper limits at $95\%$ CL on the gluino pair production cross section for a gluino GMSB model with $m_{\tilde{g}}=2400$ GeV. The one (two) standard deviation variation in the expected limit is shown in the inner green (outer yellow) band. The blue solid line shows the observed limit obtained by the CMS displaced jet search

The distribution (normalized to unity) of number of ECAL cells hit in the jet for jets in a background enriched data sample (satisfying $|\eta| < 1.48$, $p_\mathrm{T}tf > 0.08$, $\mathrm{HEF} > 0.2$, $t_{\mathrm{jet}} < -3\,\mathrm{ns}$ and $E^{\mathrm{CSC}}_\mathrm{ECAL}/E_{\mathrm{ECAL}} < 0.8$) and for signal jets satisfying signal region requirements (except those on $E_{\mathrm{ECAL}}$ and $N^{\mathrm{cell}}_{\mathrm{ECAL}}$).

The distribution (normalized to unity) of $\mathrm{HEF}$ for a data sample enriched in beam halo and noise jets (satisfying $|\eta| < 1.48$, $p_\mathrm{T} > 30 \,\mathrm{GeV}$, $p_\mathrm{T}tf < 0.08$, $t_{\mathrm{jet}} < -3\,\mathrm{ns}$, $E_{\mathrm{ECAL}} > 20\,\mathrm{GeV}$ and $t^{\mathrm{RMS}}_\mathrm{jet}/t_{\mathrm{jet}} < 0.4$) and for signal jets passing signal region selections (except on $\mathrm{HEF}$).

The distribution (normalized to unity) of $E_{\mathrm{HCAL}}$ for a data sample enriched in beam halo and noise jets (satisfying $|\eta| < 1.48$, $p_\mathrm{T} > 30 \,\mathrm{GeV}$, $p_\mathrm{T}tf < 0.08$, $t_{\mathrm{jet}} < -3\,\mathrm{ns}$, $E_{\mathrm{ECAL}} > 20\,\mathrm{GeV}$ and $t^{\mathrm{RMS}}_\mathrm{jet}/t_{\mathrm{jet}} < 0.4$) and for signal jets passing signal region selections (except on $E_{\mathrm{HCAL}}$).

The distribution (normalized to unity) of $t^{\mathrm{RMS}}_\mathrm{jet}/t_{\mathrm{jet}}$ for data sample enriched in jets from noise (satisfying $|\eta| < 1.48$, $p_\mathrm{T} > 30\,\mathrm{GeV}$, $p_\mathrm{T}tf > 0.08$, $\mathrm{HEF} > 0.2$, $t_{\mathrm{jet}} < -3\,\mathrm{ns}$, $E_{\mathrm{ECAL}} > 20 \,\mathrm{GeV}$ and $E^{\mathrm{CSC}}_\mathrm{ECAL}/E_{\mathrm{ECAL}} < 0.8$) and for signal jets passing signal region selections (except on $t^{\mathrm{RMS}}_\mathrm{jet}/t_{\mathrm{jet}}$ and $t^{\mathrm{RMS}}_\mathrm{jet}$)

The distribution (normalized to unity) of $t^{\mathrm{RMS}}_\mathrm{jet}$ for data sample enriched in jets from noise (satisfying $|\eta| < 1.48$, $p_\mathrm{T} > 30\,\mathrm{GeV}$, $p_\mathrm{T}tf > 0.08$, $\mathrm{HEF} > 0.2$, $t_{\mathrm{jet}} < -3\,\mathrm{ns}$, $E_{\mathrm{ECAL}} > 20 \,\mathrm{GeV}$ and $E^{\mathrm{CSC}}_\mathrm{ECAL}/E_{\mathrm{ECAL}} < 0.8$) and for signal jets passing signal region selections (except on $t^{\mathrm{RMS}}_\mathrm{jet}/t_{\mathrm{jet}}$ and $t^{\mathrm{RMS}}_\mathrm{jet}$)

The distribution (normalized to unity) of $p_\mathrm{T}tf$ for a data sample enriched in main bunch backgrounds (satisfying $p_\mathrm{T} > 30\,\mathrm{GeV}$, $|\eta| < 1.48$, $E^{\mathrm{CSC}}_\mathrm{ECAL}/E_{\mathrm{ECAL}} < 0.8$, $|t_{\mathrm{jet}}| < 3\,\mathrm{ns}$ and $E_{\mathrm{ECAL}} >20\,\mathrm{GeV}$) and for signal jets passing signal selections (except on $p_\mathrm{T}tf$).

The distribution (normalized to unity) of $E^{\mathrm{CSC}}_\mathrm{ECAL}/E_{\mathrm{ECAL}}$ for a data sample enriched in beam halo (satisfying $p_\mathrm{T} > 30\,\mathrm{GeV}$, $|\eta| < 1.48$, $p_\mathrm{T}tf < 0.08$, $\mathrm{HEF} < 0.2$, $t^{\mathrm{RMS}}_\mathrm{jet}/t_{\mathrm{jet}} < 0.4$, $t_{\mathrm{jet}} < -3\,\mathrm{ns}$ and $E_{\mathrm{ECAL}} > 20\,\mathrm{GeV}$) and for signal jets passing signal selections (except on $E^{\mathrm{CSC}}_\mathrm{ECAL}/E_{\mathrm{ECAL}}$)

The distribution (normalized to unity) of $\mathrm{max}(\Delta \phi_{\mathrm{DT}})$ for a data sample enriched in cosmic muons (satisfying $E^{\mathrm{CSC}}_\mathrm{ECAL}/E_{\mathrm{ECAL}} < 0.8$, $p_\mathrm{T} > 30\,\mathrm{GeV}$, $|\eta| < 1.48$, $p_\mathrm{T}tf < 0.08$, $\mathrm{HEF} > 0.2$, $t^{\mathrm{RMS}}_\mathrm{jet}/t_{\mathrm{jet}} < 0.4$, $t_{\mathrm{jet}} > 3\,\mathrm{ns}$, $E_{\mathrm{ECAL}} > 20\,\mathrm{GeV}$ and failing the HCAL noise rejection quality filters) and for signal jets passing signal selections (except on $\mathrm{max}(\Delta \phi_{\mathrm{DT}})$).

The distribution (normalized to unity) of $\mathrm{max}(\Delta \phi_{\mathrm{RPC}})$ for a data sample enriched in cosmic muons (satisfying $E^{\mathrm{CSC}}_\mathrm{ECAL}/E_{\mathrm{ECAL}} < 0.8$, $p_\mathrm{T} > 30\,\mathrm{GeV}$, $|\eta| < 1.48$, $p_\mathrm{T}tf < 0.08$, $\mathrm{HEF} > 0.2$, $t^{\mathrm{RMS}}_\mathrm{jet}/t_{\mathrm{jet}} < 0.4$, $t_{\mathrm{jet}} > 3\,\mathrm{ns}$, $E_{\mathrm{ECAL}} > 20\,\mathrm{GeV}$ and failing the HCAL noise rejection quality filters) and for signal jets passing signal selections (except on $\mathrm{max}(\Delta \phi_{\mathrm{RPC}})$).

Contribution to the delayed time of jets from the $\beta$ of the gluino is plotted against the delay contribution from the difference (assuming straight line paths) between the path taken by the gluino and the particle forming the jet from the path length for a particle travelling directly to the same position on the ECAL barrel for gluino $c\tau_{0} = 10$ m and mass of 1000 Gev (top) and 3000 GeV (bottom). The dominant contribution is shown to be the gluino $\beta$.

The observed upper limits at $95\%$ CL for the gluino pair production cross section relative to the theoretical cross section in the GMSB model, shown in the plane of $m_{\tilde{g}}$ and $c\tau_{0}$. A branching fraction of $100\%$ for the gluino decay to a gluon and a gravitino is assumed. The area below the thick black curve represents the observed exclusion region, while the dashed red lines indicate the expected limits and their $\pm 1$ standard deviation ranges from experimental uncertainties. The thin black lines show the effect of the theoretical uncertainties on the signal cross section.

The expected upper limits at $95\%$ CL for the gluino pair production cross section relative to the theoretical cross section in the GMSB model, shown in the plane of $m_{\tilde{g}}$ and $c\tau_{0}$. A branching fraction of $100\%$ for the gluino decay to a gluon and a gravitino is assumed. The area below the thick black curve represents the observed exclusion region, while the dashed red lines indicate the expected limits and their $\pm 1$ standard deviation ranges from experimental uncertainties. The thin black lines show the effect of the theoretical uncertainties on the signal cross section.

The expected plus 1 $\sigma_{\mathrm{experiment}}$ upper limits at $95\%$ CL for the gluino pair production cross section relative to the theoretical cross section in the GMSB model, shown in the plane of $m_{\tilde{g}}$ and $c\tau_{0}$. A branching fraction of $100\%$ for the gluino decay to a gluon and a gravitino is assumed. The area below the thick black curve represents the observed exclusion region, while the dashed red lines indicate the expected limits and their $\pm 1$ standard deviation ranges from experimental uncertainties. The thin black lines show the effect of the theoretical uncertainties on the signal cross section.

The expected minus 1 $\sigma_{\mathrm{experiment}}$ upper limits at $95\%$ CL for the gluino pair production cross section relative to the theoretical cross section in the GMSB model, shown in the plane of $m_{\tilde{g}}$ and $c\tau_{0}$. A branching fraction of $100\%$ for the gluino decay to a gluon and a gravitino is assumed. The area below the thick black curve represents the observed exclusion region, while the dashed red lines indicate the expected limits and their $\pm 1$ standard deviation ranges from experimental uncertainties. The thin black lines show the effect of the theoretical uncertainties on the signal cross section.

The expected upper limits at $95\%$ CL for the gluino pair production cross section in the GMSB model, shown in the plane of $m_{\tilde{g}}$ and $c\tau_{0}$. A branching fraction of $100\%$ for the gluino decay to a gluon and a gravitino is assumed. The area below the thick black curve represents the observed exclusion region, while the dashed red lines indicate the expected limits and their $\pm 1$ standard deviation ranges from experimental uncertainties. The thin black lines show the effect of the theoretical uncertainties on the signal cross section.

The expected plus 1 $\sigma_{\mathrm{experiment}}$ upper limits at $95\%$ CL for the gluino pair production cross section in the GMSB model, shown in the plane of $m_{\tilde{g}}$ and $c\tau_{0}$. A branching fraction of $100\%$ for the gluino decay to a gluon and a gravitino is assumed. The area below the thick black curve represents the observed exclusion region, while the dashed red lines indicate the expected limits and their $\pm 1$ standard deviation ranges from experimental uncertainties. The thin black lines show the effect of the theoretical uncertainties on the signal cross section.

The expected minus 1 $\sigma_{\mathrm{experiment}}$ upper limits at $95\%$ CL for the gluino pair production cross section in the GMSB model, shown in the plane of $m_{\tilde{g}}$ and $c\tau_{0}$. A branching fraction of $100\%$ for the gluino decay to a gluon and a gravitino is assumed. The area below the thick black curve represents the observed exclusion region, while the dashed red lines indicate the expected limits and their $\pm 1$ standard deviation ranges from experimental uncertainties. The thin black lines show the effect of the theoretical uncertainties on the signal cross section.

Observed upper limits at 95% CL on the fiducial cross-section times branching ratio as a function of pole mass for the zero-width, 0.5%, 1.2%, 3%, 6% and 10% relative width signals for the combined dilepton channel.

Expected upper limits at 95% CL on the fiducial cross-section times branching ratio as a function of pole mass for the zero-width, 0.5%, 1.2%, 3%, 6% and 10% relative width signals for the dielectron channel.

Observed upper limits at 95% CL on the fiducial cross-section times branching ratio as a function of pole mass for the zero-width, 0.5%, 1.2%, 3%, 6% and 10% relative width signals for the dielectron channel.

Expected upper limits at 95% CL on the fiducial cross-section times branching ratio as a function of pole mass for the zero-width, 0.5%, 1.2%, 3%, 6% and 10% relative width signals for the dimuon channel.

Observed upper limits at 95% CL on the fiducial cross-section times branching ratio as a function of pole mass for the zero-width, 0.5%, 1.2%, 3%, 6% and 10% relative width signals for the dimuon channel.

Ratio of the observed limit to the Z'$_\mathrm{SSM}$ cross-section for the combination of the dielectron and dimuon channels. A number of previous ATLAS results at $\sqrt{s}$ = 7, 8, and 13 TeV are compared to the current search.

Observed 95% exclusion contours in the HVT parameter space $\{g_h,g_f\}$ with $g_f\equiv g_l=g_q$ for a resonance mass of 3 TeV for the dilepton channel. The area outside the curves is excluded.

Observed 95% exclusion contours in the HVT parameter space $\{g_h,g_f\}$ with $g_f\equiv g_l=g_q$ for a resonance mass of 4 TeV for the dilepton channel. The area outside the curves is excluded.

Observed 95% exclusion contours in the HVT parameter space $\{g_h,g_f\}$ with $g_f\equiv g_l=g_q$ for a resonance mass of 5 TeV for the dilepton channel. The area outside the curves is excluded.

Observed 95% exclusion contours in the HVT parameter space $\{g_q,g_l\}$ with $g_h$ set to zero for a resonance mass of 3 TeV for the dilepton channel. The area outside the curves is excluded.

Observed 95% exclusion contours in the HVT parameter space $\{g_q,g_l\}$ with $g_h$ set to zero for a resonance mass of 4 TeV for the dilepton channel. The area outside the curves is excluded.

Observed 95% exclusion contours in the HVT parameter space $\{g_q,g_l\}$ with $g_h$ set to zero for a resonance mass of 5 TeV for the dilepton channel. The area outside the curves is excluded.

Expected 95% exclusion contours in the HVT parameter space $\{g_h,g_f\}$ with $g_f\equiv g_l=g_q$ for a resonance mass of 3 TeV for the dilepton channel. The area outside the curves is excluded.

Expected 95% exclusion contours in the HVT parameter space $\{g_h,g_f\}$ with $g_f\equiv g_l=g_q$ for a resonance mass of 4 TeV for the dilepton channel. The area outside the curves is excluded.

Expected 95% exclusion contours in the HVT parameter space $\{g_h,g_f\}$ with $g_f\equiv g_l=g_q$ for a resonance mass of 5 TeV for the dilepton channel. The area outside the curves is excluded.

Expected 95% exclusion contours in the HVT parameter space $\{g_q,g_l\}$ with $g_h$ set to zero for a resonance masses of 3 TeV for the dilepton channel. The area outside the curves is excluded.

Expected 95% exclusion contours in the HVT parameter space $\{g_q,g_l\}$ with $g_h$ set to zero for a resonance masses of 4 TeV for the dilepton channel. The area outside the curves is excluded.

Expected 95% exclusion contours in the HVT parameter space $\{g_q,g_l\}$ with $g_h$ set to zero for a resonance masses of 5 TeV for the dilepton channel. The area outside the curves is excluded.

Dielectron and dimuon selection efficiency as a function of pole mass for spin-0, spin-1 and spin-2 resonances. Drell-Yan MC samples are used for the spin-1 efficiencies and correspond to the zero-width signal hypothesis. Spin-0 (spin-2) efficiencies are obtained from dedicated MC samples with relative widths in the range of 0-20% (1-13%) and are averaged over all available width assumptions.

Description of the $m_{\ell\ell}^{\mathrm{true}}$-dependence for the seven relative mass resolution parameters in the dielectron channel.

Description of the $m_{\ell\ell}^{\mathrm{true}}$-dependence for the seven relative mass resolution parameters in the dimuon channel.

Upper limits at 95% CL on the cross section times branching fraction for the process $H \to 2\gamma_d + X$ with $m_H$ = 800 GeV in the hadron-hadron final state.

Upper limits at 95% CL on the cross section times branching fraction for the process $H \to 2\gamma_d + X$ with $m_H$ = 800 GeV in the muon-hadron final state.

Upper limits at 95% CL on the cross section times branching fraction for the process $H \to 4\gamma_d + X$ with $m_H$ = 800 GeV in the muon-muon final state.

Upper limits at 95% CL on the cross section times branching fraction for the process $H \to 4\gamma_d + X$ with $m_H$ = 800 GeV in the muon-hadron final state.

Upper limits at 95% CL on the cross section times branching fraction for the process $H \to 4\gamma_d + X$ with $m_H$ = 800 GeV in the hadron-hadron final state.

Extrapolated signal efficiencies as a function of proper decay length of the dark photon for the $H \to 4\gamma_d + X$ process with $m_H$ = 125 GeV in the muon-muon final state.

Extrapolated signal efficiencies as a function of proper decay length of the dark photon for the $H \to 2\gamma_d + X$ process with $m_H$ = 125 GeV in the muon-muon final state.

Extrapolated signal efficiencies as a function of proper decay length of the dark photon for the $H \to 2\gamma_d + X$ process with $m_H$ = 800 GeV in the muon-muon final state.

Extrapolated signal efficiencies as a function of proper decay length of the dark photon for the $H \to 4\gamma_d + X$ process with $m_H$ = 800 GeV in the muon-muon final state.

Extrapolated signal efficiencies as a function of proper decay length of the dark photon for the $H \to 4\gamma_d + X$ process with $m_H$ = 125 GeV in the muon-hadron final state.

Extrapolated signal efficiencies as a function of proper decay length of the dark photon for the $H \to 2\gamma_d + X$ process with $m_H$ = 125 GeV in the muon-hadron final state.

Extrapolated signal efficiencies as a function of proper decay length of the dark photon for the $H \to 2\gamma_d + X$ process with $m_H$ = 800 GeV in the muon-hadron final state.

Extrapolated signal efficiencies as a function of proper decay length of the dark photon for the $H \to 4\gamma_d + X$ process with $m_H$ = 800 GeV in the muon-hadron final state.

Extrapolated signal efficiencies as a function of proper decay length of the dark photon for the $H \to 2\gamma_d + X$ process with $m_H$ = 125 GeV in the hadron-hadron final state.

Extrapolated signal efficiencies as a function of proper decay length of the dark photon for the $H \to 2\gamma_d + X$ process with $m_H$ = 800 GeV in the hadron-hadron final state.

Extrapolated signal efficiencies as a function of proper decay length of the dark photon for the $H \to 4\gamma_d + X$ process with $m_H$ = 800 GeV in the hadron-hadron final state.

Dijet invariant mass distribution for the category with at least one b-tagged jet.

Dijet invariant mass distribution for the category with both jets b-tagged.

The 95% CL upper limits on the cross-section times acceptance times branching ratio into two jets as a function of the mass of q* signal.

The 95% CL upper limits on the cross-section times acceptance times branching ratio into two jets as a function of the mass of QBH signal.

The 95% CL upper limits on the cross-section times acceptance times branching ratio into two jets as a function of the mass of W' signal.

The 95% CL upper limits on the cross-section times acceptance times branching ratio into two jets as a function of the mass of W* signal.

The upper limits on the DM mediator Z' signal at 95% CL from the inclusive category. The 95% CL upper limits are set on the universal quark coupling $g_q$ as a function of the Z' mass.

The 95% CL upper limit on the cross-section times acceptance times b-tagging efficiency times branching ratio as a function of the mass of the b* signal.

The 95% CL upper limit on the cross-section times acceptance times b-tagging efficiency times branching ratio as a function of the signal mass in the DM mediator Z' with $g_q$ = 0.25 model.

The 95% CL upper limit on the cross-section times acceptance times b-tagging efficiency times branching ratio as a function of the signal mass in the SSM Z' model.

The 95% CL upper limit on the cross-section times acceptance times b-tagging efficiency times branching ratio as a function of the signal mass in the graviton with $k/\overline{M}_{PL}=0.2$ model.

The 95% CL upper limit on the cross-section times kinematic acceptance times branching ratio for resonances with a generic Gaussian shape, as a function of the Gaussian mean mass in the inclusive category. Different widths, from 0% up to 15% of the signal mass, are considered. Gaussian-shape signals with 0% widths correspond to signal widths smaller than the experimental resolution. For a Gaussian-shaped signal with a relative width of 15%, the limits are truncated at high mass when the broad signal starts to overlap the upper end of the mass spectrum.

The 95% CL upper limit on the cross-section times kinematic acceptance times b-tagging efficiency times branching ratio for resonances with a generic Gaussian shape, as a function of the Gaussian mean mass in the 1 b category. Different widths, from 0% up to 15% of the signal mass, are considered. Gaussian-shape signals with 0% widths correspond to signal widths smaller than the experimental resolution. For a Gaussian-shaped signal with a relative width of 15%, the limits are truncated at high mass when the broad signal starts to overlap the upper end of the mass spectrum.

The 95% CL upper limit on the cross-section times kinematic acceptance times b-tagging efficiency times branching ratio for resonances with a generic Gaussian shape, as a function of the Gaussian mean mass in the 2 b category. Different widths, from 0% up to 15% of the signal mass, are considered. Gaussian-shape signals with 0% widths correspond to signal widths smaller than the experimental resolution. For a Gaussian-shaped signal with a relative width of 15%, the limits are truncated at high mass when the broad signal starts to overlap the upper end of the mass spectrum.

The expected 95% CL upper limits on the cross-section times acceptance times b-tagging efficiency times branching ratio as a function of the DM mediator Z' mass for the current and previous iterations of the analysis. The upper limit of the previous result was obtained with the Bayesian method and is also shown scaled to the 139 fb$^{-1}$ integrated luminosity of the current result to illustrate the effect of the analysis improvements.

Acceptance for the QBH, Q* and W' benchmark signal models in the inclusive category, as a function of the simulated mass.

Acceptance for the DM Z' benchmark signal model for various $g_q$ coupling parameters in the inclusive category, as a function of the simulated mass.

Acceptance for the W* benchmark signal model in the inclusive category, as a function of the simulated mass.

Acceptance times b-tagging efficiency for the b* benchmark signal model in the 1b category, as a function of the simulated mass.

Acceptance times b-tagging efficiency for the DM Z' with $g_q$=0.25, SSM Z' and graviton with $k/\overline{M}_{PL}=0.2$ benchmark signal models in the 2b category, as a function of the simulated mass.

The expected 95% CL upper limits on the cross-section times branching ratio as a function of the DM mediator Z' mass for the current and previous iterations of the analysis. The upper limit of the previous result was obtained with the Bayesian method and is also shown scaled to the 139 fb$^{-1}$ integrated luminosity of the current result to illustrate the effect of the analysis improvements. The current b-tagging requirement is tighter than the previous one for high-$p_T$ jets, resulting in a data sample with limited size for mass above 4 TeV. The background rejection, instead, has improved significantly across the entire mass spectrum inspected by the analysis.

Upper limits at 95% CL on the cross-section of the resonant Higgs boson pair production for a spin-0 heavy scalar

Upper limits at 95% CL on the cross-section of the resonant Higgs boson pair production for a spin-0 heavy scalar

Upper limits at 95% CL on the cross-section of the resonant Higgs boson pair production for a spin-0 heavy scalar

Upper limits at 95% CL on the cross-section of the resonant Higgs boson pair production for a spin-0 heavy scalar

Upper limits at 95% CL on the cross-section of the resonant Higgs boson pair production for a spin-0 heavy scalar

Upper limits at 95% CL on the cross-section of the resonant Higgs boson pair production for a spin-0 heavy scalar

Upper limits at 95% CL on the cross-section of the resonant Higgs boson pair production for a spin-0 heavy scalar

Upper limits at 95% CL on the cross-section of the resonant Higgs boson pair production for a spin-0 heavy scalar

Upper limits at 95% CL on the cross-section of the resonant Higgs boson pair production for a spin-2 KK graviton with k/M&#772;_Pl = 1

Upper limits at 95% CL on the cross-section of the resonant Higgs boson pair production for a spin-2 KK graviton with k/M&#772;_Pl = 1

Upper limits at 95% CL on the cross-section of the resonant Higgs boson pair production for a spin-2 KK graviton with k/M&#772;_Pl = 1

Upper limits at 95% CL on the cross-section of the resonant Higgs boson pair production for a spin-2 KK graviton with k/M&#772;_Pl = 1

Upper limits at 95% CL on the cross-section of the resonant Higgs boson pair production for a spin-2 KK graviton with k/M&#772;_Pl = 1

Upper limits at 95% CL on the cross-section of the resonant Higgs boson pair production for a spin-2 KK graviton with k/M&#772;_Pl = 1

Upper limits at 95% CL on the cross-section of the resonant Higgs boson pair production for a spin-2 KK graviton with k/M&#772;_Pl = 2

Upper limits at 95% CL on the cross-section of the resonant Higgs boson pair production for a spin-2 KK graviton with k/M&#772;_Pl = 2

Upper limits at 95% CL on the cross-section of the resonant Higgs boson pair production for a spin-2 KK graviton with k/M&#772;_Pl = 2

Upper limits at 95% CL on the cross-section of the resonant Higgs boson pair production for a spin-2 KK graviton with k/M&#772;_Pl = 2

Upper limits at 95% CL on the cross-section of the resonant Higgs boson pair production for a spin-2 KK graviton with k/M&#772;_Pl = 2

Upper limits at 95% CL on the cross-section of the resonant Higgs boson pair production for a spin-2 KK graviton with k/M&#772;_Pl = 2

Expected exclusion regions for the EWK-singlet model for tan &beta; = 1. The theoretical contour points with width over mass of the resonance equal to 10% are included in the table to identify the experimental-invalid regions.

Observed exclusion regions for the EWK-singlet model for tan &beta; = 1. The theoretical contour points with width over mass of the resonance equal to 10% are included in the table to identify the experimental-invalid regions.

Expected exclusion regions for the EWK-singlet model for tan &beta; = 2. The theoretical contour points with width over mass of the resonance equal to 10% are included in the table to identify the experimental-invalid regions.

Observed exclusion regions for the EWK-singlet model for tan &beta; = 2. The theoretical contour points with width over mass of the resonance equal to 10% are included in the table to identify the experimental-invalid regions.

Expected exclusion regions in the EWK-singlet model for m_S= 260 GeV. The theoretical contour points with width over mass of the resonance equal to 10% are included in the table to identify the experimental-invalid regions.

Observed exclusion regions in the EWK-singlet model for m_S= 260 GeV. The theoretical contour points with width over mass of the resonance equal to 10% are included in the table to identify the experimental-invalid regions.

Expected exclusion regions in the hMSSM model

Expected exclusion regions (-1 sigma band) in the hMSSM model

Expected exclusion regions (+1 sigma band) in the hMSSM model

Observed exclusion regions in the hMSSM model

Upper limits at 95% CL on the signal cross-section of the ggF non-resonant Higgs boson pair production as a function of &kappa;_&lambda; with statistical uncertainties only. The limits are shown in dashed lines. The combined limit including all systematic uncertainties is also shown with a solid line.

Expected exclusion regions for the EWK-singlet model with tan &beta; = 5.0. The theoretical contour points with width over mass of the resonance equal to 10% are included in the table to identify the experimental-invalid regions.

Observed exclusion regions for the EWK-singlet model with tan &beta; = 5.0. The theoretical contour points with width over mass of the resonance equal to 10% are included in the table to identify the experimental-invalid regions.

Expected exclusion regions for the EWK-singlet model with m_S = 300 GeV. The theoretical contour points with width over mass of the resonance equal to 10% are included in the table to identify the experimental-invalid regions.

Observed exclusion regions for the EWK-singlet model with m_S = 300 GeV. The theoretical contour points with width over mass of the resonance equal to 10% are included in the table to identify the experimental-invalid regions.

Expected exclusion regions for the EWK-singlet model with m_S = 400 GeV. The theoretical contour points with width over mass of the resonance equal to 10% are included in the table to identify the experimental-invalid regions.

Observed exclusion regions for the EWK-singlet model with m_S = 400 GeV. The theoretical contour points with width over mass of the resonance equal to 10% are included in the table to identify the experimental-invalid regions.

Expected exclusion regions for the EWK-singlet model with m_S = 500 GeV. The theoretical contour points with width over mass of the resonance equal to 10% are included in the table to identify the experimental-invalid regions.

Observed exclusion regions for the EWK-singlet model with m_S = 500 GeV. The theoretical contour points with width over mass of the resonance equal to 10% are included in the table to identify the experimental-invalid regions.

The rho distributions for the mixed-flavour category events corresponding to a signal hypothesis with mH = 261 GeV and mA = 150 GeV.

The rho distributions for the same-flavour category events corresponding to a signal hypothesis with mH = 442 GeV and mA = 193 GeV. The signal is normalised to its theoretical cross section.

The rho distributions for the mixed-flavour category events corresponding to a signal hypothesis with mH = 442 GeV and mA = 193 GeV.

Expected and observed 95% CL upper limits on the product of the production cross section and branching fraction for H(A) -> ZA(H) -> bbbar as a function of mA and mH. The limits are computed using the asymptotic CLs method, combining the ee and μμ channels.

Expected and observed 95% CL upper limits on the signal strength for the Type-II 2HDM benchmark (tan(beta)=1.5, cos(beta-alpha)=0.01) as a function of mA and mH . The limits are computed using the asymptotic CLs method, combining the ee and μμ channels.

Expected and observed 95% CL upper limits on the signal strength for the Type-II 2HDM benchmark (mH = 379 GeV and m A = 172 GeV) as a function of tan(beta) and cos(beta-alpha). The limits are computed using the asymptotic CLs method, combining the ee and μμ channels.

Exclusion regions at 95% CL in the $m_{\tilde{\chi}_1^0}$ versus $m_{\tilde{g}}$ plane for the T1tttt (upper left) and T5ttbbWW (upper right) models, with off-shell third-generation squarks, and the T5tttt (lower left) and T5ttcc (lower right) models, with on-shell third-generation squarks. For the T5ttbbWW model, $m_{\tilde{\chi}_1^\pm} = m_{\tilde{\chi}_1^0} + 5 GeV$, for the T5tttt model, $m_{\tilde{t}} - m_{\tilde{\chi}_1^0} = m_t$, and for the T5ttcc model, $m_{\tilde{t}} - m_{\tilde{\chi}_1^0} = 20 GeV$ and the decay proceeds through $\tilde{t} \to c \tilde{\chi}_1^0$. The right-hand side color scale indicates the excluded cross section values for a given point in the SUSY particle mass plane. The solid black curves represent the observed exclusion limits assuming the approximate-NNLO+NNLL cross sections (thick line), or their variations of $\pm 1$ standard deviations (s.d.) (thin lines). The dashed red curves show the expected limits with the corresponding $\pm 1$ s.d. and $\pm 2$ s.d. uncertainties. Excluded regions are to the left and below the limit curves.

Exclusion regions at 95% CL in the plane of $m_{\tilde{\chi}_1^0}$ versus $m_{\tilde{g}}$ for the T5qqqqWZ model with $m_{\tilde{\chi}_1^\pm}=0.5(m_{\tilde{g}} + m_{\tilde{\chi}_1^0})$~(left) and with $m_{\tilde{\chi}_1^\pm} = m_{\tilde{\chi}_1^0} + 20 GeV$~(right). The notations are as in Fig. 7.

Exclusion regions at 95% CL in the plane of $m_{\tilde{\chi}_1^0}$ versus $m_{\tilde{g}}$ for the T5qqqqWZ model with $m_{\tilde{\chi}_1^\pm}=0.5(m_{\tilde{g}} + m_{\tilde{\chi}_1^0})$~(left) and with $m_{\tilde{\chi}_1^\pm} = m_{\tilde{\chi}_1^0} + 20 GeV$~(right). The notations are as in Fig. 7.

Exclusion regions at 95% CL in the plane of $m_{\tilde{\chi}_1^0}$ versus $m_{\tilde{g}}$ for the T5qqqqWW model with $m_{\tilde{\chi}_1^\pm}=0.5(m_{\tilde{g}} + m_{\tilde{\chi}_1^0})$~(left) and with $m_{\tilde{\chi}_1^\pm} = m_{\tilde{\chi}_1^0} + 20 GeV$~(right). The notations are as in Fig. 7.

Exclusion regions at 95% CL in the plane of $m_{\tilde{\chi}_1^0}$ versus $m_{\tilde{g}}$ for the T5qqqqWW model with $m_{\tilde{\chi}_1^\pm}=0.5(m_{\tilde{g}} + m_{\tilde{\chi}_1^0})$~(left) and with $m_{\tilde{\chi}_1^\pm} = m_{\tilde{\chi}_1^0} + 20 GeV$~(right). The notations are as in Fig. 7.

Exclusion regions at 95% CL in the plane of $m_{\tilde{\chi}_1^\pm}$ versus $m_{\tilde{b}_1}$ for the T6ttWW model with $m_{\tilde{\chi}_1^0}=50 GeV$. The notations are as in Fig. 7.

Exclusion regions at 95% CL in the plane of $m(\tilde{t}_1)$ versus $m(\tilde{t}_2)$ for the T6ttHZ model with $m(\tilde{t}_1)-m(\tilde{\chi}_1^0)=175 GeV$. The three exclusions represent $\mathcal{B}(\tilde{t}_2\to\tilde{t}_1 Z)$ of 0, 50, and 100%, respectively. The notations are as in Fig. 7.

Exclusion regions at 95% CL in the plane of $m(\tilde{t}_1)$ versus $m(\tilde{t}_2)$ for the T6ttHZ model with $m(\tilde{t}_1)-m(\tilde{\chi}_1^0)=175 GeV$. The three exclusions represent $\mathcal{B}(\tilde{t}_2\to\tilde{t}_1 Z)$ of 0, 50, and 100%, respectively. The notations are as in Fig. 7.

Exclusion regions at 95% CL in the plane of $m(\tilde{t}_1)$ versus $m(\tilde{t}_2)$ for the T6ttHZ model with $m(\tilde{t}_1)-m(\tilde{\chi}_1^0)=175 GeV$. The three exclusions represent $\mathcal{B}(\tilde{t}_2\to\tilde{t}_1 Z)$ of 0, 50, and 100%, respectively. The notations are as in Fig. 7.

Upper limits at 95% CL on the cross section for RPV gluino pair production with each gluino decaying into four quarks and one lepton (T1qqqqL, left), and each gluino decaying into a top, bottom, and strange quarks (T1tbs, right).

Upper limits at 95% CL on the cross section for RPV gluino pair production with each gluino decaying into four quarks and one lepton (T1qqqqL, left), and each gluino decaying into a top, bottom, and strange quarks (T1tbs, right).

Upper limits at 95% CL on the product of cross section, detector acceptance, and selection efficiency, $\sigma \mathcal{A} \epsilon$, for the production of an SS lepton pair with at least two jets, as a function of the minimum $p_\mathrm{T}^\mathrm{miss}$ threshold, when $H_\mathrm{T}>300 GeV$ (left), or the minimum $H_\mathrm{T}$ threshold, when $p_\mathrm{T}^\mathrm{miss}<300 GeV$ (right).

Upper limits at 95% CL on the product of cross section, detector acceptance, and selection efficiency, $\sigma \mathcal{A} \epsilon$, for the production of an SS lepton pair with at least two jets, as a function of the minimum $p_\mathrm{T}^\mathrm{miss}$ threshold, when $H_\mathrm{T}>300 GeV$ (left), or the minimum $H_\mathrm{T}$ threshold, when $p_\mathrm{T}^\mathrm{miss}<300 GeV$ (right).