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

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.

The production of $J/\psi$ mesons in proton-proton collisions at a centre-of-mass energy of $\sqrt{s}=13$ TeV is studied with the \lhcb detector. Cross-section measurements are performed as a function of the transverse momentum $p_\mathrm{T}$ and the rapidity $y$ of the $J/\psi$ meson in the region $p_\mathrm{T}<14\mathrm{GeV}/c$ and $2.0<y<4.5$, for both prompt $J/\psi$ mesons and $J/\psi$ mesons from $b$-hadron decays. The production cross-sections integrated over the kinematic coverage are $15.03\pm 0.03\pm 0.94$ $\mu$b for prompt $J/\psi$ and $2.25\pm 0.01\pm 0.14$ $\mu$b for $J/\psi$ from $b$-hadron decays, assuming zero polarization of the $J/\psi$ meson. The first uncertainties are statistical and the second systematic. The cross-section reported for $J/\psi$ mesons from $b$-hadron decays is used to extrapolate to a total $b\bar{b}$ cross-section. The ratios of the cross-sections with respect to $\sqrt{s}=8$ TeV are also determined.

The double-differential cross sections of promptly produced J/psi and psi(2S) mesons are measured in pp collisions at sqrt(s) = 7 TeV, as a function of transverse momentum pt and absolute rapidity abs(y). The analysis uses J/psi and psi(2S) dimuon samples collected by CMS, corresponding to integrated luminosities of 4.55 and 4.90 inverse femtobarns, respectively. The results are based on a two-dimensional analysis of the dimuon invariant mass and decay length, and extend to pt = 120 and 100 GeV for the J/psi and psi(2S), respectively, when integrated over the interval abs(y) < 1.2. The ratio of the psi(2S) to J/psi cross sections is also reported for abs(y) < 1.2, over the range 10 < pt < 100 GeV. These are the highest pt values for which the cross sections and ratio have been measured.

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.

The polarizations of prompt J/psi and psi(2S) mesons are measured in proton-proton collisions at sqrt(s) = 7 TeV, using a dimuon data sample collected by the CMS experiment at the LHC, corresponding to an integrated luminosity of 4.9 inverse femtobarns. The prompt J/psi and psi(2S) polarization parameters lambda[theta], lambda[phi], and lambda[theta,phi], as well as the frame-invariant quantity lambda(tilde), are measured from the dimuon decay angular distributions in three different polarization frames. The J/psi results are obtained in the transverse momentum range 14 < pt < 70 GeV, in the rapidity intervals abs(y) < 0.6 and 0.6 < abs(y) < 1.2. The corresponding psi(2S) results cover 14 < pt < 50 GeV and include a third rapidity bin, 1.2 < abs(y) < 1.5. No evidence of large transverse or longitudinal polarizations is seen in these kinematic regions, which extend much beyond those previously explored.

The ALICE collaboration reports the measurement of the inclusive J/psi yield as a function of charged particle pseudorapidity density dN_{ch}/deta in pp collisions at sqrt{s} = 7 TeV at the LHC. J/psi particles are detected for p_t > 0, in the rapidity interval |y| < 0.9 via decay into e+e-, and in the interval 2.5 < y < 4.0 via decay into mu+mu- pairs. An approximately linear increase of the J/psi yields normalized to their event average (dN_{J/psi}/dy)/<dN_{J/psi}/dy> with (dN_{ch}/deta)/<dN_{ch}/deta> is observed in both rapidity ranges, where dN_{ch}/deta is measured within |eta| < 1 and p_t > 0. In the highest multiplicity interval with <dN_{ch}/deta(bin)> = 24.1, corresponding to four times the minimum bias multiplicity density, an enhancement relative to the minimum bias J/psi yield by a factor of about 5 at 2.5 < y < 4 (8 at |y| < 0.9) is observed.