|Message: 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)|
Hello Geant4 community,
I am currently simulating a mammography x-ray tube, and to aid in computational efficiency downstream I am collecting phase space files at different locations. To test my writing and reading of those files, I'm running a 10^7 original event simulation whereby I'm collecting about 10^3 photons at my phase space plane. However, when I read the file back in as the particle generator I get, what I would consider, some unexpected results.
To read the phase space file back in I'm using the command /run/beamOn, but the number of events I have to use to get 10^3 photons out from the file is 10^7 - i.e. /run/beamOn 10000000 when reading in the phase space file results in all 10^3 from the phase space file being loaded in. This is not what I would expect, what I would expect is when using a phase space file with 10^3 photons a command of /run/beamOn 1000, and NOT /run/beamOn 10000000, would load the entire phase space file, but it doesn't.
So my question is, is this typical behavior when reading in phase space files as particle generators? or am I missing something totally obvious? or is there some other command I need to be using with a phase space file?
Another somewhat related issue, when I recycle photons using theIAEAReader->SetTimesRecycled(9) all it does is multiply the energy of the particle from the phase space file by the number of times it was recycled. i.e. if I had SetTimesRecycled(9) to recycle 10 times and I were to detect the particles using a detection volume, I would see one entry of 170keV from a 17keV phase space photon.
Thankyou in advance to anyone that may be able to help, this is one of the last hurdles I need to get over to start collecting data.
P.S. I haven't posted any code as it is essentially the same as that found in the IAEA "HowTo" guide.
|Inline Depth:||Outline Depth:||Add message:|