Message: Re: radioactive decay in geant4 dedicated examples Not Logged In (login)
 Next-in-Thread Next-in-Thread
 Next-in-Forum Next-in-Forum

More Re: radioactive decay in geant4 dedicated examples 

Forum: Documentation and Examples
Re: Question radioactive decay in geant4 dedicated examples
Re: None Re: radioactive decay in geant4 dedicated examples (michel maire)
Date: 25 Jan, 2013
From: <esther>

I get the same thing with other ions, this is rdecay01 with singleDecay.mac as an input. I get the same thing setting fullChain either to true or false.

Maybe lines regarding energy, momentum, time of life and activity (all either zero/inf/nan) might be useful. As I said in my previous post, the path to RadioactiveDecayData seems to be set correctly, and this is, as an example, the content of file for F-18:

#  18F  (109.77 m)
#  Excitation   Halflife Mode       Daughter Ex  Intensity           Q
#
# File sanitized by Vanderbilt Decay-o-matic
# $Id: repair_decay_files.py,v 1.20 2006/04/24 14:53:36 marcus Exp $ 
# Fri Oct  6 15:26:18 2006 
#
P       0.0000   6.5860e+03
W [VU repair system] Deleted zero-probability event: MshellEC
                          BetaPlus       0.0000   9.6730e-01
                          KshellEC       0.0000   3.0860e-02
                          LshellEC       0.0000   1.8380e-03
                          BetaPlus       0.0000   9.6730e+01    1655.5000
                          KshellEC       0.0000   3.0860e+00    1655.5000
                          LshellEC       0.0000   1.8380e-01    1655.5000

Thanks for your fast reply, Michel.

Esther

OUTPUT ON TERMINAL:

*************************************************************
 Geant4 version Name: geant4-09-04-patch-02    (24-June-2011)
                      Copyright : Geant4 Collaboration
                      Reference : NIM A 506 (2003), 250-303
                            WWW : http://cern.ch/geant4
*************************************************************

/run/verbose 2 /tracking/verbose 2 # /gun/particle ion /gun/ion 10 24 # /rdecay01/fullChain true # /rdecay01/histo/setFileName run1 ###/rdecay01/histo/setFileType hbook /rdecay01/histo/setFileType root /rdecay01/histo/setHisto 1 100 0. 2000 keV ----> SetHisto 1: energy spectrum: e+ e- (keV); 100 bins from 0 keV to 2000 keV /rdecay01/histo/setHisto 2 100 0. 2000 keV ----> SetHisto 2: energy spectrum: nu_e anti_nu_e (keV); 100 bins from 0 keV to 2000 keV /rdecay01/histo/setHisto 3 100 0. 1000 keV ----> SetHisto 3: energy spectrum: gamma (keV); 100 bins from 0 keV to 1000 keV /rdecay01/histo/setHisto 5 100 0. 200 eV ----> SetHisto 5: energy spectrum: ions (eV); 100 bins from 0 eV to 200 eV /rdecay01/histo/setHisto 6 100 -10. 10 keV ----> SetHisto 6: energy balance (keV); 100 bins from -10 keV to 10 keV /rdecay01/histo/setHisto 7 100 0. 1 eV ----> SetHisto 7: momentum balance (eV); 100 bins from 0 eV to 1 eV /rdecay01/histo/setHisto 8 100 0. 70 min ----> SetHisto 8: total time of life of decay chain (min); 100 bins from 0 min to 70 min /rdecay01/histo/setHisto 9 100 0. 70 min ----> SetHisto 9: activity in Bq/g (min); 100 bins from 0 min to 70 min # /rdecay01/event/printModulo 10000 # /run/beamOn 10

Region <DefaultRegionForTheWorld> -- appears in <World> world volume
 Root logical volume(s) : World 
 Pointers : G4VUserRegionInformation[0], G4UserLimits[0], G4FastSimulationManager[0], G4UserSteppingAction[0]
 Materials : G4_AIR 
 Production cuts :   gamma 1 mm      e- 1 mm      e+ 1 mm  proton 1 mm 

========= Table of registered couples ==============================

Index : 0     used in the geometry : Yes     recalculation needed : No 
 Material : G4_AIR
 Range cuts        :  gamma  1 mm     e-  1 mm     e+  1 mm  proton 1 mm 
 Energy thresholds :  gamma  990 eV     e-  990 eV     e+  990 eV  proton 100 keV
 Region(s) which use this couple : 
    DefaultRegionForTheWorld

====================================================================

Start closing geometry. G4GeometryManager::ReportVoxelStats -- Voxel Statistics

    Total memory consumed for geometry optimisation:   0 kByte
    Total CPU time elapsed for geometry optimisation: 0 seconds
Start Run processing.

*********************************************************************************************************
* G4Track Information:   Particle = Ne24[0.0],   Track ID = 1,   Parent ID = 0
*********************************************************************************************************

Step#      X         Y         Z        KineE    dEStep   StepLeng  TrakLeng    Volume     Process
    0      0 fm      0 fm      0 fm      0 eV      0 eV      0 fm      0 fm        World    initStep
    1      0 fm      0 fm      0 fm      0 eV      0 eV      0 fm      0 fm        World   UserLimit

 end of event      0 : Ne24[0.0]

*********************************************************************************************************
* G4Track Information:   Particle = Ne24[0.0],   Track ID = 1,   Parent ID = 0
*********************************************************************************************************

Step#      X         Y         Z        KineE    dEStep   StepLeng  TrakLeng    Volume     Process
    0      0 fm      0 fm      0 fm      0 eV      0 eV      0 fm      0 fm        World    initStep
    1      0 fm      0 fm      0 fm      0 eV      0 eV      0 fm      0 fm        World  RadioactiveDecay

[etc etc]

Run terminated.
Run Summary
  Number of events processed : 10
  User=0.01s Real=0s Sys=0s

 ======================== run summary ======================
 The run was 10 Ne24[0.0] of 0 eV 
 ===========================================================

 Nb of generated particles: 

      Ne24[0.0]:      10  Emean =      0 eV 	( 0 eV  --> 0 eV )

 Energy and momentum balance : final state - initial state
 (excluding gamma desexcitation from momentum balance) : 

  Energy:   mean =   -nan eV 	( 0 eV  --> 0 eV )
  Momentum: mean =   -nan eV 	( 0 eV  --> 0 eV )

 Time of life : mean =      0 ps   half-life =      0 ps    ( 0 ps  --> 0 ps )

 Activity =    inf Bq/g   (inf Ci/g) 

UserDetectorConstruction deleted. UserPhysicsList deleted. UserRunAction deleted. UserPrimaryGenerator deleted. G4 kernel has come to Quit state. EventManager deleted. UImanager deleted. Units table cleared. StateManager deleted. RunManagerKernel is deleted. RunManager is deleting.

 Add Message Add Message
to: "Re: radioactive decay in geant4 dedicated examples"

 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 ]