Geant4
9.6.p02
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
geant4_9_6_p02
examples
advanced
amsEcal
src
PhysListEmStandard.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
// $Id$
27
//
28
//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
29
//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
30
31
#include "PhysListEmStandard.hh"
32
#include "
G4SystemOfUnits.hh
"
33
#include "
G4ParticleDefinition.hh
"
34
#include "
G4ProcessManager.hh
"
35
36
#include "
G4ComptonScattering.hh
"
37
#include "
G4GammaConversion.hh
"
38
#include "
G4PhotoElectricEffect.hh
"
39
40
#include "
G4eMultipleScattering.hh
"
41
#include "
G4hMultipleScattering.hh
"
42
43
#include "
G4eIonisation.hh
"
44
#include "
G4eBremsstrahlung.hh
"
45
#include "
G4eplusAnnihilation.hh
"
46
47
#include "
G4MuIonisation.hh
"
48
#include "
G4MuBremsstrahlung.hh
"
49
#include "
G4MuPairProduction.hh
"
50
51
#include "
G4hIonisation.hh
"
52
#include "
G4hBremsstrahlung.hh
"
53
#include "
G4hPairProduction.hh
"
54
55
#include "
G4ionIonisation.hh
"
56
57
#include "
G4EmProcessOptions.hh
"
58
#include "
G4MscStepLimitType.hh
"
59
60
//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
61
62
PhysListEmStandard::PhysListEmStandard
(
const
G4String
&
name
)
63
:
G4VPhysicsConstructor
(name)
64
{}
65
66
//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
67
68
PhysListEmStandard::~PhysListEmStandard
()
69
{}
70
71
//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
72
73
void
PhysListEmStandard::ConstructProcess
()
74
{
75
// Add standard EM Processes
76
77
theParticleIterator
->
reset
();
78
while
( (*
theParticleIterator
)() ){
79
G4ParticleDefinition
* particle =
theParticleIterator
->
value
();
80
G4ProcessManager
* pmanager = particle->
GetProcessManager
();
81
G4String
particleName = particle->
GetParticleName
();
82
83
if
(particleName ==
"gamma"
) {
84
// gamma
85
pmanager->
AddDiscreteProcess
(
new
G4PhotoElectricEffect
);
86
pmanager->
AddDiscreteProcess
(
new
G4ComptonScattering
);
87
pmanager->
AddDiscreteProcess
(
new
G4GammaConversion
);
88
89
}
else
if
(particleName ==
"e-"
) {
90
//electron
91
pmanager->
AddProcess
(
new
G4eMultipleScattering
, -1, 1, 1);
92
pmanager->
AddProcess
(
new
G4eIonisation
, -1, 2, 2);
93
pmanager->
AddProcess
(
new
G4eBremsstrahlung
(), -1, 3, 3);
94
95
}
else
if
(particleName ==
"e+"
) {
96
//positron
97
pmanager->
AddProcess
(
new
G4eMultipleScattering
, -1, 1, 1);
98
pmanager->
AddProcess
(
new
G4eIonisation
, -1, 2, 2);
99
pmanager->
AddProcess
(
new
G4eBremsstrahlung
(), -1, 3, 3);
100
pmanager->
AddProcess
(
new
G4eplusAnnihilation
, 0,-1, 4);
101
102
}
else
if
( particleName ==
"mu-"
||
103
particleName ==
"mu+"
) {
104
//muon
105
pmanager->
AddProcess
(
new
G4hMultipleScattering
, -1, 1, 1);
106
pmanager->
AddProcess
(
new
G4MuIonisation
, -1, 2, 2);
107
pmanager->
AddProcess
(
new
G4MuBremsstrahlung
, -1, 3, 3);
108
pmanager->
AddProcess
(
new
G4MuPairProduction
, -1, 4, 4);
109
110
}
else
if
( particleName ==
"proton"
||
111
particleName ==
"pi-"
||
112
particleName ==
"pi+"
) {
113
//proton
114
pmanager->
AddProcess
(
new
G4hMultipleScattering
, -1, 1, 1);
115
pmanager->
AddProcess
(
new
G4hIonisation
, -1, 2, 2);
116
pmanager->
AddProcess
(
new
G4hBremsstrahlung
, -1, 3, 3);
117
pmanager->
AddProcess
(
new
G4hPairProduction
, -1, 4, 4);
118
119
}
else
if
( particleName ==
"alpha"
||
120
particleName ==
"He3"
||
121
particleName ==
"GenericIon"
) {
122
//Ions
123
pmanager->
AddProcess
(
new
G4hMultipleScattering
, -1, 1, 1);
124
pmanager->
AddProcess
(
new
G4ionIonisation
, -1, 2, 2);
125
126
}
else
if
((!particle->
IsShortLived
()) &&
127
(particle->
GetPDGCharge
() != 0.0) &&
128
(particle->
GetParticleName
() !=
"chargedgeantino"
)) {
129
//all others charged particles except geantino
130
pmanager->
AddProcess
(
new
G4hMultipleScattering
, -1, 1, 1);
131
pmanager->
AddProcess
(
new
G4hIonisation
, -1, 2, 2);
132
}
133
}
134
135
// Em options
136
//
137
// Main options and setting parameters are shown here.
138
// Several of them have default values.
139
//
140
G4EmProcessOptions
emOptions;
141
142
//physics tables
143
//
144
emOptions.
SetMinEnergy
(100*
eV
);
//default
145
emOptions.
SetMaxEnergy
(100*
TeV
);
//default
146
emOptions.
SetDEDXBinning
(12*10);
//default=12*7
147
emOptions.
SetLambdaBinning
(12*10);
//default=12*7
148
emOptions.
SetSplineFlag
(
true
);
//default
149
150
//multiple coulomb scattering
151
//
152
emOptions.
SetMscStepLimitation
(
fUseSafety
);
//default=fUseSafety
153
emOptions.
SetMscRangeFactor
(0.04);
//default
154
emOptions.
SetMscGeomFactor
(2.5);
//default
155
emOptions.
SetSkin
(3.);
//default
156
157
//energy loss
158
//
159
emOptions.
SetStepFunction
(0.2, 100*um);
//default=(0.2, 1*mm)
160
emOptions.
SetLinearLossLimit
(1.
e
-2);
//default
161
162
//ionization
163
//
164
emOptions.
SetSubCutoff
(
false
);
//default=false
165
166
//verbose
167
//
168
emOptions.
SetVerbose
(0);
169
}
170
171
//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
172
Generated on Sat May 25 2013 14:31:58 for Geant4 by
1.8.4