Observed and expected 95% credibility-level upper limits on the cross-section times branching ratio for the simplified dark-matter model. The table also shows the corresponding $1\sigma$ and $2\sigma$ bands for the expected limits. The limits are calculated using jets in events with at least one isolated lepton ($e$ or $\mu$) with $p_\text{T}^\ell \ge 60$ GeV.

The 95% CL upper limits on the cross-section times acceptance times efficiency times branching ratio, as a function of the resonance mass, $m_X$. The limits are calculated using jets in events with at least one isolated lepton ($e$ or $\mu$) with $p_\text{T}^\ell \ge 60$ GeV. The limits are shown for a generic Gaussian signal with the width $\sigma_{X}/m_{X} =0$. The table also shows the expected limits and the corresponding $1\sigma$ and $2\sigma$ bands. The tabulated data has been used in Figure 04 of the publication.

The 95% CL upper limits on the cross-section times acceptance times efficiency times branching ratio, as a function of the resonance mass, $m_X$. The limits are calculated using jets in events with at least one isolated lepton ($e$ or $\mu$) with $p_\text{T}^\ell \ge 60$ GeV. The limits are shown for a generic Gaussian signal with the width $\sigma_{X}/m_{X} =0.05$. The table also shows the expected limits and the corresponding $1\sigma$ and $2\sigma$ bands. The tabulated data has been used in Figure 04 of the publication.

The 95% CL upper limits on the cross-section times acceptance times efficiency times branching ratio, as a function of the resonance mass, $m_X$. The limits are calculated using jets in events with at least one isolated lepton ($e$ or $\mu$) with $p_\text{T}^\ell \ge 60$ GeV. The limits are shown for a generic Gaussian signal with the width $\sigma_{X}/m_{X} =0.10$. The table also shows the expected limits and the corresponding $1\sigma$ and $2\sigma$ bands. The tabulated data has been used in Figure 04 of the publication.

The 95% CL upper limits on the cross-section times acceptance times efficiency times branching ratio, as a function of the resonance mass, $m_X$. The limits are calculated using jets in events with at least one isolated lepton ($e$ or $\mu$) with $p_\text{T}^\ell \ge 60$ GeV. The limits are shown for a generic Gaussian signal with the width $\sigma_{X}/m_{X} =0.15$. The table also shows the expected limits and the corresponding $1\sigma$ and $2\sigma$ bands. The tabulated data has been used in Figure 04 of the publication.

Truth-level cross-sections ($\sigma$), acceptances ($A$), detector efficiencies ($\epsilon$) and the detector-level visible cross-sections ($\sigma^{vis}$) for different masses of the $\rho_T$ model. Only statistical uncertainties on the acceptances and efficiencies are shown. The acceptance is defined by requiring jets with $m_{jj} >216$ GeV and $|\eta|<2.4$, as well as isolated leptons ($e$ or $\mu$) with $p_\text{T}^\ell \ge 60$ GeV and $|\eta^\mu|<2.5$ ($|\eta^e|<2.47$, excluding of the gap region 1.37--1.52) for muons (electrons). The tabulated data has been used in Figure 05a of the publication.

Truth-level cross-sections ($\sigma$), acceptances ($A$), detector efficiencies ($\epsilon$) and the detector-level visible cross-sections ($\sigma^{vis}$) for different masses of the $W'/Z'$ SSM model. Only statistical uncertainties on the acceptances and efficiencies are shown. The acceptance is defined by requiring jets with $m_{jj} >216$ GeV and $|\eta|<2.4$, as well as leptons with $p_\text{T}^\ell \ge 60$ GeV and $|\eta^\mu|<2.5$ ($|\eta^e|<2.47$ with exclusion of the gap region 1.37-1.52) for muons (electrons). The tabulated data has been used in Figure 05b of the publication.

Truth-level cross-sections ($\sigma$), acceptances ($A$), detector efficiencies ($\epsilon$) and the detector-level visible cross-sections ($\sigma^{vis}$) for different masses of the $H^+$ model for $\tan \beta=1$. Only statistical uncertainties on the acceptances and efficiencies are shown. The acceptance correction was calculated assuming $\tan \beta=1$ and the narrow-width approximation. The acceptance is defined by requiring jets with $m_{jj} >216$ GeV and $|\eta|<2.4$, as well as leptons with $p_\text{T}^\ell \ge 60$ GeV and $|\eta^\mu|<2.5$ ($|\eta^e|<2.47$ with exclusion of the gap region 1.37-1.52) for muons (electrons). The tabulated data has been used in Figure 05c of the publication.

Truth-level cross-sections ($\sigma$), acceptances ($A$), detector efficiencies ($\epsilon$) and the detector-level visible cross-sections ($\sigma^{vis}$) for different masses of the $Z'$ (DM) model. Only statistical uncertainties on the acceptances and efficiencies are shown. The acceptance is defined by requiring jets with $m_{jj} > 216$ GeV and $|\eta|<2.4$, as well as leptons with $p_\text{T}^\ell \ge 60$ GeV and $|\eta^\mu|<2.5$ ($|\eta^e|<2.47$ with exclusion of the gap region 1.37-1.52) for muons (electrons). The tabulated data has been used in Figure 05d of the publication.

Lower limits at 95$\%$ CL on $\Lambda$ for the dimuon channel for different signal chiralities in the (constructive/destructive interference) SRs of the analysis.

Lower limits at 95$\%$ CL on $\Lambda$ for the dilepton combination for different signal chiralities in the (constructive/destructive interference) SRs of the analysis.

The dielectron data events passing the full selection.

The dimuon data events passing the full selection.

Distribution of the resonance mass in the c-tag Signal Region of the $ ce$ channel for the post-fit background, the observed data, and the expected signal with $m_{LQ} = 1$ TeV.

Distribution of the resonance mass in the b-tag Signal Region of the $ ce$ channel for the post-fit background, the observed data, and the expected signal with $m_{LQ} = 1$ TeV.

Distribution of the resonance mass in the untagged Signal Region of the $ c\mu$ channel for the post-fit background, the observed data, and the expected signal with $m_{LQ} = 1$ TeV.

Distribution of the resonance mass in the c-tag Signal Region of the $ c\mu$ channel for the post-fit background, the observed data, and the expected signal with $m_{LQ} = 1$ TeV.

Distribution of the resonance mass in the b-tag Signal Region of the $ c\mu$ channel for the post-fit background, the observed data, and the expected signal with $m_{LQ} = 1$ TeV.

Distribution of the resonance mass in the 0-tag Signal Region of the $ be$ channel for the post-fit background, the observed data, and the expected signal with $m_{LQ} = 1$ TeV.

Distribution of the resonance mass in the 1-tag Signal Region of the $ be$ channel for the post-fit background, the observed data, and the expected signal with $m_{LQ} = 1$ TeV.

Distribution of the resonance mass in the 2-tag Signal Region of the $ be$ channel for the post-fit background, the observed data, and the expected signal with $m_{LQ} = 1$ TeV.

Distribution of the resonance mass in the 0-tag Signal Region of the $ b\mu$ channel for the post-fit background, the observed data, and the expected signal with $m_{LQ} = 1$ TeV.

Distribution of the resonance mass in the 1-tag Signal Region of the $ b\mu$ channel for the post-fit background, the observed data, and the expected signal with $m_{LQ} = 1$ TeV.

Distribution of the resonance mass in the 2-tag Signal Region of the $ b\mu$ channel for the post-fit background, the observed data, and the expected signal with $m_{LQ} = 1$ TeV.

The observed and expected limits on the leptoquark pair production cross-section at 95% CL for $\mathcal{B}=1$ into electrons, shown as a function of $m_{LQ}$ for the $qe$ channel.

The observed and expected limits on the leptoquark pair production cross-section at 95% CL for $\mathcal{B}=1$ into muons, shown as a function of $m_{LQ}$ for the $q\mu$ channel.

The observed and expected limits on the leptoquark pair production cross-section at 95% CL for $\mathcal{B}=1$ into electrons, shown as a function of $m_{LQ}$ for the $ce$ channel.

The observed and expected limits on the leptoquark pair production cross-section at 95% CL for $\mathcal{B}=1$ into muons, shown as a function of $m_{LQ}$ for the $c\mu$ channel.

The observed and expected limits on the leptoquark pair production cross-section at 95% CL for $\mathcal{B}=1$ into electrons, shown as a function of $m_{LQ}$ for the $be$ channel.

The observed and expected limits on the leptoquark pair production cross-section at 95% CL for $\mathcal{B}=1$ into muons, shown as a function of $m_{LQ}$ for the $b\mu$ channel.

The observed and expected limits on the leptoquark branching ratio at 95% CL, shown as a function of $m_{LQ}$ for the $qe$ channel.

The observed and expected limits on the leptoquark branching ratio at 95% CL, shown as a function of $m_{LQ}$ for the $q\mu$ channel.

The observed and expected limits on the leptoquark branching ratio at 95% CL, shown as a function of $m_{LQ}$ for the $ce$ channel.

The observed and expected limits on the leptoquark branching ratio at 95% CL, shown as a function of $m_{LQ}$ for the $c\mu$ channel.

The observed and expected limits on the leptoquark branching ratio at 95% CL, shown as a function of $m_{LQ}$ for the $be$ channel.

The observed and expected limits on the leptoquark branching ratio at 95% CL, shown as a function of $m_{LQ}$ for the $b\mu$ channel.

The signal selection efficiency x acceptance summed over all signal regions, for all masses and LQ decay channels considered.

The observed and expected limits for all masses and LQ decay channels considered.

Cutflow Table in the electron channel, considering signal samples with LQ mass of 1 TeV.

Cutflow Table in the muon channel, considering signal samples with LQ mass of 1 TeV.

Summary of results for charm and bottom muon v2 as a function of pT. Uncertainties are statistical and systematic, respectively.

Ratio of unfolded width of azimuthal angular correlation distributions (P PB/ P P). Different colors correspond to different combinations of p_{T,1} and p_{T,2}

Ratio of unfolded Dijet conditional yields (P PB/ P P). Different colors correspond to different combinations of p_{T,1} and p_{T,2}

Unfolded width of azimuthal angular correlation distributions (Delta p_{T} > 3). Full markers represent p+Pb, open markers p+p

Unfolded Dijet conditional yields (Delta p_{T} > 3). Full markers represent p+Pb, open markers p+p

Ratio of unfolded width of azimuthal angular correlation distributions (P PB/ P P) (Delta p_{T} > 3). Different colors correspond to different combinations of p_{T,1} and p_{T,2}

Ratio of unfolded Dijet conditional yields (P PB/ P P) (Delta p_{T} > 3). Different colors correspond to different combinations of p_{T,1} and p_{T,2}

Unfolded azimuthal angular correlation distributions. Black markers represent p+Pb, red markers p+p

Unfolded azimuthal angular correlation distributions (Delta p_{T} > 3). Black markers represent p+Pb, red markers p+p

Observed $m_{t\bar{t}}^{reco}$ distributions in data for SR$2b$, shown together with the result of the fit with the three-shape-parameter function. The error bars indicate the effect of the fit parameter uncertainty on the background prediction. The bin width of the distributions is chosen to be the same as that used in the background parameterization.

Observed and expected upper limits on the cross-section times branching fraction of the $Z'_{TC2} \rightarrow t\bar{t}$ as a function of the $Z'_{TC2}$ mass ($\Gamma/m = 1.2$%).

Fit parameter values obtained from fits to data in the two signal regions. The fit function is $f(x) = p_{0} (1-x)^{p_{1}} x^{p_{2} + p_{3} \log(x)}$, where $x = m^{reco}_{t\bar{t}}/\sqrt{s}$.

Post-fit BDT distributions after the VBF event selection for the $\tau_{\mathrm{had}}\tau_{\mathrm{had}}$ analysis channel. The VBF signal is shown for $\mu = 0.73$ and $\tilde d = -0.01$. ''Other bkg'' denotes all background contributions not listed explicitly in the legend. The size of the combined statistical, experimental, and theoretical uncertainties is given.

Post-fit $m_{\tau\tau}^{\mathrm{MMC}}$ distributions in the low BDT score CR for the $\tau_{\mathrm{lep}}\tau_{\mathrm{lep}}$ SF analysis channel. ''Other bkg'' denotes all background contributions not listed explicitly in the legend. The size of the combined statistical, experimental, and theoretical uncertainties is given. The exact value of the $p_{\mathrm{T}}$ cut on the leptons depends on the trigger.

Post-fit $m_{\tau\tau}^{\mathrm{MMC}}$ distributions in the low BDT score CR for the $\tau_{\mathrm{lep}}\tau_{\mathrm{lep}}$ DF analysis channel. ''Other bkg'' denotes all background contributions not listed explicitly in the legend. The size of the combined statistical, experimental, and theoretical uncertainties is given.

Post-fit $m_{\tau\tau}^{\mathrm{MMC}}$ distributions in the low BDT score CR for the $\tau_{\mathrm{lep}}\tau_{\mathrm{had}}$ analysis channel. ''Other bkg'' denotes all background contributions not listed explicitly in the legend. The size of the combined statistical, experimental, and theoretical uncertainties is given. The exact value of the $p_{\mathrm{T}}$ cut on the leading lepton depends on the trigger.

Post-fit $m_{\tau\tau}^{\mathrm{MMC}}$ distributions in the low BDT score CR for the $\tau_{\mathrm{had}}\tau_{\mathrm{had}}$ analysis channel. ''Other bkg'' denotes all background contributions not listed explicitly in the legend. The size of the combined statistical, experimental, and theoretical uncertainties is given.

Post-fit Optimal Observable distributions in the low BDT score CR for the $\tau_{\mathrm{lep}}\tau_{\mathrm{lep}}$ SF analysis channel. ''Other bkg'' denotes all background contributions not listed explicitly in the legend. The size of the combined statistical, experimental, and theoretical uncertainties is given. The exact value of the $p_{\mathrm{T}}$ cut on the leptons depends on the trigger.

Post-fit Optimal Observable distributions in the low BDT score CR for the $\tau_{\mathrm{lep}}\tau_{\mathrm{lep}}$ DF analysis channel. ''Other bkg'' denotes all background contributions not listed explicitly in the legend. The size of the combined statistical, experimental, and theoretical uncertainties is given.

Post-fit Optimal Observable distributions in the low BDT score CR for the $\tau_{\mathrm{lep}}\tau_{\mathrm{had}}$ analysis channel. ''Other bkg'' denotes all background contributions not listed explicitly in the legend. The size of the combined statistical, experimental, and theoretical uncertainties is given. The exact value of the $p_{\mathrm{T}}$ cut on the leading lepton depends on the trigger.

Post-fit Optimal Observable distributions in the low BDT score CR for the $\tau_{\mathrm{had}}\tau_{\mathrm{had}}$ analysis channel. ''Other bkg'' denotes all background contributions not listed explicitly in the legend. The size of the combined statistical, experimental, and theoretical uncertainties is given.

Post-fit distributions of the event yields as a function of the Optimal Observable in the SR for the $\tau_{\mathrm{lep}}\tau_{\mathrm{lep}}$ SF analysis channel. The values of $\tilde d$, the signal strength $\mu$, the normalization of background processes, and nuisance parameters for the event yield prediction are set to those which minimize the NLL. The size of the combined statistical, experimental and theoretical uncertainties is given.

Post-fit distributions of the event yields as a function of the Optimal Observable in the SR for the $\tau_{\mathrm{lep}}\tau_{\mathrm{lep}}$ DF analysis channel. The values of $\tilde d$, the signal strength $\mu$, the normalization of background processes, and nuisance parameters for the event yield prediction are set to those which minimize the NLL. The size of the combined statistical, experimental and theoretical uncertainties is given.

Post-fit distributions of the event yields as a function of the Optimal Observable in the SR for the $\tau_{\mathrm{lep}}\tau_{\mathrm{had}}$ analysis channel. The values of $\tilde d$, the signal strength $\mu$, the normalization of background processes, and nuisance parameters for the event yield prediction are set to those which minimize the NLL. The size of the combined statistical, experimental and theoretical uncertainties is given. The exact value of the $p_{\mathrm{T}}$ cut on the leading lepton depends on the trigger.

Post-fit distributions of the event yields as a function of the Optimal Observable in the SR for the $\tau_{\mathrm{had}}\tau_{\mathrm{had}}$ analysis channel. The values of $\tilde d$, the signal strength $\mu$, the normalization of background processes, and nuisance parameters for the event yield prediction are set to those which minimize the NLL. The size of the combined statistical, experimental and theoretical uncertainties is given.

The observed $\Delta\mathrm{NLL}$ curve as a function of $\tilde d$ values. For comparison, expected $\Delta\mathrm{NLL}$ curves are also shown. The constraints on the nuisance parameters and normalization factors are first determined in a CR-only fit, and then a fit to pseudo-data corresponding to these nuisance parameters, normalization factors, and to $\tilde d=0, \mu = 1$ or $\tilde d =0, \mu = 0.73$ is performed to obtain these $\Delta\mathrm{NLL}$ curves. Moreover, a pre-fit expected $\Delta\mathrm{NLL}$ is shown, using pseudo-data corresponding to $\tilde d =0$ and $\mu = 1$ in the signal and control regions.

The expected $\Delta\mathrm{NLL}$ curves comparing the sensitivity of the fit with and without systematic uncertainties. For comparison, other curves are shown which remove the effect of jet-based systematic uncertainties, $\tau$-based systematic uncertainties, and MC statistical uncertainties.

The observed $\Delta\mathrm{NLL}$ curves for each analysis channel as a function of $\tilde d$, compared to the combined result. For the individual analysis channel $\Delta\mathrm{NLL}$ curves, only event yield information in the other SRs is used, ensuring that the Optimal Observable distributions in the other SRs do not influence the preferred value of $\tilde d$. The signal strength is constrained to be positive in these individual channel $\Delta\mathrm{NLL}$ curves. The exact value of the $p_{\mathrm{T}}$ cut on the leading lepton depends on the trigger.

Post-fit BDT distributions in the top-quark CR for the $\tau_{\mathrm{lep}}\tau_{\mathrm{lep}}$ SF channel. The size of the combined statistical, experimental, and theoretical uncertainties is given. The exact value of the $p_{\mathrm{T}}$ cut on the leptons depends on the trigger.

Post-fit BDT distributions in the top-quark CR for the $\tau_{\mathrm{lep}}\tau_{\mathrm{lep}}$ DF channel. The size of the combined statistical, experimental, and theoretical uncertainties is given.

Post-fit BDT distributions in the $Z\to \ell\ell$ CR for the $\tau_{\mathrm{lep}}\tau_{\mathrm{lep}}$ SF analysis channel. The size of the combined statistical, experimental, and theoretical uncertainties is given. The exact value of the $p_{\mathrm{T}}$ cut on the leptons depends on the trigger.

Post-fit Optimal Observable distributions in the top-quark CR for the $\tau_{\mathrm{lep}}\tau_{\mathrm{lep}}$ SF channel. The size of the combined statistical, experimental, and theoretical uncertainties is given. The exact value of the $p_{\mathrm{T}}$ cut on the leptons depends on the trigger.

Post-fit Optimal Observable distributions in the top-quark CR for the $\tau_{\mathrm{lep}}\tau_{\mathrm{lep}}$ DF channel. The size of the combined statistical, experimental, and theoretical uncertainties is given.

Post-fit Optimal Observable distributions in the $Z\to \ell\ell$ CR for the $\tau_{\mathrm{lep}}\tau_{\mathrm{lep}}$ SF analysis channel. The size of the combined statistical, experimental, and theoretical uncertainties is given. The exact value of the $p_{\mathrm{T}}$ cut on the leptons depends on the trigger.

Post-fit distribution of weighted event yields as a function of the Optimal Observable for all four SRs combined. The contributions of the different SRs are weighted by a factor of ln(1 + S/B), where S and B are the post-fit expected numbers of signal and background events in that region, respectively. The size of the combined statistical, experimental, and theoretical uncertainties is given.

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.

Distribution of $v_{3}$ from MBT events plotted as a function of the A-particle $p_\mathrm{T}$ for 0-5% centrality.

Distribution of $v_{3}$ from $p_{T}^{jet}>75$ GeV events plotted as a function of the A-particle $p_\mathrm{T}$ for 0-5% centrality.

Distribution of $v_{3}$ from $p_{T}^{jet}>100$ GeV events plotted as a function of the A-particle $p_\mathrm{T}$ for 0-5% centrality.

Distribution of $v_{2}$ from MBT events plotted as a function of the event centrality and A-particle $p_\mathrm{T}$ in (0.5-2 GeV).

Distribution of $v_{2}$ from $p_{T}^{jet}>75$ GeV events plotted as a function of the event centrality and A-particle $p_\mathrm{T}$ in (0.5-2 GeV).

Distribution of $v_{2}$ from $p_{T}^{jet}>100$ GeV events plotted as a function of the event centrality and A-particle $p_\mathrm{T}$ in (0.5-2 GeV).

Distribution of $v_{2}$ from MBT events plotted as a function of the event centrality and A-particle $p_\mathrm{T}$ in (2-9 GeV).

Distribution of $v_{2}$ from $p_{T}^{jet}>75$ GeV events plotted as a function of the event centrality and A-particle $p_\mathrm{T}$ in (2-9 GeV).

Distribution of $v_{2}$ from $p_{T}^{jet}>100$ GeV events plotted as a function of the event centrality and A-particle $p_\mathrm{T}$ in (2-9 GeV).

Distribution of $v_{2}$ from MBT events plotted as a function of the event centrality and A-particle $p_\mathrm{T}$ in (9-100 GeV).

Distribution of $v_{2}$ from $p_{T}^{jet}>75$ GeV events plotted as a function of the event centrality and A-particle $p_\mathrm{T}$ in (9-100 GeV).

Distribution of $v_{2}$ from $p_{T}^{jet}>100$ GeV events plotted as a function of the event centrality and A-particle $p_\mathrm{T}$ in (9-100 GeV).

Distribution of relative UE-UE pair fractions from MBT events plotted as a function of A-particle $p_\mathrm{T}$ for 0-5% centrality.

Distribution of relative HS-HS pair fractions from MBT events plotted as a function of A-particle $p_\mathrm{T}$ for 0-5% centrality.

Distribution of relative HS-UE pair fractions from MBT events plotted as a function of A-particle $p_\mathrm{T}$ for 0-5% centrality.

Distribution of relative UE-HS pair fractions from MBT events plotted as a function of A-particle $p_\mathrm{T}$ for 0-5% centrality.

Distribution of relative UE-UE pair fractions from $p_{T}^{jet}>75$ GeV events plotted as a function of A-particle $p_\mathrm{T}$ for 0-5% centrality.

Distribution of relative HS-HS pair fractions from $p_{T}^{jet}>75$ GeV events plotted as a function of A-particle $p_\mathrm{T}$ for 0-5% centrality.

Distribution of relative HS-UE pair fractions from $p_{T}^{jet}>75$ GeV events plotted as a function of A-particle $p_\mathrm{T}$ for 0-5% centrality.

Distribution of relative UE-HS pair fractions from $p_{T}^{jet}>75$ GeV events plotted as a function of A-particle $p_\mathrm{T}$ for 0-5% centrality.

Distribution of relative UE-UE pair fractions from $p_{T}^{jet}>100$ GeV events plotted as a function of A-particle $p_\mathrm{T}$ for 0-5% centrality.

Distribution of relative HS-HS pair fractions from $p_{T}^{jet}>100$ GeV events plotted as a function of A-particle $p_\mathrm{T}$ for 0-5% centrality.

Distribution of relative HS-UE pair fractions from $p_{T}^{jet}>100$ GeV events plotted as a function of A-particle $p_\mathrm{T}$ for 0-5% centrality.

Distribution of relative UE-HS pair fractions from $p_{T}^{jet}>100$ GeV events plotted as a function of A-particle $p_\mathrm{T}$ for 0-5% centrality.

Distribution of relative UE-UE pair fractions from MBT events plotted as a function of centrality A-particle $p_\mathrm{T}$ in (0.5-2 GeV).

Distribution of relative HS-UE pair fractions from MBT events plotted as a function of centrality A-particle $p_\mathrm{T}$ in (0.5-2 GeV).

Distribution of relative UE-UE pair fractions from $p_{T}^{jet}>75$ GeV events plotted as a function of centrality A-particle $p_\mathrm{T}$ in (0.5-2 GeV).

Distribution of relative HS-UE pair fractions from $p_{T}^{jet}>75$ GeV events plotted as a function of centrality A-particle $p_\mathrm{T}$ in (0.5-2 GeV).

Distribution of relative UE-UE pair fractions from $p_{T}^{jet}>100$ GeV events plotted as a function of centrality A-particle $p_\mathrm{T}$ in (0.5-2 GeV).

Distribution of relative HS-UE pair fractions from $p_{T}^{jet}>100$ GeV events plotted as a function of centrality A-particle $p_\mathrm{T}$ in (0.5-2 GeV).

Distribution of relative UE-UE pair fractions from MBT events plotted as a function of centrality A-particle $p_\mathrm{T}$ in (2-9 GeV).

Distribution of relative HS-UE pair fractions from MBT events plotted as a function of centrality A-particle $p_\mathrm{T}$ in (2-9 GeV).

Distribution of relative UE-UE pair fractions from $p_{T}^{jet}>75$ GeV events plotted as a function of centrality A-particle $p_\mathrm{T}$ in (2-9 GeV).

Distribution of relative HS-UE pair fractions from $p_{T}^{jet}>75$ GeV events plotted as a function of centrality A-particle $p_\mathrm{T}$ in (2-9 GeV).

Distribution of relative UE-UE pair fractions from $p_{T}^{jet}>100$ GeV events plotted as a function of centrality A-particle $p_\mathrm{T}$ in (2-9 GeV).

Distribution of relative HS-UE pair fractions from $p_{T}^{jet}>100$ GeV events plotted as a function of centrality A-particle $p_\mathrm{T}$ in (2-9 GeV).

Distribution of relative UE-UE pair fractions from MBT events plotted as a function of centrality A-particle $p_\mathrm{T}$ in (9-100 GeV).

Distribution of relative HS-UE pair fractions from MBT events plotted as a function of centrality A-particle $p_\mathrm{T}$ in (9-100 GeV).

Distribution of relative UE-UE pair fractions from $p_{T}^{jet}>75$ GeV events plotted as a function of centrality A-particle $p_\mathrm{T}$ in (9-100 GeV).

Distribution of relative HS-UE pair fractions from $p_{T}^{jet}>75$ GeV events plotted as a function of centrality A-particle $p_\mathrm{T}$ in (9-100 GeV).

Distribution of relative UE-UE pair fractions from $p_{T}^{jet}>100$ GeV events plotted as a function of centrality A-particle $p_\mathrm{T}$ in (9-100 GeV).

Distribution of relative HS-UE pair fractions from $p_{T}^{jet}>100$ GeV events plotted as a function of centrality A-particle $p_\mathrm{T}$ in (9-100 GeV).