|Message: CSDA Range: ICRU49 vs Ziegler 2000||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 have derived a simulation from the HadronTherapy advanced example. I use exactly the same physicsList as this example, but I have changed the geometry so that monoenergetic pencil beams are incident on a water cube of 40x40x40 cm3
I am running Geant4.8.0p1 and as a simple benchmark I am evaluating the CSDA range for different proton energy in water. For each primary protons that did not had any nuclear interaction (i.e. only electron ionisation) I record its total track length and average it on a large number of primaries. With the default physicsList, the simulation gives the exact same range as the NIST CSDA values (within 0.2 %).
However if I load the Ziegler 2000 stopping power instead of the default ICRU49, I get a large shift in CSDA range (more than 1 mm at 50 MeV and close to 3 cm for a 200 MeV proton beam). The change of stopping powers is made in the same way as the HadronTherapy example. In my macro file I replace the line
by the line
The statistical uncertainty on my simulation of the CSDA Range is lower than 0.2 mm
Can someone tell me if such a large shift (i.e. 3 cm at 200 MeV and 1 mm at 50 MeV)? if not, what could cause this phenomenon?
|Inline Depth:||Outline Depth:||Add message:|