|Message: Choosing a different fission model in QGSP_BIC_HP physics list||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 trying to use only the LLNL fission model and the Binary Cascade. The LLNL model has an upper energy limit of 20 MeV and that is great for me. I do not need any fissions for higher energies.
So, I removed any fission reference from G4NeutronHPBuilder, G4NeutronBuilder and G4VNeutronBuilder, and I implemented the LLNL fission model in MYPhysicsList which I registered like this:
QGSP_BIC_HP* physics = new QGSP_BIC_HP();
physics->RegisterPhysics(new MYPhysicsList()); // where the fission processes are defined, with the new model
When I start my code, I get the processes lists which I expected;
Hadronic Processes for <neutron> hadElastic Models: hElasticCHIPS: Emin(GeV)= 0.0195 Emax(GeV)= 100000 NeutronHPElastic: Emin(GeV)= 0 Emax(GeV)= 0.02 NeutronInelastic Models: QGSP: Emin(GeV)= 12 Emax(GeV)= 100000 G4LENeutronInelastic: Emin(GeV)= 9.5 Emax(GeV)= 25 Binary Cascade: Emin(GeV)= 0.0199 Emax(GeV)= 9.9 NeutronHPInelastic: Emin(GeV)= 0 Emax(GeV)= 0.02 nCapture Models: G4LCapture: Emin(GeV)= 0.0199 Emax(GeV)= 20000 NeutronHPCapture: Emin(GeV)= 0 Emax(GeV)= 0.02 nFission Models: HadronicModel: Emin(GeV)= 0 Emax(GeV)= 0.02
However, during the run , the code crushes, with the error:
G4EnergyRangeManager:GetHadronicInteraction: counter=1, Ek=34.311593, Material = Thorium, Element = Thorium *0* low=0, high=20 In src/G4EnergyRangeManager.cc, line 123: ===> GetHadronicInteraction: No Model found Target element Thorium Z= 90 A= 232 Unrecoverable error in the method ChooseHadronicInteraction of nFission TrackID= 20858 ParentID= 19440 neutron Ekin(GeV)= 0.034311593; direction= (-0.59228123,0.80261939,0.070746516) Position(mm)= (117.42113,-57.470852,46.143174); material Thorium PhysicalVolume <target>
*** G4Exception : 007 issued by : G4HadronicProcess ChooseHadronicInteraction failed. *** Fatal Exception *** core dump ***
*** G4Exception: Aborting execution *** Aborted
So it is still trying to simulate fission , even if the neutron energy is outside the 0 - 20 MeV range.
I would be very grateful if you could please suggest a way of switching off fission above 20 MeV neutron energy.