Created attachment 376 [details] patched source files There are a lot of major issues in G4ParticleHPInelasticCompFS, few examples: 1. gammas are generated from all (n,n') excited states, also those which decay alpha o p,n,d,t,h etc. 2. all continuum inelastic channels MT=91,649,699,749,799,849 made decay elastically with gamma emission. Instead these are perhaps break-up reactions involving lighter nuclei, not gammas. 3. for A(n,p,d,t,h,a)B* reactions levels refer to a different nucleus (B* with Z'=Z,Z-1,Z-2, A'=A,A-1,A-2,A-3), so we need not only the table of levels for the current nucleus A, but also for other 5 nuclei possible (not corrected!). 4. G4ParticleHPAngular recalculates hadron momentum, which is fixed in G4ParticleHPInelasticCompFS, I added a flag to tell it to keep the existing momentum, although not yet checked with other routines calling it. Perhaps G4ParticleHPAngular need further improvements. 5. G4ParticleHPEnAngCorrelation does not care about energy/momentum conservation, its a pain to see all these gammas. Perhaps reworking it in CM frame could fix that (not corrected!). 6. Be8 has not decay table! There is a lot of debugging prints in the proposed code, but it is working. I can try to improve this but I need a feedback (I am focused only on my experiment).
Thank you for reporting this problem. Now we are working on this issue. I will let you know once we have progress. Tatsumi
This will become a development item.
