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
41
#include "
G4eMultipleScattering.hh
"
42
#include "
G4eIonisation.hh
"
43
#include "
G4eBremsstrahlung.hh
"
44
#include "
G4eplusAnnihilation.hh
"
45
46
#include "
G4MuMultipleScattering.hh
"
47
#include "
G4MuIonisation.hh
"
48
#include "
G4MuBremsstrahlung.hh
"
49
#include "
G4MuPairProduction.hh
"
50
51
#include "
G4hMultipleScattering.hh
"
52
#include "
G4hIonisation.hh
"
53
#include "
G4hBremsstrahlung.hh
"
54
#include "
G4hPairProduction.hh
"
55
56
#include "
G4ionIonisation.hh
"
57
#include "
G4IonParametrisedLossModel.hh
"
58
#include "
G4NuclearStopping.hh
"
59
60
#include "
G4MscStepLimitType.hh
"
61
62
#include "
G4LossTableManager.hh
"
63
#include "
G4UAtomicDeexcitation.hh
"
64
65
#include "
G4BuilderType.hh
"
66
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
param->
SetFluo
(
true
);
81
SetPhysicsType
(
bElectromagnetic
);
82
}
83
84
//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
85
86
PhysListEmStandard::~PhysListEmStandard
()
87
{}
88
89
//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
90
91
void
PhysListEmStandard::ConstructProcess
()
92
{
93
G4PhysicsListHelper
* ph =
G4PhysicsListHelper::GetPhysicsListHelper
();
94
95
// Add standard EM Processes
96
//
97
auto
particleIterator
=
GetParticleIterator
();
98
particleIterator
->reset();
99
while
( (*
particleIterator
)() ){
100
G4ParticleDefinition
*
particle
=
particleIterator
->value();
101
G4String
particleName = particle->
GetParticleName
();
102
103
if
(particleName ==
"gamma"
) {
104
105
ph->
RegisterProcess
(
new
G4PhotoElectricEffect
, particle);
106
G4ComptonScattering
* cs =
new
G4ComptonScattering
;
107
cs->
SetEmModel
(
new
G4KleinNishinaModel
());
108
ph->
RegisterProcess
(cs, particle);
109
ph->
RegisterProcess
(
new
G4GammaConversion
, particle);
111
112
}
else
if
(particleName ==
"e-"
) {
113
114
ph->
RegisterProcess
(
new
G4eMultipleScattering
(), particle);
115
//
116
G4eIonisation
* eIoni =
new
G4eIonisation
();
117
eIoni->
SetStepFunction
(0.1, 100*
um
);
118
ph->
RegisterProcess
(eIoni, particle);
119
//
120
ph->
RegisterProcess
(
new
G4eBremsstrahlung
(), particle);
121
122
}
else
if
(particleName ==
"e+"
) {
123
124
ph->
RegisterProcess
(
new
G4eMultipleScattering
(), particle);
125
//
126
G4eIonisation
* eIoni =
new
G4eIonisation
();
127
eIoni->
SetStepFunction
(0.1, 100*
um
);
128
ph->
RegisterProcess
(eIoni, particle);
129
//
130
ph->
RegisterProcess
(
new
G4eBremsstrahlung
(), particle);
131
//
132
ph->
RegisterProcess
(
new
G4eplusAnnihilation
(), particle);
133
134
}
else
if
(particleName ==
"mu+"
||
135
particleName ==
"mu-"
) {
136
137
ph->
RegisterProcess
(
new
G4MuMultipleScattering
(), particle);
138
G4MuIonisation
* muIoni =
new
G4MuIonisation
();
139
muIoni->
SetStepFunction
(0.1, 50*
um
);
140
ph->
RegisterProcess
(muIoni, particle);
141
ph->
RegisterProcess
(
new
G4MuBremsstrahlung
(), particle);
142
ph->
RegisterProcess
(
new
G4MuPairProduction
(), particle);
143
144
}
else
if
( particleName ==
"proton"
||
145
particleName ==
"pi-"
||
146
particleName ==
"pi+"
) {
147
148
ph->
RegisterProcess
(
new
G4hMultipleScattering
(), particle);
149
G4hIonisation
* hIoni =
new
G4hIonisation
();
150
hIoni->
SetStepFunction
(0.1, 20*
um
);
151
ph->
RegisterProcess
(hIoni, particle);
152
ph->
RegisterProcess
(
new
G4hBremsstrahlung
(), particle);
153
ph->
RegisterProcess
(
new
G4hPairProduction
(), particle);
154
155
}
else
if
( particleName ==
"alpha"
||
156
particleName ==
"He3"
) {
157
158
ph->
RegisterProcess
(
new
G4hMultipleScattering
(), particle);
159
G4ionIonisation
* ionIoni =
new
G4ionIonisation
();
160
ionIoni->
SetStepFunction
(0.1, 1*
um
);
161
ph->
RegisterProcess
(ionIoni, particle);
162
ph->
RegisterProcess
(
new
G4NuclearStopping
(), particle);
163
164
}
else
if
( particleName ==
"GenericIon"
) {
165
166
ph->
RegisterProcess
(
new
G4hMultipleScattering
(), particle);
167
G4ionIonisation
* ionIoni =
new
G4ionIonisation
();
168
ionIoni->
SetEmModel
(
new
G4IonParametrisedLossModel
());
169
ionIoni->
SetStepFunction
(0.1, 1*
um
);
170
ph->
RegisterProcess
(ionIoni, particle);
171
ph->
RegisterProcess
(
new
G4NuclearStopping
(), particle);
172
173
}
else
if
((!particle->
IsShortLived
()) &&
174
(particle->
GetPDGCharge
() != 0.0) &&
175
(particle->
GetParticleName
() !=
"chargedgeantino"
)) {
176
177
//all others charged particles except geantino
178
ph->
RegisterProcess
(
new
G4hMultipleScattering
(), particle);
179
ph->
RegisterProcess
(
new
G4hIonisation
(), particle);
180
}
181
}
182
183
// Deexcitation
184
//
185
G4VAtomDeexcitation
* de =
new
G4UAtomicDeexcitation
();
186
G4LossTableManager::Instance
()->
SetAtomDeexcitation
(de);
187
188
}
189
190
//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
191
geant4
tree
geant4-10.6-release
examples
extended
electromagnetic
TestEm11
src
PhysListEmStandard.cc
Built by
Jin Huang
. updated:
Wed Jun 29 2022 17:24:51
using
1.8.2 with
ECCE GitHub integration