forked from Expander/FlexibleSUSY
/
FlexibleSUSY.m.in
128 lines (112 loc) · 3.44 KB
/
FlexibleSUSY.m.in
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
FSModelName = "@CLASSNAME@";
FSEigenstates = SARAH`EWSB;
FSDefaultSARAHModel = MSSM;
OnlyLowEnergyFlexibleSUSY = True;
(* CMSSM input parameters *)
MINPAR = {
{4, Sign[\[Mu]]}
};
FSExtraInputParameters = {
{TanBeta, TanBeta, {1}},
{Ms, Ms, {1}},
{Xtt, Xtt, {1}}
};
EWSBOutputParameters = { mHd2, mHu2 };
SUSYScale = Ms;
SUSYScaleFirstGuess = Ms;
SUSYScaleInput = {
{MassB, Ms},
{MassWB, Ms},
{MassG, Ms},
{mq2, UNITMATRIX[3] Ms^2},
{mu2, UNITMATRIX[3] Ms^2},
{md2, UNITMATRIX[3] Ms^2},
{ml2, UNITMATRIX[3] Ms^2},
{me2, UNITMATRIX[3] Ms^2},
{\[Mu], Ms},
{B[\[Mu]], Sqr[Ms]/(TanBeta + 1/TanBeta)},
{T[Yu], \[Mu]/TanBeta Yu},
{T[Yd], \[Mu] TanBeta Yd},
{T[Ye], \[Mu] TanBeta Ye},
{T[Yu][3,3], (1/TanBeta + Xtt) Yu[3,3] Ms},
{vd, vu / TanBeta}
};
InitialGuessAtLowScale = {
{vu, LowEnergyConstant[vev] Sin[ArcTan[TanBeta]]},
{vd, LowEnergyConstant[vev] Cos[ArcTan[TanBeta]]},
{Yu, Automatic},
{Yd, Automatic},
{Ye, Automatic},
{MassB, Ms},
{MassWB, Ms},
{MassG, Ms},
{mq2, UNITMATRIX[3] Ms^2},
{mu2, UNITMATRIX[3] Ms^2},
{md2, UNITMATRIX[3] Ms^2},
{ml2, UNITMATRIX[3] Ms^2},
{me2, UNITMATRIX[3] Ms^2},
{\[Mu], Ms},
{B[\[Mu]], Sqr[Ms]/(TanBeta + 1./TanBeta)},
{T[Yu], MODELPARAMETER[Yu] Ms/TanBeta},
{T[Yd], MODELPARAMETER[Yd] Ms TanBeta},
{T[Ye], MODELPARAMETER[Ye] Ms TanBeta},
{T[Yu][3,3], (1/TanBeta + Xtt) Ms }
};
LowScale = LowEnergyConstant[MZ];
LowScaleFirstGuess = LowEnergyConstant[MZ];
LowScaleInput = {
{Yu, Automatic},
{Yd, Automatic},
{Ye, Automatic},
{vu, Sqrt[4 MZDRbar^2 / (g2^2 + GUTNormalization[g1]^2 g1^2) - vd^2]}
};
UseHiggs2LoopMSSM = True;
EffectiveMu = \[Mu];
PotentialLSPParticles = { Chi, Sv, Su, Sd, Se, Cha, Glu };
DefaultPoleMassPrecision = HighPrecision;
HighPoleMassPrecision = {hh, Ah, Hpm};
MediumPoleMassPrecision = {};
LowPoleMassPrecision = {};
ExtraSLHAOutputBlocks = {
{FlexibleSUSYOutput,
{{0, Hold[HighScale]},
{1, Hold[SUSYScale]},
{2, Hold[LowScale]} } },
{ALPHA, {{ArcSin[Pole[ZH[2,2]]]}}},
{HMIX , {{1, \[Mu]},
{2, vu / vd},
{3, Sqrt[vu^2 + vd^2]},
{4, M[Ah[2]]^2},
{101, B[\[Mu]]},
{102, vd},
{103, vu} } },
{Au, {{1, 1, T[Yu][1,1] / Yu[1,1]},
{2, 2, T[Yu][2,2] / Yu[2,2]},
{3, 3, T[Yu][3,3] / Yu[3,3]} } },
{Ad, {{1, 1, T[Yd][1,1] / Yd[1,1]},
{2, 2, T[Yd][2,2] / Yd[2,2]},
{3, 3, T[Yd][3,3] / Yd[3,3]} } },
{Ae, {{1, 1, T[Ye][1,1] / Ye[1,1]},
{2, 2, T[Ye][2,2] / Ye[2,2]},
{3, 3, T[Ye][3,3] / Ye[3,3]} } },
{MSOFT, {{1, MassB},
{2, MassWB},
{3, MassG},
{21, mHd2},
{22, mHu2},
{31, SignedAbsSqrt[ml2[1,1]]},
{32, SignedAbsSqrt[ml2[2,2]]},
{33, SignedAbsSqrt[ml2[3,3]]},
{34, SignedAbsSqrt[me2[1,1]]},
{35, SignedAbsSqrt[me2[2,2]]},
{36, SignedAbsSqrt[me2[3,3]]},
{41, SignedAbsSqrt[mq2[1,1]]},
{42, SignedAbsSqrt[mq2[2,2]]},
{43, SignedAbsSqrt[mq2[3,3]]},
{44, SignedAbsSqrt[mu2[1,1]]},
{45, SignedAbsSqrt[mu2[2,2]]},
{46, SignedAbsSqrt[mu2[3,3]]},
{47, SignedAbsSqrt[md2[1,1]]},
{48, SignedAbsSqrt[md2[2,2]]},
{49, SignedAbsSqrt[md2[3,3]]} } }
};