| 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 // 26 //
27 // ------------------------------------------- 27 // --------------------------------------------------------------
28 // GEANT 4 - Underground Dark Matter Detecto 28 // GEANT 4 - Underground Dark Matter Detector Advanced Example
29 // 29 //
30 // For information related to this code c 30 // For information related to this code contact: Alex Howard
31 // e-mail: alexander.howard@cern.ch 31 // e-mail: alexander.howard@cern.ch
32 // ------------------------------------------- 32 // --------------------------------------------------------------
33 // Comments 33 // Comments
34 // 34 //
35 // Underground Advanced 35 // Underground Advanced
36 // by A. Howard and H. Araujo 36 // by A. Howard and H. Araujo
37 // (27th November 2001) 37 // (27th November 2001)
38 // 38 //
39 // PhysicsList program 39 // PhysicsList program
40 // 40 //
41 // Modified: 41 // Modified:
42 // 42 //
43 // 14-02-03 Fix bugs in msc and hIon instancia 43 // 14-02-03 Fix bugs in msc and hIon instanciation + cut per region
44 // 44 //
45 // 05-02-05 AH - changes to G4Decay - added is 45 // 05-02-05 AH - changes to G4Decay - added is not short lived protection
46 // and redefined particles to allow n 46 // and redefined particles to allow non-static creation
47 // i.e. changed construction to G4Mes 47 // i.e. changed construction to G4MesonConstructor, G4BaryonConstructor
48 // 48 //
49 // 23-10-09 LP - migrated EM physics from the 49 // 23-10-09 LP - migrated EM physics from the LowEnergy processes (not supported) to
50 // the new G4Livermore model implemen 50 // the new G4Livermore model implementation. Results unchanged.
51 // 51 //
52 // ------------------------------------------- 52 // --------------------------------------------------------------
53 53
54 #include <iomanip> 54 #include <iomanip>
55 55
56 #include "DMXPhysicsList.hh" 56 #include "DMXPhysicsList.hh"
57 57
58 #include "globals.hh" 58 #include "globals.hh"
59 #include "G4SystemOfUnits.hh" 59 #include "G4SystemOfUnits.hh"
60 #include "G4ProcessManager.hh" 60 #include "G4ProcessManager.hh"
61 #include "G4ProcessVector.hh" 61 #include "G4ProcessVector.hh"
62 62
63 #include "G4ParticleDefinition.hh" 63 #include "G4ParticleDefinition.hh"
64 #include "G4ParticleWithCuts.hh" 64 #include "G4ParticleWithCuts.hh"
65 #include "G4ParticleTypes.hh" 65 #include "G4ParticleTypes.hh"
66 #include "G4ParticleTable.hh" 66 #include "G4ParticleTable.hh"
67 67
68 #include "G4ios.hh" 68 #include "G4ios.hh"
69 #include "G4UserLimits.hh" 69 #include "G4UserLimits.hh"
70 70
71 #include "G4MesonConstructor.hh" <<
72 #include "G4BaryonConstructor.hh" <<
73 #include "G4IonConstructor.hh" <<
74 #include "G4ShortLivedConstructor.hh" <<
75 <<
76 #include "DMXMaxTimeCuts.hh" <<
77 #include "DMXMinEkineCuts.hh" <<
78 #include "G4StepLimiter.hh" <<
79 <<
80 // gamma <<
81 #include "G4PhotoElectricEffect.hh" <<
82 #include "G4LivermorePhotoElectricModel.hh" <<
83 <<
84 #include "G4ComptonScattering.hh" <<
85 #include "G4LivermoreComptonModel.hh" <<
86 <<
87 #include "G4GammaConversion.hh" <<
88 #include "G4BetheHeitler5DModel.hh" <<
89 <<
90 #include "G4RayleighScattering.hh" <<
91 #include "G4LivermoreRayleighModel.hh" <<
92 <<
93 // e- <<
94 #include "G4eMultipleScattering.hh" <<
95 <<
96 #include "G4eIonisation.hh" <<
97 #include "G4LivermoreIonisationModel.hh" <<
98 <<
99 #include "G4eBremsstrahlung.hh" <<
100 #include "G4UniversalFluctuation.hh" <<
101 <<
102 // e+ <<
103 #include "G4eIonisation.hh" <<
104 #include "G4eBremsstrahlung.hh" <<
105 #include "G4eplusAnnihilation.hh" <<
106 <<
107 // alpha and GenericIon and deuterons, triton, <<
108 //muon: <<
109 #include "G4MuIonisation.hh" <<
110 #include "G4MuBremsstrahlung.hh" <<
111 #include "G4MuPairProduction.hh" <<
112 #include "G4MuMultipleScattering.hh" <<
113 #include "G4MuonMinusCapture.hh" <<
114 <<
115 //OTHERS: <<
116 #include "G4hIonisation.hh" <<
117 #include "G4hMultipleScattering.hh" <<
118 #include "G4hBremsstrahlung.hh" <<
119 #include "G4ionIonisation.hh" <<
120 #include "G4IonParametrisedLossModel.hh" <<
121 #include "G4NuclearStopping.hh" <<
122 <<
123 //em process options to allow msc step-limitat <<
124 #include "G4EmParameters.hh" <<
125 #include "G4VAtomDeexcitation.hh" <<
126 #include "G4UAtomicDeexcitation.hh" <<
127 #include "G4LossTableManager.hh" <<
128 <<
129 #include "G4Scintillation.hh" <<
130 #include "G4OpAbsorption.hh" <<
131 #include "G4OpBoundaryProcess.hh" <<
132 #include "G4OpticalParameters.hh" <<
133 <<
134 // Elastic processes: <<
135 #include "G4HadronElasticProcess.hh" <<
136 #include "G4ChipsElasticModel.hh" <<
137 #include "G4ElasticHadrNucleusHE.hh" <<
138 <<
139 // Inelastic processes: <<
140 #include "G4HadronInelasticProcess.hh" <<
141 <<
142 // High energy FTFP model and Bertini cascade <<
143 #include "G4FTFModel.hh" <<
144 #include "G4LundStringFragmentation.hh" <<
145 #include "G4ExcitedStringDecay.hh" <<
146 #include "G4PreCompoundModel.hh" <<
147 #include "G4GeneratorPrecompoundInterface.hh" <<
148 #include "G4TheoFSGenerator.hh" <<
149 #include "G4CascadeInterface.hh" <<
150 <<
151 // Cross sections <<
152 #include "G4VCrossSectionDataSet.hh" <<
153 #include "G4CrossSectionDataSetRegistry.hh" <<
154 <<
155 #include "G4CrossSectionElastic.hh" <<
156 #include "G4CrossSectionInelastic.hh" <<
157 #include "G4BGGPionElasticXS.hh" <<
158 #include "G4BGGPionInelasticXS.hh" <<
159 #include "G4AntiNuclElastic.hh" <<
160 <<
161 #include "G4CrossSectionInelastic.hh" <<
162 #include "G4BGGNucleonInelasticXS.hh" <<
163 #include "G4BGGNucleonElasticXS.hh" <<
164 #include "G4NeutronInelasticXS.hh" <<
165 #include "G4NeutronElasticXS.hh" <<
166 #include "G4ComponentAntiNuclNuclearXS.hh" <<
167 #include "G4ComponentGGNuclNuclXsc.hh" <<
168 #include "G4ComponentGGHadronNucleusXsc.hh" <<
169 <<
170 #include "G4HadronElastic.hh" <<
171 #include "G4NeutronCaptureProcess.hh" <<
172 <<
173 // Neutron high-precision models: <20 MeV <<
174 #include "G4ParticleHPElastic.hh" <<
175 #include "G4ParticleHPElasticData.hh" <<
176 #include "G4ParticleHPCapture.hh" <<
177 #include "G4ParticleHPCaptureData.hh" <<
178 #include "G4ParticleHPInelastic.hh" <<
179 #include "G4ParticleHPInelasticData.hh" <<
180 <<
181 // Stopping processes <<
182 #include "G4HadronStoppingProcess.hh" <<
183 #include "G4HadronicAbsorptionBertini.hh" <<
184 #include "G4HadronicAbsorptionFritiof.hh" <<
185 <<
186 #include "G4HadronicParameters.hh" <<
187 <<
188 #include "G4Decay.hh" <<
189 #include "G4RadioactiveDecay.hh" <<
190 #include "G4PhysicsListHelper.hh" <<
191 #include "G4NuclideTable.hh" <<
192 #include "G4NuclearLevelData.hh" <<
193 <<
194 // Constructor /////////////////////////////// 71 // Constructor /////////////////////////////////////////////////////////////
195 DMXPhysicsList::DMXPhysicsList() : G4VUserPhys 72 DMXPhysicsList::DMXPhysicsList() : G4VUserPhysicsList()
196 { 73 {
197 74
198 defaultCutValue = 1.0*micrometer; // 75 defaultCutValue = 1.0*micrometer; //
199 cutForGamma = defaultCutValue; 76 cutForGamma = defaultCutValue;
200 cutForElectron = 1.0*nanometer; 77 cutForElectron = 1.0*nanometer;
201 cutForPositron = defaultCutValue; 78 cutForPositron = defaultCutValue;
202 79
203 VerboseLevel = 1; 80 VerboseLevel = 1;
204 OpVerbLevel = 0; 81 OpVerbLevel = 0;
205 82
206 //set a finer grid of the physic tables in o <<
207 //former LowEnergy models have 200 bins up t <<
208 G4EmParameters* param = G4EmParameters::Inst <<
209 param->SetMaxEnergy(100*GeV); <<
210 param->SetNumberOfBinsPerDecade(20); <<
211 param->SetMscStepLimitType(fMinimal); <<
212 param->SetFluo(true); <<
213 param->SetPixe(true); <<
214 param->SetAuger(true); <<
215 <<
216 G4EmParameters::Instance()->AddPhysics("Worl <<
217 G4DeexPrecoParameters* deex = G4NuclearLevel <<
218 deex->SetStoreICLevelData(true); <<
219 deex->SetMaxLifeTime(G4NuclideTable::GetInst <<
220 /std::log(2.)); <<
221 SetVerboseLevel(VerboseLevel); 83 SetVerboseLevel(VerboseLevel);
222 } 84 }
223 85
>> 86
224 // Destructor //////////////////////////////// 87 // Destructor //////////////////////////////////////////////////////////////
225 DMXPhysicsList::~DMXPhysicsList() 88 DMXPhysicsList::~DMXPhysicsList()
226 {;} 89 {;}
227 90
>> 91
228 // Construct Particles /////////////////////// 92 // Construct Particles /////////////////////////////////////////////////////
229 void DMXPhysicsList::ConstructParticle() 93 void DMXPhysicsList::ConstructParticle()
230 { 94 {
>> 95
231 // In this method, static member functions s 96 // In this method, static member functions should be called
232 // for all particles which you want to use. 97 // for all particles which you want to use.
233 // This ensures that objects of these partic 98 // This ensures that objects of these particle types will be
234 // created in the program. 99 // created in the program.
235 100
236 ConstructMyBosons(); 101 ConstructMyBosons();
237 ConstructMyLeptons(); 102 ConstructMyLeptons();
238 ConstructMyHadrons(); 103 ConstructMyHadrons();
239 ConstructMyShortLiveds(); 104 ConstructMyShortLiveds();
>> 105
240 } 106 }
241 107
>> 108
242 // construct Bosons:////////////////////////// 109 // construct Bosons://///////////////////////////////////////////////////
243 void DMXPhysicsList::ConstructMyBosons() 110 void DMXPhysicsList::ConstructMyBosons()
244 { 111 {
245 // pseudo-particles 112 // pseudo-particles
246 G4Geantino::GeantinoDefinition(); 113 G4Geantino::GeantinoDefinition();
247 G4ChargedGeantino::ChargedGeantinoDefinition 114 G4ChargedGeantino::ChargedGeantinoDefinition();
248 115
249 // gamma 116 // gamma
250 G4Gamma::GammaDefinition(); 117 G4Gamma::GammaDefinition();
251 118
252 //OpticalPhotons 119 //OpticalPhotons
253 G4OpticalPhoton::OpticalPhotonDefinition(); 120 G4OpticalPhoton::OpticalPhotonDefinition();
254 121
255 } 122 }
256 123
>> 124
257 // construct Leptons:///////////////////////// 125 // construct Leptons://///////////////////////////////////////////////////
258 void DMXPhysicsList::ConstructMyLeptons() 126 void DMXPhysicsList::ConstructMyLeptons()
259 { 127 {
260 // leptons 128 // leptons
261 G4Electron::ElectronDefinition(); 129 G4Electron::ElectronDefinition();
262 G4Positron::PositronDefinition(); 130 G4Positron::PositronDefinition();
263 G4MuonPlus::MuonPlusDefinition(); 131 G4MuonPlus::MuonPlusDefinition();
264 G4MuonMinus::MuonMinusDefinition(); 132 G4MuonMinus::MuonMinusDefinition();
265 133
266 G4NeutrinoE::NeutrinoEDefinition(); 134 G4NeutrinoE::NeutrinoEDefinition();
267 G4AntiNeutrinoE::AntiNeutrinoEDefinition(); 135 G4AntiNeutrinoE::AntiNeutrinoEDefinition();
268 G4NeutrinoMu::NeutrinoMuDefinition(); 136 G4NeutrinoMu::NeutrinoMuDefinition();
269 G4AntiNeutrinoMu::AntiNeutrinoMuDefinition() 137 G4AntiNeutrinoMu::AntiNeutrinoMuDefinition();
270 } 138 }
271 139
>> 140
>> 141 #include "G4MesonConstructor.hh"
>> 142 #include "G4BaryonConstructor.hh"
>> 143 #include "G4IonConstructor.hh"
>> 144
272 // construct Hadrons:///////////////////////// 145 // construct Hadrons://///////////////////////////////////////////////////
273 void DMXPhysicsList::ConstructMyHadrons() 146 void DMXPhysicsList::ConstructMyHadrons()
274 { 147 {
275 // mesons 148 // mesons
276 G4MesonConstructor mConstructor; 149 G4MesonConstructor mConstructor;
277 mConstructor.ConstructParticle(); 150 mConstructor.ConstructParticle();
278 151
279 // baryons 152 // baryons
280 G4BaryonConstructor bConstructor; 153 G4BaryonConstructor bConstructor;
281 bConstructor.ConstructParticle(); 154 bConstructor.ConstructParticle();
282 155
283 // ions 156 // ions
284 G4IonConstructor iConstructor; 157 G4IonConstructor iConstructor;
285 iConstructor.ConstructParticle(); 158 iConstructor.ConstructParticle();
>> 159
286 } 160 }
287 161
>> 162 #include "G4ShortLivedConstructor.hh"
>> 163
288 // construct Shortliveds:///////////////////// 164 // construct Shortliveds://///////////////////////////////////////////////////
289 void DMXPhysicsList::ConstructMyShortLiveds() 165 void DMXPhysicsList::ConstructMyShortLiveds()
290 { 166 {
291 G4ShortLivedConstructor slConstructor; 167 G4ShortLivedConstructor slConstructor;
292 slConstructor.ConstructParticle(); 168 slConstructor.ConstructParticle();
293 } 169 }
294 170
>> 171
>> 172
>> 173
295 // Construct Processes /////////////////////// 174 // Construct Processes //////////////////////////////////////////////////////
296 void DMXPhysicsList::ConstructProcess() 175 void DMXPhysicsList::ConstructProcess()
297 { 176 {
>> 177
298 AddTransportation(); 178 AddTransportation();
299 179
300 ConstructEM(); 180 ConstructEM();
301 181
302 ConstructOp(); 182 ConstructOp();
303 183
304 ConstructHad(); 184 ConstructHad();
305 185
306 ConstructGeneral(); 186 ConstructGeneral();
>> 187
307 } 188 }
308 189
>> 190
309 // Transportation //////////////////////////// 191 // Transportation ///////////////////////////////////////////////////////////
>> 192 #include "DMXMaxTimeCuts.hh"
>> 193 #include "DMXMinEkineCuts.hh"
>> 194 #include "G4StepLimiter.hh"
>> 195
310 void DMXPhysicsList::AddTransportation() { 196 void DMXPhysicsList::AddTransportation() {
311 197
312 G4VUserPhysicsList::AddTransportation(); 198 G4VUserPhysicsList::AddTransportation();
313 199
314 auto particleIterator=GetParticleIterator(); 200 auto particleIterator=GetParticleIterator();
315 particleIterator->reset(); 201 particleIterator->reset();
316 while( (*particleIterator)() ){ 202 while( (*particleIterator)() ){
317 G4ParticleDefinition* particle = particleI 203 G4ParticleDefinition* particle = particleIterator->value();
318 G4ProcessManager* pmanager = particle->Get 204 G4ProcessManager* pmanager = particle->GetProcessManager();
319 G4String particleName = particle->GetParti 205 G4String particleName = particle->GetParticleName();
320 // time cuts for ONLY neutrons: 206 // time cuts for ONLY neutrons:
321 if(particleName == "neutron") 207 if(particleName == "neutron")
322 pmanager->AddDiscreteProcess(new DMXMaxT 208 pmanager->AddDiscreteProcess(new DMXMaxTimeCuts());
323 // Energy cuts to kill charged (embedded i 209 // Energy cuts to kill charged (embedded in method) particles:
324 pmanager->AddDiscreteProcess(new DMXMinEki 210 pmanager->AddDiscreteProcess(new DMXMinEkineCuts());
325 211
326 // Step limit applied to all particles: 212 // Step limit applied to all particles:
327 pmanager->AddDiscreteProcess(new G4StepLim << 213 pmanager->AddProcess(new G4StepLimiter, -1,-1,1);
>> 214
328 } 215 }
329 } 216 }
330 217
>> 218
331 // Electromagnetic Processes ///////////////// 219 // Electromagnetic Processes ////////////////////////////////////////////////
332 // all charged particles 220 // all charged particles
>> 221
>> 222 // gamma
>> 223 #include "G4PhotoElectricEffect.hh"
>> 224 #include "G4LivermorePhotoElectricModel.hh"
>> 225
>> 226 #include "G4ComptonScattering.hh"
>> 227 #include "G4LivermoreComptonModel.hh"
>> 228
>> 229 #include "G4GammaConversion.hh"
>> 230 #include "G4LivermoreGammaConversionModel.hh"
>> 231
>> 232 #include "G4RayleighScattering.hh"
>> 233 #include "G4LivermoreRayleighModel.hh"
>> 234
>> 235
>> 236 // e-
>> 237 #include "G4eMultipleScattering.hh"
>> 238
>> 239 #include "G4eIonisation.hh"
>> 240 #include "G4LivermoreIonisationModel.hh"
>> 241
>> 242 #include "G4eBremsstrahlung.hh"
>> 243 #include "G4LivermoreBremsstrahlungModel.hh"
>> 244
>> 245
>> 246 // e+
>> 247 #include "G4eIonisation.hh"
>> 248 #include "G4eBremsstrahlung.hh"
>> 249 #include "G4eplusAnnihilation.hh"
>> 250
>> 251
>> 252 // alpha and GenericIon and deuterons, triton, He3:
>> 253
>> 254 //muon:
>> 255 #include "G4MuIonisation.hh"
>> 256 #include "G4MuBremsstrahlung.hh"
>> 257 #include "G4MuPairProduction.hh"
>> 258 #include "G4MuonMinusCapture.hh"
>> 259
>> 260 //OTHERS:
>> 261 #include "G4hIonisation.hh"
>> 262 #include "G4hMultipleScattering.hh"
>> 263 #include "G4hBremsstrahlung.hh"
>> 264 #include "G4ionIonisation.hh"
>> 265 #include "G4IonParametrisedLossModel.hh"
>> 266
>> 267 //em process options to allow msc step-limitation to be switched off
>> 268 #include "G4EmParameters.hh"
>> 269 #include "G4VAtomDeexcitation.hh"
>> 270 #include "G4UAtomicDeexcitation.hh"
>> 271 #include "G4LossTableManager.hh"
>> 272
333 void DMXPhysicsList::ConstructEM() { 273 void DMXPhysicsList::ConstructEM() {
334 274
>> 275 //set a finer grid of the physic tables in order to improve precision
>> 276 //former LowEnergy models have 200 bins up to 100 GeV
>> 277 G4EmParameters* param = G4EmParameters::Instance();
>> 278 param->SetMaxEnergy(100*GeV);
>> 279 param->SetNumberOfBinsPerDecade(20);
>> 280 param->SetMscStepLimitType(fMinimal);
>> 281 param->SetFluo(true);
>> 282 param->SetPixe(true);
>> 283 param->SetAuger(true);
335 G4LossTableManager* man = G4LossTableManager 284 G4LossTableManager* man = G4LossTableManager::Instance();
336 man->SetAtomDeexcitation(new G4UAtomicDeexci << 285 G4VAtomDeexcitation* ad = man->AtomDeexcitation();
>> 286 if(!ad) {
>> 287 man->SetAtomDeexcitation(new G4UAtomicDeexcitation());
>> 288 }
337 289
338 G4EmParameters* em_params = G4EmParameters:: <<
339 <<
340 auto particleIterator=GetParticleIterator(); 290 auto particleIterator=GetParticleIterator();
341 particleIterator->reset(); 291 particleIterator->reset();
342 while( (*particleIterator)() ){ 292 while( (*particleIterator)() ){
343 G4ParticleDefinition* particle = particleI 293 G4ParticleDefinition* particle = particleIterator->value();
344 G4ProcessManager* pmanager = particle->Get 294 G4ProcessManager* pmanager = particle->GetProcessManager();
345 G4String particleName = particle->GetParti 295 G4String particleName = particle->GetParticleName();
346 G4String particleType = particle->GetParti 296 G4String particleType = particle->GetParticleType();
347 G4double charge = particle->GetPDGCharge() 297 G4double charge = particle->GetPDGCharge();
348 298
349 if (particleName == "gamma") 299 if (particleName == "gamma")
350 { 300 {
351 //gamma 301 //gamma
352 G4RayleighScattering* theRayleigh = new G4Ra 302 G4RayleighScattering* theRayleigh = new G4RayleighScattering();
353 pmanager->AddDiscreteProcess(theRayleigh); 303 pmanager->AddDiscreteProcess(theRayleigh);
354 304
355 G4PhotoElectricEffect* thePhotoElectricEffec 305 G4PhotoElectricEffect* thePhotoElectricEffect = new G4PhotoElectricEffect();
356 thePhotoElectricEffect->SetEmModel(new G4Liv 306 thePhotoElectricEffect->SetEmModel(new G4LivermorePhotoElectricModel());
357 pmanager->AddDiscreteProcess(thePhotoElectri 307 pmanager->AddDiscreteProcess(thePhotoElectricEffect);
358 308
359 G4ComptonScattering* theComptonScattering = 309 G4ComptonScattering* theComptonScattering = new G4ComptonScattering();
360 theComptonScattering->SetEmModel(new G4Liver 310 theComptonScattering->SetEmModel(new G4LivermoreComptonModel());
361 pmanager->AddDiscreteProcess(theComptonScatt 311 pmanager->AddDiscreteProcess(theComptonScattering);
362 312
363 G4GammaConversion* theGammaConversion = new 313 G4GammaConversion* theGammaConversion = new G4GammaConversion();
364 theGammaConversion->SetEmModel(new G4BetheHe << 314 theGammaConversion->SetEmModel(new G4LivermoreGammaConversionModel());
365 pmanager->AddDiscreteProcess(theGammaConvers 315 pmanager->AddDiscreteProcess(theGammaConversion);
366 316
367 } 317 }
368 else if (particleName == "e-") 318 else if (particleName == "e-")
369 { 319 {
370 //electron 320 //electron
371 // process ordering: AddProcess(name, at res 321 // process ordering: AddProcess(name, at rest, along step, post step)
372 // Multiple scattering 322 // Multiple scattering
373 G4eMultipleScattering* msc = new G4eMultiple 323 G4eMultipleScattering* msc = new G4eMultipleScattering();
374 em_params->SetMscStepLimitType(fUseDistanceT << 324 msc->SetStepLimitType(fUseDistanceToBoundary);
375 pmanager->AddProcess(msc,-1, 1, -1); 325 pmanager->AddProcess(msc,-1, 1, -1);
376 326
377 // Ionisation 327 // Ionisation
378 G4eIonisation* eIonisation = new G4eIonisati 328 G4eIonisation* eIonisation = new G4eIonisation();
379 G4VEmModel* theIoniLiv = new G4Livermo << 329 eIonisation->SetEmModel(new G4LivermoreIonisationModel());
380 theIoniLiv->SetHighEnergyLimit(0.1*MeV << 330 eIonisation->SetStepFunction(0.2, 100*um); //improved precision in tracking
381 eIonisation->AddEmModel(0, theIoniLiv, << 331 pmanager->AddProcess(eIonisation,-1, 2, 2);
382 em_params->SetStepFunction(0.2, 100*um); //i <<
383 pmanager->AddProcess(eIonisation,-1, 2, 1); <<
384 332
385 // Bremsstrahlung 333 // Bremsstrahlung
386 G4eBremsstrahlung* eBremsstrahlung = new G4e 334 G4eBremsstrahlung* eBremsstrahlung = new G4eBremsstrahlung();
387 pmanager->AddProcess(eBremsstrahlung, -1,-3, << 335 eBremsstrahlung->SetEmModel(new G4LivermoreBremsstrahlungModel());
>> 336 pmanager->AddProcess(eBremsstrahlung, -1,-3, 3);
388 } 337 }
389 else if (particleName == "e+") 338 else if (particleName == "e+")
390 { 339 {
391 //positron 340 //positron
392 G4eMultipleScattering* msc = new G4eMultiple 341 G4eMultipleScattering* msc = new G4eMultipleScattering();
393 msc->SetStepLimitType(fUseDistanceToBoundary 342 msc->SetStepLimitType(fUseDistanceToBoundary);
394 pmanager->AddProcess(msc,-1, 1, -1); << 343 pmanager->AddProcess(msc,-1, 1, 1);
395 344
396 // Ionisation 345 // Ionisation
397 G4eIonisation* eIonisation = new G4eIonisati 346 G4eIonisation* eIonisation = new G4eIonisation();
398 // eIonisation->SetStepFunction(0.2, 100*um) << 347 eIonisation->SetStepFunction(0.2, 100*um); //
399 pmanager->AddProcess(eIonisation, << 348 pmanager->AddProcess(eIonisation, -1, 2, 2);
400 349
401 //Bremsstrahlung (use default, no low-energy 350 //Bremsstrahlung (use default, no low-energy available)
402 pmanager->AddProcess(new G4eBremsstrahlung() << 351 pmanager->AddProcess(new G4eBremsstrahlung(), -1,-1, 3);
403 352
404 //Annihilation 353 //Annihilation
405 pmanager->AddProcess(new G4eplusAnnihilation << 354 pmanager->AddProcess(new G4eplusAnnihilation(),0,-1, 4);
406 } 355 }
407 else if( particleName == "mu+" || 356 else if( particleName == "mu+" ||
408 particleName == "mu-" ) 357 particleName == "mu-" )
409 { 358 {
410 //muon 359 //muon
411 pmanager->AddProcess(new G4MuMultipleScatter << 360 pmanager->AddProcess(new G4eMultipleScattering, -1, 1, 1);
412 pmanager->AddProcess(new G4MuIonisation(), << 361 pmanager->AddProcess(new G4MuIonisation(), -1, 2, 2);
413 pmanager->AddProcess(new G4MuBremsstrahlung( << 362 pmanager->AddProcess(new G4MuBremsstrahlung(), -1,-1, 3);
414 pmanager->AddProcess(new G4MuPairProduction( << 363 pmanager->AddProcess(new G4MuPairProduction(), -1,-1, 4);
415 if( particleName == "mu-" ) 364 if( particleName == "mu-" )
416 pmanager->AddProcess(new G4MuonMinusCaptur 365 pmanager->AddProcess(new G4MuonMinusCapture(), 0,-1,-1);
417 } 366 }
418 else if (particleName == "proton" || 367 else if (particleName == "proton" ||
419 particleName == "pi+" || 368 particleName == "pi+" ||
420 particleName == "pi-") 369 particleName == "pi-")
421 { 370 {
422 //multiple scattering 371 //multiple scattering
423 pmanager->AddProcess(new G4hMultipleScatteri << 372 pmanager->AddProcess(new G4hMultipleScattering, -1, 1, 1);
424 373
425 //ionisation 374 //ionisation
426 G4hIonisation* hIonisation = new G4hIonisati 375 G4hIonisation* hIonisation = new G4hIonisation();
427 em_params->SetStepFunctionMuHad(0.2, 50*um); << 376 hIonisation->SetStepFunction(0.2, 50*um);
428 pmanager->AddProcess(hIonisation, << 377 pmanager->AddProcess(hIonisation, -1, 2, 2);
429 378
430 //bremmstrahlung 379 //bremmstrahlung
431 pmanager->AddProcess(new G4hBremsstrahlung, << 380 pmanager->AddProcess(new G4hBremsstrahlung, -1,-3, 3);
432 } 381 }
433 else if(particleName == "alpha" || 382 else if(particleName == "alpha" ||
434 particleName == "deuteron" || 383 particleName == "deuteron" ||
435 particleName == "triton" || 384 particleName == "triton" ||
436 particleName == "He3") 385 particleName == "He3")
437 { 386 {
438 //multiple scattering 387 //multiple scattering
439 pmanager->AddProcess(new G4hMultipleScatteri << 388 pmanager->AddProcess(new G4hMultipleScattering,-1,1,1);
440 389
441 //ionisation 390 //ionisation
442 G4ionIonisation* ionIoni = new G4ionIonisati 391 G4ionIonisation* ionIoni = new G4ionIonisation();
443 em_params->SetStepFunctionLightIons(0.1, 1*C << 392 ionIoni->SetStepFunction(0.1, 20*um);
444 pmanager->AddProcess(ionIoni, << 393 pmanager->AddProcess(ionIoni, -1, 2, 2);
445 pmanager->AddProcess(new G4NuclearStopping() <<
446 } 394 }
447 else if (particleName == "GenericIon") 395 else if (particleName == "GenericIon")
448 { 396 {
449 // OBJECT may be dynamically created as eith 397 // OBJECT may be dynamically created as either a GenericIon or nucleus
450 // G4Nucleus exists and therefore has partic 398 // G4Nucleus exists and therefore has particle type nucleus
451 // genericIon: 399 // genericIon:
452 400
453 //multiple scattering 401 //multiple scattering
454 pmanager->AddProcess(new G4hMultipleScatteri << 402 pmanager->AddProcess(new G4hMultipleScattering,-1,1,1);
455 403
456 //ionisation 404 //ionisation
457 G4ionIonisation* ionIoni = new G4ionIonisati 405 G4ionIonisation* ionIoni = new G4ionIonisation();
458 em_params->SetStepFunctionIons(0.1, 1*CLHEP: << 406 ionIoni->SetEmModel(new G4IonParametrisedLossModel());
459 pmanager->AddProcess(ionIoni, << 407 ionIoni->SetStepFunction(0.1, 20*um);
460 pmanager->AddProcess(new G4NuclearStopping() << 408 pmanager->AddProcess(ionIoni, -1, 2, 2);
461 } 409 }
462 410
463 else if ((!particle->IsShortLived()) && 411 else if ((!particle->IsShortLived()) &&
464 (charge != 0.0) && 412 (charge != 0.0) &&
465 (particle->GetParticleName() != "charge 413 (particle->GetParticleName() != "chargedgeantino"))
466 { 414 {
467 //all others charged particles except geanti 415 //all others charged particles except geantino
468 G4hMultipleScattering* aMultipleScatte 416 G4hMultipleScattering* aMultipleScattering = new G4hMultipleScattering();
469 G4hIonisation* ahadronIon = new G4hIon 417 G4hIonisation* ahadronIon = new G4hIonisation();
470 418
471 //multiple scattering 419 //multiple scattering
472 pmanager->AddProcess(aMultipleScattering,-1, << 420 pmanager->AddProcess(aMultipleScattering,-1,1,1);
473 421
474 //ionisation 422 //ionisation
475 pmanager->AddProcess(ahadronIon, -1,2, << 423 pmanager->AddProcess(ahadronIon, -1,2,2);
476 } 424 }
>> 425
477 } 426 }
478 } 427 }
479 428
>> 429
480 // Optical Processes ///////////////////////// 430 // Optical Processes ////////////////////////////////////////////////////////
>> 431 #include "G4Scintillation.hh"
>> 432 #include "G4OpAbsorption.hh"
>> 433 //#include "G4OpRayleigh.hh"
>> 434 #include "G4OpBoundaryProcess.hh"
>> 435
481 void DMXPhysicsList::ConstructOp() 436 void DMXPhysicsList::ConstructOp()
482 { 437 {
483 G4OpticalParameters* opParams = G4OpticalPar << 438 // default scintillation process
484 G4Scintillation* theScintProcessDef = new G4 439 G4Scintillation* theScintProcessDef = new G4Scintillation("Scintillation");
485 opParams->SetScintTrackSecondariesFirst(true << 440 // theScintProcessDef->DumpPhysicsTable();
486 opParams->SetScintByParticleType(true); << 441 theScintProcessDef->SetTrackSecondariesFirst(true);
>> 442 theScintProcessDef->SetScintillationYieldFactor(1.0); //
>> 443 theScintProcessDef->SetScintillationExcitationRatio(0.0); //
>> 444 theScintProcessDef->SetVerboseLevel(OpVerbLevel);
>> 445
>> 446 // scintillation process for alpha:
>> 447 G4Scintillation* theScintProcessAlpha = new G4Scintillation("Scintillation");
>> 448 // theScintProcessNuc->DumpPhysicsTable();
>> 449 theScintProcessAlpha->SetTrackSecondariesFirst(true);
>> 450 theScintProcessAlpha->SetScintillationYieldFactor(1.1);
>> 451 theScintProcessAlpha->SetScintillationExcitationRatio(1.0);
>> 452 theScintProcessAlpha->SetVerboseLevel(OpVerbLevel);
>> 453
>> 454 // scintillation process for heavy nuclei
>> 455 G4Scintillation* theScintProcessNuc = new G4Scintillation("Scintillation");
>> 456 // theScintProcessNuc->DumpPhysicsTable();
>> 457 theScintProcessNuc->SetTrackSecondariesFirst(true);
>> 458 theScintProcessNuc->SetScintillationYieldFactor(0.2);
>> 459 theScintProcessNuc->SetScintillationExcitationRatio(1.0);
>> 460 theScintProcessNuc->SetVerboseLevel(OpVerbLevel);
487 461
488 // optical processes 462 // optical processes
489 G4OpAbsorption* theAbsorptionProcess = new G 463 G4OpAbsorption* theAbsorptionProcess = new G4OpAbsorption();
>> 464 // G4OpRayleigh* theRayleighScatteringProcess = new G4OpRayleigh();
490 G4OpBoundaryProcess* theBoundaryProcess = ne 465 G4OpBoundaryProcess* theBoundaryProcess = new G4OpBoundaryProcess();
>> 466 // theAbsorptionProcess->DumpPhysicsTable();
>> 467 // theRayleighScatteringProcess->DumpPhysicsTable();
>> 468 theAbsorptionProcess->SetVerboseLevel(OpVerbLevel);
>> 469 // theRayleighScatteringProcess->SetVerboseLevel(OpVerbLevel);
>> 470 theBoundaryProcess->SetVerboseLevel(OpVerbLevel);
491 471
492 auto particleIterator=GetParticleIterator(); 472 auto particleIterator=GetParticleIterator();
493 particleIterator->reset(); 473 particleIterator->reset();
494 while( (*particleIterator)() ) 474 while( (*particleIterator)() )
495 { 475 {
496 G4ParticleDefinition* particle = particl 476 G4ParticleDefinition* particle = particleIterator->value();
497 G4ProcessManager* pmanager = particle->G 477 G4ProcessManager* pmanager = particle->GetProcessManager();
498 G4String particleName = particle->GetPar 478 G4String particleName = particle->GetParticleName();
499 if (theScintProcessDef->IsApplicable(*pa 479 if (theScintProcessDef->IsApplicable(*particle)) {
500 pmanager->AddProcess(theScintProcessDe << 480 // if(particle->GetPDGMass() > 5.0*GeV)
501 pmanager->SetProcessOrderingToLast(the << 481 if(particle->GetParticleName() == "GenericIon") {
502 pmanager->SetProcessOrderingToLast(the << 482 pmanager->AddProcess(theScintProcessNuc); // AtRestDiscrete
>> 483 pmanager->SetProcessOrderingToLast(theScintProcessNuc,idxAtRest);
>> 484 pmanager->SetProcessOrderingToLast(theScintProcessNuc,idxPostStep);
>> 485 }
>> 486 else if(particle->GetParticleName() == "alpha") {
>> 487 pmanager->AddProcess(theScintProcessAlpha);
>> 488 pmanager->SetProcessOrderingToLast(theScintProcessAlpha,idxAtRest);
>> 489 pmanager->SetProcessOrderingToLast(theScintProcessAlpha,idxPostStep);
>> 490 }
>> 491 else {
>> 492 pmanager->AddProcess(theScintProcessDef);
>> 493 pmanager->SetProcessOrderingToLast(theScintProcessDef,idxAtRest);
>> 494 pmanager->SetProcessOrderingToLast(theScintProcessDef,idxPostStep);
>> 495 }
503 } 496 }
504 497
505 if (particleName == "opticalphoton") { 498 if (particleName == "opticalphoton") {
506 pmanager->AddDiscreteProcess(theAbsorptionPr 499 pmanager->AddDiscreteProcess(theAbsorptionProcess);
>> 500 // pmanager->AddDiscreteProcess(theRayleighScatteringProcess);
507 pmanager->AddDiscreteProcess(theBoundaryProc 501 pmanager->AddDiscreteProcess(theBoundaryProcess);
508 } 502 }
509 } 503 }
510 } 504 }
511 505
>> 506
512 // Hadronic processes //////////////////////// 507 // Hadronic processes ////////////////////////////////////////////////////////
513 508
>> 509 // Elastic processes:
>> 510 #include "G4HadronElasticProcess.hh"
>> 511 #include "G4ChipsElasticModel.hh"
>> 512 #include "G4ElasticHadrNucleusHE.hh"
>> 513
>> 514 // Inelastic processes:
>> 515 #include "G4PionPlusInelasticProcess.hh"
>> 516 #include "G4PionMinusInelasticProcess.hh"
>> 517 #include "G4KaonPlusInelasticProcess.hh"
>> 518 #include "G4KaonZeroSInelasticProcess.hh"
>> 519 #include "G4KaonZeroLInelasticProcess.hh"
>> 520 #include "G4KaonMinusInelasticProcess.hh"
>> 521 #include "G4ProtonInelasticProcess.hh"
>> 522 #include "G4AntiProtonInelasticProcess.hh"
>> 523 #include "G4NeutronInelasticProcess.hh"
>> 524 #include "G4AntiNeutronInelasticProcess.hh"
>> 525 #include "G4DeuteronInelasticProcess.hh"
>> 526 #include "G4TritonInelasticProcess.hh"
>> 527 #include "G4AlphaInelasticProcess.hh"
>> 528
>> 529 // High energy FTFP model and Bertini cascade
>> 530 #include "G4FTFModel.hh"
>> 531 #include "G4LundStringFragmentation.hh"
>> 532 #include "G4ExcitedStringDecay.hh"
>> 533 #include "G4PreCompoundModel.hh"
>> 534 #include "G4GeneratorPrecompoundInterface.hh"
>> 535 #include "G4TheoFSGenerator.hh"
>> 536 #include "G4CascadeInterface.hh"
>> 537
>> 538 // Cross sections
>> 539 #include "G4VCrossSectionDataSet.hh"
>> 540 #include "G4CrossSectionDataSetRegistry.hh"
>> 541
>> 542 #include "G4CrossSectionElastic.hh"
>> 543 #include "G4BGGPionElasticXS.hh"
>> 544 #include "G4AntiNuclElastic.hh"
>> 545
>> 546 #include "G4CrossSectionInelastic.hh"
>> 547 #include "G4PiNuclearCrossSection.hh"
>> 548 #include "G4CrossSectionPairGG.hh"
>> 549 #include "G4BGGNucleonInelasticXS.hh"
>> 550 #include "G4ComponentAntiNuclNuclearXS.hh"
>> 551 #include "G4ComponentGGNuclNuclXsc.hh"
>> 552
>> 553 #include "G4HadronElastic.hh"
>> 554 #include "G4HadronCaptureProcess.hh"
>> 555
>> 556 // Neutron high-precision models: <20 MeV
>> 557 #include "G4ParticleHPElastic.hh"
>> 558 #include "G4ParticleHPElasticData.hh"
>> 559 #include "G4ParticleHPCapture.hh"
>> 560 #include "G4ParticleHPCaptureData.hh"
>> 561 #include "G4ParticleHPInelastic.hh"
>> 562 #include "G4ParticleHPInelasticData.hh"
>> 563
>> 564 // Stopping processes
>> 565 #include "G4PiMinusAbsorptionBertini.hh"
>> 566 #include "G4KaonMinusAbsorptionBertini.hh"
>> 567 #include "G4AntiProtonAbsorptionFritiof.hh"
>> 568 #include "G4HadronicParameters.hh"
>> 569
>> 570
514 void DMXPhysicsList::ConstructHad() 571 void DMXPhysicsList::ConstructHad()
515 { 572 {
516 //Elastic models 573 //Elastic models
>> 574 const G4double elastic_elimitPi = 1.0*GeV;
>> 575
517 G4HadronElastic* elastic_lhep0 = new G4Hadro 576 G4HadronElastic* elastic_lhep0 = new G4HadronElastic();
>> 577 G4HadronElastic* elastic_lhep1 = new G4HadronElastic();
>> 578 elastic_lhep1->SetMaxEnergy( elastic_elimitPi );
518 G4ChipsElasticModel* elastic_chip = new G4Ch 579 G4ChipsElasticModel* elastic_chip = new G4ChipsElasticModel();
519 G4ElasticHadrNucleusHE* elastic_he = new G4E 580 G4ElasticHadrNucleusHE* elastic_he = new G4ElasticHadrNucleusHE();
>> 581 elastic_he->SetMinEnergy( elastic_elimitPi );
>> 582
520 583
521 // Inelastic scattering 584 // Inelastic scattering
522 const G4double theFTFMin0 = 0.0*GeV; 585 const G4double theFTFMin0 = 0.0*GeV;
523 const G4double theFTFMin1 = 3.0*GeV; << 586 const G4double theFTFMin1 = 4.0*GeV;
524 const G4double theFTFMax = G4HadronicParamet 587 const G4double theFTFMax = G4HadronicParameters::Instance()->GetMaxEnergy();
525 const G4double theBERTMin0 = 0.0*GeV; 588 const G4double theBERTMin0 = 0.0*GeV;
526 const G4double theBERTMin1 = 19.0*MeV; 589 const G4double theBERTMin1 = 19.0*MeV;
527 const G4double theBERTMax = 6.0*GeV; << 590 const G4double theBERTMax = 5.0*GeV;
528 const G4double theHPMin = 0.0*GeV; 591 const G4double theHPMin = 0.0*GeV;
529 const G4double theHPMax = 20.0*MeV; 592 const G4double theHPMax = 20.0*MeV;
530 593
531 G4FTFModel * theStringModel = new G4FTFModel 594 G4FTFModel * theStringModel = new G4FTFModel;
532 G4ExcitedStringDecay * theStringDecay = new 595 G4ExcitedStringDecay * theStringDecay = new G4ExcitedStringDecay( new G4LundStringFragmentation );
533 theStringModel->SetFragmentationModel( theSt 596 theStringModel->SetFragmentationModel( theStringDecay );
534 G4PreCompoundModel * thePreEquilib = new G4P 597 G4PreCompoundModel * thePreEquilib = new G4PreCompoundModel( new G4ExcitationHandler );
535 G4GeneratorPrecompoundInterface * theCascade 598 G4GeneratorPrecompoundInterface * theCascade = new G4GeneratorPrecompoundInterface( thePreEquilib );
536 599
537 G4TheoFSGenerator * theFTFModel0 = new G4The 600 G4TheoFSGenerator * theFTFModel0 = new G4TheoFSGenerator( "FTFP" );
538 theFTFModel0->SetHighEnergyGenerator( theStr 601 theFTFModel0->SetHighEnergyGenerator( theStringModel );
539 theFTFModel0->SetTransport( theCascade ); 602 theFTFModel0->SetTransport( theCascade );
540 theFTFModel0->SetMinEnergy( theFTFMin0 ); 603 theFTFModel0->SetMinEnergy( theFTFMin0 );
541 theFTFModel0->SetMaxEnergy( theFTFMax ); 604 theFTFModel0->SetMaxEnergy( theFTFMax );
542 605
543 G4TheoFSGenerator * theFTFModel1 = new G4The 606 G4TheoFSGenerator * theFTFModel1 = new G4TheoFSGenerator( "FTFP" );
544 theFTFModel1->SetHighEnergyGenerator( theStr 607 theFTFModel1->SetHighEnergyGenerator( theStringModel );
545 theFTFModel1->SetTransport( theCascade ); 608 theFTFModel1->SetTransport( theCascade );
546 theFTFModel1->SetMinEnergy( theFTFMin1 ); 609 theFTFModel1->SetMinEnergy( theFTFMin1 );
547 theFTFModel1->SetMaxEnergy( theFTFMax ); 610 theFTFModel1->SetMaxEnergy( theFTFMax );
548 611
549 G4CascadeInterface * theBERTModel0 = new G4C 612 G4CascadeInterface * theBERTModel0 = new G4CascadeInterface;
550 theBERTModel0->SetMinEnergy( theBERTMin0 ); 613 theBERTModel0->SetMinEnergy( theBERTMin0 );
551 theBERTModel0->SetMaxEnergy( theBERTMax ); 614 theBERTModel0->SetMaxEnergy( theBERTMax );
552 615
553 G4CascadeInterface * theBERTModel1 = new G4C 616 G4CascadeInterface * theBERTModel1 = new G4CascadeInterface;
554 theBERTModel1->SetMinEnergy( theBERTMin1 ); 617 theBERTModel1->SetMinEnergy( theBERTMin1 );
555 theBERTModel1->SetMaxEnergy( theBERTMax ); 618 theBERTModel1->SetMaxEnergy( theBERTMax );
556 619
>> 620 G4VCrossSectionDataSet * thePiData = new G4CrossSectionPairGG( new G4PiNuclearCrossSection, 91*GeV );
557 G4VCrossSectionDataSet * theAntiNucleonData 621 G4VCrossSectionDataSet * theAntiNucleonData = new G4CrossSectionInelastic( new G4ComponentAntiNuclNuclearXS );
558 G4ComponentGGNuclNuclXsc * ggNuclNuclXsec = 622 G4ComponentGGNuclNuclXsc * ggNuclNuclXsec = new G4ComponentGGNuclNuclXsc();
559 G4VCrossSectionDataSet * theGGNuclNuclData = 623 G4VCrossSectionDataSet * theGGNuclNuclData = new G4CrossSectionInelastic(ggNuclNuclXsec);
560 G4VCrossSectionDataSet * theGGNNEl = new G4C <<
561 G4ComponentGGHadronNucleusXsc * ggHNXsec = n <<
562 G4VCrossSectionDataSet * theGGHNEl = new G4C <<
563 G4VCrossSectionDataSet * theGGHNInel = new G <<
564 624
565 auto particleIterator=GetParticleIterator(); 625 auto particleIterator=GetParticleIterator();
566 particleIterator->reset(); 626 particleIterator->reset();
567 while ((*particleIterator)()) 627 while ((*particleIterator)())
568 { 628 {
569 G4ParticleDefinition* particle = particl 629 G4ParticleDefinition* particle = particleIterator->value();
570 G4ProcessManager* pmanager = particle->G 630 G4ProcessManager* pmanager = particle->GetProcessManager();
571 G4String particleName = particle->GetPar 631 G4String particleName = particle->GetParticleName();
572 632
573 if (particleName == "pi+") 633 if (particleName == "pi+")
574 { 634 {
575 // Elastic scattering 635 // Elastic scattering
576 G4HadronElasticProcess* theElasticPr 636 G4HadronElasticProcess* theElasticProcess = new G4HadronElasticProcess;
577 theElasticProcess->AddDataSet( new G 637 theElasticProcess->AddDataSet( new G4BGGPionElasticXS( particle ) );
>> 638 theElasticProcess->RegisterMe( elastic_lhep1 );
578 theElasticProcess->RegisterMe( elast 639 theElasticProcess->RegisterMe( elastic_he );
579 pmanager->AddDiscreteProcess( theElasticPr 640 pmanager->AddDiscreteProcess( theElasticProcess );
580 //Inelastic scattering 641 //Inelastic scattering
581 G4HadronInelasticProcess* theInelasticProc << 642 G4PionPlusInelasticProcess* theInelasticProcess =
582 new G4HadronInelasticProcess( "inelastic << 643 new G4PionPlusInelasticProcess("inelastic");
583 theInelasticProcess->AddDataSet( new G4BGG << 644 theInelasticProcess->AddDataSet( thePiData );
584 theInelasticProcess->RegisterMe( theFTFMod 645 theInelasticProcess->RegisterMe( theFTFModel1 );
585 theInelasticProcess->RegisterMe( the 646 theInelasticProcess->RegisterMe( theBERTModel0 );
586 pmanager->AddDiscreteProcess( theInelastic 647 pmanager->AddDiscreteProcess( theInelasticProcess );
587 } 648 }
588 649
589 else if (particleName == "pi-") 650 else if (particleName == "pi-")
590 { 651 {
591 // Elastic scattering 652 // Elastic scattering
592 G4HadronElasticProcess* theElasticPr 653 G4HadronElasticProcess* theElasticProcess = new G4HadronElasticProcess;
593 theElasticProcess->AddDataSet( new G 654 theElasticProcess->AddDataSet( new G4BGGPionElasticXS( particle ) );
>> 655 theElasticProcess->RegisterMe( elastic_lhep1 );
594 theElasticProcess->RegisterMe( elast 656 theElasticProcess->RegisterMe( elastic_he );
595 pmanager->AddDiscreteProcess( theElasticPr 657 pmanager->AddDiscreteProcess( theElasticProcess );
596 //Inelastic scattering 658 //Inelastic scattering
597 G4HadronInelasticProcess* theInelasticProc << 659 G4PionMinusInelasticProcess* theInelasticProcess =
598 new G4HadronInelasticProcess( "inelastic << 660 new G4PionMinusInelasticProcess("inelastic");
599 theInelasticProcess->AddDataSet( new G4BGG << 661 theInelasticProcess->AddDataSet( thePiData );
600 theInelasticProcess->RegisterMe( theFTFMod 662 theInelasticProcess->RegisterMe( theFTFModel1 );
601 theInelasticProcess->RegisterMe( the 663 theInelasticProcess->RegisterMe( theBERTModel0 );
602 pmanager->AddDiscreteProcess( theInelastic 664 pmanager->AddDiscreteProcess( theInelasticProcess );
603 //Absorption 665 //Absorption
604 pmanager->AddRestProcess(new G4HadronicAbs << 666 pmanager->AddRestProcess(new G4PiMinusAbsorptionBertini, ordDefault);
605 } 667 }
>> 668
606 else if (particleName == "kaon+") 669 else if (particleName == "kaon+")
607 { 670 {
608 // Elastic scattering 671 // Elastic scattering
609 G4HadronElasticProcess* theElasticPr 672 G4HadronElasticProcess* theElasticProcess = new G4HadronElasticProcess;
610 theElasticProcess->AddDataSet( theGGHNEl ) <<
611 theElasticProcess->RegisterMe( elast 673 theElasticProcess->RegisterMe( elastic_lhep0 );
612 pmanager->AddDiscreteProcess( theElasticPr 674 pmanager->AddDiscreteProcess( theElasticProcess );
613 // Inelastic scattering 675 // Inelastic scattering
614 G4HadronInelasticProcess* theInelasticProc << 676 G4KaonPlusInelasticProcess* theInelasticProcess =
615 new G4HadronInelasticProcess( "inelastic << 677 new G4KaonPlusInelasticProcess("inelastic");
616 theInelasticProcess->AddDataSet( theGGHNIn << 678 theInelasticProcess->AddDataSet( G4CrossSectionDataSetRegistry::Instance()->
>> 679 GetCrossSectionDataSet(G4ChipsKaonPlusInelasticXS::Default_Name()));
617 theInelasticProcess->RegisterMe( theFTFMod 680 theInelasticProcess->RegisterMe( theFTFModel1 );
618 theInelasticProcess->RegisterMe( the 681 theInelasticProcess->RegisterMe( theBERTModel0 );
619 pmanager->AddDiscreteProcess( theInelastic 682 pmanager->AddDiscreteProcess( theInelasticProcess );
620 } << 683 }
>> 684
621 else if (particleName == "kaon0S") 685 else if (particleName == "kaon0S")
622 { 686 {
623 // Elastic scattering 687 // Elastic scattering
624 G4HadronElasticProcess* theElasticPr 688 G4HadronElasticProcess* theElasticProcess = new G4HadronElasticProcess;
625 theElasticProcess->AddDataSet( theGGHNEl ) <<
626 theElasticProcess->RegisterMe( elast 689 theElasticProcess->RegisterMe( elastic_lhep0 );
627 pmanager->AddDiscreteProcess( theElasticPr 690 pmanager->AddDiscreteProcess( theElasticProcess );
628 // Inelastic scattering 691 // Inelastic scattering
629 G4HadronInelasticProcess* theInelasticProc << 692 G4KaonZeroSInelasticProcess* theInelasticProcess =
630 new G4HadronInelasticProcess( "inelastic << 693 new G4KaonZeroSInelasticProcess("inelastic");
631 theInelasticProcess->AddDataSet( theGGHNIn << 694 theInelasticProcess->AddDataSet( G4CrossSectionDataSetRegistry::Instance()->
>> 695 GetCrossSectionDataSet(G4ChipsKaonZeroInelasticXS::Default_Name()));
632 theInelasticProcess->RegisterMe( theFTFMod 696 theInelasticProcess->RegisterMe( theFTFModel1 );
633 theInelasticProcess->RegisterMe( the 697 theInelasticProcess->RegisterMe( theBERTModel0 );
634 pmanager->AddDiscreteProcess( theInelastic 698 pmanager->AddDiscreteProcess( theInelasticProcess );
635 } 699 }
636 700
637 else if (particleName == "kaon0L") 701 else if (particleName == "kaon0L")
638 { 702 {
639 // Elastic scattering 703 // Elastic scattering
640 G4HadronElasticProcess* theElasticPr 704 G4HadronElasticProcess* theElasticProcess = new G4HadronElasticProcess;
641 theElasticProcess->AddDataSet( theGGHNEl ) <<
642 theElasticProcess->RegisterMe( elast 705 theElasticProcess->RegisterMe( elastic_lhep0 );
643 pmanager->AddDiscreteProcess( theElasticPr 706 pmanager->AddDiscreteProcess( theElasticProcess );
644 // Inelastic scattering 707 // Inelastic scattering
645 G4HadronInelasticProcess* theInelasticProc << 708 G4KaonZeroLInelasticProcess* theInelasticProcess =
646 new G4HadronInelasticProcess( "inelastic << 709 new G4KaonZeroLInelasticProcess("inelastic");
647 theInelasticProcess->AddDataSet( theGGHNIn << 710 theInelasticProcess->AddDataSet( G4CrossSectionDataSetRegistry::Instance()->
>> 711 GetCrossSectionDataSet(G4ChipsKaonZeroInelasticXS::Default_Name()));
648 theInelasticProcess->RegisterMe( theFTFMod 712 theInelasticProcess->RegisterMe( theFTFModel1 );
649 theInelasticProcess->RegisterMe( the 713 theInelasticProcess->RegisterMe( theBERTModel0 );
650 pmanager->AddDiscreteProcess( theInelastic 714 pmanager->AddDiscreteProcess( theInelasticProcess );
651 } 715 }
652 716
653 else if (particleName == "kaon-") 717 else if (particleName == "kaon-")
654 { 718 {
655 // Elastic scattering 719 // Elastic scattering
656 G4HadronElasticProcess* theElasticPr 720 G4HadronElasticProcess* theElasticProcess = new G4HadronElasticProcess;
657 theElasticProcess->AddDataSet( theGGHNEl ) <<
658 theElasticProcess->RegisterMe( elast 721 theElasticProcess->RegisterMe( elastic_lhep0 );
659 pmanager->AddDiscreteProcess( theElasticPr 722 pmanager->AddDiscreteProcess( theElasticProcess );
660 // Inelastic scattering 723 // Inelastic scattering
661 G4HadronInelasticProcess* theInelasticProc << 724 G4KaonMinusInelasticProcess* theInelasticProcess =
662 new G4HadronInelasticProcess( "inelastic << 725 new G4KaonMinusInelasticProcess("inelastic");
663 theInelasticProcess->AddDataSet( the << 726 theInelasticProcess->AddDataSet( G4CrossSectionDataSetRegistry::Instance()->
>> 727 GetCrossSectionDataSet(G4ChipsKaonMinusInelasticXS::Default_Name()));
664 theInelasticProcess->RegisterMe( theFTFMod 728 theInelasticProcess->RegisterMe( theFTFModel1 );
665 theInelasticProcess->RegisterMe( the 729 theInelasticProcess->RegisterMe( theBERTModel0 );
666 pmanager->AddDiscreteProcess( theInelastic 730 pmanager->AddDiscreteProcess( theInelasticProcess );
667 pmanager->AddRestProcess(new G4HadronicAbs << 731 pmanager->AddRestProcess(new G4KaonMinusAbsorptionBertini, ordDefault);
668 } 732 }
669 733
670 else if (particleName == "proton") 734 else if (particleName == "proton")
671 { 735 {
672 // Elastic scattering 736 // Elastic scattering
673 G4HadronElasticProcess* theElasticPr 737 G4HadronElasticProcess* theElasticProcess = new G4HadronElasticProcess;
674 theElasticProcess->AddDataSet( new G << 738 theElasticProcess->AddDataSet(G4CrossSectionDataSetRegistry::Instance()->
>> 739 GetCrossSectionDataSet(G4ChipsProtonElasticXS::Default_Name()));
675 theElasticProcess->RegisterMe( elast 740 theElasticProcess->RegisterMe( elastic_chip );
676 pmanager->AddDiscreteProcess( theElasticPr 741 pmanager->AddDiscreteProcess( theElasticProcess );
677 // Inelastic scattering 742 // Inelastic scattering
678 G4HadronInelasticProcess* theInelasticProc << 743 G4ProtonInelasticProcess* theInelasticProcess =
679 new G4HadronInelasticProcess( "inelastic << 744 new G4ProtonInelasticProcess("inelastic");
680 theInelasticProcess->AddDataSet( new G4BGG 745 theInelasticProcess->AddDataSet( new G4BGGNucleonInelasticXS( G4Proton::Proton() ) );
681 theInelasticProcess->RegisterMe( theFTFMod 746 theInelasticProcess->RegisterMe( theFTFModel1 );
682 theInelasticProcess->RegisterMe( the 747 theInelasticProcess->RegisterMe( theBERTModel0 );
683 pmanager->AddDiscreteProcess( theInelastic 748 pmanager->AddDiscreteProcess( theInelasticProcess );
684 } 749 }
685 else if (particleName == "anti_proton") 750 else if (particleName == "anti_proton")
686 { 751 {
687 // Elastic scattering 752 // Elastic scattering
688 const G4double elastic_elimitAntiNuc 753 const G4double elastic_elimitAntiNuc = 100.0*MeV;
689 G4AntiNuclElastic* elastic_anuc = ne 754 G4AntiNuclElastic* elastic_anuc = new G4AntiNuclElastic();
690 elastic_anuc->SetMinEnergy( elastic_ 755 elastic_anuc->SetMinEnergy( elastic_elimitAntiNuc );
691 G4CrossSectionElastic* elastic_anucx 756 G4CrossSectionElastic* elastic_anucxs = new G4CrossSectionElastic( elastic_anuc->GetComponentCrossSection() );
692 G4HadronElastic* elastic_lhep2 = new 757 G4HadronElastic* elastic_lhep2 = new G4HadronElastic();
693 elastic_lhep2->SetMaxEnergy( elastic 758 elastic_lhep2->SetMaxEnergy( elastic_elimitAntiNuc );
694 G4HadronElasticProcess* theElasticPr 759 G4HadronElasticProcess* theElasticProcess = new G4HadronElasticProcess;
695 theElasticProcess->AddDataSet( elast 760 theElasticProcess->AddDataSet( elastic_anucxs );
696 theElasticProcess->RegisterMe( elast 761 theElasticProcess->RegisterMe( elastic_lhep2 );
697 theElasticProcess->RegisterMe( elast 762 theElasticProcess->RegisterMe( elastic_anuc );
698 pmanager->AddDiscreteProcess( theElasticPr 763 pmanager->AddDiscreteProcess( theElasticProcess );
699 // Inelastic scattering 764 // Inelastic scattering
700 G4HadronInelasticProcess* theInelasticProc << 765 G4AntiProtonInelasticProcess* theInelasticProcess =
701 new G4HadronInelasticProcess( "inelastic << 766 new G4AntiProtonInelasticProcess("inelastic");
702 theInelasticProcess->AddDataSet( theAntiNu 767 theInelasticProcess->AddDataSet( theAntiNucleonData );
703 theInelasticProcess->RegisterMe( theFTFMod 768 theInelasticProcess->RegisterMe( theFTFModel0 );
704 pmanager->AddDiscreteProcess( theInelastic 769 pmanager->AddDiscreteProcess( theInelasticProcess );
705 // Absorption 770 // Absorption
706 pmanager->AddRestProcess(new G4HadronicAbs << 771 pmanager->AddRestProcess(new G4AntiProtonAbsorptionFritiof, ordDefault);
707 } 772 }
>> 773
708 else if (particleName == "neutron") { 774 else if (particleName == "neutron") {
709 // elastic scattering 775 // elastic scattering
710 G4HadronElasticProcess* theElasticProcess = 776 G4HadronElasticProcess* theElasticProcess = new G4HadronElasticProcess;
711 theElasticProcess->AddDataSet(new G4Ne << 777 theElasticProcess->AddDataSet(G4CrossSectionDataSetRegistry::Instance()->GetCrossSectionDataSet(G4ChipsNeutronElasticXS::Default_Name()));
712 G4HadronElastic* elastic_neutronChipsM 778 G4HadronElastic* elastic_neutronChipsModel = new G4ChipsElasticModel();
713 elastic_neutronChipsModel->SetMinEnergy( 19. 779 elastic_neutronChipsModel->SetMinEnergy( 19.0*MeV );
714 theElasticProcess->RegisterMe( elastic 780 theElasticProcess->RegisterMe( elastic_neutronChipsModel );
715 G4ParticleHPElastic * theElasticNeutronHP = 781 G4ParticleHPElastic * theElasticNeutronHP = new G4ParticleHPElastic;
716 theElasticNeutronHP->SetMinEnergy( the 782 theElasticNeutronHP->SetMinEnergy( theHPMin );
717 theElasticNeutronHP->SetMaxEnergy( the 783 theElasticNeutronHP->SetMaxEnergy( theHPMax );
718 theElasticProcess->RegisterMe( theElasticNeu 784 theElasticProcess->RegisterMe( theElasticNeutronHP );
719 theElasticProcess->AddDataSet( new G4Particl 785 theElasticProcess->AddDataSet( new G4ParticleHPElasticData );
720 pmanager->AddDiscreteProcess( theElasticProc 786 pmanager->AddDiscreteProcess( theElasticProcess );
721 // inelastic scattering 787 // inelastic scattering
722 G4HadronInelasticProcess* theInelasticProces << 788 G4NeutronInelasticProcess* theInelasticProcess =
723 new G4HadronInelasticProcess( "inelastic", << 789 new G4NeutronInelasticProcess("inelastic");
724 theInelasticProcess->AddDataSet( new G4Neutr << 790 theInelasticProcess->AddDataSet( new G4BGGNucleonInelasticXS( G4Neutron::Neutron() ) );
725 theInelasticProcess->RegisterMe( theFTFModel 791 theInelasticProcess->RegisterMe( theFTFModel1 );
726 theInelasticProcess->RegisterMe( theBE 792 theInelasticProcess->RegisterMe( theBERTModel1 );
727 G4ParticleHPInelastic * theNeutronInelasticH 793 G4ParticleHPInelastic * theNeutronInelasticHPModel = new G4ParticleHPInelastic;
728 theNeutronInelasticHPModel->SetMinEner 794 theNeutronInelasticHPModel->SetMinEnergy( theHPMin );
729 theNeutronInelasticHPModel->SetMaxEner 795 theNeutronInelasticHPModel->SetMaxEnergy( theHPMax );
730 theInelasticProcess->RegisterMe( theNeutronI 796 theInelasticProcess->RegisterMe( theNeutronInelasticHPModel );
731 theInelasticProcess->AddDataSet( new G4Parti 797 theInelasticProcess->AddDataSet( new G4ParticleHPInelasticData );
732 pmanager->AddDiscreteProcess(theInelasticPro 798 pmanager->AddDiscreteProcess(theInelasticProcess);
733 // capture 799 // capture
734 G4NeutronCaptureProcess* theCaptureProcess = << 800 G4HadronCaptureProcess* theCaptureProcess =
735 new G4NeutronCaptureProcess; << 801 new G4HadronCaptureProcess;
736 G4ParticleHPCapture * theLENeutronCaptureMod 802 G4ParticleHPCapture * theLENeutronCaptureModel = new G4ParticleHPCapture;
737 theLENeutronCaptureModel->SetMinEnergy(theHP 803 theLENeutronCaptureModel->SetMinEnergy(theHPMin);
738 theLENeutronCaptureModel->SetMaxEnergy(theHP 804 theLENeutronCaptureModel->SetMaxEnergy(theHPMax);
739 theCaptureProcess->RegisterMe(theLENeutronCa 805 theCaptureProcess->RegisterMe(theLENeutronCaptureModel);
740 theCaptureProcess->AddDataSet( new G4Particl 806 theCaptureProcess->AddDataSet( new G4ParticleHPCaptureData);
741 pmanager->AddDiscreteProcess(theCaptureProce 807 pmanager->AddDiscreteProcess(theCaptureProcess);
>> 808
742 } 809 }
743 else if (particleName == "anti_neutron") 810 else if (particleName == "anti_neutron")
744 { 811 {
745 // Elastic scattering 812 // Elastic scattering
746 G4HadronElasticProcess* theElasticPr 813 G4HadronElasticProcess* theElasticProcess = new G4HadronElasticProcess;
747 theElasticProcess->AddDataSet( theGGHNEl ) <<
748 theElasticProcess->RegisterMe( elast 814 theElasticProcess->RegisterMe( elastic_lhep0 );
749 pmanager->AddDiscreteProcess( theElasticPr 815 pmanager->AddDiscreteProcess( theElasticProcess );
750 // Inelastic scattering (include ann 816 // Inelastic scattering (include annihilation on-fly)
751 G4HadronInelasticProcess* theInelasticProc << 817 G4AntiNeutronInelasticProcess* theInelasticProcess =
752 new G4HadronInelasticProcess( "inelastic << 818 new G4AntiNeutronInelasticProcess("inelastic");
753 theInelasticProcess->AddDataSet( theAntiNu 819 theInelasticProcess->AddDataSet( theAntiNucleonData );
754 theInelasticProcess->RegisterMe( theFTFMod 820 theInelasticProcess->RegisterMe( theFTFModel0 );
755 pmanager->AddDiscreteProcess( theInelastic 821 pmanager->AddDiscreteProcess( theInelasticProcess );
756 } 822 }
>> 823
757 else if (particleName == "deuteron") 824 else if (particleName == "deuteron")
758 { 825 {
759 // Elastic scattering 826 // Elastic scattering
760 G4HadronElasticProcess* theElasticPr 827 G4HadronElasticProcess* theElasticProcess = new G4HadronElasticProcess;
761 theElasticProcess->AddDataSet( theGGNNEl ) <<
762 theElasticProcess->RegisterMe( elast 828 theElasticProcess->RegisterMe( elastic_lhep0 );
763 pmanager->AddDiscreteProcess( theElasticPr 829 pmanager->AddDiscreteProcess( theElasticProcess );
764 // Inelastic scattering 830 // Inelastic scattering
765 G4HadronInelasticProcess* theInelasticProc << 831 G4DeuteronInelasticProcess* theInelasticProcess =
766 new G4HadronInelasticProcess( "inelastic << 832 new G4DeuteronInelasticProcess("inelastic");
767 theInelasticProcess->AddDataSet( theGGNucl 833 theInelasticProcess->AddDataSet( theGGNuclNuclData );
768 theInelasticProcess->RegisterMe( theFTFMod 834 theInelasticProcess->RegisterMe( theFTFModel1 );
769 theInelasticProcess->RegisterMe( the 835 theInelasticProcess->RegisterMe( theBERTModel0 );
770 pmanager->AddDiscreteProcess( theInelastic 836 pmanager->AddDiscreteProcess( theInelasticProcess );
771 } 837 }
>> 838
772 else if (particleName == "triton") 839 else if (particleName == "triton")
773 { 840 {
774 // Elastic scattering 841 // Elastic scattering
775 G4HadronElasticProcess* theElasticPr 842 G4HadronElasticProcess* theElasticProcess = new G4HadronElasticProcess;
776 theElasticProcess->AddDataSet( theGGNNEl ) <<
777 theElasticProcess->RegisterMe( elast 843 theElasticProcess->RegisterMe( elastic_lhep0 );
778 pmanager->AddDiscreteProcess( theElasticPr 844 pmanager->AddDiscreteProcess( theElasticProcess );
779 // Inelastic scattering 845 // Inelastic scattering
780 G4HadronInelasticProcess* theInelasticProc << 846 G4TritonInelasticProcess* theInelasticProcess =
781 new G4HadronInelasticProcess( "inelastic << 847 new G4TritonInelasticProcess("inelastic");
782 theInelasticProcess->AddDataSet( theGGNucl 848 theInelasticProcess->AddDataSet( theGGNuclNuclData );
783 theInelasticProcess->RegisterMe( theFTFMod 849 theInelasticProcess->RegisterMe( theFTFModel1 );
784 theInelasticProcess->RegisterMe( the 850 theInelasticProcess->RegisterMe( theBERTModel0 );
785 pmanager->AddDiscreteProcess( theInelastic 851 pmanager->AddDiscreteProcess( theInelasticProcess );
786 } 852 }
787 else if (particleName == "alpha") 853 else if (particleName == "alpha")
788 { 854 {
789 // Elastic scattering 855 // Elastic scattering
790 G4HadronElasticProcess* theElasticPr 856 G4HadronElasticProcess* theElasticProcess = new G4HadronElasticProcess;
791 theElasticProcess->AddDataSet( theGGNNEl ) <<
792 theElasticProcess->RegisterMe( elast 857 theElasticProcess->RegisterMe( elastic_lhep0 );
793 pmanager->AddDiscreteProcess( theElasticPr 858 pmanager->AddDiscreteProcess( theElasticProcess );
794 // Inelastic scattering 859 // Inelastic scattering
795 G4HadronInelasticProcess* theInelasticProc << 860 G4AlphaInelasticProcess* theInelasticProcess =
796 new G4HadronInelasticProcess( "inelastic << 861 new G4AlphaInelasticProcess("inelastic");
797 theInelasticProcess->AddDataSet( the 862 theInelasticProcess->AddDataSet( theGGNuclNuclData );
798 theInelasticProcess->RegisterMe( theFTFMod 863 theInelasticProcess->RegisterMe( theFTFModel1 );
799 theInelasticProcess->RegisterMe( the 864 theInelasticProcess->RegisterMe( theBERTModel0 );
800 pmanager->AddDiscreteProcess( theInelastic 865 pmanager->AddDiscreteProcess( theInelasticProcess );
801 } 866 }
>> 867
802 } 868 }
803 } 869 }
804 870
>> 871
805 // Decays //////////////////////////////////// 872 // Decays ///////////////////////////////////////////////////////////////////
>> 873 #include "G4Decay.hh"
>> 874 #include "G4RadioactiveDecay.hh"
>> 875 #include "G4IonTable.hh"
>> 876 #include "G4Ions.hh"
>> 877
806 void DMXPhysicsList::ConstructGeneral() { 878 void DMXPhysicsList::ConstructGeneral() {
807 879
808 // Add Decay Process 880 // Add Decay Process
809 G4Decay* theDecayProcess = new G4Decay(); 881 G4Decay* theDecayProcess = new G4Decay();
810 auto particleIterator=GetParticleIterator(); 882 auto particleIterator=GetParticleIterator();
811 particleIterator->reset(); 883 particleIterator->reset();
812 while( (*particleIterator)() ) 884 while( (*particleIterator)() )
813 { 885 {
814 G4ParticleDefinition* particle = particl 886 G4ParticleDefinition* particle = particleIterator->value();
815 G4ProcessManager* pmanager = particle->G 887 G4ProcessManager* pmanager = particle->GetProcessManager();
816 888
817 if (theDecayProcess->IsApplicable(*parti 889 if (theDecayProcess->IsApplicable(*particle) && !particle->IsShortLived())
818 { 890 {
819 pmanager ->AddProcess(theDecayProcess); 891 pmanager ->AddProcess(theDecayProcess);
820 // set ordering for PostStepDoIt and AtRes 892 // set ordering for PostStepDoIt and AtRestDoIt
821 pmanager ->SetProcessOrdering(theDecayProc 893 pmanager ->SetProcessOrdering(theDecayProcess, idxPostStep);
822 pmanager ->SetProcessOrdering(theDecayProc 894 pmanager ->SetProcessOrdering(theDecayProcess, idxAtRest);
823 } 895 }
824 } 896 }
825 897
826 // Declare radioactive decay to the GenericI 898 // Declare radioactive decay to the GenericIon in the IonTable.
827 G4LossTableManager* man = G4LossTableManager << 899 const G4IonTable *theIonTable =
828 G4VAtomDeexcitation* ad = man->AtomDeexcitat << 900 G4ParticleTable::GetParticleTable()->GetIonTable();
829 if(!ad) { << 901 G4RadioactiveDecay *theRadioactiveDecay = new G4RadioactiveDecay();
830 G4EmParameters::Instance()->SetAugerCascad <<
831 ad = new G4UAtomicDeexcitation(); <<
832 man->SetAtomDeexcitation(ad); <<
833 ad->InitialiseAtomicDeexcitation(); <<
834 } <<
835 902
836 G4PhysicsListHelper::GetPhysicsListHelper()- << 903 for (G4int i=0; i<theIonTable->Entries(); i++)
837 RegisterProcess(new G4RadioactiveDecay(), << 904 {
>> 905 G4String particleName = theIonTable->GetParticle(i)->GetParticleName();
>> 906 G4String particleType = theIonTable->GetParticle(i)->GetParticleType();
>> 907
>> 908 if (particleName == "GenericIon")
>> 909 {
>> 910 G4ProcessManager* pmanager =
>> 911 theIonTable->GetParticle(i)->GetProcessManager();
>> 912 pmanager->SetVerboseLevel(VerboseLevel);
>> 913 pmanager ->AddProcess(theRadioactiveDecay);
>> 914 pmanager ->SetProcessOrdering(theRadioactiveDecay, idxPostStep);
>> 915 pmanager ->SetProcessOrdering(theRadioactiveDecay, idxAtRest);
>> 916 }
>> 917 }
838 } 918 }
839 919
840 // Cuts ////////////////////////////////////// 920 // Cuts /////////////////////////////////////////////////////////////////////
841 void DMXPhysicsList::SetCuts() 921 void DMXPhysicsList::SetCuts()
842 { 922 {
843 923
844 if (verboseLevel >1) 924 if (verboseLevel >1)
845 G4cout << "DMXPhysicsList::SetCuts:"; 925 G4cout << "DMXPhysicsList::SetCuts:";
846 926
847 if (verboseLevel>0){ 927 if (verboseLevel>0){
848 G4cout << "DMXPhysicsList::SetCuts:"; 928 G4cout << "DMXPhysicsList::SetCuts:";
849 G4cout << "CutLength : " 929 G4cout << "CutLength : "
850 << G4BestUnit(defaultCutValue,"Length") < 930 << G4BestUnit(defaultCutValue,"Length") << G4endl;
851 } 931 }
852 932
853 //special for low energy physics 933 //special for low energy physics
854 G4double lowlimit=250*eV; 934 G4double lowlimit=250*eV;
855 G4ProductionCutsTable::GetProductionCutsTabl 935 G4ProductionCutsTable::GetProductionCutsTable()->SetEnergyRange(lowlimit,100.*GeV);
856 936
857 // set cut values for gamma at first and for 937 // set cut values for gamma at first and for e- second and next for e+,
858 // because some processes for e+/e- need cut 938 // because some processes for e+/e- need cut values for gamma
859 SetCutValue(cutForGamma, "gamma"); 939 SetCutValue(cutForGamma, "gamma");
860 SetCutValue(cutForElectron, "e-"); 940 SetCutValue(cutForElectron, "e-");
861 SetCutValue(cutForPositron, "e+"); 941 SetCutValue(cutForPositron, "e+");
862 942
863 if (verboseLevel>0) DumpCutValuesTable(); 943 if (verboseLevel>0) DumpCutValuesTable();
864 } 944 }
865 945
866 946