Message: Re: Non energy conservation in a hadron elastic scattering between neutrons and protons Not Logged In (login)
 Next-in-Thread Next-in-Thread
 Next-in-Forum Next-in-Forum

Warning Re: Non energy conservation in a hadron elastic scattering between neutrons and protons 

Forum: Hadronic Processes
Re: Warning Non energy conservation in a hadron elastic scattering between neutrons and protons (Romain Roncin)
Re: Idea Re: Non energy conservation in a hadron elastic scattering between neutrons and protons (Vladimir Ivanchenko)
Date: 14 May, 2012
From: Romain Roncin <Romain Roncin>

Dear all,

To show you that there is something strange I used the QGSP_BERT_HP reference physics list instead of what I used before. I saw exactly the same results.

I join you two plots made with the QGSP_BERT_HP reference physics list : PR_10MeV.pdf corresponds to the number of protons recoils vs their vertex kinetic energy for 10MeV neutrons created in a liquid scintillator and PR_100MeV.pdf is for 100MeV neutrons created in the same liquid scintillator.

You can see that if everything is working well for the 10MeV neutrons (10000 neutrons created), there is a huge peak at 100MeV for the 100MeV neutrons (10000 neutrons created). It corresponds to what I put in my first mail : the protons are taking in only one step all the neutrons energies and even more, which is physically impossible :

- neutron  2112  E: 0 keV;  E0: 100000 keV;  Vertex Position: (0,3230,0)
  TrackID: 1;  ParentID: 0;  Step: 1;  HadronElastic
- proton  2212  E: 100054 keV;  E0: 100125 keV;  Vertex Position: (0,3058.76,0)
  TrackID: 2;  ParentID: 1;  Step: 1;  hIoni
- proton  2212  E: 99797 keV;  E0: 100125 keV;  Vertex Position: (0,3058.76,0)
  TrackID: 2;  ParentID: 1;  Step: 2;  hIoni
- proton  2212  E: 99612.8 keV;  E0: 100125 keV;  Vertex Position: (0,3058.76,0)
  TrackID: 2;  ParentID: 1;  Step: 3;  hIoni

Do you have any idea about what is happening ?

Best regards,

Romain Roncin

   Attachment:
      http://hypernews.slac.stanford.edu/HyperNews/geant4/get/AUX/2012/05/14/00.03-32114-PR_10MeV.pdf
      http://hypernews.slac.stanford.edu/HyperNews/geant4/get/AUX/2012/05/14/00.03-7671-PR_100MeV.pdf

 Add Message Add Message
to: "Re: Non energy conservation in a hadron elastic scattering between neutrons and protons"

 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 ]