Module flavio.physics.zdecays.test_zdecays
Classes
class TestAFBZ (methodName='runTest')-
Expand source code
class TestAFBZ(unittest.TestCase): def test_afbz_sm(self): for l in ['e', 'mu', 'tau']: self.assertAlmostEqual(flavio.sm_prediction('A(Z->{}{})'.format(l, l)), 0.1477, delta=0.0002, msg="Failed for {}".format(l)) self.assertAlmostEqual(flavio.sm_prediction('AFB(Z->{}{})'.format(l, l)), 0.0163, delta=0.0002, msg="Failed for {}".format(l)) self.assertAlmostEqual(flavio.sm_prediction('A(Z->bb)'), 0.935, delta=0.001) self.assertAlmostEqual(flavio.sm_prediction('A(Z->cc)'), 0.668, delta=0.001) self.assertAlmostEqual(flavio.sm_prediction('A(Z->ss)'), 0.935, delta=0.001) self.assertAlmostEqual(flavio.sm_prediction('AFB(Z->bb)'), 0.1035, delta=0.0002) self.assertAlmostEqual(flavio.sm_prediction('AFB(Z->cc)'), 0.0738, delta=0.0002)A class whose instances are single test cases.
By default, the test code itself should be placed in a method named 'runTest'.
If the fixture may be used for many test cases, create as many test methods as are needed. When instantiating such a TestCase subclass, specify in the constructor arguments the name of the test method that the instance is to execute.
Test authors should subclass TestCase for their own tests. Construction and deconstruction of the test's environment ('fixture') can be implemented by overriding the 'setUp' and 'tearDown' methods respectively.
If it is necessary to override the init method, the base class init method must always be called. It is important that subclasses should not change the signature of their init method, since instances of the classes are instantiated automatically by parts of the framework in order to be run.
When subclassing TestCase, you can set these attributes: * failureException: determines which exception will be raised when the instance's assertion methods fail; test methods raising this exception will be deemed to have 'failed' rather than 'errored'. * longMessage: determines whether long messages (including repr of objects used in assert methods) will be printed on failure in addition to any explicit message passed. * maxDiff: sets the maximum length of a diff in failure messages by assert methods using difflib. It is looked up as an instance attribute so can be configured by individual tests if required.
Create an instance of the class that will use the named test method when executed. Raises a ValueError if the instance does not have a method with the specified name.
Ancestors
- unittest.case.TestCase
Methods
def test_afbz_sm(self)-
Expand source code
def test_afbz_sm(self): for l in ['e', 'mu', 'tau']: self.assertAlmostEqual(flavio.sm_prediction('A(Z->{}{})'.format(l, l)), 0.1477, delta=0.0002, msg="Failed for {}".format(l)) self.assertAlmostEqual(flavio.sm_prediction('AFB(Z->{}{})'.format(l, l)), 0.0163, delta=0.0002, msg="Failed for {}".format(l)) self.assertAlmostEqual(flavio.sm_prediction('A(Z->bb)'), 0.935, delta=0.001) self.assertAlmostEqual(flavio.sm_prediction('A(Z->cc)'), 0.668, delta=0.001) self.assertAlmostEqual(flavio.sm_prediction('A(Z->ss)'), 0.935, delta=0.001) self.assertAlmostEqual(flavio.sm_prediction('AFB(Z->bb)'), 0.1035, delta=0.0002) self.assertAlmostEqual(flavio.sm_prediction('AFB(Z->cc)'), 0.0738, delta=0.0002)
class TestGammaZ (methodName='runTest')-
Expand source code
class TestGammaZ(unittest.TestCase): def test_obs_sm(self): # check the SM predictions self.assertAlmostEqual(flavio.sm_prediction('GammaZ'), 2.4950, delta=0.0015) self.assertAlmostEqual(flavio.sm_prediction('GammaZ'), 1 / par['tau_Z'], delta=0.0015) self.assertAlmostEqual(flavio.sm_prediction('sigma_had') / pb / 1e3, 41.488, delta=0.05) self.assertAlmostEqual(flavio.sm_prediction('Gamma(Z->ee)'), 83.966e-3, delta=0.02e-3) self.assertAlmostEqual(flavio.sm_prediction('Gamma(Z->mumu)'), 83.966e-3, delta=0.02e-3) self.assertAlmostEqual(flavio.sm_prediction('Gamma(Z->tautau)'), 83.776e-3, delta=0.02e-3) self.assertAlmostEqual(flavio.sm_prediction('Gamma(Z->uu)'), 299.936e-3, delta=0.6e-3) self.assertAlmostEqual(flavio.sm_prediction('Gamma(Z->cc)'), 299.860e-3, delta=0.6e-3) self.assertAlmostEqual(flavio.sm_prediction('Gamma(Z->dd)'), 382.770e-3, delta=0.4e-3) self.assertAlmostEqual(flavio.sm_prediction('Gamma(Z->ss)'), 382.770e-3, delta=0.4e-3) self.assertAlmostEqual(flavio.sm_prediction('Gamma(Z->bb)'), 375.724e-3, delta=0.4e-3) self.assertAlmostEqual(flavio.sm_prediction('Gamma(Z->nunu)'), 167.157e-3, delta=0.03e-3) self.assertAlmostEqual(flavio.sm_prediction('R_l'), 20750.9e-3, delta=30e-3) self.assertAlmostEqual(flavio.sm_prediction('R_c'), 172.23e-3, delta=0.09e-3) self.assertAlmostEqual(flavio.sm_prediction('R_b'), 215.80e-3, delta=0.06e-3) self.assertAlmostEqual(flavio.sm_prediction('R_e'), 20.743, delta=0.01) self.assertAlmostEqual(flavio.sm_prediction('R_mu'), 20.743, delta=0.01) self.assertAlmostEqual(flavio.sm_prediction('R_tau'), 20.743, delta=0.05) self.assertAlmostEqual(flavio.sm_prediction('R_uc'), 0.1724, delta=0.0002) self.assertAlmostEqual(flavio.sm_prediction('R_uc'), (flavio.sm_prediction('R_u') + flavio.sm_prediction('R_c')) / 2, delta=1e-15) def test_r_sm(self): # check that the Sm predictions for the Ri agree with the Gammas par = flavio.default_parameters.get_central_all() mh = par['m_h'] mt = par['m_t'] als = par['alpha_s'] Da = 0.059 mZ = par['m_Z'] arg = (mh, mt, als, Da, mZ) Rl = Zobs('Rl', *arg) Rc = Zobs('Rc', *arg) Rb = Zobs('Rb', *arg) Ge = Zobs('Gammae,mu', *arg) Gmu = Zobs('Gammae,mu', *arg) Gtau = Zobs('Gammatau', *arg) Gu = Zobs('Gammau', *arg) Gd = Zobs('Gammad,s', *arg) Gs = Zobs('Gammad,s', *arg) Gc = Zobs('Gammac', *arg) Gb = Zobs('Gammab', *arg) Ghad = Gu + Gd + Gc + Gs + Gb Gl = (Ge + Gmu + Gtau) / 3. self.assertAlmostEqual(Rl, Ghad / Gl, delta=1e-4) self.assertAlmostEqual(Rc, Gc / Ghad, delta=1e-4) self.assertAlmostEqual(Rb, Gb / Ghad, delta=1e-4) def test_obs_sm_fv(self): # check the SM predictions for LFV decays self.assertEqual(flavio.sm_prediction('BR(Z->emu)'), 0) self.assertEqual(flavio.sm_prediction('BR(Z->etau)'), 0) self.assertEqual(flavio.sm_prediction('BR(Z->mutau)'), 0) def test_Gamma_NP(self): # compare NP contributions to A.49-A.52 from 1706.08945 GF, mZ, s2w_eff = par['GF'], par['m_Z'], par['s2w']*1.0010 d_gV = 0.055 d_gA = 0.066 # A.49-A.52 from 1706.08945 dGamma_Zll = sqrt(2)*GF*mZ**3/(6*pi) * (-d_gA + (-1+4*s2w_eff)*d_gV) dGamma_Znn = sqrt(2)*GF*mZ**3/(6*pi) * (d_gA + d_gV) dGamma_Zuu = sqrt(2)*GF*mZ**3/(pi) * (d_gA -1/3*(-3+8*s2w_eff)*d_gV) /2 dGamma_Zdd = sqrt(2)*GF*mZ**3/(pi) * (-3/2*d_gA +1/2*(-3+4*s2w_eff)*d_gV) /3 # term squared in d_gV and d_gA not included in 1706.08945 d_g_squared = sqrt(2)*GF*mZ**3/(3*pi)*(abs(d_gV)**2+abs(d_gA)**2) self.assertAlmostEqual( dGamma_Zll + d_g_squared, GammaZ_NP(par, 1, smeftew.gV_SM('e', par), d_gV, smeftew.gA_SM('e', par), d_gA) ) self.assertAlmostEqual( dGamma_Znn + d_g_squared, GammaZ_NP(par, 1, smeftew.gV_SM('nue', par), d_gV, smeftew.gA_SM('nue', par), d_gA) ) self.assertAlmostEqual( dGamma_Zuu + 3*d_g_squared, GammaZ_NP(par, 3, smeftew.gV_SM('u', par), d_gV, smeftew.gA_SM('u', par), d_gA) ) self.assertAlmostEqual( dGamma_Zdd + 3*d_g_squared, GammaZ_NP(par, 3, smeftew.gV_SM('d', par), d_gV, smeftew.gA_SM('d', par), d_gA) )A class whose instances are single test cases.
By default, the test code itself should be placed in a method named 'runTest'.
If the fixture may be used for many test cases, create as many test methods as are needed. When instantiating such a TestCase subclass, specify in the constructor arguments the name of the test method that the instance is to execute.
Test authors should subclass TestCase for their own tests. Construction and deconstruction of the test's environment ('fixture') can be implemented by overriding the 'setUp' and 'tearDown' methods respectively.
If it is necessary to override the init method, the base class init method must always be called. It is important that subclasses should not change the signature of their init method, since instances of the classes are instantiated automatically by parts of the framework in order to be run.
When subclassing TestCase, you can set these attributes: * failureException: determines which exception will be raised when the instance's assertion methods fail; test methods raising this exception will be deemed to have 'failed' rather than 'errored'. * longMessage: determines whether long messages (including repr of objects used in assert methods) will be printed on failure in addition to any explicit message passed. * maxDiff: sets the maximum length of a diff in failure messages by assert methods using difflib. It is looked up as an instance attribute so can be configured by individual tests if required.
Create an instance of the class that will use the named test method when executed. Raises a ValueError if the instance does not have a method with the specified name.
Ancestors
- unittest.case.TestCase
Methods
def test_Gamma_NP(self)-
Expand source code
def test_Gamma_NP(self): # compare NP contributions to A.49-A.52 from 1706.08945 GF, mZ, s2w_eff = par['GF'], par['m_Z'], par['s2w']*1.0010 d_gV = 0.055 d_gA = 0.066 # A.49-A.52 from 1706.08945 dGamma_Zll = sqrt(2)*GF*mZ**3/(6*pi) * (-d_gA + (-1+4*s2w_eff)*d_gV) dGamma_Znn = sqrt(2)*GF*mZ**3/(6*pi) * (d_gA + d_gV) dGamma_Zuu = sqrt(2)*GF*mZ**3/(pi) * (d_gA -1/3*(-3+8*s2w_eff)*d_gV) /2 dGamma_Zdd = sqrt(2)*GF*mZ**3/(pi) * (-3/2*d_gA +1/2*(-3+4*s2w_eff)*d_gV) /3 # term squared in d_gV and d_gA not included in 1706.08945 d_g_squared = sqrt(2)*GF*mZ**3/(3*pi)*(abs(d_gV)**2+abs(d_gA)**2) self.assertAlmostEqual( dGamma_Zll + d_g_squared, GammaZ_NP(par, 1, smeftew.gV_SM('e', par), d_gV, smeftew.gA_SM('e', par), d_gA) ) self.assertAlmostEqual( dGamma_Znn + d_g_squared, GammaZ_NP(par, 1, smeftew.gV_SM('nue', par), d_gV, smeftew.gA_SM('nue', par), d_gA) ) self.assertAlmostEqual( dGamma_Zuu + 3*d_g_squared, GammaZ_NP(par, 3, smeftew.gV_SM('u', par), d_gV, smeftew.gA_SM('u', par), d_gA) ) self.assertAlmostEqual( dGamma_Zdd + 3*d_g_squared, GammaZ_NP(par, 3, smeftew.gV_SM('d', par), d_gV, smeftew.gA_SM('d', par), d_gA) ) def test_obs_sm(self)-
Expand source code
def test_obs_sm(self): # check the SM predictions self.assertAlmostEqual(flavio.sm_prediction('GammaZ'), 2.4950, delta=0.0015) self.assertAlmostEqual(flavio.sm_prediction('GammaZ'), 1 / par['tau_Z'], delta=0.0015) self.assertAlmostEqual(flavio.sm_prediction('sigma_had') / pb / 1e3, 41.488, delta=0.05) self.assertAlmostEqual(flavio.sm_prediction('Gamma(Z->ee)'), 83.966e-3, delta=0.02e-3) self.assertAlmostEqual(flavio.sm_prediction('Gamma(Z->mumu)'), 83.966e-3, delta=0.02e-3) self.assertAlmostEqual(flavio.sm_prediction('Gamma(Z->tautau)'), 83.776e-3, delta=0.02e-3) self.assertAlmostEqual(flavio.sm_prediction('Gamma(Z->uu)'), 299.936e-3, delta=0.6e-3) self.assertAlmostEqual(flavio.sm_prediction('Gamma(Z->cc)'), 299.860e-3, delta=0.6e-3) self.assertAlmostEqual(flavio.sm_prediction('Gamma(Z->dd)'), 382.770e-3, delta=0.4e-3) self.assertAlmostEqual(flavio.sm_prediction('Gamma(Z->ss)'), 382.770e-3, delta=0.4e-3) self.assertAlmostEqual(flavio.sm_prediction('Gamma(Z->bb)'), 375.724e-3, delta=0.4e-3) self.assertAlmostEqual(flavio.sm_prediction('Gamma(Z->nunu)'), 167.157e-3, delta=0.03e-3) self.assertAlmostEqual(flavio.sm_prediction('R_l'), 20750.9e-3, delta=30e-3) self.assertAlmostEqual(flavio.sm_prediction('R_c'), 172.23e-3, delta=0.09e-3) self.assertAlmostEqual(flavio.sm_prediction('R_b'), 215.80e-3, delta=0.06e-3) self.assertAlmostEqual(flavio.sm_prediction('R_e'), 20.743, delta=0.01) self.assertAlmostEqual(flavio.sm_prediction('R_mu'), 20.743, delta=0.01) self.assertAlmostEqual(flavio.sm_prediction('R_tau'), 20.743, delta=0.05) self.assertAlmostEqual(flavio.sm_prediction('R_uc'), 0.1724, delta=0.0002) self.assertAlmostEqual(flavio.sm_prediction('R_uc'), (flavio.sm_prediction('R_u') + flavio.sm_prediction('R_c')) / 2, delta=1e-15) def test_obs_sm_fv(self)-
Expand source code
def test_obs_sm_fv(self): # check the SM predictions for LFV decays self.assertEqual(flavio.sm_prediction('BR(Z->emu)'), 0) self.assertEqual(flavio.sm_prediction('BR(Z->etau)'), 0) self.assertEqual(flavio.sm_prediction('BR(Z->mutau)'), 0) def test_r_sm(self)-
Expand source code
def test_r_sm(self): # check that the Sm predictions for the Ri agree with the Gammas par = flavio.default_parameters.get_central_all() mh = par['m_h'] mt = par['m_t'] als = par['alpha_s'] Da = 0.059 mZ = par['m_Z'] arg = (mh, mt, als, Da, mZ) Rl = Zobs('Rl', *arg) Rc = Zobs('Rc', *arg) Rb = Zobs('Rb', *arg) Ge = Zobs('Gammae,mu', *arg) Gmu = Zobs('Gammae,mu', *arg) Gtau = Zobs('Gammatau', *arg) Gu = Zobs('Gammau', *arg) Gd = Zobs('Gammad,s', *arg) Gs = Zobs('Gammad,s', *arg) Gc = Zobs('Gammac', *arg) Gb = Zobs('Gammab', *arg) Ghad = Gu + Gd + Gc + Gs + Gb Gl = (Ge + Gmu + Gtau) / 3. self.assertAlmostEqual(Rl, Ghad / Gl, delta=1e-4) self.assertAlmostEqual(Rc, Gc / Ghad, delta=1e-4) self.assertAlmostEqual(Rb, Gb / Ghad, delta=1e-4)