For heavy ions, note following from TrackingAction.cxx
// note that for ions, the stored pdg value: we recreate it from the baryon number
G4ParticleDefinition* def = aTrack->GetDefinition();
int pdgid = def->GetPDGEncoding();
if ( pdgid==0 ) pdgid = 80000+def->GetBaryonNumber();
