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.

Differential cross section of J/psi mesons from b hadrons as a function of dimuon pT in pp and PbPb collisions. The PbPb cross sections are normalised by TAA for direct comparison. Global uncertainties arise from the integrated luminosity uncertainty in pp collisions, and the number of minimum bias events and TAA uncertainties for PbPb collisions.

Differential cross section of prompt J/psi mesons as a function of dimuon rapidty in pp and PbPb collisions. The PbPb cross sections are normalised by TAA for direct comparison. Global uncertainties arise from the integrated luminosity uncertainty in pp collisions, and the number of minimum bias events and TAA uncertainties for PbPb collisions.

Differential cross section of J/psi mesons from b hadrons as a function of dimuon rapidty in pp and PbPb collisions. The PbPb cross sections are normalised by TAA for direct comparison. Global uncertainties arise from the integrated luminosity uncertainty in pp collisions, and the number of minimum bias events and TAA uncertainties for PbPb collisions.

Nuclear modification factor of prompt J/psi mesons as a function of dimuon rapidity. Global uncertainties arise from the integrated luminosity uncertainty of the pp data, and the number of minimum bias events and TAA uncertainties of the PbPb data.

Nuclear modification factor of prompt J/psi mesons as a function of collision centrality. Global uncertainties arise from the full pp statistical and systematic uncertainties (including integrated luminosity uncertainty), and the number of minimum bias events uncertainty of the PbPb data.

Nuclear modification factor of prompt J/psi mesons as a function of dimuon pT. Global uncertainties arise from the integrated luminosity uncertainty of the pp data, and the number of minimum bias events and TAA uncertainties of the PbPb data.

Nuclear modification factor of prompt J/psi mesons as a function of dimuon pT in the mid-rapidity interval. Global uncertainties arise from the integrated luminosity uncertainty of the pp data, and the number of minimum bias events and TAA uncertainties of the PbPb data.

Nuclear modification factor of prompt J/psi mesons as a function of dimuon pT in the most forward rapidity interval. Global uncertainties arise from the integrated luminosity uncertainty of the pp data, and the number of minimum bias events and TAA uncertainties of the PbPb data.

Nuclear modification factor of prompt J/psi mesons as a function of collision centrality in the mid- and most forward rapidity intervals. Global uncertainties arise from the full pp statistical and systematic uncertainties (including integrated luminosity uncertainty), and the number of minimum bias events uncertainty of the PbPb data.

Nuclear modification factor of prompt J/psi mesons as a function of pT in three centrality bins. Global uncertainties arise from the integrated luminosity uncertainty of the pp data, and the number of minimum bias events and TAA uncertainties of the PbPb data.

Nuclear modification factor of prompt J/psi mesons as a function of collision centrality in the most forward rapidity range at moderate and high pT. Global uncertainties arise from the full pp statistical and systematic uncertainties (including integrated luminosity uncertainty), and the number of minimum bias events uncertainty of the PbPb data.

Nuclear modification factor of prompt J/psi mesons as a function of collision centrality in the mid-rapidity range. Global uncertainties arise from the full pp statistical and systematic uncertainties (including integrated luminosity uncertainty), and the number of minimum bias events uncertainty of the PbPb data.

Nuclear modification factor of prompt psi(2S) mesons as a function of collision centrality in the mid-rapidity range. Global uncertainties arise from the full pp statistical and systematic uncertainties (including integrated luminosity uncertainty), and the number of minimum bias events uncertainty of the PbPb data.

95% CL intervals on the nuclear modification factor of prompt psi(2S) mesons as a function of collision centrality in the mid-rapidity range. Global uncertainties arise from the full pp statistical and systematic uncertainties (including integrated luminosity uncertainty), and the number of minimum bias events uncertainty of the PbPb data.

Nuclear modification factor of prompt J/psi and psi(2S) mesons as a function of dimuon pT in the mid-rapidity range. Global uncertainties arise from the integrated luminosity uncertainty of the pp data, and the number of minimum bias events and TAA uncertainties of the PbPb data.

Nuclear modification factor of prompt J/psi mesons as a function of collision centrality in the forward rapidity range. Global uncertainties arise from the full pp statistical and systematic uncertainties (including integrated luminosity uncertainty), and the number of minimum bias events uncertainty of the PbPb data.

Nuclear modification factor of prompt psi(2S) mesons as a function of collision centrality in the forward rapidity range. Global uncertainties arise from the full pp statistical and systematic uncertainties (including integrated luminosity uncertainty), and the number of minimum bias events uncertainty of the PbPb data.

95% CL intervals on the nuclear modification factor of prompt psi(2S) mesons as a function of collision centrality in the forward rapidity range. Global uncertainties arise from the full pp statistical and systematic uncertainties (including integrated luminosity uncertainty), and the number of minimum bias events uncertainty of the PbPb data.

Nuclear modification factor of prompt J/psi mesons as a function of dimuon pT in the forward rapidity range. Global uncertainties arise from the integrated luminosity uncertainty of the pp data, and the number of minimum bias events and TAA uncertainties of the PbPb data.

Nuclear modification factor of prompt psi(2S) mesons as a function of dimuon pT in the forward rapidity range. Global uncertainties arise from the integrated luminosity uncertainty of the pp data, and the number of minimum bias events and TAA uncertainties of the PbPb data.

95% CL intervals on the nuclear modification factor of prompt psi(2S) mesons as a function of dimuon pT in the forward rapidity range. Global uncertainties arise from the integrated luminosity uncertainty of the pp data, and the number of minimum bias events and TAA uncertainties of the PbPb data.

Nuclear modification factor of nonprompt J/psi mesons as a function of dimuon rapidity. Global uncertainties arise from the integrated luminosity uncertainty of the pp data, and the number of minimum bias events and TAA uncertainties of the PbPb data.

Nuclear modification factor of nonprompt J/psi mesons as a function of collision centrality. Global uncertainties arise from the full pp statistical and systematic uncertainties (including integrated luminosity uncertainty), and the number of minimum bias events uncertainty of the PbPb data.

Nuclear modification factor of nonprompt J/psi mesons as a function of dimuon pT. Global uncertainties arise from the integrated luminosity uncertainty of the pp data, and the number of minimum bias events and TAA uncertainties of the PbPb data.

Nuclear modification factor of nonprompt J/psi mesons as a function of dimuon pT in the mid-rapidity interval. Global uncertainties arise from the integrated luminosity uncertainty of the pp data, and the number of minimum bias events and TAA uncertainties of the PbPb data.

Nuclear modification factor of nonprompt J/psi mesons as a function of dimuon pT in the most forward rapidity interval. Global uncertainties arise from the integrated luminosity uncertainty of the pp data, and the number of minimum bias events and TAA uncertainties of the PbPb data.

Nuclear modification factor of nonprompt J/psi mesons as a function of collision centrality in the mid- and most forward rapidity intervals. Global uncertainties arise from the full pp statistical and systematic uncertainties (including integrated luminosity uncertainty), and the number of minimum bias events uncertainty of the PbPb data.

Nuclear modification factor of nonprompt J/psi mesons as a function of pT in three centrality bins. Global uncertainties arise from the integrated luminosity uncertainty of the pp data, and the number of minimum bias events and TAA uncertainties of the PbPb data.

Nuclear modification factor of nonprompt J/psi mesons as a function of collision centrality in the most forward rapidity range at moderate and high pT. Global uncertainties arise from the full pp statistical and systematic uncertainties (including integrated luminosity uncertainty), and the number of minimum bias events uncertainty of the PbPb data.

95% CL intervals on the double ratio of measured yields, $(N_{\psi\mathrm{(2S)}} / N_{J/\psi})_{\mathrm{PbPb}} / (N_{\psi\mathrm{(2S)}} / N_{J/\psi})_{pp}$, as a function of pT, for the forward rapidity analysis bin.

Double ratio of measured yields, $(N_{\psi\mathrm{(2S)}} / N_{J/\psi})_{\mathrm{PbPb}} / (N_{\psi\mathrm{(2S)}} / N_{J/\psi})_{pp}$, as a function of centrality, for the midrapidity analysis bin.

95% CL intervals on the double ratio of measured yields, $(N_{\psi\mathrm{(2S)}} / N_{J/\psi})_{\mathrm{PbPb}} / (N_{\psi\mathrm{(2S)}} / N_{J/\psi})_{pp}$, as a function of centrality, for the midrapidity analysis bin.

Double ratio of measured yields, $(N_{\psi\mathrm{(2S)}} / N_{J/\psi})_{\mathrm{PbPb}} / (N_{\psi\mathrm{(2S)}} / N_{J/\psi})_{pp}$, as a function of centrality, for the forward rapidity analysis bin.

95% CL intervals on the double ratio of measured yields, $(N_{\psi\mathrm{(2S)}} / N_{J/\psi})_{\mathrm{PbPb}} / (N_{\psi\mathrm{(2S)}} / N_{J/\psi})_{pp}$, as a function of centrality, for the forward rapidity analysis bin.

Double ratio of measured yields, $(N_{\psi\mathrm{(2S)}} / N_{J/\psi})_{\mathrm{PbPb}} / (N_{\psi\mathrm{(2S)}} / N_{J/\psi})_{pp}$, integrated over centrality, for the midrapidity and forward rapidity analysis bins.

95% CL interval on the double ratio of measured yields, $(N_{\psi\mathrm{(2S)}} / N_{J/\psi})_{\mathrm{PbPb}} / (N_{\psi\mathrm{(2S)}} / N_{J/\psi})_{pp}$, integrated over centrality, for the midrapidity and forward rapidity analysis bins.

Mean PT and RMS for J/PSI production from b-hadron decay (assuming unpolarised).

DSIG/DY for prompt (unpolarised) J/PSI and from b-hadron decay integrated over PT.

Double differential cross section in PT and YRAP for prompt J/PSI production assuming no polarisation. The first systematic error is the uncorrelated systematic error and the second is the correlated systematic error.

Double differential cross section in PT and YRAP for J/PSI from b-hadron decays assuming no polarisation. The first systematic error is the uncorrelated systematic error and the second is the correlated systematic error.

Double differential cross section in PT and YRAP for prompt J/PSI production assuming fully transversely polarised J/PSIs The first systematic error is the uncorrelated systematic error and the second is the correlated systematic error.

Double differential cross section in PT and YRAP for prompt J/PSI production assuming fully longitudinally polarised J/PSIs The first systematic error is the uncorrelated systematic error and the second is the correlated systematic error.

Fraction of J/PSIs from b-hadron decay in bins of PT and YRAP assuming no polaraisation. The first systematic error is uncorrelated between bins and the second is the uncertainty due to the unknown polarisation of the prompt J/PSI.

Differential J/$\psi$ production cross section at $\sqrt{s}=2.76$ TeV. J/$\psi$ production cross section at $\sqrt{s}=2.76$ TeV. The first uncertainty is statistical, the second one is $y$-correlated, the third one is uncorrelated. Polarization-related uncertainties are not included. The value for $-0.9<y<0.9$ refers to J/$\psi\rightarrow {\rm e}^+{\rm e}^-$, the remaining values to J/$\psi\rightarrow \mu^+\mu^-$.

psi(2S) double-differential cross section times branching fraction and the corresponding scaling factors to obtain the cross sections for different polarization scenarios (azimuthal polarization parameter in the center of mass helicity frame lambda_theta^HX = +1, -1, +0.03) as a function of pT for |y| < 1.2.

Ratio of psi(2S) to J/psi differential cross section times branching fractions assuming unpolarized productions as a function of pT for |y| < 1.2.

Nuclear modification factor RAA of psi(2S) integrated over centrality, for the midrapidity and forward rapidity analysis bins.

Lambda-Phi in the CS frame for the J/psi as a function of pT for 0.6 < |y| < 1.2.

Lambda-Theta-Phi in the CS frame for the J/psi as a function of pT for 0.0 < |y| < 0.6.

Lambda-Theta-Phi in the CS frame for the J/psi as a function of pT for 0.6 < |y| < 1.2.

Lambda-Tilde in the CS frame for the J/psi as a function of pT for 0.0 < |y| < 0.6.

Lambda-Tilde in the CS frame for the J/psi as a function of pT for 0.6 < |y| < 1.2.

Lambda-Theta in the HX frame for the J/psi as a function of pT for 0.0 < |y| < 0.6.

Lambda-Theta in the HX frame for the J/psi as a function of pT for 0.6 < |y| < 1.2.

Lambda-Phi in the HX frame for the J/psi as a function of pT for 0.0 < |y| < 0.6.

Lambda-Phi in the HX frame for the J/psi as a function of pT for 0.6 < |y| < 1.2.

Lambda-Theta-Phi in the HX frame for the J/psi as a function of pT for 0.0 < |y| < 0.6.

Lambda-Theta-Phi in the HX frame for the J/psi as a function of pT for 0.6 < |y| < 1.2.

Lambda-Tilde in the HX frame for the J/psi as a function of pT for 0.0 < |y| < 0.6.

Lambda-Tilde in the HX frame for the J/psi as a function of pT for 0.6 < |y| < 1.2.

Lambda-Theta in the PX frame for the J/psi as a function of pT for 0.0 < |y| < 0.6.

Lambda-Theta in the PX frame for the J/psi as a function of pT for 0.6 < |y| < 1.2.

Lambda-Phi in the PX frame for the J/psi as a function of pT for 0.0 < |y| < 0.6.

Lambda-Phi in the PX frame for the J/psi as a function of pT for 0.6 < |y| < 1.2.

Lambda-Theta-Phi in the PX frame for the J/psi as a function of pT for 0.0 < |y| < 0.6.

Lambda-Theta-Phi in the PX frame for the J/psi as a function of pT for 0.6 < |y| < 1.2.

Lambda-Tilde in the PX frame for the J/psi as a function of pT for 0.0 < |y| < 0.6.

Lambda-Tilde in the PX frame for the J/psi as a function of pT for 0.6 < |y| < 1.2.

Lambda-Theta in the CS frame for the psi(2S) as a function of pT for 0.0 < |y| < 0.6.

Lambda-Theta in the CS frame for the psi(2S) as a function of pT for 0.6 < |y| < 1.2.

Lambda-Theta in the CS frame for the psi(2S) as a function of pT for 1.2 < |y| < 1.5.

Lambda-Phi in the CS frame for the psi(2S) as a function of pT for 0.0 < |y| < 0.6.

Lambda-Phi in the CS frame for the psi(2S) as a function of pT for 0.6 < |y| < 1.2.

Lambda-Phi in the CS frame for the psi(2S) as a function of pT for 1.2 < |y| < 1.5.

Lambda-Theta-Phi in the CS frame for the psi(2S) as a function of pT for 0.0 < |y| < 0.6.

Lambda-Theta-Phi in the CS frame for the psi(2S) as a function of pT for 0.6 < |y| < 1.2.

Lambda-Theta-Phi in the CS frame for the psi(2S) as a function of pT for 1.2 < |y| < 1.5.

Lambda-Tilde in the CS frame for the psi(2S) as a function of pT for 0.0 < |y| < 0.6.

Lambda-Tilde in the CS frame for the psi(2S) as a function of pT for 0.6 < |y| < 1.2.

Lambda-Tilde in the CS frame for the psi(2S) as a function of pT for 1.2 < |y| < 1.5.

Lambda-Theta in the HX frame for the psi(2S) as a function of pT for 0.0 < |y| < 0.6.

Lambda-Theta in the HX frame for the psi(2S) as a function of pT for 0.6 < |y| < 1.2.

Lambda-Theta in the HX frame for the psi(2S) as a function of pT for 1.2 < |y| < 1.5.

Lambda-Phi in the HX frame for the psi(2S) as a function of pT for 0.0 < |y| < 0.6.

Lambda-Phi in the HX frame for the psi(2S) as a function of pT for 0.6 < |y| < 1.2.

Lambda-Phi in the HX frame for the psi(2S) as a function of pT for 1.2 < |y| < 1.5.

Lambda-Theta-Phi in the HX frame for the psi(2S) as a function of pT for 0.0 < |y| < 0.6.

Lambda-Theta-Phi in the HX frame for the psi(2S) as a function of pT for 0.6 < |y| < 1.2.

Lambda-Theta-Phi in the HX frame for the psi(2S) as a function of pT for 1.2 < |y| < 1.5.

Lambda-Tilde in the HX frame for the psi(2S) as a function of pT for 0.0 < |y| < 0.6.

Lambda-Tilde in the HX frame for the psi(2S) as a function of pT for 0.6 < |y| < 1.2.

Lambda-Tilde in the HX frame for the psi(2S) as a function of pT for 1.2 < |y| < 1.5.

Lambda-Theta in the PX frame for the psi(2S) as a function of pT for 0.0 < |y| < 0.6.

Lambda-Theta in the PX frame for the psi(2S) as a function of pT for 0.6 < |y| < 1.2.

Lambda-Theta in the PX frame for the psi(2S) as a function of pT for 1.2 < |y| < 1.5.

Lambda-Phi in the PX frame for the psi(2S) as a function of pT for 0.0 < |y| < 0.6.

Lambda-Phi in the PX frame for the psi(2S) as a function of pT for 0.6 < |y| < 1.2.

Lambda-Phi in the PX frame for the psi(2S) as a function of pT for 1.2 < |y| < 1.5.

Lambda-Theta-Phi in the PX frame for the psi(2S) as a function of pT for 0.0 < |y| < 0.6.

Lambda-Theta-Phi in the PX frame for the psi(2S) as a function of pT for 0.6 < |y| < 1.2.

Lambda-Theta-Phi in the PX frame for the psi(2S) as a function of pT for 1.2 < |y| < 1.5.

Lambda-Tilde in the PX frame for the psi(2S) as a function of pT for 0.0 < |y| < 0.6.

Lambda-Tilde in the PX frame for the psi(2S) as a function of pT for 0.6 < |y| < 1.2.

Lambda-Tilde in the PX frame for the psi(2S) as a function of pT for 1.2 < |y| < 1.5.

The double-differential cross sections for prompt J/PSI production (assuming no polarisation) and production of J/PSI from b-hadron decays as a function of transverse momentum for the rapidity range 2.5-3.0. Also shown in the final column is the fraction (in %) of J/PSIs from the latter.

The double-differential cross sections for prompt J/PSI production (assuming no polarisation) and production of J/PSI from b-hadron decays as a function of transverse momentum for the rapidity range 3.0-3.5. Also shown in the final column is the fraction (in %) of J/PSIs from the latter.

The double-differential cross sections for prompt J/PSI production (assuming no polarisation) and production of J/PSI from b-hadron decays as a function of transverse momentum for the rapidity range 3.5-4.0. Also shown in the final column is the fraction (in %) of J/PSIs from the latter.

The double-differential cross sections for prompt J/PSI production (assuming no polarisation) and production of J/PSI from b-hadron decays as a function of transverse momentum for the rapidity range 4.0-4.5. Also shown in the final column is the fraction (in %) of J/PSIs from the latter.

The total integrated cross sections for UPSILON 1S 2S and 3S times their dimuon branching ratio in the rapidity range 2.0-4.5 and transverse momentum 0-15 GeV/c.

Differential production cross section for UPSILON 1S 2s and 3S mesons times their dimuon branching fraction as a function of transverse momentum ;.

Differential production cross section for UPSILON 1S 2s and 3S mesons times their dimuon branching fraction as a function of rapidity ;.

The double-differential cross section (times the dimuon decay probability) for the production of UPSILON 1S, 2S and 3S mesons as a function of their transverse momentum for the rapidity range 2.0-2.5, assuming no polarisation.

The double-differential cross section (times the dimuon decay probability) for the production of UPSILON 1S, 2S and 3S mesons as a function of their transverse momentum for the rapidity range 2.5-3.0, assuming no polarisation.

The double-differential cross section (times the dimuon decay probability) for the production of UPSILON 1S, 2S and 3S mesons as a function of their transverse momentum for the rapidity range 3.0-3.5, assuming no polarisation.

The double-differential cross section (times the dimuon decay probability) for the production of UPSILON 1S, 2S and 3S mesons as a function of their transverse momentum for the rapidity range 3.5-4.0, assuming no polarisation.

The double-differential cross section (times the dimuon decay probability) for the production of UPSILON 1S, 2S and 3S mesons as a function of their transverse momentum for the rapidity range 4.0-4.5, assuming no polarisation.

Ratios of the cross-sections for UPSILON(2S)->MU+MU- and UPSILON(3S)->MU+MU- with respect toe UPSILON(1S)->MU+MU- as a function of PT in the rapidity range 2.0-4.0 ;.

Ratios of the cross-sections for UPSILON(2S)->MU+MU- and UPSILON(3S)->MU+MU- with respect to UPSILON(1S)->MU+MU- as a function of rapidity in the rapidity range 2.0-4.0 ;.