59 G4ClassificationOfNewTrack classification = fUrgent;
61 if (aTrack->GetParentID() == 0) {
65 G4int pol = G4PhononPolarization::Get(aTrack->GetDefinition());
68 G4ThreeVector Ran = G4RandomDirection();
72 theKmap->SetK(aTrack, Ran);
75 G4ThreeVector momentumDir = LM->MapKtoVDir(aTrack->GetVolume(), pol, Ran);
78 G4double velocity = LM->MapKtoV(aTrack->GetVolume(), pol, Ran);
83 theTrack->SetMomentumDirection(momentumDir);
84 theTrack->SetVelocity(velocity);
85 theTrack->UseGivenVelocity(
true);
88 return classification;