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.

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 B^0_s differential production cross section is measured as functions of the transverse momentum and rapidity in pp collisions at sqrt(s) = 7 TeV, using the J/Psi phi decay, and compared with predictions based on perturbative QCD calculations at next-to-leading order. The data sample, collected by the CMS experiment at the LHC, corresponds to an integrated luminosity of 40 inverse picobarns. The B^0_s is reconstructed from the decays J/Psi to an oppositely charged muon pair and phi to K+ K-. The integrated B^0_s cross section times B^0_s to J/Psi phi branching fraction in the range 8 < pt(B) < 50 GeV/c and |y(b)| < 2.4 is measured to be 6.9 +/- 0.6 +/- 0.6 nb, where the first uncertainty is statistical and the second is systematic.

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 J/$\psi$$\to$$\mu^+\mu^-\mu^+\mu^-$ decay has been observed with a statistical significance in excess of five standard deviations. The analysis is based on an event sample of proton-proton collisions at a center-of-mass energy of 13 TeV, collected by the CMS experiment in 2018 and corresponding to an integrated luminosity of 33.6 fb${-1}$. Normalizing to the J/$\psi$$\to$$\mu^+\mu^-$ decay mode leads to a branching fraction [10.1$^{+3.3}_{-2.7}$ (stat) $\pm$ 0.4 (syst) ]$\times$ 10$^{-7}$, a value that is consistent with the standard model prediction.

Protons consist of three valence quarks, two up-quarks and one down-quark, held together by gluons and a sea of quark-antiquark pairs. Collectively, quarks and gluons are referred to as partons. In a proton-proton collision, typically only one parton of each proton undergoes a hard scattering - referred to as single-parton scattering - leaving the remainder of each proton only slightly disturbed. Here, we report the study of double- and triple-parton scatterings through the simultaneous production of three J/$\psi$ mesons, which consist of a charm quark-antiquark pair, in proton-proton collisions recorded with the CMS experiment at the Large Hadron Collider. We observed this process - reconstructed through the decays of J/$\psi$ mesons into pairs of oppositely charged muons - with a statistical significance above five standard deviations. We measured the inclusive fiducial cross section to be 272 $^{+141}_{-104}$ (stat) $\pm$ 17 (syst) fb, and compared it to theoretical expectations for triple-J/$\psi$ meson production in single-, double- and triple-parton scattering scenarios. Assuming factorization of multiple hard-scattering probabilities in terms of single-parton scattering cross sections, double- and triple-parton scattering are the dominant contributions for the measured process.

The $B_\mathrm{c}^+$ meson is observed for the first time in heavy ion collisions. Data from the CMS detector are used to study the production of the $B_\mathrm{c}^+$ meson in lead-lead (PbPb) and proton-proton (pp) collisions at a center-of-mass energy per nucleon pair of $\sqrt{s_{\mathrm{NN}}} =$ 5.02 TeV, via the $B_\mathrm{c}^+ \to (J/\psi\to\mu^+\mu^-)\mu^+\nu_\mu$ decay. The $B_\mathrm{c}^+$ nuclear modification factor, derived from the PbPb-to-pp ratio of production cross sections, is measured in two bins of the trimuon transverse momentum and of the PbPb collision centrality. The B$_\mathrm{c}^+$ meson is shown to be less suppressed than quarkonia and most of the open heavy-flavor mesons, suggesting that effects of the hot and dense nuclear matter created in heavy ion collisions contribute to its production. This measurement sets forth a promising new probe of the interplay of suppression and enhancement mechanisms in the production of heavy-flavor mesons in the quark-gluon plasma.

Yields for J/psi production in Cu+Cu collisions at sqrt (s_NN)= 200 GeV have been measured by the PHENIX experiment over the rapidity range |y| < 2.2 at transverse momenta from 0 to beyond 5 GeV/c. The invariant yield is obtained as a function of rapidity, transverse momentum and collision centrality, and compared with results in p+p and Au+Au collisions at the same energy. The Cu+Cu data provide greatly improved precision over existing Au+Au data for J/psi production in collisions with small to intermediate numbers of participants, providing a key constraint that is needed for disentangling cold and hot nuclear matter effects.