Skip to content

Commit

Permalink
Add VCMSSM model file for testing
Browse files Browse the repository at this point in the history
  • Loading branch information
Dylan Harries committed Dec 3, 2016
1 parent 676a662 commit d492b3b
Show file tree
Hide file tree
Showing 2 changed files with 208 additions and 0 deletions.
153 changes: 153 additions & 0 deletions model_files/VCMSSM/FlexibleSUSY.m.in
View file
Original file line number Diff line number Diff line change
@@ -0,0 +1,153 @@
FSModelName = "@CLASSNAME@";
FSEigenstates = SARAH`EWSB;
FSDefaultSARAHModel = MSSM;

(* VCMSSM input parameters *)

MINPAR = {
{1, m0},
{2, m12},
{4, Sign[\[Mu]]},
{5, Azero}
};

EXTPAR = {
{1, TBGuess}
};

FSExtraParameters = {
{vMSSM, {1}}
};


EWSBSubstitutions = {
{vd, vMSSM Cos[ArcTan[TanBeta]]},
{vu, vMSSM Sin[ArcTan[TanBeta]]}
};

EWSBOutputParameters = { TanBeta, \[Mu] };

EWSBInitialGuess = {
{TanBeta, vu / vd}
};

SUSYScale = Sqrt[Product[M[Su[i]]^(Abs[ZU[i,3]]^2 + Abs[ZU[i,6]]^2), {i,6}]];

SUSYScaleFirstGuess = Sqrt[m0^2 + 4 m12^2];

SUSYScaleInput = {
{vMSSM, Sqrt[vu^2 + vd^2]}
};

HighScale = g1 == g2;

HighScaleFirstGuess = 2.0 10^16;

HighScaleInput = {
{T[Ye], Azero*Ye},
{T[Yd], Azero*Yd},
{T[Yu], Azero*Yu},
{mHd2, m0^2},
{mHu2, m0^2},
{mq2, UNITMATRIX[3] m0^2},
{ml2, UNITMATRIX[3] m0^2},
{md2, UNITMATRIX[3] m0^2},
{mu2, UNITMATRIX[3] m0^2},
{me2, UNITMATRIX[3] m0^2},
{MassB, m12},
{MassWB,m12},
{MassG, m12},
{B[\[Mu]], (Azero + m0) \[Mu]}
};

LowScale = LowEnergyConstant[MZ];

LowScaleFirstGuess = LowEnergyConstant[MZ];

LowScaleInput = {
{Yu, Automatic},
{Yd, Automatic},
{Ye, Automatic},
{vd, 2 MZDRbar / Sqrt[GUTNormalization[g1]^2 g1^2 + g2^2] Cos[ArcTan[vu / vd]]},
{vu, 2 MZDRbar / Sqrt[GUTNormalization[g1]^2 g1^2 + g2^2] Sin[ArcTan[vu / vd]]}
};

InitialGuessAtLowScale = {
{vd, LowEnergyConstant[vev] Cos[ArcTan[TBGuess]]},
{vu, LowEnergyConstant[vev] Sin[ArcTan[TBGuess]]},
{Yu, Automatic},
{Yd, Automatic},
{Ye, Automatic},
{vMSSM, LowEnergyConstant[vev]},
{\[Mu], Sign[\[Mu]] LowEnergyConstant[MZ]}
};

InitialGuessAtHighScale = {};

UseHiggs2LoopMSSM = True;
EffectiveMu = \[Mu];
UseMSSM3LoopRGEs = True;

(* low-energy input to determine weak mixing angle
* Possible values: FSFermiConstant (default), FSMassW
*)
FSWeakMixingAngleOptions = FSSetOption[FSWeakMixingAngleOptions, FSWeakMixingAngleInput -> FSFermiConstant];

PotentialLSPParticles = { Chi, Sv, Su, Sd, Se, Cha, Glu };

DefaultPoleMassPrecision = MediumPrecision;
HighPoleMassPrecision = {hh, Ah, Hpm};
MediumPoleMassPrecision = {};
LowPoleMassPrecision = {};

ExtraSLHAOutputBlocks = {
{FlexibleSUSYOutput, NoScale,
{{0, Hold[HighScale]},
{1, Hold[SUSYScale]},
{2, Hold[LowScale]} } },
{FlexibleSUSYLowEnergy,
{{21, FlexibleSUSYObservable`aMuon} } },
{EFFHIGGSCOUPLINGS, NoScale,
{{1, FlexibleSUSYObservable`CpHiggsPhotonPhoton},
{2, FlexibleSUSYObservable`CpHiggsGluonGluon},
{3, FlexibleSUSYObservable`CpPseudoScalarPhotonPhoton},
{4, FlexibleSUSYObservable`CpPseudoScalarGluonGluon} } },
{ALPHA, NoScale,
{{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]]} } }
};
55 changes: 55 additions & 0 deletions model_files/VCMSSM/LesHouches.in.VCMSSM
View file
Original file line number Diff line number Diff line change
@@ -0,0 +1,55 @@
Block MODSEL # Select model
# 12 1000 # DRbar parameter output scale (GeV)
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 3 # 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 pole mass QCD corrections (0 = 1L, 1 = 2L, 2 = 3L)
14 1.000000000e-11 # beta-function zero threshold
15 0 # calculate observables (a_muon, ...)
16 0 # force positive majorana masses
17 0 # pole mass renormalization scale (0 = SUSY scale)
18 0 # pole mass renormalization scale in the EFT (0 = min(SUSY scale, Mt))
19 0 # EFT matching scale (0 = SUSY scale)
20 2 # EFT loop order for upwards matching
21 1 # EFT loop order for downwards matching
22 0 # EFT index of SM-like Higgs in the BSM model
23 1 # calculate BSM pole masses
Block FlexibleSUSYInput
0 0.00729735 # alpha_em(0)
1 125.09 # Mh pole
Block SMINPUTS # Standard Model inputs
1 1.279340000e+02 # alpha^(-1) SM MSbar(MZ)
2 1.166370000e-05 # G_Fermi
3 1.176000000e-01 # alpha_s(MZ) SM MSbar
4 9.118760000e+01 # MZ(pole)
5 4.200000000e+00 # mb(mb) SM MSbar
6 1.733000000e+02 # mtop(pole)
7 1.777000000e+00 # mtau(pole)
8 0.000000000e+00 # mnu3(pole)
9 80.404 # 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 # Input parameters
1 125. # m0
2 500. # m12
4 1 # SignMu
5 0 # Azero
Block EXTPAR
1 10. # TBGuess

0 comments on commit d492b3b

Please sign in to comment.