| 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 <<
27 /// \brief Implementation of the PhysicsList c <<
28 // <<
29 // 26 //
>> 27 // $Id: PhysicsList.cc,v 1.32 2009/11/25 19:56:36 vnivanch Exp $
>> 28 // GEANT4 tag $Name: geant4-09-03-patch-02 $
30 // 29 //
31 ////////////////////////////////////////////// 30 /////////////////////////////////////////////////////////////////////////
32 // 31 //
33 // PhysicsList 32 // PhysicsList
34 // 33 //
35 // Created: 31.04.2006 V.Ivanchenko 34 // Created: 31.04.2006 V.Ivanchenko
36 // 35 //
37 // Modified: 36 // Modified:
38 // 04.06.2006 Adoptation of Hadr01 (V.Ivanchen << 37 // 04.06.2006 Adoptation of hadr01 (V.Ivanchenko)
39 // 26.04.2007 Physics according to 8.3 Physics 38 // 26.04.2007 Physics according to 8.3 Physics List (V.Ivanchenko)
40 // 16.10.2012 Renamed used classes (A.Ribon) <<
41 // 39 //
42 ////////////////////////////////////////////// 40 ////////////////////////////////////////////////////////////////////////
43 // << 41 //
44 42
45 #include "PhysicsList.hh" 43 #include "PhysicsList.hh"
46 <<
47 #include "PhysicsListMessenger.hh" 44 #include "PhysicsListMessenger.hh"
48 45
49 #include "G4DecayPhysics.hh" 46 #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" 47 #include "G4EmStandardPhysics.hh"
56 #include "G4EmStandardPhysicsGS.hh" <<
57 #include "G4EmStandardPhysics_option1.hh" 48 #include "G4EmStandardPhysics_option1.hh"
58 #include "G4EmStandardPhysics_option2.hh" 49 #include "G4EmStandardPhysics_option2.hh"
59 #include "G4EmStandardPhysics_option3.hh" 50 #include "G4EmStandardPhysics_option3.hh"
60 #include "G4EmStandardPhysics_option4.hh" << 51 #include "G4EmLivermorePhysics.hh"
61 #include "G4Gamma.hh" << 52 #include "G4EmPenelopePhysics.hh"
62 #include "G4HadronElasticPhysics.hh" 53 #include "G4HadronElasticPhysics.hh"
63 #include "G4HadronElasticPhysicsHP.hh" << 54 #include "G4HadronDElasticPhysics.hh"
64 #include "G4HadronElasticPhysicsXS.hh" << 55 #include "G4HadronQElasticPhysics.hh"
65 #include "G4HadronHElasticPhysics.hh" 56 #include "G4HadronHElasticPhysics.hh"
66 #include "G4HadronInelasticQBBC.hh" << 57 #include "G4ChargeExchangePhysics.hh"
67 #include "G4HadronPhysicsFTFP_BERT.hh" << 58 #include "G4NeutronTrackingCut.hh"
68 #include "G4HadronPhysicsFTFP_BERT_HP.hh" << 59 #include "G4QStoppingPhysics.hh"
69 #include "G4HadronPhysicsFTFP_BERT_TRV.hh" << 60 #include "G4LHEPStoppingPhysics.hh"
70 #include "G4HadronPhysicsFTF_BIC.hh" <<
71 #include "G4HadronPhysicsQGSP_BERT.hh" <<
72 #include "G4HadronPhysicsQGSP_BERT_HP.hh" <<
73 #include "G4HadronPhysicsQGSP_BIC.hh" <<
74 #include "G4HadronPhysicsQGSP_BIC_HP.hh" <<
75 #include "G4HadronPhysicsQGSP_FTFP_BERT.hh" <<
76 #include "G4HadronPhysicsQGS_BIC.hh" <<
77 #include "G4IonBinaryCascadePhysics.hh" 61 #include "G4IonBinaryCascadePhysics.hh"
78 #include "G4IonElasticPhysics.hh" <<
79 #include "G4IonPhysics.hh" 62 #include "G4IonPhysics.hh"
80 #include "G4IonPhysicsXS.hh" << 63 #include "G4EmExtraPhysics.hh"
81 #include "G4NeutronCrossSectionXS.hh" << 64 #include "G4EmProcessOptions.hh"
82 #include "G4NeutronTrackingCut.hh" << 65
83 #include "G4ParticleTable.hh" << 66 #include "HadronPhysicsFTFP_BERT_TRV.hh"
>> 67 #include "HadronPhysicsFTFP.hh"
>> 68 #include "HadronPhysicsFTFP_BERT.hh"
>> 69 #include "HadronPhysicsFTF_BIC.hh"
>> 70 #include "HadronPhysicsLHEP.hh"
>> 71 #include "HadronPhysicsLHEP_EMV.hh"
>> 72 #include "G4HadronInelasticQBBC.hh"
>> 73 #include "HadronPhysicsQGSC.hh"
>> 74 #include "HadronPhysicsQGSC_BERT.hh"
>> 75 #include "HadronPhysicsQGSC_CHIPS.hh"
>> 76 #include "HadronPhysicsQGSC_QGSC.hh"
>> 77 #include "HadronPhysicsQGSP.hh"
>> 78 #include "HadronPhysicsQGSP_BERT.hh"
>> 79 #include "HadronPhysicsQGSP_BERT_HP.hh"
>> 80 #include "HadronPhysicsQGSP_BIC.hh"
>> 81 #include "HadronPhysicsQGSP_BIC_HP.hh"
>> 82 #include "HadronPhysicsQGSP_FTFP_BERT.hh"
>> 83 #include "HadronPhysicsQGS_BIC.hh"
>> 84
>> 85 #include "G4HadronInelasticQLHEP.hh"
>> 86 #include "G4IonPhysics.hh"
>> 87
>> 88 #include "G4LossTableManager.hh"
>> 89
>> 90 #include "G4ProcessManager.hh"
84 #include "G4ParticleTypes.hh" 91 #include "G4ParticleTypes.hh"
85 #include "G4PhysListFactoryMessenger.hh" << 92 #include "G4ParticleTable.hh"
>> 93 #include "G4Gamma.hh"
>> 94 #include "G4Electron.hh"
86 #include "G4Positron.hh" 95 #include "G4Positron.hh"
87 #include "G4ProcessManager.hh" <<
88 #include "G4Proton.hh" 96 #include "G4Proton.hh"
89 #include "G4RadioactiveDecayPhysics.hh" <<
90 #include "G4StoppingPhysics.hh" <<
91 #include "G4SystemOfUnits.hh" <<
92 97
93 //....oooOO0OOooo........oooOO0OOooo........oo 98 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo.....
94 99
95 PhysicsList::PhysicsList() : G4VModularPhysics 100 PhysicsList::PhysicsList() : G4VModularPhysicsList()
96 { 101 {
97 SetDefaultCutValue(0.7 * CLHEP::mm); << 102 G4LossTableManager::Instance();
98 verboseLevel = 1; << 103 defaultCutValue = 0.7*mm;
>> 104 cutForGamma = defaultCutValue;
>> 105 cutForElectron = defaultCutValue;
>> 106 cutForPositron = defaultCutValue;
>> 107 cutForProton = defaultCutValue;
>> 108 dump = false;
>> 109 verboseLevel = 1;
99 110
100 // Messengers << 111 pMessenger = new PhysicsListMessenger(this);
101 fMessenger = new PhysicsListMessenger(this); <<
102 fFactMessenger = new G4PhysListFactoryMessen <<
103 112
104 // Particles 113 // Particles
105 fParticleList = new G4DecayPhysics(verboseLe << 114 particleList = new G4DecayPhysics("decays");
106 115
107 // EM physics 116 // EM physics
108 fEmPhysicsList = new G4EmStandardPhysics(ver << 117 emPhysicsList = new G4EmStandardPhysics();
109 } 118 }
110 119
111 //....oooOO0OOooo........oooOO0OOooo........oo 120 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo.....
112 121
113 PhysicsList::~PhysicsList() 122 PhysicsList::~PhysicsList()
114 { 123 {
115 delete fMessenger; << 124 delete pMessenger;
116 delete fParticleList; << 125 delete particleList;
117 delete fEmPhysicsList; << 126 delete emPhysicsList;
118 for (std::size_t i = 0; i < fHadronPhys.size << 127 for(size_t i=0; i<hadronPhys.size(); i++) {
119 delete fHadronPhys[i]; << 128 delete hadronPhys[i];
120 } 129 }
121 } 130 }
122 131
123 //....oooOO0OOooo........oooOO0OOooo........oo 132 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo.....
124 133
125 void PhysicsList::ConstructParticle() 134 void PhysicsList::ConstructParticle()
126 { 135 {
127 fParticleList->ConstructParticle(); << 136 particleList->ConstructParticle();
128 } 137 }
129 138
130 //....oooOO0OOooo........oooOO0OOooo........oo 139 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo.....
131 140
132 void PhysicsList::ConstructProcess() 141 void PhysicsList::ConstructProcess()
133 { 142 {
134 AddTransportation(); 143 AddTransportation();
135 fEmPhysicsList->ConstructProcess(); << 144 emPhysicsList->ConstructProcess();
136 fParticleList->ConstructProcess(); << 145 particleList->ConstructProcess();
137 for (std::size_t i = 0; i < fHadronPhys.size << 146 for(size_t i=0; i<hadronPhys.size(); i++) {
138 fHadronPhys[i]->ConstructProcess(); << 147 hadronPhys[i]->ConstructProcess();
139 } 148 }
140 } 149 }
141 150
142 //....oooOO0OOooo........oooOO0OOooo........oo 151 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo.....
143 152
144 void PhysicsList::AddPhysicsList(const G4Strin 153 void PhysicsList::AddPhysicsList(const G4String& name)
145 { 154 {
146 if (verboseLevel > 0) { << 155 if (verboseLevel>0)
147 G4cout << "PhysicsList::AddPhysicsList: <" 156 G4cout << "PhysicsList::AddPhysicsList: <" << name << ">" << G4endl;
148 } << 157
149 if (name == "emstandard_opt0") { << 158 if (name == "emstandard_opt2") {
150 delete fEmPhysicsList; << 159
151 fEmPhysicsList = new G4EmStandardPhysics(v << 160 delete emPhysicsList;
152 } << 161 emPhysicsList = new G4EmStandardPhysics_option2();
153 else if (name == "emstandard_opt1") { << 162
154 delete fEmPhysicsList; << 163 } else if (name == "emstandard_opt3") {
155 fEmPhysicsList = new G4EmStandardPhysics_o << 164
156 } << 165 delete emPhysicsList;
157 else if (name == "emstandard_opt2") { << 166 emPhysicsList = new G4EmStandardPhysics_option3();
158 delete fEmPhysicsList; << 167
159 fEmPhysicsList = new G4EmStandardPhysics_o << 168 } else if (name == "emstandard_opt1") {
160 } << 169
161 else if (name == "emstandard_opt3") { << 170 delete emPhysicsList;
162 delete fEmPhysicsList; << 171 emPhysicsList = new G4EmStandardPhysics_option1();
163 fEmPhysicsList = new G4EmStandardPhysics_o << 172
164 } << 173 } else if (name == "emstandard_opt0") {
165 else if (name == "emstandard_opt4") { << 174
166 delete fEmPhysicsList; << 175 delete emPhysicsList;
167 fEmPhysicsList = new G4EmStandardPhysics_o << 176 emPhysicsList = new G4EmStandardPhysics();
168 } << 177
169 else if (name == "emstandardGS") { << 178 } else if (name == "FTFP_BERT_EMV") {
170 delete fEmPhysicsList; << 179
171 fEmPhysicsList = new G4EmStandardPhysicsGS <<
172 } <<
173 else if (name == "FTFP_BERT_EMV") { <<
174 AddPhysicsList("FTFP_BERT"); <<
175 AddPhysicsList("emstandard_opt1"); 180 AddPhysicsList("emstandard_opt1");
176 } <<
177 else if (name == "FTFP_BERT_EMX") { <<
178 AddPhysicsList("FTFP_BERT"); 181 AddPhysicsList("FTFP_BERT");
>> 182
>> 183 } else if (name == "FTFP_BERT_EMX") {
>> 184
179 AddPhysicsList("emstandard_opt2"); 185 AddPhysicsList("emstandard_opt2");
180 } <<
181 else if (name == "FTFP_BERT_EMY") { <<
182 AddPhysicsList("FTFP_BERT"); 186 AddPhysicsList("FTFP_BERT");
183 AddPhysicsList("emstandard_opt3"); << 187
184 } << 188 } else if (name == "FTFP_BERT") {
185 else if (name == "FTFP_BERT_EMZ") { << 189
186 AddPhysicsList("FTFP_BERT"); << 190 SetBuilderList1();
187 AddPhysicsList("emstandard_opt4"); << 191 hadronPhys.push_back( new HadronPhysicsFTFP_BERT("hadron",true));
188 } << 192 dump = true;
189 else if (name == "FTFP_BERT") { << 193
190 SetBuilderList0(false); << 194 } else if (name == "FTFP_BERT_TRV") {
191 fHadronPhys.push_back(new G4HadronPhysicsF << 195
192 } << 196 SetBuilderList1();
193 else if (name == "FTFP_BERT_TRV") { << 197 hadronPhys.push_back( new HadronPhysicsFTFP_BERT_TRV("hadron",true));
194 AddPhysicsList("emstandardGS"); << 198 dump = true;
195 SetBuilderList1(false); << 199
196 fHadronPhys.push_back(new G4HadronPhysicsF << 200 } else if (name == "FTF_BIC") {
197 } << 201
198 else if (name == "FTF_BIC") { << 202 SetBuilderList0();
199 SetBuilderList0(false); << 203 hadronPhys.push_back( new HadronPhysicsFTF_BIC("hadron",true));
200 fHadronPhys.push_back(new G4HadronPhysicsF << 204 dump = true;
201 } << 205
202 else if (name == "QBBC") { << 206 } else if (name == "LHEP") {
203 AddPhysicsList("emstandard_opt0"); << 207
204 SetBuilderList2(); 208 SetBuilderList2();
205 fHadronPhys.push_back(new G4HadronInelasti << 209 hadronPhys.push_back( new HadronPhysicsLHEP("hadron"));
206 } << 210 dump = true;
207 else if (name == "QGSP_BERT") { << 211
208 SetBuilderList0(false); << 212 } else if (name == "LHEP_EMV") {
209 fHadronPhys.push_back(new G4HadronPhysicsQ << 213
210 } <<
211 else if (name == "QGSP_FTFP_BERT") { <<
212 SetBuilderList0(false); <<
213 fHadronPhys.push_back(new G4HadronPhysicsQ <<
214 } <<
215 else if (name == "QGSP_FTFP_BERT_EMV") { <<
216 AddPhysicsList("QGSP_FTFP_BERT"); <<
217 AddPhysicsList("emstandard_opt1"); 214 AddPhysicsList("emstandard_opt1");
218 } << 215 SetBuilderList3();
219 else if (name == "QGSP_BERT_EMV") { << 216 hadronPhys.push_back( new HadronPhysicsLHEP_EMV("hadron"));
220 AddPhysicsList("QGSP_BERT"); << 217 dump = true;
>> 218
>> 219 } else if (name == "QBBC") {
>> 220
>> 221 AddPhysicsList("emstandard_opt2");
>> 222 SetBuilderList6();
>> 223 hadronPhys.push_back( new G4HadronInelasticQBBC("QBBC",verboseLevel));
>> 224
>> 225 } else if (name == "QBBC_XGG") {
>> 226
>> 227 AddPhysicsList("emstandard_opt2");
>> 228 SetBuilderList6();
>> 229 hadronPhys.push_back( new G4HadronInelasticQBBC("QBBC_XGG",verboseLevel));
>> 230
>> 231 } else if (name == "QBBC_XGGSN") {
>> 232
>> 233 AddPhysicsList("emstandard_opt2");
>> 234 SetBuilderList6();
>> 235 hadronPhys.push_back( new G4HadronInelasticQBBC("QBBC_XGG",verboseLevel));
>> 236
>> 237
>> 238 } else if (name == "QGSC_BERT") {
>> 239
>> 240 SetBuilderList4();
>> 241 hadronPhys.push_back( new HadronPhysicsQGSC_BERT("hadron",true));
>> 242 dump = true;
>> 243
>> 244 } else if (name == "QGSC_CHIPS") {
>> 245
>> 246 SetBuilderList4();
>> 247 hadronPhys.push_back( new HadronPhysicsQGSC_CHIPS("hadron",true));
>> 248 dump = true;
>> 249
>> 250 } else if (name == "QGSP_BERT") {
>> 251
>> 252 SetBuilderList1();
>> 253 hadronPhys.push_back( new HadronPhysicsQGSP_BERT("hadron",true));
>> 254 dump = true;
>> 255
>> 256 } else if (name == "QGSP_FTFP_BERT") {
>> 257
>> 258 SetBuilderList1();
>> 259 hadronPhys.push_back( new HadronPhysicsQGSP_FTFP_BERT("hadron",true));
>> 260 dump = true;
>> 261
>> 262 } else if (name == "QGSP_BERT_EMV") {
>> 263
221 AddPhysicsList("emstandard_opt1"); 264 AddPhysicsList("emstandard_opt1");
222 } <<
223 else if (name == "QGSP_BERT_EMX") { <<
224 AddPhysicsList("QGSP_BERT"); 265 AddPhysicsList("QGSP_BERT");
>> 266
>> 267 } else if (name == "QGSP_BERT_EMX") {
>> 268
225 AddPhysicsList("emstandard_opt2"); 269 AddPhysicsList("emstandard_opt2");
226 } << 270 AddPhysicsList("QGSP_BERT");
227 else if (name == "QGSP_BERT_HP") { << 271
228 SetBuilderList0(true); << 272 } else if (name == "QGSP_BERT_HP") {
229 fHadronPhys.push_back(new G4HadronPhysicsQ << 273
230 AddPhysicsList("RadioactiveDecay"); << 274 SetBuilderList1(true);
231 } << 275 hadronPhys.push_back( new HadronPhysicsQGSP_BERT_HP("hadron",true));
232 else if (name == "QGSP_BIC") { << 276
233 SetBuilderList0(false); << 277 } else if (name == "QGSP_BIC") {
234 fHadronPhys.push_back(new G4HadronPhysicsQ << 278
235 fHadronPhys.push_back(new G4IonElasticPhys << 279 SetBuilderList0();
236 } << 280 hadronPhys.push_back( new HadronPhysicsQGSP_BIC("hadron",true));
237 else if (name == "QGSP_BIC_EMY") { << 281 dump = true;
238 AddPhysicsList("QGSP_BIC"); << 282
>> 283 } else if (name == "QGSP_BIC_EMY") {
>> 284
239 AddPhysicsList("emstandard_opt3"); 285 AddPhysicsList("emstandard_opt3");
240 } << 286 SetBuilderList0();
241 else if (name == "QGS_BIC") { << 287 hadronPhys.push_back( new HadronPhysicsQGSP_BIC("hadron",true));
242 SetBuilderList0(false); << 288 dump = true;
243 fHadronPhys.push_back(new G4HadronPhysicsQ << 289
244 } << 290 } else if (name == "QGS_BIC") {
245 else if (name == "QGSP_BIC_HP") { << 291
>> 292 SetBuilderList0();
>> 293 hadronPhys.push_back( new HadronPhysicsQGS_BIC("hadron",true));
>> 294 dump = true;
>> 295
>> 296 } else if (name == "QGSP_BIC_HP") {
>> 297
246 SetBuilderList0(true); 298 SetBuilderList0(true);
247 fHadronPhys.push_back(new G4HadronPhysicsQ << 299 hadronPhys.push_back( new HadronPhysicsQGSP_BIC_HP("hadron",true));
248 AddPhysicsList("RadioactiveDecay"); << 300 dump = true;
249 } << 301
250 else if (name == "RadioactiveDecay") { << 302 } else if (name == "QGSP_BIC_EMY") {
251 fHadronPhys.push_back(new G4RadioactiveDec << 303
252 } << 304 AddPhysicsList("emstandard_opt3");
253 else { << 305 AddPhysicsList("QGSP_BIC");
>> 306
>> 307 } else {
>> 308
254 G4cout << "PhysicsList::AddPhysicsList: <" 309 G4cout << "PhysicsList::AddPhysicsList: <" << name << ">"
255 << " is not defined" << G4endl; << 310 << " is not defined"
>> 311 << G4endl;
256 } 312 }
257 } 313 }
258 314
259 //....oooOO0OOooo........oooOO0OOooo........oo 315 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo.....
260 316
261 void PhysicsList::SetBuilderList0(G4bool flagH 317 void PhysicsList::SetBuilderList0(G4bool flagHP)
262 { 318 {
263 fHadronPhys.push_back(new G4EmExtraPhysics(v << 319 hadronPhys.push_back( new G4EmExtraPhysics("extra EM"));
264 if (flagHP) { << 320 hadronPhys.push_back( new G4HadronElasticPhysics("elastic",verboseLevel,
265 fHadronPhys.push_back(new G4HadronElasticP << 321 flagHP));
266 } << 322 hadronPhys.push_back( new G4QStoppingPhysics("stopping",verboseLevel));
267 else { << 323 hadronPhys.push_back( new G4IonBinaryCascadePhysics("ionBIC"));
268 fHadronPhys.push_back(new G4HadronElasticP << 324 hadronPhys.push_back( new G4NeutronTrackingCut("nTackingCut",verboseLevel));
269 } <<
270 fHadronPhys.push_back(new G4StoppingPhysics( <<
271 fHadronPhys.push_back(new G4IonPhysics(verbo <<
272 if (!flagHP) { <<
273 fHadronPhys.push_back(new G4NeutronTrackin <<
274 } <<
275 } 325 }
276 326
277 //....oooOO0OOooo........oooOO0OOooo........oo 327 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo.....
278 328
279 void PhysicsList::SetBuilderList1(G4bool flagH 329 void PhysicsList::SetBuilderList1(G4bool flagHP)
280 { 330 {
281 fHadronPhys.push_back(new G4EmExtraPhysics(v << 331 hadronPhys.push_back( new G4EmExtraPhysics("extra EM"));
282 if (flagHP) { << 332 hadronPhys.push_back( new G4HadronHElasticPhysics(verboseLevel,flagHP));
283 fHadronPhys.push_back(new G4HadronElasticP << 333 // hadronPhys.push_back( new G4HadronElasticPhysics("elastic",verboseLevel,flagHP));
284 } << 334 hadronPhys.push_back( new G4QStoppingPhysics("stopping",verboseLevel));
285 else { << 335 hadronPhys.push_back( new G4IonPhysics("ion"));
286 fHadronPhys.push_back(new G4HadronHElastic << 336 hadronPhys.push_back( new G4NeutronTrackingCut("nTackingCut",verboseLevel));
287 } << 337 }
288 fHadronPhys.push_back(new G4StoppingPhysics( << 338
289 fHadronPhys.push_back(new G4IonPhysics(verbo << 339 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo.....
290 if (!flagHP) { << 340
291 fHadronPhys.push_back(new G4NeutronTrackin << 341 void PhysicsList::SetBuilderList2(G4bool flagHP)
292 } << 342 {
>> 343 hadronPhys.push_back( new G4EmExtraPhysics("extra EM"));
>> 344 hadronPhys.push_back( new G4HadronElasticPhysics("LElastic",verboseLevel,
>> 345 flagHP));
>> 346 hadronPhys.push_back( new G4IonPhysics("ion"));
293 } 347 }
294 348
295 //....oooOO0OOooo........oooOO0OOooo........oo 349 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo.....
296 350
297 void PhysicsList::SetBuilderList2() << 351 void PhysicsList::SetBuilderList3(G4bool flagHP)
298 { 352 {
299 fHadronPhys.push_back(new G4EmExtraPhysics(v << 353 hadronPhys.push_back( new G4EmExtraPhysics("extra EM"));
300 fHadronPhys.push_back(new G4HadronElasticPhy << 354 hadronPhys.push_back( new G4HadronElasticPhysics("LElastic",verboseLevel,
301 fHadronPhys.push_back(new G4StoppingPhysics( << 355 flagHP));
302 fHadronPhys.push_back(new G4IonPhysicsXS(ver << 356 hadronPhys.push_back( new G4QStoppingPhysics("stopping",verboseLevel));
303 fHadronPhys.push_back(new G4IonElasticPhysic << 357 hadronPhys.push_back( new G4IonPhysics("ion"));
304 fHadronPhys.push_back(new G4NeutronTrackingC <<
305 } 358 }
306 359
307 //....oooOO0OOooo........oooOO0OOooo........oo 360 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo.....
308 361
>> 362 void PhysicsList::SetBuilderList4(G4bool)
>> 363 {
>> 364 hadronPhys.push_back( new G4EmExtraPhysics("extra EM"));
>> 365 hadronPhys.push_back( new G4HadronQElasticPhysics("elastic",verboseLevel));
>> 366 hadronPhys.push_back( new G4QStoppingPhysics("stopping",verboseLevel));
>> 367 hadronPhys.push_back( new G4IonPhysics("ion"));
>> 368 hadronPhys.push_back( new G4NeutronTrackingCut("nTackingCut",verboseLevel));
>> 369 }
>> 370
>> 371 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo.....
>> 372
>> 373 void PhysicsList::SetBuilderList5(G4bool flagHP)
>> 374 {
>> 375 hadronPhys.push_back( new G4EmExtraPhysics("extra EM"));
>> 376 hadronPhys.push_back( new G4HadronDElasticPhysics(verboseLevel,flagHP));
>> 377 hadronPhys.push_back( new G4QStoppingPhysics("stopping",verboseLevel));
>> 378 hadronPhys.push_back( new G4IonBinaryCascadePhysics("ionBIC"));
>> 379 hadronPhys.push_back( new G4NeutronTrackingCut("nTackingCut",verboseLevel));
>> 380 }
>> 381
>> 382 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo.....
>> 383
>> 384 void PhysicsList::SetBuilderList6(G4bool)
>> 385 {
>> 386 hadronPhys.push_back( new G4EmExtraPhysics("extra EM"));
>> 387 hadronPhys.push_back( new G4HadronHElasticPhysics(verboseLevel,false));
>> 388 hadronPhys.push_back( new G4QStoppingPhysics("stopping",verboseLevel));
>> 389 hadronPhys.push_back( new G4IonBinaryCascadePhysics("ionBIC"));
>> 390 hadronPhys.push_back( new G4NeutronTrackingCut("nTackingCut",verboseLevel));
>> 391 }
>> 392
>> 393 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo.....
>> 394
>> 395 void PhysicsList::SetCuts()
>> 396 {
>> 397
>> 398 if (verboseLevel >0){
>> 399 G4cout << "PhysicsList::SetCuts:";
>> 400 G4cout << "CutLength : " << G4BestUnit(defaultCutValue,"Length") << G4endl;
>> 401 }
>> 402
>> 403 // set cut values for gamma at first and for e- second and next for e+,
>> 404 // because some processes for e+/e- need cut values for gamma
>> 405 SetCutValue(cutForGamma, "gamma");
>> 406 SetCutValue(cutForElectron, "e-");
>> 407 SetCutValue(cutForPositron, "e+");
>> 408 SetCutValue(cutForProton, "proton");
>> 409
>> 410 if (verboseLevel>0) DumpCutValuesTable();
>> 411 }
>> 412
>> 413 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo.....
>> 414
>> 415 void PhysicsList::SetCutForGamma(G4double cut)
>> 416 {
>> 417 cutForGamma = cut;
>> 418 SetParticleCuts(cutForGamma, G4Gamma::Gamma());
>> 419 }
>> 420
>> 421 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
>> 422
>> 423 void PhysicsList::SetCutForElectron(G4double cut)
>> 424 {
>> 425 cutForElectron = cut;
>> 426 SetParticleCuts(cutForElectron, G4Electron::Electron());
>> 427 }
>> 428
>> 429 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
>> 430
>> 431 void PhysicsList::SetCutForPositron(G4double cut)
>> 432 {
>> 433 cutForPositron = cut;
>> 434 SetParticleCuts(cutForPositron, G4Positron::Positron());
>> 435 }
>> 436
>> 437 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
>> 438
>> 439 void PhysicsList::SetCutForProton(G4double cut)
>> 440 {
>> 441 cutForProton = cut;
>> 442 SetParticleCuts(cutForProton, G4Proton::Proton());
>> 443 }
>> 444
309 void PhysicsList::List() 445 void PhysicsList::List()
310 { 446 {
311 G4cout << "### PhysicsLists available: FTFP_ << 447 G4cout << "### PhysicsLists available: FTFP_BERT FTFP_BERT_EMV FTFP_BERT_EMX FTFP_BERT_TRV FTF_BIC"
312 << "FTFP_BERT_EMX FTFP_BERT_EMZ FTFP_ << 448 << G4endl;
313 G4cout << " FTF_B << 449 G4cout << " LHEP LHEP_EMV QBBC QBBC_XGG QBBC_XGGSN"
314 << "QGSP_BERT_EMV QGSP_BERT_EMX" << G << 450 << G4endl;
315 G4cout << " QGSP_ << 451 G4cout << " QGSC_BERT QGSC_CHIPS QGSP_BERT QGSP_BERT_EMV QGSP_BIC_EMY"
316 << "QGSP_FTFP_BERT_EMV" << G4endl; << 452 << G4endl;
317 G4cout << " QGS_B << 453 G4cout << " QGSP_BERT_EMX QGSP_BERT_HP QGSP_BIC QGSP_BIC_HP"
318 << "QGSP_BIC_HP" << G4endl; << 454 << G4endl;
319 } 455 }
320 456
321 //....oooOO0OOooo........oooOO0OOooo........oo 457 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
>> 458
322 459