Skip to content

Commit

Permalink
rename m2 parameter to mu2 for compatibility with SARAH
Browse files Browse the repository at this point in the history
  • Loading branch information
Alexander Voigt authored and Alexander Voigt committed Sep 11, 2015
1 parent b90e276 commit f166196
Show file tree
Hide file tree
Showing 3 changed files with 57 additions and 57 deletions.
108 changes: 54 additions & 54 deletions model_files/SplitMSSM/FlexibleSUSY.m.in
View file
Expand Up @@ -18,15 +18,15 @@ EXTPAR = {
};

FSExtraInputParameters = {
{mq2, MSQ2IN, {3,3}}, (* 3x3 matrix *)
{mu2, MSU2IN, {3,3}}, (* 3x3 matrix *)
{md2, MSD2IN, {3,3}}, (* 3x3 matrix *)
{ml2, MSL2IN, {3,3}}, (* 3x3 matrix *)
{me2, MSE2IN, {3,3}}, (* 3x3 matrix *)
{TYu, TUIN , {3,3}} (* 3x3 matrix *)
{msq2, MSQ2IN, {3,3}}, (* 3x3 matrix *)
{msu2, MSU2IN, {3,3}}, (* 3x3 matrix *)
{msd2, MSD2IN, {3,3}}, (* 3x3 matrix *)
{msl2, MSL2IN, {3,3}}, (* 3x3 matrix *)
{mse2, MSE2IN, {3,3}}, (* 3x3 matrix *)
{TYu , TUIN , {3,3}} (* 3x3 matrix *)
};

EWSBOutputParameters = { m2 };
EWSBOutputParameters = { mu2 };

(* The high scale where we match to the MSSM *)
HighScale = MSUSY;
Expand Down Expand Up @@ -99,10 +99,10 @@ SARAH`ParameterDefinitions = Append[
(* Taken from arXiv:1407.4081 *)

(* abbreviations *)
xQU = Sqrt[Abs[mq2[3,3]/mu2[3,3]]];
xQU = Sqrt[Abs[msq2[3,3]/msu2[3,3]]];
At = TYu[3,3]/Yu[3,3];
xt = At - \[Mu]/TanBeta;
xtt = xt^2/Sqrt[Abs[mq2[3,3] mu2[3,3]]];
xtt = xt^2/Sqrt[Abs[msq2[3,3] msu2[3,3]]];

(* arXiv:1407.4081, Eq. (3) *)
lambdaTree = 1/4 (g2^2 + 3/5 g1^2) Cos[2 ArcTan[TanBeta]]^2;
Expand All @@ -115,33 +115,33 @@ lambda1LReg = 1/(4 Pi)^2 (

(* arXiv:1407.4081, Eq. (10) *)
lambda1LPhi = 1/(4 Pi)^2 (
3 Yu[3,3]^2 (Yu[3,3]^2 + 0.5 (g2^2-g1^2/5) Cos[2 ArcTan[TanBeta]]) Log[mq2[3,3]/SCALE^2]
+ 3 Yu[3,3]^2 (Yu[3,3]^2 + 0.4 g1^2 Cos[2 ArcTan[TanBeta]]) Log[mu2[3,3]/SCALE^2]
3 Yu[3,3]^2 (Yu[3,3]^2 + 0.5 (g2^2-g1^2/5) Cos[2 ArcTan[TanBeta]]) Log[msq2[3,3]/SCALE^2]
+ 3 Yu[3,3]^2 (Yu[3,3]^2 + 0.4 g1^2 Cos[2 ArcTan[TanBeta]]) Log[msu2[3,3]/SCALE^2]
+ Cos[2 ArcTan[TanBeta]]^2/300 (
3 (g1^4 + 25 g2^4) (
+ Log[mq2[1,1]/SCALE^2]
+ Log[mq2[2,2]/SCALE^2]
+ Log[mq2[3,3]/SCALE^2]
+ Log[msq2[1,1]/SCALE^2]
+ Log[msq2[2,2]/SCALE^2]
+ Log[msq2[3,3]/SCALE^2]
)
+ 24 g1^4 (
+ Log[mu2[1,1]/SCALE^2]
+ Log[mu2[2,2]/SCALE^2]
+ Log[mu2[3,3]/SCALE^2]
+ Log[msu2[1,1]/SCALE^2]
+ Log[msu2[2,2]/SCALE^2]
+ Log[msu2[3,3]/SCALE^2]
)
+ 6 g1^4 (
+ Log[md2[1,1]/SCALE^2]
+ Log[md2[2,2]/SCALE^2]
+ Log[md2[3,3]/SCALE^2]
+ Log[msd2[1,1]/SCALE^2]
+ Log[msd2[2,2]/SCALE^2]
+ Log[msd2[3,3]/SCALE^2]
)
+ (9 g1^4 + 25 g2^4) (
+ Log[ml2[1,1]/SCALE^2]
+ Log[ml2[2,2]/SCALE^2]
+ Log[ml2[3,3]/SCALE^2]
+ Log[msl2[1,1]/SCALE^2]
+ Log[msl2[2,2]/SCALE^2]
+ Log[msl2[3,3]/SCALE^2]
)
+ 18 g1^4 (
+ Log[me2[1,1]/SCALE^2]
+ Log[me2[2,2]/SCALE^2]
+ Log[me2[3,3]/SCALE^2]
+ Log[mse2[1,1]/SCALE^2]
+ Log[mse2[2,2]/SCALE^2]
+ Log[mse2[3,3]/SCALE^2]
)
)
+ 1/4800 (261 g1^4 + 630 g1^2 g2^2 + 1325 g2^4
Expand All @@ -164,17 +164,17 @@ lambda2LPhi = - (g3^2 Yu[3,3]^4)/(32 Pi^4) (
3
+ 4 Log[xQU]
+ 8 Log[xQU]^2
+ 6 (Log[mq2[3,3]/SCALE^2])^2
- 4 (1 + 3 Log[xQU]) Log[mq2[3,3]/SCALE^2]
+ 6 (Log[msq2[3,3]/SCALE^2])^2
- 4 (1 + 3 Log[xQU]) Log[msq2[3,3]/SCALE^2]
+ xtt (
+ (12 xQU Log[xQU])/(xQU^2 - 1) (2 Log[mq2[3,3]/SCALE^2] - 1)
+ (12 xQU Log[xQU])/(xQU^2 - 1) (2 Log[msq2[3,3]/SCALE^2] - 1)
- (16 xQU (xQU^2 - 2) Log[xQU]^2)/(xQU^2 - 1)^2
)
+ xtt^2 (
+ (6 xQU^2 (5 + xQU^2) Log[xQU])/(xQU^2 - 1)^3
+ (4 xQU^2 (xQU^4 - 4 xQU^2 - 5) (Log[xQU])^2)/(xQU^2 - 1)^4
- (10 xQU^2)/(xQU^2 - 1)^2
+ (12 xQU^2)/(xQU^2 - 1)^2 (1 - (xQU^2 + 1)/(xQU^2 - 1) Log[xQU]) Log[mq2[3,3]/SCALE^2]
+ (12 xQU^2)/(xQU^2 - 1)^2 (1 - (xQU^2 + 1)/(xQU^2 - 1) Log[xQU]) Log[msq2[3,3]/SCALE^2]
)
);

Expand All @@ -191,15 +191,15 @@ gYu1L = 1/(4 Pi)^2 g1 Sqrt[3/5] Sin[ArcTan[TanBeta]] (
+ (9 Yu[3,3]^2)/(4 Sin[ArcTan[TanBeta]]^2)
+ (4 g1^2 - 9 (g1^2 + 5 g2^2) Cos[ArcTan[TanBeta]]^2)/40 Log[mAInput^2/SCALE^2]
+ (g1^2/10) (
+ Log[ml2[1,1]/SCALE^2] + Log[ml2[2,2]/SCALE^2] + Log[ml2[3,3]/SCALE^2]
+ 2 Log[me2[1,1]/SCALE^2] + 2 Log[me2[2,2]/SCALE^2] + 2 Log[me2[3,3]/SCALE^2]
+ Log[msl2[1,1]/SCALE^2] + Log[msl2[2,2]/SCALE^2] + Log[msl2[3,3]/SCALE^2]
+ 2 Log[mse2[1,1]/SCALE^2] + 2 Log[mse2[2,2]/SCALE^2] + 2 Log[mse2[3,3]/SCALE^2]
)
+ (g1^2/30) (
+ Log[mq2[1,1]/SCALE^2] + Log[mq2[2,2]/SCALE^2] + Log[mq2[3,3]/SCALE^2]
+ 8 Log[mu2[1,1]/SCALE^2] + 8 Log[mu2[2,2]/SCALE^2] + 8 Log[mu2[3,3]/SCALE^2]
+ 2 Log[md2[1,1]/SCALE^2] + 2 Log[md2[2,2]/SCALE^2] + 2 Log[md2[3,3]/SCALE^2]
+ Log[msq2[1,1]/SCALE^2] + Log[msq2[2,2]/SCALE^2] + Log[msq2[3,3]/SCALE^2]
+ 8 Log[msu2[1,1]/SCALE^2] + 8 Log[msu2[2,2]/SCALE^2] + 8 Log[msu2[3,3]/SCALE^2]
+ 2 Log[msd2[1,1]/SCALE^2] + 2 Log[msd2[2,2]/SCALE^2] + 2 Log[msd2[3,3]/SCALE^2]
)
+ Yu[3,3]^2/(4 Sin[ArcTan[TanBeta]]^2) (7 Log[mq2[3,3]/SCALE^2] - 13 Log[mu2[3,3]/SCALE^2])
+ Yu[3,3]^2/(4 Sin[ArcTan[TanBeta]]^2) (7 Log[msq2[3,3]/SCALE^2] - 13 Log[msu2[3,3]/SCALE^2])
);

(* arXiv:1407.4081, Eq. (18) *)
Expand All @@ -208,13 +208,13 @@ gYd1L = 1/(4 Pi)^2 g1 Sqrt[3/5] Cos[ArcTan[TanBeta]] (
+ (3 g1^2/80) (-44 + 7 Sin[ArcTan[TanBeta]]^2)
+ (4 g1^2 - 9 (g1^2 + 5 g2^2) Sin[ArcTan[TanBeta]]^2)/40 Log[mAInput^2/SCALE^2]
+ (g1^2/10) (
+ Log[ml2[1,1]/SCALE^2] + Log[ml2[2,2]/SCALE^2] + Log[ml2[3,3]/SCALE^2]
+ 2 Log[me2[1,1]/SCALE^2] + 2 Log[me2[2,2]/SCALE^2] + 2 Log[me2[3,3]/SCALE^2]
+ Log[msl2[1,1]/SCALE^2] + Log[msl2[2,2]/SCALE^2] + Log[msl2[3,3]/SCALE^2]
+ 2 Log[mse2[1,1]/SCALE^2] + 2 Log[mse2[2,2]/SCALE^2] + 2 Log[mse2[3,3]/SCALE^2]
)
+ (g1^2/30) (
+ Log[mq2[1,1]/SCALE^2] + Log[mq2[2,2]/SCALE^2] + Log[mq2[3,3]/SCALE^2]
+ 8 Log[mu2[1,1]/SCALE^2] + 8 Log[mu2[2,2]/SCALE^2] + 8 Log[mu2[3,3]/SCALE^2]
+ 2 Log[md2[1,1]/SCALE^2] + 2 Log[md2[2,2]/SCALE^2] + 2 Log[md2[3,3]/SCALE^2]
+ Log[msq2[1,1]/SCALE^2] + Log[msq2[2,2]/SCALE^2] + Log[msq2[3,3]/SCALE^2]
+ 8 Log[msu2[1,1]/SCALE^2] + 8 Log[msu2[2,2]/SCALE^2] + 8 Log[msu2[3,3]/SCALE^2]
+ 2 Log[msd2[1,1]/SCALE^2] + 2 Log[msd2[2,2]/SCALE^2] + 2 Log[msd2[3,3]/SCALE^2]
)
);

Expand All @@ -225,27 +225,27 @@ g2u1L = 1/(4 Pi)^2 g2 Sin[ArcTan[TanBeta]] (
+ 9 Yu[3,3]^2/(4 Sin[ArcTan[TanBeta]]^2)
+ (20 g2^2 + 3 (-9 g1^2 + 35 g2^2) Cos[ArcTan[TanBeta]]^2)/120 Log[mAInput^2/SCALE^2]
+ g2^2/6 (
+ Log[ml2[1,1]/SCALE^2]
+ Log[ml2[2,2]/SCALE^2]
+ Log[ml2[3,3]/SCALE^2]
+ Log[msl2[1,1]/SCALE^2]
+ Log[msl2[2,2]/SCALE^2]
+ Log[msl2[3,3]/SCALE^2]
)
+ g2^2/2 (
+ Log[mq2[1,1]/SCALE^2]
+ Log[mq2[2,2]/SCALE^2]
+ Log[mq2[3,3]/SCALE^2]
+ Log[msq2[1,1]/SCALE^2]
+ Log[msq2[2,2]/SCALE^2]
+ Log[msq2[3,3]/SCALE^2]
)
- 0.75 Yu[3,3]^2/Sin[ArcTan[TanBeta]]^2 (3 Log[mq2[3,3]/SCALE^2] - Log[mu2[3,3]/SCALE^2])
- 0.75 Yu[3,3]^2/Sin[ArcTan[TanBeta]]^2 (3 Log[msq2[3,3]/SCALE^2] - Log[msu2[3,3]/SCALE^2])
);

(* arXiv:1407.4081, Eq. (16) *)
g2d1L = 1/(4 Pi)^2 g2 Cos[ArcTan[TanBeta]] (
- g2^2 (2/3 + 11/16 Sin[ArcTan[TanBeta]]^2)
+ (3 g1^2/80) (-2 + 7 Sin[ArcTan[TanBeta]]^2)
+ (g2^2/2) Log[mq2[1,1]/SCALE^2]
+ (g2^2/2) Log[mq2[2,2]/SCALE^2]
+ (g2^2/2) Log[mq2[3,3]/SCALE^2]
+ (g2^2/2) Log[msq2[1,1]/SCALE^2]
+ (g2^2/2) Log[msq2[2,2]/SCALE^2]
+ (g2^2/2) Log[msq2[3,3]/SCALE^2]
+ (20 g2^2 + 3 (-9 g1^2 + 35 g2^2) Sin[ArcTan[TanBeta]]^2)/120 Log[mAInput^2/SCALE^2]
+ (g2^2/6) Log[ml2[1,1]/SCALE^2]
+ (g2^2/6) Log[ml2[2,2]/SCALE^2]
+ (g2^2/6) Log[ml2[3,3]/SCALE^2]
+ (g2^2/6) Log[msl2[1,1]/SCALE^2]
+ (g2^2/6) Log[msl2[2,2]/SCALE^2]
+ (g2^2/6) Log[msl2[3,3]/SCALE^2]
);
2 changes: 1 addition & 1 deletion sarah/SplitMSSM/SplitMSSM.m
View file
Expand Up @@ -55,7 +55,7 @@
{LagSplit, {AddHC->True }}
};

LagNoHC = m2 conj[H].H - 1/2 \[Lambda] conj[H].H.conj[H].H;
LagNoHC = mu2 conj[H].H - 1/2 \[Lambda] conj[H].H.conj[H].H;

LagHC = - Yd conj[H].d.q - Ye conj[H].e.l + Yu H.u.q;

Expand Down
4 changes: 2 additions & 2 deletions sarah/SplitMSSM/parameters.m
View file
Expand Up @@ -18,9 +18,9 @@
DependenceNum -> Sqrt[2]/v* {{Mass[Fe,1],0,0},
{0, Mass[Fe,2],0},
{0, 0, Mass[Fe,3]}}}},
{m2, { Description -> "SM Mu Parameter",
{mu2, { Description -> "SM Mu Parameter",
LaTeX -> "m^2",
OutputName -> m2 }},
OutputName -> mu2 }},
{\[Lambda], { Description -> "SM Higgs Selfcouplings",
DependenceNum -> Mass[hh]^2/(2 v^2)}},
{v, { Description -> "EW-VEV",
Expand Down

0 comments on commit f166196

Please sign in to comment.