|Message: Re: Changes for heavy ions in GEmStandardPhysics.cc in G4.9.2.p02?||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)|
On Wed, 28 Oct 2009 07:10:21 GMT, michel maire wrote:
> On Tue, 27 Oct 2009 10:48:06 GMT, Daniel Matthi� wrote: > > > Hello all, I have a problem after switching from G4.9.1.p02 to > > G4.9.2.p02 in the transport calcultions of heavy ions. I recognized a > > dramatic change in the ranges of iron ions at high energies. > > Specifically I calculated the range for Fe56 at an energy of 500MeV/nuc > > (28GeV) using the G4EmStandardPhysics list. The CSDA range in water in > > this case is 97.4mm calculated with G4.9.1.p02 which is reasonable but > > 204.3mm calculated with G4.9.2.p02! What is the reason for this change? > > Do I have to use a different physics list? In the creation of my physics > > list I followed the example extended/field/field04. Regards, Daniel > > > > I have checked Fe56 range with TestEm0 and TestEm7. It is still around 96 mm > (extended/electromagnetic) > > Michel >
I found out that the problem is not in the physics list but in the material. I used the G4NistManager to obtain the definition for water:
G4NistManager *nistMan = G4NistManager::Instance(); G4Material* nWater = nistMan->FindOrBuildMaterial("G4_WATER");
If I also use that in the TestEM0 example I end up with the large range for iron.
|Inline Depth:||Outline Depth:||Add message:|