The total coherent $\mathrm{J}/\psi$ photoproduction cross section as a function of photon-nuclear center-of-mass energy per nucleon $W_{\gamma \mathrm{N}}^{\mathrm{Pb}}$, measured in PbPb ultra-peripheral collisions at $\sqrt{s_{\mathrm{NN}}}$ = 5.02 TeV. The $W_{\gamma \mathrm{N}}^{\mathrm{Pb}}$ values used correspond to the center of each rapidity range. The theoretical uncertainties is due to the uncertainties in the photon flux.

The nuclear gluon suppression factor $R_{\mathrm{g}}^{\mathrm{Pb}}$ as a function of Bjorken $x$ extracted from the CMS measurement of the coherent $\mathrm{J}/\psi$ photoproduction in PbPb ultra-peripheral collisions at $\sqrt{s_{\mathrm{NN}}}$ = 5.02 TeV. The $x$ values are evaluated at the centers of their corresponding rapidity ranges. The theoretical uncertainties are due to the uncertainties in the photon flux and the impulse approximation model.

The nuclear gluon suppression factor $R_{\mathrm{g}}^{\mathrm{Pb}}$ as a function of Bjorken $x$ extracted from the CMS measurement of the coherent $\mathrm{J}/\psi$ photoproduction in PbPb ultra-peripheral collisions at $\sqrt{s_{\mathrm{NN}}}$ = 5.02 TeV. The $x$ values are evaluated at the centers of their corresponding rapidity ranges. The theoretical uncertainties are due to the uncertainties in the photon flux and the impulse approximation model.

The total covariance matrix of the total coherent photoproduction cross section as a function of photon-nuclear center-of-mass energy per nucleon $W_{\gamma \mathrm{N}}^{\mathrm{Pb}}$. The covariance matrix includes both the experimental and theoretical (photon flux) uncertainties. The bins are ordered as increasing in $W_{\gamma \mathrm{N}}^{\mathrm{Pb}}$.

The photon flux values from STARLight as a function of rapidity, in different neutron multiplicity classes: 0n0n, 0nXn, and XnXn.

The experimental covariance matrix of the total coherent photoproduction cross section as a function of photon-nuclear center-of-mass energy per nucleon $W_{\gamma \mathrm{N}}^{\mathrm{Pb}}$. The bins are ordered as increasing in $W_{\gamma \mathrm{N}}^{\mathrm{Pb}}$.

The covariance matrix for the flux in the 0n0n neutron multiplicity class.

The theoretical (photon flux) covariance matrix of the total coherent photoproduction cross section as a function of photon-nuclear center-of-mass energy per nucleon $W_{\gamma \mathrm{N}}^{\mathrm{Pb}}$. The bins are ordered as increasing in $W_{\gamma \mathrm{N}}^{\mathrm{Pb}}$.

The covariance matrix for the flux in the 0nXn neutron multiplicity class.

The covariance matrix for the flux in the XnXn neutron multiplicity class.

The total covariance matrix of the total coherent photoproduction cross section as a function of photon-nuclear center-of-mass energy per nucleon $W_{\gamma \mathrm{N}}^{\mathrm{Pb}}$. The covariance matrix includes both the experimental and theoretical (photon flux) uncertainties. The bins are ordered as increasing in $W_{\gamma \mathrm{N}}^{\mathrm{Pb}}$.

The experimental covariance matrix of the total coherent photoproduction cross section as a function of photon-nuclear center-of-mass energy per nucleon $W_{\gamma \mathrm{N}}^{\mathrm{Pb}}$. The bins are ordered as increasing in $W_{\gamma \mathrm{N}}^{\mathrm{Pb}}$.

The theoretical (photon flux) covariance matrix of the total coherent photoproduction cross section as a function of photon-nuclear center-of-mass energy per nucleon $W_{\gamma \mathrm{N}}^{\mathrm{Pb}}$. The bins are ordered as increasing in $W_{\gamma \mathrm{N}}^{\mathrm{Pb}}$.

$M_{p\eta}$ invariant mass distribution for $E_\gamma$ = 1400-1500 MeV.

$M_{p\eta}$ invariant mass distribution for $E_\gamma$ = 1500-1600 MeV.

$M_{p\eta}$ invariant mass distribution under the conditions: $E_\gamma$ = 1400-1500 MeV; 1120 MeV $\leq M_{p\pi^0} \leq$ 1220 MeV; $\theta_p^{lab} \leq$ 25°; 25° $\leq \theta_\gamma^{lab} \leq$ 155°.

$M_{p\eta}$ invariant mass distribution under the conditions: $E_\gamma$ = 1400-1500 MeV; 1120 MeV $\leq M_{p\pi^0} \leq$ 1220 MeV; $\theta_p^{lab} \leq$ 25°; 1° $\leq \theta_\gamma^{lab} \leq$ 155°.

$M_{p\eta}$ invariant mass distribution for $E_\gamma$ = 1400-1500 MeV and $M_{p\pi^0} \leq$ 1190 MeV, allowing for all proton and photon laboratory angles.

$M_{p\eta}$ invariant mass distribution for the incident photon energy range $E_\gamma$ = 1420-1540 MeV and $M_{p\pi^0} \leq$ 1170 MeV.

$M_{p\eta}$ invariant mass distribution for the incident photon energy range $E_\gamma$ = 1420-1540 MeV and $M_{p\pi^0} \leq$ 1180 MeV.

$M_{p\eta}$ invariant mass distribution for the incident photon energy range $E_\gamma$ = 1420-1540 MeV and $M_{p\pi^0} \leq$ 1190 MeV.

$M_{p\eta}$ invariant mass distribution for the incident photon energy range $E_\gamma$ = 1420-1540 MeV and $M_{p\pi^0} \leq$ 1200 MeV.

$M_{p\eta}$ invariant mass distribution for the incident photon energy range $E_\gamma$ = 1420-1540 MeV and $M_{p\pi^0} \leq$ 1210 MeV.

$M_{p\eta}$ invariant mass distribution for the incident photon energy range $E_\gamma$ = 1420-1540 MeV and $M_{p\pi^0} \leq$ 1220 MeV.

$M_{p\eta}$ invariant mass distribution for the incident photon energy range $E_\gamma$ = 1400-1420 MeV and $M_{p\pi^0} \leq$ 1190 MeV.

$M_{p\eta}$ invariant mass distribution for the incident photon energy range $E_\gamma$ = 1420-1460 MeV and $M_{p\pi^0} \leq$ 1190 MeV.

$M_{p\eta}$ invariant mass distribution for the incident photon energy range $E_\gamma$ = 1460-1500 MeV and $M_{p\pi^0} \leq$ 1190 MeV.

$M_{p\eta}$ invariant mass distribution for the incident photon energy range $E_\gamma$ = 1500-1540 MeV and $M_{p\pi^0} \leq$ 1190 MeV.

$M_{p\eta}$ invariant mass distribution for the incident photon energy range $E_\gamma$ = 1540-1600 MeV and $M_{p\pi^0} \leq$ 1190 MeV.

$M_{p\eta}$ invariant mass distribution for the incident photon energy range $E_\gamma$ = 1420-1540 MeV and $M_{p\pi^0} \leq$ 1190 MeV.

Excess yield relative to the PWA distributions in the invariant mass distributions $M^2_{p\eta}$.

Excess yield relative to the PWA distributions in the invariant mass distributions $M^2_{\pi^0\eta}$.

Excess yield relative to the PWA distributions in the invariant mass distributions $M^2_{p\pi^0}$, with a cut on the structure in the $M^2_{p\eta}$ invariant mass distribution.

Excess yield relative to the PWA distributions in the $p-\eta$ opening angle distribution in the $\gamma p$ centre-of-mass system.

Excess yield relative to the PWA distributions in the $p-\eta$ opening angle distribution in the $\gamma p$ centre-of-mass system.

Excess yield relative to the PWA distributions in the $\pi^0-\eta$ opening angle distribution in the $\gamma p$ centre-of-mass system.

Excess yield relative to the PWA distributions in the $\pi^0-\eta$ opening angle distribution in the $\gamma p$ centre-of-mass system.

Differential cross section $d\sigma/dM_{\pi^0\eta}$ for $E_\gamma$ = 1540-1600 MeV.

$M_{p\eta}$ invariant mass distribution, taken as difference to the PWA calculations (s. Fig. 11).

$M_{\pi^0\eta}$ invariant mass distribution, taken as difference to the PWA calculations (s. Fig. 11).

$\theta_{p\eta}$ opening angle distribution, taken as difference to the PWA calculations (s. Fig. 12).

$\theta_{\pi^0\eta}$ opening angle distribution, taken as difference to the PWA calculations (s. Fig. 12).

Differential fiducial cross section (without the acceptance correction) for the DY process measured in the muon channel, as a function of $p_{\textrm{T}}$ for $15<m_{\mu\mu}<60$ GeV. The quoted error is the quadratic sum of the statistical and systematic uncertainties. The global normalisation uncertainty of 3.5% is listed separately.

Differential fiducial cross section (without the acceptance correction) for the DY process measured in the muon channel, as a function of $p_{\textrm{T}}$ for $60<m_{\mu\mu}<120$ GeV. The quoted error is the quadratic sum of the statistical and systematic uncertainties. The global normalisation uncertainty of 3.5% is listed separately.

Differential fiducial cross section (without the acceptance correction) for the DY process measured in the muon channel, as a function of $\phi^*$ for $15<m_{\mu\mu}<60$ GeV. The quoted error is the quadratic sum of the statistical and systematic uncertainties. The global normalisation uncertainty of 3.5% is listed separately.

Differential fiducial cross section (without the acceptance correction) for the DY process measured in the muon channel, as a function of $\phi^*$ for $60<m_{\mu\mu}<120$ GeV. The quoted error is the quadratic sum of the statistical and systematic uncertainties. The global normalisation uncertainty of 3.5% is listed separately.

Integrated fiducial cross section for the DY process measured in the muon channel. The quoted error is the quadratic sum of the statistical and systematic uncertainties. The global normalisation uncertainty of 3.5% is listed separately.

Differential cross section for the DY process measured in the muon channel, as a function of dimuon invariant mass. The quoted error is the quadratic sum of the statistical and systematic uncertainties. The global normalisation uncertainty of 3.5% is listed separately.

Differential cross section for the DY process measured in the muon channel, as a function of rapidity in the centre-of-mass frame for $15<m_{\mu\mu}<60$ GeV. The quoted error is the quadratic sum of the statistical and systematic uncertainties. The global normalisation uncertainty of 3.5% is listed separately.

Differential cross section for the DY process measured in the muon channel, as a function of rapidity in the centre-of-mass frame for $60<m_{\mu\mu}<120$ GeV. The quoted error is the quadratic sum of the statistical and systematic uncertainties. The global normalisation uncertainty of 3.5% is listed separately.

Differential cross section for the DY process measured in the muon channel, as a function of $p_{\textrm{T}}$ for $15<m_{\mu\mu}<60$ GeV. The quoted error is the quadratic sum of the statistical and systematic uncertainties. The global normalisation uncertainty of 3.5% is listed separately.

Differential cross section for the DY process measured in the muon channel, as a function of $p_{\textrm{T}}$ for $60<m_{\mu\mu}<120$ GeV. The quoted error is the quadratic sum of the statistical and systematic uncertainties. The global normalisation uncertainty of 3.5% is listed separately.

Differential cross section for the DY process measured in the muon channel, as a function of $\phi^*$ for $15<m_{\mu\mu}<60$ GeV. The quoted error is the quadratic sum of the statistical and systematic uncertainties. The global normalisation uncertainty of 3.5% is listed separately.

Differential cross section for the DY process measured in the muon channel, as a function of $\phi^*$ for $60<m_{\mu\mu}<120$ GeV. The quoted error is the quadratic sum of the statistical and systematic uncertainties. The global normalisation uncertainty of 3.5% is listed separately.

Integrated cross section for the DY process measured in the muon channel. The quoted error is the quadratic sum of the statistical and systematic uncertainties. The global normalisation uncertainty of 3.5% is listed separately.

Forward-backward ratios for $15<m_{\mu\mu}<60$ GeV. Quoted uncertainties are the quadratic sum of the statistical and systematic uncertainties.

Forward-backward ratios for $60<m_{\mu\mu}<120$ GeV. Quoted uncertainties are the quadratic sum of the statistical and systematic uncertainties.

Correlation matrix for the systematic uncertainties, excluding integrated luminosity, as a function of dimuon invariant mass.

Correlation matrix for the systematic uncertainties, excluding integrated luminosity, as a function of rapidity in the centre-of-mass frame for $15<m_{\mu\mu}<60$ GeV.

Correlation matrix for the systematic uncertainties, excluding integrated luminosity, as a function of rapidity in the centre-of-mass frame for $60<m_{\mu\mu}<120$ GeV.

Correlation matrix for the systematic uncertainties, excluding integrated luminosity, as a function of $p_{\textrm{T}}$ for $15<m_{\mu\mu}<60$ GeV.

Correlation matrix for the systematic uncertainties, excluding integrated luminosity, as a function of $p_{\textrm{T}}$ for $60<m_{\mu\mu}<120$ GeV.

Correlation matrix for the systematic uncertainties, excluding integrated luminosity, as a function of $\phi^*$ for $15<m_{\mu\mu}<60$ GeV.

Correlation matrix for the systematic uncertainties, excluding integrated luminosity, as a function of $\phi^*$ for $60<m_{\mu\mu}<120$ GeV.

Correlation matrix for the systematic uncertainties, excluding integrated luminosity, as a function of rapidity in the centre-of-mass frame, putting together low and high masses. Each bin is labeled lm_i or hm_i, where lm stands for $15<m_{\mu\mu}<60$ GeV, hm stands for $60<m_{\mu\mu}<120$ GeV, and i is the corresponding bin number in rapidity in the centre-of-mass frame. The global normalisation uncertainty of 3.5% is not included. Data from Figure [3, 4, 4].

Correlation matrix for the systematic uncertainties, excluding integrated luminosity, as a function of $p_{\textrm{T}}$, putting together low and high masses. Each bin is labeled lm_i or hm_i, where lm stands for $15<m_{\mu\mu}<60$ GeV, hm stands for $60<m_{\mu\mu}<120$ GeV, and i is the corresponding bin number in $p_{\textrm{T}}$. The global normalisation uncertainty of 3.5% is not included. Data from Figure [3, 4, 4].

Correlation matrix for the systematic uncertainties, excluding integrated luminosity, as a function of $\phi^*$, putting together low and high masses. Each bin is labeled lm_i or hm_i, where lm stands for $15<m_{\mu\mu}<60$ GeV, hm stands for $60<m_{\mu\mu}<120$ GeV, and i is the corresponding bin number in $\phi^*$. The global normalisation uncertainty of 3.5% is not included. Data from Figure [3, 4, 4].

Nuclear modification factor of $\Upsilon(1\mathrm{S})$ as a function of the average number of participants $\langle N_{\mathrm{part}} \rangle$ or as a function of the collision centrality.

Nuclear modification factor of $\Upsilon(2\mathrm{S})$ as a function of the average number of participants $\langle N_{\mathrm{part}} \rangle$ or as a function of the collision centrality.

Ratio of $\Upsilon(2\mathrm{S})$ and $\Upsilon(1\mathrm{S})$ yields (cf. equation 2 in the article for the definition) as a function of the average number of participants $\langle N_{\mathrm{part}} \rangle$ or as a function of the collision centrality. The global uncertainty is the quadratic sum of the branching ratio uncertainties.

Relative nuclear modification factor or double yield ratio between $\Upsilon(2\mathrm{S})$ and $\Upsilon(1\mathrm{S})$ as a function of the average number of participants $\langle N_{\mathrm{part}} \rangle$ or as a function of the collision centrality. The global uncertainty corresponds to the systematic uncertainty on the cross-section ratio in proton–proton collisions.

Nuclear modification factor of $\Upsilon(1\mathrm{S})$ as a function of transverse momentum for the 0–90% centrality interval.

Nuclear modification factor of $\Upsilon(1\mathrm{S})$ as a function of rapidity for the 0–90% centrality interval.

Nuclear modification factor of $\Upsilon(2\mathrm{S})$ as a function of rapidity for the 0–90% centrality interval.

Interpolated production cross sections and cross-section ratio for $\Upsilon(1\mathrm{S}) \rightarrow \mu^{+}\mu^{-}$ and $\Upsilon(2\mathrm{S}) \rightarrow \mu^{+}\mu^{-}$ in proton–proton collisions, integrated over the rapidity interval 2.5–4.0 and $p_{\mathrm{T}} <$ 15 GeV/$c$. Results used for the determination of the integrated nuclear modification factors and of the nuclear modification factors as a function of centrality (Figures 4 and 5 (right)).

Interpolated $p_{\mathrm{T}}$-differential cross section for $\Upsilon(1\mathrm{S}) \rightarrow \mu^{+}\mu^{-}$ in proton–proton collisions. Results used for the determination of the nuclear modification factor of $\Upsilon(1\mathrm{S})$ as a function of its transverse momentum (Figure 6).

Interpolated rapidity-differential cross section for $\Upsilon(1\mathrm{S}) \rightarrow \mu^{+}\mu^{-}$ in proton–proton collisions. Results used for the determination of the nuclear modification factor of $\Upsilon(1\mathrm{S})$ as a function of rapidity (Figure 7).

Interpolated rapidity-differential cross section for $\Upsilon(2\mathrm{S}) \rightarrow \mu^{+}\mu^{-}$ in proton–proton collisions. Results used for the determination of the nuclear modification factor of $\Upsilon(2\mathrm{S})$ as a function of rapidity (Figure 7).

The nuclear modification factor R_pPb for prompt, isolated photons with rapidity in (1.09,1.90).

The nuclear modification factor R_pPb for prompt, isolated photons with rapidity in (−1.84,0.91).

The nuclear modification factor R_pPb for prompt, isolated photons with rapidity in (−2.83,−2.02).

The ratio of R_{pPb} from rapidity (1.09,1.90) to that of rapidity (−2.83,−2.02).

pT-differential cross section of inclusive Dzero mesons in p-Pb collisions at sqrt{sNN}=5.02 TeV in the rapidity interval -0.96<y<0.04. Branching ratio of D0->Kpi : 0.0388.

pT-differential cross section of prompt Dzero mesons in p-Pb collisions at sqrt{sNN}=5.02 TeV in the rapidity interval -0.96<y<0.04. Branching ratio of D0->Kpi : 0.0388. Data points for pt<2 GeV/c from analysis "without vertexing". Data points for pt>2 GeV/c from the analysis "with vertexing" taken from PRL 113 (2014) 232301 (http://hepdata.cedar.ac.uk/view/ins1296081).

First column: production cross sections per unit of rapidity for prompt D0 mesons, inclusive D0 mesons (no feed-down subtraction) and charm quarks in -0.96 < ycm < 0.04 in p-Pb collisions at 5.02 TeV. For D0 mesons, the second (sys) error is from the luminosity uncertainty, the third (sys) error is from the branching-ratio uncertainties. For charm quarks, the second (sys) error is from the luminosity uncertainty, the third (sys) error is from the Fragmentation Function uncertainties, the fourth (sys) error is from the rapidity shapes of D0 mesons and single charm quarks. Second column: value of <pT> of prompt D0 mesons in p-Pb collisions at 5.02 TeV in the rapidity interval -0.96 < y < 0.04. The first uncertainty is statistical, the second is the systematic uncertainty.

Ratio of D+ to D0 meson yield in p-Pb collisions at sqrt{sNN}=5.02 TeV in -0.96<y< 0.04 as a function of pT. Branching ratio of D+->Kpipi : 0.0913. Branching ratio of D0->Kpi : 0.0388.

Ratio of D*+ to D0 meson yield in p-Pb collisions at sqrt{sNN}=5.02 TeV in -0.96<y< 0.04 as a function of pT. Branching ratio of D*+->D0pi->Kpipi : 0.0388*0.677. Branching ratio of D0->Kpi : 0.0388.

Ratio of Ds to D0 meson yield in p-Pb collisions at sqrt{sNN}=5.02 TeV in -0.96<y< 0.04 as a function of pT. Branching ratio of Ds->phipi->KKpi : 0.0224. Branching ratio of D0->Kpi : 0.0388.

Ratio of Ds to D+ meson yield in p-Pb collisions at sqrt{sNN}=5.02 TeV in -0.96<y< 0.04 as a function of pT. Branching ratio of Ds->phipi->KKpi : 0.0224. Branching ratio of D+->Kpipi : 0.0913.

cross section of prompt Dzero mesons in p-Pb collisions at sqrt{sNN}=5.02 TeV vs rapidity in the pt interval 2<pt<5 GeV/c. Branching ratio of D0->Kpi : 0.0388. The second (sys) error is the systematic uncertainty from the B feed-down contribution. The first (sys) error is the systematic uncertainty from the other sources.

cross section of prompt Dzero mesons in p-Pb collisions at sqrt{sNN}=5.02 TeV vs rapidity in the pt interval 5<pt<8 GeV/c. Branching ratio of D0->Kpi : 0.0388. The second (sys) error is the systematic uncertainty from the B feed-down contribution. The first (sys) error is the systematic uncertainty from the other sources.

cross section of prompt Dzero mesons in p-Pb collisions at sqrt{sNN}=5.02 TeV vs rapidity in the pt interval 8<pt<16 GeV/c. Branching ratio of D0->Kpi : 0.0388. The second (sys) error is the systematic uncertainty from the B feed-down contribution. The first (sys) error is the systematic uncertainty from the other sources.

cross section of prompt Dplus mesons in p-Pb collisions at sqrt{sNN}=5.02 TeV vs rapidity in the pt interval 2<pt<5 GeV/c. Branching ratio of D+->Kpipi : 0.0913. The second (sys) error is the systematic uncertainty from the B feed-down contribution. The first (sys) error is the systematic uncertainty from the other sources.

cross section of prompt Dplus mesons in p-Pb collisions at sqrt{sNN}=5.02 TeV vs rapidity in the pt interval 5<pt<8 GeV/c. Branching ratio of D+->Kpipi : 0.0913. The second (sys) error is the systematic uncertainty from the B feed-down contribution. The first (sys) error is the systematic uncertainty from the other sources.

cross section of prompt Dplus mesons in p-Pb collisions at sqrt{sNN}=5.02 TeV vs rapidity in the pt interval 8<pt<16 GeV/c. Branching ratio of D+->Kpipi : 0.0913. The second (sys) error is the systematic uncertainty from the B feed-down contribution. The first (sys) error is the systematic uncertainty from the other sources.

cross section of prompt Dstar mesons in p-Pb collisions at sqrt{sNN}=5.02 TeV vs rapidity in the pt interval 2<pt<5 GeV/c. Branching ratio of D*+->D0pi->Kpipi : 0.0388*0.677. The second (sys) error is the systematic uncertainty from the B feed-down contribution. The first (sys) error is the systematic uncertainty from the other sources.

cross section of prompt Dstar mesons in p-Pb collisions at sqrt{sNN}=5.02 TeV vs rapidity in the pt interval 5<pt<8 GeV/c. Branching ratio of D*+->D0pi->Kpipi : 0.0388*0.677. The second (sys) error is the systematic uncertainty from the B feed-down contribution. The first (sys) error is the systematic uncertainty from the other sources.

cross section of prompt Dstar mesons in p-Pb collisions at sqrt{sNN}=5.02 TeV vs rapidity in the pt interval 8<pt<16 GeV/c. Branching ratio of D*+->D0pi->Kpipi : 0.0388*0.677. The second (sys) error is the systematic uncertainty from the B feed-down contribution. The first (sys) error is the systematic uncertainty from the other sources.

Nuclear modification factor RpPb of D0 mesons in p-Pb collisions at sqrt{sNN}=5.02 TeV in -0.96<y<0.04 as a function of pT. Data point for 0<pt<1 GeV/c from analysis "without vertexing". Data points for pt>1 GeV/c from the analysis "with vertexing" taken from PRL 113 (2014) 232301 (http://hepdata.cedar.ac.uk/view/ins1296081).

Pt-integrated nuclear modification factor RpA of prompt D0 mesons in -0.96 < ycm < 0.04. The first error is the statistical uncertainty, the second is the systematic uncertainty.

Nuclear modicifation factor vs. N_coll $J\psi$ with|y|<1 and pT < 3 GeV/c in d Au collisions.

The invariant cross section versus transverse momentum for J/$\psi$ with $|y|\lt 1$ in $p+p$ collisions.

The invariant yield versus transverse momentum for $J/\Psi$ with |y| < 1 in 0 - 100 percent central d+Au collisions.

Differential production cross sections of $\psi$(2S) as a function of $p_{\rm T}$.

Differential production cross sections of $\psi$(2S) as a function of rapidity.

integrated production cross section of $\psi$(2S).

Differential production cross sections of $\Upsilon$(1S) as a function of $p_{\rm T}$.

Differential production cross sections of $\Upsilon$(1S) as a function of rapidity.

integrated production cross section of $\Upsilon$(1S).

Differential production cross sections of $\Upsilon$(2S) as a function of rapidity.

integrated production cross section of Upsilon (2S).

integrated production cross section of Upsilon (3S).

$\psi$(2S) over J/$\psi$ ratio as a function of $p_{\rm T}$.

$\psi$(2S) over J/$\psi$ ratio as a function of rapidity.

integrated $\psi$(2S) over J/$\psi$ production cross section ratio.

integrated $\Upsilon$(2S) over $\Upsilon$(1S) production cross section ratio.

integrated $\Upsilon$(3S) over $\Upsilon$(1S) production cross section ratio.

No description provided.

No description provided.

No description provided.

No description provided.

No description provided.

No description provided.

No description provided.

No description provided.

No description provided.

No description provided.

No description provided.

No description provided.

No description provided.

No description provided.

No description provided.

Forward cross sections.

Backward cross sections.

From S, P and D wave fit to d(sig)/domega.