Message: Re: Problem with overlapping models for nCapture Not Logged In (login)
 Next-in-Thread Next-in-Thread
 Next-in-Forum Next-in-Forum

None Re: Problem with overlapping models for nCapture 

Forum: Hadronic Processes
Re: Question Problem with overlapping models for nCapture (Arianna)
Date: 29 Jun, 2018
From: michel maire <michel maire>

On Tue, 26 Jun 2018 15:09:13 GMT, Arianna wrote:
> Hello, I'm having troubles with some simulations of neutron generation
> driven by moderate energy protons in a pitcher-catcher configuration.
> 
> I started from the example in
> geant4.10.04.p02/examples/extended/hadronic/NeutronSource and modified
> it to my needs.
> 
> From previous works I know that the data-driven model implemented in
> G4HadronPhysicsQGSP_BIC_AllHP works better than
> G4HadronPhysicsQGSP_BIC_HP. So, basically, I replaced
> G4HadronPhysicsQGSP_BIC_HP with G4HadronPhysicsQGSP_BIC_AllHP in the
> PhysicsList.cc of the example.
> 
> The problem is that using a LiF proton-neutron converter, I often get
> the following error (it didn't happen with Be or Li):
> 
> -------- EEEE ------- G4Exception-START -------- EEEE -------
> *** G4Exception : had005
>       issued by : G4HadronicProcess::PostStepDoIt
> In /geant4.10.04.p01/source/processes/hadronic/management/src/G4EnergyRangeManager.cc, line 149: 
> ===> GetHadronicInteraction: Energy ranges of two models fully overlapping
> Target element F  Z= 9  A= 19
> Unrecoverable error in the method ChooseHadronicInteraction of nCapture
> TrackID= 2  ParentID= 1  neutron
> Ekin(GeV)= 0.000875406;  direction= (0.238566,-0.631628,-0.737653)
> Position(mm)= (1.3317,-5.29241,-2.13383);  material G4_LITHIUM_FLUORIDE
> PhysicalVolume  <Target>
>  No HadronicInteraction found out
> 
> *** Fatal Exception *** core dump *** -------- EEEE --------
> G4Exception-END --------- EEEE -------
> 
> My question is: how do I know which models are overlapping and how do I
> change the min/max energies in order to remove the overlapping, but
> still have a reliable simulation?
> 

 With QGSP_BIC_HP, I read from output :

                           Hadronic Processes for neutron

  Process: nCapture
        Model:          NeutronHPCapture: 0 meV ---> 20 MeV
        Model:               nRadCapture: 19.9 MeV ---> 100 TeV
     Cr_sctns:        NeutronHPCaptureXS: 0 meV ---> 20 MeV
     Cr_sctns:        G4NeutronCaptureXS: 0 meV ---> 100 TeV

 While with AllHP it becomes :

  Process: nCapture
        Model:          NeutronHPCapture: 0 meV ---> 20 MeV
        Model:               nRadCapture: 0 meV ---> 100 TeV
     Cr_sctns:        G4NeutronCaptureXS: 0 meV ---> 100 TeV

 I believe this is a mistake. I suggest you to open a bug report.

           Michel

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

1 Idea: Re: Problem with overlapping models for nCapture   (Vladimir Ivanchenko - 03 Jul, 2018)
(_ None: Re: Problem with overlapping models for nCapture   (michel maire - 08 Jul, 2018)
 Add Message Add Message
to: "Re: Problem with overlapping models for nCapture"

 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 ]