Production of prompt J/$\psi$ meson pairs in proton-proton collisions at $\sqrt{s}$ = 7 TeV is measured with the CMS experiment at the LHC in a data sample corresponding to an integrated luminosity of about 4.7 inverse-femtobarns. The two J/$\psi$ mesons are fully reconstructed via their decays into $\mu^+\mu^-$ pairs. This observation provides for the first time access to the high-transverse-momentum region of J/$\psi$ pair production where model predictions are not yet established. The total and differential cross sections are measured in a phase space defined by the individual J/$\psi$ transverse momentum ($p_T^{J/\psi}$) and rapidity (|$y^{J/\psi}$|): |$y^{J/\psi}$| lower than 1.2 for $p_T^{J/\psi}$ greater than 6.5 GeV/c; |$y^{J/\psi}$| in [1.2,1.43] for a $p_T$ threshold that scales linearly with |$y^{J/\psi}$| from 6.5 to 4.5 GeV/c; and |$y^{J/\psi}$| in [1.43,2.2] for $p_T^{J/\psi}$ greater than 4.5 GeV/c. The total cross section, assuming unpolarized prompt J/$\psi$ pair production is 1.49 $\pm$ 0.07 (stat.) $\pm$ 0.13 (syst.) nb. Different assumptions about the J/$\psi$ polarization imply modifications to the cross section ranging from -31% to +27%.

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 ALICE experiment has measured the inclusive J/$\psi$ production in Pb-Pb collisions at $\sqrt{s_{\mathrm{NN}} } = 2.76$ TeV down to zero transverse momentum in the rapidity range $2.5 < y < 4$. A suppression of the inclusive J/$\psi$ yield in Pb-Pb is observed with respect to the one measured in pp collisions scaled by the number of binary nucleon-nucleon collisions. The nuclear modification factor, integrated over the 0-80% most central collisions, is $0.545 \pm 0.032 \rm{(stat.)} \pm 0.083 \rm{(syst.)}$ and does not exhibit a significant dependence on the collision centrality. These features appear significantly different from measurements at lower collision energies. Models including J/$\psi$ production from charm quarks in a deconfined partonic phase can describe our data.

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.

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

We have measured the polarizations of J/ψ and ψ(2S) mesons as functions of their transverse momentum pT when they are produced promptly in the rapidity range |y|<0.6 with pT≥5  GeV/c. The analysis is performed using a data sample with an integrated luminosity of about 800  pb−1 collected by the CDF II detector. For both vector mesons, we find that the polarizations become increasingly longitudinal as pT increases from 5 to 30  GeV/c. These results are compared to the predictions of nonrelativistic quantum chromodynamics and other contemporary models. The effective polarizations of J/ψ and ψ(2S) mesons from B-hadron decays are also reported.

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 ALICE experiment at the LHC has studied J/psi production at mid-rapidity in pp collisions at sqrt{s}=7 TeV through its electron pair decay on a data sample corresponding to an integrated luminosity L_int = 5.6nb-1. The fraction of J/psi from the decay of long-lived beauty hadrons was determined for J/psi candidates with transverse momentum p_t>1.3 GeV/c and rapidity |y|<0.9. The cross section for prompt J/psi mesons, i.e. directly produced J/psi and prompt decays of heavier charmonium states such as the Psi(2S) and Csi_c resonances, is sigma_prompt-J/psi(pt > 1.3 GeV/c, |y| < 0.9) = 8.3 +- 0.8(stat.) +- 1.1(syst.) + 1.5 - 1.4(syst. pol.) micro barn. The cross section for the production of b-hadrons decaying to J/psi with p_t>1.3 GeV/c and |y|<0.9 is sigma_{J/psi<-h_B} = 1.46 +- 0.38(stat.) + 0.26 -0.32(syst.) micro barn. The results are compared to QCD model predictions. The shape of the p_t and y distributions of b-quarks predicted by perturbative QCD model calculations are used to extrapolate the measured cross section to derive the b-bbar pair total cross section and dsigma/dy at mid-rapidity.

Inclusive J/$\psi$ production has been studied with the ALICE detector in p-Pb collisions at the nucleon-nucleon center of mass energy $\sqrt{s_{\rm NN}}$ = 5.02 TeV at the CERN LHC. The measurement is performed in the center of mass rapidity domains $2.03<y_{\rm cms}<3.53$ and $-4.46<y_{\rm cms}<-2.96$, down to zero transverse momentum, studying the $\mu^+\mu^-$ decay mode. In this paper, the J/$\psi$ production cross section and the nuclear modification factor $R_{\rm pPb}$ for the rapidities under study are presented. While at forward rapidity, corresponding to the proton direction, a suppression of the J/$\psi$ yield with respect to binary-scaled pp collisions is observed, in the backward region no suppression is present. The ratio of the forward and backward yields is also measured differentially in rapidity and transverse momentum. Theoretical predictions based on nuclear shadowing, as well as on models including, in addition, a contribution from partonic energy loss, are in fair agreement with the experimental results.