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
G4VContinuousDiscreteProcess.cc
Go to the documentation of this file.
Or view
the newest version in sPHENIX GitHub for file G4VContinuousDiscreteProcess.cc
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
// --------------------------------------------------------------
30
// GEANT 4 class implementation file
31
//
32
// History: first implementation, based on object model of
33
// 2nd December 1995, G.Cosmo
34
// --------------------------------------------------------------
35
// New Physics scheme 8 Jan. 1997 H.Kurahige
36
// ------------------------------------------------------------
37
38
#include "
G4VContinuousDiscreteProcess.hh
"
39
#include "
G4SystemOfUnits.hh
"
40
41
#include "
G4Step.hh
"
42
#include "
G4Track.hh
"
43
#include "
G4MaterialTable.hh
"
44
#include "
G4VParticleChange.hh
"
45
46
G4VContinuousDiscreteProcess::G4VContinuousDiscreteProcess
()
47
:
G4VProcess
(
"No Name Discrete Process"
),
48
valueGPILSelection(
CandidateForSelection
)
49
{
50
G4Exception
(
"G4VContinuousDiscreteProcess::G4VContinuousDiscreteProcess()"
,
"ProcMan102"
,
51
JustWarning
,
"Default constructor is called"
);
52
}
53
54
G4VContinuousDiscreteProcess::G4VContinuousDiscreteProcess
(
const
G4String
& aName ,
G4ProcessType
aType)
55
:
G4VProcess
(aName, aType),
56
valueGPILSelection(
CandidateForSelection
)
57
{
58
enableAtRestDoIt
=
false
;
59
}
60
61
G4VContinuousDiscreteProcess::~G4VContinuousDiscreteProcess
()
62
{
63
}
64
65
G4VContinuousDiscreteProcess::G4VContinuousDiscreteProcess
(
G4VContinuousDiscreteProcess
&
right
)
66
:
G4VProcess
(right),
67
valueGPILSelection(right.valueGPILSelection)
68
{
69
}
70
71
G4double
G4VContinuousDiscreteProcess::PostStepGetPhysicalInteractionLength
(
72
const
G4Track
&
track
,
73
G4double
previousStepSize,
74
G4ForceCondition
*
condition
75
)
76
{
77
if
( (previousStepSize <=0.0) || (
theNumberOfInteractionLengthLeft
<=0.0)) {
78
// beggining of tracking (or just after DoIt of this process)
79
ResetNumberOfInteractionLengthLeft
();
80
}
else
if
( previousStepSize > 0.0) {
81
// subtract NumberOfInteractionLengthLeft
82
SubtractNumberOfInteractionLengthLeft
(previousStepSize);
83
}
else
{
84
// zero step
85
// DO NOTHING
86
}
87
88
// condition is set to "Not Forced"
89
*condition =
NotForced
;
90
91
// get mean free path
92
currentInteractionLength
=
GetMeanFreePath
(track, previousStepSize, condition);
93
94
G4double
value
;
95
if
(
currentInteractionLength
<
DBL_MAX
) {
96
value =
theNumberOfInteractionLengthLeft
*
currentInteractionLength
;
97
}
else
{
98
value =
DBL_MAX
;
99
}
100
#ifdef G4VERBOSE
101
if
(
verboseLevel
>1){
102
G4cout
<<
"G4VContinuousDiscreteProcess::PostStepGetPhysicalInteractionLength "
;
103
G4cout
<<
"[ "
<<
GetProcessName
() <<
"]"
<<
G4endl
;
104
track.
GetDynamicParticle
()->
DumpInfo
();
105
G4cout
<<
" in Material "
<< track.
GetMaterial
()->
GetName
() <<
G4endl
;
106
G4cout
<<
"InteractionLength= "
<< value/
cm
<<
"[cm] "
<<
G4endl
;
107
}
108
#endif
109
return
value
;
110
}
111
112
113
G4VParticleChange
*
G4VContinuousDiscreteProcess::PostStepDoIt
(
114
const
G4Track
& ,
115
const
G4Step
&
116
)
117
{
118
// clear NumberOfInteractionLengthLeft
119
ClearNumberOfInteractionLengthLeft
();
120
return
pParticleChange
;
121
}
122
123
G4VParticleChange
*
G4VContinuousDiscreteProcess::AlongStepDoIt
(
124
const
G4Track
& ,
125
const
G4Step
&
126
)
127
{
128
// clear NumberOfInteractionLengthLeft
129
ClearNumberOfInteractionLengthLeft
();
130
return
pParticleChange
;
131
}
132
133
G4double
G4VContinuousDiscreteProcess::AlongStepGetPhysicalInteractionLength
(
134
const
G4Track
&
track
,
135
G4double
previousStepSize,
136
G4double
currentMinimumStep,
137
G4double
& currentSafety,
138
G4GPILSelection
* selection
139
)
140
{
141
// GPILSelection is set to defaule value of CandidateForSelection
142
valueGPILSelection
=
CandidateForSelection
;
143
144
// get Step limit proposed by the process
145
G4double
steplength =
GetContinuousStepLimit
(track,previousStepSize,currentMinimumStep, currentSafety);
146
147
// set return value for G4GPILSelection
148
*selection =
valueGPILSelection
;
149
150
#ifdef G4VERBOSE
151
if
(
verboseLevel
>1){
152
G4cout
<<
"G4VContinuousDiscreteProcess::AlongStepGetPhysicalInteractionLength "
;
153
G4cout
<<
"[ "
<<
GetProcessName
() <<
"]"
<<
G4endl
;
154
track.
GetDynamicParticle
()->
DumpInfo
();
155
G4cout
<<
" in Material "
<< track.
GetMaterial
()->
GetName
() <<
G4endl
;
156
G4cout
<<
"IntractionLength= "
<< steplength/
cm
<<
"[cm] "
<<
G4endl
;
157
}
158
#endif
159
return
steplength ;
160
}
geant4
tree
geant4-10.6-release
source
processes
management
src
G4VContinuousDiscreteProcess.cc
Built by
Jin Huang
. updated:
Wed Jun 29 2022 17:25:52
using
1.8.2 with
ECCE GitHub integration