Invariant mass distribution for $ B ^0_s \rightarrow \mu ^+\mu ^- K^+K^-$ candidates, with the mass of the $\mu^+\mu^-$ pair constrained to the nominal $J/\psi$ mass. The result of the fit is shown with signal (dashed curve) and combinatorial background (dotted curve) components and their sum (solid curve).

Background subtracted $ K ^+ K ^- $ invariant mass distribution for $ B ^0_s \rightarrow { J \mskip -3mu/\mskip -2mu\psi \mskip 2mu} K ^+ K ^- $ candidates. The vertical dotted lines separate the four intervals.

Distribution of (a) $ K ^+ K ^- $ S-wave signal events, and (b) $ K ^+ K ^- $ P-wave signal events, both in four invariant mass intervals. In (b), the distribution of simulated $ B ^0_s \rightarrow { J \mskip -3mu/\mskip -2mu\psi \mskip 2mu} \phi$ events in the four intervals assuming the same total number of P-wave events is also shown (dashed lines). Note the interference between the $ K ^+ K ^- $ S-wave and P-wave amplitudes integrated over the angular variables has vanishing contribution in these distributions.

Measured phase differences between S-wave and perpendicular P-wave amplitudes in four intervals of $m_{KK}$ for solution I (full blue circles) and solution II (full black squares). The asymmetric error bars correspond to $\Delta\ln L =-0.5$ (solid lines) and $\Delta\ln L =-2$ (dash-dotted lines).

Animated gif made out of all figures.

Numbers of signal and background events in the $m_{J/\psi KK}$ range of 5200--5550 MeV and statistical power per signal event in four intervals of $m_{KK}$.\\

Results from a simultaneous fit of the four intervals of $m_{KK}$, where the uncertainties are statistical only. Only parameters which are needed for the ambiguity resolution are shown.\\