Invariant $\mu ^+ \mu ^- $ versus $ K ^+ K ^- \mu ^+ \mu ^- $ mass. The charmonium vetoes are indicated by the solid lines. The vertical dashed lines indicate the signal region of $\pm50 {\mathrm{ Me V /}c^2} $ around the known $ B ^0_ s $ mass in which the signal decay $ B ^0_ s \rightarrow \phi\mu ^+ \mu ^- $ is visible.

Invariant mass of $ B ^0_ s \rightarrow \phi\mu ^+ \mu ^- $ candidates in six bins of invariant dimuon mass squared. The fitted signal component is denoted by the light blue shaded area, the combinatorial background component by the dark red shaded area. The solid line indicates the sum of the signal and background components.

Differential branching fraction ${\rm d}{\cal B}( B ^0_ s \rightarrow \phi\mu ^+ \mu ^- )/{\rm d}q^2$. Error bars include both statistical and systematic uncertainties added in quadrature. Shaded areas indicate the vetoed regions containing the $ { J \mskip -3mu/\mskip -2mu\psi \mskip 2mu} $ and $\psi {(2S)} $ resonances. The solid curve shows the leading order SM prediction, scaled to the fitted total branching fraction. The prediction uses the SM Wilson coefficients and leading order amplitudes given in Ref. [2], as well as the form factor calculations in Ref. [17]. $ B ^0_ s $ mixing is included as described in Ref. [1]. No error band is given for the theory prediction. The dashed curve denotes the leading order prediction scaled to a total branching fraction of $16\times 10^{-7}$ [19].

a) Longitudinal polarisation fraction $F_{\rm L}$, b) $S_3$, c) $A_6$, and d) $A_9$ in six bins of $q^2$. Error bars include statistical and systematic uncertainties added in quadrature. The solid curves are the leading order SM predictions, using the Wilson coefficients and leading order amplitudes given in Ref. [2], as well as the form factor calculations in Ref. [17]. $ B ^0_ s $ mixing is included as described in Ref. [1]. No error band is given for the theory predictions.

Animated gif made out of all figures.

Signal yield and differential branching fraction ${\rm d}{\cal B}( B ^0_ s \rightarrow \phi\mu ^+ \mu ^- )/{\rm d}q^2$ in six bins of $q^2$. Results are also quoted for the region $1<q^2<6 {\mathrm{ Ge V /}c^2} $ where theoretical predictions are most reliable. The first uncertainty is statistical, the second systematic, and the third from the branching fraction of the normalisation channel.

Results for the angular observables $F_{\rm L}$, $S_3$, $A_6$, and $A_9$ in bins of $q^2$. The first uncertainty is statistical, the second systematic.