| Summary: | Boron10 in G4NDL4.5 | ||
|---|---|---|---|
| Product: | Geant4 | Reporter: | michel.maire |
| Component: | processes/hadronic/models/neutron_hp | Assignee: | Emilio <emilio.mendoza> |
| Status: | RESOLVED FIXED | ||
| Severity: | minor | CC: | Alberto.Ribon, dennis.herbert.wright, michel.maire |
| Priority: | P4 | ||
| Version: | 10.4 | ||
| Hardware: | All | ||
| OS: | All | ||
| Attachments: |
output of Hadr03
plot of Q |
||
Created attachment 554 [details]
plot of Q
Dear Michel, I have been investigating this problem, and here are the conclusions. The neutron + B10 --> 2 alpha + triton reaction do not have final state data in G4NDL4.5. Thus, ParticleHP uses a model to generate the two alphas and the triton. The model do not try to conservate energy and momentum event-by-event, but on average. On top of this, ParticleHP tries to conservate energy and momentum event-by-event in some situations, by default. This option can be disabled by defining the environmental variable G4NEUTRONHP_DO_NOT_ADJUST_FINAL_STATE (highly recommended, in my opinion). Now, when running Geant4 with the two options you get: default --> Q = 1.85 MeV G4NEUTRONHP_DO_NOT_ADJUST_FINAL_STATE --> Q = 0 And both are wrong. I have been investigating the problem and I found a bug. When fixed the results are: default --> Q = 2.09 MeV G4NEUTRONHP_DO_NOT_ADJUST_FINAL_STATE --> Q = 327 keV So now the Q-value is OK when running when the G4NEUTRONHP_DO_NOT_ADJUST_FINAL_STATE option. In the other case Geant4 tries to force energy conservation event-by-event, unsuccessfully. For the moment, the issue is fixed when using the G4NEUTRONHP_DO_NOT_ADJUST_FINAL_STATE option. Fixing the other option can be done for this specific case but will require a much greater effort. This can be done in a future. However, it is not possible, in all the cases, to force energy conservation when using G4NDL data libraries. Only in some specific cases. Regards, Emilio HI Emilio, Thanks for your investigation. Some time ago (more than one year ...) there were discussions within hadronic group about the necessity to simplify, reduce the number of environmental variables in ParticleHP, and do not adjust final state by default. I have never seen any progress in this direction ... |
Created attachment 553 [details] output of Hadr03 See hypernews hadronic #1687 from Dongwon Lee. Neutron (5 MeV) on Boron10. In attachement, output of Hadr03 for 2 physics constructors: BinaryCascade and neutronHP (with G4NDL4.5). The second attachement is plot of Q: BIC in top, neutronHP in bottom. From which we can conclude that channel 1 of neutronHP (line 88 of output) is wrong.