| Geant4 Cross Reference |
1 // 1
2 // *******************************************
3 // * License and Disclaimer
4 // *
5 // * The Geant4 software is copyright of th
6 // * the Geant4 Collaboration. It is provided
7 // * conditions of the Geant4 Software License
8 // * LICENSE and available at http://cern.ch/
9 // * include a list of copyright holders.
10 // *
11 // * Neither the authors of this software syst
12 // * institutes,nor the agencies providing fin
13 // * work make any representation or warran
14 // * regarding this software system or assum
15 // * use. Please see the license in the file
16 // * for the full disclaimer and the limitatio
17 // *
18 // * This code implementation is the result
19 // * technical work of the GEANT4 collaboratio
20 // * By using, copying, modifying or distri
21 // * any work based on the software) you ag
22 // * use in resulting scientific publicati
23 // * acceptance of all terms of the Geant4 Sof
24 // *******************************************
25 //
26 /// \file electromagnetic/TestEm9/src/PhysicsL
27 /// \brief Implementation of the PhysicsList c
28 //
29 //
30 //--------------------------------------------
31 //
32 // ClassName: PhysicsList
33 //
34 // Author: V.Ivanchenko 14.10.2002
35 //
36 // Modified:
37 // 17.11.06 Use components from physics_lists
38 // 24.10.12 Migrated to new stopping and ion p
39 //
40 //--------------------------------------------
41 //
42
43 //....oooOO0OOooo........oooOO0OOooo........oo
44 //....oooOO0OOooo........oooOO0OOooo........oo
45
46 #include "PhysicsList.hh"
47
48 #include "PhysListEmStandard.hh"
49 #include "PhysicsListMessenger.hh"
50 #include "StepMax.hh"
51
52 #include "G4DecayPhysics.hh"
53 #include "G4Electron.hh"
54 #include "G4EmExtraPhysics.hh"
55 #include "G4EmLivermorePhysics.hh"
56 #include "G4EmLowEPPhysics.hh"
57 #include "G4EmPenelopePhysics.hh"
58 #include "G4EmStandardPhysics.hh"
59 #include "G4EmStandardPhysicsGS.hh"
60 #include "G4EmStandardPhysicsSS.hh"
61 #include "G4EmStandardPhysicsWVI.hh"
62 #include "G4EmStandardPhysics_option1.hh"
63 #include "G4EmStandardPhysics_option2.hh"
64 #include "G4EmStandardPhysics_option3.hh"
65 #include "G4EmStandardPhysics_option4.hh"
66 #include "G4Gamma.hh"
67 #include "G4HadronElasticPhysics.hh"
68 #include "G4HadronInelasticQBBC.hh"
69 #include "G4IonPhysics.hh"
70 #include "G4LossTableManager.hh"
71 #include "G4ParticleTable.hh"
72 #include "G4ParticleTypes.hh"
73 #include "G4Positron.hh"
74 #include "G4ProcessManager.hh"
75 #include "G4Proton.hh"
76 #include "G4RegionStore.hh"
77 #include "G4StoppingPhysics.hh"
78 #include "G4SystemOfUnits.hh"
79 #include "G4UnitsTable.hh"
80
81 //....oooOO0OOooo........oooOO0OOooo........oo
82
83 PhysicsList::PhysicsList()
84 : G4VModularPhysicsList(),
85 fEmPhysicsList(nullptr),
86 fDecayPhysicsList(nullptr),
87 fStepMaxProcess(nullptr),
88 fMessenger(nullptr)
89 {
90 G4LossTableManager::Instance();
91 SetDefaultCutValue(1 * mm);
92
93 fMessenger = new PhysicsListMessenger(this);
94 fStepMaxProcess = new StepMax();
95
96 // Initilise flags
97
98 SetVerboseLevel(1);
99
100 fHelIsRegisted = false;
101 fBicIsRegisted = false;
102 fGnucIsRegisted = false;
103 fStopIsRegisted = false;
104
105 // EM physics
106 fEmName = G4String("emstandard");
107 fEmPhysicsList = new G4EmStandardPhysics();
108
109 // Decay Physics is always defined
110 fDecayPhysicsList = new G4DecayPhysics();
111 }
112
113 //....oooOO0OOooo........oooOO0OOooo........oo
114
115 PhysicsList::~PhysicsList()
116 {
117 delete fMessenger;
118 delete fDecayPhysicsList;
119 delete fEmPhysicsList;
120 delete fStepMaxProcess;
121 for (size_t i = 0; i < fHadronPhys.size(); i
122 delete fHadronPhys[i];
123 }
124 }
125
126 //....oooOO0OOooo........oooOO0OOooo........oo
127
128 void PhysicsList::ConstructParticle()
129 {
130 fDecayPhysicsList->ConstructParticle();
131 }
132
133 //....oooOO0OOooo........oooOO0OOooo........oo
134
135 void PhysicsList::ConstructProcess()
136 {
137 AddTransportation();
138 fEmPhysicsList->ConstructProcess();
139 fDecayPhysicsList->ConstructProcess();
140 for (size_t i = 0; i < fHadronPhys.size(); +
141 fHadronPhys[i]->ConstructProcess();
142 }
143 AddStepMax();
144 }
145
146 //....oooOO0OOooo........oooOO0OOooo........oo
147
148 void PhysicsList::AddPhysicsList(const G4Strin
149 {
150 if (verboseLevel > 1) G4cout << "PhysicsList
151
152 if (name == fEmName) return;
153
154 if (name == "emstandard") {
155 fEmName = name;
156 delete fEmPhysicsList;
157 fEmPhysicsList = new G4EmStandardPhysics(v
158 if (verboseLevel > 0) G4cout << "PhysicsLi
159 }
160 else if (name == "emstandard_opt1") {
161 fEmName = name;
162 delete fEmPhysicsList;
163 fEmPhysicsList = new G4EmStandardPhysics_o
164 if (verboseLevel > 0) G4cout << "PhysicsLi
165 }
166 else if (name == "emstandard_opt2") {
167 fEmName = name;
168 delete fEmPhysicsList;
169 fEmPhysicsList = new G4EmStandardPhysics_o
170 if (verboseLevel > 0) G4cout << "PhysicsLi
171 }
172 else if (name == "emstandard_opt3") {
173 fEmName = name;
174 delete fEmPhysicsList;
175 fEmPhysicsList = new G4EmStandardPhysics_o
176 if (verboseLevel > 0) G4cout << "PhysicsLi
177 }
178 else if (name == "emstandard_opt4") {
179 fEmName = name;
180 delete fEmPhysicsList;
181 fEmPhysicsList = new G4EmStandardPhysics_o
182 if (verboseLevel > 0) G4cout << "PhysicsLi
183 }
184 else if (name == "local") {
185 fEmName = name;
186 delete fEmPhysicsList;
187 fEmPhysicsList = new PhysListEmStandard(ve
188 if (verboseLevel > 0) G4cout << "PhysicsLi
189 }
190 else if (name == "emlivermore") {
191 fEmName = name;
192 delete fEmPhysicsList;
193 fEmPhysicsList = new G4EmLivermorePhysics(
194 }
195 else if (name == "empenelope") {
196 fEmName = name;
197 delete fEmPhysicsList;
198 fEmPhysicsList = new G4EmPenelopePhysics(v
199 }
200 else if (name == "emlowenergy") {
201 fEmName = name;
202 delete fEmPhysicsList;
203 fEmPhysicsList = new G4EmLowEPPhysics(verb
204 }
205 else if (name == "emstandardGS") {
206 fEmName = name;
207 delete fEmPhysicsList;
208 fEmPhysicsList = new G4EmStandardPhysicsGS
209 }
210 else if (name == "emstandardSS") {
211 fEmName = name;
212 delete fEmPhysicsList;
213 fEmPhysicsList = new G4EmStandardPhysicsSS
214 }
215 else if (name == "emstandardWVI") {
216 fEmName = name;
217 delete fEmPhysicsList;
218 fEmPhysicsList = new G4EmStandardPhysicsWV
219 }
220 else if (name == "elastic" && !fHelIsRegiste
221 fHadronPhys.push_back(new G4HadronElasticP
222 fHelIsRegisted = true;
223 if (verboseLevel > 0) G4cout << "PhysicsLi
224 }
225 else if (name == "binary" && !fBicIsRegisted
226 fHadronPhys.push_back(new G4HadronInelasti
227 fHadronPhys.push_back(new G4IonPhysics());
228 fBicIsRegisted = true;
229 if (verboseLevel > 0)
230 G4cout << "PhysicsList::Add hadron inela
231 }
232 else if (name == "gamma_nuc" && !fGnucIsRegi
233 fHadronPhys.push_back(new G4EmExtraPhysics
234 fGnucIsRegisted = true;
235 if (verboseLevel > 0) G4cout << "PhysicsLi
236 }
237 else if (name == "stopping" && !fStopIsRegis
238 fHadronPhys.push_back(new G4StoppingPhysic
239 fStopIsRegisted = true;
240 if (verboseLevel > 0) G4cout << "PhysicsLi
241 }
242 else {
243 G4cout << "PhysicsList::AddPhysicsList: <"
244 << " is not defined" << G4endl;
245 }
246 }
247
248 //....oooOO0OOooo........oooOO0OOooo........oo
249
250 void PhysicsList::AddStepMax()
251 {
252 // Step limitation seen as a process
253
254 auto particleIterator = GetParticleIterator(
255 particleIterator->reset();
256 while ((*particleIterator)()) {
257 G4ParticleDefinition* particle = particleI
258 G4ProcessManager* pmanager = particle->Get
259
260 if (fStepMaxProcess->IsApplicable(*particl
261 pmanager->AddDiscreteProcess(fStepMaxPro
262 }
263 }
264 }
265
266 //....oooOO0OOooo........oooOO0OOooo........oo
267