•Renormalization

The beta functions:

Cases[amp4, CouplingConstant[__], Infinity] // Union // Renormalize

{e^(  ), L _ 4^(r  ) + λ/8, L _ 5^(r  ) + λ/4, L _ 6^(r  ) + (3 λ)/32, L _ 8^(r  ), k _ 1^(r  ) + ((12 (C^(  ))^2)/(f _ π^(ó    ))^8 + (3 C^(  ))/(f _ π^(ó    ))^4 + 3/2) λ, k _ 2^(r  ) + (2 C^(  ) λ)/(f _ π^(ó    ))^4, k _ 3^(r  ) - (3 λ)/4, k _ 4^(r  ) - (2 C^(  ) λ)/(f _ π^(ó    ))^4, k _ 7^(r  ) + ((2 C^(  ))/(f _ π^(ó    ))^4 + 1/4) λ, k _ 8^(r  ) + λ/8, k _ 9^(r  ) + (-(12 (C^(  ))^2)/(5 (f _ π^(ó    ))^8) - (3 C^(  ))/(5 (f _ π^(ó    ))^4) - 3) λ, k _ 10^(r  ) + (-C^(  )/(5 (f _ π^(ó    ))^4) - 27/20) λ, k _ 11^(r  ) + (-C^(  )/(5 (f _ π^(ó    ))^4) - 1/4) λ, k _ 14^(r  )}

The contributions:

ff2 = amp2 /. Momentum[p2] -> -Momentum[p1] /. pirule // Simplify

(2 C^(  ) (δ _ (3 i _ 1)^(2) - 1) (e^(  ))^2 + (f _ π^(ó    ))^2 (p _ 1^2 - (m _ π^0^(ó    ))^2))/(f _ π^(ó    ))^2

amploop = ((Plus @@ ampinfinities /. Momentum[p3] -> -Momentum[p1]) /. pirule /. subpar // SUNReduce[#, FullReduce -> True] & // Simplify) /. delrules // Simplify

1/(96 π^2 (f _ π^(ó    ))^4 ((m _ π^+^(ó    ))^2 - (m _ γ^(ó    ))^2)) (4 C^(  ) ((m _ π^+^(ó    ))^2 - (m _ γ^(ó    ))^2) (-2 (32 π^2 λ + log((m _ π^+^(ó    ))^2/μ^2)) (3 δ _ (3 i _ 1)^(2) - 4) (m _ π^+^(ó    ))^2 - (32 π^2 λ + log((m _ π^0^(ó    ))^2/μ^2)) (m _ π^0^(ó    ))^2 (δ _ (3 i _ 1)^(2) - 1)) (e^(  ))^2 + (f _ π^(ó    ))^2 (-2 (32 π^2 λ + log((m _ π^+^(ó    ))^2/μ^2)) (δ _ (3 i _ 1)^(2) + 1) (m _ π^+^(ó    ))^6 + 2 (3 (e^(  ))^2 (8 π^2 Overscript[J, _] _ ((m _ γ^(ó    ))^2 (m _ π^+^(ó    ))^2)(p _ 1^2) (δ _ (1 i _ 1)^(2) - 2 δ _ (2 i _ 1)^(2)) - (32 π^2 λ + log((m _ π^+^(ó    ))^2/μ^2)) (δ _ (3 i _ 1)^(2) - 1)) (f _ π^(ó    ))^2 - (32 π^2 λ + log((m _ π^+^(ó    ))^2/μ^2)) (m _ π^0^(ó    ))^2 (δ _ (3 i _ 1)^(2) - 2) + (32 π^2 λ + log((m _ π^+^(ó    ))^2/μ^2)) (m _ γ^(ó    ))^2 (δ _ (3 i _ 1)^(2) + 1)) (m _ π^+^(ó    ))^4 + ((32 π^2 λ + log((m _ π^0^(ó    ))^2/μ^2)) (4 δ _ (3 i _ 1)^(2) - 1) (m _ π^0^(ó    ))^4 + 2 (32 π^2 λ + log((m _ π^+^(ó    ))^2/μ^2)) (m _ γ^(ó    ))^2 (δ _ (3 i _ 1)^(2) - 2) (m _ π^0^(ó    ))^2 + 12 (e^(  ))^2 (f _ π^(ó    ))^2 (m _ γ^(ó    ))^2 (8 π^2 Overscript[J, _] _ ((m _ γ^(ó    ))^2 (m _ π^+^(ó    ))^2)(p _ 1^2) (δ _ (2 i _ 1)^(2) - 2 δ _ (1 i _ 1)^(2)) + (64 π^2 λ + 2 log((m _ γ^(ó    ))^2/μ^2) + 1) (δ _ (3 i _ 1)^(2) - 1))) (m _ π^+^(ó    ))^2 - (m _ γ^(ó    ))^2 ((32 π^2 λ + log((m _ π^0^(ó    ))^2/μ^2)) (4 δ _ (3 i _ 1)^(2) - 1) (m _ π^0^(ó    ))^4 + 6 (e^(  ))^2 (f _ π^(ó    ))^2 (m _ γ^(ó    ))^2 (32 π^2 δ _ (2 i _ 1)^(2) λ + 128 π^2 δ _ (3 i _ 1)^(2) λ - 128 π^2 λ - 4 log((m _ γ^(ó    ))^2/μ^2) - 24 π^2 Overscript[J, _] _ ((m _ γ^(ó    ))^2 (m _ π^+^(ó    ))^2)(p _ 1^2) δ _ (1 i _ 1)^(2) + (32 π^2 λ + log((m _ γ^(ó    ))^2/μ^2)) δ _ (1 i _ 1)^(2) + log((m _ γ^(ó    ))^2/μ^2) δ _ (2 i _ 1)^(2) + 4 log((m _ γ^(ó    ))^2/μ^2) δ _ (3 i _ 1)^(2) + 2 δ _ (3 i _ 1)^(2) - 2)) + 2 p _ 1^2 (-(32 π^2 λ + log((m _ π^+^(ó    ))^2/μ^2)) (δ _ (3 i _ 1)^(2) + 1) (m _ π^+^(ó    ))^4 + (6 (e^(  ))^2 (16 π^2 Overscript[J, _] _ ((m _ γ^(ó    ))^2 (m _ π^+^(ó    ))^2)(p _ 1^2) - 32 π^2 λ - log((m _ π^+^(ó    ))^2/μ^2)) (δ _ (3 i _ 1)^(2) - 1) (f _ π^(ó    ))^2 + (32 π^2 λ + log((m _ π^0^(ó    ))^2/μ^2)) (m _ π^0^(ó    ))^2 (δ _ (3 i _ 1)^(2) - 1) + (32 π^2 λ + log((m _ π^+^(ó    ))^2/μ^2)) (m _ γ^(ó    ))^2 (δ _ (3 i _ 1)^(2) + 1)) (m _ π^+^(ó    ))^2 - (6 (e^(  ))^2 (16 π^2 Overscript[J, _] _ ((m _ γ^(ó    ))^2 (m _ π^+^(ó    ))^2)(p _ 1^2) - 32 π^2 λ - log((m _ γ^(ó    ))^2/μ^2)) (f _ π^(ó    ))^2 + (32 π^2 λ + log((m _ π^0^(ó    ))^2/μ^2)) (m _ π^0^(ó    ))^2) (m _ γ^(ó    ))^2 (δ _ (3 i _ 1)^(2) - 1))))

ampwf4 = amp4 /. Momentum[p2] -> -Momentum[p1] /. pirule // Renormalize // Simplify

1/(9 (f _ π^(ó    ))^2) (2 (((10 k _ 1^(r  ) (f _ π^(ó    ))^8 + 5 k _ 9^(r  ) (f _ π^(ó    ))^8 + 27 C^(  ) ((f _ π^(ó    ))^4 + 4 C^(  )) λ) (δ _ (3 i _ 1)^(2) - 1) (e^(  ))^4)/(f _ π^(ó    ))^4 + (f _ π^(ó    ))^2 (-2 (23 k _ 7^(r  ) + 18 k _ 8^(r  ) + 5 k _ 11^(r  ) + k _ 14^(r  ) + (45 C^(  ) λ)/(f _ π^(ó    ))^4 + (27 λ)/4) (m _ π^0^(ó    ))^2 + 1/2 (20 k _ 2^(r  ) + 20 k _ 10^(r  ) + 9 ((4 C^(  ))/(f _ π^(ó    ))^4 - 3) λ) p _ 1^2 + 9 (4 (k _ 7^(r  ) + k _ 8^(r  ) + (2 C^(  ) λ)/(f _ π^(ó    ))^4 + (3 λ)/8) (m _ π^0^(ó    ))^2 - (2 k _ 3^(r  ) + k _ 4^(r  ) - (2 C^(  ) λ)/(f _ π^(ó    ))^4 - (3 λ)/2) p _ 1^2) δ _ (3 i _ 1)^(2)) (e^(  ))^2 - 72 (2 L _ 6^(r  ) + L _ 8^(r  ) + (3 λ)/16) (m _ π^0^(ó    ))^4 + 18 (4 L _ 4^(r  ) + 2 L _ 5^(r  ) + λ) p _ 1^2 (m _ π^0^(ó    ))^2))

coeff1 = (Coefficient[amploop /. Pair[Momentum[p1], Momentum[p1]] -> ParticleMass[Pion, SUNIndex[i1], RenormalizationState[0]]^2, LeutwylerLambda[]] // Simplify) /. delrules // Simplify

-1/(3 (f _ π^(ó    ))^4) (4 C^(  ) ((6 δ _ (3 i _ 1)^(2) - 8) (m _ π^+^(ó    ))^2 + (m _ π^0^(ó    ))^2 (δ _ (3 i _ 1)^(2) - 1)) (e^(  ))^2 + (f _ π^(ó    ))^2 (2 (δ _ (3 i _ 1)^(2) + 1) (m _ π^+^(ó    ))^4 + 2 ((δ _ (3 i _ 1)^(2) - 2) (m _ π^0^(ó    ))^2 + (m _ π^(i _ 1   ))^2 (δ _ (3 i _ 1)^(2) + 1)) (m _ π^+^(ó    ))^2 + (m _ π^0^(ó    ))^2 ((m _ π^0^(ó    ))^2 (1 - 4 δ _ (3 i _ 1)^(2)) - 2 (m _ π^(i _ 1   ))^2 (δ _ (3 i _ 1)^(2) - 1)) + 6 (e^(  ))^2 (f _ π^(ó    ))^2 ((m _ π^+^(ó    ))^2 - 3 (m _ γ^(ó    ))^2 + 2 (m _ π^(i _ 1   ))^2) (δ _ (3 i _ 1)^(2) - 1)))

coeff2 = Coefficient[ampwf4 /. Pair[Momentum[p1], Momentum[p1]] -> ParticleMass[Pion, SUNIndex[i1], RenormalizationState[0]]^2, LeutwylerLambda[]] // Simplify

1/(f _ π^(ó    ))^6 (24 (C^(  ))^2 (δ _ (3 i _ 1)^(2) - 1) (e^(  ))^4 + 2 C^(  ) (f _ π^(ó    ))^2 (3 (e^(  ))^2 (δ _ (3 i _ 1)^(2) - 1) (f _ π^(ó    ))^2 + 2 ((4 δ _ (3 i _ 1)^(2) - 5) (m _ π^0^(ó    ))^2 + (m _ π^(i _ 1   ))^2 (δ _ (3 i _ 1)^(2) + 1))) (e^(  ))^2 + (f _ π^(ó    ))^4 (-3 (m _ π^0^(ó    ))^4 + 4 (m _ π^(i _ 1   ))^2 (m _ π^0^(ó    ))^2 + 3 (e^(  ))^2 (f _ π^(ó    ))^2 ((m _ π^0^(ó    ))^2 + (m _ π^(i _ 1   ))^2) (δ _ (3 i _ 1)^(2) - 1)))

coeff1 + coeff2 /. pimassrule // SUNReduce[#, FullReduce -> True] & // Simplify

1/(3 (f _ π^(ó    ))^6) (72 (C^(  ))^2 (δ _ (3 i _ 1)^(2) - 1) (e^(  ))^4 + 2 C^(  ) (f _ π^(ó    ))^2 (9 (e^(  ))^2 (δ _ (3 i _ 1)^(2) - 1) (f _ π^(ó    ))^2 + 2 ((3 δ _ (1 i _ 1)^(2) + 3 δ _ (2 i _ 1)^(2) - 6 δ _ (3 i _ 1)^(2) + 8) (m _ π^+^(ó    ))^2 + (m _ π^0^(ó    ))^2 (17 δ _ (3 i _ 1)^(2) - 14))) (e^(  ))^2 + (f _ π^(ó    ))^4 (3 (e^(  ))^2 (f _ π^(ó    ))^2 ((δ _ (1 i _ 1)^(2) + δ _ (2 i _ 1)^(2) - 2 δ _ (3 i _ 1)^(2) + 2) (m _ π^+^(ó    ))^2 + 3 ((m _ π^0^(ó    ))^2 + 2 (m _ γ^(ó    ))^2) (δ _ (3 i _ 1)^(2) - 1)) - 2 ((δ _ (1 i _ 1)^(2) + δ _ (2 i _ 1)^(2) + δ _ (3 i _ 1)^(2) + 1) (m _ π^+^(ó    ))^4 - (m _ π^0^(ó    ))^2 (5 δ _ (1 i _ 1)^(2) + 5 δ _ (2 i _ 1)^(2) - 3 δ _ (3 i _ 1)^(2) + 2) (m _ π^+^(ó    ))^2 + (m _ π^0^(ó    ))^4 (5 - 8 δ _ (3 i _ 1)^(2)))))

% /. dmrules /. delrules // Simplify

(3 (f _ π^(ó    ))^2 ((m _ π^+^(ó    ))^2 - (m _ π^0^(ó    ))^2) (m _ γ^(ó    ))^2 (δ _ (3 i _ 1)^(2) - 1))/C^(  )

The full amplitude (to fourth order) ff4 differs from the lowest order amplitude ff2 by a factor Z, ff4 = Z^(-1) ff2.  This is equivalent to a redefinition of the pion field, π _ r= Z^(-1/2)π.

ff4 = ff2 + amploop + ampwf4 // Simplify ;

We demand that ff4 be zero on the mass shell with  p^2=m _ (π, r)^2, where  m _ (π, r)^2=  m _ π^2+Cm is the renormalized mass. Since we are only working to O(p^4), we only need Cm to first order in m _ π^2.

Cm = Collect[(-ff4 /. dmrules /. Pair[Momentum[p1], Momentum[p1]] -> ParticleMass[Pion, SUNIndex[i1], RenormalizationState[0]]^2 /. pimassrule), LeutwylerLambda[]] //. delrules // SUNReduce[#, FullReduce -> True] & // Simplify

(-80 π^2 (2 k _ 1^(r  ) + k _ 9^(r  )) ((m _ π^+^(ó    ))^2 - (m _ π^0^(ó    ))^2)^2 ((m _ π^+^(ó    ))^2 - (m _ γ^(ó    ))^2) (δ _ (3 i _ 1)^(2) - 1) (f _ π^(ó    ))^8 + C^(  ) ((m _ π^+^(ó    ))^2 - (m _ π^0^(ó    ))^2) (288 π^2 Overscript[J, _] _ ((m _ γ^(ó    ))^2 (m _ π^+^(ó    ))^2)((m _ π^0^(ó    ))^2 δ _ (3 i _ 1)^(2) - (m _ π^+^(ó    ))^2 (δ _ (3 i _ 1)^(2) - 1)) (m _ π^+^(ó    ))^4 - 27 log((m _ π^+^(ó    ))^2/μ^2) (m _ π^+^(ó    ))^4 - 72 π^2 Overscript[J, _] _ ((m _ γ^(ó    ))^2 (m _ π^+^(ó    ))^2)((m _ π^0^(ó    ))^2 δ _ (3 i _ 1)^(2) - (m _ π^+^(ó    ))^2 (δ _ (3 i _ 1)^(2) - 1)) δ _ (1 i _ 1)^(2) (m _ π^+^(ó    ))^4 + 144 π^2 Overscript[J, _] _ ((m _ γ^(ó    ))^2 (m _ π^+^(ó    ))^2)((m _ π^0^(ó    ))^2 δ _ (3 i _ 1)^(2) - (m _ π^+^(ó    ))^2 (δ _ (3 i _ 1)^(2) - 1)) δ _ (2 i _ 1)^(2) (m _ π^+^(ó    ))^4 - 288 π^2 Overscript[J, _] _ ((m _ γ^(ó    ))^2 (m _ π^+^(ó    ))^2)((m _ π^0^(ó    ))^2 δ _ (3 i _ 1)^(2) - (m _ π^+^(ó    ))^2 (δ _ (3 i _ 1)^(2) - 1)) δ _ (3 i _ 1)^(2) (m _ π^+^(ó    ))^4 + 27 log((m _ π^+^(ó    ))^2/μ^2) δ _ (3 i _ 1)^(2) (m _ π^+^(ó    ))^4 + 1472 π^2 k _ 7^(r  ) (m _ π^0^(ó    ))^2 (m _ π^+^(ó    ))^2 + 1152 π^2 k _ 8^(r  ) (m _ π^0^(ó    ))^2 (m _ π^+^(ó    ))^2 + 320 π^2 k _ 11^(r  ) (m _ π^0^(ó    ))^2 (m _ π^+^(ó    ))^2 + 64 π^2 k _ 14^(r  ) (m _ π^0^(ó    ))^2 (m _ π^+^(ó    ))^2 - 288 π^2 Overscript[J, _] _ ((m _ γ^(ó    ))^2 (m _ π^+^(ó    ))^2)((m _ π^0^(ó    ))^2 δ _ (3 i _ 1)^(2) - (m _ π^+^(ó    ))^2 (δ _ (3 i _ 1)^(2) - 1)) (m _ γ^(ó    ))^2 (m _ π^+^(ó    ))^2 + 864 π^2 λ (m _ γ^(ó    ))^2 (m _ π^+^(ó    ))^2 + 54 log((m _ γ^(ó    ))^2/μ^2) (m _ γ^(ó    ))^2 (m _ π^+^(ó    ))^2 + 18 (m _ γ^(ó    ))^2 (m _ π^+^(ó    ))^2 + 288 π^2 Overscript[J, _] _ ((m _ γ^(ó    ))^2 (m _ π^+^(ó    ))^2)((m _ π^0^(ó    ))^2 δ _ (3 i _ 1)^(2) - (m _ π^+^(ó    ))^2 (δ _ (3 i _ 1)^(2) - 1)) (m _ γ^(ó    ))^2 δ _ (1 i _ 1)^(2) (m _ π^+^(ó    ))^2 - 144 π^2 Overscript[J, _] _ ((m _ γ^(ó    ))^2 (m _ π^+^(ó    ))^2)((m _ π^0^(ó    ))^2 δ _ (3 i _ 1)^(2) - (m _ π^+^(ó    ))^2 (δ _ (3 i _ 1)^(2) - 1)) (m _ γ^(ó    ))^2 δ _ (2 i _ 1)^(2) (m _ π^+^(ó    ))^2 + 576 π^2 k _ 3^(r  ) (m _ π^0^(ó    ))^2 δ _ (3 i _ 1)^(2) (m _ π^+^(ó    ))^2 + 288 π^2 k _ 4^(r  ) (m _ π^0^(ó    ))^2 δ _ (3 i _ 1)^(2) (m _ π^+^(ó    ))^2 - 1152 π^2 k _ 7^(r  ) (m _ π^0^(ó    ))^2 δ _ (3 i _ 1)^(2) (m _ π^+^(ó    ))^2 - 1152 π^2 k _ 8^(r  ) (m _ π^0^(ó    ))^2 δ _ (3 i _ 1)^(2) (m _ π^+^(ó    ))^2 + 288 π^2 Overscript[J, _] _ ((m _ γ^(ó    ))^2 (m _ π^+^(ó    ))^2)((m _ π^0^(ó    ))^2 δ _ (3 i _ 1)^(2) - (m _ π^+^(ó    ))^2 (δ _ (3 i _ 1)^(2) - 1)) (m _ γ^(ó    ))^2 δ _ (3 i _ 1)^(2) (m _ π^+^(ó    ))^2 - 864 π^2 λ (m _ γ^(ó    ))^2 δ _ (3 i _ 1)^(2) (m _ π^+^(ó    ))^2 - 54 log((m _ γ^(ó    ))^2/μ^2) (m _ γ^(ó    ))^2 δ _ (3 i _ 1)^(2) (m _ π^+^(ó    ))^2 - 18 (m _ γ^(ó    ))^2 δ _ (3 i _ 1)^(2) (m _ π^+^(ó    ))^2 - 864 π^2 λ (m _ γ^(ó    ))^4 - 27 log((m _ γ^(ó    ))^2/μ^2) (m _ γ^(ó    ))^4 - 18 (m _ γ^(ó    ))^4 - 1472 π^2 k _ 7^(r  ) (m _ π^0^(ó    ))^2 (m _ γ^(ó    ))^2 - 1152 π^2 k _ 8^(r  ) (m _ π^0^(ó    ))^2 (m _ γ^(ó    ))^2 - 320 π^2 k _ 11^(r  ) (m _ π^0^(ó    ))^2 (m _ γ^(ó    ))^2 - 64 π^2 k _ 14^(r  ) (m _ π^0^(ó    ))^2 (m _ γ^(ó    ))^2 - 216 π^2 Overscript[J, _] _ ((m _ γ^(ó    ))^2 (m _ π^+^(ó    ))^2)((m _ π^0^(ó    ))^2 δ _ (3 i _ 1)^(2) - (m _ π^+^(ó    ))^2 (δ _ (3 i _ 1)^(2) - 1)) (m _ γ^(ó    ))^4 δ _ (1 i _ 1)^(2) + 864 π^2 λ (m _ γ^(ó    ))^4 δ _ (3 i _ 1)^(2) + 27 log((m _ γ^(ó    ))^2/μ^2) (m _ γ^(ó    ))^4 δ _ (3 i _ 1)^(2) + 18 (m _ γ^(ó    ))^4 δ _ (3 i _ 1)^(2) - 576 π^2 k _ 3^(r  ) (m _ π^0^(ó    ))^2 (m _ γ^(ó    ))^2 δ _ (3 i _ 1)^(2) - 288 π^2 k _ 4^(r  ) (m _ π^0^(ó    ))^2 (m _ γ^(ó    ))^2 δ _ (3 i _ 1)^(2) + 1152 π^2 k _ 7^(r  ) (m _ π^0^(ó    ))^2 (m _ γ^(ó    ))^2 δ _ (3 i _ 1)^(2) + 1152 π^2 k _ 8^(r  ) (m _ π^0^(ó    ))^2 (m _ γ^(ó    ))^2 δ _ (3 i _ 1)^(2) + 320 π^2 k _ 2^(r  ) ((m _ π^+^(ó    ))^2 - (m _ γ^(ó    ))^2) ((m _ π^+^(ó    ))^2 (δ _ (3 i _ 1)^(2) - 1) - (m _ π^0^(ó    ))^2 δ _ (3 i _ 1)^(2)) + 320 π^2 k _ 10^(r  ) ((m _ π^+^(ó    ))^2 - (m _ γ^(ó    ))^2) ((m _ π^+^(ó    ))^2 (δ _ (3 i _ 1)^(2) - 1) - (m _ π^0^(ó    ))^2 δ _ (3 i _ 1)^(2))) (f _ π^(ó    ))^4 + 9 (C^(  ))^2 ((m _ π^+^(ó    ))^2 - (m _ γ^(ó    ))^2) (4 log((m _ π^+^(ó    ))^2/μ^2) (δ _ (3 i _ 1)^(2) - 1) (m _ π^+^(ó    ))^4 + 2 (m _ π^0^(ó    ))^2 (-log((m _ π^+^(ó    ))^2/μ^2) (δ _ (3 i _ 1)^(2) - 2) + 256 π^2 L _ 4^(r  ) (δ _ (3 i _ 1)^(2) - 1) + 128 π^2 L _ 5^(r  ) (δ _ (3 i _ 1)^(2) - 1)) (m _ π^+^(ó    ))^2 + (m _ π^0^(ó    ))^4 (1024 π^2 L _ 6^(r  ) + 512 π^2 L _ 8^(r  ) + log((m _ π^0^(ó    ))^2/μ^2) - 512 π^2 L _ 4^(r  ) δ _ (3 i _ 1)^(2) - 256 π^2 L _ 5^(r  ) δ _ (3 i _ 1)^(2) - 2 log((m _ π^0^(ó    ))^2/μ^2) δ _ (3 i _ 1)^(2))))/(288 π^2 (C^(  ))^2 (f _ π^(ó    ))^2 ((m _ π^+^(ó    ))^2 - (m _ γ^(ó    ))^2))

Coefficient[Cm, LeutwylerLambda[]] // FullSimplify

-(3 (f _ π^(ó    ))^2 ((m _ π^+^(ó    ))^2 - (m _ π^0^(ó    ))^2) (m _ γ^(ó    ))^2 (δ _ (3 i _ 1)^(2) - 1))/C^(  )

Here follow then the mass renormalization of the charged and neutral pions in terms of scale independent coupling constants:

CmPlus = (Limit[Cm /. i1 -> 1, ParticleMass[Vector[1], RenormalizationState[0]] -> 0] /. CouplingConstant[c_[4], n_, r___] :> If[RenormalizationCoefficients[c[4]][[n]] =!= 0, RenormalizationCoefficients[c[4]][[n]]/(32 Pi^2) (CouplingConstant[c[4], n, r] + Log[ParticleMass[PionZero, RenormalizationState[0]]^2/ScaleMu^2]), 1/(16 Pi^2) (CouplingConstant[c[4], n, r])] // Simplify) /. cancelScales // FullSimplify

1/(576 π^2 (C^(  ))^2 (f _ π^(ó    ))^2) (15 (k _ 1^(r  ) - k _ 9^(r  )) ((m _ π^+^(ó    ))^2 - (m _ π^0^(ó    ))^2)^2 (f _ π^(ó    ))^8 + C^(  ) ((m _ π^+^(ó    ))^2 - (m _ π^0^(ó    ))^2) (3 (10 k _ 1^(r  ) - k _ 9^(r  ) + 9 (k _ 10^(r  ) + 16 π^2 Overscript[J, _] _ (0 (m _ π^+^(ó    ))^2)((m _ π^+^(ó    ))^2) - 2 log((m _ π^+^(ó    ))^2/(m _ π^0^(ó    ))^2))) (m _ π^+^(ó    ))^2 + (-30 k _ 1^(r  ) + 23 k _ 7^(r  ) + 9 k _ 8^(r  ) + 3 k _ 9^(r  ) - 5 k _ 11^(r  ) + 8 k _ 14^(r  )) (m _ π^0^(ó    ))^2) (f _ π^(ó    ))^4 + 2 (C^(  ))^2 (2 (30 k _ 1^(r  ) - 10 k _ 2^(r  ) - 3 k _ 9^(r  ) + k _ 10^(r  ) - 18 log((m _ π^+^(ó    ))^2/(m _ π^0^(ó    ))^2)) (m _ π^+^(ó    ))^4 - 2 (9 L _ 4^(r  ) + 9 L _ 5^(r  ) + 60 k _ 1^(r  ) - 10 k _ 2^(r  ) - 46 k _ 7^(r  ) - 6 k _ 9^(r  ) + k _ 10^(r  ) + k _ 11^(r  ) - 18 log((m _ π^+^(ó    ))^2/(m _ π^0^(ó    ))^2)) (m _ π^0^(ó    ))^2 (m _ π^+^(ó    ))^2 + (27 L _ 6^(r  ) + 288 L _ 8^(r  ) + 60 k _ 1^(r  ) - 92 k _ 7^(r  ) - 6 k _ 9^(r  ) + 2 k _ 11^(r  )) (m _ π^0^(ó    ))^4))

Cm0 = (Limit[Cm /. i1 -> 3 /. CouplingConstant[c_[4], n_, r___] :> If[RenormalizationCoefficients[c[4]][[n]] =!= 0, RenormalizationCoefficients[c[4]][[n]]/(32 Pi^2) (CouplingConstant[c[4], n, r] + Log[ParticleMass[PionZero, RenormalizationState[0]]^2/ScaleMu^2]), 1/(16 Pi^2) (CouplingConstant[c[4], n, r])], ParticleMass[Vector[1], RenormalizationState[0]] -> 0] // Simplify) /. cancelScales // FullSimplify

1/(576 π^2 C^(  ) (f _ π^(ó    ))^2) ((m _ π^0^(ó    ))^2 ((27 k _ 3^(r  ) - 5 k _ 7^(r  ) - 27 k _ 10^(r  ) + 5 k _ 11^(r  ) - 8 k _ 14^(r  )) ((m _ π^0^(ó    ))^2 - (m _ π^+^(ó    ))^2) (f _ π^(ó    ))^4 + 2 C^(  ) ((-18 L _ 4^(r  ) - 18 L _ 5^(r  ) + 27 L _ 6^(r  ) + 288 L _ 8^(r  ) + 20 k _ 2^(r  ) + 18 k _ 4^(r  ) - 20 k _ 7^(r  ) - 2 k _ 10^(r  ) + 2 k _ 11^(r  )) (m _ π^0^(ó    ))^2 - 2 (10 k _ 2^(r  ) + 9 k _ 4^(r  ) - 10 k _ 7^(r  ) - k _ 10^(r  ) + k _ 11^(r  ) - 9 log((m _ π^+^(ó    ))^2/(m _ π^0^(ó    ))^2)) (m _ π^+^(ó    ))^2)))

Change variables to compare the neutral mass renormalization with Meissner, Müller and Steininger:

Cm0Final = (Simplify /@ Collect[Expand[Cm0], CouplingConstant[_[4], __]]) /. CouplingConstant[ChPTEM2[2], RenormalizationState[0]] -> Z * DecayConstant[Pion, RenormalizationState[0]]^4 /. dmrulesinv // Simplify

1/(288 π^2 Z (f _ π^(ó    ))^4) ((m _ π^0^(ó    ))^2 (C^(  ) (e^(  ))^2 (-40 Z k _ 2^(r  ) - 27 k _ 3^(r  ) - 36 Z k _ 4^(r  ) + 40 Z k _ 7^(r  ) + 5 k _ 7^(r  ) + 4 Z k _ 10^(r  ) + 27 k _ 10^(r  ) - 4 Z k _ 11^(r  ) - 5 k _ 11^(r  ) + 8 k _ 14^(r  )) - 9 Z (f _ π^(ó    ))^2 ((2 L _ 4^(r  ) + 2 L _ 5^(r  ) - 3 L _ 6^(r  ) - 32 L _ 8^(r  )) (m _ π^0^(ó    ))^2 - 2 log((m _ π^+^(ó    ))^2/(m _ π^0^(ó    ))^2) (m _ π^+^(ó    ))^2)))

The factor Z.

Adding the mass shift to the full amputated two-point function, (p _ 1^2 - (m _ π^+^(ó    ))^2)^2 times the two-point function (compare Urech's thesis formula 1.55) gives (1-Z)(p _ 1^2 - (m _ π^+^(ó    ))^2). Thus Z is found by dividing off (p _ 1^2 - (m _ π^+^(ó    ))^2).

zz = (((ff4 + Cm) /. i1 -> I1 /. dmrules // SUNReduce[#, FullReduce -> True] &)) /. delrules ;

zzz = (zz/ff2) /. i1 -> I1 // Simplify ;

Check that the (m _ γ^(ó    ))^2 Overscript[J, _] _ ((m _ γ^(ó    ))^2 (m _ π^+^(ó    ))^2) terms go:

Limit[pvs LeutwylerJBar[Pair[Momentum[p1], Momentum[p1]], pvs, ParticleMass[PseudoScalar[3], RenormalizationState[0]]^2, ExplicitLeutwylerJ0 -> True, Dimension -> D, MassScale -> ScaleMu, ExplicitLeutwylerLambda -> True, ExplicitLeutwylerSigma -> True, ExpandGammas -> True, DimensionExpand -> False, FixPoint -> 0, TaylorOrder -> 2, B0Evaluation -> Direct1, FCIntegrate -> IntegrateHeld], pvs -> 0]

0

z = Limit[zzz /. LeutwylerJBar[a___, ParticleMass[Vector[1], RenormalizationState[0]]^2, b___] -> LeutwylerJBar[a, pv2, b] /. ParticleMass[Vector[1], RenormalizationState[0]]^n_ -> pvs^(n/2), pvs -> 0] /. LeutwylerJBar[a___, pv2, b___] -> LeutwylerJBar[a, ParticleMass[Vector[1], RenormalizationState[0]]^2, b] // Simplify

(((m _ π^+^(ó    ))^2 - (m _ π^0^(ó    ))^2) (160 π^2 k _ 2^(r  ) ((δ _ (3 I _ 1)^(2) - 1) (m _ π^+^(ó    ))^2 + p _ 1^2 - (m _ π^0^(ó    ))^2 δ _ (3 I _ 1)^(2)) + 160 π^2 k _ 10^(r  ) ((δ _ (3 I _ 1)^(2) - 1) (m _ π^+^(ó    ))^2 + p _ 1^2 - (m _ π^0^(ó    ))^2 δ _ (3 I _ 1)^(2)) - 9 (-16 π^2 Overscript[J, _] _ ((m _ γ^(ó    ))^2 (m _ π^+^(ó    ))^2)((m _ π^0^(ó    ))^2 δ _ (3 I _ 1)^(2) - (m _ π^+^(ó    ))^2 (δ _ (3 I _ 1)^(2) - 1)) (m _ π^+^(ó    ))^2 + 8 π^2 λ (m _ π^+^(ó    ))^2 + log((m _ π^+^(ó    ))^2/μ^2) (m _ π^+^(ó    ))^2 + 4 π^2 Overscript[J, _] _ ((m _ γ^(ó    ))^2 (m _ π^+^(ó    ))^2)((m _ π^0^(ó    ))^2 δ _ (3 I _ 1)^(2) - (m _ π^+^(ó    ))^2 (δ _ (3 I _ 1)^(2) - 1)) δ _ (1 I _ 1)^(2) (m _ π^+^(ó    ))^2 - 8 π^2 Overscript[J, _] _ ((m _ γ^(ó    ))^2 (m _ π^+^(ó    ))^2)((m _ π^0^(ó    ))^2 δ _ (3 I _ 1)^(2) - (m _ π^+^(ó    ))^2 (δ _ (3 I _ 1)^(2) - 1)) δ _ (2 I _ 1)^(2) (m _ π^+^(ó    ))^2 + 16 π^2 Overscript[J, _] _ ((m _ γ^(ó    ))^2 (m _ π^+^(ó    ))^2)((m _ π^0^(ó    ))^2 δ _ (3 I _ 1)^(2) - (m _ π^+^(ó    ))^2 (δ _ (3 I _ 1)^(2) - 1)) δ _ (3 I _ 1)^(2) (m _ π^+^(ó    ))^2 - 8 π^2 λ δ _ (3 I _ 1)^(2) (m _ π^+^(ó    ))^2 - log((m _ π^+^(ó    ))^2/μ^2) δ _ (3 I _ 1)^(2) (m _ π^+^(ó    ))^2 - 4 π^2 Overscript[J, _] _ ((m _ γ^(ó    ))^2 (m _ π^+^(ó    ))^2)(p _ 1^2) ((δ _ (1 I _ 1)^(2) - 2 δ _ (2 I _ 1)^(2)) (m _ π^+^(ó    ))^2 + 4 p _ 1^2 (δ _ (3 I _ 1)^(2) - 1)) - 32 π^2 k _ 3^(r  ) (m _ π^0^(ó    ))^2 δ _ (3 I _ 1)^(2) - 16 π^2 k _ 4^(r  ) (m _ π^0^(ó    ))^2 δ _ (3 I _ 1)^(2) + p _ 1^2 ((8 π^2 λ + log((m _ π^+^(ó    ))^2/μ^2)) (δ _ (3 I _ 1)^(2) - 1) + 32 π^2 k _ 3^(r  ) δ _ (3 I _ 1)^(2) + 16 π^2 k _ 4^(r  ) δ _ (3 I _ 1)^(2)))) (f _ π^(ó    ))^4 + 3 C^(  ) (-64 π^2 λ (m _ π^+^(ó    ))^4 + log((m _ π^+^(ó    ))^2/μ^2) (m _ π^+^(ó    ))^4 + 64 π^2 λ δ _ (3 I _ 1)^(2) (m _ π^+^(ó    ))^4 - log((m _ π^+^(ó    ))^2/μ^2) δ _ (3 I _ 1)^(2) (m _ π^+^(ó    ))^4 - 768 π^2 L _ 4^(r  ) (m _ π^0^(ó    ))^2 (m _ π^+^(ó    ))^2 - 384 π^2 L _ 5^(r  ) (m _ π^0^(ó    ))^2 (m _ π^+^(ó    ))^2 - 64 π^2 λ (m _ π^0^(ó    ))^2 (m _ π^+^(ó    ))^2 + log((m _ π^0^(ó    ))^2/μ^2) (m _ π^0^(ó    ))^2 (m _ π^+^(ó    ))^2 + 768 π^2 L _ 4^(r  ) (m _ π^0^(ó    ))^2 δ _ (3 I _ 1)^(2) (m _ π^+^(ó    ))^2 + 384 π^2 L _ 5^(r  ) (m _ π^0^(ó    ))^2 δ _ (3 I _ 1)^(2) (m _ π^+^(ó    ))^2 - 64 π^2 λ (m _ π^0^(ó    ))^2 δ _ (3 I _ 1)^(2) (m _ π^+^(ó    ))^2 + 2 log((m _ π^+^(ó    ))^2/μ^2) (m _ π^0^(ó    ))^2 δ _ (3 I _ 1)^(2) (m _ π^+^(ó    ))^2 - log((m _ π^0^(ó    ))^2/μ^2) (m _ π^0^(ó    ))^2 δ _ (3 I _ 1)^(2) (m _ π^+^(ó    ))^2 - 768 π^2 L _ 4^(r  ) (m _ π^0^(ó    ))^4 δ _ (3 I _ 1)^(2) - 384 π^2 L _ 5^(r  ) (m _ π^0^(ó    ))^4 δ _ (3 I _ 1)^(2) + 48 π^2 (f _ π^(ó    ))^2 ((δ _ (3 I _ 1)^(2) - 1) (m _ π^+^(ó    ))^2 + p _ 1^2 - (m _ π^0^(ó    ))^2 δ _ (3 I _ 1)^(2)) + p _ 1^2 ((64 π^2 λ - log((m _ π^+^(ó    ))^2/μ^2)) (δ _ (3 I _ 1)^(2) + 1) (m _ π^+^(ó    ))^2 + (m _ π^0^(ó    ))^2 (768 π^2 L _ 4^(r  ) + 384 π^2 L _ 5^(r  ) - (64 π^2 λ - log((m _ π^0^(ó    ))^2/μ^2)) (δ _ (3 I _ 1)^(2) - 1)))))/(144 π^2 C^(  ) (2 C^(  ) (δ _ (3 I _ 1)^(2) - 1) (e^(  ))^2 + (f _ π^(ó    ))^2 (p _ 1^2 - (m _ π^0^(ó    ))^2)))

The neutral factor:

zZero = z /. I1 -> 3 // Simplify

1/(72 π^2 C^(  ) (f _ π^(ó    ))^2) (8 π^2 (10 k _ 2^(r  ) - 18 k _ 3^(r  ) - 9 k _ 4^(r  ) + 10 k _ 10^(r  )) ((m _ π^+^(ó    ))^2 - (m _ π^0^(ó    ))^2) (f _ π^(ó    ))^4 + 3 C^(  ) (24 π^2 (f _ π^(ó    ))^2 + (64 π^2 λ - log((m _ π^+^(ó    ))^2/μ^2)) (m _ π^+^(ó    ))^2 + 192 π^2 (2 L _ 4^(r  ) + L _ 5^(r  )) (m _ π^0^(ó    ))^2))

The charged amputated two-point function with the mass shift subtracted off:

zzPlus = zz /. I1 -> 1 // Simplify

(((m _ π^+^(ó    ))^2 - (m _ π^0^(ó    ))^2) (160 π^2 k _ 2^(r  ) (p _ 1^2 - (m _ π^+^(ó    ))^2) ((m _ π^+^(ó    ))^2 - (m _ γ^(ó    ))^2) + 160 π^2 k _ 10^(r  ) (p _ 1^2 - (m _ π^+^(ó    ))^2) ((m _ π^+^(ó    ))^2 - (m _ γ^(ó    ))^2) - 9 (-12 π^2 Overscript[J, _] _ ((m _ γ^(ó    ))^2 (m _ π^+^(ó    ))^2)((m _ π^+^(ó    ))^2) (m _ π^+^(ó    ))^4 + 8 π^2 λ (m _ π^+^(ó    ))^4 + log((m _ π^+^(ó    ))^2/μ^2) (m _ π^+^(ó    ))^4 - 8 π^2 λ (m _ γ^(ó    ))^2 (m _ π^+^(ó    ))^2 - log((m _ γ^(ó    ))^2/μ^2) (m _ γ^(ó    ))^2 (m _ π^+^(ó    ))^2 + 12 π^2 Overscript[J, _] _ ((m _ γ^(ó    ))^2 (m _ π^+^(ó    ))^2)((m _ π^+^(ó    ))^2) (m _ γ^(ó    ))^4 - 4 π^2 Overscript[J, _] _ ((m _ γ^(ó    ))^2 (m _ π^+^(ó    ))^2)(p _ 1^2) ((m _ π^+^(ó    ))^2 - 3 (m _ γ^(ó    ))^2 - 4 p _ 1^2) ((m _ π^+^(ó    ))^2 - (m _ γ^(ó    ))^2) + p _ 1^2 ((8 π^2 λ + log((m _ γ^(ó    ))^2/μ^2)) (m _ γ^(ó    ))^2 - (8 π^2 λ + log((m _ π^+^(ó    ))^2/μ^2)) (m _ π^+^(ó    ))^2))) (f _ π^(ó    ))^4 + 3 C^(  ) (p _ 1^2 - (m _ π^+^(ó    ))^2) (48 π^2 (f _ π^(ó    ))^2 + (64 π^2 λ - log((m _ π^+^(ó    ))^2/μ^2)) (m _ π^+^(ó    ))^2 + (768 π^2 L _ 4^(r  ) + 384 π^2 L _ 5^(r  ) + 64 π^2 λ - log((m _ π^0^(ó    ))^2/μ^2)) (m _ π^0^(ó    ))^2) ((m _ π^+^(ó    ))^2 - (m _ γ^(ó    ))^2))/(144 π^2 C^(  ) (f _ π^(ó    ))^2 ((m _ π^+^(ó    ))^2 - (m _ γ^(ó    ))^2))

The charged factor:

zPlusFull = z /. I1 -> 1 /. dmrules // Simplify

1/(144 π^2 C^(  ) (f _ π^(ó    ))^2 (p _ 1^2 - (m _ π^+^(ó    ))^2)) ((160 π^2 k _ 2^(r  ) (p _ 1^2 - (m _ π^+^(ó    ))^2) + 160 π^2 k _ 10^(r  ) (p _ 1^2 - (m _ π^+^(ó    ))^2) - 9 ((-12 π^2 Overscript[J, _] _ ((m _ γ^(ó    ))^2 (m _ π^+^(ó    ))^2)((m _ π^+^(ó    ))^2) + 8 π^2 λ + log((m _ π^+^(ó    ))^2/μ^2)) (m _ π^+^(ó    ))^2 - (8 π^2 λ + log((m _ π^+^(ó    ))^2/μ^2)) p _ 1^2 + 4 π^2 Overscript[J, _] _ ((m _ γ^(ó    ))^2 (m _ π^+^(ó    ))^2)(p _ 1^2) (4 p _ 1^2 - (m _ π^+^(ó    ))^2))) ((m _ π^+^(ó    ))^2 - (m _ π^0^(ó    ))^2) (f _ π^(ó    ))^4 + 3 C^(  ) (p _ 1^2 - (m _ π^+^(ó    ))^2) (48 π^2 (f _ π^(ó    ))^2 + (64 π^2 λ - log((m _ π^+^(ó    ))^2/μ^2)) (m _ π^+^(ó    ))^2 + (768 π^2 L _ 4^(r  ) + 384 π^2 L _ 5^(r  ) + 64 π^2 λ - log((m _ π^0^(ó    ))^2/μ^2)) (m _ π^0^(ó    ))^2))

Check that it has a finite limit:

LeutwylerJBar[ParticleMass[PseudoScalar[3], RenormalizationState[0]]^2, ParticleMass[Vector[1], RenormalizationState[0]]^2, ParticleMass[PseudoScalar[3], RenormalizationState[0]]^2, LeutwylerJBarEvaluation -> "subthreshold", ExplicitLeutwylerSigma -> True] // Simplify

1/(32 π^2 (m _ π^+^(ó    ))^2 ((m _ π^+^(ó    ))^2 - (m _ γ^(ó    ))^2)) (2 (m _ π^+^(ó    ))^4 + ((3 log((m _ π^+^(ó    ))^2/(m _ γ^(ó    ))^2) - 2) (m _ γ^(ó    ))^2 - log(((m _ γ^(ó    ))^2 - ((m _ γ^(ó    ))^4 - 4 (m _ π^+^(ó    ))^2 (m _ γ^(ó    ))^2)^(1/2))/((m _ γ^(ó    ))^2 + ((m _ γ^(ó    ))^4 - 4 (m _ π^+^(ó    ))^2 (m _ γ^(ó    ))^2)^(1/2))) ((m _ γ^(ó    ))^4 - 4 (m _ π^+^(ó    ))^2 (m _ γ^(ó    ))^2)^(1/2)) (m _ π^+^(ó    ))^2 - log((m _ π^+^(ó    ))^2/(m _ γ^(ó    ))^2) (m _ γ^(ó    ))^4 + log(((m _ γ^(ó    ))^2 - ((m _ γ^(ó    ))^4 - 4 (m _ π^+^(ó    ))^2 (m _ γ^(ó    ))^2)^(1/2))/((m _ γ^(ó    ))^2 + ((m _ γ^(ó    ))^4 - 4 (m _ π^+^(ó    ))^2 (m _ γ^(ó    ))^2)^(1/2))) (m _ γ^(ó    ))^2 ((m _ γ^(ó    ))^4 - 4 (m _ π^+^(ó    ))^2 (m _ γ^(ó    ))^2)^(1/2))

Limit[%, ParticleMass[Photon, RenormalizationState[0]] -> 0]

1/(16 π^2)

zPlus = Limit[zPlusFull /. LeutwylerJBar[Pair[Momentum[p1], Momentum[p1]] | ParticleMass[PseudoScalar[3], RenormalizationState[0]]^2, ParticleMass[Vector[1], RenormalizationState[0]]^2, ParticleMass[PseudoScalar[3], RenormalizationState[0]]^2, ___] -> 1/(16 π^2), Pair[Momentum[p1], Momentum[p1]] -> ParticleMass[PionPlus, RenormalizationState[0]]^2] // Simplify

1/(144 π^2 C^(  ) (f _ π^(ó    ))^2) ((160 π^2 k _ 2^(r  ) + 160 π^2 k _ 10^(r  ) + 9 (8 π^2 λ + log((m _ π^+^(ó    ))^2/μ^2) - 1)) ((m _ π^+^(ó    ))^2 - (m _ π^0^(ó    ))^2) (f _ π^(ó    ))^4 + 3 C^(  ) (48 π^2 (f _ π^(ó    ))^2 + (64 π^2 λ - log((m _ π^+^(ó    ))^2/μ^2)) (m _ π^+^(ó    ))^2 + (768 π^2 L _ 4^(r  ) + 384 π^2 L _ 5^(r  ) + 64 π^2 λ - log((m _ π^0^(ó    ))^2/μ^2)) (m _ π^0^(ó    ))^2))

Check the strong part:

zZero /. {PionPlus -> Pion, PionZero -> Pion} // Simplify

(24 π^2 (f _ π^(ó    ))^2 + (384 π^2 L _ 4^(r  ) + 192 π^2 L _ 5^(r  ) + 64 π^2 λ - log((m _ π^(ó    ))^2/μ^2)) (m _ π^(ó    ))^2)/(24 π^2 (f _ π^(ó    ))^2)

zPlus /. {PionPlus -> Pion, PionZero -> Pion} // Simplify

(24 π^2 (f _ π^(ó    ))^2 + (384 π^2 L _ 4^(r  ) + 192 π^2 L _ 5^(r  ) + 64 π^2 λ - log((m _ π^(ó    ))^2/μ^2)) (m _ π^(ó    ))^2)/(24 π^2 (f _ π^(ó    ))^2)

CheckF[dum, "ChPT2P20o2.Fac"] // Renormalize // Together // Simplify

(24 π^2 (f _ π^(ó    ))^2 + (384 π^2 L _ 4^(r  ) + 192 π^2 L _ 5^(r  ) + 64 π^2 λ - log((m _ π^(ó    ))^2/μ^2)) (m _ π^(ó    ))^2)/(24 π^2 (f _ π^(ó    ))^2)

Save for later use:

$VeryVerbose = 2 ;

CheckF[z, "ChPTEM2P20o2.Fac"] ;

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

File does not exist, evaluating

Saving

CheckF[zPlus, "ChPTEM2P30o2.Fac"] ;

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

File does not exist, evaluating

Saving

CheckF[zZero, "ChPTEM2P40o2.Fac"] ;

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

File does not exist, evaluating

Saving

$VeryVerbose = 0 ;


Converted by Mathematica  (July 10, 2003)