G4FastStep Class Reference

#include <G4FastStep.hh>

Inheritance diagram for G4FastStep:

Collaboration diagram for G4FastStep:

Public Member Functions
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 ()
	G4FastStep ()
virtual	~G4FastStep ()
G4bool	operator== (const G4FastStep &right) const
G4bool	operator!= (const G4FastStep &right) const
G4Step *	UpdateStepForAtRest (G4Step *Step)
G4Step *	UpdateStepForPostStep (G4Step *Step)
void	Initialize (const G4FastTrack &)
void	DumpInfo () const
G4bool	CheckIt (const G4Track &)
Public Member Functions inherited from G4VParticleChange
	G4VParticleChange ()
virtual	~G4VParticleChange ()
G4bool	operator== (const G4VParticleChange &right) const
G4bool	operator!= (const G4VParticleChange &right) const
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)
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

Protected Member Functions
	G4FastStep (const G4FastStep &right)
G4FastStep &	operator= (const G4FastStep &right)
Protected Member Functions inherited from G4VParticleChange
	G4VParticleChange (const G4VParticleChange &right)
G4VParticleChange &	operator= (const G4VParticleChange &right)
G4Step *	UpdateStepInfo (G4Step *Step)
void	InitializeTrueStepLength (const G4Track &)
void	InitializeLocalEnergyDeposit (const G4Track &)
void	InitializeSteppingControl (const G4Track &)
void	InitializeParentWeight (const G4Track &)
void	InitializeParentGlobalTime (const G4Track &)
void	InitializeStatusChange (const G4Track &)
void	InitializeSecondaries (const G4Track &)
void	InitializeStepInVolumeFlags (const G4Track &)
G4bool	CheckSecondary (G4Track &)
G4double	GetAccuracyForWarning () const
G4double	GetAccuracyForException () const

Additional Inherited Members
Protected Attributes inherited from G4VParticleChange
G4TrackFastVector *	theListOfSecondaries
G4int	theNumberOfSecondaries
G4int	theSizeOftheListOfSecondaries
G4TrackStatus	theStatusChange
G4SteppingControl	theSteppingControlFlag
G4double	theLocalEnergyDeposit
G4double	theNonIonizingEnergyDeposit
G4double	theTrueStepLength
G4bool	theFirstStepInVolume
G4bool	theLastStepInVolume
G4double	theParentWeight
G4bool	isParentWeightProposed
G4bool	fSetSecondaryWeightByProcess
G4double	theParentGlobalTime
G4int	verboseLevel
G4bool	debugFlag
Static Protected Attributes inherited from G4VParticleChange
static const G4double	accuracyForWarning = 1.0e-9
static const G4double	accuracyForException = 0.001

Detailed Description

Definition at line 91 of file G4FastStep.hh.

Constructor & Destructor Documentation

G4FastStep::G4FastStep

(

)

Definition at line 298 of file G4FastStep.cc.

  : G4VParticleChange(),
    theEnergyChange    ( 0.0     ),
    theTimeChange      ( 0.0     ),
    theProperTimeChange( 0.0     ),
    fFastTrack         ( nullptr ),
    theWeightChange    ( 0.0     )
{
  if (verboseLevel>2)
  {
    G4cerr << "G4FastStep::G4FastStep()" << G4endl;
  }
}

G4FastStep::~G4FastStep ( )

virtual

Definition at line 312 of file G4FastStep.cc.

{
  if (verboseLevel>2)
  {
    G4cerr << "G4FastStep::~G4FastStep()" << G4endl;
  }
}

G4FastStep::G4FastStep ( const G4FastStep &  right )

protected

Definition at line 321 of file G4FastStep.cc.

  : G4VParticleChange()
{
   *this = right;
}

Member Function Documentation

G4bool G4FastStep::CheckIt ( const G4Track &  aTrack )

virtual

Reimplemented from G4VParticleChange.

Definition at line 466 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;
}

Here is the call graph for this function:

Here is the caller graph for this function:

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

Definition at line 203 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);
}

Here is the call graph for this function:

Here is the caller graph for this function:

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

Definition at line 225 of file G4FastStep.cc.

{
  // ----------------------------------------
  // Quantities in global coordinates system.
  //  
  // The allocated globalDynamics is deleted
  // by the destructor of the G4Track.
  // ----------------------------------------
  G4DynamicParticle* 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:
  //-------------------------------------
  G4Track* secondary = new G4Track(
                                   globalDynamics,
                                   time,
                                   globalPosition
                                   );
  
  //-------------------------------
  // and feed the changes:
  //-------------------------------
  AddSecondary(secondary);
  
  //--------------------------------------
  // returns the pointer on the secondary:
  //--------------------------------------
  return secondary;
}

Here is the call graph for this function:

void G4FastStep::DumpInfo ( ) const

virtual

Reimplemented from G4VParticleChange.

Definition at line 443 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;
  G4int 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);
}

Here is the call graph for this function:

Here is the caller graph for this function:

void G4FastStep::ForceSteppingHitInvocation

(

)

G4int G4FastStep::GetNumberOfSecondaryTracks

(

)

G4Track* G4FastStep::GetSecondaryTrack

(

G4int

)

G4double G4FastStep::GetTotalEnergyDeposited

(

)

const

void G4FastStep::Initialize

(

const G4FastTrack &

fastTrack

)

Definition at line 54 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();
}  

Here is the call graph for this function:

Here is the caller graph for this function:

void G4FastStep::KillPrimaryTrack

(

)

Definition at line 88 of file G4FastStep.cc.

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

Here is the call graph for this function:

Here is the caller graph for this function:

G4bool G4FastStep::operator!=

(

const G4FastStep &

right

)

const

Definition at line 361 of file G4FastStep.cc.

{
   return ((G4VParticleChange *)this != (G4VParticleChange *) &right);
}

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

protected

Definition at line 328 of file G4FastStep.cc.

{
   if (this != &right)
   {
     G4VParticleChange::operator=(right);
     theListOfSecondaries          = right.theListOfSecondaries;
     theSizeOftheListOfSecondaries = right.theSizeOftheListOfSecondaries;
     theNumberOfSecondaries        = right.theNumberOfSecondaries;
     theStatusChange               = right.theStatusChange;
     theMomentumChange             = right.theMomentumChange;
     thePolarizationChange         = right.thePolarizationChange;
     thePositionChange             = right.thePositionChange;
     theProperTimeChange           = right.theProperTimeChange;
     theTimeChange                 = right.theTimeChange;
     theEnergyChange               = right.theEnergyChange;
     theTrueStepLength             = right.theTrueStepLength;
     theLocalEnergyDeposit         = right.theLocalEnergyDeposit;
     theSteppingControlFlag        = right.theSteppingControlFlag;
     theWeightChange               = right.theWeightChange;
     fFastTrack                    = right.fFastTrack;
   }
   return *this;
}

Here is the call graph for this function:

G4bool G4FastStep::operator==

(

const G4FastStep &

right

)

const

Definition at line 356 of file G4FastStep.cc.

{
   return ((G4VParticleChange *)this == (G4VParticleChange *) &right);
}

void G4FastStep::ProposePrimaryTrackFinalEventBiasingWeight

(

G4double

)

void G4FastStep::ProposePrimaryTrackFinalKineticEnergy

(

G4double

)

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

Definition at line 151 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());
  SetPrimaryTrackFinalKineticEnergy(kineticEnergy);
}

Here is the call graph for this function:

Here is the caller graph for this function:

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

Definition at line 125 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());
}

Here is the call graph for this function:

Here is the caller graph for this function:

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

Definition at line 179 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;
}

Here is the call graph for this function:

Here is the caller graph for this function:

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

Definition at line 99 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;
}

Here is the call graph for this function:

Here is the caller graph for this function:

void G4FastStep::ProposePrimaryTrackFinalProperTime

(

G4double

)

void G4FastStep::ProposePrimaryTrackFinalTime

(

G4double

)

void G4FastStep::ProposePrimaryTrackPathLength

(

G4double

)

Here is the caller graph for this function:

void G4FastStep::ProposeTotalEnergyDeposited

(

G4double

anEnergyPart

)

Here is the caller graph for this function:

void G4FastStep::SetNumberOfSecondaryTracks

(

G4int

)

Here is the caller graph for this function:

void G4FastStep::SetPrimaryTrackFinalEventBiasingWeight

(

G4double

)

void G4FastStep::SetPrimaryTrackFinalKineticEnergy

(

G4double

)

Here is the caller graph for this function:

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

Definition at line 167 of file G4FastStep.cc.

{
  ProposePrimaryTrackFinalKineticEnergyAndDirection(kineticEnergy, direction, localCoordinates);
}

Here is the call graph for this function:

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

Definition at line 140 of file G4FastStep.cc.

{
  ProposePrimaryTrackFinalMomentumDirection(momentum, localCoordinates);
}

Here is the call graph for this function:

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

Definition at line 193 of file G4FastStep.cc.

{
  ProposePrimaryTrackFinalPolarization(polarization, localCoordinates);
}

Here is the call graph for this function:

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

Definition at line 114 of file G4FastStep.cc.

{
  ProposePrimaryTrackFinalPosition(position, localCoordinates);
}

Here is the call graph for this function:

void G4FastStep::SetPrimaryTrackFinalProperTime

(

G4double

)

void G4FastStep::SetPrimaryTrackFinalTime

(

G4double

)

void G4FastStep::SetPrimaryTrackPathLength

(

G4double

)

void G4FastStep::SetTotalEnergyDeposited

(

G4double

anEnergyPart

)

G4Step * G4FastStep::UpdateStepForAtRest ( G4Step *  Step )

virtual

Reimplemented from G4VParticleChange.

Definition at line 407 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);
}

Here is the call graph for this function:

G4Step * G4FastStep::UpdateStepForPostStep ( G4Step *  Step )

virtual

Reimplemented from G4VParticleChange.

Definition at line 370 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);
}

Here is the call graph for this function:

---

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

• source/geant4.10.03.p03/source/processes/parameterisation/include/G4FastStep.hh
• source/geant4.10.03.p03/source/processes/parameterisation/src/G4FastStep.cc