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 cross-sections of $J/\psi$ mesons in proton-proton collisions at a centre-of-mass energy of $\sqrt{s}=5$ TeV are measured using a data sample corresponding to an integrated luminosity of $9.13\pm0.18~\text{pb}^{-1}$, collected by the LHCb experiment. The cross-sections are measured differentially as a function of transverse momentum, $p_{\text{T}}$, and rapidity, $y$, and separately for $J/\psi$ mesons produced promptly and from beauty hadron decays (nonprompt). With the assumption of unpolarised $J/\psi$ mesons, the production cross-sections integrated over the kinematic range $0<p_{\text{T}}<20~\text{GeV}/c$ and $2.0<y<4.5$ are $8.154\pm0.010\pm0.283~\mu\text{b}$ for prompt $J/\psi$ mesons and $0.820\pm0.003\pm0.034~\mu\text{b}$ for nonprompt $J/\psi$ mesons, where the first uncertainties are statistical and the second systematic. These cross-sections are compared with those at $\sqrt{s}=8$ TeV and $13$ TeV, and are used to update the measurement of the nuclear modification factor in proton-lead collisions for $J/\psi$ mesons at a centre-of-mass energy per nucleon pair of $\sqrt{s_{\text{NN}}}=5$ TeV. The results are compared with theoretical predictions.

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.

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.

We report on the measurement of the inclusive J/$\psi$ polarization parameters in pp collisions at a center of mass energy $\sqrt{s} = 8$ TeV with the ALICE detector at the LHC. The analysis is based on a data sample corresponding to an integrated luminosity of 1.23 pb$^{-1}$. J/$\psi$ resonances are reconstructed in their di-muon decay channel in the rapidity interval $2.5 < y < 4.0$ and over the transverse-momentum interval $2 < p_{\rm T} < 15$ GeV/$c$. The three polarization parameters ($\lambda_\theta$, $\lambda_\varphi$, $\lambda_{\theta\varphi}$) are measured as a function of $p_{\rm T}$ both in the helicity and Collins-Soper reference frames. The measured J/$\psi$ polarization parameters are found to be compatible with zero within uncertainties, contrary to expectations from all available predictions. The results are compared with the measurement in pp collisions at $\sqrt{s} = 7$ TeV.

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.

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.

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 report the double helicity asymmetry, $A_{LL}^{J/\psi}$, in inclusive $J/\psi$ production at forward rapidity as a function of transverse momentum $p_T$ and rapidity $|y|$. The data analyzed were taken during $\sqrt{s}=510$ GeV longitudinally polarized $p

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

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.

Prompt D meson and non-prompt J/$\psi$ yields are studied as a function of the multiplicity of charged particles produced in inelastic proton-proton collisions at a centre-of-mass energy of $\sqrt{s}=7$ TeV. The results are reported as a ratio between yields in a given multiplicity interval normalised to the multiplicity-integrated ones (relative yields). They are shown as a function of the multiplicity of charged particles normalised to the average value for inelastic collisions (relative charged-particle multiplicity). D$^0$, D$^+$ and D$^{*+}$ mesons are measured in five $p_{\rm T}$ intervals from 1 to 20 GeV/$c$ and for $|y|<0.5$ via their hadronic decays. The D-meson relative yield is found to increase with increasing charged-particle multiplicity. For events with multiplicity six times higher than the average multiplicity of inelastic collisions, a yield enhancement of a factor about 15 relative to the multiplicity-integrated yield in inelastic collisions is observed. The yield enhancement is independent of transverse momentum within the uncertainties of the measurement. The D$^0$-meson relative yield is also measured as a function of the relative multiplicity at forward pseudorapidity. The non-prompt J/$\psi$, i.e. the B hadron, contribution to the inclusive J/$\psi$ production is measured in the di-electron decay channel at central rapidity. It is evaluated for $p_{\rm T}>1.3$ GeV/$c$ and $|y|<0.9$, and extrapolated to $p_{\rm T}>0$. The fraction of non-prompt J/$\psi$ in the inclusive J/$\psi$ yields shows no dependence on the charged-particle multiplicity at central rapidity. Charm and beauty hadron relative yields exhibit a similar increase with increasing charged-particle multiplicity. The measurements are compared to PYTHIA 8, EPOS 3 and percolation calculations.

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 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 prompt and non-prompt production cross-sections for the chi_c1 and chi_c2 charmonium states are measured in pp collisions at sqrt(s) = 7 TeV with the ATLAS detector at the LHC using 4.5 fb^-1 of integrated luminosity. The chi_c states are reconstructed through the radiative decay chi_c -> J/psi gamma (with J/psi -> mu+mu-) where photons are reconstructed from gamma -> e+e- conversions. The production rate of the chi_c2 state relative to the chi_c1 state is measured for prompt and non-prompt chi_c as a function of J/psi transverse momentum. The prompt chi_c cross-sections are combined with existing measurements of prompt J/psi production to derive the fraction of prompt J/psi produced in feed-down from chi_c decays. The fractions of chi_c1 and chi_c2 produced in b-hadron decays are also measured.

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.

An angular analysis of the decay $J/\psi \rightarrow \mu^{+} \mu^{-}$ is performed to measure the polarization of prompt $J/\psi$ mesons produced in $pp$ collisions at $\sqrt{s}=7\; \mathrm{TeV}$. The dataset corresponds to an integrated luminosity of 0.37 fb$^{-1}$ collected with the LHCb detector. The measurement is presented as a function of transverse momentum, $p_\mathrm{T}$, and rapidity, $y$, of the $J/\psi$ meson, in the kinematic region $2 < p_\mathrm{T} < 15\; \mathrm{GeV}/c$ and $2.0<y<4.5$.

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 production of J/psi and Upsilon mesons in pp collisions at sqrt(s) = 8 TeV is studied with the LHCb detector. The J/psi and Upsilon mesons are reconstructed in the mu+mu- decay mode and the signal yields are determined with a fit to the mu+mu- invariant mass distributions. The analysis is performed in the rapidity range 2.0<y<4.5 and transverse momentum range 0<p_T<14(15) GeV/c of the J/psi(Upsilon) mesons. The J/psi and Upsilon production cross-sections and the fraction of J/psi mesons from b-hadron decays are measured as a function of the meson p_T and y.

The production of J/psi mesons is studied with the LHCb detector using data from pp collisions at sqrt(s)=2.76 TeV corresponding to an integrated luminosity of 71 nb^-1. The differential cross-section for inclusive J/psi production is measured as a function of its transverse momentum pT. The cross-section in the fiducial region 0<pT<12 GeV/c and rapidity 2.0<y<4.5 is measured to be 5.6 +/- 0.1(stat) +/- 0.4 (syst) mub, with the assumption of unpolarised J/psi production. The fraction of J/psi production from b-hadron decays is measured to be (7.1 +/- 0.6 (stat) +/- 0.7 (syst)) %.

The prompt production of charmonium $\chi_{c}$ and $J/\psi$ states is studied in proton-proton collisions at a centre-of-mass energy of $\sqrt{s}=7$ TeV at the Large Hadron Collider. The $\chi_{c}$ and $J/\psi$ mesons are identified through their decays $\chi_{c}\rightarrow J/\psi \gamma$ and $J/\psi\rightarrow \mu^+\mu^-$ using 36 pb$^{-1}$ of data collected by the LHCb detector in 2010. The ratio of the prompt production cross-sections for $\chi_{c}$ and $J/\psi$, $\sigma (\chi_{c}\rightarrow J/\psi \gamma)/ \sigma (J/\psi)$, is determined as a function of the $J/\psi$ transverse momentum in the range $2 &lt; p_{\mathrm T}^{J/\psi} &lt; 15$ GeV/$c$. The results are in excellent agreement with next-to-leading order non-relativistic expectations and show a significant discrepancy compared with the colour singlet model prediction at leading order, especially in the low $p_{\mathrm T}^{J/\psi}$ region.

The ALICE Collaboration has measured inclusive J/psi production in pp collisions at a center of mass energy sqrt(s)=2.76 TeV at the LHC. The results presented in this Letter refer to the rapidity ranges |y|<0.9 and 2.5<y<4 and have been obtained by measuring the electron and muon pair decay channels, respectively. The integrated luminosities for the two channels are L^e_int=1.1 nb^-1 and L^mu_int=19.9 nb^-1, and the corresponding signal statistics are N_J/psi^e+e-=59 +/- 14 and N_J/psi^mu+mu-=1364 +/- 53. We present dsigma_J/psi/dy for the two rapidity regions under study and, for the forward-y range, d^2sigma_J/psi/dydp_t in the transverse momentum domain 0<p_t<8 GeV/c. The results are compared with previously published results at sqrt(s)=7 TeV and with theoretical calculations.