G4Transportation Class Reference

#include <G4Transportation.hh>

Inheritance diagram for G4Transportation:

Public Member Functions
	G4Transportation (G4int verbosityLevel=1, const G4String &aName="Transportation")
	~G4Transportation ()
G4double	AlongStepGetPhysicalInteractionLength (const G4Track &track, G4double previousStepSize, G4double currentMinimumStep, G4double &currentSafety, G4GPILSelection *selection)
G4VParticleChange *	AlongStepDoIt (const G4Track &track, const G4Step &stepData)
G4VParticleChange *	PostStepDoIt (const G4Track &track, const G4Step &stepData)
G4double	PostStepGetPhysicalInteractionLength (const G4Track &, G4double previousStepSize, G4ForceCondition *pForceCond)
G4bool	FieldExertedForce ()
G4PropagatorInField *	GetPropagatorInField ()
void	SetPropagatorInField (G4PropagatorInField *pFieldPropagator)
G4double	GetThresholdWarningEnergy () const
G4double	GetThresholdImportantEnergy () const
G4int	GetThresholdTrials () const
void	SetThresholdWarningEnergy (G4double newEnWarn)
void	SetThresholdImportantEnergy (G4double newEnImp)
void	SetThresholdTrials (G4int newMaxTrials)
void	SetHighLooperThresholds ()
void	SetLowLooperThresholds ()
void	PushThresholdsToLogger ()
void	ReportLooperThresholds ()
G4double	GetMaxEnergyKilled () const
G4double	GetSumEnergyKilled () const
void	ResetKilledStatistics (G4int report=1)
void	EnableShortStepOptimisation (G4bool optimise=true)
G4double	AtRestGetPhysicalInteractionLength (const G4Track &, G4ForceCondition *)
G4VParticleChange *	AtRestDoIt (const G4Track &, const G4Step &)
void	StartTracking (G4Track *aTrack)
virtual void	ProcessDescription (std::ostream &outFile) const
void	PrintStatistics (std::ostream &outStr) const
Public Member Functions inherited from G4VProcess
	G4VProcess (const G4String &aName="NoName", G4ProcessType aType=fNotDefined)
	G4VProcess (const G4VProcess &right)
virtual	~G4VProcess ()
G4VProcess &	operator= (const G4VProcess &)=delete
G4bool	operator== (const G4VProcess &right) const
G4bool	operator!= (const G4VProcess &right) const
G4double	GetCurrentInteractionLength () const
void	SetPILfactor (G4double value)
G4double	GetPILfactor () const
G4double	AlongStepGPIL (const G4Track &track, G4double previousStepSize, G4double currentMinimumStep, G4double &proposedSafety, G4GPILSelection *selection)
G4double	AtRestGPIL (const G4Track &track, G4ForceCondition *condition)
G4double	PostStepGPIL (const G4Track &track, G4double previousStepSize, G4ForceCondition *condition)
virtual G4bool	IsApplicable (const G4ParticleDefinition &)
virtual void	BuildPhysicsTable (const G4ParticleDefinition &)
virtual void	PreparePhysicsTable (const G4ParticleDefinition &)
virtual G4bool	StorePhysicsTable (const G4ParticleDefinition *, const G4String &, G4bool)
virtual G4bool	RetrievePhysicsTable (const G4ParticleDefinition *, const G4String &, G4bool)
const G4String &	GetPhysicsTableFileName (const G4ParticleDefinition *, const G4String &directory, const G4String &tableName, G4bool ascii=false)
const G4String &	GetProcessName () const
G4ProcessType	GetProcessType () const
void	SetProcessType (G4ProcessType)
G4int	GetProcessSubType () const
void	SetProcessSubType (G4int)
virtual const G4VProcess *	GetCreatorProcess () const
virtual void	EndTracking ()
virtual void	SetProcessManager (const G4ProcessManager *)
virtual const G4ProcessManager *	GetProcessManager ()
virtual void	ResetNumberOfInteractionLengthLeft ()
G4double	GetNumberOfInteractionLengthLeft () const
G4double	GetTotalNumberOfInteractionLengthTraversed () const
G4bool	isAtRestDoItIsEnabled () const
G4bool	isAlongStepDoItIsEnabled () const
G4bool	isPostStepDoItIsEnabled () const
virtual void	DumpInfo () const
void	SetVerboseLevel (G4int value)
G4int	GetVerboseLevel () const
virtual void	SetMasterProcess (G4VProcess *masterP)
const G4VProcess *	GetMasterProcess () const
virtual void	BuildWorkerPhysicsTable (const G4ParticleDefinition &part)
virtual void	PrepareWorkerPhysicsTable (const G4ParticleDefinition &)

Static Public Member Functions
static G4bool	EnableMagneticMoment (G4bool useMoment=true)
static G4bool	EnableGravity (G4bool useGravity=true)
static void	SetSilenceLooperWarnings (G4bool val)
static G4bool	GetSilenceLooperWarnings ()
static G4bool	EnableUseMagneticMoment (G4bool useMoment=true)
Static Public Member Functions inherited from G4VProcess
static const G4String &	GetProcessTypeName (G4ProcessType)

Protected Member Functions
void	SetTouchableInformation (const G4TouchableHandle &touchable)
void	ReportMissingLogger (const char *methodName)
Protected Member Functions inherited from G4VProcess
void	SubtractNumberOfInteractionLengthLeft (G4double prevStepSize)
void	ClearNumberOfInteractionLengthLeft ()

Protected Attributes
G4Navigator *	fLinearNavigator
G4PropagatorInField *	fFieldPropagator
G4ThreeVector	fTransportEndPosition = G4ThreeVector( 0.0, 0.0, 0.0 )
G4ThreeVector	fTransportEndMomentumDir = G4ThreeVector( 0.0, 0.0, 0.0 )
G4double	fTransportEndKineticEnergy = 0.0
G4ThreeVector	fTransportEndSpin = G4ThreeVector( 0.0, 0.0, 0.0 )
G4bool	fMomentumChanged = true
G4bool	fEndGlobalTimeComputed = false
G4double	fCandidateEndGlobalTime = 0.0
G4bool	fAnyFieldExists = false
G4bool	fParticleIsLooping = false
G4bool	fNewTrack = true
G4bool	fFirstStepInVolume = true
G4bool	fLastStepInVolume = false
G4bool	fGeometryLimitedStep = true
G4bool	fFieldExertedForce = false
G4TouchableHandle	fCurrentTouchableHandle
G4ThreeVector	fPreviousSftOrigin
G4double	fPreviousSafety
G4ParticleChangeForTransport	fParticleChange
G4double	fEndPointDistance
G4double	fThreshold_Warning_Energy = 1.0 * CLHEP::keV
G4double	fThreshold_Important_Energy = 1.0 * CLHEP::MeV
G4int	fThresholdTrials = 10
G4int	fAbandonUnstableTrials = 0
G4int	fNoLooperTrials = 0
G4double	fSumEnergyKilled = 0.0
G4double	fSumEnerSqKilled = 0.0
G4double	fMaxEnergyKilled = -1.0
G4int	fMaxEnergyKilledPDG = 0
unsigned long	fNumLoopersKilled = 0
G4double	fSumEnergyKilled_NonElectron = 0.0
G4double	fSumEnerSqKilled_NonElectron = 0.0
G4double	fMaxEnergyKilled_NonElectron = -1.0
G4int	fMaxEnergyKilled_NonElecPDG = 0
unsigned long	fNumLoopersKilled_NonElectron = 0
G4double	fSumEnergySaved = 0.0
G4double	fMaxEnergySaved = -1.0
G4double	fSumEnergyUnstableSaved = 0.0
G4bool	fShortStepOptimisation
G4SafetyHelper *	fpSafetyHelper
G4TransportationLogger *	fpLogger
Protected Attributes inherited from G4VProcess
const G4ProcessManager *	aProcessManager = nullptr
G4VParticleChange *	pParticleChange = nullptr
G4ParticleChange	aParticleChange
G4double	theNumberOfInteractionLengthLeft = -1.0
G4double	currentInteractionLength = -1.0
G4double	theInitialNumberOfInteractionLength = -1.0
G4String	theProcessName
G4String	thePhysicsTableFileName
G4ProcessType	theProcessType = fNotDefined
G4int	theProcessSubType = -1
G4double	thePILfactor = 1.0
G4int	verboseLevel = 0
G4bool	enableAtRestDoIt = true
G4bool	enableAlongStepDoIt = true
G4bool	enablePostStepDoIt = true

Static Protected Attributes
static G4bool	fUseMagneticMoment =false
static G4bool	fUseGravity = false
static G4bool	fSilenceLooperWarnings = false

Detailed Description

Definition at line 57 of file G4Transportation.hh.

Constructor & Destructor Documentation

G4Transportation()

G4Transportation::G4Transportation	(	G4int	verbosityLevel = 1,
		const G4String &	aName = "Transportation" )

Definition at line 74 of file G4Transportation.cc.

  : G4VProcess( aName, fTransportation ),
    fFieldExertedForce( false ),
    fPreviousSftOrigin( 0.,0.,0. ),
    fPreviousSafety( 0.0 ),
    fEndPointDistance( -1.0 ), 
    fShortStepOptimisation( false ) // Old default: true (=fast short steps)
{
  SetProcessSubType(static_cast<G4int>(TRANSPORTATION));
  pParticleChange= &fParticleChange;   // Required to conform to G4VProcess 
  SetVerboseLevel(verbosity);
 
  G4TransportationManager* transportMgr ; 
 
  transportMgr = G4TransportationManager::GetTransportationManager() ; 
 
  fLinearNavigator = transportMgr->GetNavigatorForTracking() ; 
 
  fFieldPropagator = transportMgr->GetPropagatorInField() ;
 
  fpSafetyHelper =   transportMgr->GetSafetyHelper();  // New 
 
  fpLogger = new G4TransportationLogger("G4Transportation", verbosity);
 
  if( G4TransportationParameters::Exists() )
  {
    auto trParams= G4TransportationParameters::Instance();
     
    SetThresholdWarningEnergy(  trParams->GetWarningEnergy() );
    SetThresholdImportantEnergy( trParams->GetImportantEnergy() );
    SetThresholdTrials( trParams->GetNumberOfTrials() );
    G4Transportation::fSilenceLooperWarnings= trParams->GetSilenceAllLooperWarnings();
  }
  else {
     SetHighLooperThresholds();
     // Use the old defaults: Warning = 100 MeV, Important = 250 MeV, No Trials = 10;
  }
 
  PushThresholdsToLogger();
  // Should be done by Set methods in SetHighLooperThresholds -- making sure
  
  static G4ThreadLocal G4TouchableHandle* pNullTouchableHandle = 0;
  if ( !pNullTouchableHandle)
  {
    pNullTouchableHandle = new G4TouchableHandle;
  }
  fCurrentTouchableHandle = *pNullTouchableHandle;
    // Points to (G4VTouchable*) 0
 
  
#ifdef G4VERBOSE
  if( verboseLevel > 0) 
  { 
     G4cout << " G4Transportation constructor> set fShortStepOptimisation to "; 
     if ( fShortStepOptimisation )  { G4cout << "true"  << G4endl; }
     else                           { G4cout << "false" << G4endl; }
  } 
#endif
}

~G4Transportation()

G4Transportation::~G4Transportation

(

)

Definition at line 136 of file G4Transportation.cc.

{
  if( fSumEnergyKilled > 0.0 )
  {
     PrintStatistics( G4cout );     
  }
  delete fpLogger;
}

Member Function Documentation

AlongStepDoIt()

G4VParticleChange * G4Transportation::AlongStepDoIt ( const G4Track & track, const G4Step & stepData )

virtual

Implements G4VProcess.

Definition at line 496 of file G4Transportation.cc.

{
#if defined(G4VERBOSE) || defined(G4DEBUG_TRANSPORT)
  static G4ThreadLocal G4long noCallsASDI=0;
  noCallsASDI++;
#else
  #define noCallsASDI 0
#endif
 
  if(fGeometryLimitedStep)
  {
    stepData.GetPostStepPoint()->SetStepStatus(fGeomBoundary);
  }
 
  fParticleChange.Initialize(track) ;
 
  //  Code for specific process 
  //
  fParticleChange.ProposePosition(fTransportEndPosition) ;
  fParticleChange.ProposeMomentumDirection(fTransportEndMomentumDir) ;
  fParticleChange.ProposeEnergy(fTransportEndKineticEnergy) ;
  fParticleChange.SetMomentumChanged(fMomentumChanged) ;
 
  fParticleChange.ProposePolarization(fTransportEndSpin);
 
  G4double deltaTime = 0.0 ;
 
  // Calculate  Lab Time of Flight (ONLY if field Equations used it!)
  // G4double endTime   = fCandidateEndGlobalTime;
  // G4double delta_time = endTime - startTime;
 
  G4double startTime = track.GetGlobalTime() ;
  
  if (!fEndGlobalTimeComputed)
  {
     // The time was not integrated .. make the best estimate possible
     //
     G4double initialVelocity = stepData.GetPreStepPoint()->GetVelocity();
     G4double stepLength      = track.GetStepLength();
 
     deltaTime= 0.0;  // in case initialVelocity = 0 
     if ( initialVelocity > 0.0 )  { deltaTime = stepLength/initialVelocity; }
 
     fCandidateEndGlobalTime   = startTime + deltaTime ;
     fParticleChange.ProposeLocalTime(  track.GetLocalTime() + deltaTime) ;
  }
  else
  {
     deltaTime = fCandidateEndGlobalTime - startTime ;
     fParticleChange.ProposeGlobalTime( fCandidateEndGlobalTime ) ;
  }
 
 
  // Now Correct by Lorentz factor to get delta "proper" Time
  
  G4double  restMass       = track.GetDynamicParticle()->GetMass() ;
  G4double deltaProperTime = deltaTime*( restMass/track.GetTotalEnergy() ) ;
 
  fParticleChange.ProposeProperTime(track.GetProperTime() + deltaProperTime) ;
  //fParticleChange.ProposeTrueStepLength( track.GetStepLength() ) ;
 
  // If the particle is caught looping or is stuck (in very difficult
  // boundaries) in a magnetic field (doing many steps) THEN can kill it ...
  //
  if ( fParticleIsLooping )
  {
      G4double endEnergy= fTransportEndKineticEnergy;
      fNoLooperTrials ++; 
      auto particleType= track.GetDynamicParticle()->GetParticleDefinition();
      
      G4bool stable = particleType->GetPDGStable();
      G4bool candidateForEnd = (endEnergy < fThreshold_Important_Energy) 
                            || (fNoLooperTrials >= fThresholdTrials) ;
      G4bool unstableAndKillable = !stable && ( fAbandonUnstableTrials != 0);
      G4bool unstableForEnd = (endEnergy < fThreshold_Important_Energy)
                              && (fNoLooperTrials >= fAbandonUnstableTrials) ;
      if( (candidateForEnd && stable) || (unstableAndKillable && unstableForEnd) )
      {
        // Kill the looping particle 
        //
        fParticleChange.ProposeTrackStatus( fStopAndKill )  ;
        G4int particlePDG= particleType->GetPDGEncoding();
        const G4int electronPDG= 11; // G4Electron::G4Electron()->GetPDGEncoding();
 
        // Simple statistics
        fSumEnergyKilled += endEnergy;
        fSumEnerSqKilled += endEnergy * endEnergy;
        fNumLoopersKilled++;
        
        if( endEnergy > fMaxEnergyKilled ) {
           fMaxEnergyKilled = endEnergy;
           fMaxEnergyKilledPDG = particlePDG; 
        }
        if(  particleType->GetPDGEncoding() != electronPDG )
        {
           fSumEnergyKilled_NonElectron += endEnergy;
           fSumEnerSqKilled_NonElectron += endEnergy * endEnergy;
           fNumLoopersKilled_NonElectron++;
           
           if( endEnergy > fMaxEnergyKilled_NonElectron )
           {
              fMaxEnergyKilled_NonElectron = endEnergy;
              fMaxEnergyKilled_NonElecPDG =  particlePDG;
           }
        }
 
        if( endEnergy > fThreshold_Warning_Energy && ! fSilenceLooperWarnings )
        {
          fpLogger->ReportLoopingTrack( track, stepData, fNoLooperTrials,
                                        noCallsASDI, __func__ );
        }
        fNoLooperTrials=0; 
      }
      else
      {
        fMaxEnergySaved = std::max( endEnergy, fMaxEnergySaved);
        if( fNoLooperTrials == 1 ) {
          fSumEnergySaved += endEnergy;
          if ( !stable )
             fSumEnergyUnstableSaved += endEnergy;
        }
#ifdef G4VERBOSE
        if( verboseLevel > 2 && ! fSilenceLooperWarnings )
        {
          G4cout << "   " << __func__
                 << " Particle is looping but is saved ..."  << G4endl             
                 << "   Number of trials = " << fNoLooperTrials << G4endl
                 << "   No of calls to  = " << noCallsASDI << G4endl;
        }
#endif
      }
  }
  else
  {
      fNoLooperTrials=0; 
  }
 
  // Another (sometimes better way) is to use a user-limit maximum Step size
  // to alleviate this problem .. 
 
  // Introduce smooth curved trajectories to particle-change
  //
  fParticleChange.SetPointerToVectorOfAuxiliaryPoints
    (fFieldPropagator->GimmeTrajectoryVectorAndForgetIt() );
 
  return &fParticleChange ;
}

AlongStepGetPhysicalInteractionLength()

G4double G4Transportation::AlongStepGetPhysicalInteractionLength ( const G4Track & track, G4double previousStepSize, G4double currentMinimumStep, G4double & currentSafety, G4GPILSelection * selection )

virtual

Implements G4VProcess.

Reimplemented in G4TransportationWithMsc.

Definition at line 190 of file G4Transportation.cc.

{
  // Initial actions moved to  StartTrack()
  // --------------------------------------
  // Note: in case another process changes touchable handle
  //    it will be necessary to add here (for all steps)
  // fCurrentTouchableHandle = aTrack->GetTouchableHandle();
 
  // GPILSelection is set to defaule value of CandidateForSelection
  // It is a return value
  //
  *selection = CandidateForSelection;
 
  // Get initial Energy/Momentum of the track
  //
  const G4ThreeVector startPosition    = track.GetPosition();
  const G4ThreeVector startMomentumDir = track.GetMomentumDirection();
 
  // The Step Point safety can be limited by other geometries and/or the
  // assumptions of any process - it's not always the geometrical safety.
  // We calculate the starting point's isotropic safety here.
  {
    const G4double MagSqShift = (startPosition - fPreviousSftOrigin).mag2();
 
    if(MagSqShift >= sqr(fPreviousSafety))
      currentSafety = 0.0;
    else
      currentSafety = fPreviousSafety - std::sqrt(MagSqShift);
  }
 
  // Is the particle charged or has it a magnetic moment?
  //
  const G4DynamicParticle* pParticle = track.GetDynamicParticle();
 
  const G4double particleMass   = pParticle->GetMass();
  const G4double particleCharge = pParticle->GetCharge();
  const G4double kineticEnergy = pParticle->GetKineticEnergy();
 
  const G4double magneticMoment    = pParticle->GetMagneticMoment();
  const G4ThreeVector particleSpin = pParticle->GetPolarization();
 
  // There is no need to locate the current volume. It is Done elsewhere:
  //   On track construction
  //   By the tracking, after all AlongStepDoIts, in "Relocation"
 
  // Check if the particle has a force, EM or gravitational, exerted on it
  //
 
  G4bool eligibleEM =
    (particleCharge != 0.0) || ((magneticMoment != 0.0) && fUseMagneticMoment);
  G4bool eligibleGrav = (particleMass != 0.0) && fUseGravity;
 
  fFieldExertedForce = false;
 
  if(eligibleEM || eligibleGrav)
  {
    if(G4FieldManager* fieldMgr =
         fFieldPropagator->FindAndSetFieldManager(track.GetVolume()))
    {
      // User can configure the field Manager for this track
      fieldMgr->ConfigureForTrack(&track);
      // Called here to allow a transition from no-field pointer
      // to finite field (non-zero pointer).
 
      // If the field manager has no field ptr, the field is zero
      //   by definition ( = there is no field ! )
      if(const G4Field* ptrField = fieldMgr->GetDetectorField())
        fFieldExertedForce =
          eligibleEM || (eligibleGrav && ptrField->IsGravityActive());
    }
  }
 
  G4double geometryStepLength = currentMinimumStep;
 
  if(currentMinimumStep == 0.0)
  {
    fEndPointDistance = 0.0;
    fGeometryLimitedStep = false;  //  Old code:  = (currentSafety == 0.0);
    // Changed to avoid problems when setting the step status (now done in AlongStepDoIt)
    fMomentumChanged           = false;
    fParticleIsLooping         = false;
    fEndGlobalTimeComputed     = false;
    fTransportEndPosition      = startPosition;
    fTransportEndMomentumDir   = startMomentumDir;
    fTransportEndKineticEnergy = kineticEnergy;
    fTransportEndSpin          = particleSpin;
  }
  else if(!fFieldExertedForce)
  {
    fGeometryLimitedStep = false;
    if(geometryStepLength > currentSafety || !fShortStepOptimisation)
    {
      const G4double linearStepLength = fLinearNavigator->ComputeStep(
        startPosition, startMomentumDir, currentMinimumStep, currentSafety);
 
      if(linearStepLength <= currentMinimumStep)
      {
        geometryStepLength = linearStepLength;
        fGeometryLimitedStep = true;
      }
      // Remember last safety origin & value.
      //
      fPreviousSftOrigin = startPosition;
      fPreviousSafety    = currentSafety;
      fpSafetyHelper->SetCurrentSafety(currentSafety, startPosition);
    }
 
    fEndPointDistance    = geometryStepLength;
 
    fMomentumChanged       = false;
    fParticleIsLooping     = false;
    fEndGlobalTimeComputed = false;
    fTransportEndPosition =
      startPosition + geometryStepLength * startMomentumDir;
    fTransportEndMomentumDir   = startMomentumDir;
    fTransportEndKineticEnergy = kineticEnergy;
    fTransportEndSpin          = particleSpin;
  }
  else  //  A field exerts force
  {
    const auto pParticleDef    = pParticle->GetDefinition();
    const auto particlePDGSpin = pParticleDef->GetPDGSpin();
    const auto particlePDGMagM = pParticleDef->GetPDGMagneticMoment();
 
    auto equationOfMotion = fFieldPropagator->GetCurrentEquationOfMotion();
 
    // The charge can change (dynamic), therefore the use of G4ChargeState
    //
    equationOfMotion->SetChargeMomentumMass(
      G4ChargeState(particleCharge, magneticMoment, particlePDGSpin),
      pParticle->GetTotalMomentum(), particleMass);
 
    G4FieldTrack aFieldTrack(startPosition,
                             track.GetGlobalTime(),  // Lab.
                             startMomentumDir, kineticEnergy, particleMass,
                             particleCharge, particleSpin, particlePDGMagM,
                             0.0,  // Length along track
                             particlePDGSpin);
 
    // Do the Transport in the field (non recti-linear)
    //
    const G4double lengthAlongCurve = fFieldPropagator->ComputeStep(
      aFieldTrack, currentMinimumStep, currentSafety, track.GetVolume(),
      kineticEnergy < fThreshold_Important_Energy);
 
    if(lengthAlongCurve < geometryStepLength)
      geometryStepLength = lengthAlongCurve;
 
    // Remember last safety origin & value.
    //
    fPreviousSftOrigin = startPosition;
    fPreviousSafety    = currentSafety;
    fpSafetyHelper->SetCurrentSafety(currentSafety, startPosition);
 
    fGeometryLimitedStep = fFieldPropagator->IsLastStepInVolume();
    //
    // It is possible that step was reduced in PropagatorInField due to
    // previous zero steps. To cope with case that reduced step is taken
    // in full, we must rely on PiF to obtain this value
 
    G4bool changesEnergy =
      fFieldPropagator->GetCurrentFieldManager()->DoesFieldChangeEnergy();
 
    fMomentumChanged   = true;
    fParticleIsLooping = fFieldPropagator->IsParticleLooping();
 
    fEndGlobalTimeComputed = changesEnergy;
    fTransportEndPosition    = aFieldTrack.GetPosition();
    fTransportEndMomentumDir = aFieldTrack.GetMomentumDir();
 
    // G4cout << " G4Transport: End of step pMag = " << aFieldTrack.GetMomentum().mag() << G4endl;
 
    fEndPointDistance  = (fTransportEndPosition - startPosition).mag();
 
    // Ignore change in energy for fields that conserve energy
    // This hides the integration error, but gives a better physical answer
    fTransportEndKineticEnergy =
      changesEnergy ? aFieldTrack.GetKineticEnergy() : kineticEnergy;
    fTransportEndSpin = aFieldTrack.GetSpin();
 
    if(fEndGlobalTimeComputed)
    {
      // If the field can change energy, then the time must be integrated
      //    - so this should have been updated
      //
      fCandidateEndGlobalTime = aFieldTrack.GetLabTimeOfFlight();
 
      // was ( fCandidateEndGlobalTime != track.GetGlobalTime() );
      // a cleaner way is to have FieldTrack knowing whether time is updated.
    }
#if defined(G4VERBOSE) || defined(G4DEBUG_TRANSPORT)
    else
    {
      // The energy should be unchanged by field transport,
      //    - so the time changed will be calculated elsewhere
      //
      // Check that the integration preserved the energy
      //     -  and if not correct this!
      G4double startEnergy = kineticEnergy;
      G4double endEnergy   = fTransportEndKineticEnergy;
 
      static G4ThreadLocal G4int no_large_ediff;
      if(verboseLevel > 1)
      {
        if(std::fabs(startEnergy - endEnergy) > perThousand * endEnergy)
        {
          static G4ThreadLocal G4int no_warnings = 0, warnModulo = 1,
                                     moduloFactor = 10;
          no_large_ediff++;
          if((no_large_ediff % warnModulo) == 0)
          {
            no_warnings++;
            std::ostringstream message;
            message << "Energy change in Step is above 1^-3 relative value. "
                    << G4endl << "     Relative change in 'tracking' step = "
                    << std::setw(15) << (endEnergy - startEnergy) / startEnergy
                    << G4endl << "     Starting E= " << std::setw(12)
                    << startEnergy / MeV << " MeV " << G4endl
                    << "     Ending   E= " << std::setw(12) << endEnergy / MeV
                    << " MeV " << G4endl
                    << "Energy has been corrected -- however, review"
                    << " field propagation parameters for accuracy." << G4endl;
            if((verboseLevel > 2) || (no_warnings < 4) ||
               (no_large_ediff == warnModulo * moduloFactor))
            {
              message << "These include EpsilonStepMax(/Min) in G4FieldManager "
                      << G4endl
                      << "which determine fractional error per step for "
                         "integrated quantities. "
                      << G4endl
                      << "Note also the influence of the permitted number of "
                         "integration steps."
                      << G4endl;
            }
            message << "Bad 'endpoint'. Energy change detected and corrected."
                    << G4endl << "Has occurred already " << no_large_ediff
                    << " times.";
            G4Exception("G4Transportation::AlongStepGetPIL()", "EnergyChange",
                        JustWarning, message);
            if(no_large_ediff == warnModulo * moduloFactor)
            {
              warnModulo *= moduloFactor;
            }
          }
        }
      }  // end of if (verboseLevel)
    }
#endif
  }
 
  // Update the safety starting from the end-point,
  // if it will become negative at the end-point.
  //
  if(currentSafety < fEndPointDistance)
  {
    if(particleCharge != 0.0)
    {
      G4double endSafety =
        fLinearNavigator->ComputeSafety(fTransportEndPosition);
      currentSafety      = endSafety;
      fPreviousSftOrigin = fTransportEndPosition;
      fPreviousSafety    = currentSafety;
      fpSafetyHelper->SetCurrentSafety(currentSafety, fTransportEndPosition);
 
      // Because the Stepping Manager assumes it is from the start point,
      //  add the StepLength
      //
      currentSafety += fEndPointDistance;
 
#ifdef G4DEBUG_TRANSPORT
      G4cout.precision(12);
      G4cout << "***G4Transportation::AlongStepGPIL ** " << G4endl;
      G4cout << "  Called Navigator->ComputeSafety at " << fTransportEndPosition
             << "    and it returned safety=  " << endSafety << G4endl;
      G4cout << "  Adding endpoint distance   " << fEndPointDistance
             << "    to obtain pseudo-safety= " << currentSafety << G4endl;
    }
    else
    {
      G4cout << "***G4Transportation::AlongStepGPIL ** " << G4endl;
      G4cout << "  Avoiding call to ComputeSafety : " << G4endl;
      G4cout << "    charge     = " << particleCharge << G4endl;
      G4cout << "    mag moment = " << magneticMoment << G4endl;
#endif
    }
  }
 
  fFirstStepInVolume = fNewTrack || fLastStepInVolume;
  fLastStepInVolume  = false;
  fNewTrack          = false;
 
  fParticleChange.ProposeFirstStepInVolume(fFirstStepInVolume);
  fParticleChange.ProposeTrueStepLength(geometryStepLength);
 
  return geometryStepLength;
}

Referenced by G4TransportationWithMsc::AlongStepGetPhysicalInteractionLength().

AtRestDoIt()

G4VParticleChange * G4Transportation::AtRestDoIt ( const G4Track & , const G4Step &  )

inlinevirtual

Implements G4VProcess.

Definition at line 144 of file G4Transportation.hh.

       { return 0; }     // No operation in AtRestDoIt

AtRestGetPhysicalInteractionLength()

G4double G4Transportation::AtRestGetPhysicalInteractionLength ( const G4Track & , G4ForceCondition *  )

inlinevirtual

Implements G4VProcess.

Definition at line 140 of file G4Transportation.hh.

       { return -1.0; }  // No operation in AtRestGPIL

EnableGravity()

G4bool G4Transportation::EnableGravity ( G4bool useGravity = true )

static

Definition at line 845 of file G4Transportation.cc.

{
  G4bool lastValue= fUseGravity;
  fUseGravity= useGravity;
  return lastValue;
}

EnableMagneticMoment()

G4bool G4Transportation::EnableMagneticMoment ( G4bool useMoment = true )

static

Definition at line 835 of file G4Transportation.cc.

{
  G4bool lastValue= fUseMagneticMoment;
  fUseMagneticMoment= useMoment;
  return lastValue;
}

Referenced by G4CoupledTransportation::EnableUseMagneticMoment(), and EnableUseMagneticMoment().

EnableShortStepOptimisation()

void G4Transportation::EnableShortStepOptimisation ( G4bool optimise = true )

inline

EnableUseMagneticMoment()

static G4bool G4Transportation::EnableUseMagneticMoment ( G4bool useMoment = true )

inlinestatic

Definition at line 135 of file G4Transportation.hh.

     { return EnableMagneticMoment(useMoment); }  // Old name - will be deprecated

FieldExertedForce()

G4bool G4Transportation::FieldExertedForce ( )

inline

Definition at line 93 of file G4Transportation.hh.

{ return fFieldExertedForce; }

GetMaxEnergyKilled()

G4double G4Transportation::GetMaxEnergyKilled ( ) const

inline

GetPropagatorInField()

G4PropagatorInField * G4Transportation::GetPropagatorInField

(

)

GetSilenceLooperWarnings()

G4bool G4Transportation::GetSilenceLooperWarnings ( )

static

Definition at line 863 of file G4Transportation.cc.

{
  return fSilenceLooperWarnings; 
}

GetSumEnergyKilled()

G4double G4Transportation::GetSumEnergyKilled ( ) const

inline

GetThresholdImportantEnergy()

G4double G4Transportation::GetThresholdImportantEnergy ( ) const

inline

GetThresholdTrials()

G4int G4Transportation::GetThresholdTrials ( ) const

inline

GetThresholdWarningEnergy()

G4double G4Transportation::GetThresholdWarningEnergy ( ) const

inline

PostStepDoIt()

G4VParticleChange * G4Transportation::PostStepDoIt ( const G4Track & track, const G4Step & stepData )

virtual

Implements G4VProcess.

Definition at line 708 of file G4Transportation.cc.

{
   G4TouchableHandle retCurrentTouchable ;   // The one to return
   G4bool isLastStep= false; 
 
  // Initialize ParticleChange  (by setting all its members equal
  //                             to corresponding members in G4Track)
  // fParticleChange.Initialize(track) ;  // To initialise TouchableChange
 
  fParticleChange.ProposeTrackStatus(track.GetTrackStatus()) ;
 
  // If the Step was determined by the volume boundary,
  // logically relocate the particle
 
  if(fGeometryLimitedStep)
  {  
    // fCurrentTouchable will now become the previous touchable, 
    // and what was the previous will be freed.
    // (Needed because the preStepPoint can point to the previous touchable)
 
    fLinearNavigator->SetGeometricallyLimitedStep() ;
    fLinearNavigator->
    LocateGlobalPointAndUpdateTouchableHandle( track.GetPosition(),
                                               track.GetMomentumDirection(),
                                               fCurrentTouchableHandle,
                                               true                      ) ;
    // Check whether the particle is out of the world volume 
    // If so it has exited and must be killed.
    //
    if( fCurrentTouchableHandle->GetVolume() == 0 )
    {
       fParticleChange.ProposeTrackStatus( fStopAndKill ) ;
    }
    retCurrentTouchable = fCurrentTouchableHandle ;
    fParticleChange.SetTouchableHandle( fCurrentTouchableHandle ) ;
 
    // Update the Step flag which identifies the Last Step in a volume
    if( !fFieldExertedForce )
       isLastStep =  fLinearNavigator->ExitedMotherVolume() 
                  || fLinearNavigator->EnteredDaughterVolume() ;
    else
       isLastStep = fFieldPropagator->IsLastStepInVolume(); 
  }
  else                 // fGeometryLimitedStep  is false
  {                    
    // This serves only to move the Navigator's location
    //
    fLinearNavigator->LocateGlobalPointWithinVolume( track.GetPosition() ) ;
 
    // The value of the track's current Touchable is retained. 
    // (and it must be correct because we must use it below to
    // overwrite the (unset) one in particle change)
    //  It must be fCurrentTouchable too ??
    //
    fParticleChange.SetTouchableHandle( track.GetTouchableHandle() ) ;
    retCurrentTouchable = track.GetTouchableHandle() ;
 
    isLastStep= false;
  }         // endif ( fGeometryLimitedStep ) 
  fLastStepInVolume= isLastStep; 
  
  fParticleChange.ProposeFirstStepInVolume(fFirstStepInVolume);
  fParticleChange.ProposeLastStepInVolume(isLastStep);    
 
  SetTouchableInformation(retCurrentTouchable);
 
  return &fParticleChange ;
}

PostStepGetPhysicalInteractionLength()

G4double G4Transportation::PostStepGetPhysicalInteractionLength ( const G4Track & , G4double previousStepSize, G4ForceCondition * pForceCond )

virtual

Implements G4VProcess.

Definition at line 651 of file G4Transportation.cc.

{
  fFieldExertedForce = false; // Not known
  *pForceCond = Forced ; 
  return DBL_MAX ;  // was kInfinity ; but convention now is DBL_MAX
}

PrintStatistics()

void G4Transportation::PrintStatistics

(

std::ostream &

outStr

)

const

Definition at line 148 of file G4Transportation.cc.

{
   outStr << " G4Transportation: Statistics for looping particles " << G4endl;   
   if( fSumEnergyKilled > 0.0 || fNumLoopersKilled > 0 )
   {
      outStr << "   Sum of energy of looping tracks killed: "
             <<  fSumEnergyKilled / CLHEP::MeV << " MeV "          
             << " from " << fNumLoopersKilled << "  tracks "  << G4endl
             <<  "  Sum of energy of non-electrons        : "
             << fSumEnergyKilled_NonElectron / CLHEP::MeV << " MeV "
             << "  from " << fNumLoopersKilled_NonElectron << " tracks "
             << G4endl;
      outStr << "   Max energy of  *any type*  looper killed: " << fMaxEnergyKilled
             << "    its PDG was " << fMaxEnergyKilledPDG  << G4endl;
      if( fMaxEnergyKilled_NonElectron > 0.0 )
      {
         outStr << "   Max energy of non-electron looper killed: " 
                << fMaxEnergyKilled_NonElectron
                << "    its PDG was " << fMaxEnergyKilled_NonElecPDG << G4endl;
      }
      if( fMaxEnergySaved > 0.0 )
      {
         outStr << "   Max energy of loopers 'saved':  " << fMaxEnergySaved << G4endl;
         outStr << "   Sum of energy of loopers 'saved': "
                <<  fSumEnergySaved << G4endl;
         outStr << "   Sum of energy of unstable loopers 'saved': "
                << fSumEnergyUnstableSaved << G4endl;
      }
   }
   else
   {
      outStr << " No looping tracks found or killed. " << G4endl;
   }
}

Referenced by ~G4Transportation().

ProcessDescription()

void G4Transportation::ProcessDescription ( std::ostream & outFile ) const

virtual

Reimplemented from G4VProcess.

Definition at line 922 of file G4Transportation.cc.

{
  G4String indent = "  "; //  : "");
  G4long oldPrec= outStr.precision(6);
  // outStr << std::setprecision(6);
  outStr << G4endl << indent << GetProcessName() << ": ";
 
  outStr << "   Parameters for looping particles: " << G4endl
         << "     warning-E = " << fThreshold_Warning_Energy / CLHEP::MeV  << " MeV "  << G4endl
         << "     important E = " << fThreshold_Important_Energy / CLHEP::MeV << " MeV " << G4endl
         << "     thresholdTrials " << fThresholdTrials << G4endl;
  outStr.precision(oldPrec);
}

PushThresholdsToLogger()

void G4Transportation::PushThresholdsToLogger

(

)

Referenced by G4Transportation(), ReportLooperThresholds(), SetHighLooperThresholds(), and SetLowLooperThresholds().

ReportLooperThresholds()

void G4Transportation::ReportLooperThresholds

(

)

Definition at line 914 of file G4Transportation.cc.

{
   PushThresholdsToLogger();  // To be absolutely certain they are in sync   
   fpLogger->ReportLooperThresholds("G4Transportation");
}

Referenced by SetHighLooperThresholds(), and SetLowLooperThresholds().

ReportMissingLogger()

void G4Transportation::ReportMissingLogger ( const char * methodName )

protected

Definition at line 903 of file G4Transportation.cc.

{
   const char* message= "Logger object missing from G4Transportation object";
   G4String classAndMethod= G4String("G4Transportation") + G4String( methodName );
   G4Exception(classAndMethod, "Missing Logger", JustWarning, message);   
}

ResetKilledStatistics()

void G4Transportation::ResetKilledStatistics ( G4int report = 1 )

inline

SetHighLooperThresholds()

void G4Transportation::SetHighLooperThresholds

(

)

Definition at line 871 of file G4Transportation.cc.

{
  // Restores the old high values -- potentially appropriate for energy-frontier
  //   HEP experiments.
  // Caution:  All tracks with E < 100 MeV that are found to loop are 
  SetThresholdWarningEnergy(    100.0 * CLHEP::MeV ); // Warn above this energy
  SetThresholdImportantEnergy(  250.0 * CLHEP::MeV ); // Extra trial above this En
 
  G4int maxTrials = 10;
  SetThresholdTrials( maxTrials );
  
  PushThresholdsToLogger();  // Again, to be sure
  if( verboseLevel )  ReportLooperThresholds();    
}

Referenced by G4Transportation().

SetLowLooperThresholds()

void G4Transportation::SetLowLooperThresholds

(

)

Definition at line 887 of file G4Transportation.cc.

{
  // These values were the default in Geant4 10.5 - beta
  SetThresholdWarningEnergy(     1.0 * CLHEP::keV ); // Warn above this En
  SetThresholdImportantEnergy(   1.0 * CLHEP::MeV ); // Extra trials above it
 
  G4int maxTrials = 30; //  A new value - was 10
  SetThresholdTrials( maxTrials );
 
  PushThresholdsToLogger();  // Again, to be sure
  if( verboseLevel )  ReportLooperThresholds();  
}

SetPropagatorInField()

void G4Transportation::SetPropagatorInField

(

G4PropagatorInField *

pFieldPropagator

)

SetSilenceLooperWarnings()

void G4Transportation::SetSilenceLooperWarnings ( G4bool val )

static

Definition at line 856 of file G4Transportation.cc.

{
  fSilenceLooperWarnings= val;  // Flag to *Supress* all 'looper' warnings  
}

SetThresholdImportantEnergy()

void G4Transportation::SetThresholdImportantEnergy ( G4double newEnImp )

inline

Referenced by G4Transportation(), SetHighLooperThresholds(), and SetLowLooperThresholds().

SetThresholdTrials()

void G4Transportation::SetThresholdTrials ( G4int newMaxTrials )

inline

Referenced by G4Transportation(), SetHighLooperThresholds(), and SetLowLooperThresholds().

SetThresholdWarningEnergy()

void G4Transportation::SetThresholdWarningEnergy ( G4double newEnWarn )

inline

Referenced by G4Transportation(), SetHighLooperThresholds(), and SetLowLooperThresholds().

SetTouchableInformation()

void G4Transportation::SetTouchableInformation ( const G4TouchableHandle & touchable )

protected

Definition at line 663 of file G4Transportation.cc.

{
  const G4VPhysicalVolume* pNewVol = touchable->GetVolume() ;
  const G4Material* pNewMaterial   = 0 ;
  G4VSensitiveDetector* pNewSensitiveDetector = 0;
 
  if( pNewVol != 0 )
  {
    pNewMaterial= pNewVol->GetLogicalVolume()->GetMaterial();
    pNewSensitiveDetector= pNewVol->GetLogicalVolume()->GetSensitiveDetector();
  }
 
  fParticleChange.SetMaterialInTouchable( (G4Material *) pNewMaterial ) ;
  fParticleChange.SetSensitiveDetectorInTouchable( pNewSensitiveDetector ) ;
  // temporarily until Get/Set Material of ParticleChange,
  // and StepPoint can be made const.
 
  const G4MaterialCutsCouple* pNewMaterialCutsCouple = 0;
  if( pNewVol != 0 )
  {
    pNewMaterialCutsCouple=pNewVol->GetLogicalVolume()->GetMaterialCutsCouple();
  }
 
  if ( pNewVol!=0 && pNewMaterialCutsCouple!=0
    && pNewMaterialCutsCouple->GetMaterial()!=pNewMaterial )
  {
    // for parametrized volume
    //
    pNewMaterialCutsCouple =
      G4ProductionCutsTable::GetProductionCutsTable()
                             ->GetMaterialCutsCouple(pNewMaterial,
                               pNewMaterialCutsCouple->GetProductionCuts());
  }
  fParticleChange.SetMaterialCutsCoupleInTouchable( pNewMaterialCutsCouple );
 
  // Set the touchable in ParticleChange
  // this must always be done because the particle change always
  // uses this value to overwrite the current touchable pointer.
  //
  fParticleChange.SetTouchableHandle(touchable) ;
}

Referenced by G4CoupledTransportation::PostStepDoIt(), and PostStepDoIt().

StartTracking()

void G4Transportation::StartTracking ( G4Track * aTrack )

virtual

Reimplemented from G4VProcess.

Reimplemented in G4TransportationWithMsc.

Definition at line 785 of file G4Transportation.cc.

{
  G4VProcess::StartTracking(aTrack);
  fNewTrack= true;
  fFirstStepInVolume= true;
  fLastStepInVolume= false;
  
  // The actions here are those that were taken in AlongStepGPIL
  // when track.GetCurrentStepNumber()==1
 
  // Whether field exists should be determined at run level -- TODO
  G4FieldManagerStore* fieldMgrStore= G4FieldManagerStore::GetInstance();
  fAnyFieldExists = fieldMgrStore->size() > 0;
  
  // reset safety value and center
  //
  fPreviousSafety    = 0.0 ; 
  fPreviousSftOrigin = G4ThreeVector(0.,0.,0.) ;
  
  // reset looping counter -- for motion in field
  fNoLooperTrials= 0; 
  // Must clear this state .. else it depends on last track's value
  //  --> a better solution would set this from state of suspended track TODO ? 
  // Was if( aTrack->GetCurrentStepNumber()==1 ) { .. }
 
  // ChordFinder reset internal state
  //
  if( fFieldPropagator && fAnyFieldExists )     
  {
     fFieldPropagator->ClearPropagatorState();   
       // Resets all state of field propagator class (ONLY) including safety
       // values (in case of overlaps and to wipe for first track).
  }
 
  // Make sure to clear the chord finders of all fields (i.e. managers)
  //
  fieldMgrStore->ClearAllChordFindersState(); 
 
  // Update the current touchable handle  (from the track's)
  //
  fCurrentTouchableHandle = aTrack->GetTouchableHandle();
 
  // Inform field propagator of new track
  //
  fFieldPropagator->PrepareNewTrack();
}

Referenced by G4CoupledTransportation::StartTracking(), and G4TransportationWithMsc::StartTracking().

Member Data Documentation

fAbandonUnstableTrials

G4int G4Transportation::fAbandonUnstableTrials = 0

protected

Definition at line 206 of file G4Transportation.hh.

Referenced by AlongStepDoIt().

fAnyFieldExists

G4bool G4Transportation::fAnyFieldExists = false

protected

Definition at line 176 of file G4Transportation.hh.

Referenced by StartTracking().

fCandidateEndGlobalTime

G4double G4Transportation::fCandidateEndGlobalTime = 0.0

protected

Definition at line 173 of file G4Transportation.hh.

Referenced by AlongStepDoIt(), G4CoupledTransportation::AlongStepGetPhysicalInteractionLength(), and AlongStepGetPhysicalInteractionLength().

fCurrentTouchableHandle

G4TouchableHandle G4Transportation::fCurrentTouchableHandle

protected

Definition at line 188 of file G4Transportation.hh.

Referenced by G4Transportation(), G4CoupledTransportation::PostStepDoIt(), PostStepDoIt(), and StartTracking().

fEndGlobalTimeComputed

G4bool G4Transportation::fEndGlobalTimeComputed = false

protected

Definition at line 172 of file G4Transportation.hh.

Referenced by AlongStepDoIt(), G4CoupledTransportation::AlongStepGetPhysicalInteractionLength(), and AlongStepGetPhysicalInteractionLength().

fEndPointDistance

G4double G4Transportation::fEndPointDistance

protected

Definition at line 197 of file G4Transportation.hh.

Referenced by G4CoupledTransportation::AlongStepGetPhysicalInteractionLength(), AlongStepGetPhysicalInteractionLength(), and G4TransportationWithMsc::AlongStepGetPhysicalInteractionLength().

fFieldExertedForce

G4bool G4Transportation::fFieldExertedForce = false

protected

Definition at line 186 of file G4Transportation.hh.

Referenced by AlongStepGetPhysicalInteractionLength(), FieldExertedForce(), PostStepDoIt(), and PostStepGetPhysicalInteractionLength().

fFieldPropagator

G4PropagatorInField* G4Transportation::fFieldPropagator

protected

Definition at line 164 of file G4Transportation.hh.

Referenced by AlongStepDoIt(), G4CoupledTransportation::AlongStepGetPhysicalInteractionLength(), AlongStepGetPhysicalInteractionLength(), G4Transportation(), PostStepDoIt(), and StartTracking().

fFirstStepInVolume

G4bool G4Transportation::fFirstStepInVolume = true

protected

Definition at line 180 of file G4Transportation.hh.

Referenced by G4CoupledTransportation::AlongStepGetPhysicalInteractionLength(), AlongStepGetPhysicalInteractionLength(), G4CoupledTransportation::IsFirstStepInAnyVolume(), G4CoupledTransportation::PostStepDoIt(), PostStepDoIt(), and StartTracking().

fGeometryLimitedStep

G4bool G4Transportation::fGeometryLimitedStep = true

protected

Definition at line 183 of file G4Transportation.hh.

Referenced by AlongStepDoIt(), G4CoupledTransportation::AlongStepGetPhysicalInteractionLength(), AlongStepGetPhysicalInteractionLength(), G4CoupledTransportation::IsLastStepInAnyVolume(), G4CoupledTransportation::PostStepDoIt(), and PostStepDoIt().

fLastStepInVolume

G4bool G4Transportation::fLastStepInVolume = false

protected

Definition at line 181 of file G4Transportation.hh.

Referenced by AlongStepGetPhysicalInteractionLength(), PostStepDoIt(), and StartTracking().

fLinearNavigator

G4Navigator* G4Transportation::fLinearNavigator

protected

Definition at line 161 of file G4Transportation.hh.

Referenced by AlongStepGetPhysicalInteractionLength(), G4TransportationWithMsc::AlongStepGetPhysicalInteractionLength(), G4Transportation(), and PostStepDoIt().

fMaxEnergyKilled

G4double G4Transportation::fMaxEnergyKilled = -1.0

protected

Definition at line 216 of file G4Transportation.hh.

Referenced by AlongStepDoIt(), and PrintStatistics().

fMaxEnergyKilled_NonElecPDG

G4int G4Transportation::fMaxEnergyKilled_NonElecPDG = 0

protected

Definition at line 222 of file G4Transportation.hh.

Referenced by AlongStepDoIt(), and PrintStatistics().

fMaxEnergyKilled_NonElectron

G4double G4Transportation::fMaxEnergyKilled_NonElectron = -1.0

protected

Definition at line 221 of file G4Transportation.hh.

Referenced by AlongStepDoIt(), and PrintStatistics().

fMaxEnergyKilledPDG

G4int G4Transportation::fMaxEnergyKilledPDG = 0

protected

Definition at line 217 of file G4Transportation.hh.

Referenced by AlongStepDoIt(), and PrintStatistics().

fMaxEnergySaved

G4double G4Transportation::fMaxEnergySaved = -1.0

protected

Definition at line 225 of file G4Transportation.hh.

Referenced by AlongStepDoIt(), and PrintStatistics().

fMomentumChanged

G4bool G4Transportation::fMomentumChanged = true

protected

Definition at line 171 of file G4Transportation.hh.

Referenced by AlongStepDoIt(), G4CoupledTransportation::AlongStepGetPhysicalInteractionLength(), AlongStepGetPhysicalInteractionLength(), and G4TransportationWithMsc::AlongStepGetPhysicalInteractionLength().

fNewTrack

G4bool G4Transportation::fNewTrack = true

protected

Definition at line 179 of file G4Transportation.hh.

Referenced by G4CoupledTransportation::AlongStepGetPhysicalInteractionLength(), AlongStepGetPhysicalInteractionLength(), and StartTracking().

fNoLooperTrials

G4int G4Transportation::fNoLooperTrials = 0

protected

Definition at line 210 of file G4Transportation.hh.

Referenced by AlongStepDoIt(), and StartTracking().

fNumLoopersKilled

unsigned long G4Transportation::fNumLoopersKilled = 0

protected

Definition at line 218 of file G4Transportation.hh.

Referenced by AlongStepDoIt(), and PrintStatistics().

fNumLoopersKilled_NonElectron

unsigned long G4Transportation::fNumLoopersKilled_NonElectron = 0

protected

Definition at line 223 of file G4Transportation.hh.

Referenced by AlongStepDoIt(), and PrintStatistics().

fParticleChange

G4ParticleChangeForTransport G4Transportation::fParticleChange

protected

Definition at line 194 of file G4Transportation.hh.

Referenced by AlongStepDoIt(), G4CoupledTransportation::AlongStepGetPhysicalInteractionLength(), AlongStepGetPhysicalInteractionLength(), G4TransportationWithMsc::AlongStepGetPhysicalInteractionLength(), G4Transportation(), G4CoupledTransportation::PostStepDoIt(), PostStepDoIt(), and SetTouchableInformation().

fParticleIsLooping

G4bool G4Transportation::fParticleIsLooping = false

protected

Definition at line 178 of file G4Transportation.hh.

Referenced by AlongStepDoIt(), G4CoupledTransportation::AlongStepGetPhysicalInteractionLength(), and AlongStepGetPhysicalInteractionLength().

fpLogger

G4TransportationLogger* G4Transportation::fpLogger

protected

Definition at line 233 of file G4Transportation.hh.

Referenced by AlongStepDoIt(), G4Transportation(), ReportLooperThresholds(), and ~G4Transportation().

fPreviousSafety

G4double G4Transportation::fPreviousSafety

protected

Definition at line 191 of file G4Transportation.hh.

Referenced by AlongStepGetPhysicalInteractionLength(), and StartTracking().

fPreviousSftOrigin

G4ThreeVector G4Transportation::fPreviousSftOrigin

protected

Definition at line 190 of file G4Transportation.hh.

Referenced by G4CoupledTransportation::AlongStepGetPhysicalInteractionLength(), AlongStepGetPhysicalInteractionLength(), G4CoupledTransportation::StartTracking(), and StartTracking().

fpSafetyHelper

G4SafetyHelper* G4Transportation::fpSafetyHelper

protected

Definition at line 232 of file G4Transportation.hh.

Referenced by G4CoupledTransportation::AlongStepGetPhysicalInteractionLength(), AlongStepGetPhysicalInteractionLength(), G4TransportationWithMsc::AlongStepGetPhysicalInteractionLength(), and G4Transportation().

fShortStepOptimisation

G4bool G4Transportation::fShortStepOptimisation

protected

Definition at line 230 of file G4Transportation.hh.

Referenced by AlongStepGetPhysicalInteractionLength(), and G4Transportation().

fSilenceLooperWarnings

G4bool G4Transportation::fSilenceLooperWarnings = false

staticprotected

Definition at line 239 of file G4Transportation.hh.

Referenced by AlongStepDoIt(), G4Transportation(), GetSilenceLooperWarnings(), and SetSilenceLooperWarnings().

fSumEnergyKilled

G4double G4Transportation::fSumEnergyKilled = 0.0

protected

Definition at line 214 of file G4Transportation.hh.

Referenced by AlongStepDoIt(), PrintStatistics(), and ~G4Transportation().

fSumEnergyKilled_NonElectron

G4double G4Transportation::fSumEnergyKilled_NonElectron = 0.0

protected

Definition at line 219 of file G4Transportation.hh.

Referenced by AlongStepDoIt(), and PrintStatistics().

fSumEnergySaved

G4double G4Transportation::fSumEnergySaved = 0.0

protected

Definition at line 224 of file G4Transportation.hh.

Referenced by AlongStepDoIt(), and PrintStatistics().

fSumEnergyUnstableSaved

G4double G4Transportation::fSumEnergyUnstableSaved = 0.0

protected

Definition at line 226 of file G4Transportation.hh.

Referenced by AlongStepDoIt(), and PrintStatistics().

fSumEnerSqKilled

G4double G4Transportation::fSumEnerSqKilled = 0.0

protected

Definition at line 215 of file G4Transportation.hh.

Referenced by AlongStepDoIt().

fSumEnerSqKilled_NonElectron

G4double G4Transportation::fSumEnerSqKilled_NonElectron = 0.0

protected

Definition at line 220 of file G4Transportation.hh.

Referenced by AlongStepDoIt().

fThreshold_Important_Energy

G4double G4Transportation::fThreshold_Important_Energy = 1.0 * CLHEP::MeV

protected

Definition at line 202 of file G4Transportation.hh.

Referenced by AlongStepDoIt(), and AlongStepGetPhysicalInteractionLength().

fThreshold_Warning_Energy

G4double G4Transportation::fThreshold_Warning_Energy = 1.0 * CLHEP::keV

protected

Definition at line 201 of file G4Transportation.hh.

Referenced by AlongStepDoIt().

fThresholdTrials

G4int G4Transportation::fThresholdTrials = 10

protected

Definition at line 203 of file G4Transportation.hh.

Referenced by AlongStepDoIt().

fTransportEndKineticEnergy

G4double G4Transportation::fTransportEndKineticEnergy = 0.0

protected

Definition at line 169 of file G4Transportation.hh.

Referenced by AlongStepDoIt(), G4CoupledTransportation::AlongStepGetPhysicalInteractionLength(), AlongStepGetPhysicalInteractionLength(), and G4TransportationWithMsc::AlongStepGetPhysicalInteractionLength().

fTransportEndMomentumDir

G4ThreeVector G4Transportation::fTransportEndMomentumDir = G4ThreeVector( 0.0, 0.0, 0.0 )

protected

Definition at line 168 of file G4Transportation.hh.

Referenced by AlongStepDoIt(), G4CoupledTransportation::AlongStepGetPhysicalInteractionLength(), AlongStepGetPhysicalInteractionLength(), and G4TransportationWithMsc::AlongStepGetPhysicalInteractionLength().

fTransportEndPosition

G4ThreeVector G4Transportation::fTransportEndPosition = G4ThreeVector( 0.0, 0.0, 0.0 )

protected

Definition at line 167 of file G4Transportation.hh.

Referenced by AlongStepDoIt(), G4CoupledTransportation::AlongStepGetPhysicalInteractionLength(), AlongStepGetPhysicalInteractionLength(), G4TransportationWithMsc::AlongStepGetPhysicalInteractionLength(), and G4CoupledTransportation::PostStepDoIt().

fTransportEndSpin

G4ThreeVector G4Transportation::fTransportEndSpin = G4ThreeVector( 0.0, 0.0, 0.0 )

protected

Definition at line 170 of file G4Transportation.hh.

Referenced by AlongStepDoIt(), G4CoupledTransportation::AlongStepGetPhysicalInteractionLength(), and AlongStepGetPhysicalInteractionLength().

fUseGravity

G4bool G4Transportation::fUseGravity = false

staticprotected

Definition at line 238 of file G4Transportation.hh.

Referenced by G4CoupledTransportation::AlongStepGetPhysicalInteractionLength(), AlongStepGetPhysicalInteractionLength(), and EnableGravity().

fUseMagneticMoment

G4bool G4Transportation::fUseMagneticMoment =false

staticprotected

Definition at line 237 of file G4Transportation.hh.

Referenced by G4CoupledTransportation::AlongStepGetPhysicalInteractionLength(), AlongStepGetPhysicalInteractionLength(), and EnableMagneticMoment().

---

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

• G4Transportation.hh
• G4Transportation.cc