Geant4 Cross Reference |
1 // 1 2 // ******************************************* 3 // * License and Disclaimer 4 // * 5 // * The Geant4 software is copyright of th 6 // * the Geant4 Collaboration. It is provided 7 // * conditions of the Geant4 Software License 8 // * LICENSE and available at http://cern.ch/ 9 // * include a list of copyright holders. 10 // * 11 // * Neither the authors of this software syst 12 // * institutes,nor the agencies providing fin 13 // * work make any representation or warran 14 // * regarding this software system or assum 15 // * use. Please see the license in the file 16 // * for the full disclaimer and the limitatio 17 // * 18 // * This code implementation is the result 19 // * technical work of the GEANT4 collaboratio 20 // * By using, copying, modifying or distri 21 // * any work based on the software) you ag 22 // * use in resulting scientific publicati 23 // * acceptance of all terms of the Geant4 Sof 24 // ******************************************* 25 // 26 // class G4SafetyHelper 27 // 28 // Class description: 29 // 30 // This class is a helper for physics processe 31 // knowledge of the safety, and the step size 32 33 // First version: J.Apostolakis, July 5th, 2 34 // ------------------------------------------- 35 #ifndef G4SAFETYHELPER_HH 36 #define G4SAFETYHELPER_HH 1 37 38 #include <vector> 39 40 #include "G4Types.hh" 41 #include "G4ThreeVector.hh" 42 #include "G4Navigator.hh" 43 44 class G4PathFinder; 45 46 class G4SafetyHelper 47 { 48 public: // with description 49 50 G4SafetyHelper(); 51 ~G4SafetyHelper(); 52 // Constructor and destructor 53 54 G4double CheckNextStep( const G4ThreeVecto 55 const G4ThreeVecto 56 const G4double cur 57 G4double& ne 58 // Return linear step for mass geometry 59 60 G4double ComputeSafety( const G4ThreeVecto 61 G4double maxRadius 62 // Return safety for all geometries. 63 // 64 // The 2nd argument is the radius of you 65 // displacement). Giving this you can re 66 // cost. If the second argument is not g 67 // isotropic safety 68 69 void Locate(const G4ThreeVector& pGlobalPo 70 const G4ThreeVector& direction 71 // Locate the point for all geometries 72 73 void ReLocateWithinVolume(const G4ThreeVec 74 // Relocate the point in the volume of i 75 76 inline void EnableParallelNavigation(G4boo 77 // To have parallel worlds considered, 78 // Alternative is to use single (mass) 79 80 void InitialiseNavigator(); 81 // Check for new navigator for tracking, 82 83 inline G4int SetVerboseLevel( G4int lev ); 84 inline G4VPhysicalVolume* GetWorldVolume() 85 inline void SetCurrentSafety(G4double val, 86 87 public: // without description 88 89 void InitialiseHelper(); 90 91 private: 92 93 G4PathFinder* fpPathFinder = nullptr; 94 G4Navigator* fpMassNavigator = nullptr; 95 96 G4bool fUseParallelGeometries = false; 97 // Flag whether to use PathFinder or sin 98 // By default, one geometry only 99 G4bool fFirstCall = true; 100 // Flag of first call 101 G4int fVerbose = 0; 102 // Whether to print warning in case of m 103 104 // State used during tracking -- for optim 105 106 G4ThreeVector fLastSafetyPosition; 107 G4double fLastSafety = 0.0; 108 109 // const G4double fRecomputeFactor = 0.0; 110 // parameter for further optimisation: 111 // if ( move < fact*safety ) do fast r 112 113 // End State (tracking) 114 }; 115 116 // ------------------------------------------- 117 // Inline definitions 118 // ------------------------------------------- 119 120 inline G4int G4SafetyHelper::SetVerboseLevel( 121 { 122 G4int oldlv = fVerbose; 123 fVerbose = lev; 124 return oldlv; 125 } 126 127 inline 128 void G4SafetyHelper::EnableParallelNavigation( 129 { 130 fUseParallelGeometries = parallel; 131 } 132 133 inline 134 G4VPhysicalVolume* G4SafetyHelper::GetWorldVol 135 { 136 return fpMassNavigator->GetWorldVolume(); 137 } 138 139 inline 140 void G4SafetyHelper::SetCurrentSafety(G4double 141 { 142 fLastSafety = val; 143 fLastSafetyPosition = pos; 144 } 145 146 #endif 147