Problem 1903

Summary:	alpha (5 MeV) + Be9 reaction
Product:	Geant4	Reporter:	michel.maire
Component:	processes/hadronic/models	Assignee:	Vladimir.Ivantchenko
Status:	RESOLVED FIXED
Severity:	minor
Priority:	P4
Version:	10.2
Hardware:	All
OS:	All
Attachments:	output of Hadr03 This is the result of Hadr03 for 10.4ref08

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