Problem 2096 - In the Bertini cascade there is a problem in interface to pre-compound
Summary: In the Bertini cascade there is a problem in interface to pre-compound
Status: RESOLVED FIXED
Alias: None
Product: Geant4
Classification: Unclassified
Component: processes/hadronic/models (show other problems)
Version: 10.4
Hardware: All All
: P4 minor
Assignee: dennis.herbert.wright
URL:
Depends on:
Blocks:
 
Reported: 2018-10-19 12:33 CEST by Vladimir.Ivantchenko
Modified: 2018-11-19 23:38 CET (History)
0 users

See Also:


Attachments
Proton spectra in n+Bi reaction at 63 MeV (test30) (15.42 KB, image/gif)
2018-10-19 12:33 CEST, Vladimir.Ivantchenko
Details

Note You need to log in before you can comment on or make changes to this problem.
Description Vladimir.Ivantchenko 2018-10-19 12:33:41 CEST
Created attachment 515 [details]
Proton spectra in n+Bi reaction at 63 MeV (test30)

If an interface to the Geant4 native pre-compound de-excitation is used by the Bertini cascade there is a problem of energy balance check. In the case of internal conversion an electron may be emitted instead of gamma. Electron mass is not taken into account in the balance computation and as a result re-sampling is forced several times. This increases CPU and damage de-excitation spectra. One of the plots demonstrating the problem is attached.
Comment 1 dennis.herbert.wright 2018-11-19 23:38:36 CET
Many problems found and fixed in interface to native G4PreCompound de-excitation. 
1) unit problems in going from MeV in G4Fragment to GeV in Bertini
2) electron charge and mass subtracted 
3) used a copy of G4Fragment from G4PreCompound instead of const ref because precompound can change the fragment
4) for momentum only, increased allowed non-conservation from 1 keV to 10 keV