adding Peter's common MSUSY E6SSMtower model file

FlexibleSUSY · Sep 2, 2016 · b0cd3f2 · b0cd3f2
1 parent e2a1098
commit b0cd3f2
Show file tree

Hide file tree

Showing 2 changed files with 147 additions and 0 deletions.
diff --git a/model_files/E6SSMtower/FlexibleSUSY.m.in b/model_files/E6SSMtower/FlexibleSUSY.m.in
@@ -0,0 +1,99 @@
+FSModelName = "@CLASSNAME@";
+FSEigenstates = SARAH`EWSB;
+FSDefaultSARAHModel = E6SSM;
+OnlyLowEnergyFlexibleSUSY = True;
+FlexibleEFTHiggs = True;
+
+EXTPAR = {
+    (* {0, MSUSY}, *)
+    (* {1, M1Input}, *)
+    (* {2, M2Input}, *)
+    (* {3, M3Input}, *)
+    (* {4, MuInput}, *)
+    (* {5, mAInput}, *)
+    (* {25, TanBeta}, *)
+    {61, LambdaInput},
+    {62, KappaInput},
+    {65, vsInput},
+    {66, Lambda12Input},
+    {70, gNInput}
+};
+
+FSExtraInputParameters = {
+    {TanBeta, TanBeta, {1}},
+    {MSUSY, MSUSY, {1}},
+    {Xtt, Xtt, {1}}
+};
+
+EWSBOutputParameters = { mHd2, mHu2, ms2 };
+
+SUSYScale = MSUSY;
+
+SUSYScaleFirstGuess = MSUSY;
+
+SUSYScaleInput = {
+    {gN, gNInput},
+    {MassB, MSUSY},
+    {MassWB, MSUSY},
+    {MassG, MSUSY},
+    {MassBp, MSUSY},
+    {me2, MSUSY^2 UNITMATRIX[3]},
+    {mq2, MSUSY^2 UNITMATRIX[3]},
+    {mu2, MSUSY^2 UNITMATRIX[3]},
+    {md2, MSUSY^2 UNITMATRIX[3]},
+    {ml2, MSUSY^2 UNITMATRIX[3]},
+    {mDx2, MSUSY^2 UNITMATRIX[3]},
+    {mDxbar2, MSUSY^2 UNITMATRIX[3]},
+    {mH1I2, MSUSY^2 UNITMATRIX[2]},
+    {mH2I2, MSUSY^2 UNITMATRIX[2]},
+    {msI2, MSUSY^2 UNITMATRIX[2]},
+    {mDx2[1,1], MSUSY^2},
+    {mDxbar2[1,1], MSUSY^2},
+    {mH1I2[1,1], MSUSY^2},
+    {mH2I2[1,1], MSUSY^2},
+    {msI2[1,1], MSUSY^2},
+    {mDx2[2,2], MSUSY^2},
+    {mDxbar2[2,2], MSUSY^2},
+    {mH1I2[2,2], MSUSY^2},
+    {mH2I2[2,2], MSUSY^2},
+    {msI2[2,2], MSUSY^2},
+    {mDx2[3,3], MSUSY^2},
+    {mDxbar2[3,3], MSUSY^2},
+    {mHp2, MSUSY^2},
+    {mHpbar2, MSUSY^2},
+    {B[\[Mu]Pr], MSUSY^2},
+    {\[Mu]Pr, MSUSY},
+    {\[Kappa], UNITMATRIX[3] KappaInput},
+    {\[Lambda]12, UNITMATRIX[2] Lambda12Input},
+    {\[Lambda], LambdaInput},
+    {T[\[Kappa]], UNITMATRIX[3] MSUSY KappaInput},
+    {T[\[Lambda]12], UNITMATRIX[2] MSUSY Lambda12Input},
+    {vs, vsInput},
+    {T[\[Lambda]], MSUSY^2 Sqrt[2] / (vs (TanBeta + 1/TanBeta)) },
+    {T[Yu], Yu \[Lambda] vs / ( Sqrt[2] TanBeta) },
+    {T[Yd], TanBeta Yd \[Lambda] vs / Sqrt[2] },
+    {T[Ye], TanBeta Ye \[Lambda] vs / Sqrt[2] },
+    {T[Yu][3,3], (\[Lambda] vs / (Sqrt[2] TanBeta) + Xtt) Yu[3,3] }
+};
+
+InitialGuessAtSUSYScale = SUSYScaleInput;
+
+SUSYScaleMatching = {
+    {vu, VEV Sin[ArcTan[TanBeta]]},
+    {vd, VEV Cos[ArcTan[TanBeta]]}
+};
+
+LowScale = LowEnergyConstant[MZ];
+
+LowScaleFirstGuess = LowEnergyConstant[MZ];
+
+UseHiggs2LoopNMSSM = True;
+EffectiveMu =  \[Lambda] vs / Sqrt[2];
+EffectiveMASqr = T[\[Lambda]] vs / Sqrt[2] (vu^2 + vd^2) / (vu vd);
+
+PotentialLSPParticles = { Chi, Sv, Su, Sd, Se, Cha, Glu };
+
+DefaultPoleMassPrecision = HighPrecision;
+HighPoleMassPrecision    = {hh, Ah, Hpm};
+MediumPoleMassPrecision  = {};
+LowPoleMassPrecision     = {};
diff --git a/model_files/E6SSMtower/LesHouches.in.E6SSMtower b/model_files/E6SSMtower/LesHouches.in.E6SSMtower
@@ -0,0 +1,48 @@
+Block MODSEL                 # Select model
+    6   0                    # flavour violation
+Block FlexibleSUSY
+    0   1.000000000e-05      # precision goal
+    1   0                    # max. iterations (0 = automatic)
+    2   0                    # algorithm (0 = two_scale, 1 = lattice)
+    3   1                    # calculate SM pole masses
+    4   2                    # pole mass loop order
+    5   2                    # EWSB loop order
+    6   3                    # beta-functions loop order
+    7   2                    # threshold corrections loop order
+    8   0                    # Higgs 2-loop corrections O(alpha_t alpha_s)
+    9   0                    # Higgs 2-loop corrections O(alpha_b alpha_s)
+   10   0                    # Higgs 2-loop corrections O((alpha_t + alpha_b)^2)
+   11   0                    # Higgs 2-loop corrections O(alpha_tau^2)
+   12   1                    # force output
+   13   1                    # Top quark 2-loop corrections QCD
+   14   1.000000000e-11      # beta-function zero threshold
+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 EXTPAR
+    61	 1.0000000000000001E-01	 # lambda(MS) 
+    62	 1.0000000000000001E-01	 # kappa(MS)
+    65   4.2426406871192862E+05       # vs(MS)
+    66	 1.0000000000000001E-01	 # lam12(MS)
+    70   8.9442719099991588E-02       # gN(MS) 
+Block MSUSY
+     3.0000000000000007E+04             # SUSY scale
+Block TanBeta
+     5.0000000000000000E+00           # tan(Beta) at the SUSY scale
+Block Xtt
+     0.0000000000000000E+00             # Xt / Ms