Geant4
10.03.p03
Main Page
Related Pages
Modules
Namespaces
Classes
Files
File List
File Members
All
Classes
Namespaces
Files
Functions
Variables
Typedefs
Enumerations
Enumerator
Friends
Macros
Groups
Pages
PhysListEmLivermore.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
//
28
//
29
//
30
// $Id: PhysListEmLivermore.cc 100279 2016-10-17 08:34:54Z gcosmo $
31
//
32
//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
33
//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
34
35
#include "PhysListEmLivermore.hh"
36
#include "
G4ParticleDefinition.hh
"
37
#include "
G4ProcessManager.hh
"
38
39
// gamma
40
41
#include "
G4PhotoElectricEffect.hh
"
42
#include "
G4LivermorePhotoElectricModel.hh
"
43
44
#include "
G4ComptonScattering.hh
"
45
#include "
G4LivermoreComptonModel.hh
"
46
47
#include "
G4GammaConversion.hh
"
48
#include "
G4LivermoreGammaConversionModel.hh
"
49
50
#include "
G4RayleighScattering.hh
"
51
#include "
G4LivermoreRayleighModel.hh
"
52
53
// e-
54
55
#include "
G4eIonisation.hh
"
56
#include "
G4LivermoreIonisationModel.hh
"
57
#include "
G4UniversalFluctuation.hh
"
58
59
#include "
G4eBremsstrahlung.hh
"
60
#include "
G4LivermoreBremsstrahlungModel.hh
"
61
62
// e+
63
64
#include "
G4eplusAnnihilation.hh
"
65
66
// mu
67
68
#include "
G4MuIonisation.hh
"
69
#include "
G4MuBremsstrahlung.hh
"
70
#include "
G4MuPairProduction.hh
"
71
72
// hadrons, ions
73
74
#include "
G4hIonisation.hh
"
75
#include "
G4ionIonisation.hh
"
76
77
#include "
G4SystemOfUnits.hh
"
78
79
//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
80
81
PhysListEmLivermore::PhysListEmLivermore
(
const
G4String
&
name
)
82
:
G4VPhysicsConstructor
(name)
83
{ }
84
85
//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
86
87
PhysListEmLivermore::~PhysListEmLivermore
()
88
{ }
89
90
//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
91
92
void
PhysListEmLivermore::ConstructProcess
()
93
{
94
// Add standard EM Processes
95
96
auto
particleIterator
=
GetParticleIterator
();
97
particleIterator
->reset();
98
while
( (*
particleIterator
)() ){
99
G4ParticleDefinition
* particle =
particleIterator
->value();
100
G4ProcessManager
* pmanager = particle->
GetProcessManager
();
101
G4String
particleName = particle->
GetParticleName
();
102
103
//Applicability range for Livermore models
104
//for higher energies, the Standard models are used
105
G4double
highEnergyLimit = 1*
GeV
;
106
107
if
(particleName ==
"gamma"
) {
108
// gamma
109
110
G4PhotoElectricEffect
* phot =
new
G4PhotoElectricEffect
();
111
G4LivermorePhotoElectricModel
*
112
photModel =
new
G4LivermorePhotoElectricModel
();
113
photModel->
SetHighEnergyLimit
(highEnergyLimit);
114
phot->
AddEmModel
(0, photModel);
115
pmanager->
AddDiscreteProcess
(phot);
116
117
G4ComptonScattering
* compt =
new
G4ComptonScattering
();
118
G4LivermoreComptonModel
*
119
comptModel =
new
G4LivermoreComptonModel
();
120
comptModel->
SetHighEnergyLimit
(highEnergyLimit);
121
compt->
AddEmModel
(0, comptModel);
122
pmanager->
AddDiscreteProcess
(compt);
123
124
G4GammaConversion
* conv =
new
G4GammaConversion
();
125
G4LivermoreGammaConversionModel
*
126
convModel =
new
G4LivermoreGammaConversionModel
();
127
convModel->
SetHighEnergyLimit
(highEnergyLimit);
128
conv->
AddEmModel
(0, convModel);
129
pmanager->
AddDiscreteProcess
(conv);
130
131
G4RayleighScattering
* rayl =
new
G4RayleighScattering
();
132
G4LivermoreRayleighModel
*
133
raylModel =
new
G4LivermoreRayleighModel
();
134
raylModel->
SetHighEnergyLimit
(highEnergyLimit);
135
rayl->
AddEmModel
(0, raylModel);
136
pmanager->
AddDiscreteProcess
(rayl);
137
138
}
else
if
(particleName ==
"e-"
) {
139
//electron
140
141
G4eIonisation
* eIoni =
new
G4eIonisation
();
142
G4LivermoreIonisationModel
*
143
eIoniModel =
new
G4LivermoreIonisationModel
();
144
eIoniModel->
SetHighEnergyLimit
(highEnergyLimit);
145
eIoni->
AddEmModel
(0, eIoniModel,
new
G4UniversalFluctuation
() );
146
pmanager->
AddProcess
(eIoni, -1,-1, 1);
147
148
G4eBremsstrahlung
* eBrem =
new
G4eBremsstrahlung
();
149
G4LivermoreBremsstrahlungModel
*
150
eBremModel =
new
G4LivermoreBremsstrahlungModel
();
151
eBremModel->
SetHighEnergyLimit
(highEnergyLimit);
152
eBrem->
AddEmModel
(0, eBremModel);
153
pmanager->
AddProcess
(eBrem, -1,-1, 2);
154
155
}
else
if
(particleName ==
"e+"
) {
156
//positron
157
pmanager->
AddProcess
(
new
G4eIonisation
, -1,-1, 1);
158
pmanager->
AddProcess
(
new
G4eBremsstrahlung
, -1,-1, 2);
159
pmanager->
AddProcess
(
new
G4eplusAnnihilation
, 0,-1, 3);
160
161
}
else
if
( particleName ==
"mu+"
||
162
particleName ==
"mu-"
) {
163
//muon
164
pmanager->
AddProcess
(
new
G4MuIonisation
, -1,-1, 1);
165
pmanager->
AddProcess
(
new
G4MuBremsstrahlung
, -1,-1, 2);
166
pmanager->
AddProcess
(
new
G4MuPairProduction
, -1,-1, 3);
167
168
}
else
if
( particleName ==
"alpha"
|| particleName ==
"GenericIon"
) {
169
pmanager->
AddProcess
(
new
G4ionIonisation
, -1,-1, 1);
170
171
}
else
if
((!particle->
IsShortLived
()) &&
172
(particle->
GetPDGCharge
() != 0.0) &&
173
(particle->
GetParticleName
() !=
"chargedgeantino"
)) {
174
//all others charged particles except geantino
175
pmanager->
AddProcess
(
new
G4hIonisation
, -1,-1, 1);
176
}
177
}
178
}
179
180
//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
181
name
const XML_Char * name
Definition:
expat.h:151
G4LivermoreBremsstrahlungModel
Definition:
G4LivermoreBremsstrahlungModel.hh:61
G4eIonisation.hh
G4UniversalFluctuation.hh
G4PhotoElectricEffect
Definition:
G4PhotoElectricEffect.hh:79
G4ComptonScattering.hh
G4hIonisation.hh
G4LivermoreIonisationModel.hh
PhysListEmLivermore::ConstructProcess
virtual void ConstructProcess()
Definition:
PhysListEmLivermore.cc:92
G4MuPairProduction.hh
G4ProcessManager::AddDiscreteProcess
G4int AddDiscreteProcess(G4VProcess *aProcess, G4int ord=ordDefault)
G4ComptonScattering
Definition:
G4ComptonScattering.hh:71
G4MuIonisation.hh
G4ParticleDefinition
Definition:
G4ParticleDefinition.hh:72
G4MuBremsstrahlung.hh
G4LivermoreComptonModel.hh
G4ParticleDefinition::GetParticleName
const G4String & GetParticleName() const
Definition:
G4ParticleDefinition.hh:120
G4VEmModel::SetHighEnergyLimit
void SetHighEnergyLimit(G4double)
Definition:
G4VEmModel.hh:732
G4LivermoreRayleighModel
Definition:
G4LivermoreRayleighModel.hh:39
G4ionIonisation.hh
G4VEnergyLossProcess::AddEmModel
void AddEmModel(G4int, G4VEmModel *, G4VEmFluctuationModel *fluc=0, const G4Region *region=nullptr)
Definition:
G4VEnergyLossProcess.cc:391
G4VPhysicsConstructor::GetParticleIterator
G4ParticleTable::G4PTblDicIterator * GetParticleIterator() const
Definition:
G4VPhysicsConstructor.cc:78
G4LivermoreIonisationModel
Definition:
G4LivermoreIonisationModel.hh:58
G4GammaConversion
Definition:
G4GammaConversion.hh:75
G4GammaConversion.hh
G4ProcessManager::AddProcess
G4int AddProcess(G4VProcess *aProcess, G4int ordAtRestDoIt=ordInActive, G4int ordAlongSteptDoIt=ordInActive, G4int ordPostStepDoIt=ordInActive)
Definition:
G4ProcessManager.cc:410
G4LivermoreBremsstrahlungModel.hh
G4ProcessManager.hh
G4ParticleDefinition.hh
G4LivermoreComptonModel
Definition:
G4LivermoreComptonModel.hh:48
G4PhotoElectricEffect.hh
G4ionIonisation
Definition:
G4ionIonisation.hh:78
G4MuBremsstrahlung
Definition:
G4MuBremsstrahlung.hh:78
G4LivermoreGammaConversionModel.hh
G4MuIonisation
Definition:
G4MuIonisation.hh:85
G4LivermoreRayleighModel.hh
G4ParticleDefinition::IsShortLived
G4bool IsShortLived() const
Definition:
G4ParticleDefinition.hh:157
G4VEmProcess::AddEmModel
void AddEmModel(G4int, G4VEmModel *, const G4Region *region=nullptr)
Definition:
G4VEmProcess.cc:203
G4RayleighScattering.hh
G4eBremsstrahlung
Definition:
G4eBremsstrahlung.hh:81
G4eIonisation
Definition:
G4eIonisation.hh:80
G4eplusAnnihilation.hh
G4LivermorePhotoElectricModel
Definition:
G4LivermorePhotoElectricModel.hh:50
G4ParticleDefinition::GetProcessManager
G4ProcessManager * GetProcessManager() const
Definition:
G4ParticleDefinition.cc:247
G4MuPairProduction
Definition:
G4MuPairProduction.hh:74
GeV
static constexpr double GeV
Definition:
G4SIunits.hh:217
G4UniversalFluctuation
Definition:
G4UniversalFluctuation.hh:64
particleIterator
std::vector< G4InuclElementaryParticle >::iterator particleIterator
Definition:
G4BigBanger.cc:65
G4eBremsstrahlung.hh
G4ProcessManager
Definition:
G4ProcessManager.hh:106
G4eplusAnnihilation
Definition:
G4eplusAnnihilation.hh:65
G4RayleighScattering
Definition:
G4RayleighScattering.hh:51
G4LivermoreGammaConversionModel
Definition:
G4LivermoreGammaConversionModel.hh:41
G4double
double G4double
Definition:
G4Types.hh:76
G4SystemOfUnits.hh
G4ParticleDefinition::GetPDGCharge
G4double GetPDGCharge() const
Definition:
G4ParticleDefinition.hh:124
PhysListEmLivermore::~PhysListEmLivermore
~PhysListEmLivermore()
Definition:
PhysListEmLivermore.cc:87
G4VPhysicsConstructor
Definition:
G4VPhysicsConstructor.hh:121
G4hIonisation
Definition:
G4hIonisation.hh:85
G4LivermorePhotoElectricModel.hh
PhysListEmLivermore::PhysListEmLivermore
PhysListEmLivermore(const G4String &name="livermore")
Definition:
PhysListEmLivermore.cc:81
G4String
Definition:
G4String.hh:45
source
geant4.10.03.p03
examples
extended
electromagnetic
TestEm13
src
PhysListEmLivermore.cc
Generated on Tue Nov 28 2017 21:43:48 for Geant4 by
1.8.5