generalize 2L Higgs contributions

FlexibleSUSY · Apr 26, 2016 · 7e6f646 · 7e6f646
1 parent 547bb97
commit 7e6f646
Showing 1 changed file with 27 additions and 10 deletions.
diff --git a/meta/WeinbergAngle.m b/meta/WeinbergAngle.m
@@ -91,31 +91,48 @@
 extPars={SINTHETAW, RHOHATRATIO, GFERMI, MW, MZ, MT, RHO2, DELTARHAT1LOOP, PIZZTMZ};
 Do[Format[extPars[[i]],CForm]=Format[ToString[extPars[[i]]],OutputForm],{i,Length[extPars]}];
 
+(*returns coefficients of Higgs-top-top vertices*)
+HiggsTopVertices[higgsName_] :=
+    Module[{indexRange, indexList, higgsVertices},
+           indexRange = TreeMasses`GetParticleIndices[higgsName][[All, 2]];
+           If[indexRange === {}, indexRange = {1}];
+           indexList = Flatten[Table @@ {Table[ToExpression["i" <> ToString[k]], {k, Length[indexRange]}], Sequence @@ Table[{ToExpression["i" <> ToString[k]], 1, indexRange[[k]]}, {k, Length[indexRange]}]}, Length[indexRange] - 1];
+           higgsVertices = Vertices`StripGroupStructure[SARAH`Vertex[{bar[TreeMasses`GetUpQuark[{3}]], TreeMasses`GetUpQuark[{3}], higgsName[#]}] & /@ indexList, SARAH`ctNr /@ Range[4]] //. SARAH`sum[idx_, start_, stop_, expr_] :> Sum[expr, {idx, start, stop}];
+           higgsVertices = Cases[higgsVertices, {{__, higgsField_}, {coeffPL_, SARAH`PL}, {coeffPR_, SARAH`PR}} /; (coeffPL/I * Susyno`LieGroups`conj[coeffPL/I] === coeffPR/I * Susyno`LieGroups`conj[coeffPR/I])
+                                 && !TreeMasses`IsGoldstone[higgsField] :> {higgsField /. List -> Sequence, coeffPL/I}];
+           Return[higgsVertices];
+          ];
+
+(*generalize Higgs dependent part of (C.5) and (C.6) in hep-ph/9606211 analogous to (C.9) and (C.10)*)
+HiggsContributions2LoopSM[] :=
+    Module[{higgsVEVlist, higgsDep},
+           higgsVEVlist = Cases[Parameters`GetDependenceSPhenoRules[], RuleDelayed[SARAH`VEVSM, repr_] :> repr];
+           If[higgsVEVlist === {}, higgsVEVlist = {SARAH`VEVSM}];
+           higgsDep = Abs[#[[2]]]^2 RHO2[FlexibleSUSY`M[#[[1]]]/MT] &;
+           Return[Simplify[3 (GFERMI MT higgsVEVlist[[1]] / (8 Pi^2 Sqrt[2]))^2 (Plus @@ (higgsDep /@ HiggsTopVertices[SARAH`HiggsBoson]) - Plus @@ (higgsDep /@ HiggsTopVertices[SARAH`PseudoScalar]))]];
+          ];
+
 (*formula according to (C.6) from hep-ph/9606211*)
 deltaRhoHat2LoopSM[]:=
-    Module[{gY, alphaDRbar, xt, hmix, expr, result},
+    Module[{gY, alphaDRbar, expr, result},
            gY = SARAH`hyperchargeCoupling FlexibleSUSY`GUTNormalization[SARAH`hyperchargeCoupling];
            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)(-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[SARAH`HiggsBoson[0]]/MT] / 3) / (1 + PIZZTMZ / MZ^2);
+           expr = (alphaDRbar SARAH`strongCoupling^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) + HiggsContributions2LoopSM[]) / (1 + PIZZTMZ / MZ^2);
            result = Parameters`CreateLocalConstRefs[expr] <> "\n";
            result = result <> "deltaRhoHat2LoopSM = ";
-           result = result <> CConversion`RValueToCFormString[expr] <> ";";
+           result = result <> CConversion`RValueToCFormString[Parameters`DecreaseIndexLiterals[expr]] <> ";";
            Return[result];
           ];
 
 (*formula according to (C.5) from hep-ph/9606211*)
 deltaRHat2LoopSM[]:=
-    Module[{gY, alphaDRbar, xt, hmix, expr, result},
+    Module[{gY, alphaDRbar, expr, result},
            gY = SARAH`hyperchargeCoupling FlexibleSUSY`GUTNormalization[SARAH`hyperchargeCoupling];
            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) RHOHATRATIO / 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) - HiggsContributions2LoopSM[] (1 - DELTARHAT1LOOP) RHOHATRATIO;
            result = Parameters`CreateLocalConstRefs[expr] <> "\n";
            result = result <> "deltaRHat2LoopSM = ";
-           result = result <> CConversion`RValueToCFormString[expr] <> ";";
+           result = result <> CConversion`RValueToCFormString[Parameters`DecreaseIndexLiterals[expr]] <> ";";
            Return[result];
           ];