change rhohat -> rhohat_ratio = rhohat / rhohat_tree

to use correct formulas; no difference in MSSM because of rhohat_tree = 1
FlexibleSUSY · Nov 17, 2015 · 93491da · 93491da
1 parent 2e0f17d
commit 93491da
Show file tree

Hide file tree

Showing 2 changed files with 15 additions and 15 deletions.
diff --git a/meta/WeinbergAngleMuonDecay.m b/meta/WeinbergAngleMuonDecay.m
@@ -20,7 +20,7 @@
            Return[massExpr[[1]] /. SARAH`Weinberg[] -> SARAH`Weinberg];
           ];
 
-extPars={SINTHETAW, RHOHAT, GFERMI, MW, MZ, MT, RHO2, DELTARHAT1LOOP, PIZZTMZ};
+extPars={SINTHETAW, RHOHATRATIO, GFERMI, MW, MZ, MT, RHO2, DELTARHAT1LOOP, PIZZTMZ};
 Do[Format[extPars[[i]],CForm]=Format[ToString[extPars[[i]]],OutputForm],{i,Length[extPars]}];
 
 (*formula according to (C.6) from hep-ph/9606211*)
@@ -42,7 +42,7 @@
                            alphaDRbar = gY^2 SARAH`leftCoupling^2 / (4 Pi (gY^2 + SARAH`leftCoupling^2));
                            xt = 3 / (8 Pi^2 Sqrt[2]) GFERMI MT^2;
                            hmix = (SARAH`HiggsMixingMatrix[0,1] / Sin[ArcTan[FlexibleSUSY`vu/FlexibleSUSY`vd]])^2;
-                           expr = alphaDRbar SARAH`strongCoupling^2/(16 Pi^3 SINTHETAW^2 (1 - SINTHETAW^2))(2.145 MT^2/MZ^2 + 0.575 Log[MT/MZ] - 0.224 - 0.144 MZ^2/MT^2) - xt^2 hmix RHO2[FlexibleSUSY`M[SARAH`HiggsBoson[0]]/MT] (1 - DELTARHAT1LOOP) RHOHAT / 3;
+                           expr = alphaDRbar SARAH`strongCoupling^2/(16 Pi^3 SINTHETAW^2 (1 - SINTHETAW^2))(2.145 MT^2/MZ^2 + 0.575 Log[MT/MZ] - 0.224 - 0.144 MZ^2/MT^2) - xt^2 hmix RHO2[FlexibleSUSY`M[SARAH`HiggsBoson[0]]/MT] (1 - DELTARHAT1LOOP) RHOHATRATIO / 3;
                            result = Parameters`CreateLocalConstRefs[expr] <> "\n";
                            result = result <> "deltaRHat2LoopSM = ";
                            result = result <> CConversion`RValueToCFormString[expr] <> ";";

diff --git a/templates/weinberg_angle.cpp.in b/templates/weinberg_angle.cpp.in
@@ -31,7 +31,7 @@ namespace flexiblesusy {
 
 #define MODELPARAMETER(p) model->get_##p()
 #define SINTHETAW sinThetaW
-#define RHOHAT rhohat
+#define RHOHATRATIO rhohat_ratio
 #define GFERMI gfermi
 #define MW mw
 #define MZ mz
@@ -149,13 +149,13 @@ double @ModelName@_weinberg_angle::calculate(double sinThetaW_start)
    while (not_converged && iteration < number_of_iterations) {
       const double deltaRhoHat = calculate_delta_rho_hat(sinThetaW_old);
 
-      double rhohat;
+      double rhohat_ratio;  // = rhohat / rhohat_tree
       if (Abs(deltaRhoHat) < 1.0)
-         rhohat = 1.0 / (1.0 - deltaRhoHat);
+         rhohat_ratio = 1.0 / (1.0 - deltaRhoHat);
       else
-         rhohat = 1.0;
+         rhohat_ratio = 1.0;
 
-      const double deltaRHat = calculate_delta_r_hat(rhohat, sinThetaW_old);
+      const double deltaRHat = calculate_delta_r_hat(rhohat_ratio, sinThetaW_old);
 
       double sin2thetasqO4 = Pi * alphaDRbar /
          (ROOT2 * Sqr(mz) * gfermi * (1.0 - deltaRHat) * rhohat_tree);
@@ -176,7 +176,7 @@ double @ModelName@_weinberg_angle::calculate(double sinThetaW_start)
       VERBOSE_MSG("Iteration step " << iteration
                   << ": prec=" << precision
                   << " dRhoHat=" << deltaRhoHat
-                  << " rhohat=" << rhohat
+                  << " rhohat_ratio=" << rhohat_ratio
                   << " dRHat=" << deltaRHat
                   << " sinThetaW_new=" << sinThetaW_new);
 
@@ -229,7 +229,7 @@ double @ModelName@_weinberg_angle::calculate_delta_rho_hat(double sinThetaW)
  *
  * @return \f$\Delta\hat{r}\f$ as defined in (C.3) and (C.5) from hep-ph/9606211
  */
-double @ModelName@_weinberg_angle::calculate_delta_r_hat(double rhohat, double sinThetaW)
+double @ModelName@_weinberg_angle::calculate_delta_r_hat(double rhohat_ratio, double sinThetaW)
 {
    const double gfermi = sm_parameters.fermi_constant;
    const double mw = sm_parameters.mw_pole;
@@ -238,11 +238,11 @@ double @ModelName@_weinberg_angle::calculate_delta_r_hat(double rhohat, double s
 
    double dvb = 0.;
    if (number_of_loops > 0)
-      dvb = calculate_delta_vb(rhohat, sinThetaW);
+      dvb = calculate_delta_vb(rhohat_ratio, sinThetaW);
 
    double deltaRHat1Loop = 0.;
    if (number_of_loops > 0)
-      deltaRHat1Loop = rhohat * piwwt_0 / Sqr(mw) - pizzt_MZ / Sqr(mz) + dvb;
+      deltaRHat1Loop = rhohat_ratio * piwwt_0 / Sqr(mw) - pizzt_MZ / Sqr(mz) + dvb;
 
    double deltaRHat2LoopSM = 0.;
    if (number_of_loops > 1) {
@@ -263,9 +263,9 @@ double @ModelName@_weinberg_angle::calculate_delta_r_hat(double rhohat, double s
  *
  * @return \f$\delta_{\text{VB}}\f$ as defined in (C.11) from hep-ph/9606211
  */
-double @ModelName@_weinberg_angle::calculate_delta_vb(double rhohat, double sinThetaW)
+double @ModelName@_weinberg_angle::calculate_delta_vb(double rhohat_ratio, double sinThetaW)
 {
-   const double deltaVbSM = calculate_delta_vb_sm(rhohat, sinThetaW);
+   const double deltaVbSM = calculate_delta_vb_sm(rhohat_ratio, sinThetaW);
 
    double deltaVbSusy = 0.;
    if (include_dvbsusy)
@@ -288,7 +288,7 @@ double @ModelName@_weinberg_angle::calculate_delta_vb(double rhohat, double sinT
  * from hep-ph/9606211
  */
 double @ModelName@_weinberg_angle::calculate_delta_vb_sm(
-   double rhohat, double sinThetaW)
+   double rhohat_ratio, double sinThetaW)
 {
    const double mz  = sm_parameters.mz_pole;
    const double mw  = sm_parameters.mw_pole;
@@ -303,7 +303,7 @@ double @ModelName@_weinberg_angle::calculate_delta_vb_sm(
    const double eDRbar     = gY * g2 / Sqrt(Sqr(gY) + Sqr(g2));
    const double alphaDRbar = Sqr(eDRbar) / (4.0 * Pi);
 
-   const double deltaVbSM = rhohat * alphaDRbar / (4.0 * Pi * sinThetaW2) *
+   const double deltaVbSM = rhohat_ratio * alphaDRbar / (4.0 * Pi * sinThetaW2) *
       (6.0 + log(cw2) / sw2 *
        (3.5 - 2.5 * sw2 - sinThetaW2 * (5.0 - 1.5 * cw2 / outcos2))
        - 4. * Log(Sqr(mz/q)));