Possible strategy of discriminating 1D new physics solutions in $b\to s\ell\ell$ decays after Moriond 2021

The recent measurements of $R_K$, $B_s\to\mu^+\mu^-$, a set of CP-averaged angular observables for the $B^0\to K^{*0}\mu^+\mu^-$ decay, and its isospin partner $B^+\to K^{*+}\mu^+\mu^-$ by the LHCb Collaboration, consistently hint at lepton universality violation in the $b\to s\ell\ell$ transitions. The so-called $B$ anamolies can be best explained in five one-dimensional scenarios, i.e, $\delta C_9^{\mu}$, $\delta C_{10}^{\mu}$, $\delta C_L^{\mu}$, $\delta C_9^{\mu}=\delta C_{10}^{\mu\prime}$, and $\delta C_9^{\mu}=-\delta C_9^{\mu\prime}$, as demonstrated in recent model independent anlayses~\cite{Alok:2019ufo,Alguero:2021anc,Geng:2021nhg,Altmannshofer:2021qrr}. In this work we explore how these scenarios can be distinguished from each other. We show that the combinations of four angular asymmetries $A_i$~$(i=3,4,5,9)$ together with the ratio $R_6$ first proposed in~\cite{Jager:2014rwa} can discriminate the five new physics scenarios in proper intervals of $q^2$ and with future high-precision measurements.