G4VParticleChange Class Reference

#include <G4VParticleChange.hh>

Inheritance diagram for G4VParticleChange:

Collaboration diagram for G4VParticleChange:

Public Member Functions
	G4VParticleChange ()
virtual	~G4VParticleChange ()
G4bool	operator== (const G4VParticleChange &right) const
G4bool	operator!= (const G4VParticleChange &right) const
virtual G4Step *	UpdateStepForAtRest (G4Step *Step)
virtual G4Step *	UpdateStepForAlongStep (G4Step *Step)
virtual G4Step *	UpdateStepForPostStep (G4Step *Step)
virtual void	Initialize (const G4Track &)
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
virtual void	DumpInfo () const
void	SetVerboseLevel (G4int vLevel)
G4int	GetVerboseLevel () const
virtual G4bool	CheckIt (const G4Track &)
void	ClearDebugFlag ()
void	SetDebugFlag ()
G4bool	GetDebugFlag () const

Protected Member Functions
	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

Protected Attributes
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
static const G4double	accuracyForWarning = 1.0e-9
static const G4double	accuracyForException = 0.001

Detailed Description

Definition at line 94 of file G4VParticleChange.hh.

Constructor & Destructor Documentation

G4VParticleChange::G4VParticleChange

(

)

Definition at line 48 of file G4VParticleChange.cc.

  :theListOfSecondaries(0),
   theNumberOfSecondaries(0),
   theSizeOftheListOfSecondaries(G4TrackFastVectorSize),
   theStatusChange(fAlive),
   theSteppingControlFlag(NormalCondition),     
   theLocalEnergyDeposit(0.0),
   theNonIonizingEnergyDeposit(0.0),
   theTrueStepLength(0.0),
   theFirstStepInVolume(false),
   theLastStepInVolume(false),
   theParentWeight(1.0),
   isParentWeightProposed(false),
   fSetSecondaryWeightByProcess(false),
   theParentGlobalTime(0.0),
   verboseLevel(1), 
   debugFlag(false)
{
#ifdef G4VERBOSE
  // activate CHeckIt if in VERBOSE mode
  debugFlag = true;
#endif
  theListOfSecondaries = new G4TrackFastVector();
}

G4VParticleChange::~G4VParticleChange ( )

virtual

Definition at line 73 of file G4VParticleChange.cc.

                                      {
  // check if tracks still exist in theListOfSecondaries
  if (theNumberOfSecondaries>0) {
#ifdef G4VERBOSE
    if (verboseLevel>0) {
      G4cout << "G4VParticleChange::~G4VParticleChange() Warning  ";
      G4cout << "theListOfSecondaries is not empty ";
    }
#endif
    for (G4int index= 0; index<theNumberOfSecondaries; index++){
      delete (*theListOfSecondaries)[index] ;
    }
  }
  delete theListOfSecondaries; 
}

G4VParticleChange::G4VParticleChange ( const G4VParticleChange &  right )

protected

Definition at line 89 of file G4VParticleChange.cc.

  :theListOfSecondaries(0),
   theNumberOfSecondaries(0),
   theSizeOftheListOfSecondaries(G4TrackFastVectorSize),
   theStatusChange( right.theStatusChange),
   theSteppingControlFlag(right.theSteppingControlFlag),
   theLocalEnergyDeposit(right.theLocalEnergyDeposit),
   theNonIonizingEnergyDeposit(right.theNonIonizingEnergyDeposit),
   theTrueStepLength(right.theTrueStepLength),
   theFirstStepInVolume( right.theFirstStepInVolume),
   theLastStepInVolume(right.theLastStepInVolume),
   theParentWeight(right.theParentWeight),
   isParentWeightProposed(false),
   fSetSecondaryWeightByProcess(right.fSetSecondaryWeightByProcess),
   theParentGlobalTime(0.0),
   verboseLevel(right.verboseLevel),
   debugFlag(right.debugFlag)
{
  theListOfSecondaries =  new G4TrackFastVector();
  theNumberOfSecondaries = right.theNumberOfSecondaries;
  for (G4int index = 0; index<theNumberOfSecondaries; index++){
    G4Track* newTrack =  new G4Track(*((*right.theListOfSecondaries)[index] ));
    theListOfSecondaries->SetElement(index, newTrack);              
  }
}

Here is the call graph for this function:

Member Function Documentation

void G4VParticleChange::AddSecondary

(

G4Track *

aSecondary

)

Definition at line 171 of file G4VParticleChange.cc.

{
  if (debugFlag) CheckSecondary(*aTrack);

  // add a secondary after size check
  if (theSizeOftheListOfSecondaries > theNumberOfSecondaries) {
    // Set weight of secondary tracks
    if (!fSetSecondaryWeightByProcess) aTrack->SetWeight(theParentWeight);
    theListOfSecondaries->SetElement(theNumberOfSecondaries, aTrack);
    theNumberOfSecondaries++;
  } else {
    delete aTrack;

#ifdef G4VERBOSE
    if (verboseLevel>0) {
      G4cout << "G4VParticleChange::AddSecondary() Warning  ";
      G4cout << "theListOfSecondaries is full !! " << G4endl;
      G4cout << " The track is deleted " << G4endl;
    }
#endif
    G4Exception("G4VParticleChange::AddSecondary",
                "TRACK101", JustWarning,
                "Secondary Bug is full. The track is deleted"); 
  }
}

Here is the call graph for this function:

Here is the caller graph for this function:

G4bool G4VParticleChange::CheckIt ( const G4Track &  )

virtual

Reimplemented in G4FastStep, G4ParticleChange, G4ParticleChangeForDecay, G4ParticleChangeForLoss, G4ParticleChangeForMSC, and G4ParticleChangeForGamma.

Definition at line 316 of file G4VParticleChange.cc.

{

  G4bool    exitWithError = false;
  G4double  accuracy;
  static G4ThreadLocal G4int nError = 0;
#ifdef G4VERBOSE
  const  G4int maxError = 30;
#endif

  // Energy deposit should not be negative
  G4bool itsOKforEnergy = true;
  accuracy = -1.0*theLocalEnergyDeposit/MeV;
  if (accuracy > accuracyForWarning) {
    itsOKforEnergy = false;
    nError += 1;
    exitWithError =  (accuracy > accuracyForException);
#ifdef G4VERBOSE
    if (nError < maxError) {
      G4cout << "  G4VParticleChange::CheckIt    : ";
      G4cout << "the energy deposit  is negative  !!" 
         << "  Difference:  " << accuracy  << "[MeV] " <<G4endl;
      G4cout << aTrack.GetDefinition()->GetParticleName()
         << " E=" << aTrack.GetKineticEnergy()/MeV
         << " pos=" << aTrack.GetPosition().x()/m
         << ", " << aTrack.GetPosition().y()/m
         << ", " << aTrack.GetPosition().z()/m
         <<G4endl;
    }
#endif
  }
 
  // true path length should not be negative
  G4bool itsOKforStepLength = true;
  accuracy = -1.0*theTrueStepLength/mm;
  if (accuracy > accuracyForWarning) {
    itsOKforStepLength = false;
    nError += 1;
    exitWithError =  (accuracy > accuracyForException);
#ifdef G4VERBOSE
    if (nError < maxError) {
      G4cout << "  G4VParticleChange::CheckIt    : ";
      G4cout << "the true step length is negative  !!"
         << "  Difference:  " << accuracy  << "[MeV] " <<G4endl;
      G4cout << aTrack.GetDefinition()->GetParticleName()
         << " E=" << aTrack.GetKineticEnergy()/MeV
         << " pos=" << aTrack.GetPosition().x()/m
         << ", " << aTrack.GetPosition().y()/m
         << ", " << aTrack.GetPosition().z()/m
         <<G4endl;
    }
#endif

  }
#ifdef G4VERBOSE
  if (!itsOKforStepLength || !itsOKforEnergy ){
  // dump out information of this particle change
    DumpInfo();
  }
#endif

  // Exit with error
  if (exitWithError) {
    G4Exception("G4VParticleChange::CheckIt",
        "TRACK001", EventMustBeAborted,
        "Step length and/or energy deposit was illegal");
  }

  // correction 
  if ( !itsOKforStepLength ) {
    theTrueStepLength =  (1.e-12)*mm;
  } 
  if ( !itsOKforEnergy ) {
    theLocalEnergyDeposit = 0.0;
  }
  return (itsOKforStepLength && itsOKforEnergy );
}

Here is the call graph for this function:

Here is the caller graph for this function:

G4bool G4VParticleChange::CheckSecondary ( G4Track &  aTrack )

protected

Definition at line 398 of file G4VParticleChange.cc.

{
  G4bool    exitWithError = false;
  G4double  accuracy;
  static G4ThreadLocal G4int nError = 0;
#ifdef G4VERBOSE
  const  G4int maxError = 30;
#endif

  // MomentumDirection should be unit vector
  G4bool itsOKforMomentum = true;  
  if (aTrack.GetKineticEnergy()>0.){
    accuracy = std::fabs((aTrack.GetMomentumDirection()).mag2()-1.0);
    if (accuracy > accuracyForWarning) {
      itsOKforMomentum = false;
      nError += 1;
      exitWithError = exitWithError || (accuracy > accuracyForException);
#ifdef G4VERBOSE
      if (nError < maxError) {
    G4cout << " G4VParticleChange::CheckSecondary  :   ";
    G4cout << "the Momentum direction is not unit vector !! " 
           << "  Difference:  " << accuracy << G4endl;
    G4cout << aTrack.GetDefinition()->GetParticleName()
           << " E=" << aTrack.GetKineticEnergy()/MeV
           << " pos=" << aTrack.GetPosition().x()/m
           << ", " << aTrack.GetPosition().y()/m
           << ", " << aTrack.GetPosition().z()/m
           <<G4endl;
      }
#endif
    }
  }
  
  // Kinetic Energy should not be negative
  G4bool itsOKforEnergy = true;
  accuracy = -1.0*(aTrack.GetKineticEnergy())/MeV;
  if (accuracy > accuracyForWarning) {
    itsOKforEnergy = false;
    nError += 1;
    exitWithError = exitWithError ||  (accuracy > accuracyForException);
#ifdef G4VERBOSE
    if (nError < maxError) {
      G4cout << " G4VParticleChange::CheckSecondary  :   ";
      G4cout << "the kinetic energy is negative  !!" 
         << "  Difference:  " << accuracy  << "[MeV] " <<G4endl;
      G4cout << " G4VParticleChange::CheckSecondary  :   ";
      G4cout << "the global time of secondary is earlier than the parent  !!" 
         << "  Difference:  " << accuracy  << "[ns] " <<G4endl;
      G4cout << aTrack.GetDefinition()->GetParticleName()
         << " E=" << aTrack.GetKineticEnergy()/MeV
         << " pos=" << aTrack.GetPosition().x()/m
         << ", " << aTrack.GetPosition().y()/m
         << ", " << aTrack.GetPosition().z()/m
       <<G4endl;
    }
#endif
  }
  // Check timing of secondaries
  G4bool itsOKforTiming = true;

  accuracy = (theParentGlobalTime - aTrack.GetGlobalTime())/ns;
  if (accuracy > accuracyForWarning){
    itsOKforTiming = false;
    nError += 1;
    exitWithError = (accuracy > accuracyForException);
#ifdef G4VERBOSE
    if (nError < maxError) {
      G4cout << " G4VParticleChange::CheckSecondary  :   ";
      G4cout << "the global time of secondary goes back comapared to the parent  !!" 
         << "  Difference:  " << accuracy  << "[ns] " <<G4endl;
      G4cout << aTrack.GetDefinition()->GetParticleName()
         << " E=" << aTrack.GetKineticEnergy()/MeV
         << " pos=" << aTrack.GetPosition().x()/m
         << ", " << aTrack.GetPosition().y()/m
           << ", " << aTrack.GetPosition().z()/m
         << " time=" << aTrack.GetGlobalTime()/ns
         << " parent time=" << theParentGlobalTime/ns
         <<G4endl;
    }
#endif
  }

  // Exit with error
  if (exitWithError) {
    G4Exception("G4VParticleChange::CheckSecondary",
        "TRACK001", EventMustBeAborted,
        "Secondary with illegal energy/momentum ");
  }

  G4bool itsOK = itsOKforMomentum && itsOKforEnergy && itsOKforTiming;

  //correction
  if (!itsOKforMomentum) {
    G4double vmag = (aTrack.GetMomentumDirection()).mag();
    aTrack.SetMomentumDirection((1./vmag)*aTrack.GetMomentumDirection());
  }
  if (!itsOKforEnergy) {
    aTrack.SetKineticEnergy(0.0);
  }
 
  if (!itsOK) {
    this->DumpInfo();
    
  }


  return itsOK;
}

Here is the call graph for this function:

Here is the caller graph for this function:

void G4VParticleChange::Clear

(

)

Here is the caller graph for this function:

void G4VParticleChange::ClearDebugFlag

(

)

Here is the caller graph for this function:

void G4VParticleChange::DumpInfo ( ) const

virtual

Reimplemented in G4FastStep, G4ParticleChange, G4ParticleChangeForLoss, G4ParticleChangeForTransport, G4ParticleChangeForDecay, G4ParticleChangeForMSC, and G4ParticleChangeForGamma.

Definition at line 252 of file G4VParticleChange.cc.

{

// Show header
  G4int olprc = G4cout.precision(3);
  G4cout << "      -----------------------------------------------" 
       << G4endl;
  G4cout << "        G4ParticleChange Information  " << std::setw(20) << G4endl;
  G4cout << "      -----------------------------------------------" 
       << G4endl;

  G4cout << "        # of 2ndaries       : " 
       << std::setw(20) << theNumberOfSecondaries
       << G4endl;

  if (theNumberOfSecondaries >0) {
    G4cout << "        Pointer to 2ndaries : " 
         << std::setw(20) << GetSecondary(0)
         << G4endl;
    G4cout << "        (Showed only 1st one)"
         << G4endl;
  }
  G4cout << "      -----------------------------------------------" 
       << G4endl;

  G4cout << "        Energy Deposit (MeV): " 
       << std::setw(20) << theLocalEnergyDeposit/MeV
       << G4endl;

  G4cout << "        Non-ionizing Energy Deposit (MeV): " 
       << std::setw(20) << theNonIonizingEnergyDeposit/MeV
       << G4endl;

  G4cout << "        Track Status        : " 
       << std::setw(20);
       if( theStatusChange == fAlive ){
         G4cout << " Alive";
       } else if( theStatusChange == fStopButAlive ){
           G4cout << " StopButAlive";
       } else if( theStatusChange == fStopAndKill ){
           G4cout << " StopAndKill";
       } else if( theStatusChange  == fKillTrackAndSecondaries ){
           G4cout << " KillTrackAndSecondaries";
       } else if( theStatusChange  == fSuspend ){
           G4cout << " Suspend";
       } else if( theStatusChange == fPostponeToNextEvent ){
           G4cout << " PostponeToNextEvent";
       }
       G4cout << G4endl;
  G4cout << "        True Path Length (mm) : " 
       << std::setw(20) << theTrueStepLength/mm
       << G4endl;
  G4cout << "        Stepping Control      : " 
       << std::setw(20) << theSteppingControlFlag
       << G4endl;   
  if (theFirstStepInVolume) {
    G4cout << "    First Step In the voulme  : "  << G4endl;
  }
  if (theLastStepInVolume) {
    G4cout << "    Last Step In the voulme  : "  << G4endl;
  }
  G4cout.precision(olprc);
}

Here is the call graph for this function:

Here is the caller graph for this function:

G4double G4VParticleChange::GetAccuracyForException ( ) const

protected

Definition at line 513 of file G4VParticleChange.cc.

{
  return accuracyForException;
}

Here is the caller graph for this function:

G4double G4VParticleChange::GetAccuracyForWarning ( ) const

protected

Definition at line 508 of file G4VParticleChange.cc.

{
  return accuracyForWarning;
}

Here is the caller graph for this function:

G4bool G4VParticleChange::GetDebugFlag

(

)

const

G4bool G4VParticleChange::GetFirstStepInVolume

(

)

const

Here is the caller graph for this function:

G4bool G4VParticleChange::GetLastStepInVolume

(

)

const

Here is the caller graph for this function:

G4double G4VParticleChange::GetLocalEnergyDeposit

(

)

const

Here is the caller graph for this function:

G4double G4VParticleChange::GetNonIonizingEnergyDeposit

(

)

const

G4int G4VParticleChange::GetNumberOfSecondaries

(

)

const

Here is the caller graph for this function:

G4double G4VParticleChange::GetParentWeight

(

)

const

Here is the caller graph for this function:

G4Track* G4VParticleChange::GetSecondary

(

G4int

anIndex

)

const

Here is the caller graph for this function:

G4SteppingControl G4VParticleChange::GetSteppingControl

(

)

const

G4TrackStatus G4VParticleChange::GetTrackStatus

(

)

const

Here is the caller graph for this function:

G4double G4VParticleChange::GetTrueStepLength

(

)

const

G4int G4VParticleChange::GetVerboseLevel

(

)

const

G4double G4VParticleChange::GetWeight

(

)

const

virtual void G4VParticleChange::Initialize ( const G4Track &  )

virtual

Reimplemented in G4ParticleChange, G4ParticleChangeForDecay, G4ParticleChangeForTransport, G4ParticleChangeForMSC, and G4ParticleChangeForNothing.

Here is the caller graph for this function:

void G4VParticleChange::InitializeLocalEnergyDeposit ( const G4Track &  )

protected

void G4VParticleChange::InitializeParentGlobalTime ( const G4Track &  )

protected

void G4VParticleChange::InitializeParentWeight ( const G4Track &  )

protected

void G4VParticleChange::InitializeSecondaries ( const G4Track &  )

protected

Here is the caller graph for this function:

void G4VParticleChange::InitializeStatusChange ( const G4Track &  )

protected

void G4VParticleChange::InitializeStepInVolumeFlags ( const G4Track &  )

protected

void G4VParticleChange::InitializeSteppingControl ( const G4Track &  )

protected

void G4VParticleChange::InitializeTrueStepLength ( const G4Track &  )

protected

G4bool G4VParticleChange::IsParentWeightSetByProcess

(

)

const

Definition at line 527 of file G4VParticleChange.cc.

{
  return  true;
}

G4bool G4VParticleChange::IsSecondaryWeightSetByProcess

(

)

const

G4bool G4VParticleChange::operator!=

(

const G4VParticleChange &

right

)

const

Definition at line 166 of file G4VParticleChange.cc.

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

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

protected

Definition at line 116 of file G4VParticleChange.cc.

{
  if (this != &right){
    if (theNumberOfSecondaries>0) {
#ifdef G4VERBOSE
      if (verboseLevel>0) {
    G4cout << "G4VParticleChange: assignment operator Warning  ";
    G4cout << "theListOfSecondaries is not empty ";
      }
#endif
      for (G4int index = 0; index<theNumberOfSecondaries; index++){
    if ( (*theListOfSecondaries)[index] ) delete ((*theListOfSecondaries)[index]) ;
      }
    }
    delete theListOfSecondaries; 
      
    theListOfSecondaries =  new G4TrackFastVector();
    theNumberOfSecondaries = right.theNumberOfSecondaries;
   for (G4int index = 0; index<theNumberOfSecondaries; index++){
    G4Track* newTrack =  new G4Track(*((*right.theListOfSecondaries)[index] ));
    theListOfSecondaries->SetElement(index, newTrack);              
  }
    theStatusChange = right.theStatusChange;
    theSteppingControlFlag = right.theSteppingControlFlag;
    theLocalEnergyDeposit = right.theLocalEnergyDeposit;
    theNonIonizingEnergyDeposit = right.theNonIonizingEnergyDeposit;
    theTrueStepLength = right.theTrueStepLength;
    
    theFirstStepInVolume = right.theFirstStepInVolume;
    theLastStepInVolume =  right.theLastStepInVolume;

    theParentWeight = right.theParentWeight;
    isParentWeightProposed = right.isParentWeightProposed;
    fSetSecondaryWeightByProcess = right.fSetSecondaryWeightByProcess;

    theParentGlobalTime = right.theParentGlobalTime;

    verboseLevel = right.verboseLevel;
    debugFlag = right.debugFlag;
    
  }
  return *this;
}

Here is the call graph for this function:

Here is the caller graph for this function:

G4bool G4VParticleChange::operator==

(

const G4VParticleChange &

right

)

const

Definition at line 160 of file G4VParticleChange.cc.

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

void G4VParticleChange::ProposeFirstStepInVolume

(

G4bool

flag

)

Here is the caller graph for this function:

void G4VParticleChange::ProposeLastStepInVolume

(

G4bool

flag

)

Here is the caller graph for this function:

void G4VParticleChange::ProposeLocalEnergyDeposit

(

G4double

anEnergyPart

)

void G4VParticleChange::ProposeNonIonizingEnergyDeposit

(

G4double

anEnergyPart

)

Here is the caller graph for this function:

void G4VParticleChange::ProposeParentWeight

(

G4double

finalWeight

)

Here is the caller graph for this function:

void G4VParticleChange::ProposeSteppingControl

(

G4SteppingControl

StepControlFlag

)

Here is the caller graph for this function:

void G4VParticleChange::ProposeTrackStatus

(

G4TrackStatus

status

)

void G4VParticleChange::ProposeTrueStepLength

(

G4double

truePathLength

)

Here is the caller graph for this function:

void G4VParticleChange::ProposeWeight

(

G4double

finalWeight

)

Here is the caller graph for this function:

void G4VParticleChange::SetDebugFlag

(

)

void G4VParticleChange::SetNumberOfSecondaries

(

G4int

totSecondaries

)

Here is the caller graph for this function:

void G4VParticleChange::SetParentWeightByProcess

(

G4bool

)

Definition at line 522 of file G4VParticleChange.cc.

{
}

Here is the caller graph for this function:

void G4VParticleChange::SetSecondaryWeightByProcess

(

G4bool

)

Here is the caller graph for this function:

void G4VParticleChange::SetVerboseLevel

(

G4int

vLevel

)

G4Step * G4VParticleChange::UpdateStepForAlongStep ( G4Step *  Step )

virtual

Reimplemented in G4ParticleChange, G4ParticleChangeForMSC, G4ParticleChangeForTransport, G4ParticleChangeForLoss, G4ParticleChangeForOccurenceBiasing, and G4ParticleChangeForNothing.

Definition at line 228 of file G4VParticleChange.cc.

{
  if (isParentWeightProposed ){ 
    G4double initialWeight = Step->GetPreStepPoint()->GetWeight();
    G4double currentWeight = Step->GetPostStepPoint()->GetWeight();
    G4double finalWeight   = (theParentWeight/initialWeight)*currentWeight;
    Step->GetPostStepPoint()->SetWeight( finalWeight );
  }   
  return UpdateStepInfo(Step);
}

Here is the call graph for this function:

Here is the caller graph for this function:

G4Step * G4VParticleChange::UpdateStepForAtRest ( G4Step *  Step )

virtual

Reimplemented in G4FastStep, G4ParticleChange, G4ParticleChangeForDecay, G4ParticleChangeForTransport, G4ParticleChangeForGamma, G4ParticleChangeForOccurenceBiasing, and G4ParticleChangeForNothing.

Definition at line 219 of file G4VParticleChange.cc.

{ 
  if (isParentWeightProposed ){
    Step->GetPostStepPoint()->SetWeight( theParentWeight );
  }
  return UpdateStepInfo(Step);
}

Here is the call graph for this function:

G4Step * G4VParticleChange::UpdateStepForPostStep ( G4Step *  Step )

virtual

Reimplemented in G4FastStep, G4ParticleChange, G4ParticleChangeForDecay, G4ParticleChangeForTransport, G4ParticleChangeForMSC, G4ParticleChangeForLoss, G4ParticleChangeForGamma, G4ParticleChangeForOccurenceBiasing, and G4ParticleChangeForNothing.

Definition at line 239 of file G4VParticleChange.cc.

{
  if (isParentWeightProposed ){ 
    Step->GetPostStepPoint()->SetWeight( theParentWeight );
  }
  return UpdateStepInfo(Step);
}

Here is the call graph for this function:

Here is the caller graph for this function:

G4Step * G4VParticleChange::UpdateStepInfo ( G4Step *  Step )

protected

Definition at line 202 of file G4VParticleChange.cc.

{
  // Update the G4Step specific attributes
  pStep->SetStepLength( theTrueStepLength );
  pStep->AddTotalEnergyDeposit( theLocalEnergyDeposit );
  pStep->AddNonIonizingEnergyDeposit( theNonIonizingEnergyDeposit );
  pStep->SetControlFlag( theSteppingControlFlag );

  if (theFirstStepInVolume) {pStep->SetFirstStepFlag();}
  else                      {pStep->ClearFirstStepFlag();} 
  if (theLastStepInVolume)  {pStep->SetLastStepFlag();}
  else                      {pStep->ClearLastStepFlag();} 

  return pStep;
}

Here is the call graph for this function:

Here is the caller graph for this function:

Member Data Documentation

const G4double G4VParticleChange::accuracyForException = 0.001

staticprotected

Definition at line 318 of file G4VParticleChange.hh.

const G4double G4VParticleChange::accuracyForWarning = 1.0e-9

staticprotected

Definition at line 317 of file G4VParticleChange.hh.

G4bool G4VParticleChange::debugFlag

protected

Definition at line 314 of file G4VParticleChange.hh.

G4bool G4VParticleChange::fSetSecondaryWeightByProcess

protected

Definition at line 286 of file G4VParticleChange.hh.

G4bool G4VParticleChange::isParentWeightProposed

protected

Definition at line 284 of file G4VParticleChange.hh.

G4bool G4VParticleChange::theFirstStepInVolume

protected

Definition at line 278 of file G4VParticleChange.hh.

G4bool G4VParticleChange::theLastStepInVolume

protected

Definition at line 279 of file G4VParticleChange.hh.

G4TrackFastVector* G4VParticleChange::theListOfSecondaries

protected

Definition at line 245 of file G4VParticleChange.hh.

G4double G4VParticleChange::theLocalEnergyDeposit

protected

Definition at line 260 of file G4VParticleChange.hh.

G4double G4VParticleChange::theNonIonizingEnergyDeposit

protected

Definition at line 269 of file G4VParticleChange.hh.

G4int G4VParticleChange::theNumberOfSecondaries

protected

Definition at line 248 of file G4VParticleChange.hh.

G4double G4VParticleChange::theParentGlobalTime

protected

Definition at line 289 of file G4VParticleChange.hh.

G4double G4VParticleChange::theParentWeight

protected

Definition at line 282 of file G4VParticleChange.hh.

G4int G4VParticleChange::theSizeOftheListOfSecondaries

protected

Definition at line 251 of file G4VParticleChange.hh.

G4TrackStatus G4VParticleChange::theStatusChange

protected

Definition at line 254 of file G4VParticleChange.hh.

G4SteppingControl G4VParticleChange::theSteppingControlFlag

protected

Definition at line 257 of file G4VParticleChange.hh.

G4double G4VParticleChange::theTrueStepLength

protected

Definition at line 274 of file G4VParticleChange.hh.

G4int G4VParticleChange::verboseLevel

protected

Definition at line 293 of file G4VParticleChange.hh.

---

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

• source/geant4.10.03.p02/source/track/include/G4VParticleChange.hh
• source/geant4.10.03.p02/source/track/src/G4VParticleChange.cc