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.

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

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.

Inclusive J/psi production cross-section as a function of J/psi pT in the J/psi rapidity (|y|) bin 1.5<|y|<2. The first uncertainty is statistical, the second is systematic and the third encapsulates any possible variation due to spin-alignment from the unpolarised central value.

Inclusive J/psi production cross-section as a function of J/psi pT in the J/psi rapidity (|y|) bin 0.75<|y|<1.5. The first uncertainty is statistical, the second is systematic and the third encapsulates any possible variation due to spin-alignment from the unpolarised central value.

Inclusive J/psi production cross-section as a function of J/psi pT in the J/psi rapidity (|y|) bin |y|<0.75. The first uncertainty is statistical, the second is systematic and the third encapsulates any possible variation due to spin-alignment from the unpolarised central value.

Non-prompt to inclusive production cross-section fraction fB as a function of J/psi pT for J/psi rapidity |y|<0.75 under the assumption that prompt and non-prompt J/psi production is unpolarised (lambda_theta = 0). The spin-alignment envelope spans the range of possible prompt cross-sections under various polarisation hypotheses, plus the range of non-prompt cross-sections within lambda_theta = +/- 0.1. The first uncertainty is statistical, the second uncertainty is systematic, the third number is the uncertainty due to spin-alignment.}.

Non-prompt to inclusive production cross-section fraction fB as a function of J/psi pT for J/psi rapidity 0.75<|y|<1.5 under the assumption that prompt and non-prompt J/psi production is unpolarised (lambda_theta = 0). The spin-alignment envelope spans the range of possible prompt cross-sections under various polarisation hypotheses, plus the range of non-prompt cross-sections within lambda_theta = +/-0.1. The first uncertainty is statistical, the second uncertainty is systematic, the third number is the uncertainty due to spin-alignment.

Non-prompt to inclusive production cross-section fraction fB as a function of J/psi pT for J/psi rapidity 1.5<|y|<2 under the assumption that prompt and non-prompt J/psi production is unpolarised (lambda_theta = 0). The spin-alignment envelope spans the range of possible prompt cross-sections under various polarisation hypotheses, plus the range of non-prompt cross-sections within lambda_theta = +/-0.1. The first uncertainty is statistical, the second uncertainty is systematic, the third number is the uncertainty due to spin-alignment.

Non-prompt to inclusive production cross-section fraction fB as a function of J/psi pT for J/psi rapidity 2<|y|<2.4 under the assumption that prompt and non-prompt J/psi production is unpolarised (lambda_theta = 0). The spin-alignment envelope spans the range of possible prompt cross-sections under various polarisation hypotheses, plus the range of non-prompt cross-sections within lambda_theta =+/-0.1. The first uncertainty is statistical, the second uncertainty is systematic, the third number is the uncertainty due to spin-alignment.

Non-prompt J/psi production cross-sections as a function of J/psi pT for J/psi rapidity |y|<0.75 under the assumption that prompt and non-prompt J/psi production is unpolarised (lambda_theta = 0), and the spin-alignment envelope spans the range of non-prompt cross-sections within lambda_theta = +/- 0.1. The first uncertainty is statistical, the second uncertainty is systematic. Comparison is made to FONLL predictions.

Non-prompt J/psi production cross-sections as a function of J/psi pT for J/psi rapidity 0.75<|y|<1.5 under the assumption that prompt and non-prompt J/psi production is unpolarised (lambda_theta = 0), and the spin-alignment envelope spans the range of non-prompt cross-sections within lambda_theta = +/- 0.1. The first uncertainty is statistical, the second uncertainty is systematic. Comparison is made to FONLL predictions.

Non-prompt J/psi production cross-sections as a function of J/psi pT for J/psi rapidity 1.5<|y|<2 under the assumption that prompt and non-prompt J/psi production is unpolarised (lambda_theta = 0), and the spin-alignment envelope spans the range of non-prompt cross-sections within lambda_theta = +/- 0.1. The first uncertainty is statistical, the second uncertainty is systematic. Comparison is made to FONLL predictions.

Non-prompt J/psi production cross-sections as a function of J/psi pT for J/psi rapidity 2<|y|<2.4 under the assumption that prompt and non-prompt J/psi production is unpolarised (lambda_theta = 0), and the spin-alignment envelope spans the range of non-prompt cross-sections within lambda_theta = +/- 0.1. The first uncertainty is statistical, the second uncertainty is systematic. Comparison is made to FONLL predictions.

Prompt J/psi production cross-sections as a function of J/psi pT for J/psi rapidity |y|<0.75. The central value assumes unpolarised (lambda_theta = 0) prompt and non-prompt production, and the spin-alignment envelope spans the range of possible prompt cross-sections under various polarisation hypotheses. The first quoted uncertainty is statistical, the second uncertainty is systematic. Comparison is made to the Colour Evaporation Model prediction.

Prompt J/psi production cross-sections as a function of J/psi pT for J/psi rapidity 0.75<|y|<1.5. The central value assumes unpolarised (lambda_theta = 0) prompt and non-prompt production, and the spin-alignment envelope spans the range of possible prompt cross-sections under various polarisation hypotheses. The first quoted uncertainty is statistical, the second uncertainty is systematic. Comparison is made to the Colour Evaporation Model prediction.

Prompt J/psi production cross-sections as a function of J/psi pT for J/psi rapidity 1.5<|y|<2. The central value assumes unpolarised (lambda_theta = 0) prompt and non-prompt production, and the spin-alignment envelope spans the range of possible prompt cross-sections under various polarisation hypotheses. The first quoted uncertainty is statistical, the second uncertainty is systematic. Comparison is made to the Colour Evaporation Model prediction.

Prompt J/psi production cross-sections as a function of J/psi pT for J/psi rapidity 2<|y|<2.4. The central value assumes unpolarised (lambda_theta = 0) prompt and non-prompt production, and the spin-alignment envelope spans the range of possible prompt cross-sections under various polarisation hypotheses. The first quoted uncertainty is statistical, the second uncertainty is systematic. Comparison is made to the Colour Evaporation Model prediction.

Unweighted J/psi candidate yields in bins of $J/psi transverse momentum and rapidity. Uncertainties are statistical only.

Summary table of all sources of considered systematic uncertainty and statistical uncertainty (as a percentage) on the corrected inclusive J/psi production cross-section, for absolute J/psi rapidities within 2<|y|<2.4. The sources of systematic error shown are, in order, Acceptance, Muon recognition, Trigger, Fitting and the Total systematics. Also shown in the last error is the possible envelope of variation the central result due to uncertainty on spin-alignment of the J\psi.

Summary table of all sources of considered systematic uncertainty and statistical uncertainty (as a percentage) on the corrected inclusive J/psi production cross-section, for absolute J/psi rapidities within 1.5<|y|<2. The sources of systematic error shown are, in order, Acceptance, Muon recognition, Trigger, Fitting and the Total systematics. Also shown in the last error is the possible envelope of variation the central result due to uncertainty on spin-alignment of the J/psi.

Summary table of all sources of considered systematic uncertainty and statistical uncertainty (as a percentage) on the corrected inclusive J/psi production cross-section, for absolute J/psi rapidities within 0.75<|y|<1.5. The sources of systematic error shown are, in order, Acceptance, Muon recognition, Trigger, Fitting and the Total systematics. Also shown in the last error is the possible envelope of variation the central result due to uncertainty on spin-alignment of the J/psi.

Summary table of all sources of considered systematic uncertainty and statistical uncertainty (as a percentage) on the corrected inclusive J/psi production cross-section, for absolute rapidities within |y|<0.75. The sources of systematic error shown are, in order, Acceptance, Muon recognition, Trigger, Fitting and the Total systematics. Also shown in the last error is the possible envelope of variation the central result due to uncertainty on spin-alignment of the J/psi.

Breakdown of sources of systematic uncertainty on the non-prompt J/psi fraction measurements, for the bin |y|<0.75, as a function of the J/psi pT.

Breakdown of sources of systematic uncertainty on the non-prompt J/psi fraction measurements, for the bin 0.75<|y|<1.5, as a function of the J/psi pT.

Breakdown of sources of systematic uncertainty on the non-prompt J/psi fraction measurements, for the bin 1.75<|y|<2.0, as a function of the J/psi pT.

Breakdown of sources of systematic uncertainty on the non-prompt J/psi fraction measurements, for the bin 2.0<|y|<2.4, as a function of the J/psi pT.

Nuclear modification factor Q_pA as function of centrality at forward (2.03<y_cms<3.53) and backward (-4.46<y_cms<-2.96) centre-of-mass rapidity. The first uncertainty is statistical, the second and third ones are the systematic uncertainties. The third uncertainty is fully correlated over centrality.

Nuclear modification factor Q_pA as function of centrality in the mid-rapidity range (-1.37 <y_cms< 0.43). The first uncertainty is statistical, the second and third one are the systematic uncertainties. The third uncertainty is fully correlated over centrality.

Values of Delta<pT^2> as a function of centrality at forward (2.03<y_cms<3.53) and backward (-4.46<y_cms<-2.96) centre-of-mass rapidity. The first uncertainty is statistical, the second and third ones are the systematic uncertainties. The third is related to the uncertainty on the <pT^2> in pp collisions and is correlated over centrality.

Nuclear modification factor Q_pA as function of pT for six centrality classes at backward (-4.46<y_cms<-2.96) centre-of-mass rapidity. The first uncertainty is statistical, the second one and the third ones are the systematic uncertainties. The third uncertainty is fully correlated over pT.

Nuclear modification factor Q_pA as a function of pT for six centrality classes at forward (2.03<y_cms<3.53) centre-of-mass rapidity. The first uncertainty is statistical, the second and the third ones are the systematic uncertainties. The third uncertainty is fully correlated over pT.

Differential production cross sections of psi(2S) as a function of pT.

Differential production cross sections of psi(2S) as a function of rapidity.

integrated production cross section of psi(2S).

psi(2S) over J/psi ratio as a function of pt.

psi(2S) over J/psi ratio as a function of rapidity.

integrated psi(2S) over J/psi production cross section ratio.

Differential production cross sections of Upsilon(1S) as a function of pT.

Differential production cross sections of Upsilon(1S) as a function of rapidity.

integrated production cross section of Upsilon (1S).

integrated production cross section of Upsilon (2S).

integrated Upsilon(1S) over Upsilon(2s) production cross section ratio.

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.