allow tree-level calculation in the tower

FlexibleSUSY · Jul 27, 2016 · 29bdad1 · 29bdad1
1 parent 71ce310
commit 29bdad1
Show file tree

Hide file tree

Showing 3 changed files with 45 additions and 8 deletions.
diff --git a/templates/standard_model_high_scale_spectrum_generator.hpp.in b/templates/standard_model_high_scale_spectrum_generator.hpp.in
@@ -199,13 +199,28 @@ void @ModelName@_standard_model_spectrum_generator<T>::run(const softsusy::QedQc
 
          eft.run_to(Q_higgs);
 
-         // 1L computation of pole mass spectrum in the SM
+         // computation of pole mass spectrum in the SM
+         const unsigned eft_pole_loops = eft.get_pole_mass_loop_order();
+         const unsigned eft_ewsb_loops = eft.get_ewsb_loop_order();
+
+         if (eft_pole_loops > 1) {
+            WARNING("Pole mass loop order " << eft_pole_loops
+                    << " in the EFT is currently not supported!  I'm using 1-loop.");
+            eft.set_pole_mass_loop_order(1);
+         }
+         if (eft_ewsb_loops > 1) {
+            WARNING("EWSB loop order " << eft_ewsb_loops
+                    << " in the EFT is currently not supported!  I'm using 1-loop.");
+            eft.set_ewsb_loop_order(1);
+         }
+
          eft.calculate_DRbar_masses();
-         eft.set_pole_mass_loop_order(1);
-         eft.set_ewsb_loop_order(1);
          eft.solve_ewsb();
          eft.calculate_spectrum();
 
+         eft.set_pole_mass_loop_order(eft_pole_loops);
+         eft.set_ewsb_loop_order(eft_ewsb_loops);
+
          this->model.get_physical().M@HiggsBoson_0@ = eft.get_physical().Mhh;
          this->model.get_physical().M@VectorZ_0@ = eft.get_physical().MVZ;
          this->model.get_physical().M@VectorW_0@ = eft.get_physical().MVWp;

diff --git a/templates/standard_model_low_scale_spectrum_generator.hpp.in b/templates/standard_model_low_scale_spectrum_generator.hpp.in
@@ -190,13 +190,28 @@ void @ModelName@_standard_model_spectrum_generator<T>::run(const softsusy::QedQc
 
          eft.run_to(Q_higgs);
 
-         // 1L computation of pole mass spectrum in the SM
+         // computation of pole mass spectrum in the SM
+         const unsigned eft_pole_loops = eft.get_pole_mass_loop_order();
+         const unsigned eft_ewsb_loops = eft.get_ewsb_loop_order();
+
+         if (eft_pole_loops > 1) {
+            WARNING("Pole mass loop order " << eft_pole_loops
+                    << " in the EFT is currently not supported!  I'm using 1-loop.");
+            eft.set_pole_mass_loop_order(1);
+         }
+         if (eft_ewsb_loops > 1) {
+            WARNING("EWSB loop order " << eft_ewsb_loops
+                    << " in the EFT is currently not supported!  I'm using 1-loop.");
+            eft.set_ewsb_loop_order(1);
+         }
+
          eft.calculate_DRbar_masses();
-         eft.set_pole_mass_loop_order(1);
-         eft.set_ewsb_loop_order(1);
          eft.solve_ewsb();
          eft.calculate_spectrum();
 
+         eft.set_pole_mass_loop_order(eft_pole_loops);
+         eft.set_ewsb_loop_order(eft_ewsb_loops);
+
          this->model.get_physical().M@HiggsBoson_0@ = eft.get_physical().Mhh;
          this->model.get_physical().M@VectorZ_0@ = eft.get_physical().MVZ;
          this->model.get_physical().M@VectorW_0@ = eft.get_physical().MVWp;

diff --git a/templates/standard_model_two_scale_matching.cpp.in b/templates/standard_model_two_scale_matching.cpp.in
@@ -42,8 +42,9 @@ CLASSNAME::@ModelName@_standard_model_Matching(
    , eft(low_)
    , constraint(constraint_)
    , scale(0.)
-   , loop_order(loop_order_)
-{}
+{
+   set_loop_order(loop_order_);
+}
 
 CLASSNAME::~@ModelName@_standard_model_Matching()
 {}
@@ -134,6 +135,12 @@ unsigned CLASSNAME::get_loop_order() const
 
 void CLASSNAME::set_loop_order(unsigned loop_order_)
 {
+   if (loop_order_ > 1) {
+      WARNING("Matching loop order " << loop_order_ << " currently not"
+              " supported!  I'm using 1-loop matching.");
+      loop_order_ = 1;
+   }
+
    loop_order = loop_order_;
 }