Measured $p_T$ cross sections in full-lepton phase space for 1.2 < |y| < 1.6.

Measured $p_T$ cross sections in full-lepton phase space for 1.6 < |y| < 2.0.

Measured $p_T$ cross sections in full-lepton phase space for 2.0 < |y| < 2.4.

Measured $p_T$ cross sections in full-lepton phase space for 2.4 < |y| < 2.8.

Measured $p_T$ cross sections in full-lepton phase space for 2.8 < |y| < 3.2.

Measured cross sections in full-lepton phase space as a function of |y|.

Normalised measured $p_T$ cross sections in full-lepton phase space for |y| < 1.6.

The total coherent $\mathrm{J}/\psi$ photoproduction cross section as a function of photon-nuclear center-of-mass energy per nucleon $W_{\gamma \mathrm{N}}^{\mathrm{Pb}}$, measured in PbPb ultra-peripheral collisions at $\sqrt{s_{\mathrm{NN}}}$ = 5.02 TeV. The $W_{\gamma \mathrm{N}}^{\mathrm{Pb}}$ values used correspond to the center of each rapidity range. The theoretical uncertainties is due to the uncertainties in the photon flux.

The nuclear gluon suppression factor $R_{\mathrm{g}}^{\mathrm{Pb}}$ as a function of Bjorken $x$ extracted from the CMS measurement of the coherent $\mathrm{J}/\psi$ photoproduction in PbPb ultra-peripheral collisions at $\sqrt{s_{\mathrm{NN}}}$ = 5.02 TeV. The $x$ values are evaluated at the centers of their corresponding rapidity ranges. The theoretical uncertainties are due to the uncertainties in the photon flux and the impulse approximation model.

The nuclear gluon suppression factor $R_{\mathrm{g}}^{\mathrm{Pb}}$ as a function of Bjorken $x$ extracted from the CMS measurement of the coherent $\mathrm{J}/\psi$ photoproduction in PbPb ultra-peripheral collisions at $\sqrt{s_{\mathrm{NN}}}$ = 5.02 TeV. The $x$ values are evaluated at the centers of their corresponding rapidity ranges. The theoretical uncertainties are due to the uncertainties in the photon flux and the impulse approximation model.

The total covariance matrix of the total coherent photoproduction cross section as a function of photon-nuclear center-of-mass energy per nucleon $W_{\gamma \mathrm{N}}^{\mathrm{Pb}}$. The covariance matrix includes both the experimental and theoretical (photon flux) uncertainties. The bins are ordered as increasing in $W_{\gamma \mathrm{N}}^{\mathrm{Pb}}$.

The photon flux values from STARLight as a function of rapidity, in different neutron multiplicity classes: 0n0n, 0nXn, and XnXn.

The experimental covariance matrix of the total coherent photoproduction cross section as a function of photon-nuclear center-of-mass energy per nucleon $W_{\gamma \mathrm{N}}^{\mathrm{Pb}}$. The bins are ordered as increasing in $W_{\gamma \mathrm{N}}^{\mathrm{Pb}}$.

The covariance matrix for the flux in the 0n0n neutron multiplicity class.

The theoretical (photon flux) covariance matrix of the total coherent photoproduction cross section as a function of photon-nuclear center-of-mass energy per nucleon $W_{\gamma \mathrm{N}}^{\mathrm{Pb}}$. The bins are ordered as increasing in $W_{\gamma \mathrm{N}}^{\mathrm{Pb}}$.

The covariance matrix for the flux in the 0nXn neutron multiplicity class.

The covariance matrix for the flux in the XnXn neutron multiplicity class.

The total covariance matrix of the total coherent photoproduction cross section as a function of photon-nuclear center-of-mass energy per nucleon $W_{\gamma \mathrm{N}}^{\mathrm{Pb}}$. The covariance matrix includes both the experimental and theoretical (photon flux) uncertainties. The bins are ordered as increasing in $W_{\gamma \mathrm{N}}^{\mathrm{Pb}}$.

The experimental covariance matrix of the total coherent photoproduction cross section as a function of photon-nuclear center-of-mass energy per nucleon $W_{\gamma \mathrm{N}}^{\mathrm{Pb}}$. The bins are ordered as increasing in $W_{\gamma \mathrm{N}}^{\mathrm{Pb}}$.

The theoretical (photon flux) covariance matrix of the total coherent photoproduction cross section as a function of photon-nuclear center-of-mass energy per nucleon $W_{\gamma \mathrm{N}}^{\mathrm{Pb}}$. The bins are ordered as increasing in $W_{\gamma \mathrm{N}}^{\mathrm{Pb}}$.

Measured cross sections for inclusive isolated-photon production as a function of $E_{\rm T}^{\gamma}$ for $1.56<|\eta^{\gamma}|<1.81$ and photon isolation cone radius $R=0.4$.

Measured cross sections for inclusive isolated-photon production as a function of $E_{\rm T}^{\gamma}$ for $1.81<|\eta^{\gamma}|<2.01$ and photon isolation cone radius $R=0.4$.

Measured cross sections for inclusive isolated-photon production as a function of $E_{\rm T}^{\gamma}$ for $2.01<|\eta^{\gamma}|<2.37$ and photon isolation cone radius $R=0.4$.

Measured cross sections for inclusive isolated-photon production as a function of $E_{\rm T}^{\gamma}$ for $|\eta^{\gamma}|<0.6$ and photon isolation cone radius $R=0.2$.

Measured cross sections for inclusive isolated-photon production as a function of $E_{\rm T}^{\gamma}$ for $0.6<|\eta^{\gamma}|<0.8$ and photon isolation cone radius $R=0.2$.

Measured cross sections for inclusive isolated-photon production as a function of $E_{\rm T}^{\gamma}$ for $0.8<|\eta^{\gamma}|<1.37$ and photon isolation cone radius $R=0.2$.

Measured cross sections for inclusive isolated-photon production as a function of $E_{\rm T}^{\gamma}$ for $1.56<|\eta^{\gamma}|<1.81$ and photon isolation cone radius $R=0.2$.

Measured cross sections for inclusive isolated-photon production as a function of $E_{\rm T}^{\gamma}$ for $1.81<|\eta^{\gamma}|<2.01$ and photon isolation cone radius $R=0.2$.

Measured cross sections for inclusive isolated-photon production as a function of $E_{\rm T}^{\gamma}$ for $2.01<|\eta^{\gamma}|<2.37$ and photon isolation cone radius $R=0.2$.

Predicted cross sections for inclusive isolated-photon production as a function of $E_{\rm T}^{\gamma}$ for $|\eta^{\gamma}|<0.6$ and isolation cone radius $0.4$ at NNLO QCD.

Predicted cross sections for inclusive isolated-photon production as a function of $E_{\rm T}^{\gamma}$ for $0.6<|\eta^{\gamma}|<0.8$ and isolation cone radius $0.4$ at NNLO QCD.

Predicted cross sections for inclusive isolated-photon production as a function of $E_{\rm T}^{\gamma}$ for $0.8<|\eta^{\gamma}|<1.37$ and isolation cone radius $0.4$ at NNLO QCD.

Predicted cross sections for inclusive isolated-photon production as a function of $E_{\rm T}^{\gamma}$ for $1.56<|\eta^{\gamma}|<1.81$ and isolation cone radius $0.4$ at NNLO QCD.

Predicted cross sections for inclusive isolated-photon production as a function of $E_{\rm T}^{\gamma}$ for $1.81<|\eta^{\gamma}|<2.01$ and isolation cone radius $0.4$ at NNLO QCD.

Predicted cross sections for inclusive isolated-photon production as a function of $E_{\rm T}^{\gamma}$ for $2.01<|\eta^{\gamma}|<2.37$ and isolation cone radius $0.4$ at NNLO QCD.

Predicted cross sections for inclusive isolated-photon production as a function of $E_{\rm T}^{\gamma}$ for $|\eta^{\gamma}|<0.6$ and isolation cone radius $0.2$ at NNLO QCD.

Predicted cross sections for inclusive isolated-photon production as a function of $E_{\rm T}^{\gamma}$ for $0.6<|\eta^{\gamma}|<0.8$ and isolation cone radius $0.2$ at NNLO QCD.

Predicted cross sections for inclusive isolated-photon production as a function of $E_{\rm T}^{\gamma}$ for $0.8<|\eta^{\gamma}|<1.37$ and isolation cone radius $0.2$ at NNLO QCD.

Predicted cross sections for inclusive isolated-photon production as a function of $E_{\rm T}^{\gamma}$ for $1.56<|\eta^{\gamma}|<1.81$ and isolation cone radius $0.2$ at NNLO QCD.

Predicted cross sections for inclusive isolated-photon production as a function of $E_{\rm T}^{\gamma}$ for $1.81<|\eta^{\gamma}|<2.01$ and isolation cone radius $0.2$ at NNLO QCD.

Predicted cross sections for inclusive isolated-photon production as a function of $E_{\rm T}^{\gamma}$ for $2.01<|\eta^{\gamma}|<2.37$ and isolation cone radius $0.2$ at NNLO QCD.

Measured ratio of the differential cross sections for inclusive isolated-photon production for $R=0.2$ and $R=0.4$ as a function of $E_{\rm T}^{\gamma}$ for $|\eta^{\gamma}|<0.6$.

Measured ratio of the differential cross sections for inclusive isolated-photon production for $R=0.2$ and $R=0.4$ as a function of $E_{\rm T}^{\gamma}$ for $0.6<|\eta^{\gamma}|<0.8$.

Measured ratio of the differential cross sections for inclusive isolated-photon production for $R=0.2$ and $R=0.4$ as a function of $E_{\rm T}^{\gamma}$ for $0.8<|\eta^{\gamma}|<1.37$.

Measured ratio of the differential cross sections for inclusive isolated-photon production for $R=0.2$ and $R=0.4$ as a function of $E_{\rm T}^{\gamma}$ for $1.56<|\eta^{\gamma}|<1.81$.

Measured ratio of the differential cross sections for inclusive isolated-photon production for $R=0.2$ and $R=0.4$ as a function of $E_{\rm T}^{\gamma}$ for $1.81<|\eta^{\gamma}|<2.01$.

Measured ratio of the differential cross sections for inclusive isolated-photon production for $R=0.2$ and $R=0.4$ as a function of $E_{\rm T}^{\gamma}$ for $2.01<|\eta^{\gamma}|<2.37$.

Predicted ratio of the differential cross sections for inclusive isolated-photon production for $R=0.2$ and $R=0.4$ as a function of $E_{\rm T}^{\gamma}$ for $|\eta^{\gamma}|<0.6$ at NNLO QCD.

Predicted ratio of the differential cross sections for inclusive isolated-photon production for $R=0.2$ and $R=0.4$ as a function of $E_{\rm T}^{\gamma}$ for $0.6<|\eta^{\gamma}|<0.8$ at NNLO QCD.

Predicted ratio of the differential cross sections for inclusive isolated-photon production for $R=0.2$ and $R=0.4$ as a function of $E_{\rm T}^{\gamma}$ for $0.8<|\eta^{\gamma}|<1.37$ at NNLO QCD.

Predicted ratio of the differential cross sections for inclusive isolated-photon production for $R=0.2$ and $R=0.4$ as a function of $E_{\rm T}^{\gamma}$ for $1.56<|\eta^{\gamma}|<1.81$ at NNLO QCD.

Predicted ratio of the differential cross sections for inclusive isolated-photon production for $R=0.2$ and $R=0.4$ as a function of $E_{\rm T}^{\gamma}$ for $1.81<|\eta^{\gamma}|<2.01$ at NNLO QCD.

Predicted ratio of the differential cross sections for inclusive isolated-photon production for $R=0.2$ and $R=0.4$ as a function of $E_{\rm T}^{\gamma}$ for $2.01<|\eta^{\gamma}|<2.37$ at NNLO QCD.

Measured fiducial integrated cross section for inclusive isolated-photon production as a function of $R$ for $|\eta^{\gamma}|<0.6$.

Measured fiducial integrated cross section for inclusive isolated-photon production as a function of $R$ for $0.6<|\eta^{\gamma}|<0.8$.

Measured fiducial integrated cross section for inclusive isolated-photon production as a function of $R$ for $0.8<|\eta^{\gamma}|<1.37$.

Measured fiducial integrated cross section for inclusive isolated-photon production as a function of $R$ for $1.56<|\eta^{\gamma}|<1.81$.

Measured fiducial integrated cross section for inclusive isolated-photon production as a function of $R$ for $1.81<|\eta^{\gamma}|<2.01$.

Measured fiducial integrated cross section for inclusive isolated-photon production as a function of $R$ for $2.01<|\eta^{\gamma}|<2.37$.

Predicted fiducial integrated cross section for inclusive isolated-photon production as a function of $R$ for $|\eta^{\gamma}|<0.6$ at NNLO QCD.

Predicted fiducial integrated cross section for inclusive isolated-photon production as a function of $R$ for $0.6<|\eta^{\gamma}|<0.8$ at NNLO QCD.

Predicted fiducial integrated cross section for inclusive isolated-photon production as a function of $R$ for $0.8<|\eta^{\gamma}|<1.37$ at NNLO QCD.

Predicted fiducial integrated cross section for inclusive isolated-photon production as a function of $R$ for $1.56<|\eta^{\gamma}|<1.81$ at NNLO QCD.

Predicted fiducial integrated cross section for inclusive isolated-photon production as a function of $R$ for $1.81<|\eta^{\gamma}|<2.01$ at NNLO QCD.

Predicted fiducial integrated cross section for inclusive isolated-photon production as a function of $R$ for $2.01<|\eta^{\gamma}|<2.37$ at NNLO QCD.

Measured unfolded differential cross-section as a function of the dilepton transverse momentum $p^{ll}_{\mathrm{T}}$. NLO predictions from Sherpa 2.2.10 and MadGraph5_aMC@NLO 2.7.3 are also shown. The uncertainty in the predictions is divided into statistical and theoretical uncertainties (scale and PDF+$\alpha_{s}$).

Measured unfolded differential cross-section as a function of the the four-body transverse momentum $p^{ll\gamma\gamma}_{\mathrm{T}}$. NLO predictions from Sherpa 2.2.10 and MadGraph5_aMC@NLO 2.7.3 are also shown. The uncertainty in the predictions is divided into statistical and theoretical uncertainties (scale and PDF+$\alpha_{s}$).

Measured unfolded differential cross-section as a function of the diphoton invariant mass $m_{\gamma\gamma}$. NLO predictions from Sherpa 2.2.10 and MadGraph5_aMC@NLO 2.7.3 are also shown. The uncertainty in the predictions is divided into statistical and theoretical uncertainties (scale and PDF+$\alpha_{s}$).

Measured unfolded differential cross-section as a function of the four-body invariant mass $m_{ll\gamma\gamma}$. NLO predictions from Sherpa 2.2.10 and MadGraph5_aMC@NLO 2.7.3 are also shown. The uncertainty in the predictions is divided into statistical and theoretical uncertainties (scale and PDF+$\alpha_{s}$).

Expected and observed $95\%$ confidence intervals for the coupling parameters $f_{T,j}/\Lambda^{4}$ of transverse dimension-8 operators. All parameter values outside of the stated range are excluded at the chosen confidence level. No unitarity constraints are applied.

Expected and observed unitarised $95\%$ confidence intervals for the coupling parameter $f_{T,8}/\Lambda^{4}$ in the clipping energy range between 1.1 and 5 TeV. The non-unitarised limits ($E_c = \infty$) are also shown. All parameter values outside of the stated range are excluded at the chosen confidence level.

Expected and observed unitarised $95\%$ confidence intervals for the coupling parameter $f_{T,0}/\Lambda^{4}$ in the clipping energy range between 1.1 and 5 TeV. The non-unitarised limits ($E_c = \infty$) are also shown. All parameter values outside of the stated range are excluded at the chosen confidence level.

Expected and observed unitarised $95\%$ confidence intervals for the coupling parameter $f_{T,1}/\Lambda^{4}$ in the clipping energy range between 1.1 and 5 TeV. The non-unitarised limits ($E_c = \infty$) are also shown. All parameter values outside of the stated range are excluded at the chosen confidence level.

Expected and observed unitarised $95\%$ confidence intervals for the coupling parameter $f_{T,2}/\Lambda^{4}$ in the clipping energy range between 1.1 and 5 TeV. The non-unitarised limits ($E_c = \infty$) are also shown. All parameter values outside of the stated range are excluded at the chosen confidence level.

Expected and observed unitarised $95\%$ confidence intervals for the coupling parameter $f_{T,5}/\Lambda^{4}$ in the clipping energy range between 1.1 and 5 TeV. The non-unitarised limits ($E_c = \infty$) are also shown. All parameter values outside of the stated range are excluded at the chosen confidence level.

Expected and observed unitarised $95\%$ confidence intervals for the coupling parameter $f_{T,6}/\Lambda^{4}$ in the clipping energy range between 1.1 and 5 TeV. The non-unitarised limits ($E_c = \infty$) are also shown. All parameter values outside of the stated range are excluded at the chosen confidence level.

Expected and observed unitarised $95\%$ confidence intervals for the coupling parameter $f_{T,7}/\Lambda^{4}$ in the clipping energy range between 1.1 and 5 TeV. The non-unitarised limits ($E_c = \infty$) are also shown. All parameter values outside of the stated range are excluded at the chosen confidence level.

Expected and observed unitarised $95\%$ confidence intervals for the coupling parameter $f_{T,9}/\Lambda^{4}$ in the clipping energy range between 1.1 and 5 TeV. The non-unitarised limits ($E_c = \infty$) are also shown. All parameter values outside of the stated range are excluded at the chosen confidence level.

Distribution of ${E}_{T}^{miss}$ after the fit of the multijet backgrounds, in the electron channel, in the $t\bar{t}$ 3-jets VR, without applying the cut on ${E}_{T}^{miss}$. Simulated events are normalised to the expected number of events given the integrated luminosity, after applying the normalisation factors obtained in the multijet fit. The last bin includes the overflow. The uncertainty band indicates the simulation's statistical uncertainty, the normalisation uncertainties for different processes ($40$ % for $W$+jets production, $30$ % for multijet background and $6$ % for top-quark processes) and the multijet background shape uncertainty in each bin, summed in quadrature. The lower panel of the figure shows the ratio of the data to the prediction.

Distribution of ${E}_{T}^{miss}$ after the fit of the multijet backgrounds, in the electron channel, in the $t\bar{t}$ 4-jets VR, without applying the cut on ${E}_{T}^{miss}$. Simulated events are normalised to the expected number of events given the integrated luminosity, after applying the normalisation factors obtained in the multijet fit. The last bin includes the overflow. The uncertainty band indicates the simulation's statistical uncertainty, the normalisation uncertainties for different processes ($40$ % for $W$+jets production, $30$ % for multijet background and $6$ % for top-quark processes) and the multijet background shape uncertainty in each bin, summed in quadrature. The lower panel of the figure shows the ratio of the data to the prediction.

Distribution of $m_{T}^{W}$ after the fit of the multijet backgrounds, in the muon channel, in the signal region, without applying the cut on $m_{T}^{W}$. Simulated events are normalised to the expected number of events given the integrated luminosity, after applying the normalisation factors obtained in the multijet fit. The last bin includes the overflow. The uncertainty band indicates the simulation's statistical uncertainty, the normalisation uncertainties for different processes ($40$ % for $W$+jets production, $30$ % for multijet background and $6$ % for top-quark processes) and the multijet background shape uncertainty in each bin, summed in quadrature. The lower panel of the figure shows the ratio of the data to the prediction.

Distribution of $m_{T}^{W}$ after the fit of the multijet backgrounds, in the muon channel, in the $W$+jets VR, without applying the cut on $m_{T}^{W}$. Simulated events are normalised to the expected number of events given the integrated luminosity, after applying the normalisation factors obtained in the multijet fit. The last bin includes the overflow. The uncertainty band indicates the simulation's statistical uncertainty, the normalisation uncertainties for different processes ($40$ % for $W$+jets production, $30$ % for multijet background and $6$ % for top-quark processes) and the multijet background shape uncertainty in each bin, summed in quadrature. The lower panel of the figure shows the ratio of the data to the prediction.

Distribution of $m_{T}^{W}$ after the fit of the multijet backgrounds, in the muon channel, in the $t\bar{t}$ 3-jets VR, without applying the cut on $m_{T}^{W}$. Simulated events are normalised to the expected number of events given the integrated luminosity, after applying the normalisation factors obtained in the multijet fit. The last bin includes the overflow. The uncertainty band indicates the simulation's statistical uncertainty, the normalisation uncertainties for different processes ($40$ % for $W$+jets production, $30$ % for multijet background and $6$ % for top-quark processes) and the multijet background shape uncertainty in each bin, summed in quadrature. The lower panel of the figure shows the ratio of the data to the prediction.

Distribution of $m_{T}^{W}$ after the fit of the multijet backgrounds, in the muon channel, in the $t\bar{t}$ 4-jets VR, without applying the cut on $m_{T}^{W}$. Simulated events are normalised to the expected number of events given the integrated luminosity, after applying the normalisation factors obtained in the multijet fit. The last bin includes the overflow. The uncertainty band indicates the simulation's statistical uncertainty, the normalisation uncertainties for different processes ($40$ % for $W$+jets production, $30$ % for multijet background and $6$ % for top-quark processes) and the multijet background shape uncertainty in each bin, summed in quadrature. The lower panel of the figure shows the ratio of the data to the prediction.

Expected distributions of the MEM discriminant $P(S|X)$ in the SR, for the s-channel single-top signal, and for the $t\bar{t}$ and $W$+jets backgrounds, for MEM discriminant values larger than $2.0\times10^{-4}$. Each distribution is normalised to unity. The binning is the same as the optimised binning used in the signal extraction fit, resulting in a non-linear horizontal scale.

Distribution of the MEM discriminant $P(S|X)$ in the $W$+jets VR. Simulated events are normalised to the expected number of events given the integrated luminosity, after applying the normalisation factors obtained in the multijet fit presented in Section 5 in the paper. The uncertainty band indicates the simulation's statistical uncertainty and the normalisation uncertainties for the various processes in each bin, summed in quadrature. The ratio of the observed number to the predicted number of events in each bin is shown in the lower panel of the figure, with different vertical axis ranges. The binning is the same as the optimised binning used in the signal extraction fit described in Section 8 in the paper, resulting in a non-linear horizontal scale.

Distribution of the MEM discriminant $P(S|X)$ in the $t\bar{t}$ 3-jets VR. Simulated events are normalised to the expected number of events given the integrated luminosity, after applying the normalisation factors obtained in the multijet fit presented in Section 5 in the paper. The uncertainty band indicates the simulation's statistical uncertainty and the normalisation uncertainties for the various processes in each bin, summed in quadrature. The ratio of the observed number to the predicted number of events in each bin is shown in the lower panel of the figure, with different vertical axis ranges. The binning is the same as the optimised binning used in the signal extraction fit described in Section 8 in the paper, resulting in a non-linear horizontal scale.

Distribution of the MEM discriminant $P(S|X)$ in the $t\bar{t}$ 4-jets VR. Simulated events are normalised to the expected number of events given the integrated luminosity, after applying the normalisation factors obtained in the multijet fit presented in Section 5 in the paper. The uncertainty bands indicate the simulation's statistical uncertainty and the normalisation uncertainties for the various processes in each bin, summed in quadrature. The ratio of the observed number to the predicted number of events in each bin is shown in the lower panel of the figure, with different vertical axis ranges. The binning is the same as the optimised binning used in the signal extraction fit described in Section 8 in the paper, resulting in a non-linear horizontal scale.

Distribution of the MEM discriminant $P(S|X)$ in the SR before the fit to data, for MEM discriminant values larger than $2.0\times10^{-4}$. The lower panel of the figure shows the ratio of the data to the prediction, with different vertical axis ranges. The uncertainty band indicates the total uncertainties and their correlations in each bin. The uncertainties in the $t\bar{t}$ and $W$+jets normalisation factors, as well as in the s-channel signal cross-section, are not defined pre-fit and therefore not included. The binning is the same as the optimised binning used in the fit, resulting in a non-linear horizontal scale.

Distribution of the MEM discriminant $P(S|X)$ in the SR after the fit to data, for MEM discriminant values larger than $2.0\times10^{-4}$. The lower panel of the figure shows the ratio of the data to the prediction, with different vertical axis ranges. The uncertainty band indicates the total uncertainties and their correlations in each bin. The binning is the same as the optimised binning used in the fit, resulting in a non-linear horizontal scale.

Distribution of the MEM discriminant $P(S|X)$ in the SR after the fit to data, for MEM discriminant values larger than $2.0\times10^{-4}$, after subtraction of all backgrounds. The fitted distribution for the simulation of the signal is shown together with the post-fit uncertainty in the backgrounds. The binning is the same as the optimised binning used in the fit, resulting in a non-linear horizontal scale.

Pre-fit and post-fit event yields in the SR, for MEM discriminant values larger than $2.0\times10^{-4}$. The central value of the event yield for each process is calculated by summing the values of the discriminant bin contents, using the nominal expected yield for the pre-fit value, and the best-fit estimate for the post-fit value. The error includes statistical and systematic uncertainties summed in quadrature. All sources of systematic uncertainties are included, taking into account correlations and anti-correlations in the post-fit case. The uncertainties in the $t\bar{t}$ and $W$+jets normalisation factors, as well as in the s-channel signal cross-section, are not defined pre-fit and therefore only included in the post-fit uncertainties.

Observed impact of the different sources of uncertainty on the measured s-channel signal cross-section, grouped by categories. The impact of each category is obtained by repeating the fit after having fixed the set of nuisance parameters corresponding to that category, subtracting the square of the resulting uncertainty from the square of the uncertainty found in the full fit, and calculating the square root. The 'Systematic uncertainties' category combines all sources of systematic uncertainties. The statistical uncertainty is obtained by repeating the fit after having fixed all nuisance parameters, including the $t\bar{t}$ and $W$+jets normalisation factors. 'Total' gives the total uncertainty on the measurement.

Observed impact of the different sources of $t\bar{t}$ modelling uncertainty on the measured s-channel signal cross-section. The impact of each category is obtained by repeating the fit after having fixed the set of nuisance parameters corresponding to that category, subtracting the square of the resulting uncertainty from the square of the uncertainty found in the full fit, and calculating the square root. 'PS & had.' refers to the parton shower and hadronisation model, and 'ME/PS matching' to the matching of the ME to the parton shower.

Observed impact of the different sources of s-channel modelling uncertainty on the measured s-channel signal cross-section. The impact of each category is obtained by repeating the fit after having fixed the set of nuisance parameters corresponding to that category, subtracting the square of the resulting uncertainty from the square of the uncertainty found in the full fit, and calculating the square root. 'PS & had.' refers to the parton shower and hadronisation model, as described in Section 7 in the paper.

Observed impact of the different sources of t-channel modelling uncertainty on the measured s-channel signal cross-section. The impact of each category is obtained by repeating the fit after having fixed the set of nuisance parameters corresponding to that category, subtracting the square of the resulting uncertainty from the square of the uncertainty found in the full fit, and calculating the square root. 'PS & had.' refers to the parton shower and hadronisation model, as described in Section 7 in the paper.

Observed impact of the different sources of $tW$ modelling uncertainty on the measured s-channel signal cross-section, grouped by categories. The impact of each category is obtained by repeating the fit after having fixed the set of nuisance parameters corresponding to that category, subtracting the square of the resulting uncertainty from the square of the uncertainty found in the full fit, and calculating the square root. 'PS & had.' refers to the parton shower and hadronisation model, and '$t\bar{t}$ overlap' to the algorithm removing the overlap between $tW$ and $t\bar{t}$ production at NLO, as described in Section 7 in the paper.

Observed impact of the different sources of PDF uncertainties on the measured s-channel signal cross-section, grouped by categories. The impact of each category is obtained by repeating the fit after having fixed the set of nuisance parameters corresponding to that category, subtracting the square of the resulting uncertainty from the square of the uncertainty found in the full fit, and calculating the square root.

Comparison between data and prediction after the fit to data in the signal region for the leading-jet $p_{T}$. The last bin includes the overflow. The uncertainty band includes all uncertainties and their correlations. The lower panel of the figure shows the ratio of the data to the prediction.

Comparison between data and prediction after the fit to data in the signal region for the leading-jet $\eta$. The uncertainty band includes all uncertainties and their correlations. The lower panel of the figure shows the ratio of the data to the prediction.

Comparison between data and prediction after the fit to data in the signal region for the subleading-jet $p_{T}$. The last bin includes the overflow. The uncertainty band includes all uncertainties and their correlations. The lower panel of the figure shows the ratio of the data to the prediction.

Comparison between data and prediction after the fit to data in the signal region for the subleading-jet $\eta$. The uncertainty band includes all uncertainties and their correlations. The lower panel of the figure shows the ratio of the data to the prediction.

Comparison between data and prediction after the fit to data in the signal region for the lepton $p_{T}$. The last bin includes the overflow. The uncertainty band includes all uncertainties and their correlations. The lower panel of the figure shows the ratio of the data to the prediction.

Comparison between data and prediction after the fit to data in the signal region for the lepton $\eta$. The uncertainty band includes all uncertainties and their correlations. The lower panel of the figure shows the ratio of the data to the prediction.

Comparison between data and prediction after the fit to data in the signal region for the ${E}_{T}^{miss}$. The last bin includes the overflow. The uncertainty band includes all uncertainties and their correlations. The lower panel of the figure shows the ratio of the data to the prediction.

Comparison between data and prediction after the fit to data in the signal region for the $m_{T}^{W}$. The last bin includes the overflow. The uncertainty band includes all uncertainties and their correlations. The lower panel of the figure shows the ratio of the data to the prediction.

Nuisance parameters ranked according to their post-fit impacts on the best-fit value of the ratio $\mu$ of the measured cross-section to the predicted cross-section. In the figure, only the 20 nuisance parameters with the largest post-fit impacts are shown. The empty (solid) blue rectangles illustrate the pre-fit (post-fit) impact on $\mu$, corresponding to the upper axis. The pre-fit (post-fit) impact of each nuisance parameter, $\Delta\mu$, is calculated as the difference in the fitted value of $\mu$ between the nominal fit and the fit when fixing the corresponding nuisance parameter to $\hat{\theta}\pm\Delta\theta$ ($\hat{\theta}\pm\Delta\hat{\theta}$), where $\hat{\theta}$ is the best-fit value of the nuisance parameter and $\Delta\theta$ ($\Delta\hat{\theta}$) is its pre-fit (post-fit) uncertainty. Several systematic uncertainties are split into different nuisance parameters, which are indicated by NP. JES (JER) indicates jet energy scale (resolution), and $\gamma$ indicates a nuisance parameter associated to the MC statistics in one of the 18 bins numbered from 0 to 17. The black points show the best-fit values of the nuisance parameters, with the error bars representing the post-fit uncertainties. Each nuisance parameter is shown wrt. its nominal value, $\theta_0$, and in units of its pre-fit uncertainty, except the free-floating normalisation factors of the $t\bar{t}$ and $W$+jets backgrounds, and the parameters associated to the MC statistics in each bin, for which the post-fit values and uncertainties are shown.

Correlation matrix of the nuisance parameters and of the ratio $\mu$ of the measured cross-section to the predicted cross-section. The correlations are given after the fit to data. In the figure, only the parameters which have a correlation of at least 0.2 with any other parameter are shown.

Distribution of the MEM discriminant $P(S|X)$ in the SR for MEM discriminant values larger than $2.0\times10^{-4}$, for the collision data used for the measurement, and for 1000 pseudo-data replicas, generated using a bootstrapping technique, in order to assess the statistical correlations between this measurement and others, for the purpose of combinations. The replicas are obtained by reweighting each observed data event by a random integer generated according to Poisson statistics, using the <a href="https://zenodo.org/record/5361038">BootstrapGenerator</a> software package , which implements a technique described in <a href="https://cds.cern.ch/record/2759945/">ATL-PHYS-PUB-2021-011</a>. The ATLAS event number and run number of each event are used as seed to uniquely but reproducibly initialise the random number generator for each event. Each pseudo-data replica is assigned an index, ranging from 0 to 999, corresponding to the random number index used consistently for each observed data event.

Measured values of the signal cross-section and of the $t\bar{t}$ and $W$+jets normalisation factors, obtained by statistical-only fits to the collision data used for the measurement, and to 1000 pseudo-data replicas, generated using a bootstrapping technique, in order to assess the statistical correlations between this measurement and others, for the purpose of combinations. The central values and their statistical uncertainties are obtained by repeating the fit after having fixed all nuisance parameters, except the $t\bar{t}$ and $W$+jets normalisation factors, which are let free-floating (unlike for the statistical uncertainty on the cross-section quoted in the paper). The replicas are obtained by reweighting each observed data event by a random integer generated according to Poisson statistics, using the <a href="https://zenodo.org/record/5361038">BootstrapGenerator</a> software package , which implements a technique described in <a href="https://cds.cern.ch/record/2759945/">ATL-PHYS-PUB-2021-011</a>. The ATLAS event number and run number of each event are used as seed to uniquely but reproducibly initialise the random number generator for each event. Each pseudo-data replica is assigned an index, ranging from 0 to 999, corresponding to the random number index used consistently for each observed data event.

The rho-parameter, i.e. the ratio of the real to imaginary part of the elastic scattering amplitude extrapolated to t=0. The systematic uncertainty includes experimental and theoretical uncerainties.

The nuclear slope parameter B from a fit of the form exp(-Bt-Ct^2-Dt^3). The systematic uncertainty includes experimental and theoretical uncerainties.

The nuclear slope parameter B from a fit of the form exp(-Bt-Ct^2-Dt^3). The systematic uncertainty includes experimental and theoretical uncerainties.

The nuclear slope parameter C from a fit of the form exp(-Bt-Ct^2-Dt^3). The systematic uncertainty includes experimental and theoretical uncerainties.

The nuclear slope parameter C from a fit of the form exp(-Bt-Ct^2-Dt^3). The systematic uncertainty includes experimental and theoretical uncerainties.

The nuclear slope parameter D from a fit of the form exp(-Bt-Ct^2-Dt^3). The systematic uncertainty includes experimental and theoretical uncerainties.

The nuclear slope parameter D from a fit of the form exp(-Bt-Ct^2-Dt^3). The systematic uncertainty includes experimental and theoretical uncerainties.

The total elastic cross section measured inside the fiducial volume. The systematic uncertainty includes experimental uncertainties.

The total elastic cross section measured inside the fiducial volume. The systematic uncertainty includes experimental uncertainties.

The total elastic cross section obtained from the fitted parameters, extrapolated to full phase space using only the nuclear amplitude.

The total elastic cross section obtained from the fitted parameters, extrapolated to full phase space using only the nuclear amplitude.

The total inelastic cross section.

The total inelastic cross section.

The ratio of elastic to total cross section.

The ratio of elastic to total cross section.

The differential elastic cross section as function of t with statistical and systematic uncertainties. The systematic uncertainties are given as signed relative change for 20 sources of experimental uncertainty associated to nuisance parameters used in the fit for the extraction of physics parameters.

The differential elastic cross section as function of t with statistical and systematic uncertainties. The systematic uncertainties are given as signed relative change for 20 sources of experimental uncertainty associated to nuisance parameters used in the fit for the extraction of physics parameters.

Statistical covariance matrix for the measurement of the differential elastic cross section as function of t.

Statistical covariance matrix for the measurement of the differential elastic cross section as function of t.

Correlation matrix for the measurement of STXS parameters in each of the 28 regions considered.

Fitted values for kappa_g and kappa_y.

Correlation matrix for the measurement of kappa_g and kappa_y.

Summary of the 68% CL confidence intervals for individual measurements of SMEFT parameters observed in data. In each case, SMEFT parameters other than the one measured are fixed to 0.

Results of the EV parameter measurement in data, in the linear and linear+quadratic parameterizations of the SMEFT. All the EVs parameters are free to vary in the fits. The ranges correspond to 68% CL confidence intervals.

Observed linear correlation coefficients of the EVs parameters in the linear parameterization.

Observed linear correlation coefficients of the EVs parameters in the linear+quadratic parameterization.

Cross-sections times H->yy branching ratio for ggF +bbH, VBF, VH, ttH, and tH production. The values are obtained from a simultaneous fit to all categories. The theory uncertainties in the predictions include uncertainties due to missing higher-order terms in the perturbative QCD calculations and choices of parton distribution functions and value of alpha_s, as well as the H->yy branching ratio uncertainty.

Best-fit values and uncertainties for STXS parameters in each of the 28 regions considered. The values for the gg->H process also include the contributions from bbH production.

Correlation matrix for the measurement of STXS parameters in each of the 33 regions considered.

Distribution of an input to the BDT classifier in the $3\ell$ category: The linear discrimminant $p_{T}^{miss,LD}$. The discriminant is defined by the relation $p_{T}^{miss,LD} = 0.6 p_{T}^{miss} + 0.4 H_{T}^{miss}$, where $H_{T}^{miss}$ corresponds to the magnitude of the vector $p_{T}$ sum of all electrons, muons, taus and small radius jets passing the selection critertia, and $p_{T}^{miss}$ represents the missing transverse momentum vector computed as the negative vector $p_{T}$ sum of all the particles reconstructed by the PF algorithm in the event.

Distributions of the transverse mass $m_{T}$ of the lepton that is not originating from the Z boson decay and the missing transverse momentum in the $3\ell$/WZ control region. The distributions expected for the WZ and ZZ as well as for other background processes are shown for the values of nuisance parameters obtained from the ML fit used in the signal extraction.

Distributions of the four lepton invariant mass $m_{4\ell}$ in the $4\ell$/ZZ control region. The distributions expected for the ZZ as well as for other background processes are shown for the values of nuisance parameters obtained from the ML fit used in the signal extraction.

Distributions of the transverse mass $m_{T}$ of the leading lepton and the missing transverse momentum in the $2\ell$(ss) control region. The distributions expected for the misidentified $\ell$ background as well as for other background processes is shown for the values of nuisance parameters obtained from the background-only ML fit, in which the HH signal is constrained to be zero.

Distributions of the reconstructed HH mass $m_{HH}$ in the $2\ell+2 au_{h}$ control region. The distributions expected for the misidentified $\ell/ au_{h}$ background as well as for other background processes is shown for the values of nuisance parameters obtained from the background-only ML fit, in which the HH signal is constrained to be zero.

Distribution in the output of the BDT trained for nonresonant HH production and evaluated for the benchmark scenario JHEP04 BM7 for the $2\ell$(ss) category. The SM HH signal is shown for a cross section amounting to $30$ times the value predicted in the SM. The distributions expected for the background processes are shown for the values of nuisance parameters obtained from the ML fit of the signal+background hypothesis to the data. The binning is chosen to approximate signal quantiles.

Distribution in the output of the BDT trained for nonresonant HH production and evaluated for the benchmark scenario JHEP04 BM7 for the $3\ell$ category. The SM HH signal is shown for a cross section amounting to $30$ times the value predicted in the SM. The distributions expected for the background processes are shown for the values of nuisance parameters obtained from the ML fit of the signal+background hypothesis to the data. The binning is chosen to approximate signal quantiles.

Distribution in the output of the BDT trained for nonresonant HH production and evaluated for the benchmark scenario JHEP04 BM7 for the $4\ell$ category. The SM HH signal is shown for a cross section amounting to $30$ times the value predicted in the SM. The distributions expected for the background processes are shown for the values of nuisance parameters obtained from the ML fit of the signal+background hypothesis to the data. The binning is chosen to approximate signal quantiles.

Distribution in the output of the BDT trained for nonresonant HH production and evaluated for the benchmark scenario JHEP04 BM7 for the $3\ell+1 au_{h}$ category. The SM HH signal is shown for a cross section amounting to $30$ times the value predicted in the SM. The distributions expected for the background processes are shown for the values of nuisance parameters obtained from the ML fit of the signal+background hypothesis to the data. The binning is chosen to approximate background quantiles.

Distribution in the output of the BDT trained for nonresonant HH production and evaluated for the benchmark scenario JHEP04 BM7 for the $2\ell+2 au_{h}$ category. The SM HH signal is shown for a cross section amounting to $30$ times the value predicted in the SM. The distributions expected for the background processes are shown for the values of nuisance parameters obtained from the ML fit of the signal+background hypothesis to the data. The binning is chosen to approximate background quantiles.

Distribution in the output of the BDT trained for nonresonant HH production and evaluated for the benchmark scenario JHEP04 BM7 for the $1\ell+3 au_{h}$ category. The SM HH signal is shown for a cross section amounting to $30$ times the value predicted in the SM. The distributions expected for the background processes are shown for the values of nuisance parameters obtained from the ML fit of the signal+background hypothesis to the data. The binning is chosen to approximate background quantiles.

Distribution in the output of the BDT trained for nonresonant HH production and evaluated for the benchmark scenario JHEP04 BM7 for the $4 au_{h}$ category. The SM HH signal is shown for a cross section amounting to $30$ times the value predicted in the SM. The distributions expected for the background processes are shown for the values of nuisance parameters obtained from the ML fit of the signal+background hypothesis to the data. The binning is chosen to approximate background quantiles.

Observed and expected $95\%$ CL upper limits on the SM HH production cross section, obtained for both individual search categories and from a simultaneous fit of all seven categories combined.

Observed and expected $95\%$ CL upper limits on the HH production cross section as a function of the Higgs boson self-coupling strength modifier $\kappa_\lambda$ for the combination of all seven categories. All Higgs boson couplings other than $\lambda$ are assumed to have the values predicted in the SM.

Observed and expected $95\%$ CL upper limits on the HH production cross section as a function of the Higgs boson self-coupling strength modifier $\kappa_\lambda$ for the seven different categories as well as their combination. All Higgs boson couplings other than $\lambda$ are assumed to have the values predicted in the SM.

Observed and expected $95\%$ CL upper limits on the HH production cross section for the twelve benchmark scenarios from doi:10.1007/JHEP04(2016)126, the additional benchmark scenario 8a from doi:10.1007/JHEP09(2018)057, the seven benchmark scenarios from doi:10.1007/JHEP03(2020)091, and for the SM.

Observed and expected $95\%$ CL upper limits on the HH production cross section for the twelve benchmark scenarios from doi:10.1007/JHEP04(2016)126, the additional benchmark scenario 8a from doi:10.1007/JHEP09(2018)057, the seven benchmark scenarios from doi:10.1007/JHEP03(2020)091, and for the SM. The limits are shown for all seven search categories as well as for their combination.

Observed and expected $95\%$ CL upper limits on the HH production cross section as a function of the effective coupling $c_{2}$ for the combination of all seven categories. All Higgs boson couplings other than $c_{2}$ are assumed to have the values predicted in the SM.

Observed and expected $95\%$ CL upper limits on the HH production cross section as a function of the effective coupling $c_{2}$ and the top Yukawa coupling modifier $\kappa_{t}$ for the combination of all seven categories. All Higgs boson couplings other than $c_{2}$ and $\kappa_{t}$ are assumed to have the values predicted in the SM. The position of the SM in the ($c_{2}-\kappa_{t}$) plane, as well as the best fit value of $(c_{2},\kappa_{t})=(1.05, 1.74)$ together with contours for the theory cross section are shown as well.

Observed and expected $95\%$ CL upper limits on the HH production cross section as a function of the Higgs boson self-coupling modifier $\kappa_{\lambda}$ and the top Yukawa coupling modifier $\kappa_{t}$ for the combination of all seven categories. All Higgs boson couplings other than $\kappa_{\lambda}$ and $\kappa_{t}$ are assumed to have the values predicted in the SM. The position of the SM in the ($\kappa_{\lambda}-\kappa_{t}$) plane, as well as the best fit value of $(\kappa_{\lambda},\kappa_{t})=(-3.54, -1.70)$ together with contours for the theory cross section are shown as well.

Observed and expected $95\%$ CL upper limits on the HH production cross section as a function of the effective coupling $c_{2}$ and the Higgs boson self-coupling modifier $\kappa_{\lambda}$ for the combination of all seven categories. All Higgs boson couplings other than $c_{2}$ and $\kappa_{\lambda}$ are assumed to have the values predicted in the SM. The position of the SM in the ($c_{2}-\kappa_{\lambda}$) plane, as well as the best fit value of $(c_{2},\kappa_{\lambda})=(-0.66, 5.33)$ together with contours for the theory cross section are shown as well.

Observed and expected $95\%$ CL upper limits on the production of new particles X of spin $0$ and mass $m_{X}$ in the range $250$-$1000$ GeV, which decay to Higgs boson pairs.

Observed and expected $95\%$ CL upper limits on the production of new particles X of spin $0$ and mass $m_{X}$ in the range $250$-$1000$ GeV, which decay to Higgs boson pairs. The Limit is shown for the seven different search categories as well as their combination.

Observed and expected $95\%$ CL upper limits on the production of new particles X of spin $2$ and mass $m_{X}$ in the range $250$-$1000$ GeV, which decay to Higgs boson pairs.

Observed and expected $95\%$ CL upper limits on the production of new particles X of spin $2$ and mass $m_{X}$ in the range $250$-$1000$ GeV, which decay to Higgs boson pairs. The Limit is shown for the seven different search categories as well as their combination.

Distribution in the output of the BDT trained for resonances of spin 2 and mass $750\,$GeV production for the $2\ell$(ss) category. The resonant HH signal is shown for a cross section amounting to $1\,$pb. The distributions expected for the background processes are shown for the values of nuisance parameters obtained from the ML fit of the signal+background hypothesis to the data. The binning is chosen to approximate signal quantiles.

Distribution in the output of the BDT trained for resonances of spin 2 and mass $750\,$GeV production for the $3\ell$ category. The resonant HH signal is shown for a cross section amounting to $1\,$pb. The distributions expected for the background processes are shown for the values of nuisance parameters obtained from the ML fit of the signal+background hypothesis to the data. The binning is chosen to approximate signal quantiles.

Bin-to-bin correlation in the measured cross section as a function of the rapidity absolute value of the first jet, $|y(\text{j}_1)|$.

Measured cross section as a function of the transverse momenum of the first jet, $p_{\text{T}}(\text{j}_1)$, and breakdown of the relative uncertainty.

Bin-to-bin correlation in the measured cross section as a function of the transverse momenum of the first jet, $p_{\text{T}}(\text{j}_1)$.

Measured cross section as a function of the rapidity absolute value of the second jet, $|y(\text{j}_2)|$, and breakdown of the relative uncertainty.

Bin-to-bin correlation in the measured cross section as a function of the rapidity absolute value of the second jet, $|y(\text{j}_2)|$.

Measured cross section as a function of the transverse momenum of the second jet, $p_{\text{T}}(\text{j}_2)$, and breakdown of the relative uncertainty.

Bin-to-bin correlation in the measured cross section as a function of the transverse momenum of the second jet, $p_{\text{T}}(\text{j}_2)$.

Measured cross section as a function of the rapidity absolute value of the third jet, $|y(\text{j}_3)|$, and breakdown of the relative uncertainty.

Bin-to-bin correlation in the measured cross section as a function of the rapidity absolute value of the third jet, $|y(\text{j}_3)|$.

Measured cross section as a function of the transverse momenum of the third jet, $p_{\text{T}}(\text{j}_3)$, and breakdown of the relative uncertainty.

Bin-to-bin correlation in the measured cross section as a function of the transverse momenum of the third jet, $p_{\text{T}}(\text{j}_3)$.

Measured cross section as a function of the $H_{\text{T}}$ observable for events with at least one jet and breakdown of the relative uncertainty.

Bin-to-bin correlation in the measured cross section as a function of the $H_{\text{T}}$ observable for events with at least one jet.

Measured cross section as a function of the $H_{\text{T}}$ observable for events with at least two jets and breakdown of the relative uncertainty.

Bin-to-bin correlation in the measured cross section as a function of the $H_{\text{T}}$ observable for events with at least two jets.

Measured cross section as a function of the $H_{\text{T}}$ observable for events with at least three jets and breakdown of the relative uncertainty.

Bin-to-bin correlation in the measured cross section as a function of the $H_{\text{T}}$ observable for events with at least three jets.

Measured differential cross section as a function of dijet mass, $M_{\text{j}_1\text{j}_2}$, for events with at least two jets and breakdown of the relative uncertainty.

Bin-to-bin correlation in the measured differential cross section as a function of dijet mass, $M_{\text{j}_1\text{j}_2}$, for events with at least two jets.

Measured differential cross section as a function of Z boson rapidity absolute value, $|y(\text{Z})|$ for events with at least one jet and breakdown of the relative uncertainty.

Bin-to-bin correlation in the measured differential cross section as a function of Z boson rapidity absolute value, $|y(\text{Z})|$ for events with at least one jet.

Measured differential cross section as a function of the leading and subleading jet rapidity difference, $y_{\text{diff}}(\text{j}_1,\text{j}_2)$, for events with at least two jets and breakdown of the relative uncertainty.

Bin-to-bin correlation in the measured differential cross section as a function of the leading and subleading jet rapidity difference, $y_{\text{diff}}(\text{j}_1,\text{j}_2)$, for events with at least two jets.

Measured differential cross section as a function of the leading and subleading jet rapidity sum, $y_{\text{sum}}(\text{j}_1,\text{j}_2)$, for events with at least two jets and breakdown of the relative uncertainty.

Bin-to-bin correlation in the measured differential cross section as a function of the leading and subleading jet rapidity sum, $y_{\text{sum}}(\text{j}_1,\text{j}_2)$, for events with at least two jets.

Measured differential cross section as a function of the Z boson and leading jet rapidity difference, $y_{\text{diff}}(\text{Z},\text{j}_1)$, for events with at least one jet and breakdown of the relative uncertainty.

Bin-to-bin correlation in the measured differential cross section as a function of the Z boson and leading jet rapidity difference, $y_{\text{diff}}(\text{Z},\text{j}_1)$, for events with at least one jet.

Measured differential cross section as a function of the Z boson and leading jet rapidity sum, $y_{\text{sum}}(\text{Z},\text{j}_1)$, for events with at least one jet and breakdown of the relative uncertainty.

Bin-to-bin correlation in the measured differential cross section as a function of the Z boson and leading jet rapidity sum, $y_{\text{sum}}(\text{Z},\text{j}_1)$, for events with at least one jet.

Measured differential cross section as a function of the Z boson and leading jet rapidity difference, $y_{\text{diff}}(\text{Z},\text{j}_1)$, for events with at least two jets and breakdown of the relative uncertainty.

Bin-to-bin correlation in the measured differential cross section as a function of the Z boson and leading jet rapidity difference, $y_{\text{diff}}(\text{Z},\text{j}_1)$, for events with at least two jets.

Measured differential cross section as a function of the Z boson and leading jet rapidity sum, $y_{\text{sum}}(\text{Z},\text{j}_1)$, for events with at least two jets and breakdown of the relative uncertainty.

Bin-to-bin correlation in the measured differential cross section as a function of the Z boson and leading jet rapidity sum, $y_{\text{sum}}(\text{Z},\text{j}_1)$, for events with at least two jets.

Measured differential cross section as a function of the Z boson and subleading jet rapidity difference, $y_{\text{diff}}(\text{Z},\text{j}_2)$, for events with at least two jets and breakdown of the relative uncertainty.

Bin-to-bin correlation in the measured differential cross section as a function of the Z boson and subleading jet rapidity difference, $y_{\text{diff}}(\text{Z},\text{j}_2)$, for events with at least two jets.

Measured differential cross section as a function of the Z boson and subleading jet rapidity sum, $y_{\text{sum}}(\text{Z},\text{j}_2)$, for events with at least two jets and breakdown of the relative uncertainty.

Bin-to-bin correlation in the measured differential cross section as a function of the Z boson and subleading jet rapidity sum, $y_{\text{sum}}(\text{Z},\text{j}_2)$, for events with at least two jets.

Measured differential cross section as a function of the Z boson and dijet rapidity difference, $y_{\text{diff}}(\text{Z},\text{dijet})$, for events with at least two jets and breakdown of the relative uncertainty.

Bin-to-bin correlation in the measured differential cross section as a function of the Z boson and dijet rapidity difference, $y_{\text{diff}}(\text{Z},\text{dijet})$, for events with at least two jets.

Measured differential cross section as a function of the Z boson and dijet rapidity sum, $y_{\text{sum}}(\text{Z},\text{dijet})$, for events with at least two jets and breakdown of the relative uncertainty.

Bin-to-bin correlation in the measured differential cross section as a function of the Z boson and dijet rapidity sum, $y_{\text{sum}}(\text{Z},\text{dijet})$, for events with at least two jets.

Measured differential cross section as a function of the Z boson and leading jet azimuthal difference, $\Delta\phi(\text{Z},\text{j}_1)$, for events with at least one jet and breakdown of the relative uncertainty.

Bin-to-bin correlation in the measured differential cross section as a function of the Z boson and leading jet azimuthal difference, $\Delta\phi(\text{Z},\text{j}_1)$, for events with at least one jet.

Measured differential cross section as a function of the Z boson and leading jet azimuthal difference, $\Delta\phi(\text{Z},\text{j}_1)$, for events with at least two jets and breakdown of the relative uncertainty.

Bin-to-bin correlation in the measured differential cross section as a function of the Z boson and leading jet azimuthal difference, $\Delta\phi(\text{Z},\text{j}_1)$, for events with at least two jets.

Measured differential cross section as a function of the Z boson and leading jet azimuthal difference for events, $\Delta\phi(\text{Z},\text{j}_1)$, with at least three jets and breakdown of the relative uncertainty.

Bin-to-bin correlation in the measured differential cross section as a function of the Z boson and leading jet azimuthal difference for events, $\Delta\phi(\text{Z},\text{j}_1)$, with at least three jets.

Measured differential cross section as a function of the Z boson and subleading jet azimuthal difference for events with at least two jets and breakdown of the relative uncertainty.

Bin-to-bin correlation in the measured differential cross section as a function of the Z boson and subleading jet azimuthal difference for events with at least two jets.

Measured differential cross section as a function of the Z boson and subleading jet azimuthal difference, $\Delta\phi(\text{Z},\text{j}_2)$, for events with at least three jets and breakdown of the relative uncertainty.

Bin-to-bin correlation in the measured differential cross section as a function of the Z boson and subleading jet azimuthal difference, $\Delta\phi(\text{Z},\text{j}_2)$, for events with at least three jets.

Measured differential cross section as a function of the Z boson and third jet azimuthal difference, $\Delta\phi(\text{Z},\text{j}_3)$, for events with at least three jets and breakdown of the relative uncertainty.

Bin-to-bin correlation in the measured differential cross section as a function of the Z boson and third jet azimuthal difference, $\Delta\phi(\text{Z},\text{j}_3)$, for events with at least three jets.

Measured differential cross section as a function of the leading and subleading jet azimuthal difference, $\Delta\phi(\text{j}_1,\text{j}_2)$, for events with at least two jets and breakdown of the relative uncertainty.

Bin-to-bin correlation in the measured differential cross section as a function of the leading and subleading jet azimuthal difference, $\Delta\phi(\text{j}_1,\text{j}_2)$, for events with at least two jets.

Measured differential cross section as a function of the leading and subleading jet azimuthal difference, $\Delta\phi(\text{j}_1,\text{j}_2)$, for events with at least three jets and breakdown of the relative uncertainty.

Bin-to-bin correlation in the measured differential cross section as a function of the leading and subleading jet azimuthal difference, $\Delta\phi(\text{j}_1,\text{j}_2)$, for events with at least three jets.

Measured differential cross section as a function of the leading and third jet azimuthal difference, $\Delta\phi(\text{j}_1,\text{j}_3)$, for events with at least three jets and breakdown of the relative uncertainty.

Bin-to-bin correlation in the measured differential cross section as a function of the leading and third jet azimuthal difference, $\Delta\phi(\text{j}_1,\text{j}_3)$, for events with at least three jets.

Measured differential cross section as a function of the subleading and third jet azimuthal difference, $\Delta\phi(\text{j}_2,\text{j}_3)$, for events with at least three jets and breakdown of the relative uncertainty.

Bin-to-bin correlation in the measured differential cross section as a function of the subleading and third jet azimuthal difference, $\Delta\phi(\text{j}_2,\text{j}_3)$, for events with at least three jets.

Measured cross section as a function of the leading jet transverse momentum, $p_{\text{T}}(\text{j}_1)$, and rapidity absolute value, $|y(\text{j}_1)|$ for events with at least one jet and breakdown of the relative uncertainty.

Bin-to-bin correlation in the measured cross section as a function of the leading jet transverse momentum, $p_{\text{T}}(\text{j}_1)$, and rapidity absolute value, $|y(\text{j}_1)|$ for events with at least one jet.

Measured cross section as a function of the leading jet and Z boson rapidity abosolute values, $|y(\text{j}_1|$ and $|y(\text(Z))|$, for events with at least one jet with $p_\text{T}>20\,\text{Gev}$ and breakdown of the relative uncertainty.

Bin-to-bin correlation in the measured cross section as a function of the leading jet and Z boson rapidity abosolute values, $|y(\text{j}_1|$ and $|y(\text(Z))|$, for events with at least one jet with $p_\text{T}>20\,\text{Gev}$.

Measured cross section as a function of the transverse momementum, $p_{\text{T}}(\text{Z})$, and rapidity, $y(\text{Z})$, of the Z boson, for events with at least one jet with $p_\text{T}>20\,\text{Gev}$ and breakdown of the relative uncertainty.

Bin-to-bin correlation in the measured cross section as a function of the transverse momementum, $p_{\text{T}}(\text{Z})$, and rapidity, $y(\text{Z})$, of the Z boson, for events with at least one jet with $p_\text{T}>20\,\text{Gev}$.

Expected $+ 1$ s.d. lower limit contour on the $m_{\mathrm{W}_{\mathrm{KK}}}$ and $m_{\mathrm{R}}$ plane after combinign with an analysis of the all-hadronic final state.

Expected - 1 s.d. lower limit contour on the $m_{\mathrm{W}_{\mathrm{KK}}}$ and $m_{\mathrm{R}}$ plane after combinign with an analysis of the all-hadronic final state.

Observed lower limit contour on the $m_{\mathrm{W}_{\mathrm{KK}}}$ and $m_{\mathrm{R}}$ plane after combinign with an analysis of the all-hadronic final state.

Post-fit distributions of the reconstructed $\ell\nu$+jets system ($m_{\mathrm{j}\ell\nu}$, $m_{\mathrm{jj}\ell\nu}$) in data and simulation for SR1.

Post-fit distributions of the reconstructed $\ell\nu$+jets system ($m_{\mathrm{j}\ell\nu}$, $m_{\mathrm{jj}\ell\nu}$) in data and simulation for SR2.

Post-fit distributions of the reconstructed $\ell\nu$+jets system ($m_{\mathrm{j}\ell\nu}$, $m_{\mathrm{jj}\ell\nu}$) in data and simulation for SR3.

Post-fit distributions of the reconstructed $\ell\nu$+jets system ($m_{\mathrm{j}\ell\nu}$, $m_{\mathrm{jj}\ell\nu}$) in data and simulation for SR5.

Post-fit distributions of the reconstructed $\ell\nu$+jets system ($m_{\mathrm{j}\ell\nu}$, $m_{\mathrm{jj}\ell\nu}$) in data and simulation for SR6.

Distribution of $m_{\mathrm{jjj}}$ for preselected events with $\mathrm{N}_{j}$ = 3

Distribution of $m_{\mathrm{j}}$ for preselected events with $\mathrm{N}_{j}$ = 3

Distribution of the deep-WH value of the highest-mass jet with 60 < $m_{\mathrm{j}}$ < 100 GeV for preselected events with $\mathrm{N}_{j}$ = 3

scale factors (SFs) for W, $t^{2}$, and q/g matched jets in the low-$m_{\mathrm{j}}$ and low-$p_{\mathrm{T}}$ (LL) bin, as functions of the deep-W discriminant value.

scale factors (SFs) for W, $t^{2}$, and q/g matched jets in the low-$m_{\mathrm{j}}$ and high-$p_{\mathrm{T}}$ (LH) bin, as functions of the deep-W discriminant value.

scale factors (SFs) for $t^{2}$, $t^{3,4}$, and q/g matched jets in the high-$m_{\mathrm{j}}$ and low-$p_{\mathrm{T}}$ (HL) bin, as functions of the deep-WH discriminant value.

scale factors (SFs) for $t^{2}$, $t^{3,4}$, and q/g matched jets in the high-$m_{\mathrm{j}}$ and high-$p_{\mathrm{T}}$ (HH) bin, as functions of the deep-WH discriminant value.

The deep-W discriminant of the jet with highest mass in the single-lepton sideband for LL samples.

The deep-W discriminant of the jet with highest mass in the single-lepton sideband for LH samples.

The deep-WH discriminant of the jet with highest mass in the single-lepton sideband for HL samples.

The deep-WH discriminant of the jet with highest mass in the single-lepton sideband for HH samples.

Comparison of the distribution of data and simulated backgrounds, as a function of the deep-W discriminant value for the highest-mass jet in CR1, after the scale factors have been applied.

Comparison of the distribution of data and simulated backgrounds, as a function of the deep-WH discriminant value for the highest-mass jet in CR2, after the scale factors have been applied.

Comparison of the distribution of data and simulated backgrounds, as a function of the deep-WH discriminant value for the highest-mass jet in CR3, after the scale factors have been applied.

Comparison of the distribution of data and simulated backgrounds, as a function of the deep-W discriminant value for the highest-mass jet in CR45, after the scale factors have been applied.

Comparison of the distribution of data and simulated backgrounds, as a function of the deep-W discriminant value for the highest-mass jet in CR6, after the scale factors have been applied.

The $m_{\mathrm{j}^{\mathrm{max}}}$ distributions for different jet types for SR1–3 events of the signal with $m_{\mathrm{W}_{\mathrm{KK}}}$ = 2.5 TeV, $m_{\mathrm{R}}$ = 0.2 TeV$ without deep-W (WH) constraints.

The deep-W distribution normalized to unity for the shown components of of the signal with $m_{\mathrm{W}_{\mathrm{KK}}}$ = 2.5 TeV, $m_{\mathrm{R}}$ = 0.2 TeV. The $t^{3,4}$ jets from the preselected sample, normalized to unity, are superimposed to compare shapes with the $R^{3q}$ and $R^{4q}$ distributions.

The deep-WH distribution normalized to unity for the shown components of of the signal with $m_{\mathrm{W}_{\mathrm{KK}}}$ = 2.5 TeV, $m_{\mathrm{R}}$ = 0.2 TeV. The $t^{3,4}$ jets from the preselected sample, normalized to unity, are superimposed to compare shapes with the $R^{3q}$ and $R^{4q}$ distributions.

The $m_{\mathrm{jj}}$ distribution for CR1 for data and simulation.

The $m_{\mathrm{jj}}$ distribution for CR2 for data and simulation.

The $m_{\mathrm{jj}}$ distribution for CR3 for data and simulation.

The $m_{\mathrm{jjj}}$ distribution for CR45 for data and simulation.

The $m_{\mathrm{jjj}}$ distribution for CR6 for data and simulation.

Post-fit distributions of the reconstructed triboson system ($m_{\mathrm{jj}}$) in data and simulation for SR1.

Post-fit distributions of the reconstructed triboson system ($m_{\mathrm{jj}}$) in data and simulation for SR2.

Post-fit distributions of the reconstructed triboson system ($m_{\mathrm{jj}}$) in data and simulation for SR3.

Post-fit distributions of the reconstructed triboson system ($m_{\mathrm{jjj}}$) in data and simulation for SR4.

Post-fit distributions of the reconstructed triboson system ($m_{\mathrm{jjj}}$) in data and simulation for SR5.

Post-fit distributions of the reconstructed triboson system ($m_{\mathrm{jjj}}$) in data and simulation for SR6.

Observed upper limits at $95\%$ CL on the signal cross section $\times$ branching fraction as functions of the $m_{\mathrm{W}_{\mathrm{KK}}}$ and $m_{\mathrm{R}}$ resonance masses.

Expected median lower limit contour on the $m_{\mathrm{W}_{\mathrm{KK}}}$ and $m_{\mathrm{R}}$ plane.

Expected $+ 1$ s.d. lower limit contour on the $m_{\mathrm{W}_{\mathrm{KK}}}$ and $m_{\mathrm{R}}$ plane.

Expected - 1 s.d. lower limit contour on the $m_{\mathrm{W}_{\mathrm{KK}}}$ and $m_{\mathrm{R}}$ plane.

Observed lower limit contour on the $m_{\mathrm{W}_{\mathrm{KK}}}$ and $m_{\mathrm{R}}$ plane.

Final state radiation correction used in the $\phi_{\eta}^{*}$ cross-section measurement. The first uncertainty is statistical and the second is systematic.

Final state radiation correction used in the $y^{Z}-p_{T}^{Z}$ cross-section measurement. The first uncertainty is statistical and the second is systematic.

Final state radiation correction used in the $y^{Z}-\phi_{\eta}^{*}$ cross-section measurement. The first uncertainty is statistical and the second is systematic.

Correlation matrix of statistical uncertainty for one-dimensional $y^Z$ measurement.

Correlation matrix of statistical uncertainty for one-dimensional $p_{T}^{Z}$ measurement.

Correlation matrix of statistical uncertainty for one-dimensional $\phi_{\eta}^{*}$ measurement.

Correlation matrix of statistical uncertainty for two-dimensional $y^Z-p_{T}^{Z}$ measurement.

Correlation matrix of statistical uncertainty for two-dimensional $y^Z-\phi_{\eta}^{*}$ measurement.

Correlation matrix of efficiency uncertainty for one-dimensional $y^Z$ measurement.

Correlation matrix of efficiency uncertainty for one-dimensional $p_{T}^{Z}$ measurement.

Correlation matrix of efficiency uncertainty for one-dimensional $\phi_{\eta}^{*}$ measurement.

Correlation matrix of efficiency uncertainty for two-dimensional $y^Z-p_{T}^{Z}$ measurement.

Correlation matrix of efficiency uncertainty for two-dimensional $y^Z-\phi_{\eta}^{*}$ measurement.

Measured total $Z$-boson cross-section for different datasets. The first uncertainty is statistical, the second systematic, and the third is due to the luminosity.

Measured single differential cross-sections in interval regions of $y^{Z}$. The first uncertainty is statistical, the second systematic, and the third is due to the luminosity.

Measured single differential cross-sections in interval regions of $p_{T}^{Z}$. The first uncertainty is statistical, the second systematic, and the third is due to the luminosity.

Measured single differential cross-sections in interval regions of $\phi_{\eta}^{*}$. The first uncertainty is statistical, the second systematic, and the third is due to the luminosity.

Measured double differential cross-sections in interval regions of $y^{Z}$ and $p_{T}^{Z}$. The first uncertainty is statistical, the second systematic, and the third is due to the luminosity.

Measured double differential cross-sections in interval regions of $y^{Z}$ and $\phi_{\eta}^{*}$. The first uncertainty is statistical, the second systematic, and the third is due to the luminosity.

Systematic uncertainties in the single differential cross-sections in interval regions of $y^{Z}$, presented in percentage. The contributions from efficiency (Eff), background (BKG), final state radiation (FSR), closure test (Closure), and alignment and calibration (Alignment) are shown.

Systematic uncertainties in the single differential cross-sections in interval regions of $p_{T}^{Z}$, presented in percentage. The contributions from efficiency (Eff), background (BKG), final state radiation (FSR), closure test (Closure), and alignment and calibration (Alignment) are shown.

Systematic uncertainties in the single differential cross-sections in interval regions of $\phi_{\eta}^{*}$, presented in percentage. The contributions from efficiency (Eff), background (BKG), final state radiation (FSR), closure test (Closure), and alignment and calibration (Alignment) are shown.

Systematic uncertainties in the double differential cross-sections in interval regions of $y^{Z}$ and $p_{T}^{Z}$, presented in percentage. The contributions from efficiency (Eff), background (BKG), final state radiation (FSR), closure test (Closure), and alignment and calibration (Alignment) are shown.

Systematic uncertainties in the double differential cross-sections in interval regions of $y^{Z}$ and $\phi_{\eta}^{*}$, presented in percentage. The contributions from efficiency (Eff), background (BKG), final state radiation (FSR), closure test (Closure), and alignment and calibration (Alignment) are shown.

Measured DPS effective cross section

pp -> J/psi J/psi J/psi X cross section

$pp \rightarrow J/\psi J/\psi J/\psi X~$ fiducial cross section

$M_{p\eta}$ invariant mass distribution for $E_\gamma$ = 1400-1500 MeV.

$M_{p\eta}$ invariant mass distribution for $E_\gamma$ = 1500-1600 MeV.

$M_{p\eta}$ invariant mass distribution under the conditions: $E_\gamma$ = 1400-1500 MeV; 1120 MeV $\leq M_{p\pi^0} \leq$ 1220 MeV; $\theta_p^{lab} \leq$ 25°; 25° $\leq \theta_\gamma^{lab} \leq$ 155°.

$M_{p\eta}$ invariant mass distribution under the conditions: $E_\gamma$ = 1400-1500 MeV; 1120 MeV $\leq M_{p\pi^0} \leq$ 1220 MeV; $\theta_p^{lab} \leq$ 25°; 1° $\leq \theta_\gamma^{lab} \leq$ 155°.

$M_{p\eta}$ invariant mass distribution for $E_\gamma$ = 1400-1500 MeV and $M_{p\pi^0} \leq$ 1190 MeV, allowing for all proton and photon laboratory angles.

$M_{p\eta}$ invariant mass distribution for the incident photon energy range $E_\gamma$ = 1420-1540 MeV and $M_{p\pi^0} \leq$ 1170 MeV.

$M_{p\eta}$ invariant mass distribution for the incident photon energy range $E_\gamma$ = 1420-1540 MeV and $M_{p\pi^0} \leq$ 1180 MeV.

$M_{p\eta}$ invariant mass distribution for the incident photon energy range $E_\gamma$ = 1420-1540 MeV and $M_{p\pi^0} \leq$ 1190 MeV.

$M_{p\eta}$ invariant mass distribution for the incident photon energy range $E_\gamma$ = 1420-1540 MeV and $M_{p\pi^0} \leq$ 1200 MeV.

$M_{p\eta}$ invariant mass distribution for the incident photon energy range $E_\gamma$ = 1420-1540 MeV and $M_{p\pi^0} \leq$ 1210 MeV.

$M_{p\eta}$ invariant mass distribution for the incident photon energy range $E_\gamma$ = 1420-1540 MeV and $M_{p\pi^0} \leq$ 1220 MeV.

$M_{p\eta}$ invariant mass distribution for the incident photon energy range $E_\gamma$ = 1400-1420 MeV and $M_{p\pi^0} \leq$ 1190 MeV.

$M_{p\eta}$ invariant mass distribution for the incident photon energy range $E_\gamma$ = 1420-1460 MeV and $M_{p\pi^0} \leq$ 1190 MeV.

$M_{p\eta}$ invariant mass distribution for the incident photon energy range $E_\gamma$ = 1460-1500 MeV and $M_{p\pi^0} \leq$ 1190 MeV.

$M_{p\eta}$ invariant mass distribution for the incident photon energy range $E_\gamma$ = 1500-1540 MeV and $M_{p\pi^0} \leq$ 1190 MeV.

$M_{p\eta}$ invariant mass distribution for the incident photon energy range $E_\gamma$ = 1540-1600 MeV and $M_{p\pi^0} \leq$ 1190 MeV.

$M_{p\eta}$ invariant mass distribution for the incident photon energy range $E_\gamma$ = 1420-1540 MeV and $M_{p\pi^0} \leq$ 1190 MeV.

Excess yield relative to the PWA distributions in the invariant mass distributions $M^2_{p\eta}$.

Excess yield relative to the PWA distributions in the invariant mass distributions $M^2_{\pi^0\eta}$.

Excess yield relative to the PWA distributions in the invariant mass distributions $M^2_{p\pi^0}$, with a cut on the structure in the $M^2_{p\eta}$ invariant mass distribution.

Excess yield relative to the PWA distributions in the $p-\eta$ opening angle distribution in the $\gamma p$ centre-of-mass system.

Excess yield relative to the PWA distributions in the $p-\eta$ opening angle distribution in the $\gamma p$ centre-of-mass system.

Excess yield relative to the PWA distributions in the $\pi^0-\eta$ opening angle distribution in the $\gamma p$ centre-of-mass system.

Excess yield relative to the PWA distributions in the $\pi^0-\eta$ opening angle distribution in the $\gamma p$ centre-of-mass system.

Differential cross section $d\sigma/dM_{\pi^0\eta}$ for $E_\gamma$ = 1540-1600 MeV.

$M_{p\eta}$ invariant mass distribution, taken as difference to the PWA calculations (s. Fig. 11).

$M_{\pi^0\eta}$ invariant mass distribution, taken as difference to the PWA calculations (s. Fig. 11).

$\theta_{p\eta}$ opening angle distribution, taken as difference to the PWA calculations (s. Fig. 12).

$\theta_{\pi^0\eta}$ opening angle distribution, taken as difference to the PWA calculations (s. Fig. 12).

The contribution from different systematic uncertainties to the measured $t\bar{t}t\bar{t}$ production cross section, grouped in categories.

The results of the fitted signal strength $\mu$ in the 1L/2LOS and 2LSS/3L combined channel.

The results of fitted inclusive ${t\bar{t}t\bar{t}}$ cross-section in the 1L/2LOS and 2LSS/3L combined channel.

Comparison between data and prediction for the distribution of the sum of the pseudo-continuous b-tagging score over the six jets with the highest score in the 1L,$\geq$9j,$\geq$3b region before the fit.

Comparison between data and prediction for the distribution of the sum of the pseudo-continuous b-tagging score over the six jets with the highest score in the 1L,$\geq$9j,$\geq$3b region after the fit.

Comparison between data and prediction for the distribution of the sum of the pseudo-continuous b-tagging score over the six jets with the highest score in the 2LOS,$\geq$7j,$\geq$3b region before the fit.

Comparison between data and prediction for the distribution of the sum of the pseudo-continuous b-tagging score over the six jets with the highest score in the 2LOS,$\geq$7j,$\geq$3b region after the fit.

Comparison between data and prediction for the distribution of the scalar sum of all jet and lepton pT in the event in the 1L,$\geq$8j,$\geq$3b region before the fit.

Comparison between data and prediction for the distribution of the scalar sum of all jet and lepton pT in the event in the 1L,$\geq$8j,$\geq$3b region after the fit.

Comparison between data and prediction for the distribution of the scalar sum of all jet and lepton pT in the event in the 2LOS,$\geq$6j,$\geq$3b region before the fit.

Comparison between data and prediction for the distribution of the scalar sum of all jet and lepton pT in the event in the 2LOS,$\geq$6j,$\geq$3b region after the fit.

Comparison between data and prediction for the distribution of number of jets in the 1L,$\geq$8j,$\geq$3b region before the fit.

Comparison between data and prediction for the distribution of number of jets in the 1L,$\geq$8j,$\geq$3b region after the fit.

Comparison between data and prediction for the distribution of number of jets in the 2LOS,$\geq$6j,$\geq$3b region before the fit.

Comparison between data and prediction for the distribution of number of jets in the 2LOS,$\geq$6j,$\geq$3b region after the fit.

Comparison between data and prediction for the distribution of b-jets multiplicity in the 1L,$\geq$8j,$\geq$3b region before the fit.

Comparison between data and prediction for the distribution of b-jets multiplicity in the 1L,$\geq$8j,$\geq$3b region after the fit.

Comparison between data and prediction for the distribution of b-jets multiplicity in the 2LOS,$\geq$6j,$\geq$3b region before the fit.

Comparison between data and prediction for the distribution of b-jets multiplicity in the 2LOS,$\geq$6j,$\geq$3b region after the fit.

Comparison between data and prediction for the distribution of the BDT score in the 1L,9j,4b signal region before the fit.

Comparison between data and prediction for the distribution of the BDT score in the 1L,9j,4b signal region after the fit.

Comparison between data and prediction for the distribution of the BDT score in the 1L,9j,$\geq$5b signal region before the fit.

Comparison between data and prediction for the distribution of the BDT score in the 1L,9j,$\geq$5b signal region after the fit.

Comparison between data and prediction for the distribution of the BDT score in the 1L,$\geq$10j,3bL signal region before the fit.

Comparison between data and prediction for the distribution of the BDT score in the 1L,$\geq$10j,3bL signal region after the fit.

Comparison between data and prediction for the distribution of the BDT score in the 1L,$\geq$10j,3bH signal region before the fit.

Comparison between data and prediction for the distribution of the BDT score in the 1L,$\geq$10j,3bH signal region after the fit.

Comparison between data and prediction for the distribution of the BDT score in the 1L,$\geq$10j,4b signal region before the fit.

Comparison between data and prediction for the distribution of the BDT score in the 1L,$\geq$10j,4b signal region after the fit.

Comparison between data and prediction for the distribution of the BDT score in the 1L,$\geq$10j,$\geq$5b signal region before the fit.

Comparison between data and prediction for the distribution of the BDT score in the 1L,$\geq$10j,$\geq$5b signal region after the fit.

Comparison between data and prediction for the distribution of the BDT score in the 2LOS,7j,$\geq$4b signal region before the fit.

Comparison between data and prediction for the distribution of the BDT score in the 2LOS,7j,$\geq$4b signal region after the fit.

Comparison between data and prediction for the distribution of the BDT score in the 2LOS,$\geq$8j,3bL signal region before the fit.

Comparison between data and prediction for the distribution of the BDT score in the 2LOS,$\geq$8j,3bL signal region after the fit.

Comparison between data and prediction for the distribution of the BDT score in the 2LOS,$\geq$8j,3bH signal region before the fit.

Comparison between data and prediction for the distribution of the BDT score in the 2LOS,$\geq$8j,3bH signal region after the fit.

Comparison between data and prediction for the distribution of the BDT score in the 2LOS,$\geq$8j,$\geq$4b signal region before the fit.

Comparison between data and prediction for the distribution of the BDT score in the 2LOS,$\geq$8j,$\geq$4b signal region after the fit.

Comparison between data and prediction for the distribution of the scalar sum of all jet and lepton pT in the event in the 1L,8j,3bV validation region before the fit.

Comparison between data and prediction for the distribution of the scalar sum of all jet and lepton pT in the event in the 1L,8j,3bV validation region after the fit.

Comparison between data and prediction for the distribution of the BDT score in the 1L,9j,3bV validation region before the fit.

Comparison between data and prediction for the distribution of the BDT score in the 1L,9j,3bV validation region after the fit.

Comparison between data and prediction for the distribution of the BDT score in the 1L,$\geq$10j,3bV validation region before the fit.

Comparison between data and prediction for the distribution of the BDT score in the 1L,$\geq$10j,3bV validation region after the fit.

Comparison between data and prediction for the distribution of the scalar sum of all jet and lepton pT in the event in the 2LOS,6j,3bV validation region before the fit.

Comparison between data and prediction for the distribution of the scalar sum of all jet and lepton pT in the event in the 2LOS,6j,3bV validation region after the fit.

Comparison between data and prediction for the distribution of the BDT score in the 2LOS,7j,3bV validation region before the fit.

Comparison between data and prediction for the distribution of the BDT score in the 2LOS,7j,3bV validation region after the fit.

Comparison between data and prediction for the distribution of the BDT score in the 2LOS,$\geq$8j,3bV validation region before the fit.

Comparison between data and prediction for the distribution of the BDT score in the 2LOS,$\geq$8j,3bV validation region after the fit.

Comparison between data and prediction for the distribution of the scalar sum of all jet and lepton pT in the event in the 1L,8j,3bL control region before the fit.

Comparison between data and prediction for the distribution of the scalar sum of all jet and lepton pT in the event in the 1L,8j,3bL control region after the fit.

Comparison between data and prediction for the distribution of the scalar sum of all jet and lepton pT in the event in the 1L,8j,3bH control region before the fit.

Comparison between data and prediction for the distribution of the scalar sum of all jet and lepton pT in the event in the 1L,8j,3bH control region after the fit.

Comparison between data and prediction for the distribution of the scalar sum of all jet and lepton pT in the event in the 1L,9j,3bL control region before the fit.

Comparison between data and prediction for the distribution of the scalar sum of all jet and lepton pT in the event in the 1L,9j,3bL control region after the fit.

Comparison between data and prediction for the distribution of the scalar sum of all jet and lepton pT in the event in the 1L,9j,3bH control region before the fit.

Comparison between data and prediction for the distribution of the scalar sum of all jet and lepton pT in the event in the 1L,9j,3bH control region after the fit.

Comparison between data and prediction for the distribution of the scalar sum of all jet and lepton pT in the event in the 1L,8j,4b control region before the fit.

Comparison between data and prediction for the distribution of the scalar sum of all jet and lepton pT in the event in the 1L,8j,4b control region after the fit.

Comparison between data and prediction for the distribution of the scalar sum of all jet and lepton pT in the event in the 1L,8j,$\geq$5b control region before the fit.

Comparison between data and prediction for the distribution of the scalar sum of all jet and lepton pT in the event in the 1L,8j,$\geq$5b control region after the fit.

Comparison between data and prediction for the distribution of the scalar sum of all jet and lepton pT in the event in the 2LOS,6j,3bL control region before the fit.

Comparison between data and prediction for the distribution of the scalar sum of all jet and lepton pT in the event in the 2LOS,6j,3bL control region after the fit.

Comparison between data and prediction for the distribution of the scalar sum of all jet and lepton pT in the event in the 2LOS,6j,3bH control region before the fit.

Comparison between data and prediction for the distribution of the scalar sum of all jet and lepton pT in the event in the 2LOS,6j,3bH control region after the fit.

Comparison between data and prediction for the distribution of the scalar sum of all jet and lepton pT in the event in the 2LOS,7j,3bL control region before the fit.

Comparison between data and prediction for the distribution of the scalar sum of all jet and lepton pT in the event in the 2LOS,7j,3bL control region after the fit.

Comparison between data and prediction for the distribution of the scalar sum of all jet and lepton pT in the event in the 2LOS,7j,3bH control region before the fit.

Comparison between data and prediction for the distribution of the scalar sum of all jet and lepton pT in the event in the 2LOS,7j,3bH control region after the fit.

Comparison between data and prediction for the distribution of the scalar sum of all jet and lepton pT in the event in the 2LOS,6j,$\geq$4b control region before the fit.

Comparison between data and prediction for the distribution of the scalar sum of all jet and lepton pT in the event in the 2LOS,6j,$\geq$4b control region after the fit.

The measured single-differential cross sections in leading jet transverse momenta for the pure electroweak Z$\gamma$jj production. The total uncertainty includes the stastical uncertianty and the systematic uncertainty, while the uncertainty of the predicted results is the theoretical uncertainty from the MadGraph5_aMC@NLO. The last bin includes overflow events.

The measured single-differential cross sections in leading lepton transverse momenta for the pure electroweak Z$\gamma$jj production. The total uncertainty includes the stastical uncertianty and the systematic uncertainty, while the uncertainty of the predicted results is the theoretical uncertainty from the MadGraph5_aMC@NLO. The last bin includes overflow events.

The measured double-differential cross sections as functions of $m_{\mathrm{jj}}$ and $|\Delta\eta_{\mathrm{jj}}|$ for the pure electroweak Z$\gamma$jj production. The total uncertainty includes the stastical uncertianty and the systematic uncertainty, while the uncertainty of the predicted results is the theoretical uncertainty from the MadGraph5_aMC@NLO. The last bin includes overflow events.

The measured single-differential cross sections in photon transverse momenta for the combined QCD-induced and electroweak Z$\gamma$jj production. The total uncertainty includes the stastical uncertianty and the systematic uncertainty, while the uncertainty of the predicted results is the theoretical uncertainty from the MadGraph5_aMC@NLO. The last bin includes overflow events.

The measured single-differential cross sections in leading jet transverse momenta for the combined QCD-induced and electroweak Z$\gamma$jj production. The total uncertainty includes the stastical uncertianty and the systematic uncertainty, while the uncertainty of the predicted results is the theoretical uncertainty from the MadGraph5_aMC@NLO. The last bin includes overflow events.

The measured single-differential cross sections in leading lepton transverse momenta for the combined QCD-induced and electroweak Z$\gamma$jj production. The total uncertainty includes the stastical uncertianty and the systematic uncertainty, while the uncertainty of the predicted results is the theoretical uncertainty from the MadGraph5_aMC@NLO. The last bin includes overflow events.

The measured double-differential cross sections as functions of $m_{\mathrm{jj}}$ and $|\Delta\eta_{\mathrm{jj}}|$ for the combined QCD-induced and electroweak Z$\gamma$jj production. The total uncertainty includes the stastical uncertianty and the systematic uncertainty, while the uncertainty of the predicted results is the theoretical uncertainty from the MadGraph5_aMC@NLO. The last bin includes overflow events.

The expected and observed limits on the aQGC parameters at 95% confidence level. The last column presents the scattering energy values for which the amplitude would violate unitarity for the observed value of the aQGC parameter. All coupling parameter limits are set in TeV$^{-4}$, whereas the unitarity bounds are in TeV.

List of relative uncertainties of the measured inclusive $t\bar{t}\text{Z}$ cross section from the combined fit. The uncertainties are symmetrised for presentation and grouped into the categories described in the text. The quadratic sum of the individual uncertainties is not equal to the total uncertainty due to correlations introduced by the fit.

The definitions of the trilepton signal regions: for the inclusive measurement, a combination of the regions with pseudo-continuous $b$-tagging 3$\ell$-Z-1$b$4$j$-PCBT and 3$\ell$-Z-2$b$3$j$-PCBT is used, whereas for the differential measurement, only the region 3$\ell$-Z-2$b$3$j$, with a fixed $b$-tagging WP is employed.

The definitions of the trilepton signal regions: for the inclusive measurement, a combination of the regions with pseudo-continuous $b$-tagging 3$\ell$-Z-1$b$4$j$-PCBT and 3$\ell$-Z-2$b$3$j$-PCBT is used, whereas for the differential measurement, only the region 3$\ell$-Z-2$b$3$j$, with a fixed $b$-tagging WP is employed.

The definitions of the four tetralepton signal regions. The regions are defined to target different $b$-jet multiplicities and flavour combinations of the non-Z leptons.

The definitions of the four tetralepton signal regions. The regions are defined to target different $b$-jet multiplicities and flavour combinations of the non-Z leptons.

The absolute particle-level differential cross-section measured in the fiducial phase-space as a function of the $p_{\textrm{T}}$ of the $Z$ boson in the 3$\ell$+4$\ell$ channel. The uncertainty is decomposed into four components which are the signal modelling uncertainty, the background modelling uncertainty, the experimental uncertainty, and the data statistical uncertainty.

The absolute particle-level differential cross-section measured in the fiducial phase-space as a function of the $p_{\textrm{T}}$ of the $Z$ boson in the 3$\ell$+4$\ell$ channel. The uncertainty is decomposed into four components which are the signal modelling uncertainty, the background modelling uncertainty, the experimental uncertainty, and the data statistical uncertainty.

The absolute parton-level differential cross-section measured in the fiducial phase-space as a function of the $p_{\textrm{T}}$ of the $Z$ boson in the 3$\ell$+4$\ell$ channel. The uncertainty is decomposed into four components which are the signal modelling uncertainty, the background modelling uncertainty, the experimental uncertainty, and the data statistical uncertainty.

The absolute parton-level differential cross-section measured in the fiducial phase-space as a function of the $p_{\textrm{T}}$ of the $Z$ boson in the 3$\ell$+4$\ell$ channel. The uncertainty is decomposed into four components which are the signal modelling uncertainty, the background modelling uncertainty, the experimental uncertainty, and the data statistical uncertainty.

The normalised particle-level differential cross-section measured in the fiducial phase-space as a function of the $p_{\textrm{T}}$ of the $Z$ boson in the 3$\ell$+4$\ell$ channel. The uncertainty is decomposed into four components which are the signal modelling uncertainty, the background modelling uncertainty, the experimental uncertainty, and the data statistical uncertainty.

The normalised particle-level differential cross-section measured in the fiducial phase-space as a function of the $p_{\textrm{T}}$ of the $Z$ boson in the 3$\ell$+4$\ell$ channel. The uncertainty is decomposed into four components which are the signal modelling uncertainty, the background modelling uncertainty, the experimental uncertainty, and the data statistical uncertainty.

The normalised parton-level differential cross-section measured in the fiducial phase-space as a function of the $p_{\textrm{T}}$ of the $Z$ boson in the 3$\ell$+4$\ell$ channel. The uncertainty is decomposed into four components which are the signal modelling uncertainty, the background modelling uncertainty, the experimental uncertainty, and the data statistical uncertainty.

The normalised parton-level differential cross-section measured in the fiducial phase-space as a function of the $p_{\textrm{T}}$ of the $Z$ boson in the 3$\ell$+4$\ell$ channel. The uncertainty is decomposed into four components which are the signal modelling uncertainty, the background modelling uncertainty, the experimental uncertainty, and the data statistical uncertainty.

The absolute parton-level differential cross-section measured in the fiducial phase-space as a function of the absolute value of rapidity of the $Z$ boson in the 3$\ell$+4$\ell$ channel. The uncertainty is decomposed into four components which are the signal modelling uncertainty, the background modelling uncertainty, the experimental uncertainty, and the data statistical uncertainty.

The absolute parton-level differential cross-section measured in the fiducial phase-space as a function of the absolute value of rapidity of the $Z$ boson in the 3$\ell$+4$\ell$ channel. The uncertainty is decomposed into four components which are the signal modelling uncertainty, the background modelling uncertainty, the experimental uncertainty, and the data statistical uncertainty.

The absolute particle-level differential cross-section measured in the fiducial phase-space as a function of the number of jets in the 3$\ell$ channel. The uncertainty is decomposed into four components which are the signal modelling uncertainty, the background modelling uncertainty, the experimental uncertainty, and the data statistical uncertainty.

The absolute particle-level differential cross-section measured in the fiducial phase-space as a function of the number of jets in the 3$\ell$ channel. The uncertainty is decomposed into four components which are the signal modelling uncertainty, the background modelling uncertainty, the experimental uncertainty, and the data statistical uncertainty.

The absolute particle-level differential cross-section measured in the fiducial phase-space as a function of the number of jets in the 4$\ell$ channel. The uncertainty is decomposed into four components which are the signal modelling uncertainty, the background modelling uncertainty, the experimental uncertainty, and the data statistical uncertainty.

The absolute particle-level differential cross-section measured in the fiducial phase-space as a function of the number of jets in the 4$\ell$ channel. The uncertainty is decomposed into four components which are the signal modelling uncertainty, the background modelling uncertainty, the experimental uncertainty, and the data statistical uncertainty.

The absolute parton-level differential cross-section measured in the fiducial phase-space as a function of the $p_{\textrm{T}}^{l \textrm{non-}Z}$ in the 3$\ell$ channel. The uncertainty is decomposed into four components which are the signal modelling uncertainty, the background modelling uncertainty, the experimental uncertainty, and the data statistical uncertainty.

The absolute parton-level differential cross-section measured in the fiducial phase-space as a function of the $p_{\textrm{T}}^{l \textrm{non-}Z}$ in the 3$\ell$ channel. The uncertainty is decomposed into four components which are the signal modelling uncertainty, the background modelling uncertainty, the experimental uncertainty, and the data statistical uncertainty.

The absolute parton-level differential cross-section measured in the fiducial phase-space as a function of the $|\Delta \phi (Z, t_{\textrm{lep}})|/\pi$ in the 3$\ell$ channel. The uncertainty is decomposed into four components which are the signal modelling uncertainty, the background modelling uncertainty, the experimental uncertainty, and the data statistical uncertainty.

The absolute parton-level differential cross-section measured in the fiducial phase-space as a function of the $|\Delta \phi (Z, t_{\textrm{lep}})|/\pi$ in the 3$\ell$ channel. The uncertainty is decomposed into four components which are the signal modelling uncertainty, the background modelling uncertainty, the experimental uncertainty, and the data statistical uncertainty.

The absolute parton-level differential cross-section measured in the fiducial phase-space as a function of the $|\Delta y (Z, t_{\textrm{lep}})|$ in the 3$\ell$ channel. The uncertainty is decomposed into four components which are the signal modelling uncertainty, the background modelling uncertainty, the experimental uncertainty, and the data statistical uncertainty.

The absolute parton-level differential cross-section measured in the fiducial phase-space as a function of the $|\Delta y (Z, t_{\textrm{lep}})|$ in the 3$\ell$ channel. The uncertainty is decomposed into four components which are the signal modelling uncertainty, the background modelling uncertainty, the experimental uncertainty, and the data statistical uncertainty.

The absolute parton-level differential cross-section measured in the fiducial phase-space as a function of the $|\Delta \phi (l_{t}^{+}, l_{\bar{t}}^{-})|/\pi$ in the 4$\ell$ channel. The uncertainty is decomposed into four components which are the signal modelling uncertainty, the background modelling uncertainty, the experimental uncertainty, and the data statistical uncertainty.

The absolute parton-level differential cross-section measured in the fiducial phase-space as a function of the $|\Delta \phi (l_{t}^{+}, l_{\bar{t}}^{-})|/\pi$ in the 4$\ell$ channel. The uncertainty is decomposed into four components which are the signal modelling uncertainty, the background modelling uncertainty, the experimental uncertainty, and the data statistical uncertainty.

The absolute parton-level differential cross-section measured in the fiducial phase-space as a function of the $|\Delta \phi (t\bar{t}, Z)|/\pi$ in the 4$\ell$ channel. The uncertainty is decomposed into four components which are the signal modelling uncertainty, the background modelling uncertainty, the experimental uncertainty, and the data statistical uncertainty.

The absolute parton-level differential cross-section measured in the fiducial phase-space as a function of the $|\Delta \phi (t\bar{t}, Z)|/\pi$ in the 4$\ell$ channel. The uncertainty is decomposed into four components which are the signal modelling uncertainty, the background modelling uncertainty, the experimental uncertainty, and the data statistical uncertainty.

The absolute parton-level differential cross-section measured in the fiducial phase-space as a function of the $p_{\textrm{T}}^{t\bar{t}}$ in the 4$\ell$ channel. The uncertainty is decomposed into four components which are the signal modelling uncertainty, the background modelling uncertainty, the experimental uncertainty, and the data statistical uncertainty.

The absolute parton-level differential cross-section measured in the fiducial phase-space as a function of the $p_{\textrm{T}}^{t\bar{t}}$ in the 4$\ell$ channel. The uncertainty is decomposed into four components which are the signal modelling uncertainty, the background modelling uncertainty, the experimental uncertainty, and the data statistical uncertainty.

The normalised parton-level differential cross-section measured in the fiducial phase-space as a function of the $p_{\textrm{T}}^{l \textrm{non-}Z}$ in the 3$\ell$ channel. The uncertainty is decomposed into four components which are the signal modelling uncertainty, the background modelling uncertainty, the experimental uncertainty, and the data statistical uncertainty.

The normalised parton-level differential cross-section measured in the fiducial phase-space as a function of the $p_{\textrm{T}}^{l \textrm{non-}Z}$ in the 3$\ell$ channel. The uncertainty is decomposed into four components which are the signal modelling uncertainty, the background modelling uncertainty, the experimental uncertainty, and the data statistical uncertainty.

The normalised parton-level differential cross-section measured in the fiducial phase-space as a function of the $|\Delta \phi (Z, t_{\textrm{lep}})|/\pi$ in the 3$\ell$ channel. The uncertainty is decomposed into four components which are the signal modelling uncertainty, the background modelling uncertainty, the experimental uncertainty, and the data statistical uncertainty.

The normalised parton-level differential cross-section measured in the fiducial phase-space as a function of the $|\Delta \phi (Z, t_{\textrm{lep}})|/\pi$ in the 3$\ell$ channel. The uncertainty is decomposed into four components which are the signal modelling uncertainty, the background modelling uncertainty, the experimental uncertainty, and the data statistical uncertainty.

The normalised parton-level differential cross-section measured in the fiducial phase-space as a function of the $|\Delta y (Z, t_{\textrm{lep}})|$ in the 3$\ell$ channel. The uncertainty is decomposed into four components which are the signal modelling uncertainty, the background modelling uncertainty, the experimental uncertainty, and the data statistical uncertainty.

The normalised parton-level differential cross-section measured in the fiducial phase-space as a function of the $|\Delta y (Z, t_{\textrm{lep}})|$ in the 3$\ell$ channel. The uncertainty is decomposed into four components which are the signal modelling uncertainty, the background modelling uncertainty, the experimental uncertainty, and the data statistical uncertainty.

The normalised parton-level differential cross-section measured in the fiducial phase-space as a function of the $|\Delta \phi (l_{t}^{+}, l_{\bar{t}}^{-})|/\pi$ in the 4$\ell$ channel. The uncertainty is decomposed into four components which are the signal modelling uncertainty, the background modelling uncertainty, the experimental uncertainty, and the data statistical uncertainty.

The normalised parton-level differential cross-section measured in the fiducial phase-space as a function of the $|\Delta \phi (l_{t}^{+}, l_{\bar{t}}^{-})|/\pi$ in the 4$\ell$ channel. The uncertainty is decomposed into four components which are the signal modelling uncertainty, the background modelling uncertainty, the experimental uncertainty, and the data statistical uncertainty.

The normalised parton-level differential cross-section measured in the fiducial phase-space as a function of the $|\Delta \phi (t\bar{t}, Z)|/\pi$ in the 4$\ell$ channel. The uncertainty is decomposed into four components which are the signal modelling uncertainty, the background modelling uncertainty, the experimental uncertainty, and the data statistical uncertainty.

The normalised parton-level differential cross-section measured in the fiducial phase-space as a function of the $|\Delta \phi (t\bar{t}, Z)|/\pi$ in the 4$\ell$ channel. The uncertainty is decomposed into four components which are the signal modelling uncertainty, the background modelling uncertainty, the experimental uncertainty, and the data statistical uncertainty.

The normalised parton-level differential cross-section measured in the fiducial phase-space as a function of the $p_{\textrm{T}}^{t\bar{t}}$ in the 4$\ell$ channel. The uncertainty is decomposed into four components which are the signal modelling uncertainty, the background modelling uncertainty, the experimental uncertainty, and the data statistical uncertainty.

The normalised parton-level differential cross-section measured in the fiducial phase-space as a function of the $p_{\textrm{T}}^{t\bar{t}}$ in the 4$\ell$ channel. The uncertainty is decomposed into four components which are the signal modelling uncertainty, the background modelling uncertainty, the experimental uncertainty, and the data statistical uncertainty.

The normalised parton-level differential cross-section measured in the fiducial phase-space as a function of the absolute value of rapidity of the $Z$ boson in the 3$\ell$+4$\ell$ channel. The uncertainty is decomposed into four components which are the signal modelling uncertainty, the background modelling uncertainty, the experimental uncertainty, and the data statistical uncertainty.

The normalised parton-level differential cross-section measured in the fiducial phase-space as a function of the absolute value of rapidity of the $Z$ boson in the 3$\ell$+4$\ell$ channel. The uncertainty is decomposed into four components which are the signal modelling uncertainty, the background modelling uncertainty, the experimental uncertainty, and the data statistical uncertainty.

The absolute particle-level differential cross-section measured in the fiducial phase-space as a function of the $p_{\textrm{T}}^{l \textrm{non-}Z}$ in the 3$\ell$ channel. The uncertainty is decomposed into four components which are the signal modelling uncertainty, the background modelling uncertainty, the experimental uncertainty, and the data statistical uncertainty.

The absolute particle-level differential cross-section measured in the fiducial phase-space as a function of the $p_{\textrm{T}}^{l \textrm{non-}Z}$ in the 3$\ell$ channel. The uncertainty is decomposed into four components which are the signal modelling uncertainty, the background modelling uncertainty, the experimental uncertainty, and the data statistical uncertainty.

The absolute particle-level differential cross-section measured in the fiducial phase-space as a function of the $|\Delta \phi (Z, t_{\textrm{lep}})|/\pi$ in the 3$\ell$ channel. The uncertainty is decomposed into four components which are the signal modelling uncertainty, the background modelling uncertainty, the experimental uncertainty, and the data statistical uncertainty.

The absolute particle-level differential cross-section measured in the fiducial phase-space as a function of the $|\Delta \phi (Z, t_{\textrm{lep}})|/\pi$ in the 3$\ell$ channel. The uncertainty is decomposed into four components which are the signal modelling uncertainty, the background modelling uncertainty, the experimental uncertainty, and the data statistical uncertainty.

The absolute particle-level differential cross-section measured in the fiducial phase-space as a function of the $|\Delta y (Z, t_{\textrm{lep}})|$ in the 3$\ell$ channel. The uncertainty is decomposed into four components which are the signal modelling uncertainty, the background modelling uncertainty, the experimental uncertainty, and the data statistical uncertainty.

The absolute particle-level differential cross-section measured in the fiducial phase-space as a function of the $|\Delta y (Z, t_{\textrm{lep}})|$ in the 3$\ell$ channel. The uncertainty is decomposed into four components which are the signal modelling uncertainty, the background modelling uncertainty, the experimental uncertainty, and the data statistical uncertainty.

The absolute particle-level differential cross-section measured in the fiducial phase-space as a function of the $|\Delta \phi (l_{t}^{+}, l_{\bar{t}}^{-})|/\pi$ in the 4$\ell$ channel. The uncertainty is decomposed into four components which are the signal modelling uncertainty, the background modelling uncertainty, the experimental uncertainty, and the data statistical uncertainty.

The absolute particle-level differential cross-section measured in the fiducial phase-space as a function of the $|\Delta \phi (l_{t}^{+}, l_{\bar{t}}^{-})|/\pi$ in the 4$\ell$ channel. The uncertainty is decomposed into four components which are the signal modelling uncertainty, the background modelling uncertainty, the experimental uncertainty, and the data statistical uncertainty.

The absolute particle-level differential cross-section measured in the fiducial phase-space as a function of the $|\Delta \phi (t\bar{t}, Z)|/\pi$ in the 4$\ell$ channel. The uncertainty is decomposed into four components which are the signal modelling uncertainty, the background modelling uncertainty, the experimental uncertainty, and the data statistical uncertainty.

The absolute particle-level differential cross-section measured in the fiducial phase-space as a function of the $|\Delta \phi (t\bar{t}, Z)|/\pi$ in the 4$\ell$ channel. The uncertainty is decomposed into four components which are the signal modelling uncertainty, the background modelling uncertainty, the experimental uncertainty, and the data statistical uncertainty.

The absolute particle-level differential cross-section measured in the fiducial phase-space as a function of the $p_{\textrm{T}}^{t\bar{t}}$ in the 4$\ell$ channel. The uncertainty is decomposed into four components which are the signal modelling uncertainty, the background modelling uncertainty, the experimental uncertainty, and the data statistical uncertainty.

The absolute particle-level differential cross-section measured in the fiducial phase-space as a function of the $p_{\textrm{T}}^{t\bar{t}}$ in the 4$\ell$ channel. The uncertainty is decomposed into four components which are the signal modelling uncertainty, the background modelling uncertainty, the experimental uncertainty, and the data statistical uncertainty.

The absolute particle-level differential cross-section measured in the fiducial phase-space as a function of the absolute value of rapidity of the $Z$ boson in the 3$\ell$+4$\ell$ channel. The uncertainty is decomposed into four components which are the signal modelling uncertainty, the background modelling uncertainty, the experimental uncertainty, and the data statistical uncertainty.

The absolute particle-level differential cross-section measured in the fiducial phase-space as a function of the absolute value of rapidity of the $Z$ boson in the 3$\ell$+4$\ell$ channel. The uncertainty is decomposed into four components which are the signal modelling uncertainty, the background modelling uncertainty, the experimental uncertainty, and the data statistical uncertainty.

The normalised particle-level differential cross-section measured in the fiducial phase-space as a function of the $p_{\textrm{T}}^{l \textrm{non-}Z}$ in the 3$\ell$ channel. The uncertainty is decomposed into four components which are the signal modelling uncertainty, the background modelling uncertainty, the experimental uncertainty, and the data statistical uncertainty.

The normalised particle-level differential cross-section measured in the fiducial phase-space as a function of the $p_{\textrm{T}}^{l \textrm{non-}Z}$ in the 3$\ell$ channel. The uncertainty is decomposed into four components which are the signal modelling uncertainty, the background modelling uncertainty, the experimental uncertainty, and the data statistical uncertainty.

The normalised particle-level differential cross-section measured in the fiducial phase-space as a function of the $|\Delta \phi (Z, t_{\textrm{lep}})|/\pi$ in the 3$\ell$ channel. The uncertainty is decomposed into four components which are the signal modelling uncertainty, the background modelling uncertainty, the experimental uncertainty, and the data statistical uncertainty.

The normalised particle-level differential cross-section measured in the fiducial phase-space as a function of the $|\Delta \phi (Z, t_{\textrm{lep}})|/\pi$ in the 3$\ell$ channel. The uncertainty is decomposed into four components which are the signal modelling uncertainty, the background modelling uncertainty, the experimental uncertainty, and the data statistical uncertainty.

The normalised particle-level differential cross-section measured in the fiducial phase-space as a function of the $|\Delta y (Z, t_{\textrm{lep}})|$ in the 3$\ell$ channel. The uncertainty is decomposed into four components which are the signal modelling uncertainty, the background modelling uncertainty, the experimental uncertainty, and the data statistical uncertainty.

The normalised particle-level differential cross-section measured in the fiducial phase-space as a function of the $|\Delta y (Z, t_{\textrm{lep}})|$ in the 3$\ell$ channel. The uncertainty is decomposed into four components which are the signal modelling uncertainty, the background modelling uncertainty, the experimental uncertainty, and the data statistical uncertainty.

The normalised particle-level differential cross-section measured in the fiducial phase-space as a function of the $|\Delta \phi (l_{t}^{+}, l_{\bar{t}}^{-})|/\pi$ in the 4$\ell$ channel. The uncertainty is decomposed into four components which are the signal modelling uncertainty, the background modelling uncertainty, the experimental uncertainty, and the data statistical uncertainty.

The normalised particle-level differential cross-section measured in the fiducial phase-space as a function of the $|\Delta \phi (l_{t}^{+}, l_{\bar{t}}^{-})|/\pi$ in the 4$\ell$ channel. The uncertainty is decomposed into four components which are the signal modelling uncertainty, the background modelling uncertainty, the experimental uncertainty, and the data statistical uncertainty.

The normalised particle-level differential cross-section measured in the fiducial phase-space as a function of the $|\Delta \phi (t\bar{t}, Z)|/\pi$ in the 4$\ell$ channel. The uncertainty is decomposed into four components which are the signal modelling uncertainty, the background modelling uncertainty, the experimental uncertainty, and the data statistical uncertainty.

The normalised particle-level differential cross-section measured in the fiducial phase-space as a function of the $|\Delta \phi (t\bar{t}, Z)|/\pi$ in the 4$\ell$ channel. The uncertainty is decomposed into four components which are the signal modelling uncertainty, the background modelling uncertainty, the experimental uncertainty, and the data statistical uncertainty.

The normalised particle-level differential cross-section measured in the fiducial phase-space as a function of the $p_{\textrm{T}}^{t\bar{t}}$ in the 4$\ell$ channel. The uncertainty is decomposed into four components which are the signal modelling uncertainty, the background modelling uncertainty, the experimental uncertainty, and the data statistical uncertainty.

The normalised particle-level differential cross-section measured in the fiducial phase-space as a function of the $p_{\textrm{T}}^{t\bar{t}}$ in the 4$\ell$ channel. The uncertainty is decomposed into four components which are the signal modelling uncertainty, the background modelling uncertainty, the experimental uncertainty, and the data statistical uncertainty.

The normalised particle-level differential cross-section measured in the fiducial phase-space as a function of the absolute value of rapidity of the $Z$ boson in the 3$\ell$+4$\ell$ channel. The uncertainty is decomposed into four components which are the signal modelling uncertainty, the background modelling uncertainty, the experimental uncertainty, and the data statistical uncertainty.

The normalised particle-level differential cross-section measured in the fiducial phase-space as a function of the absolute value of rapidity of the $Z$ boson in the 3$\ell$+4$\ell$ channel. The uncertainty is decomposed into four components which are the signal modelling uncertainty, the background modelling uncertainty, the experimental uncertainty, and the data statistical uncertainty.

The normalised particle-level differential cross-section measured in the fiducial phase-space as a function of the number of jets in the 3$\ell$ channel. The uncertainty is decomposed into four components which are the signal modelling uncertainty, the background modelling uncertainty, the experimental uncertainty, and the data statistical uncertainty.

The normalised particle-level differential cross-section measured in the fiducial phase-space as a function of the number of jets in the 3$\ell$ channel. The uncertainty is decomposed into four components which are the signal modelling uncertainty, the background modelling uncertainty, the experimental uncertainty, and the data statistical uncertainty.

The normalised particle-level differential cross-section measured in the fiducial phase-space as a function of the number of jets in the 4$\ell$ channel. The uncertainty is decomposed into four components which are the signal modelling uncertainty, the background modelling uncertainty, the experimental uncertainty, and the data statistical uncertainty.

The normalised particle-level differential cross-section measured in the fiducial phase-space as a function of the number of jets in the 4$\ell$ channel. The uncertainty is decomposed into four components which are the signal modelling uncertainty, the background modelling uncertainty, the experimental uncertainty, and the data statistical uncertainty.

The total correlation matrix of the absolute parton-level differential cross-section measured in the fiducial phase-space as a function of the $p_{\textrm{T}}$ of the $Z$ boson in the 3$\ell$+4$\ell$ channel.

The total correlation matrix of the absolute parton-level differential cross-section measured in the fiducial phase-space as a function of the $p_{\textrm{T}}$ of the $Z$ boson in the 3$\ell$+4$\ell$ channel.

The total correlation matrix of the absolute parton-level differential cross-section measured in the fiducial phase-space as a function of the absolute value of rapidity of the $Z$ boson in the 3$\ell$+4$\ell$ channel.

The total correlation matrix of the absolute parton-level differential cross-section measured in the fiducial phase-space as a function of the absolute value of rapidity of the $Z$ boson in the 3$\ell$+4$\ell$ channel.

The total correlation matrix of the absolute parton-level differential cross-section measured in the fiducial phase-space as a function of the $p_{\textrm{T}}^{l \textrm{non-}Z}$ in the 3$\ell$ channel.

The total correlation matrix of the absolute parton-level differential cross-section measured in the fiducial phase-space as a function of the $p_{\textrm{T}}^{l \textrm{non-}Z}$ in the 3$\ell$ channel.

The total correlation matrix of the absolute parton-level differential cross-section measured in the fiducial phase-space as a function of the $|\Delta \phi (Z, t_{\textrm{lep}})|/\pi$ in the 3$\ell$ channel.

The total correlation matrix of the absolute parton-level differential cross-section measured in the fiducial phase-space as a function of the $|\Delta \phi (Z, t_{\textrm{lep}})|/\pi$ in the 3$\ell$ channel.

The total correlation matrix of the absolute parton-level differential cross-section measured in the fiducial phase-space as a function of the $|\Delta y (Z, t_{\textrm{lep}})|/\pi$ in the 3$\ell$ channel.

The total correlation matrix of the absolute parton-level differential cross-section measured in the fiducial phase-space as a function of the $|\Delta y (Z, t_{\textrm{lep}})|/\pi$ in the 3$\ell$ channel.

The total correlation matrix of the absolute parton-level differential cross-section measured in the fiducial phase-space as a function of the $|\Delta \phi (l_{t}^{+}, l_{\bar{t}}^{-})|/\pi$ in the 4$\ell$ channel.

The total correlation matrix of the absolute parton-level differential cross-section measured in the fiducial phase-space as a function of the $|\Delta \phi (l_{t}^{+}, l_{\bar{t}}^{-})|/\pi$ in the 4$\ell$ channel.

The total correlation matrix of the absolute parton-level differential cross-section measured in the fiducial phase-space as a function of the $|\Delta \phi (t\bar{t}, Z)|/\pi$ in the 4$\ell$ channel.

The total correlation matrix of the absolute parton-level differential cross-section measured in the fiducial phase-space as a function of the $|\Delta \phi (t\bar{t}, Z)|/\pi$ in the 4$\ell$ channel.

The total correlation matrix of the absolute parton-level differential cross-section measured in the fiducial phase-space as a function of the $p_{\textrm{T}}^{t\bar{t}}$ in the 4$\ell$ channel.

The total correlation matrix of the absolute parton-level differential cross-section measured in the fiducial phase-space as a function of the $p_{\textrm{T}}^{t\bar{t}}$ in the 4$\ell$ channel.

The total correlation matrix of the normalised parton-level differential cross-section measured in the fiducial phase-space as a function of the $p_{\textrm{T}}$ of the $Z$ boson in the 3$\ell$+4$\ell$ channel.

The total correlation matrix of the normalised parton-level differential cross-section measured in the fiducial phase-space as a function of the $p_{\textrm{T}}$ of the $Z$ boson in the 3$\ell$+4$\ell$ channel.

The total correlation matrix of the normalised parton-level differential cross-section measured in the fiducial phase-space as a function of the absolute value of rapidity of the $Z$ boson in the 3$\ell$+4$\ell$ channel.

The total correlation matrix of the normalised parton-level differential cross-section measured in the fiducial phase-space as a function of the absolute value of rapidity of the $Z$ boson in the 3$\ell$+4$\ell$ channel.

The total correlation matrix of the normalised parton-level differential cross-section measured in the fiducial phase-space as a function of the $p_{\textrm{T}}^{l \textrm{non-}Z}$ in the 3$\ell$ channel.

The total correlation matrix of the normalised parton-level differential cross-section measured in the fiducial phase-space as a function of the $p_{\textrm{T}}^{l \textrm{non-}Z}$ in the 3$\ell$ channel.

The total correlation matrix of the normalised parton-level differential cross-section measured in the fiducial phase-space as a function of the $|\Delta \phi (Z, t_{\textrm{lep}})|/\pi$ in the 3$\ell$ channel.

The total correlation matrix of the normalised parton-level differential cross-section measured in the fiducial phase-space as a function of the $|\Delta \phi (Z, t_{\textrm{lep}})|/\pi$ in the 3$\ell$ channel.

The total correlation matrix of the normalised parton-level differential cross-section measured in the fiducial phase-space as a function of the $|\Delta y (Z, t_{\textrm{lep}})|/\pi$ in the 3$\ell$ channel.

The total correlation matrix of the normalised parton-level differential cross-section measured in the fiducial phase-space as a function of the $|\Delta y (Z, t_{\textrm{lep}})|/\pi$ in the 3$\ell$ channel.

The total correlation matrix of the normalised parton-level differential cross-section measured in the fiducial phase-space as a function of the $|\Delta \phi (l_{t}^{+}, l_{\bar{t}}^{-})|/\pi$ in the 4$\ell$ channel.

The total correlation matrix of the normalised parton-level differential cross-section measured in the fiducial phase-space as a function of the $|\Delta \phi (l_{t}^{+}, l_{\bar{t}}^{-})|/\pi$ in the 4$\ell$ channel.

The total correlation matrix of the normalised parton-level differential cross-section measured in the fiducial phase-space as a function of the $|\Delta \phi (t\bar{t}, Z)|/\pi$ in the 4$\ell$ channel.

The total correlation matrix of the normalised parton-level differential cross-section measured in the fiducial phase-space as a function of the $|\Delta \phi (t\bar{t}, Z)|/\pi$ in the 4$\ell$ channel.

The total correlation matrix of the normalised parton-level differential cross-section measured in the fiducial phase-space as a function of the $p_{\textrm{T}}^{t\bar{t}}$ in the 4$\ell$ channel.

The total correlation matrix of the normalised parton-level differential cross-section measured in the fiducial phase-space as a function of the $p_{\textrm{T}}^{t\bar{t}}$ in the 4$\ell$ channel.

The total correlation matrix of the absolute particle-level differential cross-section measured in the fiducial phase-space as a function of the $p_{\textrm{T}}$ of the $Z$ boson in the 3$\ell$+4$\ell$ channel.

The total correlation matrix of the absolute particle-level differential cross-section measured in the fiducial phase-space as a function of the $p_{\textrm{T}}$ of the $Z$ boson in the 3$\ell$+4$\ell$ channel.

The total correlation matrix of the absolute particle-level differential cross-section measured in the fiducial phase-space as a function of the absolute value of rapidity of the $Z$ boson in the 3$\ell$+4$\ell$ channel.

The total correlation matrix of the absolute particle-level differential cross-section measured in the fiducial phase-space as a function of the absolute value of rapidity of the $Z$ boson in the 3$\ell$+4$\ell$ channel.

The total correlation matrix of the absolute particle-level differential cross-section measured in the fiducial phase-space as a function of the $p_{\textrm{T}}^{l \textrm{non-}Z}$ in the 3$\ell$ channel.

The total correlation matrix of the absolute particle-level differential cross-section measured in the fiducial phase-space as a function of the $p_{\textrm{T}}^{l \textrm{non-}Z}$ in the 3$\ell$ channel.

The total correlation matrix of the absolute particle-level differential cross-section measured in the fiducial phase-space as a function of the $|\Delta \phi (Z, t_{\textrm{lep}})|/\pi$ in the 3$\ell$ channel.

The total correlation matrix of the absolute particle-level differential cross-section measured in the fiducial phase-space as a function of the $|\Delta \phi (Z, t_{\textrm{lep}})|/\pi$ in the 3$\ell$ channel.

The total correlation matrix of the absolute particle-level differential cross-section measured in the fiducial phase-space as a function of the $|\Delta y (Z, t_{\textrm{lep}})|/\pi$ in the 3$\ell$ channel.

The total correlation matrix of the absolute particle-level differential cross-section measured in the fiducial phase-space as a function of the $|\Delta y (Z, t_{\textrm{lep}})|/\pi$ in the 3$\ell$ channel.

The total correlation matrix of the absolute particle-level differential cross-section measured in the fiducial phase-space as a function of the number of jets in the 3$\ell$ channel.

The total correlation matrix of the absolute particle-level differential cross-section measured in the fiducial phase-space as a function of the number of jets in the 3$\ell$ channel.

The total correlation matrix of the absolute particle-level differential cross-section measured in the fiducial phase-space as a function of the $|\Delta \phi (l_{t}^{+}, l_{\bar{t}}^{-})|/\pi$ in the 4$\ell$ channel.

The total correlation matrix of the absolute particle-level differential cross-section measured in the fiducial phase-space as a function of the $|\Delta \phi (l_{t}^{+}, l_{\bar{t}}^{-})|/\pi$ in the 4$\ell$ channel.

The total correlation matrix of the absolute particle-level differential cross-section measured in the fiducial phase-space as a function of the $|\Delta \phi (t\bar{t}, Z)|/\pi$ in the 4$\ell$ channel.

The total correlation matrix of the absolute particle-level differential cross-section measured in the fiducial phase-space as a function of the $|\Delta \phi (t\bar{t}, Z)|/\pi$ in the 4$\ell$ channel.

The total correlation matrix of the absolute particle-level differential cross-section measured in the fiducial phase-space as a function of the $p_{\textrm{T}}^{t\bar{t}}$ in the 4$\ell$ channel.

The total correlation matrix of the absolute particle-level differential cross-section measured in the fiducial phase-space as a function of the $p_{\textrm{T}}^{t\bar{t}}$ in the 4$\ell$ channel.

The total correlation matrix of the absolute particle-level differential cross-section measured in the fiducial phase-space as a function of the number of jets in the 4$\ell$ channel.

The total correlation matrix of the absolute particle-level differential cross-section measured in the fiducial phase-space as a function of the number of jets in the 4$\ell$ channel.

The total correlation matrix of the normalised particle-level differential cross-section measured in the fiducial phase-space as a function of the $p_{\textrm{T}}$ of the $Z$ boson in the 3$\ell$+4$\ell$ channel.

The total correlation matrix of the normalised particle-level differential cross-section measured in the fiducial phase-space as a function of the $p_{\textrm{T}}$ of the $Z$ boson in the 3$\ell$+4$\ell$ channel.

The total correlation matrix of the normalised particle-level differential cross-section measured in the fiducial phase-space as a function of the absolute value of rapidity of the $Z$ boson in the 3$\ell$+4$\ell$ channel.

The total correlation matrix of the normalised particle-level differential cross-section measured in the fiducial phase-space as a function of the absolute value of rapidity of the $Z$ boson in the 3$\ell$+4$\ell$ channel.

The total correlation matrix of the normalised particle-level differential cross-section measured in the fiducial phase-space as a function of the $p_{\textrm{T}}^{l \textrm{non-}Z}$ in the 3$\ell$ channel.

The total correlation matrix of the normalised particle-level differential cross-section measured in the fiducial phase-space as a function of the $p_{\textrm{T}}^{l \textrm{non-}Z}$ in the 3$\ell$ channel.

The total correlation matrix of the normalised particle-level differential cross-section measured in the fiducial phase-space as a function of the $|\Delta \phi (Z, t_{\textrm{lep}})|/\pi$ in the 3$\ell$ channel.

The total correlation matrix of the normalised particle-level differential cross-section measured in the fiducial phase-space as a function of the $|\Delta \phi (Z, t_{\textrm{lep}})|/\pi$ in the 3$\ell$ channel.

The total correlation matrix of the normalised particle-level differential cross-section measured in the fiducial phase-space as a function of the $|\Delta y (Z, t_{\textrm{lep}})|/\pi$ in the 3$\ell$ channel.

The total correlation matrix of the normalised particle-level differential cross-section measured in the fiducial phase-space as a function of the $|\Delta y (Z, t_{\textrm{lep}})|/\pi$ in the 3$\ell$ channel.

The total correlation matrix of the normalised particle-level differential cross-section measured in the fiducial phase-space as a function of the number of jets in the 3$\ell$ channel.

The total correlation matrix of the normalised particle-level differential cross-section measured in the fiducial phase-space as a function of the number of jets in the 3$\ell$ channel.

The total correlation matrix of the normalised particle-level differential cross-section measured in the fiducial phase-space as a function of the $|\Delta \phi (l_{t}^{+}, l_{\bar{t}}^{-})|/\pi$ in the 4$\ell$ channel.

The total correlation matrix of the normalised particle-level differential cross-section measured in the fiducial phase-space as a function of the $|\Delta \phi (l_{t}^{+}, l_{\bar{t}}^{-})|/\pi$ in the 4$\ell$ channel.

The total correlation matrix of the normalised particle-level differential cross-section measured in the fiducial phase-space as a function of the $|\Delta \phi (t\bar{t}, Z)|/\pi$ in the 4$\ell$ channel.

The total correlation matrix of the normalised particle-level differential cross-section measured in the fiducial phase-space as a function of the $|\Delta \phi (t\bar{t}, Z)|/\pi$ in the 4$\ell$ channel.

The total correlation matrix of the normalised particle-level differential cross-section measured in the fiducial phase-space as a function of the $p_{\textrm{T}}^{t\bar{t}}$ in the 4$\ell$ channel.

The total correlation matrix of the normalised particle-level differential cross-section measured in the fiducial phase-space as a function of the $p_{\textrm{T}}^{t\bar{t}}$ in the 4$\ell$ channel.

The total correlation matrix of the normalised particle-level differential cross-section measured in the fiducial phase-space as a function of the number of jets in the 4$\ell$ channel.

The total correlation matrix of the normalised particle-level differential cross-section measured in the fiducial phase-space as a function of the number of jets in the 4$\ell$ channel.

Differential fiducial cross section (without the acceptance correction) for the DY process measured in the muon channel, as a function of $p_{\textrm{T}}$ for $15<m_{\mu\mu}<60$ GeV. The quoted error is the quadratic sum of the statistical and systematic uncertainties. The global normalisation uncertainty of 3.5% is listed separately.

Differential fiducial cross section (without the acceptance correction) for the DY process measured in the muon channel, as a function of $p_{\textrm{T}}$ for $60<m_{\mu\mu}<120$ GeV. The quoted error is the quadratic sum of the statistical and systematic uncertainties. The global normalisation uncertainty of 3.5% is listed separately.

Differential fiducial cross section (without the acceptance correction) for the DY process measured in the muon channel, as a function of $\phi^*$ for $15<m_{\mu\mu}<60$ GeV. The quoted error is the quadratic sum of the statistical and systematic uncertainties. The global normalisation uncertainty of 3.5% is listed separately.

Differential fiducial cross section (without the acceptance correction) for the DY process measured in the muon channel, as a function of $\phi^*$ for $60<m_{\mu\mu}<120$ GeV. The quoted error is the quadratic sum of the statistical and systematic uncertainties. The global normalisation uncertainty of 3.5% is listed separately.

Integrated fiducial cross section for the DY process measured in the muon channel. The quoted error is the quadratic sum of the statistical and systematic uncertainties. The global normalisation uncertainty of 3.5% is listed separately.

Differential cross section for the DY process measured in the muon channel, as a function of dimuon invariant mass. The quoted error is the quadratic sum of the statistical and systematic uncertainties. The global normalisation uncertainty of 3.5% is listed separately.

Differential cross section for the DY process measured in the muon channel, as a function of rapidity in the centre-of-mass frame for $15<m_{\mu\mu}<60$ GeV. The quoted error is the quadratic sum of the statistical and systematic uncertainties. The global normalisation uncertainty of 3.5% is listed separately.

Differential cross section for the DY process measured in the muon channel, as a function of rapidity in the centre-of-mass frame for $60<m_{\mu\mu}<120$ GeV. The quoted error is the quadratic sum of the statistical and systematic uncertainties. The global normalisation uncertainty of 3.5% is listed separately.

Differential cross section for the DY process measured in the muon channel, as a function of $p_{\textrm{T}}$ for $15<m_{\mu\mu}<60$ GeV. The quoted error is the quadratic sum of the statistical and systematic uncertainties. The global normalisation uncertainty of 3.5% is listed separately.

Differential cross section for the DY process measured in the muon channel, as a function of $p_{\textrm{T}}$ for $60<m_{\mu\mu}<120$ GeV. The quoted error is the quadratic sum of the statistical and systematic uncertainties. The global normalisation uncertainty of 3.5% is listed separately.

Differential cross section for the DY process measured in the muon channel, as a function of $\phi^*$ for $15<m_{\mu\mu}<60$ GeV. The quoted error is the quadratic sum of the statistical and systematic uncertainties. The global normalisation uncertainty of 3.5% is listed separately.

Differential cross section for the DY process measured in the muon channel, as a function of $\phi^*$ for $60<m_{\mu\mu}<120$ GeV. The quoted error is the quadratic sum of the statistical and systematic uncertainties. The global normalisation uncertainty of 3.5% is listed separately.

Integrated cross section for the DY process measured in the muon channel. The quoted error is the quadratic sum of the statistical and systematic uncertainties. The global normalisation uncertainty of 3.5% is listed separately.

Forward-backward ratios for $15<m_{\mu\mu}<60$ GeV. Quoted uncertainties are the quadratic sum of the statistical and systematic uncertainties.

Forward-backward ratios for $60<m_{\mu\mu}<120$ GeV. Quoted uncertainties are the quadratic sum of the statistical and systematic uncertainties.

Correlation matrix for the systematic uncertainties, excluding integrated luminosity, as a function of dimuon invariant mass.

Correlation matrix for the systematic uncertainties, excluding integrated luminosity, as a function of rapidity in the centre-of-mass frame for $15<m_{\mu\mu}<60$ GeV.

Correlation matrix for the systematic uncertainties, excluding integrated luminosity, as a function of rapidity in the centre-of-mass frame for $60<m_{\mu\mu}<120$ GeV.

Correlation matrix for the systematic uncertainties, excluding integrated luminosity, as a function of $p_{\textrm{T}}$ for $15<m_{\mu\mu}<60$ GeV.

Correlation matrix for the systematic uncertainties, excluding integrated luminosity, as a function of $p_{\textrm{T}}$ for $60<m_{\mu\mu}<120$ GeV.

Correlation matrix for the systematic uncertainties, excluding integrated luminosity, as a function of $\phi^*$ for $15<m_{\mu\mu}<60$ GeV.

Correlation matrix for the systematic uncertainties, excluding integrated luminosity, as a function of $\phi^*$ for $60<m_{\mu\mu}<120$ GeV.

Correlation matrix for the systematic uncertainties, excluding integrated luminosity, as a function of rapidity in the centre-of-mass frame, putting together low and high masses. Each bin is labeled lm_i or hm_i, where lm stands for $15<m_{\mu\mu}<60$ GeV, hm stands for $60<m_{\mu\mu}<120$ GeV, and i is the corresponding bin number in rapidity in the centre-of-mass frame. The global normalisation uncertainty of 3.5% is not included. Data from Figure [3, 4, 4].

Correlation matrix for the systematic uncertainties, excluding integrated luminosity, as a function of $p_{\textrm{T}}$, putting together low and high masses. Each bin is labeled lm_i or hm_i, where lm stands for $15<m_{\mu\mu}<60$ GeV, hm stands for $60<m_{\mu\mu}<120$ GeV, and i is the corresponding bin number in $p_{\textrm{T}}$. The global normalisation uncertainty of 3.5% is not included. Data from Figure [3, 4, 4].

Correlation matrix for the systematic uncertainties, excluding integrated luminosity, as a function of $\phi^*$, putting together low and high masses. Each bin is labeled lm_i or hm_i, where lm stands for $15<m_{\mu\mu}<60$ GeV, hm stands for $60<m_{\mu\mu}<120$ GeV, and i is the corresponding bin number in $\phi^*$. The global normalisation uncertainty of 3.5% is not included. Data from Figure [3, 4, 4].

95% CL limits on the Lagrangian approach (a.k.a. LEP parametrization) parameters in Wgamma events.

Experimental uncertainties affecting transfer factors in the analysis that is used to estimate the W background in the signal region (SR). The number of W boson events are denoted as $W_{SR}$ for the SR and in analogy as $W_{\mu\nu}$ ($W_{e\nu}$) for the single-muon (single-electron) control region (CR).

Uncertainties assigned to the simulation based processes in SR and CRs.

Uncertainties assigned to the simulation based processes in SR and CRs.

Cross sections (fb) at large Z $p_{T}$ values in the Z -> $\ell\ell$ and Z -> $\nu\nu$ channels, and their combination. The theoretical predictions from Madgraph at NLO in QCD and corrected to NLO in electroweak using the NNPDF 3.0 are also reported. With the exception of the largest Z $p_{T}$ bin, the statistical uncertainties in the measurements are much smaller than the systematic uncertainties. Both measurements and predictions correspond to $\sigma \mathcal{B}(Z -> \ell\ell)$, where $\sigma$ is the total fiducial cross section, $\mathcal{B}$ is the branching fraction, and $\ell$ is a charged lepton. The $Z -> \nu\nu$ measurement corresponds to $\sigma \mathcal{B}(Z -> \ell\ell)/\mathcal{B}(Z -> \nu\nu)$.

Cross sections (fb) at large Z $p_{T}$ values in the Z -> $\ell\ell$ and Z -> $\nu\nu$ channels, and their combination. The theoretical predictions from Madgraph at NLO in QCD and corrected to NLO in electroweak using the NNPDF 3.0 are also reported. With the exception of the largest Z $p_{T}$ bin, the statistical uncertainties in the measurements are much smaller than the systematic uncertainties. Both measurements and predictions correspond to $\sigma \mathcal{B}(Z -> \ell\ell)$, where $\sigma$ is the total fiducial cross section, $\mathcal{B}$ is the branching fraction, and $\ell$ is a charged lepton. The $Z -> \nu\nu$ measurement corresponds to $\sigma \mathcal{B}(Z -> \ell\ell)/\mathcal{B}(Z -> \nu\nu)$.

Cross sections normalized to the total cross section measurements at high Z $p_{T}$ values in the Z -> $\ell\ell$ and Z -> $\nu\nu$ channels, and in their combination. The uncertainty includes both statistical and systematic uncertainties.

Cross sections normalized to the total cross section measurements at high Z $p_{T}$ values in the Z -> $\ell\ell$ and Z -> $\nu\nu$ channels, and in their combination. The uncertainty includes both statistical and systematic uncertainties.

The relative statistical and systematic uncertainties from various sources for the absolute cross section measurements in bins of Z $p_{T}$ on charged leptons.

The relative statistical and systematic uncertainties from various sources for the absolute cross section measurements in bins of Z $p_{T}$ on charged leptons.

The relative statistical and systematic uncertainties from various sources for the relative cross section measurements in bins of Z $p_{T}$ on charged leptons.

The relative statistical and systematic uncertainties from various sources for the relative cross section measurements in bins of Z $p_{T}$ on charged leptons.

The relative statistical and systematic uncertainties from various sources for the absolute cross section measurements in bins of Z $p_{T}$ on neutrinos.

The relative statistical and systematic uncertainties from various sources for the absolute cross section measurements in bins of Z $p_{T}$ on neutrinos.

The relative statistical and systematic uncertainties from various sources for the relative cross section measurements in bins of Z $p_{T}$ on neutrinos.

The relative statistical and systematic uncertainties from various sources for the relative cross section measurements in bins of Z $p_{T}$ on neutrinos.

The relative statistical and systematic uncertainties from various sources for the absolute cross section measurements in bins of Z $p_{T}$ on neutrinos and charged leptons.

The relative statistical and systematic uncertainties from various sources for the absolute cross section measurements in bins of Z $p_{T}$ on neutrinos and charged leptons.

The relative statistical and systematic uncertainties from various sources for the relative cross section measurements in bins of Z $p_{T}$ on neutrinos and charged leptons.

The relative statistical and systematic uncertainties from various sources for the relative cross section measurements in bins of Z $p_{T}$ on neutrinos and charged leptons.

Cross section measurements in bins of Z $p_{T}$ on neutrinos and charged leptons.

Cross section measurements in bins of Z $p_{T}$ on neutrinos and charged leptons.

Data from paper Table 4. Observed and expected exclusion limits for the aQGC parameters at 95% CL, without any form factors.

Data from paper Fig.4. Observed four lepton invariant mass distribution. The last bin includes overflow.

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.

Exclusion limits at 95% CL for direct gluino pair production, where the gluinos decay to light-flavor quarks and a $\tilde{\chi}_1^\pm$ that decays to $W^\pm\tilde{\chi}_1^0$. Signal cross sections are calculated at approximately NNLO+NNLL order in $\alpha_S$, assuming unity branching fraction to $q_i\bar{q}_j W\pm\tilde{\chi}_1^0$.

Exclusion limits at 95% CL for direct gluino pair production, where the gluinos decay to bottom quarks ($\tilde{g}\to b\bar{b}\tilde{\chi}_1^0$). Signal cross sections are calculated at approximately NNLO+NNLL order in $\alpha_S$, assuming unity branching fraction for the given decay.

Exclusion limits at 95% CL for direct gluino pair production, where the gluinos decay to top quarks ($\tilde{g}\to t\bar{t}\tilde{\chi}_1^0$). Signal cross sections are calculated at approximately NNLO+NNLL order in $\alpha_S$, assuming unity branching fraction for the given decay.

Exclusion limits at 95% CL for light-flavor squark pair production, where the squarks decay to $q\tilde{\chi}_1^0$. Signal cross sections are calculated at approximately NNLO+NNLL order in $\alpha_S$, assuming unity branching fraction for the given decay and 1-fold degeneracy in the light-flavor squarks (corresponding to the inner set of curves in the limit plot). To get the theory cross section for other N-fold degeneracy assumptions (e.g. 8-fold for the outer curves in the limit plot), just multiply by N.

Exclusion limits at 95% CL for bottom squark pair production, where the squarks decay to $b\tilde{\chi}_1^0$. Signal cross sections are calculated at approximately NNLO+NNLL order in $\alpha_S$, assuming unity branching fraction for the given decay.

Exclusion limits at 95% CL for top squark pair production, where the squarks decay to $t\tilde{\chi}_1^0$. Signal cross sections are calculated at approximately NNLO+NNLL order in $\alpha_S$, assuming unity branching fraction for the given decay.

Exclusion limits at 95% CL for top squark pair production, where the squarks decay to $b\tilde{\chi}_1^\pm$ and the $\tilde{\chi}_1^0$ decay to $W^\pm\tilde{\chi}_1^0$. Signal cross sections are calculated at approximately NNLO+NNLL order in $\alpha_S$, assuming unity branching fraction for the given decay.

Exclusion limits at 95% CL for top squark pair production, where the squarks decay either to $b\tilde{\chi}_1^\pm\to bW^\pm\tilde{\chi}_1^0$ or to $t\tilde{\chi}_1^0$. Signal cross sections are calculated at approximately NNLO+NNLL order in $\alpha_S$, assuming unity branching fraction for the given decay.

Exclusion limits at 95% CL for top squark pair production, where the squarks decay to $c\tilde{\chi}_1^0$. Signal cross sections are calculated at approximately NNLO+NNLL order in $\alpha_S$, assuming unity branching fraction for the given decay.

Exclusion limits at 95% CL for the mono-$\phi$ model, in which a resonantly-produced colored scalar decays to a massive Dirac fermion and a quark. Signal cross sections are calculated at leading order in $\alpha_S$, assuming unity branching fraction for the given decay.

Cross section limits for $\mathrm{LQ}\to\mathrm{q}\nu$, where $q=u,\,d,\,s,\,\mathrm{or}\,c$. Limits are at the 95% confidence level. Theory cross sections are LO for vector LQ, and NLO for scalar LQ. Branching ratio is assumed to be 100% to $\mathrm{q}\nu$.

Cross section limits for $\mathrm{LQ}\to\mathrm{b}\nu$. Limits are at the 95% confidence level. Theory cross sections are LO for vector LQ, and NLO for scalar LQ. Branching ratio is assumed to be 100% to $\mathrm{b}\nu$.

Cross section limits for $\mathrm{LQ}\to\mathrm{t}\nu$. Limits are at the 95% confidence level. Theory cross sections are LO for vector LQ, and NLO for scalar LQ. Branching ratios are assumed to be $\mathcal{B}(\mathrm{LQ}\to\mathrm{t}\nu)=1-\beta$, and $\mathcal{B}(\mathrm{LQ}\to\mathrm{b}\tau)=\beta$.

Predictions and observations for monojet signal regions

Predictions and observations for signal regions with $250 \leq H_\mathrm{T} < 450$ GeV

Predictions and observations for signal regions with $450 \leq H_\mathrm{T} < 575$ GeV and $N_\mathrm{j}<7$

Predictions and observations for signal regions with $450 \leq H_\mathrm{T} < 575$ GeV and $N_\mathrm{j}\geq7$

Predictions and observations for signal regions with $575 \leq H_\mathrm{T} < 1200$ GeV and $N_\mathrm{j}^\mathrm{hi}<4$

Predictions and observations for signal regions with $575 \leq H_\mathrm{T} < 1200$ GeV and $4\leq N_\mathrm{j}^\mathrm{hi}<7$

Predictions and observations for signal regions with $575 \leq H_\mathrm{T} < 1200$ GeV and $N_\mathrm{j}\geq7$

Predictions and observations for signal regions with $1200 \leq H_\mathrm{T} < 1500$ GeV and $N_\mathrm{j}^\mathrm{hi}<4$

Predictions and observations for signal regions with $1200 \leq H_\mathrm{T} < 1500$ GeV and $4\leq N_\mathrm{j}^\mathrm{hi}<7$

Predictions and observations for signal regions with $1200 \leq H_\mathrm{T} < 1500$ GeV and $N_\mathrm{j}\geq7$

Predictions and observations for signal regions with $H_\mathrm{T} \geq 1500$ GeV and $N_\mathrm{j}<7$

Predictions and observations for signal regions with $H_\mathrm{T} \geq 1500$ GeV and $N_\mathrm{j}\geq7$

Covariance matrix for the 282 signal regions of the inclusive $M_\mathrm{T2}$ search

Correlation matrix for the 282 signal regions of the inclusive $M_\mathrm{T2}$ search

Bin number definitions for the $M_\mathrm{T2}$ covariance and correlation matrices

Exclusion limits at 95% CL for direct gluino pair production, where the gluinos decay to light-flavor quarks and either the lightest neutralino, or the lightest chargino, and the chargino is long-lived with $c\tau_0 = 10$ cm and mass O(100) MeV greater than the neutralino's mass. Signal cross sections are calculated at approximately NNLO+NNLL order in $\alpha_S$.

Exclusion limits at 95% CL for direct gluino pair production, where the gluinos decay to light-flavor quarks and either the lightest neutralino, or the lightest chargino, and the chargino is long-lived with $c\tau_0 = 50$ cm and mass O(100) MeV greater than the neutralino's mass. Signal cross sections are calculated at approximately NNLO+NNLL order in $\alpha_S$.

Exclusion limits at 95% CL for direct gluino pair production, where the gluinos decay to light-flavor quarks and either the lightest neutralino, or the lightest chargino, and the chargino is long-lived with $c\tau_0 = 200$ cm and mass O(100) MeV greater than the neutralino's mass. Signal cross sections are calculated at approximately NNLO+NNLL order in $\alpha_S$.

Exclusion limits at 95% CL for direct light squark pair production, where the squarks decay to light-flavor quarks and either the lightest neutralino, or the lightest chargino, and the chargino is long-lived with $c\tau_0 = 10$ cm and mass O(100) MeV greater than the neutralino's mass. Signal cross sections are calculated at approximately NNLO+NNLL order in $\alpha_S$, for a single light squark.

Exclusion limits at 95% CL for direct light squark pair production, where the squarks decay to light-flavor quarks and either the lightest neutralino, or the lightest chargino, and the chargino is long-lived with $c\tau_0 = 50$ cm and mass O(100) MeV greater than the neutralino's mass. Signal cross sections are calculated at approximately NNLO+NNLL order in $\alpha_S$, for a single light squark.

Exclusion limits at 95% CL for direct light squark pair production, where the squarks decay to light-flavor quarks and either the lightest neutralino, or the lightest chargino, and the chargino is long-lived with $c\tau_0 = 200$ cm and mass O(100) MeV greater than the neutralino's mass. Signal cross sections are calculated at approximately NNLO+NNLL order in $\alpha_S$, for a single light squark.

Exclusion limits at 95% CL for direct stop pair production, where the stops decay to either a top and the lightest neutralino, or a bottom and the lightest chargino, and the chargino is long-lived with $c\tau_0 = 10$ cm and mass O(100) MeV greater than the neutralino's mass. Signal cross sections are calculated at approximately NNLO+NNLL order in $\alpha_S$.

Exclusion limits at 95% CL for direct stop pair production, where the stops decay to either a top and the lightest neutralino, or a bottom and the lightest chargino, and the chargino is long-lived with $c\tau_0 = 50$ cm and mass O(100) MeV greater than the neutralino's mass. Signal cross sections are calculated at approximately NNLO+NNLL order in $\alpha_S$.

Exclusion limits at 95% CL for direct stop pair production, where the stops decay to either a top and the lightest neutralino, or a bottom and the lightest chargino, and the chargino is long-lived with $c\tau_0 = 200$ cm and mass O(100) MeV greater than the neutralino's mass. Signal cross sections are calculated at approximately NNLO+NNLL order in $\alpha_S$.

The maximum chargino mass excluded at 95% CL for direct gluino pair production, where the gluinos decay to light-flavor quarks and either the lightest neutralino, or the lightest chargino, and the chargino mass is O(100) MeV greater than the neutralino's mass. The chargino's lifetime is varied from $c\tau_{0} = 1$ to 2000 cm while the gluino mass is fixed to 1900 GeV. Signal cross sections are calculated at approximately NNLO+NNLL order in $\alpha_S$. If all kinematically allowed chargino masses are excluded, the curves, including 68 and 95% expected, tend to overlap. At short decay lengths, horizontal exclusion lines are obtained from the inclusive analysis, as this is not affected by track reconstruction inefficiencies, which may arise when the chargino decays before the CMS tracker, and therefore shows better sensitivity to scenarios with very small lifetime compared to the disappearing track search, based on median expected limits.

The maximum chargino mass excluded at 95% CL for direct squark pair production, where the squarks decay to light-flavor quarks and either the lightest neutralino, or the lightest chargino, and the chargino mass is O(100) MeV greater than the neutralino's mass. The chargino's lifetime is varied from $c\tau_{0} = 1$ to 2000 cm while the squark mass is fixed to 900 GeV. Signal cross sections are calculated at approximately NNLO+NNLL order in $\alpha_S$, for a single light squark. If all kinematically allowed chargino masses are excluded, the curves, including 68 and 95% expected, tend to overlap. At short decay lengths, horizontal exclusion lines are obtained from the inclusive analysis, as this is not affected by track reconstruction inefficiencies, which may arise when the chargino decays before the CMS tracker, and therefore shows better sensitivity to scenarios with very small lifetime compared to the disappearing track search, based on median expected limits.

The maximum chargino mass excluded at 95% CL for direct squark pair production, where the squarks decay to light-flavor quarks and either the lightest neutralino, or the lightest chargino, and the chargino mass is O(100) MeV greater than the neutralino's mass. The chargino's lifetime is varied from $c\tau_{0} = 1$ to 2000 cm while the squark mass is fixed to 1500 GeV. Signal cross sections are calculated at approximately NNLO+NNLL order in $\alpha_S$, and the eight light squarks' masses are assumed to be degenerate. If all kinematically allowed chargino masses are excluded, the curves, including 68 and 95% expected, tend to overlap. At short decay lengths, horizontal exclusion lines are obtained from the inclusive analysis, as this is not affected by track reconstruction inefficiencies, which may arise when the chargino decays before the CMS tracker, and therefore shows better sensitivity to scenarios with very small lifetime compared to the disappearing track search, based on median expected limits.

Exclusion limits at 95% CL for direct stop pair production, where the stops decay to either a top and the lightest neutralino, or a bottom and the lightest chargino, and the chargino mass is O(100) MeV greater than the neutralino's mass. The chargino's lifetime is varied from $c\tau_{0} = 1$ to 2000 cm while the stop mass is fixed to 1000 GeV. Signal cross sections are calculated at approximately NNLO+NNLL order in $\alpha_S$. If all kinematically allowed chargino masses are excluded, the curves, including 68 and 95% expected, tend to overlap. At short decay lengths, horizontal exclusion lines are obtained from the inclusive analysis, as this is not affected by track reconstruction inefficiencies, which may arise when the chargino decays before the CMS tracker, and therefore shows better sensitivity to scenarios with very small lifetime compared to the disappearing track search, based on median expected limits.

Exclusion limits at 95% CL for direct gluino pair production, where the gluinos decay to light-flavor quarks and either the lightest neutralino, or the lightest chargino, and the chargino mass is O(100) MeV greater than the neutralino's mass. The chargino's lifetime is varied from $c\tau_{0} = 5$ to 1000 cm while the gluino mass is fixed to 1600 GeV and the neutralino's mass is fixed to 1575 GeV. Signal cross sections are calculated at approximately NNLO+NNLL order in $\alpha_S$.

Exclusion limits at 95% CL for direct squark pair production, where the squarks decay to light-flavor quarks and either the lightest neutralino, or the lightest chargino, and the chargino mass is O(100) MeV greater than the neutralino's mass. The chargino's lifetime is varied from $c\tau_{0} = 5$ to 1000 cm while the squark mass is fixed to 2000 GeV and the neutralino's mass is fixed to 1000 GeV. Signal cross sections are calculated at approximately NNLO+NNLL order in $\alpha_S$, and the eight light squarks' masses are assumed to be degenerate.

Exclusion limits at 95% CL for direct stop pair production, where the stops decay to either a top and the lightest neutralino, or a bottom and the lightest chargino, and the chargino mass is O(100) MeV greater than the neutralino's mass. The chargino's lifetime is varied from $c\tau_{0} = 5$ to 1000 cm while the stop mass is fixed to 1100 GeV and the neutralino's mass is fixed to 1000 GeV. Signal cross sections are calculated at approximately NNLO+NNLL order in $\alpha_S$.

Predictions and observations for 2016 disappearing track signal regions

Predictions and observations for 2017-2018 pixel track signal regions

Predictions and observations for 2017-2018 medium (M) and long (L) length track signal regions

Covariance matrix for the 68 signal regions of the disappearing tracks $M_\mathrm{T2}$ search

Correlation matrix for the 68 signal regions of the disappearing tracks $M_\mathrm{T2}$ search

Summary of typical systematic uncertainties for the signal.

Numbers of observed events for all signal regions, including predicted background contributions and expected signal yields.

Limits on anomalous quartic couplings at 95% CL.

Expected and observed 95% CL upper limits on the cross section times the branching ratio sigma(P P --> W+- a)B(a --> W+- W-+) as a function of ALP mass

Expected and observed 95% CL upper limits on the photophobic ALP model parameter 1/f_a as a function of ALP mass.

Observed signal strength, and observed signal cross section.

The theoretical EWK WZ cross section in the tight EWK fiducial region. Can be multiplied by the measured EWK WZ scale factor to find the measured EWK WZ cross section.

aQGC limits on effective field theory parameters in EWK WZ events

The differential measurement of W + D*(2010) as a function of the absolute peudorapidity of the muon originating from the W boson.

The differential measurement of W+ + D*(2010)- as a function of the absolute peudorapidity of the muon originating from the W boson.

The differential measurement of W- + D*(2010)+ as a function of the absolute peudorapidity of the muon originating from the W boson.

Distributions of $\min[M(\ell,\mathrm{j})]$ in the W+jets control region, for 0 W-tagged jet category. Example signal distributions are also shown. The background distributions correspond to background-only fit to data while signal distributions are before the fit to data. Electron and muon event samples are combined. The last bin includes overflow events and its content is divided by the bin width. The distributions in each category have variable-size bins, chosen so that the statistical uncertainty in the total background in each bin is less than 30%. The lower panel in each plot shows the difference between the observed and the predicted numbers of events in that bin divided by the total uncertainty. The total uncertainty is calculated as the sum in quadrature of the statistical uncertainty in the observed measurement and the statistical and systematic uncertainties in the background-only fit to data.

Distributions of $\min[M(\ell,\mathrm{j})]$ in the W+jets control region, for 1 or more W-tagged jet category. Example signal distributions are also shown. The background distributions correspond to background-only fit to data while signal distributions are before the fit to data. Electron and muon event samples are combined. The last bin includes overflow events and its content is divided by the bin width. The distributions in each category have variable-size bins, chosen so that the statistical uncertainty in the total background in each bin is less than 30%. The lower panel in each plot shows the difference between the observed and the predicted numbers of events in that bin divided by the total uncertainty. The total uncertainty is calculated as the sum in quadrature of the statistical uncertainty in the observed measurement and the statistical and systematic uncertainties in the background-only fit to data.

Distributions of $\min[M(\ell,\mathrm{b})]$ in events with 0 t-tagged jets, 0 W-tagged jets, and 1 b-tagged jets for the combined electron and muon samples in the signal region. Example signal distributions are also shown. The background distributions correspond to the background-only fit to data, while signal distributions are before the fit to data. The last bin includes overflow events and its content is divided by the bin width. The distributions in each category have variable-size bins, chosen so that the statistical uncertainty in the total background in each bin is less than 30%. The lower panel in each plot shows the difference between the observed and the predicted numbers of events in that bin divided by the total uncertainty. The total uncertainty is calculated as the sum in quadrature of the statistical uncertainty in the observed measurement and the statistical and systematic uncertainties in the background-only fit to data.

Distributions of $\min[M(\ell,\mathrm{b})]$ in events with 0 t-tagged jets, 0 W-tagged jets, and 2 or more b-tagged jets for the combined electron and muon samples in the signal region. Example signal distributions are also shown. The background distributions correspond to the background-only fit to data, while signal distributions are before the fit to data. The last bin includes overflow events and its content is divided by the bin width. The distributions in each category have variable-size bins, chosen so that the statistical uncertainty in the total background in each bin is less than 30%. The lower panel in each plot shows the difference between the observed and the predicted numbers of events in that bin divided by the total uncertainty. The total uncertainty is calculated as the sum in quadrature of the statistical uncertainty in the observed measurement and the statistical and systematic uncertainties in the background-only fit to data.

Distributions of $\min[M(\ell,\mathrm{b})]$ in events with 0 t-tagged jets, 1 or more W-tagged jets, and 1 b-tagged jets for the combined electron and muon samples in the signal region. Example signal distributions are also shown. The background distributions correspond to the background-only fit to data, while signal distributions are before the fit to data. The last bin includes overflow events and its content is divided by the bin width. The distributions in each category have variable-size bins, chosen so that the statistical uncertainty in the total background in each bin is less than 30%. The lower panel in each plot shows the difference between the observed and the predicted numbers of events in that bin divided by the total uncertainty. The total uncertainty is calculated as the sum in quadrature of the statistical uncertainty in the observed measurement and the statistical and systematic uncertainties in the background-only fit to data.

Distributions of $\min[M(\ell,\mathrm{b})]$ in events with 0 t-tagged jets, 1 or more W-tagged jets, and 2 or more b-tagged jets for the combined electron and muon samples in the signal region. Example signal distributions are also shown. The background distributions correspond to the background-only fit to data, while signal distributions are before the fit to data. The last bin includes overflow events and its content is divided by the bin width. The distributions in each category have variable-size bins, chosen so that the statistical uncertainty in the total background in each bin is less than 30%. The lower panel in each plot shows the difference between the observed and the predicted numbers of events in that bin divided by the total uncertainty. The total uncertainty is calculated as the sum in quadrature of the statistical uncertainty in the observed measurement and the statistical and systematic uncertainties in the background-only fit to data.

Distributions of $\min[M(\ell,\mathrm{b})]$ in events with 1 or more t-tagged jets, 0 W-tagged jets, and 1 b-tagged jets for the combined electron and muon samples in the signal region. Example signal distributions are also shown. The background distributions correspond to the background-only fit to data, while signal distributions are before the fit to data. The last bin includes overflow events and its content is divided by the bin width. The distributions in each category have variable-size bins, chosen so that the statistical uncertainty in the total background in each bin is less than 30%. The lower panel in each plot shows the difference between the observed and the predicted numbers of events in that bin divided by the total uncertainty. The total uncertainty is calculated as the sum in quadrature of the statistical uncertainty in the observed measurement and the statistical and systematic uncertainties in the background-only fit to data.

Distributions of $\min[M(\ell,\mathrm{b})]$ in events with 1 or more t-tagged jets, 0 W-tagged jets, and 2 or more b-tagged jets for the combined electron and muon samples in the signal region. Example signal distributions are also shown. The background distributions correspond to the background-only fit to data, while signal distributions are before the fit to data. The last bin includes overflow events and its content is divided by the bin width. The distributions in each category have variable-size bins, chosen so that the statistical uncertainty in the total background in each bin is less than 30%. The lower panel in each plot shows the difference between the observed and the predicted numbers of events in that bin divided by the total uncertainty. The total uncertainty is calculated as the sum in quadrature of the statistical uncertainty in the observed measurement and the statistical and systematic uncertainties in the background-only fit to data.

Distributions of $\min[M(\ell,\mathrm{b})]$ in events with 1 or more t-tagged jets, 1 or more W-tagged jets, and 1 b-tagged jets for the combined electron and muon samples in the signal region. Example signal distributions are also shown. The background distributions correspond to the background-only fit to data, while signal distributions are before the fit to data. The last bin includes overflow events and its content is divided by the bin width. The distributions in each category have variable-size bins, chosen so that the statistical uncertainty in the total background in each bin is less than 30%. The lower panel in each plot shows the difference between the observed and the predicted numbers of events in that bin divided by the total uncertainty. The total uncertainty is calculated as the sum in quadrature of the statistical uncertainty in the observed measurement and the statistical and systematic uncertainties in the background-only fit to data.

Distributions of $\min[M(\ell,\mathrm{b})]$ in events with 1 or more t-tagged jets, 1 or more W-tagged jets, and 2 or more b-tagged jets for the combined electron and muon samples in the signal region. Example signal distributions are also shown. The background distributions correspond to the background-only fit to data, while signal distributions are before the fit to data. The last bin includes overflow events and its content is divided by the bin width. The distributions in each category have variable-size bins, chosen so that the statistical uncertainty in the total background in each bin is less than 30%. The lower panel in each plot shows the difference between the observed and the predicted numbers of events in that bin divided by the total uncertainty. The total uncertainty is calculated as the sum in quadrature of the statistical uncertainty in the observed measurement and the statistical and systematic uncertainties in the background-only fit to data.

Expected and observed limits at 95% CL for an LH $X_{5/3}$ after combining all categories for the same-sign dilepton final state. The theoretical uncertainty in the signal cross section is shown as a narrow band around the theoretical prediction.

Expected and observed limits at 95% CL for an RH $X_{5/3}$ after combining all categories for the same-sign dilepton final state. The theoretical uncertainty in the signal cross section is shown as a narrow band around the theoretical prediction.

Expected and observed limits at 95% CL for an LH $X_{5/3}$ after combining all categories for the single-lepton final state. The theoretical uncertainty in the signal cross section is shown as a narrow band around the theoretical prediction.

Expected and observed limits at 95% CL for an RH $X_{5/3}$ after combining all categories for the single-lepton final state. The theoretical uncertainty in the signal cross section is shown as a narrow band around the theoretical prediction.

Expected and observed limits at 95% CL for an LH $X_{5/3}$ after combining the same-sign dilepton and single-lepton final states. The theoretical uncertainty in the signal cross section is shown as a narrow band around the theoretical prediction.

Expected and observed limits at 95% CL for an RH $X_{5/3}$ after combining the same-sign dilepton and single-lepton final states. The theoretical uncertainty in the signal cross section is shown as a narrow band around the theoretical prediction.

Lower mass exclusion limits for scalar and vector LQs.

Observed upper limits on the production cross section for pair production of LQs decaying into a top quark and a muon or bottom quark and a neutrino at 95% CL in the $M_{LQ} - B(LQ \rightarrow t\mu)$ plane.

Lower mass exclusion limits for scalar and vector LQs.

Observed 95% CL upper limits on the product of the production cross section and the branching fraction for a new spin-1 Z prime resonance decaying to ZH, as a function of the resonance mass hypothesis.

Observed 95% CL upper limits on the product of the production cross section and the branching fraction for a new spin-1 V prime (Wprime/Zprime) resonance decaying to VH, as a function of the resonance mass hypothesis.

Data from Fig. 3 lower right panel. The $\textrm{pp} \to \textrm{ZZ}\to \ell\ell\ell^{\prime}\ell^{\prime}$ normalized differential cross section at $\sqrt{s} = 13$ TeV as a function of the jet multiplicity with $|\eta_{j}| < 2.4$.

Data from Fig. 4 upper right panel. The $\textrm{pp} \to \textrm{ZZ}\to \ell\ell\ell^{\prime}\ell^{\prime}$ normalized differential cross section at $\sqrt{s} = 13$ TeV as a function of the p$_T$-leading jet transverse momentum with $|\eta_{j}| < 4.7$.

Data from Fig. 5 left panel. The $\textrm{pp} \to \textrm{ZZ}\to \ell\ell\ell^{\prime}\ell^{\prime}$ normalized differential cross section at $\sqrt{s} = 13$ TeV as a function of the p$_T$-leading jet transverse momentum with $|\eta_{j}| < 4.7$.

Data from Fig. 4 lower right panel. The $\textrm{pp} \to \textrm{ZZ}\to \ell\ell\ell^{\prime}\ell^{\prime}$ normalized differential cross section at $\sqrt{s} = 13$ TeV as a function of the p$_T$-leading jet $|\eta_{j}|$ with $|\eta_{j}| < 4.7$.

Data from Fig. 5 right panel. The $\textrm{pp} \to \textrm{ZZ}\to \ell\ell\ell^{\prime}\ell^{\prime}$ normalized differential cross section at $\sqrt{s} = 13$ TeV as a function of the p$_T$-subleading jet $|\eta_{j}|$ with $|\eta_{j}| < 4.7$.

Data from Fig. 6 left panel. The $\textrm{pp} \to \textrm{ZZ}\to \ell\ell\ell^{\prime}\ell^{\prime}$ normalized differential cross section at $\sqrt{s} = 13$ TeV as a function of the invariant mass of the two p$_T$ leading jets with $|\eta_{j}| < 4.7$.

Data from Fig. 6 right panel. The $\textrm{pp} \to \textrm{ZZ}\to \ell\ell\ell^{\prime}\ell^{\prime}$ normalized differential cross section at $\sqrt{s} = 13$ TeV as a function of the invariant of pseudorapidity separation of the two $p$_T$ leading jets with $|\eta_{j}| < 4.7$.

Data from Fig. 2 upper left. The $\textrm{pp} \to \textrm{ZZ}\to \ell\ell\ell^{\prime}\ell^{\prime}$ differential cross section at $\sqrt{s} = 8$ TeV as a function of the jet multiplicity with $|\eta| < 4.7$.

Data from Fig. 3 upper left panel. The $\textrm{pp} \to \textrm{ZZ}\to \ell\ell\ell^{\prime}\ell^{\prime}$ normalized differential cross section at $\sqrt{s} = 8$ TeV as a function of the jet multiplicity with $|\eta| < 4.7$.

Data from Fig. 2 lower left panel. The $\textrm{pp} \to \textrm{ZZ}\to \ell\ell\ell^{\prime}\ell^{\prime}$ differential cross section at $\sqrt{s} = 8$ TeV as a function of the jet multiplicity with $|\eta| < 2.4$.

Data from Fig. 3 lower left panel. The $\textrm{pp} \to \textrm{ZZ}\to \ell\ell\ell^{\prime}\ell^{\prime}$ normalized differential cross section at $\sqrt{s} = 8$ TeV as a function of the jet multiplicity with $|\eta| < 2.4$.

Data from Fig. 4 upper left panel. The $\textrm{pp} \to \textrm{ZZ}\to \ell\ell\ell^{\prime}\ell^{\prime}$ normalized differential cross section at $\sqrt{s} = 8$ TeV as a function of the p$_T$-leading jet transverse momentum with $|\eta| < 4.7$.

Data from Fig. 4 lower left panel. The $\textrm{pp} \to \textrm{ZZ}\to \ell\ell\ell^{\prime}\ell^{\prime}$ normalized differential cross section at $\sqrt{s} = 8$ TeV as a function of the p$_T$-leading jet transverse $\eta$ with $|\eta| < 4.7$.

Bin-to-bin correlation in the measured cross section for Z+jets as a function of inclusive jet multiplicity, $N_{\text{jets}}^{\text{min}}$.

Measured cross section for Z+jets as a function of the transverse momementum of the Z boson, $p_{\text{T}}(\text{Z})$, for events with at least one jet and breakdown of the relative uncertainty.

Bin-to-bin correlation in the measured cross section for Z+jets as a function of the transverse momementum of the Z boson, $p_{\text{T}}(\text{Z})$, for events with at least one jet.

Measured cross section for Z+jets as a function of the transverse momenum of the first jet, $p_{\text{T}}(\text{j}_1)$, and breakdown of the relative uncertainty.

Bin-to-bin correlation in the measured cross section for Z+jets as a function of the transverse momenum of the first jet, $p_{\text{T}}(\text{j}_1)$.

Measured cross section for Z+jets as a function of the transverse momenum of the second jet, $p_{\text{T}}(\text{j}_2)$, and breakdown of the relative uncertainty.

Bin-to-bin correlation in the measured cross section for Z+jets as a function of the transverse momenum of the second jet, $p_{\text{T}}(\text{j}_2)$.

Measured cross section for Z+jets as a function of the transverse momenum of the third jet, $p_{\text{T}}(\text{j}_3)$, and breakdown of the relative uncertainty.

Bin-to-bin correlation in the measured cross section for Z+jets as a function of the transverse momenum of the third jet, $p_{\text{T}}(\text{j}_3)$.

Measured cross section for Z+jets as a function of the absolute rapidity of the first jet, $|y(\text{j}_1)|$, and breakdown of the relative uncertainty.

Bin-to-bin correlation in the measured cross section for Z+jets as a function of the absolute rapidity of the first jet, $|y(\text{j}_1)|$.

Measured cross section for Z+jets as a function of the absolute rapidity of the second jet, $|y(\text{j}_2)|$, and breakdown of the relative uncertainty.

Bin-to-bin correlation in the measured cross section for Z+jets as a function of the absolute rapidity of the second jet, $|y(\text{j}_2)|$.

Measured cross section for Z+jets as a function of the absolute rapidity of the third jet, $|y(\text{j}_3)|$, and breakdown of the relative uncertainty.

Bin-to-bin correlation in the measured cross section for Z+jets as a function of the absolute rapidity of the third jet, $|y(\text{j}_3)|$.

Measured cross section for Z+jets as a function of the $H_{\text{T}}$ observable for event with at least one jet and breakdown of the relative uncertainty.

Bin-to-bin correlation in the measured cross section for Z+jets as a function of the $H_{\text{T}}$ observable for event with at least one jet.

Measured cross section for Z+jets as a function of the $H_{\text{T}}$ observable for event with at least two jets and breakdown of the relative uncertainty.

Bin-to-bin correlation in the measured cross section for Z+jets as a function of the $H_{\text{T}}$ observable for event with at least two jets.

Measured cross section for Z+jets as a function of the $H_{\text{T}}$ observable for event with at least three jets and breakdown of the relative uncertainty.

Bin-to-bin correlation in the measured cross section for Z+jets as a function of the $H_{\text{T}}$ observable for event with at least three jets.

Measured cross section for Z+jets as a function of the transverse momentum balance between the Z boson and the accompanying jets, $p_{\text{T}}^{\text{bal}}$, for events with at least one jet and breakdown of the relative uncertainty.

Bin-to-bin correlation in the measured cross section for Z+jets as a function of the transverse momentum balance between the Z boson and the accompanying jets, $p_{\text{T}}^{\text{bal}}$, for events with at least one jet.

Measured cross section for Z+jets as a function of the transverse momentum balance between the Z boson and the accompanying jets, $p_{\text{T}}^{\text{bal}}$, for events with at least two jets and breakdown of the relative uncertainty.

Bin-to-bin correlation in the measured cross section for Z+jets as a function of the transverse momentum balance between the Z boson and the accompanying jets, $p_{\text{T}}^{\text{bal}}$, for events with at least two jets.

Measured cross section for Z+jets as a function of the transverse momentum balance between the Z boson and the accompanying jets, $p_{\text{T}}^{\text{bal}}$, for events with at least three jets and breakdown of the relative uncertainty.

Bin-to-bin correlation in the measured cross section for Z+jets as a function of the transverse momentum balance between the Z boson and the accompanying jets, $p_{\text{T}}^{\text{bal}}$, for events with at least three jets.

Measured cross section for Z+jets as a function of the JZB variable with no restriction on $p_{\text{T}}(\text{Z})$ and breakdown of the relative uncertainty.

Bin-to-bin correlation in the measured cross section for Z+jets as a function of the JZB variable with no restriction on $p_{\text{T}}(\text{Z})$.

Measured cross section for Z+jets as a function of the JZB variable for $p_{\text{T}}(\text{Z})\leq50\,$GeV) and breakdown of the relative uncertainty.

Bin-to-bin correlation in the measured cross section for Z+jets as a function of the JZB variable for $p_{\text{T}}(\text{Z})\leq50\,$GeV).

Measured cross section for Z+jets as a function of the JZB variable for $p_{\text{T}}(\text{Z})>50\,$GeV and breakdown of the relative uncertainty.

Bin-to-bin correlation in the measured cross section for Z+jets as a function of the JZB variable for $p_{\text{T}}(\text{Z})>50\,$GeV.