•Generation of coupling files

The amplitude is transformed to FeynArts notation:

mfa = FCToFA[res] ;

Terms are collected according to momenta:

mfacoll = MomentaCollect[mfa // Expand, ParticlesNumber -> 2, PerturbationOrder -> 4]

(p _ 1 ° p _ 2) ((4 c _ 2^( ) N _ 21^( ) < σ^6 . σ^I _ 1 . σ^8 >)/(3^(1/2) (f _ ϕ^(ó  ))^3) - (4 c _ 2^( ) N _ 21^( ) < σ^6 . σ^8 . σ^I _ 1 >)/(3^(1/2) (f _ ϕ^(ó  ))^3)) (m _ π^(ó  ))^2 + (p _ 1 ° p _ 2) (m _ K^0^(ó  ))^2 ((2 c _ 2^( ) N _ 21^( ) < σ^6 . σ^3 . σ^I _ 1 >)/(f _ ϕ^(ó  ))^3 + (2 c _ 2^( ) N _ 21^( ) < σ^6 . σ^8 . σ^I _ 1 >)/(3^(1/2) (f _ ϕ^(ó  ))^3) - (2 c _ 2^( ) N _ 21^( ) < σ^6 . σ^I _ 1 . σ^3 >)/(f _ ϕ^(ó  ))^3 - (2 c _ 2^( ) N _ 21^( ) < σ^6 . σ^I _ 1 . σ^8 >)/(3^(1/2) (f _ ϕ^(ó  ))^3)) + (p _ 1 ° p _ 2) (m _ K^+^(ó  ))^2 (-(2 c _ 2^( ) N _ 21^( ) < σ^6 . σ^3 . σ^I _ 1 >)/(f _ ϕ^(ó  ))^3 + (2 c _ 2^( ) N _ 21^( ) < σ^6 . σ^8 . σ^I _ 1 >)/(3^(1/2) (f _ ϕ^(ó  ))^3) + (2 c _ 2^( ) N _ 21^( ) < σ^6 . σ^I _ 1 . σ^3 >)/(f _ ϕ^(ó  ))^3 - (2 c _ 2^( ) N _ 21^( ) < σ^6 . σ^I _ 1 . σ^8 >)/(3^(1/2) (f _ ϕ^(ó  ))^3))

The coupling vectors for the FeynArts model files are generated and saved:

gencoup = GenericCoupling[mfacoll] ; gencoup

{(p _ 1 ° p _ 2) (m _ K^0^(ó  ))^2, (p _ 1 ° p _ 2) (m _ K^+^(ó  ))^2, (p _ 1 ° p _ 2) (m _ π^(ó  ))^2}

classcoup = Simplify /@ ClassesCoupling[mfacoll] ;

$VeryVerbose = 2 ;

CheckF[gencoup, XName[VertexFields -> {PhiMeson[0], Scalar[1][0]}, PerturbationOrder -> 4, PhiModel -> ChPTW3, XFileName -> Automatic] <> ".Gen"] ;

Using file name /home/fjob/.Mathematica/Applications/HighEnergyPhysics/Phi/CouplingVectors/ChPTW3P10S10o4.Gen

File exists, force evaluating

Saving

CheckF[classcoup, XName[VertexFields -> {PhiMeson[0], Scalar[1][0]}, PerturbationOrder -> 4, PhiModel -> ChPTW3, XFileName -> Automatic] <> ".Mod"] ;

Using file name /home/fjob/.Mathematica/Applications/HighEnergyPhysics/Phi/CouplingVectors/ChPTW3P10S10o4.Mod

File exists, force evaluating

Saving

$VeryVerbose = 0 ;

Converted by Mathematica  (July 10, 2003)