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.

Summary of results for cross-section of prompt psi(2S) decaying to a muon pair in pp collisions at 5.02 TeV in nb/GeV. Uncertainties are statistical and systematic, respectively.

Summary of results for cross-section of Upsilon(1S) decaying to a muon pair in pp collisions at 5.02 TeV in nb/GeV. Uncertainties are statistical and systematic, respectively.

Summary of results for cross-section of Upsilon(2S) decaying to a muon pair in pp collisions at 5.02 TeV in nb/GeV. Uncertainties are statistical and systematic, respectively.

Summary of results for cross-section of Upsilon(3S) decaying to a muon pair in pp collisions at 5.02 TeV in nb/GeV. Uncertainties are statistical and systematic, respectively.

Summary of results for cross-section of J/psi decaying to a muon pair in p+Pb collisions at 5.02 TeV in nb/GeV. Uncertainties are statistical and systematic, respectively.

Summary of results for cross-section of psi(2S) decaying to a muon pair in p+Pb collisions at 5.02 TeV in nb/GeV. Uncertainties are statistical and systematic, respectively.

Summary of results for cross-section of J/psi decaying to a muon pair in p+Pb collisions at 5.02 TeV as a function of center-of-mass rapdiity in nb/GeV. Uncertainties are statistical and systematic, respectively.

Summary of results for cross-section of psi(2S) decaying to a muon pair in p+Pb collisions at 5.02 TeV as a function of center-of-mass rapdiity in nb/GeV. Uncertainties are statistical and systematic, respectively.

Summary of results for cross-section of Upsilon(nS) decaying to a muon pair in p+Pb collisions at 5.02 TeV in nb/GeV. Uncertainties are statistical and systematic, respectively.

Summary of results for cross-section of Upsilon(nS) decaying to a muon pair in p+Pb collisions at 5.02 TeV in nb/GeV. Uncertainties are statistical and systematic, respectively.

Summary of results for RpPb of prompt J/psi in p+Pb collisions at 5.02 TeV as a function of pT. Uncertainties are statistical and local systematic and global systematic, respectively.

Summary of results for RpPb of non-prompt J/psi in p+Pb collisions at 5.02 TeV as a function of pT. Uncertainties are statistical and local systematic and global systematic, respectively.

Summary of results for RpPb of prompt J/psi in p+Pb collisions at 5.02 TeV as a function of ystar. Uncertainties are statistical and local systematic and global systematic, respectively.

Summary of results for RpPb of non-prompt J/psi in p+Pb collisions at 5.02 TeV as a function of ystar. Uncertainties are statistical and local systematic and global systematic, respectively.

Summary of results for RpPb of Upsilon(1S) in p+Pb collisions at 5.02 TeV as a function of pT. Uncertainties are statistical and local systematic and global systematic, respectively.

Summary of results for RpPb of Upsilon(1S) in p+Pb collisions at 5.02 TeV as a function of ystar. Uncertainties are statistical and local systematic and global systematic, respectively.

Summary of results for RpPb of quarkonia (prompt J/psi, non-prompt J/psi, prompt psi(2S), Upsilon(1S)) to RpPb of Z ratio in p+Pb collisions at 5.02 TeV as a function of centrality. Uncertainties are statistical and local systematic and global systematic, respectively.

Summary of results for quarkonia self-normalized yields in p+Pb collisions at 5.02 TeV as a function of self-normalized event activity. Uncertainties are statistical and systematic, respectively.

Summary of results for prompt Psi(2S) to J/psi double ratio in p+Pb collisions at 5.02 TeV as a function of center-of-mass rapidity. Uncertainties are statistical and systematic, respectively.

Summary of results for Upsilon(2S) and Upsilon(3S) to Upsilon(1S) double ratio in p+Pb collisions at 5.02 TeV. Uncertainties are statistical and systematic, respectively.

Summary of results for prompt Psi(2S) and J/psi double ratio in p+Pb collisions at 5.02 TeV as a function of centrality. Uncertainties are statistical and local systematic and global systematic, respectively.

Summary of results for Upsilon(2S) and Upsilon(3S) to Upsilon(1S) double ratio in p+Pb collisions at 5.02 TeV as a function of centrality. Uncertainties are statistical and local systematic and global systematic, respectively.

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.

Prompt J/$\psi$ $v_{2}$ as a function of rapidity.

Prompt J/$\psi$ $v_{2}$ as a function of $p_{T}$.

Prompt J/$\psi$ $v_{2}$ as a function of $p_{T}$.

Prompt J/$\psi$ $R_{\rm AA}$ as a function of centrality. The average ${N}_{\rm part}$ values correspond to events flatly distributed across centrality.

Prompt J/$\psi$ $R_{\rm AA}$ as a function of centrality. The average ${N}_{\rm part}$ values correspond to events flatly distributed across centrality.

Prompt J/$\psi$ $R_{\rm AA}$ as a function of rapidity.

Prompt J/$\psi$ $R_{\rm AA}$ as a function of rapidity.

Prompt J/$\psi$ $R_{\rm AA}$ as a function of $p_{\rm T}$.

Prompt J/$\psi$ $R_{\rm AA}$ as a function of $p_{\rm T}$.

Prompt J/$\psi$ $R_{\rm AA}$ as a function of centrality at high $p_{\rm T}$, 6.5$<p_{\rm T}<$30 GeV/c, in three different $|y|$ regions.

Prompt J/$\psi$ $R_{\rm AA}$ as a function of centrality at high $p_{\rm T}$, 6.5$<p_{\rm T}<$30 GeV/c, in three different $|y|$ regions.

Prompt J/$\psi$ $R_{\rm AA}$ as a function of centrality, at forward rapidity, 1.6$<|y|<$2.4, for two different $p_{\rm T}$ regions.

Prompt J/$\psi$ $R_{\rm AA}$ as a function of centrality, at forward rapidity, 1.6$<|y|<$2.4, for two different $p_{\rm T}$ regions.

Nonprompt J/$\psi$ $v_{2}$ as a function of $p_{\rm T}$.

Nonprompt J/$\psi$ $v_{2}$ as a function of $p_{\rm T}$.

Nonprompt J/$\psi$ $R_{\rm AA}$ as a function of centrality. The average ${N}_{\rm part}$ values correspond to events flatly distributed across centrality.

Nonprompt J/$\psi$ $R_{\rm AA}$ as a function of centrality. The average ${N}_{\rm part}$ values correspond to events flatly distributed across centrality.

Nonprompt J/$\psi$ $R_{\rm AA}$ as a function of rapidity.

Nonprompt J/$\psi$ $R_{\rm AA}$ as a function of rapidity.

Nonprompt J/$\psi$ $R_{\rm AA}$ as a function of $p_{\rm T}$.

Nonprompt J/$\psi$ $R_{\rm AA}$ as a function of $p_{\rm T}$.

Nonprompt J/$\psi$ $R_{\rm AA}$ as a function of centrality. The average ${N}_{\rm part}$ values correspond to events flatly distributed across centrality.

Nonprompt J/$\psi$ $R_{\rm AA}$ as a function of centrality. The average ${N}_{\rm part}$ values correspond to events flatly distributed across centrality.

Nonprompt J/$\psi$ $R_{\rm AA}$ as a function of centrality. The average ${N}_{\rm part}$ values correspond to events flatly distributed across centrality.

Nonprompt J/$\psi$ $R_{\rm AA}$ as a function of centrality. The average ${N}_{\rm part}$ values correspond to events flatly distributed across centrality.

Summary of results for cross-section of non-prompt $J/\psi$ decaying to a muon pair for 8 TeV data in nb/GeV. Uncertainties are statistical and systematic, respectively.

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

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

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

Summary of results for cross-section of non-prompt $\psi(2S)$ decaying to a muon pair for 8 TeV data in nb/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 7 TeV data. 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 8 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 7 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 8 TeV data. Uncertainties are statistical and systematic, respectively.

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

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

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

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

Mean weight correction factor for $J/\psi$ under the "longitudinal" spin-alignment hypothesis for 7 TeV.

Mean weight correction factor for $J/\psi$ under the "transverse zero" spin-alignment hypothesis for 7 TeV.

Mean weight correction factor for $J/\psi$ under the "transverse positive" spin-alignment transverse positive hypothesis for 7 TeV.

Mean weight correction factor for $J/\psi$ under the "transverse negative" spin-alignment hypothesis for 7 TeV.

Mean weight correction factor for $J/\psi$ under the "off-($\lambda_{\theta}$--$\lambda_{\phi}$)-plane positive" spin-alignment hypothesis for 7 TeV.

Mean weight correction factor for $J/\psi$ under the "off-($\lambda_{\theta}$--$\lambda_{\phi}$)-plane negative" spin-alignment hypothesis for 7 TeV.

Mean weight correction factor for $\psi(2S)$ under the "longitudinal" spin-alignment hypothesis for 7 TeV.

Mean weight correction factor for $\psi(2S)$ under the "transverse zero" spin-alignment hypothesis for 7 TeV.

Mean weight correction factor for $\psi(2S)$ under the "transverse positive" spin-alignment hypothesis for 7 TeV.

Mean weight correction factor for $\psi(2S)$ under the "transverse negative" spin-alignment hypothesis for 7 TeV.

Mean weight correction factor for $\psi(2S)$ under the "off-($\lambda_{\theta}$--$\lambda_{\phi}$)-plane positive" spin-alignment hypothesis for 7 TeV.

Mean weight correction factor for $\psi(2S)$ under the "off-($\lambda_{\theta}$--$\lambda_{\phi}$)-plane negative" spin-alignment hypothesis for 7 TeV.

Mean weight correction factor for $J/\psi$ under the "longitudinal" spin-alignment hypothesis for 8 TeV. Those intervals not measured in the analysis at low $p_{T}$, high rapidity are also excluded here.

Mean weight correction factor for $J/\psi$ under the "transverse zero" spin-alignment hypothesis for 8 TeV. Those intervals not measured in the analysis at low $p_{T}$, high rapidity are also excluded here.

Mean weight correction factor for $J/\psi$ under the "transverse positive" spin-alignment hypothesis for 8 TeV. Those intervals not measured in the analysis at low $p_{T}$, high rapidity are also excluded here.

Mean weight correction factor for $J/\psi$ under the "transverse negative" spin-alignment hypothesis for 8 TeV. Those intervals not measured in the analysis at low $p_{T}$, high rapidity are also excluded here.

Mean weight correction factor for $J/\psi$ under the "off-($\lambda_{\theta}$--$\lambda_{\phi}$)-plane positive" spin-alignment hypothesis for 8 TeV. Those intervals not measured in the analysis at low $p_{T}$, high rapidity are also excluded here.

Mean weight correction factor for $J/\psi$ under the "off-($\lambda_{\theta}$--$\lambda_{\phi}$)-plane negative" spin-alignment hypothesis for 8 TeV. Those intervals not measured in the analysis at low $p_{T}$, high rapidity are also excluded here.

Mean weight correction factor for $\psi(2S)$ under the "longitudinal" spin-alignment hypothesis for 8 TeV. Those intervals not measured in the analysis at low $p_{T}$, high rapidity are also excluded here.

Mean weight correction factor for $\psi(2S)$ under the "transverse zero" spin-alignment hypothesis for 8 TeV. Those intervals not measured in the analysis at low $p_{T}$, high rapidity are also excluded here.

Mean weight correction factor for $\psi(2S)$ under the "transverse positive" spin-alignment hypothesis for 8 TeV. Those intervals not measured in the analysis at low $p_{T}$, high rapidity are also excluded here.

Mean weight correction factor for $\psi(2S)$ under the "transverse negative" spin-alignment hypothesis for 8 TeV. Those intervals not measured in the analysis at low $p_{T}$, high rapidity are also excluded here.

Mean weight correction factor for $\psi(2S)$ under the "off-($\lambda_{\theta}$--$\lambda_{\phi}$)-plane positive" spin-alignment hypothesis for 8 TeV. Those intervals not measured in the analysis at low $p_{T}$, high rapidity are also excluded here.

Mean weight correction factor for $\psi(2S)$ under the "off-($\lambda_{\theta}$--$\lambda_{\phi}$)-plane negative" spin-alignment hypothesis for 8 TeV. Those intervals not measured in the analysis at low $p_{T}$, high rapidity are also excluded here.

Double differential cross-section for $J/\psi$-from-$b$ mesons as a function of $p_\perp$ in bins of $y$. The first uncertainties are statistical, the second are the correlated systematic uncertainties shared between bins and the last are the uncorrelated systematic uncertainties.

The fraction of $J/\psi$-from-$b$ mesons (in %) in bins of the $J/\psi$ $p_\perp$ and $y$. The uncertainties are statistical only. The systematic uncertainties are negligible.

The fraction of $J/\psi$-from-$b$ mesons (in %) in bins of the $J/\psi$ $p_\perp$ and $y$. The uncertainties are statistical only. The systematic uncertainties are negligible.

Differential cross-sections as a function of $p_\perp$ integrated over $y$ for prompt $J/\psi$ mesons. The first uncertainties are statistical and the second (third) are uncorrelated (correlated) systematic uncertainties amongst bins.

Differential cross-sections as a function of $p_\perp$ integrated over $y$ for prompt $J/\psi$ mesons. The first uncertainties are statistical and the second (third) are uncorrelated (correlated) systematic uncertainties amongst bins.

Differential cross-sections as a function of $p_\perp$ integrated over $y$ for $J/\psi$-from-$b$ mesons. The first uncertainties are statistical and the second (third) are uncorrelated (correlated) systematic uncertainties amongst bins.

Differential cross-sections as a function of $p_\perp$ integrated over $y$ for $J/\psi$-from-$b$ mesons. The first uncertainties are statistical and the second (third) are uncorrelated (correlated) systematic uncertainties amongst bins.

Differential cross-sections as a function of $y$ integrated over $p_\perp$ for prompt $J/\psi$ mesons. The first uncertainties are statistical and the second (third) are the correlated (uncorrelated) systematic uncertainties.

Differential cross-sections as a function of $y$ integrated over $p_\perp$ for prompt $J/\psi$ mesons. The first uncertainties are statistical and the second (third) are the correlated (uncorrelated) systematic uncertainties.

Differential cross-sections as a function of $y$ integrated over $p_\perp$ for $J/\psi$-from-$b$ mesons. The first uncertainties are statistical and the second (third) are the correlated (uncorrelated) systematic uncertainties.

Differential cross-sections as a function of $y$ integrated over $p_\perp$ for $J/\psi$-from-$b$ mesons. The first uncertainties are statistical and the second (third) are the correlated (uncorrelated) systematic uncertainties.

Ratios of differential cross-sections between measurements at $\sqrt{s} = 13$ TeV and $\sqrt{s} = 8$ TeV as a function of $p_\perp$ in bins of $y$ for prompt $J/\psi$ mesons. The systematic errors are negligible.

Ratios of differential cross-sections between measurements at $\sqrt{s} = 13$ TeV and $\sqrt{s} = 8$ TeV as a function of $p_\perp$ in bins of $y$ for prompt $J/\psi$ mesons. The systematic errors are negligible.

Ratios of differential cross-sections between measurements at $\sqrt{s} = 13$ TeV and $\sqrt{s} = 8$ TeV as a function of $p_\perp$ in bins of $y$ for $J/\psi$-from-$b$ mesons. The systematic errors are negligible.

Ratios of differential cross-sections between measurements at $\sqrt{s} = 13$ TeV and $\sqrt{s} = 8$ TeV as a function of $p_\perp$ in bins of $y$ for $J/\psi$-from-$b$ mesons. The systematic errors are negligible.

Ratios of differential cross-sections between measurements at $\sqrt{s} = 13$ TeV and $\sqrt{s} = 8$ TeV as a function of $y$ integrated over $p_\perp$ for prompt $J/\psi$ mesons. The systematic errors are negligible.

Ratios of differential cross-sections between measurements at $\sqrt{s} = 13$ TeV and $\sqrt{s} = 8$ TeV as a function of $y$ integrated over $p_\perp$ for prompt $J/\psi$ mesons. The systematic errors are negligible.

Ratios of differential cross-sections between measurements at $\sqrt{s} = 13$ TeV and $\sqrt{s} = 8$ TeV as a function of $y$ integrated over $p_\perp$ for $J/\psi$-from-$b$ mesons. The systematic errors are negligible.

Ratios of differential cross-sections between measurements at $\sqrt{s} = 13$ TeV and $\sqrt{s} = 8$ TeV as a function of $y$ integrated over $p_\perp$ for $J/\psi$-from-$b$ mesons. The systematic errors are negligible.

Ratios of differential cross-sections between measurements at $\sqrt{s} = 13$ TeV and $\sqrt{s} = 8$ TeV as a function of $p_\perp$ integrated over $y$ for prompt $J/\psi$ mesons.

Ratios of differential cross-sections between measurements at $\sqrt{s} = 13$ TeV and $\sqrt{s} = 8$ TeV as a function of $p_\perp$ integrated over $y$ for prompt $J/\psi$ mesons.

Ratios of differential cross-sections between measurements at $\sqrt{s} = 13$ TeV and $\sqrt{s} = 8$ TeV as a function of $p_\perp$ integrated over $y$ for $J/\psi$-from-$b$ mesons. The systematic errors are negligible.

Ratios of differential cross-sections between measurements at $\sqrt{s} = 13$ TeV and $\sqrt{s} = 8$ TeV as a function of $p_\perp$ integrated over $y$ for $J/\psi$-from-$b$ mesons. The systematic errors are negligible.

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.

Differential cross-section for non-prompt chi_c2 production, measured in bins of J/psi pT, assuming unpolarised chi_c production. The measurements are not corrected for the branching fractions of the decays chi_c --> J/psi + gamma and J/psi --> mu+ mu-. The uncertainty envelope associated with the unknown chi_c spin alignment is also shown.

Differential cross-section for prompt chi_c1 production, measured in bins of chi_c1 pT, assuming unpolarised chi_c production. The measurements are not corrected for the branching fractions of the decays chi_c --> J/psi + gamma and J/psi --> mu+ mu-. The uncertainty envelope associated with the unknown chi_c spin alignment is also shown.

Differential cross-section for prompt chi_c2 production, measured in bins of chi_c2 pT, assuming unpolarised chi_c production. The measurements are not corrected for the branching fractions of the decays chi_c --> J/psi + gamma and J/psi --> mu+ mu-. The uncertainty envelope associated with the unknown chi_c spin alignment is also shown.

Differential cross-section for non-prompt chi_c1 production, measured in bins of chi_c1 pT, assuming unpolarised chi_c production. The measurements are not corrected for the branching fractions of the decays chi_c --> J/psi + gamma and J/psi --> mu+ mu-. The uncertainty envelope associated with the unknown chi_c spin alignment is also shown.

Differential cross-section for non-prompt chi_c2 production, measured in bins of chi_c2 pT, assuming unpolarised chi_c production. The measurements are not corrected for the branching fractions of the decays chi_c --> J/psi + gamma and J/psi --> mu+ mu-. The uncertainty envelope associated with the unknown chi_c spin alignment is also shown.

Fraction of prompt J/psi produced in feed-down from radiative chi_c decays as a function of J/psi pT, assuming unpolarised chi_c production. The uncertainty envelope associated with the unknown chi_c spin alignment is also shown. The spin-alignment envelope assumes that prompt J/psi are produced unpolarised and represents the maximum uncertainty due to the unknown chi_c spin alignment.

Production rate of prompt chi_c2 relative to prompt chi_c1, measured in bins of J/psi pT, assuming unpolarised chi_c production. The ratio is not corrected for the branching fractions of chi_c1 and chi_c2 into J/psi gamma. The uncertainty envelope associated with the unknown chi_c spin alignment is also shown.

Production rate of non-prompt chi_c2 relative to non-prompt chi_c1, measured in bins of J/psi pT, assuming unpolarised chi_c production. The ratio is not corrected for the branching fractions of chi_c1 and chi_c2 into J/psi gamma. The uncertainty envelope associated with the unknown chi_c spin alignment is also shown.

Fraction of chi_c1 produced in b-hadron decays as a function of chi_c1 pT, assuming unpolarised chi_c production. The uncertainty envelope associated with the unknown chi_c spin alignment is also shown.

Fraction of chi_c2 produced in b-hadron decays as a function of chi_c2 pT, assuming unpolarised chi_c production. The uncertainty envelope associated with the unknown chi_c spin alignment is also shown.

Double-differential cross-section $d^2\sigma/dp_T/dy$ for prompt $J/\psi$ production in bins of $p_T$ and $y$, with statistical, systematic and polarization uncertainties.