Message: Heavy water simulation with NeutronHPThermalScattering Not Logged In (login)
 Next-in-Thread Next-in-Thread
 Next-in-Forum Next-in-Forum

Question Heavy water simulation with NeutronHPThermalScattering 

Forum: Hadronic Processes
Date: 18 Mar, 2014
From: filippo <filippo>

Hello!

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:

Process: hadElastic

Model: hElasticCHIPS: 19.5 Mev ---> 100 Tev
Model: NeutronHPElastic: 0 ev ---> 20 Mev
Cr_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!

Filippo

Inline Depth:
 1 1
 All All
Outline Depth:
 1 1
 2 2
 All All
Add message: (add)

1 Idea: Re: Heavy water simulation with NeutronHPThermalScattering   (Vladimir Ivanchenko - 18 Mar, 2014)
(_ None: Re: Heavy water simulation with NeutronHPThermalScattering   (filippo - 18 Mar, 2014)
 Add Message Add Message
to: "Heavy water simulation with NeutronHPThermalScattering"

 Subscribe Subscribe

This site runs SLAC HyperNews version 1.11-slac-98, derived from the original HyperNews


[ Geant 4 Home | Geant 4 HyperNews | Search | Request New Forum | Feedback ]