Message: Choosing a different fission model in QGSP_BIC_HP physics list Not Logged In (login)
 Next-in-Thread Next-in-Thread
 Next-in-Forum Next-in-Forum

Question Choosing a different fission model in QGSP_BIC_HP physics list 

Forum: Hadronic Processes
Date: 25 Oct, 2010
From: Cristian Bungau <Cristian Bungau>

Dear all,

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

  runManager->SetUserInitialization(physics); 

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.

Thank you,

Best regards,

Cristian

 Add Message Add Message
to: "Choosing a different fission model in QGSP_BIC_HP physics list"

 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 ]