|Message: Re: Running IAEA phase space files as particle generator||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)|
What you have observed is the expected behaviour of our IAEAphsp tool for Geant4. Let me explain you why you observe those apparently wrong results:
1) When you write an IAEAphsp file, it keeps track of what happened at each single event of the original simulation - this is the reason why you need to simulate 10^7 events in order to read *all* the phsp file (since it corresponds to 10^7 original histories). A short example: imagine you have generated a phsp from 100 original histories, and you have only scored a photon at events #44 and #76. Then, if you read back that phsp file in another Geant4 application, what you will get when applying /run/beamOn 100 (in this particular case), you will get a bunch of 'null' events, except for #44 and #76, which correspond to the events when the photons were scored.
The reason of doing this was simply to always keep track of the number of original histories you are actually shooting in a simulation (so that eventually normalization to dose per primary particle can always be done). The drawback is that you will get many 'blank' events for the sake of reproducing exactly what happened when the phsp was written. I know that other codes (such as EGSnrc) did not followed this approach in the past (I do not know if they do it now), so I understand your question.
2) As for recycling, you must know that *all* copies of the same particle are produced *at the same event*. The reason of doing this is because it is the only way that one can keep statistical correlations between particles produced by the same original history (at it applies for recycling). Simultaneously, the statistical weight of the track is divided by the number of repetitions set, so that dosimetry quantities can be normalized accordingly.
To explain your problem, I assume you are scoring the energy deposited at a certain volume, but you are not multiplying by the track's statistical weight. If you do so, you should get again 17 keV. If you are scoring a different quantity, then I will need to know more details about it.
I hope it helps,