Geant4.10
Main Page
Related Pages
Modules
Namespaces
Data Structures
Files
File List
Globals
All
Data Structures
Namespaces
Files
Functions
Variables
Typedefs
Enumerations
Enumerator
Friends
Macros
Groups
Pages
geant4.10.00.p01
examples
extended
electromagnetic
TestEm14
src
TestEm14/src/PhysListEmPenelope.cc
Go to the documentation of this file.
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
/// \file electromagnetic/TestEm14/src/PhysListEmPenelope.cc
27
/// \brief Implementation of the PhysListEmPenelope class
28
//
29
//
30
// $Id: PhysListEmPenelope.cc 68585 2013-04-01 23:35:07Z adotti $
31
//
32
//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
33
//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
34
35
#include "PhysListEmPenelope.hh"
36
#include "
G4ParticleDefinition.hh
"
37
#include "
G4ProcessManager.hh
"
38
39
// gamma
40
41
#include "
G4PhotoElectricEffect.hh
"
42
#include "
G4PenelopePhotoElectricModel.hh
"
43
44
#include "
G4ComptonScattering.hh
"
45
#include "
G4PenelopeComptonModel.hh
"
46
47
#include "
G4GammaConversion.hh
"
48
#include "
G4PenelopeGammaConversionModel.hh
"
49
50
#include "
G4RayleighScattering.hh
"
51
#include "
G4PenelopeRayleighModel.hh
"
52
53
// e-
54
55
#include "
G4eIonisation.hh
"
56
#include "
G4PenelopeIonisationModel.hh
"
57
#include "
G4UniversalFluctuation.hh
"
58
59
#include "
G4eBremsstrahlung.hh
"
60
#include "
G4PenelopeBremsstrahlungModel.hh
"
61
62
// e+
63
64
#include "
G4eplusAnnihilation.hh
"
65
#include "
G4PenelopeAnnihilationModel.hh
"
66
67
// mu
68
69
#include "
G4MuIonisation.hh
"
70
#include "
G4MuBremsstrahlung.hh
"
71
#include "
G4MuPairProduction.hh
"
72
73
// hadrons, ions
74
75
#include "
G4hIonisation.hh
"
76
#include "
G4ionIonisation.hh
"
77
78
#include "
G4EmProcessOptions.hh
"
79
#include "
G4LossTableManager.hh
"
80
#include "
G4UAtomicDeexcitation.hh
"
81
82
#include "
G4SystemOfUnits.hh
"
83
84
//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
85
86
PhysListEmPenelope::PhysListEmPenelope
(
const
G4String
&
name
)
87
:
G4VPhysicsConstructor
(name)
88
{ }
89
90
//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
91
92
PhysListEmPenelope::~PhysListEmPenelope
()
93
{ }
94
95
//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
96
97
void
PhysListEmPenelope::ConstructProcess
()
98
{
99
// Add standard EM Processes
100
101
aParticleIterator
->reset();
102
while
( (*
aParticleIterator
)() ){
103
G4ParticleDefinition
* particle =
aParticleIterator
->value();
104
G4ProcessManager
* pmanager = particle->
GetProcessManager
();
105
G4String
particleName = particle->
GetParticleName
();
106
107
//Applicability range for Penelope models
108
//for higher energies, the Standard models are used
109
G4double
highEnergyLimit = 1*
GeV
;
110
111
if
(particleName ==
"gamma"
) {
112
// gamma
113
114
G4PhotoElectricEffect
* phot =
new
G4PhotoElectricEffect
();
115
G4PenelopePhotoElectricModel
*
116
photModel =
new
G4PenelopePhotoElectricModel
();
117
photModel->
SetHighEnergyLimit
(highEnergyLimit);
118
phot->
AddEmModel
(0, photModel);
119
pmanager->
AddDiscreteProcess
(phot);
120
121
G4ComptonScattering
* compt =
new
G4ComptonScattering
();
122
G4PenelopeComptonModel
*
123
comptModel =
new
G4PenelopeComptonModel
();
124
comptModel->
SetHighEnergyLimit
(highEnergyLimit);
125
compt->
AddEmModel
(0, comptModel);
126
pmanager->
AddDiscreteProcess
(compt);
127
128
G4GammaConversion
* conv =
new
G4GammaConversion
();
129
G4PenelopeGammaConversionModel
*
130
convModel =
new
G4PenelopeGammaConversionModel
();
131
convModel->
SetHighEnergyLimit
(highEnergyLimit);
132
conv->
AddEmModel
(0, convModel);
133
pmanager->
AddDiscreteProcess
(conv);
134
135
G4RayleighScattering
* rayl =
new
G4RayleighScattering
();
136
G4PenelopeRayleighModel
*
137
raylModel =
new
G4PenelopeRayleighModel
();
138
raylModel->
SetHighEnergyLimit
(highEnergyLimit);
139
rayl->
AddEmModel
(0, raylModel);
140
pmanager->
AddDiscreteProcess
(rayl);
141
142
}
else
if
(particleName ==
"e-"
) {
143
//electron
144
145
G4eIonisation
* eIoni =
new
G4eIonisation
();
146
G4PenelopeIonisationModel
*
147
eIoniModel =
new
G4PenelopeIonisationModel
();
148
eIoniModel->
SetHighEnergyLimit
(highEnergyLimit);
149
eIoni->
AddEmModel
(0, eIoniModel,
new
G4UniversalFluctuation
() );
150
pmanager->
AddProcess
(eIoni, -1,-1, 1);
151
152
G4eBremsstrahlung
* eBrem =
new
G4eBremsstrahlung
();
153
G4PenelopeBremsstrahlungModel
*
154
eBremModel =
new
G4PenelopeBremsstrahlungModel
();
155
eBremModel->
SetHighEnergyLimit
(highEnergyLimit);
156
eBrem->
AddEmModel
(0, eBremModel);
157
pmanager->
AddProcess
(eBrem, -1,-1, 2);
158
159
}
else
if
(particleName ==
"e+"
) {
160
//positron
161
G4eIonisation
* eIoni =
new
G4eIonisation
();
162
G4PenelopeIonisationModel
*
163
eIoniModel =
new
G4PenelopeIonisationModel
();
164
eIoniModel->
SetHighEnergyLimit
(highEnergyLimit);
165
eIoni->
AddEmModel
(0, eIoniModel,
new
G4UniversalFluctuation
() );
166
pmanager->
AddProcess
(eIoni, -1,-1, 1);
167
168
G4eBremsstrahlung
* eBrem =
new
G4eBremsstrahlung
();
169
G4PenelopeBremsstrahlungModel
*
170
eBremModel =
new
G4PenelopeBremsstrahlungModel
();
171
eBremModel->
SetHighEnergyLimit
(highEnergyLimit);
172
eBrem->
AddEmModel
(0, eBremModel);
173
pmanager->
AddProcess
(eBrem, -1,-1, 2);
174
175
G4eplusAnnihilation
* eAnni =
new
G4eplusAnnihilation
();
176
G4PenelopeAnnihilationModel
*
177
eAnniModel =
new
G4PenelopeAnnihilationModel
();
178
eAnniModel->
SetHighEnergyLimit
(highEnergyLimit);
179
eAnni->
AddEmModel
(0, eAnniModel);
180
pmanager->
AddProcess
(eAnni, 0,-1, 3);
181
182
}
else
if
( particleName ==
"mu+"
||
183
particleName ==
"mu-"
) {
184
//muon
185
pmanager->
AddProcess
(
new
G4MuIonisation
, -1,-1, 1);
186
pmanager->
AddProcess
(
new
G4MuBremsstrahlung
, -1,-1, 2);
187
pmanager->
AddProcess
(
new
G4MuPairProduction
, -1,-1, 3);
188
189
}
else
if
( particleName ==
"alpha"
|| particleName ==
"GenericIon"
) {
190
pmanager->
AddProcess
(
new
G4ionIonisation
, -1,-1, 1);
191
192
}
else
if
((!particle->
IsShortLived
()) &&
193
(particle->
GetPDGCharge
() != 0.0) &&
194
(particle->
GetParticleName
() !=
"chargedgeantino"
)) {
195
//all others charged particles except geantino
196
pmanager->
AddProcess
(
new
G4hIonisation
, -1,-1, 1);
197
}
198
}
199
200
// Deexcitation
201
//
202
G4VAtomDeexcitation
* de =
new
G4UAtomicDeexcitation
();
203
de->
SetFluo
(
true
);
204
de->
SetAuger
(
false
);
205
de->
SetPIXE
(
false
);
206
G4LossTableManager::Instance
()->
SetAtomDeexcitation
(de);
207
}
208
209
//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
210
G4PenelopeAnnihilationModel.hh
G4eIonisation.hh
python.hepunit.GeV
GeV
Definition:
hepunit.py:120
G4LossTableManager::Instance
static G4LossTableManager * Instance()
Definition:
G4LossTableManager.cc:112
G4UniversalFluctuation.hh
G4PenelopeAnnihilationModel
Definition:
G4PenelopeAnnihilationModel.hh:56
G4PhotoElectricEffect
Definition:
G4PhotoElectricEffect.hh:79
G4VAtomDeexcitation::SetAuger
void SetAuger(G4bool)
Definition:
G4VAtomDeexcitation.hh:208
G4ComptonScattering.hh
G4PenelopeGammaConversionModel
Definition:
G4PenelopeGammaConversionModel.hh:57
G4PenelopeBremsstrahlungModel
Definition:
G4PenelopeBremsstrahlungModel.hh:63
G4hIonisation.hh
G4MuPairProduction.hh
G4ProcessManager::AddDiscreteProcess
G4int AddDiscreteProcess(G4VProcess *aProcess, G4int ord=ordDefault)
PhysListEmPenelope::PhysListEmPenelope
PhysListEmPenelope(const G4String &name="penelope")
Definition:
TestEm13/src/PhysListEmPenelope.cc:82
G4UAtomicDeexcitation
Definition:
G4UAtomicDeexcitation.hh:61
name
const XML_Char * name
Definition:
include/expat.h:151
G4ComptonScattering
Definition:
G4ComptonScattering.hh:71
G4LossTableManager.hh
G4PenelopeIonisationModel
Definition:
G4PenelopeIonisationModel.hh:65
G4PenelopePhotoElectricModel.hh
G4MuIonisation.hh
G4ParticleDefinition
Definition:
G4ParticleDefinition.hh:111
G4MuBremsstrahlung.hh
G4ParticleDefinition::GetProcessManager
G4ProcessManager * GetProcessManager() const
G4ParticleDefinition::GetParticleName
const G4String & GetParticleName() const
Definition:
G4ParticleDefinition.hh:159
G4VEmModel::SetHighEnergyLimit
void SetHighEnergyLimit(G4double)
Definition:
G4VEmModel.hh:683
G4ionIonisation.hh
G4PenelopeBremsstrahlungModel.hh
G4GammaConversion
Definition:
G4GammaConversion.hh:75
G4GammaConversion.hh
aParticleIterator
#define aParticleIterator
Definition:
G4VPhysicsConstructor.hh:119
G4VAtomDeexcitation
Definition:
G4VAtomDeexcitation.hh:63
G4ProcessManager::AddProcess
G4int AddProcess(G4VProcess *aProcess, G4int ordAtRestDoIt=ordInActive, G4int ordAlongSteptDoIt=ordInActive, G4int ordPostStepDoIt=ordInActive)
Definition:
G4ProcessManager.cc:410
G4ProcessManager.hh
G4ParticleDefinition.hh
G4VEnergyLossProcess::AddEmModel
void AddEmModel(G4int, G4VEmModel *, G4VEmFluctuationModel *fluc=0, const G4Region *region=0)
Definition:
G4VEnergyLossProcess.cc:344
G4PenelopeRayleighModel.hh
G4PhotoElectricEffect.hh
G4PenelopeGammaConversionModel.hh
G4ionIonisation
Definition:
G4ionIonisation.hh:78
G4MuBremsstrahlung
Definition:
G4MuBremsstrahlung.hh:78
G4UAtomicDeexcitation.hh
G4MuIonisation
Definition:
G4MuIonisation.hh:85
G4PenelopeIonisationModel.hh
PhysListEmPenelope::~PhysListEmPenelope
~PhysListEmPenelope()
Definition:
TestEm13/src/PhysListEmPenelope.cc:88
G4ParticleDefinition::IsShortLived
G4bool IsShortLived() const
Definition:
G4ParticleDefinition.hh:196
G4RayleighScattering.hh
G4eBremsstrahlung
Definition:
G4eBremsstrahlung.hh:81
G4PenelopeComptonModel
Definition:
G4PenelopeComptonModel.hh:63
PhysListEmPenelope::ConstructProcess
virtual void ConstructProcess()
Definition:
TestEm13/src/PhysListEmPenelope.cc:93
G4eIonisation
Definition:
G4eIonisation.hh:80
G4VEmProcess::AddEmModel
void AddEmModel(G4int, G4VEmModel *, const G4Region *region=0)
Definition:
G4VEmProcess.cc:188
G4eplusAnnihilation.hh
G4MuPairProduction
Definition:
G4MuPairProduction.hh:74
G4UniversalFluctuation
Definition:
G4UniversalFluctuation.hh:62
G4eBremsstrahlung.hh
G4ProcessManager
Definition:
G4ProcessManager.hh:106
G4eplusAnnihilation
Definition:
G4eplusAnnihilation.hh:65
G4RayleighScattering
Definition:
G4RayleighScattering.hh:51
G4VAtomDeexcitation::SetPIXE
void SetPIXE(G4bool)
Definition:
G4VAtomDeexcitation.hh:219
G4double
double G4double
Definition:
G4Types.hh:76
G4PenelopePhotoElectricModel
Definition:
G4PenelopePhotoElectricModel.hh:59
G4SystemOfUnits.hh
G4PenelopeComptonModel.hh
G4ParticleDefinition::GetPDGCharge
G4double GetPDGCharge() const
Definition:
G4ParticleDefinition.hh:163
G4EmProcessOptions.hh
G4LossTableManager::SetAtomDeexcitation
void SetAtomDeexcitation(G4VAtomDeexcitation *)
Definition:
G4LossTableManager.cc:1314
G4VPhysicsConstructor
Definition:
G4VPhysicsConstructor.hh:121
G4PenelopeRayleighModel
Definition:
G4PenelopeRayleighModel.hh:58
G4hIonisation
Definition:
G4hIonisation.hh:85
G4VAtomDeexcitation::SetFluo
void SetFluo(G4bool)
Definition:
G4VAtomDeexcitation.hh:198
G4String
Definition:
examples/extended/parallel/TopC/ParN02/AnnotatedFiles/G4String.hh:45
Generated on Wed Apr 30 2014 15:55:21 for Geant4.10 by
1.8.7