We report on the inclusive J/$\psi$ production cross section measured at the CERN Large Hadron Collider in proton-proton collisions at a centre-of-mass energy $\sqrt{s}~=~13$ TeV. The J/$\psi$ mesons are reconstructed in the $\rm e^{+} e^{-}$ decay channel and the measurements are performed at midrapidity ($|y|<0.9$) in the transverse-momentum interval $0< p_{\rm T} <40$ GeV/$c$, using a minimum bias data sample corresponding to an integrated luminosity $L_{\text{int}} = 32.2~\text{nb}^{-1}$ and an Electromagnetic Calorimeter triggered data sample with $L_{\text{int}} = 8.3~\mathrm{pb}^{-1}$. The $p_{\rm T}$-integrated J/$\psi$ production cross section at midrapidity, computed using the minimum bias data sample, is $\text{d}\sigma/\text{d}y|_{y=0} = 8.97\pm0.24~(\text{stat})\pm0.48~(\text{syst})\pm0.15~(\text{lumi})~\mu\text{b}$. An approximate logarithmic dependence with the collision energy is suggested by these results and available world data, in agreement with model predictions. The integrated and $p_{\rm T}$-differential measurements are compared with measurements in pp collisions at lower energies and with several recent phenomenological calculations based on the non-relativistic QCD and Color Evaporation models.

The nuclear modification factors of J/$\psi$ and $\psi$(2S) mesons are measured in PbPb collisions at a centre-of-mass energy per nucleon pair of $\sqrt{s_{\mathrm{NN}}} =$ 5.02 TeV. The analysis is based on PbPb and pp data samples collected by CMS at the LHC in 2015, corresponding to integrated luminosities of 464 $\mu$b$^{-1}$ and 28 pb$^{-1}$, respectively. The measurements are performed in the dimuon rapidity range of $|y| <$ 2.4 as a function of centrality, rapidity, and transverse momentum (p$_\mathrm{T}$) from p$_\mathrm{T}=$ 3 GeV/$c$ in the most forward region and up to 50 GeV/$c$. Both prompt and nonprompt (coming from b hadron decays) mesons are observed to be increasingly suppressed with centrality, with a magnitude similar to the one observed at $\sqrt{s_{\mathrm{NN}}}=$ 2.76 TeV for the two J/$\psi$ meson components. No dependence on rapidity is observed for either prompt or nonprompt J/$\psi$ mesons. An indication of a lower prompt J/$\psi$ meson suppression at p$_\mathrm{T} >$ 25 GeV/$c$ is seen with respect to that observed at intermediate p$_\mathrm{T}$. The prompt $\psi$(2S) meson yield is found to be more suppressed than that of the prompt J/$\psi$ mesons in the entire p$_\mathrm{T}$ range.

The modification of the production of $J/\psi$, $\psi(\mathrm{2S})$, and $\mit{\Upsilon}(n\mathrm{S})$ ($n = 1, 2, 3$) in $p$+Pb collisions with respect to their production in $pp$ collisions has been studied. The $p$+Pb and $pp$ datasets used in this paper correspond to integrated luminosities of $28$ $\mathrm{nb}^{-1}$ and $25$ $\mathrm{pb}^{-1}$ respectively, collected in 2013 and 2015 by the ATLAS detector at the LHC, both at a centre-of-mass energy per nucleon pair of 5.02 TeV. The quarkonium states are reconstructed in the dimuon decay channel. The yields of $J/\psi$ and $\psi(\mathrm{2S})$ are separated into prompt and non-prompt sources. The measured quarkonium differential cross sections are presented as a function of rapidity and transverse momentum, as is the nuclear modification factor, $R_{p\mathrm{Pb}}$ for $J/\psi$ and $\mit{\Upsilon}(\mathrm{1S})$. No significant modification of the $J/\psi$ production is observed while $\mit{\Upsilon}(\mathrm{1S})$ production is found to be suppressed at low transverse momentum in $p$+Pb collisions relative to $pp$ collisions. The production of excited charmonium and bottomonium states is found to be suppressed relative to that of the ground states in central $p$+Pb collisions.

The inclusive production cross sections at forward rapidity of J/$\psi$, $\psi$(2S), $\Upsilon$(1S) and $\Upsilon$(2S) are measured in pp collisions at $\sqrt{s} = 7$ TeV with the ALICE detector at the LHC. The analysis is based in a data sample corresponding to an integrated luminosity of 1.35 pb$^{-1}$. Quarkonia are reconstructed in the dimuon-decay channel and the signal yields are evaluated by fitting the $\mu^+\mu^-$ invariant mass distributions. The differential production cross sections are measured as a function of the transverse momentum $p_{\rm T}$ and rapidity $y$, over the ranges $0 < p_{\rm T} < 20$ GeV/$c$ for J/$\psi$, $0 < p_{\rm T} < 12$ GeV/$c$ for all other resonances and for $2.5 < y < 4$. The measured cross sections integrated over $p_{\rm T}$ and $y$, and assuming unpolarized quarkonia, are: $\sigma_{J/\psi} = 6.69 \pm 0.04 \pm 0.63$ $\mu$b, $\sigma_{\psi^{\prime}} = 1.13 \pm 0.07 \pm 0.14$ $\mu$b, $\sigma_{\Upsilon{\rm(1S)}} = 54.2 \pm 5.0 \pm 6.7$ nb and $\sigma_{\Upsilon{\rm (2S)}} = 18.4 \pm 3.7 \pm 2.2$ nb, where the first uncertainty is statistical and the second one is systematic. The results are compared to measurements performed by other LHC experiments and to theoretical models.

The ratio between the prompt psi(2S) and J/psi yields, reconstructed via their decays into muon pairs, is measured in PbPb and pp collisions at sqrt(s[NN]) = 2.76 TeV. The analysis is based on PbPb and pp data samples collected by CMS at the LHC, corresponding to integrated luminosities of 150 inverse microbarns and 5.4 inverse picobarns, respectively. The double ratio of measured yields, (N[psi(2S)]/N[J/psi])[PbPb] / (N[psi(2S)]/ N[J/psi])[pp], is computed in three PbPb collision centrality bins and two kinematic ranges: one at midrapidity, abs(y) < 1.6, covering the transverse momentum range 6.5 < pt < 30 GeV/c, and the other at forward rapidity, 1.6 < abs(y) < 2.4, extending to lower pt values, 3 < pt < 30 GeV/c. The centrality-integrated double ratio changes from 0.45 +/- 0.13 (stat) +/- 0.07 (syst) in the first range to 1.67 +/- 0.34 (stat) +/- 0.27 (syst) in the second. This difference is most pronounced in the most central collisions.

We present a measurement of inclusive J/$\psi$ production in p-Pb collisions at $\sqrt{s_{\rm NN}}$ = 5.02 TeV as a function of the centrality of the collision, as estimated from the energy deposited in the Zero Degree Calorimeters. The measurement is performed with the ALICE detector down to zero transverse momentum, $p_{\rm T}$, in the backward ($-4.46 < y_{\rm cms} < -2.96$) and forward ($2.03 < y_{\rm cms} < 3.53$) rapidity intervals in the dimuon decay channel and in the mid-rapidity region ($-1.37 < y_{\rm cms} < 0.43$) in the dielectron decay channel. The backward and forward rapidity intervals correspond to the Pb-going and p-going direction, respectively. The $p_{\rm T}$-differential J/$\psi$ production cross section at backward and forward rapidity is measured for several centrality classes, together with the corresponding average $p_{\rm T}$ and $p^2_{\rm T}$ values. The nuclear modification factor, $Q_{\rm pPb}$, is presented as a function of centrality for the three rapidity intervals, and, additionally, at backward and forward rapidity, as a function of $p_{\rm T}$ for several centrality classes. At mid- and forward rapidity, the J/$\psi$ yield is suppressed up to 40% compared to that in pp interactions scaled by the number of binary collisions. The degree of suppression increases towards central p-Pb collisions at forward rapidity, and with decreasing $p_{\rm T}$ of the J/$\psi$. At backward rapidity, the $Q_{\rm pPb}$ is compatible with unity within the total uncertainties, with an increasing trend from peripheral to central p-Pb collisions.

The nuclear modification factor RAA and the azimuthal anisotropy coefficient v[2] of prompt and nonprompt (i.e. those from decays of b hadrons) J/psi mesons, measured from PbPb and pp collisions at sqrt(s[NN]) = 2.76 TeV at the LHC, are reported. The results are presented in several event centrality intervals and several kinematic regions, for transverse momenta pt > 6.5 GeV/c and rapidity abs(y) < 2.4, extending down to pt = 3 GeV/c in the 1.6 < abs(y) < 2.4 range. The v[2] of prompt J/psi is found to be nonzero, but with no strong dependence on centrality, rapidity, or pt over the full kinematic range studied. The measured v[2] of nonprompt J/psi is consistent with zero. The RAA of prompt J/psi exhibits a suppression that increases from peripheral to central collisions but does not vary as a function of either y or pt in the fiducial range. The nonprompt J/psi RAA shows a suppression which becomes stronger as rapidity or pt increase. The v[2] and nuclear suppression of open and hidden charm, and of open charm and beauty, are compared.

We present results on transverse momentum ($p_{\rm T}$) and rapidity ($y$) differential production cross sections, mean transverse momentum and mean transverse momentum square of inclusive J/$\psi$ and $\psi(2S)$ at forward rapidity ($2.5<y<4$) as well as $\psi(2S)$-to-J/$\psi$ cross section ratios. These quantities are measured in pp collisions at center of mass energies $\sqrt{s}=5.02$ and 13 TeV with the ALICE detector. Both charmonium states are reconstructed in the dimuon decay channel, using the muon spectrometer. A comprehensive comparison to inclusive charmonium cross sections measured at $\sqrt{s}=2.76$, 7 and 8 TeV is performed. A comparison to non-relativistic quantum chromodynamics and fixed-order next-to-leading logarithm calculations, which describe prompt and non-prompt charmonium production respectively, is also presented. A good description of the data is obtained over the full $p_{\rm T}$ range, provided that both contributions are summed. In particular, it is found that for $p_{\rm T}>15$ GeV/$c$ the non-prompt contribution reaches up to 50% of the total charmonium yield.

The relative modification of the prompt psi(2S) and J/psi yields from pp to PbPb collisions, at the center of mass energy of 5.02 TeV per nucleon pair, is presented. The analysis is based on pp and PbPb data samples collected by the CMS experiment at the LHC in 2015, corresponding to integrated luminosities of 28.0 inverse picobarns and 464 inverse microbarns, respectively. The double ratio of measured yields of prompt charmonia reconstructed through their decays into muon pairs, (N[psi(2S)]/N[J/psi])[PbPb] / (N[psi(2S)]/N[J/psi])[pp], is determined as a function of PbPb collision centrality and charmonium transverse momentum pt, in two kinematic intervals: abs(y) < 1.6 covering 6.5 < pt < 30 GeV/c and 1.6 < abs(y) < 2.4 covering 3 < pt < 30 GeV/c. The centrality-integrated double ratios are 0.36 +/- 0.08 (stat) +/-0.05 (syst) in the first interval and 0.24 +/- 0.22 (stat) +/- 0.09 (syst) in the second. The double ratio is lower than unity in all the measured bins, suggesting that the psi(2S) yield is more suppressed than the J/psi yield in the explored phase space.

Inclusive J/$\psi$ production is studied in minimum-bias proton-proton collisions at a centre-of-mass energy of $\sqrt{s}$ = 5.02 TeV by ALICE at the CERN LHC. The measurement is performed at mid-rapidity ($|y| < 0.9$) in the dielectron decay channel down to zero transverse momentum $p_{\rm T}$, using a data sample corresponding to an integrated luminosity of $L_{\rm int} = 19.4 \pm$ 0.4 nb$^{-1}$. The measured $p_{\rm T}$-integrated inclusive J/$\psi$ production cross section is d$\sigma$/d$y$ = 5.64 $\pm$ 0.22 (stat.) $\pm 0.33$ (syst.) $\pm 0.12$ (lumi.) $\mu$b. The $p_{\rm T}$-differential cross section d$^{2} \sigma$/d$p_{\rm T}$d$y$ is measured in the $p_{\rm T}$ range 0$-$10 GeV/$c$ and compared with state-of-the-art perturbative QCD calculations. The J/$\psi$ $\langle p_{\rm T} \rangle$ and $\langle p_{\rm T}^{2} \rangle$ are extracted and compared with results obtained at other collision energies.