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
PhysListEmStandard.cc
Go to the documentation of this file.
Or view
the newest version in sPHENIX GitHub for file PhysListEmStandard.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
//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
28
//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
29
30
#include "PhysListEmStandard.hh"
31
#include "
G4ParticleDefinition.hh
"
32
#include "
G4ProcessManager.hh
"
33
#include "
G4PhysicsListHelper.hh
"
34
35
#include "
G4ComptonScattering.hh
"
36
#include "
G4GammaConversion.hh
"
37
#include "
G4PhotoElectricEffect.hh
"
38
#include "
G4RayleighScattering.hh
"
39
#include "
G4KleinNishinaModel.hh
"
40
#include "
G4LivermorePhotoElectricModel.hh
"
41
42
#include "
G4eMultipleScattering.hh
"
43
#include "
G4eIonisation.hh
"
44
#include "
G4eBremsstrahlung.hh
"
45
#include "
G4eplusAnnihilation.hh
"
46
47
#include "
G4MuMultipleScattering.hh
"
48
#include "
G4MuIonisation.hh
"
49
#include "
G4MuBremsstrahlung.hh
"
50
#include "
G4MuPairProduction.hh
"
51
52
#include "
G4hMultipleScattering.hh
"
53
#include "
G4hIonisation.hh
"
54
#include "
G4hBremsstrahlung.hh
"
55
#include "
G4hPairProduction.hh
"
56
57
#include "
G4ionIonisation.hh
"
58
#include "
G4IonParametrisedLossModel.hh
"
59
#include "
G4NuclearStopping.hh
"
60
61
#include "
G4MscStepLimitType.hh
"
62
63
#include "
G4LossTableManager.hh
"
64
#include "
G4UAtomicDeexcitation.hh
"
65
66
#include "
G4BuilderType.hh
"
67
#include "
G4SystemOfUnits.hh
"
68
69
//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
70
71
PhysListEmStandard::PhysListEmStandard
(
const
G4String
&
name
)
72
:
G4VPhysicsConstructor
(name)
73
{
74
G4EmParameters
* param =
G4EmParameters::Instance
();
75
param->
SetDefaults
();
76
param->
SetMinEnergy
(10*
eV
);
77
param->
SetMaxEnergy
(10*
TeV
);
78
param->
SetNumberOfBinsPerDecade
(10);
79
param->
SetMscStepLimitType
(
fUseSafetyPlus
);
80
SetPhysicsType
(
bElectromagnetic
);
81
}
82
83
//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
84
85
PhysListEmStandard::~PhysListEmStandard
()
86
{}
87
88
//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
89
90
void
PhysListEmStandard::ConstructProcess
()
91
{
92
G4PhysicsListHelper
* ph =
G4PhysicsListHelper::GetPhysicsListHelper
();
93
94
// Add standard EM Processes
95
//
96
auto
particleIterator
=
GetParticleIterator
();
97
particleIterator
->reset();
98
while
( (*
particleIterator
)() ){
99
G4ParticleDefinition
*
particle
=
particleIterator
->value();
100
G4String
particleName = particle->
GetParticleName
();
101
102
if
(particleName ==
"gamma"
) {
103
105
ph->
RegisterProcess
(
new
G4PhotoElectricEffect
, particle);
109
110
G4ComptonScattering
* cs =
new
G4ComptonScattering
;
111
cs->
SetEmModel
(
new
G4KleinNishinaModel
());
112
ph->
RegisterProcess
(cs, particle);
113
ph->
RegisterProcess
(
new
G4GammaConversion
, particle);
114
115
}
else
if
(particleName ==
"e-"
) {
116
117
ph->
RegisterProcess
(
new
G4eMultipleScattering
(), particle);
118
//
119
G4eIonisation
* eIoni =
new
G4eIonisation
();
120
eIoni->
SetStepFunction
(0.1, 100*
um
);
121
ph->
RegisterProcess
(eIoni, particle);
122
//
123
ph->
RegisterProcess
(
new
G4eBremsstrahlung
(), particle);
124
125
}
else
if
(particleName ==
"e+"
) {
126
127
ph->
RegisterProcess
(
new
G4eMultipleScattering
(), particle);
128
//
129
G4eIonisation
* eIoni =
new
G4eIonisation
();
130
eIoni->
SetStepFunction
(0.1, 100*
um
);
131
ph->
RegisterProcess
(eIoni, particle);
132
//
133
ph->
RegisterProcess
(
new
G4eBremsstrahlung
(), particle);
134
//
135
ph->
RegisterProcess
(
new
G4eplusAnnihilation
(), particle);
136
137
}
else
if
(particleName ==
"mu+"
||
138
particleName ==
"mu-"
) {
139
140
ph->
RegisterProcess
(
new
G4MuMultipleScattering
(), particle);
141
G4MuIonisation
* muIoni =
new
G4MuIonisation
();
142
muIoni->
SetStepFunction
(0.1, 50*
um
);
143
ph->
RegisterProcess
(muIoni, particle);
144
ph->
RegisterProcess
(
new
G4MuBremsstrahlung
(), particle);
145
ph->
RegisterProcess
(
new
G4MuPairProduction
(), particle);
146
147
}
else
if
( particleName ==
"proton"
||
148
particleName ==
"pi-"
||
149
particleName ==
"pi+"
) {
150
151
ph->
RegisterProcess
(
new
G4hMultipleScattering
(), particle);
152
G4hIonisation
* hIoni =
new
G4hIonisation
();
153
hIoni->
SetStepFunction
(0.1, 20*
um
);
154
ph->
RegisterProcess
(hIoni, particle);
155
ph->
RegisterProcess
(
new
G4hBremsstrahlung
(), particle);
156
ph->
RegisterProcess
(
new
G4hPairProduction
(), particle);
157
158
}
else
if
( particleName ==
"alpha"
||
159
particleName ==
"He3"
) {
160
161
ph->
RegisterProcess
(
new
G4hMultipleScattering
(), particle);
162
G4ionIonisation
* ionIoni =
new
G4ionIonisation
();
163
ionIoni->
SetStepFunction
(0.1, 1*
um
);
164
ph->
RegisterProcess
(ionIoni, particle);
165
ph->
RegisterProcess
(
new
G4NuclearStopping
(), particle);
166
167
}
else
if
( particleName ==
"GenericIon"
) {
168
169
ph->
RegisterProcess
(
new
G4hMultipleScattering
(), particle);
170
G4ionIonisation
* ionIoni =
new
G4ionIonisation
();
171
ionIoni->
SetEmModel
(
new
G4IonParametrisedLossModel
());
172
ionIoni->
SetStepFunction
(0.1, 1*
um
);
173
ph->
RegisterProcess
(ionIoni, particle);
174
ph->
RegisterProcess
(
new
G4NuclearStopping
(), particle);
175
176
}
else
if
((!particle->
IsShortLived
()) &&
177
(particle->
GetPDGCharge
() != 0.0) &&
178
(particle->
GetParticleName
() !=
"chargedgeantino"
)) {
179
180
//all others charged particles except geantino
181
ph->
RegisterProcess
(
new
G4hMultipleScattering
(), particle);
182
ph->
RegisterProcess
(
new
G4hIonisation
(), particle);
183
}
184
}
185
186
// Deexcitation
187
//
188
G4VAtomDeexcitation
* de =
new
G4UAtomicDeexcitation
();
189
de->
SetFluo
(
true
);
190
de->
SetAuger
(
false
);
191
de->
SetPIXE
(
false
);
192
G4LossTableManager::Instance
()->
SetAtomDeexcitation
(de);
193
}
194
195
//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
196
geant4
tree
geant4-10.6-release
examples
extended
electromagnetic
TestEm1
src
PhysListEmStandard.cc
Built by
Jin Huang
. updated:
Wed Jun 29 2022 17:24:51
using
1.8.2 with
ECCE GitHub integration