ECCE @ EIC Software
Reference for
ECCE @ EIC
simulation and reconstruction software on GitHub
Home page
Related Pages
Modules
Namespaces
Classes
Files
External Links
File List
File Members
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
46
class
G4SafetyHelper
47
{
48
public
:
// with description
49
50
G4SafetyHelper
();
51
~G4SafetyHelper
();
52
// Constructor and destructor
53
54
G4double
CheckNextStep
(
const
G4ThreeVector
&
position
,
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,
61
G4double
maxRadius
=
DBL_MAX
);
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
76
G4bool
RecheckDistanceToCurrentBoundary
(
const
G4ThreeVector
& pGlobalPoint,
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 );
98
inline
G4VPhysicalVolume
*
GetWorldVolume
();
99
inline
void
SetCurrentSafety
(
G4double
val,
const
G4ThreeVector
&
pos
);
100
101
public
:
// without description
102
103
void
InitialiseHelper
();
104
105
private
:
106
107
G4PathFinder
*
fpPathFinder
=
nullptr
;
108
G4Navigator
*
fpMassNavigator
=
nullptr
;
109
G4int
fMassNavigatorId
= -1;
110
111
G4bool
fUseParallelGeometries
=
false
;
112
// Flag whether to use PathFinder or single (mass) Navigator directly
113
// By default, one geometry only
114
G4bool
fFirstCall
=
true
;
115
// Flag of first call
116
G4int
fVerbose
= 0;
117
// Whether to print warning in case of move outside safety
118
119
// State used during tracking -- for optimisation
120
121
G4ThreeVector
fLastSafetyPosition
;
122
G4double
fLastSafety
= 0.0;
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
135
inline
G4int
G4SafetyHelper::SetVerboseLevel
(
G4int
lev )
136
{
137
G4int
oldlv =
fVerbose
;
138
fVerbose
= lev;
139
return
oldlv;
140
}
141
142
inline
143
void
G4SafetyHelper::EnableParallelNavigation
(
G4bool
parallel)
144
{
145
fUseParallelGeometries
= parallel;
146
}
147
148
inline
149
G4VPhysicalVolume
*
G4SafetyHelper::GetWorldVolume
()
150
{
151
return
fpMassNavigator
->
GetWorldVolume
();
152
}
153
154
inline
155
void
G4SafetyHelper::SetCurrentSafety
(
G4double
val,
const
G4ThreeVector
&
pos
)
156
{
157
fLastSafety
= val;
158
fLastSafetyPosition
=
pos
;
159
}
160
161
#endif
geant4
tree
geant4-10.6-release
source
geometry
navigation
include
G4SafetyHelper.hh
Built by
Jin Huang
. updated:
Wed Jun 29 2022 17:25:17
using
1.8.2 with
ECCE GitHub integration