Signal selection efficiency times acceptance for charged Higgs boson with masses between $60\,$GeV to $168\,$GeV. The efficiency times acceptance drops at high masses rapidly due to the low momentum of the $b$-quark produced together with the charged Higgs boson. For the $168\,$GeV mass point the efficiency times acceptance is larger than for the $160\,$GeV mass point due to having a higher probability that the top quark, which decays into the charged Higgs boson and a $b$-quark, is off-shell.

Best fit values of $\mu_{t\bar{t}}$ and $f_{\mathrm{HF}}$ parameters at charged Higgs masses between $60\,$GeV and $168\,$GeV.

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.

Expected and observed events in the 16 discovery regions along with the according control regions.

Expected signal yield and acceptance x efficiency, estimated background and observed number of events in data for the full range of simulated masses in the MS-agnostic R-hadron search.

Expected signal yield and acceptance x efficiency, estimated background and observed number of events in data for the full range of simulated masses in the full-detector R-hadron search.

p0-values and model-independent upper limits on cross-section x acceptance x efficiency for the 16 discovery regions.

Expected signal yield and acceptance x efficiency, estimated background and observed number of events in data for the full range of simulated masses in the MS-agnostic search for metastable gluino R-hadrons.

Expected signal yield and acceptance x efficiency, estimated background and observed number of events in data for the full range of simulated masses in the full-detector direct-stau search.

Expected signal yield and acceptance x efficiency, estimated background and observed number of events in data for the full range of simulated masses in the full-detector chargino search.

Upper cross-section limit in gluino R-hadron search.

Upper cross-section limit in sbottom R-hadron search.

Upper cross-section limit in stop R-hadron search.

Upper cross-section limit in stau search.

Upper cross-section limit in chargino search.

Lower mass limit as function of gluino lifetime.

Acceptance x efficiency, acceptance and efficiency for the full range of simulated masses in the MS-agnostic R-hadron search.

Upper cross-section limit in meta-stable gluino R-hadron search.

Flavor composition of 800 GeV stop R-hadrons simulated using the generic model as a function of radial distance from the interaction point.

Flavor composition of 800 GeV anti-stop R-hadrons simulated using the generic model as a function of radial distance from the interaction point.

Flavor composition of 800 GeV stop R-hadrons simulated using the Regge model as a function of radial distance from the interaction point.

Flavor composition of 800 GeV anti-stop R-hadrons simulated using the Regge model as a function of radial distance from the interaction point.

ETmiss trigger efficiency as function of true ETmiss (EtmissTurnOn).

Single-muon trigger efficiency as function of $|\eta|$ and $\beta$ (SingleMuTurnOn).

Candidate reconstruction efficiency for ID+Calo selection (IDCaloEff).

Candidate reconstruction efficiency for loose selection (LooseEff).

Efficiency for a loose candidate to be promoted to a tight candidate (TightPromotionEff).

Resolution and average of reconstructed dE/dx mass for a given simulated mass for ID+calo candidates.

Resolution and average of reconstructed ToF mass for a given simulated mass for ID+calo candidates.

Resolution and average of reconstructed ToF mass for a given simulated mass for FullDet candidates.

2$\tau$ HighMass SR eff.

2$\tau$ multibin SR eff.

2$\tau$ GMSB SR eff.

1$\tau$ Compressed SR eff.

1$\tau$ MediumMass SR eff.

2$\tau$ Compressed SR eff.

2$\tau$ HighMass SR eff.

2$\tau$ multibin SR eff.

2$\tau$ GMSB SR eff.

1$\tau$ Compressed SR acceptance.

1$\tau$ MediumMass SR acceptance.

2$\tau$ Compressed SR acceptance.

2$\tau$ HighMass SR acceptance.

2$\tau$ multibin SR acceptance.

2$\tau$ GMSB SR acceptance.

1$\tau$ Compressed SR acceptance.

1$\tau$ MediumMass SR acceptance.

2$\tau$ Compressed SR acceptance.

2$\tau$ HighMass SR acceptance.

2$\tau$ multibin SR acceptance.

2$\tau$ GMSB SR acceptance.

Cutflow table of the $1\tau$ compressed SR for the four signal benchmark scenarios of low, medium, and high mass-splitting in the simplified model as well as the GMSB model.

Cutflow table of the $1\tau$ medium-mass SR for the four signal benchmark scenarios of low, medium, and high mass-splitting in the simplified model as well as the GMSB model.

Cutflow table of the $2\tau$ compressed SR for the four signal benchmark scenarios of low, medium, and high mass-splitting in the simplified model as well as the GMSB model.

Cutflow table of the $2\tau$ high-mass SR for the four signal benchmark scenarios of low, medium, and high mass-splitting in the simplified model as well as the GMSB model.

Cutflow table of the $2\tau$ multibin SR for the four signal benchmark scenarios of low, medium, and high mass-splitting in the simplified model as well as the GMSB model.

Cutflow table of the $2\tau$ GMSB SR for the four signal benchmark scenarios of low, medium, and high mass-splitting in the simplified model as well as the GMSB model.

Best performing fit setups entering the final combination as a function of the LSP mass and the gluino mass. 'S' marks the simultaneous fit of the four simplified model single-bin SRs, 'M' denotes the simultaneous fit of the two $1\tau$ SRs and the $2\tau$ multibin SR.

Observed exclusion contour at 95% CL as a function of tanBeta and the SUSY-breaking mass scale Lambda.

Expected exclusion contour at 95% CL as a function of tanBeta and the SUSY-breaking mass scale Lambda.

Observed exclusion contour at 95% CL as a function of the LSP mass and the gluino mass.

Expected exclusion contour at 95% CL as a function of the LSP mass and the gluino mass.

Observed upper limits on the production cross section at 95% CL in pb as a function of tanBeta and SUSY breaking mass scale Lambda.

Observed upper limits on the production cross section at 95% CL in pb as a function of the LSP mass and the gluino mass.

Yields of the expected background from the SM in the bins of the multibin SR of the $2\tau$ channel with all bins being simultaneously used to constrain the background prediction. Expectation is given with the scalings computed in the combined fit applied. Uncertainties are statistial plus systematrics. Only the subsamples contributing the respective region are considered.

$m_{\mathrm{T}}^{\tau}$ in the compressed $m_{\mathrm{T}}^{\tau}$ VR of the $1\tau$ channel, illustrating the background modeling after the fit. The last bin includes overflow events.

$E_{\mathrm{T}}^{\mathrm{miss}}$ in the compressed $E_{\mathrm{T}}^{\mathrm{miss}}$ VR of the $1\tau$ channel, illustrating the background modeling after the fit. The last bin includes overflow events.

$m_{\mathrm{T}}^{\tau}$ in the medium-mass $m_{\mathrm{T}}^{\tau}$ VR of the $1\tau$ channel, illustrating the background modeling after the fit. The last bin includes overflow events.

$E_{\mathrm{T}}^{\mathrm{miss}}$ in the medium-mass $E_{\mathrm{T}}^{\mathrm{miss}}$ VR of the $1\tau$ channel, illustrating the background modeling after the fit. The last bin includes overflow events.

$H_{\mathrm{T}}$ in the medium-mass $H_{\mathrm{T}}$ VR of the $1\tau$ channel, illustrating the background modeling after the fit. The last bin includes overflow events.

$m_{\mathrm{T}}^{\tau_1}$ + $m_{\mathrm{T}}^{\tau_2}$ in the top VR of the $2\tau$ channel, illustrating the background modeling after the fit. The last bin includes overflow events.

$H_{\mathrm{T}}$ in the $W$ VR of the $2\tau$ channel, illustrating the background modeling after the fit. The last bin includes overflow events.

$m_{\mathrm{T}}^{\tau_1}$ + $m_{\mathrm{T}}^{\tau_2}$ in the $Z$ VR of the $2\tau$ channel, illustrating the background modeling after the fit. The last bin includes overflow events.

$m_{\mathrm{T}}^{\tau}$ in the compressed SR of the $1\tau$ channel before application of the $m_{\mathrm{T}}^{\tau}$ > 80 GeV requirement. The last bin includes overflow events. Signal predictions corresponding to the simplified model scenarios of low (LM), medium (MM), and high mass-splitting (HM) as well as for the GMSB benchmark are given.

$H_{\mathrm{T}}$ in the medium-mass SR of the $1\tau$ channel before application of the $H_{\mathrm{T}}$ > 1000 GeV requirement. The last bin includes overflow events. Signal predictions corresponding to the simplified model scenarios of low (LM), medium (MM), and high mass-splitting (HM) as well as for the GMSB benchmark are given.

$m_{\mathrm{T}}^{\mathrm{sum}}$ in the compressed SR of the $2\tau$ channel before application of the $m_{\mathrm{T}}^{\mathrm{sum}}$ > 1600 GeV requirement. The last bin includes overflow events. Signal predictions corresponding to the simplified model scenarios of low (LM), medium (MM), and high mass-splitting (HM) as well as for the GMSB benchmark are given.

$H_{\mathrm{T}}$ in the high-mass SR of the $2\tau$ channel before application of the $H_{\mathrm{T}}$ > 1100 GeV requirement. The last bin includes overflow events. Signal predictions corresponding to the simplified model scenarios of low (LM), medium (MM), and high mass-splitting (HM) as well as for the GMSB benchmark are given.

mT(tau_1) + mT(tau_2) in the multibin SR of the 2T channel. The last bin includes overflow events. Signal predictions corresponding to the simplified model scenarios of low (LM), medium (MM), and high mass-splitting (HM) as well as for the GMSB benchmark are given.

$H_{\mathrm{T}}$ in the GMSB SR of the $2\tau$ channel before application of the $H_{\mathrm{T}}$ > 1900 GeV requirement. The last bin includes overflow events. Signal predictions corresponding to the simplified model scenarios of low (LM), medium (MM), and high mass-splitting (HM) as well as for the GMSB benchmark are given.

Pixel-tracklet $p_{T}$ spectrum of signal ($\tau_{\tilde{\chi}_{1}^{\pm}}$ = 0.2 ns and $m_{\tilde{\chi}_{1}^{\pm}}$ = 400 GeV) in electroweak channel in the low-Emiss region.

Pixel-tracklet $p_{T}$ spectrum of total background in electroweak channel in the low-Emiss region.

Pixel-tracklet $p_{T}$ spectrum of observed data in electroweak channel in the low-Emiss region.

Pixel-tracklet $p_{T}$ spectrum of fake tracklet in strong channel in the low-Emiss region.

Pixel-tracklet $p_{T}$ spectrum of muon background in strong channel in the low-Emiss region.

Pixel-tracklet $p_{T}$ spectrum of hadron and electron background in strong channel in the low-Emiss region.

Pixel-tracklet $p_{T}$ spectrum of signal ($m_{\tilde{g}}$ = 1600 GeV, $m_{\tilde{\chi}_{1}^{\pm}} = 500 GeV$) in strong channel in the low-Emiss region.

Pixel-tracklet $p_{T}$ spectrum of total background in strong channel in the low-Emiss region.

Pixel-tracklet $p_{T}$ spectrum of observed data in strong channel in the low-Emiss region.

Pixel-tracklet $p_{T}$ spectrum of fake tracklet in electroweak channel in the high-Emiss region.

Pixel-tracklet $p_{T}$ spectrum of muon background in electroweak channel in the high-Emiss region.

Pixel-tracklet $p_{T}$ spectrum of hadron and electron background in electroweak channel in the high-Emiss region.

Pixel-tracklet $p_{T}$ spectrum of signal ($\tau_{\tilde{\chi}_{1}^{\pm}}$ = 0.2 ns and $m_{\tilde{\chi}_{1}^{\pm}}$ = 400 GeV) in electroweak channel in the high-Emiss region.

Pixel-tracklet $p_{T}$ spectrum of total background in electroweak channel in the high-Emiss region.

Pixel-tracklet $p_{T}$ spectrum of observed data in electroweak channel in the high-Emiss region.

Pixel-tracklet $p_{T}$ spectrum of fake tracklet in strong channel in the high-Emiss region.

Pixel-tracklet $p_{T}$ spectrum of muon background in strong channel in the high-Emiss region.

Pixel-tracklet $p_{T}$ spectrum of hadron and electron background in strong channel in the high-Emiss region.

Pixel-tracklet $p_{T}$ spectrum of signal ($m_{\tilde{g}}$ = 1600 GeV, $m_{\tilde{\chi}_{1}^{\pm}} = 500 GeV$) in strong channel in the high-Emiss region.

Pixel-tracklet $p_{T}$ spectrum of total background in strong channel in the high-Emiss region.

Pixel-tracklet $p_{T}$ spectrum of observed data in strong channel in the high-Emiss region.

Expected exclusion limit at 95% CL obtained in the electroweak production channel in terms of the chargino lifetime ($\tau_{\tilde{\chi}_{1}^{\pm}}$) and mass ($m_{\tilde{\chi}_{1}^{\pm}}$).

Observed exclusion limit at 95% CL obtained in the electroweak production channel in terms of the chargino lifetime ($\tau_{\tilde{\chi}_{1}^{\pm}}$) and mass ($m_{\tilde{\chi}_{1}^{\pm}}$).

Expected exclusion limit at 95% CL obtained in the strong production channel in terms of the gluino and chargino masses. The limit is shown assuming a chargino lifetime 0.2 ns.

Observed exclusion limit at 95% CL obtained in the strong production channel in terms of the gluino and chargino masses. The limit is shown assuming a chargino lifetime 0.2 ns.

Expected exclusion limit at 95% CL obtained in the strong production channel in terms of the gluino and chargino masses. The limit is shown assuming a chargino lifetime 1.0 ns.

Observed exclusion limit at 95% CL obtained in the strong production channel in terms of the gluino and chargino masses. The limit is shown assuming a chargino lifetime 1.0 ns.

Expected exclusion limit at 95% CL obtained in the electroweak production channel in terms of the chargino lifetime ($\tau_{\tilde{\chi}_{1}^{\pm}}$) and mass ($m_{\tilde{\chi}_{1}^{\pm}}$).

Observed exclusion limit at 95% CL obtained in the electroweak production channel in terms of the chargino lifetime ($\tau_{\tilde{\chi}_{1}^{\pm}}$) and mass ($m_{\tilde{\chi}_{1}^{\pm}}$).

Model dependent upper limits on cross-section (pb) for the electroweak production are shown by grey numbers in terms of the chargino lifetime ($\tau_{\tilde{\chi}_{1}^{\pm}}$) and mass ($m_{\tilde{\chi}_{1}^{\pm}}$).

Total acceptance $\times$ efficiency of the electroweak channel. The total signal acceptance $\times$ efficiency is defined as the probability of an event passing the signal region selection when an electroweak gaugino pair is produced in a pp collision.

Total acceptance $\times$ efficiency of the strong channel. In white regions, no simulation sample is available. The left-upper triangle region is not allowed kinematically in wino-LSP scenarios. The total signal acceptance $\times$ efficiency is calculated relative to events in which the gluinos decay into electroweak gaugino pairs.

The generator-level acceptance for charginos produced in the electroweak channel as a function of the chargino $eta$ and chargino decay radius (at generator level).

The acceptance $\times$ efficiency after reconstruction, for selecting and reconstructing charginos produced in the electroweak channel as a function of the chargino $\eta$ and chargino decay radius (at generator level). The acceptance $\times$ efficiency after reconstruction is the probability of a signal event, which passes all the event-level requirements, passing all the track/tracklet requirements after reconstruction.

The generator-level acceptance for charginos produced in the strong channel as a function of the chargino $\eta$ and chargino decay radius (at generator level).

The acceptance $\times$ efficiency after reconstruction, for selecting and reconstructing charginos produced in the strong channel as a function of the chargino $\eta$ and chargino decay radius (at generator level). The acceptance $\times$ efficiency after reconstruction is the probability of a signal event, which passes all the event-level requirements, passing all the track/tracklet requirements after reconstruction.

Summary of the selection criteria, and the corresponding observed number of events in data as well as the expected number of signal events in simulation for two benchmark models: a chargino produced in direct electroweak production with ($m_{\tilde{\chi}_{1}^{\pm}}$, $\tau_{\tilde{\chi}_{1}^{\pm}}$) = (400 GeV, 0.2 ns) in the high-Emiss region. The expected number of signal events is normalised to 36.1 fb${}^{-1}$. The signal selection efficiencies are also shown in parentheses. The first row shows the number of events after the application of detector and data quality conditions. Requirements below the dashed line are applied to tracks and tracklets.

Summary of the selection criteria, and the corresponding observed number of events in data as well as the expected number of signal events in simulation for two benchmark models: a chargino produced in the strong channel with ($m_{\tilde{g}}$, $m_{\tilde{\chi}_{1}^{\pm}}$, $\tau_{\tilde{\chi}_{1}^{\pm}}$) = (1600 GeV, 500 GeV, 0.2 ns) in the high-Emiss region. The expected number of signal events is normalised to 36.1 fb${}^{-1}$. The signal selection efficiencies are also shown in parentheses. The first row shows the number of events after the application of detector and data quality conditions. Requirements below the dashed line are applied to tracks and tracklets.

The event and tracklet generator-level acceptance and selection efficiency for a few electroweak signal models studied in this search. The last column shows the probability ($P$) for a reconstructed tracklet to have $p_{T}$ greater than 100 GeV.

The event and tracklet generator-level acceptance and selection efficiency for a few strong signal models studied in this search. The last column shows the probability ($P$) for a reconstructed tracklet to have $p_{T}$ greater than 100 GeV.

Systematic uncertainties in the signal event yields at $m_{\tilde{\chi}_{1}^{\pm}}$ = 400 GeV for the electroweak channel and at $m_{\tilde{g}}$ = 1600 GeV, $m_{\tilde{\chi}_{1}^{\pm}}$ = 500 GeV for the strong channel. The lifetime of the chargino is not relevant here. The uncertainty in the cross-section of the strong production is large due to the large effect from the PDF uncertainty.

Observed events, expected background for null signal, and expected signal yields for two benchmark models: electroweak channel with ($m_{\tilde{\chi}_{1}^{\pm}}$, $\tau_{\tilde{\chi}_{1}^{\pm}}$) = (400 GeV, 0.2 ns) and strong channel with ($m_{\tilde{g}}$, $m_{\tilde{\chi}_{1}^{\pm}}$, $\tau_{\tilde{\chi}_{1}^{\pm}}$) = (1600 GeV, 500 GeV, 0.2 ns) in the high-Emiss region. Also shown are the probability of a background-only experiment being more signal-like than observed ($p_0$) and the upper limit on the model-independent visible cross-section at 95\% CL. The uncertainty in the total background yield is different from the sum of uncertainties in quadrature due to anti-correlation between different backgrounds.

Effects of systematic uncertainties on the signal exclusion significance at $m_{\tilde{\chi}_{1}^{\pm}}$ = 400 GeV for the electroweak channel and at $m_{\tilde{g}}$ = 1600 GeV, $m_{\tilde{\chi}_{1}^{\pm}}$ = 500 GeV for the strong channel. The lifetime of the chargino is not relevant here. Effects of uncertainties on the fake-tracklet background is smaller in the strong channel analysis because the estimated number of the fake-tracklet background events is small.

Cross-section upper limits for the strong production, presented in unit of fb. Left-upper triangle region is unphysical because the wino mass is larger than the gluino mass.

Observed and predicted mTtot distribution in the b-tag category of the 2tau_h channel. Despite listing this as an exclusive final state (as there must be at least one b-jets), there is no explicit selection on the presence of additional light-flavour jets. Please note that the bin content is divided by the bin width in the paper figure, but not in the HepData table. The last bin includes overflows. The combined prediction for A and H bosons with masses of 300, 500 and 800 GeV and $\tan\beta$ = 10 in the hMSSM scenario are also provided.

Observed and predicted mTtot distribution for the b-inclusive selection in the 1l1tau_h channel. Please note that the bin content is divided by the bin width in the paper figure, but not in the HepData table. In the paper, the first bin is cut off at 60 GeV for aesthetics but contains underflows down to 50 GeV as in the HepData table. The last bin includes overflows. The prediction for a SSM Zprime with masses of 1500, 2000 and 2500 GeV are also provided.

Observed and predicted mTtot distribution for the b-inclusive selection in the 2tau_h channel. Please note that the bin content is divided by the bin width in the paper figure, but not in the HepData table. The last bin includes overflows. The prediction for a SSM Zprime with masses of 1500, 2000 and 2500 GeV are also provided.

Observed and expected 95% CL upper limits on the b-associated Higgs boson production cross section times ditau branching fraction as a function of the boson mass.

Observed and expected 95% CL upper limits on the Drell Yan production cross section times ditau branching fraction as a function of the Zprime boson mass.

Observed and expected 95% CL upper limits on the Higgs boson production cross section times ditau branching fraction as a function of the boson mass and the relative strength of the b-associated production.

Ratio of the 95% CL upper limits on the production cross section times branching fraction for alternate Zprime models with respect to the SSM, both observed and expected are shown.

Acceptance, acceptance times efficiency and b-tag category fraction for a scalar boson produced by gluon-gluon fusion as a function of the scalar boson mass.

Acceptance, acceptance times efficiency and b-tag category fraction for a scalar boson produced by b-associated production as a function of the scalar boson mass.

Acceptance and acceptance times efficiency for a heavy gauge boson produced by Drell Yan as a function of the gauge boson mass.

Two dimensional likelihood scan of the gluon-gluon fusion cross section times braching fraction, $\sigma(gg\phi)\times B(\phi\to\tau\tau)$, vs the b-associated production times branching fraction, $\sigma(bb\phi)\times B(\phi\to\tau\tau)$ for the Higgs boson mass ($m_\phi$) indicated in the table. For each mass, 10000 points are scanned. At each point $\Delta(\mathrm{NLL})$ is calculated, defined as the negative-log-likelihood (NLL) of the conditional fit with $\sigma(gg\phi)$ and $\sigma(bb\phi)$ fixed to their values at the point and with the minimum NLL value at any point subtracted. Vaules are provided for the fit to the observed data and to the expected data, which is the sum of Standard Model contributions not including the SM Higgs boson. The best-fit point and the preferred 68% and 95% boundaries are found at $2\Delta(\mathrm{NLL})$ values of 0.0, 2.30 and 5.90, respectively.

Two dimensional likelihood scan of the gluon-gluon fusion cross section times braching fraction, $\sigma(gg\phi)\times B(\phi\to\tau\tau)$, vs the b-associated production times branching fraction, $\sigma(bb\phi)\times B(\phi\to\tau\tau)$ for the Higgs boson mass ($m_\phi$) indicated in the table. For each mass, 10000 points are scanned. At each point $\Delta(\mathrm{NLL})$ is calculated, defined as the negative-log-likelihood (NLL) of the conditional fit with $\sigma(gg\phi)$ and $\sigma(bb\phi)$ fixed to their values at the point and with the minimum NLL value at any point subtracted. Vaules are provided for the fit to the observed data and to the expected data, which is the sum of Standard Model contributions not including the SM Higgs boson. The best-fit point and the preferred 68% and 95% boundaries are found at $2\Delta(\mathrm{NLL})$ values of 0.0, 2.30 and 5.90, respectively.

Two dimensional likelihood scan of the gluon-gluon fusion cross section times braching fraction, $\sigma(gg\phi)\times B(\phi\to\tau\tau)$, vs the b-associated production times branching fraction, $\sigma(bb\phi)\times B(\phi\to\tau\tau)$ for the Higgs boson mass ($m_\phi$) indicated in the table. For each mass, 10000 points are scanned. At each point $\Delta(\mathrm{NLL})$ is calculated, defined as the negative-log-likelihood (NLL) of the conditional fit with $\sigma(gg\phi)$ and $\sigma(bb\phi)$ fixed to their values at the point and with the minimum NLL value at any point subtracted. Vaules are provided for the fit to the observed data and to the expected data, which is the sum of Standard Model contributions not including the SM Higgs boson. The best-fit point and the preferred 68% and 95% boundaries are found at $2\Delta(\mathrm{NLL})$ values of 0.0, 2.30 and 5.90, respectively.

Two dimensional likelihood scan of the gluon-gluon fusion cross section times braching fraction, $\sigma(gg\phi)\times B(\phi\to\tau\tau)$, vs the b-associated production times branching fraction, $\sigma(bb\phi)\times B(\phi\to\tau\tau)$ for the Higgs boson mass ($m_\phi$) indicated in the table. For each mass, 10000 points are scanned. At each point $\Delta(\mathrm{NLL})$ is calculated, defined as the negative-log-likelihood (NLL) of the conditional fit with $\sigma(gg\phi)$ and $\sigma(bb\phi)$ fixed to their values at the point and with the minimum NLL value at any point subtracted. Vaules are provided for the fit to the observed data and to the expected data, which is the sum of Standard Model contributions not including the SM Higgs boson. The best-fit point and the preferred 68% and 95% boundaries are found at $2\Delta(\mathrm{NLL})$ values of 0.0, 2.30 and 5.90, respectively.

Two dimensional likelihood scan of the gluon-gluon fusion cross section times braching fraction, $\sigma(gg\phi)\times B(\phi\to\tau\tau)$, vs the b-associated production times branching fraction, $\sigma(bb\phi)\times B(\phi\to\tau\tau)$ for the Higgs boson mass ($m_\phi$) indicated in the table. For each mass, 10000 points are scanned. At each point $\Delta(\mathrm{NLL})$ is calculated, defined as the negative-log-likelihood (NLL) of the conditional fit with $\sigma(gg\phi)$ and $\sigma(bb\phi)$ fixed to their values at the point and with the minimum NLL value at any point subtracted. Vaules are provided for the fit to the observed data and to the expected data, which is the sum of Standard Model contributions not including the SM Higgs boson. The best-fit point and the preferred 68% and 95% boundaries are found at $2\Delta(\mathrm{NLL})$ values of 0.0, 2.30 and 5.90, respectively.

Two dimensional likelihood scan of the gluon-gluon fusion cross section times braching fraction, $\sigma(gg\phi)\times B(\phi\to\tau\tau)$, vs the b-associated production times branching fraction, $\sigma(bb\phi)\times B(\phi\to\tau\tau)$ for the Higgs boson mass ($m_\phi$) indicated in the table. For each mass, 10000 points are scanned. At each point $\Delta(\mathrm{NLL})$ is calculated, defined as the negative-log-likelihood (NLL) of the conditional fit with $\sigma(gg\phi)$ and $\sigma(bb\phi)$ fixed to their values at the point and with the minimum NLL value at any point subtracted. Vaules are provided for the fit to the observed data and to the expected data, which is the sum of Standard Model contributions not including the SM Higgs boson. The best-fit point and the preferred 68% and 95% boundaries are found at $2\Delta(\mathrm{NLL})$ values of 0.0, 2.30 and 5.90, respectively.

Two dimensional likelihood scan of the gluon-gluon fusion cross section times braching fraction, $\sigma(gg\phi)\times B(\phi\to\tau\tau)$, vs the b-associated production times branching fraction, $\sigma(bb\phi)\times B(\phi\to\tau\tau)$ for the Higgs boson mass ($m_\phi$) indicated in the table. For each mass, 10000 points are scanned. At each point $\Delta(\mathrm{NLL})$ is calculated, defined as the negative-log-likelihood (NLL) of the conditional fit with $\sigma(gg\phi)$ and $\sigma(bb\phi)$ fixed to their values at the point and with the minimum NLL value at any point subtracted. Vaules are provided for the fit to the observed data and to the expected data, which is the sum of Standard Model contributions not including the SM Higgs boson. The best-fit point and the preferred 68% and 95% boundaries are found at $2\Delta(\mathrm{NLL})$ values of 0.0, 2.30 and 5.90, respectively.

Two dimensional likelihood scan of the gluon-gluon fusion cross section times braching fraction, $\sigma(gg\phi)\times B(\phi\to\tau\tau)$, vs the b-associated production times branching fraction, $\sigma(bb\phi)\times B(\phi\to\tau\tau)$ for the Higgs boson mass ($m_\phi$) indicated in the table. For each mass, 10000 points are scanned. At each point $\Delta(\mathrm{NLL})$ is calculated, defined as the negative-log-likelihood (NLL) of the conditional fit with $\sigma(gg\phi)$ and $\sigma(bb\phi)$ fixed to their values at the point and with the minimum NLL value at any point subtracted. Vaules are provided for the fit to the observed data and to the expected data, which is the sum of Standard Model contributions not including the SM Higgs boson. The best-fit point and the preferred 68% and 95% boundaries are found at $2\Delta(\mathrm{NLL})$ values of 0.0, 2.30 and 5.90, respectively.

Two dimensional likelihood scan of the gluon-gluon fusion cross section times braching fraction, $\sigma(gg\phi)\times B(\phi\to\tau\tau)$, vs the b-associated production times branching fraction, $\sigma(bb\phi)\times B(\phi\to\tau\tau)$ for the Higgs boson mass ($m_\phi$) indicated in the table. For each mass, 10000 points are scanned. At each point $\Delta(\mathrm{NLL})$ is calculated, defined as the negative-log-likelihood (NLL) of the conditional fit with $\sigma(gg\phi)$ and $\sigma(bb\phi)$ fixed to their values at the point and with the minimum NLL value at any point subtracted. Vaules are provided for the fit to the observed data and to the expected data, which is the sum of Standard Model contributions not including the SM Higgs boson. The best-fit point and the preferred 68% and 95% boundaries are found at $2\Delta(\mathrm{NLL})$ values of 0.0, 2.30 and 5.90, respectively.

Two dimensional likelihood scan of the gluon-gluon fusion cross section times braching fraction, $\sigma(gg\phi)\times B(\phi\to\tau\tau)$, vs the b-associated production times branching fraction, $\sigma(bb\phi)\times B(\phi\to\tau\tau)$ for the Higgs boson mass ($m_\phi$) indicated in the table. For each mass, 10000 points are scanned. At each point $\Delta(\mathrm{NLL})$ is calculated, defined as the negative-log-likelihood (NLL) of the conditional fit with $\sigma(gg\phi)$ and $\sigma(bb\phi)$ fixed to their values at the point and with the minimum NLL value at any point subtracted. Vaules are provided for the fit to the observed data and to the expected data, which is the sum of Standard Model contributions not including the SM Higgs boson. The best-fit point and the preferred 68% and 95% boundaries are found at $2\Delta(\mathrm{NLL})$ values of 0.0, 2.30 and 5.90, respectively.

Two dimensional likelihood scan of the gluon-gluon fusion cross section times braching fraction, $\sigma(gg\phi)\times B(\phi\to\tau\tau)$, vs the b-associated production times branching fraction, $\sigma(bb\phi)\times B(\phi\to\tau\tau)$ for the Higgs boson mass ($m_\phi$) indicated in the table. For each mass, 10000 points are scanned. At each point $\Delta(\mathrm{NLL})$ is calculated, defined as the negative-log-likelihood (NLL) of the conditional fit with $\sigma(gg\phi)$ and $\sigma(bb\phi)$ fixed to their values at the point and with the minimum NLL value at any point subtracted. Vaules are provided for the fit to the observed data and to the expected data, which is the sum of Standard Model contributions not including the SM Higgs boson. The best-fit point and the preferred 68% and 95% boundaries are found at $2\Delta(\mathrm{NLL})$ values of 0.0, 2.30 and 5.90, respectively.

Two dimensional likelihood scan of the gluon-gluon fusion cross section times braching fraction, $\sigma(gg\phi)\times B(\phi\to\tau\tau)$, vs the b-associated production times branching fraction, $\sigma(bb\phi)\times B(\phi\to\tau\tau)$ for the Higgs boson mass ($m_\phi$) indicated in the table. For each mass, 10000 points are scanned. At each point $\Delta(\mathrm{NLL})$ is calculated, defined as the negative-log-likelihood (NLL) of the conditional fit with $\sigma(gg\phi)$ and $\sigma(bb\phi)$ fixed to their values at the point and with the minimum NLL value at any point subtracted. Vaules are provided for the fit to the observed data and to the expected data, which is the sum of Standard Model contributions not including the SM Higgs boson. The best-fit point and the preferred 68% and 95% boundaries are found at $2\Delta(\mathrm{NLL})$ values of 0.0, 2.30 and 5.90, respectively.

Two dimensional likelihood scan of the gluon-gluon fusion cross section times braching fraction, $\sigma(gg\phi)\times B(\phi\to\tau\tau)$, vs the b-associated production times branching fraction, $\sigma(bb\phi)\times B(\phi\to\tau\tau)$ for the Higgs boson mass ($m_\phi$) indicated in the table. For each mass, 10000 points are scanned. At each point $\Delta(\mathrm{NLL})$ is calculated, defined as the negative-log-likelihood (NLL) of the conditional fit with $\sigma(gg\phi)$ and $\sigma(bb\phi)$ fixed to their values at the point and with the minimum NLL value at any point subtracted. Vaules are provided for the fit to the observed data and to the expected data, which is the sum of Standard Model contributions not including the SM Higgs boson. The best-fit point and the preferred 68% and 95% boundaries are found at $2\Delta(\mathrm{NLL})$ values of 0.0, 2.30 and 5.90, respectively.

Two dimensional likelihood scan of the gluon-gluon fusion cross section times braching fraction, $\sigma(gg\phi)\times B(\phi\to\tau\tau)$, vs the b-associated production times branching fraction, $\sigma(bb\phi)\times B(\phi\to\tau\tau)$ for the Higgs boson mass ($m_\phi$) indicated in the table. For each mass, 10000 points are scanned. At each point $\Delta(\mathrm{NLL})$ is calculated, defined as the negative-log-likelihood (NLL) of the conditional fit with $\sigma(gg\phi)$ and $\sigma(bb\phi)$ fixed to their values at the point and with the minimum NLL value at any point subtracted. Vaules are provided for the fit to the observed data and to the expected data, which is the sum of Standard Model contributions not including the SM Higgs boson. The best-fit point and the preferred 68% and 95% boundaries are found at $2\Delta(\mathrm{NLL})$ values of 0.0, 2.30 and 5.90, respectively.

Two dimensional likelihood scan of the gluon-gluon fusion cross section times braching fraction, $\sigma(gg\phi)\times B(\phi\to\tau\tau)$, vs the b-associated production times branching fraction, $\sigma(bb\phi)\times B(\phi\to\tau\tau)$ for the Higgs boson mass ($m_\phi$) indicated in the table. For each mass, 10000 points are scanned. At each point $\Delta(\mathrm{NLL})$ is calculated, defined as the negative-log-likelihood (NLL) of the conditional fit with $\sigma(gg\phi)$ and $\sigma(bb\phi)$ fixed to their values at the point and with the minimum NLL value at any point subtracted. Vaules are provided for the fit to the observed data and to the expected data, which is the sum of Standard Model contributions not including the SM Higgs boson. The best-fit point and the preferred 68% and 95% boundaries are found at $2\Delta(\mathrm{NLL})$ values of 0.0, 2.30 and 5.90, respectively.

Observed and expected 95% CL upper limits on the gluon-gluon fusion Higgs boson production cross section times ditau branching fraction as a function of the boson mass.

Distribution of $R_\text{ISR}$ for SRC1-5 after the likelihood fit. The stacked histograms show the SM expectation and the hatched uncertainty band around the SM expectation shows the MC statistical and detector-related systematic uncertainties. A representative signal point is shown for each distribution.

Distribution of $m_\text{T}^{b,\text{max}}$ for SRD-high after the likelihood fit. The stacked histograms show the SM expectation and the hatched uncertainty band around the SM expectation shows the MC statistical and detector-related systematic uncertainties. A representative signal point is shown for each distribution.

Distribution of $H_\text{T}$ for SRE after the likelihood fit. The stacked histograms show the SM expectation and the hatched uncertainty band around the SM expectation shows the MC statistical and detector-related systematic uncertainties. A representative signal point is shown for each distribution.

Expected (blue solid line) exclusion limits at 95% CL as a function of stop and LSP masses in the scenario where both top squarks decay via $\tilde{t}\to t^{(*)} \tilde{\chi^{0}_{1}}$.

Observed (red solid line) exclusion limits at 95% CL as a function of stop and LSP masses in the scenario where both top squarks decay via $\tilde{t}\to t^{(*)} \tilde{\chi^{0}_{1}}$.

Expected exclusion limits at 95% CL as a function of $\tilde{t}$ and $\tilde{\chi^{0}_{1}}$ masses and branching fraction to $\tilde{t}\to t^{(*)} \tilde{\chi^{0}_{1}}$ in the Natural SUSY-inspired mixed grid scenario where $m_{\tilde{\chi^{\pm}_{1}}}=m_{\tilde{\chi^{0}_{1}}}$+1 GeV.

Observed exclusion limits at 95% CL as a function of $\tilde{t}$ and $\tilde{\chi^{0}_{1}}$ masses and branching fraction to $\tilde{t}\to t^{(*)} \tilde{\chi^{0}_{1}}$ in the Natural SUSY-inspired mixed grid scenario where $m_{\tilde{\chi^{\pm}_{1}}}=m_{\tilde{\chi^{0}_{1}}}$+1 GeV.

Expected exclusion limits at 95% CL as a function of $\tilde{t}$ and $\tilde{\chi^{0}_{1}}$ masses and branching fraction to $\tilde{t}\to t^{(*)} \tilde{\chi^{0}_{1}}$ in the Natural SUSY-inspired mixed grid scenario where $m_{\tilde{\chi^{\pm}_{1}}}=m_{\tilde{\chi^{0}_{1}}}$+1 GeV.

Observed exclusion limits at 95% CL as a function of $\tilde{t}$ and $\tilde{\chi^{0}_{1}}$ masses and branching fraction to $\tilde{t}\to t^{(*)} \tilde{\chi^{0}_{1}}$ in the Natural SUSY-inspired mixed grid scenario where $m_{\tilde{\chi^{\pm}_{1}}}=m_{\tilde{\chi^{0}_{1}}}$+1 GeV.

Expected exclusion limits at 95% CL as a function of $\tilde{t}$ and $\tilde{\chi^{0}_{1}}$ masses and branching fraction to $\tilde{t}\to t^{(*)} \tilde{\chi^{0}_{1}}$ in the Natural SUSY-inspired mixed grid scenario where $m_{\tilde{\chi^{\pm}_{1}}}=m_{\tilde{\chi^{0}_{1}}}$+1 GeV.

Observed exclusion limits at 95% CL as a function of $\tilde{t}$ and $\tilde{\chi^{0}_{1}}$ masses and branching fraction to $\tilde{t}\to t^{(*)} \tilde{\chi^{0}_{1}}$ in the Natural SUSY-inspired mixed grid scenario where $m_{\tilde{\chi^{\pm}_{1}}}=m_{\tilde{\chi^{0}_{1}}}$+1 GeV.

Expected exclusion limits at 95% CL as a function of $\tilde{t}$ and $\tilde{\chi^{0}_{1}}$ masses and branching fraction to $\tilde{t}\to t^{(*)} \tilde{\chi^{0}_{1}}$ in the Natural SUSY-inspired mixed grid scenario where $m_{\tilde{\chi^{\pm}_{1}}}=m_{\tilde{\chi^{0}_{1}}}$+1 GeV.

Observed exclusion limits at 95% CL as a function of $\tilde{t}$ and $\tilde{\chi^{0}_{1}}$ masses and branching fraction to $\tilde{t}\to t^{(*)} \tilde{\chi^{0}_{1}}$ in the Natural SUSY-inspired mixed grid scenario where $m_{\tilde{\chi^{\pm}_{1}}}=m_{\tilde{\chi^{0}_{1}}}$+1 GeV.

Expected exclusion limits at 95% CL as a function of $m_{\tilde{t}}$ and $m_{\tilde{\chi^{0}_{1}}}$ for the pMSSM-inspired non-asymptotic Higgsino simplified model for a small tan$\beta$ assumption.

Observed exclusion limits at 95% CL as a function of $m_{\tilde{t}}$ and $m_{\tilde{\chi^{0}_{1}}}$ for the pMSSM-inspired non-asymptotic Higgsino simplified model for a small tan$\beta$ assumption.

Expected exclusion limits at 95% CL as a function of $m_{\tilde{t}}$ and $m_{\tilde{\chi^{0}_{1}}}$ for the pMSSM-inspired non-asymptotic Higgsino simplified model for a large tan$\beta$ assumption.

Observed exclusion limits at 95% CL as a function of $m_{\tilde{t}}$ and $m_{\tilde{\chi^{0}_{1}}}$ for the pMSSM-inspired non-asymptotic Higgsino simplified model for a large tan$\beta$ assumption.

Expected exclusion limits at 95% CL as a function of $m_{\tilde{t}}$ and $m_{\tilde{\chi^{0}_{1}}}$ for the pMSSM-inspired non-asymptotic Higgsino simplified model for a small right-handed top-squark mass parameter assumption.

Observed exclusion limits at 95% CL as a function of $m_{\tilde{t}}$ and $m_{\tilde{\chi^{0}_{1}}}$ for the pMSSM-inspired non-asymptotic Higgsino simplified model for a small right-handed top-squark mass parameter assumption.

Expected exclusion limits at 95% CL as a function of $\tilde{t}$ and $\tilde{\chi^{0}_{1}}$ masses for the Wino NLSP pMSSM model for a negative value of $\mu$.

Observed exclusion limits at 95% CL as a function of $\tilde{t}$ and $\tilde{\chi^{0}_{1}}$ masses for the Wino NLSP pMSSM model for a negative value of $\mu$.

Expected exclusion limits at 95% CL as a function of $\tilde{t}$ and $\tilde{\chi^{0}_{1}}$ masses for the Wino NLSP pMSSM model for a positive value of $\mu$.

Observed exclusion limits at 95% CL as a function of $\tilde{t}$ and $\tilde{\chi^{0}_{1}}$ masses for the Wino NLSP pMSSM model for a positive value of $\mu$.

Expected exclusion limits at 95% CL as a function of $\tilde{t}$ and $\tilde{\chi^{0}_{1}}$ masses for for the left-handed top-squark mass parameter scan in the well-tempered pMSSM model.

Observed exclusion limits at 95% CL as a function of $\tilde{t}$ and $\tilde{\chi^{0}_{1}}$ masses for for the left-handed top-squark mass parameter scan in the well-tempered pMSSM model.

Expected exclusion limits at 95% CL as a function of $\tilde{t}$ and $\tilde{\chi^{0}_{1}}$ masses for for the right-handed top-squark mass parameter scan in the well-tempered pMSSM model.

Observed exclusion limits at 95% CL as a function of $\tilde{t}$ and $\tilde{\chi^{0}_{1}}$ masses for for the right-handed top-squark mass parameter scan in the well-tempered pMSSM model.

Expected exclusion limits at 95% CL exclusion as a function of $\tilde{g}$ and $\tilde{t}$ masses in the scenario where both gluinos decay via $\tilde{g}\to t \tilde{t}\to t\tilde{\chi^{0}_{1}}+$soft and $\Delta m(\tilde{t},\tilde{\chi^{0}_{1}})=5$ GeV.

Observed exclusion limits at 95% CL exclusion as a function of $\tilde{g}$ and $\tilde{t}$ masses in the scenario where both gluinos decay via $\tilde{g}\to t \tilde{t}\to t\tilde{\chi^{0}_{1}}+$soft and $\Delta m(\tilde{t},\tilde{\chi^{0}_{1}})=5$ GeV.

Results of the exclusion fits for top squark pair production in the case where both top squarks decay via $\tilde{t}\to t^{(*)} \tilde{\chi^{0}_{1}}$ from the combination of SRA, SRB, SRC, SRD and SRE, based on the best expected $CL_s$. The numbers centered on the grid points indicate which of the signal regions gave the best expected $CL_s$ (with 1, 2, 3, 4, 5, 6 corresponding to SRA, SRB, SRC, SRD-low,SR D-high, SRE respectively).

Results of the exclusion fits in the grid with two stop decay channels: $\tilde{t}\to t^{(*)} \tilde{\chi^{0}_{1}}$ and $\tilde{t}\to b \tilde{\chi^{\pm}_{1}} \to b W^{(*)} \tilde{\chi^{0}_{1}}$, with $m(\tilde{\chi^{\pm}_{1}})-m(\tilde{\chi^{0}_{1}}) = 1$ GeV. The results are shown as a function of the branching ratio to $\tilde{t}\to t^{(*)} \tilde{\chi^{0}_{1}}$: 0%. The results are based on taking the signal region with the best expected $CL_s$, using SRA, SRB, SRC, SRD-low and SRD-high, The numbers centered on the grid points indicate the signal region used (1: SRA, 2: SRB, 3: SRD-low, 4: SRD-high, 5: SRC).

Results of the exclusion fits in the grid with two stop decay channels: $\tilde{t}\to t^{(*)} \tilde{\chi^{0}_{1}}$ and $\tilde{t}\to b \tilde{\chi^{\pm}_{1}} \to b W^{(*)} \tilde{\chi^{0}_{1}}$, with $m(\tilde{\chi^{\pm}_{1}})-m(\tilde{\chi^{0}_{1}}) = 1$ GeV. The results are shown as a function of the branching ratio to $\tilde{t}\to t^{(*)} \tilde{\chi^{0}_{1}}$: 25%. The results are based on taking the signal region with the best expected $CL_s$, using SRA, SRB, SRC, SRD-low and SRD-high, The numbers centered on the grid points indicate the signal region used (1: SRA, 2: SRB, 3: SRD-low, 4: SRD-high, 5: SRC).

Results of the exclusion fits in the grid with two stop decay channels: $\tilde{t}\to t^{(*)} \tilde{\chi^{0}_{1}}$ and $\tilde{t}\to b \tilde{\chi^{\pm}_{1}} \to b W^{(*)} \tilde{\chi^{0}_{1}}$, with $m(\tilde{\chi^{\pm}_{1}})-m(\tilde{\chi^{0}_{1}}) = 1$ GeV. The results are shown as a function of the branching ratio to $\tilde{t}\to t^{(*)} \tilde{\chi^{0}_{1}}$: 50%. The results are based on taking the signal region with the best expected $CL_s$, using SRA, SRB, SRC, SRD-low and SRD-high, The numbers centered on the grid points indicate the signal region used (1: SRA, 2: SRB, 3: SRD-low, 4: SRD-high, 5: SRC).

Results of the exclusion fits in the grid with two stop decay channels: $\tilde{t}\to t^{(*)} \tilde{\chi^{0}_{1}}$ and $\tilde{t}\to b \tilde{\chi^{\pm}_{1}} \to b W^{(*)} \tilde{\chi^{0}_{1}}$, with $m(\tilde{\chi^{\pm}_{1}})-m(\tilde{\chi^{0}_{1}}) = 1$ GeV. The results are shown as a function of the branching ratio to $\tilde{t}\to t^{(*)} \tilde{\chi^{0}_{1}}$: 0 (top left), 25% (top right), 50% (middle left), 75% middle right) and 100% (bottom). The results are based on taking the signal region with the best expected $CL_s$, using SRA, SRB, SRC, SRD-low and SRD-high, The numbers centered on the grid points indicate the signal region used (1: SRA, 2: SRB, 3: SRD-low, 4: SRD-high, 5: SRC).

Results of the exclusion fits in the wino NLSP grid for negative values of $\mu$. The results are based on taking the signal region with the best expected $CL_s$, using SRA, SRB, SRD-low and SRD-high, where SRA and SRB are the statistical combinations of their respective regions. The numbers centered on the grid points indicate the signal region used (1: SRA, 2: SRB, 4: SRD-low, 5: SRD-high).

Results of the exclusion fits in the wino NLSP grid for positive values of $\mu$. The results are based on taking the signal region with the best expected $CL_s$, using SRA, SRB, SRD-low and SRD-high, where SRA and SRB are the statistical combinations of their respective regions. The numbers centered on the grid points indicate the signal region used (1: SRA, 2: SRB, 4: SRD-low, 5: SRD-high).

Results of the exclusion fits in the well-tempered neutralino grid for the $m_{q3L}$ scenario. The results are based on taking the signal region with the best expected $CL_s$, using SRA, SRB, SRD-low and SRD-high, where SRA and SRB are the statistical combinations of their respective regions. The numbers centered on the grid points indicate the signal region used (1: SRA, 2: SRB, 4: SRD-low, 5: SRD-high).

Results of the exclusion fits in the well-tempered neutralino grid for the $m_{tR}$ scenario. The results are based on taking the signal region with the best expected $CL_s$, using SRA, SRB, SRD-low and SRD-high, where SRA and SRB are the statistical combinations of their respective regions. The numbers centered on the grid points indicate the signal region used (1: SRA, 2: SRB, 4: SRD-low, 5: SRD-high).

Results of the exclusion fits in the non-asymptotic higgsino grid with $m(\tilde{\chi^{\pm}_{1}}) - m(\tilde{\chi^{0}_{1}}) = 5$ GeV. A scenarios with large tan$\beta$ (top left) is shown. The results are based on taking the signal region with the best expected $CL_s$, using SRA, SRB, SRD-low and SRD-high, where SRA and SRB are the statistical combinations of their respective regions. The numbers centered on the grid points indicate the signal region used (1: SRA, 2: SRB, 3: SRD-low, 4: SRD-high).

Results of the exclusion fits in the non-asymptotic higgsino grid with $m(\tilde{\chi^{\pm}_{1}}) - m(\tilde{\chi^{0}_{1}}) = 5$ GeV. A scenarios with small tan$\beta$ are shown. The results are based on taking the signal region with the best expected $CL_s$, using SRA, SRB, SRD-low and SRD-high, where SRA and SRB are the statistical combinations of their respective regions. The numbers centered on the grid points indicate the signal region used (1: SRA, 2: SRB, 3: SRD-low, 4: SRD-high).

Results of the exclusion fits in the non-asymptotic higgsino grid with $m(\tilde{\chi^{\pm}_{1}}) - m(\tilde{\chi^{0}_{1}}) = 5$ GeV. A scenario with a mostly right-handed top squark partner is shown. The results are based on taking the signal region with the best expected $CL_s$, using SRA, SRB, SRD-low and SRD-high, where SRA and SRB are the statistical combinations of their respective regions. The numbers centered on the grid points indicate the signal region used (1: SRA, 2: SRB, 3: SRD-low, 4: SRD-high).

Acceptance for SRA-TT for top squark pair production in the case where both top squarks decay via $\tilde{t}\to t^{(*)} \tilde{\chi^{0}_{1}}$.

Acceptance for SRB-TT for top squark pair production in the case where both top squarks decay via $\tilde{t}\to t^{(*)} \tilde{\chi^{0}_{1}}$.

Acceptance for SRE for top squark pair production in the case where both top squarks decay via $\tilde{t}\to t^{(*)} \tilde{\chi^{0}_{1}}$.

Acceptance for SRA-TW for top squark pair production in the case where both top squarks decay via $\tilde{t}\to t^{(*)} \tilde{\chi^{0}_{1}}$.

Acceptance for SRB-TW for top squark pair production in the case where both top squarks decay via $\tilde{t}\to t^{(*)} \tilde{\chi^{0}_{1}}$.

Acceptance for SRA-T0 for top squark pair production in the case where both top squarks decay via $\tilde{t}\to t^{(*)} \tilde{\chi^{0}_{1}}$.

Acceptance for SRB-T0 for top squark pair production in the case where both top squarks decay via $\tilde{t}\to t^{(*)} \tilde{\chi^{0}_{1}}$.

Acceptance for SRC1 for top squark pair production in the case where both top squarks decay via $\tilde{t}\to t^{(*)} \tilde{\chi^{0}_{1}}$.

Acceptance for SRC2 for top squark pair production in the case where both top squarks decay via $\tilde{t}\to t^{(*)} \tilde{\chi^{0}_{1}}$.

Acceptance for SRC3 for top squark pair production in the case where both top squarks decay via $\tilde{t}\to t^{(*)} \tilde{\chi^{0}_{1}}$.

Acceptance for SRC4 for top squark pair production in the case where both top squarks decay via $\tilde{t}\to t^{(*)} \tilde{\chi^{0}_{1}}$.

Acceptance for SRC5 for top squark pair production in the case where both top squarks decay via $\tilde{t}\to t^{(*)} \tilde{\chi^{0}_{1}}$.

Acceptance for SRD-low for the natural SUSY-inspired mixed grid in which two decay modes are considered, the $\tilde{t}\to t^{(*)} \tilde{\chi^{0}_{1}}$ and $\tilde{t}\to b \tilde{\chi^{\pm}_{1}} \to b W^{(*)} \tilde{\chi^{0}_{1}}$, with $m(\tilde{\chi^{\pm}_{1}})-m(\tilde{\chi^{0}_{1}}) = 1$ GeV.

Acceptance for SRD-high for the natural SUSY-inspired mixed grid in which two decay modes are considered, the $\tilde{t}\to t^{(*)} \tilde{\chi^{0}_{1}}$ and $\tilde{t}\to b \tilde{\chi^{\pm}_{1}} \to b W^{(*)} \tilde{\chi^{0}_{1}}$, with $m(\tilde{\chi^{\pm}_{1}})-m(\tilde{\chi^{0}_{1}}) = 1$ GeV.

Acceptance for SRE for gluino pair production in the case where both gluinos decay via $\tilde{g}\to t \tilde{t} \to t\tilde{\chi^{0}_{1}}+$soft and $\Delta m(\tilde{t},\tilde{\chi^{0}_{1}})=5$ GeV.

Efficiencies for SRA-TT for top squark pair production in the case where both top squarks decay via $\tilde{t}\to t^{(*)} \tilde{\chi^{0}_{1}}$.

Efficiencies for SRB-TT for top squark pair production in the case where both top squarks decay via $\tilde{t}\to t^{(*)} \tilde{\chi^{0}_{1}}$.

Efficiencies for SRE for top squark pair production in the case where both top squarks decay via $\tilde{t}\to t^{(*)} \tilde{\chi^{0}_{1}}$.

Efficiencies for SRA-TW for top squark pair production in the case where both top squarks decay via $\tilde{t}\to t^{(*)} \tilde{\chi^{0}_{1}}$.

Efficiencies for SRB-TW for top squark pair production in the case where both top squarks decay via $\tilde{t}\to t^{(*)} \tilde{\chi^{0}_{1}}$.

Efficiencies for SRA-T0 for top squark pair production in the case where both top squarks decay via $\tilde{t}\to t^{(*)} \tilde{\chi^{0}_{1}}$.

Efficiencies for SRB-T0 for top squark pair production in the case where both top squarks decay via $\tilde{t}\to t^{(*)} \tilde{\chi^{0}_{1}}$.

Efficiencies for SRC1 for top squark pair production in the case where both top squarks decay via $\tilde{t}\to t^{(*)} \tilde{\chi^{0}_{1}}$.

Efficiencies for SRC2 for top squark pair production in the case where both top squarks decay via $\tilde{t}\to t^{(*)} \tilde{\chi^{0}_{1}}$.

Efficiencies for SRC3 for top squark pair production in the case where both top squarks decay via $\tilde{t}\to t^{(*)} \tilde{\chi^{0}_{1}}$.

Efficiencies for SRC4 for top squark pair production in the case where both top squarks decay via $\tilde{t}\to t^{(*)} \tilde{\chi^{0}_{1}}$.

Efficiencies for SRC5 for top squark pair production in the case where both top squarks decay via $\tilde{t}\to t^{(*)} \tilde{\chi^{0}_{1}}$.

Efficiencies for SRD-low for the natural SUSY-inspired mixed grid in which two decay modes are considered, the $\tilde{t}\to t^{(*)} \tilde{\chi^{0}_{1}}$ and $\tilde{t}\to b \tilde{\chi^{\pm}_{1}} \to b W^{(*)} \tilde{\chi^{0}_{1}}$, with $m(\tilde{\chi^{\pm}_{1}})-m(\tilde{\chi^{0}_{1}}) = 1$ GeV.

Efficiencies for SRD-high for the natural SUSY-inspired mixed grid in which two decay modes are considered, the $\tilde{t}\to t^{(*)} \tilde{\chi^{0}_{1}}$ and $\tilde{t}\to b \tilde{\chi^{\pm}_{1}} \to b W^{(*)} \tilde{\chi^{0}_{1}}$, with $m(\tilde{\chi^{\pm}_{1}})-m(\tilde{\chi^{0}_{1}}) = 1$ GeV.

Efficiencies for SRE for gluino pair production in the case where both gluinos decay via $\tilde{g}\to t \tilde{t} \to t\tilde{\chi^{0}_{1}}+$soft and $\Delta m(\tilde{t},\tilde{\chi^{0}_{1}})=5$ GeV.

Upper limit cross-section, in femtobarn, for the $\tilde{t}\to t^{(*)} \tilde{\chi^{0}_{1}}$ grid.

Upper limit cross-section, in femtobarn, for the $\tilde{g}\to t \tilde{t} \to t\tilde{\chi^{0}_{1}}+$soft grid.

Cutflow for SRA for a signal model with top squark pair production in the case where both top squarks decay via $\tilde{t}_1\to t^{(*)} \tilde\chi^0_1$ with $m(\tilde{t}_1,\tilde\chi^0_1)=$ (800,1) GeV. An integrated luminosity of 36.1 $\mathrm{fb}^{-1}$ is assumed when calculating the weighted yields. For the derivation skim at least one of the following four criteria is required: $H_{\mathrm{T}}$ $>$ 150 GeV; at least one loose electron with $p_{\mathrm{T}}$ $>$ 100 GeV or at least two loose electrons with $p_{\mathrm{T}}$ $>$ 20 GeV; at least one muon with $p_{\mathrm{T}}$ $>$ 100 GeV or at least two muons with $p_{\mathrm{T}}$ $>$ 20 GeV; or at least one photon with $p_{\mathrm{T}}$ $>$ 100 GeV or at least two photons with $p_{\mathrm{T}}$ $>$ 50 GeV.

Cutflow SRB for a signal model with top squark pair production in the case where both top squarks decay via $\tilde{t}_1\to t^{(*)} \tilde\chi^0_1$ with $m(\tilde{t}_1,\tilde\chi^0_1)=$ (600,300) GeV. An integrated luminosity of 36.1 $\mathrm{fb}^{-1}$ is assumed when calculating the weighted yields. For the derivation skim at least one of the following four criteria is required: $H_{\mathrm{T}}$ $>$ 150 GeV; at least one loose electron with $p_{\mathrm{T}}$ $>$ 100 GeV or at least two loose electrons with $p_{\mathrm{T}}$ $>$ 20 GeV; at least one muon with $p_{\mathrm{T}}$ $>$ 100 GeV or at least two muons with $p_{\mathrm{T}}$ $>$ 20 GeV; or at least one photon with $p_{\mathrm{T}}$ $>$ 100 GeV or at least two photons with $p_{\mathrm{T}}$ $>$ 50 GeV.

Cutflow for SRC1 and SRC2 for a signal model with top squark pair production in the case where both top squarks decay via $\tilde{t}_1\to t^{(*)} \tilde\chi^0_1$ with $m(\tilde{t}_1,\tilde\chi^0_1)=$ (250,77) GeV. An integrated luminosity of 36.1 $\mathrm{fb}^{-1}$ is assumed when calculating the weighted yields. For the derivation skim at least one of the following four criteria is required: $H_{\mathrm{T}}$ $>$ 150 GeV; at least one loose electron with $p_{\mathrm{T}}$ $>$ 100 GeV or at least two loose electrons with $p_{\mathrm{T}}$ $>$ 20 GeV; at least one muon with $p_{\mathrm{T}}$ $>$ 100 GeV or at least two muons with $p_{\mathrm{T}}$ $>$ 20 GeV; or at least one photon with $p_{\mathrm{T}}$ $>$ 100 GeV or at least two photons with $p_{\mathrm{T}}$ $>$ 50 GeV.

Cutflow for SRC3, SRC4, and SRC5 for a signal model with top squark pair production in the case where both top squarks decay via $\tilde{t}_1\to t^{(*)} \tilde\chi^0_1$ with $m(\tilde{t}_1,\tilde\chi^0_1)=$ (500,327) GeV. An integrated luminosity of 36.1 $\mathrm{fb}^{-1}$ is assumed when calculating the weighted yields. For the derivation skim at least one of the following four criteria is required: $H_{\mathrm{T}}$ $>$ 150 GeV; at least one loose electron with $p_{\mathrm{T}}$ $>$ 100 GeV or at least two loose electrons with $p_{\mathrm{T}}$ $>$ 20 GeV; at least one muon with $p_{\mathrm{T}}$ $>$ 100 GeV or at least two muons with $p_{\mathrm{T}}$ $>$ 20 GeV; or at least one photon with $p_{\mathrm{T}}$ $>$ 100 GeV or at least two photons with $p_{\mathrm{T}}$ $>$ 50 GeV.

Cutflow for SRD-high for a signal model with bottom squark pair production in the case where both bottom squarks decay via $b\tilde\chi^{\pm}_1\to bW^{(*)} \tilde\chi^0_1$, with $m(\tilde\chi^{\pm}_1)-m(\tilde\chi^0_1) = 1$ GeV with $m(\tilde{b}_1,\tilde\chi^{\pm}_1)=$ (750,200) GeV. An integrated luminosity of 36.1 $\mathrm{fb}^{-1}$ is assumed when calculating the weighted yields. For the derivation skim at least one of the following four criteria is required: $H_{\mathrm{T}}$ $>$ 150 GeV; at least one loose electron with $p_{\mathrm{T}}$ $>$ 100 GeV or at least two loose electrons with $p_{\mathrm{T}}$ $>$ 20 GeV; at least one muon with $p_{\mathrm{T}}$ $>$ 100 GeV or at least two muons with $p_{\mathrm{T}}$ $>$ 20 GeV; or at least one photon with $p_{\mathrm{T}}$ $>$ 100 GeV or at least two photons with $p_{\mathrm{T}}$ $>$ 50 GeV.

Cutflow for SRD-low for a signal model with bottom squark pair production in the case where both bottom squarks decay via $b\tilde\chi^{\pm}_1\to bW^{(*)} \tilde\chi^0_1$, with $m(\tilde\chi^{\pm}_1)-m(\tilde\chi^0_1) = 1$ GeV witht $m(\tilde{b}_1,\tilde\chi^{\pm}_1)=$ (400,200) GeV. An integrated luminosity of 36.1 $\mathrm{fb}^{-1}$ is assumed when calculating the weighted yields. For the derivation skim at least one of the following four criteria is required: $H_{\mathrm{T}}$ $>$ 150 GeV; at least one loose electron with $p_{\mathrm{T}}$ $>$ 100 GeV or at least two loose electrons with $p_{\mathrm{T}}$ $>$ 20 GeV; at least one muon with $p_{\mathrm{T}}$ $>$ 100 GeV or at least two muons with $p_{\mathrm{T}}$ $>$ 20 GeV; or at least one photon with $p_{\mathrm{T}}$ $>$ 100 GeV or at least two photons with $p_{\mathrm{T}}$ $>$ 50 GeV.

Cutflow for SRE for a signal model with gluino pair production in the case where both gluinos decay via $\tilde{g}\to t\tilde{t}_1\to t\tilde\chi^0_1+$soft and $\Delta m(\tilde{t}_1, \tilde\chi^0_1)=5$ GeV with $m(\tilde{g},\tilde{t}_1)=$ (1700,400) GeV. An integrated luminosity of 36.1 $\mathrm{fb}^{-1}$ is assumed when calculating the weighted yields. For the derivation skim at least one of the following four criteria is required: $H_{\mathrm{T}}$ $>$ 150 GeV; at least one loose electron with $p_{\mathrm{T}}$ $>$ 100 GeV or at least two loose electrons with $p_{\mathrm{T}}$ $>$ 20 GeV; at least one muon with $p_{\mathrm{T}}$ $>$ 100 GeV or at least two muons with $p_{\mathrm{T}}$ $>$ 20 GeV; or at least one photon with $p_{\mathrm{T}}$ $>$ 100 GeV or at least two photons with $p_{\mathrm{T}}$ $>$ 50 GeV.

Upper limits at the 95% CL on the cross section for SSM W' production and decay to the muon+neutrino channel as a function of the W' pole mass.

Combined (electron and muon channels) upper limits at the 95% CL on the cross section for SSM W' production and decay to a single lepton generation as a function of the W' pole mass.

Product of acceptance and efficiency for the electron and muon selections as a function of the W' pole mass.

Event yields for the different processes estimated with the fit to the $O_\mathrm{NN}$ distribution compared to the numbers of observed events. Only the statistical uncertainties are quoted. The $Z,VV+\mathrm{jets}$ contributions and the multijet background are fixed in the fit; therefore no uncertainty is quoted for these processes.

Detailed list of the contribution from each source of uncertainty to the total uncertainty in the measured values of $\sigma_{\mathrm{fid}}(tq)$ and $\sigma_{\mathrm{fid}}(\bar tq)$. The estimation of the systematic uncertainties has a statistical uncertainty of $0.3\%$. Uncertainties contributing less than $0.5\%$ are marked with ‘<0.5’.

Significant contributions to the total relative uncertainty in the measured value of $R_{t}$. The estimation of the systematic uncertainties has a statistical uncertainty of $0.3~\%$. Uncertainties contributing less than $0.5~\%$ are not shown.

Slopes $a$ of the mass dependence of the measured cross$-$sections.

Predicted (post-fit) and observed event yields for the signal region (SR), after the requirement on the neural network discriminant, $O_{\mathrm{NN}}~>~0.8$. The multijet background prediction is obtained from the fit to the $E_{\mathrm{T}}^{\mathrm{miss}}$ distribution described in Section 6, while all the other predictions and uncertainties are derived from the total cross$-$section measurement. In some cases there is no uncertainty quoted. In these cases the uncertainty is < 0.5.

Predicted (post-fit) and observed event yields for the signal region (SR), after the requirement on the second neural network discriminant, $O_{\mathrm{NN2}}~>~0.8$. The multijet background prediction is obtained from the fit to the $E_{\mathrm{T}}^{\mathrm{miss}}$ distribution described in Section 6, while all the other predictions and uncertainties are derived from the total cross$-$section measurement. In some cases there is no uncertainty quoted. In these cases the uncertainty is < 0.5.

Migration matrix for $p_{\mathrm{T}}(\hat{t\hspace{-0.2mm}})$ at the particle level. The pseudo top quark is shown on the $y$-axis and the reconstructed variable is shown on the $x$-axis.

Migration matrix for $p_{\mathrm{T}}(t)$ at the parton level. The parton-level quark is shown on the $y$-axis and the reconstructed variable is shown on the $x$-axis.

Migration matrix for $|y(\hat{t\hspace{-0.2mm}})|$ at the particle level. The pseudo top quark is shown on the $y$-axis and the reconstructed variable is shown on the $x$-axis.

Migration matrix for $|y(t)|$ at the parton level. The parton-level quark is shown on the $y$-axis and the reconstructed variable is shown on the $x$-axis.

Uncertainties in the normalisations of the different backgrounds for all processes, as derived from the total cross-section measurement.

Absolute and normalised unfolded differential $tq$ production cross$-$section as a function of $p_{\mathrm{T}}(\hat{t\hspace{-0.2mm}})$ at particle level.

Absolute and normalised unfolded differential $\bar tq$ production cross$-$section as a function of $p_{\mathrm{T}}(\hat{t\hspace{-0.2mm}})$ at particle level.

Absolute and normalised unfolded differential $tq$ production cross$-$section as a function of $|y(\hat{t\hspace{-0.2mm}})|$ at particle level.

Absolute and normalised unfolded differential $\bar tq$ production cross$-$section as a function of $|y(\hat{t\hspace{-0.2mm}})|$ at particle level.

Absolute and normalised unfolded differential $tq$ production cross$-$section as a function of $p_{\mathrm{T}}(\hat{j\hspace{-0.2mm}})$ at particle level.

Absolute and normalised unfolded differential $\bar tq$ production cross$-$section as a function of $p_{\mathrm{T}}(\hat{j\hspace{-0.2mm}})$ at particle level.

Absolute and normalised unfolded differential $tq$ production cross$-$section as a function of $|y(\hat{j\hspace{-0.2mm}})|$ at particle level.

Absolute and normalised unfolded differential $\bar tq$ production cross$-$section as a function of $|y(\hat{j\hspace{-0.2mm}})|$ at particle level.

Absolute and normalised unfolded differential $tq$ production cross$-$section as a function of $p_{\mathrm{T}}(t)$ at parton level.

Absolute and normalised unfolded differential $\bar tq$ production cross$-$section as a function of $p_{\mathrm{T}}(t)$ at parton level.

Absolute and normalised unfolded differential $tq$ production cross$-$section as a function of $|y(t)|$ at parton level.

Absolute and normalised unfolded differential $\bar tq$ production cross$-$section as a function of $|y(t)|$ at parton level.

Statistical correlation matrix for the absolute differential cross-section as a function of $p_{\mathrm{T}}(\hat{t\hspace{-0.2mm}})$ for $tq$ events(at the particle level). It includes the statistical uncertainty due to the number of data events and MC statistics.

Statistical correlation matrix for the absolute differential cross-section as a function of $p_{\mathrm{T}}(\hat{t\hspace{-0.2mm}})$ for $ \bar tq$ events (at the particle level). It includes the statistical uncertainty due to the number of data events and MC statistics.

Statistical correlation matrix for the normalised differential cross-section as a function of $p_{\mathrm{T}}(\hat{t\hspace{-0.2mm}})$ for $tq$ events (at the particle level). It includes the statistical uncertainty due to the number of data events and MC statistics.

Statistical correlation matrix for the normalised differential cross-section as a function of $p_{\mathrm{T}}(\hat{t\hspace{-0.2mm}})$ for $\bar tq$ events (at the particle level). It includes the statistical uncertainty due to the number of data events and MC statistics.

Statistical correlation matrix for the absolute differential cross-section as a function of $|y(\hat{t\hspace{-0.2mm}})|$ for $tq$ events (at the particle level). It includes the statistical uncertainty due to the number of data events and MC statistics.

Statistical correlation matrix for the absolute differential cross-section as a function of $|y(\hat{t\hspace{-0.2mm}})|$ for $ \bar tq$ events (at the particle level). It includes the statistical uncertainty due to the number of data events and MC statistics.

Statistical correlation matrix for the normalised differential cross-section as a function of $|y(\hat{t\hspace{-0.2mm}})|$ for $tq$ events (at the particle level). It includes the statistical uncertainty due to the number of data events and MC statistics.

Statistical correlation matrix for the normalised differential cross-section as a function of $|y(\hat{t\hspace{-0.2mm}})|$ for $\bar tq$ events (at the particle level). It includes the statistical uncertainty due to the number of data events and MC statistics.

Statistical correlation matrix for the absolute differential cross-section as a function of $p_{\mathrm{T}}(\hat{j\hspace{-0.2mm}})$ for $tq$ events (at the particle level). It includes the statistical uncertainty due to the number of data events and MC statistics.

Statistical correlation matrix for the absolute differential cross-section as a function of $p_{\mathrm{T}}(\hat{j\hspace{-0.2mm}})$ for $\bar tq$ events (at the particle level). It includes the statistical uncertainty due to the number of data events and MC statistics.

Statistical correlation matrix for the normalised differential cross-section as a function of $p_{\mathrm{T}}(\hat{j\hspace{-0.2mm}})$ for $tq$ events (at the particle level). It includes the statistical uncertainty due to the number of data events and MC statistics.

Statistical correlation matrix for the normalised differential cross-section as a function of $p_{\mathrm{T}}(\hat{j\hspace{-0.2mm}})$ for $\bar tq$ events (at the particle level). It includes the statistical uncertainty due to the number of data events and MC statistics.

Statistical correlation matrix for the absolute differential cross-section as a function of $|y(\hat{j\hspace{-0.2mm}})|$ for $tq$ events (at the particle level). It includes the statistical uncertainty due to the number of data events and MC statistics.

Statistical correlation matrix for the absolute differential cross-section as a function of $|y(\hat{j\hspace{-0.2mm}})|$ for $\bar tq$ events (at the particle level). It includes the statistical uncertainty due to the number of data events and MC statistics.

Statistical correlation matrix for the normalised differential cross-section as a function of $|y(\hat{j\hspace{-0.2mm}})|$ for $tq$ events (at the particle level). It includes the statistical uncertainty due to the number of data events and MC statistics.

Statistical correlation matrix for the normalised differential cross-section as a function of $|y(\hat{j\hspace{-0.2mm}})|$ for $\bar tq$ events (at the particle level). It includes the statistical uncertainty due to the number of data events and MC statistics.

Statistical correlation matrix for the absolute differential cross-section as a function of $p_{\mathrm{T}}(t)$ for $tq$ events (at the parton level). It includes the statistical uncertainty due to the number of data events and MC statistics.

Statistical correlation matrix for the absolute differential cross-section as a function of $p_{\mathrm{T}}(t)$ for $ \bar tq$ events (at the parton level). It includes the statistical uncertainty due to the number of data events and MC statistics.

Statistical correlation matrix for the normalised differential cross-section as a function of $p_{\mathrm{T}}(t)$ for $tq$ events (at the parton level). It includes the statistical uncertainty due to the number of data events and MC statistics.

Statistical correlation matrix for the normalised differential cross-section as a function of $p_{\mathrm{T}}(t)$ for $ \bar tq$ events (at the parton level). It includes the statistical uncertainty due to the number of data events and MC statistics.

Statistical correlation matrix for the absolute differential cross-section as a function of $|y(t)|$ for $tq$ events (at the parton level). It includes the statistical uncertainty due to the number of data events and MC statistics.

Statistical correlation matrix for the absolute differential cross-section as a function of $|y(t)|$ for $\bar tq$ events (at the parton level). It includes the statistical uncertainty due to the number of data events and MC statistics.

Statistical correlation matrix for the normalised differential cross-section as a function of $|y(t)|$ for $tq$ events (at the parton level). It includes the statistical uncertainty due to the number of data events and MC statistics.

Statistical correlation matrix for the normalised differential cross-section as a function of $|y(t)|$ for $\bar tq$ events (at the parton level). It includes the statistical uncertainty due to the number of data events and MC statistics.

Fiducial acceptance $A_{\mathrm{fid}}$ for different $t$-channel single top-quark MC samples. $^{\mathrm{(a)}}$ Calculation taken from AcerMC $+$ $\mathrm{P{\scriptsize YTHIA}6}$. $^{\mathrm{(b)}}$ Calculation taken from $\mathrm{P{\scriptsize OWHEG}}$-$\mathrm{B{\scriptsize OX}}$ $+$ $\mathrm{P{\scriptsize YTHIA}6}$.

Fiducial acceptance $A_{\mathrm{fid}}$ for different $t$-channel single top-antiquark MC samples. $^{\mathrm{(a)}}$ Calculation taken from AcerMC $+$ $\mathrm{P{\scriptsize YTHIA}6}$. $^{\mathrm{(b)}}$ Calculation taken from $\mathrm{P{\scriptsize OWHEG}}$-$\mathrm{B{\scriptsize OX}}$ $+$ $\mathrm{P{\scriptsize YTHIA}6}$.

Uncertainties for the absolute differential $tq$ cross-section as a function of $p_{\mathrm{T}}(\hat{t\hspace{-0.2mm}})$ at particle level per bin ([0,35,50,75,100,150,200,300] GeV) in percent of $\dfrac{\mathrm{d}\sigma(tq)}{\mathrm{d}p_{\mathrm{T}}(\hat{t\hspace{-0.2mm}})}$. If the uncertainty reported in the paper is "0.0" for both the $\textit{plus}$ and $\textit{minus}$ variation, the value "+0.01" is assigned to the $\textit{plus}$ variation for technical reasons.

Uncertainties for the normalised differential $tq$ cross-section as a function of $p_{\mathrm{T}}(\hat{t\hspace{-0.2mm}})$ at particle level per bin ([0,35,50,75,100,150,200,300] GeV) in percent of $\left( \dfrac{1}{\sigma}\right)\dfrac{\mathrm{d}\sigma(tq)}{\mathrm{d}p_{\mathrm{T}}(\hat{t\hspace{-0.2mm}})}$. If the uncertainty reported in the paper is "0.0" for both the $\textit{plus}$ and $\textit{minus}$ variation, the value "+0.01" is assigned to the $\textit{plus}$ variation for technical reasons.

Uncertainties for the absolute differential $\bar tq$ cross-section as a function of $p_{\mathrm{T}}(\hat{t\hspace{-0.2mm}})$ at particle level per bin ([0,35,50,75,100,150,300] GeV) in percent of $\dfrac{\mathrm{d}\sigma(\bar tq)}{\mathrm{d}p_{\mathrm{T}}(\hat{t\hspace{-0.2mm}})}$.

Uncertainties for the normalised differential $\bar tq$ cross-section as a function of $p_{\mathrm{T}}(\hat{t\hspace{-0.2mm}})$ at particle level per bin ([0,35,50,75,100,150,300] GeV) in percent of $\left( \dfrac{1}{\sigma}\right)\dfrac{\mathrm{d}\sigma(\bar tq)}{\mathrm{d}p_{\mathrm{T}}(\hat{t\hspace{-0.2mm}})}$.

Uncertainties for the absolute differential $tq$ cross-section as a function of $|y(\hat{t\hspace{-0.2mm}})|$ at particle level per bin ([0,0.15,0.3,0.45,0.7,1.0,1.3,2.2]) in percent of $\dfrac{\mathrm{d}\sigma(tq)}{\mathrm{d}|y(\hat{t\hspace{-0.2mm}})|}$. If the uncertainty reported in the paper is "0.0" for both the $\textit{plus}$ and $\textit{minus}$ variation, the value "+0.01" is assigned to the $\textit{plus}$ variation for technical reasons.

Uncertainties for the normalised differential $tq$ cross-section as a function of $|y(\hat{t\hspace{-0.2mm}})|$ at particle level per bin ([0,0.15,0.3,0.45,0.7,1.0,1.3,2.2]) in percent of $\left(\dfrac{1}{\sigma}\right)\dfrac{\mathrm{d}\sigma(tq)}{\mathrm{d}|y(\hat{t\hspace{-0.2mm}})|}$. If the uncertainty reported in the paper is "0.0" for both the $\textit{plus}$ and $\textit{minus}$ variation, the value "+0.01" is assigned to the $\textit{plus}$ variation for technical reasons.

Uncertainties for the absolute differential $\bar tq$ cross-section as a function of $|y(\hat{t\hspace{-0.2mm}})|$ at particle level per bin ([0,0.15,0.3,0.45,0.7,1.0,1.3,2.2]) in percent of $\dfrac{\mathrm{d}\sigma(\bar tq)}{\mathrm{d}|y(\hat{t\hspace{-0.2mm}})|}$. If the uncertainty reported in the paper is "0.0" for both the $\textit{plus}$ and $\textit{minus}$ variation, the value "+0.01" is assigned to the $\textit{plus}$ variation for technical reasons.

Uncertainties for the normalised differential $\bar tq$ cross-section as a function of $|y(\hat{t\hspace{-0.2mm}})|$ at particle level per bin ([0,0.15,0.3,0.45,0.7,1.0,1.3,2.2]) in percent of $\left(\dfrac{1}{\sigma}\right)\dfrac{\mathrm{d}\sigma(\bar tq)}{\mathrm{d}|y(\hat{t\hspace{-0.2mm}})|}$. If the uncertainty reported in the paper is "0.0" for both the $\textit{plus}$ and $\textit{minus}$ variation, the value "+0.01" is assigned to the $\textit{plus}$ variation for technical reasons.

Uncertainties for the absolute differential $tq$ cross-section as a function of $p_{\mathrm{T}}(\hat{j\hspace{-0.2mm}})$ at particle level per bin ([30,45,60,75,100,150,300] GeV) in percent of $\dfrac{\mathrm{d}\sigma(tq)}{\mathrm{d}p_{\mathrm{T}}(\hat{j\hspace{-0.2mm}})}$. If the uncertainty reported in the paper is "0.0" for both the $\textit{plus}$ and $\textit{minus}$ variation, the value "+0.01" is assigned to the $\textit{plus}$ variation for technical reasons.

Uncertainties for the normalised differential $tq$ cross-section as a function of $p_{\mathrm{T}}(\hat{j\hspace{-0.2mm}})$ at particle level per bin ([30,45,60,75,100,150,300] GeV) in percent of $\left(\dfrac{1}{\sigma}\right)\dfrac{\mathrm{d}\sigma(tq)}{\mathrm{d}p_{\mathrm{T}}(\hat{j\hspace{-0.2mm}})}$. If the uncertainty reported in the paper is "0.0" for both the $\textit{plus}$ and $\textit{minus}$ variation, the value "+0.01" is assigned to the $\textit{plus}$ variation for technical reasons.

Uncertainties for the absolute differential $\bar tq$ cross-section as a function of $p_{\mathrm{T}}(\hat{j\hspace{-0.2mm}})$ at particle level per bin ([30,45,60,75,100,150,300] GeV) in percent of $\dfrac{\mathrm{d}\sigma(\bar tq)}{\mathrm{d}p_{\mathrm{T}}(\hat{j\hspace{-0.2mm}})}$. If the uncertainty reported in the paper is "0.0" for both the $\textit{plus}$ and $\textit{minus}$ variation, the value "+0.01" is assigned to the $\textit{plus}$ variation for technical reasons.

Uncertainties for the normalised differential $\bar tq$ cross-section as a function of $p_{\mathrm{T}}(\hat{j\hspace{-0.2mm}})$ at particle level per bin ([30,45,60,75,100,150,300] GeV) in percent of $\left(\dfrac{1}{\sigma}\right)\dfrac{\mathrm{d}\sigma(\bar tq)}{\mathrm{d}p_{\mathrm{T}}(\hat{j\hspace{-0.2mm}})}$. If the uncertainty reported in the paper is "0.0" for both the $\textit{plus}$ and $\textit{minus}$ variation, the value "+0.01" is assigned to the $\textit{plus}$ variation for technical reasons.

Uncertainties for the absolute differential $tq$ cross-section as a function of $|y(\hat{j\hspace{-0.2mm}})|$ at particle level per bin ([0.0, 1.2, 1.7, 2.2, 2.7, 3.3, 4.5]) in percent of $\dfrac{\mathrm{d}\sigma(tq)}{\mathrm{d}|y(\hat{j\hspace{-0.2mm}})|}$. If the uncertainty reported in the paper is "0.0" for both the $\textit{plus}$ and $\textit{minus}$ variation, the value "+0.01" is assigned to the $\textit{plus}$ variation for technical reasons.

Uncertainties for the normalised differential $tq$ cross-section as a function of $|y(\hat{j\hspace{-0.2mm}})|$ at particle level per bin ([0.0, 1.2, 1.7, 2.2, 2.7, 3.3, 4.5]) in percent of $\left(\dfrac{1}{\sigma}\right)\dfrac{\mathrm{d}\sigma(tq)}{\mathrm{d}|y(\hat{j\hspace{-0.2mm}})|}$. If the uncertainty reported in the paper is "0.0" for both the $\textit{plus}$ and $\textit{minus}$ variation, the value "+0.01" is assigned to the $\textit{plus}$ variation for technical reasons.

Uncertainties for the absolute differential $\bar tq$ cross-section as a function of $|y(\hat{j\hspace{-0.2mm}})|$ at particle level per bin ([0.0, 1.2, 1.7, 2.2, 2.7, 3.3, 4.5]) in percent of $\dfrac{\mathrm{d}\sigma(\bar tq)}{\mathrm{d}|y(\hat{j\hspace{-0.2mm}})|}$. If the uncertainty reported in the paper is "0.0" for both the $\textit{plus}$ and $\textit{minus}$ variation, the value "+0.01" is assigned to the $\textit{plus}$ variation for technical reasons.

Uncertainties for the normalised differential $\bar tq$ cross-section as a function of $|y(\hat{j\hspace{-0.2mm}})|$ at particle level per bin ([0.0, 1.2, 1.7, 2.2, 2.7, 3.3, 4.5]) in percent of $\left(\dfrac{1}{\sigma}\right)\dfrac{\mathrm{d}\sigma(\bar tq)}{\mathrm{d}|y(\hat{j\hspace{-0.2mm}})|}$. If the uncertainty reported in the paper is "0.0" for both the $\textit{plus}$ and $\textit{minus}$ variation, the value "+0.01" is assigned to the $\textit{plus}$ variation for technical reasons.

Uncertainties for the absolute differential $tq$ cross-section as a function of $p_{\mathrm{T}}(t)$ at parton level per bin ([0,50,100,150,200,300] GeV) in percent of $\dfrac{\mathrm{d}\sigma(tq)}{\mathrm{d}p_{\mathrm{T}}(t)}$.

Uncertainties for the normalised differential $tq$ cross-section as a function of $p_{\mathrm{T}}(t)$ at parton level per bin ([0,50,100,150,200,300] GeV) in percent of $\left(\dfrac{1}{\sigma}\right)\dfrac{\mathrm{d}\sigma(tq)}{\mathrm{d}p_{\mathrm{T}}(t)}$. If the uncertainty reported in the paper is "0.0" for both the $\textit{plus}$ and $\textit{minus}$ variation, the value "+0.01" is assigned to the $\textit{plus}$ variation for technical reasons.

Uncertainties for the absolute differential $\bar tq $ cross-section as a function of $p_{\mathrm{T}}(t)$ at parton level per bin ([0,50,100,150,300] GeV) in percent of $\dfrac{\mathrm{d}\sigma(\bar tq)}{\mathrm{d}p_{\mathrm{T}}(t)}$.

Uncertainties for the normalised differential $\bar tq $ cross-section as a function of $p_{\mathrm{T}}(t)$ at parton level per bin ([0,50,100,150,300] GeV) in percent of $\left(\dfrac{1}{\sigma}\right)\dfrac{\mathrm{d}\sigma(\bar tq)}{\mathrm{d}p_{\mathrm{T}}(t)}$.

Uncertainties for the absolute differential $ tq $ cross-section as a function of $|y(t)|$ at parton level per bin ([0,0.3,0.7,1.3,2.2]) in percent of $\dfrac{\mathrm{d}\sigma(tq)}{\mathrm{d}|y(t)|}$. If the uncertainty reported in the paper is "0.0" for both the $\textit{plus}$ and $\textit{minus}$ variation, the value "+0.01" is assigned to the $\textit{plus}$ variation for technical reasons.

Uncertainties for the normalised differential $ tq $ cross-section as a function of $|y(t)|$ at parton level per bin ([0,0.3,0.7,1.3,2.2]) in percent of $\left(\dfrac{1}{\sigma}\right)\dfrac{\mathrm{d}\sigma(tq)}{\mathrm{d}|y(t)|}$. If the uncertainty reported in the paper is "0.0" for both the $\textit{plus}$ and $\textit{minus}$ variation, the value "+0.01" is assigned to the $\textit{plus}$ variation for technical reasons.

Uncertainties for the absolute differential $ \bar tq $ cross-section as a function of $|y(t)|$ at parton level per bin ([0,0.3,0.7,1.3,2.2]) in percent of $\dfrac{\mathrm{d}\sigma(\bar tq)}{\mathrm{d}|y(t)|}$.

Uncertainties for the normalised differential $ \bar tq $ cross-section as a function of $|y(t)|$ at parton level per bin ([0,0.3,0.7,1.3,2.2]) in percent of $\left(\dfrac{1}{\sigma}\right)\dfrac{\mathrm{d}\sigma(\bar tq)}{\mathrm{d}|y(t)|}$. If the uncertainty reported in the paper is "0.0" for both the $\textit{plus}$ and $\textit{minus}$ variation, the value "+0.01" is assigned to the $\textit{plus}$ variation for technical reasons.