|Message: Re: Thermal neutron scattering on plexiglass||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 have been messing around with G4NeutronHPThermalScattering, and I can't tell you what is the best model for your situation, but I can at least tell you how to turn on G4NeutronHPThermalScattering so that you can see the difference yourself.
> As far as I > know, this is not enabled by default in Geant4.9.3.p2, which is the > version I am using.
No physics processes are enabled by default in any version of Geant4, it always depends on your physics list. However, none of the reference physics lists that I have found include thermal scattering. You will need to create your own physics list and register G4NeutronHPThermalScattering. The easiest way to do this is to create a custom G4VPhysicsConstructor and use it in your own modular physics list. I simply modified G4HadronElasticPhysicsHP to include thermal scattering below 4 eV (note: elastic hadronic interactions are always in a separate module from the rest of the hadronic interactions). Attached are the modified files (ThermalNeutronScattering.hh/.cc). You should be able to just add these to your source and then add the following line to your modular physics list: this->RegisterPhysics( new ThermalNeutronScattering(verboseLevel) );
> Would "G4NeutronHPThermalScattering" actually > apply to G4_PLEXIGALSS?
No, the list of materials applicable to G4NeutronHPThermalScattering can be found in G4NeutronHPThermalScatteringNames.cc. The NIST materials are listed at bottom, i.e. G4_something. The group at the top are elements, i.e. if you define an element with this exact name and then use it in a material, G4NeutronHPThermalScattering will be applied. A run time message will list all materials/elements that are recognized by G4NeutronHPThermalScattering.
Now, even though "PLEXIGLASS" is not one the relevant materials, you might try using polyethylene as an approximation. If you use G4_PLEXIGLASS, then G4NeutronHPElastic will be used, which contains no S(α,β) information and is based on a free gas approximation. What model is best is a question for a neutronics expert. I did notice that T. Koi uses polyethlyene as an approximation for PMMA in a presentation that I found on the web.
Attachment: http://hypernews.slac.stanford.edu/HyperNews/geant4/get/AUX/2012/09/06/10.16-21769-lNeutronScattering.cc.txt http://hypernews.slac.stanford.edu/HyperNews/geant4/get/AUX/2012/09/06/10.16-41759-lNeutronScattering.hh.txt