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.

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.

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.

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.