Skip to content

Commit

Permalink
adapting 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 21, 2018
1 parent 3184bf5 commit 0ad0c39
Show file tree
Hide file tree
Showing 6 changed files with 65 additions and 65 deletions.
4 changes: 2 additions & 2 deletions test/test_CMSSMCKM_low_scale_constraint.cpp
View file
Expand Up @@ -256,7 +256,7 @@ BOOST_AUTO_TEST_CASE( test_low_energy_constraint_with_flavour_mixing )
const double fs_mt = m.calculate_MFu_DRbar(qedqcd.displayPoleMt(), 2);
const double fs_mb = m.calculate_MFd_DRbar(qedqcd.displayMass(mBottom), 2);
const double fs_me = m.calculate_MFe_DRbar(qedqcd.displayMass(mTau), 2);
const double fs_MZ = m.calculate_MVZ_DRbar(Electroweak_constants::MZ);
const double fs_MZ = AbsSqrt(m.calculate_M2VZ_DRbar(Electroweak_constants::MZ));
const double fs_old_vd = m.get_vd();
const double fs_old_vu = m.get_vu();
// const double fs_old_vev = Sqrt(Sqr(fs_old_vu) + Sqr(fs_old_vd));
Expand All @@ -265,7 +265,7 @@ BOOST_AUTO_TEST_CASE( test_low_energy_constraint_with_flavour_mixing )
const double fs_new_vev = Sqrt(Sqr(fs_new_vu) + Sqr(fs_new_vd));

BOOST_CHECK_CLOSE_FRACTION(fs_mt, ss_mt, 9.5e-5);
BOOST_CHECK_CLOSE_FRACTION(fs_mb, ss_mb, 1.0e-12);
BOOST_CHECK_CLOSE_FRACTION(fs_mb, ss_mb, 2.0e-12);
BOOST_CHECK_CLOSE_FRACTION(fs_me, ss_me, 7.0e-4);
BOOST_CHECK_CLOSE_FRACTION(fs_MZ, ss_MZ, 5.0e-7);
BOOST_CHECK_CLOSE_FRACTION(fs_new_vev, ss_new_vev, 5.0e-7);
Expand Down
34 changes: 17 additions & 17 deletions test/test_CMSSMMassWInput_spectrum.cpp
View file
Expand Up @@ -246,8 +246,8 @@ void CMSSMMassWInput_weinberg_angle_low_scale_constraint::calculate_DRbar_gauge_
mW_run = qedqcd.displayPoleMW();

if (model->get_thresholds()) {
mZ_run = model->calculate_MVZ_DRbar(qedqcd.displayPoleMZ());
mW_run = model->calculate_MVWm_DRbar(qedqcd.displayPoleMW());
mZ_run = AbsSqrt(model->calculate_M2VZ_DRbar(qedqcd.displayPoleMZ()));
mW_run = AbsSqrt(model->calculate_M2VWm_DRbar(qedqcd.displayPoleMW()));
}

const double AlphaS = alpha_s_drbar;
Expand Down Expand Up @@ -324,9 +324,9 @@ void CMSSMMassWInput_weinberg_angle_low_scale_constraint::fill_data(
BOOST_REQUIRE(mz_pole > 0.);
BOOST_CHECK_EQUAL(scale, mz_pole);

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

Weinberg_angle::Self_energy_data se_data;
se_data.scale = scale;
Expand All @@ -348,9 +348,9 @@ void CMSSMMassWInput_weinberg_angle_low_scale_constraint::fill_data(
// check corrected self-energies
{
model->MFu(2) = mt_pole;
const double pizztMZ_check = Re(model->self_energy_VZ_1loop(mz_pole));
const double pizztMZ_check = Re(model->self_energy_VZ_1loop(Sqr(mz_pole)));
const double piwwt0_check = Re(model->self_energy_VWm_1loop(0));
const double piwwtMW_check = Re(model->self_energy_VWm_1loop(mw_pole));
const double piwwtMW_check = Re(model->self_energy_VWm_1loop(Sqr(mw_pole)));
model->MFu(2) = mt_drbar;

BOOST_CHECK_CLOSE_FRACTION(pizztMZ_check, pizztMZ_corrected, 1.0e-10);
Expand Down Expand Up @@ -790,9 +790,9 @@ BOOST_AUTO_TEST_CASE( test_CMSSMMassWInput_spectrum )
const DoubleVector
MCha_1l(ToDoubleVector(fs.get_physical().MCha)),
MChi_1l(ToDoubleVector(fs.get_physical().MChi)),
MHpm_1l(ToDoubleVector(fs.get_physical().MHpm)),
MAh_1l(ToDoubleVector(fs.get_physical().MAh)),
Mhh_1l(ToDoubleVector(fs.get_physical().Mhh));
MHpm_1l(ToDoubleVector(fs.get_physical().get_MHpm())),
MAh_1l(ToDoubleVector(fs.get_physical().get_MAh())),
Mhh_1l(ToDoubleVector(fs.get_physical().get_Mhh()));
const DoubleVector mch_1l(ss.displayPhys().mch),
mn_1l(ss.displayPhys().mneut.apply(fabs));
const double mHpm_1l = ss.displayPhys().mHpm;
Expand All @@ -816,7 +816,7 @@ BOOST_AUTO_TEST_CASE( test_CMSSMMassWInput_spectrum )
BOOST_CHECK_CLOSE_FRACTION(Mhh_1l(2), mH0_1l, 0.004);

// down-type squarks
const DoubleVector Sd_1l(ToDoubleVector(fs.get_physical().MSd));
const DoubleVector Sd_1l(ToDoubleVector(fs.get_physical().get_MSd()));
const DoubleVector md_1l(ss.displayPhys().md.flatten().sort());
BOOST_CHECK_CLOSE_FRACTION(Sd_1l(1), md_1l(1), 0.006);
BOOST_CHECK_CLOSE_FRACTION(Sd_1l(2), md_1l(2), 0.006);
Expand All @@ -826,7 +826,7 @@ BOOST_AUTO_TEST_CASE( test_CMSSMMassWInput_spectrum )
BOOST_CHECK_CLOSE_FRACTION(Sd_1l(6), md_1l(6), 0.006);

// up-type squarks
const DoubleVector Su_1l(ToDoubleVector(fs.get_physical().MSu));
const DoubleVector Su_1l(ToDoubleVector(fs.get_physical().get_MSu()));
const DoubleVector mu_1l(ss.displayPhys().mu.flatten().sort());
BOOST_CHECK_CLOSE_FRACTION(Su_1l(1), mu_1l(1), 0.006);
BOOST_CHECK_CLOSE_FRACTION(Su_1l(2), mu_1l(2), 0.006);
Expand All @@ -836,7 +836,7 @@ BOOST_AUTO_TEST_CASE( test_CMSSMMassWInput_spectrum )
BOOST_CHECK_CLOSE_FRACTION(Su_1l(6), mu_1l(6), 0.006);

// sleptons
const DoubleVector Se_1l(ToDoubleVector(fs.get_physical().MSe));
const DoubleVector Se_1l(ToDoubleVector(fs.get_physical().get_MSe()));
const DoubleVector me_1l(ss.displayPhys().me.flatten().sort());
BOOST_CHECK_CLOSE_FRACTION(Se_1l(1), me_1l(1), 0.0005);
BOOST_CHECK_CLOSE_FRACTION(Se_1l(2), me_1l(2), 0.0008);
Expand All @@ -847,7 +847,7 @@ BOOST_AUTO_TEST_CASE( test_CMSSMMassWInput_spectrum )

// sneutrinos
const DoubleVector msnu_1l(ss.displayPhys().msnu);
const DoubleVector Snu_1l(ToDoubleVector(fs.get_physical().MSv));
const DoubleVector Snu_1l(ToDoubleVector(fs.get_physical().get_MSv()));
BOOST_CHECK_CLOSE_FRACTION(Snu_1l(1), msnu_1l(1), 0.0056);
BOOST_CHECK_CLOSE_FRACTION(Snu_1l(2), msnu_1l(2), 0.0056);
BOOST_CHECK_CLOSE_FRACTION(Snu_1l(3), msnu_1l(3), 0.0090);
Expand Down Expand Up @@ -1035,11 +1035,11 @@ void CMSSMMassWInput_iterative_low_scale_constraint::apply()
model->set_vu(vu);

model->calculate_DRbar_masses();
model->calculate_Mhh_pole();
model->calculate_MVZ_pole();
model->calculate_M2hh_pole();
model->calculate_M2VZ_pole();

const double mH = model->get_physical().Mhh(0);
const double mZ = model->get_physical().MVZ;
const double mH = AbsSqrt(model->get_physical().M2hh(0));
const double mZ = AbsSqrt(model->get_physical().M2VZ);

#define LowEnergyConstant(p) Electroweak_constants::p
#define STANDARD_DEVIATION(p) Electroweak_constants::Error_##p
Expand Down
18 changes: 9 additions & 9 deletions test/test_NMSSM_one_loop_spectrum.cpp
View file
Expand Up @@ -125,8 +125,8 @@ BOOST_AUTO_TEST_CASE( test_NMSSM_pole_masses )
BOOST_CHECK_CLOSE(mch(2), MCha(2), 1.0e-12);

// photon, gluon mass
const double vp = m.get_physical().MVP;
const double vg = m.get_physical().MVG;
const double vp = m.get_physical().M2VP;
const double vg = m.get_physical().M2VG;
BOOST_CHECK_EQUAL(vp, 0.0);
BOOST_CHECK_EQUAL(vg, 0.0);

Expand All @@ -136,7 +136,7 @@ BOOST_AUTO_TEST_CASE( test_NMSSM_pole_masses )
BOOST_CHECK_CLOSE(MGlu, mGluino, 1.0e-04);

// down-type squarks
const DoubleVector Sd(ToDoubleVector(m.get_physical().MSd));
const DoubleVector Sd(ToDoubleVector(m.get_physical().get_MSd()));
const DoubleVector md(s.displayPhys().md.flatten().sort());
BOOST_CHECK_CLOSE(Sd(1), md(1), 1.0e-10);
BOOST_CHECK_CLOSE(Sd(2), md(2), 1.0e-10);
Expand All @@ -146,7 +146,7 @@ BOOST_AUTO_TEST_CASE( test_NMSSM_pole_masses )
BOOST_CHECK_CLOSE(Sd(6), md(6), 1.0e-10);

// up-type squarks
const DoubleVector Su(ToDoubleVector(m.get_physical().MSu));
const DoubleVector Su(ToDoubleVector(m.get_physical().get_MSu()));
const DoubleVector mu(s.displayPhys().mu.flatten().sort());
BOOST_CHECK_CLOSE(Su(1), mu(1), 1.0e-10);
BOOST_CHECK_CLOSE(Su(2), mu(2), 1.0e-10);
Expand All @@ -156,7 +156,7 @@ BOOST_AUTO_TEST_CASE( test_NMSSM_pole_masses )
BOOST_CHECK_CLOSE(Su(6), mu(6), 1.0e-10);

// down-type sleptons
const DoubleVector Se(ToDoubleVector(m.get_physical().MSe));
const DoubleVector Se(ToDoubleVector(m.get_physical().get_MSe()));
const DoubleVector me(s.displayPhys().me.flatten().sort());
BOOST_CHECK_CLOSE(Se(1), me(1), 1.0e-10);
BOOST_CHECK_CLOSE(Se(2), me(2), 1.0e-10);
Expand All @@ -166,7 +166,7 @@ BOOST_AUTO_TEST_CASE( test_NMSSM_pole_masses )
BOOST_CHECK_CLOSE(Se(6), me(6), 1.0e-10);

// up-type sleptons
const DoubleVector Sv(ToDoubleVector(m.get_physical().MSv));
const DoubleVector Sv(ToDoubleVector(m.get_physical().get_MSv()));
const DoubleVector msnu(s.displayPhys().msnu.sort());
BOOST_CHECK_CLOSE(Sv(1), msnu(1), 1.0e-10);
BOOST_CHECK_CLOSE(Sv(2), msnu(2), 1.0e-10);
Expand Down Expand Up @@ -207,20 +207,20 @@ BOOST_AUTO_TEST_CASE( test_NMSSM_pole_masses )
}

// neutral CP even Higgs
const DoubleVector hh(ToDoubleVector(m.get_physical().Mhh));
const DoubleVector hh(ToDoubleVector(m.get_physical().get_Mhh()));
const DoubleVector mh0(s.displayPhys().mh0);
BOOST_CHECK_CLOSE(hh(1), mh0(1), 0.002);
BOOST_CHECK_CLOSE(hh(2), mh0(2), 0.0001);
BOOST_CHECK_CLOSE(hh(3), mh0(3), 0.1);

// neutral CP odd Higgs
const DoubleVector Ah(ToDoubleVector(m.get_physical().MAh));
const DoubleVector Ah(ToDoubleVector(m.get_physical().get_MAh()));
const DoubleVector mA0(s.displayPhys().mA0);
BOOST_CHECK_CLOSE(Ah(2), mA0(1), 0.0002);
BOOST_CHECK_CLOSE(Ah(3), mA0(2), 0.0006);

// charged Higgs
const DoubleVector Hpm(ToDoubleVector(m.get_physical().MHpm));
const DoubleVector Hpm(ToDoubleVector(m.get_physical().get_MHpm()));
const double mHpm = s.displayPhys().mHpm;
BOOST_CHECK_CLOSE(Hpm(2), mHpm, 0.003);
}
8 changes: 4 additions & 4 deletions test/test_NMSSM_self_energies.cpp
View file
Expand Up @@ -94,13 +94,13 @@ BOOST_AUTO_TEST_CASE( test_NMSSM_self_energy_neutral_higgs )
s.setPhys(physical); // now modified

m.set_precision(std::pow(10,-0.1)); // performs 1 iteration
m.calculate_Mhh_pole();
m.calculate_MAh_pole();
m.calculate_M2hh_pole();
m.calculate_M2Ah_pole();

hh_ss = s.displayPhys().mh0;
hh_fs = m.get_physical().Mhh;
hh_fs = m.get_physical().get_Mhh();
Ah_ss = s.displayPhys().mA0;
Ah_fs = m.get_physical().MAh;
Ah_fs = m.get_physical().get_MAh();

BOOST_CHECK_CLOSE_FRACTION(hh_ss(1), hh_fs(0), 2.0e-04);
BOOST_CHECK_CLOSE_FRACTION(hh_ss(2), hh_fs(1), 4.0e-04);
Expand Down
20 changes: 10 additions & 10 deletions test/test_NMSSM_spectrum.cpp
View file
Expand Up @@ -384,9 +384,9 @@ BOOST_AUTO_TEST_CASE( test_NMSSM_spectrum )
const DoubleVector
MCha_1l(ToDoubleVector(fs.get_physical().MCha)),
MChi_1l(ToDoubleVector(fs.get_physical().MChi)),
MHpm_1l(ToDoubleVector(fs.get_physical().MHpm)),
MAh_1l(ToDoubleVector(fs.get_physical().MAh)),
Mhh_1l(ToDoubleVector(fs.get_physical().Mhh));
MHpm_1l(ToDoubleVector(fs.get_physical().get_MHpm())),
MAh_1l(ToDoubleVector(fs.get_physical().get_MAh())),
Mhh_1l(ToDoubleVector(fs.get_physical().get_Mhh()));
const DoubleVector mch_1l(ss.displayPhys().mch),
mn_1l(ss.displayPhys().mneut.apply(fabs));
const double mHpm_1l = ss.displayPhys().mHpm;
Expand Down Expand Up @@ -415,7 +415,7 @@ BOOST_AUTO_TEST_CASE( test_NMSSM_spectrum )
BOOST_TEST_MESSAGE("FlexibleSUSY:\n mh_1l = " << Mhh_1l << " mA_1l = " << MAh_1l);

// down-type squarks
const DoubleVector Sd_1l(ToDoubleVector(fs.get_physical().MSd));
const DoubleVector Sd_1l(ToDoubleVector(fs.get_physical().get_MSd()));
const DoubleVector md_1l(ss.displayPhys().md.flatten().sort());
BOOST_CHECK_CLOSE_FRACTION(Sd_1l(1), md_1l(1), 0.006);
BOOST_CHECK_CLOSE_FRACTION(Sd_1l(2), md_1l(2), 0.006);
Expand All @@ -425,7 +425,7 @@ BOOST_AUTO_TEST_CASE( test_NMSSM_spectrum )
BOOST_CHECK_CLOSE_FRACTION(Sd_1l(6), md_1l(6), 0.006);

// up-type squarks
const DoubleVector Su_1l(ToDoubleVector(fs.get_physical().MSu));
const DoubleVector Su_1l(ToDoubleVector(fs.get_physical().get_MSu()));
const DoubleVector mu_1l(ss.displayPhys().mu.flatten().sort());
BOOST_CHECK_CLOSE_FRACTION(Su_1l(1), mu_1l(1), 0.006);
BOOST_CHECK_CLOSE_FRACTION(Su_1l(2), mu_1l(2), 0.006);
Expand All @@ -435,7 +435,7 @@ BOOST_AUTO_TEST_CASE( test_NMSSM_spectrum )
BOOST_CHECK_CLOSE_FRACTION(Su_1l(6), mu_1l(6), 0.006);

// sleptons
const DoubleVector Se_1l(ToDoubleVector(fs.get_physical().MSe));
const DoubleVector Se_1l(ToDoubleVector(fs.get_physical().get_MSe()));
const DoubleVector me_1l(ss.displayPhys().me.flatten().sort());
BOOST_CHECK_CLOSE_FRACTION(Se_1l(1), me_1l(1), 0.0005);
BOOST_CHECK_CLOSE_FRACTION(Se_1l(2), me_1l(2), 0.0008);
Expand All @@ -446,7 +446,7 @@ BOOST_AUTO_TEST_CASE( test_NMSSM_spectrum )

// sneutrinos
const DoubleVector msnu_1l(ss.displayPhys().msnu.sort());
const DoubleVector Snu_1l(ToDoubleVector(fs.get_physical().MSv));
const DoubleVector Snu_1l(ToDoubleVector(fs.get_physical().get_MSv()));
BOOST_CHECK_CLOSE_FRACTION(Snu_1l(1), msnu_1l(1), 0.0056);
BOOST_CHECK_CLOSE_FRACTION(Snu_1l(2), msnu_1l(2), 0.0056);
BOOST_CHECK_CLOSE_FRACTION(Snu_1l(3), msnu_1l(3), 0.0090);
Expand All @@ -465,9 +465,9 @@ BOOST_AUTO_TEST_CASE( test_NMSSM_spectrum )
const NMSSM<Two_scale> fs_2l(nmssm_tester.get_model());

const DoubleVector
MHpm_2l(ToDoubleVector(fs_2l.get_physical().MHpm)),
MAh_2l(ToDoubleVector(fs_2l.get_physical().MAh)),
Mhh_2l(ToDoubleVector(fs_2l.get_physical().Mhh));
MHpm_2l(ToDoubleVector(fs_2l.get_physical().get_MHpm())),
MAh_2l(ToDoubleVector(fs_2l.get_physical().get_MAh())),
Mhh_2l(ToDoubleVector(fs_2l.get_physical().get_Mhh()));
const double mHpm_2l = ss_2l.displayPhys().mHpm;
const DoubleVector mA_2l(ss_2l.displayPhys().mA0);
const DoubleVector mh_2l(ss_2l.displayPhys().mh0);
Expand Down
46 changes: 23 additions & 23 deletions test/test_SMSSM_one_loop_spectrum.cpp
View file
Expand Up @@ -136,8 +136,8 @@ BOOST_AUTO_TEST_CASE( test_SMSSM_pole_masses )
BOOST_CHECK_CLOSE(mch(2), MCha(2), 1.0e-12);

// photon, gluon mass
const double vp = m.get_physical().MVP;
const double vg = m.get_physical().MVG;
const double vp = m.get_physical().M2VP;
const double vg = m.get_physical().M2VG;
BOOST_CHECK_EQUAL(vp, 0.0);
BOOST_CHECK_EQUAL(vg, 0.0);

Expand All @@ -147,17 +147,17 @@ BOOST_AUTO_TEST_CASE( test_SMSSM_pole_masses )
BOOST_CHECK_CLOSE(MGlu, mGluino, 4.0e-5);

// down-type squarks
const DoubleVector Sd(ToDoubleVector(m.get_physical().MSd));
const DoubleVector Sd(ToDoubleVector(m.get_physical().get_MSd()));
const DoubleVector md(s.displayPhys().md.flatten().sort());
BOOST_CHECK_CLOSE(Sd(1), md(1), 1.0e-10);
BOOST_CHECK_CLOSE(Sd(2), md(2), 1.0e-10);
BOOST_CHECK_CLOSE(Sd(3), md(3), 1.0e-10);
BOOST_CHECK_CLOSE(Sd(3), md(3), 1.0e-09);
BOOST_CHECK_CLOSE(Sd(4), md(4), 1.0e-10);
BOOST_CHECK_CLOSE(Sd(5), md(5), 1.0e-10);
BOOST_CHECK_CLOSE(Sd(6), md(6), 1.0e-10);

// up-type squarks
const DoubleVector Su(ToDoubleVector(m.get_physical().MSu));
const DoubleVector Su(ToDoubleVector(m.get_physical().get_MSu()));
const DoubleVector mu(s.displayPhys().mu.flatten().sort());
BOOST_CHECK_CLOSE(Su(1), mu(1), 1.0e-10);
BOOST_CHECK_CLOSE(Su(2), mu(2), 1.0e-10);
Expand All @@ -167,7 +167,7 @@ BOOST_AUTO_TEST_CASE( test_SMSSM_pole_masses )
BOOST_CHECK_CLOSE(Su(6), mu(6), 1.0e-10);

// down-type sleptons
const DoubleVector Se(ToDoubleVector(m.get_physical().MSe));
const DoubleVector Se(ToDoubleVector(m.get_physical().get_MSe()));
const DoubleVector me(s.displayPhys().me.flatten().sort());
BOOST_CHECK_CLOSE(Se(1), me(1), 1.0e-10);
BOOST_CHECK_CLOSE(Se(2), me(2), 1.0e-10);
Expand All @@ -177,7 +177,7 @@ BOOST_AUTO_TEST_CASE( test_SMSSM_pole_masses )
BOOST_CHECK_CLOSE(Se(6), me(6), 1.0e-10);

// up-type sleptons
const DoubleVector Sv(ToDoubleVector(m.get_physical().MSv));
const DoubleVector Sv(ToDoubleVector(m.get_physical().get_MSv()));
const DoubleVector msnu(s.displayPhys().msnu.sort());
BOOST_CHECK_CLOSE(Sv(1), msnu(1), 1.0e-10);
BOOST_CHECK_CLOSE(Sv(2), msnu(2), 1.0e-10);
Expand Down Expand Up @@ -218,19 +218,19 @@ BOOST_AUTO_TEST_CASE( test_SMSSM_pole_masses )
}

// neutral CP even Higgs
const DoubleVector hh(ToDoubleVector(m.get_physical().Mhh));
const DoubleVector hh(ToDoubleVector(m.get_physical().get_Mhh()));
const DoubleVector mh0(s.displayPhys().mh0);
BOOST_CHECK_CLOSE(hh(1), mh0(1), 0.0015);
BOOST_CHECK_CLOSE(hh(2), mh0(2), 2.0e-5);
BOOST_CHECK_CLOSE(hh(3), mh0(3), 0.008);

// neutral CP odd Higgs
const DoubleVector Ah(ToDoubleVector(m.get_physical().MAh));
const DoubleVector Ah(ToDoubleVector(m.get_physical().get_MAh()));
const DoubleVector mA0(s.displayPhys().mA0);
BOOST_CHECK_CLOSE(Ah(2), mA0(1), 3.0e-5);
BOOST_CHECK_CLOSE(Ah(3), mA0(2), 7.0e-4);
// charged Higgs
const DoubleVector Hpm(ToDoubleVector(m.get_physical().MHpm));
const DoubleVector Hpm(ToDoubleVector(m.get_physical().get_MHpm()));
const double mHpm = s.displayPhys().mHpm;
BOOST_CHECK_CLOSE(Hpm(2), mHpm, 0.0009);
}
Expand Down Expand Up @@ -286,19 +286,19 @@ BOOST_AUTO_TEST_CASE( test_self_energies )
BOOST_CHECK_CLOSE(nmTS, fsTS, 2.0e-12);

/// CP-even Higgs self-energies
const double momentum = 100.;
const double momentum = 100., momentum2 = Sqr(momentum);
const double nmHSE11 = s.pis1s1(momentum,s.displayMu());
const double fsHSE11 = m.self_energy_hh_1loop(momentum, 0, 0).real();
const double fsHSE11 = m.self_energy_hh_1loop(momentum2, 0, 0).real();
const double nmHSE12 = s.pis1s2(momentum,s.displayMu());
const double fsHSE12 = m.self_energy_hh_1loop(momentum, 0, 1).real();
const double fsHSE12 = m.self_energy_hh_1loop(momentum2, 0, 1).real();
const double nmHSE22 = s.pis2s2(momentum,s.displayMu());
const double fsHSE22 = m.self_energy_hh_1loop(momentum, 1, 1).real();
const double fsHSE22 = m.self_energy_hh_1loop(momentum2, 1, 1).real();
const double nmHSE13 = s.pis1s3(momentum,s.displayMu());
const double fsHSE13 = m.self_energy_hh_1loop(momentum, 0, 2).real();
const double fsHSE13 = m.self_energy_hh_1loop(momentum2, 0, 2).real();
const double nmHSE23 = s.pis2s3(momentum,s.displayMu());
const double fsHSE23 = m.self_energy_hh_1loop(momentum, 1, 2).real();
const double fsHSE23 = m.self_energy_hh_1loop(momentum2, 1, 2).real();
const double nmHSE33 = s.pis3s3(momentum,s.displayMu());
const double fsHSE33 = m.self_energy_hh_1loop(momentum, 2, 2).real();
const double fsHSE33 = m.self_energy_hh_1loop(momentum2, 2, 2).real();

BOOST_CHECK_CLOSE(nmHSE11, fsHSE11, 1.0e-12);
BOOST_CHECK_CLOSE(nmHSE12, fsHSE12, 1.0e-12);
Expand All @@ -309,17 +309,17 @@ BOOST_AUTO_TEST_CASE( test_self_energies )

/// CP-odd self-energies
const double nmPHSE11 = s.pip1p1(momentum,s.displayMu());
const double fsPHSE11 = m.self_energy_Ah_1loop(momentum, 0, 0).real();
const double fsPHSE11 = m.self_energy_Ah_1loop(momentum2, 0, 0).real();
const double nmPHSE12 = s.pip1p2(momentum,s.displayMu());
const double fsPHSE12 = m.self_energy_Ah_1loop(momentum, 0, 1).real();
const double fsPHSE12 = m.self_energy_Ah_1loop(momentum2, 0, 1).real();
const double nmPHSE22 = s.pip2p2(momentum,s.displayMu());
const double fsPHSE22 = m.self_energy_Ah_1loop(momentum, 1, 1).real();
const double fsPHSE22 = m.self_energy_Ah_1loop(momentum2, 1, 1).real();
const double nmPHSE13 = s.pip1p3(momentum,s.displayMu());
const double fsPHSE13 = m.self_energy_Ah_1loop(momentum, 0, 2).real();
const double fsPHSE13 = m.self_energy_Ah_1loop(momentum2, 0, 2).real();
const double nmPHSE23 = s.pip2p3(momentum,s.displayMu());
const double fsPHSE23 = m.self_energy_Ah_1loop(momentum, 1, 2).real();
const double fsPHSE23 = m.self_energy_Ah_1loop(momentum2, 1, 2).real();
const double nmPHSE33 = s.pip3p3(momentum,s.displayMu());
const double fsPHSE33 = m.self_energy_Ah_1loop(momentum, 2, 2).real();
const double fsPHSE33 = m.self_energy_Ah_1loop(momentum2, 2, 2).real();

BOOST_CHECK_CLOSE(nmPHSE11, fsPHSE11, 1.0e-12);
BOOST_CHECK_CLOSE(nmPHSE12, fsPHSE12, 1.0e-12);
Expand Down

0 comments on commit 0ad0c39

Please sign in to comment.