•Generation of coupling files

The amplitude is transformed to FeynArts notation:

mfa = FCToFA[melsimplified] ;

Terms are collected according to momenta:

mfacoll = MomentaCollect[mfa, ParticlesNumber -> 2, PerturbationOrder -> 1]

-1/(3 (f _ ϕ^(ó  ))^2) (i c _ 2^( ) p _ 3 _ μ _ 1 (-6 d _ (6 I _ 2 k1)^(3) d _ (I _ 1 I _ 3 k1)^(3) - 6 i f _ (6 I _ 2 k1)^(3) d _ (I _ 1 I _ 3 k1)^(3) + 3 d _ (6 I _ 1 k1)^(3) d _ (I _ 2 I _ 3 k1)^(3) + 3 i d _ (I _ 2 I _ 3 k1)^(3) f _ (6 I _ 1 k1)^(3) + 3 i d _ (I _ 1 I _ 2 k1)^(3) f _ (6 I _ 3 k1)^(3) + 3 d _ (6 I _ 3 k1)^(3) (d _ (I _ 1 I _ 2 k1)^(3) + i f _ (I _ 1 I _ 2 k1)^(3)) - 3 f _ (6 I _ 3 k1)^(3) f _ (I _ 1 I _ 2 k1)^(3) - 3 i d _ (6 I _ 1 k1)^(3) f _ (I _ 2 I _ 3 k1)^(3) + 3 f _ (6 I _ 1 k1)^(3) f _ (I _ 2 I _ 3 k1)^(3) + 2 δ _ (6 I _ 3) δ _ (I _ 1 I _ 2) - 4 δ _ (6 I _ 2) δ _ (I _ 1 I _ 3) + 2 δ _ (6 I _ 1) δ _ (I _ 2 I _ 3))) - 1/(3 (f _ ϕ^(ó  ))^2) (i c _ 2^( ) p _ 2 _ μ _ 1 (-6 d _ (6 I _ 3 k1)^(3) d _ (I _ 1 I _ 2 k1)^(3) - 6 i f _ (6 I _ 3 k1)^(3) d _ (I _ 1 I _ 2 k1)^(3) + 3 d _ (6 I _ 1 k1)^(3) d _ (I _ 2 I _ 3 k1)^(3) + 3 i d _ (I _ 2 I _ 3 k1)^(3) f _ (6 I _ 1 k1)^(3) + 3 i d _ (I _ 1 I _ 3 k1)^(3) f _ (6 I _ 2 k1)^(3) + 3 d _ (6 I _ 2 k1)^(3) (d _ (I _ 1 I _ 3 k1)^(3) + i f _ (I _ 1 I _ 3 k1)^(3)) - 3 f _ (6 I _ 2 k1)^(3) f _ (I _ 1 I _ 3 k1)^(3) + 3 i d _ (6 I _ 1 k1)^(3) f _ (I _ 2 I _ 3 k1)^(3) - 3 f _ (6 I _ 1 k1)^(3) f _ (I _ 2 I _ 3 k1)^(3) - 4 δ _ (6 I _ 3) δ _ (I _ 1 I _ 2) + 2 δ _ (6 I _ 2) δ _ (I _ 1 I _ 3) + 2 δ _ (6 I _ 1) δ _ (I _ 2 I _ 3)))

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

gencoup = GenericCoupling[mfacoll] ; gencoup

{p _ 3 _ μ _ 1, p _ 2 _ μ _ 1}

classcoup = ClassesCoupling[mfacoll] // Together ; classcoup // StandardForm

{{(CouplingConstant[ChPTW3[2], 1] (-3 i SU3D[6, I3, k1] SU3D[I1, I2, k1] + 6 i SU3D[6, I2, k1] SU3D[I1, I3, k1] - 3 i SU3D[6, I1, k1] SU3D[I2, I3, k1] + 3 SU3D[I2, I3, k1] SU3F[6, I1, k1] - 6 SU3D[I1, I3, k1] SU3F[6, I2, k1] + 3 SU3D[I1, I2, k1] SU3F[6, I3, k1] + 3 SU3D[6, I3, k1] SU3F[I1, I2, k1] + 3 i SU3F[6, I3, k1] SU3F[I1, I2, k1] - 3 SU3D[6, I1, k1] SU3F[I2, I3, k1] - 3 i SU3F[6, I1, k1] SU3F[I2, I3, k1] - 2 i SUNDelta[6, I3] SUNDelta[I1, I2] + 4 i SUNDelta[6, I2] SUNDelta[I1, I3] - 2 i SUNDelta[6, I1] SUNDelta[I2, I3]))/(3 DecayConstant[PseudoScalar[1], RenormalizationState[0]]^2)}, {(CouplingConstant[ChPTW3[2], 1] (6 i SU3D[6, I3, k1] SU3D[I1, I2, k1] - 3 i SU3D[6, I2, k1] SU3D[I1, I3, k1] - 3 i SU3D[6, I1, k1] SU3D[I2, I3, k1] + 3 SU3D[I2, I3, k1] SU3F[6, I1, k1] + 3 SU3D[I1, I3, k1] SU3F[6, I2, k1] - 6 SU3D[I1, I2, k1] SU3F[6, I3, k1] + 3 SU3D[6, I2, k1] SU3F[I1, I3, k1] + 3 i SU3F[6, I2, k1] SU3F[I1, I3, k1] + 3 SU3D[6, I1, k1] SU3F[I2, I3, k1] + 3 i SU3F[6, I1, k1] SU3F[I2, I3, k1] + 4 i SUNDelta[6, I3] SUNDelta[I1, I2] - 2 i SUNDelta[6, I2] SUNDelta[I1, I3] - 2 i SUNDelta[6, I1] SUNDelta[I2, I3]))/(3 DecayConstant[PseudoScalar[1], RenormalizationState[0]]^2)}}

$VeryVerbose = 2 ;

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

Using file name D:\\Program Files\\Wolfram Research\\Mathematica\\4.1\\AddOns\\Applications\\HighEnergyPhysics\\Phi\\CouplingVectors\\ChPTW3A00P10P10o2.Gen

File does not exist, evaluating

Saving

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

Using file name D:\\Program Files\\Wolfram Research\\Mathematica\\4.1\\AddOns\\Applications\\HighEnergyPhysics\\Phi\\CouplingVectors\\ChPTW3A00P10P10o2.Mod

File does not exist, evaluating

Saving

$VeryVerbose = 0 ;

Converted by Mathematica  (July 10, 2003)