Model-independent unbinned analysis of B to K^*(to K^+π^-)μ^+μ^-: zeroes, bounds, Wilson coefficients and symmetries

We present a model-independent method to study the four-body decay $B\to K^*(\to K^+\pi^-)\mu^+\mu^-$, based on extracting continuous observables with a moments approach. The method allows the observables to be determined unbinned in both the dilepton and $K^+\pi^-$ invariant masses on which the decay dynamics depend. This will allow the method to shed new light on how the observables depend on the P- and S-wave contributions to the $K^+\pi^-$ system. This approach contrasts with the state-of-the-art analyses, which bin in dilepton and $K^+\pi^-$ mass, or use a model for the dependence of the underlying decay amplitudes on these masses. The method does not require making a statistical fit, and so avoids problems of biases and poor uncertainty estimation when dealing with small samples or a large number of fit parameters. We provide the Standard Model predictions for the unbinned optimised observables, derive new geometrical bounds on their values and study the robustness of these bounds in the presence of a scalar new physics contribution. We explore the zero-crossing points of $P_2$ and $P_{4,5}^\prime$ observables as a function of a new physics contribution to the dominant vector Wilson coefficient, $C_9^{\rm NP}$. We also discuss the conditions that can be used to test the theoretical model of the amplitudes needed for an experimental amplitude analysis. Finally, as an illustration, we show how the proposed method might be used to extract the zero-crossing points, make a comparison with the bounds and test a non-trivial relation between the observable values.