Differential cross section times dimuon branching fraction of Y(1S) as a function of $y^{Y}_{CM}$ in pPb collisions. The global uncertainty arises from the integrated luminosity uncertainty in pPb collisions.

Differential cross section times dimuon branching fraction of Y(2S) as a function of $y^{Y}_{CM}$ in pPb collisions. The global uncertainty arises from the integrated luminosity uncertainty in pPb collisions.

Differential cross section times dimuon branching fraction of Y(3S) as a function of $y^{Y}_{CM}$ in pPb collisions. The global uncertainty arises from the integrated luminosity uncertainty in pPb collisions.

Differential cross section times dimuon branching fraction of Y(1S) as a function of pT in pp collisions. The global uncertainty arises from the integrated luminosity uncertainty in pp collisions.

Differential cross section times dimuon branching fraction of Y(2S) as a function of pT in pp collisions. The global uncertainty arises from the integrated luminosity uncertainty in pp collisions.

Differential cross section times dimuon branching fraction of Y(3S) as a function of pT in pp collisions. The global uncertainty arises from the integrated luminosity uncertainty in pp collisions.

Differential cross section times dimuon branching fraction of Y(1S) as a function of $|y^{Y}_{CM}|$ in pp collisions. The global uncertainty arises from the integrated luminosity uncertainty in pp collisions.

Differential cross section times dimuon branching fraction of Y(2S) as a function of $|y^{Y}_{CM}|$ in pp collisions. The global uncertainty arises from the integrated luminosity uncertainty in pp collisions.

Differential cross section times dimuon branching fraction of Y(3S) as a function of $|y^{Y}_{CM}|$ in pp collisions. The global uncertainty arises from the integrated luminosity uncertainty in pp collisions.

Nuclear modification factor of Y(1S) as a function of pT. The global uncertainty arises from the integrated luminosity uncertainties in pPb and pp collisions.

Nuclear modification factor of Y(2S) as a function of pT. The global uncertainty arises from the integrated luminosity uncertainties in pPb and pp collisions.

Nuclear modification factor of Y(3S) as a function of pT. The global uncertainty arises from the integrated luminosity uncertainties in pPb and pp collisions.

Nuclear modification factor of Y(1S) as a function of $y^{Y}_{CM}$. The global uncertainty arises from the integrated luminosity uncertainties in pPb and pp collisions.

Nuclear modification factor of Y(2S) as a function of $y^{Y}_{CM}$. The global uncertainty arises from the integrated luminosity uncertainties in pPb and pp collisions.

Nuclear modification factor of Y(3S) as a function of $y^{Y}_{CM}$. The global uncertainty arises from the integrated luminosity uncertainties in pPb and pp collisions.

Nuclear modification factor of Y(1S) at forward and backward $y^{Y}_{CM}$ for pT < 6 GeV/c. The global uncertainty arises from the integrated luminosity uncertainties in pPb and pp collisions.

Nuclear modification factor of Y(2S) at forward and backward $y^{Y}_{CM}$ for pT < 6 GeV/c. The global uncertainty arises from the integrated luminosity uncertainties in pPb and pp collisions.

Nuclear modification factor of Y(3S) at forward and backward $y^{Y}_{CM}$ for pT < 6 GeV/c. The global uncertainty arises from the integrated luminosity uncertainties in pPb and pp collisions.

Nuclear modification factor of Y(1S) at forward and backward $y^{Y}_{CM}$ for 6 < pT < 30 GeV/c. The global uncertainty arises from the integrated luminosity uncertainties in pPb and pp collisions.

Nuclear modification factor of Y(2S) at forward and backward $y^{Y}_{CM}$ for 6 < pT < 30 GeV/c. The global uncertainty arises from the integrated luminosity uncertainties in pPb and pp collisions.

Nuclear modification factor of Y(3S) at forward and backward $y^{Y}_{CM}$ for 6 < pT < 30 GeV/c. The global uncertainty arises from the integrated luminosity uncertainties in pPb and pp collisions.

RFB of Y(1S) versus $N^{|\eta_{lab}|<2.4}_{tracks}$.

RFB of Y(2S) versus $N^{|\eta_{lab}|<2.4}_{tracks}$.

RFB of Y(3S) versus $N^{|\eta_{lab}|<2.4}_{tracks}$.

RFB of Y(1S) versus $E^{|\eta_{lab}|>4}_{T}$.

RFB of Y(2S) versus $E^{|\eta_{lab}|>4}_{T}$.

RFB of Y(3S) versus $E^{|\eta_{lab}|>4}_{T}$.

Nuclear modification factor of Y(1S), Y(2S), and Y(3S) integrated over pT and $y^{Y}_{CM}$. The global uncertainty arises from the integrated luminosity uncertainties in pPb and pp collisions.

$D^+ D^0$~mass distributions for selected candidates with the $D^0$ background subtracted. Uncertainties on the data points are statistical only and represent one standard deviation, calculated as a sum in quadrature of the assigned weights from the background-subtraction procedure.

Mass distributions for selected $D^+ D^+$ candidates with the $D^0$ background subtracted. Uncertainties on the data points are statistical only and represent one standard deviation, calculated as a sum in quadrature of the assigned weights from the background-subtraction procedure.

Mass distributions for selected $D^+ D^0 \pi^+$ candidates with the $D^0$ background subtracted. Uncertainties on the data points are statistical only and represent one standard deviation, calculated as a sum in quadrature of the assigned weights from the background-subtraction procedure.

Background-subtracted distributions for the multiplicity of tracks reconstructed in the vertex detector for $T_{cc}^+\to D^0 D^0 \pi^+$ signal, low-mass $D^0\bar{D}^0$ and $D^0D^0$ pairs. The binning scheme is chosen to have an approximately uniform distribution for $D^0\bar{D}^0$ pairs. The distributions for the $D^0\bar{D}^0$ and $D^0D^0$ pairs are normalised to the same yields as the $T_{cc}^+\to D^0 D^0 \pi^+$ signal. For better visualisation, the points are slightly displaced from the bin centres. Uncertainties on the data points are statistical only and represent one standard deviation, calculated as a sum in quadrature of the assigned weights from the background-subtraction procedure.

Background-subtracted transverse momentum spectra for $T_{cc}^+\to D^0 D^0 \pi^+$ signal, low-mass $D^0\bar{D}^0$ and $D^0D^0$ pairs. The binning scheme is chosen to have an approximately uniform distribution for $D^0\bar{D}^0$ pairs. The distributions for the $D^0\bar{D}^0$ and $D^0D^0$ pairs are normalised to the same yields as the $T_{cc}^+\to D^0 D^0 \pi^+$ signal. For better visualisation, the points are slightly displaced from the bin centres. Uncertainties on the data points are statistical only and represent one standard deviation, calculated as a sum in quadrature of the assigned weights from the background-subtraction procedure.

Mass distributions for selected $D^0 \bar{D}^0$ candidates with the $D^0$ background subtracted. Uncertainties on the data points are statistical only and represent one standard deviation, calculated as a sum in quadrature of the assigned weights from the background-subtraction procedure.

Mass distributions for selected $D^0 D^0$ candidates with the $D^0$ background subtracted. Uncertainties on the data points are statistical only and represent one standard deviation, calculated as a sum in quadrature of the assigned weights from the background-subtraction procedure.

Mass distributions for selected $\bar{D}^0D^0\pi^+$ candidates with the $D^0$ background subtracted. Uncertainties on the data points are statistical only and represent one standard deviation, calculated as a sum in quadrature of the assigned weights from the background-subtraction procedure.

Mass distributions for selected $D^0D^0\pi^+$ candidates with the $D^0$ background subtracted. Uncertainties on the data points are statistical only and represent one standard deviation, calculated as a sum in quadrature of the assigned weights from the background-subtraction procedure.

Mass distributions for selected $D^0D^-$ candidates with the $D^0$ background subtracted. Uncertainties on the data points are statistical only and represent one standard deviation, calculated as a sum in quadrature of the assigned weights from the background-subtraction procedure.

Mass distributions for selected $D^0D^+$ candidates with the $D^0$ background subtracted. Uncertainties on the data points are statistical only and represent one standard deviation, calculated as a sum in quadrature of the assigned weights from the background-subtraction procedure.

Distribution of $D^0 D^0 \pi^+$ mass where the contribution of the non-$D^0$ background has been statistically subtracted by assigning the a weight to every candidate.

Mass distribution for $D^0 \pi^+$ pairs from selected $D^0 D^0 \pi^+$ candidates with a mass below the $D^{*+}D^0$ mass threshold with non-$D^0$ background subtracted by assigning the a weight to every candidate.

$D^0 D^0$~mass distributions for selected candidates with the $D^0$ background subtracted by assigning the a weight to every candidate.

$D^+ D^0$~mass distributions for selected candidates with the $D^0$ background subtracted by assigning the a weight to every candidate.

Mass distributions for selected $D^+ D^+$ candidates with the $D^0$ background subtracted by assigning the a weight to every candidate.

Mass distributions for selected $D^+ D^0 \pi^+$ candidates with the $D^0$ background subtracted by assigning the a weight to every candidate.

Nuclear modification factor for prompt charged particles at 5TeV with respect to pseudorapidity and transverse momentum in the forward region region. The pseudorapidity is expressed in the nucleon-nucleon center-of-mass system. First uncertainty is statistical, the second is systematic and the third is from the luminosity. Luminosity uncertainty is fully correlated among the different kinematic ranges.

Nuclear modification factor for prompt charged particles at 5TeV with respect to pseudorapidity and transverse momentum in the backward region region. The pseudorapidity is expressed in the nucleon-nucleon center-of-mass system. First uncertainty is statistical, the second is systematic and the third is from the luminosity. Luminosity uncertainty is fully correlated among the different kinematic ranges.

Distribution of $m(l,\gamma)$ in the $N_{jet}\geq 3$ signal region.

Distribution of $\Delta R(l,\gamma)$ in the $N_{jet}\geq 3$ signal region.

Distribution of $\Delta R(j,\gamma)$ in the $N_{jet}\geq 3$ signal region.

Fit result of the multijet template obtained with loosely isolated leptons and the electroweak background to the measured $m_{T}(W)$ distribution with isolated leptons in the $N_{jet}=2$, $N_{b jet}=0$ selection for electrons.

Fit result of the multijet template obtained with loosely isolated leptons and the electroweak background to the measured $m_{T}(W)$ distribution with isolated leptons in the $N_{jet}=2$, $N_{b jet}=0$ selection for muons.

Distribution of the invariant mass of the lepton and the photon ($m(l,\gamma)$) in the $N_{jet}\geq 3$, $N_{b jet}=0$ selection for the e channel.

Distribution of the invariant mass of the lepton and the photon ($m(l,\gamma)$) in the $N_{jet}\geq 3$, $N_{b jet}=0$ selection for the $\mu$ channel.

Extracted scale factors for the contribution from misidentified electrons for the three data-taking periods, and the Z$\gamma$, W$\gamma$ simulations.

Predicted and observed yields in the control regions in the $N_{jet}= 3$ and $\geq 4$ seletions using the post-fit values of the nuisance parameters.

Predicted and observed yields in the signal regions in the $N_{jet}= 3$ and $\geq 4$ seletions using the post-fit values of the nuisance parameters.

The measured inclusive ttgamma cross section in the fiducial phase space compared to the prediction from simulation using Madgraph_aMC@NLO at a center-of-mass energy of 13 TeV.

Summary of the measured cross section ratios with respect to the NLO cross section prediction for signal regions binned in the electron channel, muon channel and the combined single lepton measurement.

The unfolded differential cross sections for $p_{T}(\gamma)$ and the comparison to simulations.

The unfolded differential cross sections for $|\eta(\gamma)|$ and the comparison to simulations.

The unfolded differential cross sections for $\Delta R(l,\gamma)$ and the comparison to simulations.

The covariance matrix of systematic uncertainties for the unfolded differential measurement for $p_{T}(\gamma)$.

The covariance matrix of systematic uncertainties for the unfolded differential measurement for $|\eta(\gamma)|$.

The covariance matrix of systematic uncertainties for the unfolded differential measurement for $\Delta R(l,\gamma)$.

The covariance matrix of statistic uncertainties for the unfolded differential measurement for $p_{T}(\gamma)$.

The covariance matrix of statistic uncertainties for the unfolded differential measurement for $|\eta(\gamma)|$.

The covariance matrix of statistic uncertainties for the unfolded differential measurement for $\Delta R(l,\gamma)$.

The correlation matrix of statistical uncertainties for the unfolded differential measurement for $p_{T}(\gamma)$.

The correlation matrix of statistical uncertainties for the unfolded differential measurement for $|\eta(\gamma)|$.

The correlation matrix of statistical uncertainties for the unfolded differential measurement for $\Delta R(l,\gamma)$.

The correlation matrix of systematic uncertainties for the unfolded differential measurement for $p_{T}(\gamma)$.

The correlation matrix of systematic uncertainties for the unfolded differential measurement for $|\eta(\gamma)|$.

The correlation matrix of systematic uncertainties for the unfolded differential measurement for $\Delta R(l,\gamma)$.

Summary of the one-dimensional intervals at 68 and 95% CL.

The observed and predicted post-fit yields for the combined Run 2 data set in the SR3 signal region for the electron channel.

The observed and predicted post-fit yields for the combined Run 2 data set in the SR3 signal region for the muon channel.

The observed and predicted post-fit yields for the combined Run 2 data set in the SR4p signal region for the electron channel.

The observed and predicted post-fit yields for the combined Run 2 data set in the SR4p signal region for the muon channel.

Negative log-likelihood ratio values with respect to the best fit value of the one-dimensional profiled scan for the Wilson coefficient $c_{tZ}$.

Negative log-likelihood ratio values with respect to the best fit value of the one-dimensional profiled scan for the Wilson coefficient $c^{I}_{tZ}$.

Negative log-likelihood ratio values with respect to the best fit value of the one-dimensional scan for the Wilson coefficient $c_{tZ}$.

Negative log-likelihood ratio values with respect to the best fit value of the one-dimensional scan for the Wilson coefficient $c^{I}_{tZ}$.

Negative log-likelihood ratio values with respect to the best fit value of the two-dimensional scan for the Wilson coefficients $c_{tZ}$ and $c^{I}_{tZ}$.

Differential cross sections as function of transverse momentum imbalance between Z boson and dijet for Z+ ≥ 2 jets events.

Differential cross sections as function of transverse momentum imbalance between leading and subleading jet for Z+ ≥ 2 jets events.

Area normalized distribution as function of Delta Phi between Z boson and the leading jet for Z+ ≥ 1 jet events.

Area normalized distribution as function of transverse momentum imbalance between Z boson and the leading jet for Z+ ≥ 1 jet events.

Area normalized distribution as function of Delta Phi between Z boson and dijet for Z+ ≥ 2 jets events.

Area normalized distribution as function of transverse momentum imbalance between Z boson and dijet for Z+ ≥ 2 jets events.

Area normalized distribution as function of transverse momentum imbalance between leading and subleading jet for Z+ ≥ 2 jets events

Response matrix for Delta Phi between Z boson and the leading jet for Z+ ≥ 1 jet events.

Response matrix for transverse momentum imbalance between Z boson and the leading jet for Z+ ≥ 1 jet events.

Response matrix for Delta Phi between Z boson and dijet for Z+ ≥ 2 jets events.

Response matrix for transverse momentum imbalance between Z boson and dijet for Z+ ≥ 2 jets events.

Response matrix for transverse momentum imbalance between the leading jet and subleading jet for Z+ ≥ 2 jets events.

Correlation matrix for Delta Phi between Z boson and the leading jet for Z+ ≥ 1 jet events (for absolute differential cross section measurements).

Correlation matrix for transverse momentum imbalance between Z boson and the leading jet for Z+ ≥ 1 jet events (for absolute differential cross section measurements).

Correlation matrix for Delta Phi between Z boson and dijet for Z+ ≥ 2 jets events (for absolute differential cross section measurements).

Correlation matrix for transverse momentum imbalance between Z boson and dijet for Z+ ≥ 2 jets events (for absolute differential cross section measurements).

Correlation matrix for transverse momentum imbalance between the leading jet and subleading jet for Z+ ≥ 2 jets events (for absolute differential cross section measurements).

Correlation matrix for Delta Phi between Z boson and the leading jet for Z+ ≥ 1 jet events (for normalized differential cross section measurements).

Correlation matrix for transverse momentum imbalance between Z boson and the leading jet for Z+ ≥ 1 jet events (for normalized differential cross section measurements).

Correlation matrix for Delta Phi between Z boson and dijet for Z+ ≥ 2 jets events (for normalized differential cross section measurements).

Correlation matrix for transverse momentum imbalance between Z boson and dijet for Z+ ≥ 2 jets events (for normalized differential cross section measurements).

Correlation matrix for transverse momentum imbalance between the leading jet and subleading jet for Z+ ≥ 2 jets events (for normalized differential cross section measurements).

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].

Numerical results of $b \bar{b}$- and $c \bar{c}$-dijet cross-sections, $c \bar{c}$/$b \bar{b}$ dijet cross-section ratios and their total uncertainties as a function of the leading jet $p_T$.

Numerical results of $b \bar{b}$- and $c \bar{c}$-dijet cross-sections, $c\bar{c}$/$b \bar{b}$ dijet cross-section ratios and their total uncertainties as a function of $m_{jj}$ (dijet invariant mass).

Covariance matrix, corresponding to the total uncertainties, obtained between the leading jet eta intervals of the $b \bar{b}$-dijet differential cross sections. The unit of all the elements of the matrix is nb$^2$.

Covariance matrix, corresponding to the total uncertainties, obtained between the leading jet eta intervals of the $c \bar{c}$-dijet differential cross sections. The unit of all the elements of the matrix is nb$^2$.

Covariance matrix, corresponding to the total uncertainties, obtained between the leading jet $\eta$ intervals of the $b \bar{b}$ (horizontal) and $c \bar{c}$ (vertical) differential cross sections. The unit of all the elements of the matrix is nb$^2$.

Covariance matrix, corresponding to the total uncertainties, obtained between the $\Delta y^*$ intervals of the $b \bar{b}$-dijet differential cross sections. The unit of all the elements of the matrix is nb$^2$.

Covariance matrix, corresponding to the total uncertainties, obtained between the $\Delta y^*$ intervals of the $c \bar{c}$-dijet differential cross sections. The unit of all the elements of the matrix is nb$^2$.

Covariance matrix, corresponding to the total uncertainties, obtained between the $\Delta y^*$ intervals of the $b \bar{b}$ (horizontal) and $c \bar{c}$ (vertical) differential cross sections. The unit of all the elements of the matrix is nb$^2$.

Covariance matrix, corresponding to the total uncertainties, obtained between the leading jet $p_T$ intervals of the $b \bar{b}$-dijet differential cross sections. The unit of all the elements of the matrix is (nb GeV$/c)^2$ and the $p_T$ intervals are given in GeV$/c$.

Covariance matrix, corresponding to the total uncertainties, obtained between the leading jet $p_T$ intervals of the $c \bar{c}$-dijet differential cross sections. The unit of all the elements of the matrix is (nb GeV$/c)^2$ and the $p_T$ intervals are given in GeV$/$c .

Covariance matrix, corresponding to the total uncertainties, obtained between the leading jet $p_T$ intervals of the $b \bar{b}$ (horizontal) and $c \bar{c}$ (vertical) differential cross sections. The unit of all the elements of the matrix is (nb GeV$/c)^2$ and the $p_T$ intervals are given in GeV$/c$.

Covariance matrix, corresponding to the total uncertainties, obtained between the $m_{jj}$ intervals of the $b \bar{b}$-dijet differential cross sections. The unit of all the elements of the matrix is (nb GeV$/c^2)^2$ and the mass intervals are given in GeV$/c^2$.

Covariance matrix, corresponding to the total uncertainties, obtained between the $m_{jj}$ intervals of the $c \bar{c}$-dijet differential cross sections. The unit of all the elements of the matrix is (nb GeV$/c^2)^2$ and the mass intervals are given in GeV$/c^2$.

Covariance matrix, corresponding to the total uncertainties, obtained between the $m_{jj}$ intervals of the $b \bar{b}$ (horizontal) and $c \bar{c}$ (vertical) differential cross sections. The unit of all the elements of the matrix is $($nb GeV$/c^2)^2$ and the mass intervals are given in GeV$/c^2$.