|Message: fission above 20 MeV||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 simulate neutron production in an uranium target. I use standard QGSP_INCL_ABLA physics list and I am interested in energy range below 1 GeV. INCL_ABLA simulates neutron induced fission reactions (I can see fission fragments) but it calls them "NeutronInelastic". At the same time "Low Energy model" also simulates fission reactions and calls them "nFission" (no fission fragments and no energy deposition). Therefore fission reaction are simulated twice: by INCL/ABLA and by LowEnergy model... And LowEnergy model does it not accurately.
QGSP_INCL_ABLA physics list:
... theNeutrons=new G4NeutronBuilder; theNeutrons->RegisterMe(theQGSPNeutron=new G4QGSPNeutronBuilder(QuasiElastic, ProjectileDiffraction)); theNeutrons->RegisterMe(theLEPNeutron=new G4LEPNeutronBuilder); theLEPNeutron->SetMinInelasticEnergy(9.5*GeV); theLEPNeutron->SetMaxInelasticEnergy(25*GeV);
theNeutrons->RegisterMe(theBertiniNeutron=new G4BertiniNeutronBuilder); theBertiniNeutron->SetMinEnergy(2.9*GeV); theBertiniNeutron->SetMaxEnergy(9.9*GeV);
theNeutrons->RegisterMe(theInclAblaNeutron=new G4InclAblaNeutronBuilder); theInclAblaNeutron->SetMinEnergy(0.0*GeV); theInclAblaNeutron->SetMaxEnergy(3.0*GeV); ...
Is it possible to disable Low Energy Fission for neutron energy < 3 GeV (where INCL_ABLA works)? If I add line
I get an exception:
G4EnergyRangeManager:GetHadronicInteraction: counter=2, Ek=25.8938, Material = TargetMat, Element = el_U *0* low=3000, high=2e+07 *1* low=0, high=20 In src/G4EnergyRangeManager.cc, line 127: ===> GetHadronicInteraction: No Model found Target element el Z= 92 A= 238 Unrecoverable error in the method ChooseHadronicInteraction of nFission TrackID= 1655 ParentID= 781 neutron Ekin(GeV)= 0.0258938; direction= (-0.640663,-0.730955,0.235066) Position(mm)= (-4.60825,0.0549675,-303.604); material TargetMat PhysicalVolume <Target>
Thanks for any help in advance.
|Inline Depth:||Outline Depth:||Add message:|