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 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 prompt and non-prompt production cross-sections for $\psi(2S)$ mesons are measured using 2.1 fb$^{-1}$ of $pp$ collision data at a centre-of-mass energy of 7 TeV recorded by the ATLAS experiment at the LHC. The measurement exploits the $\psi(2S)\to J/\psi(\to\mu^+\mu^-)\pi^+\pi^-$ decay mode, and probes $\psi(2S)$ mesons with transverse momenta in the range $10\leq p_T<100$ GeV and rapidity $|y|<2.0$. The results are compared to other measurements of $\psi(2S)$ production at the LHC and to various theoretical models for prompt and non-prompt quarkonium production.

Production of inclusive charmonia in pp collisions at center-of-mass energy of $\sqrt{s}$ = 13 TeV and p-Pb collisions at center-of-mass energy per nucleon pair of $\sqrt{s_{\rm NN}}$ = 8.16 TeV is studied as a function of charged-particle pseudorapidity density with ALICE. Ground and excited charmonium states (J/$\psi$, $\psi$(2S)) are measured from their dimuon decays in the interval of rapidity in the center-of-mass frame $2.5 < y_{\rm cms} < 4.0$ for pp collisions, and $2.03 < y_{\rm cms} < 3.53$ and $-4.46 < y_{\rm cms} < -2.96$ for p-Pb collisions. The charged-particle pseudorapidity density is measured around midrapidity ($|\eta|<1.0$). In pp collisions, the measured charged-particle multiplicity extends to about six times the average value, while in p-Pb collisions at forward (backward) rapidity a multiplicity corresponding to about three (four) times the average is reached. The $\psi$(2S) yield increases with the charged-particle pseudorapidity density. The ratio of $\psi$(2S) over J/$\psi$ yield does not show a significant multiplicity dependence in either colliding system, suggesting a similar behavior of J/$\psi$ and $\psi$(2S) yields with respect to charged-particle pseudorapidity density. Results for the $\psi$(2S) yield and its ratio with respect to J/$\psi$ agree with available model calculations.

The production rates of prompt and non-prompt $J/\psi$ and $\psi(2\mathrm{S})$ mesons are measured using 2.1 $fb^{-1}$ and 11.4 $fb^{-1}$ of data collected with the ATLAS experiment at the LHC, in proton-proton collisions at $\sqrt{s}=7$ and 8 TeV respectively. Production cross-sections for both prompt and non-prompt production sources, ratios of $\psi(2\mathrm{S})$ to $J/\psi$ production, and fractions of non-prompt to inclusive production for $J/\psi$ and $\psi(2\mathrm{S})$ are measured double-differentially as a function of meson $p_{T}$ and rapidity. These measurements are made in a restricted fiducial volume and also corrected for geometrical acceptance after which they are compared to a variety of theoretical predictions.

The differential cross-section for the inclusive production of $\psi(2S)$ mesons in $pp$ collisions at $\sqrt{s}$ =7 TeV has been measured with the LHCb detector. The data sample corresponds to an integrated luminosity of 36 pb$^{-1}$ collected during 2010. The $\psi(2S)$ mesons are reconstructed in the decay channels $\psi(2S) \rightarrow \mu^{+} \mu^{-}$ and $\psi(2S)$ $\rightarrow$ $J$/$\psi \pi^{+}\pi^{-}$, with the $J$/$\psi$ meson decaying into two muons. Results are presented both for promptly produced $\psi(2S)$ mesons and for those originating from $b$-hadron decays. In the kinematic range $p_{\rm T}(\psi(2S))$ $\leq$ 16 GeV/$c$ and 2 < $y$($\psi(2S))$ $\leq$ 4.5 we measure \begin{eqnarray*} \sigma_{\rm prompt}(\psi(2S)) &=& 1.44 \pm 0.01~(\text{stat})\pm 0.12~(\text{syst})^{+0.20}_{-0.40}~(\text{pol})~{\rm \mu b}, \\ \sigma_{b}(\psi(2S)) &=& 0.25 \pm 0.01~(\text{stat}) \pm 0.02~(\text{syst})~{\rm \mu b}, \end{eqnarray*} where the last uncertainty on the prompt cross-section is due to the unknown $\psi(2S)$ polarization. Recent QCD calculations are found to be in good agreement with our measurements. Combining the present result with the LHCb $J$/$\psi$ measurements we determine the inclusive branching fraction \begin{equation*} \mathcal{B}(b \rightarrow \psi(2S) X) = (2.73 \pm 0.06~(\text{stat}) \pm 0.16~(\text{syst}) \pm 0.24~(\text{BF})) \times 10^{-3}, \end{equation*} where the last uncertainty is due to the $\mathcal{B}(b \rightarrow J/\psi X)$, $\mathcal{B}(J/\psi \rightarrow \mu^+ \mu^-)$ and $\mathcal{B}(\psi(2S) \rightarrow e^+ e^-)$ branching fraction uncertainties. All above results are corrected by an erratum included as an appendix.

Differential production cross sections of J/$\psi$ and $\psi$(2S) charmonium and $\Upsilon$(nS) (n = 1, 2, 3) bottomonium states are measured in proton-proton collisions at $\sqrt{s} =$ 13 TeV, with data collected by the CMS detector at the LHC, corresponding to an integrated luminosity of 2.3 fb$^{-1}$ for the J/$\psi$ and 2.7 fb$^{-1}$ for the other mesons. The five quarkonium states are reconstructed in the dimuon decay channel, for dimuon rapidity $|y| <$ 1.2. The double-differential cross sections for each state are measured as a function of $y$ and transverse momentum, and compared to theoretical expectations. In addition, ratios are presented of cross sections for prompt $\psi$(2S) to J/$\psi$, $\Upsilon$(2S) to $\Upsilon$(1S), and $\Upsilon$(3S) to $\Upsilon$(1S) production.

Measurements of the differential production cross-sections of prompt and non-prompt $J/\psi$ and $\psi(2$S$)$ mesons with transverse momenta between 8 and 360 GeV and rapidity in the range $|y|<2$ are reported. Furthermore, measurements of the non-prompt fractions of $J/\psi$ and $\psi(2$S$)$, and the prompt and non-prompt $\psi(2$S$)$-to-$J/\psi$ production ratios, are presented. The analysis is performed using 140 fb$^{-1}$ of $\sqrt{s}=13$ TeV $pp$ collision data recorded by the ATLAS detector at the LHC during the years 2015-2018.

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.