Problem 1000

Summary: Energy balance for binary cascade
Product: Geant4 Reporter: gennady.vasilyev
Component: processes/hadronic/models/binary_cascadeAssignee: Gunter.Folger
Status: RESOLVED FIXED    
Severity: trivial    
Priority: P5    
Version: 9.1   
Hardware: All   
OS: All   

Description gennady.vasilyev 2008-02-07 08:59:20 CET
Wrong energy balance for alpha + He

At use of QGSP_BIC_HP (4.9.1 p01) I have received a deviation from conservation of energy for alpha + He:

* G4Track Information:   Particle = alpha,   Track ID = 1,   Parent ID = 0
*********************************************************************************************************

Step#    X(mm)    Y(mm)    Z(mm) KinE(MeV)  dE(MeV) StepLeng TrackLeng  NextVolume ProcName
    0        0        0       10       353        0        0         0       World initStep
    1        0        0     6.34       333     20.3     3.66      3.66       World ionIoni
    2 -0.00376  0.00289        5       325     7.77     1.34         5    Absorber Transportation
    3 -0.00993   0.0096     1.84       306     19.1     3.16      8.16    Absorber ionIoni
    4    -0.01  0.00982    -1.01       287     19.2     2.85        11    Absorber ionIoni
    5 -0.00756  0.00926    -2.03         0     7.48     1.02        12    Absorber alphaInelastic
	:----- List of 2ndaries - #SpawnInStep=  7(Rest= 0,Along= 0,Post= 7), #SpawnTotal=  7 ---------------
	:  -0.00756   0.00926     -2.03       440            neutron
	:  -0.00756   0.00926     -2.03       170             proton
	:  -0.00756   0.00926     -2.03       139            neutron
	:  -0.00756   0.00926     -2.03      92.4             proton
	:  -0.00756   0.00926     -2.03      21.2            neutron
	:  -0.00756   0.00926     -2.03       103           He2[0.0]
	:  -0.00756   0.00926     -2.03       160              gamma

Best regards, G.I.Vasilyev.
Comment 1 Gunter.Folger 2008-02-12 12:26:39 CET
Thanks for reporting this. 
Can you confirm that in your geometry, He is used as absorber. What density do you use?

Gunter
Comment 2 Koichi Murakami 2008-02-21 05:44:58 CET
Hello.

1 cm layer 0.1 He +0.9 H. 1 g/cm2

Gennady.


Comment 3 Gunter.Folger 2011-03-18 15:14:00 CET
Energy/momentum checking is now in place.