•Generation of coupling files

The amplitude calculated is transformed to FeynArts notation:

mfa = FCToFA[resul] ;

Terms are collected according to momenta. Here are many more than in a previous calculation CHECK!!!

mfacoll = (Collect[mfa // Expand, {FourVector[__] * ParticleMass[__]^2, FourVector[__] * MomentaScalarProduct[__]}] /. {FourVector[a__] * ParticleMass[b__]^2 * c_ :> FourVector[a] * ParticleMass[b]^2 * Collect[c, {_DecayConstant, _UTrace1}], FourVector[a__] * MomentaScalarProduct[b__] * c_ :> FourVector[a] * MomentaScalarProduct[b] * Collect[c, {_DecayConstant, _UTrace1}]}) ;

mfacoll

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

gencoup = GenericCoupling[mfacoll] ; gencoup

{p _ 3 _ μ _ 1 (p _ 2 ° p _ 4), p _ 2 _ μ _ 1 (p _ 3 ° p _ 4), p _ 4 _ μ _ 1 (p _ 2 ° p _ 3), p _ 3 _ μ _ 1 (p _ 3 ° p _ 4), p _ 3 _ μ _ 1 (p _ 1 ° p _ 4), p _ 4 _ μ _ 1 (p _ 1 ° p _ 3), p _ 2 _ μ _ 1 (p _ 2 ° p _ 4), p _ 2 _ μ _ 1 (p _ 1 ° p _ 4), p _ 4 _ μ _ 1 (p _ 1 ° p _ 2), p _ 1 _ σ p _ 3 _ ν p _ 4 _ μ g^(μ _ 1 ρ), p _ 1 _ ρ p _ 3 _ ν p _ 4 _ μ g^(μ _ 1 σ), p _ 1 _ σ p _ 2 _ ν p _ 4 _ μ g^(μ _ 1 ρ), p _ 1 _ ρ p _ 2 _ ν p _ 4 _ μ g^(μ _ 1 σ), p _ 1 _ σ p _ 2 _ ν p _ 4 _ μ g^(μ _ 1 ρ), p _ 1 _ ρ p _ 2 _ ν p _ 4 _ μ g^(μ _ 1 σ), p _ 1 _ σ p _ 2 _ ν p _ 4 _ μ g^(μ _ 1 ρ), p _ 1 _ σ p _ 3 _ ν p _ 4 _ μ g^(μ _ 1 ρ), p _ 1 _ ρ p _ 2 _ ν p _ 4 _ μ g^(μ _ 1 σ), p _ 1 _ ρ p _ 3 _ ν p _ 4 _ μ g^(μ _ 1 σ), p _ 1 _ σ p _ 3 _ ν p _ 4 _ μ g^(μ _ 1 ρ), p _ 1 _ ρ p _ 3 _ ν p _ 4 _ μ g^(μ _ 1 σ), p _ 1 _ σ p _ 2 _ ν p _ 4 _ μ g^(μ _ 1 ρ), p _ 1 _ σ p _ 3 _ ν p _ 4 _ μ g^(μ _ 1 ρ), p _ 1 _ ρ p _ 2 _ ν p _ 4 _ μ g^(μ _ 1 σ), p _ 1 _ ρ p _ 3 _ ν p _ 4 _ μ g^(μ _ 1 σ), p _ 1 _ σ p _ 3 _ ν p _ 4 _ μ g^(μ _ 1 ρ), p _ 1 _ ρ p _ 3 _ ν p _ 4 _ μ g^(μ _ 1 σ), p _ 1 _ σ p _ 3 _ ν p _ 4 _ μ g^(μ _ 1 ρ), p _ 1 _ ρ p _ 3 _ ν p _ 4 _ μ g^(μ _ 1 σ), p _ 1 _ σ p _ 2 _ ν p _ 4 _ μ g^(μ _ 1 ρ), p _ 1 _ ρ p _ 2 _ ν p _ 4 _ μ g^(μ _ 1 σ), p _ 1 _ σ p _ 2 _ ν p _ 4 _ μ g^(μ _ 1 ρ), p _ 1 _ ρ p _ 2 _ ν p _ 4 _ μ g^(μ _ 1 σ), p _ 4 _ μ _ 1 (m _ K^0^(ó  ))^2, p _ 4 _ μ _ 1 (m _ K^+^(ó  ))^2, p _ 4 _ μ _ 1 (m _ π^(ó  ))^2}

classcoup = Simplify /@ ClassesCoupling[mfacoll] ;

classcoup // LeafCount

12257

classcoup // Length

36

$VeryVerbose = 2 ;

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

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

File exists, force evaluating

Saving

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

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

File exists, force evaluating

Saving

$VeryVerbose = 0 ;

Converted by Mathematica  (July 10, 2003)