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
G4LowECapture.cc
Go to the documentation of this file.
Or view
the newest version in sPHENIX GitHub for file G4LowECapture.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
// ClassName: G4LowECapture
30
//
31
// Author: V.Ivanchenko 31 August 2010
32
//
33
//----------------------------------------------------------------------------
34
//
35
//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
36
//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
37
38
#include "
G4LowECapture.hh
"
39
#include "
G4SystemOfUnits.hh
"
40
#include "
G4ParticleDefinition.hh
"
41
#include "
G4Step.hh
"
42
#include "
G4Track.hh
"
43
#include "
G4Region.hh
"
44
#include "
G4RegionStore.hh
"
45
46
//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
47
48
G4LowECapture::G4LowECapture
(
G4double
ekinlim)
49
:
G4VDiscreteProcess
(
"Capture"
,
fElectromagnetic
),
50
kinEnergyThreshold(ekinlim), isIon(
false
),
51
nRegions(0)
52
{}
53
54
//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
55
56
G4LowECapture::~G4LowECapture
()
57
{}
58
59
//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
60
61
void
G4LowECapture::SetKinEnergyLimit
(
G4double
val)
62
{
63
kinEnergyThreshold
= val;
64
if
(
verboseLevel
> 0) {
65
G4cout
<<
"### G4LowECapture: Tracking cut E(MeV) = "
66
<<
kinEnergyThreshold
/
MeV
<<
G4endl
;
67
}
68
}
69
70
//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
71
72
void
G4LowECapture::AddRegion
(
const
G4String
& nam)
73
{
74
G4String
r
= nam;
75
if
(r ==
""
|| r ==
"world"
|| r ==
"World"
) r =
"DefaultRegionForTheWorld"
;
76
for
(
G4int
i=0; i<
nRegions
; ++i) {
77
if
(
regionName
[i] == r) {
return
; }
78
}
79
regionName
.push_back(r);
80
++
nRegions
;
81
if
(
verboseLevel
> 1) {
82
G4cout
<<
"### G4LowECapture: new G4Region <"
<< r <<
">"
<<
G4endl
;
83
}
84
}
85
86
//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
87
88
void
G4LowECapture::BuildPhysicsTable
(
const
G4ParticleDefinition
&
part
)
89
{
90
G4RegionStore
* store =
G4RegionStore::GetInstance
();
91
for
(
G4int
i=0; i<
nRegions
; ++i) {
92
const
G4Region
*
r
= store->
GetRegion
(
regionName
[i]);
93
if
(r &&
verboseLevel
> 0) {
94
G4cout
<<
"### G4LowECapture: new G4Region <"
95
<<
regionName
[i] <<
"> with tracking cut "
96
<<
kinEnergyThreshold
/
keV
<<
" keV"
<<
G4endl
;
97
}
98
if
(r) {
region
.push_back(r); }
99
}
100
nRegions =
region
.size();
101
102
// ions reusing G4GenericIon parameters
103
if
(part.
GetParticleType
() ==
"nucleus"
) {
104
G4String
pname
= part.
GetParticleName
();
105
if
(pname !=
"deuteron"
&& pname !=
"triton"
&&
106
pname !=
"alpha"
&& pname !=
"He3"
&&
107
pname !=
"alpha+"
&& pname !=
"helium"
&&
108
pname !=
"hydrogen"
) {
isIon
=
true
; }
109
}
110
}
111
112
//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
113
114
G4bool
G4LowECapture::IsApplicable
(
const
G4ParticleDefinition
&)
115
{
116
return
true
;
117
}
118
119
G4double
G4LowECapture::PostStepGetPhysicalInteractionLength
(
120
const
G4Track
& aTrack,
G4double
,
G4ForceCondition
*
condition
)
121
{
122
*condition =
NotForced
;
123
G4double
limit =
DBL_MAX
;
124
G4double
eLimit =
kinEnergyThreshold
;
125
if
(
isIon
) {
126
eLimit *= aTrack.
GetDefinition
()->
GetPDGMass
()/
CLHEP::proton_mass_c2
;
127
}
128
if
(aTrack.
GetKineticEnergy
() < eLimit) {
129
for
(
G4int
i=0; i<
nRegions
; ++i) {
130
if
(aTrack.
GetVolume
()->
GetLogicalVolume
()->
GetRegion
() ==
region
[i]) {
131
limit = 0.0;
132
break
;
133
}
134
}
135
}
136
return
limit;
137
}
138
139
//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
140
141
G4VParticleChange
*
G4LowECapture::PostStepDoIt
(
const
G4Track
& aTrack,
142
const
G4Step
&)
143
{
144
pParticleChange
->
Initialize
(aTrack);
145
pParticleChange
->
ProposeTrackStatus
(
fStopAndKill
);
146
pParticleChange
->
ProposeLocalEnergyDeposit
(aTrack.
GetKineticEnergy
());
147
return
pParticleChange
;
148
}
149
150
//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
151
152
G4double
G4LowECapture::GetMeanFreePath
(
const
G4Track
&,
G4double
,
153
G4ForceCondition
*)
154
{
155
return
DBL_MAX
;
156
}
157
158
//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
159
160
geant4
tree
geant4-10.6-release
source
processes
electromagnetic
lowenergy
src
G4LowECapture.cc
Built by
Jin Huang
. updated:
Wed Jun 29 2022 17:25:32
using
1.8.2 with
ECCE GitHub integration