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.

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.

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.

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

A first measurement of the top quark mass using the decay channel t to (W to l nu) (b to J/psi + X to mu+ mu- + X) is presented. The analysis uses events selected from the proton-proton collisions recorded with the CMS detector at the LHC at a center-of-mass energy of 8 TeV. The data correspond to an integrated luminosity of 19.7 inverse femtobarns, with 666 t t-bar and single top quark candidate events containing a reconstructed J/psi candidate decaying into an oppositely-charged muon pair. The mass of the (J/psi + l) system, where l is an electron or a muon from W boson decay, is used to extract a top quark mass of 173.5 +/- 3.0 (stat) +/- 0.9 (syst) GeV.

The production of the X(3872) is studied in pp collisions at sqrt(s) = 7 TeV, using decays to J/psi pi pi, where the J/psi decays to two muons. The data were recorded by the CMS experiment and correspond to an integrated luminosity of 4.8 inverse femtobarns. The measurements are performed in a kinematic range in which the X(3872) candidates have a transverse momentum 10 < pt < 50 GeV and rapidity abs(y) < 1.2. The ratio of the X(3872) and psi(2S) cross sections times their branching fractions into J/psi pi pi is measured as a function of pt. In addition, the fraction of X(3872) originating from B decays is determined. From these measurements the prompt X(3872) differential cross section times branching fraction as a function of pt is extracted. The pi pi mass spectrum of the J/psi pi pi system in the X(3872) decays is also investigated.

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.

The Lambda(b) differential production cross section and the cross section ratio anti-Lambda(b)/Lambda(b) are measured as functions of transverse momentum pt(Lambda(b)) and rapidity abs(y(Lambda(b))) in pp collisions at sqrt(s) = 7 TeV using data collected by the CMS experiment at the LHC. The measurements are based on Lambda(b) decays reconstructed in the exclusive final state J/Psi Lambda, with the subsequent decays J/Psi to an opposite-sign muon pair and Lambda to proton pion, using a data sample corresponding to an integrated luminosity of 1.9 inverse femtobarns. The product of the cross section times the branching ratio for Lambda(b) to J/Psi Lambda versus pt(Lambda(b)) falls faster than that of b mesons. The measured value of the cross section times the branching ratio for pt(Lambda(b)) > 10 GeV and abs(y(Lambda(b))) < 2.0 is 1.06 +/- 0.06 +/- 0.12 nb, and the integrated cross section ratio for anti-Lambda(b)/Lambda(b) is 1.02 +/- 0.07 +/- 0.09, where the uncertainties are statistical and systematic, respectively.