From 0c2a9c0a8e9297822988dcf884964e3a844f1b4e Mon Sep 17 00:00:00 2001 From: Alexander Voigt Date: Tue, 26 Jul 2016 16:36:35 +0200 Subject: [PATCH] calculate full EFT spectrum at the 1L level --- ...d_model_high_scale_spectrum_generator.hpp.in | 17 ++++------------- ...rd_model_low_scale_spectrum_generator.hpp.in | 12 ++++-------- 2 files changed, 8 insertions(+), 21 deletions(-) diff --git a/templates/standard_model_high_scale_spectrum_generator.hpp.in b/templates/standard_model_high_scale_spectrum_generator.hpp.in index d427c710c..09fc9182b 100644 --- a/templates/standard_model_high_scale_spectrum_generator.hpp.in +++ b/templates/standard_model_high_scale_spectrum_generator.hpp.in @@ -196,21 +196,12 @@ void @ModelName@_standard_model_spectrum_generator::run(const softsusy::QedQc eft.run_to(Q_higgs); + // 1L computation of pole mass spectrum in the SM eft.calculate_DRbar_masses(); - eft.solve_ewsb(); - // eft.calculate_spectrum(); - // does not work properly, find out why! - - eft.calculate_MVWp_pole(); - eft.calculate_MVZ_pole(); - eft.calculate_MFu_pole(); - eft.calculate_MFd_pole(); - eft.calculate_MFe_pole(); - - // since Higgs mass matching is locked to 1-loop, so is this calculation eft.set_pole_mass_loop_order(1); - eft.solve_ewsb_one_loop(); - eft.calculate_Mhh_pole(); + eft.set_ewsb_loop_order(1); + eft.solve_ewsb(); + eft.calculate_spectrum(); this->model.get_physical().M@HiggsBoson_0@ = eft.get_physical().Mhh; this->model.get_physical().M@VectorZ_0@ = eft.get_physical().MVZ; diff --git a/templates/standard_model_low_scale_spectrum_generator.hpp.in b/templates/standard_model_low_scale_spectrum_generator.hpp.in index fa5c72cf4..658864a5b 100644 --- a/templates/standard_model_low_scale_spectrum_generator.hpp.in +++ b/templates/standard_model_low_scale_spectrum_generator.hpp.in @@ -186,16 +186,12 @@ void @ModelName@_standard_model_spectrum_generator::run(const softsusy::QedQc eft.run_to(Q_higgs); + // 1L computation of pole mass spectrum in the SM eft.calculate_DRbar_masses(); + eft.set_pole_mass_loop_order(1); + eft.set_ewsb_loop_order(1); eft.solve_ewsb(); - // eft.calculate_spectrum(); - // does not work properly, find out why! - - eft.calculate_MVWp_pole(); - eft.calculate_MVZ_pole(); - eft.calculate_MFu_pole(); - eft.calculate_MFd_pole(); - eft.calculate_MFe_pole(); + eft.calculate_spectrum(); // since Higgs mass matching is locked to 1-loop, so is this calculation eft.set_pole_mass_loop_order(1);