•Preliminaries

udrules = {PionZero -> PionPlus}

{π^0 -> π^+}

subpar = Table[(ParticleMass[PseudoScalar[2], i | SUNIndex[i], r___] -> ParticleMass[Select[$IsoSpinProjectionRules, (! FreeQ[#, {i}] &)][[1]][[1]], r]), {i, 3}]

{ParticleMass(π, 1 | 1, r___) -> ParticleMass(π^+, r), ParticleMass(π, 2 | 2, r___) -> ParticleMass(π^+, r), ParticleMass(π, 3 | 3, r___) -> ParticleMass(π^0, r)}

{ParticleMass(π, 1, r___) -> ParticleMass (π^+, r), ParticleMass(π, 2, r___) -> ParticleMass (π^+, r), ParticleMass(π, 3, r___) -> ParticleMass (π^0, r)}

pimassrule = ParticleMass[Pion, i_ ? (! IntegerQ[#] &) | SUNIndex[i_ ? (! IntegerQ[#] &)], r___]^2 -> (SU2Delta[SUNIndex[i], SUNIndex[1]] + SU2Delta[SUNIndex[i], SUNIndex[2]]) ParticleMass[PionPlus, r]^2 + SU2Delta[SUNIndex[i], SUNIndex[3]] ParticleMass[PionZero, r]^2

ParticleMass(π, i_ ? (¬ IntegerQ[#1] &) | i_ ? (¬ IntegerQ[#1] &), r___)^2 -> (δ _ (1 i)^(2) + δ _ (2 i)^(2)) ParticleMass(π^+, r)^2 + ParticleMass(π^0, r)^2 δ _ (3 i)^(2)

pirule = ParticleMass[Pion, RenormalizationState[0]] -> ParticleMass[PionZero, RenormalizationState[0]]

m _ π^(ó  ) -> m _ π^0^(ó  )

delrules = {SU2Delta[SUNIndex[1], SUNIndex[i1]] + SU2Delta[SUNIndex[2], SUNIndex[i1]] -> (1 - SU2Delta[SUNIndex[3], SUNIndex[i1]]), a_ * SU2Delta[SUNIndex[1], SUNIndex[i1]] + a_ * SU2Delta[SUNIndex[2], SUNIndex[i1]] -> a * (1 - SU2Delta[SUNIndex[3], SUNIndex[i1]])}

{δ _ (1 i _ 1)^(2) + δ _ (2 i _ 1)^(2) -> 1 - δ _ (3 i _ 1)^(2), a_ δ _ (1 i _ 1)^(2) + a_ δ _ (2 i _ 1)^(2) -> a (1 - δ _ (3 i _ 1)^(2))}

cancelScales = Log[ParticleMass[PionPlus, RenormalizationState[0]]^2/ScaleMu^2] -> Log[ParticleMass[PionZero, RenormalizationState[0]]^2/ScaleMu^2] + Log[ParticleMass[PionPlus, RenormalizationState[0]]^2/ParticleMass[PionZero, RenormalizationState[0]]^2]

log((m _ π^+^(ó  ))^2/μ^2) -> log((m _ π^+^(ó  ))^2/(m _ π^0^(ó  ))^2) + log((m _ π^0^(ó  ))^2/μ^2)

dmrules = CouplingConstant[QED[1], RenormalizationState[0]]^n_ :> DecayConstant[PseudoScalar[2], RenormalizationState[0]]^n/(2 CouplingConstant[ChPTEM2[2], RenormalizationState[0]])^(n/2) (ParticleMass[PionPlus, RenormalizationState[0]]^2 - ParticleMass[PionZero, RenormalizationState[0]]^2)^(n/2)

(e^( ))^n_ :> ((f _ π^(ó  ))^n ((m _ PionPlus^(ó  ))^2 - (m _ PionZero^(ó  ))^2)^(n/2))/(2 C^( ))^(n/2)

dmrulesinv = {(ParticleMass[PionPlus, RenormalizationState[0]]^2 - ParticleMass[PionZero, RenormalizationState[0]]^2) -> 2 CouplingConstant[QED[1], RenormalizationState[0]]^2 CouplingConstant[ChPTEM2[2], RenormalizationState[0]]/DecayConstant[Pion, RenormalizationState[0]]^2, (ParticleMass[PionZero, RenormalizationState[0]]^2 - ParticleMass[PionPlus, RenormalizationState[0]]^2) -> -2 CouplingConstant[QED[1], RenormalizationState[0]]^2 CouplingConstant[ChPTEM2[2], RenormalizationState[0]]/DecayConstant[Pion, RenormalizationState[0]]^2}

{(m _ π^+^(ó  ))^2 - (m _ π^0^(ó  ))^2 -> (2 C^( ) (e^( ))^2)/(f _ π^(ó  ))^2, (m _ π^0^(ó  ))^2 - (m _ π^+^(ó  ))^2 -> -(2 C^( ) (e^( ))^2)/(f _ π^(ó  ))^2}

Converted by Mathematica  (July 10, 2003)