The dimuon invariant-mass distribution in the 5-10% multiplicity class. Each dimuon is weighted with the inverse of the detector acceptance and efficiency as a function of the dimuon $p_{\rm T}$ and $y$. The data are fitted to a combination of double-sided Crystal Ball (CB2) functions for the J/$\psi$ and $\psi(2S)$ signals, and a Variable-Width Gaussian function for the background. Only statistical uncertainties are shown. The bottom panel shows the ratio between the data and the fit.