| 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 //
27 //
28 //
29 // class G4ITSafetyHelper
30 //
31 // Class description:
32 //
33 // G4ITSafetyHelper is a duplicated version of
34 // This class is a helper for physics processe
35 // knowledge of the safety, and the step size
36
37 // First version: J. Apostolakis, July 5th,
38 // Modified:
39 // 10.04.07 V.Ivanchenko Use unique G4ITSafe
40 // -------------------------------------------
41
42 #ifndef G4ITSAFETYHELPER_HH
43 #define G4ITSAFETYHELPER_HH 1
44
45 #include <vector>
46
47 #include "G4Types.hh"
48 #include "G4ThreeVector.hh"
49 #include "G4ITNavigator.hh"
50 #include "G4TrackState.hh"
51
52 class G4PathFinder;
53
54 class G4ITSafetyHelper : public G4TrackStateDe
55 {
56 public:
57 // with description
58 G4ITSafetyHelper();
59 ~G4ITSafetyHelper() override;
60 //
61 // Constructor and destructor
62
63 G4double CheckNextStep(const G4ThreeVector&
64 const G4ThreeVector&
65 const G4double curren
66 G4double& newSafety);
67 //
68 // Return linear step for mass geometry
69
70 G4double ComputeSafety(const G4ThreeVector&
71 G4double maxRadius =
72 //
73 // Return safety for all geometries.
74 //
75 // The 2nd argument is the radius of your i
76 // Giving this you can reduce the average
77 // If the second argument is not given, thi
78
79 void Locate(const G4ThreeVector& pGlobalPoin
80 const G4ThreeVector& direction);
81 //
82 // Locate the point for all geometries
83
84 void ReLocateWithinVolume(const G4ThreeVecto
85 //
86 // Relocate the point in the volume of inter
87
88 inline void EnableParallelNavigation(G4bool
89 //
90 // To have parallel worlds considered, must
91 // Alternative is to use single (mass) Navi
92
93 void InitialiseNavigator();
94 //
95 // Check for new navigator for tracking, and
96
97 G4int SetVerboseLevel( G4int lev )
98 { G4int oldlv= fVerbose; fVerbose= lev; retu
99
100 inline G4VPhysicalVolume* GetWorldVolume();
101 inline void SetCurrentSafety(G4double val, c
102
103 public: // without description
104
105 void InitialiseHelper();
106
107 private:
108
109 G4PathFinder* fpPathFinder;
110 G4ITNavigator* fpMassNavigator;
111 G4int fMassNavigatorId;
112
113 G4bool fUseParallelGeometries{false};
114 // Flag whether to use PathFinder or single
115 G4bool fFirstCall{true};
116 // Flag of first call
117 G4int fVerbose{0};
118 // Whether to print warning in case of move
119
120 public:
121 // State used during tracking -- for optimis
122 class State
123 {
124 friend class G4ITSafetyHelper;
125 G4ThreeVector fLastSafetyPosition;
126 G4double fLastSafety{0.0};
127
128 public:
129 State() :
130 fLastSafetyPosition(0.0,0.0,0.0)
131 {}
132
133 virtual ~State()
134 = default;
135 };
136
137 // const G4double fRecomputeFactor;
138 // parameter for further optimisation:
139 // if ( move < fact*safety ) do fast recomp
140 // End State (tracking)
141 };
142
143 RegisterTrackState(G4ITSafetyHelper, State)
144
145 // Inline definitions
146
147 inline
148 void G4ITSafetyHelper::EnableParallelNavigatio
149 {
150 fUseParallelGeometries = parallel;
151 }
152
153 inline G4VPhysicalVolume* G4ITSafetyHelper::Ge
154 {
155 return fpMassNavigator->GetWorldVolume();
156 }
157
158 inline
159 void G4ITSafetyHelper::SetCurrentSafety(G4doub
160 {
161 fpTrackState->fLastSafety = val;
162 fpTrackState->fLastSafetyPosition = pos;
163 }
164
165 #endif
166