A series of measurements from different experiments at the LHC along with the Belle experiment have reported deviations from their SM predictions of the angular observable \(P^{\prime }_5\)  in \(B \rightarrow K^* \, \mu ^+\,\mu ^-\)  decay at the level of 3\(\sigma \)  [1]}, [2]}, [3]}. In addition, the measured value of the branching ratio of \(B_s \rightarrow \phi \, \mu ^+\,\mu ^-\)  decay shows tension with the SM prediction at the level of 3.5\(\sigma \)  [4]} . These atypical measurements can be analysed in a model-independent framework using the language of effective field theory [5]}, [6]}, [7]}, [8]}, [9]}, [10]}, [11]}, [12]}, [13]}, [14]}, [15]}. Irrespective of the adopted methodology for the treatment of hadronic uncertainties, the global analysis of data favours new physics scenarios with non-zero value(s) of the muonic Wilson coefficients.
