ECCE @ EIC Software
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Macros Groups Pages
G4ITSafetyHelper.hh
Go to the documentation of this file. Or view the newest version in sPHENIX GitHub for file G4ITSafetyHelper.hh
1 //
2 // ********************************************************************
3 // * License and Disclaimer *
4 // * *
5 // * The Geant4 software is copyright of the Copyright Holders of *
6 // * the Geant4 Collaboration. It is provided under the terms and *
7 // * conditions of the Geant4 Software License, included in the file *
8 // * LICENSE and available at http://cern.ch/geant4/license . These *
9 // * include a list of copyright holders. *
10 // * *
11 // * Neither the authors of this software system, nor their employing *
12 // * institutes,nor the agencies providing financial support for this *
13 // * work make any representation or warranty, express or implied, *
14 // * regarding this software system or assume any liability for its *
15 // * use. Please see the license in the file LICENSE and URL above *
16 // * for the full disclaimer and the limitation of liability. *
17 // * *
18 // * This code implementation is the result of the scientific and *
19 // * technical work of the GEANT4 collaboration. *
20 // * By using, copying, modifying or distributing the software (or *
21 // * any work based on the software) you agree to acknowledge its *
22 // * use in resulting scientific publications, and indicate your *
23 // * acceptance of all terms of the Geant4 Software license. *
24 // ********************************************************************
25 //
26 //
27 //
28 //
29 // class G4ITSafetyHelper
30 //
31 // Class description:
32 //
33 // G4ITSafetyHelper is a duplicated version of G4SafetyHelper
34 // This class is a helper for physics processes which require
35 // knowledge of the safety, and the step size for the 'mass' geometry
36 
37 // First version: J. Apostolakis, July 5th, 2006
38 // Modified:
39 // 10.04.07 V.Ivanchenko Use unique G4ITSafetyHelper
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 G4TrackStateDependent<G4ITSafetyHelper>
55 {
56 public:
57  // with description
60  //
61  // Constructor and destructor
62 
64  const G4ThreeVector& direction,
65  const G4double currentMaxStep,
66  G4double& newSafety);
67  //
68  // Return linear step for mass geometry
69 
70  G4double ComputeSafety(const G4ThreeVector& pGlobalPoint,
71  G4double maxRadius = DBL_MAX); // Radius of interest
72  //
73  // Return safety for all geometries.
74  //
75  // The 2nd argument is the radius of your interest (e.g. maximum displacement )
76  // Giving this you can reduce the average computational cost.
77  // If the second argument is not given, this is the real isotropic safety
78 
79  void Locate(const G4ThreeVector& pGlobalPoint,
80  const G4ThreeVector& direction);
81  //
82  // Locate the point for all geometries
83 
84  void ReLocateWithinVolume(const G4ThreeVector& pGlobalPoint );
85  //
86  // Relocate the point in the volume of interest
87 
88  inline void EnableParallelNavigation(G4bool parallel);
89  //
90  // To have parallel worlds considered, must be true.
91  // Alternative is to use single (mass) Navigator directly
92 
93  void InitialiseNavigator();
94  //
95  // Check for new navigator for tracking, and reinitialise pointer
96 
98  { G4int oldlv= fVerbose; fVerbose= lev; return oldlv;}
99 
101  inline void SetCurrentSafety(G4double val, const G4ThreeVector& pos);
102 
103 public: // without description
104 
105  void InitialiseHelper();
106 
107 private:
108 
110  G4ITNavigator* fpMassNavigator;
112 
114  // Flag whether to use PathFinder or single (mass) Navigator directly
116  // Flag of first call
118  // Whether to print warning in case of move outside safety
119 
120 public:
121  // State used during tracking -- for optimisation
122  class State
123  {
124  friend class G4ITSafetyHelper;
127 
128  public:
129  State() :
130  fLastSafetyPosition(0.0,0.0,0.0),
131  fLastSafety(0.0)
132  {}
133 
134  virtual ~State()
135  {}
136  };
137 
138  // const G4double fRecomputeFactor;
139  // parameter for further optimisation:
140  // if ( move < fact*safety ) do fast recomputation of safety
141  // End State (tracking)
142 };
143 
145 
146 // Inline definitions
147 
148 inline
150 {
151  fUseParallelGeometries = parallel;
152 }
153 
155 {
156  return fpMassNavigator->GetWorldVolume();
157 }
158 
159 inline
161 {
162  fpTrackState->fLastSafety = val;
163  fpTrackState->fLastSafetyPosition = pos;
164 }
165 
166 #endif