|Message: looking for a Variance Reduction techique to use with Pencil Beam simulations||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 need to find an adequate VR technique to decrease the CPU time of Pencil Beam simulations. The Pencil Beam model couldn't be simpler. It consist of a 6MeV mono-energetic photon point source placed on the top surface of a water cylinder. The goal is to score dose2water in the water cylinder.
I tried to use Topas VR technique named "Secondary Biasing. It expects, among the other parameters,the maximum energy, and the name for each process generating the particles to be split.
I wish to split ONLY "primary" 6MeV photons into 100 fictitious photons. Therefore I set the MaximumEnergy to 6.18 MeV (nomimal Linac energy). Unluckily this VR technique only filters by an upper energy bound. Thus all photons whose energy is lower than the chosen MaximumEnergy will be split producing replicas of low energy photons which we do not care to track.
Furthermore, this VR technique expects the names of the processes that generate the particles to be split (gammas). The problem is to single out the true "primary" photons. In fact the PencilBeam model does NOT include an electron source generating the photons. Geant names the generating process for such 6MeV photons "Primary". To prove it I have added to the PencilBeam model an N-Tuple scorer recording all gammas that hit the water cylinder upper surface, just below the photon source. I am attaching the .header and the .phsp files where you can see (use Excel or MatLab) that most photons are regarded as "Primary" even if their energy is by far lower than 6MeV. The reason is that Geant4 considers photons produced by Compton events as "Primary" regardless of the number of Compton interactions they have had.
Is there an appropriate VR technique for the PencilBeam simulations ? In general, is there a VR technique to use whenever the particle source generates gammas rather than electrons ?
Any suggestion / comment is greatly appreciated. Thank you in advance. Best regards, Mauede
Attachment: http://hypernews.slac.stanford.edu/HyperNews/geant4/get/AUX/2016/10/24/13.23-58559-ROOTOutput.header http://hypernews.slac.stanford.edu/HyperNews/geant4/get/AUX/2016/10/24/13.23-82529-ROOTOutput.phsp