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.

Measured DPS effective cross section

pp -> J/psi J/psi J/psi X cross section

$pp \rightarrow J/\psi J/\psi J/\psi X~$ fiducial cross section

$v_{2}$ of $\Upsilon(\mathrm{1S})$ as a function of $p_{\mathrm{T}}$ in three centrality ranges, 10-30%, 30-50% and 50-90%.

$v_{2}$ of $\Upsilon(\mathrm{1S})$ as a function of $p_{\mathrm{T}}$ in 5-60% centrality range.

Summary of the measured values of $ d\sigma / d{m}$ (pb/GeV) in the dielectron channel with the statistical ($\delta_{\text{stat}}$), experimental ($\delta_{\text{exp}}$) and theoretical ($\delta_{\text{theo}}$) uncertainties, respectively. Here, $\delta_{\text{tot}}$ is the quadratic sum of the three components.

Summary of the measured values of fiducial $ d\sigma / d{m}$ (pb/GeV) (with no FSR correction applied) in the dimuon channel with the statistical ($\delta_{\text{stat}}$) and experimental ($\delta_{\text{exp}}$) uncertainties shown separately. Here, $\delta_{\text{tot}}$ is the quadratic sum of the two components.

of the measured values of fiducial $ d\sigma / d{m}$ (pb/GeV) (with no FSR correction applied) in the dielectron channel with the statistical ($\delta_{\text{stat}}$) and experimental ($\delta_{\text{exp}}$) uncertainties shown separately. Here, $\delta_{\text{tot}}$ is the quadratic sum of the two components.

Summary of the combined values of $ d\sigma / d{m}$ (pb/GeV) using the results from both the dimuon and dielectron channels. Here, $\delta_{\text{tot}}$ is the quadratic sum of the statistical, experimental and theoretical uncertainties.

Covariance matrix of the Drell-Yan differential cross section as a function of the dimuon invariant mass, measured in the full phase space, with FSR correction is applied. Luminosity uncertainty is not included in the covariance.

Covariance matrix of the Drell-Yan differential cross section as a function of the dielectron invariant mass, measured in the full phase space, with FSR correction is applied. Luminosity uncertainty is not included in the covariance.

Covariance matrix of the differential DY cross section measured for the combination of the two channels in the full phase space.

Differential cross sections of the Y(2S) meson as a function of pT for pp collisions. The global uncertainty arises from the integrated luminosity uncertainty in pp collisions.

Differential cross sections of the Y(3S) meson as a function of pT for pp collisions. The global uncertainty arises from the integrated luminosity uncertainty in pp collisions.

Differential cross sections of the Y(3S) meson as a function of pT for pp collisions. The global uncertainty arises from the integrated luminosity uncertainty in pp collisions.

Differential cross sections of the Y(1S) meson as a function of pT for PbPb collisions. Global uncertainties arise from the number of minimum bias events and TAA uncertainties for PbPb collisions.

Differential cross sections of the Y(1S) meson as a function of pT for PbPb collisions. Global uncertainties arise from the number of minimum bias events and TAA uncertainties for PbPb collisions.

Differential cross sections of the Y(2S) meson as a function of pT for PbPb collisions. Global uncertainties arise from the number of minimum bias events and TAA uncertainties for PbPb collisions.

Differential cross sections of the Y(2S) meson as a function of pT for PbPb collisions. Global uncertainties arise from the number of minimum bias events and TAA uncertainties for PbPb collisions.

Differential cross sections of the Y(3S) meson as a function of pT for PbPb collisions.

Differential cross sections of the Y(3S) meson as a function of pT for PbPb collisions.

Differential cross sections of the Y(1S) meson as a function of rapidity for pp collisions. The global uncertainty arises from the integrated luminosity uncertainty in pp collisions.

Differential cross sections of the Y(1S) meson as a function of rapidity for pp collisions. The global uncertainty arises from the integrated luminosity uncertainty in pp collisions.

Differential cross sections of the Y(2S) meson as a function of rapidity for pp collisions. The global uncertainty arises from the integrated luminosity uncertainty in pp collisions.

Differential cross sections of the Y(2S) meson as a function of rapidity for pp collisions. The global uncertainty arises from the integrated luminosity uncertainty in pp collisions.

Differential cross sections of the Y(3S) meson as a function of rapidity for pp collisions. The global uncertainty arises from the integrated luminosity uncertainty in pp collisions.

Differential cross sections of the Y(3S) meson as a function of rapidity for pp collisions. The global uncertainty arises from the integrated luminosity uncertainty in pp collisions.

Differential cross sections of the Y(1S) meson as a function of pT for PbPb collisions. Global uncertainties arise from the number of minimum bias events and TAA uncertainties for PbPb collisions.

Differential cross sections of the Y(1S) meson as a function of pT for PbPb collisions. Global uncertainties arise from the number of minimum bias events and TAA uncertainties for PbPb collisions.

Differential cross sections of the Y(2S) meson as a function of pT for PbPb collisions. Global uncertainties arise from the number of minimum bias events and TAA uncertainties for PbPb collisions.

Differential cross sections of the Y(2S) meson as a function of pT for PbPb collisions. Global uncertainties arise from the number of minimum bias events and TAA uncertainties for PbPb collisions.

Differential cross sections of the Y(3S) meson as a function of pT for PbPb collisions.

Differential cross sections of the Y(3S) meson as a function of pT for PbPb collisions.

Nuclear modification factor of the Y(1S) meson as a function of pT. Global uncertainty combines the uncertainties from TAA, pp luminosity, and number of minimum bias events in PbPb collisions.

Nuclear modification factor of the Y(1S) meson as a function of pT. Global uncertainty combines the uncertainties from TAA, pp luminosity, and number of minimum bias events in PbPb collisions.

Nuclear modification factor of the Y(2S) meson as a function of pT. Global uncertainty combines the uncertainties from TAA, pp luminosity, and number of minimum bias events in PbPb collisions.

Nuclear modification factor of the Y(2S) meson as a function of pT. Global uncertainty combines the uncertainties from TAA, pp luminosity, and number of minimum bias events in PbPb collisions.

Nuclear modification factor of the Y(3S) meson as a function of pT.

Nuclear modification factor of the Y(3S) meson as a function of pT.

Nuclear modification factor of the Y(1S) meson as a function of rapidity. Global uncertainty combines the uncertainties from TAA, pp luminosity, and number of minimum bias events in PbPb collisions.

Nuclear modification factor of the Y(1S) meson as a function of rapidity. Global uncertainty combines the uncertainties from TAA, pp luminosity, and number of minimum bias events in PbPb collisions.

Nuclear modification factor of the Y(2S) meson as a function of rapidity. Global uncertainty combines the uncertainties from TAA, pp luminosity, and number of minimum bias events in PbPb collisions.

Nuclear modification factor of the Y(2S) meson as a function of rapidity. Global uncertainty combines the uncertainties from TAA, pp luminosity, and number of minimum bias events in PbPb collisions.

Nuclear modification factor of the Y(3S) meson as a function of rapidity.

Nuclear modification factor of the Y(3S) meson as a function of rapidity.

Nuclear modification factor of the Y(1S) meson as a function of collision centrality. Global uncertainty combines the pp luminosity and number of minimum bias events in PbPb collisions. The PP Y(1S) global uncertainty consists of the statistical and systematic uncertainties of pp yield for Y(1S) state.

Nuclear modification factor of the Y(1S) meson as a function of collision centrality. Global uncertainty combines the pp luminosity and number of minimum bias events in PbPb collisions. The PP Y(1S) global uncertainty consists of the statistical and systematic uncertainties of pp yield for Y(1S) state.

Nuclear modification factor of the Y(2S) meson as a function of collision centrality. Global uncertainty combines the pp luminosity and number of minimum bias events in PbPb collisions. The PP Y(2S) global uncertainty consists of the statistical and systematic uncertainties of pp yield for Y(2S) state.

Nuclear modification factor of the Y(2S) meson as a function of collision centrality. Global uncertainty combines the pp luminosity and number of minimum bias events in PbPb collisions. The PP Y(2S) global uncertainty consists of the statistical and systematic uncertainties of pp yield for Y(2S) state.

Confidence intervals on the nuclear modification factors of the Y(3S) mesons as a function of collision centrality.

Confidence intervals on the nuclear modification factors of the Y(3S) mesons as a function of collision centrality.

Confidence intervals on the nuclear modification factors of the Y(2S) mesons as a function of collision centrality.

Confidence intervals on the nuclear modification factors of the Y(2S) mesons as a function of collision centrality.

Nuclear modification factors of the Y(1S) and Y(2S) mesons in the integrated pT, rapidity and centrality.

Nuclear modification factors of the Y(1S) and Y(2S) mesons in the integrated pT, rapidity and centrality.

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.

Values of $g_1^{NS}$ for the COMPASS data in ($x$, $Q^2$) bins.

Parameters of two solutions of the QCD fit. The assumed functional form at $Q^2_0 = 1 ($GeV$/c)^2$ is $q(x) = \eta x^\alpha (1-x)^\beta (1 + \gamma x)/N$ where N is the normalisation. The parametrisation is used for the singlet and gluon helicity distribution and the two non-singlet distributions $\Delta q_3$ and $\Delta q_8$.

Normalisations of the different data sets used in the QCD fit.

Differential cross section for Y(2S) states as a function of their transverse momentum and per unit of rapidity in pp collisions. Statistical (systematic) uncertainties are displayed as error bars (boxes). Global relative uncertainties of 3.7% are not displayed.

Differential cross section for Y(3S) states as a function of their transverse momentum and per unit of rapidity in pp collisions. Statistical (systematic) uncertainties are displayed as error bars (boxes). Global relative uncertainties of 3.7% are not displayed.

Differential cross section for Y(3S) states as a function of their transverse momentum and per unit of rapidity in pp collisions. Statistical (systematic) uncertainties are displayed as error bars (boxes). Global relative uncertainties of 3.7% are not displayed.

Differential cross section for Y(1S) states as a function of their transverse momentum and per unit of rapidity in PbPb collisions. Statistical (systematic) uncertainties are displayed as error bars (boxes). Global relative uncertainties of 6.5% are not displayed.

Differential cross section for Y(1S) states as a function of their transverse momentum and per unit of rapidity in PbPb collisions. Statistical (systematic) uncertainties are displayed as error bars (boxes). Global relative uncertainties of 6.5% are not displayed.

Differential cross section for Y(2S) states as a function of their transverse momentum and per unit of rapidity in PbPb collisions. Statistical (systematic) uncertainties are displayed as error bars (boxes). Global relative uncertainties of 6.5% are not displayed.

Differential cross section for Y(2S) states as a function of their transverse momentum and per unit of rapidity in PbPb collisions. Statistical (systematic) uncertainties are displayed as error bars (boxes). Global relative uncertainties of 6.5% are not displayed.

Differential cross section for Y(1S) states as a function of their rapidity and integrated over transverse momentum in pp collisions. Statistical (systematic) uncertainties are displayed as error bars (boxes). Global relative uncertainties of 3.7% are not displayed.

Differential cross section for Y(1S) states as a function of their rapidity and integrated over transverse momentum in pp collisions. Statistical (systematic) uncertainties are displayed as error bars (boxes). Global relative uncertainties of 3.7% are not displayed.

Differential cross section for Y(2S) states as a function of their rapidity and integrated over transverse momentum in pp collisions. Statistical (systematic) uncertainties are displayed as error bars (boxes). Global relative uncertainties of 3.7% are not displayed.

Differential cross section for Y(2S) states as a function of their rapidity and integrated over transverse momentum in pp collisions. Statistical (systematic) uncertainties are displayed as error bars (boxes). Global relative uncertainties of 3.7% are not displayed.

Differential cross section for Y(3S) states as a function of their rapidity and integrated over transverse momentum in pp collisions. Statistical (systematic) uncertainties are displayed as error bars (boxes). Global relative uncertainties of 3.7% are not displayed.

Differential cross section for Y(3S) states as a function of their rapidity and integrated over transverse momentum in pp collisions. Statistical (systematic) uncertainties are displayed as error bars (boxes). Global relative uncertainties of 3.7% are not displayed.

Differential cross section for Y(1S) states as a function of their rapidity and integrated over transverse momentum in PbPb collisions. Statistical (systematic) uncertainties are displayed as error bars (boxes). Global relative uncertainties of 6.5% are not displayed.

Differential cross section for Y(1S) states as a function of their rapidity and integrated over transverse momentum in PbPb collisions. Statistical (systematic) uncertainties are displayed as error bars (boxes). Global relative uncertainties of 6.5% are not displayed.

Differential cross section for Y(2S) states as a function of their rapidity and integrated over transverse momentum in PbPb collisions. Statistical (systematic) uncertainties are displayed as error bars (boxes). Global relative uncertainties of 6.5% are not displayed.

Differential cross section for Y(2S) states as a function of their rapidity and integrated over transverse momentum in PbPb collisions. Statistical (systematic) uncertainties are displayed as error bars (boxes). Global relative uncertainties of 6.5% are not displayed.

Nuclear modification factor for Υ(1S) states in PbPb collisions as a function of pT. Statistical (systematic) uncertainties are displayed as error bars (boxes), while the global (fully correlated) uncertainty (7.5%) is displayed as a grey box at unity.

Nuclear modification factor for Υ(1S) states in PbPb collisions as a function of pT. Statistical (systematic) uncertainties are displayed as error bars (boxes), while the global (fully correlated) uncertainty (7.5%) is displayed as a grey box at unity.

Nuclear modification factor for Υ(2S) states in PbPb collisions as a function of pT. Statistical (systematic) uncertainties are displayed as error bars (boxes), while the global (fully correlated) uncertainty (7.5%) is displayed as a grey box at unity.

Nuclear modification factor for Υ(2S) states in PbPb collisions as a function of pT. Statistical (systematic) uncertainties are displayed as error bars (boxes), while the global (fully correlated) uncertainty (7.5%) is displayed as a grey box at unity.

Nuclear modification factor for Υ(1S) states in PbPb collisions as a function of |y|. Statistical (systematic) uncertainties are displayed as error bars (boxes), while the global (fully correlated) uncertainty (7.5%) is displayed as a grey box at unity.

Nuclear modification factor for Υ(1S) states in PbPb collisions as a function of |y|. Statistical (systematic) uncertainties are displayed as error bars (boxes), while the global (fully correlated) uncertainty (7.5%) is displayed as a grey box at unity.

Nuclear modification factor for Υ(2S) states in PbPb collisions as a function of |y|. Statistical (systematic) uncertainties are displayed as error bars (boxes), while the global (fully correlated) uncertainty (7.5%) is displayed as a grey box at unity.

Nuclear modification factor for Υ(2S) states in PbPb collisions as a function of |y|. Statistical (systematic) uncertainties are displayed as error bars (boxes), while the global (fully correlated) uncertainty (7.5%) is displayed as a grey box at unity.

Nuclear modification factors for Υ(1S) meson production in PbPb collisions, as a function of centrality, displayed as the average number of participating nucleons. Statistical (systematic) uncertainties are displayed as error bars (boxes), while the global (fully correlated) uncertainties from the PbPb data (3.2%) or from the pp reference (6.3%) are displayed at unity as empty and filled red black boxes, respectively.

Nuclear modification factors for Υ(1S) meson production in PbPb collisions, as a function of centrality, displayed as the average number of participating nucleons. Statistical (systematic) uncertainties are displayed as error bars (boxes), while the global (fully correlated) uncertainties from the PbPb data (3.2%) or from the pp reference (6.3%) are displayed at unity as empty and filled red black boxes, respectively.

Nuclear modification factors for Υ(2S) meson production in PbPb collisions, as a function of centrality, displayed as the average number of participating nucleons. Statistical (systematic) uncertainties are displayed as error bars (boxes), while the global (fully correlated) uncertainties from the PbPb data (3.2%) or from the pp reference (6.9%) are displayed at unity as empty and filled black black boxes, respectively.

Nuclear modification factors for Υ(2S) meson production in PbPb collisions, as a function of centrality, displayed as the average number of participating nucleons. Statistical (systematic) uncertainties are displayed as error bars (boxes), while the global (fully correlated) uncertainties from the PbPb data (3.2%) or from the pp reference (6.9%) are displayed at unity as empty and filled black black boxes, respectively.

Nuclear modification factors for Υ(1S), Υ(2S) and Υ(3S) meson production in PbPb collisions, integrated over centrality. For the Υ(3S), the upper limit derived on the nuclear modification factor is given.

Nuclear modification factors for Υ(1S), Υ(2S) and Υ(3S) meson production in PbPb collisions, integrated over centrality. For the Υ(3S), the upper limit derived on the nuclear modification factor is given.

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.

Covariance matrix of total experimental uncertainties (statistical and systematic uncertainties added in quadrature) of absolute double differential fiducial cross section. The bin index is PT_i + 10*y_j.

Measured differential cross section dsigma/dET in pb/GeV given q_T < 10 GeV. For the data values, the first uncertainty is statistical and the second is systematic. For the theory values, the uncertainty combines statistical, PDF, and renormalization/factorization scale components.

Measured differential cross section dsigma/dET in pb/GeV given q_T > 50 GeV. For the data values, the first uncertainty is statistical and the second is systematic. For the theory values, the uncertainty combines statistical, PDF, and renormalization/factorization scale components.

Measured differential cross section dsigma/dDR in pb/radians. For the data values, the first uncertainty is statistical and the second is systematic. For the theory values, the uncertainty combines statistical, PDF, and renormalization/factorization scale components.

Measured differential cross section dsigma/dDR in pb/radians given q_T < 10 GeV. For the data values, the first uncertainty is statistical and the second is systematic. For the theory values, the uncertainty combines statistical, PDF, and renormalization/factorization scale components.

Measured differential cross section dsigma/dDR in pb/radians given q_T > 50 GeV. For the data values, the first uncertainty is statistical and the second is systematic. For the theory values, the uncertainty combines statistical, PDF, and renormalization/factorization scale components.

Fiducial cross sections for pT(gamma) > 15 GeV with statistical, systematic, and luminosity uncertainties, respectively. Combined result of electron and muon channels.

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.

The $p_{\rm T}$ bin width and corrected yield ratios R21 and R31 for $|y| < 0.6$. The statistical and systematic uncertainties in the corrected yield ratios are given as the percentage of the ratio.

The $p_{\rm T}$ bin width and corrected yield ratios R21 and R31 for $0.6 < |y| < 1.2$. The statistical and systematic uncertainties in the corrected yield ratios are given as the percentage of the ratio.

The $p_{\rm T}$ bin width and corrected yield ratios R21 and R31 for $|y| < 1.2$. The statistical and systematic uncertainties in the corrected yield ratios are given as the percentage of the ratio.

The dimuon acceptance $\mathcal{A}$ calculated using the CMS measured polarization and its positive and negative uncertainties $\mathcal{A}(\sigma^{\pm})$, and the values of the acceptance assuming no polarization $\mathcal{A}$(unpol), transverse polarization $\mathcal{A}(T)$, and longitudinal polarization $\mathcal{A}(L)$ for the $\Upsilon$(1S) state and $|y| <0.6$.

The dimuon acceptance $\mathcal{A}$ calculated using the CMS measured polarization and its positive and negative uncertainties $\mathcal{A}(\sigma^{\pm})$, and the values of the acceptance assuming no polarization $\mathcal{A}$(unpol), transverse polarization $\mathcal{A}(T)$, and longitudinal polarization $\mathcal{A}(L)$ for the $\Upsilon $(1S) state and $0.6 < |y| < 1.2$.

The dimuon acceptance $\mathcal{A}$ calculated using the CMS measured polarization and its positive and negative uncertainties $\mathcal{A}(\sigma^{\pm})$, and the values of the acceptance assuming no polarization $\mathcal{A}$(unpol), transverse polarization $\mathcal{A}(T)$, and longitudinal polarization $\mathcal{A}(L)$ for the $\Upsilon $(1S) state and $|y| < 1.2$.

The dimuon acceptance $\mathcal{A}$ calculated using the CMS measured polarization and its positive and negative uncertainties $\mathcal{A}(\sigma^{\pm})$, and the values of the acceptance assuming no polarization $\mathcal{A}$(unpol), transverse polarization $\mathcal{A}(T)$, and longitudinal polarization $\mathcal{A}(L)$ for the $\Upsilon $(2S) state and $|y| < 0.6$.

The dimuon acceptance $\mathcal{A}$ calculated using the CMS measured polarization and its positive and negative uncertainties $\mathcal{A}(\sigma^{\pm})$, and the values of the acceptance assuming no polarization $\mathcal{A}$(unpol), transverse polarization $\mathcal{A}(T)$, and longitudinal polarization $\mathcal{A}(L)$ for the $\Upsilon $(2S) state and $0.6 < |y| < 1.2$.

The dimuon acceptance $\mathcal{A}$ calculated using the CMS measured polarization and its positive and negative uncertainties $\mathcal{A}(\sigma^{\pm})$, and the values of the acceptance assuming no polarization $\mathcal{A}$(unpol), transverse polarization $\mathcal{A}(T)$, and longitudinal polarization $\mathcal{A}(L)$ for the $\Upsilon $(2S) state and $|y| < 1.2$.

The dimuon acceptance $\mathcal{A}$ calculated using the CMS measured polarization and its positive and negative uncertainties $\mathcal{A}(\sigma^{\pm})$, and the values of the acceptance assuming no polarization $\mathcal{A}$(unpol), transverse polarization $\mathcal{A}(T)$, and longitudinal polarization $\mathcal{A}(L)$ for the $\Upsilon $(3S) state and $|y| < 0.6$.

The dimuon acceptance $\mathcal{A}$ calculated using the CMS measured polarization and its positive and negative uncertainties $\mathcal{A}(\sigma^{\pm})$, and the values of the acceptance assuming no polarization $\mathcal{A}$(unpol), transverse polarization $\mathcal{A}(T)$, and longitudinal polarization $\mathcal{A}(L)$ for the $\Upsilon $(3S) state and $0.6 < |y| < 1.2$.

The dimuon acceptance $\mathcal{A}$ calculated using the CMS measured polarization and its positive and negative uncertainties $\mathcal{A}(\sigma^{\pm})$, and the values of the acceptance assuming no polarization $\mathcal{A}$(unpol), transverse polarization $\mathcal{A}(T)$, and longitudinal polarization $\mathcal{A}(L)$ for the $\Upsilon $(3S) state and $|y| < 1.2$.

The Drell-Yan pre-FSR dilepton rapidity distribution D(SIG)/DABS(YRAP) within the detector acceptance, for the mass bin 30-45 GeV, as measured in the combined dilepton channel.

The Drell-Yan pre-FSR dilepton rapidity distribution D(SIG)/DABS(YRAP) within the detector acceptance, for the mass bin 45-60 GeV, as measured in the combined dilepton channel.

The Drell-Yan pre-FSR dilepton rapidity distribution D(SIG)/DABS(YRAP) within the detector acceptance, for the mass bin 60-120 GeV, as measured in the combined dilepton channel.

The Drell-Yan pre-FSR dilepton rapidity distribution D(SIG)/DABS(YRAP) within the detector acceptance, for the mass bin 120-200 GeV, as measured in the combined dilepton channel.

The Drell-Yan pre-FSR dilepton rapidity distribution D(SIG)/DABS(YRAP) within the detector acceptance, for the mass bin 200-1500 GeV, as measured in the combined dilepton channel.

Measured ratio of the Drell-Yan normalized differential cross sections at center-of-mass energies of 7 and 8 TeV in the combined dilepton channel for the full phase space.

Measured ratio of the Drell-Yan normalized differential cross sections as a function of the absolute dilepton rapidity within the detector acceptance, at center-of-mass energies of 7 and 8 TeV in the combined dilepton channel. The mass bin 20-30 GeV.

Measured ratio of the Drell-Yan normalized differential cross sections as a function of the absolute dilepton rapidity within the detector acceptance, at center-of-mass energies of 7 and 8 TeV in the combined dilepton channel. The mass bin 30-45 GeV.

Measured ratio of the Drell-Yan normalized differential cross sections as a function of the absolute dilepton rapidity within the detector acceptance, at center-of-mass energies of 7 and 8 TeV in the combined dilepton channel. The mass bin 45-60 GeV.

Measured ratio of the Drell-Yan normalized differential cross sections as a function of the absolute dilepton rapidity within the detector acceptance, at center-of-mass energies of 7 and 8 TeV in the combined dilepton channel. The mass bin 60-120 GeV.

Measured ratio of the Drell-Yan normalized differential cross sections as a function of the absolute dilepton rapidity within the detector acceptance, at center-of-mass energies of 7 and 8 TeV in the combined dilepton channel. The mass bin 120-200 GeV.

Measured ratio of the Drell-Yan normalized differential cross sections as a function of the absolute dilepton rapidity within the detector acceptance, at center-of-mass energies of 7 and 8 TeV in the combined dilepton channel. The mass bin 200-1500 GeV.

Table of results for the dimuon channel for $1<|y|<2$ region $30 < M_{\ell\ell} < 60$ GeV. The first quoted uncertainty is statistical and the second is the total experimental systematic uncertainty.

Table of results for the dimuon channel for $160 < M_{\ell\ell} < 300$ GeV region. The first quoted uncertainty is statistical and the second is the total experimental systematic uncertainty.

Table of results for the dimuon channel for $300 < M_{\ell\ell} < 500$ GeV region. The first quoted uncertainty is statistical and the second is the total experimental systematic uncertainty.

The measured Z boson production cross section in pp collisions as a function of the Z boson rapidity for the dielectron decay channel.

Event centrality dependence of the $\mbox{Z}\rightarrow\mu^{+}\mu^{-}$ yields per MB event in PbPb collisions, divided by the expected average nuclear overlap function, $T_{\rm AA}$. On the horizontal axis, event centrality is depicted as the average number of participating nucleons, $N_{\rm part}$.

Event centrality dependence of the $\mbox{Z}\rightarrow e^{+}e^{-}$ yields per MB event in PbPb collisions, divided by the expected average nuclear overlap function, $T_{\rm AA}$. On the horizontal axis, event centrality is depicted as the average number of participating nucleons, $N_{\rm part}$.

The measured Z boson yields per MB event in PbPb collisions as a function of the Z boson $p_{\rm T}$ for the dimuon decay channel.

The measured Z boson yields per MB event in PbPb collisions as a function of the Z boson $p_{\rm T}$ for the dielectron decay channel.

The measured Z boson yields per MB event in PbPb collisions as a function of the Z boson rapidity for the dimuon decay channel.

The measured Z boson yields per MB event in PbPb collisions as a function of the Z boson rapidity for the dielectron decay channel.

The measured Z boson yields per MB event in PbPb collisions as a function of the Z boson rapidity for the combined leptonic decay channel.

The measured Z boson yields per MB event in PbPb collisions as a function of the Z boson $p_{\rm T}$ for the combined leptonic decay channel.

Event centrality dependence of the $\mbox{Z}\rightarrow l^{+}l^{-}$ yields per MB event in PbPb collisions, divided by the expected average nuclear overlap function, $T_{\rm AA}$. On the horizontal axis, event centrality is depicted as the average number of participating nucleons, $N_{\rm part}$.

The measured Z boson production cross section in pp collisions as a function of the Z boson rapidity for the combined leptonic decay channel.

The measured Z boson production cross section in pp collisions as a function of the Z boson pT for the combined leptonic decay channel in |y|<1.44.

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

The differential cross section measurement as a function of the transverse momentum of the second leading jet.

The differential cross section measurement as a function of the transverse momentum of the third leading jet.

The differential cross section measurement as a function of the transverse momentum of the fourth leading jet.

The differential cross section measurement as a function of the scalar sum of the transverse momenta of jets (HT) for jet multiplicities of one or more.

The differential cross section measurement as a function of the scalar sum of the transverse momenta of jets (HT) for jet multiplicities of two or more.

The differential cross section measurement as a function of the scalar sum of the transverse momenta of jets (HT) for jet multiplicities of three or more.

The differential cross section measurement as a function of the scalar sum of the transverse momenta of jets (HT) for jet multiplicities of four or more.

The differential cross section measurement as a function of the absolute value of the pseudorapidity of the first leading jet.

The differential cross section measurement as a function of the absolute value of the pseudorapidity of the first second jet.

The differential cross section measurement as a function of the absolute value of the pseudorapidity of the third leading jet.

The differential cross section measurement as a function of the absolute value of the pseudorapidity of the fourth leading jet.

The differential cross section measurement as a function of the absolute value of the azimuthal separation between the first leading jet and the muon.

The differential cross section measurement as a function of the absolute value of the azimuthal separation between the second leading jet and the muon.

The differential cross section measurement as a function of the absolute value of the azimuthal separation between the third leading jet and the muon.

The differential cross section measurement as a function of the absolute value of the azimuthal separation between the fourth leading jet and the muon.

Normalised differential TOP TOPBAR production cross section as a function of the number of additional jets with PT(JET) > 30 GeV in the lepton+jets channel. The total uncertainties are shown and include the statistical and systematic uncertainties.

Measured gap fraction in the dilepton channel as a function of the additional jet transverse momentum, GAPFRAC(PT) = N(PT)/N(TOTAL) where N(PT) is the number of events that do not contain additional jets above the given PT threshold.

Measured gap fraction in the dilepton channel as a function of the scalar sum of the transverse momentum of the additional jets, HT = SUM(PT(ADD. JETS)), GAPFRAC(HT) = N(HT)/N(TOTAL) where N(HT) is the number of events in which HT (with PT > 30 GeV) is less than the given threshold.

The product SIG( Zb+-(10610) PI-+ ) * BR( UPSILON(NS) PI-+ ) for N = 1, 2, 3.

The product SIG( Zb+-(10650) PI-+ ) * BR( UPSILON(NS) PI-+ ) for N = 1, 2, 3.

The Zb+-(10610) mass measured in each decay channel.

The Zb+-(10610) width measured in each decay channel.

The Zb+-(10650) mass measured in each decay channel.

The Zb+-(10650) width measured in each decay channel.

The fully corrected measurement of ASYM(FB) as a function of the muon-pair invariant mass.

The cross section ratio for a Z boson produced with at least one b-jet, relative to a Z boson produced with at least one jet of any flavour.

ASYMUU(COS(PHI(HADRON))) asymmetries for positive and negative hadrons as a function of XB. The errors are statistical and systematic.

ASYMUU(COS(PHI(HADRON))) asymmetries for positive and negative hadrons as a function of Z. The errors are statistical and systematic.

ASYMUU(COS(PHI(HADRON))) asymmetries for positive and negative hadrons as a function of PT(HADRON). The errors are statistical and systematic.

ASYMUU(COS(2*PHI(HADRON))) asymmetries for positive and negative hadrons as a function of XB. The errors are statistical and systematic.

ASYMUU(COS(2*PHI(HADRON))) asymmetries for positive and negative hadrons as a function of Z. The errors are statistical and systematic.

ASYMUU(COS(2*PHI(HADRON))) asymmetries for positive and negative hadrons as a function of PT(HADRON). The errors are statistical and systematic.

ASYMUU(COS(PHI(HADRON))) 3d asymmetries for positive and negative hadrons in the first XB bin. The errors are statistical and systematic.

ASYMUU(COS(PHI(HADRON))) 3d asymmetries for positive and negative hadrons in the second XB bin. The errors are statistical and systematic.

ASYMUU(COS(PHI(HADRON))) 3d asymmetries for positive and negative hadrons in the third XB bin. The errors are statistical and systematic.

ASYMUU(COS(PHI(HADRON))) 3d asymmetries for positive and negative hadrons in the last XB bin. The errors are statistical and systematic.

ASYMUU(COS(2*PHI(HADRON))) 3d asymmetries for positive and negative hadrons in the first XB bin. The errors are statistical and systematic.

ASYMUU(COS(2*PHI(HADRON))) 3d asymmetries for positive and negative hadrons in the second XB bin. The errors are statistical and systematic.

ASYMUU(COS(2*PHI(HADRON))) 3d asymmetries for positive and negative hadrons in the third XB bin. The errors are statistical and systematic.

ASYMUU(COS(2*PHI(HADRON))) 3d asymmetries for positive and negative hadrons in the last XB bin. The errors are statistical and systematic.

Covariance matrix (statistical and systematic uncertainties combined) with the muon $p_{T}>35$ GeV. The units are in $10^{-4}$.

The DY pre-FSR cross section measurements for the combined dimuon and dielectron channels normalized to the Z-peak region in the full acceptance.

The DY dimuon rapidity spectrum within the detector acceptance, normalized to the Z-peak region, for the mass bin 20-30 GeV. The two columns correspond to pre-FSR and post-FSR. The uncertainties are the total experimental uncertainties.

The DY dimuon rapidity spectrum within the detector acceptance, normalized to the Z-peak region, for the mass bin 30-45 GeV. The two columns correspond to pre-FSR and post-FSR. The uncertainties are the total experimental uncertainties.

The DY dimuon rapidity spectrum within the detector acceptance, normalized to the Z-peak region, for the mass bin 45-60 GeV. The two columns correspond to pre-FSR and post-FSR. The uncertainties are the total experimental uncertainties.

The DY dimuon rapidity spectrum within the detector acceptance, normalized to the Z-peak region, for the mass bin 60-120 GeV. The two columns correspond to pre-FSR and post-FSR. The uncertainties are the total experimental uncertainties.

The DY dimuon rapidity spectrum within the detector acceptance, normalized to the Z-peak region, for the mass bin 120-200 GeV. The two columns correspond to pre-FSR and post-FSR. The uncertainties are the total experimental uncertainties.

The DY dimuon rapidity spectrum within the detector acceptance, normalized to the Z-peak region, for the mass bin 200-1500 GeV. The two columns correspond to pre-FSR and post-FSR. The uncertainties are the total experimental uncertainties.

Contributions from individual sources of systematic uncertainty for the ($p^{\mu}_{T} > 35$, $E_T^{missing} > 35$) GeV kinematic region. All uncertainty values are multiplied by 100. The columns (1-7) correspond to: 1.0 = Electro-Weak background 2.0 = Multi-Jet background 3.0 = Charge mis-identification 4.0 = Relative charge efficiency 5.0 = Magnet polarity weighting 6.0 = Momentum/$E_T^{missing}$ resolution 7.0 = Trigger isolation.

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.

PT dependences of the differential multiplicities for 0.0080 < x_Bjorken < 0.0120 and 1.00 < Q^2 < 1.50 GeV^2 for Positive hadrons.

PT dependences of the differential multiplicities for 0.0080 < x_Bjorken < 0.0120 and 1.50 < Q^2 < 2.10 GeV^2 for Positive hadrons.

PT dependences of the differential multiplicities for 0.0120 < x_Bjorken < 0.0180 and 1.00 < Q^2 < 1.50 GeV^2 for Positive hadrons.

PT dependences of the differential multiplicities for 0.0120 < x_Bjorken < 0.0180 and 1.50 < Q^2 < 2.50 GeV^2 for Positive hadrons.

PT dependences of the differential multiplicities for 0.0120 < x_Bjorken < 0.0180 and 2.50 < Q^2 < 3.50 GeV^2 for Positive hadrons.

PT dependences of the differential multiplicities for 0.0180 < x_Bjorken < 0.0250 and 1.00 < Q^2 < 1.50 GeV^2 for Positive hadrons.

PT dependences of the differential multiplicities for 0.0180 < x_Bjorken < 0.0250 and 1.50 < Q^2 < 2.50 GeV^2 for Positive hadrons.

PT dependences of the differential multiplicities for 0.0180 < x_Bjorken < 0.0250 and 2.50 < Q^2 < 3.50 GeV^2 for Positive hadrons.

PT dependences of the differential multiplicities for 0.0180 < x_Bjorken < 0.0250 and 3.50 < Q^2 < 5.00 GeV^2 for Positive hadrons.

PT dependences of the differential multiplicities for 0.0250 < x_Bjorken < 0.0350 and 1.00 < Q^2 < 1.20 GeV^2 for Positive hadrons.

PT dependences of the differential multiplicities for 0.0250 < x_Bjorken < 0.0400 and 1.20 < Q^2 < 1.50 GeV^2 for Positive hadrons.

PT dependences of the differential multiplicities for 0.0250 < x_Bjorken < 0.0400 and 1.50 < Q^2 < 2.50 GeV^2 for Positive hadrons.

PT dependences of the differential multiplicities for 0.0250 < x_Bjorken < 0.0400 and 2.50 < Q^2 < 3.50 GeV^2 for Positive hadrons.

PT dependences of the differential multiplicities for 0.0250 < x_Bjorken < 0.0400 and 3.50 < Q^2 < 6.00 GeV^2 for Positive hadrons.

PT dependences of the differential multiplicities for 0.0400 < x_Bjorken < 0.0500 and 1.50 < Q^2 < 2.50 GeV^2 for Positive hadrons.

PT dependences of the differential multiplicities for 0.0400 < x_Bjorken < 0.0700 and 2.50 < Q^2 < 3.50 GeV^2 for Positive hadrons.

PT dependences of the differential multiplicities for 0.0400 < x_Bjorken < 0.0700 and 3.50 < Q^2 < 6.00 GeV^2 for Positive hadrons.

PT dependences of the differential multiplicities for 0.0400 < x_Bjorken < 0.0700 and 6.00 < Q^2 < 10.00 GeV^2 for Positive hadrons.

PT dependences of the differential multiplicities for 0.0700 < x_Bjorken < 0.1200 and 3.50 < Q^2 < 6.00 GeV^2 for Positive hadrons.

PT dependences of the differential multiplicities for 0.0700 < x_Bjorken < 0.1200 and 6.00 < Q^2 < 10.00 GeV^2 for Positive hadrons.

PT dependences of the differential multiplicities for 0.0045 < x_Bjorken < 0.0060 and 1.00 < Q^2 < 1.25 GeV^2 for Negative hadrons.

PT dependences of the differential multiplicities for 0.0060 < x_Bjorken < 0.0080 and 1.00 < Q^2 < 1.30 GeV^2 for Negative hadrons.

PT dependences of the differential multiplicities for 0.0060 < x_Bjorken < 0.0080 and 1.30 < Q^2 < 1.70 GeV^2 for Negative hadrons.

PT dependences of the differential multiplicities for 0.0080 < x_Bjorken < 0.0120 and 1.00 < Q^2 < 1.50 GeV^2 for Negative hadrons.

PT dependences of the differential multiplicities for 0.0080 < x_Bjorken < 0.0120 and 1.50 < Q^2 < 2.10 GeV^2 for Negative hadrons.

PT dependences of the differential multiplicities for 0.0120 < x_Bjorken < 0.0180 and 1.00 < Q^2 < 1.50 GeV^2 for Negative hadrons.

PT dependences of the differential multiplicities for 0.0120 < x_Bjorken < 0.0180 and 1.50 < Q^2 < 2.50 GeV^2 for Negative hadrons.

PT dependences of the differential multiplicities for 0.0120 < x_Bjorken < 0.0180 and 2.50 < Q^2 < 3.50 GeV^2 for Negative hadrons.

PT dependences of the differential multiplicities for 0.0180 < x_Bjorken < 0.0250 and 1.00 < Q^2 < 1.50 GeV^2 for Negative hadrons.

PT dependences of the differential multiplicities for 0.0180 < x_Bjorken < 0.0250 and 1.50 < Q^2 < 2.50 GeV^2 for Negative hadrons.

PT dependences of the differential multiplicities for 0.0180 < x_Bjorken < 0.0250 and 2.50 < Q^2 < 3.50 GeV^2 for Negative hadrons.

PT dependences of the differential multiplicities for 0.0180 < x_Bjorken < 0.0250 and 3.50 < Q^2 < 5.00 GeV^2 for Negative hadrons.

PT dependences of the differential multiplicities for 0.0250 < x_Bjorken < 0.0350 and 1.00 < Q^2 < 1.20 GeV^2 for Negative hadrons.

PT dependences of the differential multiplicities for 0.0250 < x_Bjorken < 0.0400 and 1.20 < Q^2 < 1.50 GeV^2 for Negative hadrons.

PT dependences of the differential multiplicities for 0.0250 < x_Bjorken < 0.0400 and 1.50 < Q^2 < 2.50 GeV^2 for Negative hadrons.

PT dependences of the differential multiplicities for 0.0250 < x_Bjorken < 0.0400 and 2.50 < Q^2 < 3.50 GeV^2 for Negative hadrons.

PT dependences of the differential multiplicities for 0.0250 < x_Bjorken < 0.0400 and 3.50 < Q^2 < 6.00 GeV^2 for Negative hadrons.

PT dependences of the differential multiplicities for 0.0400 < x_Bjorken < 0.0500 and 1.50 < Q^2 < 2.50 GeV^2 for Negative hadrons.

PT dependences of the differential multiplicities for 0.0400 < x_Bjorken < 0.0700 and 2.50 < Q^2 < 3.50 GeV^2 for Negative hadrons.

PT dependences of the differential multiplicities for 0.0400 < x_Bjorken < 0.0700 and 3.50 < Q^2 < 6.00 GeV^2 for Negative hadrons.

PT dependences of the differential multiplicities for 0.0400 < x_Bjorken < 0.0700 and 6.00 < Q^2 < 10.00 GeV^2 for Negative hadrons.

PT dependences of the differential multiplicities for 0.0700 < x_Bjorken < 0.1200 and 3.50 < Q^2 < 6.00 GeV^2 for Negative hadrons.

PT dependences of the differential multiplicities for 0.0700 < x_Bjorken < 0.1200 and 6.00 < Q^2 < 10.00 GeV^2 for Negative hadrons.

Fitted average PTs for positive hadrons.

Fitted average PTs for negative hadrons.

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

The fiducial and acceptance corrected UPSI(3S) production cross sections (times di-muon branching ratio) as a function of PT for the |rapidity| range < 2.4. Note these are integrated cross sections and the acceptance-corrected cross sections assume the UPSI(3S) are unpolarized with the variations due to the 4 extreme polarization scenarios shown in the last 4 columns. The fiducial cross sections do not need to make any assumptions on the polarizations scenarios. The luminosity uncertainty of 4% is not included in the systematic errors.

The fiducial and acceptance corrected UPSI(1S) production cross sections (times di-muon branching ratio) as a function of PT for the |rapidity| range < 0.4. Note these are integrated cross sections and the acceptance-corrected cross sections assume the UPSI(1S) are unpolarized with the variations due to the 4 extreme polarization scenarios shown in the last 4 columns. The fiducial cross sections do not need to make any assumptions on the polarizations scenarios. The luminosity uncertainty of 4% is not included in the systematic errors.

The fiducial and acceptance corrected UPSI(2S) production cross sections (times di-muon branching ratio) as a function of PT for the |rapidity| range < 0.4. Note these are integrated cross sections and the acceptance-corrected cross sections assume the UPSI(2S) are unpolarized with the variations due to the 4 extreme polarization scenarios shown in the last 4 columns. The fiducial cross sections do not need to make any assumptions on the polarizations scenarios. The luminosity uncertainty of 4% is not included in the systematic errors.

The fiducial and acceptance corrected UPSI(3S) production cross sections (times di-muon branching ratio) as a function of PT for the |rapidity| range < 0.4. Note these are integrated cross sections and the acceptance-corrected cross sections assume the UPSI(3S) are unpolarized with the variations due to the 4 extreme polarization scenarios shown in the last 4 columns. The fiducial cross sections do not need to make any assumptions on the polarizations scenarios. The luminosity uncertainty of 4% is not included in the systematic errors.

The fiducial and acceptance corrected UPSI(1S) production cross sections (times di-muon branching ratio) as a function of PT for the |rapidity| range 0.4 to 0.8. Note these are integrated cross sections and the acceptance-corrected cross sections assume the UPSI(1S) are unpolarized with the variations due to the 4 extreme polarization scenarios shown in the last 4 columns. The fiducial cross sections do not need to make any assumptions on the polarizations scenarios. The luminosity uncertainty of 4% is not included in the systematic errors.

The fiducial and acceptance corrected UPSI(2S) production cross sections (times di-muon branching ratio) as a function of PT for the |rapidity| range 0.4 to 0.8. Note these are integrated cross sections and the acceptance-corrected cross sections assume the UPSI(2S) are unpolarized with the variations due to the 4 extreme polarization scenarios shown in the last 4 columns. The fiducial cross sections do not need to make any assumptions on the polarizations scenarios. The luminosity uncertainty of 4% is not included in the systematic errors.

The fiducial and acceptance corrected UPSI(3S) production cross sections (times di-muon branching ratio) as a function of PT for the |rapidity| range 0.4 to 0.8. Note these are integrated cross sections and the acceptance-corrected cross sections assume the UPSI(3S) are unpolarized with the variations due to the 4 extreme polarization scenarios shown in the last 4 columns. The fiducial cross sections do not need to make any assumptions on the polarizations scenarios. The luminosity uncertainty of 4% is not included in the systematic errors.

The fiducial and acceptance corrected UPSI(1S) production cross sections (times di-muon branching ratio) as a function of PT for the |rapidity| range 0.8 to 1.2. Note these are integrated cross sections and the acceptance-corrected cross sections assume the UPSI(1S) are unpolarized with the variations due to the 4 extreme polarization scenarios shown in the last 4 columns. The fiducial cross sections do not need to make any assumptions on the polarizations scenarios. The luminosity uncertainty of 4% is not included in the systematic errors.

The fiducial and acceptance corrected UPSI(2S) production cross sections (times di-muon branching ratio) as a function of PT for the |rapidity| range 0.8 to 1.2. Note these are integrated cross sections and the acceptance-corrected cross sections assume the UPSI(2S) are unpolarized with the variations due to the 4 extreme polarization scenarios shown in the last 4 columns. The fiducial cross sections do not need to make any assumptions on the polarizations scenarios. The luminosity uncertainty of 4% is not included in the systematic errors.

The fiducial and acceptance corrected UPSI(3S) production cross sections (times di-muon branching ratio) as a function of PT for the |rapidity| range 0.8 to 1.2. Note these are integrated cross sections and the acceptance-corrected cross sections assume the UPSI(3S) are unpolarized with the variations due to the 4 extreme polarization scenarios shown in the last 4 columns. The fiducial cross sections do not need to make any assumptions on the polarizations scenarios. The luminosity uncertainty of 4% is not included in the systematic errors.

The fiducial and acceptance corrected UPSI(1S) production cross sections (times di-muon branching ratio) as a function of PT for the |rapidity| range 1.2 to 1.6. Note these are integrated cross sections and the acceptance-corrected cross sections assume the UPSI(1S) are unpolarized with the variations due to the 4 extreme polarization scenarios shown in the last 4 columns. The fiducial cross sections do not need to make any assumptions on the polarizations scenarios. The luminosity uncertainty of 4% is not included in the systematic errors.

The fiducial and acceptance corrected UPSI(2S) production cross sections (times di-muon branching ratio) as a function of PT for the |rapidity| range 1.2 to 1.6. Note these are integrated cross sections and the acceptance-corrected cross sections assume the UPSI(2S) are unpolarized with the variations due to the 4 extreme polarization scenarios shown in the last 4 columns. The fiducial cross sections do not need to make any assumptions on the polarizations scenarios. The luminosity uncertainty of 4% is not included in the systematic errors.

The fiducial and acceptance corrected UPSI(3S) production cross sections (times di-muon branching ratio) as a function of PT for the |rapidity| range 1.2 to 1.6. Note these are integrated cross sections and the acceptance-corrected cross sections assume the UPSI(3S) are unpolarized with the variations due to the 4 extreme polarization scenarios shown in the last 4 columns. The fiducial cross sections do not need to make any assumptions on the polarizations scenarios. The luminosity uncertainty of 4% is not included in the systematic errors.

The fiducial and acceptance corrected UPSI(1S) production cross sections (times di-muon branching ratio) as a function of PT for the |rapidity| range 1.6 to 2.0. Note these are integrated cross sections and the acceptance-corrected cross sections assume the UPSI(1S) are unpolarized with the variations due to the 4 extreme polarization scenarios shown in the last 4 columns. The fiducial cross sections do not need to make any assumptions on the polarizations scenarios. The luminosity uncertainty of 4% is not included in the systematic errors.

The fiducial and acceptance corrected UPSI(2S) production cross sections (times di-muon branching ratio) as a function of PT for the |rapidity| range 1.6 to 2.0. Note these are integrated cross sections and the acceptance-corrected cross sections assume the UPSI(2S) are unpolarized with the variations due to the 4 extreme polarization scenarios shown in the last 4 columns. The fiducial cross sections do not need to make any assumptions on the polarizations scenarios. The luminosity uncertainty of 4% is not included in the systematic errors.

The fiducial and acceptance corrected UPSI(3S) production cross sections (times di-muon branching ratio) as a function of PT for the |rapidity| range 1.6 to 2.0. Note these are integrated cross sections and the acceptance-corrected cross sections assume the UPSI(3S) are unpolarized with the variations due to the 4 extreme polarization scenarios shown in the last 4 columns. The fiducial cross sections do not need to make any assumptions on the polarizations scenarios. The luminosity uncertainty of 4% is not included in the systematic errors.

The fiducial and acceptance corrected UPSI(1S) production cross sections (times di-muon branching ratio) as a function of PT for the |rapidity| range 2.0 to 2.4. Note these are integrated cross sections and the acceptance-corrected cross sections assume the UPSI(1S) are unpolarized with the variations due to the 4 extreme polarization scenarios shown in the last 4 columns. The fiducial cross sections do not need to make any assumptions on the polarizations scenarios. The luminosity uncertainty of 4% is not included in the systematic errors.

The fiducial and acceptance corrected UPSI(2S) production cross sections (times di-muon branching ratio) as a function of PT for the |rapidity| range 2.0 to 2.4. Note these are integrated cross sections and the acceptance-corrected cross sections assume the UPSI(2S) are unpolarized with the variations due to the 4 extreme polarization scenarios shown in the last 4 columns. The fiducial cross sections do not need to make any assumptions on the polarizations scenarios. The luminosity uncertainty of 4% is not included in the systematic errors.

The fiducial and acceptance corrected UPSI(3S) production cross sections (times di-muon branching ratio) as a function of PT for the |rapidity| range 2.0 to 2.4. Note these are integrated cross sections and the acceptance-corrected cross sections assume the UPSI(3S) are unpolarized with the variations due to the 4 extreme polarization scenarios shown in the last 4 columns. The fiducial cross sections do not need to make any assumptions on the polarizations scenarios. The luminosity uncertainty of 4% is not included in the systematic errors.

The fiducial and acceptance corrected UPSI(1S) production cross sections (times di-muon branching ratio) as a function of rapidity for PT < 50 GeV/c. Note these are integrated cross sections and the acceptance-corrected cross sections assume the UPSI(1S) are unpolarized with the variations due to the 4 extreme polarization scenarios shown in the last 4 columns. The fiducial cross sections do not need to make any assumptions on the polarizations scenarios. The luminosity uncertainty of 4% is not included in the systematic errors.

The fiducial and acceptance corrected UPSI(2S) production cross sections (times di-muon branching ratio) as a function of rapidity for PT < 50 GeV/c. Note these are integrated cross sections and the acceptance-corrected cross sections assume the UPSI(2S) are unpolarized with the variations due to the 4 extreme polarization scenarios shown in the last 4 columns. The fiducial cross sections do not need to make any assumptions on the polarizations scenarios. The luminosity uncertainty of 4% is not included in the systematic errors.

The fiducial and acceptance corrected UPSI(3S) production cross sections (times di-muon branching ratio) as a function of rapidity for PT < 50 GeV/c. Note these are integrated cross sections and the acceptance-corrected cross sections assume the UPSI(3S) are unpolarized with the variations due to the 4 extreme polarization scenarios shown in the last 4 columns. The fiducial cross sections do not need to make any assumptions on the polarizations scenarios. The luminosity uncertainty of 4% is not included in the systematic errors.

The ratios of the UPSI(3S) to UPSI(1S) fiducial and acceptance cross sections as a function of PT for the |rapidiy| range < 2.4.

The ratios of the UPSI(2S) to UPSI(1S) fiducial and acceptance cross sections as a function of PT for the |rapidiy| range < 2.4.

The ratios of the UPSI(3S) to UPSI(2S) fiducial and acceptance cross sections as a function of PT for the |rapidiy| range < 2.4.

The acceptance-corrected UPSI(1S) production cross section, integrated and differential in PT, times the respective dimuon branching fraction, integrated over the rapidity range |y| < 1.2. The first (sys) error is the quadratic sum of all systematic uncertainties excluding that associated with the polarization of the UPSI(1S) which is shown separately as the second (sys) value. The latter has been obtained by calculating the acceptance with the +68 and -68% CL as given by the measured UPSI(1S) polarization [see. PRL 110 (2013) 081802].

The acceptance-corrected UPSI(2S) production cross section, integrated and differential in PT, times the respective dimuon branching fraction, integrated over the rapidity range |y| < 1.2. The first (sys) error is the quadratic sum of all systematic uncertainties excluding that associated with the polarization of the UPSI(2S) which is shown separately as the second (sys) value. The latter has been obtained by calculating the acceptance with the +68 and -68% CL as given by the measured UPSI(2S) polarization [see. PRL 110 (2013) 081802].

The acceptance-corrected UPSI(3S) production cross section, integrated and differential in PT, times the respective dimuon branching fraction, integrated over the rapidity range |y| < 1.2. The first (sys) error is the quadratic sum of all systematic uncertainties excluding that associated with the polarization of the UPSI(3S) which is shown separately as the second (sys) value. The latter has been obtained by calculating the acceptance with the +68 and -68% CL as given by the measured UPSI(3S) polarization [see. PRL 110 (2013) 081802].