import flavio
from wilson import Wilson
with open('SN.yaml', 'r') as f:
    my_w = Wilson.load_wc(f)
print(
    "R_K : ", flavio.np_prediction('<Rmue>(B0->Kll)', my_w, q2min=1, q2max=6), "\n",
    "BR(B -> K mu tau) : ", flavio.np_prediction('BR(B+->Kmutau,taumu)', my_w), "\n",
    "Delta M_s / (Delta M_s)_SM : ", flavio.np_prediction('DeltaM_s', my_w)/flavio.sm_prediction('DeltaM_s'), "\n",
    "delta a_mu : ", flavio.np_prediction('a_mu', my_w)-flavio.sm_prediction('a_mu'), "\n",
    "BR(mu -> e gamma) : ", flavio.np_prediction('BR(mu->egamma)', my_w), "\n"
)