|Message: Adding hadronic models||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)|
I am learning to simulate deuteron beam on beryllium target at about 20MeV. For hadronic process I use physics list QGSP_BIC_HP.
I noticed G4LEDeuteronInelastic is already included in G4IonBinaryCascadePhysics.cc in QGSP_BIC_HP. But if I register again the model G4LEDeuteronInelastic in my physics list, there is a big difference in neutron production. It gives ~50% more neutrons in one test simulation of 5M deuteron beams, compared with using QGSP_BIC_HP alone (for hadronic process).
I am a little confused why is this so. Will the order of registering hadronic processes in the code makes any difference?
My physics list is similar to hadronic example h01. This is what I add in end of PhysicsList::ConstructProcess():
// for low-energy deuteron physics G4ParticleDefinition *theDeuteron = G4Deuteron::DeuteronDefinition(); G4ProcessManager *deuteronProcMan = theDeuteron->GetProcessManager(); G4DeuteronInelasticProcess* deuteronInelasticProcess = new G4DeuteronInelasticProcess(); G4LEDeuteronInelastic* model = new G4LEDeuteronInelastic(); deuteronInelasticProcess->RegisterMe(model); deuteronProcMan->AddDiscreteProcess(deuteronInelasticProcess);
Thanks for your advice.
|Inline Depth:||Outline Depth:||Add message:|