#include <G4MonopoleTransportation.hh>

Inheritance diagram for G4MonopoleTransportation:

Public Member Functions
	G4MonopoleTransportation (const G4Monopole *p, G4int verbosityLevel=1)

	~G4MonopoleTransportation ()

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

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

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

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

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)

G4double	GetMaxEnergyKilled () const

G4double	GetSumEnergyKilled () const

void	ResetKilledStatistics (G4int report=1)

void	EnableShortStepOptimisation (G4bool optimise=true)

virtual G4double	AtRestGetPhysicalInteractionLength (const G4Track &, G4ForceCondition *)

virtual G4VParticleChange *	AtRestDoIt (const G4Track &, const G4Step &)

virtual void	StartTracking (G4Track *aTrack)

Public Member Functions inherited from G4VProcess
	G4VProcess (const G4String &aName="NoName", G4ProcessType aType=fNotDefined)

	G4VProcess (const G4VProcess &right)

virtual	~G4VProcess ()

G4int	operator== (const G4VProcess &right) const

G4int	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 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 &)

Protected Member Functions
G4bool	DoesGlobalFieldExist ()

Protected Member Functions inherited from G4VProcess
void	SubtractNumberOfInteractionLengthLeft (G4double previousStepSize)

void	ClearNumberOfInteractionLengthLeft ()

Additional Inherited Members
Static Public Member Functions inherited from G4VProcess
static const G4String &	GetProcessTypeName (G4ProcessType)

Protected Attributes inherited from G4VProcess
const G4ProcessManager *	aProcessManager

G4VParticleChange *	pParticleChange

G4ParticleChange	aParticleChange

G4double	theNumberOfInteractionLengthLeft

G4double	currentInteractionLength

G4double	theInitialNumberOfInteractionLength

G4String	theProcessName

G4String	thePhysicsTableFileName

G4ProcessType	theProcessType

G4int	theProcessSubType

G4double	thePILfactor

G4bool	enableAtRestDoIt

G4bool	enableAlongStepDoIt

G4bool	enablePostStepDoIt

G4int	verboseLevel

Detailed Description

Definition at line 66 of file G4MonopoleTransportation.hh.

Constructor & Destructor Documentation

G4MonopoleTransportation::G4MonopoleTransportation	(	const G4Monopole *	p,
		G4int	verbosityLevel = `1`
	)

Definition at line 60 of file G4MonopoleTransportation.cc.

References G4MonopoleFieldSetup::GetMonopoleFieldSetup(), G4TransportationManager::GetNavigatorForTracking(), G4TransportationManager::GetPropagatorInField(), G4TransportationManager::GetSafetyHelper(), G4TransportationManager::GetTransportationManager(), G4VProcess::SetProcessSubType(), TRANSPORTATION, and G4VProcess::verboseLevel.

   : G4VProcess( G4String("MonopoleTransportation"), fTransportation ),
     fParticleDef(mpl),
     fMagSetup(0),
     fLinearNavigator(0),
     fFieldPropagator(0),
     fParticleIsLooping( false ),
     fPreviousSftOrigin (0.,0.,0.),
     fPreviousSafety    ( 0.0 ),
     fThreshold_Warning_Energy( 100 * MeV ),  
     fThreshold_Important_Energy( 250 * MeV ), 
     fThresholdTrials( 10 ), 
     // fUnimportant_Energy( 1 * MeV ), 
     fNoLooperTrials(0),
     fSumEnergyKilled( 0.0 ), fMaxEnergyKilled( 0.0 ), 
     fShortStepOptimisation(false),    // Old default: true (=fast short steps)
     fpSafetyHelper(0)
 {
   verboseLevel = verb;
 
   // set Process Sub Type
   SetProcessSubType(TRANSPORTATION);  
 
   fMagSetup = G4MonopoleFieldSetup::GetMonopoleFieldSetup();
   
   G4TransportationManager* transportMgr ; 
 
   transportMgr = G4TransportationManager::GetTransportationManager() ; 
 
   fLinearNavigator = transportMgr->GetNavigatorForTracking() ; 
 
   // fGlobalFieldMgr = transportMgr->GetFieldManager() ;
 
   fFieldPropagator = transportMgr->GetPropagatorInField() ;
 
   fpSafetyHelper =   transportMgr->GetSafetyHelper();  // New 
 
   // Cannot determine whether a field exists here,
   //  because it would only work if the field manager has informed 
   //  about the detector's field before this transportation process 
   //  is constructed.
   // Instead later the method DoesGlobalFieldExist() is called
 
   static G4TouchableHandle nullTouchableHandle;  // Points to (G4VTouchable*) 0
   fCurrentTouchableHandle = nullTouchableHandle; 
 
   fEndGlobalTimeComputed  = false;
   fCandidateEndGlobalTime = 0;
 }

G4MonopoleTransportation::~G4MonopoleTransportation ( )

Definition at line 113 of file G4MonopoleTransportation.cc.

References G4cout, G4endl, and G4VProcess::verboseLevel.

 {
   if( (verboseLevel > 0) && (fSumEnergyKilled > 0.0 ) ){ 
     G4cout << " G4MonopoleTransportation: Statistics for looping particles " 
            << G4endl;
     G4cout << "   Sum of energy of loopers killed: " <<  fSumEnergyKilled << G4endl;
     G4cout << "   Max energy of loopers killed: " <<  fMaxEnergyKilled << G4endl;
   } 
 }

Member Function Documentation

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

virtual

Implements G4VProcess.

Definition at line 401 of file G4MonopoleTransportation.cc.

References G4ParticleTable::FindParticle(), fStopAndKill, G4cout, G4endl, G4DynamicParticle::GetDefinition(), G4Track::GetDynamicParticle(), G4Track::GetGlobalTime(), G4Track::GetKineticEnergy(), G4DynamicParticle::GetMass(), G4ParticleTable::GetParticleTable(), G4Step::GetPreStepPoint(), G4Track::GetProperTime(), G4Track::GetStepLength(), G4Track::GetTotalEnergy(), G4StepPoint::GetVelocity(), G4Track::GetVelocity(), G4PropagatorInField::GimmeTrajectoryVectorAndForgetIt(), G4ParticleChangeForTransport::Initialize(), python.hepunit::MeV, G4ParticleChange::ProposeEnergy(), G4ParticleChange::ProposeGlobalTime(), G4ParticleChange::ProposeMomentumDirection(), G4ParticleChange::ProposePolarization(), G4ParticleChange::ProposePosition(), G4ParticleChange::ProposeProperTime(), G4VParticleChange::ProposeTrackStatus(), G4ParticleChangeForTransport::SetMomentumChanged(), G4ParticleChangeForTransport::SetPointerToVectorOfAuxiliaryPoints(), and G4VProcess::verboseLevel.

 {
   static G4int noCalls=0;
   static const G4ParticleDefinition* fOpticalPhoton =
            G4ParticleTable::GetParticleTable()->FindParticle("opticalphoton");
 
   noCalls++;
 
   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 finalVelocity   = track.GetVelocity() ;
      G4double initialVelocity = stepData.GetPreStepPoint()->GetVelocity() ;
      G4double stepLength      = track.GetStepLength() ;
 
      deltaTime= 0.0;  // in case initialVelocity = 0 
      const G4DynamicParticle* fpDynamicParticle = track.GetDynamicParticle();
      if (fpDynamicParticle->GetDefinition()== fOpticalPhoton)
      {
         //  A photon is in the medium of the final point
         //  during the step, so it has the final velocity.
         deltaTime = stepLength/finalVelocity ;
      }
      else if (finalVelocity > 0.0)
      {
         G4double meanInverseVelocity ;
         // deltaTime = stepLength/finalVelocity ;
         meanInverseVelocity = 0.5
                             * ( 1.0 / initialVelocity + 1.0 / finalVelocity ) ;
         deltaTime = stepLength * meanInverseVelocity ;
      }
      else if( initialVelocity > 0.0 )
      {
         deltaTime = stepLength/initialVelocity ;
      }
      fCandidateEndGlobalTime   = startTime + deltaTime ;
   }
   else
   {
      deltaTime = fCandidateEndGlobalTime - startTime ;
   }
 
   fParticleChange.ProposeGlobalTime( fCandidateEndGlobalTime ) ;
 
   // Now Correct by Lorentz factor to get "proper" deltaTime
   
   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 this kills it ...
   //
   if ( fParticleIsLooping )
   {
       G4double endEnergy= fTransportEndKineticEnergy;
 
       if( (endEnergy < fThreshold_Important_Energy) 
           || (fNoLooperTrials >= fThresholdTrials ) ){
         // Kill the looping particle 
         //
         fParticleChange.ProposeTrackStatus( fStopAndKill )  ;
 
         // 'Bare' statistics
         fSumEnergyKilled += endEnergy; 
         if( endEnergy > fMaxEnergyKilled) { fMaxEnergyKilled= endEnergy; }
 
 #ifdef G4VERBOSE
         if( (verboseLevel > 1) || 
             ( endEnergy > fThreshold_Warning_Energy )  ) { 
           G4cout << " G4MonopoleTransportation is killing track that is looping or stuck "
                  << G4endl
                  << "   This track has " << track.GetKineticEnergy() / MeV
                  << " MeV energy." << G4endl;
           G4cout << "   Number of trials = " << fNoLooperTrials 
                  << "   No of calls to AlongStepDoIt = " << noCalls 
                  << G4endl;
         }
 #endif
         fNoLooperTrials=0; 
       }
       else{
         fNoLooperTrials ++; 
 #ifdef G4VERBOSE
         if( (verboseLevel > 2) ){
           G4cout << "   G4MonopoleTransportation::AlongStepDoIt(): Particle looping -  "
                  << "   Number of trials = " << fNoLooperTrials 
                  << "   No of calls to  = " << noCalls 
                  << 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 ;
 }

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

virtual

Implements G4VProcess.

Definition at line 131 of file G4MonopoleTransportation.cc.

References CandidateForSelection, G4Navigator::ComputeSafety(), G4PropagatorInField::ComputeStep(), G4Navigator::ComputeStep(), G4FieldManager::ConfigureForTrack(), G4PropagatorInField::FindAndSetFieldManager(), G4cout, G4endl, G4DynamicParticle::GetCharge(), G4PropagatorInField::GetChordFinder(), G4FieldManager::GetDetectorField(), G4Track::GetDynamicParticle(), G4Track::GetGlobalTime(), G4ChordFinder::GetIntegrationDriver(), G4FieldTrack::GetKineticEnergy(), G4Track::GetKineticEnergy(), G4FieldTrack::GetLabTimeOfFlight(), G4FieldTrack::GetMomentumDir(), G4DynamicParticle::GetMomentumDirection(), G4Track::GetMomentumDirection(), G4ParticleDefinition::GetPDGMass(), G4ParticleDefinition::GetPDGSpin(), G4Track::GetPolarization(), G4FieldTrack::GetPosition(), G4Track::GetPosition(), G4Track::GetProperTime(), G4FieldTrack::GetSpin(), G4MagInt_Driver::GetStepper(), G4DynamicParticle::GetTotalMomentum(), G4Track::GetVelocity(), G4Track::GetVolume(), G4PropagatorInField::IsParticleLooping(), CLHEP::Hep3Vector::mag2(), G4Monopole::MagneticCharge(), G4VParticleChange::ProposeTrueStepLength(), G4EquationOfMotion::SetChargeMomentumMass(), G4SafetyHelper::SetCurrentSafety(), G4MonopoleFieldSetup::SetStepperAndChordFinder(), and sqr().

 {  
   fMagSetup->SetStepperAndChordFinder(1); 
   // change to monopole equation
 
   G4double geometryStepLength, newSafety ; 
   fParticleIsLooping = false ;
 
   // 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 G4DynamicParticle* pParticle      = track.GetDynamicParticle() ;
   G4ThreeVector startMomentumDir          = pParticle->GetMomentumDirection() ;
   G4ThreeVector startPosition             = track.GetPosition() ;
 
   // G4double   theTime        = track.GetGlobalTime() ;
 
   // 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.
   //
   G4ThreeVector OriginShift = startPosition - fPreviousSftOrigin ;
   G4double      MagSqShift  = OriginShift.mag2() ;
   if( MagSqShift >= sqr(fPreviousSafety) )
   {
      currentSafety = 0.0 ;
   }
   else
   {
      currentSafety = fPreviousSafety - std::sqrt(MagSqShift) ;
   }
 
   // Is the monopole charged ?
   //
   G4double particleMagneticCharge = fParticleDef->MagneticCharge() ; 
   G4double particleElectricCharge = pParticle->GetCharge();
 
   fGeometryLimitedStep = false ;
   // fEndGlobalTimeComputed = false ;
 
   // 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 whether the particle have an (EM) field force exerting upon it
   //
   G4FieldManager* fieldMgr=0;
   G4bool          fieldExertsForce = false ;
     
   if( (particleMagneticCharge != 0.0) )
   {      
      fieldMgr= fFieldPropagator->FindAndSetFieldManager( track.GetVolume() ); 
      if (fieldMgr != 0) {
        // Message the field Manager, to configure it for this track
        fieldMgr->ConfigureForTrack( &track );
        // Moved here, in order to allow a transition
        //   from a zero-field  status (with fieldMgr->(field)0
        //   to a finite field  status
 
        // If the field manager has no field, there is no field !
        fieldExertsForce = (fieldMgr->GetDetectorField() != 0);
      }      
   }
 
   // G4cout << " G4Transport:  field exerts force= " << fieldExertsForce
   //          << "  fieldMgr= " << fieldMgr << G4endl;
 
   // Choose the calculation of the transportation: Field or not 
   //
   if( !fieldExertsForce ) 
   {
      G4double linearStepLength ;
      if( fShortStepOptimisation && (currentMinimumStep <= currentSafety) )
      {
        // The Step is guaranteed to be taken
        //
        geometryStepLength   = currentMinimumStep ;
        fGeometryLimitedStep = false ;
      }
      else
      {
        //  Find whether the straight path intersects a volume
        //
        linearStepLength = fLinearNavigator->ComputeStep( startPosition, 
                                                          startMomentumDir,
                                                          currentMinimumStep, 
                                                          newSafety) ;
        // Remember last safety origin & value.
        //
        fPreviousSftOrigin = startPosition ;
        fPreviousSafety    = newSafety ; 
        // fpSafetyHelper->SetCurrentSafety( newSafety, startPosition);
 
        // The safety at the initial point has been re-calculated:
        //
        currentSafety = newSafety ;
           
        fGeometryLimitedStep= (linearStepLength <= currentMinimumStep); 
        if( fGeometryLimitedStep )
        {
          // The geometry limits the Step size (an intersection was found.)
          geometryStepLength   = linearStepLength ;
        } 
        else
        {
          // The full Step is taken.
          geometryStepLength   = currentMinimumStep ;
        }
      }
      endpointDistance = geometryStepLength ;
 
      // Calculate final position
      //
      fTransportEndPosition = startPosition+geometryStepLength*startMomentumDir ;
 
      // Momentum direction, energy and polarisation are unchanged by transport
      //
      fTransportEndMomentumDir   = startMomentumDir ; 
      fTransportEndKineticEnergy = track.GetKineticEnergy() ;
      fTransportEndSpin          = track.GetPolarization();
      fParticleIsLooping         = false ;
      fMomentumChanged           = false ; 
      fEndGlobalTimeComputed     = false ;
   }
   else   //  A field exerts force
   {
      G4double       momentumMagnitude = pParticle->GetTotalMomentum() ;
      G4ThreeVector  EndUnitMomentum ;
      G4double       lengthAlongCurve ;
      G4double       restMass = fParticleDef->GetPDGMass() ;
 
      G4ChargeState chargeState(particleElectricCharge,             // The charge can change (dynamic)
                                fParticleDef->GetPDGSpin(),
                                0,   //  Magnetic moment:  pParticleDef->GetMagneticMoment(),
                                0,   //  Electric Dipole moment - not in Particle Definition 
                                particleMagneticCharge );   // in Mev/c 
 
      G4EquationOfMotion* equationOfMotion = 
      (fFieldPropagator->GetChordFinder()->GetIntegrationDriver()->GetStepper())
      ->GetEquationOfMotion();
 
      equationOfMotion->SetChargeMomentumMass( chargeState,       //  Was particleMagneticCharge - in Mev/c 
                                               momentumMagnitude, //  Was particleElectricCharge 
                                               restMass ) ;  
      // SetChargeMomentumMass now passes both the electric and magnetic charge - in chargeState
 
      G4ThreeVector spin        = track.GetPolarization() ;
      G4FieldTrack  aFieldTrack = G4FieldTrack( startPosition, 
                                                track.GetMomentumDirection(),
                                                0.0, 
                                                track.GetKineticEnergy(),
                                                restMass,
                                                track.GetVelocity(),
                                                track.GetGlobalTime(), // Lab.
                                                track.GetProperTime(), // Part.
                                                &spin                  ) ;
      if( currentMinimumStep > 0 ) 
      {
         // Do the Transport in the field (non recti-linear)
         //
         lengthAlongCurve = fFieldPropagator->ComputeStep( aFieldTrack,
                                                           currentMinimumStep, 
                                                           currentSafety,
                                                           track.GetVolume() ) ;
         fGeometryLimitedStep= lengthAlongCurve < currentMinimumStep; 
         if( fGeometryLimitedStep ) {
            geometryStepLength   = lengthAlongCurve ;
         } else {
            geometryStepLength   = currentMinimumStep ;
         }
      }
      else
      {
         geometryStepLength   = lengthAlongCurve= 0.0 ;
         fGeometryLimitedStep = false ;
      }
 
      // Remember last safety origin & value.
      //
      fPreviousSftOrigin = startPosition ;
      fPreviousSafety    = currentSafety ;         
      // fpSafetyHelper->SetCurrentSafety( newSafety, startPosition);
        
      // Get the End-Position and End-Momentum (Dir-ection)
      //
      fTransportEndPosition = aFieldTrack.GetPosition() ;
 
      // Momentum:  Magnitude and direction can be changed too now ...
      //
      fMomentumChanged         = true ; 
      fTransportEndMomentumDir = aFieldTrack.GetMomentumDir() ;
 
      fTransportEndKineticEnergy  = aFieldTrack.GetKineticEnergy() ; 
 
      fCandidateEndGlobalTime   = aFieldTrack.GetLabTimeOfFlight();
      fEndGlobalTimeComputed    = true;
 
      fTransportEndSpin = aFieldTrack.GetSpin();
      fParticleIsLooping = fFieldPropagator->IsParticleLooping() ;
      endpointDistance   = (fTransportEndPosition - startPosition).mag() ;
   }
 
   // If we are asked to go a step length of 0, and we are on a boundary
   // then a boundary will also limit the step -> we must flag this.
   //
   if( currentMinimumStep == 0.0 ) 
   {
       if( currentSafety == 0.0 )  fGeometryLimitedStep = true ;
   }
 
   // Update the safety starting from the end-point,
   // if it will become negative at the end-point.
   //
   if( currentSafety < endpointDistance ) 
   {
       // if( particleMagneticCharge == 0.0 ) 
       //    G4cout  << "  Avoiding call to ComputeSafety : charge = 0.0 " << G4endl;
  
       if( particleMagneticCharge != 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     += endpointDistance ;
 
 #ifdef G4DEBUG_TRANSPORT 
          G4cout.precision(12) ;
          G4cout << "***G4MonopoleTransportation::AlongStepGPIL ** " << G4endl  ;
          G4cout << "  Called Navigator->ComputeSafety at " << fTransportEndPosition
                 << "    and it returned safety= " << endSafety << G4endl ; 
          G4cout << "  Adding endpoint distance " << endpointDistance 
                 << "   to obtain pseudo-safety= " << currentSafety << G4endl ; 
 #endif
       }
   }            
 
   fParticleChange.ProposeTrueStepLength(geometryStepLength) ;
 
   fMagSetup->SetStepperAndChordFinder(0);
   // change back to usual equation
 
   return geometryStepLength ;
 }

virtual G4VParticleChange* G4MonopoleTransportation::AtRestDoIt	(	const G4Track &	,
		const G4Step &
	)

inlinevirtual

Implements G4VProcess.

Definition at line 135 of file G4MonopoleTransportation.hh.

138 {return 0;};

virtual G4double G4MonopoleTransportation::AtRestGetPhysicalInteractionLength	(	const G4Track &	,
		G4ForceCondition *
	)

inlinevirtual

Implements G4VProcess.

Definition at line 129 of file G4MonopoleTransportation.hh.

132 { return -1.0; };

G4bool G4MonopoleTransportation::DoesGlobalFieldExist ( )

protected

Referenced by StartTracking().

void G4MonopoleTransportation::EnableShortStepOptimisation ( G4bool optimise = true )

inline

G4double G4MonopoleTransportation::GetMaxEnergyKilled ( ) const

inline

G4PropagatorInField* G4MonopoleTransportation::GetPropagatorInField ( )

G4double G4MonopoleTransportation::GetSumEnergyKilled ( ) const

inline

G4double G4MonopoleTransportation::GetThresholdImportantEnergy ( ) const

inline

G4int G4MonopoleTransportation::GetThresholdTrials ( ) const

inline

G4double G4MonopoleTransportation::GetThresholdWarningEnergy ( ) const

inline

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

virtual

Implements G4VProcess.

Definition at line 549 of file G4MonopoleTransportation.cc.

References fStopAndKill, G4VPhysicalVolume::GetLogicalVolume(), G4MaterialCutsCouple::GetMaterial(), G4LogicalVolume::GetMaterial(), G4ProductionCutsTable::GetMaterialCutsCouple(), G4LogicalVolume::GetMaterialCutsCouple(), G4Track::GetMomentumDirection(), G4Track::GetPosition(), G4MaterialCutsCouple::GetProductionCuts(), G4ProductionCutsTable::GetProductionCutsTable(), G4LogicalVolume::GetSensitiveDetector(), G4Track::GetTouchableHandle(), G4Track::GetTrackStatus(), G4VTouchable::GetVolume(), G4Navigator::LocateGlobalPointWithinVolume(), G4VParticleChange::ProposeTrackStatus(), G4Navigator::SetGeometricallyLimitedStep(), G4ParticleChangeForTransport::SetMaterialCutsCoupleInTouchable(), G4ParticleChangeForTransport::SetMaterialInTouchable(), G4ParticleChangeForTransport::SetSensitiveDetectorInTouchable(), and G4ParticleChangeForTransport::SetTouchableHandle().

 {
   G4TouchableHandle retCurrentTouchable ;   // The one to return
 
   // 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 ) ;
   }
   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() ;
   }         // endif ( fGeometryLimitedStep ) 
 
   const G4VPhysicalVolume* pNewVol = retCurrentTouchable->GetVolume() ;
   const G4Material* pNewMaterial   = 0 ;
   const G4VSensitiveDetector* pNewSensitiveDetector   = 0 ;
                                                                                        
   if( pNewVol != 0 )
   {
     pNewMaterial= pNewVol->GetLogicalVolume()->GetMaterial();
     pNewSensitiveDetector= pNewVol->GetLogicalVolume()->GetSensitiveDetector();
   }
 
   // ( <const_cast> pNewMaterial ) ;
   // ( <const_cast> pNewSensitiveDetector) ;
 
   fParticleChange.SetMaterialInTouchable( (G4Material *) pNewMaterial ) ;
   fParticleChange.SetSensitiveDetectorInTouchable( (G4VSensitiveDetector *) pNewSensitiveDetector ) ;
 
   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 );
 
   // temporarily until Get/Set Material of ParticleChange, 
   // and StepPoint can be made const. 
   // 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(retCurrentTouchable) ;
   
   return &fParticleChange ;
 }

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

virtual

Implements G4VProcess.

Definition at line 539 of file G4MonopoleTransportation.cc.

References DBL_MAX, and Forced.

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

void G4MonopoleTransportation::ResetKilledStatistics ( G4int report = 1 )

inline

void G4MonopoleTransportation::SetPropagatorInField ( G4PropagatorInField * pFieldPropagator )

void G4MonopoleTransportation::SetThresholdImportantEnergy ( G4double newEnImp )

inline

void G4MonopoleTransportation::SetThresholdTrials ( G4int newMaxTrials )

inline

void G4MonopoleTransportation::SetThresholdWarningEnergy ( G4double newEnWarn )

inline

void G4MonopoleTransportation::StartTracking ( G4Track * aTrack )

virtual

Reimplemented from G4VProcess.

Definition at line 651 of file G4MonopoleTransportation.cc.

References G4FieldManagerStore::ClearAllChordFindersState(), G4PropagatorInField::ClearPropagatorState(), DoesGlobalFieldExist(), G4FieldManagerStore::GetInstance(), G4Track::GetTouchableHandle(), and G4VProcess::StartTracking().

 {
   G4VProcess::StartTracking(aTrack);
 
 // The actions here are those that were taken in AlongStepGPIL
 //   when track.GetCurrentStepNumber()==1
 
   // 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( DoesGlobalFieldExist() ) {
      fFieldPropagator->ClearPropagatorState();   
        // Resets all state of field propagator class (ONLY)
        //  including safety values (in case of overlaps and to wipe for first track).
 
      // G4ChordFinder* chordF= fFieldPropagator->GetChordFinder();
      // if( chordF ) chordF->ResetStepEstimate();
   }
 
   // Make sure to clear the chord finders of all fields (ie managers)
   static G4FieldManagerStore* fieldMgrStore= G4FieldManagerStore::GetInstance();
   fieldMgrStore->ClearAllChordFindersState(); 
 
   // Update the current touchable handle  (from the track's)
   //
   fCurrentTouchableHandle = aTrack->GetTouchableHandle();
 }

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

Public Member Functions

Protected Member Functions

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation