Tutorial

Quit[]

Loading FeynRules

$FeynRulesPath = SetDirectory["~/feynrules/branch/feynrules-current"]

<<FeynRules`

Loading the model

SetDirectory["~/feynrules/branch/feynrules-current/Models/Tutorial"]

LoadModel["SM.fr","Tutorial.fr"]

Restriction files, to make the quarks massless, and the CKM matrix diagonal:

LoadRestriction["DiagonalCKM.rst","Massless.rst"]

The Lagrangian

LScalarKin=1/2del[pi1,mu]del[pi1,mu]+1/2del[pi2,mu]del[pi2,mu]-MM1^2/2pi1^2-MM2^2/2pi2^2-MM12^2pi1 pi2

LFermionKin=I uvbar.Ga[mu].DC[uv,mu]-Muv uvbar.uv+I evbar.Ga[mu].DC[ev,mu]-Mev evbar.ev

LDecays=lam1 pi1 uvbar[ss1,ii].uq[ss2,3,ii]ProjP[ss1,ss2]+lam2 pi2 uvbar[ss1,ii].uq[ss2,3,ii]ProjP[ss1,ss2]+lam1p pi1 evbar[ss1].l[ss2,1]ProjP[ss1,ss2]+lam2p pi2 evbar[ss1].l[ss2,1]ProjP[ss1,ss2]
LDecays=LDecays+HC[LDecays]

LNew=LScalarKin+LFermionKin+LDecays

Checking the hermiticity of the model

CheckMassSpectrum[LNew]

The Feynman rules

vertices=FeynmanRules[LNew]

Computing decays

decays=ComputeWidths[vertices];

x1=PartialWidth[{uv,t,p1}]
x2=PartialWidth[{uv,t,p2}]
x3=PartialWidth[{p2,e,evbar}]
x4=PartialWidth[{ev,e,p1}]

NumericalValue[x1]
NumericalValue[x2]
NumericalValue[x3]
NumericalValue[x4]

UFO

SetDirectory["~/mg5amcnlo/models"]
WriteUFO[LSM+LNew]