Observed and expected limits on $\tan\beta$ as a function of $m_{H^+}$ in the $M_h^{125}(\tilde{\chi})$ scenario. Limits are shown for $\tan\beta$ values in the range of 0.5-60 due to the availability of the model prediction. The bands surrounding the expected limits show the 68% and 95% confidence intervals. Uncertainties in the predicted $H^+$ cross-sections or branching ratios are not considered.

Observed and expected limits on $\tan\beta$ as a function of $m_{H^+}$ in the $M_h^{125}(\tilde{\tau})$ scenario. Limits are shown for $\tan\beta$ values in the range of 0.5-60 due to the availability of the model prediction. The bands surrounding the expected limits show the 68% and 95% confidence intervals. Uncertainties in the predicted $H^+$ cross-sections or branching ratios are not considered.

Observed and expected limits on $\tan\beta$ as a function of $m_{H^+}$ in the $M_h^{125}$(alignment) scenario. Limits are shown for $\tan\beta$ values in the range of 0.5-60 due to the availability of the model prediction. The bands surrounding the expected limits show the 68% and 95% confidence intervals. Uncertainties in the predicted $H^+$ cross-sections or branching ratios are not considered.

Observed and expected limits on $\tan\beta$ as a function of $m_{H^+}$ in the $M_{h_1}^{125}$(CPV) scenario. Limits are shown for $\tan\beta$ values in the range of 0.5-60 due to the availability of the model prediction. The bands surrounding the expected limits show the 68% and 95% confidence intervals. Uncertainties in the predicted $H^+$ cross-sections or branching ratios are not considered.

Event yields of the SM background processes and the 800 GeV $H^{+}$ sample in the four analysis regions before the fit to the data. Uncertainties include both statistical and systematic uncertainties. The yields of the $H^{+}$ signal sample correspond to a cross-section times branching fraction of 10 pb.

Event acceptance for the different $H^+$ mass signal samples.

Distribution of $S$, which is one of the discriminating variables used to define the $2\ell^{sc}$ SRs. The events are selected with the preselection requirements listed in Table 4 in the paper. The data (dots) are compared with the expected contributions from the relevant background sources (histograms). The expected signal distributions for $m_{H^{\pm\pm}} = 300~GeV$ are also shown, scaled to the observed number of events. The last bin includes overflows.

Distribution of $E_{T}^{miss}$, which is one of the discriminating variables used to define the $3\ell$ SRs. The events are selected with the preselection requirements listed in Table 4 in the paper. The data (dots) are compared with the expected contributions from the relevant background sources (histograms). The expected signal distributions for $m_{H^{\pm\pm}} = 300~GeV$ are also shown, scaled to the observed number of events. The last bin includes overflows.

Distribution of $\Delta R_{\ell^{\pm}\ell^{\pm}}$, which is one of the discriminating variables used to define the $3\ell$ SRs. The events are selected with the preselection requirements listed in Table 4 in the paper. The data (dots) are compared with the expected contributions from the relevant background sources (histograms). The expected signal distributions for $m_{H^{\pm\pm}} = 300~GeV$ are also shown, scaled to the observed number of events. The last bin includes overflows.

Distribution of $m_{x\ell}$ ($x$=3), which is one of the discriminating variables used to define the $3\ell$ SRs. The events are selected with the preselection requirements listed in Table 4 in the paper. The data (dots) are compared with the expected contributions from the relevant background sources (histograms). The expected signal distributions for $m_{H^{\pm\pm}} = 300~GeV$ are also shown, scaled to the observed number of events. The last bin includes overflows.

Distribution of $p_{T}^{leading jet}$, which is one of the discriminating variables used to define the $3\ell$ SRs. The events are selected with the preselection requirements listed in Table 4 in the paper. The data (dots) are compared with the expected contributions from the relevant background sources (histograms). The expected signal distributions for $m_{H^{\pm\pm}} = 300~GeV$ are also shown, scaled to the observed number of events. The last bin includes overflows.

Distribution of $E_{T}^{miss}$, which is one of the discriminating variables used to define the $4\ell$ SRs. The events are selected with the preselection requirements listed in Table 4 in the paper. The data (dots) are compared with the expected contributions from the relevant background sources (histograms). The expected signal distributions for $m_{H^{\pm\pm}} = 300~GeV$ are also shown, scaled to the observed number of events. The last bin includes overflows.

Distribution of $\Delta R_{\ell^{\pm}\ell^{\pm}}^{min}$, which is one of the discriminating variables used to define the $4\ell$ SRs. The events are selected with the preselection requirements listed in Table 4 in the paper. The data (dots) are compared with the expected contributions from the relevant background sources (histograms). The expected signal distributions for $m_{H^{\pm\pm}} = 300~GeV$ are also shown, scaled to the observed number of events. The last bin includes overflows.

Distribution of $m_{x\ell}$ ($x$=4), which is one of the discriminating variables used to define the $4\ell$ SRs. The events are selected with the preselection requirements listed in Table 4 in the paper. The data (dots) are compared with the expected contributions from the relevant background sources (histograms). The expected signal distributions for $m_{H^{\pm\pm}} = 300~GeV$ are also shown, scaled to the observed number of events. The last bin includes overflows.

Distribution of $p_{T}^{\ell_{1}}$, which is one of the discriminating variables used to define the $4\ell$ SRs. The events are selected with the preselection requirements listed in Table 4 in the paper. The data (dots) are compared with the expected contributions from the relevant background sources (histograms). The expected signal distributions for $m_{H^{\pm\pm}} = 300~GeV$ are also shown, scaled to the observed number of events. The last bin includes overflows.

Contributions from different categories of uncertainties relative to the expected background yields in the defined SRs, as obtained after performing the likelihood ratio test discussed in Section 9 in the paper. The uncertainties are shown for the combination of the individual channels of the $2\ell^{sc}$, $3\ell$ and $4\ell$ SRs. The SRs are indicated along the horizontal axis. In the HEPData entry, the x-axis is simplified for easier visualisation. The first number indicates the sub channel (2:$2\ell^{sc}$, 3:$3\ell$, 4:$4\ell$), while the second number indicates the mass point (2:200, 3:300, 4:400, 5:500).

Data event yields compared with the expected contributions from relevant background sources, for the combination of the individual channels of the $2\ell^{sc}$, $3\ell$ and $4\ell$ SRs. The total uncertainties in the expected event yields are shown as the hatched bands. The SRs are indicated along the horizontal axis. In the HEPData entry, the x-axis is simplified for easier visualisation. The first number indicates the sub channel (2:$2\ell^{sc}$, 3:$3\ell$, 4:$4\ell$), while the second number indicates the mass point (2:200, 3:300, 4:400, 5:500).

The $E_{T}^{miss}$ distribution for the SRs of the $m_{H^{\pm\pm}} = 300~GeV$ signal mass hypothesis, where the selection requirement on $E_{T}^{miss}$ has been removed. In the attached plot, the signals are stacked on top of the backgrounds while individuals contributions of the $2\ell^{sc}$ channel are shown in HEPData. The last bin, isolated by a vertical red dashed line, is inclusive and corresponds to the SR.

The $E_{T}^{miss}$ distribution for the SRs of the $m_{H^{\pm\pm}} = 300~GeV$ signal mass hypothesis, where the selection requirement on $E_{T}^{miss}$ has been removed. In the attached plot, the signals are stacked on top of the backgrounds while individuals contributions of the $3\ell$ channel are shown in HEPData. The last bin, isolated by a vertical red dashed line, is inclusive and corresponds to the SR.

The $E_{T}^{miss}$ distribution for the SRs of the $m_{H^{\pm\pm}} = 300~GeV$ signal mass hypothesis, where the selection requirement on $E_{T}^{miss}$ has been removed. In the attached plot, the signals are stacked on top of the backgrounds while individuals contributions of the $4\ell$ channel are shown in HEPData. The last bin, isolated by a vertical red dashed line, is inclusive and corresponds to the SR.

Observed and expected upper limits of the $H^{\pm\pm}$ pair production cross section times branching fraction at 95% CL obtained from the combination of 2$\ell^{sc}$, 3$\ell$ and 4$\ell$ channels. The region above the observed limit is excluded by the measurement. The bands represent the expected exclusion curves within one and two standard deviations.

The theoretical prediction of Figure 9(a) in the paper.

Observed and expected upper limits of the $H^{\pm\pm}$ and $H^{\pm}$ production cross section times branching fraction at 95% CL obtained from the combination of 2$\ell^{sc}$, 3$\ell$ and 4$\ell$ channels. The region above the observed limit is excluded by the measurement. The bands represent the expected exclusion curves within one and two standard deviations.

The theoretical prediction of Figure 9(b) in the paper.

Data event yields compared with the estimated background in the $m_{H^{\pm\pm}} = 200~GeV$ or $m_{H^{\pm\pm}} = 220~GeV$ SRs. SFOC 0 and SFOC 1,2 refer to the number of same-flavour opposite charge lepton pairs. The total uncertainties in the estimated background yields are shown as the hashed bands. In the HEPData entry, the x-axis is simplified for easier visualisation (1:$e^{\pm}e^{\pm}$, 2:$e^{\pm}\mu^{\pm}$, 3:$\mu^{\pm}\mu^{\pm}$, 4:SFOC 0, 5:SFOC 1,2, 6:$4\ell$).

Data event yields compared with the estimated background in the $m_{H^{\pm\pm}} = 300~GeV$ or $m_{H^{\pm\pm}} = 350~GeV$ SRs. SFOC 0 and SFOC 1,2 refer to the number of same-flavour opposite charge lepton pairs. The total uncertainties in the estimated background yields are shown as the hashed bands. In the HEPData entry, the x-axis is simplified for easier visualisation (1:$e^{\pm}e^{\pm}$, 2:$e^{\pm}\mu^{\pm}$, 3:$\mu^{\pm}\mu^{\pm}$, 4:SFOC 0, 5:SFOC 1,2, 6:$4\ell$).

Data event yields compared with the estimated background in the $m_{H^{\pm\pm}} = 400~GeV$ or $m_{H^{\pm\pm}} = 450~GeV$ SRs. SFOC 0 and SFOC 1,2 refer to the number of same-flavour opposite charge lepton pairs. The total uncertainties in the estimated background yields are shown as the hashed bands. In the HEPData entry, the x-axis is simplified for easier visualisation (1:$e^{\pm}e^{\pm}$, 2:$e^{\pm}\mu^{\pm}$, 3:$\mu^{\pm}\mu^{\pm}$, 4:SFOC 0, 5:SFOC 1,2, 6:$4\ell$).

Data event yields compared with the estimated background in the $m_{H^{\pm\pm}} = 500~GeV$ or $m_{H^{\pm\pm}} = 550~GeV$ or $m_{H^{\pm\pm}} = 600~GeV$ SRs. SFOC 0 and SFOC 1,2 refer to the number of same-flavour opposite charge lepton pairs. The total uncertainties in the estimated background yields are shown as the hashed bands. In the HEPData entry, the x-axis is simplified for easier visualisation (1:$e^{\pm}e^{\pm}$, 2:$e^{\pm}\mu^{\pm}$, 3:$\mu^{\pm}\mu^{\pm}$, 4:SFOC 0, 5:SFOC 1,2, 6:$4\ell$).

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

Summary of the observed and expected 95% CL upper limits on B as a function of $m_{\mathrm{LQ}_{3}^{\mathrm{d}}}$.

Summary of the observed and expected 95% CL upper limits on B as a function of $m_{\mathrm{LQ}_{3}^{\mathrm{d}}}$.

Cutflow of the preselection requirements (see Section 5) for $\mathrm{LQ}_{3}^{\mathrm{d}}$ signals with $m_{\mathrm{LQ}_{3}^{\mathrm{d}}}$=0.9, 1.1, and 1.3 TeV, assuming B=1. The yields correspond to an integrated luminosity of 139 fb$^{-1}$.

Cutflow of the preselection requirements (see Section 5) for $\mathrm{LQ}_{3}^{\mathrm{d}}$ signals with $m_{\mathrm{LQ}_{3}^{\mathrm{d}}}$=0.9, 1.1, and 1.3 TeV, assuming B=1. The yields correspond to an integrated luminosity of 139 fb$^{-1}$.

Cutflow of the signal region requirements in 1$\ell+\geq 1\tau$ channel (see Table 3) for $\mathrm{LQ}_{3}^{\mathrm{d}}$ signals with $m_{\mathrm{LQ}_{3}^{\mathrm{d}}}$=0.9, 1.1, and 1.3 TeV, assuming B=1. Events that satisfy the preselection requirements are considered. The yields correspond to an integrated luminosity of 139 fb$^{-1}$.

Cutflow of the signal region requirements in 1$\ell+\geq 1\tau$ channel (see Table 3) for $\mathrm{LQ}_{3}^{\mathrm{d}}$ signals with $m_{\mathrm{LQ}_{3}^{\mathrm{d}}}$=0.9, 1.1, and 1.3 TeV, assuming B=1. Events that satisfy the preselection requirements are considered. The yields correspond to an integrated luminosity of 139 fb$^{-1}$.

Cutflow of the signal region requirements in the 2$\ell$OS+$\geq 1\tau$ channel (see Table 4) for $\mathrm{LQ}_{3}^{\mathrm{d}}$ signals with $m_{\mathrm{LQ}_{3}^{\mathrm{d}}}$=0.9, 1.1, and 1.3 TeV, assuming B=1. Events that satisfy the preselection requirements are considered. The yields correspond to an integrated luminosity of 139 fb$^{-1}$.

Cutflow of the signal region requirements in the 2$\ell$OS+$\geq 1\tau$ channel (see Table 4) for $\mathrm{LQ}_{3}^{\mathrm{d}}$ signals with $m_{\mathrm{LQ}_{3}^{\mathrm{d}}}$=0.9, 1.1, and 1.3 TeV, assuming B=1. Events that satisfy the preselection requirements are considered. The yields correspond to an integrated luminosity of 139 fb$^{-1}$.

Cutflow of the signal region requirements in the $2\ell$SS/$3\ell+\geq 1\tau$ channel (see Table 5) for $\mathrm{LQ}_{3}^{\mathrm{d}}$ signals with $m_{\mathrm{LQ}_{3}^{\mathrm{d}}}$=0.9, 1.1, and 1.3 TeV, assuming B=1. Events that satisfy the preselection requirements are considered. In this channel, two signal regions (SR-L and SR-H) are defined based on $p_{\mathrm{T}, 1}^{\tau}$, with SR-L and SR-H requiring $125< p_{\mathrm{T}, 1}^{\tau} < 225$ GeV and $p_{\mathrm{T}, 1}^{\tau}>225$ GeV, respectively. The yields correspond to an integrated luminosity of 139 fb$^{-1}$.

Cutflow of the signal region requirements in the $2\ell$SS/$3\ell+\geq 1\tau$ channel (see Table 5) for $\mathrm{LQ}_{3}^{\mathrm{d}}$ signals with $m_{\mathrm{LQ}_{3}^{\mathrm{d}}}$=0.9, 1.1, and 1.3 TeV, assuming B=1. Events that satisfy the preselection requirements are considered. In this channel, two signal regions (SR-L and SR-H) are defined based on $p_{\mathrm{T}, 1}^{\tau}$, with SR-L and SR-H requiring $125< p_{\mathrm{T}, 1}^{\tau} < 225$ GeV and $p_{\mathrm{T}, 1}^{\tau}>225$ GeV, respectively. The yields correspond to an integrated luminosity of 139 fb$^{-1}$.

Acceptances on TRUTH level of the DMt samples in the tW1L channel for all bins in the SR. The acceptance is defined as the number of weighted TRUTH events in the SR over the number of expected events without any selections. The maps include all samples in the $m_H - tan\beta$ plane with $m_a = 250~GeV$.

Efficiencies of the DMt samples in the tW2L SR. The efficiency is defined as the number of weighted reconstructed events over the number of weighted TRUTH events in the SR. The maps include all samples in the $m_a - m_H$ plane with $tan\beta = 1$.

Acceptances on TRUTH level of the DMt samples in the tW2L SR. The acceptance is defined as the number of weighted TRUTH events in the SR over the number of expected events without any selections. The maps include all samples in the $m_a - m_H$ plane with $tan\beta = 1$.

Efficiencies of the DMt samples in the tW2L SR. The efficiency is defined as the number of weighted reconstructed events over the number of weighted TRUTH events in the SR. The maps include all samples in the $m_H - tan\beta$ plane with $m_a = 250~GeV$.

Acceptances on TRUTH level of the DMt samples in the tW2L SR. The acceptance is defined as the number of weighted TRUTH events in the SR over the number of expected events without any selections. The maps include all samples in the $m_H - tan\beta$ plane with $m_a = 250~GeV$.

Efficiencies of the DMt samples in the tj1L channel for all bins in the SR. The efficiency is defined as the number of weighted reconstructed events over the number of weighted TRUTH events in the SR. The map includes all used samples in the $m_H - tan\beta$ plane with $m_a = 250~GeV$.

Acceptances on TRUTH level of the DMt samples in the tj1L channel for all bins in the SR. The acceptance is defined as the number of weighted TRUTH events in the SR over the number of expected events without any selections. The map includes all used samples in the $m_H - tan\beta$ plane with $m_a = 250~GeV$.

Upper limits on signal strength (excluded cross section over theoretical cross section) of the tW1L analysis considering only DMt signal.

Upper limits on excluded cross sections of the tW1L analysis considering only the DMt signal.

The expected exclusion contours as a function of $(m_a, m_{H^{\pm}})$, assuming only $tW$+DM contributions, for the tW1L analysis channel.

The observed exclusion contours as a function of $(m_a, m_{H^{\pm}})$, assuming only $tW$+DM contributions, for the tW1L analysis channel.

Upper limits on signal strength (excluded cross section over theoretical cross section) of the tW1L analysis considering only DMt signal.

Upper limits on excluded cross sections of the tW1L analysis considering only the DMt signal.

The expected exclusion contours as a function of $(m_{H^{\pm}}, \tan\beta)$, assuming only $tW$+DM contributions, for the tW1L analysis channel.

The observed exclusion contours as a function of $(m_{H^{\pm}}, \tan\beta)$, assuming only $tW$+DM contributions, for the tW1L analysis channel.

Upper limits on signal strength (excluded cross section over theoretical cross section) of the tW2L analysis considering only DMt signal.

Upper limits on excluded cross sections of the tW2L analysis considering only the DMt signal.

The expected exclusion contours as a function of $(m_a, m_{H^{\pm}})$, assuming only $tW$+DM contributions, for the tW2L analysis channel.

The observed exclusion contours as a function of $(m_a, m_{H^{\pm}})$, assuming only $tW$+DM contributions, for the tW2L analysis channel.

Upper limits on signal strength (excluded cross section over theoretical cross section) of the tW2L analysis considering only DMt signal.

Upper limits on excluded cross sections of the tW2L analysis considering only the DMt signal.

The expected exclusion contours as a function of $(m_{H^{\pm}}, \tan\beta)$, assuming only $tW$+DM contributions, for the tW2L analysis channel.

The observed exclusion contours as a function of $(m_{H^{\pm}}, \tan\beta)$, assuming only $tW$+DM contributions, for the tW2L analysis channel.

Upper limits on signal strength (excluded cross section over theoretical cross section) of the combined tW1L and tW2L analyses considering only the DMt signal.

Upper limits on excluded cross sections of the combined tW1L and tW2L analyses considering only the DMt signal.

The expected exclusion contours as a function of $(m_a, m_{H^{\pm}})$, assuming only $tW$+DM contributions, for the statistical combination of the tW1L and tW2L analysis channel.

The observed exclusion contours as a function of $(m_a, m_{H^{\pm}})$, assuming only $tW$+DM contributions, for the statistical combination of the tW1L and tW2L analysis channel.

Upper limits on signal strength (excluded cross section over theoretical cross section) of the combined tW1L and tW2L analyses considering only the DMt signal.

Upper limits on excluded cross sections of the combined tW1L and tW2L analyses considering only the DMt signal.

The expected exclusion contours as a function of $(m_{H^{\pm}}, \tan\beta)$, assuming only $tW$+DM contributions, for the statistical combination of the tW1L and tW2L analysis channel.

The observed exclusion contours as a function of $(m_{H^{\pm}}, \tan\beta)$, assuming only $tW$+DM contributions, for the statistical combination of the tW1L and tW2L analysis channel.

Upper limits on signal strength (excluded cross section over theoretical cross section) of the tW1L analysis considering the DMt$\bar{t}$+DMt signal.

The expected exclusion contours as a function of $(m_a, m_{H^{\pm}})$, assuming DM$t\bar{t}$ and DM$t$ contributions, for the tW1L analysis channel.

The observed exclusion contours as a function of $(m_a, m_{H^{\pm}})$, assuming DM$t\bar{t}$ and DM$t$ contributions, for the tW1L analysis channel.

Upper limits on signal strength (excluded cross section over theoretical cross section) of the tW1L analysis considering the DMt$\bar{t}$+DMt signal.

The expected exclusion contours as a function of $(m_{H^{\pm}}, \tan\beta)$, assuming DM$t\bar{t}$ and DM$t$ contributions, for the tW1L analysis channel.

The observed exclusion contours as a function of $(m_{H^{\pm}}, \tan\beta)$, assuming DM$t\bar{t}$ and DM$t$ contributions, for the tW1L analysis channel.

Upper limits on signal strength (excluded cross section over theoretical cross section) of the tW2L analysis considering the DMt$\bar{t}$+DMt signal.

The expected exclusion contours as a function of $(m_a, m_{H^{\pm}})$, assuming DM$t\bar{t}$ and DM$t$ contributions, for the tW2L analysis channel.

The observed exclusion contours as a function of $(m_a, m_{H^{\pm}})$, assuming DM$t\bar{t}$ and DM$t$ contributions, for the tW2L analysis channel.

Upper limits on signal strength (excluded cross section over theoretical cross section) of the tW2L analysis considering the DMt$\bar{t}$+DMt signal.

The expected exclusion contours as a function of $(m_{H^{\pm}}, \tan\beta)$, assuming DM$t\bar{t}$ and DM$t$ contributions, for the tW2L analysis channel.

The observed exclusion contours as a function of $(m_{H^{\pm}}, \tan\beta)$, assuming DM$t\bar{t}$ and DM$t$ contributions, for the tW2L analysis channel.

Upper limits on signal strength (excluded cross section over theoretical cross section) of the combined tW1L and tW2L analyses considering the DMt$\bar{t}$+DMt signal.

The expected exclusion contours as a function of $(m_a, m_{H^{\pm}})$, assuming DM$t\bar{t}$ and DM$t$ contributions, for the statistical combination of the tW1L and tW2L analysis channel.

The observed exclusion contours as a function of $(m_a, m_{H^{\pm}})$, assuming DM$t\bar{t}$ and DM$t$ contributions, for the statistical combination of the tW1L and tW2L analysis channel.

Upper limits on signal strength (excluded cross section over theoretical cross section) of the combined tW1L and tW2L analyses considering the DMt$\bar{t}$+DMt signal.

The expected exclusion contours as a function of $(m_{H^{\pm}}, \tan\beta)$, assuming DM$t\bar{t}$ and DM$t$ contributions, for the statistical combination of the tW1L and tW2L analysis channel.

The observed exclusion contours as a function of $(m_{H^{\pm}}, \tan\beta)$, assuming DM$t\bar{t}$ and DM$t$ contributions, for the statistical combination of the tW1L and tW2L analysis channel.

Upper limits on signal strength (excluded cross section over theoretical cross section) of the tj1L analysis considering only the DMt signal.

Upper limits on upper limits on excluded cross sections of the tj1L analysis considering only the DMt signal.

The expected and observed cross section exclusion limits as a function of $m_{H^{\pm}}$ in the tj1L analysis channel for signal models with $m_a = 250~GeV$, and $\tan\beta=0.3$. The $\sigma^{}_\mathrm{BSM}$ is the cross section of the $t$-channel DM production process.

The expected and observed cross section exclusion limits as a function of $m_{H^{\pm}}$ in the tj1L analysis channel for signal models with $m_a = 250~GeV$, and $\tan\beta=0.5$. The $\sigma^{}_\mathrm{BSM}$ is the cross section of the $t$-channel DM production process.

Cross sections of the DMt samples in the tW1L channel. The maps include all samples in the $m_a - m_H$ plane with $tan\beta = 1$.

Cross sections of the DMt samples in the tW1L channel. The maps include all samples in the $m_H - tan\beta$ plane with $m_a = 250~GeV$.

Cross sections times branching ratio of the DMt samples in the tW2L channel. The maps include all samples in the $m_a - m_H$ plane with $tan\beta = 1$.

Cross sections times branching ratio of the DMt samples in the tW2L channel. The maps include all samples in the $m_H - tan\beta$ plane with $m_a = 250~GeV$.

Cross sections of the DMt samples in the tj1L channel. The map includes all samples in the $m_H - tan\beta$ plane with $m_a = 250~GeV$.

MC generator filter efficiencies of the DMt samples in the tW1L channel. The maps include all samples in the $m_a - m_H$ plane with $tan\beta = 1$.

MC generator filter efficiencies of the DMt samples in the tW1L channel. The maps include all samples in the $m_H - tan\beta$ plane with $m_a = 250~GeV$.

MC generator filter efficiencies of the DMt samples in the tW2L channel. The maps include all samples in the $m_a - m_H$ plane with $tan\beta = 1$.

MC generator filter efficiencies of the DMt samples in the tW2L channel. The maps include all samples in the $m_H - tan\beta$ plane with $m_a = 250~GeV$.

MC generator filter efficiencies of the DMt samples in the tj1L channel. The map includes all samples in the $m_H - tan\beta$ plane with $m_a = 250~GeV$.

Background-only fit results for the tW1L and tW2L signal regions. The backgrounds which contribute only a small amount (rare processes such as triboson, Higgs boson production processes, $t\bar{t}t\bar{t}$, $t\bar{t}WW$ and non-prompt or misidentified leptons background) are grouped and labelled as ``Others´´. The quoted uncertainties on the fitted SM background include both the statistical and systematic uncertainties.

Background-only fit results for the tj1L signal regions. The backgrounds which contribute only a small amount ($Z$+jets, rare processes such as $tWZ$, triboson, Higgs boson production processes, ,$t\bar{t}t\bar{t}$, $t\bar{t}WW$) are grouped and labelled as ``Others´´. The quoted uncertainties on the fitted SM background include both the statistical and systematic uncertainties.

Cutflow of the weighted events with statistical uncertainties for two DMt samples in all bins of the tW1L channel. The PreSelection includes at least 1 lepton in the event, at least 1 $b$-jet with $p_{\mathrm{T}} > 50~GeV$, $m\mathrm{_{T}^{lep}} > 30~GeV$, $\Delta\phi\mathrm{_{4jets, MET}^{min}} > 0.5$ and $E\mathrm{_{T}^{miss}} > 200~GeV$.

Cutflow of the weighted events with statistical uncertainties for two DMt samples in the tW2L channel. The PreSelection includes at least 2 leptons in the event, at least 1 $b$-jet with $p_{\mathrm{T}} > 40~GeV$, $m_{ll} > 40~GeV$, $m\mathrm{_{T2}} > 40~GeV$, $\Delta\phi\mathrm{_{4jets, MET}^{min}} > 0.5$ and $E\mathrm{_{T}^{miss}} > 200~GeV$.

Cutflow of the weighted events with the statistical uncertainties (except for the first cuts) for two DMt samples in all bins off the tj1L channel. The PreSelection includes at least 1 lepton in the event, at least 1 $b$-jet with $p_{\mathrm{T}} > 50~GeV$, $m\mathrm{_{T}^{lep}} > 30~GeV$, $\Delta\phi\mathrm{_{4jets, MET}^{min}} > 0.5$ and $E\mathrm{_{T}^{miss}} > 200~GeV$.

The best-fit expected and observed distributions of the combined NN output in the VRSS for the $e\tau$ channel for events with 3-prong $\tau_\text{had-vis}$ candidates. The last bin in each plot includes overflow events.

The best-fit expected and observed distributions of the combined NN output in the SR for the $e\tau$ channel for events with 1-prong $\tau_\text{had-vis}$ candidates. The last bin in each plot includes overflow events.

The best-fit expected and observed distributions of the combined NN output in the SR for the $e\tau$ channel for events with 3-prong $\tau_\text{had-vis}$ candidates. The last bin in each plot includes overflow events.

The best-fit expected and observed distributions of the combined NN output in the SR for the $\mu\tau$ channel for events with 1-prong $\tau_\text{had-vis}$ candidates. The last bin in each plot includes overflow events.

The best-fit expected and observed distributions of the combined NN output in the SR for the $\mu\tau$ channel for events with 3-prong $\tau_\text{had-vis}$ candidates. The last bin in each plot includes overflow events.

Observed and expected upper limits on $\mathcal{B}(Z\rightarrow\ell\tau)$ at 95% confidence level.

Measured normalised differential fiducial cross sections of $\gamma\gamma \rightarrow \gamma\gamma$ production in Pb+Pb collisions at $\sqrt{s_{NN}}$ = 5.02 TeV for diphoton $|cos(\theta*)|$ are shown as points with error bars giving the statistical uncertainty and grey bands indicating the size of the total uncertainty. The results are compared with the prediction from the SuperChic v3.0 MC generator (solid line).

Measured normalised differential fiducial cross sections of $\gamma\gamma \rightarrow \gamma\gamma$ production in Pb+Pb collisions at $\sqrt{s_{NN}}$ = 5.02 TeV for average photon transverse momentum are shown as points with error bars giving the statistical uncertainty and grey bands indicating the size of the total uncertainty. The results are compared with the prediction from the SuperChic v3.0 MC generator (solid line) with bands denoting the theoretical uncertainty.

Measured normalised differential fiducial cross sections of $\gamma\gamma \rightarrow \gamma\gamma$ production in Pb+Pb collisions at $\sqrt{s_{NN}}$ = 5.02 TeV for average photon transverse momentum are shown as points with error bars giving the statistical uncertainty and grey bands indicating the size of the total uncertainty. The results are compared with the prediction from the SuperChic v3.0 MC generator (solid line).

Measured differential fiducial cross sections of $\gamma\gamma \rightarrow \gamma\gamma$ production in Pb+Pb collisions at $\sqrt{s_{NN}}$ = 5.02 TeV for diphoton rapidity are shown as points with error bars giving the statistical uncertainty and grey bands indicating the size of the total uncertainty. The results are compared with the prediction from the SuperChic v3.0 MC generator (solid line) with bands denoting the theoretical uncertainty.

Measured normalised differential fiducial cross sections of $\gamma\gamma \rightarrow \gamma\gamma$ production in Pb+Pb collisions at $\sqrt{s_{NN}}$ = 5.02 TeV for diphoton rapidity are shown as points with error bars giving the statistical uncertainty and grey bands indicating the size of the total uncertainty. The results are compared with the prediction from the SuperChic v3.0 MC generator (solid line).

The 95% CL upper limit on the ALP cross section $\sigma_{\gamma\gamma\rightarrow a \rightarrow\gamma\gamma}$ for the $\gamma\gamma\rightarrow a \rightarrow \gamma\gamma$ process as a function of ALP mass m$_{a}$. The observed upper limit is shown as a solid black line and the expected upper limit is shown by the dashed black line with its $\pm1$ and $\pm2$ standard deviation bands. The discontinuity at $m_{a}=70 GeV$ is caused by the increase of the mass-bin width which brings an increase in signal acceptance.

The 95% CL upper limit on the ALP coupling $1/\Lambda_{a}$ for the $\gamma\gamma\rightarrow a \rightarrow \gamma\gamma$ process as a function of ALP mass m$_{a}$. The observed upper limit is shown as a solid black line and the expected upper limit is shown by the dashed black line with its $\pm1$ and $\pm2$ standard deviation bands. The discontinuity at $m_{a} = 70 GeV$ is caused by the increase of the mass-bin width which brings an increase in signal acceptance.

Fiducial cross section for light-by-light scattering

The number of data and fitted simulated events in each bin of the $|d_{0}^{\mu}|$ distribution in the $10<p_{\textrm{T}}^{\mu}<20$ GeV selection in the $\mu-\mu$ channel.

The number of data and fitted simulated events in each bin of the $|d_{0}^{\mu}|$ distribution in the $20<p_{\textrm{T}}^{\mu}<250$ GeV selection in the $e-\mu$ channel.

The number of data and fitted simulated events in each bin of the $|d_{0}^{\mu}|$ distribution in the $20<p_{\textrm{T}}^{\mu}<250$ GeV selection in the $\mu-\mu$ channel.

The measurement of the ratio of the rate of decay of W bosons to τ-leptons and muons, $R(\tau/\mu)=B(W\rightarrow\tau\nu_\tau)/B(W\rightarrow \mu\nu_\mu)$.

Measured and predicted fiducial cross-sections in both the lllljj and ll$\nu\nu$jj channels for the inclusive ZZjj processes. Uncertainties due to different sources are presented.

Measured and predicted fiducial cross-sections in both the lllljj and ll$\nu\nu$jj channels for the inclusive ZZjj processes. Uncertainties due to different sources are presented.

Observed and expected multivariate discriminant distribution in the $\ell\ell\ell\ell jj$ QCD CR.

Observed and expected multivariate discriminant distribution in the $\ell\ell\ell\ell jj$ QCD CR.

Observed and expected multivariate discriminant distribution in the $\ell\ell\ell\ell jj$ SR.

Observed and expected multivariate discriminant distribution in the $\ell\ell\ell\ell jj$ SR.

Observed and expected multivariate discriminant distribution in the $\ell\ell\nu\nu jj$ SR.

Observed and expected multivariate discriminant distribution in the $\ell\ell\nu\nu jj$ SR.

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.