Cutflow for signal samples in the muon-tau channel for 2016 signal samples. Each entry other than the total is the relative efficiency with respect to the previous selection.

Cutflow for signal samples in the muon-tau channel for 2017 signal samples. Each entry other than the total is the relative efficiency with respect to the previous selection.

Cutflow for signal samples in the muon-tau channel for 2018 signal samples. Each entry other than the total is the relative efficiency with respect to the previous selection.

Cutflow for signal samples in the electron-tau channel for 2016 signal samples. Each entry other than the total is the relative efficiency with respect to the previous selection.

Cutflow for signal samples in the electron-tau channel for 2017 signal samples. Each entry other than the total is the relative efficiency with respect to the previous selection.

Cutflow for signal samples in the electron-tau channel for 2018 signal samples. Each entry other than the total is the relative efficiency with respect to the previous selection.

Yields for the muon-tau channel in 2016 as a function of reconstructed mass. Data and each background are listed with their respective uncertainties. The last bin contains the overflow bin.

Yields for the muon-tau channel in 2017 as a function of reconstructed mass. Data and each background are listed with their respective uncertainties. The last bin contains the overflow bin.

Yields for the muon-tau channel in 2018 as a function of reconstructed mass. Data and each background are listed with their respective uncertainties. The last bin contains the overflow bin.

Yields for the electron-tau channel in 2016 as a function of reconstructed mass. Data and each background are listed with their respective uncertainties. The last bin contains the overflow bin.

Yields for the electron-tau channel in 2017 as a function of reconstructed mass. Data and each background are listed with their respective uncertainties. The last bin contains the overflow bin.

Yields for the electron-tau channel in 2018 as a function of reconstructed mass. Data and each background are listed with their respective uncertainties. The last bin contains the overflow bin.

Yields for the di-tau channel in 2016 as a function of reconstructed mass. Data and each background are listed with their respective uncertainties. The last bin contains the overflow bin.

Yields for the di-tau channel in 2017 as a function of reconstructed mass. Data and each background are listed with their respective uncertainties. The last bin contains the overflow bin.

Yields for the di-tau channel in 2018 as a function of reconstructed mass. Data and each background are listed with their respective uncertainties. The last bin contains the overflow bin.

Covariance matrix for all bins and all channels used in the background-only fit of this analysis. There are the e-tau, the mu-tau, and the tau-tau analysis channels, each present for 2016, 2017, and 2018. Each has 6 bins in the reconstructed mass (mrec[GeV]) indicated with their borders in the bin naming. The bin labels have channel name, year name, then bin name.

95% confidence level expected limit with 1 and 2 sigma uncertainties, as well as observed limits on Z' model cross section times branching ratio to di-tau in the e-tau channel.

95% confidence level expected limit with 1 and 2 sigma uncertainties, as well as observed limits on Z' model cross section times branching ratio to di-tau in the mu-tau channel.

95% confidence level expected limit with 1 and 2 sigma uncertainties, as well as observed limits on Z' model cross section times branching ratio to ditau in the ditau channel.

95% confidence level expected limit with 1 and 2 sigma uncertainties, as well as observed limits on Z' model cross section times branching ratio to ditau for the combination of all channels.

Exclusion limits on the production cross section of the Wkk -> WR model from the different anomaly detection methods

Exclusion limits on the production cross section of the G -> HH model from the different anomaly detection methods

Exclusion limits on the production cross section times acceptance times efficiency for generic Double Crystal Ball signal shape from the VAE-QR method

Exclusion limits on the production cross section times acceptance times efficiency for generic Double Crystal Ball signal shape from the CWoLa Hunting method

Exclusion limits on the production cross section times acceptance times efficiency for generic Double Crystal Ball signal shape from the TNT method

Exclusion limits on the production cross section times acceptance times efficiency for generic Double Crystal Ball signal shape from the CATHODE method

Exclusion limits on the production cross section times acceptance times efficiency for generic Double Crystal Ball signal shape from the CATHODE-b method

Exclusion limits on the production cross section times acceptance times efficiency for generic Double Crystal Ball signal shape from the QUAK - generic method

Number of events for the different lepton flavour categories in the signal region accounting for the fit to data. The hatched band includes all systematic uncertainties in the MC prediction. The vertical bars of the data account for the statistical uncertainty. The ratio panels show the ratio between data (black markers) with respect to the total prediction after the fit to data. Processes with a small contribution to this region are grouped in the ``Other" category

Sum of the charge of the final-state leptons in the signal region accounting for the fit to data. This distribution is not included in the fit. The hatched band includes all systematic uncertainties in the MC prediction. The vertical bars of the data account for the statistical uncertainty. The ratio panels show the ratio between data (black markers) with respect to the total prediction after the fit to data. Processes with a small contribution to this region are grouped in the ``Other" category.

Invariant mass of the trilepton system in the signal region accounting for the fit to data. This distribution is not included in the fit. The hatched band includes all systematic uncertainties in the MC prediction. The vertical bars of the data account for the statistical uncertainty. The ratio panels show the ratio between data (black markers) with respect to the total prediction after the fit to data. Processes with a small contribution to this region are grouped in the ``Other" category

Measured fiducial cross sections and their corresponding uncertainties for the flavourexclusive and flavour-inclusive categories. The predictions from both POWHEG at NLO in QCD and LO EW as well as several ones obtained from MATRIX (NLO QCD, NNLO QCD, NNLO QCD × NLO EW) are also included

Total WZ production cross section for each of the flavour-exclusive and for the flavour-inclusive categories. The vertical bands show different theoretical predictions for the WZ cross section at NLO in QCD (red dashed line) and NNLO QCD × NLO EW (blue solid line), as well as their corresponding scale uncertainties. For each measurement, the best fit value is denoted with a purple point, with two delimiters on the error bars that account for the total and statistical uncertainties

Observed 95% CL upper limits on the production cross section times branching ratio of the A → ZH → Ztt process in the (mA , mH) plane.

Median expected (dashed lines) and observed (filled contours) 95% CL exclusion regions in the (mH, mA) parameter space of the type-II 2HDM for tan β = 0.5 (blue), 1 (orange), and 1.5 (red), derived assuming narrow resonances. The enclosed regions are excluded. The dotted lines indicate points of constant natural line width of the A boson relative to its mass in the 2HDM.

Median expected (dashed lines) and observed (filled contours) 95% CL exclusion regions in the (mH, mA) parameter space of the type-II 2HDM for tan β = 0.5 (blue), 1 (orange), and 1.5 (red), derived assuming narrow resonances. The enclosed regions are excluded. The dotted lines indicate points of constant natural line width of the A boson relative to its mass in the 2HDM.

Median expected (dashed lines) and observed (filled contours) 95% CL exclusion regions in the (mH, mA) parameter space of the type-II 2HDM for tan β = 0.5 (blue), 1 (orange), and 1.5 (red), derived assuming narrow resonances. The enclosed regions are excluded. The dotted lines indicate points of constant natural line width of the A boson relative to its mass in the 2HDM.

Median expected (dashed lines) and observed (filled contours) 95% CL exclusion regions in the (mH, mA) parameter space of the type-II 2HDM for tan β = 0.5 (blue), 1 (orange), and 1.5 (red), derived assuming narrow resonances. The enclosed regions are excluded. The dotted lines indicate points of constant natural line width of the A boson relative to its mass in the 2HDM.

Median expected (dashed lines) and observed (filled contours) 95% CL exclusion regions in the (mH, mA) parameter space of the type-II 2HDM for tan β = 0.5 (blue), 1 (orange), and 1.5 (red), derived assuming narrow resonances. The enclosed regions are excluded. The dotted lines indicate points of constant natural line width of the A boson relative to its mass in the 2HDM.

Median expected (dashed lines) and observed (filled contours) 95% CL exclusion regions in the (mH, mA) parameter space of the type-II 2HDM for tan β = 0.5 (blue), 1 (orange), and 1.5 (red), derived assuming narrow resonances. The enclosed regions are excluded. The dotted lines indicate points of constant natural line width of the A boson relative to its mass in the 2HDM.

Median expected (dashed lines) and observed (filled contours) 95% CL exclusion regions in the type-II 2HDM (tan β, mA) parameter space at mH = 400 GeV, derived assuming narrow resonances. The enclosed regions are excluded. The dotted lines indicate points of constant natural line width of the A boson relative to its mass in the 2HDM.

Median expected (dashed lines) and observed (filled contours) 95% CL exclusion regions in the type-II 2HDM (tan β, mA) parameter space at mH = 400 GeV, derived assuming narrow resonances. The enclosed regions are excluded. The dotted lines indicate points of constant natural line width of the A boson relative to its mass in the 2HDM.

Median expected (dashed lines) and observed (filled contours) 95% CL exclusion regions in the type-II 2HDM (tanβ, cos(β−α)) parameter space at mA = 600 GeV and mH = 400 GeV, derived assuming narrow resonances. The enclosed regions are excluded. The dotted lines indicate points of constant natural line width of the A boson relative to its mass in the 2HDM.

Median expected (dashed lines) and observed (filled contours) 95% CL exclusion regions in the type-II 2HDM (tanβ, cos(β−α)) parameter space at mA = 600 GeV and mH = 400 GeV, derived assuming narrow resonances. The enclosed regions are excluded. The dotted lines indicate points of constant natural line width of the A boson relative to its mass in the 2HDM.

Invariant mass of the lepton pair plus jet system, for data, background (pre- and post-reweighting) and three $H\rightarrow Za$ signal hypotheses (for $a\rightarrow q\bar{q}/gg$ inclusively). Events are required to pass the complete event selection but not the classification NN requirement. The background normalization is set equal to that of the data for events passing the preselection and being in the $m_{\ell\ell j}$ 100-180 GeV region. The signal normalization assumes the SM Higgs boson inclusive production cross-section, $\mathcal{B}(H\to Za)=100\%$, and it is scaled up by a factor of 100. The error bars (hatched regions) represent the data (MC) sample's statistical uncertainty in the histograms and the ratio plots. Vertical arrows indicate data points that fall outside the displayed $y$-axis range.

Regression NN output variable, for data, reweighted background, and three $H\rightarrow Za$ signal hypotheses ($a\rightarrow q\bar{q}/gg$ inclusively). Events are required to pass the complete event selection, including a $120<m_{\ell\ell j}<140$ GeV requirement, but not the classification NN output variable requirement. The background normalization is set equal to that of the data for events passing the preselection and being in the $m_{\ell\ell j}$ 100-180 GeV region. The signal normalization assumes the SM Higgs boson inclusive production cross-section, $\mathcal{B}(H\to Za)=100\%$, and it is scaled up by a factor of 100. The error bars (hatched regions) represent the data (MC) sample's statistical uncertainty in the histograms and the ratio plots. Vertical arrows indicate data points that fall outside the displayed $y$-axis range.

Classification NN output variable, for data, reweighted background, and three $H\rightarrow Za$ signal hypotheses ($a\rightarrow q\bar{q}/gg$ inclusively). Events are required to pass the complete event selection, including a $120<m_{\ell\ell j}<140$ GeV requirement, but not the classification NN output variable requirement. The background normalization is set equal to that of the data for events passing the preselection and being in the $m_{\ell\ell j}$ 100-180 GeV region. The signal normalization assumes the SM Higgs boson inclusive production cross-section, $\mathcal{B}(H\to Za)=100\%$, and it is scaled up by a factor of 100. The error bars (hatched regions) represent the data (MC) sample's statistical uncertainty in the histograms and the ratio plots. Vertical arrows indicate data points that fall outside the displayed $y$-axis range.

Invariant mass of the lepton pair plus jet system for data, reweighted background, and various signal hypotheses. Events are required to pass the complete event selection (preselection and classification NN requirement). The pre-fit background distribution is shown along with three $H\rightarrow Za$ signal hypotheses coming from the pure MC sample ($a\rightarrow q\bar{q}/gg$ inclusively). The background normalization is set equal to that of the data for events passing the preselection and being in the region 100-180 GeV. The error bars represent the data sample's statistical uncertainty and the grey bands represent the full background model uncertainty before the fit (assuming all the uncertainties are uncorrelated). The lower panel shows the ratio of the data to the pre-fit background prediction.

Invariant mass of the lepton pair plus jet system for data, reweighted background, and 2.0 GeV signal hypotheses. Events are required to pass the complete event selection (preselection and classification NN requirement). The background shape and normalization come from the fit to the data. The signal includes only $gg$ decays, its shape comes from the fit, and it is normalized to $\mathcal{B}(H\rightarrow Za)=100\%$. The error bars represent the data sample's statistical uncertainty and the grey bands represent the full background model uncertainty after the fit. The lower panel shows the per-bin signal significance, $(N_\text{data}-N_\text{MC})\,\text{\Large{/}}\sqrt{\sigma_\text{data}^2+\sigma_\text{MC}^2}$. The post-fit background uncertainty is calculated from a fit under the background-only hypothesis. Vertical arrows indicate data points that fall outside the displayed $y$-axis range.

Observed 95$\%$ CL upper limits on $\mathcal{B}(H \to Za)=100\%$ for $\mathcal{B}(a \to gg)=100\%$, and the expectation under the background-only hypothesis together with its $\pm1\sigma$ and $\pm2\sigma$ intervals. The weaker limits from the previous version of the analysis (PRL 125(2020) 221802) are also shown. Linear interpolation is used to set limits between the mass points for which MC signal samples were generated.

Observed 95$\%$ confidence level upper limits on $\mathcal{B}(H \to Za)=100\%$ for $\mathcal{B}(a\rightarrow q\bar{q})=100\%$, and the expectation under the background-only hypothesis together with its $\pm1\sigma$ and $\pm2\sigma$ intervals. The weaker limits from the previous version of the analysis (PRL 125(2020) 221802) are also shown. Linear interpolation is used to set limits between the mass points for which MC signal samples were generated.

Observed two-dimensional likelihood scans for $c_{Hq}^{(1)}$ vs. $c_{Hq}^{(3)}$ while other coefficients are fixed at their SM values.

Observed two-dimensional likelihood scans for $c_{Hq}^{(1)}$ vs. $c_{Hu}$ while allowing the other coefficients to float freely at each point of the sca.

Observed two-dimensional likelihood scans for $c_{Hq}^{(1)}$ vs. $c_{Hu}$ while other coefficients are fixed at their SM values.

Observed two-dimensional likelihood scans for $c_{Hq}^{(1)}$ vs. $c_{Hd}$ while allowing the other coefficients to float freely at each point of the sca.

Observed two-dimensional likelihood scans for $c_{Hq}^{(1)}$ vs. $c_{Hd}$ while other coefficients are fixed at their SM values.

Observed two-dimensional likelihood scans for $c_{Hq}^{(3)}$ vs. $c_{Hu}$ while allowing the other coefficients to float freely at each point of the sca.

Observed two-dimensional likelihood scans for $c_{Hq}^{(3)}$ vs. $c_{Hu}$ while other coefficients are fixed at their SM values.

Observed two-dimensional likelihood scans for $c_{Hq}^{(3)}$ vs. $c_{Hd}$ while allowing the other coefficients to float freely at each point of the sca.

Observed two-dimensional likelihood scans for $c_{Hq}^{(3)}$ vs. $c_{Hd}$ while other coefficients are fixed at their SM values.

Observed two-dimensional likelihood scans for $c_{Hu}$ vs. $c_{Hd}$ while allowing the other coefficients to float freely at each point of the sca.

Observed two-dimensional likelihood scans for $c_{Hu}$ vs. $c_{Hd}$ while other coefficients are fixed at their SM values.

Observed two-dimensional likelihood scans for $c_{Hq}^{(1)}$ vs. $g^{(2)}_{ZZ}$ while allowing the other coefficients to float freely at each point of the sca.

Observed two-dimensional likelihood scans for $c_{Hq}^{(1)}$ vs. $g^{(2)}_{ZZ}$ while other coefficients are fixed at their SM values.

Observed two-dimensional likelihood scans for $c_{Hq}^{(3)}$ vs. $g^{(2)}_{ZZ}$ while allowing the other coefficients to float freely at each point of the sca.

Observed two-dimensional likelihood scans for $c_{Hq}^{(3)}$ vs. $g^{(2)}_{ZZ}$ while other coefficients are fixed at their SM values.

Observed two-dimensional likelihood scans for $c_{Hu}$ vs. $g^{(2)}_{ZZ}$ while allowing the other coefficients to float freely at each point of the sca.

Observed two-dimensional likelihood scans for $c_{Hu}$ vs. $g^{(2)}_{ZZ}$ while other coefficients are fixed at their SM values.

Observed two-dimensional likelihood scans for $c_{Hq}^{(1)}$ vs. $g^{(4)}_{ZZ}$ while allowing the other coefficients to float freely at each point of the sca.

Observed two-dimensional likelihood scans for $c_{Hq}^{(1)}$ vs. $g^{(4)}_{ZZ}$ while other coefficients are fixed at their SM values.

Observed two-dimensional likelihood scans for $c_{Hq}^{(3)}$ vs. $g^{(4)}_{ZZ}$ while allowing the other coefficients to float freely at each point of the sca.

Observed two-dimensional likelihood scans for $c_{Hq}^{(3)}$ vs. $g^{(4)}_{ZZ}$ while other coefficients are fixed at their SM values.

Observed two-dimensional likelihood scans for $c_{Hu}$ vs. $g^{(4)}_{ZZ}$ while allowing the other coefficients to float freely at each point of the sca.

Observed two-dimensional likelihood scans for $c_{Hu}$ vs. $g^{(4)}_{ZZ}$ while other coefficients are fixed at their SM values.

Observed two-dimensional likelihood scans for $c_{Hd}$ vs. $g^{(2)}_{ZZ}$ while allowing the other coefficients to float freely at each point of the sca.

Observed two-dimensional likelihood scans for $c_{Hd}$ vs. $g^{(2)}_{ZZ}$ while other coefficients are fixed at their SM values.

Observed two-dimensional likelihood scans for $c_{Hd}$ vs. $g^{(4)}_{ZZ}$ while allowing the other coefficients to float freely at each point of the sca.

Observed two-dimensional likelihood scans for $c_{Hd}$ vs. $g^{(4)}_{ZZ}$ while other coefficients are fixed at their SM values.

Observed two-dimensional likelihood scans for $g^{(2)}_{ZZ}$ vs. $g^{(4)}_{ZZ}$ while allowing the other coefficients to float freely at each point of the sca.

Observed two-dimensional likelihood scans for $g^{(2)}_{ZZ}$ vs. $g^{(4)}_{ZZ}$ while other coefficients are fixed at their SM values.

Invariant mass distribution of final state particles in SR1 HF category ($\text{H}\to\text{J}/\psi\gamma$ signal)

Invariant mass distribution of final state particles in SR1 Z-LP category ($\text{Z}\to\text{J}/\psi\gamma$ signal)

Invariant mass distribution of final state particles in SR1 Z-HP category ($\text{Z}\to\text{J}/\psi\gamma$ signal)

Invariant mass distribution of final state particles in SR2 category ($\text{H}\to\psi(\text{2S})\gamma$ signal)

Invariant mass distribution of final state particles in SR2 category ($\text{H}\to\psi(\text{2S})\gamma$ signal)

Invariant mass distribution of final state particles in control region category ($\text{H}\to\mu\mu\gamma$ background)

Observed and expected (with $\pm 1$, $\pm 2$ standard deviation bands) exclusion limits at 95% C.L. on the branching fraction $\mathcal{B}$ of the $(\text{H,Z}) \to \psi(\text{nS})\gamma$ decays

Observed (expected) upper limits at 95% C.L. on the normalized values with respect to the SM expectation, denoted as the signal strength parameter $\mu$, of the product of the cross section $\sigma$ and the branching fraction $\mathcal{B}$ of the $(\text{H,Z}) \to \psi(\text{nS})\gamma$ decays, and the branching fraction, assuming a SM Z and H boson cross section.

Results for the $P_3$ angular observable. The first uncertainties are statistical and the second systematic.

Results for the $P_4^\prime$ angular observable. The first uncertainties are statistical and the second systematic.

Results for the $P_5^\prime$ angular observable. The first uncertainties are statistical and the second systematic.

Results for the $P_6^\prime$ angular observable. The first uncertainties are statistical and the second systematic.

Results for the $P_8^\prime$ angular observable. The first uncertainties are statistical and the second systematic.

Results for the CP averaged observables $F_L$ and $P_1$–$P_8^\prime$. The first uncertainties are statistical and the second systematic.

Correlation matrix of angular observables Fl and P1–P8', considering only the statistical uncertainties from the maximum-likelihood fit in the interval 1.1 < $q^2$ < 2 GeV$^2$

Correlation matrix of angular observables Fl and P1–P8', considering only the statistical uncertainties from the maximum-likelihood fit in the interval 2 < $q^2$ < 4.3 GeV$^2$

Correlation matrix of angular observables Fl and P1–P8', considering only the statistical uncertainties from the maximum-likelihood fit in the interval 4.3 < $q^2$ < 6 GeV$^2$

Correlation matrix of angular observables Fl and P1–P8', considering only the statistical uncertainties from the maximum-likelihood fit in the interval 6 < $q^2$ < 8.68 GeV$^2$

Correlation matrix of angular observables Fl and P1–P8', considering only the statistical uncertainties from the maximum-likelihood fit in the interval 10.09 < $q^2$ < 12.86 GeV$^2$

Correlation matrix of angular observables Fl and P1–P8', considering only the statistical uncertainties from the maximum-likelihood fit in the interval 14.18 < $q^2$ < 16 GeV$^2$

Distributions in $m_{Z'}^{rec}$ for data in the SRs, together with fits of the background functions under the background-only hypothesis for the electron channel. The number of observed events in each bin is divided by the bin width. The signal predictions are shown for different Z' boson masses, normalized to an arbitrary cross section of 1 fb. In the panels below the distributions, the ratios of data to the background function are displayed. The shaded green areas represent the statistical uncertainty from the fit. The $\chi^2$ values per number of degrees of freedom ($\chi^2$/n.d.f.) and the corresponding $p$-values are provided for each fit.

Expected upper limits at 95% CL on the product of the production cross section and the branching fraction as a function of the Z' boson mass.Each line corresponds to the three different final states and their combined result.

Expected and observed upper limits at 95% CL on the product of the production cross section and the branching fraction as functions of the Z' boson mass from the combination of all final states. The limits are compared with the predictions of the HVT model and the expected limits (shown by red curves) from a previous analysis.

Prediction from the HVT model.

Observed upper limits at 95% CL on gF for different Z' boson masses as functions of the product of $g_{H}$ with the sign of $g_{F}$. The two benchmark scenarios of the HVT model are shown by the black markers.

Differential cross section of tZq as a function of the transverse momentum of the lepton coming from the W boson decay. The overflow is included in the last bin.

Differential cross section of ttZ+tWZ as a function of the azimuthal angle between the two leptons coming from the Z boson. The overflow is included in the last bin.

Differential cross section of tZq as a function of the azimuthal angle between the two leptons coming from the Z boson. The overflow is included in the last bin.

Differential cross section of ttZ+tWZ as a function of the $\Delta$R between the Z boson and the lepton coming from the W boson. The overflow is included in the last bin.

Differential cross section of tZq as a function of the $\Delta$R between the Z boson and the lepton coming from the W boson. The overflow is included in the last bin.

Differential cross section of ttZ+tWZ as a function of the cosine of the polar angle between the Z boson and the negatively charged lepton coming from its decay, boosted into the Z boson restframe.

Differential cross section of tZq as a function of the cosine of the polar angle between the Z boson and the negatively charged lepton coming from its decay, boosted into the Z boson restframe.

Normalized differential cross section of ttZ+tWZ as a function of the transverse momentum of the Z boson. The overflow is included in the last bin.

Normalized differential cross section of tZq as a function of the transverse momentum of the Z boson. The overflow is included in the last bin.

Normalized differential cross section of ttZ+tWZ as a function of the transverse momentum of the lepton coming from the W boson decay. The overflow is included in the last bin.

Normalized differential cross section of tZq as a function of the transverse momentum of the lepton coming from the W boson decay. The overflow is included in the last bin.

Normalized differential cross section of ttZ+tWZ as a function of the azimuthal angle between the two leptons coming from the Z boson. The overflow is included in the last bin.

Normalized differential cross section of tZq as a function of the azimuthal angle between the two leptons coming from the Z boson. The overflow is included in the last bin.

Normalized differential cross section of ttZ+tWZ as a function of the $\Delta$R between the Z boson and the lepton coming from the W boson. The overflow is included in the last bin.

Normalized differential cross section of tZq as a function of the $\Delta$R between the Z boson and the lepton coming from the W boson. The overflow is included in the last bin.

Normalized differential cross section of ttZ+tWZ as a function of the cosine of the polar angle between the Z boson and the negatively charged lepton coming from its decay, boosted into the Z boson restframe.

Normalized differential cross section of tZq as a function of the cosine of the polar angle between the Z boson and the negatively charged lepton coming from its decay, boosted into the Z boson restframe.

Covariance matrix for the differential cross sections as a function of the transverse momentum of the Z boson.

Covariance matrix for the differential cross sections as a function of the transverse momentum of the lepton coming from the W boson decay.

Covariance matrix for the differential cross sections as a function of the $\Delta$R between the Z boson and the lepton coming from the W boson.

Covariance matrix for the differential cross sections as a function of the azimuthal angle between the two leptons coming from the Z boson.

Covariance matrix for the differential cross sections as a function of the cosine of the polar angle between the Z boson and the negatively charged lepton coming from its decay, boosted into the Z boson restframe.