#include <G4ErrorPropagationNavigator.hh>

Inheritance diagram for G4ErrorPropagationNavigator:

Collaboration diagram for G4ErrorPropagationNavigator:

Public Member Functions
	G4ErrorPropagationNavigator ()

	~G4ErrorPropagationNavigator ()

G4double	ComputeStep (const G4ThreeVector &pGlobalPoint, const G4ThreeVector &pDirection, const G4double pCurrentProposedStepLength, G4double &pNewSafety)

G4double	ComputeSafety (const G4ThreeVector &globalpoint, const G4double pProposedMaxLength=DBL_MAX, const G4bool keepState=true)

G4ThreeVector	GetGlobalExitNormal (const G4ThreeVector &point, G4bool *valid)

G4double	TargetSafetyFromPoint (const G4ThreeVector &pGlobalpoint)

Public Member Functions inherited from G4Navigator
	G4Navigator ()

virtual	~G4Navigator ()

G4double	CheckNextStep (const G4ThreeVector &pGlobalPoint, const G4ThreeVector &pDirection, const G4double pCurrentProposedStepLength, G4double &pNewSafety)

virtual G4VPhysicalVolume *	ResetHierarchyAndLocate (const G4ThreeVector &point, const G4ThreeVector &direction, const G4TouchableHistory &h)

virtual G4VPhysicalVolume *	LocateGlobalPointAndSetup (const G4ThreeVector &point, const G4ThreeVector *direction=0, const G4bool pRelativeSearch=true, const G4bool ignoreDirection=true)

virtual void	LocateGlobalPointWithinVolume (const G4ThreeVector &position)

void	LocateGlobalPointAndUpdateTouchableHandle (const G4ThreeVector &position, const G4ThreeVector &direction, G4TouchableHandle &oldTouchableToUpdate, const G4bool RelativeSearch=true)

void	LocateGlobalPointAndUpdateTouchable (const G4ThreeVector &position, const G4ThreeVector &direction, G4VTouchable *touchableToUpdate, const G4bool RelativeSearch=true)

void	LocateGlobalPointAndUpdateTouchable (const G4ThreeVector &position, G4VTouchable *touchableToUpdate, const G4bool RelativeSearch=true)

void	SetGeometricallyLimitedStep ()

virtual G4bool	RecheckDistanceToCurrentBoundary (const G4ThreeVector &pGlobalPoint, const G4ThreeVector &pDirection, const G4double CurrentProposedStepLength, G4double prDistance, G4double prNewSafety=0) const

G4VPhysicalVolume *	GetWorldVolume () const

void	SetWorldVolume (G4VPhysicalVolume *pWorld)

G4GRSVolume *	CreateGRSVolume () const

G4GRSSolid *	CreateGRSSolid () const

G4TouchableHistory *	CreateTouchableHistory () const

G4TouchableHistory *	CreateTouchableHistory (const G4NavigationHistory *) const

virtual G4TouchableHistoryHandle	CreateTouchableHistoryHandle () const

virtual G4ThreeVector	GetLocalExitNormal (G4bool *valid)

virtual G4ThreeVector	GetLocalExitNormalAndCheck (const G4ThreeVector &point, G4bool *valid)

G4int	GetVerboseLevel () const

void	SetVerboseLevel (G4int level)

G4bool	IsActive () const

void	Activate (G4bool flag)

G4bool	EnteredDaughterVolume () const

G4bool	ExitedMotherVolume () const

void	CheckMode (G4bool mode)

G4bool	IsCheckModeActive () const

void	SetPushVerbosity (G4bool mode)

void	PrintState () const

const G4AffineTransform &	GetGlobalToLocalTransform () const

const G4AffineTransform	GetLocalToGlobalTransform () const

G4AffineTransform	GetMotherToDaughterTransform (G4VPhysicalVolume *dVolume, G4int dReplicaNo, EVolume dVolumeType)

void	ResetStackAndState ()

G4int	SeverityOfZeroStepping (G4int *noZeroSteps) const

G4ThreeVector	GetCurrentLocalCoordinate () const

G4ThreeVector	NetTranslation () const

G4RotationMatrix	NetRotation () const

void	EnableBestSafety (G4bool value=false)

Additional Inherited Members
Protected Member Functions inherited from G4Navigator
void	SetSavedState ()

void	RestoreSavedState ()

virtual void	ResetState ()

G4ThreeVector	ComputeLocalPoint (const G4ThreeVector &rGlobPoint) const

G4ThreeVector	ComputeLocalAxis (const G4ThreeVector &pVec) const

EVolume	VolumeType (const G4VPhysicalVolume *pVol) const

EVolume	CharacteriseDaughters (const G4LogicalVolume *pLog) const

G4int	GetDaughtersRegularStructureId (const G4LogicalVolume *pLog) const

virtual void	SetupHierarchy ()

Protected Attributes inherited from G4Navigator
G4double	kCarTolerance

G4double	fMinStep

G4double	fSqTol

G4NavigationHistory	fHistory

G4bool	fEnteredDaughter

G4bool	fExitedMother

G4bool	fWasLimitedByGeometry

G4ThreeVector	fStepEndPoint

G4ThreeVector	fLastStepEndPointLocal

G4int	fVerbose

Detailed Description

Definition at line 50 of file G4ErrorPropagationNavigator.hh.

Constructor & Destructor Documentation

G4ErrorPropagationNavigator::G4ErrorPropagationNavigator ( )

Definition at line 47 of file G4ErrorPropagationNavigator.cc.

   : G4Navigator()
 {
 }

G4ErrorPropagationNavigator::~G4ErrorPropagationNavigator ( )

Definition at line 54 of file G4ErrorPropagationNavigator.cc.

55 {

56 }

Member Function Documentation

G4double G4ErrorPropagationNavigator::ComputeSafety	(	const G4ThreeVector &	globalpoint,
		const G4double	pProposedMaxLength = `DBL_MAX`,
		const G4bool	keepState = `true`
	)

virtual

Reimplemented from G4Navigator.

Definition at line 155 of file G4ErrorPropagationNavigator.cc.

 {
   G4double safetyGeom = G4Navigator::ComputeSafety(pGlobalPoint,
                                                   pMaxLength, keepState);
 
   G4double safetyTarget = TargetSafetyFromPoint( pGlobalPoint ); 
 
   return std::min(safetyGeom, safetyTarget); 
 }

Here is the call graph for this function:

G4double G4ErrorPropagationNavigator::ComputeStep	(	const G4ThreeVector &	pGlobalPoint,
		const G4ThreeVector &	pDirection,
		const G4double	pCurrentProposedStepLength,
		G4double &	pNewSafety
	)

virtual

Reimplemented from G4Navigator.

Definition at line 61 of file G4ErrorPropagationNavigator.cc.

 {
   G4double safetyGeom= DBL_MAX;
    
   G4double Step = G4Navigator::ComputeStep(pGlobalPoint, pDirection,
                                            pCurrentProposedStepLength,
                                            safetyGeom);
   
   G4ErrorPropagatorData * g4edata
     = G4ErrorPropagatorData::GetErrorPropagatorData();
 
   if (g4edata !=0)
   {
     const G4ErrorTarget* target = g4edata->GetTarget();
     if( target != 0 )
     {
       G4double StepPlane= target->GetDistanceFromPoint(pGlobalPoint,pDirection);
 
       if( StepPlane < 0. ) // Negative means target is crossed, will not be found
       {
         StepPlane = DBL_MAX; 
       }
 #ifdef G4VERBOSE
       if( G4ErrorPropagatorData::verbose() >= 4 )
       {
         G4cout << "G4ErrorPropagationNavigator::ComputeStep()" << G4endl
                << "  Target step: " << StepPlane
                << ", Transportation step: " << Step << G4endl;
         target->Dump( "G4ErrorPropagationNavigator::ComputeStep Target " );
       }
 #endif
 
       if(StepPlane<Step)
       {
 #ifdef G4VERBOSE
         if( G4ErrorPropagatorData::verbose() >= 2 )
         {
           G4cout << "G4ErrorPropagationNavigator::ComputeStep()" << G4endl
                  << "  TargetCloserThanBoundary: " << StepPlane << " < "
                  << Step << G4endl;
         }
 #endif
         Step = StepPlane;
         g4edata->SetState(G4ErrorState_TargetCloserThanBoundary);
       }
       else
       {
         g4edata->SetState(G4ErrorState_Propagating);
       }
     }
   }
   G4double safetyTarget = TargetSafetyFromPoint(pGlobalPoint);
    // Avoid call to G4Navigator::ComputeSafety - which could have side effects
   pNewSafety= std::min(safetyGeom, safetyTarget); 
 
 #ifdef G4VERBOSE
   if( G4ErrorPropagatorData::verbose() >= 3 )
   {
     G4cout << "G4ErrorPropagationNavigator::ComputeStep()" << G4endl
            << "  Step: " << Step << ", ComputeSafety: " << pNewSafety
            << G4endl;
   }
 #endif
 
   return Step;
 }

Here is the call graph for this function:

G4ThreeVector G4ErrorPropagationNavigator::GetGlobalExitNormal	(	const G4ThreeVector &	point,
		G4bool *	valid
	)

virtual

Reimplemented from G4Navigator.

Definition at line 170 of file G4ErrorPropagationNavigator.cc.

 {
   G4ErrorPropagatorData *g4edata
         = G4ErrorPropagatorData::GetErrorPropagatorData();
   const G4ErrorTarget* target = 0;
 
   G4ThreeVector normal(0.0, 0.0, 0.0);
   G4double      distance= 0;
   
   // Determine which 'geometry' limited the step
   if (g4edata)
   {
     target = g4edata->GetTarget();
     if(target)
     {
       distance = target->GetDistanceFromPoint(point);
     }
   }
   
   if( distance > kCarTolerance   // Not reached the target.
      || (!target) )
             //  If a target does not exist, this seems the best we can do
   {
     normal= G4Navigator::GetGlobalExitNormal(point, valid);
   }
   else
   {
     switch( target->GetType() )
     {
       case G4ErrorTarget_GeomVolume:
         // The volume is in the 'real' mass geometry
         normal= G4Navigator::GetGlobalExitNormal(point, valid);
         break;
       case G4ErrorTarget_TrkL:
         normal= G4ThreeVector( 0.0, 0.0, 0.0);
         *valid= false;
         G4Exception("G4ErrorPropagationNavigator::GetGlobalExitNormal",
                     "Geometry1003",
                     JustWarning, "Unexpected value of Target type");
         break;
       case G4ErrorTarget_PlaneSurface:
       case G4ErrorTarget_CylindricalSurface:
         const G4ErrorSurfaceTarget* surfaceTarget=
           static_cast<const G4ErrorSurfaceTarget*>(target);
         normal= surfaceTarget->GetTangentPlane(point).normal().unit();
         *valid= true;
         break;
 
 //      default:
 //        normal= G4ThreeVector( 0.0, 0.0, 0.0);
 //        *valid= false;
 //        G4Exception("G4ErrorPropagationNavigator::GetGlobalExitNormal",
 //                    "Geometry:003",
 //                    FatalException, "Impossible value of Target type");
 //        exit(1);
 //        break;
     }
   }
   return normal;
 }

Here is the call graph for this function:

G4double G4ErrorPropagationNavigator::TargetSafetyFromPoint ( const G4ThreeVector & pGlobalpoint )

Definition at line 134 of file G4ErrorPropagationNavigator.cc.

 {
   G4double safety= DBL_MAX;
 
   G4ErrorPropagatorData *g4edata
     = G4ErrorPropagatorData::GetErrorPropagatorData();
 
   if (g4edata !=0)
   {
     const G4ErrorTarget* target = g4edata->GetTarget();
     if( target != 0 )
     {
        safety = target->GetDistanceFromPoint(pGlobalpoint);
     }
   }
   return safety;
 }

Here is the call graph for this function:

Here is the caller graph for this function:

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

source/geant4.10.03.p03/source/geometry/navigation/include/G4ErrorPropagationNavigator.hh
source/geant4.10.03.p03/source/geometry/navigation/src/G4ErrorPropagationNavigator.cc

Public Member Functions

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation