|Message: bug? in Be(p,n) reaction gamma spectrums using QGSP_BERT_HP and QGSP_BIC_HP has big difference||Not Logged In (login)|
Click on the Forum title, e.g. on the "Forums by Category" page, to read a sequence of postings to the Forum and its threads all in one page. If you are only interested in one thread or the thread following a specific posting, click the thread or the posting, which takes you to a smaller page, which contains only the part you are interested in and may be easier to navigate.
Messages are "chained" if there are only replies at the first level, i.e. 1/1.html, 1/1/1.html etc. In case of "chained" messages the message number is replaced by the icon and there is no indentation.
Inline: Display the subject line only or also the text of the posting(s); for the choice "All" the "Outline" choices are switched off.
|1||0||1||no text / full text of posting|
|2||1||All||text for level 1 only / text for All postings|
Outline: Choose the depth of the posting thread, successive toggle controls provide increasing detail.
|1||2||1||2 levels / 1 level (original posting)|
|2||3||2||3 levels / 2 levels|
|3||3||All||3 levels / all levels (all postings)|
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");
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:||Outline Depth:||Add message:|