The expected (dotted black) and observed (solid black) upper limits at 95% CL on $\sigma_{T^\prime bq(T^\prime \rightarrow tH)}$ in the hadronic channel displayed as a function of $\rm{M}_{T^\prime}$. The green (yellow) band represents the 68% (95%) of the limit values expected under the background-only hypothesis. The theoretical cross sections for the singlet $T^\prime$ production with representative $\kappa_T$-values fixed at 0.1, 0.15, 0.2 and 0.25 (for $\Gamma/\rm{M}_{T^\prime} < 5\%$) are shown as solid red lines.

Selection of the hadronic signal region for the three BDTs and $\rm{M}_{T^\prime}$ windows.

Selection of the leptonic signal region for the three BDTs and $\rm{M}_{T^\prime}$ windows.

The expected yields of different processes in each signal window for events with a $T^{\prime}$ with mass in the range $\rm{M}_{T^\prime} \in [600, 1200]$ GeV, and the observed number of events in the signal region $m_{\gamma\gamma} \in [115, 135]$. Here, the yields for the $T^\prime$ are for $\kappa_{T}$ fixed at 0.2.

Signal efficiency for the hadronic channel (left) and the leptonic channel (right) in the signal regions as optimized for each of the three different $\rm{M}_{T^\prime}$ ranges, [600, 700], [700, 1000] and [1000, 1200] GeV, represented as red, green and blue curves respectively. The efficiency is defined as the ratio of the events after the final selection to the total expected events.

Correlation matrix between the fiducial cross sections for the four individual decay final states and the $ZZ^{*}$ normalisation factor.

Differential fiducial cross section for the transverse momentum $p_{T}^{4l}$ of the Higgs boson. The measured cross sections are compared to predictions provided by NNLOPS + XH. NNLOPS is normalised to the N3LO total cross section with a K-factor = 1.1 . Measured value in the last bin is un upper limit at 95% CL.

Correlation matrix between the measured cross sections and the $ZZ^{*}$ background normalization corresponding to the differential fiducial cross section for the transverse momentum $p_{T}^{4l}$ of the Higgs boson.

Differential fiducial cross section for the invariant mass $m_{12}$ of the leading Z boson. The measured cross sections are compared to predictions provided by NNLOPS + XH. NNLOPS is normalised to the N3LO total cross section with a K-factor = 1.1 .

Correlation matrix between the measured cross sections and the $ZZ^{*}$ background normalization corresponding to the differential fiducial cross section for the invariant mass $m_{12}$ of the leading Z boson.

Differential fiducial cross section for the invariant mass $m_{34}$ of the subleading Z boson. The measured cross sections are compared to predictions provided by NNLOPS + XH. NNLOPS is normalised to the N3LO total cross section with a K-factor = 1.1 .

Correlation matrix between the measured cross sections and the $ZZ^{*}$ background normalization corresponding to the differential fiducial cross section for the invariant mass $m_{34}$ of the subleading Z boson.

Differential fiducial cross section for the rapidity $|y_{4l}|$ of the Higgs boson. The measured cross sections are compared to predictions provided by NNLOPS + XH. NNLOPS is normalised to the N3LO total cross section with a K-factor = 1.1 .

Correlation matrix between the measured cross sections and the $ZZ^{*}$ background normalization corresponding to the differential fiducial cross section for the rapidity $|y_{4l}|$ of the Higgs boson.

Differential fiducial cross section for the production angle $|\cos\theta^{*}|$ of the leading Z boson. The measured cross sections are compared to predictions provided by NNLOPS + XH. NNLOPS is normalised to the N3LO total cross section with a K-factor = 1.1 .

Correlation matrix between the measured cross sections and the $ZZ^{*}$ background normalization corresponding to the differential fiducial cross section for the production angle $|\cos\theta^{*}|$ of the leading Z boson.

Differential fiducial cross section for the production angle $\cos\theta_{1}$ of the anti-lepton from the leading Z boson. The measured cross sections are compared to predictions provided by NNLOPS + XH. NNLOPS is normalised to the N3LO total cross section with a K-factor = 1.1 .

Correlation matrix between the measured cross sections and the $ZZ^{*}$ background normalization corresponding to the differential fiducial cross section for the production angle $\cos\theta_{1}$ of the anti-lepton from the leading Z boson.

Differential fiducial cross section for the production angle $\cos\theta_{2}$ of the anti-lepton from the subleading Z boson. The measured cross sections are compared to predictions provided by NNLOPS + XH. NNLOPS is normalised to the N3LO total cross section with a K-factor = 1.1 .

Correlation matrix between the measured cross sections and the $ZZ^{*}$ background normalization corresponding to the differential fiducial cross section for the production angle $\cos\theta_{2}$ of the anti-lepton from the subleading Z boson.

Differential fiducial cross section for the azimuthal angle $\phi$ of the decay planes of the two reconstructed Z bosons. The measured cross sections are compared to predictions provided by NNLOPS + XH. NNLOPS is normalised to the N3LO total cross section with a K-factor = 1.1 .

Correlation matrix between the measured cross sections and the $ZZ^{*}$ background normalization corresponding to the differential fiducial cross section for the azimuthal angle $\phi$ of the decay planes of the two reconstructed Z bosons.

Differential fiducial cross section for the azimuthal angle $\phi_{1}$ of the decay plane of the leading Z boson and the plane formed between its four-momentum and the z-axis. The measured cross sections are compared to predictions provided by NNLOPS + XH. NNLOPS is normalised to the N3LO total cross section with a K-factor = 1.1 .

Correlation matrix between the measured cross sections and the $ZZ^{*}$ background normalization corresponding to the differential fiducial cross section for the azimuthal angle $\phi_{1}$ of the decay plane of the leading Z boson and the plane formed between its four-momentum and the z-axis.

Differential fiducial cross section for the jet multiplicity $N_{jets}$. The measured cross sections are compared to predictions provided by NNLOPS + XH. NNLOPS is normalised to the N3LO total cross section with a K-factor = 1.1 .

Correlation matrix between the measured cross sections and the $ZZ^{*}$ background normalization corresponding to the differential fiducial cross section for the jet multiplicity $N_{jets}$.

Differential fiducial cross section for the inclusive jet multiplicity $N_{jets}$. The measured cross sections are compared to predictions provided by NNLOPS + XH. NNLOPS is normalised to the N3LO total cross section with a K-factor = 1.1 .

Differential fiducial cross section for the number of b-quark initiated jets $N_{b-jets}$. The measured cross sections are compared to predictions provided by NNLOPS + XH. NNLOPS is normalised to the N3LO total cross section with a K-factor = 1.1 .

Correlation matrix between the measured cross sections and the $ZZ^{*}$ background normalization corresponding to the differential fiducial cross section for the number of b-quark initiated jets $N_{b-jets}$.

Differential fiducial cross section for the transverse momentum of the leading jet $p_{T}^{lead.jet}$ in events with at least one jet. The measured cross sections are compared to predictions provided by NNLOPS + XH. NNLOPS is normalised to the N3LO total cross section with a K-factor = 1.1 .

Correlation matrix between the measured cross sections and the $ZZ^{*}$ background normalization corresponding to the differential fiducial cross section for the transverse momentum of the leading jet $p_{T}^{lead.jet}$ in events with at least one jet.

Differential fiducial cross section for the transverse momentum of the subleading jet $p_{T}^{sublead.jet}$ in events with at least two jets. The measured cross sections are compared to predictions provided by NNLOPS + XH. NNLOPS is normalised to the N3LO total cross section with a K-factor = 1.1 .

Correlation matrix between the measured cross sections and the $ZZ^{*}$ background normalization corresponding to the differential fiducial cross section for the transverse momentum of the subleading jet $p_{T}^{sublead.jet}$ in events with at least two jets.

Differential fiducial cross section for the invariant mass of the two highest-pT jets $m_{jj}$ in events with at least two jets. The measured cross sections are compared to predictions provided by NNLOPS + XH. NNLOPS is normalised to the N3LO total cross section with a K-factor = 1.1 .

Correlation matrix between the measured cross sections and the $ZZ^{*}$ background normalization corresponding to the differential fiducial cross section for the invariant mass of the two highest-pT jets $m_{jj}$ in events with at least two jets.

Differential fiducial cross section for the distance between the two highest-pT jets in pseudorapidity $\Delta\eta_{jj}$. The measured cross sections are compared to predictions provided by NNLOPS + XH. NNLOPS is normalised to the N3LO total cross section with a K-factor = 1.1 .

Correlation matrix between the measured cross sections and the $ZZ^{*}$ background normalization corresponding to the differential fiducial cross section for the distance between the two highest-pT jets in pseudorapidity $\Delta\eta_{jj}$.

Differential fiducial cross section for the distance between the two highest-pT jets in $\phi$ $\Delta\phi_{jj}$. The measured cross sections are compared to predictions provided by NNLOPS + XH. NNLOPS is normalised to the N3LO total cross section with a K-factor = 1.1 .

Correlation matrix between the measured cross sections and the $ZZ^{*}$ background normalization corresponding to the differential fiducial cross section for the distance between the two highest-pT jets in $\phi$ $\Delta\phi_{jj}$.

Differential fiducial cross section for the transverse momentum of the four lepton plus jet system, in events with at least one jet $p_{T}^{4lj}$. The measured cross sections are compared to predictions provided by NNLOPS + XH. NNLOPS is normalised to the N3LO total cross section with a K-factor = 1.1 .

Correlation matrix between the measured cross sections and the $ZZ^{*}$ background normalization corresponding to the differential fiducial cross section for the transverse momentum of the four lepton plus jet system, in events with at least one jet $p_{T}^{4lj}$.

Differential fiducial cross section for the transverse momentum of the four lepton plus di-jet system, in events with at least two jets $p_{T}^{4ljj}$. The measured cross sections are compared to predictions provided by NNLOPS + XH. NNLOPS is normalised to the N3LO total cross section with a K-factor = 1.1 . Measured value in the last bin is un upper limit at 95% CL.

Correlation matrix between the measured cross sections and the $ZZ^{*}$ background normalization corresponding to the differential fiducial cross section for the transverse momentum of the four lepton plus di-jet system, in events with at least two jets $p_{T}^{4ljj}$.

Differential fiducial cross section for the invariant mass of the four lepton plus jet system in events with at least one jet $m_{4lj}$. The measured cross sections are compared to predictions provided by NNLOPS + XH. NNLOPS is normalised to the N3LO total cross section with a K-factor = 1.1 .

Correlation matrix between the measured cross sections and the $ZZ^{*}$ background normalization corresponding to the differential fiducial cross section for the invariant mass of the four lepton plus jet system in events with at least one jet $m_{4lj}$.

Differential fiducial cross section for the invariant mass of the four lepton plus di-jet system in events with at least two jets $m_{4ljj}$. The measured cross sections are compared to predictions provided by NNLOPS + XH. NNLOPS is normalised to the N3LO total cross section with a K-factor = 1.1 .

Correlation matrix between the measured cross sections and the $ZZ^{*}$ background normalization corresponding to the differential fiducial cross section for the invariant mass of the four lepton plus di-jet system in events with at least two jets $m_{4ljj}$.

Differential fiducial cross section for the leading vs. subleading Z boson mass $m_{12}$vs.$m_{34}$. The measured cross sections are compared to predictions provided by NNLOPS + XH. NNLOPS is normalised to the N3LO total cross section with a K-factor = 1.1 .

Correlation matrix between the measured cross sections and the $ZZ^{*}$ background normalization corresponding to the differential fiducial cross section for the leading vs. subleading Z boson mass $m_{12}$vs.$m_{34}$.

Differential fiducial cross section for the leading vs. subleading Z boson mass $m_{12}$vs.$m_{34}$ in $ll\mu\mu$ final states. The measured cross sections are compared to predictions provided by NNLOPS + XH. NNLOPS is normalised to the N3LO total cross section with a K-factor = 1.1 .

Differential fiducial cross section for the leading vs. subleading Z boson mass $m_{12}$vs.$m_{34}$ in $llee$ final states. The measured cross sections are compared to predictions provided by NNLOPS + XH. NNLOPS is normalised to the N3LO total cross section with a K-factor = 1.1 .

Correlation matrix between the measured cross sections and the $ZZ^{*}$ background normalization corresponding to the differential fiducial cross section for the leading vs. subleading Z boson mass m12 vs. m34 in $ll\mu\mu$ and $llee$ final states.

Differential fiducial cross section of the $p_{T}^{4l}$ distribution in $|y_{4l}|$ bins. The measured cross sections are compared to predictions provided by NNLOPS + XH. NNLOPS is normalised to the N3LO total cross section with a K-factor = 1.1 .

Correlation matrix between the measured cross sections and the $ZZ^{*}$ background normalization corresponding to the differential fiducial cross section of the $p_{T}^{4l}$ distribution in $|y_{4l}|$ bins.

Differential fiducial cross section of the $p_{T}^{4l}$ distribution in $N_{jets}$ bins. The measured cross sections are compared to predictions provided by NNLOPS + XH. NNLOPS is normalised to the N3LO total cross section with a K-factor = 1.1 .

Correlation matrix between the measured cross sections and the $ZZ^{*}$ background normalization corresponding to the differential fiducial cross section of the $p_{T}^{4l}$ distribution in $N_{jets}$ bins.

Differential fiducial cross section for transverse momentum of the four lepton system vs. the transverse momentum of the four lepton plus jet system $p_{T}^{4l}$vs.$p_{T}^{4lj}$. The measured cross sections are compared to predictions provided by NNLOPS + XH. NNLOPS is normalised to the N3LO total cross section with a K-factor = 1.1 .

Correlation matrix between the measured cross sections and the $ZZ^{*}$ background normalization corresponding to the differential fiducial cross section for transverse momentum of the four lepton system vs. the transverse momentum of the four lepton plus jet system $p_{T}^{4l}$vs.$p_{T}^{4lj}$.

Differential fiducial cross section for the transverse momentum of the four lepton plus jet system vs the invariant mass of the four lepton plus jet system $p_{T}^{4l}$vs.$m_{4lj}$. The measured cross sections are compared to predictions provided by NNLOPS + XH. NNLOPS is normalised to the N3LO total cross section with a K-factor = 1.1 .

Correlation matrix between the measured cross sections and the $ZZ^{*}$ background normalization corresponding to the differential fiducial cross section for the transverse momentum of the four lepton plus jet system vs the invariant mass of the four lepton plus jet system $p_{T}^{4l}$vs.$m_{4lj}$.

Differential fiducial cross section for the transverse momentum of the four lepton vs the transverse momentum of the leading jet $p_{T}^{4l}$vs.$p_{T}^{l.jet}$. The measured cross sections are compared to predictions provided by NNLOPS + XH. NNLOPS is normalised to the N3LO total cross section with a K-factor = 1.1 .

Correlation matrix between the measured cross sections and the $ZZ^{*}$ background normalization corresponding to the differential fiducial cross section for the transverse momentum of the four lepton vs the transverse momentum of the leading jet $p_{T}^{4l}$vs.$p_{T}^{lead.jet}$.

Differential fiducial cross section for the transverse momentum of the leading jet vs the rapidity of the leading jet $p_{T}^{lead.jet}$vs.$|y^{lead.jet}|$. The measured cross sections are compared to predictions provided by NNLOPS + XH. NNLOPS is normalised to the N3LO total cross section with a K-factor = 1.1 .

Correlation matrix between the measured cross sections and the $ZZ^{*}$ background normalization corresponding to the differential fiducial cross section for the transverse momentum of the leading jet vs the rapidity of the leading jet $p_{T}^{lead.jet}$vs.$|y^{lead.jet}|$.

Differential fiducial cross section for the transverse momentum of the leading jet vs the transverse momentum of the subleading jet $p_{T}^{lead.jet}$vs.$p_{T}^{sublead.jet}$. The measured cross sections are compared to predictions provided by NNLOPS + XH. NNLOPS is normalised to the N3LO total cross section with a K-factor = 1.1 .

Correlation matrix between the measured cross sections and the $ZZ^{*}$ background normalization corresponding to the differential fiducial cross section for the transverse momentum of the leading jet vs the transverse momentum of the subleading jet $p_{T}^{lead.jet}$vs.$p_{T}^{sublead.jet}$.

Differential fiducial cross section for the leading Z boson mass $m_{12}$ in $4\mu$ and $4e$ final states. The measured cross sections are compared to predictions provided by NNLOPS + XH. NNLOPS is normalised to the N3LO total cross section with a K-factor = 1.1 .

Differential fiducial cross section for the leading Z boson mass $m_{12}$ in $2e2\mu$ and $2\mu2e$ final states. The measured cross sections are compared to predictions provided by NNLOPS + XH. NNLOPS is normalised to the N3LO total cross section with a K-factor = 1.1 .

Correlation matrix between the measured cross sections and the $ZZ^{*}$ background normalization corresponding to the differential fiducial cross section for the leading Z boson mass $m_{12}$ in $4l$ and $2l2l$ final states.

Differential fiducial cross section for the subleading Z boson mass $m_{34}$ in $4\mu$ and $4e$ final states. The measured cross sections are compared to predictions provided by NNLOPS + XH. NNLOPS is normalised to the N3LO total cross section with a K-factor = 1.1 .

Differential fiducial cross section for the subleading Z boson mass $m_{34}$ in $2e2\mu$ and $2\mu2e$ final states. The measured cross sections are compared to predictions provided by NNLOPS + XH. NNLOPS is normalised to the N3LO total cross section with a K-factor = 1.1 .

Correlation matrix between the measured cross sections and the $ZZ^{*}$ background normalization corresponding to the differential fiducial cross section for the subleading Z boson mass $m_{34}$ in $4l$ and $2l2l$ final states.

Differential fiducial cross section for the azimuthal angle $\phi$ of the decay planes of the two reconstructed Z bosons in $4\mu$ and $4e$ final states. The measured cross sections are compared to predictions provided by NNLOPS + XH. NNLOPS is normalised to the N3LO total cross section with a K-factor = 1.1 .

Differential fiducial cross section for the azimuthal angle $\phi$ of the decay planes of the two reconstructed Z bosons in $2e2\mu$ and $2\mu2e$ final states. The measured cross sections are compared to predictions provided by NNLOPS + XH. NNLOPS is normalised to the N3LO total cross section with a K-factor = 1.1 .

Correlation matrix between the measured cross sections and the $ZZ^{*}$ background normalization corresponding to the differential fiducial cross section for the azimuthal angle $\phi$ of the decay planes of the two reconstructed Z bosons in $4l$ and $2l2l$ final states.

Differential fiducial cross section for the leading vs. subleading Z boson mass $m_{12}$vs.$m_{34}$ in $4\mu$ and $4e$ final states. The measured cross sections are compared to predictions provided by NNLOPS + XH. NNLOPS is normalised to the N3LO total cross section with a K-factor = 1.1 .

Differential fiducial cross section for the leading vs. subleading Z boson mass $m_{12}$vs.$m_{34}$ in $2\mu2e$ and $2e2\mu$ final states. The measured cross sections are compared to predictions provided by NNLOPS + XH. NNLOPS is normalised to the N3LO total cross section with a K-factor = 1.1 .

Correlation matrix between the measured cross sections and the $ZZ^{*}$ background normalization corresponding to the differential fiducial cross section for the leading vs. subleading Z boson mass $m_{12}$vs.$m_{34}$ in $4l$ and $2l2l$ final states.

The expected (pre-fit) and observed numbers of events for an integrated luminosity of 139 fb$^{-1}$ at $\sqrt{s}$=13 TeV in the signal region 115 < $m_{4l}$< 130 GeV and sideband region $m_{4l}$ in 105-115 GeV or 130-160 GeV (350 GeV for $tXX$-enriched) in each reconstructed event category assuming the SM Higgs boson signal with a mass $m_{H}$= 125 GeV. Combined statistical and systematic uncertainties are included for the predictions. Expected contributions that are below 0.2% of the total yield in each reconstructed event category are not shown and replaced by -.

The expected and the observed (post-fit) the jet multiplicity distribution after the inclusive event selection for an integrated luminosity of 139fb$^{-1}$ at $\sqrt{s}$=13TeV. The SM Higgs boson signal is assumed to have a mass $m_{H}$=125 GeV.

The expected and the observed (post-fit) the four-lepton transverse momentum distribution for events with zero jets, after the inclusive event selection for an integrated luminosity of 139fb$^{-1}$ at $\sqrt{s}$=13TeV. The SM Higgs boson signal is assumed to have a mass $m_{H}$=125 GeV.

The expected and the observed (post-fit) the four-lepton transverse momentum distribution for events with one jet, after the inclusive event selection for an integrated luminosity of 139fb$^{-1}$ at $\sqrt{s}$=13TeV. The SM Higgs boson signal is assumed to have a mass $m_{H}$=125 GeV.

The expected and the observed (post-fit) the four-lepton transverse momentum distribution for events with at least two jets, after the inclusive event selection for an integrated luminosity of 139fb$^{-1}$ at $\sqrt{s}$=13TeV. The SM Higgs boson signal is assumed to have a mass $m_{H}$=125 GeV.

The expected and the observed (post-fit) the dijet invariant mass distribution for events with at least two jets, after the inclusive event selection for an integrated luminosity of 139fb$^{-1}$ at $\sqrt{s}$=13TeV. The SM Higgs boson signal is assumed to have a mass $m_{H}$=125 GeV.

The expected and the observed $NN_{ggF}$ output (post-fit) distribution in 0$j$-$p_{T}^{4l}$-Low category for an integrated luminosity of 139fb$^{-1}$ at $\sqrt{s}$=13TeV. The SM Higgs boson signal is assumed to have a mass $m_{H}$=125 GeV.

The expected and the observed $NN_{ggF}$ output (post-fit) distribution in 0$j$-$p_{T}^{4l}$-Med category for an integrated luminosity of 139fb$^{-1}$ at $\sqrt{s}$=13TeV. The SM Higgs boson signal is assumed to have a mass $m_{H}$=125 GeV.

The expected and the observed $NN_{VBF}$ output (post-fit) distribution in 1$j$-$p_{T}^{4l}$-Low category with $NN_{ZZ}$<0.25 for an integrated luminosity of 139fb$^{-1}$ at $\sqrt{s}$=13TeV. The SM Higgs boson signal is assumed to have a mass $m_{H}$=125 GeV.

The expected and the observed $NN_{ZZ}$ output (post-fit) distribution in 1$j$-$p_{T}^{4l}$-Low category with $NN_{ZZ}$>0.25 for an integrated luminosity of 139fb$^{-1}$ at $\sqrt{s}$=13TeV. The SM Higgs boson signal is assumed to have a mass $m_{H}$=125 GeV.

The expected and the observed $NN_{VBF}$ output (post-fit) distribution in 1$j$-$p_{T}^{4l}$-Med category with $NN_{ZZ}$<0.25 for an integrated luminosity of 139fb$^{-1}$ at $\sqrt{s}$=13TeV. The SM Higgs boson signal is assumed to have a mass $m_{H}$=125 GeV.

The expected and the observed $NN_{ZZ}$ output (post-fit) distribution in 1$j$-$p_{T}^{4l}$-Med category with $NN_{ZZ}$>0.25 for an integrated luminosity of 139fb$^{-1}$ at $\sqrt{s}$=13TeV. The SM Higgs boson signal is assumed to have a mass $m_{H}$=125 GeV.

The expected and the observed $NN_{VBF}$ output (post-fit) distribution in 1$j$-$p_{T}^{4l}$-High category for an integrated luminosity of 139fb$^{-1}$ at $\sqrt{s}$=13TeV. The SM Higgs boson signal is assumed to have a mass $m_{H}$=125 GeV.

The expected and the observed $NN_{VBF}$ output (post-fit) distribution in 2$j$ category with $NN_{VH}$<0.2 for an integrated luminosity of 139fb$^{-1}$ at $\sqrt{s}$=13TeV. The SM Higgs boson signal is assumed to have a mass $m_{H}$=125 GeV.

The expected and the observed $NN_{VH}$ output (post-fit) distribution in 2$j$ category with $NN_{VH}$>0.2 for an integrated luminosity of 139fb$^{-1}$ at $\sqrt{s}$=13TeV. The SM Higgs boson signal is assumed to have a mass $m_{H}$=125 GeV.

The expected and the observed $NN_{VBF}$ output (post-fit) distribution in 2$j$-BSM-like category for an integrated luminosity of 139fb$^{-1}$ at $\sqrt{s}$=13TeV. The SM Higgs boson signal is assumed to have a mass $m_{H}$=125 GeV.

The expected and the observed $NN_{ttH}$ output (post-fit) distribution in $ttH$-Had-enriched category with $NN_{tXX}$<0.4 for an integrated luminosity of 139fb$^{-1}$ at $\sqrt{s}$=13TeV. The SM Higgs boson signal is assumed to have a mass $m_{H}$=125 GeV.

The expected and the observed $NN_{tXX}$ output (post-fit) distribution in $ttH$-Had-enriched category with $NN_{tXX}$>0.4 for an integrated luminosity of 139fb$^{-1}$ at $\sqrt{s}$=13TeV. The SM Higgs boson signal is assumed to have a mass $m_{H}$=125 GeV.

The expected and the observed $NN_{ttH}$ output (post-fit) distribution in $VH$-Lep-enriched category for an integrated luminosity of 139fb$^{-1}$ at $\sqrt{s}$=13TeV. The SM Higgs boson signal is assumed to have a mass $m_{H}$=125 GeV.

The expected and the observed events in the categories where no NN discriminant is used for an integrated luminosity of 139fb$^{-1}$ at $\sqrt{s}$=13TeV. The SM Higgs boson signal is assumed to have a mass $m_{H}$=125 GeV.

The expected and the observed events in the side bands used to constraint the $ZZ^{*}$ and $tXX$ backgrounds for an integrated luminosity of 139fb$^{-1}$ at $\sqrt{s}$=13TeV. The SM Higgs boson signal is assumed to have a mass $m_{H}$=125 GeV.

The expected SM cross section $(\sigma \cdot BR)_{SM}$, the observed value of $(\sigma \cdot BR)$, and their ratio $(\sigma \cdot BR)/(\sigma \cdot BR)_{SM}$ for the inclusive production and for each Production Mode and Reduced Stage-1.1 production bin for an integrated luminosity of 139fb$^{-1}$ at $\sqrt{s}$=13 TeV.

The correlation matrices between the measured cross-sections and the $ZZ$ and $tXX$ normalisation factors for Production Mode Stage.

The correlation matrices between the measured cross-sections and the $ZZ$ and $tXX$ normalisation factors for Reduced Stage 1.1.

Likelihood contours at 68% CLs in the (ggF,VBF) plane. The observed best fit value is the first value of the table.

Likelihood contours at 95% CLs in the (ggF,VBF) plane. The observed best fit value is the first value of the table.

Likelihood contours at 68% CLs in the (ggF,VH) plane. The observed best fit value is the first value of the table.

Likelihood contours at 95% CLs in the (ggF,VH) plane. The observed best fit value is the first value of the table.

Likelihood contours at 68% CLs in the (VBF,VH) plane. The observed best fit value is the first value of the table.

Likelihood contours at 95% CLs in the (VBF,VH) plane. The observed best fit value is the first value of the table.

Likelihood contours at 68% CLs in the ($gg2H-0j-p_T^H-Low$,$gg2H-0j-p_T^H-High$) plane. The observed best fit value is the first value of the table.

Likelihood contours at 95% CLs in the ($gg2H-0j-p_T^H-Low$,$gg2H-0j-p_T^H-High$) plane. The observed best fit value is the first value of the table.

Likelihood contours at 68% CLs in the ($\kappa_{V}$vs.$\kappa_{F}$) plane. The observed best fit value is the first value of the table.

Likelihood contours at 95% CLs in the ($\kappa_{V}$vs.$\kappa_{F}$) plane. The observed best fit value is the first value of the table.

The expected signal yield ratio for chosen CP-even and CP-odd EFT parameter values together with the corresponding cross-section measurement in each production bin of Reduced Stage 1.1. The parameter values correspond approximately to the expected confidence intervals at the 68% CLs obtained from the statistical interpretation of data.

The expected and observed confidence intervals at 68% and 95% CLs of the SMEFT Wilson coefficients for an integrated luminosity of 139fb$^{-1}$ at $\sqrt{s}$=13 TeV. For the coupling with two best fit values, the two values have been splitted in the two different columns 'Two Best-fit value - first' and 'Two Best-fit value - second', and the symbol '-' has been inserted in the missing fields.

Expected 2D-fit likelihood curves at the 95% CLs for the $c_{HW}$ versus $c_{HB}$ coupling parameters at an integrated luminosity of 139 fb$^{-1}$ and $\sqrt{s}$=13 TeV. Except for the two fitted Wilson coefficients, all others are set to zero.

Observed 2D-fit likelihood curves at the 95% CLs for the $c_{HW}$ versus $c_{HB}$ coupling parameters at an integrated luminosity of 139 fb$^{-1}$ and $\sqrt{s}$=13 TeV. Except for the two fitted Wilson coefficients, all others are set to zero.

Expected 2D-fit likelihood curves at the 95% CLs for the $c_{HB}$ versus $c_{HG}$ coupling parameters at an integrated luminosity of 139 fb$^{-1}$ and $\sqrt{s}$=13 TeV. Except for the two fitted Wilson coefficients, all others are set to zero.

Observed 2D-fit likelihood curves at the 95% CLs for the $c_{HB}$ versus $c_{HG}$ coupling parameters at an integrated luminosity of 139 fb$^{-1}$ and $\sqrt{s}$=13 TeV. Except for the two fitted Wilson coefficients, all others are set to zero.

Expected 2D-fit likelihood curves at the 95% CLs for the $c_{uH}$ versus $c_{HG}$ coupling parameters at an integrated luminosity of 139 fb$^{-1}$ and $\sqrt{s}$=13 TeV. Except for the two fitted Wilson coefficients, all others are set to zero.

Observed 2D-fit likelihood curves at the 95% CLs for the $c_{uH}$ versus $c_{HG}$ coupling parameters at an integrated luminosity of 139 fb$^{-1}$ and $\sqrt{s}$=13 TeV. Except for the two fitted Wilson coefficients, all others are set to zero.

The best-fit values and the corresponding deviation from the SM prediction obtained from two-dimensional likelihood scans of the CP-odd BSM coupling parameters with 139fb$^{-1}$ data at $\sqrt{s}$=13 TeV. For the coupling with two best fit values, the two values have been splitted in the two different columns 'Two Observed best-fit first/second_POI - first' and 'Two Observed best-fit first/second_POI - second', and the symbol '-' has been inserted in the missing fields.

Expected 2D-fit likelihood curves at the 95% CLs for the $c_{H\tilde{W}}$ versus $c_{H\tilde{B}}$ coupling parameters at an integrated luminosity of 139 fb$^{-1}$ and $\sqrt{s}$=13 TeV. Except for the two fitted Wilson coefficients, all others are set to zero.

Observed 2D-fit likelihood curves at the 95% CLs for the $c_{H\tilde{W}}$ versus $c_{H\tilde{B}}$ coupling parameters at an integrated luminosity of 139 fb$^{-1}$ and $\sqrt{s}$=13 TeV. Except for the two fitted Wilson coefficients, all others are set to zero.

Expected 2D-fit likelihood curves at the 95% CLs for the $c_{H\tilde{B}}$ versus $c_{H\tilde{G}}$ coupling parameters at an integrated luminosity of 139 fb$^{-1}$ and $\sqrt{s}$=13 TeV. Except for the two fitted Wilson coefficients, all others are set to zero.

Observed 2D-fit likelihood curves at the 95% CLs for the $c_{H\tilde{B}}$ versus $c_{H\tilde{G}}$ coupling parameters at an integrated luminosity of 139 fb$^{-1}$ and $\sqrt{s}$=13 TeV. Except for the two fitted Wilson coefficients, all others are set to zero.

Expected 2D-fit likelihood curves at the 95% CLs for the $c_{\tilde{u}H}$ versus $c_{H\tilde{G}}$ coupling parameters at an integrated luminosity of 139 fb$^{-1}$ and $\sqrt{s}$=13 TeV. Except for the two fitted Wilson coefficients, all others are set to zero.

Observed 2D-fit likelihood curves at the 95% CLs for the $c_{\tilde{u}H}$ versus $c_{H\tilde{G}}$ coupling parameters at an integrated luminosity of 139 fb$^{-1}$ and $\sqrt{s}$=13 TeV. Except for the two fitted Wilson coefficients, all others are set to zero.

Signal acceptance times efficiency, expressed as a percentage, obtained from the ratio of the number of simulated signal events satisfying the event selection criteria in each reconstructed category to the total number of events generated in the phase space specified by a given Production model bin.

Signal acceptance times efficiency, expressed as a percentage, obtained from the ratio of the number of simulated signal events satisfying the event selection criteria in each reconstructed category to the total number of events generated in the phase space specified by a given Reduced Stage-1.1 production bin.

Observed covariance matrix for the different parameters of interest when fitting the for Production Mode Stage.

Observed covariance matrix for the different parameters of interest when fitting the for Reduced Stage 1.1.

The expected SM cross section $(\sigma\cdot BR)_{SM}$, the observed value of $(\sigma\cdot BR)$, and their ratio $(\sigma\cdot BR)/(\sigma\cdot BR)_{SM}$ for the alternate Reduced Stage-1.1 scheme for an integrated luminosity of 139fb$^{-1}$ at $\sqrt{s}$=13 TeV.

The correlation matrices between the measured cross-sections and the ZZ and $tXX$ normalisation factors for alternative Reduced Stage 1.1.

Expected 2D-fit likelihood curves at the 95% CLs for the $c_{HW}$ versus $c_{HWB}$ coupling parameters at an integrated luminosity of 139 fb$^{-1}$ and $\sqrt{s}$=13 TeV. Except for the two fitted Wilson coefficients, all others are set to zero.

Observed 2D-fit likelihood curves at the 95% CLs for the $c_{HW}$ versus $c_{HWB}$ coupling parameters at an integrated luminosity of 139 fb$^{-1}$ and $\sqrt{s}$=13 TeV. Except for the two fitted Wilson coefficients, all others are set to zero.

Expected 2D-fit likelihood curves at the 95% CLs for the $c_{HB}$ versus $c_{HWB}$ coupling parameters at an integrated luminosity of 139 fb$^{-1}$ and $\sqrt{s}$=13 TeV. Except for the two fitted Wilson coefficients, all others are set to zero.

Observed 2D-fit likelihood curves at the 95% CLs for the $c_{HB}$ versus $c_{HWB}$ coupling parameters at an integrated luminosity of 139 fb$^{-1}$ and $\sqrt{s}$=13 TeV. Except for the two fitted Wilson coefficients, all others are set to zero.

Expected 2D-fit likelihood curves at the 95% CLs for the $c_{HW}$ versus $c_{HG}$ coupling parameters at an integrated luminosity of 139 fb$^{-1}$ and $\sqrt{s}$=13 TeV. Except for the two fitted Wilson coefficients, all others are set to zero.

Observed 2D-fit likelihood curves at the 95% CLs for the $c_{HW}$ versus $c_{HG}$ coupling parameters at an integrated luminosity of 139 fb$^{-1}$ and $\sqrt{s}$=13 TeV. Except for the two fitted Wilson coefficients, all others are set to zero.

Expected 2D-fit likelihood curves at the 95% CLs for the $c_{HWB}$ versus $c_{HG}$ coupling parameters at an integrated luminosity of 139 fb$^{-1}$ and $\sqrt{s}$=13 TeV. Except for the two fitted Wilson coefficients, all others are set to zero.

Observed 2D-fit likelihood curves at the 95% CLs for the $c_{HWB}$ versus $c_{HG}$ coupling parameters at an integrated luminosity of 139 fb$^{-1}$ and $\sqrt{s}$=13 TeV. Except for the two fitted Wilson coefficients, all others are set to zero.

Expected 2D-fit likelihood curves at the 95% CLs for the $c_{H\tilde{W}}$ versus $c_{H\tilde{W}B}$ coupling parameters at an integrated luminosity of 139 fb$^{-1}$ and $\sqrt{s}$=13 TeV. Except for the two fitted Wilson coefficients, all others are set to zero.

Observed 2D-fit likelihood curves at the 95% CLs for the $c_{H\tilde{W}}$ versus $c_{H\tilde{W}B}$ coupling parameters at an integrated luminosity of 139 fb$^{-1}$ and $\sqrt{s}$=13 TeV. Except for the two fitted Wilson coefficients, all others are set to zero.

Expected 2D-fit likelihood curves at the 95% CLs for the $c_{H\tilde{B}}$ versus $c_{H\tilde{W}B}$ coupling parameters at an integrated luminosity of 139 fb$^{-1}$ and $\sqrt{s}$=13 TeV. Except for the two fitted Wilson coefficients, all others are set to zero.

Observed 2D-fit likelihood curves at the 95% CLs for the $c_{H\tilde{B}}$ versus $c_{H\tilde{W}B}$ coupling parameters at an integrated luminosity of 139 fb$^{-1}$ and $\sqrt{s}$=13 TeV. Except for the two fitted Wilson coefficients, all others are set to zero.

Expected 2D-fit likelihood curves at the 95% CLs for the $c_{H\tilde{W}}$ versus $c_{H\tilde{G}}$ coupling parameters at an integrated luminosity of 139 fb$^{-1}$ and $\sqrt{s}$=13 TeV. Except for the two fitted Wilson coefficients, all others are set to zero.

Observed 2D-fit likelihood curves at the 95% CLs for the $c_{H\tilde{W}}$ versus $c_{H\tilde{G}}$ coupling parameters at an integrated luminosity of 139 fb$^{-1}$ and $\sqrt{s}$=13 TeV. Except for the two fitted Wilson coefficients, all others are set to zero.

Expected 2D-fit likelihood curves at the 95% CLs for the $c_{H\tilde{W}B}$ versus $c_{H\tilde{G}}$ coupling parameters at an integrated luminosity of 139 fb$^{-1}$ and $\sqrt{s}$=13 TeV. Except for the two fitted Wilson coefficients, all others are set to zero.

Observed 2D-fit likelihood curves at the 95% CLs for the $c_{H\tilde{W}B}$ versus $c_{H\tilde{G}}$ coupling parameters at an integrated luminosity of 139 fb$^{-1}$ and $\sqrt{s}$=13 TeV. Except for the two fitted Wilson coefficients, all others are set to zero.

The observed values of Sigma*BR(H->ZZ*), the SM expected cross section sBRsm and their ratio Sigma*BR/(Sigma*BR)_SM for the inclusive production and in each Stage-0 and reduced Stage-1 production bin for an integrated luminosity of 36.1/fb and at sqrt(s)=13 TeV. The bbH contribution is considered as a part of the ggF production bins. The upper limits correspond to the 95% CL obtained with pseudo-experiments using the CL_s method. The uncertainties are given as (stat.)+(exp.)+(th.) for Stage 0 and as (stat.)+(syst.) for reduced Stage 1. Values without uncertainity are 95% CL upper limits.

Signal acceptance obtained as the ratio of the number of simulated signal events satisfying the event selection criteria in each reconstructed event category over the total number of events generated in the phase space specified by a given reduced Stage-1 ggF production bin. Results are obtained in bins of BDT discriminants using coarse binning with several bins merged into one. Values less than 0.0001 are set to 0.

Signal acceptance obtained as the ratio of the number of simulated signal events satisfying the event selection criteria in each reconstructed event category over the total number of events generated in the phase space specified by the given reduced Stage-1 VBF and VH production bins. Results are obtained in bins of BDT discriminants using coarse binning with several bins merged into one. Values less than 0.0001 are set to 0.

The signal strengths mu for the inclusive production and in each Stage-0 and reduced Stage-1 production bin for an integrated luminosity of 36.1/fb and at sqrt(s)=13 TeV. The bbH contribution is considered as a part of the ggF production bins. The upper limits correspond to the 95% CL obtained with pseudo-experiments using the CL_s method. The uncertainties are given as (stat.)+(exp.)+(th.) for Stage 0 and as (stat.)+(syst.) for reduced Stage 1. Values without uncertainity are 95% CL upper limits.

Signal acceptance (in percent) obtained as the ratio of the number of simulated signal events satisfying the event selection criteria in each reconstructed event category to the total number of generated events, as predicted by the MadGraph5_aMC@NLO generator assuming the SM coupling tensor structure or the BSM tensor structure with ($\kappa_{SM}$ = 1, | $\kappa_{AVV}$ | $\neq$ 0).

Number of expected ggF Higgs boson events for an integrated luminosity of $\mathcal L=36.1 \text{fb}^{-1}$ and at $\sqrt{\mathrm{s}}=13$ TeV, as predicted by the MadGraph5_aMC@NLO generator assuming the SM coupling tensor structure or the BSM tensor structure with ($\kappa_{SM}=1$, $|\kappa_{Avv}|=6$). The highest-order SM predicition for the sum of the ggF, ttH and bbH contributions is also shown for comparison.

Number of expected VBF and VH Higgs boson events for an integrated luminosity of $\mathcal L=36.1 \text{fb}^{-1}$ and at $\sqrt{\mathrm{s}}=13$ TeV, as predicted by the MadGraph5_aMC@NLO generator assuming the SM coupling tensor structure or the BSM tensor structure with ($\kappa_{SM}=1$, $|\kappa_{Avv}|=5$). The highest-order SM predicition for the sum of the VBF and VH contributions is also shown for comparison.

Expected Correlation Matrix for Stage 0

Observed Correlation Matrix for Stage 0. As upper limits are derived for ttH and VH POIs using the observed data, the corresponding terms inside the matrix are set to zero.

Expected Correlation Matrix for Reduced Stage 1

Observed Correlation Matrix for Reduced Stage 1. As upper limits are derived for ttH and VH POIs using the observed data, the corresponding terms inside the matrix are set to zero.

Expected Covariance Matrix for Stage 0

Observed Covariance Matrix for Stage 0. As upper limits are derived for ttH and VH POIs using the observed data, the corresponding terms inside the matrix are set to zero.

Expected Covariance Matrix for Reduced Stage 1

Observed Covariance Matrix for Reduced Stage 1. As upper limits are derived for ttH and VH POIs using the observed data, the corresponding terms inside the matrix are set to zero.

Likelihood contours at 68% CL in the (Sigma_ggF*B , Sigma_VBF*B ) plane

Likelihood contours at 95% CL in the (Sigma_ggF*B , Sigma_VBF*B ) plane

Expected two-dimensional negative log-likelihood scans for $\kappa_{HVV}$ versus $\kappa_{AVV}$ coupling parameters using $\mathcal L=36.1 \text{fb}^{-1}$ of data and at $\sqrt{\mathrm{s}}=13$ TeV. The couplings $\kappa_{Hgg}$ and $\kappa_{SM}$ are fixed to the SM value of one in the fit. The 95% CL exclusion limits are shown.

Observed two-dimensional negative log-likelihood scans for $\kappa_{HVV}$ versus $\kappa_{AVV}$ coupling parameters using $\mathcal L=36.1 \text{fb}^{-1}$ of data and at $\sqrt{\mathrm{s}}=13$ TeV. The couplings $\kappa_{Hgg}$ and $\kappa_{SM}$ are fixed to the SM value of one in the fit. The 95% CL exclusion limits are shown.

Expected two-dimensional negative log-likelihood scans for $\kappa_{HVV}$ versus $\kappa_{AVV}$ coupling parameters using $\mathcal L=36.1 \text{fb}^{-1}$ of data and at $\sqrt{\mathrm{s}}=13$ TeV. The coupling $\kappa_{Hgg}$ is fixed to the SM value of one in the fit. The coupling $\kappa_{SM}$ is left as a free parameter of the fit. The 95% CL exclusion limits are shown.

Observed two-dimensional negative log-likelihood scans for $\kappa_{HVV}$ versus $\kappa_{AVV}$ coupling parameters using $\mathcal L=36.1 \text{fb}^{-1}$ of data and at $\sqrt{\mathrm{s}}=13$ TeV. The coupling $\kappa_{Hgg}$ is fixed to the SM value of one in the fit. The coupling $\kappa_{SM}$ is left as a free parameter of the fit. The 95% CL exclusion limits are shown.

Expected two-dimensional negative log-likelihood scans for $\kappa_{HVV}$ versus $\kappa_{SM}$ coupling parameters using $\mathcal L=36.1 \text{fb}^{-1}$ of data and at $\sqrt{\mathrm{s}}=13$ TeV. The 95% CL exclusion limits are shown.

Observed two-dimensional negative log-likelihood scans for $\kappa_{HVV}$ versus $\kappa_{SM}$ coupling parameters using $\mathcal L=36.1 \text{fb}^{-1}$ of data and at $\sqrt{\mathrm{s}}=13$ TeV. The 95% CL exclusion limits are shown.

Expected two-dimensional negative log-likelihood scans for $\kappa_{AVV}$ versus $\kappa_{SM}$ coupling parameters using $\mathcal L=36.1 \text{fb}^{-1}$ of data and at $\sqrt{\mathrm{s}}=13$ TeV. The 95% CL exclusion limits are shown.

Observed two-dimensional negative log-likelihood scans for $\kappa_{AVV}$ versus $\kappa_{SM}$ coupling parameters using $\mathcal L=36.1 \text{fb}^{-1}$ of data and at $\sqrt{\mathrm{s}}=13$ TeV. The 95% CL exclusion limits are shown.