53 #define G4DEBUG_NAVIGATION 1
62 fTopPhysical(0), fCheck(false),
fPushed(false), fWarnPush(true)
162 const G4bool relativeSearch,
163 const G4bool ignoreDirection )
165 G4bool notKnownContained=
true, noResult;
177 if( considerDirection && pGlobalDirection != 0 )
179 globalDirection=*pGlobalDirection;
187 G4cout <<
"*** G4ITNavigator1::LocateGlobalPointAndSetup: ***" <<
G4endl;
188 G4cout <<
" Called with arguments: " << G4endl
189 <<
" Globalpoint = " << globalPoint << G4endl
190 <<
" RelativeSearch = " << relativeSearch <<
G4endl;
196 G4cout.precision(oldcoutPrec);
200 if ( !relativeSearch )
283 notKnownContained =
false;
304 G4int noLevelsExited=0 ;
306 while (notKnownContained)
312 insideCode = targetSolid->
Inside(localPoint);
316 G4String solidResponse =
"-kInside-";
318 solidResponse =
"-kOutside-";
320 solidResponse =
"-kSurface-";
321 G4cout <<
"*** G4ITNavigator1::LocateGlobalPointAndSetup(): ***" <<
G4endl
322 <<
" Invoked Inside() for solid: " << targetSolid->
GetName()
323 <<
". Solid replied: " << solidResponse <<
G4endl
324 <<
" For local point p: " << localPoint <<
G4endl;
349 if( noLevelsExited > 1 )
373 if( (!
fExiting)&&considerDirection )
378 G4bool directionExiting =
false;
389 directionExiting = normal.dot(localDirection) > 0.0;
390 isExiting = isExiting || directionExiting;
406 if( noLevelsExited > 1 )
428 notKnownContained=
false;
433 notKnownContained=
false;
454 if (!targetPhysical) {
break; }
541 #ifdef G4DEBUG_NAVIGATION
545 G4cout <<
"*** G4ITNavigator1::LocateGlobalPointAndSetup() ***" <<
G4endl;
560 if (targetPhysical) { curPhysVol_Name = targetPhysical->
GetName(); }
561 G4cout <<
" Return value = new volume = " << curPhysVol_Name <<
G4endl;
566 G4cout <<
"Upon exiting LocateGlobalPointAndSetup():" <<
G4endl;
569 G4cout.precision(oldcoutPrec);
575 return targetPhysical;
598 #ifdef G4DEBUG_NAVIGATION
601 G4cout <<
"Entering LocateGlobalWithinVolume(): History = " <<
G4endl;
635 G4Exception(
"G4ITNavigator1::LocateGlobalPointWithinVolume()",
637 "Not applicable for replicated volumes.");
789 const G4double pCurrentProposedStepLength,
815 G4cout <<
"*** G4ITNavigator1::ComputeStep: ***" <<
G4endl;
817 <<
" - Proposed step length = " << pCurrentProposedStepLength
819 #ifdef G4DEBUG_NAVIGATION
822 G4cout <<
" Called with the arguments: " << G4endl
823 <<
" Globalpoint = " << std::setw(25) << pGlobalpoint << G4endl
824 <<
" Direction = " << std::setw(25) << pDirection <<
G4endl;
841 G4double moveLenSq = (newLocalPoint-oldLocalPoint).mag2();
863 pCurrentProposedStepLength,
880 pCurrentProposedStepLength,
908 "Point is relocated in voxels, while it should be outside!");
911 pCurrentProposedStepLength,
926 pCurrentProposedStepLength,
945 pCurrentProposedStepLength,
959 pCurrentProposedStepLength,
971 G4Exception(
"G4ITNavigator1::ComputeStep()",
"GeomNav0001",
982 G4bool calculatedExitNormal=
false;
988 pCurrentProposedStepLength,
992 calculatedExitNormal,
1026 #ifdef G4DEBUG_NAVIGATION
1029 G4cout <<
"G4ITNavigator1::ComputeStep(): another zero step, # "
1031 <<
" at " << pGlobalpoint
1032 <<
" in volume " << motherPhysical->
GetName()
1033 <<
" nav-comp-step calls # " << sNavCScalls
1045 std::ostringstream message;
1046 message <<
"Track stuck or not moving." <<
G4endl
1047 <<
" Track stuck, not moving for "
1049 <<
" in volume -" << motherPhysical->
GetName()
1050 <<
"- at point " << pGlobalpoint <<
G4endl
1051 <<
" direction: " << pDirection <<
"." <<
G4endl
1052 <<
" Potential geometry or navigation problem !"
1054 <<
" Trying pushing it of " << Step <<
" mm ...";
1055 G4Exception(
"G4ITNavigator1::ComputeStep()",
"GeomNav1002",
1056 JustWarning, message,
"Potential overlap in geometry!");
1065 std::ostringstream message;
1066 message <<
"Stuck Track: potential geometry or navigation problem."
1068 <<
" Track stuck, not moving for "
1070 <<
" in volume -" << motherPhysical->
GetName()
1071 <<
"- at point " << pGlobalpoint <<
G4endl
1072 <<
" direction: " << pDirection <<
".";
1074 G4Exception(
"G4ITNavigator1::ComputeStep()",
"GeomNav0003",
1091 #ifdef G4DEBUG_NAVIGATION
1094 G4cout <<
" At G4Nav CompStep End - if(exiting) - fExiting= " <<
fExiting
1155 #ifdef G4DEBUG_NAVIGATION
1158 desc <<
"Problem in ComputeStep: Replica Navigation did not provide"
1159 <<
" valid exit Normal. " <<
G4endl;
1160 desc <<
" Do not know how calculate it in this case." <<
G4endl;
1161 desc <<
" Location = " << finalLocalPoint <<
G4endl;
1162 desc <<
" Volume name = " << motherPhysical->
GetName()
1164 G4Exception(
"G4ITNavigator1::ComputeStep()",
"GeomNav0003",
1165 JustWarning, desc,
"Normal not available for exiting.");
1210 G4cout <<
" Returned step= " << Step;
1214 G4cout <<
" Requested step= " << pCurrentProposedStepLength ;
1232 const G4double pCurrentProposedStepLength,
1243 pCurrentProposedStepLength,
1306 for ( i=1; i<=cdepth; i++ )
1335 ComputeMaterial(replicaNo, current, &touchable) );
1362 if( candidateLogical )
1384 currentSolid= candidateLogical->
GetSolid();
1385 inSideIt = currentSolid->
Inside(daughterPointOwnLocal);
1386 onSurface = (inSideIt ==
kSurface);
1391 safety = (currentSolid->
DistanceToIn(daughterPointOwnLocal));
1394 else if (inSideIt ==
kInside )
1396 safety = (currentSolid->
DistanceToOut(daughterPointOwnLocal));
1403 nextSolidExitNormal =
1408 ExitNormal = -nextSolidExitNormal;
1416 std::ostringstream message;
1417 message <<
"Point not on surface ! " <<
G4endl
1419 << daughterPointOwnLocal <<
G4endl
1420 <<
" Physical volume = "
1422 <<
" Logical volume = "
1424 <<
" Solid = " << currentSolid->
GetName()
1427 << *currentSolid <<
G4endl;
1430 message <<
"Point is Outside. " << G4endl
1431 <<
" Safety (from outside) = " << safety <<
G4endl;
1435 message <<
"Point is Inside. " << G4endl
1436 <<
" Safety (from inside) = " << safety <<
G4endl;
1438 G4Exception(
"G4ITNavigator1::GetLocalExitNormal()",
"GeomNav1001",
1456 G4Exception(
"G4ITNavigator1::GetLocalExitNormal()",
1458 "Incorrect call to GetLocalSurfaceNormal." );
1471 desc <<
" Parameters of solid: " << *daughterSolid
1473 G4Exception(
"G4ITNavigator1::GetLocalExitNormal()",
1475 "Surface Normal returned by Solid is not a Unit Vector." );
1493 message <<
"Function called when *NOT* at a Boundary." <<
G4endl;
1494 G4Exception(
"G4ITNavigator1::GetLocalExitNormal()",
1510 G4int enteringReplicaNo,
1513 switch (enteringVolumeType)
1518 G4Exception(
"G4ITNavigator1::GetMotherToDaughterTransform()",
1520 "Method NOT Implemented yet for replica volumes.");
1528 pParam->
ComputeSolid(enteringReplicaNo, pEnteringPhysVol);
1563 #ifdef G4DEBUG_NAVIGATION
1568 G4ThreeVector ExpectedBoundaryPointLocal;
1571 ExpectedBoundaryPointLocal =
1591 G4bool* pNormalCalculated)
1601 *pNormalCalculated =
true;
1609 #ifdef G4DEBUG_NAVIGATION
1617 edN <<
" State of Navigator: " <<
G4endl;
1621 G4Exception(
"G4ITNavigator1::GetGlobalExitNormal()",
1623 "LocalExitNormalAndCheck() did not calculate Normal.");
1627 G4double localMag2= localNormal.mag2();
1632 edN <<
"G4ITNavigator1::GetGlobalExitNormal: "
1633 <<
" Using Local Normal - from call to GetLocalExitNormalAndCheck. "
1635 <<
" Local Exit Normal = " << localNormal <<
" || = "
1636 << std::sqrt(localMag2) <<
G4endl
1637 <<
" Global Exit Normal = " << globalNormal <<
" || = "
1638 << globalNormal.mag() <<
G4endl;
1641 G4Exception(
"G4ITNavigator1::GetGlobalExitNormal()",
1643 "Value obtained from new local *solid* is incorrect.");
1644 localNormal = localNormal.unit();
1650 #ifdef G4DEBUG_NAVIGATION
1665 edDfn <<
"Found difference in normals in case of exiting mother "
1666 <<
"- when Get is called after ComputingStep " <<
G4endl;
1667 edDfn <<
" Magnitude of diff = " << diffNorm.mag() <<
G4endl;
1668 edDfn <<
" Normal stored (Global) = " << fExitNormalGlobalFrame
1670 edDfn <<
" Global Computed from Local = " << globalNormal <<
G4endl;
1671 G4Exception(
"G4ITNavigator1::GetGlobalExitNormal()",
"GeomNav0003",
1677 return globalNormal;
1681 #define G4NEW_SAFETY 1
1698 #ifdef G4DEBUG_NAVIGATION
1702 G4cout <<
"*** G4ITNavigator1::ComputeSafety: ***" <<
G4endl
1703 <<
" Called at point: " << pGlobalpoint <<
G4endl;
1707 <<
" - Maximum length = " << pMaxLength <<
G4endl;
1710 G4cout <<
" ----- Upon entering Compute Safety:" <<
G4endl;
1722 if( !(endpointOnSurface && stayedOnEndpoint) )
1735 #ifdef G4DEBUG_NAVIGATION
1738 G4cout <<
" G4ITNavigator1::ComputeSafety() relocates-in-volume to point: "
1739 << pGlobalpoint <<
G4endl;
1756 *motherPhysical, pMaxLength);
1757 newSafety= safetyTwo;
1762 newSafety= safetyOldVoxel;
1781 G4Exception(
"G4ITNavigator1::ComputeSafety()",
"GeomNav0001",
1794 #ifdef G4DEBUG_NAVIGATION
1797 G4cout <<
" G4ITNavigator1::ComputeSafety() finds that point - "
1798 << pGlobalpoint <<
" - is on surface " <<
G4endl;
1815 #ifdef G4DEBUG_NAVIGATION
1820 G4cout <<
" Returned value of Safety = " << newSafety <<
G4endl;
1822 G4cout.precision(oldcoutPrec);
1846 G4cout <<
"The current state of G4ITNavigator1 is: " <<
G4endl;
1849 <<
" Exiting = " <<
fExiting << G4endl
1850 <<
" Entering = " <<
fEntering << G4endl
1851 <<
" BlockedPhysicalVolume= " ;
1864 G4cout << std::setw(30) <<
" ExitNormal " <<
" "
1865 << std::setw( 5) <<
" Valid " <<
" "
1866 << std::setw( 9) <<
" Exiting " <<
" "
1867 << std::setw( 9) <<
" Entering" <<
" "
1868 << std::setw(15) <<
" Blocked:Volume " <<
" "
1869 << std::setw( 9) <<
" ReplicaNo" <<
" "
1870 << std::setw( 8) <<
" LastStepZero " <<
" "
1874 <<
", " << std::setw(7) <<
fExitNormal.z() <<
" ) "
1876 << std::setw( 9) <<
fExiting <<
" "
1879 G4cout << std::setw(15) <<
"None";
1893 G4cout.precision(oldcoutPrec);
1922 G4double shiftOrigin = std::sqrt(shiftOriginSafSq);
1925 if( diffShiftSaf > fAccuracyForWarning )
1929 std::ostringstream message, suggestion;
1930 message <<
"Accuracy error or slightly inaccurate position shift."
1932 <<
" The Step's starting point has moved "
1933 << std::sqrt(moveLenSq)/
mm <<
" mm " <<
G4endl
1934 <<
" since the last call to a Locate method." <<
G4endl
1935 <<
" This has resulted in moving "
1936 << shiftOrigin/
mm <<
" mm "
1937 <<
" from the last point at which the safety "
1938 <<
" was calculated " <<
G4endl
1939 <<
" which is more than the computed safety= "
1940 << fPreviousSafety/
mm <<
" mm at that point." <<
G4endl
1941 <<
" This difference is "
1942 << diffShiftSaf/
mm <<
" mm." <<
G4endl
1943 <<
" The tolerated accuracy is "
1944 << fAccuracyForException/
mm <<
" mm.";
1948 if( ((++warnNow % 100) == 1) )
1951 <<
" This problem can be due to either " <<
G4endl
1952 <<
" - a process that has proposed a displacement"
1953 <<
" larger than the current safety , or" <<
G4endl
1954 <<
" - inaccuracy in the computation of the safety";
1955 suggestion <<
"We suggest that you " << G4endl
1956 <<
" - find i) what particle is being tracked, and "
1957 <<
" ii) through what part of your geometry " << G4endl
1958 <<
" for example by re-running this event with "
1960 <<
" /tracking/verbose 1 " << G4endl
1961 <<
" - check which processes you declare for"
1962 <<
" this particle (and look at non-standard ones)"
1964 <<
" - in case, create a detailed logfile"
1965 <<
" of this event using:" << G4endl
1966 <<
" /tracking/verbose 6 ";
1970 message,
G4String(suggestion.str()));
1971 G4cout.precision(oldcoutPrec);
1972 G4cerr.precision(oldcerrPrec);
1974 #ifdef G4DEBUG_NAVIGATION
1977 G4cerr <<
"WARNING - G4ITNavigator1::ComputeStep()" <<
G4endl
1978 <<
" The Step's starting point has moved "
1979 << std::sqrt(moveLenSq) <<
"," <<
G4endl
1980 <<
" which has taken it to the limit of"
1981 <<
" the current safety. " <<
G4endl;
1986 if ( shiftOriginSafSq >
sqr(safetyPlus) )
1988 std::ostringstream message;
1989 message <<
"May lead to a crash or unreliable results." <<
G4endl
1990 <<
" Position has shifted considerably without"
1991 <<
" notifying the navigator !" <<
G4endl
1992 <<
" Tolerated safety: " << safetyPlus <<
G4endl
1993 <<
" Computed shift : " << shiftOriginSafSq;
1994 G4Exception(
"G4ITNavigator1::ComputeStep()",
"GeomNav1002",
2011 G4int oldcoutPrec = os.precision(4);
2014 os <<
"The current state of G4ITNavigator1 is: " <<
G4endl;
2017 <<
" Exiting = " << n.
fExiting << G4endl
2018 <<
" Entering = " << n.
fEntering << G4endl
2019 <<
" BlockedPhysicalVolume= " ;
2032 os << std::setw(30) <<
" ExitNormal " <<
" "
2033 << std::setw( 5) <<
" Valid " <<
" "
2034 << std::setw( 9) <<
" Exiting " <<
" "
2035 << std::setw( 9) <<
" Entering" <<
" "
2036 << std::setw(15) <<
" Blocked:Volume " <<
" "
2037 << std::setw( 9) <<
" ReplicaNo" <<
" "
2038 << std::setw( 8) <<
" LastStepZero " <<
" "
2042 <<
", " << std::setw(7) << n.
fExitNormal.z() <<
" ) "
2044 << std::setw( 9) << n.
fExiting <<
" "
2047 { os << std::setw(15) <<
"None"; }
2066 os.precision(oldcoutPrec);
G4int fActionThreshold_NoZeroSteps
G4SmartVoxelHeader * GetVoxelHeader() const
friend std::ostream & operator<<(std::ostream &os, const G4ITNavigator1 &n)
virtual void ResetState()
G4VPhysicalVolume * GetTopVolume() const
virtual void LocateGlobalPointWithinVolume(const G4ThreeVector &position)
G4TouchableHistory * CreateTouchableHistory() const
virtual G4bool LevelLocate(G4NavigationHistory &history, const G4VPhysicalVolume *blockedVol, const G4int blockedNum, const G4ThreeVector &globalPoint, const G4ThreeVector *globalDirection, const G4bool pLocatedOnEdge, G4ThreeVector &localPoint)
G4VoxelNavigation fvoxelNav
G4double ComputeStep(const G4ThreeVector &globalPoint, const G4ThreeVector &globalDirection, const G4double currentProposedStepLength, G4double &newSafety, G4NavigationHistory &history, G4bool &validExitNormal, G4ThreeVector &exitNormal, G4bool &exiting, G4bool &entering, G4VPhysicalVolume *(*pBlockedPhysical), G4int &blockedReplicaNo)
G4bool fChangedGrandMotherRefFrame
const G4ThreeVector & GetTranslation() const
static const G4double kInfinity
G4SaveNavigatorState * fpSaveState
#define fWasLimitedByGeometry
std::ostringstream G4ExceptionDescription
void UpdateMaterial(G4Material *pMaterial)
G4int GetDaughtersRegularStructureId(const G4LogicalVolume *pLog) const
CLHEP::Hep3Vector G4ThreeVector
CLHEP::HepRotation G4RotationMatrix
G4bool fLastTriedStepComputation
virtual void ComputeDimensions(G4VPVParameterisation *p, const G4int n, const G4VPhysicalVolume *pRep)
virtual G4VPhysicalVolume * LocateGlobalPointAndSetup(const G4ThreeVector &point, const G4ThreeVector *direction=0, const G4bool pRelativeSearch=true, const G4bool ignoreDirection=true)
G4ThreeVector fExitNormal
virtual G4VSolid * ComputeSolid(const G4int, G4VPhysicalVolume *)
G4double GetSurfaceTolerance() const
G4VoxelSafety * fpVoxelSafety
G4ThreeVector ComputeLocalPoint(const G4ThreeVector &rGlobPoint) const
G4double ComputeSafety(const G4ThreeVector &localPoint, const G4VPhysicalVolume ¤tPhysical, G4double maxLength=DBL_MAX)
EVolume VolumeType(const G4VPhysicalVolume *pVol) const
virtual G4ThreeVector GetLocalExitNormal(G4bool *valid)
G4bool sLocatedOutsideWorld
void SetSolid(G4VSolid *pSolid)
G4SmartVoxelNode * ParamVoxelLocate(G4SmartVoxelHeader *pHead, const G4ThreeVector &localPoint)
EVolume GetVolumeType(G4int n) const
G4ThreeVector fLastStepEndPointLocal
void SetVerboseLevel(G4int level)
G4double ComputeStep(const G4ThreeVector &globalPoint, const G4ThreeVector &globalDirection, const G4ThreeVector &localPoint, const G4ThreeVector &localDirection, const G4double currentProposedStepLength, G4double &newSafety, G4NavigationHistory &history, G4bool &validExitNormal, G4bool &calculatedExitNormal, G4ThreeVector &exitNormal, G4bool &exiting, G4bool &entering, G4VPhysicalVolume *(*pBlockedPhysical), G4int &blockedReplicaNo)
virtual G4GeometryType GetEntityType() const =0
G4ReferenceCountedHandle< G4TouchableHistory > G4TouchableHistoryHandle
G4bool LevelLocate(G4NavigationHistory &history, const G4VPhysicalVolume *blockedVol, const G4int blockedNum, const G4ThreeVector &globalPoint, const G4ThreeVector *globalDirection, const G4bool pLocatedOnEdge, G4ThreeVector &localPoint)
void ComputeTransformation(const G4int replicaNo, G4VPhysicalVolume *pVol, G4ThreeVector &point) const
G4NormalNavigation fnormalNav
void ResetStackAndState()
G4int GetTopReplicaNo() const
virtual G4double ComputeStep(const G4ThreeVector &pGlobalPoint, const G4ThreeVector &pDirection, const G4double pCurrentProposedStepLength, G4double &pNewSafety)
virtual void SetupHierarchy()
G4double ComputeSafety(const G4ThreeVector &globalpoint, const G4NavigationHistory &history, const G4double pMaxLength=DBL_MAX)
static double normal(HepRandomEngine *eptr)
G4int GetVerboseLevel() const
virtual G4VPhysicalVolume * ResetHierarchyAndLocate(const G4ThreeVector &point, const G4ThreeVector &direction, const G4TouchableHistory &h)
G4int fAbandonThreshold_NoZeroSteps
void SetNormalNavigation(G4NormalNavigation *fnormnav)
G4double ComputeSafety(const G4ThreeVector &localPoint, const G4NavigationHistory &history, const G4double pProposedMaxLength=DBL_MAX)
G4GLOB_DLL std::ostream G4cout
G4ThreeVector ComputeLocalAxis(const G4ThreeVector &pVec) const
G4RegularNavigation fregularNav
const G4String & GetName() const
virtual EInside Inside(const G4ThreeVector &p) const =0
G4double ComputeStep(const G4ThreeVector &globalPoint, const G4ThreeVector &globalDirection, const G4double currentProposedStepLength, G4double &newSafety, G4NavigationHistory &history, G4bool &validExitNormal, G4ThreeVector &exitNormal, G4bool &exiting, G4bool &entering, G4VPhysicalVolume *(*pBlockedPhysical), G4int &blockedReplicaNo)
EVolume GetTopVolumeType() const
virtual G4int GetRegularStructureId() const =0
virtual G4double ComputeSafety(const G4ThreeVector &globalpoint, const G4NavigationHistory &history, const G4double pMaxLength=DBL_MAX)
virtual G4ThreeVector SurfaceNormal(const G4ThreeVector &p) const =0
virtual G4VPVParameterisation * GetParameterisation() const =0
virtual void SetCopyNo(G4int CopyNo)=0
const G4AffineTransform & GetGlobalToLocalTransform() const
G4ThreeVector fStepEndPoint
virtual G4ThreeVector GetLocalExitNormalAndCheck(const G4ThreeVector &point, G4bool *valid)
G4int GetReplicaNo(G4int n) const
G4NavigationHistory fHistory
G4bool LevelLocate(G4NavigationHistory &history, const G4VPhysicalVolume *blockedVol, const G4int blockedNum, const G4ThreeVector &globalPoint, const G4ThreeVector *globalDirection, const G4bool pLocatedOnEdge, G4ThreeVector &localPoint)
void ComputeStepLog(const G4ThreeVector &pGlobalpoint, G4double moveLenSq) const
virtual G4double DistanceToIn(const G4ThreeVector &p, const G4ThreeVector &v) const =0
G4ITNavigatorState_Lock1 * GetNavigatorState()
G4double CheckNextStep(const G4ThreeVector &pGlobalPoint, const G4ThreeVector &pDirection, const G4double pCurrentProposedStepLength, G4double &pNewSafety)
G4ThreeVector fPreviousSftOrigin
#define G4DEBUG_NAVIGATION
void SetNavigatorState(G4ITNavigatorState_Lock1 *)
void G4Exception(const char *originOfException, const char *exceptionCode, G4ExceptionSeverity severity, const char *comments)
G4bool sWasLimitedByGeometry
G4ThreeVector sExitNormal
G4ParameterisedNavigation fparamNav
static const double perMillion
void NewLevel(G4VPhysicalVolume *pNewMother, EVolume vType=kNormal, G4int nReplica=-1)
G4bool fLocatedOutsideWorld
G4double ComputeStep(const G4ThreeVector &localPoint, const G4ThreeVector &localDirection, const G4double currentProposedStepLength, G4double &newSafety, G4NavigationHistory &history, G4bool &validExitNormal, G4ThreeVector &exitNormal, G4bool &exiting, G4bool &entering, G4VPhysicalVolume *(*pBlockedPhysical), G4int &blockedReplicaNo)
G4LogicalVolume * GetLogicalVolume() const
const G4AffineTransform GetLocalToGlobalTransform() const
const G4AffineTransform & GetTransform(G4int n) const
G4ThreeVector fLastLocatedPointLocal
EVolume CharacteriseDaughters(const G4LogicalVolume *pLog) const
virtual void ComputeTransformation(const G4int, G4VPhysicalVolume *) const =0
G4ThreeVector sPreviousSftOrigin
virtual G4int GetCopyNo() const =0
G4bool fCalculatedExitNormal
const G4RotationMatrix * GetRotation() const
G4bool EnteredDaughterVolume() const
G4SmartVoxelNode * VoxelLocate(G4SmartVoxelHeader *pHead, const G4ThreeVector &localPoint)
G4double ComputeSafety(const G4ThreeVector &localPoint, const G4NavigationHistory &history, const G4double pProposedMaxLength=DBL_MAX)
G4ThreeVector sLastLocatedPointLocal
const G4AffineTransform & GetTopTransform() const
const G4NavigationHistory * GetHistory() const
G4VPhysicalVolume * spBlockedPhysicalVolume
G4VPhysicalVolume * fBlockedPhysicalVolume
G4int MoveUpHistory(G4int num_levels=1)
const G4String & GetName() const
G4double ComputeSafety(const G4ThreeVector &globalPoint, const G4ThreeVector &localPoint, G4NavigationHistory &history, const G4double pProposedMaxLength=DBL_MAX)
G4bool LevelLocate(G4NavigationHistory &history, const G4VPhysicalVolume *blockedVol, const G4int blockedNum, const G4ThreeVector &globalPoint, const G4ThreeVector *globalDirection, const G4bool pLocatedOnEdge, G4ThreeVector &localPoint)
G4VPhysicalVolume * GetVolume(G4int n) const
virtual G4ThreeVector GetGlobalExitNormal(const G4ThreeVector &point, G4bool *valid)
G4ReplicaNavigation freplicaNav
G4bool LevelLocate(G4NavigationHistory &history, const G4VPhysicalVolume *blockedVol, const G4int blockedNum, const G4ThreeVector &globalPoint, const G4ThreeVector *globalDirection, const G4bool pLocatedOnEdge, G4ThreeVector &localPoint)
virtual G4double DistanceToOut(const G4ThreeVector &p, const G4ThreeVector &v, const G4bool calcNorm=false, G4bool *validNorm=0, G4ThreeVector *n=0) const =0
virtual G4double ComputeStep(const G4ThreeVector &globalPoint, const G4ThreeVector &globalDirection, const G4double currentProposedStepLength, G4double &newSafety, G4NavigationHistory &history, G4bool &validExitNormal, G4ThreeVector &exitNormal, G4bool &exiting, G4bool &entering, G4VPhysicalVolume *(*pBlockedPhysical), G4int &blockedReplicaNo)
virtual G4double ComputeSafety(const G4ThreeVector &globalpoint, const G4double pProposedMaxLength=DBL_MAX, const G4bool keepState=true)
EInside BackLocate(G4NavigationHistory &history, const G4ThreeVector &globalPoint, G4ThreeVector &localPoint, const G4bool &exiting, G4bool ¬KnownInside) const
G4AffineTransform GetMotherToDaughterTransform(G4VPhysicalVolume *dVolume, G4int dReplicaNo, EVolume dVolumeType)
virtual G4bool CheckOverlaps(G4int res=1000, G4double tol=0., G4bool verbose=true, G4int errMax=1)
G4ThreeVector fGrandMotherExitNormal
G4ThreeVector fExitNormalGlobalFrame
static G4GeometryTolerance * GetInstance()
G4VSolid * GetSolid() const
G4GLOB_DLL std::ostream G4cerr
G4bool fWasLimitedByGeometry
virtual G4TouchableHistoryHandle CreateTouchableHistoryHandle() const
virtual ~G4ITNavigator1()