|Message: Heavy water simulation with NeutronHPThermalScattering||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'm simulating a simple geometry with geant4.10:
the setup is:
- a cube of air (NistAir) as experimental room volume - a cube of heavy water as target
I use GPS to launch neutrons of 5 MeV going through the cube of heavy water.
the questions are:
1) I define the heavy water as follows, using the names provided in G4NeutronHPThermalScatteringNames.hh:
G4Element* H2 = new G4Element("TS_D_of_Heavy_Water" ,"d_heavy_water" , 1., 2*1.0079*g/mole); G4Element* O = new G4Element("Oxygen","O" , 8., 16.00*g/mole);
//G4int ncomponents, natoms;
G4Material* H2O =new G4Material("Heavy_Water_ts", 1.100*g/cm3, ncomponents=2, kStateLiquid, 293*kelvin, 1*atmosphere); H2O->AddElement(H2, natoms=2); H2O->AddElement(O, natoms=1);
in order to obtain the thermal scattering for neutrons. Is this a correct definition of the material???
2) In my physicslist.cc I added:
// Neutron Physics RegisterPhysics( new NeutronHPphysics("neutronHP"));
where neutronHP refers to the file called NeutronHPphysics.cc of the example/extended/hadronic/Hdr04/src
is this the correct way to activate the thermal scattering process?
I'm asking because if I use the QGSP_BERT_HP only geant when running states in the terminal output:
Hadronic Process for neutron:
Model: hElasticCHIPS: 19.5 Mev ---> 100 Tev Model: NeutronHPElastic: 0 ev ---> 20 MevCr_sctns: NeutronHPElasticsXS: 0 ev ---> 20 Mev
Cr_sctns: ChipsNeutronElasticXS: 0 ev ---> 20 Mev Cr_sctns: GheishElastic: 0 ev ---> 20 Mev
[follow others processes]
while using my definition of Physicslist.cc written above I obtain:
Process: hadElastic Model: NeutronHPElastic: 4 eV ---> 20 MeV Model: NeutronHPThermalScattering: 0 meV ---> 4 eV Cr_sctns: NeutronHPThermalScatteringData: 0 meV ---> 4 eV Cr_sctns: NeutronHPElasticXS: 0 meV ---> 20 MeV Cr_sctns: GheishaElastic: 0 meV ---> 100 TeV
[follow others processes]
where Thermal scattering appears in the range 0 meV ---> 4 eV
Thank you so much!
|Inline Depth:||Outline Depth:||Add message:|