ECCE @ EIC Software
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Macros Groups Pages
G4SafetyHelper.hh
Go to the documentation of this file. Or view the newest version in sPHENIX GitHub for file G4SafetyHelper.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 // class G4SafetyHelper
27 //
28 // Class description:
29 //
30 // This class is a helper for physics processes which require
31 // knowledge of the safety, and the step size for the 'mass' geometry
32 
33 // First version: J.Apostolakis, July 5th, 2006
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 
47 {
48  public: // with description
49 
50  G4SafetyHelper();
52  // Constructor and destructor
53 
55  const G4ThreeVector& direction,
56  const G4double currentMaxStep,
57  G4double& newSafety );
58  // Return linear step for mass geometry
59 
60  G4double ComputeSafety( const G4ThreeVector& pGlobalPoint,
62  // Return safety for all geometries.
63  //
64  // The 2nd argument is the radius of your interest (e.g. maximum
65  // displacement). Giving this you can reduce the average computational
66  // cost. If the second argument is not given, this is the real
67  // isotropic safety
68 
69  void Locate(const G4ThreeVector& pGlobalPoint,
70  const G4ThreeVector& direction);
71  // Locate the point for all geometries
72 
73  void ReLocateWithinVolume(const G4ThreeVector& pGlobalPoint );
74  // Relocate the point in the volume of interest
75 
77  const G4ThreeVector& pDirection,
78  const G4double pCurrentProposedStepLength,
79  G4double* prDistance,
80  G4double* prNewSafety = nullptr) const;
81  // Trial method for checking potential displacement for MS
82  // Check new Globalpoint, to see whether it is in current volume
83  // (mother) and not in potential entering daughter.
84  // If in mother, check distance to boundary along pDirection.
85  // If in entering daughter, check distance back to boundary.
86  // NOTE:
87  // Can be called only after ComputeStep is called - before ReLocation
88  // Deals only with current volume (and potentially entered)
89 
90  inline void EnableParallelNavigation(G4bool parallel);
91  // To have parallel worlds considered, must be true.
92  // Alternative is to use single (mass) Navigator directly
93 
94  void InitialiseNavigator();
95  // Check for new navigator for tracking, and reinitialise pointer
96 
97  inline G4int SetVerboseLevel( G4int lev );
99  inline void SetCurrentSafety(G4double val, const G4ThreeVector& pos);
100 
101  public: // without description
102 
103  void InitialiseHelper();
104 
105  private:
106 
110 
112  // Flag whether to use PathFinder or single (mass) Navigator directly
113  // By default, one geometry only
115  // Flag of first call
117  // Whether to print warning in case of move outside safety
118 
119  // State used during tracking -- for optimisation
120 
123 
124  // const G4double fRecomputeFactor = 0.0;
125  // parameter for further optimisation:
126  // if ( move < fact*safety ) do fast recomputation of safety
127 
128  // End State (tracking)
129 };
130 
131 // --------------------------------------------------------------------
132 // Inline definitions
133 // --------------------------------------------------------------------
134 
136 {
137  G4int oldlv = fVerbose;
138  fVerbose = lev;
139  return oldlv;
140 }
141 
142 inline
144 {
145  fUseParallelGeometries = parallel;
146 }
147 
148 inline
150 {
152 }
153 
154 inline
156 {
157  fLastSafety = val;
159 }
160 
161 #endif