Examples of fit to the OS dimuon invariant mass distribution in the mass region $2 < m_{\mu^{+}\mu^{-}} < 5$ GeV/$c^{2}$ for $\pt < 20$ GeV/$c$ (left), and $7< m_{\mu^{+}\mu^{-}} < 13$ GeV/$c^{2}$ for $\pt < 15$ GeV/$c$ (right). The J/$\psi$, $\psi$(2S) and $\Upsilon$(nS) signals are modelled with extended Crystal Ball functions, while the background is described by a pseudo Gaussian with a width increasing linearly with the invariant mass. The fit is performed on the full data sample. The widths of the $\psi$(2S) , $\Upsilon$(2S) and $\Upsilon$(3S), for these examples, are fixed to 73 MeV/$c^2$, 156 MeV/$c^2$ and 161 MeV/$c^2$, respectively.

Transverse momentum dependence of the inclusive J/$\psi$ cross section. The measurements are compared to theoretical calculations from Refs. [44–46] (left) and Refs. [47–49] (right). The calculations of the non-prompt contribution [49] are also shown separately. See text for details.

Rapidity dependence of the inclusive J/$\psi$ cross section. The measurements are compared to theoretical calculations from Refs. [46, 47].

Rapidity dependence of the inclusive J/$\psi$ cross section for various $p_{\rm T}$ ranges, compared to theoretical calculations [44–49]. The theoretical calculations are scaled as for the data for $0.3<\pt<2$ GeV/$c$. See text for details.

The left and right panels show the $p_{\rm T}$ dependence for the inclusive $\psi$(2S) production cross section in pp collisions at $\s = 5.02$ TeV. The results are compared with the theory predictions based on NRQCD [44–46] (left) and CEM [47,48] (right) models. The calculation of the non-prompt contribution from FONLL calculations [49] are also shown separately. See text for details.

Rapidity dependence for the inclusive $\psi$(2S) production cross section in pp collisions at $\s = 5.02$ TeV. The results are compared with the theory predictions based on NRQCD+CGC  and ICEM [47] models. See text for details.

The inclusive $\psi$(2S)-to-J$\psi$ cross section ratio as a function of $p_{\rm T}$ (left) and $y$ (right), compared with theoretical calculations [44, 47, 49]. See text for details.

Transverse momentum dependence of the $\Upsilon$(1S) cross section (left) and $y$ dependence of the $\Upsilon$(1S), $\Upsilon$(2S), and $\Upsilon$(3S) (right) measured by ALICE and CMS. The two panels also show theoretical calculations [48,51]. See text for details.

Transverse momentum dependence of the inclusive J/$\psi$ cross section, at forward $y$, measured in pp collisions at $\sqrt{s} = 5.02$, 7 [12], 8 [13], and 13 [9] TeV (top panels), and ratio of the measurements at 5.02, 7, and 8 TeV to the 13 TeV data (bottom panels). The data are compared with the NRQCD theoretical calculations from Butenschön et al + FONLL (left panels) [44, 49] and with theoretical calculations from ICEM + FONLL (right panels) [47,49].

Rapidity dependence of the inclusive J/$\psi$ (left) and $\psi$(2S) (right) cross section, at forward $y$, measured in pp collisions at $\s = 5.02$, 7 [12], 8 [13], and 13 [9] TeV (top panels), and ratio of the measurements at 5.02, 7, and 8 TeV to the 13 TeV data (bottom panels). The data are compared with theoretical calculations from ICEM + FONLL [47, 49].

Transverse momentum dependence of the inclusive $\psi$(2S) cross section, at forward $y$, measured in pp collisions at $\s = 5.02$, 7 [12], 8 [13], and 13 [9] TeV (top panels), and ratio of the measurements at 5.02, 7, 8 TeV to the 13 TeV data (bottom panels). The data are compared with the NRQCD theoretical calculations from Butenschön et al + FONLL (left panels) [44,49] and with theoretical calculations from ICEM + FONLL (right panels) [47,49].

Inclusive $\psi$(2S)-to-J/$\psi$ cross section ratio as a function of $p_{\rm T}$, at forward $y$, in pp collisions at $\s = 5.02$, 7 [12], 8 [13], and 13 [9] TeV (left panel). The data at $\s = 5.02$ TeV are compared with NRQCD theoretical calculations from Butenschön et al + FONLL [44,49] and with theoretical calculations from ICEM + FONLL [47,49] (right panel).

Inclusive $\psi$(2S)-to-J/$\psi$ cross section ratio as a function of $y$ in pp collisions at $\s = 5.02$, 7 [12], 8 [13], and 13 [9] TeV (left panel). The data at $\s = 5.02$ TeV are compared with theoretical calculations from ICEM + FONLL [47,49] (right panel).

Inclusive $\psi$(2S)-to-J/$\psi$ cross section ratio (left) and J/$\psi$, $\psi$(2S), $\Upsilon$(1S), $\Upsilon$(2S), and $\Upsilon$(3S) $p_{\rm T}$-integrated cross section per unit of rapidity (right) as a function of the collision energy in pp collisions [9, 11–13]. In the left panel, the systematic boxes include the BR uncertainties from both resonances, on top of the MC input and signal extraction systematic uncertainties. The 13 TeV data point is computed from the published individual J/$\psi$ and $\psi$(2S) $p_{\rm T}$-integrated cross sections. The statistical and systematic uncertainties are assumed to be uncorrelated between the resonances when computing the ratio. In the right panel, the luminosity and branching ratio uncertainties are included in the systematic boxes. The data are compared with theoretical calculations from ICEM + FONLL[47,49].

Inclusive $\psi$(2S)-to-J/$\psi$ cross section ratio as a function of $p_{\rm T}$, at forward $y$, in pp collisions at $\s$ = 7 [12] (top left), 8 [13] (top right), and 13 TeV [9] (bottom left). The data are compared with NRQCD theoretical calculations from Butenschön et al + FONLL [44,49] and with theoretical calculations from ICEM + FONLL [47,49].

Inclusive $\psi$(2S)-to-J/$\psi$ cross section ratio as a function of $y$ in pp collisions at $\s$ = 7 [12] (top left), 8 [13] (top right), and 13 TeV [9] (bottom left). The data are compared with theoretical calculations from ICEM + FONLL [47,49].