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 $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 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$ $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 $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.

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.

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.

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 ;.