Message: bug in geant4 9.2 ? helium counter, energies of secondaries. Not Logged In (login)

Question bug in geant4 9.2 ? helium counter, energies of secondaries.  

Forum: Hadronic Processes
Date: 29 Jan, 2009
From: Marcin Kasztelan <Marcin Kasztelan>

Hello.

Probably I found a bug in new geant4-09-02 (19-December-2008)

*************************************************************

My simulations: I simulate neutron helium counters. In the process n+3He -> p + tritium + 764 keV It worked ok to version 9.1 of geant but now i have wrong energies of secondaries in that reaction. Off course I use the same source code of my program.

-------------------------------------------
My material:
Material: Licznik_gas     density:  2.236 mg/cm3  RadL:  62.954 m    Nucl.Int.Length: 468.405 m    Imean: 180.463 eV   temperature: 294.00 K  pressure:   5.50 atm
   --->  Element: Hel3 (He3)   Z =  2.0   N =   3.0   A =   3.02 g/mole
   --->  Isotope:   He3   Z =  2   N =   3   A =   3.02 g/mole   abundance: 100.00 %  ElmMassFraction:  22.89 %  ElmAbundance  89.19 %
   --->  Element: Krypton (Kr)   Z = 36.0   N =  83.8   A =  83.80 g/mole  ElmMassFraction:  77.11 %  ElmAbundance  10.81 %
--------------------------------------------

And some output from geant4 9.1: ********* BEGIN ******************

*************************************************************
 Geant4 version Name: geant4-09-01-patch-03    (12-September-2008)
                      Copyright : Geant4 Collaboration
                      Reference : NIM A 506 (2003), 250-303
                            WWW : http://cern.ch/geant4
*************************************************************

<<< Geant4 Physics List engine packaging library: PACK 5.4
<<< Geant4 Physics List simulation engine: LHEP_PRECO_HP 3.2

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

Step#    X(mm)    Y(mm)    Z(mm) KinE(MeV)  dE(MeV) StepLeng TrackLeng  NextVolume ProcName
    0     -115      350        0   2.4e-08        0        0         0       World initStep
    1     -100      372    -87.5   2.4e-08        0     91.4      91.4  Licznik_ob Transportation
    2     -100      372    -88.5   2.4e-08        0     1.04      92.5      CENTER Transportation
    3    -98.8      374    -97.1         0        0     9.02       101      CENTER NeutronInelastic
    :----- List of 2ndaries - #SpawnInStep=  2(Rest= 0,Along= 0,Post= 2), #SpawnTotal=  2 ---------------
    :     -98.8       374     -97.1     0.573             proton
    :     -98.8       374     -97.1     0.191             triton
    :----------------------------------------------------------------- EndOf2ndaries Info ---------------

*********************************************************************************************************
* G4Track Information:   Particle = triton,   Track ID = 3,   Parent ID = 1
*********************************************************************************************************

Step#    X(mm)    Y(mm)    Z(mm) KinE(MeV)  dE(MeV) StepLeng TrackLeng  NextVolume ProcName
    0    -98.8      374    -97.1     0.191        0        0         0      CENTER initStep
    1      -99      374    -98.6    0.0212     0.17     1.79      1.79      CENTER hIoni
    2    -98.9      373    -99.2         0   0.0212    0.779      2.57      CENTER hIoni

*********************************************************************************************************
* G4Track Information:   Particle = proton,   Track ID = 2,   Parent ID = 1
*********************************************************************************************************

Step#    X(mm)    Y(mm)    Z(mm) KinE(MeV)  dE(MeV) StepLeng TrackLeng  NextVolume ProcName
    0    -98.8      374    -97.1     0.573        0        0         0      CENTER initStep
    1    -98.7      374    -97.1     0.568  0.00345   0.0525    0.0525      CENTER hIoni
    :----- List of 2ndaries - #SpawnInStep=  1(Rest= 0,Along= 0,Post= 1), #SpawnTotal=  1 ---------------
    :     -98.7       374     -97.1   0.00102                 e-
    :----------------------------------------------------------------- EndOf2ndaries Info ---------------

********* END ******************

and output from newest geant4-09-02 (19-December-2008)

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

Step#    X(mm)    Y(mm)    Z(mm) KinE(MeV)  dE(MeV) StepLeng TrackLeng  NextVolume ProcName
   24     -159      328     -124  1.08e-08        0        0       591     grafit2 initStep
   25     -125      319     -104  1.08e-08        0     40.7       632       World Transportation
   26     -112      315    -95.3  1.08e-08        0     16.2       648  Licznik_ob Transportation
   27     -110      315    -94.3  1.08e-08        0     1.91       650      CENTER Transportation
   28     -105      313    -91.5         0        0     5.45       655      CENTER NeutronInelastic
    :----- List of 2ndaries - #SpawnInStep=  2(Rest= 0,Along= 0,Post= 2), #SpawnTotal=  2 ---------------
    :      -105       313     -91.5     0.954             proton
    :      -105       313     -91.5     0.319             triton
    :----------------------------------------------------------------- EndOf2ndaries Info ---------------

*********************************************************************************************************
* G4Track Information:   Particle = triton,   Track ID = 22,   Parent ID = 1
*********************************************************************************************************

Step#    X(mm)    Y(mm)    Z(mm) KinE(MeV)  dE(MeV) StepLeng TrackLeng  NextVolume ProcName
    0     -105      313    -91.5     0.319        0        0         0      CENTER initStep
    1     -105      312    -91.1     0.132    0.187     2.09      2.09      CENTER hIoni
    2     -104      310    -90.8     0.044   0.0881     1.61       3.7      CENTER hIoni
    3     -104      309    -90.8         0    0.044     1.11      4.81      CENTER hIoni

*********************************************************************************************************
* G4Track Information:   Particle = proton,   Track ID = 21,   Parent ID = 1
*********************************************************************************************************

Step#    X(mm)    Y(mm)    Z(mm) KinE(MeV)  dE(MeV) StepLeng TrackLeng  NextVolume ProcName
    0     -105      313    -91.5     0.954        0        0         0      CENTER initStep
    1     -105      314    -91.5      0.95  0.00246      0.1       0.1      CENTER hIoni
    :----- List of 2ndaries - #SpawnInStep=  1(Rest= 0,Along= 0,Post= 1), #SpawnTotal=  1 ---------------
    :      -105       314     -91.5   0.00198                 e-
    :----------------------------------------------------------------- EndOf2ndaries Info ---------------

********* END ******************

As we see a secondaries have wrong energies:

proton 0.954 MeV - should have 0.573 MeV

tritium 0.319 MeV - should have 0.191 MeV

Can somebody help me with that?

Greetings

Marcin Kasztelan

--------------------------------------------------------------------

Inline Depth:
 0 0
 All All
Outline Depth:
 1 1
 2 2
 All All
Add message: (add)

None Re: bug in geant4 9.2 ? helium counter, energies of secondaries.   by Koi, Tatsumi <Koi, Tatsumi>,   06 Feb, 2009
Re: Question bug in geant4 9.2 ? helium counter, energies of secondaries. (Marcin Kasztelan)
Hi Marcin, Thank you for reporting this problem. The reaction(n-3He->p+T) Q value from Geant4 nuclear masses become 1273keV. If we use masses of AME2003, then the Q value becomes your value 763keV. I did not expect this level difference between them. The reason is clear, however more investigation required for fix. Does some one do full comparison nuclear mass (or atomic mass) of Geant4 to AME2003?

Tatsumi

1 None: Re: bug in geant4 9.2 ? helium counter, energies of secondaries.   (Koi, Tatsumi - 05 Mar, 2009)
(_ None: Re: bug in geant4 9.2 ? helium counter, energies of secondaries.   (Marcin Kasztelan - 05 Mar, 2009)
 Add Message Add Message
to: "bug in geant4 9.2 ? helium counter, energies of secondaries. "

 Subscribe Subscribe

This site runs SLAC HyperNews version 1.11-slac-98, derived from the original HyperNews


[ Geant 4 Home | Geant 4 HyperNews | Search | Request New Forum | Feedback ]