|Message: Re: Running IAEA Phase-Spaces in Example N01||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)|
Exactly, you were using the reader of IAEA phsp files correctly. Just let me explain with more details the way it works reading independent histories.
When a phase-space file has been recorded in IAEA format including the variable called "n_stat", then one event of the simulation is corresponded with one independent history of the original simulation which recorded such a phase-space file (a.k.a. original histories).
"n_stat" is an extra integer variable which stores, for each particle, the number of original histories employed to get such a particle after that stored previously in the phsp file. It is a variable often used in EGSnrc and Penelope.
Thus, when "n_stat" is defined, the number of original histories read from a phsp file can be clearly identified with the number of events of the simulation which reads the phsp file. So, if an original history did not hit the phsp plane, it is reproduced as a "blank" event in the Geant4 simulation in order to keep the correspondence with original histories.
Otherwise, if "n_stat" is not defined, each independent history of the IAEAphsp file is just simulated in separated events, with no blank events, since the book-keeping of original histories cannot be performed.
I hope this helps to understand this behaviour. You may also take a look at this recent paper, which explains the main features of these classes, or contact me as well:
MA Cortes-Giraldo et al., Int. J. Rad. Biol. 88: 200-208 (2012) [doi:10.3109/09553002.2011.627977]
Best regads, Miguel