From b0cd3f24c912fe0bfdc63b04d6a9a6be351cbffd Mon Sep 17 00:00:00 2001 From: Alexander Voigt Date: Fri, 2 Sep 2016 23:28:16 +0200 Subject: [PATCH] adding Peter's common MSUSY E6SSMtower model file --- model_files/E6SSMtower/FlexibleSUSY.m.in | 99 +++++++++++++++++++ .../E6SSMtower/LesHouches.in.E6SSMtower | 48 +++++++++ 2 files changed, 147 insertions(+) create mode 100644 model_files/E6SSMtower/FlexibleSUSY.m.in create mode 100644 model_files/E6SSMtower/LesHouches.in.E6SSMtower diff --git a/model_files/E6SSMtower/FlexibleSUSY.m.in b/model_files/E6SSMtower/FlexibleSUSY.m.in new file mode 100644 index 000000000..730b74e42 --- /dev/null +++ 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 new file mode 100644 index 000000000..69064cd76 --- /dev/null +++ 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