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.

IRing2 for HT2>=500 GeV, NJets>=5

IRing2 for HT2>=1000 GeV, NJets>=2

IRing2 for HT2>=1000 GeV, NJets>=3

IRing2 for HT2>=1000 GeV, NJets>=4

IRing2 for HT2>=1000 GeV, NJets>=5

IRing2 for HT2>=1500 GeV, NJets>=2

IRing2 for HT2>=1500 GeV, NJets>=3

IRing2 for HT2>=1500 GeV, NJets>=4

IRing2 for HT2>=1500 GeV, NJets>=5

IRing128 for HT2>=500 GeV, NJets>=2

IRing128 for HT2>=500 GeV, NJets>=3

IRing128 for HT2>=500 GeV, NJets>=4

IRing128 for HT2>=500 GeV, NJets>=5

IRing128 for HT2>=1000 GeV, NJets>=2

IRing128 for HT2>=1000 GeV, NJets>=3

IRing128 for HT2>=1000 GeV, NJets>=4

IRing128 for HT2>=1000 GeV, NJets>=5

IRing128 for HT2>=1500 GeV, NJets>=2

IRing128 for HT2>=1500 GeV, NJets>=3

IRing128 for HT2>=1500 GeV, NJets>=4

IRing128 for HT2>=1500 GeV, NJets>=5

ICyl16 for HT2>=500 GeV, NJets>=2

ICyl16 for HT2>=500 GeV, NJets>=3

ICyl16 for HT2>=500 GeV, NJets>=4

ICyl16 for HT2>=500 GeV, NJets>=5

ICyl16 for HT2>=1000 GeV, NJets>=2

ICyl16 for HT2>=1000 GeV, NJets>=3

ICyl16 for HT2>=1000 GeV, NJets>=4

ICyl16 for HT2>=1000 GeV, NJets>=5

ICyl16 for HT2>=1500 GeV, NJets>=2

ICyl16 for HT2>=1500 GeV, NJets>=3

ICyl16 for HT2>=1500 GeV, NJets>=4

ICyl16 for HT2>=1500 GeV, NJets>=5

IRing2 covariance for HT2>=500 GeV, NJets>=2 (Table 1)

IRing2 covariance for HT2>=500 GeV, NJets>=3 (Table 2)

IRing2 covariance for HT2>=500 GeV, NJets>=4 (Table 3)

IRing2 covariance for HT2>=500 GeV, NJets>=5 (Table 4)

IRing2 covariance for HT2>=1000 GeV, NJets>=2 (Table 5)

IRing2 covariance for HT2>=1000 GeV, NJets>=3 (Table 6)

IRing2 covariance for HT2>=1000 GeV, NJets>=4 (Table 7)

IRing2 covariance for HT2>=1000 GeV, NJets>=5 (Table 8)

IRing2 covariance for HT2>=1500 GeV, NJets>=2 (Table 9)

IRing2 covariance for HT2>=1500 GeV, NJets>=3 (Table 10)

IRing2 covariance for HT2>=1500 GeV, NJets>=4 (Table 11)

IRing2 covariance for HT2>=1500 GeV, NJets>=5 (Table 12)

IRing128 covariance for HT2>=500 GeV, NJets>=2 (Table 13)

IRing128 covariance for HT2>=500 GeV, NJets>=3 (Table 14)

IRing128 covariance for HT2>=500 GeV, NJets>=4 (Table 15)

IRing128 covariance for HT2>=500 GeV, NJets>=5 (Table 16)

IRing128 covariance for HT2>=1000 GeV, NJets>=2 (Table 17)

IRing128 covariance for HT2>=1000 GeV, NJets>=3 (Table 18)

IRing128 covariance for HT2>=1000 GeV, NJets>=4 (Table 19)

IRing128 covariance for HT2>=1000 GeV, NJets>=5 (Table 20)

IRing128 covariance for HT2>=1500 GeV, NJets>=2 (Table 21)

IRing128 covariance for HT2>=1500 GeV, NJets>=3 (Table 22)

IRing128 covariance for HT2>=1500 GeV, NJets>=4 (Table 23)

IRing128 covariance for HT2>=1500 GeV, NJets>=5 (Table 24)

ICyl16 covariance for HT2>=500 GeV, NJets>=2 (Table 25)

ICyl16 covariance for HT2>=500 GeV, NJets>=3 (Table 26)

ICyl16 covariance for HT2>=500 GeV, NJets>=4 (Table 27)

ICyl16 covariance for HT2>=500 GeV, NJets>=5 (Table 28)

ICyl16 covariance for HT2>=1000 GeV, NJets>=2 (Table 29)

ICyl16 covariance for HT2>=1000 GeV, NJets>=3 (Table 30)

ICyl16 covariance for HT2>=1000 GeV, NJets>=4 (Table 31)

ICyl16 covariance for HT2>=1000 GeV, NJets>=5 (Table 32)

ICyl16 covariance for HT2>=1500 GeV, NJets>=2 (Table 33)

ICyl16 covariance for HT2>=1500 GeV, NJets>=3 (Table 34)

ICyl16 covariance for HT2>=1500 GeV, NJets>=4 (Table 35)

ICyl16 covariance for HT2>=1500 GeV, NJets>=5 (Table 36)

IRing2 covariance, complete

1-IRing128 covariance, complete

1-ICyl16 covariance, complete

Observed (solid line) and expected (dashed line) 95% CL upper limits on the cross-section of single plus non-resonant plus pair vector LQ production from the combination of the high b-jet $p_{T}$ signal region for the $\tau_\text{lep}\tau_\text{had}$ and $\tau_\text{had}\tau_\text{had}$ channels. [$U_1^{MIN}$ model ($\kappa$ = 1) with $\lambda$ = 1.0]

Observed (solid line) and expected (dashed line) 95% CL upper limits on the cross-section of single plus non-resonant plus pair vector LQ production from the combination of the high b-jet $p_{T}$ signal region for the $\tau_\text{lep}\tau_\text{had}$ and $\tau_\text{had}\tau_\text{had}$ channels. [$U_1^{MIN}$ model ($\kappa$ = 1) with $\lambda$ = 1.7]

Observed (solid line) and expected (dashed line) 95% CL upper limits on the cross-section of single plus non-resonant plus pair vector LQ production from the combination of the high b-jet $p_{T}$ signal region for the $\tau_\text{lep}\tau_\text{had}$ and $\tau_\text{had}\tau_\text{had}$ channels. [$U_1^{MIN}$ model ($\kappa$ = 1) with $\lambda$ = 2.5]

The two-dimensional 95% CL exclusion limits in the $\lambda$-$m_{LQ}$ plane for singly plus non-resonant produced vector LQ (green lines) and for the sum, referred as Total, of single plus non-resonant plus pair vector LQ production (blue lines). Regions to the left of the lines are excluded. The dotted area shows the preferred region where the chosen LQ model can explain observed B anomalies. [$\kappa$ = 0]

The two-dimensional 95% CL exclusion limits in the $\lambda$-$m_{LQ}$ plane for singly plus non-resonant produced vector LQ (green lines) and for the sum, referred as Total, of single plus non-resonant plus pair vector LQ production (blue lines). Regions to the left of the lines are excluded. The dotted area shows the preferred region where the chosen LQ model can explain observed B anomalies. [$\kappa$ = 1]

Observed (solid line) and expected (dashed line) 95% CL upper limits for $\lambda$ = 1.0 on the cross-section of single plus non-resonant plus pair $\widetilde{S_{1}}$ production hypotheses from the combination of the high b-jet $p_{T}$ category for the $\tau_\text{lep}\tau_\text{had}$ and $\tau_\text{had}\tau_\text{had}$ channels.

Observed (solid line) and expected (dashed line) 95% CL upper limits for $\lambda$ = 1.7 on the cross-section of single plus non-resonant plus pair $\widetilde{S_{1}}$ production hypotheses from the combination of the high b-jet $p_{T}$ category for the $\tau_\text{lep}\tau_\text{had}$ and $\tau_\text{had}\tau_\text{had}$ channels.

Observed (solid line) and expected (dashed line) 95% CL upper limits for $\lambda$ = 2.5 on the cross-section of single plus non-resonant plus pair $\widetilde{S_{1}}$ production hypotheses from the combination of the high b-jet $p_{T}$ category for the $\tau_\text{lep}\tau_\text{had}$ and $\tau_\text{had}\tau_\text{had}$ channels.

The two-dimensional 95% CL exclusion limits in the $\lambda$-$m_{LQ}$ plane for singly plus non-resonant produced $\widetilde{S_{1}}$ (green lines) and for the sum, referred as Total, of single plus non-resonant plus pair vector LQ production (blue lines). Regions to the left of the lines are excluded.

Observed (solid line) and expected (dashed line) 95% upper limits on the visible cross-section, $\sigma_{\mathrm{vis}}$, obtained from model-independent search by a signal-plus-background fit in the high and low b-jet $p_{T}$ category for the $\tau_\text{lep}\tau_\text{had}$ and $\tau_\text{had}\tau_\text{had}$ channels.

Observed (solid line) and expected (dashed line) 95% CL upper limits on the cross-section of singly plus non-resonant produced vector LQ signals from the combination of the high b-jet $p_{T}$ category for the $\tau_\text{lep}\tau_\text{had}$ and $\tau_\text{had}\tau_\text{had}$ channels. [$U_1^{YM}$ model ($\kappa$ = 0) with $\lambda$ = 1.0]

Observed (solid line) and expected (dashed line) 95% CL upper limits on the cross-section of singly plus non-resonant produced vector LQ signals from the combination of the high b-jet $p_{T}$ category for the $\tau_\text{lep}\tau_\text{had}$ and $\tau_\text{had}\tau_\text{had}$ channels. [$U_1^{YM}$ model ($\kappa$ = 0) with $\lambda$ = 1.7]

Observed (solid line) and expected (dashed line) 95% CL upper limits on the cross-section of singly plus non-resonant produced vector LQ signals from the combination of the high b-jet $p_{T}$ category for the $\tau_\text{lep}\tau_\text{had}$ and $\tau_\text{had}\tau_\text{had}$ channels. [$U_1^{YM}$ model ($\kappa$ = 0) with $\lambda$ = 2.5]

Observed (solid line) and expected (dashed line) 95% CL upper limits on the cross-section of singly plus non-resonant produced vector LQ signals from the combination of the high b-jet $p_{T}$ category for the $\tau_\text{lep}\tau_\text{had}$ and $\tau_\text{had}\tau_\text{had}$ channels. [$U_1^{MIN}$ model ($\kappa$ = 1) with $\lambda$ = 1.0]

Observed (solid line) and expected (dashed line) 95% CL upper limits on the cross-section of singly plus non-resonant produced vector LQ signals from the combination of the high b-jet $p_{T}$ category for the $\tau_\text{lep}\tau_\text{had}$ and $\tau_\text{had}\tau_\text{had}$ channels. [$U_1^{MIN}$ model ($\kappa$ = 1) with $\lambda$ = 1.7]

Observed (solid line) and expected (dashed line) 95% CL upper limits on the cross-section of singly plus non-resonant produced vector LQ signals from the combination of the high b-jet $p_{T}$ category for the $\tau_\text{lep}\tau_\text{had}$ and $\tau_\text{had}\tau_\text{had}$ channels. [$U_1^{MIN}$ model ($\kappa$ = 1) with $\lambda$ = 2.5]

Observed (solid line) and expected (dashed line) 95% CL upper limits for $\lambda$ = 1.0 on the cross-section of singly produced $\widetilde{S_{1}}$ signal hypotheses from the combination of the high b-jet $p_{T}$ category for the $\tau_\text{lep}\tau_\text{had}$ and $\tau_\text{had}\tau_\text{had}$ channels.

Observed (solid line) and expected (dashed line) 95% CL upper limits for $\lambda$ = 1.7 on the cross-section of singly produced $\widetilde{S_{1}}$ signal hypotheses from the combination of the high b-jet $p_{T}$ category for the $\tau_\text{lep}\tau_\text{had}$ and $\tau_\text{had}\tau_\text{had}$ channels.

Observed (solid line) and expected (dashed line) 95% CL upper limits for $\lambda$ = 2.5 on the cross-section of singly produced $\widetilde{S_{1}}$ signal hypotheses from the combination of the high b-jet $p_{T}$ category for the $\tau_\text{lep}\tau_\text{had}$ and $\tau_\text{had}\tau_\text{had}$ channels.

The signal acceptance times efficiency for single plus non-resonant produced LQ in the $\widetilde{S_{1}}$ model from the combination of the high and low b-jet $p_{T}$ signal regions for the $\tau_\text{lep}\tau_\text{had}$ and $\tau_\text{had}\tau_\text{had}$ channels. [$\lambda$ = 1.0]

The signal acceptance times efficiency for single plus non-resonant produced LQ in the $\widetilde{S_{1}}$ model from the combination of the high and low b-jet $p_{T}$ signal regions for the $\tau_\text{lep}\tau_\text{had}$ and $\tau_\text{had}\tau_\text{had}$ channels. [$\lambda$ = 1.7]

The signal acceptance times efficiency for single plus non-resonant produced LQ in the $\widetilde{S_{1}}$ model from the combination of the high and low b-jet $p_{T}$ signal regions for the $\tau_\text{lep}\tau_\text{had}$ and $\tau_\text{had}\tau_\text{had}$ channels. [$\lambda$ = 2.5]

The signal acceptance times efficiency for single plus non-resonant produced LQ in the $U_1^{MIN}$ model from the combination of the high and low b-jet $p_{T}$ signal regions for the $\tau_\text{lep}\tau_\text{had}$ and $\tau_\text{had}\tau_\text{had}$ channels. [$\lambda$ = 1.0]

The signal acceptance times efficiency for single plus non-resonant produced LQ in the $U_1^{MIN}$ model from the combination of the high and low b-jet $p_{T}$ signal regions for the $\tau_\text{lep}\tau_\text{had}$ and $\tau_\text{had}\tau_\text{had}$ channels. [$\lambda$ = 1.7]

The signal acceptance times efficiency for single plus non-resonant produced LQ in the $U_1^{MIN}$ model from the combination of the high and low b-jet $p_{T}$ signal regions for the $\tau_\text{lep}\tau_\text{had}$ and $\tau_\text{had}\tau_\text{had}$ channels. [$\lambda$ = 2.5]

The signal acceptance times efficiency for single plus non-resonant produced LQ in the $U_1^{YM}$ model from the combination of the high and low b-jet $p_{T}$ signal regions for the $\tau_\text{lep}\tau_\text{had}$ and $\tau_\text{had}\tau_\text{had}$ channels. [$\lambda$ = 1.0]

The signal acceptance times efficiency for single plus non-resonant produced LQ in the $U_1^{YM}$ model from the combination of the high and low b-jet $p_{T}$ signal regions for the $\tau_\text{lep}\tau_\text{had}$ and $\tau_\text{had}\tau_\text{had}$ channels. [$\lambda$ = 1.7]

The signal acceptance times efficiency for single plus non-resonant produced LQ in the $U_1^{YM}$ model from the combination of the high and low b-jet $p_{T}$ signal regions for the $\tau_\text{lep}\tau_\text{had}$ and $\tau_\text{had}\tau_\text{had}$ channels. [$\lambda$ = 2.5]

Observed (solid line) and expected (dashed line) 95% CL upper limits on the cross-section of singly produced vector LQ signals from the combination of the high and low b-jet $p_{T}$ category for the $\tau_\text{lep}\tau_\text{had}$ and $\tau_\text{had}\tau_\text{had}$ channels, neglecting the interference between the LQ non-resonant production processes and the SM processes. However, the interference between the non-resonant LQ production and SM diagrams, such as those from Z+jets, could be non-negligible for events in the low b-jet $p_{T}$ signal region. [$U_1^{YM}$ model ($\kappa$ = 0) with $\lambda$ = 1.0]

Observed (solid line) and expected (dashed line) 95% CL upper limits on the cross-section of singly produced vector LQ signals from the combination of the high and low b-jet $p_{T}$ category for the $\tau_\text{lep}\tau_\text{had}$ and $\tau_\text{had}\tau_\text{had}$ channels, neglecting the interference between the LQ non-resonant production processes and the SM processes. However, the interference between the non-resonant LQ production and SM diagrams, such as those from Z+jets, could be non-negligible for events in the low b-jet $p_{T}$ signal region. [$U_1^{YM}$ model ($\kappa$ = 0) with $\lambda$ = 1.7]

Observed (solid line) and expected (dashed line) 95% CL upper limits on the cross-section of singly produced vector LQ signals from the combination of the high and low b-jet $p_{T}$ category for the $\tau_\text{lep}\tau_\text{had}$ and $\tau_\text{had}\tau_\text{had}$ channels, neglecting the interference between the LQ non-resonant production processes and the SM processes. However, the interference between the non-resonant LQ production and SM diagrams, such as those from Z+jets, could be non-negligible for events in the low b-jet $p_{T}$ signal region. [$U_1^{YM}$ model ($\kappa$ = 0) with $\lambda$ = 2.5]

Observed (solid line) and expected (dashed line) 95% CL upper limits on the cross-section of singly produced vector LQ signals from the combination of the high and low b-jet $p_{T}$ category for the $\tau_\text{lep}\tau_\text{had}$ and $\tau_\text{had}\tau_\text{had}$ channels, neglecting the interference between the LQ non-resonant production processes and the SM processes. However, the interference between the non-resonant LQ production and SM diagrams, such as those from Z+jets, could be non-negligible for events in the low b-jet $p_{T}$ signal region. [$U_1^{MIN}$ model ($\kappa$ = 1) with $\lambda$ = 1.0]

Observed (solid line) and expected (dashed line) 95% CL upper limits on the cross-section of singly produced vector LQ signals from the combination of the high and low b-jet $p_{T}$ category for the $\tau_\text{lep}\tau_\text{had}$ and $\tau_\text{had}\tau_\text{had}$ channels, neglecting the interference between the LQ non-resonant production processes and the SM processes. However, the interference between the non-resonant LQ production and SM diagrams, such as those from Z+jets, could be non-negligible for events in the low b-jet $p_{T}$ signal region. [$U_1^{MIN}$ model ($\kappa$ = 1) with $\lambda$ = 1.7]

Observed (solid line) and expected (dashed line) 95% CL upper limits on the cross-section of singly produced vector LQ signals from the combination of the high and low b-jet $p_{T}$ category for the $\tau_\text{lep}\tau_\text{had}$ and $\tau_\text{had}\tau_\text{had}$ channels, neglecting the interference between the LQ non-resonant production processes and the SM processes. However, the interference between the non-resonant LQ production and SM diagrams, such as those from Z+jets, could be non-negligible for events in the low b-jet $p_{T}$ signal region. [$U_1^{MIN}$ model ($\kappa$ = 1) with $\lambda$ = 2.5]

Observed (solid line) and expected (dashed line) 95% CL upper limits on the cross-section of the single plus non-resonant plus pair production of vector LQ from the combination of the high and low b-jet $p_{T}$ category for the $\tau_\text{lep}\tau_\text{had}$ and $\tau_\text{had}\tau_\text{had}$ channels, neglecting the interference between the LQ non-resonant production processes and the SM processes. However, the interference between the non-resonant LQ production and SM diagrams, such as those from Z+jets, could be non-negligible for events in the low b-jet $p_{T}$ signal region. [$U_1^{YM}$ model ($\kappa$ = 0) with $\lambda$ = 1.0]

Observed (solid line) and expected (dashed line) 95% CL upper limits on the cross-section of the single plus non-resonant plus pair production of vector LQ from the combination of the high and low b-jet $p_{T}$ category for the $\tau_\text{lep}\tau_\text{had}$ and $\tau_\text{had}\tau_\text{had}$ channels, neglecting the interference between the LQ non-resonant production processes and the SM processes. However, the interference between the non-resonant LQ production and SM diagrams, such as those from Z+jets, could be non-negligible for events in the low b-jet $p_{T}$ signal region. [$U_1^{YM}$ model ($\kappa$ = 0) with $\lambda$ = 1.7]

Observed (solid line) and expected (dashed line) 95% CL upper limits on the cross-section of the single plus non-resonant plus pair production of vector LQ from the combination of the high and low b-jet $p_{T}$ category for the $\tau_\text{lep}\tau_\text{had}$ and $\tau_\text{had}\tau_\text{had}$ channels, neglecting the interference between the LQ non-resonant production processes and the SM processes. However, the interference between the non-resonant LQ production and SM diagrams, such as those from Z+jets, could be non-negligible for events in the low b-jet $p_{T}$ signal region. [$U_1^{YM}$ model ($\kappa$ = 0) with $\lambda$ = 2.5]

Observed (solid line) and expected (dashed line) 95% CL upper limits on the cross-section of the single plus non-resonant plus pair production of vector LQ from the combination of the high and low b-jet $p_{T}$ category for the $\tau_\text{lep}\tau_\text{had}$ and $\tau_\text{had}\tau_\text{had}$ channels, neglecting the interference between the LQ non-resonant production processes and the SM processes. However, the interference between the non-resonant LQ production and SM diagrams, such as those from Z+jets, could be non-negligible for events in the low b-jet $p_{T}$ signal region. [$U_1^{MIN}$ model ($\kappa$ = 1) with $\lambda$ = 1.0]

Observed (solid line) and expected (dashed line) 95% CL upper limits on the cross-section of the single plus non-resonant plus pair production of vector LQ from the combination of the high and low b-jet $p_{T}$ category for the $\tau_\text{lep}\tau_\text{had}$ and $\tau_\text{had}\tau_\text{had}$ channels, neglecting the interference between the LQ non-resonant production processes and the SM processes. However, the interference between the non-resonant LQ production and SM diagrams, such as those from Z+jets, could be non-negligible for events in the low b-jet $p_{T}$ signal region. [$U_1^{MIN}$ model ($\kappa$ = 1) with $\lambda$ = 1.7]

Observed (solid line) and expected (dashed line) 95% CL upper limits on the cross-section of the single plus non-resonant plus pair production of vector LQ from the combination of the high and low b-jet $p_{T}$ category for the $\tau_\text{lep}\tau_\text{had}$ and $\tau_\text{had}\tau_\text{had}$ channels, neglecting the interference between the LQ non-resonant production processes and the SM processes. However, the interference between the non-resonant LQ production and SM diagrams, such as those from Z+jets, could be non-negligible for events in the low b-jet $p_{T}$ signal region. [$U_1^{MIN}$ model ($\kappa$ = 1) with $\lambda$ = 2.5]

The two-dimensional 95% CL exclusion limits in the $\lambda$-$m_{LQ}$ plane for singly plus non-resonant produced vector LQ (green lines) and for the sum, referred as Total, of single plus non-resonant plus pair vector LQ production (blue lines).Regions to the left of the lines are excluded. The dotted area shows the preferred region where the chosen LQ model can explain observed B anomalies. Results are extracted from the combination of the high and low b-jet $p_{T}$ signal region for the $\tau_\text{lep}\tau_\text{had}$ and $\tau_\text{had}\tau_\text{had}$ channels, neglecting the interference between the LQ non-resonant production processes and the SM processes. However, the interference between the non-resonant LQ production and SM diagrams, such as those from Z+jets, could be non-negligible for events in the low b-jet $p_{T}$ signal region. [$\kappa$ = 0]

The two-dimensional 95% CL exclusion limits in the $\lambda$-$m_{LQ}$ plane for singly plus non-resonant produced vector LQ (green lines) and for the sum, referred as Total, of single plus non-resonant plus pair vector LQ production (blue lines).Regions to the left of the lines are excluded. The dotted area shows the preferred region where the chosen LQ model can explain observed B anomalies. Results are extracted from the combination of the high and low b-jet $p_{T}$ signal region for the $\tau_\text{lep}\tau_\text{had}$ and $\tau_\text{had}\tau_\text{had}$ channels, neglecting the interference between the LQ non-resonant production processes and the SM processes. However, the interference between the non-resonant LQ production and SM diagrams, such as those from Z+jets, could be non-negligible for events in the low b-jet $p_{T}$ signal region. [$\kappa$ = 1]

Observed (solid line) and expected (dashed line) 95% CL upper limits for $\lambda$ = 1.0 on the cross-section of singly-produced $\widetilde{S_{1}}$ signals from the combination of the high and low b-jet $p_{T}$ categories for the $\tau_\text{lep}\tau_\text{had}$ and $\tau_\text{had}\tau_\text{had}$ channels, neglecting the interference between the LQ non-resonant production processes and the SM processes. However, the interference between the non-resonant LQ production and SM diagrams, such as those from Z+jets, could be non-negligible for events in the low b-jet $p_{T}$ signal region.

Observed (solid line) and expected (dashed line) 95% CL upper limits for $\lambda$ = 1.7 on the cross-section of singly-produced $\widetilde{S_{1}}$ signals from the combination of the high and low b-jet $p_{T}$ categories for the $\tau_\text{lep}\tau_\text{had}$ and $\tau_\text{had}\tau_\text{had}$ channels, neglecting the interference between the LQ non-resonant production processes and the SM processes. However, the interference between the non-resonant LQ production and SM diagrams, such as those from Z+jets, could be non-negligible for events in the low b-jet $p_{T}$ signal region.

Observed (solid line) and expected (dashed line) 95% CL upper limits for $\lambda$ = 2.5 on the cross-section of singly-produced $\widetilde{S_{1}}$ signals from the combination of the high and low b-jet $p_{T}$ categories for the $\tau_\text{lep}\tau_\text{had}$ and $\tau_\text{had}\tau_\text{had}$ channels, neglecting the interference between the LQ non-resonant production processes and the SM processes. However, the interference between the non-resonant LQ production and SM diagrams, such as those from Z+jets, could be non-negligible for events in the low b-jet $p_{T}$ signal region.

Observed (solid line) and expected (dashed line) 95% CL upper limits for $\lambda$ = 1.0 on the cross-section of single plus non-resonant plus pair $\widetilde{S_{1}}$ production signal from the combination of the high and low b-jet $p_{T}$ categories for the $\tau_\text{lep}\tau_\text{had}$ and $\tau_\text{had}\tau_\text{had}$ channels, neglecting the interference between the LQ non-resonant production processes and the SM processes. However, the interference between the non-resonant LQ production and SM diagrams, such as those from Z+jets, could be non-negligible for events in the low b-jet $p_{T}$ signal region.

Observed (solid line) and expected (dashed line) 95% CL upper limits for $\lambda$ = 1.7 on the cross-section of single plus non-resonant plus pair $\widetilde{S_{1}}$ production signal from the combination of the high and low b-jet $p_{T}$ categories for the $\tau_\text{lep}\tau_\text{had}$ and $\tau_\text{had}\tau_\text{had}$ channels, neglecting the interference between the LQ non-resonant production processes and the SM processes. However, the interference between the non-resonant LQ production and SM diagrams, such as those from Z+jets, could be non-negligible for events in the low b-jet $p_{T}$ signal region.

Observed (solid line) and expected (dashed line) 95% CL upper limits for $\lambda$ = 2.5 on the cross-section of single plus non-resonant plus pair $\widetilde{S_{1}}$ production signal from the combination of the high and low b-jet $p_{T}$ categories for the $\tau_\text{lep}\tau_\text{had}$ and $\tau_\text{had}\tau_\text{had}$ channels, neglecting the interference between the LQ non-resonant production processes and the SM processes. However, the interference between the non-resonant LQ production and SM diagrams, such as those from Z+jets, could be non-negligible for events in the low b-jet $p_{T}$ signal region.

The two-dimensional 95% CL exclusion limits in the $\lambda$-$m_{LQ}$ plane for singly plus non-resonant produced $\widetilde{S_{1}}$ (green lines) and for the sum, referred as Total, of single plus non-resonant plus pair $\widetilde{S_{1}}$ production (blue lines). Regions to the left of the lines are excluded. Results are extracted from the combination of the high and low b-jet $p_{T}$ signal region for the $\tau_\text{lep}\tau_\text{had}$ and $\tau_\text{had}\tau_\text{had}$ channels, neglecting the interference between the LQ non-resonant production processes and the SM processes. However, the interference between the non-resonant LQ production and SM diagrams, such as those from Z+jets, could be non-negligible for events in the low b-jet $p_{T}$ signal region.

The average tetrajet invariant mass distributions in data, along with the fitted background estimates for 0.16 < $\alpha$ < 0.18.

The average tetrajet invariant mass distributions in data, along with the fitted background estimates for 0.18 < $\alpha$ < 0.20.

The average tetrajet invariant mass distributions in data, along with the fitted background estimates for 0.20 < $\alpha$ < 0.22.

The average tetrajet invariant mass distributions in data, along with the fitted background estimates for 0.22 < $\alpha$ < 0.24.

The average tetrajet invariant mass distributions in data, along with the fitted background estimates for 0.24 < $\alpha$ < 0.26.

The average tetrajet invariant mass distributions in data, along with the fitted background estimates for 0.26 < $\alpha$ < 0.28.

The average tetrajet invariant mass distributions in data, along with the fitted background estimates for 0.28 < $\alpha$ < 0.30.

The average tetrajet invariant mass distributions in data, along with the fitted background estimates for 0.30 < $\alpha$ < 0.32.

The average tetrajet invariant mass distributions in data, along with the fitted background estimates for 0.32 < $\alpha$ < 0.34.

The average dijet invariant mass distributions in data, along with the fitted background estimates for 0.10 < $\alpha$ < 0.12.

The average dijet invariant mass distributions in data, along with the fitted background estimates for 0.12 < $\alpha$ < 0.14.

The average dijet invariant mass distributions in data, along with the fitted background estimates for 0.14 < $\alpha$ < 0.16.

The average dijet invariant mass distributions in data, along with the fitted background estimates for 0.16 < $\alpha$ < 0.18.

The average dijet invariant mass distributions in data, along with the fitted background estimates for 0.18 < $\alpha$ < 0.20.

The average dijet invariant mass distributions in data, along with the fitted background estimates for 0.20 < $\alpha$ < 0.22.

The average dijet invariant mass distributions in data, along with the fitted background estimates for 0.22 < $\alpha$ < 0.24.

The average dijet invariant mass distributions in data, along with the fitted background estimates for 0.24 < $\alpha$ < 0.26.

The average dijet invariant mass distributions in data, along with the fitted background estimates for 0.26 < $\alpha$ < 0.28.

The average dijet invariant mass distributions in data, along with the fitted background estimates for 0.28 < $\alpha$ < 0.30.

The average dijet invariant mass distributions in data, along with the fitted background estimates for 0.30 < $\alpha$ < 0.32.

The average dijet invariant mass distributions in data, along with the fitted background estimates for 0.32 < $\alpha$ < 0.34.

The expected and observed limits on $m_{4j}$ for signal templates using a 4-parameter fit function for 0.10 < $\alpha$ < 0.12.

The expected and observed limits on $m_{4j}$ for signal templates using a 4-parameter fit function for 0.12 < $\alpha$ < 0.14.

The expected and observed limits on $m_{4j}$ for signal templates using a 4-parameter fit function for 0.14 < $\alpha$ < 0.16.

The expected and observed limits on $m_{4j}$ for signal templates using a 4-parameter fit function for 0.16 < $\alpha$ < 0.18.

The expected and observed limits on $m_{4j}$ for signal templates using a 4-parameter fit function for 0.18 < $\alpha$ < 0.20.

The expected and observed limits on $m_{4j}$ for signal templates using a 4-parameter fit function for 0.20 < $\alpha$ < 0.22.

The expected and observed limits on $m_{4j}$ for signal templates using a 4-parameter fit function for 0.22 < $\alpha$ < 0.24.

The expected and observed limits on $m_{4j}$ for signal templates using a 4-parameter fit function for 0.24 < $\alpha$ < 0.26.

The expected and observed limits on $m_{4j}$ for signal templates using a 4-parameter fit function for 0.26 < $\alpha$ < 0.28.

The expected and observed limits on $m_{4j}$ for signal templates using a 4-parameter fit function for 0.28 < $\alpha$ < 0.30.

The expected and observed limits on $m_{4j}$ for signal templates using a 4-parameter fit function for 0.30 < $\alpha$ < 0.32.

The expected and observed limits on $m_{4j}$ for signal templates using a 4-parameter fit function for 0.32 < $\alpha$ < 0.34.

The expected and observed limits on $m_{\langle 2j \rangle}$ for signal templates using a 4-parameter fit function for 0.10 < $\alpha$ < 0.12.

The expected and observed limits on $m_{\langle 2j \rangle}$ for signal templates using a 4-parameter fit function for 0.12 < $\alpha$ < 0.14.

The expected and observed limits on $m_{\langle 2j \rangle}$ for signal templates using a 4-parameter fit function for 0.14 < $\alpha$ < 0.16.

The expected and observed limits on $m_{\langle 2j \rangle}$ for signal templates using a 4-parameter fit function for 0.16 < $\alpha$ < 0.18.

The expected and observed limits on $m_{\langle 2j \rangle}$ for signal templates using a 4-parameter fit function for 0.18 < $\alpha$ < 0.20.

The expected and observed limits on $m_{\langle 2j \rangle}$ for signal templates using a 4-parameter fit function for 0.20 < $\alpha$ < 0.22.

The expected and observed limits on $m_{\langle 2j \rangle}$ for signal templates using a 4-parameter fit function for 0.22 < $\alpha$ < 0.24.

The expected and observed limits on $m_{\langle 2j \rangle}$ for signal templates using a 4-parameter fit function for 0.24 < $\alpha$ < 0.26.

The expected and observed limits on $m_{\langle 2j \rangle}$ for signal templates using a 4-parameter fit function for 0.26 < $\alpha$ < 0.28.

The expected and observed limits on $m_{\langle 2j \rangle}$ for signal templates using a 4-parameter fit function for 0.28 < $\alpha$ < 0.30.

The expected and observed limits on $m_{\langle 2j \rangle}$ for signal templates using a 4-parameter fit function for 0.30 < $\alpha$ < 0.32.

The expected and observed limits on $m_{\langle 2j \rangle}$ for signal templates using a 4-parameter fit function for 0.32 < $\alpha$ < 0.34.

The expected and observed limits on $m_{4j}$ for Gaussian signal templates using a 4-parameter fit function for 0.10 < $\alpha$ < 0.12.

The expected and observed limits on $m_{4j}$ for Gaussian signal templates using a 4-parameter fit function for 0.12 < $\alpha$ < 0.14.

The expected and observed limits on $m_{4j}$ for Gaussian signal templates using a 4-parameter fit function for 0.14 < $\alpha$ < 0.16.

The expected and observed limits on $m_{4j}$ for Gaussian signal templates using a 4-parameter fit function for 0.16 < $\alpha$ < 0.18.

The expected and observed limits on $m_{4j}$ for Gaussian signal templates using a 4-parameter fit function for 0.18 < $\alpha$ < 0.20.

The expected and observed limits on $m_{4j}$ for Gaussian signal templates using a 4-parameter fit function for 0.20 < $\alpha$ < 0.22.

The expected and observed limits on $m_{4j}$ for Gaussian signal templates using a 4-parameter fit function for 0.22 < $\alpha$ < 0.24.

The expected and observed limits on $m_{4j}$ for Gaussian signal templates using a 4-parameter fit function for 0.24 < $\alpha$ < 0.26.

The expected and observed limits on $m_{4j}$ for Gaussian signal templates using a 4-parameter fit function for 0.26 < $\alpha$ < 0.28.

The expected and observed limits on $m_{4j}$ for Gaussian signal templates using a 4-parameter fit function for 0.28 < $\alpha$ < 0.30.

The expected and observed limits on $m_{4j}$ for Gaussian signal templates using a 4-parameter fit function for 0.30 < $\alpha$ < 0.32.

The expected and observed limits on $m_{4j}$ for Gaussian signal templates using a 4-parameter fit function for 0.32 < $\alpha$ < 0.34.

The expected and observed limits on $m_{\langle 2j \rangle}$ for Gaussian signal templates using a 4-parameter fit function for 0.10 < $\alpha$ < 0.12.

The expected and observed limits on $m_{\langle 2j \rangle}$ for Gaussian signal templates using a 4-parameter fit function for 0.12 < $\alpha$ < 0.14.

The expected and observed limits on $m_{\langle 2j \rangle}$ for Gaussian signal templates using a 4-parameter fit function for 0.14 < $\alpha$ < 0.16.

The expected and observed limits on $m_{\langle 2j \rangle}$ for Gaussian signal templates using a 4-parameter fit function for 0.16 < $\alpha$ < 0.18.

The expected and observed limits on $m_{\langle 2j \rangle}$ for Gaussian signal templates using a 4-parameter fit function for 0.18 < $\alpha$ < 0.20.

The expected and observed limits on $m_{\langle 2j \rangle}$ for Gaussian signal templates using a 4-parameter fit function for 0.20 < $\alpha$ < 0.22.

The expected and observed limits on $m_{\langle 2j \rangle}$ for Gaussian signal templates using a 4-parameter fit function for 0.22 < $\alpha$ < 0.24.

The expected and observed limits on $m_{\langle 2j \rangle}$ for Gaussian signal templates using a 4-parameter fit function for 0.24 < $\alpha$ < 0.26.

The expected and observed limits on $m_{\langle 2j \rangle}$ for Gaussian signal templates using a 4-parameter fit function for 0.26 < $\alpha$ < 0.28.

The expected and observed limits on $m_{\langle 2j \rangle}$ for Gaussian signal templates using a 4-parameter fit function for 0.28 < $\alpha$ < 0.30.

The expected and observed limits on $m_{\langle 2j \rangle}$ for Gaussian signal templates using a 4-parameter fit function for 0.30 < $\alpha$ < 0.32.

The expected and observed limits on $m_{\langle 2j \rangle}$ for Gaussian signal templates using a 4-parameter fit function for 0.32 < $\alpha$ < 0.34.

The product of the analysis acceptance and selection efficiency for all analysis selection criteria is shown as a function of $m_Y$ and $m_{X}/m_{Y}$.

The cross-section of the signal samples is shown as a function of $m_Y$ and $m_{X}/m_{Y}$.

The expected minus one standard deviation exclusion limit at 95% CL for the clockwork gravity model projected in the $k–M_{5}$ parameter space for the $ee$ channel for the case with mass thresholds.

The expected plus two standard deviation exclusion limit at 95% CL for the clockwork gravity model projected in the $k–M_{5}$ parameter space for the $ee$ channel for the case with mass thresholds.

The expected minus two standard deviation exclusion limit at 95% CL for the clockwork gravity model projected in the $k–M_{5}$ parameter space for the $ee$ channel for the case with mass thresholds.

The observed exclusion limit at 95% CL for the clockwork gravity model projected in the $k–M_{5}$ parameter space for the $\gamma\gamma$ channel for the case with mass thresholds.

The median expected exclusion limit at 95% CL for the clockwork gravity model projected in the $k–M_{5}$ parameter space for the $\gamma\gamma$ channel for the case with mass thresholds.

The expected plus one standard deviation exclusion limit at 95% CL for the clockwork gravity model projected in the $k–M_{5}$ parameter space for the $\gamma\gamma$ channel for the case with mass thresholds.

The expected minus one standard deviation exclusion limit at 95% CL for the clockwork gravity model projected in the $k–M_{5}$ parameter space for the $\gamma\gamma$ channel for the case with mass thresholds.

The expected plus two standard deviation exclusion limit at 95% CL for the clockwork gravity model projected in the $k–M_{5}$ parameter space for the $\gamma\gamma$ channel for the case with mass thresholds.

The expected minus two standard deviation exclusion limit at 95% CL for the clockwork gravity model projected in the $k–M_{5}$ parameter space for the $\gamma\gamma$ channel for the case with mass thresholds.

The observed exclusion limit at 95% CL for the clockwork gravity model projected in the $k–M_{5}$ parameter space for the $ee$ channel for the case without mass thresholds.

The median expected exclusion limit at 95% CL for the clockwork gravity model projected in the $k–M_{5}$ parameter space for the $ee$ channel for the case without mass thresholds.

The expected plus one standard deviation exclusion limit at 95% CL for the clockwork gravity model projected in the $k–M_{5}$ parameter space for the $ee$ channel for the case without mass thresholds.

The expected minus one standard deviation exclusion limit at 95% CL for the clockwork gravity model projected in the $k–M_{5}$ parameter space for the $ee$ channel for the case without mass thresholds.

The expected plus two standard deviation exclusion limit at 95% CL for the clockwork gravity model projected in the $k–M_{5}$ parameter space for the $ee$ channel for the case without mass thresholds.

The expected minus two standard deviation exclusion limit at 95% CL for the clockwork gravity model projected in the $k–M_{5}$ parameter space for the $ee$ channel for the case without mass thresholds.

The observed exclusion limit at 95% CL for the clockwork gravity model projected in the $k–M_{5}$ parameter space for the $\gamma\gamma$ channel for the case without mass thresholds.

The median expected exclusion limit at 95% CL for the clockwork gravity model projected in the $k–M_{5}$ parameter space for the $\gamma\gamma$ channel for the case without mass thresholds.

The expected plus one standard deviation exclusion limit at 95% CL for the clockwork gravity model projected in the $k–M_{5}$ parameter space for the $\gamma\gamma$ channel for the case without mass thresholds.

The expected minus one standard deviation exclusion limit at 95% CL for the clockwork gravity model projected in the $k–M_{5}$ parameter space for the $\gamma\gamma$ channel for the case without mass thresholds.

The expected plus two standard deviation exclusion limit at 95% CL for the clockwork gravity model projected in the $k–M_{5}$ parameter space for the $\gamma\gamma$ channel for the case without mass thresholds.

The expected minus two standard deviation exclusion limit at 95% CL for the clockwork gravity model projected in the $k–M_{5}$ parameter space for the $\gamma\gamma$ channel for the case without mass thresholds.

Background-only fit to the observed numbers of events (black points) as a function of the diphoton invariant mass $m_{\gamma\gamma}$.

Expected and observed 95% CL upper limits on the signal fiducial cross section $\sigma_{\textrm{fid}}$, assuming 100% branching ratio for ALP decay into two photons, as a function of the hypothetical ALP mass $m_{\textrm{X}}$. The 1$\sigma$ and 2$\sigma$ confidence intervals are shown by the coloured bands.

Expected and observed 95% CL upper limits on the the ALP coupling constant, assuming 100% branching ratio for ALP decay into two photons, as a function of the hypothetical ALP mass $m_{\textrm{X}}$. The 1$\sigma$ and 2$\sigma$ confidence intervals are shown by the coloured bands. Contours of the ALP natural width $\Gamma$ are illustrated by the smooth blue solid lines.

Fiducial region acceptance of signal events for coupling constant f$^{-1}=0.05$ TeV$^{-1}$. EL, SD, and DD stand for the exclusive, single-dissociative, and double-dissociative processes. The acceptances are parameterised as $p_0 + p^{i}_{1} e^{p_{2}^{i} m_{\textrm{X}}} + p_{3}^{i} e^{p_{4}^{i} m_{\textrm{X}}} $ , where $i$ specifies the process type, exclusive, single-dissociative, or double-dissociative.

Correction factors $A_{\textrm{X}}/A_{0}$ and $C_{\textrm{X}}$ for the exclusive signal events with ALP coupling constant f$^{-1}$=0.05 TeV$^{-1}$, where acceptance $A_{\textrm{X}}$ and efficiency $C_{\textrm{X}}$ are defined in arXiv:2102.13405. $A_{\textrm{X}}$/$A_{0}$ and $C_{\textrm{X}}$ are parameterised as $p_0 + p^{i}_{1} e^{p_{2}^{i} m_{\textrm{X}}} + p_{3}^{i} e^{p_{4}^{i} m_{\textrm{X}}}$, where $i$ specifies the process type.

Simulated ALP signal cutflow for $m_{\rm{X}}$=300 GeV and $m_{\rm{X}}$=1200 GeV. The yields are normalized to a luminosity of 14.6 fb$^{-1}$. Selections are applied for each side after the selection on acoplanarity. The union of the sets of selected events is taken finally.

Observed (solid line) and expected (dashed line) 95\% CL limits on the product of the production cross section times the branching ratio of a narrow width spin-0 resonance $X$ produced via gluon-gluon initial states decaying to a $Z$ boson and a photon, $\sigma(pp\to X)\cdot\mathcal{B}(X\to Z\gamma)$, as a function of the resonance mass $m_{X}$. Observed (expected) results are derived from ensemble tests for $m_{X}>1850 (900)$~\GeV. The green and yellow solid bands correspond to the $\pm1\sigma$ and $\pm2\sigma$ intervals for the expected upper limits. The limits in the $m_{X}$ ranges from 220 GeV to 3400 GeV and are obtained from the combined $ee\gamma$ and $\mu\mu\gamma$ channels.

Observed (solid line) and expected (dashed line) 95\% CL limits on the product of the production cross section times the branching ratio of a narrow width spin-2 resonance $X$ produced via gluon-gluon initial states decaying to a $Z$ boson and a photon, $\sigma(pp\to X)\cdot\mathcal{B}(X\to Z\gamma)$, as a function of the resonance mass $m_{X}$. Observed (expected) results are derived from ensemble tests for $m_{X}>1850 (900)$~\GeV. The green and yellow solid bands correspond to the $\pm1\sigma$ and $\pm2\sigma$ intervals for the expected upper limit respectively. The limits in the $m_{X}$ ranges from 220 GeV to 3400 GeV and are obtained from the combined $ee\gamma$ and $\mu\mu\gamma$ channels..

Observed (solid line) and expected (dashed line) 95\% CL limits on the product of the production cross section times the branching ratio of a narrow width spin-2 resonance $X$ produced via $q\bar{q}$ initial states decaying to a $Z$ boson and a photon, $\sigma(pp\to X)\cdot\mathcal{B}(X\to Z\gamma)$, as a function of the resonance mass $m_{X}$. Observed (expected) results are derived from ensemble tests for $m_{X}>1850 (900)$~\GeV. The green and yellow solid bands correspond to the $\pm1\sigma$ and $\pm2\sigma$ intervals for the expected upper limit respectively. The limits in the $m_{X}$ ranges from 220 GeV to 3400 GeV and are obtained from the combined $ee\gamma$ and $\mu\mu\gamma$ channels.

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

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

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

Distribution (events/100 GeV) of the reconstructed $m_{tb}$ for data and backgrounds in the 0-lepton channel's validation region before the fit to data. The systematics uncertainty is shown for the pre-fit background sum, including the background statistical uncertainty. The individual background components are obtained before the fit, too. There is also the expected signal distribution. The distribution of the $W^{\prime}$ boson signal for a mass of 3 TeV is normalised to the predicted cross-section. The last bin in each distribution includes overflow.

Reconstructed $m_{tb}$ distributions (events/50 GeV) for data and backgrounds in the 1-lepton channel's 2j1b signal region after the background-only fit to data. The systematics uncertainty is shown for the post-fit background sum, including the background statistical uncertainty. The individual background components are obtained after the fit, too. There are also the pre-fit background sum and the expected signal distribution. The distribution of the $W^{\prime}$ boson signal for a mass of 3 TeV is normalised to the predicted cross-section. The last bin in each distribution includes overflow.

Reconstructed $m_{tb}$ distributions (events/50 GeV) for data and backgrounds in the 1-lepton channel's 3j1b signal region after the background-only fit to data. The systematics uncertainty is shown for the post-fit background sum, including the background statistical uncertainty. The individual background components are obtained after the fit, too. There are also the pre-fit background sum and the expected signal distribution. The distribution of the $W^{\prime}$ boson signal for a mass of 3 TeV is normalised to the predicted cross-section. The last bin in each distribution includes overflow.

Reconstructed $m_{tb}$ distributions (events/50 GeV) for data and backgrounds in the 1-lepton channel's 2j2b signal region after the background-only fit to data. The systematics uncertainty is shown for the post-fit background sum, including the background statistical uncertainty. The individual background components are obtained after the fit, too. There are also the pre-fit background sum and the expected signal distribution. The distribution of the $W^{\prime}$ boson signal for a mass of 3 TeV is normalised to the predicted cross-section. The last bin in each distribution includes overflow.

Reconstructed $m_{tb}$ distributions (events/50 GeV) for data and backgrounds in the 1-lepton channel's 3j2b signal region after the background-only fit to data. The systematics uncertainty is shown for the post-fit background sum, including the background statistical uncertainty. The individual background components are obtained after the fit, too. There are also the pre-fit background sum and the expected signal distribution. The distribution of the $W^{\prime}$ boson signal for a mass of 3 TeV is normalised to the predicted cross-section. The last bin in each distribution includes overflow.

Reconstructed $m_{tb}$ distributions (events/50 GeV) for data and backgrounds in the 1-lepton channel's 2j1b control region after the background-only fit to data. The systematics uncertainty is shown for the post-fit background sum, including the background statistical uncertainty. The individual background components are obtained after the fit, too. There are also the pre-fit background sum and the expected signal distribution. The distribution of the $W^{\prime}$ boson signal for a mass of 3 TeV is normalised to the predicted cross-section. The last bin in each distribution includes overflow.

Reconstructed $m_{tb}$ distributions (events/50 GeV) for data and backgrounds in the 1-lepton channel's 3j1b control region after the background-only fit to data. The systematics uncertainty is shown for the post-fit background sum, including the background statistical uncertainty. The individual background components are obtained after the fit, too. There are also the pre-fit background sum and the expected signal distribution. The distribution of the $W^{\prime}$ boson signal for a mass of 3 TeV is normalised to the predicted cross-section. The last bin in each distribution includes overflow.

Reconstructed $m_{tb}$ distributions (events/50 GeV) for data and backgrounds in the 1-lepton channel's 2j1b validation region after the background-only fit to data. The systematics uncertainty is shown for the post-fit background sum, including the background statistical uncertainty. The individual background components are obtained after the fit, too. There are also the pre-fit background sum and the expected signal distribution. The distribution of the $W^{\prime}$ boson signal for a mass of 3 TeV is normalised to the predicted cross-section. The last bin in each distribution includes overflow.

Reconstructed $m_{tb}$ distributions (events/50 GeV) for data and backgrounds in the 1-lepton channel's 3j1b validation region after the background-only fit to data. The systematics uncertainty is shown for the post-fit background sum, including the background statistical uncertainty. The individual background components are obtained after the fit, too. There are also the pre-fit background sum and the expected signal distribution. The distribution of the $W^{\prime}$ boson signal for a mass of 3 TeV is normalised to the predicted cross-section. The last bin in each distribution includes overflow.

Observed and expected 95% CL limits on the cross-section times branching ratio for the production of a $W^{\prime}$ boson with decay to tb and left-handed couplings as a function of the mass of the $W^{\prime}$ boson and a coupling value of $g^{\prime}/g$=0.5. The observed and expected limits are derived using a linear interpolation between simulated signal mass hypotheses. The uncertainty on the theory prediction includes components from factorisation and renormalisation scales, PDFs, the strong coupling constant, and the top-quark mass.

Observed and expected 95% CL limits on the cross-section times branching ratio for the production of a $W^{\prime}$ boson with decay to tb and left-handed couplings as a function of the mass of the $W^{\prime}$ boson and a coupling value of $g^{\prime}/g$=1.0. The observed and expected limits are derived using a linear interpolation between simulated signal mass hypotheses. The uncertainty on the theory prediction includes components from factorisation and renormalisation scales, PDFs, the strong coupling constant, and the top-quark mass.

Observed and expected 95% CL limits on the cross-section times branching ratio for the production of a $W^{\prime}$ boson with decay to tb and left-handed couplings as a function of the mass of the $W^{\prime}$ boson and a coupling value of $g^{\prime}/g$=2.0. The observed and expected limits are derived using a linear interpolation between simulated signal mass hypotheses. The uncertainty on the theory prediction includes components from factorisation and renormalisation scales, PDFs, the strong coupling constant, and the top-quark mass.

Observed and expected 95% CL limits on the cross-section times branching ratio for the production of a $W^{\prime}$ boson with decay to tb and right-handed couplings as a function of the mass of the $W^{\prime}$ boson and a coupling value of $g^{\prime}/g$=0.5. The observed and expected limits are derived using a linear interpolation between simulated signal mass hypotheses. The uncertainty on the theory prediction includes components from factorisation and renormalisation scales, PDFs, the strong coupling constant, and the top-quark mass.

Observed and expected 95% CL limits on the cross-section times branching ratio for the production of a $W^{\prime}$ boson with decay to tb and right-handed couplings as a function of the mass of the $W^{\prime}$ boson and a coupling value of $g^{\prime}/g$=1.0. The observed and expected limits are derived using a linear interpolation between simulated signal mass hypotheses. The uncertainty on the theory prediction includes components from factorisation and renormalisation scales, PDFs, the strong coupling constant, and the top-quark mass.

Observed and expected 95% CL limits on the cross-section times branching ratio for the production of a $W^{\prime}$ boson with decay to tb and right-handed couplings as a function of the mass of the $W^{\prime}$ boson and a coupling value of $g^{\prime}/g$=2.0. The observed and expected limits are derived using a linear interpolation between simulated signal mass hypotheses. The uncertainty on the theory prediction includes components from factorisation and renormalisation scales, PDFs, the strong coupling constant, and the top-quark mass.

95% CL limits on the coupling value ($g^{\prime}/g$) and the $W^{\prime}$-boson mass for left-handed $W^{\prime}$-boson couplings. The area above the line is excluded.

95% CL limits on the coupling value ($g^{\prime}/g$) and the $W^{\prime}$-boson mass for right-handed $W^{\prime}$-boson couplings. The area above the line is excluded.

The product of fiducial acceptance and selection efficiency for the three signal regions of the 0-lepton channel as a function of the mass of the $W^{\prime}$ boson with left-handed chirality. The $W^{\prime}$ boson coupling strength is set to $g^{\prime}/g=1$.

The product of fiducial acceptance and selection efficiency for the three signal regions of the 0-lepton channel as a function of the mass of the $W^{\prime}$ boson with left-handed chirality. The $W^{\prime}$ boson coupling strength is set to $g^{\prime}/g=0.5$.

The product of fiducial acceptance and selection efficiency for the three signal regions of the 0-lepton channel as a function of the mass of the $W^{\prime}$ boson with left-handed chirality. The $W^{\prime}$ boson coupling strength is set to $g^{\prime}/g=2$.

The product of fiducial acceptance and selection efficiency for the four signal regions of the 1-lepton channel as a function of the mass of the $W^{\prime}$ boson with left-handed chirality. The $W^{\prime}$ boson coupling strength is set to $g^{\prime}/g=1$.

The product of fiducial acceptance and selection efficiency for the four signal regions of the 1-lepton channel as a function of the mass of the $W^{\prime}$ boson with left-handed chirality. The $W^{\prime}$ boson coupling strength is set to $g^{\prime}/g=0.5$.

The product of fiducial acceptance and selection efficiency for the four signal regions of the 1-lepton channel as a function of the mass of the $W^{\prime}$ boson with left-handed chirality. The $W^{\prime}$ boson coupling strength is set to $g^{\prime}/g=2$.

The product of fiducial acceptance and selection efficiency for the three signal regions of the 0-lepton channel as a function of the mass of the $W^{\prime}$ boson with right-handed chirality. The $W^{\prime}$ boson coupling strength is set to $g^{\prime}/g=1$.

The product of fiducial acceptance and selection efficiency for the three signal regions of the 0-lepton channel as a function of the mass of the $W^{\prime}$ boson with right-handed chirality. The $W^{\prime}$ boson coupling strength is set to $g^{\prime}/g=0.5$.

The product of fiducial acceptance and selection efficiency for the three signal regions of the 0-lepton channel as a function of the mass of the $W^{\prime}$ boson with right-handed chirality. The $W^{\prime}$ boson coupling strength is set to $g^{\prime}/g=2$.

The product of fiducial acceptance and selection efficiency for the four signal regions of the 1-lepton channel as a function of the mass of the $W^{\prime}$ boson with right-handed chirality. The $W^{\prime}$ boson coupling strength is set to $g^{\prime}/g=1$.

The product of fiducial acceptance and selection efficiency for the four signal regions of the 1-lepton channel as a function of the mass of the $W^{\prime}$ boson with right-handed chirality. The $W^{\prime}$ boson coupling strength is set to $g^{\prime}/g=0.5$.

The product of fiducial acceptance and selection efficiency for the four signal regions of the 1-lepton channel as a function of the mass of the $W^{\prime}$ boson with right-handed chirality. The $W^{\prime}$ boson coupling strength is set to $g^{\prime}/g=2$.