| 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.17 2006/11/24 16:48:57 vnivanch Exp $
>> 28 // GEANT4 tag $Name: geant4-08-02-patch-01-ref $
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 <<
40 // 16.10.2012 Renamed used classes (A.Ribon) <<
41 // 38 //
42 ////////////////////////////////////////////// 39 ////////////////////////////////////////////////////////////////////////
43 // << 40 //
44 41
45 #include "PhysicsList.hh" 42 #include "PhysicsList.hh"
46 <<
47 #include "PhysicsListMessenger.hh" 43 #include "PhysicsListMessenger.hh"
48 44
49 #include "G4DecayPhysics.hh" 45 #include "G4DecayPhysics.hh"
50 #include "G4Electron.hh" << 46 #include "G4EmStandardPhysics72.hh"
51 #include "G4EmExtraPhysics.hh" << 47 #include "G4EmStandardPhysics71.hh"
52 #include "G4EmLivermorePhysics.hh" <<
53 #include "G4EmParameters.hh" <<
54 #include "G4EmPenelopePhysics.hh" <<
55 #include "G4EmStandardPhysics.hh" 48 #include "G4EmStandardPhysics.hh"
56 #include "G4EmStandardPhysicsGS.hh" <<
57 #include "G4EmStandardPhysics_option1.hh" <<
58 #include "G4EmStandardPhysics_option2.hh" <<
59 #include "G4EmStandardPhysics_option3.hh" <<
60 #include "G4EmStandardPhysics_option4.hh" <<
61 #include "G4Gamma.hh" <<
62 #include "G4HadronElasticPhysics.hh" 49 #include "G4HadronElasticPhysics.hh"
63 #include "G4HadronElasticPhysicsHP.hh" << 50 #include "G4HadronQElasticPhysics.hh"
64 #include "G4HadronElasticPhysicsXS.hh" <<
65 #include "G4HadronHElasticPhysics.hh" 51 #include "G4HadronHElasticPhysics.hh"
66 #include "G4HadronInelasticQBBC.hh" << 52 #include "G4NeutronTrackingCut.hh"
67 #include "G4HadronPhysicsFTFP_BERT.hh" << 53 #include "G4QStoppingPhysics.hh"
68 #include "G4HadronPhysicsFTFP_BERT_HP.hh" << 54 #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" 55 #include "G4IonBinaryCascadePhysics.hh"
78 #include "G4IonElasticPhysics.hh" <<
79 #include "G4IonPhysics.hh" 56 #include "G4IonPhysics.hh"
80 #include "G4IonPhysicsXS.hh" << 57 #include "G4EmExtraPhysics.hh"
81 #include "G4NeutronCrossSectionXS.hh" << 58 #include "G4EmProcessOptions.hh"
82 #include "G4NeutronTrackingCut.hh" << 59
83 #include "G4ParticleTable.hh" << 60 #include "HadronPhysicsQGSP_BIC.hh"
>> 61 #include "HadronPhysicsLHEP.hh"
>> 62 #include "HadronPhysicsLHEP_BIC.hh"
>> 63 #include "HadronPhysicsLHEP_BIC_HP.hh"
>> 64 #include "HadronPhysicsLHEP_BERT.hh"
>> 65 #include "HadronPhysicsLHEP_BERT_HP.hh"
>> 66 #include "HadronPhysicsQGSP_BERT_HP.hh"
>> 67 #include "HadronPhysicsQGSP_BERT.hh"
>> 68 #include "HadronPhysicsQGSC.hh"
>> 69 #include "HadronPhysicsQGSP.hh"
>> 70 #include "G4HadronInelasticQBBC.hh"
>> 71 #include "G4HadronInelasticQLHEP.hh"
>> 72 #include "G4IonPhysics.hh"
>> 73 #include "G4HadronProcessStore.hh"
>> 74
>> 75 #include "G4LossTableManager.hh"
>> 76
>> 77 #include "G4ProcessManager.hh"
84 #include "G4ParticleTypes.hh" 78 #include "G4ParticleTypes.hh"
85 #include "G4PhysListFactoryMessenger.hh" << 79 #include "G4ParticleTable.hh"
>> 80 #include "G4Gamma.hh"
>> 81 #include "G4Electron.hh"
86 #include "G4Positron.hh" 82 #include "G4Positron.hh"
87 #include "G4ProcessManager.hh" <<
88 #include "G4Proton.hh" <<
89 #include "G4RadioactiveDecayPhysics.hh" <<
90 #include "G4StoppingPhysics.hh" <<
91 #include "G4SystemOfUnits.hh" <<
92 83
93 //....oooOO0OOooo........oooOO0OOooo........oo 84 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo.....
94 85
95 PhysicsList::PhysicsList() : G4VModularPhysics 86 PhysicsList::PhysicsList() : G4VModularPhysicsList()
96 { 87 {
97 SetDefaultCutValue(0.7 * CLHEP::mm); << 88 G4LossTableManager::Instance();
98 verboseLevel = 1; << 89 defaultCutValue = 1.*mm;
>> 90 cutForGamma = defaultCutValue;
>> 91 cutForElectron = defaultCutValue;
>> 92 cutForPositron = defaultCutValue;
>> 93 dump = false;
99 94
100 // Messengers << 95 pMessenger = new PhysicsListMessenger(this);
101 fMessenger = new PhysicsListMessenger(this); <<
102 fFactMessenger = new G4PhysListFactoryMessen <<
103 96
104 // Particles 97 // Particles
105 fParticleList = new G4DecayPhysics(verboseLe << 98 particleList = new G4DecayPhysics("decays");
106 99
107 // EM physics 100 // EM physics
108 fEmPhysicsList = new G4EmStandardPhysics(ver << 101 emPhysicsList = new G4EmStandardPhysics("standard");
109 } 102 }
110 103
111 //....oooOO0OOooo........oooOO0OOooo........oo 104 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo.....
112 105
113 PhysicsList::~PhysicsList() 106 PhysicsList::~PhysicsList()
114 { 107 {
115 delete fMessenger; << 108 delete pMessenger;
116 delete fParticleList; << 109 delete particleList;
117 delete fEmPhysicsList; << 110 delete emPhysicsList;
118 for (std::size_t i = 0; i < fHadronPhys.size << 111 for(size_t i=0; i<hadronPhys.size(); i++) {
119 delete fHadronPhys[i]; << 112 delete hadronPhys[i];
120 } 113 }
121 } 114 }
122 115
123 //....oooOO0OOooo........oooOO0OOooo........oo 116 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo.....
124 117
125 void PhysicsList::ConstructParticle() 118 void PhysicsList::ConstructParticle()
126 { 119 {
127 fParticleList->ConstructParticle(); << 120 particleList->ConstructParticle();
128 } 121 }
129 122
130 //....oooOO0OOooo........oooOO0OOooo........oo 123 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo.....
131 124
132 void PhysicsList::ConstructProcess() 125 void PhysicsList::ConstructProcess()
133 { 126 {
134 AddTransportation(); 127 AddTransportation();
135 fEmPhysicsList->ConstructProcess(); << 128 emPhysicsList->ConstructProcess();
136 fParticleList->ConstructProcess(); << 129 particleList->ConstructProcess();
137 for (std::size_t i = 0; i < fHadronPhys.size << 130 for(size_t i=0; i<hadronPhys.size(); i++) {
138 fHadronPhys[i]->ConstructProcess(); << 131 hadronPhys[i]->ConstructProcess();
139 } 132 }
>> 133
>> 134 if(dump) G4HadronProcessStore::Instance()->Dump(1);
140 } 135 }
141 136
142 //....oooOO0OOooo........oooOO0OOooo........oo 137 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo.....
143 138
144 void PhysicsList::AddPhysicsList(const G4Strin 139 void PhysicsList::AddPhysicsList(const G4String& name)
145 { 140 {
146 if (verboseLevel > 0) { << 141 if (verboseLevel>0)
147 G4cout << "PhysicsList::AddPhysicsList: <" 142 G4cout << "PhysicsList::AddPhysicsList: <" << name << ">" << G4endl;
148 } << 143
149 if (name == "emstandard_opt0") { << 144 if (name == "LHEP") {
150 delete fEmPhysicsList; << 145
151 fEmPhysicsList = new G4EmStandardPhysics(v << 146 hadronPhys.push_back( new G4EmExtraPhysics("gamma_nuc"));
152 } << 147 hadronPhys.push_back( new G4HadronElasticPhysics("LElastic",
153 else if (name == "emstandard_opt1") { << 148 verboseLevel,false));
154 delete fEmPhysicsList; << 149 hadronPhys.push_back( new HadronPhysicsLHEP());
155 fEmPhysicsList = new G4EmStandardPhysics_o << 150 hadronPhys.push_back( new G4IonPhysics("ion"));
156 } << 151
157 else if (name == "emstandard_opt2") { << 152 } else if (name == "G4standard_exp") {
158 delete fEmPhysicsList; << 153
159 fEmPhysicsList = new G4EmStandardPhysics_o << 154 delete emPhysicsList;
160 } << 155 emPhysicsList = new G4EmStandardPhysics72("standard");
161 else if (name == "emstandard_opt3") { << 156
162 delete fEmPhysicsList; << 157 } else if (name == "G4standard_fast") {
163 fEmPhysicsList = new G4EmStandardPhysics_o << 158
164 } << 159 delete emPhysicsList;
165 else if (name == "emstandard_opt4") { << 160 emPhysicsList = new G4EmStandardPhysics71("standard");
166 delete fEmPhysicsList; << 161
167 fEmPhysicsList = new G4EmStandardPhysics_o << 162 } else if (name == "QGSC") {
168 } << 163
169 else if (name == "emstandardGS") { << 164 hadronPhys.push_back( new G4EmExtraPhysics("gamma_nuc"));
170 delete fEmPhysicsList; << 165 hadronPhys.push_back( new G4HadronQElasticPhysics("QElastic",verboseLevel));
171 fEmPhysicsList = new G4EmStandardPhysicsGS << 166 hadronPhys.push_back( new HadronPhysicsQGSC());
172 } << 167 hadronPhys.push_back( new G4QStoppingPhysics("stopping"));
173 else if (name == "FTFP_BERT_EMV") { << 168 hadronPhys.push_back( new G4IonPhysics("ion"));
174 AddPhysicsList("FTFP_BERT"); << 169 hadronPhys.push_back( new G4NeutronTrackingCut());
175 AddPhysicsList("emstandard_opt1"); << 170 dump = true;
176 } << 171
177 else if (name == "FTFP_BERT_EMX") { << 172 } else if (name == "QGSP") {
178 AddPhysicsList("FTFP_BERT"); << 173
179 AddPhysicsList("emstandard_opt2"); << 174 hadronPhys.push_back( new G4EmExtraPhysics("gamma_nuc"));
180 } << 175 hadronPhys.push_back( new G4HadronElasticPhysics("elastic",
181 else if (name == "FTFP_BERT_EMY") { << 176 verboseLevel,false));
182 AddPhysicsList("FTFP_BERT"); << 177 hadronPhys.push_back( new HadronPhysicsQGSP());
183 AddPhysicsList("emstandard_opt3"); << 178 hadronPhys.push_back( new G4QStoppingPhysics("stopping"));
184 } << 179 hadronPhys.push_back( new G4IonPhysics("ion"));
185 else if (name == "FTFP_BERT_EMZ") { << 180 hadronPhys.push_back( new G4NeutronTrackingCut());
186 AddPhysicsList("FTFP_BERT"); << 181 dump = true;
187 AddPhysicsList("emstandard_opt4"); << 182
188 } << 183 } else if (name == "LHEP_EMV") {
189 else if (name == "FTFP_BERT") { << 184
190 SetBuilderList0(false); << 185 AddPhysicsList("G4standard_fast");
191 fHadronPhys.push_back(new G4HadronPhysicsF << 186 AddPhysicsList("LHEP");
192 } << 187 dump = true;
193 else if (name == "FTFP_BERT_TRV") { << 188
194 AddPhysicsList("emstandardGS"); << 189 } else if (name == "QGSP_EMV") {
195 SetBuilderList1(false); << 190
196 fHadronPhys.push_back(new G4HadronPhysicsF << 191 AddPhysicsList("G4standard_fast");
197 } << 192 AddPhysicsList("QGSP");
198 else if (name == "FTF_BIC") { << 193 dump = true;
199 SetBuilderList0(false); << 194
200 fHadronPhys.push_back(new G4HadronPhysicsF << 195 } else if (name == "QGSP_EMX") {
201 } << 196
202 else if (name == "QBBC") { << 197 AddPhysicsList("G4standard_exp");
203 AddPhysicsList("emstandard_opt0"); << 198 AddPhysicsList("QGSP");
204 SetBuilderList2(); << 199 dump = true;
205 fHadronPhys.push_back(new G4HadronInelasti << 200
206 } << 201 } else if (name == "QGSP_BERT") {
207 else if (name == "QGSP_BERT") { << 202
208 SetBuilderList0(false); << 203 hadronPhys.push_back( new G4EmExtraPhysics("gamma_nuc"));
209 fHadronPhys.push_back(new G4HadronPhysicsQ << 204 hadronPhys.push_back( new G4HadronElasticPhysics("elastic",
210 } << 205 verboseLevel,false));
211 else if (name == "QGSP_FTFP_BERT") { << 206 hadronPhys.push_back( new HadronPhysicsQGSP_BERT());
212 SetBuilderList0(false); << 207 hadronPhys.push_back( new G4QStoppingPhysics("stopping"));
213 fHadronPhys.push_back(new G4HadronPhysicsQ << 208 hadronPhys.push_back( new G4IonPhysics("ion"));
214 } << 209 hadronPhys.push_back( new G4NeutronTrackingCut());
215 else if (name == "QGSP_FTFP_BERT_EMV") { << 210 dump = true;
216 AddPhysicsList("QGSP_FTFP_BERT"); << 211
217 AddPhysicsList("emstandard_opt1"); << 212 } else if (name == "QGSP_BERT_HP") {
218 } << 213
219 else if (name == "QGSP_BERT_EMV") { << 214 hadronPhys.push_back( new G4EmExtraPhysics("gamma_nuc"));
220 AddPhysicsList("QGSP_BERT"); << 215 hadronPhys.push_back( new G4HadronElasticPhysics("elastic",
221 AddPhysicsList("emstandard_opt1"); << 216 verboseLevel,true));
222 } << 217 hadronPhys.push_back( new HadronPhysicsQGSP_BERT_HP());
223 else if (name == "QGSP_BERT_EMX") { << 218 hadronPhys.push_back( new G4QStoppingPhysics("stopping"));
224 AddPhysicsList("QGSP_BERT"); << 219 hadronPhys.push_back( new G4IonPhysics("ion"));
225 AddPhysicsList("emstandard_opt2"); << 220 hadronPhys.push_back( new G4NeutronTrackingCut());
226 } << 221 dump = true;
227 else if (name == "QGSP_BERT_HP") { << 222
228 SetBuilderList0(true); << 223 } else if (name == "QGSP_BIC") {
229 fHadronPhys.push_back(new G4HadronPhysicsQ << 224
230 AddPhysicsList("RadioactiveDecay"); << 225 SetStandardList(false);
231 } << 226 hadronPhys.push_back( new HadronPhysicsQGSP_BIC());
232 else if (name == "QGSP_BIC") { << 227 dump = true;
233 SetBuilderList0(false); << 228
234 fHadronPhys.push_back(new G4HadronPhysicsQ << 229 } else if (name == "QGSP_BIC_HP") {
235 fHadronPhys.push_back(new G4IonElasticPhys << 230
236 } << 231 SetStandardList(true);
237 else if (name == "QGSP_BIC_EMY") { << 232 hadronPhys.push_back( new HadronPhysicsQGSP_BIC());
238 AddPhysicsList("QGSP_BIC"); << 233 dump = true;
239 AddPhysicsList("emstandard_opt3"); << 234
240 } << 235 } else if (name == "QBBC") {
241 else if (name == "QGS_BIC") { << 236
242 SetBuilderList0(false); << 237 SetStandardList(false);
243 fHadronPhys.push_back(new G4HadronPhysicsQ << 238 hadronPhys.push_back( new G4HadronInelasticQBBC("QBBC",verboseLevel));
244 } << 239
245 else if (name == "QGSP_BIC_HP") { << 240 } else if (name == "QBBCG") {
246 SetBuilderList0(true); << 241
247 fHadronPhys.push_back(new G4HadronPhysicsQ << 242 SetStandardList(false, true);
248 AddPhysicsList("RadioactiveDecay"); << 243 hadronPhys.push_back( new G4HadronInelasticQBBC("QBBC",verboseLevel,
249 } << 244 false,false,false,false,true));
250 else if (name == "RadioactiveDecay") { << 245
251 fHadronPhys.push_back(new G4RadioactiveDec << 246 } else if (name == "QBEC") {
252 } << 247
253 else { << 248 SetStandardList(false);
>> 249 hadronPhys.push_back( new G4HadronInelasticQBBC("QBBC",verboseLevel,
>> 250 false,true,false,false));
>> 251
>> 252 } else if (name == "QBBC_HP") {
>> 253
>> 254 SetStandardList(true);
>> 255 hadronPhys.push_back( new G4HadronInelasticQBBC("QBBC",verboseLevel,
>> 256 false,false,false,true));
>> 257
>> 258 } else if (name == "QBEC_HP") {
>> 259
>> 260 SetStandardList(true);
>> 261 hadronPhys.push_back( new G4HadronInelasticQBBC("QBBC",verboseLevel,
>> 262 false,true,false,true));
>> 263
>> 264 } else if (name == "QBBC_CHIPS") {
>> 265
>> 266 SetStandardList(false);
>> 267 hadronPhys.push_back( new G4HadronInelasticQBBC("QBBC",verboseLevel,
>> 268 false,false,true,false));
>> 269
>> 270 } else if (name == "QBBC_FTF") {
>> 271
>> 272 SetStandardList(false);
>> 273 hadronPhys.push_back( new G4HadronInelasticQBBC("QBBC",verboseLevel,
>> 274 true,false,false,false));
>> 275
>> 276 } else {
>> 277
254 G4cout << "PhysicsList::AddPhysicsList: <" 278 G4cout << "PhysicsList::AddPhysicsList: <" << name << ">"
255 << " is not defined" << G4endl; << 279 << " is not defined"
>> 280 << G4endl;
256 } 281 }
257 } 282 }
258 283
259 //....oooOO0OOooo........oooOO0OOooo........oo 284 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo.....
260 285
261 void PhysicsList::SetBuilderList0(G4bool flagH << 286 void PhysicsList::SetStandardList(G4bool flagHP, G4bool glauber)
262 { 287 {
263 fHadronPhys.push_back(new G4EmExtraPhysics(v << 288 hadronPhys.push_back( new G4EmExtraPhysics("gamma_nuc"));
264 if (flagHP) { << 289 hadronPhys.push_back( new G4HadronHElasticPhysics("elastic",verboseLevel,
265 fHadronPhys.push_back(new G4HadronElasticP << 290 flagHP,glauber));
266 } << 291 hadronPhys.push_back( new G4QStoppingPhysics("stopping"));
267 else { << 292 hadronPhys.push_back( new G4IonBinaryCascadePhysics("binary_ion"));
268 fHadronPhys.push_back(new G4HadronElasticP << 293 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 } 294 }
276 295
277 //....oooOO0OOooo........oooOO0OOooo........oo 296 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo.....
278 297
279 void PhysicsList::SetBuilderList1(G4bool flagH << 298 void PhysicsList::SetCuts()
280 { 299 {
281 fHadronPhys.push_back(new G4EmExtraPhysics(v << 300
282 if (flagHP) { << 301 if (verboseLevel >0){
283 fHadronPhys.push_back(new G4HadronElasticP << 302 G4cout << "PhysicsList::SetCuts:";
284 } << 303 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 } 304 }
>> 305
>> 306 // set cut values for gamma at first and for e- second and next for e+,
>> 307 // because some processes for e+/e- need cut values for gamma
>> 308 SetCutValue(cutForGamma, "gamma");
>> 309 SetCutValue(cutForElectron, "e-");
>> 310 SetCutValue(cutForPositron, "e+");
>> 311
>> 312 if (verboseLevel>0) DumpCutValuesTable();
293 } 313 }
294 314
295 //....oooOO0OOooo........oooOO0OOooo........oo 315 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo.....
296 316
297 void PhysicsList::SetBuilderList2() << 317 void PhysicsList::SetCutForGamma(G4double cut)
298 { 318 {
299 fHadronPhys.push_back(new G4EmExtraPhysics(v << 319 cutForGamma = cut;
300 fHadronPhys.push_back(new G4HadronElasticPhy << 320 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 } 321 }
306 322
307 //....oooOO0OOooo........oooOO0OOooo........oo << 323 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
>> 324
>> 325 void PhysicsList::SetCutForElectron(G4double cut)
>> 326 {
>> 327 cutForElectron = cut;
>> 328 SetParticleCuts(cutForElectron, G4Electron::Electron());
>> 329 }
>> 330
>> 331 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
>> 332
>> 333 void PhysicsList::SetCutForPositron(G4double cut)
>> 334 {
>> 335 cutForPositron = cut;
>> 336 SetParticleCuts(cutForPositron, G4Positron::Positron());
>> 337 }
308 338
309 void PhysicsList::List() 339 void PhysicsList::List()
310 { 340 {
311 G4cout << "### PhysicsLists available: FTFP_ << 341 G4cout << "### PhysicsLists available: LHEP LHEP_EMV QGSP QGSP_BERT QGSP_BERT_HP QGSP_BIC QGSP_BIC_HP"
312 << "FTFP_BERT_EMX FTFP_BERT_EMZ FTFP_ << 342 << G4endl;
313 G4cout << " FTF_B << 343 G4cout << " QGSP_EMV QGSP_EMX QGSC QBBC QBEC QBBC_HP QBEC_HP QBBC_CHIPS QBBC_FTF"
314 << "QGSP_BERT_EMV QGSP_BERT_EMX" << G << 344 << G4endl;
315 G4cout << " QGSP_ <<
316 << "QGSP_FTFP_BERT_EMV" << G4endl; <<
317 G4cout << " QGS_B <<
318 << "QGSP_BIC_HP" << G4endl; <<
319 } 345 }
320 346
321 //....oooOO0OOooo........oooOO0OOooo........oo 347 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
>> 348
322 349