Skip to content

Commit

Permalink
adapt SM and CMSSM tests that do now pass
Browse files Browse the repository at this point in the history
  • Loading branch information
Alexander Voigt authored and Alexander Voigt committed Mar 20, 2018
1 parent 3ccdb54 commit a40b557
Show file tree
Hide file tree
Showing 4 changed files with 20 additions and 20 deletions.
8 changes: 4 additions & 4 deletions test/test_CMSSM_weinberg_angle.cpp
View file
Expand Up @@ -450,9 +450,9 @@ void setup_data(const CMSSM_input_parameters& input,
const double ss_pizztMZ = ss.piZZT(mz_pole, scale, false);
const double ss_piwwt0 = ss.piWWT(0., scale, false);
const double ss_piwwtMW = ss.piWWT(mw_pole, scale, false);
const double fs_pizztMZ = Re(fs.self_energy_VZ_1loop(mz_pole));
const double fs_pizztMZ = Re(fs.self_energy_VZ_1loop(Sqr(mz_pole)));
const double fs_piwwt0 = Re(fs.self_energy_VWm_1loop(0));
const double fs_piwwtMW = Re(fs.self_energy_VWm_1loop(mw_pole));
const double fs_piwwtMW = Re(fs.self_energy_VWm_1loop(Sqr(mw_pole)));

BOOST_CHECK_CLOSE_FRACTION(ss_pizztMZ, fs_pizztMZ, 5.0e-08);
BOOST_CHECK_CLOSE_FRACTION(ss_piwwtMW, fs_piwwtMW, 3.0e-05);
Expand Down Expand Up @@ -595,8 +595,8 @@ BOOST_AUTO_TEST_CASE( test_self_energy_top_correction )
const double ss_pizztMZ = ss.piZZT(mz_pole, scale, false);
const double ss_piwwtMW = ss.piWWT(mw_pole, scale, false);
const double ss_piwwt0 = ss.piWWT(0. , scale, false);
const double fs_pizztMZ = Re(fs.self_energy_VZ_1loop(mz_pole));
const double fs_piwwtMW = Re(fs.self_energy_VWm_1loop(mw_pole));
const double fs_pizztMZ = Re(fs.self_energy_VZ_1loop(Sqr(mz_pole)));
const double fs_piwwtMW = Re(fs.self_energy_VWm_1loop(Sqr(mw_pole)));
const double fs_piwwt0 = Re(fs.self_energy_VWm_1loop(0));

BOOST_CHECK_CLOSE_FRACTION(ss_pizztMZ, fs_pizztMZ, 5.0e-08);
Expand Down
8 changes: 4 additions & 4 deletions test/test_CMSSM_weinberg_angle_meta.cpp
View file
Expand Up @@ -85,8 +85,8 @@ void setup_data(const CMSSM_mass_eigenstates& model,
data.msve_drbar = msve;
data.msvm_drbar = msvm;

const double pizztMZ = Re(model.self_energy_VZ_1loop(mz_pole));
const double piwwtMW = Re(model.self_energy_VWm_1loop(mw_pole));
const double pizztMZ = Re(model.self_energy_VZ_1loop(Sqr(mz_pole)));
const double piwwtMW = Re(model.self_energy_VWm_1loop(Sqr(mw_pole)));
const double piwwt0 = Re(model.self_energy_VWm_1loop(0.));

Weinberg_angle::Self_energy_data se_data;
Expand Down Expand Up @@ -192,8 +192,8 @@ BOOST_AUTO_TEST_CASE( test_delta_r )
sm_parameters.alpha_s = 0.1176;
CMSSM_weinberg_angle wein(&model, sm_parameters);
// initialize self-energies
wein.pizzt_MZ = wein.calculate_self_energy_VZ(Electroweak_constants::MZ);
wein.piwwt_MW = wein.calculate_self_energy_VWm(Electroweak_constants::MW);
wein.pizzt_MZ = wein.calculate_self_energy_VZ(Sqr(Electroweak_constants::MZ));
wein.piwwt_MW = wein.calculate_self_energy_VWm(Sqr(Electroweak_constants::MW));
wein.piwwt_0 = wein.calculate_self_energy_VWm(0.);
const double delta_r_2 = wein.calculate_delta_r_hat(outrho, outsin);

Expand Down
16 changes: 8 additions & 8 deletions test/test_SM_weinberg_angle.cpp
View file
Expand Up @@ -197,9 +197,9 @@ BOOST_AUTO_TEST_CASE( test_delta_vb )
const double fs_g3 = fs.get_g3();
const double fs_hmu = fs.get_Ye(1,1);

const double pizztMZ = Re(fs.self_energy_VZ_1loop(mz_pole));
const double pizztMZ = Re(fs.self_energy_VZ_1loop(Sqr(mz_pole)));
const double piwwt0 = Re(fs.self_energy_VWp_1loop(0.));
const double piwwtMW = Re(fs.self_energy_VWp_1loop(mw_pole));
const double piwwtMW = Re(fs.self_energy_VWp_1loop(Sqr(mw_pole)));

Weinberg_angle::Self_energy_data se_data;
se_data.scale = scale;
Expand Down Expand Up @@ -270,9 +270,9 @@ BOOST_AUTO_TEST_CASE( test_delta_r )
const double fs_g3 = fs.get_g3();
const double fs_hmu = fs.get_Ye(1,1);

const double pizztMZ = Re(fs.self_energy_VZ_1loop(mz_pole));
const double pizztMZ = Re(fs.self_energy_VZ_1loop(Sqr(mz_pole)));
const double piwwt0 = Re(fs.self_energy_VWp_1loop(0.));
const double piwwtMW = Re(fs.self_energy_VWp_1loop(mw_pole));
const double piwwtMW = Re(fs.self_energy_VWp_1loop(Sqr(mw_pole)));

Weinberg_angle::Self_energy_data se_data;
se_data.scale = scale;
Expand Down Expand Up @@ -354,9 +354,9 @@ BOOST_AUTO_TEST_CASE( test_delta_rho )
const double fs_g3 = fs.get_g3();
const double fs_hmu = fs.get_Ye(1,1);

const double pizztMZ = Re(fs.self_energy_VZ_1loop(mz_pole));
const double pizztMZ = Re(fs.self_energy_VZ_1loop(Sqr(mz_pole)));
const double piwwt0 = Re(fs.self_energy_VWp_1loop(0.));
const double piwwtMW = Re(fs.self_energy_VWp_1loop(mw_pole));
const double piwwtMW = Re(fs.self_energy_VWp_1loop(Sqr(mw_pole)));

Weinberg_angle::Self_energy_data se_data;
se_data.scale = scale;
Expand Down Expand Up @@ -438,9 +438,9 @@ BOOST_AUTO_TEST_CASE( test_rho_sinTheta )
const double fs_g3 = fs.get_g3();
const double fs_hmu = fs.get_Ye(1,1);

const double pizztMZ = Re(fs.self_energy_VZ_1loop(mz_pole));
const double pizztMZ = Re(fs.self_energy_VZ_1loop(Sqr(mz_pole)));
const double piwwt0 = Re(fs.self_energy_VWp_1loop(0.));
const double piwwtMW = Re(fs.self_energy_VWp_1loop(mw_pole));
const double piwwtMW = Re(fs.self_energy_VWp_1loop(Sqr(mw_pole)));

Weinberg_angle::Self_energy_data se_data;
se_data.scale = scale;
Expand Down
8 changes: 4 additions & 4 deletions test/test_SM_weinberg_angle_meta.cpp
View file
Expand Up @@ -40,8 +40,8 @@ void setup_data(const SM_mass_eigenstates& model,
data.alpha_em_drbar = Sqr(e_drbar) / (4.0 * Pi);
data.mh_drbar = model.get_Mhh();

const double pizztMZ = Re(model.self_energy_VZ_1loop(mz_pole));
const double piwwtMW = Re(model.self_energy_VWp_1loop(mw_pole));
const double pizztMZ = Re(model.self_energy_VZ_1loop(Sqr(mz_pole)));
const double piwwtMW = Re(model.self_energy_VWp_1loop(Sqr(mw_pole)));
const double piwwt0 = Re(model.self_energy_VWp_1loop(0.));

Weinberg_angle::Self_energy_data se_data;
Expand Down Expand Up @@ -145,8 +145,8 @@ BOOST_AUTO_TEST_CASE( test_delta_r )
sm_parameters.alpha_s = 0.1176;
SM_weinberg_angle wein(&model, sm_parameters);
// initialize self-energies
wein.pizzt_MZ = wein.calculate_self_energy_VZ(Electroweak_constants::MZ);
wein.piwwt_MW = wein.calculate_self_energy_VWp(Electroweak_constants::MW);
wein.pizzt_MZ = wein.calculate_self_energy_VZ(Sqr(Electroweak_constants::MZ));
wein.piwwt_MW = wein.calculate_self_energy_VWp(Sqr(Electroweak_constants::MW));
wein.piwwt_0 = wein.calculate_self_energy_VWp(0.);
const double delta_r_2 = wein.calculate_delta_r_hat(outrho, outsin);

Expand Down

0 comments on commit a40b557

Please sign in to comment.