def test_smpar_roundtrip(self):
wc = get_random_wc('SMEFT', 'Warsaw', cmax=1e-6)
smeft = SMEFT(wc=None)
smeft.scale_in = 160
smeft._set_initial_wcxf(wc, get_smpar=False)
CSM = smpar.smeftpar(smeft.scale_in, smeft.C_in, 'Warsaw')
Cboth = CSM.copy()
Cboth.update(smeft.C_in)
Cback = smpar.smpar(Cboth)
for k in smpar.p:
if k in ['m_Z']:
self.assertAlmostEqual(smpar.p[k] / Cback[k],
1,
msg="Failed for {}".format(k),
delta=0.05)
elif k in ['delta']:
self.assertAlmostEqual(smpar.p[k] / Cback[k],
1,
msg="Failed for {}".format(k),
delta=1e-3)
else:
self.assertAlmostEqual(smpar.p[k] / Cback[k],
1,
msg="Failed for {}".format(k),
delta=1e-5)
def test_wcxf_smpar(self):
wc = get_random_wc('SMEFT', 'Warsaw', 1e5, 1e-11)
smeft = SMEFT(wc)
C_out = smeft._rgevolve(91.1876)
p_out = smpar.smpar(C_out)
for k in p_out:
self.assertAlmostEqual(p_out[k] / smpar.p[k], 1,
delta=0.05,
msg="Failed for {}".format(k))
def test_smpar_small(self):
wc = get_random_wc('SMEFT', 'Warsaw', cmax=1e-24)
smeft = SMEFT(wc=None)
smeft.scale_in = 160
smeft._set_initial_wcxf(wc, get_smpar=False)
CSM = smpar.smeftpar(smeft.scale_in, smeft.C_in, 'Warsaw')
Cboth = CSM.copy()
Cboth.update(smeft.C_in)
Cback = smpar.smpar(Cboth)
for k in smpar.p:
if k not in ['m_Z', 'gamma']:
self.assertAlmostEqual(smpar.p[k], Cback[k],
msg="Failed for {}".format(k))
