Prompt dimuon mass spectra showing the (black) inclusive and (red) $X+b$ candidates with all fiducial and selection requirements applied. The grey boxes show the regions vetoed due to large contributions from QCD resonances.

Displaced dimuon mass spectra showing the (black) inclusive and (red) promptly produced candidates with all fiducial and selection requirements applied. The grey box shows the region vetoed due to the large doubly misidentified $ K ^0_{\mathrm{S}}$ background, whose low-mass tail extends into the search region. A dedicated study of this region is presented in a search for $ K ^0_{\mathrm{S}} \rightarrow \mu^+\mu^-$ decays [44].

Signed local significances in the $ m( X ) < 20\text{ Ge V} $ region for the (top) inclusive and (bottom) associated beauty searches for prompt $ X \rightarrow \mu^+\mu^-$ decays. If the best-fit signal-yield estimator is negative, the signed significance is negative and { vice versa}. The grey regions are excluded either due to a nearby large QCD resonance contribution, or because the overlap of the bin with the fiducial region in Table 1 is small.

Signed local significances in the $ m( X ) > 20\text{ Ge V} $ region for the (top) inclusive and (bottom) associated beauty searches for prompt $ X \rightarrow \mu^+\mu^-$ decays. The lower limit on the vertical axis of $\log_{10}[\Gamma( X ) / \text{ Me V} ] = -\infty$ corresponds to $\Gamma( X ) = 0$.

Fit to the $ m(\mu^+\mu^-)$ spectrum in events with at least two $b$-tagged jets. The ${27 < m( X ) < 30\text{ Ge V} }$ search region is marked by the vertical dashed lines.

Signed local significances for the (top) promptly produced and (bottom) inclusive searches for displaced $ X \rightarrow \mu^+\mu^-$ decays. The black points show the individual candidates.

Upper limits at 90% confidence level on the cross section $\sigma( X \rightarrow \mu^+\mu^- )$ in the $ m( X ) < 20\text{ Ge V} $ region for the (top) inclusive and (bottom) associated beauty searches for prompt $ X \rightarrow \mu^+\mu^-$ decays.

Upper limits at 90% confidence level on the cross section $\sigma( X \rightarrow \mu^+\mu^- )$ in the $ m( X ) > 20\text{ Ge V} $ region for the (top) inclusive and (bottom) associated beauty searches for prompt $ X \rightarrow \mu^+\mu^-$ decays.

Upper limits at 90% confidence level on the cross section $\sigma( X \rightarrow \mu^+\mu^- )$ for the (top) promptly produced and (bottom) inclusive searches for displaced $ X \rightarrow \mu^+\mu^-$ decays.

Upper limits at 90% confidence level on the $X$--$H$ mixing angle, $\theta_H$, for the 2HDM scenario discussed in the text (blue) from this analysis compared with existing limits from (red) {\sc BaBar} [58], (green) CMS Run 1 [59], (magenta) CMS Run 2 [60] and (yellow) LHCb Run 1 [61].

Upper limits at 90% confidence level on the $\gamma$--$Z_{\rm HV}$ kinetic mixing strength for the HV scenario discussed in the text. Here, $X$ denotes a composite HV vector boson.

Animated gif made out of all figures.

Fiducial regions of the searches for prompt and displaced $ X \rightarrow \mu^+\mu^-$ decays.

Summary of systematic uncertainties. The luminosity and efficiency uncertainties are highly correlated, which is accounted for when obtaining the total uncertainties.