Skip to content

Commit

Permalink
put model setup into separate function
Browse files Browse the repository at this point in the history
  • Loading branch information
Alexander Voigt committed Dec 2, 2015
1 parent 6458f4c commit 9499486
Showing 1 changed file with 17 additions and 10 deletions.
27 changes: 17 additions & 10 deletions test/test_CMSSM_model.cpp
Expand Up @@ -1487,15 +1487,9 @@ void compare_self_energy_CP_odd_higgs(CMSSM<Two_scale> model,
TEST_CLOSE(Ah_ss(1), Ah_fs(1), 4.0e-3);
}

void compare_models(int loopLevel)
void setup_models(CMSSM<Two_scale>& m, MssmSoftsusy& softSusy,
CMSSM_input_parameters input, int loopLevel)
{
CMSSM_input_parameters input;
input.m0 = 125.;
input.m12 = 500.;
input.TanBeta = 10.;
input.SignMu = 1;
input.Azero = 0.;

const double ALPHASMZ = 0.1176;
const double ALPHAMZ = 1.0 / 127.918;
const double sinthWsq = 0.23122;
Expand Down Expand Up @@ -1533,7 +1527,6 @@ void compare_models(int loopLevel)
Ye(2,2) = 1.77699 * root2 / (vev * cosBeta);
mm0 = sqr(m0) * Eigen::Matrix<double,3,3>::Identity();

CMSSM<Two_scale> m(input);
m.set_scale(91);
m.set_loops(loopLevel);
m.set_g1(g1);
Expand All @@ -1560,7 +1553,6 @@ void compare_models(int loopLevel)
m.set_vu(vu);
m.set_vd(vd);

MssmSoftsusy softSusy;
softSusy.setMu(91);
softSusy.setLoops(loopLevel);
softSusy.setGaugeCoupling(1, g1);
Expand All @@ -1586,6 +1578,21 @@ void compare_models(int loopLevel)
softSusy.setM3Squared(BMu);
softSusy.setHvev(vev);
softSusy.setTanb(tanBeta);
}

void compare_models(int loopLevel)
{
CMSSM_input_parameters input;
input.m0 = 125.;
input.m12 = 500.;
input.TanBeta = 10.;
input.SignMu = 1;
input.Azero = 0.;

CMSSM<Two_scale> m(input);
MssmSoftsusy softSusy;

setup_models(m, softSusy, input, loopLevel);

std::cout << "comparing parameters ... ";
test_parameter_equality(softSusy, m);
Expand Down

0 comments on commit 9499486

Please sign in to comment.