Skip to content

Commit

Permalink
adding SM with extra singlet (for testing FlexibleBSM)
Browse files Browse the repository at this point in the history
  • Loading branch information
Alexander Voigt authored and Alexander Voigt committed Nov 30, 2014
1 parent b51e887 commit eec603c
Show file tree
Hide file tree
Showing 5 changed files with 319 additions and 0 deletions.
49 changes: 49 additions & 0 deletions model_files/NSM/FlexibleSUSY.m.in
View file
@@ -0,0 +1,49 @@

FSModelName = "@CLASSNAME@";
FSEigenstates = SARAH`EWSB;

(* input parameters *)

MINPAR = {
{1, LambdaInput1},
{2, LambdaInput2},
{3, LambdaInput3},
{4, LambdaInput4},
{5, LambdaInput5}
};

EWSBOutputParameters = { mH2, mS2 };

SUSYScale = vS;

SUSYScaleFirstGuess = 1000;

SUSYScaleInput = {
{\[Lambda]1, LambdaInput1},
{\[Lambda]2, LambdaInput2},
{\[Lambda]3, LambdaInput3},
{\[Lambda]4, LambdaInput4},
{\[Lambda]5, LambdaInput5}
};

LowScale = SM[MZ];

LowScaleFirstGuess = SM[MZ];

LowScaleInput = {
{Yu, Automatic},
{Yd, Automatic},
{Ye, Automatic},
{vH, 2 MZDRbar / Sqrt[GUTNormalization[g1]^2 g1^2 + g2^2]},
{vS, 2 MZDRbar / Sqrt[GUTNormalization[g1]^2 g1^2 + g2^2]}
};

InitialGuessAtLowScale = {
{Yu, Automatic},
{Yd, Automatic},
{Ye, Automatic},
{vH, SM[vev]},
{vS, SM[vev]}
};

OnlyLowEnergyFlexibleSUSY = True;
38 changes: 38 additions & 0 deletions model_files/NSM/LesHouches.in.NSM
View file
@@ -0,0 +1,38 @@
Block MODSEL # Select model
# 12 1000 # 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 1 # pole mass loop order
5 1 # EWSB loop order
6 2 # beta-functions loop order
7 1 # 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^2 + alpha_t alpha_b + alpha_b^2)
11 0 # Higgs 2-loop corrections O(alpha_tau^2)
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)
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 0.1 # LambdaInput1
2 0.1 # LambdaInput2
3 0.1 # LambdaInput3
4 0 # LambdaInput4
5 0 # LambdaInput5
98 changes: 98 additions & 0 deletions sarah/NSM/NSM.m
View file
@@ -0,0 +1,98 @@
Off[General::spell]

Model`Name = "NSM";
Model`NameLaTeX ="Next-to Standard Model";
Model`Authors = "F.Staub, A.Voigt";
Model`Date = "2014-11-27";

(*-------------------------------------------*)
(* Particle Content*)
(*-------------------------------------------*)

(* Gauge Groups *)

Gauge[[1]] = {B, U[1], hypercharge, g1, False};
Gauge[[2]] = {WB, SU[2], left, g2, True};
Gauge[[3]] = {G, SU[3], color, g3, False};

(* Matter Fields *)

FermionFields[[1]] = {q, 3, {uL, dL}, 1/6, 2, 3};
FermionFields[[2]] = {l, 3, {vL, eL}, -1/2, 2, 1};
FermionFields[[3]] = {d, 3, conj[dR], 1/3, 1, -3};
FermionFields[[4]] = {u, 3, conj[uR], -2/3, 1, -3};
FermionFields[[5]] = {e, 3, conj[eR], 1, 1, 1};

ScalarFields[[1]] = {H, 1, {Hp, H0}, 1/2, 2, 1};
ScalarFields[[2]] = {S, 1, S0 , 0, 1, 1};

RealScalars = {S};

(*----------------------------------------------*)
(* DEFINITION *)
(*----------------------------------------------*)

NameOfStates = {GaugeES, EWSB};

(* ----- Before EWSB ----- *)

DEFINITION[GaugeES][LagrangianInput] = {
{LagHC , {AddHC->True }},
{LagNoHC, {AddHC->False}}
};

LagZ2 = \
- mH2 conj[H].H \
- mS2 S.S \
+ \[Lambda]1 conj[H].H.conj[H].H \
+ \[Lambda]2 S.S.S.S \
+ \[Lambda]3 S.S.conj[H].H;

LagZ2v = \
+ \[Lambda]4 S.conj[H].H \
+ \[Lambda]5 S.S.S;

LagNoHC = LagZ2 + LagZ2v;

LagHC = Yd conj[H].d.q + Ye conj[H].e.l + Yu H.u.q;

(* Gauge Sector *)

DEFINITION[EWSB][GaugeSector] = {
{{VB,VWB[3]}, {VP,VZ}, ZZ},
{{VWB[1],VWB[2]}, {VWp,conj[VWp]}, ZW}
};

(* ----- VEVs ---- *)

DEFINITION[EWSB][VEVs] = {
{H0, {vH, 1/Sqrt[2]}, {Ah, \[ImaginaryI]/Sqrt[2]}, {phiH, 1/Sqrt[2]}},
{S0, {vS, 1}, {0, 0}, {phiS, 1}}
};

DEFINITION[EWSB][MatterSector] = {
{{{dL}, {conj[dR]}}, {{DL,Vd}, {DR,Ud}}},
{{{uL}, {conj[uR]}}, {{UL,Vu}, {UR,Uu}}},
{{{eL}, {conj[eR]}}, {{EL,Ve}, {ER,Ue}}},
{{phiH, phiS} , {hh, ZH}}
};

(*------------------------------------------------------*)
(* Dirac-Spinors *)
(*------------------------------------------------------*)

DEFINITION[EWSB][DiracSpinors] = {
Fd -> {DL, conj[DR]},
Fe -> {EL, conj[ER]},
Fu -> {UL, conj[UR]},
Fv -> {vL, 0}
};

DEFINITION[EWSB][GaugeES] = {
Fd1 -> {FdL, 0},
Fd2 -> {0, FdR},
Fu1 -> {Fu1, 0},
Fu2 -> {0, Fu2},
Fe1 -> {Fe1, 0},
Fe2 -> {0, Fe2}
};
36 changes: 36 additions & 0 deletions sarah/NSM/parameters.m
View file
@@ -0,0 +1,36 @@
ParameterDefinitions = {
{g1, { Description -> "Hypercharge-Coupling"}},
{g2, { Description -> "Left-Coupling"}},
{g3, { Description -> "Strong-Coupling"}},
{AlphaS, { Description -> "Alpha Strong"}},
{e, { Description -> "electric charge"}},
{Gf, { Description -> "Fermi's constant"}},
{aEWinv, { Description -> "inverse weak coupling constant at mZ"}},
{Yu, { Description -> "Up-Yukawa-Coupling"}},
{Yd, { Description -> "Down-Yukawa-Coupling"}},
{Ye, { Description -> "Lepton-Yukawa-Coupling"}},
{mH2, { Description -> "SM Mu Parameter"}},
{mS2, { Description -> "Singlet mass parameter"}},
{\[Lambda]1,{ Description -> "SM Higgs Selfcouplings"}},
{\[Lambda]2,{ Description -> "Singlet 4-selfcoupling"}},
{\[Lambda]3,{ Description -> "Singlet 2S-2H coupling"}},
{\[Lambda]4,{ Description -> "Singlet S-2H coupling"}},
{\[Lambda]5,{ Description -> "Singlet 3-selfcoupling"}},
{vH, { Description -> "EW-VEV"}},
{vS, { Description -> "Singlet-VEV",
DependenceSPheno -> None }},
{mH2, { Description -> "SM Higgs Mass Parameter"}},
{ThetaW, { Description -> "Weinberg-Angle",
DependenceNum -> ArcSin[Sqrt[1 - Mass[VWp]^2/Mass[VZ]^2]] }},
{ZZ, { Description -> "Photon-Z Mixing Matrix"}},
{ZH, { Description -> "Higgs Mixing Matrix"}},
{ZW, { Description -> "W Mixing Matrix",
Dependence -> 1/Sqrt[2] {{1, 1},
{\[ImaginaryI],-\[ImaginaryI]}} }},
{Vu, { Description ->"Left-Up-Mixing-Matrix"}},
{Vd, { Description ->"Left-Down-Mixing-Matrix"}},
{Uu, { Description ->"Right-Up-Mixing-Matrix"}},
{Ud, { Description ->"Right-Down-Mixing-Matrix"}},
{Ve, { Description ->"Left-Lepton-Mixing-Matrix"}},
{Ue, { Description ->"Right-Lepton-Mixing-Matrix"}}
};
98 changes: 98 additions & 0 deletions sarah/NSM/particles.m
View file
@@ -0,0 +1,98 @@
ParticleDefinitions[GaugeES] = {
{H0, { PDG -> {0},
Width -> 0,
Mass -> Automatic,
FeynArtsNr -> 1,
LaTeX -> "H^0",
OutputName -> "H0" }},

{Hp, { PDG -> {0},
Width -> 0,
Mass -> Automatic,
FeynArtsNr -> 2,
LaTeX -> "H^+",
OutputName -> "Hp" }},

{S0, { PDG -> {0},
Width -> 0,
Mass -> Automatic,
FeynArtsNr -> 3,
LaTeX -> "S^0",
OutputName -> "S0" }},

{VB, { Description -> "B-Boson"}},
{VG, { Description -> "Gluon"}},
{VWB, { Description -> "W-Bosons"}},
{gB, { Description -> "B-Boson Ghost"}},
{gG, { Description -> "Gluon Ghost" }},
{gWB, { Description -> "W-Boson Ghost"}}

};

ParticleDefinitions[EWSB] = {
{hh , { Description -> "Higgs",
PDG -> {25,35},
PDG.IX -> {101000001, 101000002} }},

{Ah , { Description -> "Pseudo-Scalar Higgs",
PDG -> {0},
PDG.IX ->{0},
Mass -> {0},
Width -> {0} }},

{Hp, { Description -> "Charged Higgs",
PDG -> {0},
PDG.IX ->{0},
Width -> {0},
Mass -> {0},
LaTeX -> {"H^+","H^-"},
OutputName -> {"Hp","Hm"}
}},

{VP, { Description -> "Photon"}},
{VZ, { Description -> "Z-Boson",
Goldstone -> Ah }},
{VG, { Description -> "Gluon" }},
{VWp, { Description -> "W+ - Boson",
Goldstone -> Hp }},
{gP, { Description -> "Photon Ghost"}},
{gWp, { Description -> "Positive W+ - Boson Ghost"}},
{gWpC, { Description -> "Negative W+ - Boson Ghost" }},
{gZ, { Description -> "Z-Boson Ghost" }},
{gG, { Description -> "Gluon Ghost" }},

{Fd, { Description -> "Down-Quarks"}},
{Fu, { Description -> "Up-Quarks"}},
{Fe, { Description -> "Leptons" }},
{Fv, { Description -> "Neutrinos" }}
};

WeylFermionAndIndermediate = {
{H, { PDG -> {0},
Width -> 0,
Mass -> Automatic,
LaTeX -> "H",
OutputName -> "" }},
{S, { PDG -> {0},
Width -> 0,
Mass -> Automatic,
LaTeX -> "S",
OutputName -> "" }},

{dR, {LaTeX -> "d_R" }},
{eR, {LaTeX -> "e_R" }},
{lep, {LaTeX -> "l" }},
{uR, {LaTeX -> "u_R" }},
{q, {LaTeX -> "q" }},
{eL, {LaTeX -> "e_L" }},
{dL, {LaTeX -> "d_L" }},
{uL, {LaTeX -> "u_L" }},
{vL, {LaTeX -> "\\nu_L" }},

{DR, {LaTeX -> "D_R" }},
{ER, {LaTeX -> "E_R" }},
{UR, {LaTeX -> "U_R" }},
{EL, {LaTeX -> "E_L" }},
{DL, {LaTeX -> "D_L" }},
{UL, {LaTeX -> "U_L" }}
};

0 comments on commit eec603c

Please sign in to comment.