•ϕϕϕϕϕ

ll = ArgumentsSupply[Lagrangian[ChPTW3[2]], x, RenormalizationState[0], ExpansionOrder -> 5, DropOrder -> 5, DiagonalToU -> True] ;

Redundant terms are discarded:

$VeryVerbose = 2 ;

lll = DiscardTerms[ll, Retain -> {Particle[PhiMeson , RenormalizationState[0]] -> 5}, CommutatorReduce -> True, Method -> Expand] // Simplify

DeclareUScalar[QuantumField[aa___, Particle[bb__], cc___][x_]] ;

$VeryVerbose = 0 ;

lala = ($IsoIndicesCounter = 0 ; PhiToFC[IsoIndicesSupply[#]]) & /@ Expand[lll] ;

DeclareNonCommutative[UMatrix[UGenerator[SUNIndex[j_], opt___], opts___]] ; DeclareNonCommutative[UMatrix[UGenerator[ExplicitSUNIndex[j_], opt___], opts___]] ;

Calculation of the Feynman rule:

fields = {QuantumField[Particle[PhiMeson, RenormalizationState[0]], SUNIndex[I1]][p1], QuantumField[Particle[PhiMeson, RenormalizationState[0]], SUNIndex[I2]][p2], QuantumField[Particle[PhiMeson, RenormalizationState[0]], SUNIndex[I3]][p3], QuantumField[Particle[PhiMeson, RenormalizationState[0]], SUNIndex[I4]][p4], QuantumField[Particle[PhiMeson, RenormalizationState[0]], SUNIndex[I5]][p5]}

{ϕ^( )^I _ 1, ϕ^( )^I _ 2, ϕ^( )^I _ 3, ϕ^( )^I _ 4, ϕ^( )^I _ 5}

lala // Length

18

ii = 0 ; res = CheckF[(++ ii ; WriteString["stdout", ii, " "] ; FeynRule[#, fields]) & /@ lala, "5MesonVertex"] ;

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18

resu = CheckF[Collect[res, _UTrace1], "5resu"] ;

resul = CheckF[Collect[resu // IndicesCleanup, HoldPattern[Plus[__ ? ((! FreeQ[{##}, Momentum | ParticleMass, Infinity, Heads -> True]) &)]]], "5resul"] ;

Share[]

13405608

Converted by Mathematica  (July 10, 2003)