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

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:
 1 1
 All All
Outline Depth:
 1 1
 2 2
 All All
Add message: (add)

1 None: Re: bug in geant4 9.2 ? helium counter, energies of secondaries.   (Koi, Tatsumi - 06 Feb, 2009)
(_ 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 ]