Distribution of the diphoton invariant mass in validation region VR2. The solid histograms are stacked to show the SM expectations after the 2×2D background estimation technique is applied. Background and signal predictions are normalised to the luminosity. The background category "h (other)" includes events originating from VBF, Vh, ggF, thq, thW and bb̄h, all subdominant in this signature. Statistical and systematic uncertainties are indicated by the shaded area. The lower panel of each plot shows the ratio of the data to the SM prediction for the respective bin. The first and last bins include the underflows and overflows respectively.

Distribution of the missing transverse momentum in validation region VR2. The solid histograms are stacked to show the SM expectations after the 2×2D background estimation technique is applied. Background and signal predictions are normalised to the luminosity. The background category "h (other)" includes events originating from VBF, Vh, ggF, thq, thW and bb̄h, all subdominant in this signature. Statistical and systematic uncertainties are indicated by the shaded area. The lower panel of each plot shows the ratio of the data to the SM prediction for the respective bin. The first and last bins include the underflows and overflows respectively.

Distribution of the diphoton invariant mass with all selections of the signal regions applied, except on m_{γγ} itself, for signal region SR1h. The background estimation techniques described in the text are applied. The different backgrounds are stacked to add up to the total SM prediction in each bin. The predicted yields for signal benchmark models of varying χ̃⁰₁ mass are also overlaid (not stacked) assuming B(χ̃⁰₁ → hG̃ ) to equal 100%. Background and signal predictions are normalised to the luminosity. The background category "h (other)" includes events originating from VBF, Vh, ggF, thq, thW and bb̄h, all subdominant in this signature. The sizes of the statistical and systematic uncertainties are indicated by the shaded areas. The lower panels show the ratio of the data to the SM prediction. Arrows indicate the borders of the signal region (|m_{γγ}-125 GeV|<5 GeV). The first and last bins include the underflows and overflows respectively.

Distribution of the diphoton invariant mass with all selections of the signal regions applied, except on m_{γγ} itself, for signal region SR1Z. The background estimation techniques described in the text are applied. The different backgrounds are stacked to add up to the total SM prediction in each bin. The predicted yields for signal benchmark models of varying χ̃⁰₁ mass are also overlaid (not stacked) assuming B(χ̃⁰₁ → hG̃ ) to equal 50%. Background and signal predictions are normalised to the luminosity. The background category "h (other)" includes events originating from VBF, Vh, ggF, thq, thW and bb̄h, all subdominant in this signature. The sizes of the statistical and systematic uncertainties are indicated by the shaded areas. The lower panels show the ratio of the data to the SM prediction. Arrows indicate the borders of the signal region (|m_{γγ}-125 GeV|<5 GeV). The first and last bins include the underflows and overflows respectively.

Distribution of the diphoton invariant mass with all selections of the signal regions applied, except on m_{γγ} itself, for signal region SR2. The background estimation techniques described in the text are applied. The different backgrounds are stacked to add up to the total SM prediction in each bin. The predicted yields for signal benchmark models of varying χ̃⁰₁ mass are also overlaid (not stacked) assuming B(χ̃⁰₁ → hG̃ ) to equal 100%. Background and signal predictions are normalised to the luminosity. The background category "h (other)" includes events originating from VBF, Vh, ggF, thq, thW and bb̄h, all subdominant in this signature. The sizes of the statistical and systematic uncertainties are indicated by the shaded areas. The lower panels show the ratio of the data to the SM prediction. Arrows indicate the borders of the signal region (|m_{γγ}-125 GeV|<5 GeV). The first and last bins include the underflows and overflows respectively.

Observed and expected limits on the pure higgsino cross-section at 95% CL assuming B(χ̃⁰₁ → hG̃ )=100% for different χ̃⁰₁ masses, obtained by a statistical combination of the three signal regions SR1h, SR1Z and SR2. The inner and outer bands indicate the 1σ and 2σ variation on the expected limit respectively.

Observed and expected 95% CL limits on the pure-higgsino branching fraction to B(χ̃⁰₁ → hG̃ ) as a function of the higgsino mass m(χ̃⁰₁) assuming it decays via either χ̃⁰₁→ hG̃ or χ̃⁰₁→ ZG̃. Limits are obtained by performing a statistical combination of the three signal regions SR1h, SR1Z and SR2. The ± 1σ variation on the expected limit is shown. The dotted lines indicate the observed limit obtained by a variation of theoretical prediction for the neutralino production cross-section by ±1 σ. Values of B(χ̃⁰₁ → hG̃ ) larger than the observed 95% CL limit are excluded, as indicated by the hatched area.

Observed and expected 95% CL limits on the pure-higgsino branching fraction to B(χ̃⁰₁ → hG̃ ) as a function of the higgsino mass m(χ̃⁰₁) assuming it decays via either χ̃⁰₁→ hG̃ or χ̃⁰₁→ ZG̃. Limits are obtained by performing a statistical combination of the three signal regions SR1h, SR1Z and SR2. The ± 1σ variation on the expected limit is shown. The dotted lines indicate the observed limit obtained by a variation of theoretical prediction for the neutralino production cross-section by ±1 σ. Values of B(χ̃⁰₁ → hG̃ ) larger than the observed 95% CL limit are excluded, as indicated by the hatched area.

Observed and expected 95% CL limits on the pure-higgsino branching fraction to B(χ̃⁰₁ → hG̃ ) as a function of the higgsino mass m(χ̃⁰₁) assuming it decays via either χ̃⁰₁→ hG̃ or χ̃⁰₁→ ZG̃. Limits are obtained by performing a statistical combination of the three signal regions SR1h, SR1Z and SR2. The ± 1σ variation on the expected limit is shown. The dotted lines indicate the observed limit obtained by a variation of theoretical prediction for the neutralino production cross-section by ±1 σ. Values of B(χ̃⁰₁ → hG̃ ) larger than the observed 95% CL limit are excluded, as indicated by the hatched area.

Observed and expected 95% CL limits on the pure-higgsino branching fraction to B(χ̃⁰₁ → hG̃ ) as a function of the higgsino mass m(χ̃⁰₁) assuming it decays via either χ̃⁰₁→ hG̃ or χ̃⁰₁→ ZG̃. Limits are obtained by performing a statistical combination of the three signal regions SR1h, SR1Z and SR2. The ± 1σ variation on the expected limit is shown. The dotted lines indicate the observed limit obtained by a variation of theoretical prediction for the neutralino production cross-section by ±1 σ. Values of B(χ̃⁰₁ → hG̃ ) larger than the observed 95% CL limit are excluded, as indicated by the hatched area.

Observed and expected 95% CL limits on the pure-higgsino branching fraction to B(χ̃⁰₁ → hG̃ ) as a function of the higgsino mass m(χ̃⁰₁) assuming it decays via either χ̃⁰₁→ hG̃ or χ̃⁰₁→ ZG̃. Limits are obtained by performing a statistical combination of the three signal regions SR1h, SR1Z and SR2. The ± 1σ variation on the expected limit is shown. The dotted lines indicate the observed limit obtained by a variation of theoretical prediction for the neutralino production cross-section by ±1 σ. Values of B(χ̃⁰₁ → hG̃ ) larger than the observed 95% CL limit are excluded, as indicated by the hatched area.

Observed and expected 95% CL limits on the pure-higgsino branching fraction to B(χ̃⁰₁ → hG̃ ) as a function of the higgsino mass m(χ̃⁰₁) assuming it decays via either χ̃⁰₁→ hG̃ or χ̃⁰₁→ ZG̃. Limits are obtained by performing a statistical combination of the three signal regions SR1h, SR1Z and SR2. The ± 1σ variation on the expected limit is shown. The dotted lines indicate the observed limit obtained by a variation of theoretical prediction for the neutralino production cross-section by ±1 σ. Values of B(χ̃⁰₁ → hG̃ ) larger than the observed 95% CL limit are excluded, as indicated by the hatched area.

Numbers of signal and background events in the signal regions. The respective background estimation techniques are applied. The background category "h (other)" includes events originating from VBF, Vh, ggF, thq, thW and bb&#772;h, all subdominant in this signature. The different backgrounds are stacked to add up to the total Standard Model prediction in each bin. The predicted yields for signal benchmark models of varying &chi;&#771;⁰₁ mass are also plotted (not stacked), assuming B(&chi;&#771;⁰₁ &rarr; hG&#771; )=100% and a &chi;&#771;⁰₁ mass of 130 or 200 GeV. The statistical and systematic uncertainties are indicated by the shaded areas in the top plot. The bottom panel shows the statistical significance&nbsp;<a href="https://atlas.web.cern.ch/Atlas/GROUPS/PHYSICS/PUBNOTES/ATL-PHYS-PUB-2020-025/">[Ref]</a> of the difference between the SM prediction and the observed data in each region.

Observed 95% CL limits in pb on the pure higgsino cross-section, shown in the m(&chi;&#771;⁰₁)-B(&chi;&#771;⁰₁ &rarr; hG&#771; ) plane. Limits are obtained by a statistical combination of the three signal regions SR1h, SR1Z and SR2, assuming the neutralino to decay via either &chi;&#771;⁰₁&rarr; hG&#771; or &chi;&#771;⁰₁&rarr; ZG&#771;.

Acceptances for all signal model points considered in the analysis, shown in the m(&chi;&#771;⁰₁)-B(&chi;&#771;⁰₁ &rarr; hG&#771; ) plane. Acceptances are provided separately for signal region SR1h, assuming the neutralino to decay via either &chi;&#771;⁰₁&rarr; hG&#771; or &chi;&#771;⁰₁&rarr; ZG&#771;.

Acceptances for all signal model points considered in the analysis, shown in the m(&chi;&#771;⁰₁)-B(&chi;&#771;⁰₁ &rarr; hG&#771; ) plane. Acceptances are provided separately for signal region SR1Z, assuming the neutralino to decay via either &chi;&#771;⁰₁&rarr; hG&#771; or &chi;&#771;⁰₁&rarr; ZG&#771;.

Acceptances for all signal model points considered in the analysis, shown in the m(&chi;&#771;⁰₁)-B(&chi;&#771;⁰₁ &rarr; hG&#771; ) plane. Acceptances are provided separately for signal region SR2, assuming the neutralino to decay via either &chi;&#771;⁰₁&rarr; hG&#771; or &chi;&#771;⁰₁&rarr; ZG&#771;.

Efficiencies for all signal model points considered in the analysis, shown in the m(&chi;&#771;⁰₁)-B(&chi;&#771;⁰₁ &rarr; hG&#771; ) plane. Efficiencies are provided separately for signal region SR1h, assuming the neutralino to decay via either &chi;&#771;⁰₁&rarr; hG&#771; or &chi;&#771;⁰₁&rarr; ZG&#771;.

Efficiencies for all signal model points considered in the analysis, shown in the m(&chi;&#771;⁰₁)-B(&chi;&#771;⁰₁ &rarr; hG&#771; ) plane. Efficiencies are provided separately for signal region SR1Z, assuming the neutralino to decay via either &chi;&#771;⁰₁&rarr; hG&#771; or &chi;&#771;⁰₁&rarr; ZG&#771;.

Efficiencies for all signal model points considered in the analysis, shown in the m(&chi;&#771;⁰₁)-B(&chi;&#771;⁰₁ &rarr; hG&#771; ) plane. Efficiencies are provided separately for signal region SR2, assuming the neutralino to decay via either &chi;&#771;⁰₁&rarr; hG&#771; or &chi;&#771;⁰₁&rarr; ZG&#771;.

Observed and expected numbers of events in the three signal regions. The background category "h (other)" includes events originating from VBF, Vh, ggF, thq, thW and bb&#772;h, all subdominant in this signature. The table also includes model-independent 95% CL upper limits on the visible number of BSM events (S⁹⁵_obs), the number of BSM events given the expected number of background events (S⁹⁵_exp) and the visible BSM cross-section (&lang;&epsilon; &sigma;&rang;_obs⁹⁵), all calculated from pseudo-experiments. The discovery p-value (p₀) is also shown and its value is capped at 0.5 if the observed number of events is below the expected number of events.

Cutflows of two benchmark signal points assuming B(&chi;&#771;⁰₁ &rarr; hG&#771; )=100% for all three discovery signal regions. The initial selection includes the leptons veto. Only statistical uncertainties are included. Expected yields are normalised to a luminosity of 139 fb^-1.

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.

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

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

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

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

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

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

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

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

Expected CLs values in the High mA SR, mA vs ma signal grid.

Observed CLs values in the High mA SR, mA vs ma signal grid.

Expected CLs values in the Low mA SR, mA vs ma signal grid.

Observed CLs values in the High mA SR, mA vs ma signal grid.

CLs+b values in the Low mA SR, mA vs tanB signal grid.

CLs+b values in the High mA SR, mA vs ma signal grid.

CLs+b values in the Low mA SR, mA vs ma signal grid.

Cut flow of the 2HDM+a signal points, gluon–gluon fusion production, Low mA SR. tanB = 1, $sin\theta$ = 0.35. The first two entries in the tables are number of raw MC events, third entry is theoretical prediction, and all other lines include the correct weights. Note that during the generation Higgs boson’s branching ratio to taus has been set to 1. An additional factor of 0.0627 is used to account for that, starting from the ‘Initial’ entry.

Cut flow of the 2HDM+a signal points, gluon–gluon fusion production, High mA SR. tanB = 1, $sin\theta$ = 0.35. The first two entries in the tables are number of raw MC events, third entry is theoretical prediction, and all other lines include the correct weights. Note that during the generation Higgs boson’s branching ratio to taus has been set to 1. An additional factor of 0.0627 is used to account for that, starting from the ‘Initial’ entry.

Cut flow of the 2HDM+a signal points, bb annihilation production, Low mA SR. tanB = 1, $sin\theta$ = 0.35. The first two entries in the tables are number of raw MC events, third entry is theoretical prediction, and all other lines include the correct weights. Note that during the generation Higgs boson’s branching ratio to taus has been set to 1. An additional factor of 0.0627 is used to account for that, starting from the ‘Initial’ entry.

Cut flow of the 2HDM+a signal points, bb annihilation production, High mA SR. tanB = 1, $sin\theta$ = 0.35. The first two entries in the tables are number of raw MC events, third entry is theoretical prediction, and all other lines include the correct weights. Note that during the generation Higgs boson’s branching ratio to taus has been set to 1. An additional factor of 0.0627 is used to account for that, starting from the ‘Initial’ entry.

A comparison of the observed and expected yields in the four bins of the Low mA SR.

A comparison of the observed and expected yields in the two bins of the High mA SR.

Expected exclusion contours at 95% CL as a function of mA and ma.

Observed exclusion contours at 95% CL as a function of mA and ma.

Expected +- 1sigma exclusion contours at 95% CL as a function of mA and ma.

Expected exclusion contours at 95% CL as a function of mA and ma.

Observed exclusion contours at 95% CL as a function of mA and ma.

Expected +- 1sigma exclusion contours at 95% CL as a function of mA and ma.

Acceptance, High mA SR, mA vs tanB grid, 400-750 GeV, bb prod

Acceptance, High mA SR, mA vs tanB grid, >750 GeV, bb prod

Acceptance, Low mA SR, mA vs tanB grid, 100-250 GeV, bb prod

Acceptance, Low mA SR, mA vs tanB grid, 250-400 GeV, bb prod

Acceptance, Low mA SR, mA vs tanB grid, 400-550 GeV, bb prod

Acceptance, Low mA SR, mA vs tanB grid, >550 GeV, bb prod

Acceptance, High mA SR, mA vs ma grid, 400-750 GeV, bb prod

Acceptance, High mA SR, mA vs ma grid, >750 GeV, bb prod

Acceptance, Low mA SR, mA vs ma grid, 100-250 GeV, bb prod

Acceptance, Low mA SR, mA vs ma grid, 250-400 GeV, bb prod

Acceptance, Low mA SR, mA vs ma grid, 400-550 GeV, bb prod

Acceptance, Low mA SR, mA vs ma grid, >550 GeV, bb prod

Acceptance, High mA SR, mA vs tanB grid, 400-750 GeV, ggF prod

Acceptance, High mA SR, mA vs tanB grid, >750 GeV, ggF prod

Acceptance, Low mA SR, mA vs tanB grid, 100-250 GeV, ggF prod

Acceptance, Low mA SR, mA vs tanB grid, 250-400 GeV, ggF prod

Acceptance, Low mA SR, mA vs tanB grid, 400-550 GeV, ggF prod

Acceptance, Low mA SR, mA vs tanB grid, >550 GeV, ggF prod

Acceptance, High mA SR, mA vs ma grid, 400-750 GeV, ggF prod

Acceptance, High mA SR, mA vs ma grid, >750 GeV, ggF prod

Acceptance, Low mA SR, mA vs ma grid, 100-250 GeV, ggF prod

Acceptance, Low mA SR, mA vs ma grid, 250-400 GeV, ggF prod

Acceptance, Low mA SR, mA vs ma grid, 400-550 GeV, ggF prod

Acceptance, Low mA SR, mA vs ma grid, >550 GeV, ggF prod

Efficiency, High mA SR, mA vs tanB grid, 400-750 GeV, bb prod

Efficiency, High mA SR, mA vs tanB grid, >750 GeV, bb prod

Efficiency, Low mA SR, mA vs tanB grid, 100-250 GeV, bb prod

Efficiency, Low mA SR, mA vs tanB grid, 250-400 GeV, bb prod

Efficiency, Low mA SR, mA vs tanB grid, 400-550 GeV, bb prod

Efficiency, Low mA SR, mA vs tanB grid, >550 GeV, bb prod

Efficiency, High mA SR, mA vs ma grid, 400-750 GeV, bb prod

Efficiency, High mA SR, mA vs ma grid, >750 GeV, bb prod

Efficiency, Low mA SR, mA vs ma grid, 100-250 GeV, bb prod

Efficiency, Low mA SR, mA vs ma grid, 250-400 GeV, bb prod

Efficiency, Low mA SR, mA vs ma grid, 400-550 GeV, bb prod

Efficiency, Low mA SR, mA vs ma grid, >550 GeV, bb prod

Efficiency, High mA SR, mA vs tanB grid, 400-750 GeV, ggF prod

Efficiency, High mA SR, mA vs tanB grid, >750 GeV, ggF prod

Efficiency, Low mA SR, mA vs tanB grid, 100-250 GeV, ggF prod

Efficiency, Low mA SR, mA vs tanB grid, 250-400 GeV, ggF prod

Efficiency, Low mA SR, mA vs tanB grid, 400-550 GeV, ggF prod

Efficiency, Low mA SR, mA vs tanB grid, >550 GeV, ggF prod

Efficiency, High mA SR, mA vs ma grid, 400-750 GeV, ggF prod

Efficiency, High mA SR, mA vs ma grid, >750 GeV, ggF prod

Efficiency, Low mA SR, mA vs ma grid, 100-250 GeV, ggF prod

Efficiency, Low mA SR, mA vs ma grid, 250-400 GeV, ggF prod

Efficiency, Low mA SR, mA vs ma grid, 400-550 GeV, ggF prod

Efficiency, Low mA SR, mA vs ma grid, >550 GeV, ggF prod

Fit results of the simultaneous measurements of the $H\to e\tau$ and $H\to \mu\tau$ signals (2POI) showing best-fit values of the LFV branching ratios of the Higgs boson $\hat{B}$($H\to e\tau$). The results from standalone channel/categories fits are compared with the results of the combined fit.

Fit results of the simultaneous measurements of the $H\to e\tau$ and $H\to \mu\tau$ signals (2POI) showing upper limits at 95% C.L. on the LFV branching ratios of the Higgs boson $H\to \mu\tau$. The results from standalone channel/categories fits are compared with the results of the combined fit.

Fit results of the simultaneous measurements of the $H\to e\tau$ and $H\to \mu\tau$ signals (2POI) showing upper limits at 95% C.L. on the LFV branching ratios of the Higgs boson $H\to \mu\tau$. The results from standalone channel/categories fits are compared with the results of the combined fit.

Fit results of the simultaneous measurements of the $H\to e\tau$ and $H\to \mu\tau$ signals (2POI) showing best-fit values of the LFV branching ratios of the Higgs boson $\hat{B}$($H\to \mu\tau$). The results from standalone channel/categories fits are compared with the results of the combined fit.

Fit results of the simultaneous measurements of the $H\to e\tau$ and $H\to \mu\tau$ signals (2POI) showing best-fit values of the LFV branching ratios of the Higgs boson $\hat{B}$($H\to \mu\tau$). The results from standalone channel/categories fits are compared with the results of the combined fit.

Fit results of the independent searches (1 POI) showing upper limits at 95% C.L. on the LFV branching ratios of the Higgs boson $H\to e\tau$. The results from standalone channel/categories fits are compared with the results of the combined fit.

Fit results of the independent searches (1 POI) showing upper limits at 95% C.L. on the LFV branching ratios of the Higgs boson $H\to e\tau$. The results from standalone channel/categories fits are compared with the results of the combined fit.

Fit results of the independent searches (1 POI) showing best-fit values of the LFV branching ratios of the Higgs boson $\hat{B}$($H\to e\tau$). The results from standalone channel/categories fits are compared with the results of the combined fit.

Fit results of the independent searches (1 POI) showing best-fit values of the LFV branching ratios of the Higgs boson $\hat{B}$($H\to e\tau$). The results from standalone channel/categories fits are compared with the results of the combined fit.

Fit results of the independent searches (1 POI) showing upper limits at 95% C.L. on the LFV branching ratios of the Higgs boson $H\to \mu\tau$. The results from standalone channel/categories fits are compared with the results of the combined fit.

Fit results of the independent searches (1 POI) showing upper limits at 95% C.L. on the LFV branching ratios of the Higgs boson $H\to \mu\tau$. The results from standalone channel/categories fits are compared with the results of the combined fit.

Fit results of the independent searches (1 POI) showing best-fit values of the LFV branching ratios of the Higgs boson $\hat{B}$($H\to \mu\tau$). The results from standalone channel/categories fits are compared with the results of the combined fit.

Fit results of the independent searches (1 POI) showing best-fit values of the LFV branching ratios of the Higgs boson $\hat{B}$($H\to \mu\tau$). The results from standalone channel/categories fits are compared with the results of the combined fit.

BDT score distribution, after the simultaneous fit of the $H\to e\tau$ and $H\to \mu\tau$ signals, obtained by fitting the data of the MC-template $\ell\tau_{\ell'}$ channel for $e\tau_{\mu}$ final state in the non-VBF category.

BDT score distribution, after the simultaneous fit of the $H\to e\tau$ and $H\to \mu\tau$ signals, obtained by fitting the data of the MC-template $\ell\tau_{\ell'}$ channel for $e\tau_{\mu}$ final state in the non-VBF category.

BDT score distribution, after the simultaneous fit of the $H\to e\tau$ and $H\to \mu\tau$ signals, obtained by fitting the data of the MC-template $\ell\tau_{\ell'}$ channel for $\mu\tau_{e}$ final state in the non-VBF category.

BDT score distribution, after the simultaneous fit of the $H\to e\tau$ and $H\to \mu\tau$ signals, obtained by fitting the data of the MC-template $\ell\tau_{\ell'}$ channel for $\mu\tau_{e}$ final state in the non-VBF category.

BDT score distribution, after the simultaneous fit of the $H\to e\tau$ and $H\to \mu\tau$ signals, obtained by fitting the data of the MC-template $\ell\tau_{\ell'}$ channel for $e\tau_{\mu}$ final state in the VBF category.

BDT score distribution, after the simultaneous fit of the $H\to e\tau$ and $H\to \mu\tau$ signals, obtained by fitting the data of the MC-template $\ell\tau_{\ell'}$ channel for $e\tau_{\mu}$ final state in the VBF category.

BDT score distribution, after the simultaneous fit of the $H\to e\tau$ and $H\to \mu\tau$ signals, obtained by fitting the data of the MC-template $\ell\tau_{\ell'}$ channel for $\mu\tau_{e}$ final state in the VBF category.

BDT score distribution, after the simultaneous fit of the $H\to e\tau$ and $H\to \mu\tau$ signals, obtained by fitting the data of the MC-template $\ell\tau_{\ell'}$ channel for $\mu\tau_{e}$ final state in the VBF category.

BDT score distribution, after the simultaneous fit of the $H\to e\tau$ and $H\to \mu\tau$ signals, obtained by fitting the data of the MC-template $\ell\tau_{had}$ channel for $e\tau_{had}$ final state in the non-VBF category.

BDT score distribution, after the simultaneous fit of the $H\to e\tau$ and $H\to \mu\tau$ signals, obtained by fitting the data of the MC-template $\ell\tau_{had}$ channel for $e\tau_{had}$ final state in the non-VBF category.

BDT score distribution, after the simultaneous fit of the $H\to e\tau$ and $H\to \mu\tau$ signals, obtained by fitting the data of the MC-template $\ell\tau_{had}$ channel for $\mu\tau_{had}$ final state in the non-VBF category.

BDT score distribution, after the simultaneous fit of the $H\to e\tau$ and $H\to \mu\tau$ signals, obtained by fitting the data of the MC-template $\ell\tau_{had}$ channel for $\mu\tau_{had}$ final state in the non-VBF category.

BDT score distribution, after the simultaneous fit of the $H\to e\tau$ and $H\to \mu\tau$ signals, obtained by fitting the data of the MC-template $\ell\tau_{had}$ channel for $e\tau_{had}$ final state in the VBF category.

BDT score distribution, after the simultaneous fit of the $H\to e\tau$ and $H\to \mu\tau$ signals, obtained by fitting the data of the MC-template $\ell\tau_{had}$ channel for $e\tau_{had}$ final state in the VBF category.

BDT score distribution, after the simultaneous fit of the $H\to e\tau$ and $H\to \mu\tau$ signals, obtained by fitting the data of the MC-template $\ell\tau_{had}$ channel for $\mu\tau_{had}$ final state in the VBF category.

BDT score distribution, after the simultaneous fit of the $H\to e\tau$ and $H\to \mu\tau$ signals, obtained by fitting the data of the MC-template $\ell\tau_{had}$ channel for $\mu\tau_{had}$ final state in the VBF category.

NN score distribution, after an independent fit of the $H\to \mu\tau$ signal, obtained by fitting the data of the Symmetry $\ell\tau_{\ell'}$ channel, for the $\mu\tau_{e}$ final state in the non-VBF category.

NN score distribution, after an independent fit of the $H\to \mu\tau$ signal, obtained by fitting the data of the Symmetry $\ell\tau_{\ell'}$ channel, for the $\mu\tau_{e}$ final state in the non-VBF category.

NN score distribution, after an independent fit of the $H\to \mu\tau$ signal, obtained by fitting the data of the Symmetry $\ell\tau_{\ell'}$ channel, for the $\mu\tau_{e}$ final state in the VBF category.

NN score distribution, after an independent fit of the $H\to \mu\tau$ signal, obtained by fitting the data of the Symmetry $\ell\tau_{\ell'}$ channel, for the $\mu\tau_{e}$ final state in the VBF category.

Best-fit $\mathcal{B}(H\to \mu\tau)- \mathcal{B}(H\to e\tau)$ values, given in %, obtained from the standalone fits of the $\ell\tau_\ell'$ final state with either background estimation method. For the MC-template method, the correlated uncertainties are fixed to their best-fit values and the uncertainties are re-derived considering only the uncorrelated sources. For the Symmetry method, the full uncertainty is considered.

Best-fit $\mathcal{B}(H\to \mu\tau)- \mathcal{B}(H\to e\tau)$ values, given in %, obtained from the standalone fits of the $\ell\tau_\ell'$ final state with either background estimation method. For the MC-template method, the correlated uncertainties are fixed to their best-fit values and the uncertainties are re-derived considering only the uncorrelated sources. For the Symmetry method, the full uncertainty is considered.

Best-fit value of the branching ratios $\mathcal{B}(H\to e\tau)$ and $\mathcal{B}(H\to \mu\tau)$, given in %, and likelihood contours at 68% and 95% C.L. Obtained from the simultaneous fit of $H\to e\tau$ and $H\to \mu\tau$ signals based on the MC-template method, compared to the SM expectation.

Best-fit value of the branching ratios $\mathcal{B}(H\to e\tau)$ and $\mathcal{B}(H\to \mu\tau)$, given in %, and likelihood contours at 68% and 95% C.L. Obtained from the simultaneous fit of $H\to e\tau$ and $H\to \mu\tau$ signals based on the MC-template method, compared to the SM expectation.

Expected 2-D likelihood scan of the signal strength $\mu_{\tau\tau}$ versus the $CP$-mixing angle $\phi_{\tau}$.

Free-floating parameters in the measurement.

Leading tau Pt distributions obtained before the fit to data (Pre-Fit) showing the expected background and tuH signals after applying fake factors in the $t_{h}\tau_{lep}\tau_{had}$-2j region. Other MC includes single top, V+jets, and other small backgrounds. The tuH signals with nominal branching ratio of 0.1% are scaled using normalization factors of 2 to 50. Statistical and systematic uncertainties are included in the "Total background".

Leading tau Pt distributions obtained before the fit to data (Pre-Fit) showing the expected background and tuH signals after applying fake factors in the $t_{h}\tau_{lep}\tau_{had}$-3j region. Other MC includes single top, V+jets, and other small backgrounds. The tuH signals with nominal branching ratio of 0.1% are scaled using normalization factors of 2 to 50. Statistical and systematic uncertainties are included in the "Total background".

Leading tau Pt distributions obtained before the fit to data (Pre-Fit) showing the expected background and tuH signals after applying fake factors in the $t_{\ell}\tau_{had}\tau_{had}$SS region. Other MC includes single top, V+jets, and other small backgrounds. The tuH signals with nominal branching ratio of 0.1% are scaled using normalization factors of 2 to 50. Statistical and systematic uncertainties are included in the "Total background".

Leading tau Pt distributions obtained before the fit to data (Pre-Fit) showing the expected background and tuH signals after applying fake factors in the $t_{h}\tau_{had}\tau_{had}$-2j region. Rare includes single top, V+jets, and other small backgrounds. $\tau_{sub}$ real includes the contribution of fakes for which the sub-leading tau is real. The tuH signals with nominal branching ratio of 0.1% are scaled using normalization factors of 2 to 50. Statistical and systematic uncertainties are included in the "Total background".

Leading tau Pt distributions obtained before the fit to data (Pre-Fit) showing the expected background and tuH signals after applying fake factors in the $t_{h}\tau_{had}\tau_{had}$-3j region. Rare includes single top, V+jets, and other small backgrounds. $\tau_{sub}$ real includes the contribution of fakes for which the sub-leading tau is real. The tuH signals with nominal branching ratio of 0.1% are scaled using normalization factors of 2 to 50. Statistical and systematic uncertainties are included in the "Total background".

Leading tau Pt distributions obtained before the fit to data (Pre-Fit) showing the expected background and tuH signals after applying fake factors in the $t_{h}\tau_{had}\tau_{had}$-3jSS region. Rare includes single top, V+jets, and other small backgrounds. $\tau_{sub}$ real includes the contribution of fakes for which the sub-leading tau is real. The tuH signals with nominal branching ratio of 0.1% are scaled using normalization factors of 2 to 50. Statistical and systematic uncertainties are included in the "Total background".

Di-tau mass distributions obtained before the fit to data (Pre-Fit) showing the expected background and tuH signals after applying fake factors in the $t_{\ell}\tau_{had}\tau_{had}$ region. Other MC includes single top, V+jets, and other small backgrounds. The tuH signals with nominal branching ratio of 0.1% are scaled using normalization factors of 2 to 50. Statistical and systematic uncertainties are included in the "Total background".

Di-tau mass distributions obtained before the fit to data (Pre-Fit) showing the expected background and tuH signals after applying fake factors in the $t_{\ell}\tau_{had}$-1j region. Other MC includes single top, V+jets, and other small backgrounds. The tuH signals with nominal branching ratio of 0.1% are scaled using normalization factors of 2 to 50. Statistical and systematic uncertainties are included in the "Total background".

Di-tau mass distributions obtained before the fit to data (Pre-Fit) showing the expected background and tuH signals after applying fake factors in the $t_{\ell}\tau_{had}$-2j region. Other MC includes single top, V+jets, and other small backgrounds. The tuH signals with nominal branching ratio of 0.1% are scaled using normalization factors of 2 to 50. Statistical and systematic uncertainties are included in the "Total background".

Di-tau mass distributions obtained before the fit to data (Pre-Fit) showing the expected background and tuH signals after applying fake factors in the $t_{h}\tau_{lep}\tau_{had}$-2j region. Other MC includes single top, V+jets, and other small backgrounds. The tuH signals with nominal branching ratio of 0.1% are scaled using normalization factors of 2 to 50. Statistical and systematic uncertainties are included in the "Total background".

Di-tau mass distributions obtained before the fit to data (Pre-Fit) showing the expected background and tuH signals after applying fake factors in the $t_{h}\tau_{lep}\tau_{had}$-3j region. Other MC includes single top, V+jets, and other small backgrounds. The tuH signals with nominal branching ratio of 0.1% are scaled using normalization factors of 2 to 50. Statistical and systematic uncertainties are included in the "Total background".

Di-tau mass distributions obtained before the fit to data (Pre-Fit) showing the expected background and tuH signals after applying fake factors in the $t_{\ell}\tau_{had}\tau_{had}$SS region. Other MC includes single top, V+jets, and other small backgrounds. The tuH signals with nominal branching ratio of 0.1% are scaled using normalization factors of 2 to 50. Statistical and systematic uncertainties are included in the "Total background".

Di-tau mass distributions obtained before the fit to data (Pre-Fit) showing the expected background and tuH signals after applying fake factors in the $t_{h}\tau_{had}\tau_{had}$-2j region. Rare includes single top, V+jets, and other small backgrounds. $\tau_{sub}$ real includes the contribution of fakes for which the sub-leading tau is real. The tuH signals with nominal branching ratio of 0.1% are scaled using normalization factors of 2 to 50. Statistical and systematic uncertainties are included in the "Total background".

Di-tau mass distributions obtained before the fit to data (Pre-Fit) showing the expected background and tuH signals after applying fake factors in the $t_{h}\tau_{had}\tau_{had}$-3j region. Rare includes single top, V+jets, and other small backgrounds. $\tau_{sub}$ real includes the contribution of fakes for which the sub-leading tau is real. The tuH signals with nominal branching ratio of 0.1% are scaled using normalization factors of 2 to 50. Statistical and systematic uncertainties are included in the "Total background".

Di-tau mass distributions obtained before the fit to data (Pre-Fit) showing the expected background and tuH signals after applying fake factors in the $t_{h}\tau_{had}\tau_{had}$-3jSS region. Rare includes single top, V+jets, and other small backgrounds. $\tau_{sub}$ real includes the contribution of fakes for which the sub-leading tau is real. The tuH signals with nominal branching ratio of 0.1% are scaled using normalization factors of 2 to 50. Statistical and systematic uncertainties are included in the "Total background".

BDT output distributions obtained from a signal+background fit to the data for the tuH search in the $t_{\ell}\tau_{had}\tau_{had}$ region, Other MC includes single top, V+jets, and other small backgrounds. Statistical and systematic uncertainties are included in the "Total background". The signal shapes of tt(uH), tH, and their sum are also shown using a normalisation of 2 x $\mathcal{B}(\mathrm{t}\to\mathrm{uH})$ of 0.1%.

BDT output distributions obtained from a signal+background fit to the data for the tuH search in the $t_{\ell}\tau_{had}$-1j region, Other MC includes single top, V+jets, and other small backgrounds. Statistical and systematic uncertainties are included in the "Total background". The signal shapes of tt(uH), tH, and their sum are also shown using a normalisation of 2 x $\mathcal{B}(\mathrm{t}\to\mathrm{uH})$ of 0.1%.

BDT output distributions obtained from a signal+background fit to the data for the tuH search in the $t_{\ell}\tau_{had}$-2j region, Other MC includes single top, V+jets, and other small backgrounds. Statistical and systematic uncertainties are included in the "Total background". The signal shapes of tt(uH), tH, and their sum are also shown using a normalisation of 2 x $\mathcal{B}(\mathrm{t}\to\mathrm{uH})$ of 0.1%.

BDT output distributions obtained from a signal+background fit to the data for the tuH search in the $t_{h}\tau_{lep}\tau_{had}$-2j region, Other MC includes single top, V+jets, and other small backgrounds. Statistical and systematic uncertainties are included in the "Total background". The signal shapes of tt(uH), tH, and their sum are also shown using a normalisation of 2 x $\mathcal{B}(\mathrm{t}\to\mathrm{uH})$ of 0.1%.

BDT output distributions obtained from a signal+background fit to the data for the tuH search in the $t_{h}\tau_{lep}\tau_{had}$-3j region, Other MC includes single top, V+jets, and other small backgrounds. Statistical and systematic uncertainties are included in the "Total background". The signal shapes of tt(uH), tH, and their sum are also shown using a normalisation of 2 x $\mathcal{B}(\mathrm{t}\to\mathrm{uH})$ of 0.1%.

BDT output distributions obtained from a signal+background fit to the data for the tuH search in the $t_{\ell}\tau_{had}\tau_{had}$SS region, Other MC includes single top, V+jets, and other small backgrounds. Statistical and systematic uncertainties are included in the "Total background". The signal shapes of tt(uH), tH, and their sum are also shown using a normalisation of 2 x $\mathcal{B}(\mathrm{t}\to\mathrm{uH})$ of 0.1%.

BDT output distributions obtained from a signal+background fit to the data for the tuH search in the $t_{h}\tau_{had}\tau_{had}$-2j region, Rare includes single top, V+jets, and other small backgrounds. $\tau_{sub}$ real includes the contribution of fakes for which the sub-leading tau is real. Statistical and systematic uncertainties are included in the "Total background". The signal shapes of tt(uH), tH, and their sum are also shown using a normalisation of 2 x $\mathcal{B}(\mathrm{t}\to\mathrm{uH})$ of 0.1%.

BDT output distributions obtained from a signal+background fit to the data for the tuH search in the $t_{h}\tau_{had}\tau_{had}$-3j region, Rare includes single top, V+jets, and other small backgrounds. $\tau_{sub}$ real includes the contribution of fakes for which the sub-leading tau is real. Statistical and systematic uncertainties are included in the "Total background". The signal shapes of tt(uH), tH, and their sum are also shown using a normalisation of 2 x $\mathcal{B}(\mathrm{t}\to\mathrm{uH})$ of 0.1%.

BDT output distributions obtained from a signal+background fit to the data for the tuH search in the $t_{h}\tau_{had}\tau_{had}$-3jSS region, Rare includes single top, V+jets, and other small backgrounds. $\tau_{sub}$ real includes the contribution of fakes for which the sub-leading tau is real. Statistical and systematic uncertainties are included in the "Total background". The signal shapes of tt(uH), tH, and their sum are also shown using a normalisation of 2 x $\mathcal{B}(\mathrm{t}\to\mathrm{uH})$ of 0.1%.

BDT output distributions obtained from a signal+background fit to the data for the tcH search in the $t_{\ell}\tau_{had}\tau_{had}$ region, Other MC includes single top, V+jets, and other small backgrounds. Statistical and systematic uncertainties are included in the "Total background". The signal shapes of tt(cH), tH, and their sum are also shown using a normalisation of 2 x $\mathcal{B}(\mathrm{t}\to\mathrm{cH})$ of 0.1%.

BDT output distributions obtained from a signal+background fit to the data for the tcH search in the $t_{\ell}\tau_{had}$-1j region, Other MC includes single top, V+jets, and other small backgrounds. Statistical and systematic uncertainties are included in the "Total background". The signal shapes of tt(cH), tH, and their sum are also shown using a normalisation of 2 x $\mathcal{B}(\mathrm{t}\to\mathrm{cH})$ of 0.1%.

BDT output distributions obtained from a signal+background fit to the data for the tcH search in the $t_{\ell}\tau_{had}$-2j region, Other MC includes single top, V+jets, and other small backgrounds. Statistical and systematic uncertainties are included in the "Total background". The signal shapes of tt(cH), tH, and their sum are also shown using a normalisation of 2 x $\mathcal{B}(\mathrm{t}\to\mathrm{cH})$ of 0.1%.

BDT output distributions obtained from a signal+background fit to the data for the tcH search in the $t_{h}\tau_{lep}\tau_{had}$-2j region, Other MC includes single top, V+jets, and other small backgrounds. Statistical and systematic uncertainties are included in the "Total background". The signal shapes of tt(cH), tH, and their sum are also shown using a normalisation of 2 x $\mathcal{B}(\mathrm{t}\to\mathrm{cH})$ of 0.1%.

BDT output distributions obtained from a signal+background fit to the data for the tcH search in the $t_{h}\tau_{lep}\tau_{had}$-3j region,Other MC includes single top, V+jets, and other small backgrounds. Statistical and systematic uncertainties are included in the "Total background". The signal shapes of tt(cH), tH, and their sum are also shown using a normalisation of 2 x $\mathcal{B}(\mathrm{t}\to\mathrm{cH})$ of 0.1%.

BDT output distributions obtained from a signal+background fit to the data for the tcH search in the $t_{\ell}\tau_{had}\tau_{had}$SS region, Other MC includes single top, V+jets, and other small backgrounds. Statistical and systematic uncertainties are included in the "Total background". The signal shapes of tt(cH), tH, and their sum are also shown using a normalisation of 2 x $\mathcal{B}(\mathrm{t}\to\mathrm{cH})$ of 0.1%.

BDT output distributions obtained from a signal+background fit to the data for the tcH search in the $t_{h}\tau_{had}\tau_{had}$-2j region, Rare includes single top, V+jets, and other small backgrounds. $\tau_{sub}$ real includes the contribution of fakes for which the sub-leading tau is real. Statistical and systematic uncertainties are included in the "Total background". The signal shapes of tt(cH), tH, and their sum are also shown using a normalisation of 2 x $\mathcal{B}(\mathrm{t}\to\mathrm{cH})$ of 0.1%.

BDT output distributions obtained from a signal+background fit to the data for the tcH search in the $t_{h}\tau_{had}\tau_{had}$-3j region, Rare includes single top, V+jets, and other small backgrounds. $\tau_{sub}$ real includes the contribution of fakes for which the sub-leading tau is real. Statistical and systematic uncertainties are included in the "Total background". The signal shapes of tt(cH), tH, and their sum are also shown using a normalisation of 2 x $\mathcal{B}(\mathrm{t}\to\mathrm{cH})$ of 0.1%.

BDT output distributions obtained from a signal+background fit to the data for the tcH search in the $t_{h}\tau_{had}\tau_{had}$-3jSS region, Rare includes single top, V+jets, and other small backgrounds. $\tau_{sub}$ real includes the contribution of fakes for which the sub-leading tau is real. Statistical and systematic uncertainties are included in the "Total background". The signal shapes of tt(cH), tH, and their sum are also shown using a normalisation of 2 x $\mathcal{B}(\mathrm{t}\to\mathrm{cH})$ of 0.1%.

95% CL upper limits on $\mathcal{B}(\mathrm{t}\to\mathrm{cH})$ for the individual searches as well as their combination, assuming $\mathcal{B}(\mathrm{t}\to\mathrm{uH}) = 0$. The observed limits are compared with the expected (median) limits under the background-only hypothesis. The surrounding shaded bands correspond to the 68% and 95% CL intervals around the expected limits, denoted by $\pm 1\sigma$ and $\pm 2\sigma$, respectively.

95% CL upper limits on $\mathcal{B}(\mathrm{t}\to\mathrm{uH})$ for the individual searches as well as their combination, assuming $\mathcal{B}(\mathrm{t}\to\mathrm{cH}) = 0$. The observed limits are compared with the expected (median) limits under the background-only hypothesis. The surrounding shaded bands correspond to the 68% and 95% CL intervals around the expected limits, denoted by $\pm 1\sigma$ and $\pm 2\sigma$, respectively.

Observed upper limits at 95% CL on the branching fractions in the plane of $\mathcal{B}(\mathrm{t}\to\mathrm{uH})$ and $\mathcal{B}(\mathrm{t}\to\mathrm{cH})$.

Expected upper limits at 95% CL on the branching fractions in the plane of $\mathcal{B}(\mathrm{t}\to\mathrm{uH})$ and $\mathcal{B}(\mathrm{t}\to\mathrm{cH})$.

Expected $+2\sigma$ upper limits at 95% CL on the branching fractions in the plane of $\mathcal{B}(\mathrm{t}\to\mathrm{uH})$ and $\mathcal{B}(\mathrm{t}\to\mathrm{cH})$.

Expected $+1\sigma$ upper limits at 95% CL on the branching fractions in the plane of $\mathcal{B}(\mathrm{t}\to\mathrm{uH})$ and $\mathcal{B}(\mathrm{t}\to\mathrm{cH})$.

Expected $-1\sigma$ upper limits at 95% CL on the branching fractions in the plane of $\mathcal{B}(\mathrm{t}\to\mathrm{uH})$ and $\mathcal{B}(\mathrm{t}\to\mathrm{cH})$.

Expected $-2\sigma$ upper limits at 95% CL on the branching fractions in the plane of $\mathcal{B}(\mathrm{t}\to\mathrm{uH})$ and $\mathcal{B}(\mathrm{t}\to\mathrm{cH})$.

Observed upper limits at 95% CL on the anomalous couplings in the plane of $C_{\mathrm{c\phi}}$ and $C_{\mathrm{u\phi}}$.

Expected upper limits at 95% CL on the anomalous couplings in the plane of $C_{\mathrm{c\phi}}$ and $C_{\mathrm{u\phi}}$.

Expected $+2\sigma$ upper limits at 95% CL on the anomalous couplings in the plane of $C_{\mathrm{c\phi}}$ and $C_{\mathrm{u\phi}}$.

Expected $+1\sigma$ upper limits at 95% CL on the anomalous couplings in the plane of $C_{\mathrm{c\phi}}$ and $C_{\mathrm{u\phi}}$.

Expected $-1\sigma$ upper limits at 95% CL on the anomalous couplings in the plane of $C_{\mathrm{c\phi}}$ and $C_{\mathrm{u\phi}}$.

Expected $-2\sigma$ upper limits at 95% CL on the anomalous couplings in the plane of $C_{\mathrm{c\phi}}$ and $C_{\mathrm{u\phi}}$.

Predicted and observed yields in each of the analysis regions considered in leptonic channel.

Predicted and observed yields in each of the analysis regions considered in hadronic channel.

Absolute uncertainties on $\mathcal{B}(\mathrm{t}\to\mathrm{qH})$ obtained from the combined fit to the data.

Summary of 95% CL upper limits on $\mathcal{B}(\mathrm{t}\to\mathrm{qH})$, significance and best-fit branching ratio in the signal regions. The values in the tables are in the form of observed(expected).

Expected, with the corresponding $\pm 1\sigma$ intervals, and observed 95% CL branching fraction upper limits for the Higgs and $Z$ boson decays into a quarkonium state and a photon. Standard Model production of the Higgs boson is assumed. The corresponding upper limits on the production cross section times branching fraction $\sigma\times\mathcal{B}$ are also shown.

Observed yields in the boost_0_ge2J signal region category of the hh channel.

Observed yields in the boost_1_1J signal region category of the hh channel.

Observed yields in the boost_1_1J signal region category of the hh channel.

Observed yields in the boost_1_ge2J signal region category of the hh channel.

Observed yields in the boost_1_ge2J signal region category of the hh channel.

Observed yields in the boost_2 signal region category of the hh channel.

Observed yields in the boost_2 signal region category of the hh channel.

Observed yields in the boost_3 signal region category of the hh channel.

Observed yields in the boost_3 signal region category of the hh channel.

Observed yields in the ttH_0 signal region category of the hh channel.

Observed yields in the ttH_0 signal region category of the hh channel.

Observed yields in the ttH_1 signal region category of the hh channel.

Observed yields in the ttH_1 signal region category of the hh channel.

Observed yields in the VH_0 signal region category of the hh channel.

Observed yields in the VH_0 signal region category of the hh channel.

Observed yields in the VH_1 signal region category of the hh channel.

Observed yields in the VH_1 signal region category of the hh channel.

Observed yields in the VBF_0 signal region category of the hh channel.

Observed yields in the VBF_0 signal region category of the hh channel.

Observed yields in the VBF_1 signal region category of the hh channel.

Observed yields in the VBF_1 signal region category of the hh channel.

Observed yields in the boost_0_1J signal region category of the lh channel.

Observed yields in the boost_0_1J signal region category of the lh channel.

Observed yields in the boost_0_ge2J signal region category of the lh channel.

Observed yields in the boost_0_ge2J signal region category of the lh channel.

Observed yields in the boost_1_1J signal region category of the lh channel.

Observed yields in the boost_1_1J signal region category of the lh channel.

Observed yields in the boost_1_ge2J signal region category of the lh channel.

Observed yields in the boost_1_ge2J signal region category of the lh channel.

Observed yields in the boost_2 signal region category of the lh channel.

Observed yields in the boost_2 signal region category of the lh channel.

Observed yields in the boost_3 signal region category of the lh channel.

Observed yields in the boost_3 signal region category of the lh channel.

Observed yields in the VH_0 signal region category of the lh channel.

Observed yields in the VH_0 signal region category of the lh channel.

Observed yields in the VH_1 signal region category of the lh channel.

Observed yields in the VH_1 signal region category of the lh channel.

Observed yields in the VBF_0 signal region category of the lh channel.

Observed yields in the VBF_0 signal region category of the lh channel.

Observed yields in the VBF_1 signal region category of the lh channel.

Observed yields in the VBF_1 signal region category of the lh channel.

Observed yields in the boost_0_1J signal region category of the ll channel.

Observed yields in the boost_0_1J signal region category of the ll channel.

Observed yields in the boost_0_ge2J signal region category of the ll channel.

Observed yields in the boost_0_ge2J signal region category of the ll channel.

Observed yields in the boost_1_1J signal region category of the ll channel.

Observed yields in the boost_1_1J signal region category of the ll channel.

Observed yields in the boost_1_ge2J signal region category of the ll channel.

Observed yields in the boost_1_ge2J signal region category of the ll channel.

Observed yields in the boost_2 signal region category of the ll channel.

Observed yields in the boost_2 signal region category of the ll channel.

Observed yields in the boost_3 signal region category of the ll channel.

Observed yields in the boost_3 signal region category of the ll channel.

Observed yields in the VH_0 signal region category of the ll channel.

Observed yields in the VH_0 signal region category of the ll channel.

Observed yields in the VH_1 signal region category of the ll channel.

Observed yields in the VH_1 signal region category of the ll channel.

Observed yields in the VBF_0 signal region category of the ll channel.

Observed yields in the VBF_0 signal region category of the ll channel.

Observed yields in the VBF_1 signal region category of the ll channel.

Observed yields in the VBF_1 signal region category of the ll channel.

Best-fit values and uncertainties for the pp $\rightarrow H\rightarrow\tau\tau$ cross-section measurement and the measurement in the four dominant production modes. All measurements include the branching ratio of $H\rightarrow\tau\tau$ and are performed with true Higgs boson rapidity $|y_{H}|<2.5$. The SM predictions for each region, computed using the inclusive cross-section calculations and the simulated event samples are also shown. The contributions to the total uncertainty in the measurements from statistical (Stat. unc.) or systematic uncertainties (Syst. unc.) are given separately.

Best-fit values and uncertainties for the pp $\rightarrow H\rightarrow\tau\tau$ cross-section measurement and the measurement in the four dominant production modes. All measurements include the branching ratio of $H\rightarrow\tau\tau$ and are performed with true Higgs boson rapidity $|y_{H}|<2.5$. The SM predictions for each region, computed using the inclusive cross-section calculations and the simulated event samples are also shown. The contributions to the total uncertainty in the measurements from statistical (Stat. unc.) or systematic uncertainties (Syst. unc.) are given separately.

The measured correlations between each parameter of interest in the measurement of the cross-sections per production mode

The measured correlations between each parameter of interest in the measurement of the cross-sections per production mode

Best-fit values and uncertainties for the $H\rightarrow\tau\tau$ cross-sections, in the reduced stage 1.2 STXS scheme. The EW production mode includes vector-boson fusion and qq$\rightarrow$V($\rightarrow$qq)H processes. All measurements include the branching ratio of $H\rightarrow\tau\tau$ and are performed with true Higgs boson rapidity $|y_{H}|<2.5$. The SM predictions for each region, computed using the inclusive cross-section calculations and the simulated event samples are also shown. The contributions to the total uncertainty in the measurements from statistical (Stat. unc.) or systematic uncertainties (Syst. unc.) are given separately. The asterisk symbol (*) indicates that the criteria for m$_{jj}$ only apply to events with at least two reconstructed jets.

Best-fit values and uncertainties for the $H\rightarrow\tau\tau$ cross-sections, in the reduced stage 1.2 STXS scheme. The EW production mode includes vector-boson fusion and qq$\rightarrow$V($\rightarrow$qq)H processes. All measurements include the branching ratio of $H\rightarrow\tau\tau$ and are performed with true Higgs boson rapidity $|y_{H}|<2.5$. The SM predictions for each region, computed using the inclusive cross-section calculations and the simulated event samples are also shown. The contributions to the total uncertainty in the measurements from statistical (Stat. unc.) or systematic uncertainties (Syst. unc.) are given separately. The asterisk symbol (*) indicates that the criteria for m$_{jj}$ only apply to events with at least two reconstructed jets.

The measured correlations between each pair of parameters of interest in the STXS measurement. The asterisk symbol (*) indicates that the criteria for mjj only apply to events with at least two reconstructed jets.

The measured correlations between each pair of parameters of interest in the STXS measurement. The asterisk symbol (*) indicates that the criteria for mjj only apply to events with at least two reconstructed jets.