•Generation of coupling files

mfa = FCToFA[melsimplified1] ;

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

-1/(36 f _ ϕ^(ó  )) (p _ 2 _ μ _ 1 (15 d _ (I _ 1 I _ 4 k1)^(3) d _ (I _ 2 I _ 3 k1)^(3) + 15 d _ (I _ 1 I _ 3 k1)^(3) d _ (I _ 2 I _ 4 k1)^(3) - 30 d _ (I _ 1 I _ 2 k1)^(3) d _ (I _ 3 I _ 4 k1)^(3) - 9 f _ (I _ 1 I _ 4 k1)^(3) f _ (I _ 2 I _ 3 k1)^(3) - 9 f _ (I _ 1 I _ 3 k1)^(3) f _ (I _ 2 I _ 4 k1)^(3) + 10 δ _ (I _ 1 I _ 4) δ _ (I _ 2 I _ 3) + 10 δ _ (I _ 1 I _ 3) δ _ (I _ 2 I _ 4) - 20 δ _ (I _ 1 I _ 2) δ _ (I _ 3 I _ 4))) - 1/(36 f _ ϕ^(ó  )) (p _ 3 _ μ _ 1 (15 d _ (I _ 1 I _ 4 k1)^(3) d _ (I _ 2 I _ 3 k1)^(3) - 30 d _ (I _ 1 I _ 3 k1)^(3) d _ (I _ 2 I _ 4 k1)^(3) + 15 d _ (I _ 1 I _ 2 k1)^(3) d _ (I _ 3 I _ 4 k1)^(3) + 9 f _ (I _ 1 I _ 4 k1)^(3) f _ (I _ 2 I _ 3 k1)^(3) - 9 f _ (I _ 1 I _ 2 k1)^(3) f _ (I _ 3 I _ 4 k1)^(3) + 10 δ _ (I _ 1 I _ 4) δ _ (I _ 2 I _ 3) - 20 δ _ (I _ 1 I _ 3) δ _ (I _ 2 I _ 4) + 10 δ _ (I _ 1 I _ 2) δ _ (I _ 3 I _ 4))) - 1/(36 f _ ϕ^(ó  )) (p _ 4 _ μ _ 1 (-30 d _ (I _ 1 I _ 4 k1)^(3) d _ (I _ 2 I _ 3 k1)^(3) + 15 d _ (I _ 1 I _ 3 k1)^(3) d _ (I _ 2 I _ 4 k1)^(3) + 15 d _ (I _ 1 I _ 2 k1)^(3) d _ (I _ 3 I _ 4 k1)^(3) + 9 f _ (I _ 1 I _ 3 k1)^(3) f _ (I _ 2 I _ 4 k1)^(3) + 9 f _ (I _ 1 I _ 2 k1)^(3) f _ (I _ 3 I _ 4 k1)^(3) - 20 δ _ (I _ 1 I _ 4) δ _ (I _ 2 I _ 3) + 10 δ _ (I _ 1 I _ 3) δ _ (I _ 2 I _ 4) + 10 δ _ (I _ 1 I _ 2) δ _ (I _ 3 I _ 4)))

gencoup = GenericCoupling[mfacoll] ; gencoup

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

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

{{(-15*SU3D[I1, I4, k1]*SU3D[I2, I3, k1] - 15*SU3D[I1, I3, k1]*SU3D[I2, I4, k1] +
    30*SU3D[I1, I2, k1]*SU3D[I3, I4, k1] + 9*SU3F[I1, I4, k1]*SU3F[I2, I3, k1] +
    9*SU3F[I1, I3, k1]*SU3F[I2, I4, k1] - 10*SUNDelta[I1, I4]*SUNDelta[I2, I3] -
    10*SUNDelta[I1, I3]*SUNDelta[I2, I4] + 20*SUNDelta[I1, I2]*SUNDelta[I3, I4])/
   (36*DecayConstant[PseudoScalar[1], RenormalizationState[0]])},
{(-15*SU3D[I1, I4, k1]*SU3D[I2, I3, k1] + 30*SU3D[I1, I3, k1]*SU3D[I2, I4, k1] -
    15*SU3D[I1, I2, k1]*SU3D[I3, I4, k1] - 9*SU3F[I1, I4, k1]*SU3F[I2, I3, k1] +
    9*SU3F[I1, I2, k1]*SU3F[I3, I4, k1] - 10*SUNDelta[I1, I4]*SUNDelta[I2, I3] +
    20*SUNDelta[I1, I3]*SUNDelta[I2, I4] - 10*SUNDelta[I1, I2]*SUNDelta[I3, I4])/
   (36*DecayConstant[PseudoScalar[1], RenormalizationState[0]])},
{(30*SU3D[I1, I4, k1]*SU3D[I2, I3, k1] - 15*SU3D[I1, I3, k1]*SU3D[I2, I4, k1] -
    15*SU3D[I1, I2, k1]*SU3D[I3, I4, k1] - 9*SU3F[I1, I3, k1]*SU3F[I2, I4, k1] -
    9*SU3F[I1, I2, k1]*SU3F[I3, I4, k1] + 20*SUNDelta[I1, I4]*SUNDelta[I2, I3] -
    10*SUNDelta[I1, I3]*SUNDelta[I2, I4] - 10*SUNDelta[I1, I2]*SUNDelta[I3, I4])/
   (36*DecayConstant[PseudoScalar[1], RenormalizationState[0]])}}

$VeryVerbose = 2 ;

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

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

File exists, force evaluating

Saving

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

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

File exists, force evaluating

Saving

$VeryVerbose = 0 ;

Converted by Mathematica  (July 10, 2003)