We present a fully automated procedure providing an easy way to perform, systematically, phenomenological analyses in flavor physics for general BSM scenarios. This procedure relies on MARTY-1.4, is model independent and requires as input only the Lagrangian of the theory. Once the Lagrangian has been defined, tree-level and one-loop Wilson coefficients can be calculated symbolically by MARTY, from which flavor observables can be computed numerically by available software programs. We focus in particular on $b\to s\gamma$ and the recently updated $b\to s\ell^+\ell^-$ observables which are in tension with the SM, and present a general procedure to extract the relevant one-loop coefficients, such as $C_7$, $C_7^\prime$, $C_9$ and $C_{10}$.