G4PropagatorInField Class Reference

#include <G4PropagatorInField.hh>

Collaboration diagram for G4PropagatorInField:

Public Member Functions
	G4PropagatorInField (G4Navigator *theNavigator, G4FieldManager *detectorFieldMgr, G4VIntersectionLocator *vLocator=0)
	~G4PropagatorInField ()
G4double	ComputeStep (G4FieldTrack &pFieldTrack, G4double pCurrentProposedStepLength, G4double &pNewSafety, G4VPhysicalVolume *pPhysVol=0)
G4ThreeVector	EndPosition () const
G4ThreeVector	EndMomentumDir () const
G4bool	IsParticleLooping () const
G4double	GetEpsilonStep () const
void	SetEpsilonStep (G4double newEps)
G4FieldManager *	FindAndSetFieldManager (G4VPhysicalVolume *pCurrentPhysVol)
G4ChordFinder *	GetChordFinder ()
G4int	SetVerboseLevel (G4int verbose)
G4int	GetVerboseLevel () const
G4int	Verbose () const
void	SetVerboseTrace (G4bool enable)
G4bool	GetVerboseTrace ()
G4int	GetMaxLoopCount () const
void	SetMaxLoopCount (G4int new_max)
void	printStatus (const G4FieldTrack &startFT, const G4FieldTrack &currentFT, G4double requestStep, G4double safety, G4int step, G4VPhysicalVolume *startVolume)
G4FieldTrack	GetEndState () const
G4double	GetMinimumEpsilonStep () const
void	SetMinimumEpsilonStep (G4double newEpsMin)
G4double	GetMaximumEpsilonStep () const
void	SetMaximumEpsilonStep (G4double newEpsMax)
void	SetLargestAcceptableStep (G4double newBigDist)
G4double	GetLargestAcceptableStep ()
void	SetTrajectoryFilter (G4VCurvedTrajectoryFilter *filter)
std::vector< G4ThreeVector > *	GimmeTrajectoryVectorAndForgetIt () const
void	ClearPropagatorState ()
void	SetDetectorFieldManager (G4FieldManager *newGlobalFieldManager)
void	SetUseSafetyForOptimization (G4bool)
G4bool	GetUseSafetyForOptimization ()
G4bool	IntersectChord (const G4ThreeVector &StartPointA, const G4ThreeVector &EndPointB, G4double &NewSafety, G4double &LinearStepLength, G4ThreeVector &IntersectionPoint)
G4bool	IsFirstStepInVolume ()
G4bool	IsLastStepInVolume ()
void	PrepareNewTrack ()
G4VIntersectionLocator *	GetIntersectionLocator ()
void	SetIntersectionLocator (G4VIntersectionLocator *pLocator)
G4double	GetDeltaIntersection () const
G4double	GetDeltaOneStep () const
G4FieldManager *	GetCurrentFieldManager ()
G4EquationOfMotion *	GetCurrentEquationOfMotion ()
void	SetNavigatorForPropagating (G4Navigator *SimpleOrMultiNavigator)
G4Navigator *	GetNavigatorForPropagating ()
void	SetThresholdNoZeroStep (G4int noAct, G4int noHarsh, G4int noAbandon)
G4int	GetThresholdNoZeroSteps (G4int i)
G4double	GetZeroStepThreshold ()
void	SetZeroStepThreshold (G4double newLength)
void	RefreshIntersectionLocator ()

Protected Member Functions
void	PrintStepLengthDiagnostic (G4double currentProposedStepLength, G4double decreaseFactor, G4double stepTrial, const G4FieldTrack &aFieldTrack)
void	ReportLoopingParticle (G4int count, double StepTaken, G4VPhysicalVolume *pPhysVol)
void	ReportStuckParticle (G4int noZeroSteps, G4double proposedStep, G4double lastTriedStep, G4VPhysicalVolume *physVol)

Private Attributes
G4int	fMax_loop_count
G4bool	fUseSafetyForOptimisation
G4int	fActionThreshold_NoZeroSteps
G4int	fSevereActionThreshold_NoZeroSteps
G4int	fAbandonThreshold_NoZeroSteps
G4double	fZeroStepThreshold
G4double	fLargestAcceptableStep
G4double	kCarTolerance
G4bool	fAllocatedLocator
G4FieldManager *	fDetectorFieldMgr
G4VIntersectionLocator *	fIntersectionLocator
G4VCurvedTrajectoryFilter *	fpTrajectoryFilter
G4Navigator *	fNavigator
G4FieldManager *	fCurrentFieldMgr
G4bool	fSetFieldMgr
G4double	fEpsilonStep
G4FieldTrack	End_PointAndTangent
G4bool	fParticleIsLooping
G4int	fNoZeroStep
G4double	fFull_CurveLen_of_LastAttempt
G4double	fLast_ProposedStepLength
G4ThreeVector	fPreviousSftOrigin
G4double	fPreviousSafety
G4int	fVerboseLevel
G4bool	fVerbTracePiF
G4bool	fFirstStepInVolume
G4bool	fLastStepInVolume
G4bool	fNewTrack

Detailed Description

Definition at line 64 of file G4PropagatorInField.hh.

Constructor & Destructor Documentation

G4PropagatorInField()

G4PropagatorInField::G4PropagatorInField	(	G4Navigator *	theNavigator,
		G4FieldManager *	detectorFieldMgr,
		G4VIntersectionLocator *	vLocator = 0
	)

Definition at line 60 of file G4PropagatorInField.cc.

  : 
    fMax_loop_count(1000),
    fUseSafetyForOptimisation(true),   // (false) is less sensitive to incorrect safety
    fZeroStepThreshold( 0.0 ),         // length of what is recognised as 'zero' step
    fDetectorFieldMgr(detectorFieldMgr), 
    fpTrajectoryFilter( 0 ),
    fNavigator(theNavigator),
    fCurrentFieldMgr(detectorFieldMgr),
    fSetFieldMgr(false),
    End_PointAndTangent(G4ThreeVector(0.,0.,0.),
                        G4ThreeVector(0.,0.,0.),0.0,0.0,0.0,0.0,0.0),
    fParticleIsLooping(false),
    fNoZeroStep(0), 
    fVerboseLevel(0),
    fVerbTracePiF(false),
    fFirstStepInVolume(true),
    fLastStepInVolume(true),
    fNewTrack(true)
{
  if(fDetectorFieldMgr) { fEpsilonStep = fDetectorFieldMgr->GetMaximumEpsilonStep();}
  else                  { fEpsilonStep= 1.0e-5; } 
  fActionThreshold_NoZeroSteps = 2; 
  fSevereActionThreshold_NoZeroSteps = 10; 
  fAbandonThreshold_NoZeroSteps = 50; 
  fFull_CurveLen_of_LastAttempt = -1; 
  fLast_ProposedStepLength = -1;
  fLargestAcceptableStep = 1000.0 * meter;

  fPreviousSftOrigin= G4ThreeVector(0.,0.,0.);
  fPreviousSafety= 0.0;
  kCarTolerance = G4GeometryTolerance::GetInstance()->GetSurfaceTolerance();
  fZeroStepThreshold= std::max( 1.0e5 * kCarTolerance, 1.0e-1 * micrometer );

#ifdef G4DEBUG_FIELD
  G4cout << " PiF: Zero Step Threshold set to "
         << fZeroStepThreshold / millimeter
     << " mm." << G4endl;
  G4cout << " PiF:   Value of kCarTolerance = "
         << kCarTolerance / millimeter 
     << " mm. " << G4endl;
  fVerboseLevel = 3;   
  fVerbTracePiF = true;   
#endif 

  // Defining Intersection Locator and his parameters
  if (vLocator==0)
  {
    fIntersectionLocator= new G4MultiLevelLocator(theNavigator);
    fAllocatedLocator= true;
  }
  else
  {
    fIntersectionLocator= vLocator;
    fAllocatedLocator= false;
  }
  RefreshIntersectionLocator();  //  Copy all relevant parameters
}

Here is the call graph for this function:

~G4PropagatorInField()

G4PropagatorInField::~G4PropagatorInField

(

)

Definition at line 123 of file G4PropagatorInField.cc.

{
  if(fAllocatedLocator)  { delete  fIntersectionLocator; }
}

Member Function Documentation

ClearPropagatorState()

void G4PropagatorInField::ClearPropagatorState

(

)

Definition at line 669 of file G4PropagatorInField.cc.

{
  // Goal: Clear all memory of previous steps,  cached information

  fParticleIsLooping= false;
  fNoZeroStep= 0;

  End_PointAndTangent= G4FieldTrack( G4ThreeVector(0.,0.,0.),
                                     G4ThreeVector(0.,0.,0.),
                                     0.0,0.0,0.0,0.0,0.0); 
  fFull_CurveLen_of_LastAttempt = -1; 
  fLast_ProposedStepLength = -1;

  fPreviousSftOrigin= G4ThreeVector(0.,0.,0.);
  fPreviousSafety= 0.0;
}

Here is the call graph for this function:

Here is the caller graph for this function:

ComputeStep()

G4double G4PropagatorInField::ComputeStep	(	G4FieldTrack &	pFieldTrack,
		G4double	pCurrentProposedStepLength,
		G4double &	pNewSafety,
		G4VPhysicalVolume *	pPhysVol = 0
	)

Definition at line 145 of file G4PropagatorInField.cc.

{  
  // If CurrentProposedStepLength is too small for finding Chords
  // then return with no action (for now - TODO: some action)
  //
  if(CurrentProposedStepLength<kCarTolerance)
  {
    return kInfinity;
  }

  // Introducing smooth trajectory display (jacek 01/11/2002)
  //
  if (fpTrajectoryFilter)
  {
    fpTrajectoryFilter->CreateNewTrajectorySegment();
  }

  fFirstStepInVolume = fNewTrack ? true : fLastStepInVolume;
  fLastStepInVolume= false;
  fNewTrack= false; 

  if( fVerboseLevel > 2 )
  {
    G4cout << "G4PropagatorInField::ComputeStep() called" << G4endl;
    G4cout << "   Starting FT: " << pFieldTrack;
    G4cout << "   Requested length = " << CurrentProposedStepLength << G4endl;
    G4cout << "   PhysVol = ";
    if( pPhysVol )
       G4cout << pPhysVol->GetName() << G4endl;
    else
       G4cout << " N/A ";
    G4cout << G4endl;
  }
  
  // Parameters for adaptive Runge-Kutta integration
  
  G4double      h_TrialStepSize;        // 1st Step Size 
  G4double      TruePathLength = CurrentProposedStepLength;
  G4double      StepTaken = 0.0; 
  G4double      s_length_taken, epsilon ; 
  G4bool        intersects;
  G4bool        first_substep = true;

  G4double      NewSafety;
  fParticleIsLooping = false;

  // If not yet done, 
  //   Set the field manager to the local  one if the volume has one, 
  //                      or to the global one if not
  //
  if( !fSetFieldMgr ) fCurrentFieldMgr= FindAndSetFieldManager( pPhysVol ); 
  // For the next call, the field manager must again be set
  fSetFieldMgr= false;

  G4FieldTrack  CurrentState(pFieldTrack);
  G4FieldTrack  OriginalState = CurrentState;

  // If the Step length is "infinite", then an approximate-maximum Step
  // length (used to calculate the relative accuracy) must be guessed.
  //
  if( CurrentProposedStepLength >= fLargestAcceptableStep )
  {
    G4ThreeVector StartPointA, VelocityUnit;
    StartPointA  = pFieldTrack.GetPosition();
    VelocityUnit = pFieldTrack.GetMomentumDir();

    G4double trialProposedStep = 1.e2 * ( 10.0 * cm + 
      fNavigator->GetWorldVolume()->GetLogicalVolume()->
                  GetSolid()->DistanceToOut(StartPointA, VelocityUnit) );
    CurrentProposedStepLength= std::min( trialProposedStep,
                                           fLargestAcceptableStep ); 
  }
  epsilon = fCurrentFieldMgr->GetDeltaOneStep() / CurrentProposedStepLength;
  // G4double raw_epsilon= epsilon;
  G4double epsilonMin= fCurrentFieldMgr->GetMinimumEpsilonStep();
  G4double epsilonMax= fCurrentFieldMgr->GetMaximumEpsilonStep();
  if( epsilon < epsilonMin ) epsilon = epsilonMin;
  if( epsilon > epsilonMax ) epsilon = epsilonMax;
  SetEpsilonStep( epsilon );


  // Values for Intersection Locator has to be updated on each call for the
  // case that CurrentFieldManager has changed from the one of previous step
  RefreshIntersectionLocator();

  // G4cout << "G4PiF: Epsilon of current step - raw= " << raw_epsilon
  //        << " final= " << epsilon << G4endl;

  //  Shorten the proposed step in case of earlier problems (zero steps)
  // 
  if( fNoZeroStep > fActionThreshold_NoZeroSteps )
  {
    G4double stepTrial;

    stepTrial= fFull_CurveLen_of_LastAttempt; 
    if( (stepTrial <= 0.0) && (fLast_ProposedStepLength > 0.0) ) 
      stepTrial= fLast_ProposedStepLength; 

    G4double decreaseFactor = 0.9; // Unused default
    if(   (fNoZeroStep < fSevereActionThreshold_NoZeroSteps)
       && (stepTrial > 100.0*fZeroStepThreshold) )
    {
      // Attempt quick convergence
      //
      decreaseFactor= 0.25;
    } 
    else
    {
      // We are in significant difficulties, probably at a boundary that
      // is either geometrically sharp or between very different materials.
      // Careful decreases to cope with tolerance are required.
      //
      if( stepTrial > 100.0*fZeroStepThreshold )
        decreaseFactor = 0.35;     // Try decreasing slower
      else if( stepTrial > 30.0*fZeroStepThreshold )
        decreaseFactor= 0.5;       // Try yet slower decrease
      else if( stepTrial > 10.0*fZeroStepThreshold )
        decreaseFactor= 0.75;      // Try even slower decreases
      else
        decreaseFactor= 0.9;       // Try very slow decreases
     }
     stepTrial *= decreaseFactor;

#ifdef G4DEBUG_FIELD
     G4cerr << " G4PropagatorInField::ComputeStep(): " << G4endl
        << "  Decreasing step -  in volume " << pPhysVol;
     if( pPhysVol )
        G4cerr << "    with name " << pPhysVol->GetName(); 
     G4cerr << G4endl;
     PrintStepLengthDiagnostic(CurrentProposedStepLength, decreaseFactor,
                               stepTrial, pFieldTrack);
#endif
     if( stepTrial == 0.0 )  //  Change to make it < 0.1 * kCarTolerance ??
     {
       std::ostringstream message;
       message << "Particle abandoned due to lack of progress in field."
               << G4endl
               << "  Properties : " << pFieldTrack << G4endl
               << "  Attempting a zero step = " << stepTrial << G4endl
               << "  while attempting to progress after " << fNoZeroStep
               << " trial steps. Will abandon step.";
       G4Exception("G4PropagatorInField::ComputeStep()", "GeomNav1002",
                   JustWarning, message);
       fParticleIsLooping= true;
       return 0;  // = stepTrial;
     }
     if( stepTrial < CurrentProposedStepLength )
       CurrentProposedStepLength = stepTrial;
  }
  fLast_ProposedStepLength = CurrentProposedStepLength;

  G4int do_loop_count = 0; 
  do  // Loop checking, 07.10.2016, J.Apostolakis
  { 
    G4FieldTrack SubStepStartState = CurrentState;
    G4ThreeVector SubStartPoint = CurrentState.GetPosition(); 
    
    if(!first_substep)
    {
      if( fVerboseLevel > 4 )
      {
        G4cout << " PiF: Calling Nav/Locate Global Point within-Volume "
               << G4endl;
      }
      fNavigator->LocateGlobalPointWithinVolume( SubStartPoint );
    }

    // How far to attempt to move the particle !
    //
    h_TrialStepSize = CurrentProposedStepLength - StepTaken;

    // Integrate as far as "chord miss" rule allows.
    //
    s_length_taken = GetChordFinder()->AdvanceChordLimited( 
                             CurrentState,    // Position & velocity
                             h_TrialStepSize,
                             fEpsilonStep,
                             fPreviousSftOrigin,
                             fPreviousSafety
                             );
    //  CurrentState is now updated with the final position and velocity. 

    fFull_CurveLen_of_LastAttempt = s_length_taken;

    G4ThreeVector  EndPointB = CurrentState.GetPosition(); 
    G4ThreeVector  InterSectionPointE;
    G4double       LinearStepLength;
 
    // Intersect chord AB with geometry
    intersects= IntersectChord( SubStartPoint, EndPointB, 
                                NewSafety,     LinearStepLength, 
                                InterSectionPointE );
    // E <- Intersection Point of chord AB and either volume A's surface 
    //                                  or a daughter volume's surface ..

    if( first_substep )
    { 
       currentSafety = NewSafety;
    } // Updating safety in other steps is potential future extention

    if( intersects )
    {
       G4FieldTrack IntersectPointVelct_G(CurrentState);  // FT-Def-Construct

       // Find the intersection point of AB true path with the surface
       //   of vol(A), if it exists. Start with point E as first "estimate".
       G4bool recalculatedEndPt= false;
       
       G4bool found_intersection = fIntersectionLocator->
         EstimateIntersectionPoint( SubStepStartState, CurrentState, 
                                    InterSectionPointE, IntersectPointVelct_G,
                                    recalculatedEndPt, fPreviousSafety,
                                    fPreviousSftOrigin);
       intersects = found_intersection;
       if( found_intersection )
       {        
          End_PointAndTangent= IntersectPointVelct_G;  // G is our EndPoint ...
          StepTaken = TruePathLength = IntersectPointVelct_G.GetCurveLength()
                                      - OriginalState.GetCurveLength();
       }
       else
       {
          // Either "minor" chords do not intersect
          // or else stopped (due to too many steps)
          //
          if( recalculatedEndPt )
          {
             G4double endAchieved = IntersectPointVelct_G.GetCurveLength();
             G4double endExpected = CurrentState.GetCurveLength(); 

             // Detect failure - due to too many steps
             G4bool shortEnd = endAchieved
                             < (endExpected*(1.0-CLHEP::perMillion));

             G4double stepAchieved = endAchieved
                                   - SubStepStartState.GetCurveLength();

             // Update remaining state - must work for 'full' step or
             // abandonned intersection
             //
             CurrentState= IntersectPointVelct_G;
             s_length_taken = stepAchieved;
             if( shortEnd )
             {
                fParticleIsLooping = true;
             } 
          }
       }
    }
    if( !intersects )
    {
      StepTaken += s_length_taken; 
      // For smooth trajectory display (jacek 01/11/2002)
      if (fpTrajectoryFilter) {
        fpTrajectoryFilter->TakeIntermediatePoint(CurrentState.GetPosition());
      }
    }
    first_substep = false;

#ifdef G4DEBUG_FIELD
    if( fNoZeroStep > fActionThreshold_NoZeroSteps )
    {
      printStatus( SubStepStartState,  // or OriginalState,
                   CurrentState,  CurrentProposedStepLength, 
                   NewSafety,     do_loop_count,  pPhysVol );
    }
    if( (fVerboseLevel > 1) && (do_loop_count > fMax_loop_count-10 ))
    {
      if( do_loop_count == fMax_loop_count-9 )
      {
        G4cout << " G4PropagatorInField::ComputeStep(): " << G4endl
               << "  Difficult track - taking many sub steps." << G4endl;
      }
      printStatus( SubStepStartState, CurrentState, CurrentProposedStepLength, 
                   NewSafety, do_loop_count, pPhysVol );
    }
#endif

    do_loop_count++;

  } while( (!intersects )
        && (!fParticleIsLooping)
        && (StepTaken + kCarTolerance < CurrentProposedStepLength)  
        && ( do_loop_count < fMax_loop_count ) );

  if( do_loop_count >= fMax_loop_count  )
  {
    fParticleIsLooping = true;
  }
  if ( fParticleIsLooping && (fVerboseLevel > 0) )
  {
    ReportLoopingParticle( do_loop_count, StepTaken, pPhysVol );
  }
    
  if( !intersects )
  {
    // Chord AB or "minor chords" do not intersect
    // B is the endpoint Step of the current Step.
    //
    End_PointAndTangent = CurrentState; 
    TruePathLength = StepTaken;   //  Original code
    // Tried the following to avoid potential issue with round-off error
    // - but has issues... Suppressing this change JA 2015/05/02
    // TruePathLength = CurrentProposedStepLength;
  }
  fLastStepInVolume = intersects;
  
  // Set pFieldTrack to the return value
  //
  pFieldTrack = End_PointAndTangent;

#ifdef G4VERBOSE
  // Check that "s" is correct
  //
  if( std::fabs(OriginalState.GetCurveLength() + TruePathLength 
      - End_PointAndTangent.GetCurveLength()) > 3.e-4 * TruePathLength )
  {
    std::ostringstream message;
    message << "Curve length mis-match between original state "
            << "and proposed endpoint of propagation." << G4endl
            << "  The curve length of the endpoint should be: " 
            << OriginalState.GetCurveLength() + TruePathLength << G4endl
            << "  and it is instead: "
            << End_PointAndTangent.GetCurveLength() << "." << G4endl
            << "  A difference of: "
            << OriginalState.GetCurveLength() + TruePathLength 
               - End_PointAndTangent.GetCurveLength() << G4endl
            << "  Original state = " << OriginalState   << G4endl
            << "  Proposed state = " << End_PointAndTangent;
    G4Exception("G4PropagatorInField::ComputeStep()",
                "GeomNav0003", FatalException, message);
  }
#endif

  if( TruePathLength+kCarTolerance >= CurrentProposedStepLength )
  {
     fNoZeroStep = 0;     
  }
  else
  {     
     // In particular anomalous cases, we can get repeated zero steps
     // We identify these cases and take corrective action when they occur.
     // 
     if( TruePathLength < std::max( fZeroStepThreshold, 0.5*kCarTolerance ) )
     {
        fNoZeroStep++;
     }
     else{
        fNoZeroStep = 0;
     }
  }
  if( fNoZeroStep > fAbandonThreshold_NoZeroSteps )
  { 
     fParticleIsLooping = true;
     ReportStuckParticle( fNoZeroStep, CurrentProposedStepLength, fFull_CurveLen_of_LastAttempt,
                               pPhysVol );
     fNoZeroStep = 0; 
  }
 
  return TruePathLength;
}

Here is the call graph for this function:

Here is the caller graph for this function:

EndMomentumDir()

G4ThreeVector G4PropagatorInField::EndMomentumDir ( ) const

inline

EndPosition()

G4ThreeVector G4PropagatorInField::EndPosition ( ) const

inline

FindAndSetFieldManager()

G4FieldManager * G4PropagatorInField::FindAndSetFieldManager

(

G4VPhysicalVolume *

pCurrentPhysVol

)

Definition at line 687 of file G4PropagatorInField.cc.

{
  G4FieldManager* currentFieldMgr;

  currentFieldMgr = fDetectorFieldMgr;
  if( pCurrentPhysicalVolume)
  {
     G4FieldManager *pRegionFieldMgr= 0, *localFieldMgr = 0;
     G4LogicalVolume* pLogicalVol= pCurrentPhysicalVolume->GetLogicalVolume();

     if( pLogicalVol ) { 
    // Value for Region, if any, Overrides 
    G4Region*  pRegion= pLogicalVol->GetRegion();
    if( pRegion ) { 
       pRegionFieldMgr= pRegion->GetFieldManager();
       if( pRegionFieldMgr ) 
         currentFieldMgr= pRegionFieldMgr;
    }

    // 'Local' Value from logical volume, if any, Overrides 
    localFieldMgr= pLogicalVol->GetFieldManager();
    if ( localFieldMgr ) 
       currentFieldMgr = localFieldMgr;
     }
  }
  fCurrentFieldMgr= currentFieldMgr;

  // Flag that field manager has been set
  //
  fSetFieldMgr= true;

  return currentFieldMgr;
}

Here is the call graph for this function:

Here is the caller graph for this function:

GetChordFinder()

G4ChordFinder* G4PropagatorInField::GetChordFinder ( )

inline

Here is the caller graph for this function:

GetCurrentEquationOfMotion()

G4EquationOfMotion* G4PropagatorInField::GetCurrentEquationOfMotion ( )

inline

Here is the caller graph for this function:

GetCurrentFieldManager()

G4FieldManager* G4PropagatorInField::GetCurrentFieldManager ( )

inline

Here is the caller graph for this function:

GetDeltaIntersection()

G4double G4PropagatorInField::GetDeltaIntersection ( ) const

inline

Here is the caller graph for this function:

GetDeltaOneStep()

G4double G4PropagatorInField::GetDeltaOneStep ( ) const

inline

Here is the caller graph for this function:

GetEndState()

G4FieldTrack G4PropagatorInField::GetEndState ( ) const

inline

Here is the caller graph for this function:

GetEpsilonStep()

G4double G4PropagatorInField::GetEpsilonStep ( ) const

inline

GetIntersectionLocator()

G4VIntersectionLocator* G4PropagatorInField::GetIntersectionLocator ( )

inline

Here is the caller graph for this function:

GetLargestAcceptableStep()

G4double G4PropagatorInField::GetLargestAcceptableStep ( )

inline

Here is the caller graph for this function:

GetMaximumEpsilonStep()

G4double G4PropagatorInField::GetMaximumEpsilonStep ( ) const

inline

Here is the caller graph for this function:

GetMaxLoopCount()

G4int G4PropagatorInField::GetMaxLoopCount ( ) const

inline

Here is the caller graph for this function:

GetMinimumEpsilonStep()

G4double G4PropagatorInField::GetMinimumEpsilonStep ( ) const

inline

Here is the caller graph for this function:

GetNavigatorForPropagating()

G4Navigator* G4PropagatorInField::GetNavigatorForPropagating ( )

inline

Here is the caller graph for this function:

GetThresholdNoZeroSteps()

G4int G4PropagatorInField::GetThresholdNoZeroSteps ( G4int  i )

inline

Here is the caller graph for this function:

GetUseSafetyForOptimization()

G4bool G4PropagatorInField::GetUseSafetyForOptimization ( )

inline

Here is the caller graph for this function:

GetVerboseLevel()

G4int G4PropagatorInField::GetVerboseLevel ( ) const

inline

GetVerboseTrace()

G4bool G4PropagatorInField::GetVerboseTrace ( )

inline

Definition at line 103 of file G4PropagatorInField.hh.

{ return fVerbTracePiF; }

Here is the call graph for this function:

GetZeroStepThreshold()

G4double G4PropagatorInField::GetZeroStepThreshold ( )

inline

Here is the caller graph for this function:

GimmeTrajectoryVectorAndForgetIt()

std::vector< G4ThreeVector > * G4PropagatorInField::GimmeTrajectoryVectorAndForgetIt

(

)

const

Definition at line 646 of file G4PropagatorInField.cc.

{
  // NB, GimmeThePointsAndForgetThem really forgets them, so it can
  // only be called (exactly) once for each step.

  if (fpTrajectoryFilter)
  {
    return fpTrajectoryFilter->GimmeThePointsAndForgetThem();
  }
  else
  {
    return 0;
  }
}

Here is the call graph for this function:

Here is the caller graph for this function:

IntersectChord()

G4bool G4PropagatorInField::IntersectChord ( const G4ThreeVector &  StartPointA, const G4ThreeVector &  EndPointB, G4double &  NewSafety, G4double &  LinearStepLength, G4ThreeVector &  IntersectionPoint  )

inline

Here is the caller graph for this function:

IsFirstStepInVolume()

G4bool G4PropagatorInField::IsFirstStepInVolume ( )

inline

Definition at line 162 of file G4PropagatorInField.hh.

{ return fFirstStepInVolume; }

IsLastStepInVolume()

G4bool G4PropagatorInField::IsLastStepInVolume ( )

inline

Definition at line 163 of file G4PropagatorInField.hh.

{ return fLastStepInVolume;  }

Here is the caller graph for this function:

IsParticleLooping()

G4bool G4PropagatorInField::IsParticleLooping ( ) const

inline

Here is the caller graph for this function:

PrepareNewTrack()

void G4PropagatorInField::PrepareNewTrack ( )

inline

Definition at line 164 of file G4PropagatorInField.hh.

{ fNewTrack = true; fFirstStepInVolume=false; fLastStepInVolume=false; } 

Here is the call graph for this function:

Here is the caller graph for this function:

printStatus()

void G4PropagatorInField::printStatus	(	const G4FieldTrack &	startFT,
		const G4FieldTrack &	currentFT,
		G4double	requestStep,
		G4double	safety,
		G4int	step,
		G4VPhysicalVolume *	startVolume
	)

Definition at line 516 of file G4PropagatorInField.cc.

{
  const G4int verboseLevel=fVerboseLevel;
  const G4ThreeVector StartPosition       = StartFT.GetPosition();
  const G4ThreeVector StartUnitVelocity   = StartFT.GetMomentumDir();
  const G4ThreeVector CurrentPosition     = CurrentFT.GetPosition();
  const G4ThreeVector CurrentUnitVelocity = CurrentFT.GetMomentumDir();

  G4double step_len = CurrentFT.GetCurveLength() - StartFT.GetCurveLength();

  G4int oldprec;   // cout/cerr precision settings
      
  if( ((stepNo == 0) && (verboseLevel <3)) || (verboseLevel >= 3) )
  {
    oldprec = G4cout.precision(4);
    G4cout << std::setw( 6)  << " " 
           << std::setw( 25) << " Current Position  and  Direction" << " "
           << G4endl; 
    G4cout << std::setw( 5) << "Step#" 
           << std::setw(10) << "  s  " << " "
           << std::setw(10) << "X(mm)" << " "
           << std::setw(10) << "Y(mm)" << " "  
           << std::setw(10) << "Z(mm)" << " "
           << std::setw( 7) << " N_x " << " "
           << std::setw( 7) << " N_y " << " "
           << std::setw( 7) << " N_z " << " " ;
    G4cout << std::setw( 7) << " Delta|N|" << " "
           << std::setw( 9) << "StepLen" << " "  
           << std::setw(12) << "StartSafety" << " "  
           << std::setw( 9) << "PhsStep" << " ";  
    if( startVolume )
      { G4cout << std::setw(18) << "NextVolume" << " "; }
    G4cout.precision(oldprec);
    G4cout << G4endl;
  }
  if((stepNo == 0) && (verboseLevel <=3))
  {
    // Recurse to print the start values
    //
    printStatus( StartFT, StartFT, -1.0, safety, -1, startVolume);
  }
  if( verboseLevel <= 3 )
  {
    if( stepNo >= 0)
      { G4cout << std::setw( 4) << stepNo << " "; }
    else
      { G4cout << std::setw( 5) << "Start" ; }
    oldprec = G4cout.precision(8);
    G4cout << std::setw(10) << CurrentFT.GetCurveLength() << " "; 
    G4cout.precision(8);
    G4cout << std::setw(10) << CurrentPosition.x() << " "
           << std::setw(10) << CurrentPosition.y() << " "
           << std::setw(10) << CurrentPosition.z() << " ";
    G4cout.precision(4);
    G4cout << std::setw( 7) << CurrentUnitVelocity.x() << " "
           << std::setw( 7) << CurrentUnitVelocity.y() << " "
           << std::setw( 7) << CurrentUnitVelocity.z() << " ";
    G4cout.precision(3); 
    G4cout << std::setw( 7)
           << CurrentFT.GetMomentum().mag()-StartFT.GetMomentum().mag() << " "; 
    G4cout << std::setw( 9) << step_len << " "; 
    G4cout << std::setw(12) << safety << " ";
    if( requestStep != -1.0 ) 
      { G4cout << std::setw( 9) << requestStep << " "; }
    else
      { G4cout << std::setw( 9) << "Init/NotKnown" << " "; }
    if( startVolume != 0)
      { G4cout << std::setw(12) << startVolume->GetName() << " "; }
    G4cout.precision(oldprec);
    G4cout << G4endl;
  }
  else // if( verboseLevel > 3 )
  {
    //  Multi-line output
      
    G4cout << "Step taken was " << step_len  
           << " out of PhysicalStep = " <<  requestStep << G4endl;
    G4cout << "Final safety is: " << safety << G4endl;
    G4cout << "Chord length = " << (CurrentPosition-StartPosition).mag()
           << G4endl;
    G4cout << G4endl; 
  }
}

Here is the call graph for this function:

Here is the caller graph for this function:

PrintStepLengthDiagnostic()

void G4PropagatorInField::PrintStepLengthDiagnostic ( G4double  currentProposedStepLength, G4double  decreaseFactor, G4double  stepTrial, const G4FieldTrack &  aFieldTrack  )

protected

Definition at line 610 of file G4PropagatorInField.cc.

{
  G4int  iprec= G4cout.precision(8); 
  G4cout << " " << std::setw(12) << " PiF: NoZeroStep " 
         << " " << std::setw(20) << " CurrentProposed len " 
         << " " << std::setw(18) << " Full_curvelen_last" 
         << " " << std::setw(18) << " last proposed len " 
         << " " << std::setw(18) << " decrease factor   " 
         << " " << std::setw(15) << " step trial  " 
         << G4endl;

  G4cout << " " << std::setw(10) << fNoZeroStep << "  "
         << " " << std::setw(20) << CurrentProposedStepLength
         << " " << std::setw(18) << fFull_CurveLen_of_LastAttempt
         << " " << std::setw(18) << fLast_ProposedStepLength 
         << " " << std::setw(18) << decreaseFactor
         << " " << std::setw(15) << stepTrial
         << G4endl;
  G4cout.precision( iprec ); 

}

Here is the caller graph for this function:

RefreshIntersectionLocator()

void G4PropagatorInField::RefreshIntersectionLocator

(

)

Definition at line 132 of file G4PropagatorInField.cc.

{
  fIntersectionLocator->SetEpsilonStepFor(fEpsilonStep);
  fIntersectionLocator->SetDeltaIntersectionFor(fCurrentFieldMgr->GetDeltaIntersection());
  fIntersectionLocator->SetChordFinderFor(GetChordFinder());
  fIntersectionLocator->SetSafetyParametersFor( fUseSafetyForOptimisation);
}

Here is the call graph for this function:

Here is the caller graph for this function:

ReportLoopingParticle()

void G4PropagatorInField::ReportLoopingParticle ( G4int  count, double  StepTaken, G4VPhysicalVolume *  pPhysVol  )

protected

Definition at line 736 of file G4PropagatorInField.cc.

{
   std::ostringstream message;
   message << "  Killing looping particle " 
           << " after " << count << " field substeps "
           << " totaling " << StepTaken / mm << " mm " ;
   if( pPhysVol )
   {
      message << " in *volume* " << pPhysVol->GetName() ;
   }
   else
   {
      message << " in unknown or null volume. " ;
   }
   G4Exception("G4PropagatorInField::ComputeStep()", "GeomNav1002",
               JustWarning, message);   
}

Here is the call graph for this function:

Here is the caller graph for this function:

ReportStuckParticle()

void G4PropagatorInField::ReportStuckParticle ( G4int  noZeroSteps, G4double  proposedStep, G4double  lastTriedStep, G4VPhysicalVolume *  physVol  )

protected

Definition at line 756 of file G4PropagatorInField.cc.

{
   std::ostringstream message;
   message << "Particle is stuck; it will be killed." << G4endl
           << "  Zero progress for "  << noZeroSteps << " attempted steps." 
           << G4endl
           << "  Proposed Step is " << proposedStep
           << " but Step Taken is "<< lastTriedStep << G4endl;
   if( physVol )
      message << " in volume " << physVol->GetName() ; 
   else
      message << " in unknown or null volume. " ; 
   G4Exception("G4PropagatorInField::ComputeStep()",
               "GeomNav1002", JustWarning, message);
}

Here is the call graph for this function:

Here is the caller graph for this function:

SetDetectorFieldManager()

void G4PropagatorInField::SetDetectorFieldManager ( G4FieldManager *  newGlobalFieldManager )

inline

Here is the caller graph for this function:

SetEpsilonStep()

void G4PropagatorInField::SetEpsilonStep ( G4double  newEps )

inline

Here is the caller graph for this function:

SetIntersectionLocator()

void G4PropagatorInField::SetIntersectionLocator ( G4VIntersectionLocator *  pLocator )

inline

Here is the caller graph for this function:

SetLargestAcceptableStep()

void G4PropagatorInField::SetLargestAcceptableStep ( G4double  newBigDist )

inline

Here is the caller graph for this function:

SetMaximumEpsilonStep()

void G4PropagatorInField::SetMaximumEpsilonStep ( G4double  newEpsMax )

inline

Here is the caller graph for this function:

SetMaxLoopCount()

void G4PropagatorInField::SetMaxLoopCount ( G4int  new_max )

inline

Here is the caller graph for this function:

SetMinimumEpsilonStep()

void G4PropagatorInField::SetMinimumEpsilonStep ( G4double  newEpsMin )

inline

Here is the caller graph for this function:

SetNavigatorForPropagating()

void G4PropagatorInField::SetNavigatorForPropagating ( G4Navigator *  SimpleOrMultiNavigator )

inline

Here is the caller graph for this function:

SetThresholdNoZeroStep()

void G4PropagatorInField::SetThresholdNoZeroStep ( G4int  noAct, G4int  noHarsh, G4int  noAbandon  )

inline

Here is the caller graph for this function:

SetTrajectoryFilter()

void G4PropagatorInField::SetTrajectoryFilter

(

G4VCurvedTrajectoryFilter *

filter

)

Definition at line 664 of file G4PropagatorInField.cc.

{
  fpTrajectoryFilter = filter;
}

Here is the caller graph for this function:

SetUseSafetyForOptimization()

void G4PropagatorInField::SetUseSafetyForOptimization ( G4bool  )

inline

Here is the caller graph for this function:

SetVerboseLevel()

G4int G4PropagatorInField::SetVerboseLevel

(

G4int

verbose

)

Definition at line 721 of file G4PropagatorInField.cc.

{
  G4int oldval= fVerboseLevel;
  fVerboseLevel= level;

  // Forward the verbose level 'reduced' to ChordFinder,
  // MagIntegratorDriver ... ? 
  //
  G4MagInt_Driver* integrDriver= GetChordFinder()->GetIntegrationDriver(); 
  integrDriver->SetVerboseLevel( fVerboseLevel - 2 );
  G4cout << "Set Driver verbosity to " << fVerboseLevel - 2 << G4endl;

  return oldval;
}

Here is the call graph for this function:

SetVerboseTrace()

void G4PropagatorInField::SetVerboseTrace ( G4bool  enable )

inline

Definition at line 102 of file G4PropagatorInField.hh.

{ fVerbTracePiF = enable; }

SetZeroStepThreshold()

void G4PropagatorInField::SetZeroStepThreshold ( G4double  newLength )

inline

Here is the caller graph for this function:

Verbose()

G4int G4PropagatorInField::Verbose ( ) const

inline

Member Data Documentation

End_PointAndTangent

G4FieldTrack G4PropagatorInField::End_PointAndTangent

private

Definition at line 266 of file G4PropagatorInField.hh.

fAbandonThreshold_NoZeroSteps

G4int G4PropagatorInField::fAbandonThreshold_NoZeroSteps

private

Definition at line 221 of file G4PropagatorInField.hh.

fActionThreshold_NoZeroSteps

G4int G4PropagatorInField::fActionThreshold_NoZeroSteps

private

Definition at line 219 of file G4PropagatorInField.hh.

fAllocatedLocator

G4bool G4PropagatorInField::fAllocatedLocator

private

Definition at line 232 of file G4PropagatorInField.hh.

fCurrentFieldMgr

G4FieldManager* G4PropagatorInField::fCurrentFieldMgr

private

Definition at line 260 of file G4PropagatorInField.hh.

fDetectorFieldMgr

G4FieldManager* G4PropagatorInField::fDetectorFieldMgr

private

Definition at line 237 of file G4PropagatorInField.hh.

fEpsilonStep

G4double G4PropagatorInField::fEpsilonStep

private

Definition at line 265 of file G4PropagatorInField.hh.

fFirstStepInVolume

G4bool G4PropagatorInField::fFirstStepInVolume

private

Definition at line 282 of file G4PropagatorInField.hh.

fFull_CurveLen_of_LastAttempt

G4double G4PropagatorInField::fFull_CurveLen_of_LastAttempt

private

Definition at line 271 of file G4PropagatorInField.hh.

fIntersectionLocator

G4VIntersectionLocator* G4PropagatorInField::fIntersectionLocator

private

Definition at line 240 of file G4PropagatorInField.hh.

fLargestAcceptableStep

G4double G4PropagatorInField::fLargestAcceptableStep

private

Definition at line 225 of file G4PropagatorInField.hh.

fLast_ProposedStepLength

G4double G4PropagatorInField::fLast_ProposedStepLength

private

Definition at line 272 of file G4PropagatorInField.hh.

fLastStepInVolume

G4bool G4PropagatorInField::fLastStepInVolume

private

Definition at line 283 of file G4PropagatorInField.hh.

fMax_loop_count

G4int G4PropagatorInField::fMax_loop_count

private

Definition at line 214 of file G4PropagatorInField.hh.

fNavigator

G4Navigator* G4PropagatorInField::fNavigator

private

Definition at line 250 of file G4PropagatorInField.hh.

fNewTrack

G4bool G4PropagatorInField::fNewTrack

private

Definition at line 284 of file G4PropagatorInField.hh.

fNoZeroStep

G4int G4PropagatorInField::fNoZeroStep

private

Definition at line 268 of file G4PropagatorInField.hh.

fParticleIsLooping

G4bool G4PropagatorInField::fParticleIsLooping

private

Definition at line 267 of file G4PropagatorInField.hh.

fPreviousSafety

G4double G4PropagatorInField::fPreviousSafety

private

Definition at line 275 of file G4PropagatorInField.hh.

fPreviousSftOrigin

G4ThreeVector G4PropagatorInField::fPreviousSftOrigin

private

Definition at line 274 of file G4PropagatorInField.hh.

fpTrajectoryFilter

G4VCurvedTrajectoryFilter* G4PropagatorInField::fpTrajectoryFilter

private

Definition at line 243 of file G4PropagatorInField.hh.

fSetFieldMgr

G4bool G4PropagatorInField::fSetFieldMgr

private

Definition at line 262 of file G4PropagatorInField.hh.

fSevereActionThreshold_NoZeroSteps

G4int G4PropagatorInField::fSevereActionThreshold_NoZeroSteps

private

Definition at line 220 of file G4PropagatorInField.hh.

fUseSafetyForOptimisation

G4bool G4PropagatorInField::fUseSafetyForOptimisation

private

Definition at line 216 of file G4PropagatorInField.hh.

fVerboseLevel

G4int G4PropagatorInField::fVerboseLevel

private

Definition at line 278 of file G4PropagatorInField.hh.

fVerbTracePiF

G4bool G4PropagatorInField::fVerbTracePiF

private

Definition at line 279 of file G4PropagatorInField.hh.

fZeroStepThreshold

G4double G4PropagatorInField::fZeroStepThreshold

private

Definition at line 222 of file G4PropagatorInField.hh.

kCarTolerance

G4double G4PropagatorInField::kCarTolerance

private

Definition at line 229 of file G4PropagatorInField.hh.

---

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

• G4PropagatorInField.hh
• G4PropagatorInField.cc