•ππππ

IsoVector[QuantumField[Particle[AxialVector[0], ___], ___], ___][_] := 0 ; IsoVector[QuantumField[Particle[Vector[0], ___], ___], ___][_] := 0 ; IsoVector[QuantumField[Particle[Vector[1], ___], ___], ___][_] := 0 ; QuantumField[Particle[Vector[1], ___], ___][_] := 0 ; IsoVector[QuantumField[Particle[Scalar[1 | 2], ___], ___], ___][_] := 0 ; QuantumField[Particle[Scalar[1 | 2], ___], ___][_] := 0 ; IsoVector[QuantumField[Particle[PseudoScalar[0], ___], ___], ___][_] := 0 ; QuantumField[Particle[PseudoScalar[0], ___], ___][_] := 0 ;

LoadLagrangian[ChPT2[4]] ;

lag = Lagrangian[ChPT2[4]] + Lagrangian[ChPTEM2[4]] /. {CouplingConstant[ChPT2[4], _ ? (FreeQ[#, 1 | 2 | 3 | 4 | 5 | 6 | 8] &), ___] -> 0, CouplingConstant[ChPTEM2[4], _ ? (FreeQ[#, 2 | 3 | 4 | 7 | 10 | 11] &), ___] -> 0}

k _ 11^( ) (< Q '6 Q > '6 (< ÷„ '6 χ^† > + < χ '6 ÷„^† >)) (f _ π^(ó  ))^2 + k _ 10^( ) (< Q '6 Q > '6 < ÷s _ μ(÷„) '6 ÷s _ μ(÷„)^† >) (f _ π^(ó  ))^2 + k _ 4^( ) (< ÷„^† '6 ÷s _ μ(÷„) '6 Q > '6 < ÷s _ μ(÷„) '6 ÷„^† '6 Q >) (f _ π^(ó  ))^2 + k _ 3^( ) (< ÷„^† '6 ÷s _ μ(÷„) '6 Q > '6 < ÷„^† '6 ÷s _ μ(÷„) '6 Q > + < ÷s _ μ(÷„) '6 ÷„^† '6 Q > '6 < ÷s _ μ(÷„) '6 ÷„^† '6 Q >) (f _ π^(ó  ))^2 + k _ 7^( ) (< Q '6 ÷„ '6 Q '6 ÷„^† > '6 (< χ '6 ÷„^† > + < χ^† '6 ÷„ >)) (f _ π^(ó  ))^2 + k _ 2^( ) (< Q '6 ÷„ '6 Q '6 ÷„^† > '6 < ÷s _ μ(÷„) '6 ÷s _ μ(÷„)^† >) (f _ π^(ó  ))^2 + L _ 6^( ) ((< ÷„ '6 χ^† > + < χ '6 ÷„^† >) '6 (< ÷„ '6 χ^† > + < χ '6 ÷„^† >)) + L _ 4^( ) (< ÷s _ μ(÷„) '6 ÷s _ μ(÷„)^† > '6 (< ÷„ '6 χ^† > + < χ '6 ÷„^† >)) + L _ 1^( ) (< ÷s _ μ(÷„) '6 ÷s _ μ(÷„)^† > '6 < ÷s _ ν(÷„) '6 ÷s _ ν(÷„)^† >) + L _ 2^( ) (< ÷s _ μ(÷„) '6 ÷s _ ν(÷„)^† > '6 < ÷s _ μ(÷„) '6 ÷s _ ν(÷„)^† >) + L _ 5^( ) < ÷s _ μ(÷„) '6 ÷s _ μ(÷„)^† '6 (÷„ '6 χ^† + χ '6 ÷„^†) > + L _ 8^( ) (< ÷„ '6 χ^† '6 ÷„ '6 χ^† > + < χ '6 ÷„^† '6 χ '6 ÷„^† >) + L _ 3^( ) < ÷s _ μ(÷„) '6 ÷s _ μ(÷„)^† '6 ÷s _ ν(÷„) '6 ÷s _ ν(÷„)^† >

ll = (WriteString["stdout", "."] ; UNMSplit[#, x, DropOrder -> 4]) & /@ lag ;

.............

lll = ArgumentsSupply[ll /. {UMatrix[UChiralSpurionRight] -> UQuarkChargeMatrix[RenormalizationState[0], DiagonalToU -> True], UMatrix[UChiralSpurionLeft] -> UQuarkChargeMatrix[RenormalizationState[0], DiagonalToU -> True], UMatrix[UChiralSpurion] -> UQuarkChargeMatrix[RenormalizationState[0], DiagonalToU -> True]}, x, RenormalizationState[0], DiagonalToU -> True, ExpansionOrder -> 4, DropOrder -> 4] ;

ArgumentsSupply :: argxpr : Warning : The argument x is already in the expression.

llle = (WriteString["stdout", "."] ; ExpandU[#] // CommutatorReduce // Simplify) & /@ lll ;

.............

llll = ($IsoIndicesCounter = 0 ; WriteString["stdout", "."] ; # // IsoIndicesSupply // SUNReduce[#, FullReduce -> True] & // IndicesCleanup // CommutatorReduce[#, FullReduce -> True] & // Simplify) & /@ llle ;

.............

lala = Simplify /@ Collect[llll /. NM -> Times, CouplingConstant[_[4], __]] ;

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

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

mel = ((WriteString["stdout", "."] ; I * FunctionalD[PhiToFC[#], fields]) & /@ lala) ;

.............

melsimplified = Collect[mel // Contract // SUNReduce[#, FullReduce -> True] &, {_DecayConstant, _CouplingConstant, _Pair}] ;

mfacoll = MomentaCollect[melsimplified // Expand, ParticlesNumber -> 4, PerturbationOrder -> 4, ScalarProductForm -> Pair] ;

gencoup = GenericCoupling[mfacoll] ; gencoup

{(p _ 1 ° p _ 4) (p _ 2 ° p _ 3), (p _ 1 ° p _ 3) (p _ 2 ° p _ 4), (p _ 1 ° p _ 2) (p _ 3 ° p _ 4), (p _ 1 ° p _ 4) (m _ π^(ó  ))^2, (p _ 2 ° p _ 3) (m _ π^(ó  ))^2, (p _ 1 ° p _ 3) (m _ π^(ó  ))^2, (p _ 2 ° p _ 4) (m _ π^(ó  ))^2, (p _ 1 ° p _ 2) (m _ π^(ó  ))^2, (p _ 3 ° p _ 4) (m _ π^(ó  ))^2, (m _ π^(ó  ))^4, (e^( ))^2 (p _ 2 ° p _ 3), (e^( ))^2 (p _ 1 ° p _ 4), (e^( ))^2 (p _ 2 ° p _ 4), (e^( ))^2 (p _ 1 ° p _ 3), (e^( ))^2 (p _ 3 ° p _ 4), (e^( ))^2 (p _ 1 ° p _ 2), (e^( ))^2 (m _ π^(ó  ))^2}

classcoup = ClassesCoupling[mfacoll] // Together ;

$VeryVerbose = 2 ;

CheckF[gencoup, XName[VertexFields -> {PseudoScalar[2][0], PseudoScalar[2][0], PseudoScalar[2][0], PseudoScalar[2][0]}, PerturbationOrder -> 4, PhiModel -> ChPTEM2, XFileName -> Automatic] <> ".Gen"] ;

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

File exists, loading

CheckF[classcoup, XName[VertexFields -> {PseudoScalar[2][0], PseudoScalar[2][0], PseudoScalar[2][0], PseudoScalar[2][0]}, PerturbationOrder -> 4, PhiModel -> ChPTEM2, XFileName -> Automatic] <> ".Mod"] ;

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

File exists, loading

$VeryVerbose = 0 ;

Converted by Mathematica  (July 10, 2003)