|Message: np-Elastic Cross-Section Discontinuity||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 am working on characterizing a plastic neutron scintillation detector's efficiency. To do this, I throw neutrons into the detector and see at what rate they interact and are detected. Thus, my result is mostly dependent on np-elastic cross-section and light production. I'm getting a very strange result, and it seems to be tied to the np-elastic.
I am using FTFP_BERT_HP and am attempting to look at the efficiency for neutrons ranging from a 5 MeV up to 200 MeV (with a detection threshold of 4 MeV). When I do this, I get a discontinuity at 32.7 MeV (see attached picture). If I look at the code here: $GEANTINSTALL/source/process/hadronic/cross_sections/G4HadronCrossSections.cc it seems that 32.7 MeV is the point where the cross-section tables stop and, I'm assuming, a new model is selected. This discontinuity is fatal to my calculation, so I need to find a way to overcome it.
I see that there is a SAID database result for np-elastic that comes pre-packaged in the GEANT Data Files I installed upon installing GEANT... however, I'm unsure how I tell my physics list to use only this data instead of whatever it is using now. Is there a way for me to tell GEANT that I want it to use only that data (or perhaps the ENDF result for nH)?
Currently, I just register my PhysicsList as FTFP_BERT_HP in the Main() and register extra processes as necessary (OpticalPhysics or what-have-you). So, I'm not even sure what cross-sections GEANT is truly using. I have tried using G4EmCalculator, but it always returns 0 pbarns as the cross-section for "hadElastic" process, so I can't tell what it is selecting. I'm hoping to make the EmCalculator and will update with those values if I can figure out why it's returning 0.
Thanks for any and all help, Zach
|Inline Depth:||Outline Depth:||Add message:|