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_WVI.cc
Go to the documentation of this file.
Or view
the newest version in sPHENIX GitHub for file PhysListEmStandard_WVI.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
//
28
//
29
//
30
//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
31
//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
32
33
#include "PhysListEmStandard_WVI.hh"
34
#include "DetectorConstruction.hh"
35
36
#include "
G4ParticleDefinition.hh
"
37
#include "
G4ProcessManager.hh
"
38
39
#include "
G4ComptonScattering.hh
"
40
#include "
MyKleinNishinaCompton.hh
"
41
#include "
G4GammaConversion.hh
"
42
#include "
G4PhotoElectricEffect.hh
"
43
44
#include "
G4eMultipleScattering.hh
"
45
#include "
G4MuMultipleScattering.hh
"
46
#include "
G4WentzelVIModel.hh
"
47
#include "
G4CoulombScattering.hh
"
48
#include "
G4eCoulombScatteringModel.hh
"
49
50
#include "
G4eIonisation.hh
"
51
#include "MyMollerBhabhaModel.hh"
52
#include "
G4eBremsstrahlung.hh
"
53
#include "
G4eplusAnnihilation.hh
"
54
55
#include "
G4hIonisation.hh
"
56
#include "
G4hMultipleScattering.hh
"
57
58
#include "
G4EmParameters.hh
"
59
#include "
G4MscStepLimitType.hh
"
60
61
#include "
G4BuilderType.hh
"
62
#include "
G4SystemOfUnits.hh
"
63
64
//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
65
66
PhysListEmStandard_WVI::PhysListEmStandard_WVI
(
const
G4String
&
name
,
67
DetectorConstruction
* det)
68
:
G4VPhysicsConstructor
(name), fDetector(det)
69
{
70
G4EmParameters
* param =
G4EmParameters::Instance
();
71
param->
SetDefaults
();
72
param->
SetVerbose
(1);
73
param->
SetMinEnergy
(10*
eV
);
74
param->
SetMaxEnergy
(10*
GeV
);
75
param->
SetNumberOfBinsPerDecade
(20);
76
param->
SetLowestElectronEnergy
(10*
eV
);
77
param->
SetStepFunction
(0.2, 100*
um
);
78
param->
SetBuildCSDARange
(
true
);
79
param->
SetMaxEnergyForCSDARange
(10*
GeV
);
80
param->
SetMscThetaLimit
(0.15);
81
SetPhysicsType
(
bElectromagnetic
);
82
}
83
84
//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
85
86
PhysListEmStandard_WVI::~PhysListEmStandard_WVI
()
87
{}
88
89
//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
90
91
void
PhysListEmStandard_WVI::ConstructProcess
()
92
{
93
// Add standard EM Processes
94
//
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
if
(particleName ==
"gamma"
) {
104
// gamma
105
106
G4ComptonScattering
* compton =
new
G4ComptonScattering
();
107
MyKleinNishinaCompton
* comptonModel =
108
new
MyKleinNishinaCompton
(
fDetector
);
109
comptonModel->
SetCSFactor
(1000.);
110
compton->
SetEmModel
(comptonModel );
111
112
pmanager->
AddDiscreteProcess
(
new
G4PhotoElectricEffect
);
113
pmanager->
AddDiscreteProcess
(compton);
114
pmanager->
AddDiscreteProcess
(
new
G4GammaConversion
);
115
116
}
else
if
(particleName ==
"e-"
) {
117
//electron
118
119
G4MuMultipleScattering
* eMsc =
new
G4MuMultipleScattering
();
120
G4WentzelVIModel
* wvi =
new
G4WentzelVIModel
();
121
eMsc->
SetEmModel
(wvi);
122
123
G4eIonisation
* eIoni =
new
G4eIonisation
();
124
eIoni->
SetEmModel
(
new
MyMollerBhabhaModel
());
125
126
G4CoulombScattering
* cs =
new
G4CoulombScattering
();
127
G4eCoulombScatteringModel
* single =
new
G4eCoulombScatteringModel
();
128
cs->
SetEmModel
(single, 1);
129
130
pmanager->
AddProcess
(eMsc, -1, 1,-1);
131
pmanager->
AddProcess
(eIoni,-1, 2, 1);
132
pmanager->
AddProcess
(cs, -1,-1, 2);
133
134
}
else
if
(particleName ==
"e+"
) {
135
//positron
136
137
G4MuMultipleScattering
* pMsc =
new
G4MuMultipleScattering
();
138
pMsc->
SetEmModel
(
new
G4WentzelVIModel
());
139
G4eIonisation
* pIoni =
new
G4eIonisation
();
140
pIoni->
SetEmModel
(
new
MyMollerBhabhaModel
);
141
142
pmanager->
AddProcess
(pMsc, -1, 1,-1);
143
pmanager->
AddProcess
(pIoni, -1, 2, 1);
144
pmanager->
AddProcess
(
new
G4eplusAnnihilation
, 0,-1, 2);
145
pmanager->
AddProcess
(
new
G4CoulombScattering
, -1,-1, 3);
146
147
}
else
if
( particleName ==
"proton"
) {
148
//proton
149
G4hMultipleScattering
* msc =
new
G4hMultipleScattering
();
150
msc->
AddEmModel
(1,
new
G4WentzelVIModel
());
151
pmanager->
AddProcess
(msc, -1, 1,-1);
152
pmanager->
AddProcess
(
new
G4hIonisation
, -1, 2, 1);
153
pmanager->
AddProcess
(
new
G4CoulombScattering
, -1,-1, 2);
154
}
155
}
156
}
157
158
//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
159
geant4
tree
geant4-10.6-release
examples
extended
medical
fanoCavity
src
PhysListEmStandard_WVI.cc
Built by
Jin Huang
. updated:
Wed Jun 29 2022 17:25:08
using
1.8.2 with
ECCE GitHub integration