Summary of results for cross-section of non-prompt $\psi(2S)$ decaying to a muon pair for 13 TeV data in fb/GeV. Uncertainties are statistical and systematic, respectively.

Summary of results for non-prompt fraction of $J/\psi$ decaying to a muon pair as a percentage for 13 TeV data. Uncertainties are statistical and systematic, respectively.

Summary of results for non-prompt fraction of $\psi(2S)$ decaying to a muon pair as a percentage for 13 TeV data. Uncertainties are statistical and systematic, respectively.

Summary of results for prompt $\psi(2S)$ over $J/\psi$ ratio as a percentage for 13 TeV data. Uncertainties are statistical and systematic, respectively.

Summary of results for non-prompt $\psi(2S)$ over $J/\psi$ ratio as a percentage for 13 TeV data. Uncertainties are statistical and systematic, respectively.

Correction factors for an assumed angular dependence $\propto 1+\lambda_{\theta} \cos^2\theta $ in the helicity frame, for several values of the parameter $\lambda_{\theta}$. The correction factors were found to be the same (within 1% - 2%) for $J\psi$ and $\psi(2S)$ mesons, for prompt and non-prompt production mechanisms, and for the three rapidity intervals considered in this paper.

The measured cross section as a function of exclusive jet multiplicity, $N_{\text{jets}}$, when $50<p_T<100$ GeV

The measured cross section as a function of exclusive jet multiplicity, $N_{\text{jets}}$, when $p_T>100$ GeV

The measured cross section as a function of $\Delta\phi_{Z,jet1}$, when $p_T<10$ GeV

The measured cross section as a function of $\Delta\phi_{Z,jet1}$, when $10<p_T<30$ GeV

The measured cross section as a function of $\Delta\phi_{Z,jet1}$, when $30<p_T<50$ GeV

The measured cross section as a function of $\Delta\phi_{Z,jet1}$, when $50<p_T<100$ GeV

The measured cross section as a function of $\Delta\phi_{Z,jet1}$, when $p_T>100$ GeV

The measured cross section as a function of $\Delta\phi_{jet1,jet2}$, when $p_T<10$ GeV

The measured cross section as a function of $\Delta\phi_{jet1,jet2}$, when $10<p_T<30$ GeV

The measured cross section as a function of $\Delta\phi_{jet1,jet2}$, when $30<p_T<50$ GeV

The measured cross section as a function of $\Delta\phi_{jet1,jet2}$, when $50<p_T<100$ GeV

The measured cross section as a function of $\Delta\phi_{jet1,jet2}$, when $p_T>100$ GeV

Ratio of $p_{\mathrm{T,ch\ jet}}$-differential cross section of charm jets tagged with $\mathrm{D^{0}}$ mesons for different $R$: $\sigma(R=0.2)/\sigma(R=0.4)$ and $\sigma(R=0.2)/\sigma(R=0.6)$ in pp collisions at $\sqrt{s}=13$ TeV.

Ratio of $p_{\mathrm{T,ch\ jet}}$-differential cross section of charm jets tagged with $\mathrm{D^{0}}$ mesons for different $R$: $\sigma(R=0.2)/\sigma(R=0.4)$ and $\sigma(R=0.2)/\sigma(R=0.6)$ in pp collisions at $\sqrt{s}=5.02$ TeV.

The fraction of $\mathrm{D^{0}}$ jets over inclusive charged-particle jets in pp collisions at $\sqrt{s}=5.02$ TeV for $R=0.2$, $0.4$, and $0.6$.

Distributions of $z^{\mathrm{ch}}_{||}$-differential yield of charm jets tagged with $\mathrm{D^{0}}$ mesons normalised by the number of $\mathrm{D^{0}}$ jets within each distribution in pp collisions at $\sqrt{s}=13$ TeV in four jet-$p_{\mathrm{T}}$ intervals $5<p_{\mathrm{T,ch\ jet}}<7$ GeV/$c$, $7<p_{\mathrm{T,ch\ jet}}<10$ GeV/$c$, $10<p_{\mathrm{T,ch\ jet}}<15$ GeV/$c$, and $15<p_{\mathrm{T,ch\ jet}}<50$ GeV/$c$ for jet radii $R=0.2$, $0.4$, and $0.6$.

Distributions of $z^{\mathrm{ch}}_{||}$-differential yield of charm jets tagged with $\mathrm{D^{0}}$ mesons normalised by the number of $\mathrm{D^{0}}$ jets within each distribution in pp collisions at $\sqrt{s}=5.02$ TeV in four jet-$p_{\mathrm{T}}$ intervals $5<p_{\mathrm{T,ch\ jet}}<7$ GeV/$c$, $7<p_{\mathrm{T,ch\ jet}}<10$ GeV/$c$, $10<p_{\mathrm{T,ch\ jet}}<15$ GeV/$c$, and $15<p_{\mathrm{T,ch\ jet}}<50$ GeV/$c$ for jet radii $R=0.2$, $0.4$, and $0.6$.

$p_{\mathrm{T,ch\ jet}}$-differential cross section of charm jets tagged with $\mathrm{D^{0}}$ mesons for $R=0.3$ in pp collisions at $\sqrt{s}=5.02$ TeV.

The fraction of $\mathrm{D^{0}}$ jets over inclusive charged-particle jets in pp collisions at $\sqrt{s}=5.02$ TeV for $R=0.3$.

Distributions of $z^{\mathrm{ch}}_{||}$-differential yield of charm jets tagged with $\mathrm{D^{0}}$ mesons normalised by the number of $\mathrm{D^{0}}$ jets within each distribution in pp collisions at $\sqrt{s}=5.02$ TeV in four jet-$p_{\mathrm{T}}$ intervals $5<p_{\mathrm{T,ch\ jet}}<7$ GeV/$c$, $7<p_{\mathrm{T,ch\ jet}}<10$ GeV/$c$, $10<p_{\mathrm{T,ch\ jet}}<15$ GeV/$c$, and $15<p_{\mathrm{T,ch\ jet}}<50$ GeV/$c$ for jet radius $R=0.3$.

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.

The measured differential cross section in the $76 \le M_{ll} < 106$ GeV mass window, in bins of the dilepton transverse momentum. The values are normalized by the bin width. This entry contains the covariance matrix of the results.

The measured differential cross section in the $106 \le M_{ll} < 170$ GeV mass window, in bins of the dilepton transverse momentum. The values are normalized by the bin width.

The measured differential cross section in the $106 \le M_{ll} < 170$ GeV mass window, in bins of the dilepton transverse momentum. The values are normalized by the bin width. This entry contains the covariance matrix of the results.

The measured differential cross section in the $170 \le M_{ll} < 350$ GeV mass window, in bins of the dilepton transverse momentum. The values are normalized by the bin width.

The measured differential cross section in the $170 \le M_{ll} < 350$ GeV mass window, in bins of the dilepton transverse momentum. The values are normalized by the bin width. This entry contains the covariance matrix of the results.

The measured differential cross section in the $350 \le M_{ll} < 1000$ GeV mass window, in bins of the dilepton transverse momentum. The values are normalized by the bin width.

The measured differential cross section in the $350 \le M_{ll} < 1000$ GeV mass window, in bins of the dilepton transverse momentum. The values are normalized by the bin width. This entry contains the covariance matrix of the results.

The measured differential cross section in the $50 \le M_{ll} < 76$ GeV mass window, in bins of the dilepton transverse momentum. At least one jet is required. The values are normalized by the bin width.

The measured differential cross section in the $50 \le M_{ll} < 76$ GeV mass window, in bins of the dilepton transverse momentum. At least one jet is required. The values are normalized by the bin width. This entry contains the covariance matrix of the results.

The measured differential cross section in the $76 \le M_{ll} < 106$ GeV mass window, in bins of the dilepton transverse momentum. At least one jet is required. The values are normalized by the bin width.

The measured differential cross section in the $76 \le M_{ll} < 106$ GeV mass window, in bins of the dilepton transverse momentum. At least one jet is required. The values are normalized by the bin width. This entry contains the covariance matrix of the results.

The measured differential cross section in the $106 \le M_{ll} < 170$ GeV mass window, in bins of the dilepton transverse momentum. At least one jet is required. The values are normalized by the bin width.

The measured differential cross section in the $106 \le M_{ll} < 170$ GeV mass window, in bins of the dilepton transverse momentum. At least one jet is required. The values are normalized by the bin width. This entry contains the covariance matrix of the results.

The measured differential cross section in the $170 \le M_{ll} < 350$ GeV mass window, in bins of the dilepton transverse momentum. At least one jet is required. The values are normalized by the bin width.

The measured differential cross section in the $170 \le M_{ll} < 350$ GeV mass window, in bins of the dilepton transverse momentum. At least one jet is required. The values are normalized by the bin width. This entry contains the covariance matrix of the results.

The measured differential cross section in the $50 \le M_{ll} < 76$ GeV mass window, in bins of the $\varphi^*$ variable. The values are normalized by the bin width.

The measured differential cross section in the $50 \le M_{ll} < 76$ GeV mass window, in bins of the $\varphi^*$ variable. The values are normalized by the bin width. This entry contains the covariance matrix of the results.

The measured differential cross section in the $76 \le M_{ll} < 106$ GeV mass window, in bins of the $\varphi^*$ variable. The values are normalized by the bin width.

The measured differential cross section in the $76 \le M_{ll} < 106$ GeV mass window, in bins of the $\varphi^*$ variable. The values are normalized by the bin width. This entry contains the covariance matrix of the results.

The measured differential cross section in the $106 \le M_{ll} < 170$ GeV mass window, in bins of the $\varphi^*$ variable. The values are normalized by the bin width.

The measured differential cross section in the $106 \le M_{ll} < 170$ GeV mass window, in bins of the $\varphi^*$ variable. The values are normalized by the bin width. This entry contains the covariance matrix of the results.

The measured differential cross section in the $170 \le M_{ll} < 350$ GeV mass window, in bins of the $\varphi^*$ variable. The values are normalized by the bin width.

The measured differential cross section in the $170 \le M_{ll} < 350$ GeV mass window, in bins of the $\varphi^*$ variable. The values are normalized by the bin width. This entry contains the covariance matrix of the results.

The measured differential cross section in the $350 \le M_{ll} < 1000$ GeV mass window, in bins of the $\varphi^*$ variable. The values are normalized by the bin width.

The measured differential cross section in the $350 \le M_{ll} < 1000$ GeV mass window, in bins of the $\varphi^*$ variable. The values are normalized by the bin width. This entry contains the covariance matrix of the results.

The measured differential cross section in the $50 \le M_{ll} < 76$ GeV mass window, in bins of the dilepton transverse momentum, divided by the measured differential cross section in the $76 \le M_{ll} < 106$ GeV mass window. The values are not normalized by the bin width.

The measured differential cross section in the $50 \le M_{ll} < 76$ GeV mass window, in bins of the dilepton transverse momentum, divided by the measured differential cross section in the $76 \le M_{ll} < 106$ GeV mass window. The values are not normalized by the bin width. This entry contains the covariance matrix of the results.

The measured differential cross section in the $106 \le M_{ll} < 170$ GeV mass window, in bins of the dilepton transverse momentum, divided by the measured differential cross section in the $76 \le M_{ll} < 106$ GeV mass window. The values are not normalized by the bin width.

The measured differential cross section in the $106 \le M_{ll} < 170$ GeV mass window, in bins of the dilepton transverse momentum, divided by the measured differential cross section in the $76 \le M_{ll} < 106$ GeV mass window. The values are not normalized by the bin width. This entry contains the covariance matrix of the results.

The measured differential cross section in the $170 \le M_{ll} < 350$ GeV mass window, in bins of the dilepton transverse momentum, divided by the measured differential cross section in the $76 \le M_{ll} < 106$ GeV mass window. The values are not normalized by the bin width.

The measured differential cross section in the $170 \le M_{ll} < 350$ GeV mass window, in bins of the dilepton transverse momentum, divided by the measured differential cross section in the $76 \le M_{ll} < 106$ GeV mass window. The values are not normalized by the bin width. This entry contains the covariance matrix of the results.

The measured differential cross section in the $350 \le M_{ll} < 1000$ GeV mass window, in bins of the dilepton transverse momentum, divided by the measured differential cross section in the $76 \le M_{ll} < 106$ GeV mass window. The values are not normalized by the bin width.

The measured differential cross section in the $350 \le M_{ll} < 1000$ GeV mass window, in bins of the dilepton transverse momentum, divided by the measured differential cross section in the $76 \le M_{ll} < 106$ GeV mass window. The values are not normalized by the bin width. This entry contains the covariance matrix of the results.

The measured differential cross section in the $50 \le M_{ll} < 76$ GeV mass window, in bins of the dilepton transverse momentum, divided by the measured differential cross section in the $76 \le M_{ll} < 106$ GeV mass window. At least one jet is required. The values are not normalized by the bin width.

The measured differential cross section in the $50 \le M_{ll} < 76$ GeV mass window, in bins of the dilepton transverse momentum, divided by the measured differential cross section in the $76 \le M_{ll} < 106$ GeV mass window. At least one jet is required. The values are not normalized by the bin width. This entry contains the covariance matrix of the results.

The measured differential cross section in the $106 \le M_{ll} < 170$ GeV mass window, in bins of the dilepton transverse momentum, divided by the measured differential cross section in the $76 \le M_{ll} < 106$ GeV mass window. At least one jet is required. The values are not normalized by the bin width.

The measured differential cross section in the $106 \le M_{ll} < 170$ GeV mass window, in bins of the dilepton transverse momentum, divided by the measured differential cross section in the $76 \le M_{ll} < 106$ GeV mass window. At least one jet is required. The values are not normalized by the bin width. This entry contains the covariance matrix of the results.

The measured differential cross section in the $170 \le M_{ll} < 350$ GeV mass window, in bins of the dilepton transverse momentum, divided by the measured differential cross section in the $76 \le M_{ll} < 106$ GeV mass window. At least one jet is required. The values are not normalized by the bin width.

The measured differential cross section in the $170 \le M_{ll} < 350$ GeV mass window, in bins of the dilepton transverse momentum, divided by the measured differential cross section in the $76 \le M_{ll} < 106$ GeV mass window. At least one jet is required. The values are not normalized by the bin width. This entry contains the covariance matrix of the results.

The measured differential cross section in the $50 \le M_{ll} < 76$ GeV mass window, in bins of the $\varphi^*$ variable, divided by the measured differential cross section in the $76 \le M_{ll} < 106$ GeV mass window. At least one jet is required. The values are not normalized by the bin width.

The measured differential cross section in the $50 \le M_{ll} < 76$ GeV mass window, in bins of the $\varphi^*$ variable, divided by the measured differential cross section in the $76 \le M_{ll} < 106$ GeV mass window. At least one jet is required. The values are not normalized by the bin width. This entry contains the covariance matrix of the results.

The measured differential cross section in the $106 \le M_{ll} < 170$ GeV mass window, in bins of the $\varphi^*$ variable, divided by the measured differential cross section in the $76 \le M_{ll} < 106$ GeV mass window. At least one jet is required. The values are not normalized by the bin width.

The measured differential cross section in the $106 \le M_{ll} < 170$ GeV mass window, in bins of the $\varphi^*$ variable, divided by the measured differential cross section in the $76 \le M_{ll} < 106$ GeV mass window. At least one jet is required. The values are not normalized by the bin width. This entry contains the covariance matrix of the results.

The measured differential cross section in the $170 \le M_{ll} < 350$ GeV mass window, in bins of the $\varphi^*$ variable, divided by the measured differential cross section in the $76 \le M_{ll} < 106$ GeV mass window. At least one jet is required. The values are not normalized by the bin width.

The measured differential cross section in the $170 \le M_{ll} < 350$ GeV mass window, in bins of the $\varphi^*$ variable, divided by the measured differential cross section in the $76 \le M_{ll} < 106$ GeV mass window. At least one jet is required. The values are not normalized by the bin width. This entry contains the covariance matrix of the results.

The measured differential cross section in the $350 \le M_{ll} < 1000$ GeV mass window, in bins of the $\varphi^*$ variable, divided by the measured differential cross section in the $76 \le M_{ll} < 106$ GeV mass window. At least one jet is required. The values are not normalized by the bin width.

The measured differential cross section in the $350 \le M_{ll} < 1000$ GeV mass window, in bins of the $\varphi^*$ variable, divided by the measured differential cross section in the $76 \le M_{ll} < 106$ GeV mass window. At least one jet is required. The values are not normalized by the bin width. This entry contains the covariance matrix of the results.

Response matrix for pT ll mass 50-76 (electron channel)

Response matrix for pT ll mass 50-76 (muon channel)

Response matrix for pT ll mass 76-106 (electron channel)

Response matrix for pT ll mass 76-106 (muon channel)

Response matrix for pT ll mass 106-170 (electron channel)

Response matrix for pT ll mass 106-170 (muon channel)

Response matrix for pT ll mass 170-350 (electron channel)

Response matrix for pT ll mass 170-350 (muon channel)

Response matrix for pT ll mass 350-1000 (electron channel)

Response matrix for pT ll mass 50-76 jet (electron channel)

Response matrix for pT ll mass 50-76 jet (muon channel)

Response matrix for pT ll mass 76-106 jet (electron channel)

Response matrix for pT ll mass 76-106 jet (muon channel)

Response matrix for pT ll mass 106-170 jet (electron channel)

Response matrix for pT ll mass 106-170 jet (muon channel)

Response matrix for pT ll mass 170-350 jet (electron channel)

Response matrix for pT ll mass 170-350 jet (muon channel)

Response matrix for phistar mass 50-76 (electron channel)

Response matrix for phistar mass 50-76 (muon channel)

Response matrix for phistar mass 76-106 (electron channel)

Response matrix for phistar mass 76-106 (muon channel)

Response matrix for phistar mass 106-170 (electron channel)

Response matrix for phistar mass 106-170 (muon channel)

Response matrix for phistar mass 170-350 (electron channel)

Response matrix for phistar mass 170-350 (muon channel)

Response matrix for phistar mass 350-1000 (electron channel)

Response matrix for phistar mass 350-1000 (muon channel)

Single-differential production cross-sections for nonprompt $J/\psi$ mesons as a function of $p_\text{T}$. The first uncertainties are statistical, the second are correlated systematic uncertainties shared between intervals, and the last are uncorrelated systematic uncertainties.

Single-differential production cross-sections for prompt $J/\psi$ mesons as a function of $y$. The first uncertainties are statistical, the second are correlated systematic uncertainties shared between intervals, and the last are uncorrelated systematic uncertainties.

Single-differential production cross-sections for nonprompt $J/\psi$ mesons as a function of $y$. The first uncertainties are statistical, the second are correlated systematic uncertainties shared between intervals, and the last are uncorrelated systematic uncertainties.

Fraction of nonprompt $J/\psi$ mesons (in \%) in ($p_\text{T},y$) intervals. The first uncertainty is statistical and the second is systematic.

Nuclear modification factor $R_{p\text{Pb}}$ for prompt $J/\psi$ mesons as a function of $y$. The first uncertainty is statistical and the second is systematic.

Nuclear modification factor $R_{p\text{Pb}}$ for nonprompt $J/\psi$ mesons as a function of $y$. The first uncertainty is statistical and the second is systematic.

Cross-section ratios between 8 TeV and 5 TeV measurements for prompt $J/\psi$ mesons as a function of $p_\text{T}$. The first uncertainty is statistical and the second is systematic.

Cross-section ratios between 8 TeV and 5 TeV measurements for prompt $J/\psi$ mesons as a function of $y$. The first uncertainty is statistical and the second is systematic.

Cross-section ratios between 13 TeV and 5 TeV measurements for prompt $J/\psi$ mesons as a function of $p_\text{T}$. The first uncertainty is statistical and the second is systematic.

Cross-section ratios between 13 TeV and 5 TeV measurements for prompt $J/\psi$ mesons as a function of $y$. The first uncertainty is statistical and the second is systematic.

Cross-section ratios between 8 TeV and 5 TeV measurements for nonprompt $J/\psi$ mesons as a function of $p_\text{T}$. The first uncertainty is statistical and the second is systematic.

Cross-section ratios between 8 TeV and 5 TeV measurements for nonprompt $J/\psi$ mesons as a function of $y$. The first uncertainty is statistical and the second is systematic.

Cross-section ratios between 13 TeV and 5 TeV measurements for nonprompt $J/\psi$ mesons as a function of $p_\text{T}$. The first uncertainty is statistical and the second is systematic.

Cross-section ratios between 13 TeV and 5 TeV measurements for nonprompt $J/\psi$ mesons as a function of $y$. The first uncertainty is statistical and the second is systematic.

Relative changes of cross-sections (in \%), for a polarisation of $\lambda_{\theta}=-0.2$ rather than zero, in ($p_\text{T},y$) intervals.

Relative changes of cross-sections (in \%), for a polarisation of $\lambda_{\theta}=-1$ rather than zero, in ($p_\text{T},y$) intervals.

Relative changes of cross-sections (in \%), for a polarisation of $\lambda_{\theta}=+1$ rather than zero, in ($p_\text{T},y$) intervals.

Ratios of transverse-momentum distributions of charged particles in three intervals of multiplicities to the respective one for inclusive (INEL>0) collisions. The spectra were normalised by the integral prior to division.