| Geant4 Cross Reference |
1 // 1 //
2 // ******************************************* 2 // ********************************************************************
3 // * License and Disclaimer 3 // * License and Disclaimer *
4 // * 4 // * *
5 // * The Geant4 software is copyright of th 5 // * The Geant4 software is copyright of the Copyright Holders of *
6 // * the Geant4 Collaboration. It is provided 6 // * the Geant4 Collaboration. It is provided under the terms and *
7 // * conditions of the Geant4 Software License 7 // * conditions of the Geant4 Software License, included in the file *
8 // * LICENSE and available at http://cern.ch/ 8 // * LICENSE and available at http://cern.ch/geant4/license . These *
9 // * include a list of copyright holders. 9 // * include a list of copyright holders. *
10 // * 10 // * *
11 // * Neither the authors of this software syst 11 // * Neither the authors of this software system, nor their employing *
12 // * institutes,nor the agencies providing fin 12 // * institutes,nor the agencies providing financial support for this *
13 // * work make any representation or warran 13 // * work make any representation or warranty, express or implied, *
14 // * regarding this software system or assum 14 // * regarding this software system or assume any liability for its *
15 // * use. Please see the license in the file 15 // * use. Please see the license in the file LICENSE and URL above *
16 // * for the full disclaimer and the limitatio 16 // * for the full disclaimer and the limitation of liability. *
17 // * 17 // * *
18 // * This code implementation is the result 18 // * This code implementation is the result of the scientific and *
19 // * technical work of the GEANT4 collaboratio 19 // * technical work of the GEANT4 collaboration. *
20 // * By using, copying, modifying or distri 20 // * By using, copying, modifying or distributing the software (or *
21 // * any work based on the software) you ag 21 // * any work based on the software) you agree to acknowledge its *
22 // * use in resulting scientific publicati 22 // * use in resulting scientific publications, and indicate your *
23 // * acceptance of all terms of the Geant4 Sof 23 // * acceptance of all terms of the Geant4 Software license. *
24 // ******************************************* 24 // ********************************************************************
25 // 25 //
26 /// \file hadronic/Hadr01/src/PhysicsList.cc 26 /// \file hadronic/Hadr01/src/PhysicsList.cc
27 /// \brief Implementation of the PhysicsList c 27 /// \brief Implementation of the PhysicsList class
28 // 28 //
29 // 29 //
30 // 30 //
31 ////////////////////////////////////////////// 31 /////////////////////////////////////////////////////////////////////////
32 // 32 //
33 // PhysicsList 33 // PhysicsList
34 // 34 //
35 // Created: 31.04.2006 V.Ivanchenko 35 // Created: 31.04.2006 V.Ivanchenko
36 // 36 //
37 // Modified: 37 // Modified:
38 // 04.06.2006 Adoptation of Hadr01 (V.Ivanchen 38 // 04.06.2006 Adoptation of Hadr01 (V.Ivanchenko)
39 // 26.04.2007 Physics according to 8.3 Physics 39 // 26.04.2007 Physics according to 8.3 Physics List (V.Ivanchenko)
40 // 16.10.2012 Renamed used classes (A.Ribon) 40 // 16.10.2012 Renamed used classes (A.Ribon)
41 // 41 //
42 ////////////////////////////////////////////// 42 ////////////////////////////////////////////////////////////////////////
43 // << 43 //
44 44
45 #include "PhysicsList.hh" 45 #include "PhysicsList.hh"
46 <<
47 #include "PhysicsListMessenger.hh" 46 #include "PhysicsListMessenger.hh"
48 47
49 #include "G4DecayPhysics.hh" 48 #include "G4DecayPhysics.hh"
50 #include "G4Electron.hh" <<
51 #include "G4EmExtraPhysics.hh" <<
52 #include "G4EmLivermorePhysics.hh" <<
53 #include "G4EmParameters.hh" <<
54 #include "G4EmPenelopePhysics.hh" <<
55 #include "G4EmStandardPhysics.hh" 49 #include "G4EmStandardPhysics.hh"
56 #include "G4EmStandardPhysicsGS.hh" <<
57 #include "G4EmStandardPhysics_option1.hh" 50 #include "G4EmStandardPhysics_option1.hh"
58 #include "G4EmStandardPhysics_option2.hh" 51 #include "G4EmStandardPhysics_option2.hh"
59 #include "G4EmStandardPhysics_option3.hh" 52 #include "G4EmStandardPhysics_option3.hh"
60 #include "G4EmStandardPhysics_option4.hh" 53 #include "G4EmStandardPhysics_option4.hh"
61 #include "G4Gamma.hh" << 54 #include "G4EmStandardPhysicsGS.hh"
>> 55 #include "G4EmLivermorePhysics.hh"
>> 56 #include "G4EmPenelopePhysics.hh"
62 #include "G4HadronElasticPhysics.hh" 57 #include "G4HadronElasticPhysics.hh"
63 #include "G4HadronElasticPhysicsHP.hh" <<
64 #include "G4HadronElasticPhysicsXS.hh" 58 #include "G4HadronElasticPhysicsXS.hh"
>> 59 #include "G4HadronElasticPhysicsHP.hh"
65 #include "G4HadronHElasticPhysics.hh" 60 #include "G4HadronHElasticPhysics.hh"
66 #include "G4HadronInelasticQBBC.hh" << 61 #include "G4NeutronTrackingCut.hh"
>> 62 #include "G4NeutronCrossSectionXS.hh"
>> 63 #include "G4StoppingPhysics.hh"
>> 64 #include "G4IonBinaryCascadePhysics.hh"
>> 65 #include "G4IonPhysics.hh"
>> 66 #include "G4IonPhysicsXS.hh"
>> 67 #include "G4IonElasticPhysics.hh"
>> 68 #include "G4EmExtraPhysics.hh"
>> 69 #include "G4EmParameters.hh"
>> 70 #include "G4PhysListFactoryMessenger.hh"
>> 71
67 #include "G4HadronPhysicsFTFP_BERT.hh" 72 #include "G4HadronPhysicsFTFP_BERT.hh"
68 #include "G4HadronPhysicsFTFP_BERT_HP.hh" 73 #include "G4HadronPhysicsFTFP_BERT_HP.hh"
69 #include "G4HadronPhysicsFTFP_BERT_TRV.hh" 74 #include "G4HadronPhysicsFTFP_BERT_TRV.hh"
70 #include "G4HadronPhysicsFTF_BIC.hh" 75 #include "G4HadronPhysicsFTF_BIC.hh"
>> 76 #include "G4HadronInelasticQBBC.hh"
71 #include "G4HadronPhysicsQGSP_BERT.hh" 77 #include "G4HadronPhysicsQGSP_BERT.hh"
72 #include "G4HadronPhysicsQGSP_BERT_HP.hh" 78 #include "G4HadronPhysicsQGSP_BERT_HP.hh"
73 #include "G4HadronPhysicsQGSP_BIC.hh" 79 #include "G4HadronPhysicsQGSP_BIC.hh"
74 #include "G4HadronPhysicsQGSP_BIC_HP.hh" 80 #include "G4HadronPhysicsQGSP_BIC_HP.hh"
75 #include "G4HadronPhysicsQGSP_FTFP_BERT.hh" 81 #include "G4HadronPhysicsQGSP_FTFP_BERT.hh"
76 #include "G4HadronPhysicsQGS_BIC.hh" 82 #include "G4HadronPhysicsQGS_BIC.hh"
77 #include "G4IonBinaryCascadePhysics.hh" << 83 #include "G4RadioactiveDecayPhysics.hh"
78 #include "G4IonElasticPhysics.hh" << 84
79 #include "G4IonPhysics.hh" << 85 #include "G4ProcessManager.hh"
80 #include "G4IonPhysicsXS.hh" <<
81 #include "G4NeutronCrossSectionXS.hh" <<
82 #include "G4NeutronTrackingCut.hh" <<
83 #include "G4ParticleTable.hh" <<
84 #include "G4ParticleTypes.hh" 86 #include "G4ParticleTypes.hh"
85 #include "G4PhysListFactoryMessenger.hh" << 87 #include "G4ParticleTable.hh"
>> 88 #include "G4Gamma.hh"
>> 89 #include "G4Electron.hh"
86 #include "G4Positron.hh" 90 #include "G4Positron.hh"
87 #include "G4ProcessManager.hh" <<
88 #include "G4Proton.hh" 91 #include "G4Proton.hh"
89 #include "G4RadioactiveDecayPhysics.hh" << 92
90 #include "G4StoppingPhysics.hh" <<
91 #include "G4SystemOfUnits.hh" 93 #include "G4SystemOfUnits.hh"
92 94
93 //....oooOO0OOooo........oooOO0OOooo........oo 95 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo.....
94 96
95 PhysicsList::PhysicsList() : G4VModularPhysics 97 PhysicsList::PhysicsList() : G4VModularPhysicsList()
96 { 98 {
97 SetDefaultCutValue(0.7 * CLHEP::mm); << 99 SetDefaultCutValue(0.7*CLHEP::mm);
>> 100
98 verboseLevel = 1; 101 verboseLevel = 1;
99 102
100 // Messengers <<
101 fMessenger = new PhysicsListMessenger(this); 103 fMessenger = new PhysicsListMessenger(this);
102 fFactMessenger = new G4PhysListFactoryMessen 104 fFactMessenger = new G4PhysListFactoryMessenger(this);
103 105
104 // Particles 106 // Particles
105 fParticleList = new G4DecayPhysics(verboseLe 107 fParticleList = new G4DecayPhysics(verboseLevel);
106 108
107 // EM physics 109 // EM physics
108 fEmPhysicsList = new G4EmStandardPhysics(ver 110 fEmPhysicsList = new G4EmStandardPhysics(verboseLevel);
109 } 111 }
110 112
111 //....oooOO0OOooo........oooOO0OOooo........oo 113 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo.....
112 114
113 PhysicsList::~PhysicsList() 115 PhysicsList::~PhysicsList()
114 { 116 {
115 delete fMessenger; 117 delete fMessenger;
116 delete fParticleList; 118 delete fParticleList;
117 delete fEmPhysicsList; 119 delete fEmPhysicsList;
118 for (std::size_t i = 0; i < fHadronPhys.size << 120 for(size_t i=0; i<fHadronPhys.size(); i++) {
119 delete fHadronPhys[i]; 121 delete fHadronPhys[i];
120 } 122 }
121 } 123 }
122 124
123 //....oooOO0OOooo........oooOO0OOooo........oo 125 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo.....
124 126
125 void PhysicsList::ConstructParticle() 127 void PhysicsList::ConstructParticle()
126 { 128 {
127 fParticleList->ConstructParticle(); 129 fParticleList->ConstructParticle();
128 } 130 }
129 131
130 //....oooOO0OOooo........oooOO0OOooo........oo 132 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo.....
131 133
132 void PhysicsList::ConstructProcess() 134 void PhysicsList::ConstructProcess()
133 { 135 {
134 AddTransportation(); 136 AddTransportation();
135 fEmPhysicsList->ConstructProcess(); 137 fEmPhysicsList->ConstructProcess();
136 fParticleList->ConstructProcess(); 138 fParticleList->ConstructProcess();
137 for (std::size_t i = 0; i < fHadronPhys.size << 139 for(size_t i=0; i<fHadronPhys.size(); i++) {
138 fHadronPhys[i]->ConstructProcess(); 140 fHadronPhys[i]->ConstructProcess();
139 } 141 }
140 } 142 }
141 143
142 //....oooOO0OOooo........oooOO0OOooo........oo 144 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo.....
143 145
144 void PhysicsList::AddPhysicsList(const G4Strin 146 void PhysicsList::AddPhysicsList(const G4String& name)
145 { 147 {
146 if (verboseLevel > 0) { << 148 if (verboseLevel>0) {
147 G4cout << "PhysicsList::AddPhysicsList: <" 149 G4cout << "PhysicsList::AddPhysicsList: <" << name << ">" << G4endl;
148 } 150 }
149 if (name == "emstandard_opt0") { 151 if (name == "emstandard_opt0") {
>> 152
150 delete fEmPhysicsList; 153 delete fEmPhysicsList;
151 fEmPhysicsList = new G4EmStandardPhysics(v 154 fEmPhysicsList = new G4EmStandardPhysics(verboseLevel);
152 } << 155
153 else if (name == "emstandard_opt1") { << 156 } else if (name == "emstandard_opt1") {
>> 157
154 delete fEmPhysicsList; 158 delete fEmPhysicsList;
155 fEmPhysicsList = new G4EmStandardPhysics_o 159 fEmPhysicsList = new G4EmStandardPhysics_option1(verboseLevel);
156 } << 160
157 else if (name == "emstandard_opt2") { << 161 } else if (name == "emstandard_opt2") {
>> 162
158 delete fEmPhysicsList; 163 delete fEmPhysicsList;
159 fEmPhysicsList = new G4EmStandardPhysics_o 164 fEmPhysicsList = new G4EmStandardPhysics_option2(verboseLevel);
160 } << 165
161 else if (name == "emstandard_opt3") { << 166 } else if (name == "emstandard_opt3") {
>> 167
162 delete fEmPhysicsList; 168 delete fEmPhysicsList;
163 fEmPhysicsList = new G4EmStandardPhysics_o 169 fEmPhysicsList = new G4EmStandardPhysics_option3(verboseLevel);
164 } << 170
165 else if (name == "emstandard_opt4") { << 171 } else if (name == "emstandard_opt4") {
>> 172
166 delete fEmPhysicsList; 173 delete fEmPhysicsList;
167 fEmPhysicsList = new G4EmStandardPhysics_o 174 fEmPhysicsList = new G4EmStandardPhysics_option4(verboseLevel);
168 } << 175
169 else if (name == "emstandardGS") { << 176 } else if (name == "emstandardGS") {
>> 177
170 delete fEmPhysicsList; 178 delete fEmPhysicsList;
171 fEmPhysicsList = new G4EmStandardPhysicsGS 179 fEmPhysicsList = new G4EmStandardPhysicsGS(verboseLevel);
172 } << 180
173 else if (name == "FTFP_BERT_EMV") { << 181 } else if (name == "FTFP_BERT_EMV") {
>> 182
174 AddPhysicsList("FTFP_BERT"); 183 AddPhysicsList("FTFP_BERT");
175 AddPhysicsList("emstandard_opt1"); 184 AddPhysicsList("emstandard_opt1");
176 } << 185
177 else if (name == "FTFP_BERT_EMX") { << 186 } else if (name == "FTFP_BERT_EMX") {
>> 187
178 AddPhysicsList("FTFP_BERT"); 188 AddPhysicsList("FTFP_BERT");
179 AddPhysicsList("emstandard_opt2"); 189 AddPhysicsList("emstandard_opt2");
180 } << 190
181 else if (name == "FTFP_BERT_EMY") { << 191 } else if (name == "FTFP_BERT_EMY") {
>> 192
182 AddPhysicsList("FTFP_BERT"); 193 AddPhysicsList("FTFP_BERT");
183 AddPhysicsList("emstandard_opt3"); 194 AddPhysicsList("emstandard_opt3");
184 } << 195
185 else if (name == "FTFP_BERT_EMZ") { << 196 } else if (name == "FTFP_BERT_EMZ") {
>> 197
186 AddPhysicsList("FTFP_BERT"); 198 AddPhysicsList("FTFP_BERT");
187 AddPhysicsList("emstandard_opt4"); 199 AddPhysicsList("emstandard_opt4");
188 } << 200
189 else if (name == "FTFP_BERT") { << 201 } else if (name == "FTFP_BERT") {
>> 202
190 SetBuilderList0(false); 203 SetBuilderList0(false);
191 fHadronPhys.push_back(new G4HadronPhysicsF << 204 fHadronPhys.push_back( new G4HadronPhysicsFTFP_BERT(verboseLevel));
192 } << 205
193 else if (name == "FTFP_BERT_TRV") { << 206 } else if (name == "FTFP_BERT_TRV") {
>> 207
194 AddPhysicsList("emstandardGS"); 208 AddPhysicsList("emstandardGS");
195 SetBuilderList1(false); 209 SetBuilderList1(false);
196 fHadronPhys.push_back(new G4HadronPhysicsF << 210 fHadronPhys.push_back( new G4HadronPhysicsFTFP_BERT_TRV(verboseLevel));
197 } << 211
198 else if (name == "FTF_BIC") { << 212 } else if (name == "FTF_BIC") {
>> 213
199 SetBuilderList0(false); 214 SetBuilderList0(false);
200 fHadronPhys.push_back(new G4HadronPhysicsF << 215 fHadronPhys.push_back( new G4HadronPhysicsFTF_BIC(verboseLevel));
201 } << 216
202 else if (name == "QBBC") { << 217 } else if (name == "QBBC") {
>> 218
203 AddPhysicsList("emstandard_opt0"); 219 AddPhysicsList("emstandard_opt0");
204 SetBuilderList2(); 220 SetBuilderList2();
205 fHadronPhys.push_back(new G4HadronInelasti << 221 fHadronPhys.push_back( new G4HadronInelasticQBBC(verboseLevel));
206 } << 222
207 else if (name == "QGSP_BERT") { << 223 } else if (name == "QGSP_BERT") {
>> 224
208 SetBuilderList0(false); 225 SetBuilderList0(false);
209 fHadronPhys.push_back(new G4HadronPhysicsQ << 226 fHadronPhys.push_back( new G4HadronPhysicsQGSP_BERT(verboseLevel));
210 } << 227
211 else if (name == "QGSP_FTFP_BERT") { << 228 } else if (name == "QGSP_FTFP_BERT") {
>> 229
212 SetBuilderList0(false); 230 SetBuilderList0(false);
213 fHadronPhys.push_back(new G4HadronPhysicsQ << 231 fHadronPhys.push_back( new G4HadronPhysicsQGSP_FTFP_BERT(verboseLevel));
214 } << 232
215 else if (name == "QGSP_FTFP_BERT_EMV") { << 233 } else if (name == "QGSP_FTFP_BERT_EMV") {
>> 234
216 AddPhysicsList("QGSP_FTFP_BERT"); 235 AddPhysicsList("QGSP_FTFP_BERT");
217 AddPhysicsList("emstandard_opt1"); 236 AddPhysicsList("emstandard_opt1");
218 } << 237
219 else if (name == "QGSP_BERT_EMV") { << 238 } else if (name == "QGSP_BERT_EMV") {
>> 239
220 AddPhysicsList("QGSP_BERT"); 240 AddPhysicsList("QGSP_BERT");
221 AddPhysicsList("emstandard_opt1"); 241 AddPhysicsList("emstandard_opt1");
222 } << 242
223 else if (name == "QGSP_BERT_EMX") { << 243 } else if (name == "QGSP_BERT_EMX") {
>> 244
224 AddPhysicsList("QGSP_BERT"); 245 AddPhysicsList("QGSP_BERT");
225 AddPhysicsList("emstandard_opt2"); 246 AddPhysicsList("emstandard_opt2");
226 } << 247
227 else if (name == "QGSP_BERT_HP") { << 248 } else if (name == "QGSP_BERT_HP") {
>> 249
228 SetBuilderList0(true); 250 SetBuilderList0(true);
229 fHadronPhys.push_back(new G4HadronPhysicsQ << 251 fHadronPhys.push_back( new G4HadronPhysicsQGSP_BERT_HP(verboseLevel));
230 AddPhysicsList("RadioactiveDecay"); 252 AddPhysicsList("RadioactiveDecay");
231 } << 253
232 else if (name == "QGSP_BIC") { << 254 } else if (name == "QGSP_BIC") {
>> 255
233 SetBuilderList0(false); 256 SetBuilderList0(false);
234 fHadronPhys.push_back(new G4HadronPhysicsQ << 257 fHadronPhys.push_back( new G4HadronPhysicsQGSP_BIC(verboseLevel));
235 fHadronPhys.push_back(new G4IonElasticPhys << 258 fHadronPhys.push_back( new G4IonElasticPhysics(verboseLevel));
236 } << 259
237 else if (name == "QGSP_BIC_EMY") { << 260 } else if (name == "QGSP_BIC_EMY") {
>> 261
238 AddPhysicsList("QGSP_BIC"); 262 AddPhysicsList("QGSP_BIC");
239 AddPhysicsList("emstandard_opt3"); 263 AddPhysicsList("emstandard_opt3");
240 } << 264
241 else if (name == "QGS_BIC") { << 265 } else if (name == "QGS_BIC") {
>> 266
242 SetBuilderList0(false); 267 SetBuilderList0(false);
243 fHadronPhys.push_back(new G4HadronPhysicsQ << 268 fHadronPhys.push_back( new G4HadronPhysicsQGS_BIC(verboseLevel));
244 } << 269
245 else if (name == "QGSP_BIC_HP") { << 270 } else if (name == "QGSP_BIC_HP") {
>> 271
246 SetBuilderList0(true); 272 SetBuilderList0(true);
247 fHadronPhys.push_back(new G4HadronPhysicsQ << 273 fHadronPhys.push_back( new G4HadronPhysicsQGSP_BIC_HP(verboseLevel));
248 AddPhysicsList("RadioactiveDecay"); 274 AddPhysicsList("RadioactiveDecay");
249 } << 275
250 else if (name == "RadioactiveDecay") { << 276 } else if (name == "RadioactiveDecay") {
251 fHadronPhys.push_back(new G4RadioactiveDec << 277
252 } << 278 fHadronPhys.push_back( new G4RadioactiveDecayPhysics(verboseLevel));
253 else { << 279
>> 280 } else {
>> 281
254 G4cout << "PhysicsList::AddPhysicsList: <" 282 G4cout << "PhysicsList::AddPhysicsList: <" << name << ">"
255 << " is not defined" << G4endl; << 283 << " is not defined"
>> 284 << G4endl;
256 } 285 }
257 } 286 }
258 287
259 //....oooOO0OOooo........oooOO0OOooo........oo 288 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo.....
260 289
261 void PhysicsList::SetBuilderList0(G4bool flagH 290 void PhysicsList::SetBuilderList0(G4bool flagHP)
262 { 291 {
263 fHadronPhys.push_back(new G4EmExtraPhysics(v << 292 fHadronPhys.push_back( new G4EmExtraPhysics(verboseLevel));
264 if (flagHP) { << 293 if(flagHP) {
265 fHadronPhys.push_back(new G4HadronElasticP << 294 fHadronPhys.push_back( new G4HadronElasticPhysicsHP(verboseLevel) );
266 } << 295 } else {
267 else { << 296 fHadronPhys.push_back( new G4HadronElasticPhysics(verboseLevel) );
268 fHadronPhys.push_back(new G4HadronElasticP << 297 }
269 } << 298 fHadronPhys.push_back( new G4StoppingPhysics(verboseLevel));
270 fHadronPhys.push_back(new G4StoppingPhysics( << 299 fHadronPhys.push_back( new G4IonPhysics(verboseLevel));
271 fHadronPhys.push_back(new G4IonPhysics(verbo << 300 if(!flagHP) {
272 if (!flagHP) { << 301 fHadronPhys.push_back( new G4NeutronTrackingCut(verboseLevel));
273 fHadronPhys.push_back(new G4NeutronTrackin <<
274 } 302 }
275 } 303 }
276 304
277 //....oooOO0OOooo........oooOO0OOooo........oo 305 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo.....
278 306
279 void PhysicsList::SetBuilderList1(G4bool flagH 307 void PhysicsList::SetBuilderList1(G4bool flagHP)
280 { 308 {
281 fHadronPhys.push_back(new G4EmExtraPhysics(v << 309 fHadronPhys.push_back( new G4EmExtraPhysics(verboseLevel));
282 if (flagHP) { << 310 if(flagHP) {
283 fHadronPhys.push_back(new G4HadronElasticP << 311 fHadronPhys.push_back( new G4HadronElasticPhysicsHP(verboseLevel) );
284 } << 312 } else {
285 else { << 313 fHadronPhys.push_back( new G4HadronHElasticPhysics(verboseLevel) );
286 fHadronPhys.push_back(new G4HadronHElastic << 314 }
287 } << 315 fHadronPhys.push_back( new G4StoppingPhysics(verboseLevel));
288 fHadronPhys.push_back(new G4StoppingPhysics( << 316 fHadronPhys.push_back( new G4IonPhysics(verboseLevel));
289 fHadronPhys.push_back(new G4IonPhysics(verbo << 317 if(!flagHP) {
290 if (!flagHP) { << 318 fHadronPhys.push_back( new G4NeutronTrackingCut(verboseLevel));
291 fHadronPhys.push_back(new G4NeutronTrackin <<
292 } 319 }
293 } 320 }
294 321
295 //....oooOO0OOooo........oooOO0OOooo........oo 322 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo.....
296 323
297 void PhysicsList::SetBuilderList2() 324 void PhysicsList::SetBuilderList2()
298 { 325 {
299 fHadronPhys.push_back(new G4EmExtraPhysics(v << 326 fHadronPhys.push_back( new G4EmExtraPhysics(verboseLevel));
300 fHadronPhys.push_back(new G4HadronElasticPhy << 327 fHadronPhys.push_back( new G4HadronElasticPhysicsXS(verboseLevel) );
301 fHadronPhys.push_back(new G4StoppingPhysics( << 328 fHadronPhys.push_back( new G4StoppingPhysics(verboseLevel));
302 fHadronPhys.push_back(new G4IonPhysicsXS(ver << 329 fHadronPhys.push_back( new G4IonPhysicsXS(verboseLevel));
303 fHadronPhys.push_back(new G4IonElasticPhysic << 330 fHadronPhys.push_back( new G4IonElasticPhysics(verboseLevel));
304 fHadronPhys.push_back(new G4NeutronTrackingC << 331 fHadronPhys.push_back( new G4NeutronTrackingCut(verboseLevel));
305 } 332 }
306 333
307 //....oooOO0OOooo........oooOO0OOooo........oo 334 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo.....
308 335
309 void PhysicsList::List() 336 void PhysicsList::List()
310 { 337 {
311 G4cout << "### PhysicsLists available: FTFP_ 338 G4cout << "### PhysicsLists available: FTFP_BERT FTFP_BERT_EMV "
312 << "FTFP_BERT_EMX FTFP_BERT_EMZ FTFP_ << 339 << "FTFP_BERT_EMX FTFP_BERT_EMZ FTFP_BERT_TRV"
>> 340 << G4endl;
313 G4cout << " FTF_B 341 G4cout << " FTF_BIC QBBC QGSP_BERT "
314 << "QGSP_BERT_EMV QGSP_BERT_EMX" << G << 342 << "QGSP_BERT_EMV QGSP_BERT_EMX"
>> 343 << G4endl;
315 G4cout << " QGSP_ 344 G4cout << " QGSP_BERT_HP QGSP_FTFP_BERT "
316 << "QGSP_FTFP_BERT_EMV" << G4endl; << 345 << "QGSP_FTFP_BERT_EMV"
>> 346 << G4endl;
317 G4cout << " QGS_B 347 G4cout << " QGS_BIC QGSP_BIC QGSP_BIC_EMY "
318 << "QGSP_BIC_HP" << G4endl; << 348 << "QGSP_BIC_HP"
>> 349 << G4endl;
319 } 350 }
320 351
321 //....oooOO0OOooo........oooOO0OOooo........oo 352 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
>> 353
322 354