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

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

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

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

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

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

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

Product of acceptance and efficiency for pp->tbH(->Wh) as function of the charged Higgs boson mass for the resolved qqbb signal region.

Product of acceptance and efficiency for pp->tbH(->Wh) as function of the charged Higgs boson mass for the resolved lvbb signal region.

Product of acceptance and efficiency for pp->tbH(->Wh) as function of the charged Higgs boson mass for the merged lvbb signal region

Product of acceptance and efficiency for pp->tbH(->Wh) as function of the charged Higgs boson mass for the merged qqbb signal region

Distributions of the mWh observable in the low-purity signal regions of the resolved qqbb 5jex3bex event categories. The term ‘Others’ summarises events from tHjb, tWh, tttt, and SM Vh production. The distributions are presented after a background-only maximum-likelihood fit to data. The individual background uncertainty does not take into account the possible correlations between the nuisance parameter. The expected signal contribution assuming $m_{H^{\pm}}$ = 700 GeV, normalised to the expected limit of the cross-section times branching ratio ($\sigma_{sig}$ × B) of 0.064 pb, is shown as a dashed histogram.

Distributions of the mWh observable in the low-purity signal regions of the resolved qqbb 5jex4bin event categories. The term ‘Others’ summarises events from tHjb, tWh, tttt, and SM Vh production. The distributions are presented after a background-only maximum-likelihood fit to data. The individual background uncertainty does not take into account the possible correlations between the nuisance parameter. The expected signal contribution assuming $m_{H^{\pm}}$ = 700 GeV, normalised to the expected limit of the cross-section times branching ratio ($\sigma_{sig}$ × B) of 0.064 pb, is shown as a dashed histogram.

Distributions of the mWh observable in the low-purity signal regions of the resolved qqbb 6jin3bex event categories. The term ‘Others’ summarises events from tHjb, tWh, tttt, and SM Vh production. The distributions are presented after a background-only maximum-likelihood fit to data. The individual background uncertainty does not take into account the possible correlations between the nuisance parameter. The expected signal contribution assuming $m_{H^{\pm}}$ = 700 GeV, normalised to the expected limit of the cross-section times branching ratio ($\sigma_{sig}$ × B) of 0.064 pb, is shown as a dashed histogram.

Distributions of the mWh observable in the low-purity signal regions of the resolved qqbb 6jin4bin event categories. The term ‘Others’ summarises events from tHjb, tWh, tttt, and SM Vh production. The distributions are presented after a background-only maximum-likelihood fit to data. The individual background uncertainty does not take into account the possible correlations between the nuisance parameter. The expected signal contribution assuming $m_{H^{\pm}}$ = 700 GeV, normalised to the expected limit of the cross-section times branching ratio ($\sigma_{sig}$ × B) of 0.064 pb, is shown as a dashed histogram.

Distributions of the mWh observable in the high-purity signal regions of the resolved qqbb 5jex3bex event categories. The term ‘Others’ summarises events from tHjb, tWh, tttt, and SM Vh production. The distributions are presented after a background-only maximum-likelihood fit to data. The individual background uncertainty does not take into account the possible correlations between the nuisance parameter. The expected signal contribution assuming $m_{H^{\pm}}$ = 700 GeV, normalised to the expected limit of the cross-section times branching ratio ($\sigma_{sig}$ × B) of 0.064 pb, is shown as a dashed histogram.

Distributions of the mWh observable in the high-purity signal regions of the resolved qqbb 5jex4bin event categories. The term ‘Others’ summarises events from tHjb, tWh, tttt, and SM Vh production. The distributions are presented after a background-only maximum-likelihood fit to data. The individual background uncertainty does not take into account the possible correlations between the nuisance parameter. The expected signal contribution assuming $m_{H^{\pm}}$ = 700 GeV, normalised to the expected limit of the cross-section times branching ratio ($\sigma_{sig}$ × B) of 0.064 pb, is shown as a dashed histogram.

Distributions of the mWh observable in the high-purity signal regions of the resolved qqbb 6jin3bex event categories. The term ‘Others’ summarises events from tHjb, tWh, tttt, and SM Vh production. The distributions are presented after a background-only maximum-likelihood fit to data. The individual background uncertainty does not take into account the possible correlations between the nuisance parameter. The expected signal contribution assuming $m_{H^{\pm}}$ = 700 GeV, normalised to the expected limit of the cross-section times branching ratio ($\sigma_{sig}$ × B) of 0.064 pb, is shown as a dashed histogram.

Distributions of the mWh observable in the high-purity signal regions of the resolved qqbb 6jin4bin event categories. The term ‘Others’ summarises events from tHjb, tWh, tttt, and SM Vh production. The distributions are presented after a background-only maximum-likelihood fit to data. The individual background uncertainty does not take into account the possible correlations between the nuisance parameter. The expected signal contribution assuming $m_{H^{\pm}}$ = 700 GeV, normalised to the expected limit of the cross-section times branching ratio ($\sigma_{sig}$ × B) of 0.064 pb, is shown as a dashed histogram.

Distributions of the mWh observable in the signal regions of the resolved lvbb 5jex3bex event categories. The term ‘Others’ summarises events from tHjb, tWh, tttt, and SM Vh production. The distributions are presented after a background-only maximum-likelihood fit to data. The individual background uncertainty does not take into account the possible correlations between the nuisance parameter. The expected signal contribution assuming $m_{H^{\pm}}$ = 700 GeV, normalised to the expected limit of the cross-section times branching ratio ($\sigma_{sig}$ × B) of 0.064 pb, is shown as a dashed histogram.

Distributions of the mWh observable in the signal regions of the resolved lvbb 5jex4bin event categories. The term ‘Others’ summarises events from tHjb, tWh, tttt, and SM Vh production. The distributions are presented after a background-only maximum-likelihood fit to data. The individual background uncertainty does not take into account the possible correlations between the nuisance parameter. The expected signal contribution assuming $m_{H^{\pm}}$ = 700 GeV, normalised to the expected limit of the cross-section times branching ratio ($\sigma_{sig}$ × B) of 0.064 pb, is shown as a dashed histogram.

Distributions of the mWh observable in the signal regions of the resolved lvbb 6jin3bex event categories. The term ‘Others’ summarises events from tHjb, tWh, tttt, and SM Vh production. The distributions are presented after a background-only maximum-likelihood fit to data. The individual background uncertainty does not take into account the possible correlations between the nuisance parameter. The expected signal contribution assuming $m_{H^{\pm}}$ = 700 GeV, normalised to the expected limit of the cross-section times branching ratio ($\sigma_{sig}$ × B) of 0.064 pb, is shown as a dashed histogram.

Distributions of the mWh observable in the signal regions of the resolved lvbb 6jin4bin event categories. The term ‘Others’ summarises events from tHjb, tWh, tttt, and SM Vh production. The distributions are presented after a background-only maximum-likelihood fit to data. The individual background uncertainty does not take into account the possible correlations between the nuisance parameter. The expected signal contribution assuming $m_{H^{\pm}}$ = 700 GeV, normalised to the expected limit of the cross-section times branching ratio ($\sigma_{sig}$ × B) of 0.064 pb, is shown as a dashed histogram.

The mWh distributions in the Low-NN-score SRs of the merged qqbb 0b event categories. The term ‘Others’ summarises events from tHjb, tWh, tttt, and SM Vh production. The background prediction is shown after a background-only maximum-likelihood fit to data. The individual background uncertainty does not take into account the possible correlations between the nuisance parameter. The expected signal contributions assuming $m_{H^{\pm}}$ = 900 GeV and $m_{H^{\pm}}$ = 2000 GeV, normalised the expected limit of the cross-section times branching ratio ($\sigma_{sig}$ × B) values of 0.036 pb and 2.7 fb respectively, are shown as dashed histograms.

The mWh distributions in the Low-NN-score SRs of the merged qqbb 1b event categories. The term ‘Others’ summarises events from tHjb, tWh, tttt, and SM Vh production. The background prediction is shown after a background-only maximum-likelihood fit to data. The individual background uncertainty does not take into account the possible correlations between the nuisance parameter. The expected signal contributions assuming $m_{H^{\pm}}$ = 900 GeV and $m_{H^{\pm}}$ = 2000 GeV, normalised the expected limit of the cross-section times branching ratio ($\sigma_{sig}$ × B) values of 0.036 pb and 2.7 fb respectively, are shown as dashed histograms.

The mWh distributions in the High-NN-score SRs of the merged qqbb 0b event categories. The term ‘Others’ summarises events from tHjb, tWh, tttt, and SM Vh production. The background prediction is shown after a background-only maximum-likelihood fit to data. The individual background uncertainty does not take into account the possible correlations between the nuisance parameter. The expected signal contributions assuming $m_{H^{\pm}}$ = 900 GeV and $m_{H^{\pm}}$ = 2000 GeV, normalised the expected limit of the cross-section times branching ratio ($\sigma_{sig}$ × B) values of 0.036 pb and 2.7 fb respectively, are shown as dashed histograms.

The mWh distributions in the High-NN-score SRs of the merged qqbb 1b event categories. The term ‘Others’ summarises events from tHjb, tWh, tttt, and SM Vh production. The background prediction is shown after a background-only maximum-likelihood fit to data. The individual background uncertainty does not take into account the possible correlations between the nuisance parameter. The expected signal contributions assuming $m_{H^{\pm}}$ = 900 GeV and $m_{H^{\pm}}$ = 2000 GeV, normalised the expected limit of the cross-section times branching ratio ($\sigma_{sig}$ × B) values of 0.036 pb and 2.7 fb respectively, are shown as dashed histograms.

The mWh distributions in the Low-NN-score SRs of the merged lvbb 0b event categories. The term ‘Others’ summarises events from tHjb, tWh, tttt, and SM Vh production. The background prediction is shown after a background-only maximum-likelihood fit to data. The individual background uncertainty does not take into account the possible correlations between the nuisance parameter. The expected signal contributions assuming $m_{H^{\pm}}$ = 900 GeV and $m_{H^{\pm}}$ = 2000 GeV, normalised the expected limit of the cross-section times branching ratio ($\sigma_{sig}$ × B) values of 0.036 pb and 2.7 fb respectively, are shown as dashed histograms.

The mWh distributions in the Low-NN-score SRs of the merged lvbb 1b event categories. The term ‘Others’ summarises events from tHjb, tWh, tttt, and SM Vh production. The background prediction is shown after a background-only maximum-likelihood fit to data. The individual background uncertainty does not take into account the possible correlations between the nuisance parameter. The expected signal contributions assuming $m_{H^{\pm}}$ = 900 GeV and $m_{H^{\pm}}$ = 2000 GeV, normalised the expected limit of the cross-section times branching ratio ($\sigma_{sig}$ × B) values of 0.036 pb and 2.7 fb respectively, are shown as dashed histograms.

The mWh distributions in the Medium-NN-score SRs of the merged lvbb 0b event categories. The term ‘Others’ summarises events from tHjb, tWh, tttt, and SM Vh production. The background prediction is shown after a background-only maximum-likelihood fit to data. The individual background uncertainty does not take into account the possible correlations between the nuisance parameter. The expected signal contributions assuming $m_{H^{\pm}}$ = 900 GeV and $m_{H^{\pm}}$ = 2000 GeV, normalised the expected limit of the cross-section times branching ratio ($\sigma_{sig}$ × B) values of 0.036 pb and 2.7 fb respectively, are shown as dashed histograms.

The mWh distributions in the Medium-NN-score SRs of the merged lvbb 1b event categories. The term ‘Others’ summarises events from tHjb, tWh, tttt, and SM Vh production. The background prediction is shown after a background-only maximum-likelihood fit to data. The individual background uncertainty does not take into account the possible correlations between the nuisance parameter. The expected signal contributions assuming $m_{H^{\pm}}$ = 900 GeV and $m_{H^{\pm}}$ = 2000 GeV, normalised the expected limit of the cross-section times branching ratio ($\sigma_{sig}$ × B) values of 0.036 pb and 2.7 fb respectively, are shown as dashed histograms.

The mWh distributions in the High-NN-score SRs of the merged lvbb 0b event categories. The term ‘Others’ summarises events from tHjb, tWh, tttt, and SM Vh production. The background prediction is shown after a background-only maximum-likelihood fit to data. The individual background uncertainty does not take into account the possible correlations between the nuisance parameter. The expected signal contributions assuming $m_{H^{\pm}}$ = 900 GeV and $m_{H^{\pm}}$ = 2000 GeV, normalised the expected limit of the cross-section times branching ratio ($\sigma_{sig}$ × B) values of 0.036 pb and 2.7 fb respectively, are shown as dashed histograms.

The mWh distributions in the High-NN-score SRs of the merged lvbb 1b event categories. The term ‘Others’ summarises events from tHjb, tWh, tttt, and SM Vh production. The background prediction is shown after a background-only maximum-likelihood fit to data. The individual background uncertainty does not take into account the possible correlations between the nuisance parameter. The expected signal contributions assuming $m_{H^{\pm}}$ = 900 GeV and $m_{H^{\pm}}$ = 2000 GeV, normalised the expected limit of the cross-section times branching ratio ($\sigma_{sig}$ × B) values of 0.036 pb and 2.7 fb respectively, are shown as dashed histograms.

Cutflow table for a few selected signal samples after applying the selection requirements of the resolved analysis. At the initial cut stage, all events from the tbH ± → W ± h(→ b b) production process (including those from the fully hadronic decay channel) are considered.

Cutflow table for a few selected signal samples after applying the selection requirements of the merged analysis. At the initial cut stage, all events from the tbH ± → W ± h(→ bb) production process (including those from the fully hadronic decay channel) are considered.

The impact of the eight most important systematic uncertainties on \(R_t\), in %.

The post-fit yields in the two SRs. All uncertainties applied in the analysis are included.

Limits on EFT parameters and Vtb extracted from the analysis. The upper and lower limits correspond to 95% CL.

Results of the main analysis.

Pre-fit distribution of the NN discriminant \(D_{\text{nn}}\) in SR plus.

Post-fit distribution of the NN discriminant \(D_{\text{nn}}\) in SR plus.

Post-fit distribution of the \(\Delta\phi(\vec{p}_{T}^{miss},\mu)\) variable in the muon-plus CR.

Pre-fit distribution of the NN discriminant \(D_{\text{nn}}\) in SR minus.

Post-fit distribution of the NN discriminant \(D_{\text{nn}}\) in SR minus.

Post-fit distribution of the \(\Delta\phi(\vec{p}_{T}^{miss},\mu)\) variable in the muon-minus CR.

Pre-fit distribution of the invariant mass of the untagged jet and the b-tagged jet, \(m(jb)\), in SR plus.

Pre-fit distribution of the absolute value of the pseudorapidity of the untagged jet, \(|\eta(j)|\), in SR plus.

Pre-fit distribution of the absolute value of the difference in \(p_{\text{T}}\) between the reconstructed W boson and the jet pair, \(|\Delta p_{\text{T}}(W, jb)|\), in SR plus.

Pre-fit distribution of the difference in azimuth angle between the reconstructed W boson and the jet pair, \(|\Delta\phi(W, jb)|\), in SR plus.

Pre-fit distribution of the invariant mass of the reconstructed top quark, \(m(t)\), in SR plus.

Pre-fit distribution of the absolute value of the difference in pseudorapidity between the charged lepton and the untagged jet, \(|\Delta\eta(\ell ,j)|\), SR plus.

Pre-fit distribution of the angular distance of the charged lepton and the untagged jet, \(\Delta R(\ell ,j)\), in SR plus.

Pre-fit distribution of the absolute value of the difference in pseudorapidity between the b-tagged jet and the charged lepton, \(|\Delta\eta(b,\ell )|\), in SR plus.

Observed exclusion limits at the 95% CL in the $\tilde{t}_1 - \tilde\chi^0_1$ mass plane, assuming BR$(\tilde{t}_1 \rightarrow t + \tilde\chi^0_1) = 0.5$ and a $-1 \sigma$ deviation of the NNLO+NNLL theoretical cross-section of a $\tilde{t}_1$ pair-production.

Expected exclusion limits at the 95% CL in the $\tilde{t}_1 - \tilde\chi^0_1$ mass plane, assuming BR$(\tilde{t}_1 \rightarrow t + \tilde\chi^0_1) = 0.5$.

The $1\sigma$ variation of expected 95% CL exclusion contour obtained by varying MC statistical uncertainties, detector-related systematic uncertainties, and theoretical uncertainties in the $\tilde{t}_1 - \tilde\chi^0_1$ mass plane, assuming BR$(\tilde{t}_1 \rightarrow t + \tilde\chi^0_1) = 0.5$.

The $2\sigma$ variation of expected 95%CL exclusion contour obtained by varying MC statistical uncertainties, detector-related systematic uncertainties, and theoretical uncertainties in the $\tilde{t}_1 - \tilde\chi^0_1$ mass plane, assuming BR$(\tilde{t}_1 \rightarrow t + \tilde\chi^0_1) = 0.5$.

Observed exclusion limits at the 95% CL in the $m_{\tilde{t}_1} - $BR$(\tilde{t}_1 \rightarrow t + \tilde\chi^0_1)$ plane, assuming $m_{\tilde\chi^0_1} = 1 \mathrm{GeV}$.

Observed exclusion limits at the 95% CL in the $\tilde{t}_1 - \tilde\chi^0_1$ mass plane, assuming $m_{\tilde\chi^0_1} = 1 \mathrm{ GeV}$ and a $+1 \sigma$ deviation of the NNLO+NNLL theoretical cross-section of a $\tilde{t}_1$ pair-production.

Observed exclusion limits at the 95% CL in the $\tilde{t}_1 - \tilde\chi^0_1$ mass plane, assuming $m_{\tilde\chi^0_1} = 1 \mathrm{ GeV}$ and a $-1 \sigma$ deviation of the NNLO+NNLL theoretical cross-section of a $\tilde{t}_1$ pair-production.

Expected exclusion limits at the 95% CL in the $m_{\tilde{t}_1} - $BR$(\tilde{t}_1 \rightarrow t + \tilde\chi^0_1)$ plane, assuming $m_{\tilde\chi^0_1} = 1 \mathrm{GeV}$.

The $1\sigma$ variation of expected 95% CL exclusion contour obtained by varying MC statistical uncertainties, detector-related systematic uncertainties, and theoretical uncertainties in the $m_{\tilde{t}_1} - $BR$(\tilde{t}_1 \rightarrow t + \tilde\chi^0_1)$ plane, assuming $m_{\tilde\chi^0_1} = 1 \mathrm{GeV}$.

The $2\sigma$ variation of expected 95%CL exclusion contour obtained by varying MC statistical uncertainties, detector-related systematic uncertainties, and theoretical uncertainties in the $m_{\tilde{t}_1} - $BR$(\tilde{t}_1 \rightarrow t + \tilde\chi^0_1)$ plane, assuming $m_{\tilde\chi^0_1} = 1 \mathrm{GeV}$.

Observed upper limits on the top-spartner pair production cross-section in fb at the 95% CL in the $\tilde{t}_1 - \tilde\chi^0_1$ mass plane, assuming BR$(\tilde{t}_1 \rightarrow t + \tilde\chi^0_1) = 0.5$

Expected upper limits on the top-spartner pair production cross-section in fb at the 95% CL in the $\tilde{t}_1 - \tilde\chi^0_1$ mass plane, assuming BR$(\tilde{t}_1 \rightarrow t + \tilde\chi^0_1) = 0.5$

Observed upper limits on the top-spartner pair production cross-section in fb at the 95% CL in the $m_{\tilde{t}_1} - $BR$(\tilde{t}_1 \rightarrow t + \tilde\chi^0_1)$ plane, assuming $m_{\tilde\chi^0_1} = 1 \mathrm{GeV}$.

Expected upper limits on the top-spartner pair production cross-section in fb at the 95% CL in the $m_{\tilde{t}_1} - $BR$(\tilde{t}_1 \rightarrow t + \tilde\chi^0_1)$ plane, assuming $m_{\tilde\chi^0_1} = 1 \mathrm{GeV}$.

Post-fit distribution of $N_{\mathrm{tops}}^{\mathrm{DNN}} (\mathrm{R=1.0})$ in the SRA signal region presented without the associated SRA applied to the variable. For each bin yields for the data and total SM prediction are provided. The SM prediction is provided with the total uncertainty, including the MC statistical uncertainty, detector-related systematic uncertainties and theoretical uncertainties. The rightmost bin includes overflow events.

Post-fit distribution of $m_{\mathrm{T2}}(j^{b}_{R=1.0}, c)$ in the SRA signal region presented without the associated SRA applied to the variable. For each bin yields for the data and total SM prediction are provided. The SM prediction is provided with the total uncertainty, including the MC statistical uncertainty, detector-related systematic uncertainties and theoretical uncertainties. The rightmost bin includes overflow events.

Post-fit distribution of $p_{\mathrm{T}}(c_{1})$ in the SRB signal region presented without the associated SRB applied to the variable. For each bin yields for the data and total SM prediction are provided. The SM prediction is provided with the total uncertainty, including the MC statistical uncertainty, detector-related systematic uncertainties and theoretical uncertainties. The rightmost bin includes overflow events.

Post-fit distribution of $m_{\mathrm{T}}(j, \mathrm{E}_{\mathrm{T}}^{\mathrm{miss}})_{\mathrm{close}}$ in the SRB signal region presented without the associated SRB applied to the variable. For each bin yields for the data and total SM prediction are provided. The SM prediction is provided with the total uncertainty, including the MC statistical uncertainty, detector-related systematic uncertainties and theoretical uncertainties. The rightmost bin includes overflow events.

Post-fit distribution of $\mathrm{E}_{\mathrm{T}}^{\mathrm{miss}} \textrm{Significance}$ in the SRC signal region presented without the associated SRC applied to the variable. For each bin yields for the data and total SM prediction are provided. The SM prediction is provided with the total uncertainty, including the MC statistical uncertainty, detector-related systematic uncertainties and theoretical uncertainties. The rightmost bin includes overflow events.

Post-fit distribution of $m_{\mathrm{T}}(j, \mathrm{E}_{\mathrm{T}}^{\mathrm{miss}})_{\mathrm{close}}$ in the SRC signal region presented without the associated SRC applied to the variable. For each bin yields for the data and total SM prediction are provided. The SM prediction is provided with the total uncertainty, including the MC statistical uncertainty, detector-related systematic uncertainties and theoretical uncertainties. The rightmost bin includes overflow events.

Post-fit distribution of NN signal score in the SRD signal region presented without the associated SRD applied to the variable. For each bin yields for the data and total SM prediction are provided. The SM prediction is provided with the total uncertainty, including the MC statistical uncertainty, detector-related systematic uncertainties and theoretical uncertainties. The rightmost bin includes overflow events.

Post-fit distribution of $m_{\mathrm{eff}}$ in the SRD signal region. For each bin yields for the data and total SM prediction are provided. The SM prediction is provided with the total uncertainty, including the MC statistical uncertainty, detector-related systematic uncertainties and theoretical uncertainties. The rightmost bin includes overflow events.

Pull plots showing the SRA, SRB and SRC post-fit data and SM agreement using the background-only fit configuration

Pull plots showing the SRD post-fit data and SM agreement using the background-only fit configuration

Cutflow for the reference point $(m_{\tilde{t}_1}, m_{\tilde\chi^0_1})=$ (900,1) , (700,300), (550,375) in the Signal region A. Results are shown including all correction factors applied to simulation, and is normalised to 139 fb$^{-1}$.

Cutflow for the reference point $(m_{\tilde{t}_1}, m_{\tilde\chi^0_1})=$ (900,1) , (700,300), (550,375) in the Signal region B. Results are shown including all correction factors applied to simulation, and is normalised to 139 fb$^{-1}$.

Cutflow for the reference point $(m_{\tilde{t}_1}, m_{\tilde\chi^0_1})=$ (900,1) , (700,300), (550,375) in the Signal region C. Results are shown including all correction factors applied to simulation, and is normalised to 139 fb$^{-1}$.

Cutflow for the reference point $(m_{\tilde{t}_1}, m_{\tilde\chi^0_1})=$ (900,1) , (700,300), (550,375) in the Signal region D 750. Results are shown including all correction factors applied to simulation, and is normalised to 139 fb$^{-1}$.

Cutflow for the reference point $(m_{\tilde{t}_1}, m_{\tilde\chi^0_1})=$ (900,1) , (700,300), (550,375) in the Signal region D 1000. Results are shown including all correction factors applied to simulation, and is normalised to 139 fb$^{-1}$.

Cutflow for the reference point $(m_{\tilde{t}_1}, m_{\tilde\chi^0_1})=$ (900,1) , (700,300), (550,375) in the Signal region D 1250. Results are shown including all correction factors applied to simulation, and is normalised to 139 fb$^{-1}$.

Cutflow for the reference point $(m_{\tilde{t}_1}, m_{\tilde\chi^0_1})=$ (900,1) , (700,300), (550,375) in the Signal region D 1500. Results are shown including all correction factors applied to simulation, and is normalised to 139 fb$^{-1}$.

Cutflow for the reference point $(m_{\tilde{t}_1}, m_{\tilde\chi^0_1})=$ (900,1) , (700,300), (550,375) in the Signal region D 1750. Results are shown including all correction factors applied to simulation, and is normalised to 139 fb$^{-1}$.

Cutflow for the reference point $(m_{\tilde{t}_1}, m_{\tilde\chi^0_1})=$ (900,1) , (700,300), (550,375) in the Signal region D 2000. Results are shown including all correction factors applied to simulation, and is normalised to 139 fb$^{-1}$.

Signal acceptance in the SRA$^{[450,575]}$ region in the $\tilde{t}_1 - \tilde\chi^0_1$ mass plane, assuming BR$(\tilde{t}_1 \rightarrow t + \tilde\chi^0_1) = 0.5$. Please mind that the acceptance given in the table is multiplied by factor of $10^{3}$ (a factor of 10 larger than in the plots).

Signal efficiency in the SRA$^{[450,575]}$ region in the $\tilde{t}_1 - \tilde\chi^0_1$ mass plane, assuming BR$(\tilde{t}_1 \rightarrow t + \tilde\chi^0_1) = 0.5$.

Signal acceptance in the SRA$^{\geq 575}$ region in the $\tilde{t}_1 - \tilde\chi^0_1$ mass plane, assuming BR$(\tilde{t}_1 \rightarrow t + \tilde\chi^0_1) = 0.5$. Please mind that the acceptance given in the table is multiplied by factor of $10^{3}$ (a factor of 10 larger than in the plots).

Signal efficiency in the SRA$^{\geq 575}$ region in the $\tilde{t}_1 - \tilde\chi^0_1$ mass plane, assuming BR$(\tilde{t}_1 \rightarrow t + \tilde\chi^0_1) = 0.5$.

Signal acceptance in the SRB$^{[100,150]}$ region in the $\tilde{t}_1 - \tilde\chi^0_1$ mass plane, assuming BR$(\tilde{t}_1 \rightarrow t + \tilde\chi^0_1) = 0.5$. Please mind that the acceptance given in the table is multiplied by factor of $10^{3}$ (a factor of 10 larger than in the plots).

Signal efficiency in the SRB$^{[100,150]}$ region in the $\tilde{t}_1 - \tilde\chi^0_1$ mass plane, assuming BR$(\tilde{t}_1 \rightarrow t + \tilde\chi^0_1) = 0.5$.

Signal acceptance in the SRB$^{[150,400]}$ region in the $\tilde{t}_1 - \tilde\chi^0_1$ mass plane, assuming BR$(\tilde{t}_1 \rightarrow t + \tilde\chi^0_1) = 0.5$. Please mind that the acceptance given in the table is multiplied by factor of $10^{3}$ (a factor of 10 larger than in the plots).

Signal efficiency in the SRB$^{[150,400]}$ region in the $\tilde{t}_1 - \tilde\chi^0_1$ mass plane, assuming BR$(\tilde{t}_1 \rightarrow t + \tilde\chi^0_1) = 0.5$.

Signal acceptance in the SRB$^{\geq 400}$ region in the $\tilde{t}_1 - \tilde\chi^0_1$ mass plane, assuming BR$(\tilde{t}_1 \rightarrow t + \tilde\chi^0_1) = 0.5$. Please mind that the acceptance given in the table is multiplied by factor of $10^{3}$ (a factor of 10 larger than in the plots).

Signal efficiency in the SRB$^{\geq 400}$ region in the $\tilde{t}_1 - \tilde\chi^0_1$ mass plane, assuming BR$(\tilde{t}_1 \rightarrow t + \tilde\chi^0_1) = 0.5$.

Signal acceptance in the SRC$^{[100,150]}$ region in the $\tilde{t}_1 - \tilde\chi^0_1$ mass plane, assuming BR$(\tilde{t}_1 \rightarrow t + \tilde\chi^0_1) = 0.5$. Please mind that the acceptance given in the table is multiplied by factor of $10^{3}$ (a factor of 10 larger than in the plots).

Signal efficiency in the SRC$^{[100,150]}$ region in the $\tilde{t}_1 - \tilde\chi^0_1$ mass plane, assuming BR$(\tilde{t}_1 \rightarrow t + \tilde\chi^0_1) = 0.5$.

Signal acceptance in the SRC$^{[150,300]}$ region in the $\tilde{t}_1 - \tilde\chi^0_1$ mass plane, assuming BR$(\tilde{t}_1 \rightarrow t + \tilde\chi^0_1) = 0.5$. Please mind that the acceptance given in the table is multiplied by factor of $10^{3}$ (a factor of 10 larger than in the plots).

Signal efficiency in the SRC$^{[150,300]}$ region in the $\tilde{t}_1 - \tilde\chi^0_1$ mass plane, assuming BR$(\tilde{t}_1 \rightarrow t + \tilde\chi^0_1) = 0.5$.

Signal acceptance in the SRC$^{[300,500]}$ region in the $\tilde{t}_1 - \tilde\chi^0_1$ mass plane, assuming BR$(\tilde{t}_1 \rightarrow t + \tilde\chi^0_1) = 0.5$. Please mind that the acceptance given in the table is multiplied by factor of $10^{3}$ (a factor of 10 larger than in the plots).

Signal efficiency in the SRC$^{[300,500]}$ region in the $\tilde{t}_1 - \tilde\chi^0_1$ mass plane, assuming BR$(\tilde{t}_1 \rightarrow t + \tilde\chi^0_1) = 0.5$.

Signal acceptance in the SRC$^{\geq 500}$ region in the $\tilde{t}_1 - \tilde\chi^0_1$ mass plane, assuming BR$(\tilde{t}_1 \rightarrow t + \tilde\chi^0_1) = 0.5$. Please mind that the acceptance given in the table is multiplied by factor of $10^{3}$ (a factor of 10 larger than in the plots).

Signal efficiency in the SRC$^{\geq 500}$ region in the $\tilde{t}_1 - \tilde\chi^0_1$ mass plane, assuming BR$(\tilde{t}_1 \rightarrow t + \tilde\chi^0_1) = 0.5$.

Signal acceptance in the SRD750 region in the $\tilde{t}_1 - \tilde\chi^0_1$ mass plane, assuming BR$(\tilde{t}_1 \rightarrow t + \tilde\chi^0_1) = 0.5$. Please mind that the acceptance given in the table is multiplied by factor of $10^{3}$ (a factor of 10 larger than in the plots).

Signal efficiency in the SRD750 region in the $\tilde{t}_1 - \tilde\chi^0_1$ mass plane, assuming BR$(\tilde{t}_1 \rightarrow t + \tilde\chi^0_1) = 0.5$.

Signal acceptance in the SRD1000 region in the $\tilde{t}_1 - \tilde\chi^0_1$ mass plane, assuming BR$(\tilde{t}_1 \rightarrow t + \tilde\chi^0_1) = 0.5$. Please mind that the acceptance given in the table is multiplied by factor of $10^{3}$ (a factor of 10 larger than in the plots).

Signal efficiency in the SRD1000 region in the $\tilde{t}_1 - \tilde\chi^0_1$ mass plane, assuming BR$(\tilde{t}_1 \rightarrow t + \tilde\chi^0_1) = 0.5$.

Signal acceptance in the SRD1250 region in the $\tilde{t}_1 - \tilde\chi^0_1$ mass plane, assuming BR$(\tilde{t}_1 \rightarrow t + \tilde\chi^0_1) = 0.5$. Please mind that the acceptance given in the table is multiplied by factor of $10^{3}$ (a factor of 10 larger than in the plots).

Signal efficiency in the SRD1250 region in the $\tilde{t}_1 - \tilde\chi^0_1$ mass plane, assuming BR$(\tilde{t}_1 \rightarrow t + \tilde\chi^0_1) = 0.5$.

Signal acceptance in the SRD1500 region in the $\tilde{t}_1 - \tilde\chi^0_1$ mass plane, assuming BR$(\tilde{t}_1 \rightarrow t + \tilde\chi^0_1) = 0.5$. Please mind that the acceptance given in the table is multiplied by factor of $10^{3}$ (a factor of 10 larger than in the plots).

Signal efficiency in the SRD1500 region in the $\tilde{t}_1 - \tilde\chi^0_1$ mass plane, assuming BR$(\tilde{t}_1 \rightarrow t + \tilde\chi^0_1) = 0.5$.

Signal acceptance in the SRD1750 region in the $\tilde{t}_1 - \tilde\chi^0_1$ mass plane, assuming BR$(\tilde{t}_1 \rightarrow t + \tilde\chi^0_1) = 0.5$. Please mind that the acceptance given in the table is multiplied by factor of $10^{3}$ (a factor of 10 larger than in the plots).

Signal efficiency in the SRD1750 region in the $\tilde{t}_1 - \tilde\chi^0_1$ mass plane, assuming BR$(\tilde{t}_1 \rightarrow t + \tilde\chi^0_1) = 0.5$.

Signal acceptance in the SRD2000 region in the $\tilde{t}_1 - \tilde\chi^0_1$ mass plane, assuming BR$(\tilde{t}_1 \rightarrow t + \tilde\chi^0_1) = 0.5$. Please mind that the acceptance given in the table is multiplied by factor of $10^{3}$ (a factor of 10 larger than in the plots).

Signal efficiency in the SRD2000 region in the $\tilde{t}_1 - \tilde\chi^0_1$ mass plane, assuming BR$(\tilde{t}_1 \rightarrow t + \tilde\chi^0_1) = 0.5$.

Distribution of the reconstructed $\tau\tau$ invariant mass ($m_{\tau\tau}$) for all events in the VBF_1 signal region for $p_{\text{T}}^{H}>200$ GeV. The observed Higgs boson signal corresponds to $(\sigma\times B)/(\sigma\times B)_{\text{SM}}\,=\,0.99$. Entries with values above the $x$-axis range are shown in the last bin of each distribution. The prediction for each sample is determined from the likelihood fit performed to measure the total $pp\rightarrow H\rightarrow\tau\tau$ cross-section.

Distribution of the reconstructed $\tau\tau$ invariant mass ($m_{\tau\tau}$) for all events in the VBF_0 signal region for $m_{jj}<1$ TeV. The observed Higgs boson signal corresponds to $(\sigma\times B)/(\sigma\times B)_{\text{SM}}\,=\,0.99$. Entries with values above the $x$-axis range are shown in the last bin of each distribution. The prediction for each sample is determined from the likelihood fit performed to measure the total $pp\rightarrow H\rightarrow \tau\tau$ cross-section.

Distribution of the reconstructed $\tau\tau$ invariant mass ($m_{\tau\tau}$) for all events in the VBF_0 signal region for $m_{jj}>1$ TeV. The observed Higgs boson signal corresponds to $(\sigma\times B)/(\sigma\times B)_{\text{SM}}\,=\,0.99$. Entries with values above the $x$-axis range are shown in the last bin of each distribution. The prediction for each sample is determined from the likelihood fit performed to measure the total $pp\rightarrow H \rightarrow \tau\tau$ cross-section.

Distribution of the reconstructed $\tau\tau$ invariant mass ($m_{\tau\tau}$) for all events in the VBF_1 signal region for $m_{jj}<1$ TeV. The observed Higgs boson signal corresponds to $(\sigma\times B)/(\sigma\times B)_{\text{SM}}\,=\,0.99$. Entries with values above the $x$-axis range are shown in the last bin of each distribution. The prediction for each sample is determined from the likelihood fit performed to measure the total $pp\rightarrow \rightarrow \tau\tau$ cross-section.

Distribution of the reconstructed $\tau\tau$ invariant mass ($m_{\tau\tau}$) for all events in the VBF_1 signal region for $m_{jj}>1$ TeV. The observed Higgs boson signal corresponds to $(\sigma\times B)/(\sigma\times B)_{\text{SM}}\,=\,0.99$. Entries with values above the $x$-axis range are shown in the last bin of each distribution. The prediction for each sample is determined from the likelihood fit performed to measure the total $pp\rightarrow H\rightarrow \tau\tau$ cross-section.

Comparison of the observed and expected event yields in the ttH categories. The window category contains events with $m_{\tau\tau} \in [100, 150]$ GeV, while events outside of this mass region are included in the sideband region. The observed Higgs boson signal corresponds to $(\sigma\times B)/(\sigma\times B)_{\text{SM}}\,=\,0.99$. The prediction for each sample is determined from the likelihood fit performed to measure the total $pp\rightarrow H\rightarrow\tau\tau$ cross-section.

The measured values of $\sigma_{H}\times B(H\rightarrow\tau\tau)$ relative to the SM expectations in the total (labelled as 'Combined') and per-production-mode fit. The total $\pm1\sigma$ uncertainty in each measurement is indicated by the error bar. The total systematic uncertainty in each measurement is also indicated.

The measured values for $\sigma_{H}\times B(H\rightarrow\tau\tau)$ relative to the SM expectations in the simplified template cross-section measurement. The total $\pm1\sigma$ uncertainty in each measurement is indicated by the error bar. The total systematic uncertainty in each measurement is also indicated.

Upper limits at 95% CL on simplified template cross-sections in the individual ttH $p_{\text{T}}^{H}$ bins relative to their SM expectation in the $H\rightarrow\tau\tau$ decay channel, derived using the CL$_\text{s}$ method. The observed upper limits, together with the expected limits both in the background-only hypothesis and in the SM hypothesis (dotted red lines) are included. In the case of the expected limits in the background-only hypothesis, one- and two-standard-deviation uncertainty bands are also shown.

Measured fiducial differential cross-sections as a function of $p_\text{T}(j_0)$. The measurements are compared with particle-level SM predictions from the POWHEG+PYTHIA8, POWHEG+HERWIG7 and MadGraph5_aMC@NLO+Pythia8 generators for the combined VBF, ggF, $VH$, and $t\bar{t}H$ Higgs boson production modes. The shaded box around each data point shows the statistical uncertainty, while the total uncertainty is indicated by the error bar. The contribution from the ggF, $VH$, and $t\bar{t}H$ production modes as predicted by POWHEG+PYTHIA8 is also shown.

Measured fiducial differential cross-sections as a function of $p_\text{T}^{H}$. The measurements are compared with particle-level SM predictions from the POWHEG+PYTHIA8, POWHEG+HERWIG7 and MadGraph5_aMC@NLO+Pythia8 generators for the combined VBF, ggF, $VH$, and $t\bar{t} H$ Higgs boson production modes. The shaded box around each data point shows the statistical uncertainty, while the total uncertainty is indicated by the error bar. The contribution from the ggF, $VH$, and $t\bar{t} H$ production modes as predicted by POWHEG+PYTHIA8 is also shown.

Measured fiducial differential cross-sections as a function of $\Delta\phi_{jj}^{\text{signed}}$. The measurements are compared with particle-level SM predictions from the POWHEG+PYTHIA8, POWHEG+HERWIG7 and MadGraph5_aMC@NLO+Pythia8 generators for the combined VBF, ggF, $VH$, and $t\bar{t}H$ Higgs boson production modes. The shaded box around each data point shows the statistical uncertainty, while the total uncertainty is indicated by the error bar. The contribution from the ggF, $VH$, and $t\bar{t} H$ production modes as predicted by POWHEG+PYTHIA8 is also shown.

Measured fiducial differential cross-sections as a function of $\Delta\phi_{jj}^{\text{signed}}$ vs $p_\text{T}^\text{H}$. The measurements are compared with particle-level SM predictions from the POWHEG+PYTHIA8, POWHEG+HERWIG7 and MadGraph5_aMC@NLO+Pythia8 generators for the combined VBF, ggF, $VH$, and $t\bar{t}H$ Higgs boson production modes. The shaded box around each data point shows the statistical uncertainty, while the total uncertainty is indicated by the error bar. The contribution from the ggF, $VH$, and $t\bar{t}H$ production modes as predicted by POWHEG+PYTHIA8 is also shown.

The unfolded data compared with the nominal SM predictions from POWHEG+PYTHIA8, as well as the predictions for several values of Wilson coefficients in the SMEFT, for $\Delta\phi_{jj}^{\text{signed}}$.

The unfolded data compared with the nominal SM predictions from POWHEG+PYTHIA8, as well as the predictions for several values of Wilson coefficients in the SMEFT, for $\Delta\phi_{jj}^{\text{signed}}$.

The unfolded data compared with the nominal SM predictions from POWHEG+PYTHIA8, as well as the predictions for several values of Wilson coefficients in the SMEFT, for $\Delta\phi_{jj}^{\text{signed}}$ vs $p_\text{T}^{H}$.

The unfolded data compared with the nominal SM predictions from POWHEG+PYTHIA8, as well as the predictions for several values of Wilson coefficients in the SMEFT, for $p_\text{T}^{H}$.

The 95% expected and observed confidence intervals for each of the six considered Wilson coefficients. Each coefficient is treated individually and both the linear and linear + quadratic models are considered when evaluating the sensitivity to the terms that enter at $\mathcal{O}(\Lambda^{-4})$. For the CP-even operators the $\Delta\phi_{jj}^{\text{signed}}$ distribution is used to extract the confidence interval. The limits are computed at a new-physics scale $\Lambda=1$ TeV.

The 95% expected and observed confidence intervals for each of the six considered Wilson coefficients. Each coefficient is treated individually and both the linear and linear + quadratic models are considered when evaluating the sensitivity to the terms that enter at $\mathcal{O}(\Lambda^{-4})$. For the CP-odd operators the $\Delta\phi_{jj}^{\text{signed}}$ vs $p_\text{T}^{H}$ distribution is used to extract the confidence interval. The limits are computed at a new-physics scale $\Lambda=1$ TeV.

Two-dimensional observed and expected limits at 68% and 95% confidence level obtained from $c_{HW}$ and $c_{HB}$ using only the SM-dimension-6 interference in the SMEFT basis. The limits are computed at a new-physics scale $\Lambda$ = 1 TeV. In this case, the $\Delta\phi_{jj}^{signed}$ distribution is used. IMPORTANT: HEPData requires the same number of rows for each column. However, each curve of on the plot has a different number of points. For that reason, the curves are padded with data of "0.0, 0.0". The "0.0, 0.0" points should be omitted in analysis.

Two-dimensional observed and expected limits at 68% and 95% confidence level obtained from $c_{H\tilde{W}}$ and $c_{H\tilde{B}}$ using only the SM-dimension-6 interference in the SMEFT basis. The limits are computed at a new-physics scale $\Lambda$ = 1 TeV. In this purely CP-odd case, the $\Delta\phi_{jj}^{signed}$ vs $p_{T}^{H}$ distribution is used. IMPORTANT: HEPData requires the same number of rows for each column. However, each curve of on the plot has a different number of points. For that reason, the curves are padded with data of "0.0, 0.0". The "0.0, 0.0" points should be omitted in analysis.

Two-dimensional observed and expected limits at 68% and 95% confidence level obtained from $c_{HW}$ and $c_{H\tilde{W}}$ using only the SM-dimension-6 interference in the SMEFT basis. The limits are computed at a new-physics scale $\Lambda$ = 1 TeV. In this case, the $\Delta\phi_{jj}^{signed}$ distribution is used. IMPORTANT: HEPData requires the same number of rows for each column. However, each curve of on the plot has a different number of points. For that reason, the curves are padded with data of "0.0, 0.0". The "0.0, 0.0" points should be omitted in analysis.

The expected values of $\sigma_{H}\times B(H\rightarrow\tau\tau)$ relative to the SM expectations in the total (labelled as 'Combined') and per-production-mode fit. The total $\pm1\sigma$ uncertainty in each measurement is indicated by the error bar. The total systematic uncertainty in each measurement is also indicated.

The expected values for $\sigma_{H}\times B(H\rightarrow\tau\tau)$ relative to the SM expectations in the simplified template cross-section measurement. The total $\pm1\sigma$ uncertainty in each measurement is indicated by the error bar. The total systematic uncertainty in each measurement is also indicated.

The measured values for $\sigma_{H}\times B(H\rightarrow\tau\tau)$ in the simplified template cross-section measurement. The total $\pm1\sigma$ uncertainty in the measurement is indicated by the black error bars, with the individual contribution from the statistical uncertainty in blue.

The measured correlations between differential fiducial cross-sections, $\sigma$, and the normalisation factors, $\alpha$, for the $Z(\to\tau\tau) + \text{jets}$ and top backgrounds for $p_\text{T}(j_0)$. The bins denoted as $\sigma_{1,2,3,4}$ correspond to the measurements in the four bins listed elsewhere for each observable.

The measured correlations between differential fiducial cross-sections, $\sigma$, and the normalisation factors, $\alpha$, for the $Z(\to\tau\tau) + \text{jets}$ and top backgrounds for $p_\text{T}^{H}}$. The bins denoted as $\sigma_{1,2,3,4}$ correspond to the measurements in the four bins listed elsewhere for each observable.

The measured correlations between differential fiducial cross-sections, $\sigma$, and the normalisation factors, $\alpha$, for the $Z(\to\tau\tau) + \text{jets}$ and top backgrounds for $\Delta\phi_{jj}^{\text{signed}}$. The bins denoted as $\sigma_{1,2,3,4}$ correspond to the measurements in the four bins listed elsewhere for each observable.

The measured correlations between differential fiducial cross-sections, $\sigma$, and the normalisation factors, $\alpha$, for the $Z(\to\tau\tau) + \text{jets}$ and top backgrounds for $\Delta\phi_{jj}^{\text{signed}}$ vs $p_\text{T}^{H}$. The bins denoted as $\sigma_{1,2,3,4}$ correspond to the measurements in the four bins listed elsewhere for each observable.

Figure 8(left) of the article. Differential fiducial cross-section of the Z boson $p_T$ in events with $Z \left( \rightarrow \ell\ell \right) + 1 b$-jet. The thin inner band corresponds to the statistical uncertainty of the data, and the outer band to statistical and systematic uncertainties of the data, added in quadrature.

Figure 8(right) of the article. Differential fiducial cross-section of the leading $b$-jet $p_T$ in events with $Z \left( \rightarrow \ell\ell \right) + 1 b$-jet. The thin inner band corresponds to the statistical uncertainty of the data, and the outer band to statistical and systematic uncertainties of the data, added in quadrature.

Figure 9 of the article. Differential fiducial cross-section of the $\Delta R$ between the $Z$ boson and the leading $b$-jet in events with $Z \left( \rightarrow \ell\ell \right) + 1 b$-jet. The thin inner band corresponds to the statistical uncertainty of the data, and the outer band to statistical and systematic uncertainties of the data, added in quadrature.

Figure 10(left) of the article. Differential fiducial cross-section of the $\Delta \phi$ between the leading and sub-leading $b$-jets in events with $Z \left( \rightarrow \ell\ell \right) + 2 b$-jets. The thin inner band corresponds to the statistical uncertainty of the data, and the outer band to statistical and systematic uncertainties of the data, added in quadrature.

Figure 10(right) of the article. Differential fiducial cross-section of the invariant mass of the leading and sub-leading $b$-jets in events with $Z \left( \rightarrow \ell\ell \right) + 2 b$-jets. The thin inner band corresponds to the statistical uncertainty of the data, and the outer band to statistical and systematic uncertainties of the data, added in quadrature.

Figure 11(left) of the article. Differential fiducial cross-section of the Z boson $p_T$ in events with $Z \left( \rightarrow \ell\ell \right) + 1 c$-jet. The thin inner band corresponds to the statistical uncertainty of the data, and the outer band to statistical and systematic uncertainties of the data, added in quadrature.

Figure 11(right) of the article. Differential fiducial cross-section of the leading $c$-jet $p_T$ in events with $Z \left( \rightarrow \ell\ell ight) + 1 c$-jet. The thin inner band corresponds to the statistical uncertainty of the data, and the outer band to statistical and systematic uncertainties of the data, added in quadrature.

Figure 12 of the article. Differential fiducial cross-section of the leading $c$-jet $x_F$ in events with $Z \left( \rightarrow \ell\ell \right) + 1 c$-jet. The thin inner band corresponds to the statistical uncertainty of the data, and the outer band to statistical and systematic uncertainties of the data, added in quadrature.

Correction factors from parton level done with flavour dressing algorithms as seen in Ref. [2] of the paper to hadron level with R 0.3 cone labelling (see Tables 18 & 24) for the differential cross-section of the Z boson $p_T$ in events with $Z (\rightarrow ll) \ge 1 $ b-jets. Systematic errors include Sherpa 2.2.11 and MGaMC+Py8 FxFx uncertainties for what concerns the passage from hadron level flavour dressing to hadron level cone labelling, plus statistical uncertainty of Pythia8 simulations used for parton to hadron corrections. See Table 4.

Correction factors from parton level done with flavour dressing algorithms as seen in Ref. [2] of the paper to hadron level with R 0.3 cone labelling (see Tables 19 & 25) for the differential cross-section of the Z boson $p_T$ in events with $Z (\rightarrow ll) \ge 1 $ c-jets. Systematic errors include Sherpa 2.2.11 and MGaMC+Py8 FxFx uncertainties for what concerns the passage from hadron level flavour dressing to hadron level cone labelling, plus statistical uncertainty of Pythia8 simulations used for parton to hadron corrections. See Table 9.

Correction factors from parton level done with flavour dressing algorithms as seen in Ref. [2] of the paper to hadron level with R 0.3 cone labelling (see Tables 20 & 26) for the differential cross-section of the leading $b$-jet $p_T$ in events with $Z (\rightarrow ll) \ge 1 $ b-jets. Systematic errors include Sherpa 2.2.11 and MGaMC+Py8 FxFx uncertainties for what concerns the passage from hadron level flavour dressing to hadron level cone labelling, plus statistical uncertainty of Pythia8 simulations used for parton to hadron corrections. See Table 5.

Correction factors from parton level done with flavour dressing algorithms as seen in Ref. [2] of the paper to hadron level with R 0.3 cone labelling (see Tables 21 & 27) for the differential cross-section of the leading $c$-jet $p_T$ in events with $Z (\rightarrow ll) \ge 1 $ c-jets. Systematic errors include Sherpa 2.2.11 and MGaMC+Py8 FxFx uncertainties for what concerns the passage from hadron level flavour dressing to hadron level cone labelling, plus statistical uncertainty of Pythia8 simulations used for parton to hadron corrections. See Table 10.

Correction factors from parton level done with flavour dressing algorithms as seen in Ref. [2] of the paper to hadron level with R 0.3 cone labelling (see Tables 22 & 28) for the differential cross-section of the leading $c$-jet $x_F$ in events with $Z (\rightarrow ll) \ge 1 $ c-jets. Systematic errors include Sherpa 2.2.11 and MGaMC+Py8 FxFx uncertainties for what concerns the passage from hadron level flavour dressing to hadron level cone labelling, plus statistical uncertainty of Pythia8 simulations used for parton to hadron corrections. See Table 11.

Correction factors from parton level done with flavour dressing algorithms as seen in Ref. [2] of the paper of the paper to hadron level with R 0.3 cone labelling (see Tables 23 & 29) for the differential cross-section of $\Delta R$ between the $Z$ boson and the leading $b$-jet in events with $Z (\rightarrow ll) \ge 1 $ b-jets. Systematic errors include Sherpa 2.2.11 and MGaMC+Py8 FxFx uncertainties for what concerns the passage from hadron level flavour dressing to hadron level cone labelling, plus statistical uncertainty of Pythia8 simulations used for parton to hadron corrections. See Table 6.

Correction factors from hadron level done with flavour dressing algorithms as seen in Ref. [2] of the paper to hadron level R 0.3 cone labelling for the differential cross-section of the Z boson $p_T$ in events with $Z (\rightarrow ll) \ge 1 $ b-jets. Systematic errors include Sherpa 2.2.11 and MGaMC+Py8 FxFx uncertainties. See Table 12.

Correction factors from hadron level done with flavour dressing algorithms as seen in Ref. [2] of the paper to hadron level with R 0.3 cone labelling for the differential cross-section of the Z boson $p_T$ in events with $Z (\rightarrow ll) \ge 1 $ c-jets. Systematic errors include Sherpa 2.2.11 and MGaMC+Py8 FxFx uncertainties. See Table 13.

Correction factors from hadron level done with flavour dressing algorithms as seen in Ref. [2] of the paper to hadron level with R 0.3 cone labelling for the differential cross-section of the leading $b$-jet $p_T$ in events with $Z (\rightarrow ll) \ge 1 $ b-jets. Systematic errors include Sherpa 2.2.11 and MGaMC+Py8 FxFx uncertainties. See Table 14.

Correction factors from hadron level done with flavour dressing algorithms as seen in Ref. [2] of the paper to hadron level with R 0.3 cone labelling for the differential cross-section of the leading $c$-jet $p_T$ in events with $Z (\rightarrow ll) \ge 1 $ c-jets. Systematic errors include Sherpa 2.2.11 and MGaMC+Py8 FxFx uncertainties. See Table 15.

Correction factors from hadron level done with flavour dressing algorithms as seen in Ref. [2] of the paper to hadron level with R 0.3 cone labelling for the differential cross-section of the leading $c$-jet $x_F$ in events with $Z (\rightarrow ll) \ge 1 $ c-jets. Systematic errors include Sherpa 2.2.11 and MGaMC+Py8 FxFx uncertainties. See Table 16.

Correction factors from hadron level done with flavour dressing algorithms as seen in Ref. [2] of the paper of the paper to hadron level with R 0.3 cone labelling for the differential cross-section of $\Delta R$ between the $Z$ boson and the leading $b$-jet in events with $Z (\rightarrow ll) \ge 1 $ b-jets. Systematic errors include Sherpa 2.2.11 and MGaMC+Py8 FxFx uncertainties. See Table 17.

Correction factors from parton level to hadron level for the differential cross-section of the Z boson $p_T$ in events with $Z (\rightarrow ll) \ge 1 $ b-jets. Systematic errors statistical uncertainty of Pythia8 simulations. See Table 12.

Correction factors from parton level to hadron level for the differential cross-section of the Z boson $p_T$ in events with $Z (\rightarrow ll) \ge 1 $ c-jets. Systematic errors include statistical uncertainty of Pythia8 simulations. See Table 13.

Correction factors from parton level to hadron level for the differential cross-section of the leading $b$-jet $p_T$ in events with $Z (\rightarrow ll) \ge 1 $ b-jets. Systematic errors include statistical uncertainty of Pythia8 simulations. See Table 14.

Correction factors from parton level to hadron level for the differential cross-section of the leading $c$-jet $p_T$ in events with $Z (\rightarrow ll) \ge 1 $ c-jets. Systematic errors include statistical uncertainty of Pythia8 simulations. See Table 15.

Correction factors from parton level to hadron level for the differential cross-section of the leading $c$-jet $x_F$ in events with $Z (\rightarrow ll) \ge 1 $ c-jets. Systematic errors include statistical uncertainty of Pythia8 simulations. See Table 16.

Correction factors from parton level to hadron level for the differential cross-section of $\Delta R$ between the $Z$ boson and the leading $b$-jet in events with $Z (\rightarrow ll) \ge 1 $ b-jets. Systematic errors include statistical uncertainty of Pythia8 simulations. See Table 17.

Post-fit $m_\mathrm{MMC}$ distribution in the high-mass SR for the $m_A = 90\,\mathrm{GeV}$ signal mass hypothesis. $m_\mathrm{MMC}$ is the mass reconstructed by the Missing Mass Calculator. Processes contributing to the background Others are $Z/\gamma^* \rightarrow ee/\mu\mu$ and SM Higgs. The subscript on the $A\to\tau\tau$ process indicates the mass of the $A$ boson. otal includes all backgrounds and the signal process. The high-mass Signal Region is defined as: - 1 electron and 1 muon with opposite charge - $p_\mathrm{T}$ requirements of the leptons are a combination of the following: - $p_\mathrm{T}^e > 18\,\mathrm{GeV}$ and $p_\mathrm{T}^\mu > 15\,\mathrm{GeV}$ or - $p_\mathrm{T}^e > 10\,\mathrm{GeV}$ and $p_\mathrm{T}^\mu > 25\,\mathrm{GeV}$ or - $p_\mathrm{T}^e > 27\,\mathrm{GeV}$ and $p_\mathrm{T}^\mu > 10\,\mathrm{GeV}$ - $\vert \eta_e \vert < 2.47$, excluding $1.37 < \vert \eta_e \vert < 1.52$ - $\vert \eta_\mu \vert < 2.7$ - no jets with $b$-quarks - $\Delta R_{\ell\ell} < 1.0$ - $E_\mathrm{T}^\mathrm{miss} > 30\,\mathrm{GeV}$ - $m_\mathrm{T}^\mathrm{tot} = \sqrt{\left(p_\mathrm{T}^e+p_\mathrm{T}^\mu+E_\mathrm{T}^\mathrm{miss}\right)^2-\left(\vec{p}_\mathrm{T}^{\,e}+\vec{p}_\mathrm{T}^{\,\mu}+\vec{E}_\mathrm{T}^{\,\mathrm{miss}}\right)^2} < 65\,\mathrm{GeV}$ - $35\,\mathrm{GeV} < m_\mathrm{MMC} < 130\,\mathrm{GeV}$

Expected and observed $95\%$ CL limits on the production cross-section for $gg\rightarrow A$ times the branching ratio for $A$ decaying into two $\tau$-leptons for $A$ boson masses ranging from $20$ to $90\,\mathrm{GeV}$.

Expected and observed $95\%$ CL limits on the production cross-section for $gg\rightarrow A$ times the branching ratio for $A$ decaying into two $\tau$-leptons for $A$ boson masses ranging from $20$ to $90\,\mathrm{GeV}$.

Expected and observed $95\%$ CL limits on the absolute value of the up-type quark coupling parameter $|\zeta_\mathrm{u}|$ for masses of the $A$ boson ranging from $20$ to $90\,\mathrm{GeV}$.

Expected and observed $95\%$ CL limits on the absolute value of the up-type quark coupling parameter $|\zeta_\mathrm{u}|$ for masses of the $A$ boson ranging from $20$ to $90\,\mathrm{GeV}$.

Observed $p_0$ values for the fits on the production cross-section for $gg\rightarrow A$ times the branching ratio for $A$ decaying into two $\tau$-leptons for $A$ boson masses from $20$ to $90\,\mathrm{GeV}$. For the fit on the up-type quark coupling parameter the $p_0$ values agree up to a precision of $O(10^{-5})$.

Observed $p_0$ values for the fits on the production cross-section for $gg\rightarrow A$ times the branching ratio for $A$ decaying into two $\tau$-leptons for $A$ boson masses from $20$ to $90\,\mathrm{GeV}$. For the fit on the up-type quark coupling parameter the $p_0$ values agree up to a precision of $O(10^{-5})$.

Pre-fit background composition of the SR$3\ell$ Prod. The table shows the event yields as opposed to just the percentages of the relevant background processes.

Post-fit plot of $H_\text{T}(\text{jets})$ in the SR$2\ell$ Dec from a signal-blinded background-only fit.

Post-fit plot of $m(t_\text{SM}, b\text{-jet}_0)$ in the SR$2\ell$ Prod from a signal-blinded background-only fit.

Post-fit plot of $m(\ell_\text{OS},\ell_\text{SS,1})$ in the SR$3\ell$ Dec from a signal-blinded background-only fit.

Post-fit plot of $m(\ell_\text{OS},\ell_\text{SS,1})$ in the SR$3\ell$ Prod from a signal-blinded background-only fit.

Post-fit plot of $D_\text{NN}(tHc)$ in the SR$2\ell$ Dec from the full fit to data.

Post-fit plot of $D_\text{NN}(tHc)$ in the SR$2\ell$ Prod from the full fit to data.

Post-fit plot of $D_\text{NN}(tHc)$ in the SR$3\ell$ Dec from the full fit to data.

Post-fit plot of $D_\text{NN}(tHc)$ in the SR$3\ell$ Prod from the full fit to data.

Post-fit plot of $p_\text{T}(\ell_1)$ in the CR$2\ell$ HF$e$ from the full fit to data.

Post-fit plot of $p_\text{T}(\ell_1)$ in the CR$2\ell$ HF$\mu$ from the full fit to data.

Post-fit plot of $p_\text{T}(\ell_1)$ in the CR$2\ell$ $t\bar{t}V$ from the full fit to data.

Post-fit plot of $p_\text{T}(\ell_2)$ in the CR$3\ell$ HF$e$ from the full fit to data.

Post-fit plot of $p_\text{T}(\ell_2)$ in the CR$3\ell$ HF$\mu$ from the full fit to data.

Post-fit plot of $p_\text{T}(b\text{-jet}_0)$ in the CR$3\ell$ $t\bar{t}W$ from the full fit to data.

Post-fit plot of $p_\text{T}(b\text{-jet}_0)$ in the CR$3\ell$ $t\bar{t}Z$ from the full fit to data.

Observed and expected upper exclusion limits on the branching ratio $\mathcal{B}(t\to Hu)$ for different analyses and their statistical combination.

Observed and expected upper exclusion limits on the branching ratio $\mathcal{B}(t\to Hc)$ for different analyses and their statistical combination.

Post-fit normalisation factors of free-floating background processes and the signal normalisation.

Post-fit predicted and observed yields in all $2\ell$SS signal and control regions. Pre-fit signal contributions for a signal cross section equivalent to $\mathcal{B}(t\to Hq)=0.1\,\%$ are given as well.

Post-fit predicted and observed yields in all $3\ell$ signal and control regions. Pre-fit signal contributions for a signal cross section equivalent to $\mathcal{B}(t\to Hq)=0.1\,\%$ are given as well.

Expected upper limits on $\mathcal{B}(t\to Hq)$ for the nominal fit and alternative fit configurations. One contains the full phase space but only considers statistical uncertainties. Two other configurations consider the full set of systematic uncertainties, but only encompass one final state.

Expected and observed upper limits on $\mathcal{B}(t\to Hq)$ and $|C_{u\phi}^{qt,tq}|$ for the full fit containing all systematic uncertainties.

Pre-fit plot of $H_\text{T}(\text{jets})$ in the SR$2\ell$ Dec from a signal-blinded background-only fit.

Pre-fit plot of $m(t_\text{SM}, b\text{-jet}_0)$ in the SR$2\ell$ Prod from a signal-blinded background-only fit.

Pre-fit plot of $m(\ell_\text{OS},\ell_\text{SS,1})$ in the SR$3\ell$ Dec from a signal-blinded background-only fit.

Pre-fit plot of $m(\ell_\text{OS},\ell_\text{SS,1})$ in the SR$3\ell$ Prod from a signal-blinded background-only fit.

Post-fit plot of $D_\text{NN}(tHu)$ in the SR$2\ell$ Dec from the full fit to data.

Post-fit plot of $D_\text{NN}(tHu)$ in the SR$2\ell$ Prod from the full fit to data.

Post-fit plot of $D_\text{NN}(tHu)$ in the SR$3\ell$ Dec from the full fit to data.

Post-fit plot of $D_\text{NN}(tHu)$ in the SR$3\ell$ Prod from the full fit to data.

Pre-fit plot of $D_\text{NN}(tHc)$ in the SR$2\ell$ Dec from the full fit to data.

Pre-fit plot of $D_\text{NN}(tHc)$ in the SR$2\ell$ Prod from the full fit to data.

Pre-fit plot of $D_\text{NN}(tHc)$ in the SR$3\ell$ Dec from the full fit to data.

Pre-fit plot of $D_\text{NN}(tHc)$ in the SR$3\ell$ Prod from the full fit to data.

Pre-fit plot of $D_\text{NN}(tHu)$ in the SR$2\ell$ Dec from the full fit to data.

Pre-fit plot of $D_\text{NN}(tHu)$ in the SR$2\ell$ Prod from the full fit to data.

Pre-fit plot of $D_\text{NN}(tHu)$ in the SR$3\ell$ Dec from the full fit to data.

Pre-fit plot of $D_\text{NN}(tHu)$ in the SR$3\ell$ Prod from the full fit to data.

Pre-fit plot of $p_\text{T}(\ell_1)$ in the CR$2\ell$ HF$e$ from the full fit to data.

Pre-fit plot of $p_\text{T}(\ell_1)$ in the CR$2\ell$ HF$\mu$ from the full fit to data.

Pre-fit plot of $p_\text{T}(\ell_1)$ in the CR$2\ell$ $t\bar{t}V$ from the full fit to data.

Pre-fit plot of $p_\text{T}(\ell_2)$ in the CR$3\ell$ HF$e$ from the full fit to data.

Pre-fit plot of $p_\text{T}(\ell_2)$ in the CR$3\ell$ HF$\mu$ from the full fit to data.

Pre-fit plot of $p_\text{T}(b\text{-jet}_0)$ in the CR$3\ell$ $t\bar{t}W$ from the full fit to data.

Pre-fit plot of $p_\text{T}(b\text{-jet}_0)$ in the CR$3\ell$ $t\bar{t}Z$ from the full fit to data.

Ranking of fit nuisance parameters according to their impact on the post-fit $tHu$ signal normalisation when fixed to $\pm1\sigma$

Ranking of fit nuisance parameters according to their impact on the post-fit $tHc$ signal normalisation when fixed to $\pm1\sigma$

Expected upper exclusion limits on the branching ratio $\mathcal{B}(t\to Hu)$ for each individual final state and the full analysis.

Expected upper exclusion limits on the branching ratio $\mathcal{B}(t\to Hc)$ for each individual final state and the full analysis.

The observed and expected ($\pm2\sigma$) upper limits at 95% CL on the branching ratio for $H\rightarrow aa\rightarrow \gamma\gamma\tau\tau$ as a function of the resonance mass hypothesis $m_{a}$.

The polynomial parameterisation of the total signal selection efficiency as a function of $m_{a}$, defined as the fraction of the reconstructed and selected events out of the total number of events in the signal simulated samples.

Local observed significance for signal discovery at different (mX, mS). The points show where the significance was evaluated. The band at mS = 125 GeV is not shown as those points are equivalent to those already probed in Phys. Rev. D 106 (2022) 052001.

Signal region yield for each of the simulated signal samples and interpolated signal points that were analyzed in the 1 b-jet Signal Region, normalized to a cross section of 1 fb. The signal region yields for one of the signals is found in Table 2 of the paper.

Signal region yield for each of the simulated signal samples and interpolated signal points that were analyzed in the 2 b-jet Signal Region, normalized to a cross section of 1 fb. The signal region yields for one of the signals is found in Table 2 of the paper.