| Geant4 Cross Reference |
1 // 1
2 // *******************************************
3 // * License and Disclaimer
4 // *
5 // * The Geant4 software is copyright of th
6 // * the Geant4 Collaboration. It is provided
7 // * conditions of the Geant4 Software License
8 // * LICENSE and available at http://cern.ch/
9 // * include a list of copyright holders.
10 // *
11 // * Neither the authors of this software syst
12 // * institutes,nor the agencies providing fin
13 // * work make any representation or warran
14 // * regarding this software system or assum
15 // * use. Please see the license in the file
16 // * for the full disclaimer and the limitatio
17 // *
18 // * This code implementation is the result
19 // * technical work of the GEANT4 collaboratio
20 // * By using, copying, modifying or distri
21 // * any work based on the software) you ag
22 // * use in resulting scientific publicati
23 // * acceptance of all terms of the Geant4 Sof
24 // *******************************************
25 //
26 //
27 //
28 // Authors: Elena Guardincerri (Elena.Guardinc
29 // Alfonso Mantero (Alfonso.Mantero@g
30 //
31 // History:
32 // -----------
33 //
34 // 16 Sept 2001 First committed to cvs
35 // 12 Sep 2003 Bug in auger production fixe
36 //
37 // -------------------------------------------
38
39 #include "G4RDAtomicDeexcitation.hh"
40 #include "Randomize.hh"
41 #include "G4PhysicalConstants.hh"
42 #include "G4SystemOfUnits.hh"
43 #include "G4Gamma.hh"
44 #include "G4Electron.hh"
45 #include "G4RDAtomicTransitionManager.hh"
46 #include "G4RDFluoTransition.hh"
47
48 G4RDAtomicDeexcitation::G4RDAtomicDeexcitation
49 minGammaEnergy(100.*eV),
50 minElectronEnergy(100.*eV),
51 fAuger(false)
52 {}
53
54 G4RDAtomicDeexcitation::~G4RDAtomicDeexcitatio
55 {}
56
57 std::vector<G4DynamicParticle*>* G4RDAtomicDee
58 {
59
60 std::vector<G4DynamicParticle*>* vectorOfPar
61
62 vectorOfParticles = new std::vector<G4Dynami
63 G4DynamicParticle* aParticle;
64 G4int provShellId = 0;
65 G4int counter = 0;
66
67 // The aim of this loop is to generate more
68 // from the same ionizing event
69 do
70 {
71 if (counter == 0)
72 // First call to GenerateParticles(...):
73 // givenShellId is given by the process
74 {
75 provShellId = SelectTypeOfTransition(Z, gi
76 //std::cout << "AtomicDeexcitation::Genera
77 //<< provShellId << std::endl;
78
79 if ( provShellId >0)
80 {
81 aParticle = GenerateFluorescence(Z,giv
82 //std::cout << "AtomicDeexcitation::Ge
83 }
84 else if ( provShellId == -1)
85 {
86 aParticle = GenerateAuger(Z, givenShel
87 //std::cout << "AtomicDeexcitation::Ge
88 }
89 else
90 {
91 G4Exception("G4RDAtomicDeexcitation::G
92 "InvalidSetup", Fata
93 "Starting shell unco
94 }
95 }
96 else
97 // Following calls to GenerateParticles(...)
98 // newShellId is given by GenerateFluorescen
99 {
100 provShellId = SelectTypeOfTransition(Z,new
101 //std::cout << "AtomicDeexcitation::Genera
102 //<< provShellId << ", new ShellId = "<< n
103 //<< std::endl;
104
105
106 if (provShellId >0)
107 {
108 aParticle = GenerateFluorescence(Z,new
109 //std::cout << "AtomicDeexcitation::Ge
110 }
111 else if ( provShellId == -1)
112 {
113 aParticle = GenerateAuger(Z, newShellI
114 //std::cout << "AtomicDeexcitation::Ge
115 }
116 else
117 {
118 G4Exception("G4RDAtomicDeexcitation::G
119 "InvalidSetup", Fata
120 "Starting shell unco
121 }
122 }
123 counter++;
124 if (aParticle != 0) {vectorOfParticles->
125 else {provShellId = -2;}
126 }
127
128 // Look this in a particular way: only one a
129 while (provShellId > -2);
130
131 return vectorOfParticles;
132 }
133
134 G4int G4RDAtomicDeexcitation::SelectTypeOfTran
135 {
136 if (shellId <=0 )
137 {
138 G4Exception("G4RDAtomicDeexcitation::Sel
139 "InvalidCondition", FatalExc
140 "Zero or negative shellId!")
141 }
142
143 const G4RDAtomicTransitionManager* transiti
144 G4RDAtomicTransitionManager::Instance(
145 G4int provShellId = -1;
146 G4int shellNum = 0;
147 G4int maxNumOfShells = transitionManager->Nu
148
149 //std::cout << "AtomicDeexcitation::SelectTy
150 //<< maxNumOfShells<< std::endl;
151
152 const G4RDFluoTransition* refShell = transit
153
154 // This loop gives shellNum the value of the
155 // in the vector storing the list of the she
156 // a radiative transition
157 if ( shellId <= refShell->FinalShellId())
158 {
159 while (shellId != transitionManager->Rea
160 {
161 if(shellNum ==maxNumOfShells-1)
162 {
163 break;
164 }
165 shellNum++;
166 }
167 G4int transProb = 0; //AM change 29/6/07
168
169 G4double partialProb = G4UniformRand();
170 G4double partSum = 0;
171 const G4RDFluoTransition* aShell = trans
172 G4int trSize = (aShell->TransitionProba
173
174 // Loop over the shells wich can provide
175 // radiative transition towards shellId:
176 // in every loop the partial sum of the
177 // is calculated and compared with a ran
178 // If the partial sum is greater, the sh
179 // is chosen as the starting shell for a
180 // and its identity is returned
181 // Else, terminateded the loop, -1 is re
182 while(transProb < trSize){
183
184 partSum += aShell->TransitionProbability(tr
185
186 if(partialProb <= partSum)
187 {
188 provShellId = aShell->OriginatingShellI
189 break;
190 }
191 transProb++;
192 }
193
194 // here provShellId is the right one or
195 // if -1, the control is passed to the A
196 }
197
198
199
200 else
201 {
202
203 provShellId = -1;
204
205 }
206 return provShellId;
207 }
208
209 G4DynamicParticle* G4RDAtomicDeexcitation::Gen
210 G4int shellId,
211 G4int provShellId )
212 {
213
214
215 const G4RDAtomicTransitionManager* transiti
216 // G4int provenienceShell = provShellId;
217
218 //isotropic angular distribution for the out
219 G4double newcosTh = 1.-2.*G4UniformRand();
220 G4double newsinTh = std::sqrt(1.-newcosTh*n
221 G4double newPhi = twopi*G4UniformRand();
222
223 G4double xDir = newsinTh*std::sin(newPhi);
224 G4double yDir = newsinTh*std::cos(newPhi);
225 G4double zDir = newcosTh;
226
227 G4ThreeVector newGammaDirection(xDir,yDir,zD
228
229 G4int shellNum = 0;
230 G4int maxNumOfShells = transitionManager->Nu
231
232 // find the index of the shell named shellId
233 while (shellId != transitionManager->
234 ReachableShell(Z,shellNum)->FinalShellId())
235 {
236 if(shellNum == maxNumOfShells-1)
237 {
238 break;
239 }
240 shellNum++;
241 }
242 // number of shell from wich an electron can
243 size_t transitionSize = transitionManager->
244 ReachableShell(Z,shellNum)->OriginatingShe
245
246 size_t index = 0;
247
248 // find the index of the shell named provShe
249 // storing the shells from which shellId can
250 while (provShellId != transitionManager->
251 ReachableShell(Z,shellNum)->OriginatingShel
252 {
253 if(index == transitionSize-1)
254 {
255 break;
256 }
257 index++;
258 }
259 // energy of the gamma leaving provShellId f
260 G4double transitionEnergy = transitionManage
261 ReachableShell(Z,shellNum)->TransitionEner
262
263 // This is the shell where the new vacancy i
264 // shell where the electron came from
265 newShellId = transitionManager->
266 ReachableShell(Z,shellNum)->OriginatingShe
267
268
269 G4DynamicParticle* newPart = new G4DynamicPa
270 newGammaDirection,
271 transitionEnergy);
272 return newPart;
273 }
274
275 G4DynamicParticle* G4RDAtomicDeexcitation::Gen
276 {
277 if(!fAuger) return 0;
278
279
280 const G4RDAtomicTransitionManager* transiti
281 G4RDAtomicTransitionManager::Instance(
282
283
284
285 if (shellId <=0 )
286 {
287 G4Exception("G4RDAtomicDeexcitation::Gen
288 "InvalidCondition", FatalExc
289 "Zero or negative shellId!")
290 }
291
292 // G4int provShellId = -1;
293 G4int maxNumOfShells = transitionManager->Nu
294
295 const G4RDAugerTransition* refAugerTransitio
296 transitionManager->ReachableAugerShell
297
298
299 // This loop gives to shellNum the value of
300 // in the vector storing the list of the vac
301 // that can originate a NON-radiative transi
302
303 // ---- MGP ---- Next line commented out to
304 // G4int p = refAugerTransition->FinalShellI
305
306 G4int shellNum = 0;
307
308
309 if ( shellId <= refAugerTransition->FinalShe
310 //"FinalShellId" is final from the point o
311 // being the Id of the shell in which ther
312 {
313 G4int pippo = transitionManager->Reachab
314 if (shellId != pippo ) {
315 do {
316 shellNum++;
317 if(shellNum == maxNumOfShells)
318 {
319 // G4cout << "G4RDAtomicDeexcitation
320 // G4cout << "Absorbed enrgy deposited
321 return 0;
322 // // G4Exception("G4RDAtomicDeexcitat
323 }
324 }
325 while (shellId != (transitionManager->Reacha
326 }
327 /* {
328
329 if(shellNum == maxNumOfShells-1)
330 {
331 G4Exception("G4RDAtomicDeexcitation: N
332 }
333 shellNum++;
334 }*/
335
336
337
338
339 // Now we have that shellnum is the shel
340
341 // G4cout << " the index of the she
342
343 // But we have now to select two shells:
344 // and another for the auger emission.
345
346 G4int transitionLoopShellIndex = 0;
347 G4double partSum = 0;
348 const G4RDAugerTransition* anAugerTransi
349 transitionManager->ReachableAugerS
350
351 // G4cout << " corresponding to the
352
353
354 G4int transitionSize =
355 (anAugerTransition->TransitionOrig
356 while (transitionLoopShellIndex < transi
357
358 std::vector<G4int>::const_iterator pos
359 anAugerTransition->TransitionOr
360
361 G4int transitionLoopShellId = *(pos+tr
362 G4int numberOfPossibleAuger =
363 (anAugerTransition->AugerTransit
364 G4int augerIndex = 0;
365 // G4int partSum2 = 0;
366
367
368 if (augerIndex < numberOfPossibleAuger) {
369
370 do
371 {
372 G4double thisProb = anAugerTransition-
373 transitionLoopShellId);
374 partSum += thisProb;
375 augerIndex++;
376
377 } while (augerIndex < numberOfPossibleAu
378 }
379 transitionLoopShellIndex++;
380 }
381
382
383
384 // Now we have the entire probability of
385 // located in shellNum (index of shellId
386
387 // AM *********************** F I X E D
388 // Here we duplicate the previous loop,
389 // to be under the random number shoot b
390 // previuos loop, while integrating the
391 // 5 minutes from now: a line:
392 // G4int numberOfPossibleAuger = (anAuge
393 // AugerTransitionProbabilities(transiti
394 // to be inserted.
395 // AM *********************** F I X E D
396
397 // Remains to get the same result with a
398
399 // AM *********************** F I X E D
400 // Another Bug: in EADL Auger Transition
401 // a vacancy in one shell, but not all o
402 // doesn't occur in the main one a local
403 // generating the last transition presen
404 // AM *********************** F I X E D
405
406
407 G4double totalVacancyAugerProbability =
408
409
410 //And now we start to select the right a
411 G4int transitionRandomShellIndex = 0;
412 G4int transitionRandomShellId = 1;
413 G4int augerIndex = 0;
414 partSum = 0;
415 G4double partialProb = G4UniformRand();
416 // G4int augerOriginatingShellId = 0;
417
418 G4int numberOfPossibleAuger =
419 (anAugerTransition->AugerTransitionProbabi
420 G4bool foundFlag = false;
421
422 while (transitionRandomShellIndex < tran
423
424 std::vector<G4int>::const_iterator pos
425 anAugerTransition->TransitionOr
426
427 transitionRandomShellId = *(pos+transi
428
429 augerIndex = 0;
430 numberOfPossibleAuger = (anAugerTransition->
431 AugerTransitionProbabilities(transiti
432
433 while (augerIndex < numberOfPossibleAu
434 G4double thisProb =anAugerTransition->Auge
435 transitionRandomShellId);
436
437 partSum += thisProb;
438
439 if (partSum >= (partialProb*totalVac
440 foundFlag = true;
441 break;
442 }
443 augerIndex++;
444 }
445 if (partSum >= (partialProb*totalVacan
446 transitionRandomShellIndex++;
447 }
448
449 // Now we have the index of the shell fr
450 // and the id of the shell, from which t
451 // If no Transition has been found, 0 is
452
453 if (!foundFlag) {return 0;}
454
455 // Isotropic angular distribution for th
456 G4double newcosTh = 1.-2.*G4UniformRand(
457 G4double newsinTh = std::sqrt(1.-newcos
458 G4double newPhi = twopi*G4UniformRand();
459
460 G4double xDir = newsinTh*std::sin(newPh
461 G4double yDir = newsinTh*std::cos(newPhi
462 G4double zDir = newcosTh;
463
464 G4ThreeVector newElectronDirection(xDir,
465
466 // energy of the auger electron emitted
467
468
469 G4double transitionEnergy = anAugerTrans
470 /*
471 G4cout << "AUger TransitionId " << anAugerTr
472 G4cout << "augerIndex: " << augerIndex << G4
473 G4cout << "transitionShellId: " << transitio
474 */
475
476 // This is the shell where the new vacan
477 // shell where the electron came from
478 newShellId = transitionRandomShellId;
479
480
481 G4DynamicParticle* newPart = new G4Dynam
482 newElectronDirection,
483 transitionEnergy);
484 return newPart;
485
486 }
487 else
488 {
489 //G4Exception("G4RDAtomicDeexcitation: n
490 return 0;
491 }
492
493 }
494
495 void G4RDAtomicDeexcitation::SetCutForSecondar
496 {
497 minGammaEnergy = cut;
498 }
499
500 void G4RDAtomicDeexcitation::SetCutForAugerEle
501 {
502 minElectronEnergy = cut;
503 }
504
505 void G4RDAtomicDeexcitation::ActivateAugerElec
506 {
507 fAuger = val;
508 }
509
510
511
512
513
514
515
516