|Message: cannot reproduce rare events with Random Status saved||Not Logged In (login)|
Click on the Forum title, e.g. on the "Forums by Category" page, to read a sequence of postings to the Forum and its threads all in one page. If you are only interested in one thread or the thread following a specific posting, click the thread or the posting, which takes you to a smaller page, which contains only the part you are interested in and may be easier to navigate.
Messages are "chained" if there are only replies at the first level, i.e. 1/1.html, 1/1/1.html etc. In case of "chained" messages the message number is replaced by the icon and there is no indentation.
Inline: Display the subject line only or also the text of the posting(s); for the choice "All" the "Outline" choices are switched off.
|1||0||1||no text / full text of posting|
|2||1||All||text for level 1 only / text for All postings|
Outline: Choose the depth of the posting thread, successive toggle controls provide increasing detail.
|1||2||1||2 levels / 1 level (original posting)|
|2||3||2||3 levels / 2 levels|
|3||3||All||3 levels / all levels (all postings)|
I have a problem when trying to replay selected rare events by Random Status saving/restoring. My physics list consists of default G4DecayPhysics, G4RadioactiveDecayPhysics and G4EmPenelopePhysics. Radioactive decay of a radioisotope is used as Primary Generator (and random numbers are directly invoked at this stage for the position selection). I'm interested in events when emitted electrons reach detector with a specified energy. Two schemes of Random Status storing were tried (both with no success):
1. Simple and usual one: add
in main.cc and then in Stepping Action call:
2. A little bit more complicated (as I wasn't sure that the previous one saves Random Status before Primory Generator is invoked): add
in main.cc and then in Stepping Action execute:
G4String rnd = G4RunManager::GetRunManager()->GetCurrentEvent()->GetRandomNumberStatus();
std::ofstream of; of.open(filename); of << rnd.substr(rnd.find("\n")+1,rnd.npos); of.close();
Doesn't matter which of the above methods was utilized, it didn't provide the selected event reconstruction (after UI command /random/restoreEngineFrom filename). That was clear since no electron penetrated the detector.
What is the reason for such behavior?
Thanks in advance, Ilya Kudinov
P.S. Geant4 version is 9.6, the latest release.
|Inline Depth:||Outline Depth:||Add message:|