Observed two-dimensional likelihood scans for $c_{Hq}^{(1)}$ vs. $c_{Hq}^{(3)}$ while other coefficients are fixed at their SM values.

Observed two-dimensional likelihood scans for $c_{Hq}^{(1)}$ vs. $c_{Hu}$ while allowing the other coefficients to float freely at each point of the sca.

Observed two-dimensional likelihood scans for $c_{Hq}^{(1)}$ vs. $c_{Hu}$ while other coefficients are fixed at their SM values.

Observed two-dimensional likelihood scans for $c_{Hq}^{(1)}$ vs. $c_{Hd}$ while allowing the other coefficients to float freely at each point of the sca.

Observed two-dimensional likelihood scans for $c_{Hq}^{(1)}$ vs. $c_{Hd}$ while other coefficients are fixed at their SM values.

Observed two-dimensional likelihood scans for $c_{Hq}^{(3)}$ vs. $c_{Hu}$ while allowing the other coefficients to float freely at each point of the sca.

Observed two-dimensional likelihood scans for $c_{Hq}^{(3)}$ vs. $c_{Hu}$ while other coefficients are fixed at their SM values.

Observed two-dimensional likelihood scans for $c_{Hq}^{(3)}$ vs. $c_{Hd}$ while allowing the other coefficients to float freely at each point of the sca.

Observed two-dimensional likelihood scans for $c_{Hq}^{(3)}$ vs. $c_{Hd}$ while other coefficients are fixed at their SM values.

Observed two-dimensional likelihood scans for $c_{Hu}$ vs. $c_{Hd}$ while allowing the other coefficients to float freely at each point of the sca.

Observed two-dimensional likelihood scans for $c_{Hu}$ vs. $c_{Hd}$ while other coefficients are fixed at their SM values.

Observed two-dimensional likelihood scans for $c_{Hq}^{(1)}$ vs. $g^{(2)}_{ZZ}$ while allowing the other coefficients to float freely at each point of the sca.

Observed two-dimensional likelihood scans for $c_{Hq}^{(1)}$ vs. $g^{(2)}_{ZZ}$ while other coefficients are fixed at their SM values.

Observed two-dimensional likelihood scans for $c_{Hq}^{(3)}$ vs. $g^{(2)}_{ZZ}$ while allowing the other coefficients to float freely at each point of the sca.

Observed two-dimensional likelihood scans for $c_{Hq}^{(3)}$ vs. $g^{(2)}_{ZZ}$ while other coefficients are fixed at their SM values.

Observed two-dimensional likelihood scans for $c_{Hu}$ vs. $g^{(2)}_{ZZ}$ while allowing the other coefficients to float freely at each point of the sca.

Observed two-dimensional likelihood scans for $c_{Hu}$ vs. $g^{(2)}_{ZZ}$ while other coefficients are fixed at their SM values.

Observed two-dimensional likelihood scans for $c_{Hq}^{(1)}$ vs. $g^{(4)}_{ZZ}$ while allowing the other coefficients to float freely at each point of the sca.

Observed two-dimensional likelihood scans for $c_{Hq}^{(1)}$ vs. $g^{(4)}_{ZZ}$ while other coefficients are fixed at their SM values.

Observed two-dimensional likelihood scans for $c_{Hq}^{(3)}$ vs. $g^{(4)}_{ZZ}$ while allowing the other coefficients to float freely at each point of the sca.

Observed two-dimensional likelihood scans for $c_{Hq}^{(3)}$ vs. $g^{(4)}_{ZZ}$ while other coefficients are fixed at their SM values.

Observed two-dimensional likelihood scans for $c_{Hu}$ vs. $g^{(4)}_{ZZ}$ while allowing the other coefficients to float freely at each point of the sca.

Observed two-dimensional likelihood scans for $c_{Hu}$ vs. $g^{(4)}_{ZZ}$ while other coefficients are fixed at their SM values.

Observed two-dimensional likelihood scans for $c_{Hd}$ vs. $g^{(2)}_{ZZ}$ while allowing the other coefficients to float freely at each point of the sca.

Observed two-dimensional likelihood scans for $c_{Hd}$ vs. $g^{(2)}_{ZZ}$ while other coefficients are fixed at their SM values.

Observed two-dimensional likelihood scans for $c_{Hd}$ vs. $g^{(4)}_{ZZ}$ while allowing the other coefficients to float freely at each point of the sca.

Observed two-dimensional likelihood scans for $c_{Hd}$ vs. $g^{(4)}_{ZZ}$ while other coefficients are fixed at their SM values.

Observed two-dimensional likelihood scans for $g^{(2)}_{ZZ}$ vs. $g^{(4)}_{ZZ}$ while allowing the other coefficients to float freely at each point of the sca.

Observed two-dimensional likelihood scans for $g^{(2)}_{ZZ}$ vs. $g^{(4)}_{ZZ}$ while other coefficients are fixed at their SM values.

Distributions in $m_{Z'}^{rec}$ for data in the SRs, together with fits of the background functions under the background-only hypothesis for the muon channel. The number of observed events in each bin is divided by the bin width. The signal predictions are shown for different Z' boson masses, normalized to an arbitrary cross section of 1 fb. In the panels below the distributions, the ratios of data to the background function are displayed. The shaded green areas represent the statistical uncertainty from the fit. The $\chi^2$ values per number of degrees of freedom ($\chi^2$/n.d.f.) and the corresponding $p$-values are provided for each fit.

Distributions in $m_{Z'}^{T}$ for data in the SRs, together with fits of the background functions under the background-only hypothesis for the invisible channel. The number of observed events in each bin is divided by the bin width. The signal predictions are shown for different Z' boson masses, normalized to an arbitrary cross section of 1 fb. In the panels below the distributions, the ratios of data to the background function are displayed. The shaded green areas represent the statistical uncertainty from the fit. The $\chi^2$ values per number of degrees of freedom ($\chi^2$/n.d.f.) and the corresponding $p$-values are provided for each fit.

Distributions in $m_{Z'}^{T}$ for data in the SRs, together with fits of the background functions under the background-only hypothesis for the invisible channel. The number of observed events in each bin is divided by the bin width. The signal predictions are shown for different Z' boson masses, normalized to an arbitrary cross section of 1 fb. In the panels below the distributions, the ratios of data to the background function are displayed. The shaded green areas represent the statistical uncertainty from the fit. The $\chi^2$ values per number of degrees of freedom ($\chi^2$/n.d.f.) and the corresponding $p$-values are provided for each fit.

Distributions in $m_{Z'}^{rec}$ for data in the SRs, together with fits of the background functions under the background-only hypothesis for the electron channel. The number of observed events in each bin is divided by the bin width. The signal predictions are shown for different Z' boson masses, normalized to an arbitrary cross section of 1 fb. In the panels below the distributions, the ratios of data to the background function are displayed. The shaded green areas represent the statistical uncertainty from the fit. The $\chi^2$ values per number of degrees of freedom ($\chi^2$/n.d.f.) and the corresponding $p$-values are provided for each fit.

Distributions in $m_{Z'}^{rec}$ for data in the SRs, together with fits of the background functions under the background-only hypothesis for the electron channel. The number of observed events in each bin is divided by the bin width. The signal predictions are shown for different Z' boson masses, normalized to an arbitrary cross section of 1 fb. In the panels below the distributions, the ratios of data to the background function are displayed. The shaded green areas represent the statistical uncertainty from the fit. The $\chi^2$ values per number of degrees of freedom ($\chi^2$/n.d.f.) and the corresponding $p$-values are provided for each fit.

Expected upper limits at 95% CL on the product of the production cross section and the branching fraction as a function of the Z' boson mass.Each line corresponds to the three different final states and their combined result.

Expected upper limits at 95% CL on the product of the production cross section and the branching fraction as a function of the Z' boson mass.Each line corresponds to the three different final states and their combined result.

Expected and observed upper limits at 95% CL on the product of the production cross section and the branching fraction as functions of the Z' boson mass from the combination of all final states. The limits are compared with the predictions of the HVT model and the expected limits (shown by red curves) from a previous analysis.

Expected and observed upper limits at 95% CL on the product of the production cross section and the branching fraction as functions of the Z' boson mass from the combination of all final states. The limits are compared with the predictions of the HVT model and the expected limits (shown by red curves) from a previous analysis.

Prediction from the HVT model.

Prediction from the HVT model.

Observed upper limits at 95% CL on gF for different Z' boson masses as functions of the product of $g_{H}$ with the sign of $g_{F}$. The two benchmark scenarios of the HVT model are shown by the black markers.

Observed upper limits at 95% CL on gF for different Z' boson masses as functions of the product of $g_{H}$ with the sign of $g_{F}$. The two benchmark scenarios of the HVT model are shown by the black markers.

Correlations of the ttH signal-strength modifiers in the different Higgs pT bins of the STXS measurement.

Observed and expected upper 95% CL limit on the tH signal strength modifier $\mu_{tH}$ with $\mu_{ttH}=1$ for different channels and years, where the uncertainties are uncorrelated between the channels and years, and in their combination. The one and two standard deviation confidence intervals on the expected limit are also listed (sigma).

Likelihood-ratio test statistic as a function of the ttH and tH signal strength modifiers $\mu_{ttH}$ and $\mu_{tH}$. The observed best fit point is $(\mu_{tH}, \mu_{ttH}) = (-3.83, +0.35)$.

Observed and expected likelihood ratio test statistic as a function of $\kappa_{t}$ and $\kappa_{V}$. The observed best fit point is $(\kappa_{t}, \kappa_{V}) = (+0.59, +1.40)$.

Observed and expected likelihood ratio test statistic as a function of $\kappa_{t}$ for $\kappa_{V}=1$. The observed best fit point is $\kappa_{t} = +0.54$.

Observed and expected likelihood ratio test statistic as a function of $\kappa_{t}$ and $\tilde{\kappa}_{t}$ with $\kappa_{V}=1$. The observed best fit point is $(\kappa_{t}, \tilde{\kappa}_{t}) = (+0.53, 0.00)$.

Observed and expected likelihood ratio test statistic as a function of $f_{CP}$. The overall modifier of the top-Higgs coupling strength is profiled such that the likelihood attains its minimum at each point in the plane, and $\kappa_{V}$ is set to 1. The observed best fit point is $f_{CP} = 0.00$.

Observed and expected likelihood ratio test statistic as a function of $\cos(\alpha)$. The overall modifier of the top-Higgs coupling strength is profiled such that the likelihood attains its minimum at each point in the plane, and $\kappa_{V}$ is set to 1. The observed best fit point is $\cos(\alpha) = +1.00$.

Observed and expected likelihood ratio test statistic as a function of $\kappa_{t}$ and $\tilde{\kappa}_{t}$ with $\kappa_{V}=1$, in the $H \rightarrow bb$, $H\rightarrow\gamma\gamma$, $H\rightarrow ZZ$, $H\rightarrow WW$, and $H\rightarrow\tau\tau$ decay channels and in their combination.

Observed and expected likelihood ratio test statistic as a function of $f_{CP}$, in the $H \rightarrow bb$, $H\rightarrow\gamma\gamma$, $H\rightarrow ZZ$, $H\rightarrow WW$, and $H\rightarrow\tau\tau$ decay channel and in their combination. The overall modifier of the top-Higgs coupling strength is profiled such that the likelihood attains its minimum at each point in the plane, and $\kappa_{V}$ is set to 1.

Observed and expected likelihood ratio test statistic as a function of $\cos(\alpha)$, in the $H \rightarrow bb$, $H\rightarrow\gamma\gamma$, $H\rightarrow ZZ$, $H\rightarrow WW$, and $H\rightarrow\tau\tau$ decay channels and in their combination. The overall modifier of the top-Higgs coupling strength is profiled such that the likelihood attains its minimum at each point in the plane, and $\kappa_{V}$ is set to 1.

Example description

Example description

Example description

Observed profiled likelihood on $f_{\Lambda1}$ using Approach 1. The signal strength modifiers are treated as free parameters. Axis scales are varied to improve the visibility of important features.

Expected profiled likelihood on $f_{a3}$ using Approach 1. The signal strength modifiers are treated as free parameters.

Observed profiled likelihood on $f_{a3}$ using Approach 1. The signal strength modifiers are treated as free parameters.

Expected profiled likelihood on $f_{\Lambda1}^{Z\gamma}$ using Approach 1. The signal strength modifiers are treated as free parameters.

Observed profiled likelihood on $f_{\Lambda1}^{Z\gamma}$ using Approach 1. The signal strength modifiers are treated as free parameters.

Expected profiled likelihood on $f_{a2}$ using Approach 2. The other two anomalous coupling cross section fractions are fixed to zero. The signal strength modifiers are treated as free parameters. Axis scales are varied to improve the visibility of important features.

Observed profiled likelihood on $f_{a2}$ using Approach 2. The other two anomalous coupling cross section fractions are fixed to zero. The signal strength modifiers are treated as free parameters. Axis scales are varied to improve the visibility of important features.

Expected profiled likelihood on $f_{\Lambda1}$ using Approach 2. The other two anomalous coupling cross section fractions are fixed to zero. The signal strength modifiers are treated as free parameters. Axis scales are varied to improve the visibility of important features.

Observed profiled likelihood on $f_{\Lambda1}$ using Approach 2. The other two anomalous coupling cross section fractions are fixed to zero. The signal strength modifiers are treated as free parameters. Axis scales are varied to improve the visibility of important features.

Expected profiled likelihood on $f_{a3}$ using Approach 2. The other two anomalous coupling cross section fractions are fixed to zero. The signal strength modifiers are treated as free parameters.

Observed profiled likelihood on $f_{a3}$ using Approach 2. The other two anomalous coupling cross section fractions are fixed to zero. The signal strength modifiers are treated as free parameters.

Expected profiled likelihood on $f_{a2}$ using Approach 2. The other two anomalous coupling cross section fractions are left floating in the fit. The signal strength modifiers are treated as free parameters. Axis scales are varied to improve the visibility of important features.

Observed profiled likelihood on $f_{a2}$ using Approach 2. The other two anomalous coupling cross section fractions are left floating in the fit. The signal strength modifiers are treated as free parameters. Axis scales are varied to improve the visibility of important features.

Expected profiled likelihood on $f_{\Lambda1}$ using Approach 2. The other two anomalous coupling cross section fractions are left floating in the fit. The signal strength modifiers are treated as free parameters. Axis scales are varied to improve the visibility of important features.

Observed profiled likelihood on $f_{\Lambda1}$ using Approach 2. The other two anomalous coupling cross section fractions are left floating in the fit. The signal strength modifiers are treated as free parameters. Axis scales are varied to improve the visibility of important features.

Expected profiled likelihood on $f_{a3}$ using Approach 2. The other two anomalous coupling cross section fractions are left floating in the fit. The signal strength modifiers are treated as free parameters.

Observed profiled likelihood on $f_{a3}$ using Approach 2. The other two anomalous coupling cross section fractions are left floating in the fit. The signal strength modifiers are treated as free parameters.

Expected profiled likelihood on the $\delta c_{Z}$ coupling of the SMEFT Higgs basis.

Observed profiled likelihood on the $\delta c_{Z}$ coupling of the SMEFT Higgs basis.

Expected profiled likelihood on the $c_{Z\Box}$ coupling of the SMEFT Higgs basis.

Observed profiled likelihood on the $c_{Z\Box}$ coupling of the SMEFT Higgs basis.

Expected profiled likelihood on the $c_{ZZ}$ coupling of the SMEFT Higgs basis.

Observed profiled likelihood on the $c_{ZZ}$ coupling of the SMEFT Higgs basis.

Expected profiled likelihood on the $\tilde{c}_{ZZ}$ coupling of the SMEFT Higgs basis.

Observed profiled likelihood on the $\tilde{c}_{ZZ}$ coupling of the SMEFT Higgs basis.

Expected profiled likelihood on $f_{a3}^{ggH}$. The signal strength modifiers and the CP-odd HVV anomalous coupling cross section fraction are treated as free parameters.

Observed profiled likelihood on $f_{a3}^{ggH}$. The signal strength modifiers and the CP-odd HVV anomalous coupling cross section fraction are treated as free parameters.

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 \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 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 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 $\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 $\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 $\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 $\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 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 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).