Likelihood profiles of the $\kappa_\lambda$ parameter.

Likelihood profiles of the $\kappa_{2V}$ parameter.

The 95% CL upper limits on the $HH$ cross-section for different values of $\kappa_\lambda$. The expected limits assume no $HH$ production. The red line shows the theory prediction for the combined ggF $HH$ and VBF $HH$ cross-section as a function of $\kappa_\lambda$ where all parameters and couplings are set to their SM values except $\kappa_\lambda$. The bands surrounding the red cross-section line indicate the theoretical uncertainty on the predicted cross-section. The $\kappa_\lambda$ values in the interval $[-6.7, 15.3]$ cannot be excluded by this analysis.

The 95% CL upper limits on the VBF $HH$ cross section for different values of $\kappa_{2V}$. The expected limits assume no VBF $HH$ production. SM ggF $HH$ production is considered as background. The red line shows the theory prediction for the VBF $HH$ cross-section as a function of $\kappa_{2V}$ where all parameters and couplings are set to their SM values except $\kappa_{2V}$. The uncertainty band is smaller than the width of the plotted line. The $\kappa_{2V}$ values in the interval $[-0.3, 2.4]$ cannot be excluded by this analysis.

Correlation matrix of all the parameters including nuisance parameters used in the profile likelihood unfolding fit on $A_{C}$.

Observed and expected (background-only fitted) invariant mass spectra of Wgamma events. The fitted function is ${ d N}/{ d m} = p_{0} * (m/\sqrt{s})^{p_{1} + p_{2} * \log(m/\sqrt{s}) + p_{3} * \log^{2}(m/\sqrt{s})}$

Expected and observed 95% CL upper limits on the product of the cross section and branching fraction for narrow scalar Wgamma resonances. Limits are compared to predicted cross sections for the heavy scalar triplet model described in arXiv:1912.08234

Expected and observed 95% CL upper limits on the product of the cross section and branching fraction for broad scalar Wgamma resonances.

Expected and observed 95% CL upper limits on the product of the cross section and branching fraction for narrow vector Wgamma resonances. Limits are compared to predicted cross sections for the heavy vector triplet model described in arXiv:1912.08234

Expected and observed 95% CL upper limits on the product of the cross section and branching fraction for broad vector Wgamma resonances.

Expected and observed model-independent 95% CL upper limits on the product of the cross section, branching fraction and signal acceptance for general Wgamma resonances.

Expected and observed model-independent 95% CL upper limits on the product of the cross section, branching fraction, signal acceptance and W tagging efficiency for general Jgamma resonances.

Definition of the tetralepton signal regions.

Definition of the fiducial volumes at particle- and parton-level. Leptons refer exclusively to electrons and muons - they are dressed with additional radiation at particle-level, but not at parton-level.

Definition of the dilepton $t\bar{t}$ validation regions.

Pre-fit distribution of the number of $b$-jets in 2L-$e\mu$-6j2b, this distribution is not used in the fit.

Pre-fit distribution of the DNN output 2L-$e\mu$-6j1b, this distribution is not used in the fit.

Pre-fit distribution of the DNN output 2L-$e\mu$-5j2b, this distribution is not used in the fit.

Pre-fit distribution of the DNN output 2L-$e\mu$-6j2b, this distribution is not used in the fit.

Definition of the tetralepton control region.

Definition of the trilepton fakes control regions.

Pre-fit distribution of jet multiplicity in CR-$t\bar{t}$-e region.

Pre-fit distribution of loose lepton transverse momentum in CR-$t\bar{t}$-$\mu$ region.

Pre-fit distribution of the transverse mass of the trailing lepton and the missing transverse momentum in CR-Z-e region.

Post-fit distribution of jet multiplicity in CR-$t\bar{t}$-e region

Post-fit distribution of loose lepton transverse momentum in CR-$t\bar{t}$-$\mu$ region

Post-fit distribution of the transverse mass of the trailing lepton and the missing transverse momentum in CR-Z-e region

Post-fit distribution of NN output in SR-2L-5j2b region.

Post-fit distribution of NN output in SR-2L-6j1b region.

Post-fit distribution of NN output in SR-2L-6j2b region.

Post-fit distribution of DNN-$t\bar{t}Z$ output in 3L-SR-ttZ region.

Post-fit distribution of DNN-$t\bar{t}Z$ outputt in 3L-SR-tZq region.

Post fit events yields in 3L-SR-WZ region.

Post-fit distribution of NN output in 4L-SR-SF region.

Post-fit distribution of NN output in 4L-SR-DF region.

Post-fit distribution of b-tagger output for leading b-jet in 4L-CR-ZZ region.

Measured values of the background normalizations obtained from the combined fit. The uncertainties include statistical and systematic sources.

Measured $\sigma_{t\bar{t}\text{Z}}$ cross sections obtained from the fits in the different lepton channels. The uncertainties include statistical and systematic sources.

Grouped impact of systematic uncertainties in the combined inclusive fit to data.

Unfolded absolute cross section as a function of $p^{Z}_{T}$ in the combination of $3\ell$ and $4\ell$ channels at particle-level (Figure 8 top-left).

Unfolded absolute cross section as a function of $p^{Z}_{T}$ in the combination of $3\ell$ and $4\ell$ channels at parton-level (Figure 8 top-right).

Unfolded normalized cross section as a function of $p^{Z}_{T}$ in the combination of $3\ell$ and $4\ell$ channels at particle-level (Figure 8 bottom-left).

Unfolded normalized cross section as a function of $p^{Z}_{T}$ in the combination of $3\ell$ and $4\ell$ channels at parton-level (Figure 8 bottom-right).

Unfolded absolute cross section as a function of $|y^{Z}$| in the combination of $3\ell$ and $4\ell$ channels at particle-level (Figure 17 top-left and Figure 11 top-left).

Unfolded absolute cross section as a function of $|y^{Z}$| in the combination of $3\ell$ and $4\ell$ channels at parton-level (Figure 17 top-right).

Unfolded normalized cross section as a function of $|y^{Z}$| in the combination of $3\ell$ and $4\ell$ channels at particle-level (Figure 17 bottom-left).

Unfolded normalized cross section as a function of $|y^{Z}$| in the combination of $3\ell$ and $4\ell$ channels at parton-level (Figure 17 bottom-right).

Unfolded absolute cross section as a function of cos $\theta_{Z}^{*}$ in the combination of $3\ell$ and $4\ell$ channels at particle-level (Figure 18 top-left and Figure 11 top-right).

Unfolded absolute cross section as a function of cos $\theta_{Z}^{*}$ in the combination of $3\ell$ and $4\ell$ channels at parton-level (Figure 18 top-right).

Unfolded normalized cross section as a function of cos $\theta_{Z}^{*}$ in the combination of $3\ell$ and $4\ell$ channels at particle-level (Figure 18 bottom-left).

Unfolded normalized cross section as a function of cos $\theta_{Z}^{*}$ in the combination of $3\ell$ and $4\ell$ channels at parton-level (Figure 18 bottom-right).

Unfolded absolute cross section as a function of $p_{T}^{\mathrm{top}}$ in the combination of $3\ell$ and $4\ell$ channels at particle-level (Figure 19 top-left and Figure 11 bottom-left).

Unfolded absolute cross section as a function of $p_{T}^{\mathrm{top}}$ in the combination of $3\ell$ and $4\ell$ channels at parton-level (Figure 19, top-right).

Unfolded normalized cross section as a function of $p_{T}^{\mathrm{top}}$ in the combination of $3\ell$ and $4\ell$ channels at particle-level (Figure 19, bottom-left).

Unfolded normalized cross section as a function of $p_{T}^{\mathrm{top}}$ in the combination of $3\ell$ and $4\ell$ channels at parton-level (Figure 19, bottom-right).

Unfolded absolute cross section as a function of $p_{T}^{t\bar{t}}$ in the combination of $3\ell$ and $4\ell$ channels at particle-level (Figure 20 top-left and Figure 11 bottom-right).

Unfolded absolute cross section as a function of $p_{T}^{t\bar{t}}$ in the combination of $3\ell$ and $4\ell$ channels at parton-level (Figure 20, top-right).

Unfolded normalized cross section as a function of $p_{T}^{t\bar{t}}$ in the combination of $3\ell$ and $4\ell$ channels at particle-level (Figure 20, bottom-left)

Unfolded normalized cross section as a function of $p_{T}^{t\bar{t}}$ in the combination of $3\ell$ and $4\ell$ channels at parton-level (Figure 20, bottom-right)

Unfolded absolute cross section as a function of $|\Delta\Phi(t\bar{t}, Z)|/\pi$ in the combination of $3\ell$ and $4\ell$ channels at particle-level (Figure 21 top-left and Figure 12 top-left).

Unfolded absolute cross section as a function of $|\Delta\Phi(t\bar{t}, Z)|/\pi$ in the combination of $3\ell$ and $4\ell$ channels at parton-level (Figure 21, top-right).

Unfolded normalized cross section as a function of $|\Delta\Phi(t\bar{t}, Z)|/\pi$ in the combination of $3\ell$ and $4\ell$ channels at particle-level (Figure 21, bottom-left).

Unfolded normalized cross section as a function of $|\Delta\Phi(t\bar{t}, Z)|/\pi$ in the combination of $3\ell$ and $4\ell$ channels at parton-level (Figure 21, top-right).

Unfolded absolute cross section as a function of $m^{t\bar{t}Z}$ in the combination of $3\ell$ and $4\ell$ channels at particle-level (Figure 22 top-left and Figure 12 bottom-left).

Unfolded absolute cross section as a function of $m^{t\bar{t}Z}$ in the combination of $3\ell$ and $4\ell$ channels at parton-level (Figure 22, top-right).

Unfolded normalized cross section as a function of $m^{t\bar{t}Z}$ in the combination of $3\ell$ and $4\ell$ channels at particle-level (Figure 22, bottom-left).

Unfolded normalized cross section as a function of $m^{t\bar{t}Z}$ in the combination of $3\ell$ and $4\ell$ channels at parton-level (Figure 22, bottom-right).

Unfolded absolute cross section as a function of $m^{t\bar{t}}$ in the combination of $3\ell$ and $4\ell$ channels at particle-level (Figure 23 top-left and Figure 12 bottom-right).

Unfolded absolute cross section as a function of $m^{t\bar{t}}$ in the combination of $3\ell$ and $4\ell$ channels at parton-level (Figure 23, top-right).

Unfolded normalized cross section as a function of $m^{t\bar{t}}$ in the combination of $3\ell$ and $4\ell$ channels at particle-level (Figure 23, bottom-left).

Unfolded normalized cross section as a function of $m^{t\bar{t}}$ in the combination of $3\ell$ and $4\ell$ channels at parton-level (Figure 23, bottom-right).

Unfolded absolute cross section as a function of $|y^{t\bar{t}Z}|$ in the combination of $3\ell$ and $4\ell$ channels at particle-level (Figure 24 top-left and Figure 12 top-right).

Unfolded absolute cross section as a function of $|y^{t\bar{t}Z}|$ in the combination of $3\ell$ and $4\ell$ channels at parton-level (Figure 24, top-right).

Unfolded normalized cross section as a function of $|y^{t\bar{t}Z}|$ in the combination of $3\ell$ and $4\ell$ channels at particle-level (Figure 24, bottom-left).

Unfolded normalized cross section as a function of $|y^{t\bar{t}Z}|$ in the combination of $3\ell$ and $4\ell$ channels at parton-level (Figure 24, bottom-right).

Unfolded absolute cross section as a function of $H_{\text{T}}^{\text{l}}$ in the trilepton channel at particle-level (Figure 25 top-left and Figure 9 top-left).

Unfolded absolute cross section as a function of $H_{\text{T}}^{\text{l}}$ in the trilepton channel at parton-level (Figure 25 top-right).

Unfolded normalized cross section as a function of $H_{\text{T}}^{\text{l}}$ in the trilepton channel at particle-level (Figure 25 bottom-left).

Unfolded normalized cross section as a function of $H_{\text{T}}^{\text{l}}$ in the trilepton channel at parton-level (Figure 25 bottom-right).

Unfolded absolute cross section as a function of $|\Delta\Phi(Z, t_{lep})|/\pi$ in the trilepton channel at particle-level (Figure 26 top-left and Figure 10 bottom-left).

Unfolded absolute cross section as a function of $|\Delta\Phi(Z, t_{lep})|/\pi$ in the trilepton channel at parton-level (Figure 26 top-right).

Unfolded normalized cross section as a function of $|\Delta\Phi(Z, t_{lep})|/\pi$ in the trilepton channel at particle-level (Figure 26 bottom-left).

Unfolded normalized cross section as a function of $|\Delta\Phi(Z, t_{lep})|/\pi$ in the trilepton channel at parton-level (Figure 26 bottom-right).

Unfolded absolute cross section as a function of $|\Delta y(Z, t_{lep})|$ in the trilepton channel at particle-level (Figure 27 top-left and Figure 10 bottom-right).

Unfolded absolute cross section as a function of $|\Delta y(Z, t_{lep})|$ in the trilepton channel at parton-level (Figure 27 top-right).

Unfolded normalized cross section as a function of $|\Delta y(Z, t_{lep})|$ in the trilepton channel at particle-level (Figure 27 bottom-left).

Unfolded normalized cross section as a function of $|\Delta y(Z, t_{lep})|$ in the trilepton channel at parton-level (Figure 27 bottom-right).

Unfolded absolute cross section as a function of $p_{\text{T}}^{\ell, non-Z}$ in the trilepton channel at particle-level (Figure 28 top-left and Figure 10 top-left).

Unfolded absolute cross section as a function of $p_{\text{T}}^{\ell, non-Z}$ in the trilepton channel at parton-level (Figure 28 top-right).

Unfolded normalized cross section as a function of $p_{\text{T}}^{\ell, non-Z}$ in the trilepton channel at particle-level (Figure 28 bottom-left).

Unfolded normalized cross section as a function of $p_{\text{T}}^{\ell, non-Z}$ in the trilepton channel at parton-level (Figure 28 bottom-right).

Unfolded absolute cross section as a function of $N_{\text{jets}}$ in the trilepton channel at particle-level (Figure 29 left and Figure 9 bottom-left).

Unfolded normalized cross section as a function of $N_{\text{jets}}$ in the trilepton channel at particle-level (Figure 29 right).

Unfolded absolute cross section as a function of $H_{\text{T}}^{\text{l}}$ in the tetralepton channel at particle-level (Figure 30 top-left and Figure 9 top-right).

Unfolded absolute cross section as a function of $H_{\text{T}}^{\text{l}}$ in the tetralepton channel at parton-level (Figure 30 top-right).

Unfolded normalized cross section as a function of $H_{\text{T}}^{\text{l}}$ in the tetralepton channel at particle-level (Figure 30 bottom-left).

Unfolded normalized cross section as a function of $H_{\text{T}}^{\text{l}}$ in the tetralepton channel at parton-level (Figure 30 bottom-right).

Unfolded absolute cross section as a function of $|\Delta\Phi(l_{t}^{+}, l_{\bar{t}}^{-})|/\pi$ in the tetralepton channel at particle-level (Figure 31 top-left and Figure 10 top-right).

Unfolded absolute cross section as a function of $|\Delta\Phi(l_{t}^{+}, l_{\bar{t}}^{-})|/\pi$ in the tetralepton channel at parton-level (Figure 31 top-right).

Unfolded normalized cross section as a function of $|\Delta\Phi(l_{t}^{+}, l_{\bar{t}}^{-})|/\pi$ in the tetralepton channel at particle-level (Figure 31 bottom-left).

Unfolded normalized cross section as a function of $|\Delta\Phi(l_{t}^{+}, l_{\bar{t}}^{-})|/\pi$ in the tetralepton channel at parton-level (Figure 31 bottom-right).

Unfolded absolute cross section as a function of $N_{\text{jets}}$ in the tetralepton channel at particle-level (Figure 32 left and Figure 9 bottom-right).

Unfolded normalized cross section as a function of $N_{\text{jets}}$ in the tetralepton channel at particle-level (Figure 32 right).

Bootstrap replicas (0-499) for data in all regions used in inclusive cross section measurement. The used bootstrap method is described in ATL-PHYS-PUB-2021-011 (https://atlas.web.cern.ch/Atlas/GROUPS/PHYSICS/PUBNOTES/ATL-PHYS-PUB-2021-011/).

Bootstrap replicas (500-999) for data in all regions used in inclusive cross section measurement. The used bootstrap method is described in ATL-PHYS-PUB-2021-011 (https://atlas.web.cern.ch/Atlas/GROUPS/PHYSICS/PUBNOTES/ATL-PHYS-PUB-2021-011/).

Bootstrap replicas (0-499) for data, variable $|\Delta\Phi(t\bar{t}, Z)|/\pi$. The used bootstrap method is described in ATL-PHYS-PUB-2021-011 (https://atlas.web.cern.ch/Atlas/GROUPS/PHYSICS/PUBNOTES/ATL-PHYS-PUB-2021-011/).

Bootstrap replicas (500-999) for data, variable $|\Delta\Phi(t\bar{t}, Z)|/\pi$. The used bootstrap method is described in ATL-PHYS-PUB-2021-011 (https://atlas.web.cern.ch/Atlas/GROUPS/PHYSICS/PUBNOTES/ATL-PHYS-PUB-2021-011/).

Bootstrap replicas (0-499) for data, variable $|\Delta\Phi(Z, t_{lep})|/\pi$. The used bootstrap method is described in ATL-PHYS-PUB-2021-011 (https://atlas.web.cern.ch/Atlas/GROUPS/PHYSICS/PUBNOTES/ATL-PHYS-PUB-2021-011/).

Bootstrap replicas (500-999) for data, variable $|\Delta\Phi(Z, t_{lep})|/\pi$. The used bootstrap method is described in ATL-PHYS-PUB-2021-011 (https://atlas.web.cern.ch/Atlas/GROUPS/PHYSICS/PUBNOTES/ATL-PHYS-PUB-2021-011/).

Bootstrap replicas (0-499) for data, variable $m^{t\bar{t}}$. The used bootstrap method is described in ATL-PHYS-PUB-2021-011 (https://atlas.web.cern.ch/Atlas/GROUPS/PHYSICS/PUBNOTES/ATL-PHYS-PUB-2021-011/).

Bootstrap replicas (500-999) for data, variable $m^{t\bar{t}}$. The used bootstrap method is described in ATL-PHYS-PUB-2021-011 (https://atlas.web.cern.ch/Atlas/GROUPS/PHYSICS/PUBNOTES/ATL-PHYS-PUB-2021-011/).

Bootstrap replicas (0-499) for data, variable $N_{\text{jets}}$ in $3\ell$ channel. The used bootstrap method is described in ATL-PHYS-PUB-2021-011 (https://atlas.web.cern.ch/Atlas/GROUPS/PHYSICS/PUBNOTES/ATL-PHYS-PUB-2021-011/).

Bootstrap replicas (500-999) for data, variable $N_{\text{jets}}$ in $3\ell$ channel. The used bootstrap method is described in ATL-PHYS-PUB-2021-011 (https://atlas.web.cern.ch/Atlas/GROUPS/PHYSICS/PUBNOTES/ATL-PHYS-PUB-2021-011/).

Bootstrap replicas (0-499) for data, variable $|y^{t\bar{t}Z}|$. The used bootstrap method is described in ATL-PHYS-PUB-2021-011 (https://atlas.web.cern.ch/Atlas/GROUPS/PHYSICS/PUBNOTES/ATL-PHYS-PUB-2021-011/).

Bootstrap replicas (500-999) for data, variable $|y^{t\bar{t}Z}|$. The used bootstrap method is described in ATL-PHYS-PUB-2021-011 (https://atlas.web.cern.ch/Atlas/GROUPS/PHYSICS/PUBNOTES/ATL-PHYS-PUB-2021-011/).

Bootstrap replicas (0-499) for data, variable $H_{\text{T}}^{\text{l}}$ in $3\ell$ channel. The used bootstrap method is described in ATL-PHYS-PUB-2021-011 (https://atlas.web.cern.ch/Atlas/GROUPS/PHYSICS/PUBNOTES/ATL-PHYS-PUB-2021-011/).

Bootstrap replicas (500-999) for data, variable $H_{\text{T}}^{\text{l}}$ in $3\ell$ channel. The used bootstrap method is described in ATL-PHYS-PUB-2021-011 (https://atlas.web.cern.ch/Atlas/GROUPS/PHYSICS/PUBNOTES/ATL-PHYS-PUB-2021-011/).

Bootstrap replicas (0-499) for data, variable $y^{Z}$. The used bootstrap method is described in ATL-PHYS-PUB-2021-011 (https://atlas.web.cern.ch/Atlas/GROUPS/PHYSICS/PUBNOTES/ATL-PHYS-PUB-2021-011/).

Bootstrap replicas (500-999) for data, variable $y^{Z}$. The used bootstrap method is described in ATL-PHYS-PUB-2021-011 (https://atlas.web.cern.ch/Atlas/GROUPS/PHYSICS/PUBNOTES/ATL-PHYS-PUB-2021-011/).

Bootstrap replicas (0-499) for data, variable $p_{T}^{\mathrm{top}}$. The used bootstrap method is described in ATL-PHYS-PUB-2021-011 (https://atlas.web.cern.ch/Atlas/GROUPS/PHYSICS/PUBNOTES/ATL-PHYS-PUB-2021-011/).

Bootstrap replicas (500-999) for data, variable $p_{T}^{\mathrm{top}}$. The used bootstrap method is described in ATL-PHYS-PUB-2021-011 (https://atlas.web.cern.ch/Atlas/GROUPS/PHYSICS/PUBNOTES/ATL-PHYS-PUB-2021-011/).

Bootstrap replicas (0-499) for data, variable cos $\theta^{*}_{Z}$. The used bootstrap method is described in ATL-PHYS-PUB-2021-011 (https://atlas.web.cern.ch/Atlas/GROUPS/PHYSICS/PUBNOTES/ATL-PHYS-PUB-2021-011/).

Bootstrap replicas (500-999) for data, variable cos $\theta^{*}_{Z}$. The used bootstrap method is described in ATL-PHYS-PUB-2021-011 (https://atlas.web.cern.ch/Atlas/GROUPS/PHYSICS/PUBNOTES/ATL-PHYS-PUB-2021-011/).

Bootstrap replicas (0-499) for data, variable $p_{\text{T}}^{\ell, non-Z}$. The used bootstrap method is described in ATL-PHYS-PUB-2021-011 (https://atlas.web.cern.ch/Atlas/GROUPS/PHYSICS/PUBNOTES/ATL-PHYS-PUB-2021-011/).

Bootstrap replicas (500-999) for data, variable $p_{\text{T}}^{\ell, non-Z}$. The used bootstrap method is described in ATL-PHYS-PUB-2021-011 (https://atlas.web.cern.ch/Atlas/GROUPS/PHYSICS/PUBNOTES/ATL-PHYS-PUB-2021-011/).

Bootstrap replicas (0-499) for data, variable $|\Delta\Phi(l_{t}^{+}, l_{\bar{t}}^{-})|/\pi$. The used bootstrap method is described in ATL-PHYS-PUB-2021-011 (https://atlas.web.cern.ch/Atlas/GROUPS/PHYSICS/PUBNOTES/ATL-PHYS-PUB-2021-011/).

Bootstrap replicas (500-999) for data, variable $|\Delta\Phi(l_{t}^{+}, l_{\bar{t}}^{-})|/\pi$. The used bootstrap method is described in ATL-PHYS-PUB-2021-011 (https://atlas.web.cern.ch/Atlas/GROUPS/PHYSICS/PUBNOTES/ATL-PHYS-PUB-2021-011/).

Bootstrap replicas (0-499) for data, variable $H_{\text{T}}^{\text{l}}$ in $4\ell$ channel. The used bootstrap method is described in ATL-PHYS-PUB-2021-011 (https://atlas.web.cern.ch/Atlas/GROUPS/PHYSICS/PUBNOTES/ATL-PHYS-PUB-2021-011/).

Bootstrap replicas (500-999) for data, variable $H_{\text{T}}^{\text{l}}$ in $4\ell$ channel. The used bootstrap method is described in ATL-PHYS-PUB-2021-011 (https://atlas.web.cern.ch/Atlas/GROUPS/PHYSICS/PUBNOTES/ATL-PHYS-PUB-2021-011/).

Bootstrap replicas (0-499) for data, variable $m^{t\bar{t}Z}$. The used bootstrap method is described in ATL-PHYS-PUB-2021-011 (https://atlas.web.cern.ch/Atlas/GROUPS/PHYSICS/PUBNOTES/ATL-PHYS-PUB-2021-011/).

Bootstrap replicas (500-999) for data, variable $m^{t\bar{t}Z}$. The used bootstrap method is described in ATL-PHYS-PUB-2021-011 (https://atlas.web.cern.ch/Atlas/GROUPS/PHYSICS/PUBNOTES/ATL-PHYS-PUB-2021-011/).

Bootstrap replicas (0-499) for data, variable $N_{\text{jets}}$ in $4\ell$ channel. The used bootstrap method is described in ATL-PHYS-PUB-2021-011 (https://atlas.web.cern.ch/Atlas/GROUPS/PHYSICS/PUBNOTES/ATL-PHYS-PUB-2021-011/).

Bootstrap replicas (500-999) for data, variable $N_{\text{jets}}$ in $4\ell$ channel. The used bootstrap method is described in ATL-PHYS-PUB-2021-011 (https://atlas.web.cern.ch/Atlas/GROUPS/PHYSICS/PUBNOTES/ATL-PHYS-PUB-2021-011/).

Bootstrap replicas (0-499) for data, variable $|\Delta y(Z, t_{lep})|$. The used bootstrap method is described in ATL-PHYS-PUB-2021-011 (https://atlas.web.cern.ch/Atlas/GROUPS/PHYSICS/PUBNOTES/ATL-PHYS-PUB-2021-011/).

Bootstrap replicas (500-999) for data, variable $|\Delta y(Z, t_{lep})|$. The used bootstrap method is described in ATL-PHYS-PUB-2021-011 (https://atlas.web.cern.ch/Atlas/GROUPS/PHYSICS/PUBNOTES/ATL-PHYS-PUB-2021-011/).

Bootstrap replicas (0-499) for data, variable $p^{Z}_{T}$. The used bootstrap method is described in ATL-PHYS-PUB-2021-011 (https://atlas.web.cern.ch/Atlas/GROUPS/PHYSICS/PUBNOTES/ATL-PHYS-PUB-2021-011/).

Bootstrap replicas (500-999) for data, variable $p^{Z}_{T}$. The used bootstrap method is described in ATL-PHYS-PUB-2021-011 (https://atlas.web.cern.ch/Atlas/GROUPS/PHYSICS/PUBNOTES/ATL-PHYS-PUB-2021-011/).

Bootstrap replicas (0-499) for data, variable $p_{T}^{t\bar{t}}$. The used bootstrap method is described in ATL-PHYS-PUB-2021-011 (https://atlas.web.cern.ch/Atlas/GROUPS/PHYSICS/PUBNOTES/ATL-PHYS-PUB-2021-011/).

Bootstrap replicas (500-999) for data, variable $p_{T}^{t\bar{t}}$. The used bootstrap method is described in ATL-PHYS-PUB-2021-011 (https://atlas.web.cern.ch/Atlas/GROUPS/PHYSICS/PUBNOTES/ATL-PHYS-PUB-2021-011/).

Parton-level acceptance and selection efficiency histograms for $|\Delta\Phi(Z, t_{lep})|/\pi$ variable.

Parton-level acceptance and selection efficiency histograms for $|\Delta y(Z, t_{lep})|$ variable.

Parton-level acceptance and selection efficiency histograms for $H_{\text{T}}^{\text{ l}}$ variable.

Parton-level acceptance and selection efficiency histograms for $p_{\text{T}}^{\ell, non-Z}$ variable.

Parton-level acceptance and selection efficiency histograms for $|\Delta\Phi(l_{t}^{+}, l_{\bar{t}}^{-})|/\pi$ variable.

Parton-level acceptance and selection efficiency histograms for $H_{\text{T}}^{\text{ l}}$ variable.

Parton-level acceptance and selection efficiency histograms for cos $\theta_{Z}^{*}$ variable.

Parton-level acceptance and selection efficiency histograms for $p^{Z}_{T}$ variable.

Parton-level acceptance and selection efficiency histograms for $|y^{Z}$| variable.

Parton-level acceptance and selection efficiency histograms for $|\Delta\Phi(t\bar{t}, Z)|/\pi$ variable.

Parton-level acceptance and selection efficiency histograms for $m^{t\bar{t}}$ variable.

Parton-level acceptance and selection efficiency histograms for $m^{t\bar{t}Z}$ variable.

Parton-level acceptance and selection efficiency histograms for $p_{T}^{\mathrm{top}}$ variable.

Parton-level acceptance and selection efficiency histograms for $p_{T}^{t\bar{t}}$ variable.

Parton-level acceptance and selection efficiency histograms for $|y^{t\bar{t}Z}|$ variable.

Particle-level acceptance and selection efficiency histograms for $|\Delta\Phi(Z, t_{lep})|/\pi$ variable.

Particle-level acceptance and selection efficiency histograms for $|\Delta y(Z, t_{lep})|$ variable.

Particle-level acceptance and selection efficiency histograms for $H_{\text{T}}^{\text{ l}}$ variable.

Particle-level acceptance and selection efficiency histograms for $N_{\text{jets}}$ variable.

Particle-level acceptance and selection efficiency histograms for $p_{\text{T}}^{\ell, non-Z}$ variable.

Particle-level acceptance and selection efficiency histograms for $|\Delta\Phi(l_{t}^{+}, l_{\bar{t}}^{-})|/\pi$ variable.

Particle-level acceptance and selection efficiency histograms for $H_{\text{T}}^{\text{ l}}$ variable.

Particle-level acceptance and selection efficiency histograms for $N_{\text{jets}}$ variable.

Particle-level acceptance and selection efficiency histograms for cos $\theta_{Z}^{*}$ variable.

Particle-level acceptance and selection efficiency histograms for $p^{Z}_{T}$ variable.

Particle-level acceptance and selection efficiency histograms for $|y^{Z}$| variable.

Particle-level acceptance and selection efficiency histograms for $|\Delta\Phi(t\bar{t}, Z)|/\pi$ variable.

Particle-level acceptance and selection efficiency histograms for $m^{t\bar{t}}$ variable.

Particle-level acceptance and selection efficiency histograms for $m^{t\bar{t}Z}$ variable.

Particle-level acceptance and selection efficiency histograms for $p_{T}^{\mathrm{top}}$ variable.

Particle-level acceptance and selection efficiency histograms for $p_{T}^{t\bar{t}}$ variable.

Particle-level acceptance and selection efficiency histograms for $|y^{t\bar{t}Z}|$ variable.

Migration matrix for cos $\theta_{Z}^{*}$ variable at particle-level in region SR-3L-ttZ.

Migration matrix for cos $\theta_{Z}^{*}$ variable at particle-level in region SR-3L-tZq.

Migration matrix for cos $\theta_{Z}^{*}$ variable at particle-level in region SR-3L-WZ.

Migration matrix for cos $\theta_{Z}^{*}$ variable at particle-level in region SR-4L-DF.

Migration matrix for cos $\theta_{Z}^{*}$ variable at particle-level in region SR-4L-SF.

Migration matrix for cos $\theta_{Z}^{*}$ variable at particle-level in region CR-4L-ZZ.

Migration matrix for cos $\theta_{Z}^{*}$ variable at parton-level in region SR-3L-ttZ.

Migration matrix for cos $\theta_{Z}^{*}$ variable at parton-level in region SR-3L-tZq.

Migration matrix for cos $\theta_{Z}^{*}$ variable at parton-level in region SR-3L-WZ.

Migration matrix for cos $\theta_{Z}^{*}$ variable at parton-level in region SR-4L-DF.

Migration matrix for cos $\theta_{Z}^{*}$ variable at parton-level in region SR-4L-SF.

Migration matrix for cos $\theta_{Z}^{*}$ variable at parton-level in region CR-4L-ZZ.

Migration matrix for $|\Delta\Phi(t\bar{t}, Z)|/\pi$ variable at particle-level in region SR-3L-ttZ.

Migration matrix for $|\Delta\Phi(t\bar{t}, Z)|/\pi$ variable at particle-level in region SR-3L-tZq.

Migration matrix for $|\Delta\Phi(t\bar{t}, Z)|/\pi$ variable at particle-level in region SR-3L-WZ.

Migration matrix for $|\Delta\Phi(t\bar{t}, Z)|/\pi$ variable at particle-level in region SR-4L-DF.

Migration matrix for $|\Delta\Phi(t\bar{t}, Z)|/\pi$ variable at particle-level in region SR-4L-SF.

Migration matrix for $|\Delta\Phi(t\bar{t}, Z)|/\pi$ variable at particle-level in region CR-4L-ZZ.

Migration matrix for $|\Delta\Phi(t\bar{t}, Z)|/\pi$ variable at parton-level in region SR-3L-ttZ.

Migration matrix for $|\Delta\Phi(t\bar{t}, Z)|/\pi$ variable at parton-level in region SR-3L-tZq.

Migration matrix for $|\Delta\Phi(t\bar{t}, Z)|/\pi$ variable at parton-level in region SR-3L-WZ.

Migration matrix for $|\Delta\Phi(t\bar{t}, Z)|/\pi$ variable at parton-level in region SR-4L-DF.

Migration matrix for $|\Delta\Phi(t\bar{t}, Z)|/\pi$ variable at parton-level in region SR-4L-SF.

Migration matrix for $|\Delta\Phi(t\bar{t}, Z)|/\pi$ variable at parton-level in region CR-4L-ZZ.

Migration matrix for $|\Delta\Phi(l_{t}^{+}, l_{\bar{t}}^{-})|/\pi$ variable at particle-level in region SR-4L-DF.

Migration matrix for $|\Delta\Phi(l_{t}^{+}, l_{\bar{t}}^{-})|/\pi$ variable at particle-level in region SR-4L-SF.

Migration matrix for $|\Delta\Phi(l_{t}^{+}, l_{\bar{t}}^{-})|/\pi$ variable at particle-level in region CR-4L-ZZ.

Migration matrix for $|\Delta\Phi(l_{t}^{+}, l_{\bar{t}}^{-})|/\pi$ variable at parton-level in region SR-4L-DF.

Migration matrix for $|\Delta\Phi(l_{t}^{+}, l_{\bar{t}}^{-})|/\pi$ variable at parton-level in region SR-4L-SF.

Migration matrix for $|\Delta\Phi(l_{t}^{+}, l_{\bar{t}}^{-})|/\pi$ variable at parton-level in region CR-4L-ZZ.

Migration matrix for $|\Delta\Phi(Z, t_{lep})|/\pi$ variable at particle-level in region SR-3L-ttZ.

Migration matrix for $|\Delta\Phi(Z, t_{lep})|/\pi$ variable at particle-level in region SR-3L-tZq.

Migration matrix for $|\Delta\Phi(Z, t_{lep})|/\pi$ variable at particle-level in region SR-3L-WZ.

Migration matrix for $|\Delta\Phi(Z, t_{lep})|/\pi$ variable at parton-level in region SR-3L-ttZ.

Migration matrix for $|\Delta\Phi(Z, t_{lep})|/\pi$ variable at parton-level in region SR-3L-tZq.

Migration matrix for $|\Delta\Phi(Z, t_{lep})|/\pi$ variable at parton-level in region SR-3L-WZ.

Migration matrix for $|\Delta y(Z, t_{lep})|$ variable at particle-level in region SR-3L-ttZ.

Migration matrix for $|\Delta y(Z, t_{lep})|$ variable at particle-level in region SR-3L-tZq.

Migration matrix for $|\Delta y(Z, t_{lep})|$ variable at particle-level in region SR-3L-WZ.

Migration matrix for $|\Delta y(Z, t_{lep})|$ variable at parton-level in region SR-3L-ttZ.

Migration matrix for $|\Delta y(Z, t_{lep})|$ variable at parton-level in region SR-3L-tZq.

Migration matrix for $|\Delta y(Z, t_{lep})|$ variable at parton-level in region SR-3L-WZ.

Migration matrix for $H_{\text{T}}^{\text{ l}}$ variable at particle-level in region SR-4L-DF.

Migration matrix for $H_{\text{T}}^{\text{ l}}$ variable at particle-level in region SR-4L-SF.

Migration matrix for $H_{\text{T}}^{\text{ l}}$ variable at particle-level in region CR-4L-ZZ.

Migration matrix for $H_{\text{T}}^{\text{ l}}$ variable at parton-level in region SR-4L-DF.

Migration matrix for $H_{\text{T}}^{\text{ l}}$ variable at parton-level in region SR-4L-SF.

Migration matrix for $H_{\text{T}}^{\text{ l}}$ variable at parton-level in region CR-4L-ZZ.

Migration matrix for $H_{\text{T}}^{\text{ l}}$ variable at particle-level in region SR-3L-ttZ.

Migration matrix for $H_{\text{T}}^{\text{ l}}$ variable at particle-level in region SR-3L-tZq.

Migration matrix for $H_{\text{T}}^{\text{ l}}$ variable at particle-level in region SR-3L-WZ.

Migration matrix for $H_{\text{T}}^{\text{ l}}$ variable at parton-level in region SR-3L-ttZ.

Migration matrix for $H_{\text{T}}^{\text{ l}}$ variable at parton-level in region SR-3L-tZq.

Migration matrix for $H_{\text{T}}^{\text{ l}}$ variable at parton-level in region SR-3L-WZ.

Migration matrix for $m^{t\bar{t}Z}$ variable at particle-level in region SR-3L-ttZ.

Migration matrix for $m^{t\bar{t}Z}$ variable at particle-level in region SR-3L-tZq.

Migration matrix for $m^{t\bar{t}Z}$ variable at particle-level in region SR-3L-WZ.

Migration matrix for $m^{t\bar{t}Z}$ variable at particle-level in region SR-4L-DF.

Migration matrix for $m^{t\bar{t}Z}$ variable at particle-level in region SR-4L-SF.

Migration matrix for $m^{t\bar{t}Z}$ variable at particle-level in region CR-4L-ZZ.

Migration matrix for $m^{t\bar{t}Z}$ variable at parton-level in region SR-3L-ttZ.

Migration matrix for $m^{t\bar{t}Z}$ variable at parton-level in region SR-3L-tZq.

Migration matrix for $m^{t\bar{t}Z}$ variable at parton-level in region SR-3L-WZ.

Migration matrix for $m^{t\bar{t}Z}$ variable at parton-level in region SR-4L-DF.

Migration matrix for $m^{t\bar{t}Z}$ variable at parton-level in region SR-4L-SF.

Migration matrix for $m^{t\bar{t}Z}$ variable at parton-level in region CR-4L-ZZ.

Migration matrix for $m^{t\bar{t}}$ variable at particle-level in region SR-3L-ttZ.

Migration matrix for $m^{t\bar{t}}$ variable at particle-level in region SR-3L-tZq.

Migration matrix for $m^{t\bar{t}}$ variable at particle-level in region SR-3L-WZ.

Migration matrix for $m^{t\bar{t}}$ variable at particle-level in region SR-4L-DF.

Migration matrix for $m^{t\bar{t}}$ variable at particle-level in region SR-4L-SF.

Migration matrix for $m^{t\bar{t}}$ variable at particle-level in region CR-4L-ZZ.

Migration matrix for $m^{t\bar{t}}$ variable at parton-level in region SR-3L-ttZ.

Migration matrix for $m^{t\bar{t}}$ variable at parton-level in region SR-3L-tZq.

Migration matrix for $m^{t\bar{t}}$ variable at parton-level in region SR-3L-WZ.

Migration matrix for $m^{t\bar{t}}$ variable at parton-level in region SR-4L-DF.

Migration matrix for $m^{t\bar{t}}$ variable at parton-level in region SR-4L-SF.

Migration matrix for $m^{t\bar{t}}$ variable at parton-level in region CR-4L-ZZ.

Migration matrix for $N_{\text{jets}}$ variable at particle-level in region SR-4L-DF.

Migration matrix for $N_{\text{jets}}$ variable at particle-level in region SR-4L-SF.

Migration matrix for $N_{\text{jets}}$ variable at particle-level in region CR-4L-ZZ.

Migration matrix for $N_{\text{jets}}$ variable at particle-level in region SR-3L-ttZ.

Migration matrix for $N_{\text{jets}}$ variable at particle-level in region SR-3L-tZq.

Migration matrix for $N_{\text{jets}}$ variable at particle-level in region SR-3L-WZ.

Migration matrix for $p^{Z}_{T}$ variable at particle-level in region SR-3L-ttZ.

Migration matrix for $p^{Z}_{T}$ variable at particle-level in region SR-3L-tZq.

Migration matrix for $p^{Z}_{T}$ variable at particle-level in region SR-3L-WZ.

Migration matrix for $p^{Z}_{T}$ variable at particle-level in region SR-4L-DF.

Migration matrix for $p^{Z}_{T}$ variable at particle-level in region SR-4L-SF.

Migration matrix for $p^{Z}_{T}$ variable at particle-level in region CR-4L-ZZ.

Migration matrix for $p^{Z}_{T}$ variable at parton-level in region SR-3L-ttZ.

Migration matrix for $p^{Z}_{T}$ variable at parton-level in region SR-3L-tZq.

Migration matrix for $p^{Z}_{T}$ variable at parton-level in region SR-3L-WZ.

Migration matrix for $p^{Z}_{T}$ variable at parton-level in region SR-4L-DF.

Migration matrix for $p^{Z}_{T}$ variable at parton-level in region SR-4L-SF.

Migration matrix for $p^{Z}_{T}$ variable at parton-level in region CR-4L-ZZ.

Migration matrix for $p_{T}^{\mathrm{top}}$ variable at particle-level in region SR-3L-ttZ.

Migration matrix for $p_{T}^{\mathrm{top}}$ variable at particle-level in region SR-3L-tZq.

Migration matrix for $p_{T}^{\mathrm{top}}$ variable at particle-level in region SR-3L-WZ.

Migration matrix for $p_{T}^{\mathrm{top}}$ variable at particle-level in region SR-4L-DF.

Migration matrix for $p_{T}^{\mathrm{top}}$ variable at particle-level in region SR-4L-SF.

Migration matrix for $p_{T}^{\mathrm{top}}$ variable at particle-level in region CR-4L-ZZ.

Migration matrix for $p_{T}^{\mathrm{top}}$ variable at parton-level in region SR-3L-ttZ.

Migration matrix for $p_{T}^{\mathrm{top}}$ variable at parton-level in region SR-3L-tZq.

Migration matrix for $p_{T}^{\mathrm{top}}$ variable at parton-level in region SR-3L-WZ.

Migration matrix for $p_{T}^{\mathrm{top}}$ variable at parton-level in region SR-4L-DF.

Migration matrix for $p_{T}^{\mathrm{top}}$ variable at parton-level in region SR-4L-SF.

Migration matrix for $p_{T}^{\mathrm{top}}$ variable at parton-level in region CR-4L-ZZ.

Migration matrix for $p_{T}^{t\bar{t}}$ variable at particle-level in region SR-3L-ttZ.

Migration matrix for $p_{T}^{t\bar{t}}$ variable at particle-level in region SR-3L-tZq.

Migration matrix for $p_{T}^{t\bar{t}}$ variable at particle-level in region SR-3L-WZ.

Migration matrix for $p_{T}^{t\bar{t}}$ variable at particle-level in region SR-4L-DF.

Migration matrix for $p_{T}^{t\bar{t}}$ variable at particle-level in region SR-4L-SF.

Migration matrix for $p_{T}^{t\bar{t}}$ variable at particle-level in region CR-4L-ZZ.

Migration matrix for $p_{T}^{t\bar{t}}$ variable at parton-level in region SR-3L-ttZ.

Migration matrix for $p_{T}^{t\bar{t}}$ variable at parton-level in region SR-3L-tZq.

Migration matrix for $p_{T}^{t\bar{t}}$ variable at parton-level in region SR-3L-WZ.

Migration matrix for $p_{T}^{t\bar{t}}$ variable at parton-level in region SR-4L-DF.

Migration matrix for $p_{T}^{t\bar{t}}$ variable at parton-level in region SR-4L-SF.

Migration matrix for $p_{T}^{t\bar{t}}$ variable at parton-level in region CR-4L-ZZ.

Migration matrix for $p_{\text{T}}^{\ell, non-Z}$ variable at particle-level in region SR-3L-ttZ.

Migration matrix for $p_{\text{T}}^{\ell, non-Z}$ variable at particle-level in region SR-3L-tZq.

Migration matrix for $p_{\text{T}}^{\ell, non-Z}$ variable at particle-level in region SR-3L-WZ.

Migration matrix for $p_{\text{T}}^{\ell, non-Z}$ variable at parton-level in region SR-3L-ttZ.

Migration matrix for $p_{\text{T}}^{\ell, non-Z}$ variable at parton-level in region SR-3L-tZq.

Migration matrix for $p_{\text{T}}^{\ell, non-Z}$ variable at parton-level in region SR-3L-WZ.

Migration matrix for $|y^{Z}$| variable at particle-level in region SR-3L-ttZ.

Migration matrix for $|y^{Z}$| variable at particle-level in region SR-3L-tZq.

Migration matrix for $|y^{Z}$| variable at particle-level in region SR-3L-WZ.

Migration matrix for $|y^{Z}$| variable at particle-level in region SR-4L-DF.

Migration matrix for $|y^{Z}$| variable at particle-level in region SR-4L-SF.

Migration matrix for $|y^{Z}$| variable at particle-level in region CR-4L-ZZ.

Migration matrix for $|y^{Z}$| variable at parton-level in region SR-3L-ttZ.

Migration matrix for $|y^{Z}$| variable at parton-level in region SR-3L-tZq.

Migration matrix for $|y^{Z}$| variable at parton-level in region SR-3L-WZ.

Migration matrix for $|y^{Z}$| variable at parton-level in region SR-4L-DF.

Migration matrix for $|y^{Z}$| variable at parton-level in region SR-4L-SF.

Migration matrix for $|y^{Z}$| variable at parton-level in region CR-4L-ZZ.

Migration matrix for $|y^{t\bar{t}Z}|$ variable at particle-level in region SR-3L-ttZ.

Migration matrix for $|y^{t\bar{t}Z}|$ variable at particle-level in region SR-3L-tZq.

Migration matrix for $|y^{t\bar{t}Z}|$ variable at particle-level in region SR-3L-WZ.

Migration matrix for $|y^{t\bar{t}Z}|$ variable at particle-level in region SR-4L-DF.

Migration matrix for $|y^{t\bar{t}Z}|$ variable at particle-level in region SR-4L-SF.

Migration matrix for $|y^{t\bar{t}Z}|$ variable at particle-level in region CR-4L-ZZ.

Migration matrix for $|y^{t\bar{t}Z}|$ variable at parton-level in region SR-3L-ttZ.

Migration matrix for $|y^{t\bar{t}Z}|$ variable at parton-level in region SR-3L-tZq.

Migration matrix for $|y^{t\bar{t}Z}|$ variable at parton-level in region SR-3L-WZ.

Migration matrix for $|y^{t\bar{t}Z}|$ variable at parton-level in region SR-4L-DF.

Migration matrix for $|y^{t\bar{t}Z}|$ variable at parton-level in region SR-4L-SF.

Migration matrix for $|y^{t\bar{t}Z}|$ variable at parton-level in region CR-4L-ZZ.

Covariance matrix for absolute cross section as a function of $p_{T}^{\mathrm{top}}$ at particle-level.

Covariance matrix for normalized cross section as a function of $p_{T}^{\mathrm{top}}$ at particle-level.

Covariance matrix for absolute cross section as a function of $p_{T}^{\mathrm{top}}$ at parton-level.

Covariance matrix for normalized cross section as a function of $p_{T}^{\mathrm{top}}$ at parton-level.

Covariance matrix for absolute cross section as a function of $p_{T}^{t\bar{t}}$ at particle-level.

Covariance matrix for normalized cross section as a function of $p_{T}^{t\bar{t}}$ at particle-level.

Covariance matrix for absolute cross section as a function of $p_{T}^{t\bar{t}}$ at parton-level.

Covariance matrix for normalized cross section as a function of $p_{T}^{t\bar{t}}$ at parton-level.

Covariance matrix for absolute cross section as a function of $|\Delta\Phi(t\bar{t}, Z)|/\pi$ at particle-level.

Covariance matrix for normalized cross section as a function of $|\Delta\Phi(t\bar{t}, Z)|/\pi$ at particle-level.

Covariance matrix for absolute cross section as a function of $|\Delta\Phi(t\bar{t}, Z)|/\pi$ at parton-level.

Covariance matrix for normalized cross section as a function of $|\Delta\Phi(t\bar{t}, Z)|/\pi$ at parton-level.

Covariance matrix for absolute cross section as a function of $m^{t\bar{t}Z}$ at particle-level.

Covariance matrix for normalized cross section as a function of $m^{t\bar{t}Z}$ at particle-level.

Covariance matrix for absolute cross section as a function of $m^{t\bar{t}Z}$ at parton-level.

Covariance matrix for normalized cross section as a function of $m^{t\bar{t}Z}$ at parton-level.

Covariance matrix for absolute cross section as a function of $m^{t\bar{t}}$ at particle-level.

Covariance matrix for normalized cross section as a function of $m^{t\bar{t}}$ at particle-level.

Covariance matrix for absolute cross section as a function of $m^{t\bar{t}}$ at parton-level.

Covariance matrix for normalized cross section as a function of $m^{t\bar{t}}$ at parton-level.

Covariance matrix for absolute cross section as a function of $|y^{t\bar{t}Z}|$ at particle-level.

Covariance matrix for normalized cross section as a function of $|y^{t\bar{t}Z}|$ at particle-level.

Covariance matrix for absolute cross section as a function of $|y^{t\bar{t}Z}|$ at parton-level.

Covariance matrix for normalized cross section as a function of $|y^{t\bar{t}Z}|$ at parton-level.

Covariance matrix for absolute cross section as a function of cos $\theta_{Z}^{*}$ at particle-level.

Covariance matrix for normalized cross section as a function of cos $\theta_{Z}^{*}$ at particle-level.

Covariance matrix for absolute cross section as a function of cos $\theta_{Z}^{*}$ at parton-level.

Covariance matrix for normalized cross section as a function of cos $\theta_{Z}^{*}$ at parton-level.

Covariance matrix for absolute cross section as a function of $|\Delta\Phi(l_{t}^{+}, l_{\bar{t}}^{-})|/\pi$ at particle-level.

Covariance matrix for normalized cross section as a function of $|\Delta\Phi(l_{t}^{+}, l_{\bar{t}}^{-})|/\pi$ at particle-level.

Covariance matrix for absolute cross section as a function of $|\Delta\Phi(l_{t}^{+}, l_{\bar{t}}^{-})|/\pi$ at parton-level.

Covariance matrix for normalized cross section as a function of $|\Delta\Phi(l_{t}^{+}, l_{\bar{t}}^{-})|/\pi$ at parton-level.

Covariance matrix for absolute cross section as a function of $|\Delta\Phi(Z, t_{lep})|/\pi$ at particle-level.

Covariance matrix for normalized cross section as a function of $|\Delta\Phi(Z, t_{lep})|/\pi$ at particle-level.

Covariance matrix for absolute cross section as a function of $|\Delta\Phi(Z, t_{lep})|/\pi$ at parton-level.

Covariance matrix for normalized cross section as a function of $|\Delta\Phi(Z, t_{lep})|/\pi$ at parton-level.

Covariance matrix for absolute cross section as a function of $|\Delta y(Z, t_{lep})|$ at particle-level.

Covariance matrix for normalized cross section as a function of $|\Delta y(Z, t_{lep})|$ at particle-level.

Covariance matrix for absolute cross section as a function of $|\Delta y(Z, t_{lep})|$ at parton-level.

Covariance matrix for normalized cross section as a function of $|\Delta y(Z, t_{lep})|$ at parton-level.

Covariance matrix for absolute cross section as a function of $H_{\text{T}}^{\text{l}}$ at in the tetralepton channel particle-level.

Covariance matrix for normalized cross section as a function of $H_{\text{T}}^{\text{l}}$ at in the tetralepton channel particle-level.

Covariance matrix for absolute cross section as a function of $H_{\text{T}}^{\text{l}}$ at in the tetralepton channel parton-level.

Covariance matrix for normalized cross section as a function of $H_{\text{T}}^{\text{l}}$ in the tetralepton channel at parton-level.

Covariance matrix for absolute cross section as a function of $H_{\text{T}}^{\text{l}}$ in the trilepton channel at particle-level.

Covariance matrix for normalized cross section as a function of $H_{\text{T}}^{\text{l}}$ in the trilepton channel at particle-level.

Covariance matrix for absolute cross section as a function of $H_{\text{T}}^{\text{l}}$ in the trilepton channel at parton-level.

Covariance matrix for normalized cross section as a function of $H_{\text{T}}^{\text{l}}$ in the trilepton channel at parton-level.

Covariance matrix for absolute cross section as a function of $N_{\text{jets}}$ in the tetralepton channel at particle-level.

Covariance matrix for normalized cross section as a function of $N_{\text{jets}}$ in the tetralepton channel at particle-level.

Covariance matrix for absolute cross section as a function of $N_{\text{jets}}$ in the trilepton channel at particle-level.

Covariance matrix for normalized cross section as a function of $N_{\text{jets}}$ in the trilepton channel at particle-level.

Covariance matrix for absolute cross section as a function of $p^{Z}_{T}$ at particle-level.

Covariance matrix for normalized cross section as a function of $p^{Z}_{T}$ at particle-level.

Covariance matrix for absolute cross section as a function of $p^{Z}_{T}$ at parton-level.

Covariance matrix for normalized cross section as a function of $p^{Z}_{T}$ at parton-level.

Covariance matrix for absolute cross section as a function of $p_{\text{T}}^{\ell, non-Z}$ at particle-level.

Covariance matrix for normalized cross section as a function of $p_{\text{T}}^{\ell, non-Z}$ at particle-level.

Covariance matrix for absolute cross section as a function of $p_{\text{T}}^{\ell, non-Z}$ at parton-level.

Covariance matrix for normalized cross section as a function of $p_{\text{T}}^{\ell, non-Z}$ at parton-level.

Covariance matrix for absolute cross section as a function of $|y^{Z}$| at particle-level.

Covariance matrix for normalized cross section as a function of $|y^{Z}$| at particle-level.

Covariance matrix for absolute cross section as a function of $|y^{Z}$| at parton-level.

Covariance matrix for normalized cross section as a function of $|y^{Z}$| at parton-level.

Ranking of nuisance parameters and background normalizations on signal strength for inclusive cross section measurement in combination of all channels

Observed and expected 68% and 95% credible intervals for the top-boson operators, in the marginalised linear fit.

Observed and expected 68% and 95% credible intervals for the top-boson operators, in the marginalised quadratic fit.

Observed and expected 68% and 95% credible intervals for the top-boson operators, in the independent quadratic fits (allowing only one Wilson Coefficient to be non-zero).

Observed and expected 68% and 95% credible intervals for the four-quark operators, in the marginalised linear fit.

Observed and expected 68% and 95% credible intervals for the four-quark operators, in the marginalised quadratic fit.

Observed and expected 68% and 95% credible intervals for the four-quark operators, in the independent quadratic fits (allowing only one Wilson Coefficient to be non-zero).

Observed and expected 68% and 95% credible intervals for Fisher-rotated directions of EFT sensitivity, in the marginalised linear fit.

Correlation matrix of the input particle-level observables used in the EFT fit.

Table with post-fit yields in the four regions used in the measurement. The correlations of the nuisance parameters, as obtained by the fit, are considered for the calculation of the uncertainties.

Ratio of the $t\bar{t}$ to the $Z$-boson cross section compared to the predictions for several sets of parton distribution functions. The total cross section for $t\bar{t}$ production and the fiducial cross section for $Z$-boson production are used. The uncertainty on the prediction represents the PDF+$\alpha_{s}$ variations.

Nuisance parameter ranking plot for the $t\bar{t}$ signal-strength. Only the 10 most important nuisance parameters are shown. The impact of each nuisance parameter on the measured cross section, is defined as the shift induced in measured cross-section as the nuisance parameter is shifted by its pre-fit and post-fit uncertainties from its nominal best-fit value and then fixed while all other nuisance parameters profiled.

Nuisance parameter ranking plot for the ratio of the $t\bar{t}$ and $Z$-boson signal-strengths. Only the 10 most important nuisance parameters are shown. The impact of each nuisance parameter on the measured ratio, is defined as the shift induced in measured ratio as the nuisance parameter is shifted by its pre-fit and post-fit uncertainties from its nominal best-fit value and then fixed while all other nuisance parameters profiled.

Nuisance parameter ranking plot for the $Z$-boson signal-strength. Only the 10 most important nuisance parameters are shown. The impact of each nuisance parameter on the measured cross section, is defined as the shift induced in measured cross-section as the nuisance parameter is shifted by its pre-fit and post-fit uncertainties from its nominal best-fit value and then fixed while all other nuisance parameters profiled.

Fiducial differential cross section of the electroweak $W^\pm W^\pm jj$ production as a function of $N_{\mathrm{gap}\,\mathrm{jets}}$. The correlation of uncertainties of the measured cross section across bins is presented in Table 14.

Fiducial differential cross section of the electroweak $W^\pm W^\pm jj$ production as a function of $\xi_{\mathrm{j}3}$. The correlation of uncertainties of the measured cross section across bins is presented in Table 15.

Fiducial differential cross section of the inclusive $W^\pm W^\pm jj$ production as a function of $m_{\ell\ell}$. The correlation of uncertainties of the measured cross section across bins is presented in Table 16.

Fiducial differential cross section of the inclusive $W^\pm W^\pm jj$ production as a function of $m_{\mathrm{T}}$. The correlation of uncertainties of the measured cross section across bins is presented in Table 17.

Fiducial differential cross section of the inclusive $W^\pm W^\pm jj$ production as a function of $m_{\mathrm{jj}}$. The correlation of uncertainties of the measured cross section across bins is presented in Table 18.

Fiducial differential cross section of the inclusive $W^\pm W^\pm jj$ production as a function of $N_{\mathrm{gap}\,\mathrm{jets}}$. The correlation of uncertainties of the measured cross section across bins is presented in Table 19.

Fiducial differential cross section of the inclusive $W^\pm W^\pm jj$ production as a function of $\xi_{\mathrm{j}3}$. The correlation of uncertainties of the measured cross section across bins is presented in Table 20.

Observed correlations between the bins of the LH-unfolded cross section of the electroweak $W^\pm W^\pm jj$ production as a function of $m_{\ell\ell}$. See Table 1.

Observed correlations between the bins of the LH-unfolded cross section of the electroweak $W^\pm W^\pm jj$ production as a function of $m_{\mathrm{T}}$. See Table 2.

Observed correlations between the bins of the LH-unfolded cross section of the electroweak $W^\pm W^\pm jj$ production as a function of $m_{\mathrm{jj}}$. See Table 3.

Observed correlations between the bins of the LH-unfolded cross section of the electroweak $W^\pm W^\pm jj$ production as a function of $N_{\mathrm{gap}\,\mathrm{jets}}$. See Table 4.

Observed correlations between the bins of the LH-unfolded cross section of the electroweak $W^\pm W^\pm jj$ production as a function of $\xi_{\mathrm{j}3}$. See Table 5.

Observed correlations between the bins of the LH-unfolded cross section of the inclusive $W^\pm W^\pm jj$ production as a function of $m_{\ell\ell}$. See Table 6.

Observed correlations between the bins of the LH-unfolded cross section of the inclusive $W^\pm W^\pm jj$ production as a function of $m_{\mathrm{T}}$. See Table 7.

Observed correlations between the bins of the LH-unfolded cross section of the inclusive $W^\pm W^\pm jj$ production as a function of $m_{\mathrm{jj}}$. See Table 8.

Observed correlations between the bins of the LH-unfolded cross section of the inclusive $W^\pm W^\pm jj$ production as a function of $N_{\mathrm{gap}\,\mathrm{jets}}$. See Table 9.

Observed correlations between the bins of the LH-unfolded cross section of the inclusive $W^\pm W^\pm jj$ production as a function of $\xi_{\mathrm{j}3}$. See Table 10.

Evolution of the one-dimensional expected and observed limits at 95% CL on the parameters corresponding to the quartic operators with label M0 as a function of the cut-off scale. The unitarity bounds as a function of the cut-off scale are defined for one non-zero Wilson coefficient.

Evolution of the one-dimensional expected and observed limits at 95% CL on the parameters corresponding to the quartic operators with label M1 as a function of the cut-off scale. The unitarity bounds as a function of the cut-off scale are defined for one non-zero Wilson coefficient.

Evolution of the one-dimensional expected and observed limits at 95% CL on the parameters corresponding to the quartic operators with label M7 as a function of the cut-off scale. The unitarity bounds as a function of the cut-off scale are defined for one non-zero Wilson coefficient. The limits on M7 were obtained without taking into account the SM-EFT interference for the EW W$^\pm$Zjj final state.

Evolution of the one-dimensional expected and observed limits at 95% CL on the parameters corresponding to the quartic operators with label S02 as a function of the cut-off scale. The unitarity bounds as a function of the cut-off scale are defined for one non-zero Wilson coefficient.

Evolution of the one-dimensional expected and observed limits at 95% CL on the parameters corresponding to the quartic operators with label S1 as a function of the cut-off scale. The unitarity bounds as a function of the cut-off scale are defined for one non-zero Wilson coefficient.

Evolution of the one-dimensional expected and observed limits at 95% CL on the parameters corresponding to the quartic operators with label T0 as a function of the cut-off scale. The unitarity bounds as a function of the cut-off scale are defined for one non-zero Wilson coefficient.

Evolution of the one-dimensional expected and observed limits at 95% CL on the parameters corresponding to the quartic operators with label T1 as a function of the cut-off scale. The unitarity bounds as a function of the cut-off scale are defined for one non-zero Wilson coefficient.

Evolution of the one-dimensional expected and observed limits at 95% CL on the parameters corresponding to the quartic operators with label T2 as a function of the cut-off scale. The unitarity bounds as a function of the cut-off scale are defined for one non-zero Wilson coefficient.

Expected and observed exclusion limits at 95% CL for $\sigma_{\mathrm{VBF}}(\mathrm{H}^{\pm\pm}_5) \times \mathcal{B}(\mathrm{H}^{\pm\pm}_5 \rightarrow W^{\pm}W^{\pm})$ as a function of the doubly-charged Higgs mass.

Expected and observed exclusion limits at 95% CL for $\sin\theta_{\mathrm{H}}$ as a function of the doubly-charged Higgs mass in the Georgi-Machacek model.

The mass planes of the reconstructed Higgs boson candidates for the 2Pass selections of the analysis, shown for the VBF $\kappa_{2V} = 0$ HH samples.

The mass planes of the reconstructed Higgs boson candidates for the 2Pass selections of the analysis, shown for the $m_{X} = 1$ TeV spin-0 narrow-width resonance HH samples.

Observed values of $-2\ln\Lambda$ as a function of $\kappa_{2V}$ for the resolved (dotted green) and boosted (dashed blue) analyses, and their combination (solid black), with all other coupling modifiers fixed to their SM predictions.

Expected values of $-2\ln\Lambda$ as a function of $\kappa_{2V}$ for the resolved (dotted green) and boosted (dashed blue) analyses, and their combination (solid black), with all other coupling modifiers fixed to their SM predictions.

Observed $-2\ln\Lambda$ values of boosted-only result in the $\kappa_{\Lambda}-\kappa_{2V}$ plane.

Observed $-2\ln\Lambda$ values of resolved-only result in the $\kappa_{\Lambda}-\kappa_{2V}$ plane.

Observed $-2\ln\Lambda$ values of combination result in the $\kappa_{\Lambda}-\kappa_{2V}$ plane.

Expected $-2\ln\Lambda$ values of boosted-only result in the $\kappa_{\Lambda}-\kappa_{2V}$ plane.

Expected $-2\ln\Lambda$ values of resolved-only result in the $\kappa_{\Lambda}-\kappa_{2V}$ plane.

Expected $-2\ln\Lambda$ values of combination result in the $\kappa_{\Lambda}-\kappa_{2V}$ plane.

Expected (dashed black lines) and observed (solid black lines) 95% CL upper limits on the cross-section of spin-0 heavy resonances with narrow width assumptions. The SM H → bb branching ratio is assumed in both cases. The ±1$\sigma$ and ±2$\sigma$ uncertainty ranges for the expected limits are shown as coloured bands.

Expected (dashed black lines) and observed (solid black lines) 95% CL upper limits on the cross-section of spin-0 heavy resonances with broad width assumptions. The SM H → bb branching ratio is assumed in both cases. The ±1$\sigma$ and ±2$\sigma$ uncertainty ranges for the expected limits are shown as coloured bands.

Observed and expected 95% CL exclusion limits on the production cross-sections of the combined VBF HH process as a function of $\kappa_{2V}$ . The expected exclusion limits assume no HH production and are obtained using a fit to the data with the background-only hypothesis. The red line shows the theory prediction for the VBF HH cross-section as a function of $\kappa_{2V}$ where all other parameters and couplings are set to their SM values.

Observed and expected 95% CL exclusion limits on the production cross-sections of the combined HH process as a function of $\kappa_{\lambda}$ . The expected exclusion limits assume no HH production and are obtained using a fit to the data with the background-only hypothesis. The red line shows the theory prediction for the HH cross-section as a function of $\kappa_{\lambda}$ where all other parameters and couplings are set to their SM values.

Observed −2 ln Λ as a function of $\kappa_{\lambda}$ with the $\kappa_{2V}$ and $\kappa_{V}$ modifiers fixed to unity.

Expected −2 ln Λ as a function of $\kappa_{\lambda}$ with the $\kappa_{2V}$ and $\kappa_{V}$ modifiers fixed to unity.

Cumulative acceptance times efficiency as a function of resonance mass for each event selection step for the narrow-width signal

Cumulative acceptance times efficiency as a function of resonance mass for each event selection step for the broad-width signal

Cumulative acceptance times efficiency as a function of $\kappa_{2V}$ for each event selection step for the nonresonant signal.

The distributions of the invariant mass of the HH system for various values of $\kappa_{2V}$ .

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.

Summary of the observed and expected 95% CL upper limits on the cross section for $\mathrm{LQ}_{\mathrm{mix}}^{\mathrm{d}}$ pair production as a function of $m_{\mathrm{LQ}_{\mathrm{mix}}^{\mathrm{d}}}$ under the assumptions of B(LQ$\rightarrow t\mu$)=1.

Summary of observed and predicted yields in the four signal region categories. The background prediction is shown after the combined likelihood fit to data under the background-only hypothesis across all control region and signal region categories. The expected signal yields that are obtained by using their theoretical cross sections are also shown with their pre-fit uncertainties, assuming B=1 and $\mu$=1. The "Other" contribution is dominated by $t\bar{t}t\bar{t}$ and $t\bar{t}WW$ in the $3\ell$ SRs, whereas it is dominated by $tWZ$ and $t\bar{t}WW$ in the $4\ell$ SRs. Dashes refer to components that are negligible or not applicable.

Summary of the data, background and signal yields in the four signal region categories. The background prediction is shown before the combined likelihood fit to data. Dashes (--) refer to components that are negligible or not applicable.

Summary of the data and background yields in the seven control region categories. The background prediction is shown before the combined likelihood fit to data. Dashes (--) refer to components that are negligible or not applicable.

Summary of the data and background yields in the seven control region categories. The background prediction is shown after the combined likelihood fit to data under the background-only hypothesis across all control region categories. Dashes (--) refer to components that are negligible or not applicable.

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.