As for the flavour sector, there are numerous observables that have to be computed and checked for each parameter space point that satisfies the scalar sector constraints. The inverted procedure outlined in sec:Yukawa already ensures that our choice of the parameter-space points, even before a numerical scan, automatically reproduces the correct measured values of the quark masses and of the CKM matrix elements. But the presence of charged scalars, as well as neutral ones
with tree-level CKM-suppressed FCNC interactions, means that we must verify particularly sensitive quantities, such as the \(b\rightarrow s \gamma \)  width, the neutral \(K\) - and \(B\) -meson mass differences and the CP-violating \(\epsilon _K\)  phase, among others. To this end, we employ the FlavourKit package [1]} as part of the SPheno tool in order to compute the full list of Wilson coefficients for a given parameter-space point. The Wilson coefficients
are then passed over to the Flavio python package [2]} enabling to compute all the relevant flavour physics observables for a given point and to compare them to the corresponding SM values.
