Message: bug? in Be(p,n) reaction gamma spectrums using QGSP_BERT_HP and QGSP_BIC_HP has big difference Not Logged In (login)
 Next-in-Thread Next-in-Thread
 Next-in-Forum Next-in-Forum

Question bug? in Be(p,n) reaction gamma spectrums using QGSP_BERT_HP and QGSP_BIC_HP has big difference  

Forum: Hadronic Processes
Date: 08 Dec, 2013
From: ja <ja>

Hi!
   In my simulation of BNCT, energy of proton 9 MeV using the following code in main.cc

  G4VModularPhysicsList* physics = new QGSP_BIC_HP; 
or
  G4VModularPhysicsList* physics = new QGSP_BERT_HP; 

  physics->RegisterPhysics(new G4RadioactiveDecayPhysics());
  physics->RegisterPhysics(new G4EmStandardPhysics());
physics->SetCutValue(0.001*mm,"proton");
 physics->SetCutValue(0.01*mm,"gamma");
 physics->SetCutValue(0.01*mm,"electron");
 physics->SetCutValue(0.01*mm,"positron");

  runManager->SetUserInitialization(physics);

the gamma ray spectrum given by QGSP_BIC_HP has a strong peak at 3.5 MeV due to 9 Be(p,a)6 Li, while the spectrum obtained using QGSP_BERT_HP is a continues one without peak.

Both QGSP_BIC_HP and QGSP_BERT_HP are suggested in the G4 manual for simulation of this kind of issues.

So can any one tell me if there is a bug in my code or in G4 source code. And show me how to make the correct physics list in the main.cc

Any helpful response is deeply appreciated!

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

1 None: Re: bug? in Be(p,n) reaction gamma spectrums using QGSP_BERT_HP and QGSP_BIC_HP has big difference   (michel maire - 11 Dec, 2013)
1 None: Re: bug? in Be(p,n) reaction gamma spectrums using QGSP_BERT_HP and QGSP_BIC_HP has big difference   (=?ISO-8859-1?Q?Jos=E9_Manuel_Quesada_Molina?= - 02 Jan, 2014)
2 Idea: Re: bug? in Be(p,n) reaction gamma spectrums using QGSP_BERT_HP and QGSP_BIC_HP has big difference   (Vladimir Ivanchenko - 13 Dec, 2013)
1 Feedback: Re: bug? in Be(p,n) reaction gamma spectrums using QGSP_BERT_HP and QGSP_BIC_HP has big difference   (ja - 15 Dec, 2013)
 Add Message Add Message
to: "bug? in Be(p,n) reaction gamma spectrums using QGSP_BERT_HP and QGSP_BIC_HP has big difference "

 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 ]