import flavio
from wilson import Wilson
w = Wilson({ 'egamma_12': 5.5e-13, 'egamma_21': 2.6e-15 }, scale=5, eft='WET', basis='JMS')
print(flavio.np_prediction('BR(mu->egamma)', w))