Problem 1903 - alpha (5 MeV) + Be9 reaction
Summary: alpha (5 MeV) + Be9 reaction
Status: RESOLVED FIXED
Alias: None
Product: Geant4
Classification: Unclassified
Component: processes/hadronic/models (show other problems)
Version: 10.2
Hardware: All All
: P4 minor
Assignee: Vladimir.Ivantchenko
URL:
Depends on:
Blocks:
 
Reported: 2016-10-25 16:00 CEST by michel.maire
Modified: 2018-10-01 17:40 CEST (History)
0 users

See Also:


Attachments
output of Hadr03 (2.06 KB, text/plain)
2017-09-29 11:52 CEST, michel.maire
Details
This is the result of Hadr03 for 10.4ref08 (14.35 KB, text/plain)
2018-09-26 19:18 CEST, Vladimir.Ivantchenko
Details

Note You need to log in before you can comment on or make changes to this problem.
Description michel.maire 2016-10-25 16:00:33 CEST
I am still on the reaction alpha(5.4 MeV) on Be9 !
I noticed differences between ref06 and ref07.
In ref06, there are important quasi-elastic interactions and few 3 alphas production.
In ref07, these 2 channels have disappeared. 
As a consequence, the neutron production increases from 50% to 100% … which is not a small effect !
Is this expected ?
Comment 1 dennis.herbert.wright 2016-11-29 18:13:38 CET
Which physics list are you using?
Comment 2 michel.maire 2016-11-29 19:02:30 CET
In Hadr03, constructor G4IonPhysics

# Macro file for "Hadr03.cc"
#
/control/verbose 2
/run/verbose 1
#
/testhadr/det/setIsotopeMat Be9 4 9 1.85 g/cm3
#
/run/initialize
#
/process/list
/process/inactivate hadElastic
#
/gun/particle alpha
/gun/energy 5 MeV
#
/run/printProgress 10000
/run/beamOn 100000
Comment 3 dennis.herbert.wright 2016-12-07 00:59:23 CET
The hadronic excitation code as used by G4Precompound is likely responsible and there have been changes there recently. I am re-assigning this to Vladimir.
Comment 4 Vladimir.Ivantchenko 2017-09-28 19:48:57 CEST
Hi Michel,

can you, please,  run it with 10.3ref09 ?

Vladimir
Comment 5 michel.maire 2017-09-29 11:52:51 CEST
Created attachment 471 [details]
output of Hadr03

HI Vladimir,

  output in attachement
Comment 6 Vladimir.Ivantchenko 2018-09-26 19:12:40 CEST
Hi Michel,

I have tried to check if we have a problem or not with this reaction now. In the attachment there is the result. I obtain it with QBBC hadron inelastic and G4IonPhysicsXS ion inelastic. There are a significant fraction of 3 alpha final state but for some reason in the printout there is no gamma, however, tracking verbose shows  gamma in each event.

Can we close this problem report or not?

Cheers,
Vladimir
Comment 7 Vladimir.Ivantchenko 2018-09-26 19:18:37 CEST
Created attachment 512 [details]
This is the result of Hadr03 for 10.4ref08
Comment 8 michel.maire 2018-09-27 15:27:56 CEST
HI Vladimir,

 The bug report was created because I noticed signifiant differences between 10.2 ref-06 and ref-07. Since this time, the situation is rather stable, including 10.4 ref-08. But I do not know if the present neutron yield is reasonable or not.

 Concerning the number of gamma per event, I think it is just a problem of merging in MT mode. I will have a look.

     Michel
Comment 9 Vladimir.Ivantchenko 2018-10-01 17:39:59 CEST
Because Hadr03 shows stable results, we can close this problem.

VI