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