Message: fission above 20 MeV Not Logged In (login)
 Next-in-Thread Next-in-Thread
 Next-in-Forum Next-in-Forum

None fission above 20 MeV 

Forum: Hadronic Processes
Date: 02 Dec, 2011
From: Yury Malyshkin <Yury Malyshkin>

Hi all,

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

  theLEPNeutron->SetMinEnergy(3*GeV);

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.

Yury

Inline Depth:
 1 1
 All All
Outline Depth:
 1 1
 2 2
 All All
Add message: (add)

1 Idea: Re: fission above 20 MeV   (Vladimir Ivanchenko - 09 Jan, 2012)
(_ None: Re: fission above 20 MeV   (Yury Malyshkin - 25 Jan, 2012)
 Add Message Add Message
to: "fission above 20 MeV"

 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 ]