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.
$\mathrm{J}\mspace{-2mu}/\mspace{-2mu}\psi \to \mu\mu\mu\mu$ branching fraction
$\mathcal{B}(\mathrm{J}\mspace{-2mu}/\mspace{-2mu}\psi \to \mu\mu\mu\mu)$ / $\mathcal{B}(\mathrm{J}\mspace{-2mu}/\mspace{-2mu}\psi \to \mu\mu)$ ratio
Using proton-proton collision data corresponding to an integrated luminosity of 140 fb$^{-1}$ collected by the CMS experiment at $\sqrt{s}$ = 13 TeV, the $\Lambda_\text{b}^0$$\to$ J/$\psi\Xi^-$K$^+$ decay is observed for the first time, with a statistical significance exceeding 5 standard deviations. The relative branching fraction, with respect to the $\Lambda_\text{b}^0$$\to$$\psi$(2S)$\Lambda$ decay, is measured to be $\mathcal{B}$($\Lambda_\text{b}^0$$\to$ J/$\psi\Xi^-$K$^+$) / $\mathcal{B}$( $\Lambda_\text{b}^0$$\to$$\psi$(2S)$\Lambda$) = [3.38 $\pm$ 1.02 $\pm$ 0.61 $\pm$ 0.03]%, where the first uncertainty is statistical, the second is systematic, and the third is related to the uncertainties in $\mathcal{B}$($\psi$(2S) $\to$ J/$\psi\pi^+\pi^-$) and $\mathcal{B}$($\Xi^-$ $\to$ $\Lambda\pi^-$).
The measured branching fraction ratio
A search for the rare $\eta$$\to$$\mu^+\mu^-\mu^+\mu^-$ double-Dalitz decay is performed using a sample of proton-proton collisions, collected by the CMS experiment at the CERN LHC with high-rate muon triggers in 2017-2018 and corresponding to an integrated luminosity of 101 fb$^{-1}$. A signal having a statistical significance well in excess of 5 standard deviations is observed. Using the $\eta$$\to$$\mu^+ \mu^-$ decay as normalization, the branching fraction $\mathcal{B}(\eta$$\to$$\mu^+\mu^-\mu^+\mu^-)$ = [5.0 $\pm$ 0.8 (stat) $\pm$ 0.7 (syst) $\pm$ 0.7 ($\mathcal{B}_{2\mu}$)] $\times$ 10$^{-9}$ is measured, where the last term is the uncertainty in the normalization channel branching fraction. This work achieves an improved precision of over five orders of magnitude compared to previous results, leading to the first measurement of this branching fraction, which is found to agree with theoretical predictions.
The total efficiencies for the four-muon ($A_{4\mu}^{i,j}$, red and blue points) and two-muon ($A_{2\mu}^{i,j}$, orange and green points) decay channels, as functions of the generated meson's $p_{\mathrm{T}}$ and $y$, evaluated through MC simulation.
Measured ratio of $\mathcal{B}_{4\mu}/\mathcal{B}_{2\mu}$
Measured branching fraction $\mathcal{B}_{4\mu}$