Expected (dashed black line) and observed (solid red line) 95% CL exclusion limits on the higgsino simplified model being considered. These are shown with $\pm 1\sigma_{\mathrm{exp}}$ (yellow band) from experimental systematic and statistical uncertainties, and with $\pm 1\sigma_{\mathrm{theory}}^{\mathrm{SUSY}}$ (red dotted lines) from signal cross-section uncertainties, respectively. The limits set by the latest ATLAS searches using the soft lepton and disappearing track signatures are illustrated by the blue and green regions, respectively, while the limit imposed by the LEP experiments is shown in gray. The dot-dashed gray line indicates the predicted mass-splitting for the pure higgsino scenario.

Expected (dashed black line) and observed (solid red line) 95% CL exclusion limits on the higgsino simplified model being considered. These are shown with $\pm 1\sigma_{\mathrm{exp}}$ (yellow band) from experimental systematic and statistical uncertainties, and with $\pm 1\sigma_{\mathrm{theory}}^{\mathrm{SUSY}}$ (red dotted lines) from signal cross-section uncertainties, respectively. The limits set by the latest ATLAS searches using the soft lepton and disappearing track signatures are illustrated by the blue and green regions, respectively, while the limit imposed by the LEP experiments is shown in gray. The dot-dashed gray line indicates the predicted mass-splitting for the pure higgsino scenario.

Expected (dashed black line) and observed (solid red line) 95% CL exclusion limits on the higgsino simplified model being considered. These are shown with $\pm 1\sigma_{\mathrm{exp}}$ (yellow band) from experimental systematic and statistical uncertainties, and with $\pm 1\sigma_{\mathrm{theory}}^{\mathrm{SUSY}}$ (red dotted lines) from signal cross-section uncertainties, respectively. The limits set by the latest ATLAS searches using the soft lepton and disappearing track signatures are illustrated by the blue and green regions, respectively, while the limit imposed by the LEP experiments is shown in gray. The dot-dashed gray line indicates the predicted mass-splitting for the pure higgsino scenario.

Expected (dashed black line) and observed (solid red line) 95% CL exclusion limits on the higgsino simplified model being considered. These are shown with $\pm 1\sigma_{\mathrm{exp}}$ (yellow band) from experimental systematic and statistical uncertainties, and with $\pm 1\sigma_{\mathrm{theory}}^{\mathrm{SUSY}}$ (red dotted lines) from signal cross-section uncertainties, respectively. The limits set by the latest ATLAS searches using the soft lepton and disappearing track signatures are illustrated by the blue and green regions, respectively, while the limit imposed by the LEP experiments is shown in gray. The dot-dashed gray line indicates the predicted mass-splitting for the pure higgsino scenario.

Expected and observed CLs values per signal point represented by the grey numbers. The expected (dashed) and observed (solid) 95% CL exclusion limits are overlaid along with $\pm 1\sigma_{\mathrm{exp}}$ (yellow band) from experimental systematic and statistical uncertainties, and with $\pm 1\sigma_{\mathrm{theory}}^{\mathrm{SUSY}}$ (red dotted lines) from signal cross-section uncertainties, respectively.

Expected and observed CLs values per signal point represented by the grey numbers. The expected (dashed) and observed (solid) 95% CL exclusion limits are overlaid along with $\pm 1\sigma_{\mathrm{exp}}$ (yellow band) from experimental systematic and statistical uncertainties, and with $\pm 1\sigma_{\mathrm{theory}}^{\mathrm{SUSY}}$ (red dotted lines) from signal cross-section uncertainties, respectively.

Expected and observed cross-section upper-limit per signal point represented by the grey numbers. The expected (dashed) and observed (solid) 95% CL exclusion limits are overlaid along with $\pm 1\sigma_{\mathrm{exp}}$ (yellow band) from experimental systematic and statistical uncertainties, and with $\pm 1\sigma_{\mathrm{theory}}^{\mathrm{SUSY}}$ (red dotted lines) from signal cross-section uncertainties, respectively.

Expected and observed cross-section upper-limit per signal point represented by the grey numbers. The expected (dashed) and observed (solid) 95% CL exclusion limits are overlaid along with $\pm 1\sigma_{\mathrm{exp}}$ (yellow band) from experimental systematic and statistical uncertainties, and with $\pm 1\sigma_{\mathrm{theory}}^{\mathrm{SUSY}}$ (red dotted lines) from signal cross-section uncertainties, respectively.

Truth-level signal acceptances for each production process ($\tilde{\chi}_1^\pm \tilde{\chi}_1^0$, $\tilde{\chi}_1^\pm \tilde{\chi}_2^0$, $\tilde{\chi}_1^+ \tilde{\chi}_1^-$, and $\tilde{\chi}_2^0 \tilde{\chi}_1^0$) in a SR with the $S(d_0)$ requirement removed. The acceptance is defined as the fraction of accepted events divided by the total number of events in the generator-level signal Monte Carlo simulation, where the signal candidate track is identified as the charged particle with the largest distance between the interaction vertex and the secondary vertex of the higgsino decays.

Truth-level signal acceptances for each production process ($\tilde{\chi}_1^\pm \tilde{\chi}_1^0$, $\tilde{\chi}_1^\pm \tilde{\chi}_2^0$, $\tilde{\chi}_1^+ \tilde{\chi}_1^-$, and $\tilde{\chi}_2^0 \tilde{\chi}_1^0$) in a SR with the $S(d_0)$ requirement removed. The acceptance is defined as the fraction of accepted events divided by the total number of events in the generator-level signal Monte Carlo simulation, where the signal candidate track is identified as the charged particle with the largest distance between the interaction vertex and the secondary vertex of the higgsino decays.

Truth-level signal acceptances for each production process ($\tilde{\chi}_1^\pm \tilde{\chi}_1^0$, $\tilde{\chi}_1^\pm \tilde{\chi}_2^0$, $\tilde{\chi}_1^+ \tilde{\chi}_1^-$, and $\tilde{\chi}_2^0 \tilde{\chi}_1^0$) in a SR with the $S(d_0)$ requirement removed. The acceptance is defined as the fraction of accepted events divided by the total number of events in the generator-level signal Monte Carlo simulation, where the signal candidate track is identified as the charged particle with the largest distance between the interaction vertex and the secondary vertex of the higgsino decays.

Truth-level signal acceptances for each production process ($\tilde{\chi}_1^\pm \tilde{\chi}_1^0$, $\tilde{\chi}_1^\pm \tilde{\chi}_2^0$, $\tilde{\chi}_1^+ \tilde{\chi}_1^-$, and $\tilde{\chi}_2^0 \tilde{\chi}_1^0$) in a SR with the $S(d_0)$ requirement removed. The acceptance is defined as the fraction of accepted events divided by the total number of events in the generator-level signal Monte Carlo simulation, where the signal candidate track is identified as the charged particle with the largest distance between the interaction vertex and the secondary vertex of the higgsino decays.

Truth-level signal acceptances for each production process ($\tilde{\chi}_1^\pm \tilde{\chi}_1^0$, $\tilde{\chi}_1^\pm \tilde{\chi}_2^0$, $\tilde{\chi}_1^+ \tilde{\chi}_1^-$, and $\tilde{\chi}_2^0 \tilde{\chi}_1^0$) in a SR with the $S(d_0)$ requirement removed. The acceptance is defined as the fraction of accepted events divided by the total number of events in the generator-level signal Monte Carlo simulation, where the signal candidate track is identified as the charged particle with the largest distance between the interaction vertex and the secondary vertex of the higgsino decays.

Truth-level signal acceptances for each production process ($\tilde{\chi}_1^\pm \tilde{\chi}_1^0$, $\tilde{\chi}_1^\pm \tilde{\chi}_2^0$, $\tilde{\chi}_1^+ \tilde{\chi}_1^-$, and $\tilde{\chi}_2^0 \tilde{\chi}_1^0$) in a SR with the $S(d_0)$ requirement removed. The acceptance is defined as the fraction of accepted events divided by the total number of events in the generator-level signal Monte Carlo simulation, where the signal candidate track is identified as the charged particle with the largest distance between the interaction vertex and the secondary vertex of the higgsino decays.

Signal efficiencies in SR-Low for each production process ($\tilde{\chi}_1^\pm \tilde{\chi}_1^0$, $\tilde{\chi}_1^\pm \tilde{\chi}_2^0$, $\tilde{\chi}_1^+ \tilde{\chi}_1^-$, and $\tilde{\chi}_2^0 \tilde{\chi}_1^0$), defined by the number of events of reconstructed-level signal simulation divided by the number of events obtained at generator level, where the $S(d_0)$ selecton efficiency has the largest impact. The higgsino decay products from $\Delta \mathrm{m}(\tilde{\chi}_1^\pm,\tilde{\chi}_1^0) < 0.4$ GeV signal have $p_{\mathrm{T}}$ too low to be reconstructed as the signal candidate tracks, and therefore the identified signal candidate tracks are typically from pile-up collisions or underlying events similar to the QCD track background, causing a low $S(d_0)$ selection efficiency in these plots.

Signal efficiencies in SR-Low for each production process ($\tilde{\chi}_1^\pm \tilde{\chi}_1^0$, $\tilde{\chi}_1^\pm \tilde{\chi}_2^0$, $\tilde{\chi}_1^+ \tilde{\chi}_1^-$, and $\tilde{\chi}_2^0 \tilde{\chi}_1^0$), defined by the number of events of reconstructed-level signal simulation divided by the number of events obtained at generator level, where the $S(d_0)$ selecton efficiency has the largest impact. The higgsino decay products from $\Delta \mathrm{m}(\tilde{\chi}_1^\pm,\tilde{\chi}_1^0) < 0.4$ GeV signal have $p_{\mathrm{T}}$ too low to be reconstructed as the signal candidate tracks, and therefore the identified signal candidate tracks are typically from pile-up collisions or underlying events similar to the QCD track background, causing a low $S(d_0)$ selection efficiency in these plots.

Signal efficiencies in SR-Low for each production process ($\tilde{\chi}_1^\pm \tilde{\chi}_1^0$, $\tilde{\chi}_1^\pm \tilde{\chi}_2^0$, $\tilde{\chi}_1^+ \tilde{\chi}_1^-$, and $\tilde{\chi}_2^0 \tilde{\chi}_1^0$), defined by the number of events of reconstructed-level signal simulation divided by the number of events obtained at generator level, where the $S(d_0)$ selecton efficiency has the largest impact. The higgsino decay products from $\Delta \mathrm{m}(\tilde{\chi}_1^\pm,\tilde{\chi}_1^0) < 0.4$ GeV signal have $p_{\mathrm{T}}$ too low to be reconstructed as the signal candidate tracks, and therefore the identified signal candidate tracks are typically from pile-up collisions or underlying events similar to the QCD track background, causing a low $S(d_0)$ selection efficiency in these plots.

Signal efficiencies in SR-Low for each production process ($\tilde{\chi}_1^\pm \tilde{\chi}_1^0$, $\tilde{\chi}_1^\pm \tilde{\chi}_2^0$, $\tilde{\chi}_1^+ \tilde{\chi}_1^-$, and $\tilde{\chi}_2^0 \tilde{\chi}_1^0$), defined by the number of events of reconstructed-level signal simulation divided by the number of events obtained at generator level, where the $S(d_0)$ selecton efficiency has the largest impact. The higgsino decay products from $\Delta \mathrm{m}(\tilde{\chi}_1^\pm,\tilde{\chi}_1^0) < 0.4$ GeV signal have $p_{\mathrm{T}}$ too low to be reconstructed as the signal candidate tracks, and therefore the identified signal candidate tracks are typically from pile-up collisions or underlying events similar to the QCD track background, causing a low $S(d_0)$ selection efficiency in these plots.

Signal efficiencies in SR-Low for each production process ($\tilde{\chi}_1^\pm \tilde{\chi}_1^0$, $\tilde{\chi}_1^\pm \tilde{\chi}_2^0$, $\tilde{\chi}_1^+ \tilde{\chi}_1^-$, and $\tilde{\chi}_2^0 \tilde{\chi}_1^0$), defined by the number of events of reconstructed-level signal simulation divided by the number of events obtained at generator level, where the $S(d_0)$ selecton efficiency has the largest impact. The higgsino decay products from $\Delta \mathrm{m}(\tilde{\chi}_1^\pm,\tilde{\chi}_1^0) < 0.4$ GeV signal have $p_{\mathrm{T}}$ too low to be reconstructed as the signal candidate tracks, and therefore the identified signal candidate tracks are typically from pile-up collisions or underlying events similar to the QCD track background, causing a low $S(d_0)$ selection efficiency in these plots.

Signal efficiencies in SR-Low for each production process ($\tilde{\chi}_1^\pm \tilde{\chi}_1^0$, $\tilde{\chi}_1^\pm \tilde{\chi}_2^0$, $\tilde{\chi}_1^+ \tilde{\chi}_1^-$, and $\tilde{\chi}_2^0 \tilde{\chi}_1^0$), defined by the number of events of reconstructed-level signal simulation divided by the number of events obtained at generator level, where the $S(d_0)$ selecton efficiency has the largest impact. The higgsino decay products from $\Delta \mathrm{m}(\tilde{\chi}_1^\pm,\tilde{\chi}_1^0) < 0.4$ GeV signal have $p_{\mathrm{T}}$ too low to be reconstructed as the signal candidate tracks, and therefore the identified signal candidate tracks are typically from pile-up collisions or underlying events similar to the QCD track background, causing a low $S(d_0)$ selection efficiency in these plots.

Signal efficiencies in SR-High for each production process ($\tilde{\chi}_1^\pm \tilde{\chi}_1^0$, $\tilde{\chi}_1^\pm \tilde{\chi}_2^0$, $\tilde{\chi}_1^+ \tilde{\chi}_1^-$, and $\tilde{\chi}_2^0 \tilde{\chi}_1^0$), defined by the number of events of reconstructed-level signal simulation divided by the number of events obtained at generator level, where the $S(d_0)$ selecton efficiency has the largest impact. The higgsino decay products from $\Delta \mathrm{m}(\tilde{\chi}_1^\pm,\tilde{\chi}_1^0) < 0.4$ GeV signal have $p_{\mathrm{T}}$ too low to be reconstructed as the signal candidate tracks, and therefore the identified signal candidate tracks are typically from pile-up collisions or underlying events similar to the QCD track background, causing a low $S(d_0)$ selection efficiency in these plots.

Signal efficiencies in SR-High for each production process ($\tilde{\chi}_1^\pm \tilde{\chi}_1^0$, $\tilde{\chi}_1^\pm \tilde{\chi}_2^0$, $\tilde{\chi}_1^+ \tilde{\chi}_1^-$, and $\tilde{\chi}_2^0 \tilde{\chi}_1^0$), defined by the number of events of reconstructed-level signal simulation divided by the number of events obtained at generator level, where the $S(d_0)$ selecton efficiency has the largest impact. The higgsino decay products from $\Delta \mathrm{m}(\tilde{\chi}_1^\pm,\tilde{\chi}_1^0) < 0.4$ GeV signal have $p_{\mathrm{T}}$ too low to be reconstructed as the signal candidate tracks, and therefore the identified signal candidate tracks are typically from pile-up collisions or underlying events similar to the QCD track background, causing a low $S(d_0)$ selection efficiency in these plots.

Signal efficiencies in SR-High for each production process ($\tilde{\chi}_1^\pm \tilde{\chi}_1^0$, $\tilde{\chi}_1^\pm \tilde{\chi}_2^0$, $\tilde{\chi}_1^+ \tilde{\chi}_1^-$, and $\tilde{\chi}_2^0 \tilde{\chi}_1^0$), defined by the number of events of reconstructed-level signal simulation divided by the number of events obtained at generator level, where the $S(d_0)$ selecton efficiency has the largest impact. The higgsino decay products from $\Delta \mathrm{m}(\tilde{\chi}_1^\pm,\tilde{\chi}_1^0) < 0.4$ GeV signal have $p_{\mathrm{T}}$ too low to be reconstructed as the signal candidate tracks, and therefore the identified signal candidate tracks are typically from pile-up collisions or underlying events similar to the QCD track background, causing a low $S(d_0)$ selection efficiency in these plots.

Signal efficiencies in SR-High for each production process ($\tilde{\chi}_1^\pm \tilde{\chi}_1^0$, $\tilde{\chi}_1^\pm \tilde{\chi}_2^0$, $\tilde{\chi}_1^+ \tilde{\chi}_1^-$, and $\tilde{\chi}_2^0 \tilde{\chi}_1^0$), defined by the number of events of reconstructed-level signal simulation divided by the number of events obtained at generator level, where the $S(d_0)$ selecton efficiency has the largest impact. The higgsino decay products from $\Delta \mathrm{m}(\tilde{\chi}_1^\pm,\tilde{\chi}_1^0) < 0.4$ GeV signal have $p_{\mathrm{T}}$ too low to be reconstructed as the signal candidate tracks, and therefore the identified signal candidate tracks are typically from pile-up collisions or underlying events similar to the QCD track background, causing a low $S(d_0)$ selection efficiency in these plots.

Signal efficiencies in SR-High for each production process ($\tilde{\chi}_1^\pm \tilde{\chi}_1^0$, $\tilde{\chi}_1^\pm \tilde{\chi}_2^0$, $\tilde{\chi}_1^+ \tilde{\chi}_1^-$, and $\tilde{\chi}_2^0 \tilde{\chi}_1^0$), defined by the number of events of reconstructed-level signal simulation divided by the number of events obtained at generator level, where the $S(d_0)$ selecton efficiency has the largest impact. The higgsino decay products from $\Delta \mathrm{m}(\tilde{\chi}_1^\pm,\tilde{\chi}_1^0) < 0.4$ GeV signal have $p_{\mathrm{T}}$ too low to be reconstructed as the signal candidate tracks, and therefore the identified signal candidate tracks are typically from pile-up collisions or underlying events similar to the QCD track background, causing a low $S(d_0)$ selection efficiency in these plots.

Signal efficiencies in SR-High for each production process ($\tilde{\chi}_1^\pm \tilde{\chi}_1^0$, $\tilde{\chi}_1^\pm \tilde{\chi}_2^0$, $\tilde{\chi}_1^+ \tilde{\chi}_1^-$, and $\tilde{\chi}_2^0 \tilde{\chi}_1^0$), defined by the number of events of reconstructed-level signal simulation divided by the number of events obtained at generator level, where the $S(d_0)$ selecton efficiency has the largest impact. The higgsino decay products from $\Delta \mathrm{m}(\tilde{\chi}_1^\pm,\tilde{\chi}_1^0) < 0.4$ GeV signal have $p_{\mathrm{T}}$ too low to be reconstructed as the signal candidate tracks, and therefore the identified signal candidate tracks are typically from pile-up collisions or underlying events similar to the QCD track background, causing a low $S(d_0)$ selection efficiency in these plots.

Event selection cutflows for signal samples with $m(\tilde{\chi}_{1}^0)$ = 150 GeV and $\Delta m(\tilde{\chi}_{1}^\pm, \tilde{\chi}_{1}^0)$ = 1.5, 1.0, and 0.75 GeV, including all six production processes ($\tilde{\chi}_1^\pm \tilde{\chi}_1^0$, $\tilde{\chi}_1^\pm \tilde{\chi}_2^0$, $\tilde{\chi}_1^+ \tilde{\chi}_1^-$, and $\tilde{\chi}_2^0 \tilde{\chi}_1^0$). The cross-section used to obtain the initial number of events ($\sigma(\mathrm{n}_{\mathrm{jets}}) \geq 1$) refers to an emission of at least one gluon or quark with $p_{\mathrm{T}} > 50$ GeV at the parton level.

Event selection cutflows for signal samples with $m(\tilde{\chi}_{1}^0)$ = 150 GeV and $\Delta m(\tilde{\chi}_{1}^\pm, \tilde{\chi}_{1}^0)$ = 0.5, 0.35, and 0.25 GeV, including all six production processes ($\tilde{\chi}_1^\pm \tilde{\chi}_1^0$, $\tilde{\chi}_1^\pm \tilde{\chi}_2^0$, $\tilde{\chi}_1^+ \tilde{\chi}_1^-$, and $\tilde{\chi}_2^0 \tilde{\chi}_1^0$). The cross-section used to obtain the initial number of events ($\sigma(\mathrm{n}_{\mathrm{jets}}) \geq 1$) refers to an emission of at least one gluon or quark with $p_{\mathrm{T}} > 50$ GeV at the parton level.

The double ratio of particle yield of f0((980) to K*(892)0

Transverse-momentum-differential particle yield ratio of f$_{0}$(980) to charged pions and its double ratio

Transverse-momentum-differential particle yield ratio of f$_{0}$(980) to K$^{*}(892)$^{0}$ and its double ratio

Nuclear modification factor of f0(980) in p-Pb collisions at 5.02 TeV

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 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.

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.

Summary of observed and expected background yields in the SR after the likelihood fit under the background-only hypothesis.

Intensive skewness of $\langle p_\mathrm{T}\rangle$ as a function of $\langle\mathrm{d}N_\mathrm{ch}/\mathrm{d}\eta\rangle^{1/3}_{|\eta|<0.5}$ in pp collisions at $\sqrt{s}$ = 5.02 TeV.

Intensive skewness of $\langle p_\mathrm{T}\rangle$ as a function of $\langle\mathrm{d}N_\mathrm{ch}/\mathrm{d}\eta\rangle^{1/3}_{|\eta|<0.5}$ in Xe$-$Xe collisions at $\sqrt{s_\mathrm{NN}}$ = 5.02 TeV.

Intensive skewness of $\langle p_\mathrm{T}\rangle$ as a function of $\langle\mathrm{d}N_\mathrm{ch}/\mathrm{d}\eta\rangle^{1/3}_{|\eta|<0.5}$ in Pb$-$Pb collisions at $\sqrt{s_\mathrm{NN}}$ = 5.02 TeV.

Kurtosis of $\langle p_\mathrm{T}\rangle$ as a function of $\langle\mathrm{d}N_\mathrm{ch}/\mathrm{d}\eta\rangle^{1/3}_{|\eta|<0.5}$ in pp collisions at $\sqrt{s}$ = 5.02 TeV.

Kurtosis of $\langle p_\mathrm{T}\rangle$ as a function of $\langle\mathrm{d}N_\mathrm{ch}/\mathrm{d}\eta\rangle^{1/3}_{|\eta|<0.5}$ in Xe$-$Xe collisions at $\sqrt{s_\mathrm{NN}}$ = 5.02 TeV.

Kurtosis of $\langle p_\mathrm{T}\rangle$ as a function of $\langle\mathrm{d}N_\mathrm{ch}/\mathrm{d}\eta\rangle^{1/3}_{|\eta|<0.5}$ in Pb$-$Pb collisions at $\sqrt{s_\mathrm{NN}}$ = 5.02 TeV.

Observed and expected distributions of the invariant mass of the four-lepton system in the $A\to ZH\to 4\ell+X$ search for SR1 under a background-only fit to data. The total background includes the $q\overline{q}\to ZZ$, $gg\to ZZ$, $q\overline{q}\to ZZ$ (EW), $VVV$, $t\overline{t}V$, $t\overline{t}$, $Z$+jets and $WZ$ processes. The distribution of the $(m_{A}, m_{H}) = (510, 380)$ GeV signal is normalised to the observed upper limit on the cross-section (29.9 fb).

Observed and expected distributions of the invariant mass of the four-lepton system in the $A\to ZH\to 4\ell+X$ search for SR2 under a background-only fit to data. The total background includes the $q\overline{q}\to ZZ$, $gg\to ZZ$, $q\overline{q}\to ZZ$ (EW), $VVV$, $t\overline{t}V$, $t\overline{t}$, $Z$+jets and $WZ$ processes. The distribution of the $(m_{A}, m_{H}) = (510, 380)$ GeV signal is normalised to the observed upper limit on the cross-section (29.9 fb).

Observed and expected distributions of the invariant mass of the four-lepton system in the $A\to ZH\to 4\ell+X$ search for SR3 under a background-only fit to data. The total background includes the $q\overline{q}\to ZZ$, $gg\to ZZ$, $q\overline{q}\to ZZ$ (EW), $VVV$, $t\overline{t}V$, $t\overline{t}$, $Z$+jets and $WZ$ processes. The distribution of the $(m_{A}, m_{H}) = (510, 380)$ GeV signal is normalised to the observed upper limit on the cross-section (29.9 fb).

Observed and expected distributions of the invariant mass of the four-lepton system in the $A\to ZH\to 4\ell+X$ search for SR4 under a background-only fit to data. The total background includes the $q\overline{q}\to ZZ$, $gg\to ZZ$, $q\overline{q}\to ZZ$ (EW), $VVV$, $t\overline{t}V$, $t\overline{t}$, $Z$+jets and $WZ$ processes. The distribution of the $(m_{A}, m_{H}) = (510, 380)$ GeV signal is normalised to the observed upper limit on the cross-section (29.9 fb).

Observed and expected distributions of the invariant mass of the four-lepton system in the $A\to ZH\to 4\ell+X$ search for SR5 under a background-only fit to data. The total background includes the $q\overline{q}\to ZZ$, $gg\to ZZ$, $q\overline{q}\to ZZ$ (EW), $VVV$, $t\overline{t}V$, $t\overline{t}$, $Z$+jets and $WZ$ processes. The distribution of the $(m_{A}, m_{H}) = (510, 380)$ GeV signal is normalised to the observed upper limit on the cross-section (29.9 fb).

Observed and expected distributions of the invariant mass of the four-lepton system in the $A\to ZH\to 4\ell+X$ search for SR6 under a background-only fit to data. The total background includes the $q\overline{q}\to ZZ$, $gg\to ZZ$, $q\overline{q}\to ZZ$ (EW), $VVV$, $t\overline{t}V$, $t\overline{t}$, $Z$+jets and $WZ$ processes. The distribution of the $(m_{A}, m_{H}) = (510, 380)$ GeV signal is normalised to the observed upper limit on the cross-section (29.9 fb).

Observed and expected distributions of the invariant mass of the four-lepton system in the $A\to ZH\to 4\ell+X$ search for SR7 under a background-only fit to data. The total background includes the $q\overline{q}\to ZZ$, $gg\to ZZ$, $q\overline{q}\to ZZ$ (EW), $VVV$, $t\overline{t}V$, $t\overline{t}$, $Z$+jets and $WZ$ processes. The distribution of the $(m_{A}, m_{H}) = (510, 380)$ GeV signal is normalised to the observed upper limit on the cross-section (29.9 fb).

Local p0-values in the $(m_{H}, m_{R})$ plane for the $R\to SH\to 4\ell+E^{\textrm{miss}}_{\textrm{T}}$ search with $m_{S} = 160$ GeV.

Local p0-values in the $(m_{H}, m_{A})$ plane for the $A\to ZH\to 4\ell+X$ search.

The observed upper limits at 95% confidence level on $\sigma(gg\to R)\times \mathcal{B}(R\to SH)\times (H\to ZZ)$ across the $(m_{H}, m_{R})$ plane with $m_{S} = 160$ GeV for the $R\to SH\to 4\ell+E^{\textrm{miss}}_{\textrm{T}}$ search.

The expected upper limits at 95% confidence level on $\sigma(gg\to R)\times \mathcal{B}(R\to SH)\times (H\to ZZ)$ across the $(m_{H}, m_{R})$ plane with $m_{S} = 160$ GeV for the $R\to SH\to 4\ell+E^{\textrm{miss}}_{\textrm{T}}$ search.

The observed upper limits at 95% confidence level on $\sigma(gg\to A)\times \mathcal{B}(A\to ZH)\times (H\to ZZ)$ across the $(m_{H}, m_{A})$ plane for the $A\to ZH\to 4\ell+X$ search.

The expected upper limits at 95% confidence level on $\sigma(gg\to A)\times \mathcal{B}(A\to ZH)\times (H\to ZZ)$ across the $(m_{H}, m_{A})$ plane for the $A\to ZH\to 4\ell+X$ search.

Non-prompt J/$\psi$ fraction as a function of average number of participants in Pb-Pb at 5.02 TeV

Non-prompt J/$\psi$ yields as a function of transverse momentum in Pb-Pb at 5.02 TeV, centrality 0-10%

Prompt J/$\psi$ yields as a function of transverse momentum in Pb-Pb at 5.02 TeV, centrality 0-10%

Prompt J/$\psi$ nuclear modification factor as a function of transverse momentum in Pb-Pb at 5.02 TeV, centrality 0-10%. Global uncertainty 2.3%

Prompt J/$\psi$ nuclear modification factor as a function of transverse momentum in Pb-Pb at 5.02 TeV, centrality 10-30%. Global uncertainty 2.3%

Prompt J/$\psi$ nuclear modification factor as a function of transverse momentum in Pb-Pb at 5.02 TeV, centrality 30-50%. Global uncertainty 2.7%

Prompt J/$\psi$ nuclear modification factor as a function of average number of participants in Pb-Pb at 5.02 TeV. Global uncertainty 2.2%

Non-prompt J/$\psi$ nuclear modification factor as a function of transverse momentum in Pb-Pb at 5.02 TeV, centrality 0-10%. Global uncertainty 2.3%

Non-prompt J/$\psi$ nuclear modification factor as a function of transverse momentum in Pb-Pb at 5.02 TeV,, centrality 10-30%. Global uncertainty 2.3%

Non-prompt J/$\psi$ nuclear modification factor as a function of transverse momentum in Pb-Pb at 5.02 TeV, centrality 30-50%. Global uncertainty 2.7%

Non-prompt J/$\psi$ nuclear modification factor as a function of average number of participants in Pb-Pb at 5.02 TeV. Global uncertainty 2.7%

Prompt J/$\psi$ yields as a function of transverse momentum in Pb-Pb at 5.02 TeV, centrality 30-50%

Non-prompt J/$\psi$ yields as a function of transverse momentum in Pb-Pb at 5.02 TeV, centrality 30-50%

$R_\mathrm{T}$ distribution using events with trigger particles $5<p_\mathrm{T}^\mathrm{trig}<40~\mathrm{GeV}/c$ in the pseudorapidity range of $|\eta|<0.8$ and with $p_\mathrm{T}>0.5~\mathrm{GeV}/c$ in pp collisions at $\sqrt{s}=13~\mathrm{TeV}$

$p_\mathrm{T}$ as a function of $R_\mathrm{T}$ in the toward region for different $R_\mathrm{T}$ intervals using events with trigger particles $8<p_\mathrm{T}^\mathrm{trig}<15~\mathrm{GeV}/c$ in the pseudorapidity range of $|\eta|<0.8$ and with $p_\mathrm{T}>0.5~\mathrm{GeV}/c$ for pp collisions at $\sqrt{s}=5.02~\mathrm{TeV}$

$p_\mathrm{T}$ as a function of $R_\mathrm{T}$ in the away region for different $R_\mathrm{T}$ intervals using events with trigger particles $8<p_\mathrm{T}^\mathrm{trig}<15~\mathrm{GeV}/c$ in the pseudorapidity range of $|\eta|<0.8$ and with $p_\mathrm{T}>0.5~\mathrm{GeV}/c$ for pp collisions at $\sqrt{s}=5.02~\mathrm{TeV}$

$p_\mathrm{T}$ as a function of $R_\mathrm{T}$ in the transverse region for different $R_\mathrm{T}$ intervals using events with trigger particles $8<p_\mathrm{T}^\mathrm{trig}<15~\mathrm{GeV}/c$ in the pseudorapidity range of $|\eta|<0.8$ and with $p_\mathrm{T}>0.5~\mathrm{GeV}/c$ for pp collisions at $\sqrt{s}=5.02~\mathrm{TeV}$

$p_\mathrm{T}$ as a function of $R_\mathrm{T}$ in the toward region for different $R_\mathrm{T}$ intervals using events with trigger particles $8<p_\mathrm{T}^\mathrm{trig}<15~\mathrm{GeV}/c$ in the pseudorapidity range of $|\eta|<0.8$ and with $p_\mathrm{T}>0.5~\mathrm{GeV}/c$ for p-Pb collisions at $\sqrt{s_\mathrm{NN}}=5.02~\mathrm{TeV}$

$p_\mathrm{T}$ as a function of $R_\mathrm{T}$ in the away region for different $R_\mathrm{T}$ intervals using events with trigger particles $8<p_\mathrm{T}^\mathrm{trig}<15~\mathrm{GeV}/c$ in the pseudorapidity range of $|\eta|<0.8$ and with $p_\mathrm{T}>0.5~\mathrm{GeV}/c$ for p-Pb collisions at $\sqrt{s_\mathrm{NN}}=5.02~\mathrm{TeV}$

$p_\mathrm{T}$ as a function of $R_\mathrm{T}$ in the transverse region for different $R_\mathrm{T}$ intervals using events with trigger particles $8<p_\mathrm{T}^\mathrm{trig}<15~\mathrm{GeV}/c$ in the pseudorapidity range of $|\eta|<0.8$ and with $p_\mathrm{T}>0.5~\mathrm{GeV}/c$ for p-Pb collisions at $\sqrt{s_\mathrm{NN}}=5.02~\mathrm{TeV}$

$p_\mathrm{T}$ as a function of $R_\mathrm{T}$ in the toward region for different $R_\mathrm{T}$ intervals using events with trigger particles $8<p_\mathrm{T}^\mathrm{trig}<15~\mathrm{GeV}/c$ in the pseudorapidity range of $|\eta|<0.8$ and with $p_\mathrm{T}>0.5~\mathrm{GeV}/c$ for Pb-Pb collisions at $\sqrt{s_\mathrm{NN}}=5.02~\mathrm{TeV}$

$p_\mathrm{T}$ as a function of $R_\mathrm{T}$ in the away region for different $R_\mathrm{T}$ intervals using events with trigger particles $8<p_\mathrm{T}^\mathrm{trig}<15~\mathrm{GeV}/c$ in the pseudorapidity range of $|\eta|<0.8$ and with $p_\mathrm{T}>0.5~\mathrm{GeV}/c$ for Pb-Pb collisions at $\sqrt{s_\mathrm{NN}}=5.02~\mathrm{TeV}$

$p_\mathrm{T}$ as a function of $R_\mathrm{T}$ in the transverse region for different $R_\mathrm{T}$ intervals using events with trigger particles $8<p_\mathrm{T}^\mathrm{trig}<15~\mathrm{GeV}/c$ in the pseudorapidity range of $|\eta|<0.8$ and with $p_\mathrm{T}>0.5~\mathrm{GeV}/c$ for Pb-Pb collisions at $\sqrt{s_\mathrm{NN}}=5.02~\mathrm{TeV}$

Average $p_\mathrm{T}$ as a function of $R_\mathrm{T}$ in the toward region using events with trigger particles $8<p_\mathrm{T}^\mathrm{trig}<15~\mathrm{GeV}/c$ in the pseudorapidity range of $|\eta|<0.8$ and with $p_\mathrm{T}>0.5~\mathrm{GeV}/c$ for pp collisions at $\sqrt{s_\mathrm{NN}}=5.02~\mathrm{TeV}$

Average $p_\mathrm{T}$ as a function of $R_\mathrm{T}$ in the toward region using events with trigger particles $8<p_\mathrm{T}^\mathrm{trig}<15~\mathrm{GeV}/c$ in the pseudorapidity range of $|\eta|<0.8$ and with $p_\mathrm{T}>0.5~\mathrm{GeV}/c$ for p-Pb collisions at $\sqrt{s_\mathrm{NN}}=5.02~\mathrm{TeV}$

Average $p_\mathrm{T}$ as a function of $R_\mathrm{T}$ in the toward region using events with trigger particles $8<p_\mathrm{T}^\mathrm{trig}<15~\mathrm{GeV}/c$ in the pseudorapidity range of $|\eta|<0.8$ and with $p_\mathrm{T}>0.5~\mathrm{GeV}/c$ for Pb-Pb collisions at $\sqrt{s_\mathrm{NN}}=5.02~\mathrm{TeV}$

Average $p_\mathrm{T}$ as a function of $R_\mathrm{T}$ in the away region using events with trigger particles $8<p_\mathrm{T}^\mathrm{trig}<15~\mathrm{GeV}/c$ in the pseudorapidity range of $|\eta|<0.8$ and with $p_\mathrm{T}>0.5~\mathrm{GeV}/c$ for pp collisions at $\sqrt{s_\mathrm{NN}}=5.02~\mathrm{TeV}$

Average $p_\mathrm{T}$ as a function of $R_\mathrm{T}$ in the away region using events with trigger particles $8<p_\mathrm{T}^\mathrm{trig}<15~\mathrm{GeV}/c$ in the pseudorapidity range of $|\eta|<0.8$ and with $p_\mathrm{T}>0.5~\mathrm{GeV}/c$ for p-Pb collisions at $\sqrt{s_\mathrm{NN}}=5.02~\mathrm{TeV}$

Average $p_\mathrm{T}$ as a function of $R_\mathrm{T}$ in the away region using events with trigger particles $8<p_\mathrm{T}^\mathrm{trig}<15~\mathrm{GeV}/c$ in the pseudorapidity range of $|\eta|<0.8$ and with $p_\mathrm{T}>0.5~\mathrm{GeV}/c$ for Pb-Pb collisions at $\sqrt{s_\mathrm{NN}}=5.02~\mathrm{TeV}$

Average $p_\mathrm{T}$ as a function of $R_\mathrm{T}$ in the transverse region using events with trigger particles $8<p_\mathrm{T}^\mathrm{trig}<15~\mathrm{GeV}/c$ in the pseudorapidity range of $|\eta|<0.8$ and with $p_\mathrm{T}>0.5~\mathrm{GeV}/c$ for pp collisions at $\sqrt{s_\mathrm{NN}}=5.02~\mathrm{TeV}$

Average $p_\mathrm{T}$ as a function of $R_\mathrm{T}$ in the transverse region using events with trigger particles $8<p_\mathrm{T}^\mathrm{trig}<15~\mathrm{GeV}/c$ in the pseudorapidity range of $|\eta|<0.8$ and with $p_\mathrm{T}>0.5~\mathrm{GeV}/c$ for p-Pb collisions at $\sqrt{s_\mathrm{NN}}=5.02~\mathrm{TeV}$

Average $p_\mathrm{T}$ as a function of $R_\mathrm{T}$ in the transverse region using events with trigger particles $8<p_\mathrm{T}^\mathrm{trig}<15~\mathrm{GeV}/c$ in the pseudorapidity range of $|\eta|<0.8$ and with $p_\mathrm{T}>0.5~\mathrm{GeV}/c$ for Pb-Pb collisions at $\sqrt{s_\mathrm{NN}}=5.02~\mathrm{TeV}$

$p_\mathrm{T}$-integrated yield normalised to $<N_\mathrm{ch}^\mathrm{T}>$ as a function of $R_\mathrm{T}$ in the toward region using events with trigger particles $8<p_\mathrm{T}^\mathrm{trig}<15~\mathrm{GeV}/c$ in the pseudorapidity range of $|\eta|<0.8$ and with $p_\mathrm{T}>0.5~\mathrm{GeV}/c$ for pp collisions at $\sqrt{s}=5.02~\mathrm{TeV}$

$p_\mathrm{T}$-integrated yield normalised to $<N_\mathrm{ch}^\mathrm{T}>$ as a function of $R_\mathrm{T}$ in the away region using events with trigger particles $8<p_\mathrm{T}^\mathrm{trig}<15~\mathrm{GeV}/c$ in the pseudorapidity range of $|\eta|<0.8$ and with $p_\mathrm{T}>0.5~\mathrm{GeV}/c$ for pp collisions at $\sqrt{s}=5.02~\mathrm{TeV}$

$p_\mathrm{T}$-integrated yield normalised to $<N_\mathrm{ch}^\mathrm{T}>$ as a function of $R_\mathrm{T}$ in the toward region using events with trigger particles $8<p_\mathrm{T}^\mathrm{trig}<15~\mathrm{GeV}/c$ in the pseudorapidity range of $|\eta|<0.8$ and with $p_\mathrm{T}>0.5~\mathrm{GeV}/c$ for p-Pb collisions at $\sqrt{s_\mathrm{NN}}=5.02~\mathrm{TeV}$

$p_\mathrm{T}$-integrated yield normalised to $<N_\mathrm{ch}^\mathrm{T}>$ as a function of $R_\mathrm{T}$ in the away region using events with trigger particles $8<p_\mathrm{T}^\mathrm{trig}<15~\mathrm{GeV}/c$ in the pseudorapidity range of $|\eta|<0.8$ and with $p_\mathrm{T}>0.5~\mathrm{GeV}/c$ for p-Pb collisions at $\sqrt{s_\mathrm{NN}}=5.02~\mathrm{TeV}$

$p_\mathrm{T}$-integrated yield normalised to $<N_\mathrm{ch}^\mathrm{T}>$ as a function of $R_\mathrm{T}$ in the toward region using events with trigger particles $8<p_\mathrm{T}^\mathrm{trig}<15~\mathrm{GeV}/c$ in the pseudorapidity range of $|\eta|<0.8$ and with $p_\mathrm{T}>0.5~\mathrm{GeV}/c$ for Pb-Pb collisions at $\sqrt{s_\mathrm{NN}}=5.02~\mathrm{TeV}$

$p_\mathrm{T}$-integrated yield normalised to $<N_\mathrm{ch}^\mathrm{T}>$ as a function of $R_\mathrm{T}$ in the away region using events with trigger particles $8<p_\mathrm{T}^\mathrm{trig}<15~\mathrm{GeV}/c$ in the pseudorapidity range of $|\eta|<0.8$ and with $p_\mathrm{T}>0.5~\mathrm{GeV}/c$ for Pb-Pb collisions at $\sqrt{s_\mathrm{NN}}=5.02~\mathrm{TeV}$