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.

Ratio of measured PSI(2S) over J/PSI(1S) cross section in charged-particle multiplicity intervals and integrated in multiplicity.

Ratio of measured PSI(2S) over J/PSI(1S) cross section in charged-particle multiplicity intervals and integrated in multiplicity.

Ratio of measured PSI(2S) over J/PSI(1S) cross section in charged-particle multiplicity intervals and integrated in multiplicity.

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.

Double-differential cross section times the dimuon branching fraction of the Y(2S) meson for different ranges of pT in bins of |y| and for the full |y| < 1.2 range, for the unpolarized decay hypothesis. The global uncertainty in the integrated luminosity of 2.3% is not included in the systematic uncertainties.

Double-differential cross section times the dimuon branching fraction of the Y(3S) meson for different ranges of pT in bins of |y| and for the full |y| < 1.2 range, for the unpolarized decay hypothesis. The global uncertainty in the integrated luminosity of 2.3% is not included in the systematic uncertainties.

Multiplicative scaling factors to obtain the J/psi differential cross sections for different polarization scenarios (Lambda_{Theta} = +1, +0.10, -1) from the unpolarized cross section measurements given in Table 1.

Multiplicative scaling factors to obtain the psi(2S) differential cross sections for different polarization scenarios (Lambda_{Theta} = +1, +0.03, -1) from the unpolarized cross section measurements given in Table 2.

Multiplicative scaling factors to obtain the Y(1S) differential cross sections for different polarization scenarios (Lambda_{Theta} = +1, k, -1) from the unpolarized cross section measurements given in Table 3. The parameter k corresponds to a linear interpolation of the CMS measured value as a function of pT for pT < 50 GeV. For pT > 50 GeV, where no measurements exist, k is taken as the average of all the measured values for pT < 50 GeV.

Multiplicative scaling factors to obtain the Y(2S) differential cross sections for different polarization scenarios (Lambda_{Theta} = +1, k, -1) from the unpolarized cross section measurements given in Table 3. The parameter k corresponds to a linear interpolation of the CMS measured value as a function of pT for pT < 50 GeV. For pT > 50 GeV, where no measurements exist, k is taken as the average of all the measured values for pT < 50 GeV.

Multiplicative scaling factors to obtain the Y(3S) differential cross sections for different polarization scenarios (Lambda_{Theta} = +1, k, -1) from the unpolarized cross section measurements given in Table 3. The parameter k corresponds to a linear interpolation of the CMS measured value as a function of pT for pT < 50 GeV. For pT > 50 GeV, where no measurements exist, k is taken as the average of all the measured values for pT < 50 GeV.

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

Ratios of differential cross section times branching fractions of the prompt Y(2S) to Y(1S) and Y(3S) to Y(1S) mesons assuming unpolarized productions as a function of pT for |y| < 1.2.

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.

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.

The measured prompt PSI(2S) polarisation parameter LAMBDA(INV) in bins of YRAP and PT in the helicity frame.

The measured prompt PSI(2S) polarisation parameter LAMBDA(THETA) in bins of YRAP and PT in the Collins-Soper frame.

The measured prompt PSI(2S) polarisation parameter LAMBDA(THETA PHI) in bins of YRAP and PT in the Collins-Soper frame.

The measured prompt PSI(2S) polarisation parameter LAMBDA(PHI) in bins of YRAP and PT in the Collins-Soper frame.

The measured prompt PSI(2S) polarisation parameter LAMBDA(INV) in bins of YRAP and PT in the Collins-Soper frame.

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.

Total nonprompt X(3872) fraction without acceptance corrections.

Differential cross section for prompt P P --> X + anything times branching fraction X-> J/PSI pi+ pi-.

Cross section for prompt P P --> X + anything times branching fraction X-> J/PSI pi+ pi-.

Dipion invariant-mass spectrum for X(3872) ->J/PSI pi+pi- decays corrected for acceptance and efficiency, dividing by the total cross section for 0.5 < M(PI+PI-) < 0.78 GEV.