Message: Re: Energy Conservation violation warning while using INCLXX physics list Not Logged In (login)
 Next-in-Thread Next-in-Thread
 Next-in-Forum Next-in-Forum

None Re: Energy Conservation violation warning while using INCLXX physics list 

Forum: Hadronic Processes
Re: Question Energy Conservation violation warning while using INCLXX physics list (Dipak Rimal)
Re: None Re: Energy Conservation violation warning while using INCLXX physics list (Davide Mancusi)
Date: 18 Dec, 2014
From: Dipak Rimal <Dipak Rimal>

Hi Davide,

Thank you for your response.

1) I figured out that the crash was unrelated to these warning but some other issues with G4AtomicShell. The error message using the latest G4 release is :

/geant4.10.01_src/source/materials/src/G4AtomicShells.cc:761: static G4double G4AtomicShells::GetBindingEnergy(G4int, G4int): Assertion `Z>0 && Z<101 && ShellNb<fNumberOfShells[Z]' failed. leading to run termination.

2) We are simulating 1 GeV proton + Hg interaction producing charged pions, which ultimately decay to produce neutrinos. We are studying neutrino flux and the low energy neutron background.

We are studying different physics lists that best describes these processes. We found that number of neutrinos produced per flavor per proton on target (also the hadron production) was entirely different while using the physics lists with Bertini cascade (BERT) and the Binary Intranuclear cascade (BIC) models. So we are trying INCLXX as an alternative.

> On Mon, 15 Dec 2014 22:29:12 GMT, Davide Mancusi wrote:
> Hello Dipak,
> 
> On Mon, 15 Dec 2014 19:27:49 GMT, Dipak Rimal wrote:
> 
> > Hello G4 experts I am trying to use INCLXX physics to simulate
> > interaction of 1 GeV proton on Hg target. I am getting frequent energy
> > conservation violation warnings such as:
> > 
> > [INCL++] Warning: energy conservation violated by 6.73985 MeV in
> > 47.8781-MeV neutron + Be9 inelastic reaction, in direct kinematics. Will
> > resample. [INCL++] Warning: maximum number of tries exceeded for the
> > proposed 5.14186-MeV proton + Hg201 inelastic reaction, in direct
> > kinematics.
> 
> These warnings are harmless. The first one is typically generated in
> corner cases of reactions on light nuclei. Note that events violating
> energy conservation are resampled, so your calculation will always
> satisfy energy conservation. The second warning is emitted because
> Geant4 requested an inelastic event for the mentioned projectile-target
> combination (because the parametrized proton reaction cross section does
> not vanish), but INCL++ cannot generate one (because e.g. the proton is
> below INCL++'s internal Coulomb barrier).
> 
> Both cases are corner cases. Similar things happen in all models; INCL++
> is just a more chatty about it. You should not worry too much about
> them, unless you are specifically interested in e.g. sub-barrier
> reactions.
> 
> > The simulation even quits sometimes after exceeding maximum number of
> > tries. I tried the latest release of GEANT4 ( 10.1 ) too but I still get
> > same warning messages.
> 
> If the simulation quits or crashes unexpectedly, it cannot be due to
> these warnings.
> 
> > This happens for low energies interactions. Since I am working on low
> > energy regime, I am a bit concerned about these warnings. Can you please
> > comment on the severity of these warnings and possible resolutions?
> 
> If you are interested in a precise description of reactions below ~50
> MeV, then intranuclear-cascade models are not your best bet. How low is
> your low energy?
> 
> Cheers, Davide
> 

On Mon, 15 Dec 2014 22:29:12 GMT, Davide Mancusi wrote:

> Hello Dipak,
> 
> On Mon, 15 Dec 2014 19:27:49 GMT, Dipak Rimal wrote:
> 
> > Hello G4 experts I am trying to use INCLXX physics to simulate
> > interaction of 1 GeV proton on Hg target. I am getting frequent energy
> > conservation violation warnings such as:
> > 
> > [INCL++] Warning: energy conservation violated by 6.73985 MeV in
> > 47.8781-MeV neutron + Be9 inelastic reaction, in direct kinematics. Will
> > resample. [INCL++] Warning: maximum number of tries exceeded for the
> > proposed 5.14186-MeV proton + Hg201 inelastic reaction, in direct
> > kinematics.
> 
> These warnings are harmless. The first one is typically generated in
> corner cases of reactions on light nuclei. Note that events violating
> energy conservation are resampled, so your calculation will always
> satisfy energy conservation. The second warning is emitted because
> Geant4 requested an inelastic event for the mentioned projectile-target
> combination (because the parametrized proton reaction cross section does
> not vanish), but INCL++ cannot generate one (because e.g. the proton is
> below INCL++'s internal Coulomb barrier).
> 
> Both cases are corner cases. Similar things happen in all models; INCL++
> is just a more chatty about it. You should not worry too much about
> them, unless you are specifically interested in e.g. sub-barrier
> reactions.
> 
> > The simulation even quits sometimes after exceeding maximum number of
> > tries. I tried the latest release of GEANT4 ( 10.1 ) too but I still get
> > same warning messages.
> 
> If the simulation quits or crashes unexpectedly, it cannot be due to
> these warnings.
> 
> > This happens for low energies interactions. Since I am working on low
> > energy regime, I am a bit concerned about these warnings. Can you please
> > comment on the severity of these warnings and possible resolutions?
> 
> If you are interested in a precise description of reactions below ~50
> MeV, then intranuclear-cascade models are not your best bet. How low is
> your low energy?
> 
> Cheers, Davide
> 

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

1 Idea: Re: Energy Conservation violation warning while using INCLXX physics list   (Vladimir Ivanchenko - 18 Dec, 2014)
 Add Message Add Message
to: "Re: Energy Conservation violation warning while using INCLXX 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 ]