Module flavio.physics.running.test_running
Classes
class TestRunning (methodName='runTest')-
Expand source code
class TestRunning(unittest.TestCase): def test_alphae(self): # compare to alpha_em at m_tau as given on p. 4 of # http://pdg.lbl.gov/2015/reviews/rpp2014-rev-standard-model.pdf alpha_tau = get_alpha_e(par_pdg, 1.777) self.assertAlmostEqual(1/alpha_tau/133.465,1.,places=3) # check at thresholds mc = config['RGE thresholds']['mc'] self.assertEqual(get_alpha_e(par_pdg, mc), get_alpha_e(par_pdg, mc, nf_out=4)) self.assertEqual(get_alpha_e(par_pdg, mc, nf_out=3), get_alpha_e(par_pdg, mc, nf_out=4)) mb = config['RGE thresholds']['mb'] self.assertEqual(get_alpha_e(par_pdg, mb), get_alpha_e(par_pdg, mb, nf_out=5)) self.assertEqual(get_alpha_e(par_pdg, mb, nf_out=4), get_alpha_e(par_pdg, mb, nf_out=5)) mt = config['RGE thresholds']['mt'] self.assertEqual(get_alpha_e(par_pdg, mt), get_alpha_e(par_pdg, mt, nf_out=6)) self.assertEqual(get_alpha_e(par_pdg, mt, nf_out=5), get_alpha_e(par_pdg, mt, nf_out=6)) def test_alphas(self): alpha_b = get_alpha(par_noqed, 4.2) # compare to 3-loop alpha_s at 4.2 GeV according to RunDec self.assertAlmostEqual(alpha_b['alpha_s']/0.225911,1.,places=4) def test_runningmasses(self): # compare to RunDec np.testing.assert_almost_equal(get_mb(par, 120.)/2.79211, 1,decimal=2) np.testing.assert_almost_equal(get_mt(par, 120.)/167.225, 1,decimal=2) def test_polemasses(self): np.testing.assert_almost_equal(get_mb_pole(par, nl=2)/4.78248, 1,decimal=2) np.testing.assert_almost_equal(get_mc_pole(par, nl=2)/1.68375, 1,decimal=2) np.testing.assert_almost_equal(get_mb_pole(par, nl=3)/4.92987, 1,decimal=2) def test_ksmass(self): # compare to 1107.3100 # KS -> MSbar conversion self.assertAlmostEqual( flavio.physics.running.masses.mKS2mMS(4.55, 4, get_alpha(par_ks, 4.55)['alpha_s'], Mu=1, nl=2), 4.20, delta=0.01) self.assertAlmostEqual( flavio.physics.running.masses.mKS2mMS(1.15, 3, get_alpha(par_ks, 1.15)['alpha_s'], Mu=1, nl=2), 1.329, delta=0.01) # MSbar -> KS conversion self.assertAlmostEqual(get_mb_KS(par_ks, 1.), 4.553, delta=0.01) self.assertAlmostEqual(get_mc_KS(par_ks, 1.), 1.091, delta=0.13) # this is satisfied poorly! def test_mb1S(self): par = flavio.default_parameters.get_central_all() par['m_b'] = 4.18 par['alpha_s'] = 0.1185 alpha_s = get_alpha(par, par['m_b'], nf_out=5)['alpha_s'] mb1S = get_mb_1S(par, nl=3) self.assertAlmostEqual(mb1S, 4.67, delta=0.005) def test_f_perp(self): par = flavio.default_parameters.get_central_all() self.assertEqual(par['f_perp_K*0'], get_f_perp(par, 'K*0', 2)) self.assertEqual(par['f_perp_K*0'], get_f_perp(par, 'K*0', 2, nf_out=3)) self.assertEqual(par['f_perp_K*0'], get_f_perp(par, 'K*0', 2, nf_out=4)) # rough comparison of 1503.05534 and 0804.0473 self.assertAlmostEqual(get_f_perp(par, 'K*0', 1) / par['f_perp_K*0'], 0.159 / (0.712 * 0.204), delta=0.025) self.assertAlmostEqual(get_f_perp(par, 'K*0', 1) / par['f_perp_K*0'], 0.159 / (0.712 * 0.204), delta=0.025) # just check this doesn't raise get_f_perp(par, 'phi', 100, nf_out=6)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_alphae(self)-
Expand source code
def test_alphae(self): # compare to alpha_em at m_tau as given on p. 4 of # http://pdg.lbl.gov/2015/reviews/rpp2014-rev-standard-model.pdf alpha_tau = get_alpha_e(par_pdg, 1.777) self.assertAlmostEqual(1/alpha_tau/133.465,1.,places=3) # check at thresholds mc = config['RGE thresholds']['mc'] self.assertEqual(get_alpha_e(par_pdg, mc), get_alpha_e(par_pdg, mc, nf_out=4)) self.assertEqual(get_alpha_e(par_pdg, mc, nf_out=3), get_alpha_e(par_pdg, mc, nf_out=4)) mb = config['RGE thresholds']['mb'] self.assertEqual(get_alpha_e(par_pdg, mb), get_alpha_e(par_pdg, mb, nf_out=5)) self.assertEqual(get_alpha_e(par_pdg, mb, nf_out=4), get_alpha_e(par_pdg, mb, nf_out=5)) mt = config['RGE thresholds']['mt'] self.assertEqual(get_alpha_e(par_pdg, mt), get_alpha_e(par_pdg, mt, nf_out=6)) self.assertEqual(get_alpha_e(par_pdg, mt, nf_out=5), get_alpha_e(par_pdg, mt, nf_out=6)) def test_alphas(self)-
Expand source code
def test_alphas(self): alpha_b = get_alpha(par_noqed, 4.2) # compare to 3-loop alpha_s at 4.2 GeV according to RunDec self.assertAlmostEqual(alpha_b['alpha_s']/0.225911,1.,places=4) def test_f_perp(self)-
Expand source code
def test_f_perp(self): par = flavio.default_parameters.get_central_all() self.assertEqual(par['f_perp_K*0'], get_f_perp(par, 'K*0', 2)) self.assertEqual(par['f_perp_K*0'], get_f_perp(par, 'K*0', 2, nf_out=3)) self.assertEqual(par['f_perp_K*0'], get_f_perp(par, 'K*0', 2, nf_out=4)) # rough comparison of 1503.05534 and 0804.0473 self.assertAlmostEqual(get_f_perp(par, 'K*0', 1) / par['f_perp_K*0'], 0.159 / (0.712 * 0.204), delta=0.025) self.assertAlmostEqual(get_f_perp(par, 'K*0', 1) / par['f_perp_K*0'], 0.159 / (0.712 * 0.204), delta=0.025) # just check this doesn't raise get_f_perp(par, 'phi', 100, nf_out=6) def test_ksmass(self)-
Expand source code
def test_ksmass(self): # compare to 1107.3100 # KS -> MSbar conversion self.assertAlmostEqual( flavio.physics.running.masses.mKS2mMS(4.55, 4, get_alpha(par_ks, 4.55)['alpha_s'], Mu=1, nl=2), 4.20, delta=0.01) self.assertAlmostEqual( flavio.physics.running.masses.mKS2mMS(1.15, 3, get_alpha(par_ks, 1.15)['alpha_s'], Mu=1, nl=2), 1.329, delta=0.01) # MSbar -> KS conversion self.assertAlmostEqual(get_mb_KS(par_ks, 1.), 4.553, delta=0.01) self.assertAlmostEqual(get_mc_KS(par_ks, 1.), 1.091, delta=0.13) # this is satisfied poorly! def test_mb1S(self)-
Expand source code
def test_mb1S(self): par = flavio.default_parameters.get_central_all() par['m_b'] = 4.18 par['alpha_s'] = 0.1185 alpha_s = get_alpha(par, par['m_b'], nf_out=5)['alpha_s'] mb1S = get_mb_1S(par, nl=3) self.assertAlmostEqual(mb1S, 4.67, delta=0.005) def test_polemasses(self)-
Expand source code
def test_polemasses(self): np.testing.assert_almost_equal(get_mb_pole(par, nl=2)/4.78248, 1,decimal=2) np.testing.assert_almost_equal(get_mc_pole(par, nl=2)/1.68375, 1,decimal=2) np.testing.assert_almost_equal(get_mb_pole(par, nl=3)/4.92987, 1,decimal=2) def test_runningmasses(self)-
Expand source code
def test_runningmasses(self): # compare to RunDec np.testing.assert_almost_equal(get_mb(par, 120.)/2.79211, 1,decimal=2) np.testing.assert_almost_equal(get_mt(par, 120.)/167.225, 1,decimal=2)