use expr arXiv:1407.4081, Eq. (13) in degenerate case

FlexibleSUSY · Sep 11, 2015 · 94492a8 · 94492a8
1 parent 4d960bc
commit 94492a8
Show file tree

Hide file tree

Showing 2 changed files with 11 additions and 7 deletions.
diff --git a/model_files/SuperSplitMSSM/FlexibleSUSY.m.in b/model_files/SuperSplitMSSM/FlexibleSUSY.m.in
@@ -34,8 +34,12 @@ HighScale = MSUSY;
 HighScaleFirstGuess = MSUSY;
 
 HighScaleInput = {
-    {\[Lambda], lambdaTree + lambda1LReg + lambda1LPhi + lambda1LChi1
-                + lambda1LChi2 + lambda2LPhiHSS}
+    {\[Lambda], lambdaTree
+                + lambda1LReg + lambda1LPhi + lambda1LChi1 + lambda1LChi2
+                + If[IsClose[SCALE^2,msq2[3,3]] &&
+                     IsClose[SCALE^2,msu2[3,3]] &&
+                     IsClose[SCALE,M3Input],
+                     lambda2LPhiHSSDegenerate, lambda2LPhiHSSFull]}
 };
 
 (* The scale where we impose the EWSB conditions
@@ -179,7 +183,7 @@ lambda1LChi2 = 1/(4 Pi)^2 (
 (* lambda 2-loop threshold correction O(alpha_s alpha_t^2),
    arXiv:1407.4081, Eq. (36)
    Valid in the limit M3 = MQ[3,3] = MU[3,3] = MSUSY *)
-lambda2LPhiHSS = With[{ r = xt / SCALE },
+lambda2LPhiHSSDegenerate = With[{ r = xt / SCALE },
     (g3^2 Yu[3,3]^4)/(96 Pi^4) (
         -12 r - 6 r^2 + 14 r^3 + 1/2 r^4 - r^5
     )
@@ -188,7 +192,7 @@ lambda2LPhiHSS = With[{ r = xt / SCALE },
 (* lambda 2-loop threshold correction O(alpha_s alpha_t^2),
    arXiv:1407.4081, provided by Pietro Slavich,
    General expression. *)
-lambda2LPhiHSS = With[{
+lambda2LPhiHSSFull = With[{
     htSM = Yu[3,3] Sin[ArcTan[TanBeta]],
     CF = 4/3,
     Nc = 3,

diff --git a/model_files/SuperSplitMSSM/LesHouches.in.SuperSplitMSSM b/model_files/SuperSplitMSSM/LesHouches.in.SuperSplitMSSM
@@ -7,7 +7,7 @@ Block FlexibleSUSY
     3   1                    # calculate SM pole masses
     4   2                    # pole mass loop order
     5   2                    # EWSB loop order
-    6   2                    # beta-functions loop order
+    6   3                    # beta-functions loop order
     7   2                    # threshold corrections loop order
     8   1                    # Higgs 2-loop corrections O(alpha_t alpha_s)
     9   1                    # Higgs 2-loop corrections O(alpha_b alpha_s)
@@ -37,7 +37,7 @@ Block EXTPAR                 # Input parameters
     0   1000                 # MSUSY
     1   1000                 # M1(MSUSY)
     2   1000                 # M2(MSUSY)
-    3   1001                 # M3(MSUSY)
+    3   1000                 # M3(MSUSY)
     4   1000                 # Mu(MSUSY)
     5   1000                 # mA(MSUSY)
     6   173.34               # MEWSB
@@ -63,7 +63,7 @@ Block MSL2IN
 Block MSU2IN
   1  1     1.00000000E+06   # mu2(1,1)
   2  2     1.00000000E+06   # mu2(2,2)
-  3  3     1.00100000E+06   # mu2(3,3)
+  3  3     1.00000000E+06   # mu2(3,3)
 Block MSD2IN
   1  1     1.00000000E+06   # md2(1,1)
   2  2     1.00000000E+06   # md2(2,2)