G4FastStep Class Reference

#include <G4FastStep.hh>

Inheritance diagram for G4FastStep:

Public Member Functions
	G4FastStep ()=default
	~G4FastStep () override=default
	G4FastStep (const G4FastStep &right)=delete
G4FastStep &	operator= (const G4FastStep &right)=delete
void	KillPrimaryTrack ()
void	ProposePrimaryTrackFinalPosition (const G4ThreeVector &, G4bool localCoordinates=true)
void	SetPrimaryTrackFinalPosition (const G4ThreeVector &, G4bool localCoordinates=true)
void	ProposePrimaryTrackFinalTime (G4double)
void	SetPrimaryTrackFinalTime (G4double)
void	ProposePrimaryTrackFinalProperTime (G4double)
void	SetPrimaryTrackFinalProperTime (G4double)
void	ProposePrimaryTrackFinalMomentumDirection (const G4ThreeVector &, G4bool localCoordinates=true)
void	SetPrimaryTrackFinalMomentum (const G4ThreeVector &, G4bool localCoordinates=true)
void	ProposePrimaryTrackFinalKineticEnergy (G4double)
void	SetPrimaryTrackFinalKineticEnergy (G4double)
void	ProposePrimaryTrackFinalKineticEnergyAndDirection (G4double, const G4ThreeVector &, G4bool localCoordinates=true)
void	SetPrimaryTrackFinalKineticEnergyAndDirection (G4double, const G4ThreeVector &, G4bool localCoordinates=true)
void	ProposePrimaryTrackFinalPolarization (const G4ThreeVector &, G4bool localCoordinates=true)
void	SetPrimaryTrackFinalPolarization (const G4ThreeVector &, G4bool localCoordinates=true)
void	ProposePrimaryTrackPathLength (G4double)
void	SetPrimaryTrackPathLength (G4double)
void	ProposePrimaryTrackFinalEventBiasingWeight (G4double)
void	SetPrimaryTrackFinalEventBiasingWeight (G4double)
void	SetNumberOfSecondaryTracks (G4int)
G4int	GetNumberOfSecondaryTracks ()
G4Track *	CreateSecondaryTrack (const G4DynamicParticle &, G4ThreeVector, G4ThreeVector, G4double, G4bool localCoordinates=true)
G4Track *	CreateSecondaryTrack (const G4DynamicParticle &, G4ThreeVector, G4double, G4bool localCoordinates=true)
G4Track *	GetSecondaryTrack (G4int)
void	ProposeTotalEnergyDeposited (G4double anEnergyPart)
void	SetTotalEnergyDeposited (G4double anEnergyPart)
G4double	GetTotalEnergyDeposited () const
void	ForceSteppingHitInvocation ()
G4Step *	UpdateStepForAtRest (G4Step *Step) override
G4Step *	UpdateStepForPostStep (G4Step *Step) override
void	Initialize (const G4FastTrack &)
void	DumpInfo () const override
G4bool	CheckIt (const G4Track &) override
Public Member Functions inherited from G4VParticleChange
	G4VParticleChange ()
virtual	~G4VParticleChange ()=default
	G4VParticleChange (const G4VParticleChange &right)=delete
G4VParticleChange &	operator= (const G4VParticleChange &right)=delete
virtual G4Step *	UpdateStepForAlongStep (G4Step *Step)
G4double	GetTrueStepLength () const
void	ProposeTrueStepLength (G4double truePathLength)
G4double	GetLocalEnergyDeposit () const
void	ProposeLocalEnergyDeposit (G4double anEnergyPart)
G4double	GetNonIonizingEnergyDeposit () const
void	ProposeNonIonizingEnergyDeposit (G4double anEnergyPart)
G4TrackStatus	GetTrackStatus () const
void	ProposeTrackStatus (G4TrackStatus status)
const G4Track *	GetCurrentTrack () const
G4SteppingControl	GetSteppingControl () const
void	ProposeSteppingControl (G4SteppingControl StepControlFlag)
G4bool	GetFirstStepInVolume () const
G4bool	GetLastStepInVolume () const
void	ProposeFirstStepInVolume (G4bool flag)
void	ProposeLastStepInVolume (G4bool flag)
void	Clear ()
void	SetNumberOfSecondaries (G4int totSecondaries)
G4int	GetNumberOfSecondaries () const
G4Track *	GetSecondary (G4int anIndex) const
void	AddSecondary (G4Track *aSecondary)
G4double	GetWeight () const
G4double	GetParentWeight () const
void	ProposeWeight (G4double finalWeight)
void	ProposeParentWeight (G4double finalWeight)
void	SetSecondaryWeightByProcess (G4bool)
G4bool	IsSecondaryWeightSetByProcess () const
void	SetParentWeightByProcess (G4bool)
G4bool	IsParentWeightSetByProcess () const
void	SetVerboseLevel (G4int vLevel)
G4int	GetVerboseLevel () const
void	ClearDebugFlag ()
void	SetDebugFlag ()
G4bool	GetDebugFlag () const

Additional Inherited Members
Protected Member Functions inherited from G4VParticleChange
G4Step *	UpdateStepInfo (G4Step *Step)
void	InitializeLocalEnergyDeposit ()
void	InitializeSteppingControl ()
void	InitializeParentWeight (const G4Track &)
void	InitializeStatusChange (const G4Track &)
void	InitializeSecondaries ()
void	InitializeFromStep (const G4Step *)
G4double	ComputeBeta (G4double kinEnergy)
G4bool	CheckSecondary (G4Track &)
G4double	GetAccuracyForWarning () const
G4double	GetAccuracyForException () const
Protected Attributes inherited from G4VParticleChange
const G4Track *	theCurrentTrack = nullptr
std::vector< G4Track * >	theListOfSecondaries
G4TrackStatus	theStatusChange = fAlive
G4SteppingControl	theSteppingControlFlag = NormalCondition
G4double	theLocalEnergyDeposit = 0.0
G4double	theNonIonizingEnergyDeposit = 0.0
G4double	theTrueStepLength = 0.0
G4double	theParentWeight = 1.0
G4double	theParentGlobalTime = 0.0
G4int	theNumberOfSecondaries = 0
G4int	theSizeOftheListOfSecondaries = 0
G4int	verboseLevel = 1
G4int	nError = 0
G4bool	theFirstStepInVolume = false
G4bool	theLastStepInVolume = false
G4bool	isParentWeightProposed = false
G4bool	fSetSecondaryWeightByProcess = false
G4bool	debugFlag = false
Static Protected Attributes inherited from G4VParticleChange
static const G4double	accuracyForWarning = 1.0e-9
static const G4double	accuracyForException = 0.001
static const G4int	maxError = 10

Detailed Description

Definition at line 89 of file G4FastStep.hh.

Constructor & Destructor Documentation

G4FastStep() [1/2]

G4FastStep::G4FastStep ( )

default

~G4FastStep()

G4FastStep::~G4FastStep ( )

overridedefault

G4FastStep() [2/2]

G4FastStep::G4FastStep ( const G4FastStep & right )

delete

Member Function Documentation

CheckIt()

G4bool G4FastStep::CheckIt ( const G4Track & aTrack )

overridevirtual

Reimplemented from G4VParticleChange.

Definition at line 359 of file G4FastStep.cc.

{
  //
  //      In the G4FastStep::CheckIt
  //      We only check a bit
  //
  //      If the user violates the energy,
  //      We don't care, we agree.
  //
  //      But for theMomentumDirectionChange,
  //      We do pay attention.
  //      And if too large is its range,
  //      We issue an Exception.
  //
  //
  // It means, the G4FastStep::CheckIt issues an exception only for the
  // theMomentumDirectionChange which should be an unit vector
  // and it corrects it because it could cause problems for the ulterior
  // tracking.For the rest, only warning are issued.
 
  G4bool itsOK = true;
  G4bool exitWithError = false;
  G4double accuracy;
 
  // Energy should not be larger than the initial value
  accuracy = (theEnergyChange - aTrack.GetKineticEnergy()) / MeV;
  if (accuracy > GetAccuracyForWarning()) {
    G4ExceptionDescription ed;
    ed << "The energy becomes larger than the initial value, difference = " << accuracy << " MeV"
       << G4endl;
    G4Exception("G4FastStep::CheckIt(const G4Track& aTrack)", "FastSim006", JustWarning, ed);
    itsOK = false;
    if (accuracy > GetAccuracyForException()) {
      exitWithError = true;
    }
  }
 
  G4bool itsOKforMomentum = true;
  if (theEnergyChange > 0.) {
    accuracy = std::abs(theMomentumChange.mag2() - 1.0);
    if (accuracy > GetAccuracyForWarning()) {
      G4ExceptionDescription ed;
      ed << "The Momentum Change is not a unit vector, difference = " << accuracy << G4endl;
      G4Exception("G4FastStep::CheckIt(const G4Track& aTrack)", "FastSim007", JustWarning, ed);
      itsOK = itsOKforMomentum = false;
      if (accuracy > GetAccuracyForException()) {
        exitWithError = true;
      }
    }
  }
 
  accuracy = (aTrack.GetGlobalTime() - theTimeChange) / ns;
  if (accuracy > GetAccuracyForWarning()) {
    G4ExceptionDescription ed;
    ed << "The global time is getting backward, difference = " << accuracy << " ns" << G4endl;
    G4Exception("G4FastStep::CheckIt(const G4Track& aTrack)", "FastSim008", JustWarning, ed);
    itsOK = false;
  }
 
  accuracy = (aTrack.GetProperTime() - theProperTimeChange) / ns;
  if (accuracy > GetAccuracyForWarning()) {
    G4ExceptionDescription ed;
    ed << "The proper time is getting backward, difference = " << accuracy << " ns" << G4endl;
    G4Exception("G4FastStep::CheckIt(const G4Track& aTrack)", "FastSim009", JustWarning, ed);
    itsOK = false;
  }
 
  if (!itsOK) {
    G4cout << "ERROR - G4FastStep::CheckIt() " << G4endl;
    G4cout << "        Pointer : " << this << G4endl;
    DumpInfo();
  }
 
  // Exit with error
  if (exitWithError) {
    G4ExceptionDescription ed;
    ed << "An inaccuracy in G4FastStep is beyond tolerance." << G4endl;
    G4Exception("G4FastStep::CheckIt(const G4Track& aTrack)", "FastSim010", FatalException, ed);
  }
 
  // correction for Momentum only.
  if (!itsOKforMomentum) {
    G4double vmag = theMomentumChange.mag();
    theMomentumChange = (1. / vmag) * theMomentumChange;
  }
 
  itsOK = (itsOK) && G4VParticleChange::CheckIt(aTrack);
  return itsOK;
}

Referenced by UpdateStepForAtRest(), and UpdateStepForPostStep().

CreateSecondaryTrack() [1/2]

G4Track * G4FastStep::CreateSecondaryTrack	(	const G4DynamicParticle &	dynamics,
		G4ThreeVector	position,
		G4double	time,
		G4bool	localCoordinates = true )

Definition at line 199 of file G4FastStep.cc.

{
  // ----------------------------------------
  // Quantities in global coordinates system.
  //
  // The allocated globalDynamics is deleted
  // by the destructor of the G4Track.
  // ----------------------------------------
  auto globalDynamics = new G4DynamicParticle(dynamics);
  G4ThreeVector globalPosition(position);
 
  // -----------------------------------
  // Convert to global system if needed:
  // -----------------------------------
  if (localCoordinates) {
    // -- Momentum Direction:
    globalDynamics->SetMomentumDirection(
      fFastTrack->GetInverseAffineTransformation()->TransformAxis(
        globalDynamics->GetMomentumDirection()));
    // -- Polarization:
    G4ThreeVector globalPolarization;
    globalPolarization = fFastTrack->GetInverseAffineTransformation()->TransformAxis(
      globalDynamics->GetPolarization());
    globalDynamics->SetPolarization(globalPolarization.x(), globalPolarization.y(),
                                    globalPolarization.z());
 
    // -- Position:
    globalPosition = fFastTrack->GetInverseAffineTransformation()->TransformPoint(globalPosition);
  }
 
  //-------------------------------------
  // Create the G4Track of the secondary:
  //-------------------------------------
  auto secondary = new G4Track(globalDynamics, time, globalPosition);
 
  //-------------------------------
  // and feed the changes:
  //-------------------------------
  AddSecondary(secondary);
 
  //--------------------------------------
  // returns the pointer on the secondary:
  //--------------------------------------
  return secondary;
}

CreateSecondaryTrack() [2/2]

G4Track * G4FastStep::CreateSecondaryTrack	(	const G4DynamicParticle &	dynamics,
		G4ThreeVector	polarization,
		G4ThreeVector	position,
		G4double	time,
		G4bool	localCoordinates = true )

Definition at line 182 of file G4FastStep.cc.

{
  G4DynamicParticle dummyDynamics(dynamics);
 
  // ------------------------------------------
  // Add the polarization to the dummyDynamics:
  // ------------------------------------------
  dummyDynamics.SetPolarization(polarization.x(), polarization.y(), polarization.z());
 
  return CreateSecondaryTrack(dummyDynamics, position, time, localCoordinates);
}

Referenced by CreateSecondaryTrack(), and G4BaierKatkov::GeneratePhoton().

DumpInfo()

void G4FastStep::DumpInfo ( ) const

overridevirtual

Reimplemented from G4VParticleChange.

Definition at line 330 of file G4FastStep.cc.

{
  // use base-class DumpInfo
  G4VParticleChange::DumpInfo();
 
  G4cout << "        Position - x (mm)   : " << G4BestUnit(thePositionChange.x(), "Length")
         << G4endl;
  G4cout << "        Position - y (mm)   : " << G4BestUnit(thePositionChange.y(), "Length")
         << G4endl;
  G4cout << "        Position - z (mm)   : " << G4BestUnit(thePositionChange.z(), "Length")
         << G4endl;
  G4cout << "        Time (ns)           : " << G4BestUnit(theTimeChange, "Time") << G4endl;
  G4cout << "        Proper Time (ns)    : " << G4BestUnit(theProperTimeChange, "Time") << G4endl;
  G4long olprc = G4cout.precision(3);
  G4cout << "        Momentum Direct - x : " << std::setw(20) << theMomentumChange.x() << G4endl;
  G4cout << "        Momentum Direct - y : " << std::setw(20) << theMomentumChange.y() << G4endl;
  G4cout << "        Momentum Direct - z : " << std::setw(20) << theMomentumChange.z() << G4endl;
  G4cout.precision(olprc);
  G4cout << "        Kinetic Energy (MeV): " << G4BestUnit(theEnergyChange, "Energy") << G4endl;
  G4cout.precision(3);
  G4cout << "        Polarization - x    : " << std::setw(20) << thePolarizationChange.x()
         << G4endl;
  G4cout << "        Polarization - y    : " << std::setw(20) << thePolarizationChange.y()
         << G4endl;
  G4cout << "        Polarization - z    : " << std::setw(20) << thePolarizationChange.z()
         << G4endl;
  G4cout.precision(olprc);
}

Referenced by CheckIt().

ForceSteppingHitInvocation()

void G4FastStep::ForceSteppingHitInvocation

(

)

GetNumberOfSecondaryTracks()

G4int G4FastStep::GetNumberOfSecondaryTracks

(

)

GetSecondaryTrack()

G4Track * G4FastStep::GetSecondaryTrack

(

G4int

)

GetTotalEnergyDeposited()

G4double G4FastStep::GetTotalEnergyDeposited

(

)

const

Initialize()

void G4FastStep::Initialize

(

const G4FastTrack &

fastTrack

)

Definition at line 53 of file G4FastStep.cc.

{
  // keeps the fastTrack reference
  fFastTrack = &fastTrack;
 
  // currentTrack will be used to Initialize the other data members
  const G4Track& currentTrack = *(fFastTrack->GetPrimaryTrack());
 
  // use base class's method at first
  G4VParticleChange::Initialize(currentTrack);
 
  // set Energy/Momentum etc. equal to those of the parent particle
  const G4DynamicParticle* pParticle = currentTrack.GetDynamicParticle();
  theEnergyChange = pParticle->GetKineticEnergy();
  theMomentumChange = pParticle->GetMomentumDirection();
  thePolarizationChange = pParticle->GetPolarization();
  theProperTimeChange = pParticle->GetProperTime();
 
  // set Position/Time etc. equal to those of the parent track
  thePositionChange = currentTrack.GetPosition();
  theTimeChange = currentTrack.GetGlobalTime();
 
  // switch off stepping hit invokation by default:
  theSteppingControlFlag = AvoidHitInvocation;
 
  // event biasing weigth:
  theWeightChange = currentTrack.GetWeight();
}

Referenced by G4FastSimulationManager::AtRestGetFastSimulationManagerTrigger(), and G4FastSimulationManager::PostStepGetFastSimulationManagerTrigger().

KillPrimaryTrack()

void G4FastStep::KillPrimaryTrack

(

)

Definition at line 87 of file G4FastStep.cc.

{
  ProposePrimaryTrackFinalKineticEnergy(0.);
  ProposeTrackStatus(fStopAndKill);
}

operator=()

G4FastStep & G4FastStep::operator= ( const G4FastStep & right )

delete

ProposePrimaryTrackFinalEventBiasingWeight()

void G4FastStep::ProposePrimaryTrackFinalEventBiasingWeight

(

G4double

)

ProposePrimaryTrackFinalKineticEnergy()

void G4FastStep::ProposePrimaryTrackFinalKineticEnergy

(

G4double

)

Referenced by G4ChannelingFastSimModel::DoIt(), KillPrimaryTrack(), and ProposePrimaryTrackFinalKineticEnergyAndDirection().

ProposePrimaryTrackFinalKineticEnergyAndDirection()

void G4FastStep::ProposePrimaryTrackFinalKineticEnergyAndDirection	(	G4double	kineticEnergy,
		const G4ThreeVector &	direction,
		G4bool	localCoordinates = true )

Definition at line 138 of file G4FastStep.cc.

{
  // Compute global direction if needed...
  G4ThreeVector globalDirection = direction;
  if (localCoordinates)
    globalDirection = fFastTrack->GetInverseAffineTransformation()->TransformAxis(direction);
  // ...and feed the globalMomentum (ensuring unitarity)
  SetMomentumChange(globalDirection.unit());
  ProposePrimaryTrackFinalKineticEnergy(kineticEnergy);
}

Referenced by SetPrimaryTrackFinalKineticEnergyAndDirection().

ProposePrimaryTrackFinalMomentumDirection()

void G4FastStep::ProposePrimaryTrackFinalMomentumDirection	(	const G4ThreeVector &	momentum,
		G4bool	localCoordinates = true )

Definition at line 117 of file G4FastStep.cc.

{
  // Compute the momentum in global reference
  // system if needed ...
  G4ThreeVector globalMomentum = momentum;
  if (localCoordinates)
    globalMomentum = fFastTrack->GetInverseAffineTransformation()->TransformAxis(momentum);
  // ...and feed the globalMomentum (ensuring unitarity)
  SetMomentumChange(globalMomentum.unit());
}

Referenced by G4ChannelingFastSimModel::DoIt(), and SetPrimaryTrackFinalMomentum().

ProposePrimaryTrackFinalPolarization()

void G4FastStep::ProposePrimaryTrackFinalPolarization	(	const G4ThreeVector &	polarization,
		G4bool	localCoordinates = true )

Definition at line 161 of file G4FastStep.cc.

{
  // Compute polarization in global system if needed:
  G4ThreeVector globalPolarization(polarization);
  if (localCoordinates)
    globalPolarization =
      fFastTrack->GetInverseAffineTransformation()->TransformAxis(globalPolarization);
  // Feed the particle globalPolarization:
  thePolarizationChange = globalPolarization;
}

Referenced by SetPrimaryTrackFinalPolarization().

ProposePrimaryTrackFinalPosition()

void G4FastStep::ProposePrimaryTrackFinalPosition	(	const G4ThreeVector &	position,
		G4bool	localCoordinates = true )

Definition at line 96 of file G4FastStep.cc.

{
  // Compute the position coordinate in global
  // reference system if needed ...
  G4ThreeVector globalPosition = position;
  if (localCoordinates)
    globalPosition = fFastTrack->GetInverseAffineTransformation()->TransformPoint(position);
  // ...and feed the globalPosition:
  thePositionChange = globalPosition;
}

Referenced by G4ChannelingFastSimModel::DoIt(), and SetPrimaryTrackFinalPosition().

ProposePrimaryTrackFinalProperTime()

void G4FastStep::ProposePrimaryTrackFinalProperTime

(

G4double

)

ProposePrimaryTrackFinalTime()

void G4FastStep::ProposePrimaryTrackFinalTime

(

G4double

)

Referenced by G4ChannelingFastSimModel::DoIt().

ProposePrimaryTrackPathLength()

void G4FastStep::ProposePrimaryTrackPathLength

(

G4double

)

ProposeTotalEnergyDeposited()

void G4FastStep::ProposeTotalEnergyDeposited

(

G4double

anEnergyPart

)

SetNumberOfSecondaryTracks()

void G4FastStep::SetNumberOfSecondaryTracks

(

G4int

)

Referenced by G4ChannelingFastSimModel::DoIt().

SetPrimaryTrackFinalEventBiasingWeight()

void G4FastStep::SetPrimaryTrackFinalEventBiasingWeight

(

G4double

)

SetPrimaryTrackFinalKineticEnergy()

void G4FastStep::SetPrimaryTrackFinalKineticEnergy

(

G4double

)

SetPrimaryTrackFinalKineticEnergyAndDirection()

void G4FastStep::SetPrimaryTrackFinalKineticEnergyAndDirection	(	G4double	kineticEnergy,
		const G4ThreeVector &	direction,
		G4bool	localCoordinates = true )

Definition at line 151 of file G4FastStep.cc.

{
  ProposePrimaryTrackFinalKineticEnergyAndDirection(kineticEnergy, direction, localCoordinates);
}

SetPrimaryTrackFinalMomentum()

void G4FastStep::SetPrimaryTrackFinalMomentum	(	const G4ThreeVector &	momentum,
		G4bool	localCoordinates = true )

Definition at line 129 of file G4FastStep.cc.

{
  ProposePrimaryTrackFinalMomentumDirection(momentum, localCoordinates);
}

SetPrimaryTrackFinalPolarization()

void G4FastStep::SetPrimaryTrackFinalPolarization	(	const G4ThreeVector &	polarization,
		G4bool	localCoordinates = true )

Definition at line 173 of file G4FastStep.cc.

{
  ProposePrimaryTrackFinalPolarization(polarization, localCoordinates);
}

SetPrimaryTrackFinalPosition()

void G4FastStep::SetPrimaryTrackFinalPosition	(	const G4ThreeVector &	position,
		G4bool	localCoordinates = true )

Definition at line 108 of file G4FastStep.cc.

{
  ProposePrimaryTrackFinalPosition(position, localCoordinates);
}

SetPrimaryTrackFinalProperTime()

void G4FastStep::SetPrimaryTrackFinalProperTime

(

G4double

)

SetPrimaryTrackFinalTime()

void G4FastStep::SetPrimaryTrackFinalTime

(

G4double

)

SetPrimaryTrackPathLength()

void G4FastStep::SetPrimaryTrackPathLength

(

G4double

)

SetTotalEnergyDeposited()

void G4FastStep::SetTotalEnergyDeposited

(

G4double

anEnergyPart

)

UpdateStepForAtRest()

G4Step * G4FastStep::UpdateStepForAtRest ( G4Step * Step )

overridevirtual

Reimplemented from G4VParticleChange.

Definition at line 295 of file G4FastStep.cc.

{
  // A physics process always calculates the final state of the particle
 
  // G4StepPoint* pPreStepPoint  = pStep->GetPreStepPoint();
  G4StepPoint* pPostStepPoint = pStep->GetPostStepPoint();
  G4Track* aTrack = pStep->GetTrack();
  // G4double     mass = aTrack->GetDynamicParticle()->GetMass();
 
  // update kinetic energy and momentum direction
  pPostStepPoint->SetMomentumDirection(theMomentumChange);
  pPostStepPoint->SetKineticEnergy(theEnergyChange);
 
  // update polarization
  pPostStepPoint->SetPolarization(thePolarizationChange);
 
  // update position and time
  pPostStepPoint->SetPosition(thePositionChange);
  pPostStepPoint->SetGlobalTime(theTimeChange);
  pPostStepPoint->AddLocalTime(theTimeChange - aTrack->GetGlobalTime());
  pPostStepPoint->SetProperTime(theProperTimeChange);
 
  // update weight
  pPostStepPoint->SetWeight(theWeightChange);
 
  if (debugFlag) CheckIt(*aTrack);
 
  //  Update the G4Step specific attributes
  return UpdateStepInfo(pStep);
}

UpdateStepForPostStep()

G4Step * G4FastStep::UpdateStepForPostStep ( G4Step * Step )

overridevirtual

Reimplemented from G4VParticleChange.

Definition at line 260 of file G4FastStep.cc.

{
  // A physics process always calculates the final state of the particle
 
  // Take note that the return type of GetMomentumChange is a
  // pointer to G4ParticleMometum. Also it is a normalized
  // momentum vector.
 
  //  G4StepPoint* pPreStepPoint  = pStep->GetPreStepPoint();
  G4StepPoint* pPostStepPoint = pStep->GetPostStepPoint();
  G4Track* aTrack = pStep->GetTrack();
  //  G4double     mass = aTrack->GetDynamicParticle()->GetMass();
 
  // update kinetic energy and momentum direction
  pPostStepPoint->SetMomentumDirection(theMomentumChange);
  pPostStepPoint->SetKineticEnergy(theEnergyChange);
 
  // update polarization
  pPostStepPoint->SetPolarization(thePolarizationChange);
 
  // update position and time
  pPostStepPoint->SetPosition(thePositionChange);
  pPostStepPoint->SetGlobalTime(theTimeChange);
  pPostStepPoint->AddLocalTime(theTimeChange - aTrack->GetGlobalTime());
  pPostStepPoint->SetProperTime(theProperTimeChange);
 
  // update weight
  pPostStepPoint->SetWeight(theWeightChange);
 
  if (debugFlag) CheckIt(*aTrack);
 
  //  Update the G4Step specific attributes
  return UpdateStepInfo(pStep);
}

---

The documentation for this class was generated from the following files:

• G4FastStep.hh
• G4FastStep.cc