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
TestEm12
src
electromagnetic/TestEm12/src/PhysListEmStandardSS.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/TestEm12/src/PhysListEmStandardSS.cc
27
/// \brief Implementation of the PhysListEmStandardSS class
28
//
29
// $Id: PhysListEmStandardSS.cc 77092 2013-11-21 10:50:40Z gcosmo $
30
//
31
//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
32
//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
33
34
#include "PhysListEmStandardSS.hh"
35
36
#include "
G4ParticleDefinition.hh
"
37
#include "
G4ProcessManager.hh
"
38
39
#include "
G4ComptonScattering.hh
"
40
#include "
G4GammaConversion.hh
"
41
#include "
G4PhotoElectricEffect.hh
"
42
43
#include "
G4CoulombScattering.hh
"
44
#include "
G4eCoulombScatteringModel.hh
"
45
#include "
G4IonCoulombScatteringModel.hh
"
46
47
#include "
G4eIonisation.hh
"
48
#include "
G4eBremsstrahlung.hh
"
49
#include "
G4eplusAnnihilation.hh
"
50
51
#include "
G4MuIonisation.hh
"
52
#include "
G4MuBremsstrahlung.hh
"
53
#include "
G4MuPairProduction.hh
"
54
55
#include "
G4hIonisation.hh
"
56
#include "
G4ionIonisation.hh
"
57
58
#include "
G4EmProcessOptions.hh
"
59
60
#include "
G4SystemOfUnits.hh
"
61
62
//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
63
64
PhysListEmStandardSS::PhysListEmStandardSS
(
const
G4String
&
name
)
65
:
G4VPhysicsConstructor
(name)
66
{}
67
68
//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
69
70
PhysListEmStandardSS::~PhysListEmStandardSS
()
71
{}
72
73
//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
74
75
void
PhysListEmStandardSS::ConstructProcess
()
76
{
77
// Add standard EM Processes
78
79
aParticleIterator
->reset();
80
while
( (*
aParticleIterator
)() ){
81
G4ParticleDefinition
* particle =
aParticleIterator
->value();
82
G4ProcessManager
* pmanager = particle->
GetProcessManager
();
83
G4String
particleName = particle->
GetParticleName
();
84
85
if
(particleName ==
"gamma"
) {
86
// gamma
87
pmanager->
AddDiscreteProcess
(
new
G4PhotoElectricEffect
);
88
pmanager->
AddDiscreteProcess
(
new
G4ComptonScattering
);
89
pmanager->
AddDiscreteProcess
(
new
G4GammaConversion
);
90
91
}
else
if
(particleName ==
"e-"
) {
92
//electron
93
pmanager->
AddProcess
(
new
G4eIonisation
, -1, 1, 1);
94
pmanager->
AddProcess
(
new
G4eBremsstrahlung
, -1, 2, 2);
95
96
G4CoulombScattering
* cs =
new
G4CoulombScattering
();
97
G4eCoulombScatteringModel
*
model
=
98
new
G4eCoulombScatteringModel
();
99
model->
SetLowEnergyThreshold
(10*
eV
);
100
cs->
SetEmModel
(model, 1);
101
pmanager->
AddDiscreteProcess
(cs);
102
103
}
else
if
(particleName ==
"e+"
) {
104
//positron
105
pmanager->
AddProcess
(
new
G4eIonisation
, -1, 1, 1);
106
pmanager->
AddProcess
(
new
G4eBremsstrahlung
, -1, 2, 2);
107
pmanager->
AddProcess
(
new
G4eplusAnnihilation
, 0,-1, 3);
108
pmanager->
AddDiscreteProcess
(
new
G4CoulombScattering
);
109
110
}
else
if
(particleName ==
"mu+"
||
111
particleName ==
"mu-"
) {
112
//muon
113
pmanager->
AddProcess
(
new
G4MuIonisation
, -1, 1, 1);
114
pmanager->
AddProcess
(
new
G4MuBremsstrahlung
, -1, 2, 2);
115
pmanager->
AddProcess
(
new
G4MuPairProduction
, -1, 3, 3);
116
pmanager->
AddDiscreteProcess
(
new
G4CoulombScattering
);
117
118
}
else
if
(particleName ==
"alpha"
|| particleName ==
"He3"
) {
119
pmanager->
AddProcess
(
new
G4ionIonisation
, -1, 1, 1);
120
G4CoulombScattering
* cs =
new
G4CoulombScattering
();
121
//cs->AddEmModel(0, new G4IonCoulombScatteringModel());
122
cs->
SetBuildTableFlag
(
false
);
123
pmanager->
AddDiscreteProcess
(cs);
124
125
}
else
if
(particleName ==
"GenericIon"
) {
126
pmanager->
AddProcess
(
new
G4ionIonisation
, -1, 1, 1);
127
G4CoulombScattering
* cs =
new
G4CoulombScattering
();
128
//cs->AddEmModel(0, new G4IonCoulombScatteringModel());
129
cs->
SetBuildTableFlag
(
false
);
130
pmanager->
AddDiscreteProcess
(cs);
131
132
}
else
if
((!particle->
IsShortLived
()) &&
133
(particle->
GetPDGCharge
() != 0.0) &&
134
(particle->
GetParticleName
() !=
"chargedgeantino"
)) {
135
//all others charged particles except geantino
136
G4CoulombScattering
* cs =
new
G4CoulombScattering
();
137
pmanager->
AddProcess
(
new
G4hIonisation
, -1, 1, 1);
138
cs->
SetBuildTableFlag
(
false
);
139
pmanager->
AddDiscreteProcess
(cs);
140
}
141
}
142
143
// Em options
144
//
145
// Main options and setting parameters are shown here.
146
// Several of them have default values.
147
//
148
G4EmProcessOptions
emOptions;
149
150
//physics tables
151
//
152
emOptions.
SetMinEnergy
(100*
eV
);
//default
153
emOptions.
SetMaxEnergy
(100*
TeV
);
//default
154
emOptions.
SetDEDXBinning
(12*20);
//default=12*7
155
emOptions.
SetLambdaBinning
(12*20);
//default=12*7
156
// emOptions.SetSplineFlag(true); //default
157
158
//energy loss
159
//
160
//emOptions.SetStepFunction(0.2, 100*um); //default=(0.2, 1*mm)
161
//emOptions.SetLinearLossLimit(1.e-2); //default
162
163
//ionization
164
//
165
emOptions.
SetSubCutoff
(
false
);
//default
166
167
// scattering
168
emOptions.
SetPolarAngleLimit
(0.0);
169
}
170
171
//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
172
G4eCoulombScatteringModel
Definition:
G4eCoulombScatteringModel.hh:78
G4eCoulombScatteringModel.hh
G4eIonisation.hh
G4PhotoElectricEffect
Definition:
G4PhotoElectricEffect.hh:79
PhysListEmStandardSS::PhysListEmStandardSS
PhysListEmStandardSS(const G4String &name="standardSS")
Definition:
electromagnetic/TestEm1/src/PhysListEmStandardSS.cc:63
G4VEmProcess::SetBuildTableFlag
void SetBuildTableFlag(G4bool val)
Definition:
G4VEmProcess.hh:658
G4EmProcessOptions
Definition:
G4EmProcessOptions.hh:63
G4ComptonScattering.hh
G4hIonisation.hh
G4EmProcessOptions::SetMinEnergy
void SetMinEnergy(G4double val)
Definition:
G4EmProcessOptions.cc:109
G4MuPairProduction.hh
G4ProcessManager::AddDiscreteProcess
G4int AddDiscreteProcess(G4VProcess *aProcess, G4int ord=ordDefault)
name
const XML_Char * name
Definition:
include/expat.h:151
G4ComptonScattering
Definition:
G4ComptonScattering.hh:71
G4MuIonisation.hh
G4ParticleDefinition
Definition:
G4ParticleDefinition.hh:111
G4MuBremsstrahlung.hh
G4ParticleDefinition::GetProcessManager
G4ProcessManager * GetProcessManager() const
G4ParticleDefinition::GetParticleName
const G4String & GetParticleName() const
Definition:
G4ParticleDefinition.hh:159
G4CoulombScattering
Definition:
G4CoulombScattering.hh:55
G4ionIonisation.hh
G4EmProcessOptions::SetDEDXBinning
void SetDEDXBinning(G4int val)
Definition:
G4EmProcessOptions.cc:137
G4CoulombScattering.hh
G4VEmProcess::SetEmModel
void SetEmModel(G4VEmModel *, G4int index=1)
Definition:
G4VEmProcess.cc:198
G4EmProcessOptions::SetLambdaBinning
void SetLambdaBinning(G4int val)
Definition:
G4EmProcessOptions.cc:151
G4GammaConversion
Definition:
G4GammaConversion.hh:75
G4GammaConversion.hh
aParticleIterator
#define aParticleIterator
Definition:
G4VPhysicsConstructor.hh:119
G4ProcessManager::AddProcess
G4int AddProcess(G4VProcess *aProcess, G4int ordAtRestDoIt=ordInActive, G4int ordAlongSteptDoIt=ordInActive, G4int ordPostStepDoIt=ordInActive)
Definition:
G4ProcessManager.cc:410
G4eCoulombScatteringModel::SetLowEnergyThreshold
void SetLowEnergyThreshold(G4double val)
Definition:
G4eCoulombScatteringModel.hh:195
model
const XML_Char XML_Content * model
Definition:
include/expat.h:151
python.hepunit.eV
eV
Definition:
hepunit.py:118
G4IonCoulombScatteringModel.hh
G4ProcessManager.hh
G4ParticleDefinition.hh
G4PhotoElectricEffect.hh
G4ionIonisation
Definition:
G4ionIonisation.hh:78
G4MuBremsstrahlung
Definition:
G4MuBremsstrahlung.hh:78
G4EmProcessOptions::SetMaxEnergy
void SetMaxEnergy(G4double val)
Definition:
G4EmProcessOptions.cc:116
PhysListEmStandardSS::ConstructProcess
virtual void ConstructProcess()
Definition:
electromagnetic/TestEm1/src/PhysListEmStandardSS.cc:74
G4MuIonisation
Definition:
G4MuIonisation.hh:85
G4ParticleDefinition::IsShortLived
G4bool IsShortLived() const
Definition:
G4ParticleDefinition.hh:196
G4eBremsstrahlung
Definition:
G4eBremsstrahlung.hh:81
G4eIonisation
Definition:
G4eIonisation.hh:80
G4eplusAnnihilation.hh
G4MuPairProduction
Definition:
G4MuPairProduction.hh:74
G4eBremsstrahlung.hh
G4ProcessManager
Definition:
G4ProcessManager.hh:106
G4eplusAnnihilation
Definition:
G4eplusAnnihilation.hh:65
G4SystemOfUnits.hh
G4ParticleDefinition::GetPDGCharge
G4double GetPDGCharge() const
Definition:
G4ParticleDefinition.hh:163
PhysListEmStandardSS::~PhysListEmStandardSS
virtual ~PhysListEmStandardSS()
Definition:
electromagnetic/TestEm1/src/PhysListEmStandardSS.cc:69
G4EmProcessOptions::SetSubCutoff
void SetSubCutoff(G4bool val, const G4Region *r=0)
Definition:
G4EmProcessOptions.cc:88
G4EmProcessOptions.hh
G4VPhysicsConstructor
Definition:
G4VPhysicsConstructor.hh:121
G4hIonisation
Definition:
G4hIonisation.hh:85
python.hepunit.TeV
TeV
Definition:
hepunit.py:121
G4String
Definition:
examples/extended/parallel/TopC/ParN02/AnnotatedFiles/G4String.hh:45
G4EmProcessOptions::SetPolarAngleLimit
void SetPolarAngleLimit(G4double val)
Definition:
G4EmProcessOptions.cc:369
Generated on Wed Apr 30 2014 15:55:21 for Geant4.10 by
1.8.7