Message: bug in HP for gammas? Not Logged In (login)
 Next-in-Thread Next-in-Thread
 Next-in-Forum Next-in-Forum

Question bug in HP for gammas? 

Forum: Hadronic Processes
Date: 08 Jun, 2012
From: Zoltan Elekes <Zoltan Elekes>

Dear Developers,

When I was comparing experimental gamma-ray spectra induced by low-energy neutrons (~2 MeV) on different materials (including e.g., copper) taken with a HPGe detector, I have observed that simulations performed with version 9.2 matches quite well with the experiment using QGSP_BIC_HP physics list. However, version 9.4 gave very different result and no match with experimental spectrum. I discussed this problem with experts Vladimir Ivantchenko and Laurent Desorgher and they suggested that the problem might be with HP package.

Therefore, I have performed a simplified simulation. I shoot 1.4 MeV neutrons on a copper cube and detect the gamma-rays with a Ge crystal. I simulated the spectra with version 9.2, 9.4 and 9.5. The results are attached, the x axis is in keV. As is seen, version 9.2 can nicely reproduce the expected gammas and their relative strengths, however version 9.4 and 9.5 give completely wrong gamma-ray energies. Somehow, they are shifted and the relative strengths are also not OK.

Minor problem: using an iron cube, I have also found that the the gamma-ray energy of 1238 keV in 56Fe (from level 2085 keV to level 847 keV) appears at 1233 keV even with version 9.2.

The datafiles, I have used are the following:

version 9.2: PhotonEvaporation2.0 G4EMLOW6.2 G4NDL3.13 RadioactiveDecay3.2 G4ABLA3.0

version 9.4: PhotonEvaporation2.1 G4EMLOW6.19 G4NDL3.14 RadioactiveDecay3.3 G4ABLA3.0 G4NEUTRONXS1.0 G4PII1.2 RealSurface1.0

version 9.5: PhotonEvaporation2.2 G4EMLOW6.23 G4NDL4.0 RadioactiveDecay3.4 G4ABLA3.0 G4NEUTRONXS1.1 G4PII1.3 RealSurface1.0

Regards, Zoltan

   Attachment:
      http://hypernews.slac.stanford.edu/HyperNews/geant4/get/AUX/2012/06/08/02.33-34613-picture9.2.png
      http://hypernews.slac.stanford.edu/HyperNews/geant4/get/AUX/2012/06/08/02.33-33586-picture9.4.png
      http://hypernews.slac.stanford.edu/HyperNews/geant4/get/AUX/2012/06/08/02.33-84445-picture9.5.png

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

1 Feedback: Re: bug in HP for gammas?   (Malachi Schram - 08 Jun, 2012)
1 More: Re: bug in HP for gammas?   (Malachi Schram - 11 Jun, 2012)
1 Note: Re: bug in HP for gammas?   (Zoltan Elekes - 11 Jun, 2012)
... 1 Message(s)
2 None: Re: bug in HP for gammas?   (Koi, Tatsumi - 16 Jul, 2012)
1 Feedback: Re: bug in HP for gammas?   (Malachi Schram - 17 Jul, 2012)
2 None: Re: bug in HP for gammas?   (Zoltan Elekes - 18 Jul, 2012)
... 3 Message(s)
 Add Message Add Message
to: "bug in HP for gammas?"

 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 ]