Method Update

source/track/src/G4FieldTrackUpdator.cc:52–85

---------------------------------------------------------------------

//---------------------------------------------------------------------
void G4FieldTrackUpdator::Update(G4FieldTrack* ftrk, const G4Track* trk)
{
  const G4DynamicParticle* ptDynamicParticle = trk->GetDynamicParticle();

  // The following properties must be updated 1) for each new track, and
  ftrk->SetRestMass(ptDynamicParticle->GetMass());
  // 2) Since ion can lose/gain electrons, this must be done at every step

  ftrk->UpdateState(trk->GetPosition(), trk->GetGlobalTime(),
                    trk->GetMomentumDirection(), trk->GetKineticEnergy());

#ifdef G4CHECK
  if((trk->GetMomentum() - ftrk->GetMomentum()).mag2() >
     1.e-16 * trk->GetMomentum().mag2())
  {
    G4cerr << "ERROR> G4FieldTrackUpdator sees *Disagreement* in momentum "
           << G4endl;
    G4cout << "  FTupdator: Tracking Momentum= " << trk->GetMomentum()
           << G4endl;
    G4cout << "  FTupdator: FldTrack Momentum= " << ftrk->GetMomentum()
           << G4endl;
    G4cout << "  FTupdator: FldTrack-Tracking= "
           << ftrk->GetMomentum() - trk->GetMomentum() << G4endl;
  }
#endif

  ftrk->SetProperTimeOfFlight(trk->GetProperTime());

  ftrk->SetChargeAndMoments(ptDynamicParticle->GetCharge(),
                            ptDynamicParticle->GetMagneticMoment());
  ftrk->SetPDGSpin(ptDynamicParticle->GetParticleDefinition()->GetPDGSpin());
  // The charge can change during tracking
  ftrk->SetSpin(ptDynamicParticle->GetPolarization());
}