Geant4 Cross Reference

Cross-Referencing   Geant4
Geant4/examples/extended/radioactivedecay/rdecay01/README

Version: [ ReleaseNotes ] [ 1.0 ] [ 1.1 ] [ 2.0 ] [ 3.0 ] [ 3.1 ] [ 3.2 ] [ 4.0 ] [ 4.0.p1 ] [ 4.0.p2 ] [ 4.1 ] [ 4.1.p1 ] [ 5.0 ] [ 5.0.p1 ] [ 5.1 ] [ 5.1.p1 ] [ 5.2 ] [ 5.2.p1 ] [ 5.2.p2 ] [ 6.0 ] [ 6.0.p1 ] [ 6.1 ] [ 6.2 ] [ 6.2.p1 ] [ 6.2.p2 ] [ 7.0 ] [ 7.0.p1 ] [ 7.1 ] [ 7.1.p1 ] [ 8.0 ] [ 8.0.p1 ] [ 8.1 ] [ 8.1.p1 ] [ 8.1.p2 ] [ 8.2 ] [ 8.2.p1 ] [ 8.3 ] [ 8.3.p1 ] [ 8.3.p2 ] [ 9.0 ] [ 9.0.p1 ] [ 9.0.p2 ] [ 9.1 ] [ 9.1.p1 ] [ 9.1.p2 ] [ 9.1.p3 ] [ 9.2 ] [ 9.2.p1 ] [ 9.2.p2 ] [ 9.2.p3 ] [ 9.2.p4 ] [ 9.3 ] [ 9.3.p1 ] [ 9.3.p2 ] [ 9.4 ] [ 9.4.p1 ] [ 9.4.p2 ] [ 9.4.p3 ] [ 9.4.p4 ] [ 9.5 ] [ 9.5.p1 ] [ 9.5.p2 ] [ 9.6 ] [ 9.6.p1 ] [ 9.6.p2 ] [ 9.6.p3 ] [ 9.6.p4 ] [ 10.0 ] [ 10.0.p1 ] [ 10.0.p2 ] [ 10.0.p3 ] [ 10.0.p4 ] [ 10.1 ] [ 10.1.p1 ] [ 10.1.p2 ] [ 10.1.p3 ] [ 10.2 ] [ 10.2.p1 ] [ 10.2.p2 ] [ 10.2.p3 ] [ 10.3 ] [ 10.3.p1 ] [ 10.3.p2 ] [ 10.3.p3 ] [ 10.4 ] [ 10.4.p1 ] [ 10.4.p2 ] [ 10.4.p3 ] [ 10.5 ] [ 10.5.p1 ] [ 10.6 ] [ 10.6.p1 ] [ 10.6.p2 ] [ 10.6.p3 ] [ 10.7 ] [ 10.7.p1 ] [ 10.7.p2 ] [ 10.7.p3 ] [ 10.7.p4 ] [ 11.0 ] [ 11.0.p1 ] [ 11.0.p2 ] [ 11.0.p3, ] [ 11.0.p4 ] [ 11.1 ] [ 11.1.1 ] [ 11.1.2 ] [ 11.1.3 ] [ 11.2 ] [ 11.2.1 ] [ 11.2.2 ] [ 11.3.0 ] 

  1 -------------------------------------------------------------------
  2 
  3      =========================================================
  4      Geant4 - an Object-Oriented Toolkit for Simulation in HEP
  5      =========================================================
  6 
  7                             rdecay01 
  8                             --------
  9           
 10  Survey G4RadioactiveDecay process. See Physics Reference Manual, chapter 36.
 11  See also http://ie.lbl.gov/decay.html
 12   
 13  1- Geometry construction
 14     ---------------------
 15 
 16  It is a simple box which represente an 'infinite' homogeneous medium.
 17   
 18  2- Physics list
 19     ------------
 20 
 21  PhysicsList.cc defines only G4RadioactiveDecay, G4Transportation processes,
 22  and relevant particle definitions.
 23  Therefore, once created, particles or ions travel as geantino.
 24           
 25  3- Primary generator
 26     -----------------
 27  
 28  Default kinematic is an ion (Ne24), at rest, at coordinate origin.
 29  Can be changed with particleGun commands.
 30           
 31  4- Physics
 32     -------
 33  As said above, all particles and ions behave as geantino, eg. no energy loss.
 34      
 35  A flag, /rdecay01/fullChain (true or false), allows to limit to
 36  single decay or full decay chain (default). 
 37  In case of full decay chain, G4TrackStatus of ions is set to fStopButAlive
 38  in order to force decay at rest.   
 39  In case of single decay, G4TrackStatus of secondary ion is set to fStopAndKill.
 40    
 41  At each decay, one counts and plots energy spectrum of created particles and
 42  ions, and energy-momentum balance of that decay.
 43  
 44  Total time of life of decay chain is plotted. Activity is computed.
 45  
 46  The command /rdecay01/timeWindow allows to survey activity of each nuclide in a specified
 47  time window [t1,t2] : population at t1 and t2, nb of decays within [t1,t2], mean activity.
 48  See timeWindow.mac
 49  
 50  Few macros are given in example. Debug.mac is to be run in interactive mode. 
 51  
 52  4-a User data files
 53  -------------------
 54  Users can redefine RadioactiveDecay and PhotonEvaporation data, via commands:
 55  /grdm/setRadioactiveDecayFile
 56  /grdm/setPhotoEvaporationFile 
 57  Examples of such files are given in subdirectory UserData. 
 58  Formats are described in readme
 59  Examples in macros Cf238.mac and No252.mac
 60  
 61  4-b example of biasing
 62  ----------------------
 63  macro timeWindowBiased.mac illustrates one of the biasing capabilities of the 
 64  radioactiveDecay package. Ca47 is "forced" to decay within 20 days
 65  (eg. 1728000 seconds in data file).
 66  It is instructive to plot time of life (histo 8) with and without the weight of the track:
 67  see lines 189/190 of TrackingAction.cc
 68  And also to compare with analog decay mode : comment out /grdm/ commands in the macro. 
 69     
 70  5- Visualisation
 71     -------------
 72  
 73  Visualization Manager is set in the main().
 74  Initialisation of the drawing is done via the commands
 75  /vis/.. in the macro vis.mac. This macro is automatically read from the main 
 76  in case of interactive running mode.
 77  
 78  e- red
 79  e+ blue
 80  nu_e white
 81  anti_nu_e white
 82  gamma green
 83  alpha yellow
 84  GenericIon grey
 85  
 86  6- How to start ?
 87     --------------
 88  
 89   - Execute rdecay01 in 'batch' mode from macro files
 90   % rdecay01  singleDecay.mac
 91   % rdecay01  rdecay01.in > rdecay01.out  
 92     
 93   - Execute rdecay01 in 'interactive mode' with visualization
 94   % rdecay01
 95   ....
 96   Idle>   ---> type your commands. For instance:
 97   Idle> /control/execute debug.mac 
 98   ....
 99   Idle> /run/beamOn 1
100   ....        
101   Idle> exit
102   
103  Macros provided in this example:
104   - Co60.mac:  Co60
105   - Gd158.mac: Gd158 excited state
106   - No158.mac: read user data file
107   - Po212.mac: Po212 excited state
108   - Ra228.mac: Ra228 excited state
109   - alpha.mac: Po212 alpha decay
110   - atomicDeexcitation.mac: plot Auger cascade
111   - fullChain.mac: U238
112   - neutron.mac: Li10 neutron emission
113   - proton.mac:  Co53 proton  emission
114   - singleDecay.mac: Ne10
115   - timeWindow.mac: print activity within a given time window
116   - timeWindowBiased.mac: force decay within a given time window  
117     
118  Macros to be run interactively:
119   - Cf238.mac: read user data file
120   - debug.mac: Pb210
121   - electronicCapture.mac: Fe55 electronic capture
122   - vis.mac: To activate visualization
123   
124  7- Histograms
125     ----------
126  
127   rdecay01 produces several 1D histograms which are saved as
128   rdecay01.root by default.
129 
130       1 : energy spectrum: e+ e-
131       2 : energy spectrum: nu_e anti_nu_ev
132       3 : energy spectrum: gamma
133       4 : energy spectrum: alpha
134       5 : energy spectrum: ions
135       6 : total kinetic energy (Q)
136       7 : momentum balance
137       8 : total time of life of decay chain
138       9 : total visible energy
139                             
140    The histograms are managed by G4AnalysisManager and its Messenger. 
141    The histos can be individually activated with the command :
142    /analysis/h1/set id nbBins  valMin valMax unit 
143    where unit is the desired unit for the histo (MeV or keV, deg or mrad, etc..)
144    
145    One can control the name of the histograms file with the command:
146    /analysis/setFileName  name  (default rdecay1)
147    
148    It is possible to choose the format of the histogram file : root (default),
149    xml, csv, by using namespace in HistoManager.hh
150    
151    It is also possible to print selected histograms on an ascii file:
152    /analysis/h1/setAscii id
153    All selected histos will be written on a file name.ascii (default rdecay1)
154