Step 1

Load FeynRules and the model

Quit[]

$FeynRulesPath =SetDirectory["C:\\Users\\Céline\\Documents\\UCL\\Doctorat\\FeynRules-current\\FeynRules-current-4"];
<< FeynRules`

LoadModel["C:/users/Céline/Documents/UCL/Doctorat\\FeynRules-current\\FeynRules-current-4\Models\SM\SM.fr","C:/users/Céline/Documents/UCL/Doctorat/Generic//tutorial.fr"]

LNP

Masses

FullSimplify[GetMassTerms[LNP+LSM,FlavorExpand->SU2W,SelectParticles->{{Z,Z},{Z,Zprime},{Zprime,Zprime}}]]

FullSimplify[GetMassTerms[LNP+LSM,FlavorExpand->SU2W,SelectParticles->{{Z,Z},{Z,Zprime},{Zprime,Zprime}}]/.{ee->gw*sw}/.{sw->g1/gw cw}/.{cw->gw/Sqrt[gw^2+g1^2]}/.{v->2MZ /Sqrt[gw^2+g1^2]}]

CheckDiagonalMassTerms[LNP+LSM,FlavorExpand->SU2W,SelectParticles->{{Z,Z},{Z,Zprime},{Zprime,Zprime},{Z,A},{Zprime,A},{A,A}}]

CheckMassSpectrum[LNP+LSM,FlavorExpand->SU2W,SelectParticles->{{Z,Z},{Z,Zprime},{Zprime,Zprime},{Z,A},{Zprime,A},{A,A}}]

Mass diagonalization

M2={{MZ^2,2grho MZ^2},{2 grho MZ^2,mrho^2}}

Eigenvalues[M2]

EV=Eigenvectors[M2]

Step 2

Load FeynRules and the model

Quit[]

$FeynRulesPath =SetDirectory["C:\\Users\\Céline\\Documents\\UCL\\Doctorat\\FeynRules-current\\FeynRules-current-4"];
<< FeynRules`

LoadModel["C:/users/Céline/Documents/UCL/Doctorat/Generic//SMUnphysZ.fr","C:/users/Céline/Documents/UCL/Doctorat/Generic//tutorial2.fr"]

LNP

Checks

Kinetic

CheckDiagonalKineticTerms[LNP+LSM,FlavorExpand->SU2W,SelectParticles->{{A,A},{Z1,Z1},{Z2,Z2},{Z1,A},{Z2,A},{Z1,Z2}}]

CheckKineticTermNormalisation[LNP+LSM,FlavorExpand->SU2W,SelectParticles->{{A,A},{Z1,Z1},{Z2,Z2}}]

Hermiticity

CheckHermiticity[LSM+LNP]

Mass

CheckDiagonalMassTerms[LNP+LSM,FlavorExpand->SU2W]

CheckMassSpectrum[LNP+LSM,FlavorExpand->SU2W,SelectParticles->{A,A}]

MassTest=CheckMassSpectrum[LNP+LSM,FlavorExpand->SU2W,SelectParticles->{{Z2,Z2},{Z1,Z1}}]

Feynman rules

FeynmanRules[LSM+LNP];

To MadGraph

SetDirectory["C:/users/Céline/Documents/UCL/Doctorat/Generic"];
WriteMGOutput[LSM+LNP,Output->"MC4BSMZ"]

Step 3 : SILH

Load FeynRules and the model

Quit[]

$FeynRulesPath =SetDirectory["C:\\Users\\Céline\\Documents\\UCL\\Doctorat\\FeynRules-current\\FeynRules-current-4"];
<< FeynRules`

LoadModel["C:/users/Céline/Documents/UCL/Doctorat/Generic//SILHUnphysZ.fr","C:/users/Céline/Documents/UCL/Doctorat/Generic//tutorial2.fr"]

Lagrangian and redefinitions

LH := cH/2/frho^2del[Phibar.Phi,mu]^2

GetKineticTerms[LHiggs+LH,SelectParticles->{H,H}]

Simplify[GetKineticTerms[LHiggs+LH,SelectParticles->{H,H}]/.H->Hren /Sqrt[1+cH v^2/frho^2]]

L6:=c6  λ/frho^2 (Phibar .Phi)^3

SelectFieldContent[LHiggs+L6,{H}]

muhsol=Solve[SelectFieldContent[LHiggs+L6,{H}]==0,muH]

GetMassTerms[LHiggs+L6,SelectParticles->{H,H}]

lsol=Solve[GetMassTerms[LHiggs+L6,SelectParticles->{H,H}]==mH^2/2 H^2/.muhsol[[2]],λ]/.{frho->v/ξ}

Series[λ/.lsol,{ξ,0,2}]

Series[muH^2/.muhsol[[2]]/.{frho->v/ξ},{ξ,0,2}]

Ly:=cy /frho^2 Phibar.Phi LYukawa

ttmass=CheckMassSpectrum[Ly+LYukawa,SelectParticles->{tt,ttbar},FlavorExpand->Generation]

yl[3]->Sqrt[2]ymtau/v(1-cy v^2/frho^2/2)

Step 4 : SILH

Load FeynRules and the model

Quit[]

$FeynRulesPath =SetDirectory["C:\\Users\\Céline\\Documents\\UCL\\Doctorat\\FeynRules-current\\FeynRules-current-4"];
<< FeynRules`

LoadModel["C:/users/Céline/Documents/UCL/Doctorat/Generic//SILHUnphysZ2.fr","C:/users/Céline/Documents/UCL/Doctorat/Generic//tutorial2.fr"]

LH := cH/2/frho^2del[Phibar.Phi,mu]^2
L6:=c6  λ/frho^2 (Phibar .Phi)^3
Ly:=cy /frho^2 Phibar.Phi LYukawa

Checks

Hermiticity

CheckHermiticity[LH+L6+Ly]

Kinetic

CheckKineticTermNormalisation[LHiggs+LH,SelectParticles->{Hren,Hren}]

Masses

CheckMassSpectrum[LHiggs+L6+LH,SelectParticles->{Hren,Hren}]

tmass=CheckMassSpectrum[Ly+LYukawa+L6+LH,SelectParticles->{t,tbar},FlavorExpand->{Generation,SU2W}]

bmass=CheckMassSpectrum[Ly+LYukawa+L6+LH,SelectParticles->{b,bbar},FlavorExpand->{Generation,SU2W}]

cmass=CheckMassSpectrum[Ly+LYukawa+L6+LH,SelectParticles->{c,cbar},FlavorExpand->{Generation,SU2W}]

taumass=CheckMassSpectrum[Ly+LYukawa+L6+LH,SelectParticles->{tt,ttbar},FlavorExpand->{Generation,SU2W}]

Feynman rules

LTot = Normal[Series[ExpandIndices[LSM+LNP+L6+LH+Ly,FlavorExpand->{SU2W}]/.{frho->v/ξ},{ξ,0,2}]]

FeynmanRules[LTot,FlavorExpand->{Generation,SU2W}]

To MadGraph

SetDirectory["C:/users/Céline/Documents/UCL/Doctorat/Generic"];
WriteMGOutput[LTot,MaxParticles->4]

To CalcHep

SetDirectory["C:/users/Céline/Documents/UCL/Doctorat/Generic"];
WriteCHOutput[LTot,MaxParticles->4]

SetDirectory["C:/users/Céline/Documents/UCL/Doctorat/Generic"];
WriteCHOutput[LTot,MaxParticles->4,MaxExpressionLength->122]