adapt examples to new QedQcd class

FlexibleSUSY · Dec 21, 2016 · e692638 · e692638
1 parent f5fe626
commit e692638
Show file tree

Hide file tree

Showing 3 changed files with 9 additions and 9 deletions.
diff --git a/examples/customized-betas/MSSMcbs_two_scale_initial_guesser.cpp b/examples/customized-betas/MSSMcbs_two_scale_initial_guesser.cpp
@@ -92,11 +92,11 @@ void MSSMcbs_initial_guesser<Two_scale>::guess_susy_parameters()
    mtau_guess = leAtMt.displayMass(mTau);
 
    // guess gauge couplings at mt
-   const DoubleVector alpha_sm(leAtMt.getGaugeMu(mtpole, sinThetaW2));
+   const auto alpha_sm(leAtMt.getGaugeMu(mtpole, sinThetaW2));
 
-   model->set_g1(sqrt(4.0 * M_PI * alpha_sm(1)));
-   model->set_g2(sqrt(4.0 * M_PI * alpha_sm(2)));
-   model->set_g3(sqrt(4.0 * M_PI * alpha_sm(3)));
+   model->set_g1(sqrt(4.0 * M_PI * alpha_sm(0)));
+   model->set_g2(sqrt(4.0 * M_PI * alpha_sm(1)));
+   model->set_g3(sqrt(4.0 * M_PI * alpha_sm(2)));
    model->set_scale(mtpole);
 
    // apply user-defined initial guess at the low scale

diff --git a/examples/customized-betas/MSSMcbs_two_scale_low_scale_constraint.cpp b/examples/customized-betas/MSSMcbs_two_scale_low_scale_constraint.cpp
@@ -151,7 +151,7 @@ void MSSMcbs_low_scale_constraint<Two_scale>::update_scale()
 
 void MSSMcbs_low_scale_constraint<Two_scale>::calculate_DRbar_gauge_couplings()
 {
-   assert(qedqcd.displayMu() == get_scale() && "Error: low-energy data"
+   assert(qedqcd.get_scale() == get_scale() && "Error: low-energy data"
           " set is not defined at the same scale as the low-energy"
           " constraint.  You need to run the low-energy data set to this"
           " scale!");

diff --git a/examples/tower/MSSMD5O_MSSMRHN_two_scale_initial_guesser.cpp b/examples/tower/MSSMD5O_MSSMRHN_two_scale_initial_guesser.cpp
@@ -79,11 +79,11 @@ void MSSMD5O_MSSMRHN_initial_guesser<Two_scale>::guess_susy_parameters()
    const double mtpole = leAtMt.displayPoleMt();
 
    // guess gauge couplings at mt
-   const DoubleVector alpha_sm(leAtMt.getGaugeMu(mtpole, sinThetaW2));
+   const auto alpha_sm(leAtMt.getGaugeMu(mtpole, sinThetaW2));
 
-   model_1->set_g1(sqrt(4.0 * M_PI * alpha_sm(1)));
-   model_1->set_g2(sqrt(4.0 * M_PI * alpha_sm(2)));
-   model_1->set_g3(sqrt(4.0 * M_PI * alpha_sm(3)));
+   model_1->set_g1(sqrt(4.0 * M_PI * alpha_sm(0)));
+   model_1->set_g2(sqrt(4.0 * M_PI * alpha_sm(1)));
+   model_1->set_g3(sqrt(4.0 * M_PI * alpha_sm(2)));
    model_1->set_scale(mtpole);
    model_1->set_loops(2);