add new meta file for generation of 2-loop SM contributions to deltaR…

… and deltaRho
FlexibleSUSY · Jul 1, 2015 · c2a53dc · c2a53dc
1 parent 1cb841f
commit c2a53dc
Show file tree

Hide file tree

Showing 5 changed files with 32 additions and 5 deletions.
diff --git a/meta/FlexibleSUSY.m b/meta/FlexibleSUSY.m
@@ -1,5 +1,5 @@
 
-BeginPackage["FlexibleSUSY`", {"SARAH`", "AnomalousDimension`", "BetaFunction`", "TextFormatting`", "CConversion`", "TreeMasses`", "EWSB`", "Traces`", "SelfEnergies`", "Vertices`", "Phases`", "LoopMasses`", "WriteOut`", "Constraint`", "ThresholdCorrections`", "ConvergenceTester`", "Utils`", "ThreeLoopSM`"}];
+BeginPackage["FlexibleSUSY`", {"SARAH`", "AnomalousDimension`", "BetaFunction`", "TextFormatting`", "CConversion`", "TreeMasses`", "EWSB`", "Traces`", "SelfEnergies`", "Vertices`", "Phases`", "LoopMasses`", "WriteOut`", "Constraint`", "ThresholdCorrections`", "ConvergenceTester`", "Utils`", "ThreeLoopSM`", "WeinbergAngleMuonDecay`"}];
 
 FS`Version = StringTrim[FSImportString[FileNameJoin[{Global`$flexiblesusyConfigDir,"version"}]]];
 FS`GitCommit = StringTrim[FSImportString[FileNameJoin[{Global`$flexiblesusyConfigDir,"git_commit"}]]];
@@ -610,7 +610,7 @@
 WriteWeinbergAngleClass[files_List] :=
    Module[{},
           WriteOut`ReplaceInFiles[files,
-                 { "@HalloOma@" -> IndentText["IchSageOmaHallo"],
+                 { "@deltaRho2LoopSM@" -> IndentText[WeinbergAngleMuonDecay`deltaRho2LoopSM[]],
                    Sequence @@ GeneralReplacementRules[]
                  } ];
           ];

diff --git a/meta/WeinbergAngleMuonDecay.m b/meta/WeinbergAngleMuonDecay.m
@@ -0,0 +1,23 @@
+(* ::Package:: *)
+
+BeginPackage["WeinbergAngleMuonDecay`", {"SARAH`", "CConversion`", "Parameters`", "TextFormatting`"}];
+
+deltaRho2LoopSM::usage="";
+
+Begin["`Private`"];
+
+deltaRho2LoopSM[]:=Module[{gY, alphaDRbar, xt, gfermi, mt, hmix, expr, sinThetaW, mw, mz, RHO2, result},
+                          gY = g1 Sqrt[0.6];
+                          alphaDRbar = gY^2 Susyno`LieGroups`g2^2 / (4 Pi (gY^2 + Susyno`LieGroups`g2^2));
+                          xt = 3.0 / (8 Pi^2 Sqrt[2]) gfermi mt^2;
+                          hmix = (FlexibleSUSY`ZH[0,1] / Sin[ArcTan[FlexibleSUSY`vu/FlexibleSUSY`vd]])^2;
+                          expr = alphaDRbar g3^2/(16 Pi^3 sinThetaW^2)(-2.145 mt^2/mw^2 + 1.262 Log[mt/mz] - 2.24 - 0.85 mz^2/mt^2) + xt^2 hmix RHO2[FlexibleSUSY`M[hh[0]]/mt] / 3;
+                          result = Parameters`CreateLocalConstRefs[expr];
+                          result = result <> "\n" <> "\n";
+                          result = result <> CConversion`RValueToCFormString[expr];
+                          Return[result];
+                          ];
+
+End[];
+
+EndPackage[];
diff --git a/meta/module.mk b/meta/module.mk
@@ -24,7 +24,8 @@ META_SRC     := \
 		$(DIR)/Utils.m \
 		$(DIR)/Vertices.m \
 		$(DIR)/WriteOut.m \
-		$(DIR)/WeinbergAngle.m
+		$(DIR)/WeinbergAngle.m \
+		$(DIR)/WeinbergAngleMuonDecay.m
 
 .PHONY:         all-$(MODNAME) clean-$(MODNAME) distclean-$(MODNAME)
 

diff --git a/templates/weinberg_angle.cpp.in b/templates/weinberg_angle.cpp.in
@@ -200,10 +200,11 @@ double @ModelName@_weinberg_angle::calculate_delta_rho(double rhohat, double sin
       deltaRho1Loop = pizztMZ_corrected / (rhohat * Sqr(mz)) -
          piwwtMW_corrected / Sqr(mw);
 
+   const double gfermi = sm_parameters.fermi_constant;
+
    const double mh = model->get_M@HiggsBoson_0@;
    const double g3 = model->get_@strongCoupling@();
-   const double xt = 3.0 * sm_parameters.fermi_constant * Sqr(mt) *
-                     ROOT2 * oneOver16PiSqr;
+   const double xt = 3.0 * gfermi * Sqr(mt) * ROOT2 * oneOver16PiSqr;
    const double sinb   = Sin(ArcTan(model->get_vu() / model->get_vd()));
    const double hmix_r = Sqr(model->get_ZH(0,1) / sinb);
 

diff --git a/templates/weinberg_angle.hpp.in b/templates/weinberg_angle.hpp.in
@@ -77,3 +77,5 @@ private:
 } // namespace flexiblesusy
 
 #endif
+
+@deltaRho2LoopSM@