Skip to content

Commit

Permalink
use input parameter to specify high-scale matching conndition loop
Browse files Browse the repository at this point in the history
order
  • Loading branch information
Alexander Voigt authored and Alexander Voigt committed Jan 27, 2017
1 parent 2dc07b9 commit 3e6a7a3
Show file tree
Hide file tree
Showing 10 changed files with 109 additions and 107 deletions.
28 changes: 14 additions & 14 deletions model_files/HGTHDMIIMSSMBC/FlexibleSUSY.m.in
Expand Up @@ -36,26 +36,26 @@ HighScaleInput = {
{g2u , g2 },
{g2up , GUTNormalization[g1] g1 (* gY *)},
{Lambda1, 1/2 (1/4 ( (GUTNormalization[g1] g1)^2 + g2^2)
+ UnitStep[THRESHOLD-1] UnitStep[LambdaLoopOrder-1] (deltaLambda1th1L + deltaLambda1Phi1L)
+ UnitStep[THRESHOLD-2] UnitStep[LambdaLoopOrder-2] deltaLambda1th2L)},
+ UnitStep[LambdaLoopOrder-1] (deltaLambda1th1L + deltaLambda1Phi1L)
+ UnitStep[LambdaLoopOrder-2] deltaLambda1th2L)},
{Lambda2, 1/2 (1/4 ( (GUTNormalization[g1] g1)^2 + g2^2)
+ UnitStep[THRESHOLD-1] UnitStep[LambdaLoopOrder-1] (deltaLambda2th1L + deltaLambda2Phi1L)
+ UnitStep[THRESHOLD-2] UnitStep[LambdaLoopOrder-2] deltaLambda2th2L)},
+ UnitStep[LambdaLoopOrder-1] (deltaLambda2th1L + deltaLambda2Phi1L)
+ UnitStep[LambdaLoopOrder-2] deltaLambda2th2L)},
{Lambda3, 1/4 (-(GUTNormalization[g1] g1)^2 + g2^2)
+ UnitStep[THRESHOLD-1] UnitStep[LambdaLoopOrder-1] (deltaLambda3th1L + deltaLambda3Phi1L)
+ UnitStep[THRESHOLD-2] UnitStep[LambdaLoopOrder-2] deltaLambda3th2L},
+ UnitStep[LambdaLoopOrder-1] (deltaLambda3th1L + deltaLambda3Phi1L)
+ UnitStep[LambdaLoopOrder-2] deltaLambda3th2L},
{Lambda4, -1/2 g2^2
+ UnitStep[THRESHOLD-1] UnitStep[LambdaLoopOrder-1] (deltaLambda4th1L + deltaLambda4Phi1L)
+ UnitStep[THRESHOLD-2] UnitStep[LambdaLoopOrder-2] deltaLambda4th2L},
+ UnitStep[LambdaLoopOrder-1] (deltaLambda4th1L + deltaLambda4Phi1L)
+ UnitStep[LambdaLoopOrder-2] deltaLambda4th2L},
{Lambda5, 0
+ UnitStep[THRESHOLD-1] UnitStep[LambdaLoopOrder-1] (deltaLambda5th1L + deltaLambda5Phi1L)
+ UnitStep[THRESHOLD-2] UnitStep[LambdaLoopOrder-2] deltaLambda5th2L},
+ UnitStep[LambdaLoopOrder-1] (deltaLambda5th1L + deltaLambda5Phi1L)
+ UnitStep[LambdaLoopOrder-2] deltaLambda5th2L},
{Lambda6, 0
+ UnitStep[THRESHOLD-1] UnitStep[LambdaLoopOrder-1] (deltaLambda6th1L + deltaLambda6Phi1L)
+ UnitStep[THRESHOLD-2] UnitStep[LambdaLoopOrder-2] deltaLambda6th2L},
+ UnitStep[LambdaLoopOrder-1] (deltaLambda6th1L + deltaLambda6Phi1L)
+ UnitStep[LambdaLoopOrder-2] deltaLambda6th2L},
{Lambda7, 0
+ UnitStep[THRESHOLD-1] UnitStep[LambdaLoopOrder-1] (deltaLambda7th1L + deltaLambda7Phi1L)
+ UnitStep[THRESHOLD-2] UnitStep[LambdaLoopOrder-2] deltaLambda7th2L}
+ UnitStep[LambdaLoopOrder-1] (deltaLambda7th1L + deltaLambda7Phi1L)
+ UnitStep[LambdaLoopOrder-2] deltaLambda7th2L}
};

(* The scale where we impose the EWSB conditions
Expand Down
28 changes: 14 additions & 14 deletions model_files/HGTHDMIIMSSMBCApprox/FlexibleSUSY.m.in
Expand Up @@ -96,26 +96,26 @@ HighScaleInput = {
{g2u , g2 },
{g2up , GUTNormalization[g1] g1 (* gY *)},
{Lambda1, 1/2 (deltaLambda0L[[1]]
+ UnitStep[THRESHOLD-1] UnitStep[LambdaLoopOrder-1] deltaLambda1L[[1]]
+ UnitStep[THRESHOLD-2] UnitStep[LambdaLoopOrder-2] deltaLambda1th2L)},
+ UnitStep[LambdaLoopOrder-1] deltaLambda1L[[1]]
+ UnitStep[LambdaLoopOrder-2] deltaLambda1th2L)},
{Lambda2, 1/2 (deltaLambda0L[[2]]
+ UnitStep[THRESHOLD-1] UnitStep[LambdaLoopOrder-1] deltaLambda1L[[2]]
+ UnitStep[THRESHOLD-2] UnitStep[LambdaLoopOrder-2] deltaLambda2th2L)},
+ UnitStep[LambdaLoopOrder-1] deltaLambda1L[[2]]
+ UnitStep[LambdaLoopOrder-2] deltaLambda2th2L)},
{Lambda3, deltaLambda0L[[3]]
+ UnitStep[THRESHOLD-1] UnitStep[LambdaLoopOrder-1] deltaLambda1L[[3]]
+ UnitStep[THRESHOLD-2] UnitStep[LambdaLoopOrder-2] deltaLambda3th2L},
+ UnitStep[LambdaLoopOrder-1] deltaLambda1L[[3]]
+ UnitStep[LambdaLoopOrder-2] deltaLambda3th2L},
{Lambda4, deltaLambda0L[[4]]
+ UnitStep[THRESHOLD-1] UnitStep[LambdaLoopOrder-1] deltaLambda1L[[4]]
+ UnitStep[THRESHOLD-2] UnitStep[LambdaLoopOrder-2] deltaLambda4th2L},
+ UnitStep[LambdaLoopOrder-1] deltaLambda1L[[4]]
+ UnitStep[LambdaLoopOrder-2] deltaLambda4th2L},
{Lambda5, deltaLambda0L[[5]]
+ UnitStep[THRESHOLD-1] UnitStep[LambdaLoopOrder-1] deltaLambda1L[[5]]
+ UnitStep[THRESHOLD-2] UnitStep[LambdaLoopOrder-2] deltaLambda5th2L},
+ UnitStep[LambdaLoopOrder-1] deltaLambda1L[[5]]
+ UnitStep[LambdaLoopOrder-2] deltaLambda5th2L},
{Lambda6, deltaLambda0L[[6]]
+ UnitStep[THRESHOLD-1] UnitStep[LambdaLoopOrder-1] deltaLambda1L[[6]]
+ UnitStep[THRESHOLD-2] UnitStep[LambdaLoopOrder-2] deltaLambda6th2L},
+ UnitStep[LambdaLoopOrder-1] deltaLambda1L[[6]]
+ UnitStep[LambdaLoopOrder-2] deltaLambda6th2L},
{Lambda7, deltaLambda0L[[7]]
+ UnitStep[THRESHOLD-1] UnitStep[LambdaLoopOrder-1] deltaLambda1L[[7]]
+ UnitStep[THRESHOLD-2] UnitStep[LambdaLoopOrder-2] deltaLambda7th2L}
+ UnitStep[LambdaLoopOrder-1] deltaLambda1L[[7]]
+ UnitStep[LambdaLoopOrder-2] deltaLambda7th2L}
};

(* The scale where we impose the EWSB conditions
Expand Down
28 changes: 14 additions & 14 deletions model_files/HGTHDMIIMSSMBCFull/FlexibleSUSY.m.in
Expand Up @@ -86,26 +86,26 @@ HighScaleInput = {
{g2u , g2 },
{g2up , GUTNormalization[g1] g1 (* gY *)},
{Lambda1, 1/2 (deltaLambda0L[[1]]
+ UnitStep[THRESHOLD-1] UnitStep[LambdaLoopOrder-1] deltaLambda1L[[1]]
+ UnitStep[THRESHOLD-2] UnitStep[LambdaLoopOrder-2] deltaLambda1th2L)},
+ UnitStep[LambdaLoopOrder-1] deltaLambda1L[[1]]
+ UnitStep[LambdaLoopOrder-2] deltaLambda1th2L)},
{Lambda2, 1/2 (deltaLambda0L[[2]]
+ UnitStep[THRESHOLD-1] UnitStep[LambdaLoopOrder-1] deltaLambda1L[[2]]
+ UnitStep[THRESHOLD-2] UnitStep[LambdaLoopOrder-2] deltaLambda2th2L)},
+ UnitStep[LambdaLoopOrder-1] deltaLambda1L[[2]]
+ UnitStep[LambdaLoopOrder-2] deltaLambda2th2L)},
{Lambda3, deltaLambda0L[[3]]
+ UnitStep[THRESHOLD-1] UnitStep[LambdaLoopOrder-1] deltaLambda1L[[3]]
+ UnitStep[THRESHOLD-2] UnitStep[LambdaLoopOrder-2] deltaLambda3th2L},
+ UnitStep[LambdaLoopOrder-1] deltaLambda1L[[3]]
+ UnitStep[LambdaLoopOrder-2] deltaLambda3th2L},
{Lambda4, deltaLambda0L[[4]]
+ UnitStep[THRESHOLD-1] UnitStep[LambdaLoopOrder-1] deltaLambda1L[[4]]
+ UnitStep[THRESHOLD-2] UnitStep[LambdaLoopOrder-2] deltaLambda4th2L},
+ UnitStep[LambdaLoopOrder-1] deltaLambda1L[[4]]
+ UnitStep[LambdaLoopOrder-2] deltaLambda4th2L},
{Lambda5, deltaLambda0L[[5]]
+ UnitStep[THRESHOLD-1] UnitStep[LambdaLoopOrder-1] deltaLambda1L[[5]]
+ UnitStep[THRESHOLD-2] UnitStep[LambdaLoopOrder-2] deltaLambda5th2L},
+ UnitStep[LambdaLoopOrder-1] deltaLambda1L[[5]]
+ UnitStep[LambdaLoopOrder-2] deltaLambda5th2L},
{Lambda6, deltaLambda0L[[6]]
+ UnitStep[THRESHOLD-1] UnitStep[LambdaLoopOrder-1] deltaLambda1L[[6]]
+ UnitStep[THRESHOLD-2] UnitStep[LambdaLoopOrder-2] deltaLambda6th2L},
+ UnitStep[LambdaLoopOrder-1] deltaLambda1L[[6]]
+ UnitStep[LambdaLoopOrder-2] deltaLambda6th2L},
{Lambda7, deltaLambda0L[[7]]
+ UnitStep[THRESHOLD-1] UnitStep[LambdaLoopOrder-1] deltaLambda1L[[7]]
+ UnitStep[THRESHOLD-2] UnitStep[LambdaLoopOrder-2] deltaLambda7th2L}
+ UnitStep[LambdaLoopOrder-1] deltaLambda1L[[7]]
+ UnitStep[LambdaLoopOrder-2] deltaLambda7th2L}
};

(* The scale where we impose the EWSB conditions
Expand Down
4 changes: 2 additions & 2 deletions model_files/HSSUSY/FlexibleSUSY.m.in
Expand Up @@ -36,11 +36,11 @@ HighScaleFirstGuess = MSUSY;

HighScaleInput = {
{\[Lambda], lambdaTree
+ UnitStep[THRESHOLD-1] UnitStep[LambdaLoopOrder-1] (
+ UnitStep[LambdaLoopOrder-1] (
+ lambda1LReg + lambda1LPhi
+ lambda1LChi1 + lambda1LChi2
+ lambda1Lbottom + lambda1Ltau)
+ UnitStep[THRESHOLD-2] UnitStep[LambdaLoopOrder-2] (
+ UnitStep[LambdaLoopOrder-2] (
Which[IsCloseRel[SCALE^2,msq2[3,3],0.01] &&
IsCloseRel[SCALE^2,msu2[3,3],0.01] &&
IsCloseRel[SCALE,M3Input,0.01],
Expand Down
28 changes: 14 additions & 14 deletions model_files/HTHDMIIMSSMBC/FlexibleSUSY.m.in
Expand Up @@ -29,26 +29,26 @@ HighScaleFirstGuess = MSUSY;

HighScaleInput = {
{Lambda1, 1/2 (1/4 ( (GUTNormalization[g1] g1)^2 + g2^2)
+ UnitStep[THRESHOLD-1] UnitStep[LambdaLoopOrder-1] (deltaLambda1th1L + deltaLambda1Phi1L)
+ UnitStep[THRESHOLD-2] UnitStep[LambdaLoopOrder-2] deltaLambda1th2L)},
+ UnitStep[LambdaLoopOrder-1] (deltaLambda1th1L + deltaLambda1Phi1L)
+ UnitStep[LambdaLoopOrder-2] deltaLambda1th2L)},
{Lambda2, 1/2 (1/4 ( (GUTNormalization[g1] g1)^2 + g2^2)
+ UnitStep[THRESHOLD-1] UnitStep[LambdaLoopOrder-1] (deltaLambda2th1L + deltaLambda2Phi1L)
+ UnitStep[THRESHOLD-2] UnitStep[LambdaLoopOrder-2] deltaLambda2th2L)},
+ UnitStep[LambdaLoopOrder-1] (deltaLambda2th1L + deltaLambda2Phi1L)
+ UnitStep[LambdaLoopOrder-2] deltaLambda2th2L)},
{Lambda3, 1/4 (-(GUTNormalization[g1] g1)^2 + g2^2)
+ UnitStep[THRESHOLD-1] UnitStep[LambdaLoopOrder-1] (deltaLambda3th1L + deltaLambda3Phi1L)
+ UnitStep[THRESHOLD-2] UnitStep[LambdaLoopOrder-2] deltaLambda3th2L},
+ UnitStep[LambdaLoopOrder-1] (deltaLambda3th1L + deltaLambda3Phi1L)
+ UnitStep[LambdaLoopOrder-2] deltaLambda3th2L},
{Lambda4, -1/2 g2^2
+ UnitStep[THRESHOLD-1] UnitStep[LambdaLoopOrder-1] (deltaLambda4th1L + deltaLambda4Phi1L)
+ UnitStep[THRESHOLD-2] UnitStep[LambdaLoopOrder-2] deltaLambda4th2L},
+ UnitStep[LambdaLoopOrder-1] (deltaLambda4th1L + deltaLambda4Phi1L)
+ UnitStep[LambdaLoopOrder-2] deltaLambda4th2L},
{Lambda5, 0
+ UnitStep[THRESHOLD-1] UnitStep[LambdaLoopOrder-1] (deltaLambda5th1L + deltaLambda5Phi1L)
+ UnitStep[THRESHOLD-2] UnitStep[LambdaLoopOrder-2] deltaLambda5th2L},
+ UnitStep[LambdaLoopOrder-1] (deltaLambda5th1L + deltaLambda5Phi1L)
+ UnitStep[LambdaLoopOrder-2] deltaLambda5th2L},
{Lambda6, 0
+ UnitStep[THRESHOLD-1] UnitStep[LambdaLoopOrder-1] (deltaLambda6th1L + deltaLambda6Phi1L)
+ UnitStep[THRESHOLD-2] UnitStep[LambdaLoopOrder-2] deltaLambda6th2L},
+ UnitStep[LambdaLoopOrder-1] (deltaLambda6th1L + deltaLambda6Phi1L)
+ UnitStep[LambdaLoopOrder-2] deltaLambda6th2L},
{Lambda7, 0
+ UnitStep[THRESHOLD-1] UnitStep[LambdaLoopOrder-1] (deltaLambda7th1L + deltaLambda7Phi1L)
+ UnitStep[THRESHOLD-2] UnitStep[LambdaLoopOrder-2] deltaLambda7th2L}
+ UnitStep[LambdaLoopOrder-1] (deltaLambda7th1L + deltaLambda7Phi1L)
+ UnitStep[LambdaLoopOrder-2] deltaLambda7th2L}
};

(* The scale where we impose the EWSB conditions
Expand Down
15 changes: 8 additions & 7 deletions model_files/SplitMSSM/FlexibleSUSY.m.in
Expand Up @@ -15,7 +15,8 @@ EXTPAR = {
{5, mAInput},
{6, MEWSB},
{7, AtInput},
{25, TanBeta}
{25, TanBeta},
{100, LambdaLoopOrder}
};

FSExtraInputParameters = {
Expand All @@ -35,14 +36,14 @@ HighScaleFirstGuess = MSUSY;

HighScaleInput = {
{\[Lambda], lambdaTree
+ UnitStep[THRESHOLD-1] (lambda1LReg + lambda1LPhi)
+ UnitStep[THRESHOLD-2] (
+ UnitStep[LambdaLoopOrder-1] (lambda1LReg + lambda1LPhi)
+ UnitStep[LambdaLoopOrder-2] (
If[IsCloseRel[msq2[3,3],msu2[3,3],0.01],
lambda2LPhiMSqMSuDegenerate, lambda2LPhiFull])},
{gYu , gYuTree + UnitStep[THRESHOLD-1] gYu1L},
{gYd , gYdTree + UnitStep[THRESHOLD-1] gYd1L},
{g2u , g2uTree + UnitStep[THRESHOLD-1] g2u1L},
{g2d , g2dTree + UnitStep[THRESHOLD-1] g2d1L}
{gYu , gYuTree + UnitStep[LambdaLoopOrder-1] gYu1L},
{gYd , gYdTree + UnitStep[LambdaLoopOrder-1] gYd1L},
{g2u , g2uTree + UnitStep[LambdaLoopOrder-1] g2u1L},
{g2d , g2dTree + UnitStep[LambdaLoopOrder-1] g2d1L}
};

(* The scale where we impose the EWSB conditions
Expand Down
1 change: 1 addition & 0 deletions model_files/SplitMSSM/LesHouches.in.SplitMSSM
Expand Up @@ -43,6 +43,7 @@ Block EXTPAR # Input parameters
6 173.34 # MEWSB
7 0 # At(MSUSY)
25 10 # TanBeta(MSUSY)
100 2 # LambdaLoopOrder
Block MSQ2IN
1 1 1.00000000E+06 # mq2(1,1)
2 2 1.00000000E+06 # mq2(2,2)
Expand Down
28 changes: 14 additions & 14 deletions model_files/THDMIIMSSMBC/FlexibleSUSY.m.in
Expand Up @@ -29,26 +29,26 @@ HighScaleFirstGuess = MSUSY;

HighScaleInput = {
{Lambda1, 1/2 (1/4 ( (GUTNormalization[g1] g1)^2 + g2^2)
+ UnitStep[THRESHOLD-1] UnitStep[LambdaLoopOrder-1] (deltaLambda1th1L + deltaLambda1Phi1L)
+ UnitStep[THRESHOLD-2] UnitStep[LambdaLoopOrder-2] deltaLambda1th2L)},
+ UnitStep[LambdaLoopOrder-1] (deltaLambda1th1L + deltaLambda1Phi1L)
+ UnitStep[LambdaLoopOrder-2] deltaLambda1th2L)},
{Lambda2, 1/2 (1/4 ( (GUTNormalization[g1] g1)^2 + g2^2)
+ UnitStep[THRESHOLD-1] UnitStep[LambdaLoopOrder-1] (deltaLambda2th1L + deltaLambda2Phi1L)
+ UnitStep[THRESHOLD-2] UnitStep[LambdaLoopOrder-2] deltaLambda2th2L)},
+ UnitStep[LambdaLoopOrder-1] (deltaLambda2th1L + deltaLambda2Phi1L)
+ UnitStep[LambdaLoopOrder-2] deltaLambda2th2L)},
{Lambda3, 1/4 (-(GUTNormalization[g1] g1)^2 + g2^2)
+ UnitStep[THRESHOLD-1] UnitStep[LambdaLoopOrder-1] (deltaLambda3th1L + deltaLambda3Phi1L)
+ UnitStep[THRESHOLD-2] UnitStep[LambdaLoopOrder-2] deltaLambda3th2L},
+ UnitStep[LambdaLoopOrder-1] (deltaLambda3th1L + deltaLambda3Phi1L)
+ UnitStep[LambdaLoopOrder-2] deltaLambda3th2L},
{Lambda4, -1/2 g2^2
+ UnitStep[THRESHOLD-1] UnitStep[LambdaLoopOrder-1] (deltaLambda4th1L + deltaLambda4Phi1L)
+ UnitStep[THRESHOLD-2] UnitStep[LambdaLoopOrder-2] deltaLambda4th2L},
+ UnitStep[LambdaLoopOrder-1] (deltaLambda4th1L + deltaLambda4Phi1L)
+ UnitStep[LambdaLoopOrder-2] deltaLambda4th2L},
{Lambda5, 0
+ UnitStep[THRESHOLD-1] UnitStep[LambdaLoopOrder-1] (deltaLambda5th1L + deltaLambda5Phi1L)
+ UnitStep[THRESHOLD-2] UnitStep[LambdaLoopOrder-2] deltaLambda5th2L},
+ UnitStep[LambdaLoopOrder-1] (deltaLambda5th1L + deltaLambda5Phi1L)
+ UnitStep[LambdaLoopOrder-2] deltaLambda5th2L},
{Lambda6, 0
+ UnitStep[THRESHOLD-1] UnitStep[LambdaLoopOrder-1] (deltaLambda6th1L + deltaLambda6Phi1L)
+ UnitStep[THRESHOLD-2] UnitStep[LambdaLoopOrder-2] deltaLambda6th2L},
+ UnitStep[LambdaLoopOrder-1] (deltaLambda6th1L + deltaLambda6Phi1L)
+ UnitStep[LambdaLoopOrder-2] deltaLambda6th2L},
{Lambda7, 0
+ UnitStep[THRESHOLD-1] UnitStep[LambdaLoopOrder-1] (deltaLambda7th1L + deltaLambda7Phi1L)
+ UnitStep[THRESHOLD-2] UnitStep[LambdaLoopOrder-2] deltaLambda7th2L}
+ UnitStep[LambdaLoopOrder-1] (deltaLambda7th1L + deltaLambda7Phi1L)
+ UnitStep[LambdaLoopOrder-2] deltaLambda7th2L}
};

(* The scale where we impose the EWSB conditions
Expand Down
28 changes: 14 additions & 14 deletions model_files/THDMIIMSSMBCApprox/FlexibleSUSY.m.in
Expand Up @@ -87,26 +87,26 @@ deltaLambda1L = Simplify[Plus @@ Select[List @@ Expand[#],

HighScaleInput = {
{Lambda1, 1/2 (deltaLambda0L[[1]]
+ UnitStep[THRESHOLD-1] UnitStep[LambdaLoopOrder-1] deltaLambda1L[[1]]
+ UnitStep[THRESHOLD-2] UnitStep[LambdaLoopOrder-2] deltaLambda1th2L)},
+ UnitStep[LambdaLoopOrder-1] deltaLambda1L[[1]]
+ UnitStep[LambdaLoopOrder-2] deltaLambda1th2L)},
{Lambda2, 1/2 (deltaLambda0L[[2]]
+ UnitStep[THRESHOLD-1] UnitStep[LambdaLoopOrder-1] deltaLambda1L[[2]]
+ UnitStep[THRESHOLD-2] UnitStep[LambdaLoopOrder-2] deltaLambda2th2L)},
+ UnitStep[LambdaLoopOrder-1] deltaLambda1L[[2]]
+ UnitStep[LambdaLoopOrder-2] deltaLambda2th2L)},
{Lambda3, deltaLambda0L[[3]]
+ UnitStep[THRESHOLD-1] UnitStep[LambdaLoopOrder-1] deltaLambda1L[[3]]
+ UnitStep[THRESHOLD-2] UnitStep[LambdaLoopOrder-2] deltaLambda3th2L},
+ UnitStep[LambdaLoopOrder-1] deltaLambda1L[[3]]
+ UnitStep[LambdaLoopOrder-2] deltaLambda3th2L},
{Lambda4, deltaLambda0L[[4]]
+ UnitStep[THRESHOLD-1] UnitStep[LambdaLoopOrder-1] deltaLambda1L[[4]]
+ UnitStep[THRESHOLD-2] UnitStep[LambdaLoopOrder-2] deltaLambda4th2L},
+ UnitStep[LambdaLoopOrder-1] deltaLambda1L[[4]]
+ UnitStep[LambdaLoopOrder-2] deltaLambda4th2L},
{Lambda5, deltaLambda0L[[5]]
+ UnitStep[THRESHOLD-1] UnitStep[LambdaLoopOrder-1] deltaLambda1L[[5]]
+ UnitStep[THRESHOLD-2] UnitStep[LambdaLoopOrder-2] deltaLambda5th2L},
+ UnitStep[LambdaLoopOrder-1] deltaLambda1L[[5]]
+ UnitStep[LambdaLoopOrder-2] deltaLambda5th2L},
{Lambda6, deltaLambda0L[[6]]
+ UnitStep[THRESHOLD-1] UnitStep[LambdaLoopOrder-1] deltaLambda1L[[6]]
+ UnitStep[THRESHOLD-2] UnitStep[LambdaLoopOrder-2] deltaLambda6th2L},
+ UnitStep[LambdaLoopOrder-1] deltaLambda1L[[6]]
+ UnitStep[LambdaLoopOrder-2] deltaLambda6th2L},
{Lambda7, deltaLambda0L[[7]]
+ UnitStep[THRESHOLD-1] UnitStep[LambdaLoopOrder-1] deltaLambda1L[[7]]
+ UnitStep[THRESHOLD-2] UnitStep[LambdaLoopOrder-2] deltaLambda7th2L}
+ UnitStep[LambdaLoopOrder-1] deltaLambda1L[[7]]
+ UnitStep[LambdaLoopOrder-2] deltaLambda7th2L}
};

(* The scale where we impose the EWSB conditions
Expand Down

0 comments on commit 3e6a7a3

Please sign in to comment.