Warning, /geant4/examples/extended/radioactivedecay/rdecay01/README.md is written in an unsupported language. File is not indexed.
0001 \page Examplerdecay01 Example rdecay01
0002
0003 Survey G4RadioactiveDecay process. See Physics Reference Manual, chapter 36.
0004 See also http://ie.lbl.gov/decay.html
0005
0006 ## Geometry construction
0007
0008 It is a simple box which represente an 'infinite' homogeneous medium.
0009
0010 ## Physics list
0011
0012 PhysicsList defines only G4RadioactiveDecay, G4Transportation processes,
0013 and relevant particle definitions.
0014 Therefore, once created, particles or ions travel as geantino.
0015
0016 ## Primary generator
0017
0018 Default kinematic is an ion (Ne24), at rest, at coordinate origin.
0019 Can be changed with particleGun commands.
0020
0021 ## Physics
0022
0023 As said above, all particles and ions behave as geantino, eg. no energy loss.
0024
0025 A flag:
0026 ```
0027 /rdecay01/fullChain (true or false)
0028 ```
0029 allows to limit to single decay or full decay chain (default).
0030 In case of full decay chain, G4TrackStatus of ions is set to fStopButAlive
0031 in order to force decay at rest.
0032 In case of single decay, G4TrackStatus of secondary ion is set to fStopAndKill.
0033
0034 At each decay, one counts and plots energy spectrum of created particles and
0035 ions, and energy-momentum balance of that decay.
0036
0037 Total time of life of decay chain is plotted. Activity is computed.
0038
0039 The command
0040 ```
0041 /rdecay01/timeWindow
0042 ```
0043 allows to survey activity of each nuclide in a specified
0044 time window [t1,t2] : population at t1 and t2, nb of decays within [t1,t2], mean activity.
0045 See timeWindow.mac
0046
0047 Few macros are given in example. Debug.mac is to be run in interactive mode.
0048
0049 ## User data files
0050
0051 Users can redefine RadioactiveDecay and PhotonEvaporation data, via commands:
0052 ```
0053 /grdm/setRadioactiveDecayFile
0054 /grdm/setPhotoEvaporationFile
0055 ```
0056 Examples of such files are given in subdirectory UserData. \n
0057 Formats are described in readme \n
0058 Examples in macros Cf238.mac and No252.mac
0059
0060 ## example of biasing
0061
0062 macro timeWindowBiased.mac illustrates one of the biasing capabilities of the
0063 radioactiveDecay package. Ca47 is forced to decay within 20 days
0064 (eg. 1728000 seconds in data file).
0065 It is instructive to plot time of life (histo 8) with and without the weight of the track:
0066 see line 189/190 of TrackingAction.cc
0067 And also to compare with analog decay mode : comment out /grdm/ commands in the macro.
0068 ## Visualisation
0069
0070 Visualization Manager is set in main().
0071 Initialisation of the drawing is done via the commands
0072 /vis/.. in the macro vis.mac. This macro is automatically read from the main
0073 in case of interactive running mode.
0074
0075 - e- red
0076 - e+ blue
0077 - nu_e white
0078 - anti_nu_e white
0079 - gamma green
0080 - alpha yellow
0081 - GenericIon grey
0082
0083 ## How to start ?
0084
0085 - Execute rdecay01 in 'batch' mode from macro files
0086 ```
0087 % ./rdecay01 singleDecay.mac
0088 % ./rdecay01 rdecay01.in > rdecay01.out
0089 ```
0090
0091 - Execute rdecay01 in 'interactive mode' with visualization
0092 ```
0093 % ./rdecay01
0094 ....
0095 Idle> ---> type your commands. For instance:
0096 Idle> /control/execute debug.mac
0097 ....
0098 Idle> /run/beamOn 1
0099 ....
0100 Idle> exit
0101 ```
0102
0103 Macros provided in this example:
0104 - Co60.mac: Co60
0105 - Gd158.mac: Gd158 excited state
0106 - No158.mac: read user data file
0107 - Po212.mac: Po212 excited state
0108 - Ra228.mac: Ra228 excited state
0109 - alpha.mac: Po212 alpha decay
0110 - atomicDeexcitation.mac: plot Auger cascade
0111 - fullChain.mac: U238
0112 - neutron.mac: Li10 neutron emission
0113 - proton.mac: Co53 proton emission
0114 - singleDecay.mac: Ne10
0115 - timeWindow.mac: print activity within a given time window
0116 - timeWindowBiased.mac: force decay within a given time window
0117
0118 Macros to be run interactively:
0119 - Cf238.mac: read user data file
0120 - debug.mac: Pb210
0121 - electronicCapture.mac: Fe55 electronic capture
0122 - vis.mac: To activate visualization
0123
0124 ## Histograms
0125
0126 rdecay01 produces several 1D histograms which are saved as
0127 rdecay01.root by default.
0128
0129 - 1 : energy spectrum: e+ e-
0130 - 2 : energy spectrum: nu_e anti_nu_ev
0131 - 3 : energy spectrum: gamma
0132 - 4 : energy spectrum: alpha
0133 - 5 : energy spectrum: ions
0134 - 6 : total kinetic energy (Q)
0135 - 7 : momentum balance
0136 - 8 : total time of life of decay chain
0137 - 9 : total visible energy
0138
0139 The histograms are managed by G4AnalysiManager and its Messenger.
0140 The histos can be individually activated with the command :
0141 ```
0142 /analysis/h1/set id nbBins valMin valMax unit
0143 ```
0144 where unit is the desired unit for the histo (MeV or keV, deg or mrad, etc..)
0145
0146 One can control the name of the histograms file with the command:
0147 ```
0148 /analysis/setFileName name (default rdecay1)
0149 ```
0150
0151 It is possible to choose the format of the histogram file : root (default),
0152 xml, csv, by using namespace in HistoManager.hh
0153
0154 It is also possible to print selected histograms on an ascii file:
0155 ```
0156 /analysis/h1/setAscii id
0157 ```
0158 All selected histos will be written on a file name.ascii (default rdecay1)
