Message: Re: Running IAEA Phase-Spaces in Example N01 Not Logged In (login)
 Next-in-Thread Next-in-Thread
 Next-in-Forum Next-in-Forum

None Re: Running IAEA Phase-Spaces in Example N01 

Forum: Medical Applications
Re: None Running IAEA Phase-Spaces in Example N01 (Landon Clark)
Re: None Re: Running IAEA Phase-Spaces in Example N01 (Landon Clark)
Date: 16 Jul, 2012
From: Miguel Cortes-Giraldo <Miguel Cortes-Giraldo>

Dear Landon,

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

 Add Message Add Message
to: "Re: Running IAEA Phase-Spaces in Example N01"

 Subscribe Subscribe

This site runs SLAC HyperNews version 1.11-slac-98, derived from the original HyperNews


[ Geant 4 Home | Geant 4 HyperNews | Search | Request New Forum | Feedback ]