•Generation of coupling files

Terms are collected according to momenta:

mfacoll = MomentaCollect[melsimplified, ParticlesNumber -> 4, PerturbationOrder -> 2, ScalarProductForm -> Pair]

(i (δ _ (I _ 1 I _ 4) δ _ (I _ 2 I _ 3) + δ _ (I _ 1 I _ 3) δ _ (I _ 2 I _ 4) + δ _ (I _ 1 I _ 2) δ _ (I _ 3 I _ 4)) (m _ π^(ó  ))^2)/(3 (f _ π^(ó  ))^2) - (i p _ 1 · p _ 2 (δ _ (I _ 1 I _ 4) δ _ (I _ 2 I _ 3) + δ _ (I _ 1 I _ 3) δ _ (I _ 2 I _ 4) - 2 δ _ (I _ 1 I _ 2) δ _ (I _ 3 I _ 4)))/(3 (f _ π^(ó  ))^2) - (i p _ 3 · p _ 4 (δ _ (I _ 1 I _ 4) δ _ (I _ 2 I _ 3) + δ _ (I _ 1 I _ 3) δ _ (I _ 2 I _ 4) - 2 δ _ (I _ 1 I _ 2) δ _ (I _ 3 I _ 4)))/(3 (f _ π^(ó  ))^2) - (i p _ 1 · p _ 3 (δ _ (I _ 1 I _ 4) δ _ (I _ 2 I _ 3) - 2 δ _ (I _ 1 I _ 3) δ _ (I _ 2 I _ 4) + δ _ (I _ 1 I _ 2) δ _ (I _ 3 I _ 4)))/(3 (f _ π^(ó  ))^2) - (i p _ 2 · p _ 4 (δ _ (I _ 1 I _ 4) δ _ (I _ 2 I _ 3) - 2 δ _ (I _ 1 I _ 3) δ _ (I _ 2 I _ 4) + δ _ (I _ 1 I _ 2) δ _ (I _ 3 I _ 4)))/(3 (f _ π^(ó  ))^2) - (i p _ 1 · p _ 4 (-2 δ _ (I _ 1 I _ 4) δ _ (I _ 2 I _ 3) + δ _ (I _ 1 I _ 3) δ _ (I _ 2 I _ 4) + δ _ (I _ 1 I _ 2) δ _ (I _ 3 I _ 4)))/(3 (f _ π^(ó  ))^2) - (i p _ 2 · p _ 3 (-2 δ _ (I _ 1 I _ 4) δ _ (I _ 2 I _ 3) + δ _ (I _ 1 I _ 3) δ _ (I _ 2 I _ 4) + δ _ (I _ 1 I _ 2) δ _ (I _ 3 I _ 4)))/(3 (f _ π^(ó  ))^2)

The coupling vectors are generated and saved:

gencoup = GenericCoupling[mfacoll] ; gencoup

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

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

{{-(i (SUNDelta[I1, I4] SUNDelta[I2, I3] + SUNDelta[I1, I3] SUNDelta[I2, I4] - 2 SUNDelta[I1, I2] SUNDelta[I3, I4]))/(3 DecayConstant[PseudoScalar[2], RenormalizationState[0]]^2)}, {-(i (SUNDelta[I1, I4] SUNDelta[I2, I3] + SUNDelta[I1, I3] SUNDelta[I2, I4] - 2 SUNDelta[I1, I2] SUNDelta[I3, I4]))/(3 DecayConstant[PseudoScalar[2], RenormalizationState[0]]^2)}, {-(i (SUNDelta[I1, I4] SUNDelta[I2, I3] - 2 SUNDelta[I1, I3] SUNDelta[I2, I4] + SUNDelta[I1, I2] SUNDelta[I3, I4]))/(3 DecayConstant[PseudoScalar[2], RenormalizationState[0]]^2)}, {-(i (SUNDelta[I1, I4] SUNDelta[I2, I3] - 2 SUNDelta[I1, I3] SUNDelta[I2, I4] + SUNDelta[I1, I2] SUNDelta[I3, I4]))/(3 DecayConstant[PseudoScalar[2], RenormalizationState[0]]^2)}, {(i (2 SUNDelta[I1, I4] SUNDelta[I2, I3] - SUNDelta[I1, I3] SUNDelta[I2, I4] - SUNDelta[I1, I2] SUNDelta[I3, I4]))/(3 DecayConstant[PseudoScalar[2], RenormalizationState[0]]^2)}, {(i (2 SUNDelta[I1, I4] SUNDelta[I2, I3] - SUNDelta[I1, I3] SUNDelta[I2, I4] - SUNDelta[I1, I2] SUNDelta[I3, I4]))/(3 DecayConstant[PseudoScalar[2], RenormalizationState[0]]^2)}, {(i (SUNDelta[I1, I4] SUNDelta[I2, I3] + SUNDelta[I1, I3] SUNDelta[I2, I4] + SUNDelta[I1, I2] SUNDelta[I3, I4]))/(3 DecayConstant[PseudoScalar[2], RenormalizationState[0]]^2)}}

With XFileName->Automatic, the coupling vectors are saved under names generated from the options specified or defaults.

$VeryVerbose = 2 ;

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

Utilities |

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

File exists, loading

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

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

File exists, loading

$VeryVerbose = 0 ;

Converted by Mathematica  (July 10, 2003)