Message: Re: Polyethylene / He3 counter with QGSP_BERT_HP physics list Not Logged In (login)
 Next-in-Thread Next-in-Thread
 Next-in-Forum Next-in-Forum

Ok Re: Polyethylene / He3 counter with QGSP_BERT_HP physics list 

Forum: Hadronic Processes
Re: Question Polyethylene / He3 counter with QGSP_BERT_HP physics list (Maduka Kaluarachchi)
Re: None Re: Polyethylene / He3 counter with QGSP_BERT_HP physics list (Koi, Tatsumi)
Date: 13 Sep, 2016
From: Maduka Kaluarachchi <Maduka Kaluarachchi>

Hi Tatsumi,

Thank you very much for the quick response and it worked. Now I have very good agreement with standard responses. I thought just using QGSP_BERT_HP will automatically activate thermal scattering model.

I'm briefly explaining the changes I did for this to work if any one has the same issue.

To activate the thermal neutron scattering model, I referred the post by Michael Kelsey (http://hypernews.slac.stanford.edu/HyperNews/geant4/get/hadronprocess/1471.html). I have attached the physics list I used based on this post.

Apart from this one more important change is the material definition. Materials relevant to thermal processes must be defined with their corresponding specific names defined in the class G4NeutronHPThermalScatteringNames. I have defined Polyethylene as follows,

G4Element* elTSH = new G4Element( "TS_H_of_Polyethylene", "H_POLYETHYLENE" , 1.0 , 1.0079*g/mole );

G4Material* HDpoly_TS = new G4Material( "Polyethylene_TS" , density = 0.963*g/cm3 , ncomponents = 2 );

HDpoly_TS -> AddElement(elTSH,natoms=2);

HDpoly_TS -> AddElement(elC,natoms=1);

Refer page 21 of this presentation by Tatsumi (http://geant4.slac.stanford.edu/SLACTutorial14/HadronicPhysics2.pdf).

Thank you again Tatsumi.

   Attachment:
      http://hypernews.slac.stanford.edu/HyperNews/geant4/get/AUX/2016/09/13/22.02-99241-PhysList.cc
      http://hypernews.slac.stanford.edu/HyperNews/geant4/get/AUX/2016/09/13/22.02-92818-PhysList.hh

 Add Message Add Message
to: "Re: Polyethylene / He3 counter with QGSP_BERT_HP physics list"

 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 ]