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.

Observed event yields in $\textrm{SR}$ compared with post-fit expectations from Monte Carlo simulations, as a function of the scalar sum of lepton and jet transverse momenta, $H_{\mathrm{T}}$. The last bin includes overflow events. `Signal (prod.)' and `Signal (dec.)' refer to the single-top-quark production and top-quark pair decay signal contributions, respectively.

Expected and observed 95$\%$ CL upper limits on the branching ratio corresponding to the decay of a top quark to a muon and a tau lepton through a cLFV process for scalar, vector and tensor couplings. In the vector case, all vector operators are assumed to contribute simultaneously with the same effective coupling strength. The table shows the endpoint values of $\mathcal{B}(\mathrm{t}\to\mu\tau\mathrm{u})$ and $\mathcal{B}(\mathrm{t}\to\mu\tau\mathrm{c})$, which are interpolated linearly for each limit.

Expected and observed 95$\%$ CL upper limits on the Wilson coefficients for scalar, vector and tensor couplings of a top quark to a muon and a tau lepton through a cLFV process. In the vector case, all vector operators are assumed to contribute simultaneously with the same effective coupling strength. The table shows the endpoint values of the $|c^{\mu\tau\mathrm{ut}}|$ and $|c^{\mu\tau\mathrm{ct}}|$ Wilson coefficients, which are interpolated assuming a linear relationship between $\mathcal{B}(\mathrm{t}\to\mu\tau\mathrm{u})$ and $\mathcal{B}(\mathrm{t}\to\mu\tau\mathrm{c})$.

Observed event yields in $\textrm{SR}$ compared with pre-fit expectations from Monte Carlo simulations, as a function of the scalar sum of lepton and jet transverse momenta, $H_{\mathrm{T}}$. The last bin includes overflow events. The signal yields represent a leptoquark mass of $m_{S_{1}} = 1$ TeV and a coupling strength of $\lambda^{\textrm{LQ}} = 2.0$.

Observed event yields in $\textrm{SR}$ compared with post-fit expectations from Monte Carlo simulations, as a function of the scalar sum of lepton and jet transverse momenta, $H_{\mathrm{T}}$. The last bin includes overflow events. The signal yields represent a leptoquark mass of $m_{S_{1}} = 1$ TeV and a coupling strength of $\lambda^{\textrm{LQ}} = 2.0$.

Observed and expected exclusion 95$\%$ CL upper limits on the $S_{1}$ leptoquark coupling strength $\lambda^{\textrm{LQ}}$ as a function of LQ mass, $m_{S_{1}}$.

Theoretical cross-sections for single-top-quark production and top-quark decays through cLFV interactions involving vector, scalar and tensor EFT Wilson coefficients. The column titled as $c^{(ijk3)}_{\textrm{vector}}$ represents the individual cross-section contributions from each of $c_{\mathsf{lq}}^{-(ijk3)}$, $c_{\mathsf{eq}}^{(ijk3)}$, $c_{\mathsf{lu}}^{(ijk3)}$ and $c_{\mathsf{eu}}^{(ijk3)}$. The coefficient indices represent the lepton flavour generations ($i,j = 1,2,3$ where $i \neq j$) and light quark flavour generations ($k = 1,2$). The single-top-quark production cross-sections are quoted for $u$- and $c$-quark couplings separately, while they are combined for the $t\bar{t}$ decay process ($q_k = u,c$). The scale and PDF uncertainties are given. The value of each Wilson coefficient is set to 1.0 for the calculation of the cross-section.

Requirements for each analysis region. The symbol $\ell_{3}$ denotes the lowest $p_{\mathrm{T}}$ lepton. In $\mathrm{CR}t\bar{t}\mu$ an additional requirement is placed on the leading muon $p_{\mathrm{T}}$ ($p^{\mu_1}_\mathrm{T}$) and dilepton invariant masses in order to reject signal events.

Post-fit event yields for each analysis region entering the fit, with correlations on the full systematic uncertainties taken into account as determined in the fit under a signal+background hypothesis. The 'fake electron' category in the CR$t\bar{t}\mu$ control region collects small contributions primarily from $t\bar{t}\gamma$ and $Z\gamma$ which enter the event selection due to photon conversions.

Expected and observed 95$\%$ CL upper limits on Wilson coefficients corresponding to 2Q2L EFT operators which could introduce a cLFV top decay in the $\mu\tau$ channel, and existing limits from JHEP04 (2019) 014 (previous). The previous limits shown for $c^{1(ijk3)}_{\mathsf{lequ}}$ and $c^{3(ijk3)}_{\mathsf{lequ}}$ are tightened by a factor of $\sqrt{2}$, as JHEP04 (2019) 014 does not assume that these operators are Hermitian. Results are shown separately for the $\mu\tau ut$ and $\mu\tau ct$ interactions. The lepton generations are denoted by $i,j=2,3$ for $\mu$ and $\tau$ (where $i \neq j$) and the quark generations are denoted by $k=1,2$ for $u$ and $c$, respectively.

Expected and observed 95$\%$ CL upper limits on the branching ratio ($\mathcal{B}$) corresponding to the decay of a top quark to a muon and a $\tau$ lepton through a cLFV process using specific Wilson coefficients corresponding to 2Q2L EFT operators. Results are shown separately for $\mu\tau ut$ and $\mu\tau ct$ interactions. The lepton generations are denoted by $i,j=2,3$ for $\mu$ and $\tau$ (where $i \neq j$) and the quark generations are denoted by $k=1,2$ for $u$ and $c$, respectively.

Expected and observed 95% CL upper limits on the inclusive branching ratio ($\mathcal{B}$) corresponding to the decay of a top quark to a muon and a $\tau$-lepton through a cLFV process. Limits are shown for the statistical uncertainty only and for the full set of statistical and systematic uncertainties.

Expected and observed 95$\%$ CL upper limits on the Wilson coefficients corresponding to EFT operators inducing the decay of a top quark to a muon, a tau lepton and an up quark through a cLFV process ($\mu\tau ut$ interaction).

Expected and observed 95$\%$ CL upper limits on the Wilson coefficients corresponding to EFT operators inducing the decay of a top quark to a muon, a tau lepton and a charm quark through a cLFV process ($\mu\tau ct$ interaction).

Ranking of nuisance parameters by post-fit impact on cLFV signal strength $\mu$ (the parameter of interest, POI). The pre-fit (post-fit) impact of each nuisance parameter is calculated as the difference to the fitted value of cLFV signal strength between the nominal fit and a fit when fixing the corresponding nuisance parameter to $\hat{\theta}\pm\Delta\theta$ ($\hat{\theta}\pm\Delta\hat{\theta}$), where $\hat{\theta}$ is the best-fit value of the nuisance parameter and $\Delta\theta$ ($\Delta\hat{\theta}$) is its pre-fit (post-fit) uncertainty. The pulls of the nuisance parameters are calculated relative to zero. These pulls and their post-fit errors, $\Delta\hat{\theta} / \Delta\theta$, are given in the second column of the table.