adding first version of MRSSM-tower model file

FlexibleSUSY · Apr 6, 2016 · 8ebd6f4 · 8ebd6f4
1 parent 416f573
commit 8ebd6f4
Show file tree

Hide file tree

Showing 2 changed files with 168 additions and 0 deletions.
diff --git a/model_files/MRSSMtower/FlexibleSUSY.m.in b/model_files/MRSSMtower/FlexibleSUSY.m.in
@@ -0,0 +1,103 @@
+
+FSModelName = "MRSSMtower";
+FSEigenstates = SARAH`EWSB;
+FSDefaultSARAHModel = MRSSM;
+AutomaticInputAtMSUSY = False;
+
+SMTower = True;
+
+OnlyLowEnergyFlexibleSUSY = True;
+
+(* CMSSM input parameters *)
+
+MINPAR = {
+    {3, TanBeta}
+};
+
+EXTPAR = {
+    {0, Ms}
+};
+
+EWSBOutputParameters = { mHd2, mHu2, mS2, mT2 };
+
+SUSYScale = Ms;
+
+SUSYScaleFirstGuess = Ms;
+
+SUSYScaleInput = {
+    {mq2, UNITMATRIX[3] Sqr[Ms]},
+    {mu2, UNITMATRIX[3] Sqr[Ms]},
+    {md2, UNITMATRIX[3] Sqr[Ms]},
+    {ml2, UNITMATRIX[3] Sqr[Ms]},
+    {me2, UNITMATRIX[3] Sqr[Ms]},
+    {\[Mu], Ms},
+    {B[\[Mu]], Sqr[Ms]/(TanBeta + 1/TanBeta)},
+    {LamSD , SARAH`LHInput[LamSD]},
+    {LamSU , SARAH`LHInput[LamSU]},
+    {LamTD , SARAH`LHInput[LamTD]},
+    {LamTU , SARAH`LHInput[LamTU]},
+    {MDBS  , SARAH`LHInput[MDBS]},
+    {MDGoc , SARAH`LHInput[MDGoc]},
+    {MDWBT , SARAH`LHInput[MDWBT]},
+    {moc2  , Sqr[Ms]},
+    {mRd2  , SARAH`LHInput[mRd2]},
+    {mRu2  , SARAH`LHInput[mRu2]},
+    {MuD   , SARAH`LHInput[MuD]}, 
+    {MuU   , SARAH`LHInput[MuU]},
+    {vS    , SARAH`LHInput[vS]},
+    {vT    , SARAH`LHInput[vT]},
+    {B[MuD], SARAH`LHInput[B[MuD]]},
+    {B[MuU], SARAH`LHInput[B[MuU]]}
+};
+
+InitialGuessAtSUSYScale = {
+    {mq2, UNITMATRIX[3] Sqr[Ms]},
+    {mu2, UNITMATRIX[3] Sqr[Ms]},
+    {md2, UNITMATRIX[3] Sqr[Ms]},
+    {ml2, UNITMATRIX[3] Sqr[Ms]},
+    {me2, UNITMATRIX[3] Sqr[Ms]},
+    {\[Mu], Ms},
+    {B[\[Mu]], Sqr[Ms]/(TanBeta + 1/TanBeta)}
+};
+
+(* InitialGuessAtLowScale is ignored if SMTower is true *)
+
+InitialGuessAtLowScale = {
+    {vd, LowEnergyConstant[vev] Cos[ArcTan[TanBeta]]},
+    {vu, LowEnergyConstant[vev] Sin[ArcTan[TanBeta]]},
+    {Yu, Automatic},
+    {Yd, Automatic},
+    {Ye, Automatic}
+};
+
+(* vev is a keyword in SUSYScaleUserMatching, describing SM VEV matched into full theory *)
+
+SUSYScaleUserMatching = {
+    {vu, vev Sin[ArcTan[INPUTPARAMETER[TanBeta]]]},
+    {vd, vev Cos[ArcTan[INPUTPARAMETER[TanBeta]]]}
+};
+
+LowScale = LowEnergyConstant[MZ];
+
+LowScaleFirstGuess = LowEnergyConstant[MZ];
+
+DefaultPoleMassPrecision = HighPrecision;
+HighPoleMassPrecision    = {hh, Ah, Hpm};
+MediumPoleMassPrecision  = {};
+LowPoleMassPrecision     = {};
+
+FSWeakMixingAngleOptions = FSSetOption[
+    FSWeakMixingAngleOptions,
+    FSWeakMixingAngleInput -> FSMassW
+];
+FSWeakMixingAngleOptions = FSSetOption[
+    FSWeakMixingAngleOptions,
+    FSWeakMixingAngleExpr  -> ArcSin[Sqrt[1 - (Mass[VWm]^2 - g2^2*vT^2)/Mass[VZ]^2]]
+];
+
+SARAH`ParameterDefinitions = Append[
+    Cases[SARAH`ParameterDefinitions, {p_ /; p =!= SARAH`ThetaW, ___}],
+    {SARAH`ThetaW,
+     { Description -> "Weinberg-Angle",
+       DependenceNum -> ArcSin[Sqrt[1 - (Mass[VWm]^2 - g2^2*vT^2)/Mass[VZ]^2]]}}
+];
diff --git a/model_files/MRSSMtower/LesHouches.in.MRSSMtower b/model_files/MRSSMtower/LesHouches.in.MRSSMtower
@@ -0,0 +1,65 @@
+Block MODSEL     #
+    1 0          #  1/0: High/low scale input
+    2 1          # Boundary Condition
+    6 1          # Generation Mixing
+#   12 1000       # output scale
+Block FlexibleSUSY
+    0   1.000000000e-04      # precision goal
+    1   0                    # max. iterations (0 = automatic)
+    2   0                    # algorithm (0 = two_scale, 1 = lattice)
+    3   0                    # calculate SM pole masses
+    4   2                    # pole mass loop order
+    5   2                    # EWSB loop order
+    6   2                    # 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)
+   10   1                    # Higgs 2-loop corrections O((alpha_t + alpha_b)^2)
+   11   1                    # Higgs 2-loop corrections O(alpha_tau^2)
+   12   0                    # force output
+   13   1                    # Top quark 2-loop corrections QCD
+   14   1.000000000e-11      # beta-function zero threshold
+   15   0                    # calculate observables (a_muon, ...)
+   16   0                    # force positive majorana masses
+Block SMINPUTS               # Standard Model inputs
+    1   1.279440000e+02      # alpha^(-1) SM MSbar(MZ)
+    2   1.166380000e-05      # G_Fermi
+    3   1.184000000e-01      # alpha_s(MZ) SM MSbar
+    4   9.118760000e+01      # MZ(pole)
+    5   4.180000000e+00      # mb(mb) SM MSbar
+    6   1.733400000e+02      # mtop(pole)
+    7   1.777000000e+00      # mtau(pole)
+    8   0.000000000e+00      # mnu3(pole)
+    9   80.384               # MW pole
+   11   5.109989020e-04      # melectron(pole)
+   12   0.000000000e+00      # mnu1(pole)
+   13   1.056583570e-01      # mmuon(pole)
+   14   0.000000000e+00      # mnu2(pole)
+   21   4.750000000e-03      # md(2 GeV) MS-bar
+   22   2.400000000e-03      # mu(2 GeV) MS-bar
+   23   1.040000000e-01      # ms(2 GeV) MS-bar
+   24   1.270000000e+00      # mc(mc) MS-bar
+Block MINPAR
+    3    5                # TanBeta
+Block EXTPAR
+    0    1e3              # Ms
+Block NMSSMRUNIN
+    5    1                # vS
+Block MSOFTIN
+    300  10               # MDB
+    302  2000             # MDO
+    301  100              # MDW
+    50   5e4              # mRd2
+    51   5e4              # mRu2
+Block HMIXIN
+    301 -0.1              # LSD
+    302 -0.1              # LSU
+    303 -0.1              # LTD
+    304 -0.1              # LTU
+    201  250              # MuD
+    202  250              # MuU
+    310  1                # vT
+    1    0                # Mu
+    203  0                # BmuD
+    204  0                # BmuU
+    101  1e4              # Bmu