|Message: Missing Implementation of Deuteron Elastic Cross Sections||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 want to study deuteron and antideuteron interaction with material at low energies (below 10 GeV), and by printing out the elastic and inelastic cross-sections implemented in Geant4, I noticed that the deuteron elastic cross section is actually not implemented at all: hardcoded to 0 in line 1320 of G4HadronCrossSections.cc. This was very surprising, since for antideuterons it is implemented (in G4ComponentAntiNuclNuclearXS.cc). Searching in this forum I found an old thread (http://hypernews.slac.stanford.edu/HyperNews/geant4/get/hadronprocess/608.html) which already pointed out this missing implementation. The answer was that it should be fixed in the G4UHadronElasticProcess class since version 8.1, but this class, as well as G4WHadronElasticProcess which followed in later versions, do not exist anymore. So my questions are:
1) Is there a reason for the missing implementation of the deuteron elastic cross-sections? Especially, why is it implemented for antideuterons but not for deuterons?
2)Can you suggest me any model that I could implement myself in order to fix this problem?
Note that this problem is valid not only for deuterons, but also for tritons and alpha particles.
Thank you and Best Regards,
Technical University Munich
|Inline Depth:||Outline Depth:||Add message:|