Hello all! I use G4 under windows and I try to simulate space resolution of a gas detector filled with Xe for gammas with energy 20keV. When I start simutation without Auger electrons all works stable, but when I add the line in phisics list: theLEPhotoElectric->ActivateAuger(true) all crashed!!! The crash happens in: G4AtomicDeexitation.cc(line 395): G4int numberOfPossibleAuger = (anAugerTransition->AugerTransitionProbabilities(transitionRandomShellId))->size (); More interesting the fact, that value of "numberOfPossibleAuger" variable is never used!!! and it change the value just immediately in the cycle. I changed the code, to simple create the variable (see below) and all start to work stable!!! G4int numberOfPossibleAuger; // = (anAugerTransition->AugerTransitionProbabilities(transitionRandomShellId))->size(); Please fix this bug.... Additionally.. in the file g4atomicdeexcitation.cc( line107) : warning C4701: potentially uninitialized local variable 'aParticle' used.... Please initiate the aParticle variable too... Sincerely Yours, Viacheslav Porosev
The problem report is reassigned to the developer responsible for atomic relaxation.
I cannot reproduce this problem. Can you post some informations about your system? Thanx, alf
Created attachment 11 [details] printscreens while crash...
Dear colleague... Recently, I rebuild G4.9.1 and problem when Auger activated still presents at the same place.. I made an attachment, where computer screen while crash is shown ... bug-1.jpg illustrates it. G4 crash when it try to GenerateAuger for Z=36 and shellID=1 and it could't to find correct "shellID" in "AugerTransitionProbabilities" method... see details on bug-2.jpg. The previously suggested method lets to fix the problem…. Sincerely Yours, Viacheslav Porosev
Hello all! I would like to say that G492_patch01 suffer from the same problem in the same place. The "AugerTransitionProbabilities(transitionRandomShellId)" could not get shellId when transitionRandomShellId is equals to 1. The previuosly mentioned solution helps!! Sincerely Yours, Viacheslav Porosev G4DynamicParticle* G4AtomicDeexcitation::GenerateAuger(G4int Z, G4int shellId) { ... Lines 390 etc... G4int transitionRandomShellId = 1; G4int augerIndex = 0; partSum = 0; G4double partialProb = G4UniformRand(); // G4int augerOriginatingShellId = 0; G4int numberOfPossibleAuger;// = // (anAugerTransition->AugerTransitionProbabilities(transitionRandomShellId))->size();
Please use Geant4 version 9.4
