•The kaon

zkaon = CheckF[dum, "ChPT3P60o2.Fac"]

Using file name D:\\Program Files\\Wolfram Research\\Mathematica\\4.1\\AddOns\\Applications\\HighEnergyPhysics\\Phi\\Factors\\ChPT3P60o2.Fac

File exists, loading

1 - 1/(64 π^2 (f _ ϕ^(ó  ))^2) ((32 π^2 λ + log((m _ π^(ó  ))^2/μ^2)) (m _ π^(ó  ))^2 + 2 (32 π^2 (λ - 8 L _ 5^( )) + log((m _ K^(ó  ))^2/μ^2)) (m _ K^(ó  ))^2 + (32 π^2 λ + log((m _ η^(ó  ))^2/μ^2)) (m _ η^(ó  ))^2 - 512 π^2 L _ 4^( ) ((m _ π^(ó  ))^2 + 2 (m _ K^(ó  ))^2))

ff2 = amp2 /. p2 -> -p1

-i f _ ϕ^(ó  ) p _ 1^μ _ 1

amploop = ampinfinities /. ρ1 -> μ1 /. i1 -> 4 // SUNReduce // Simplify

(i p _ 1^μ _ 1 ((32 π^2 λ + log((m _ π^(ó  ))^2/μ^2)) (m _ π^(ó  ))^2 + 2 (32 π^2 λ + log((m _ K^(ó  ))^2/μ^2)) (m _ K^(ó  ))^2 + (32 π^2 λ + log((m _ η^(ó  ))^2/μ^2)) (m _ η^(ó  ))^2))/(32 π^2 f _ ϕ^(ó  ))

ampwf4 = amp4 /. p2 -> -p1 /. I1 -> 4 // SUNReduce // Simplify

-(8 i p _ 1^μ _ 1 (L _ 5^( ) (m _ K^(ó  ))^2 + L _ 4^( ) ((m _ π^(ó  ))^2 + 2 (m _ K^(ó  ))^2)))/f _ ϕ^(ó  )

amppf4 = ff2 (1 + (2 - zkaon))/2 // Simplify

-1/2 i f _ ϕ^(ó  ) p _ 1^μ _ 1 (1/(64 π^2 (f _ ϕ^(ó  ))^2) ((32 π^2 λ + log((m _ π^(ó  ))^2/μ^2)) (m _ π^(ó  ))^2 + 2 (32 π^2 (λ - 8 L _ 5^( )) + log((m _ K^(ó  ))^2/μ^2)) (m _ K^(ó  ))^2 + (32 π^2 λ + log((m _ η^(ó  ))^2/μ^2)) (m _ η^(ó  ))^2 - 512 π^2 L _ 4^( ) ((m _ π^(ó  ))^2 + 2 (m _ K^(ó  ))^2)) + 2)

ff4 = (amploop + ampwf4 + amppf4 /. gellmannOkubo // ExpandAll // Simplify) /. etalogs

-1/(128 π^2 f _ ϕ^(ó  )) (i p _ 1^μ _ 1 (128 π^2 (f _ ϕ^(ó  ))^2 - 64 π^2 λ (m _ π^(ó  ))^2 - 3 log((m _ π^(ó  ))^2/μ^2) (m _ π^(ó  ))^2 + log((m _ η^(ó  ))^2/μ^2) (m _ π^(ó  ))^2 + 512 π^2 L _ 5^( ) (m _ K^(ó  ))^2 - 320 π^2 λ (m _ K^(ó  ))^2 - 6 log((m _ K^(ó  ))^2/μ^2) (m _ K^(ó  ))^2 - 4 log((m _ η^(ó  ))^2/μ^2) (m _ K^(ó  ))^2 + 512 π^2 L _ 4^( ) ((m _ π^(ó  ))^2 + 2 (m _ K^(ó  ))^2)))

ff0 = Renormalize[ff4] // Simplify

-1/(128 π^2 f _ ϕ^(ó  )) (i p _ 1^μ _ 1 (128 π^2 (f _ ϕ^(ó  ))^2 - 3 log((m _ π^(ó  ))^2/μ^2) (m _ π^(ó  ))^2 + log((m _ η^(ó  ))^2/μ^2) (m _ π^(ó  ))^2 + 512 π^2 L _ 5^(r ) (m _ K^(ó  ))^2 - 6 log((m _ K^(ó  ))^2/μ^2) (m _ K^(ó  ))^2 - 4 log((m _ η^(ó  ))^2/μ^2) (m _ K^(ó  ))^2 + 512 π^2 L _ 4^(r ) ((m _ π^(ó  ))^2 + 2 (m _ K^(ó  ))^2)))

c = Collect[Coefficient[-ff0/I/DecayConstant[PhiMeson, RenormalizationState[0]], Pair[LorentzIndex[μ1], Momentum[p1]]], _DecayConstant]

1/(128 π^2 (f _ ϕ^(ó  ))^2) (-3 log((m _ π^(ó  ))^2/μ^2) (m _ π^(ó  ))^2 + log((m _ η^(ó  ))^2/μ^2) (m _ π^(ó  ))^2 + 512 π^2 L _ 5^(r ) (m _ K^(ó  ))^2 - 6 log((m _ K^(ó  ))^2/μ^2) (m _ K^(ó  ))^2 - 4 log((m _ η^(ó  ))^2/μ^2) (m _ K^(ó  ))^2 + 512 π^2 L _ 4^(r ) ((m _ π^(ó  ))^2 + 2 (m _ K^(ó  ))^2)) + 1

CheckF[c, "ChPT3A00P60o2.Fac"] ;

Using file name D:\\Program Files\\Wolfram Research\\Mathematica\\4.1\\AddOns\\Applications\\HighEnergyPhysics\\Phi\\Factors\\ChPT3A00P60o2.Fac

File exists, loading

Converted by Mathematica  (July 10, 2003)