|Message: Re: HP Physics List||Not Logged In (login)|
|1||0||1||no text / full text of posting|
|2||1||All||text for level 1 only / text for All postings|
|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)|
Ok, I'm willing to upgrade to version 8, but first I'd like to explain the particular problem I'm having with this physics list and hopefully you can tell me whether or not this problem is fixed in version 8.
I am shooting thermal neutrons through Lithium 6. I have set up my material using G4Isotope objects so that I know my material is in fact pure lithium 6 and not natural lithium. When the step is decided which process to do it sees my material as lithium 6 and I can tell this because it chooses the NeutronInelastic process and creates an alpha and a triton, which is good since the cross section for that process at that energy is around 10^5 barns. However, in looking through the code where it actually performs the process, I find that it calls a function ChooseAandZ(const G4DynamicParticle *aParticle, const G4Material *aMaterial ) to determine what element it is shooting into. In this method is uses data tables to find the natural abundancy of each isotope of lithium and then uses those as probabilities to determine which isotope it is shooting into, instead of using the type of lithium that I set up for the material. Of course it chooses lithium 7 most of the time since that is much more abundant than lithium 6. This is giving me incorrect initial energies for my alpha and triton. Honestly, I'm not sure how it is giving me an alpha and triton at all if it is using lithium 7 since there should be an extra particle left over.
If you could tell me whether or not this problem has been fixed in version 8, before I switch over to it, that would be incredibly helpful, Thank you.
|Inline Depth:||Outline Depth:||Add message:|