|Message: Re: G4DNA Carbon Ions||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)|
after verifying the cross sections nothing comes out meaning all seems smooth. However, the issue is most probably caused by the switch between normal and relativistic formulas in the G4DNARuddIonisationExtendedModel Class, and specifically in functions RejectionFunction and ProposedSampledEnergy, the value at which we use relativistic formulas is around ~100 MeV/amu and for carbon ions it corresponds to ~ 1200 MeV.
So if you definitely need to eliminate this you'll have to decrease all the energy limits down to let's say zero (in the file G4DNARuddIonisationExtendedModel.cc) and therefore you'll have only relativistic formulas in action, however since we never tried this you will have to validate what you obtain before using it, to make sure the proper energy loss is still there.
Anyway we will be addressing this issue for our next release, but until then you may look for the magic numbers ((tau/MeV)<5.447761194e-2) and ((k/MeV)/(particle->GetPDGMass()/MeV) <= 0.1 ) and change this into ((tau/MeV)<0) and ((k/MeV)/(particle->GetPDGMass()/MeV) <= 0 ), recompile Geant4 and your application and check the rest out.
Hope this helps a bit Ziad
|Inline Depth:||Outline Depth:||Add message:|