Quit[];

$FeynRulesPath=SetDirectory["/home/alba/Desktop/archivos_para_descarga/feynrules-current"]

<<FeynRules`

SetDirectory[NotebookDirectory[]]

(*En nuestro caso LoadModel[“SM.fr”,”Simple.fr”]*)

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

(*Las siguientes 3 instrucciones checan que el modelo no tenga problemas*)

(*En nuestro caso,Lagrangianos=LSM+Lax*)

CheckHermiticity[LSM+Lax]

CheckDiagonalQuadraticTerms[LSM+Lax]

CheckDiagonalMassTerms[LSM+Lax]

CheckMassSpectrum[LSM+Lax]

(*Ejemplos de Reglas de Feynman del Lagrangiano o una parte de él*)

FeynmanRules[Lax]//Simplify

DMFR=FeynmanRules[Lax,MaxParticles->4]//Simplify;

SMFR=FeynmanRules[LSM+Lax,SelectParticles->{H,H,H}]

FeynmanGauge=True;

(*Genera los archivos que necesita CalcHEP y micromegas para correr nuestro modelo*)

WriteCHOutput[LSM+Lax]

(*Genera los archivos que necesita MadGraph y maddm para correr nuestro modelo*)

FeynmanGauge=False;

WriteUFO[Lagrangiano]

(*Otros outputs*)

WriteFeynArtsOutput[Lagrangiano]