Problem 1790

Summary: some issues with continuum in G4NeutronHPInelasticCompFS
Product: Geant4 Reporter: Artem Zontikov <zontikov.a>
Component: processes/hadronic/models/neutron_hpAssignee: tkoi
Status: RESOLVED FIXED    
Severity: trivial CC: dennis.herbert.wright
Priority: P5    
Version: 10.1   
Hardware: All   
OS: All   

Description Artem Zontikov 2015-10-02 00:21:31 CEST
1. Here http://www-geant4.kek.jp/lxr/source/processes/hadronic/models/neutron_hp/src/G4NeutronHPInelasticCompFS.cc#L327 maximum energy for scattered neutron is defined as E'=E+QI (where QI < 0) but ENDF-6 Formats Manual states for MT=91 that 0 <= E' <= E-U (where U is the threshold energy). See chapters 5.2 and 5.3 for details.

2. Here http://www-geant4.kek.jp/lxr/source/processes/hadronic/models/neutron_hp/src/G4NeutronHPInelasticCompFS.cc#L340 level increment could lead to further deexcitation via photon with energy greater than incident neutron energy. Is it a feature?

3. ENDF standard say that photons from the continuum should be present in MF=12 - MF=15. In case there is no such data it may be an idea worth considering to add at http://www-geant4.kek.jp/lxr/source/processes/hadronic/models/neutron_hp/src/G4NeutronHPInelasticCompFS.cc#L416 something like:
if ( !(G4NeutronHPManager::GetInstance()->GetDoNotAdjustFinalState()) ) thePhotons->push_back(theRestEnergy);
This will implement the treatment of continuum in the same manner as adopted for discrete levels.

Or even better add a specific function in G4NeutronHPManager which allows user to skip gammas from deexcitation of continuum - all of them or only gammas which are not from discrete levels.
Comment 1 tkoi 2015-11-14 04:53:58 CET

*** This problem has been marked as a duplicate of problem 1789 ***