Warning, /geant4/examples/extended/radioactivedecay/rdecay02/README.md is written in an unsupported language. File is not indexed.
0001 \page Examplerdecay02 Example rdecay02
0002
0003 Rdecay02 is created to show how to use the G4RadioactiveDecay process to
0004 simulate the decays of radioactive isotopes as well as the induced
0005 radioactivity resulted from nuclear interactions.
0006
0007 In this example a simple geometry consists of a cylindric target placed
0008 in the centre of a tube detector. Various primary event generation and
0009 tallying options are available.
0010
0011 ## GEOMETRY
0012
0013 The world is filled with "Air" and there are two components in it:
0014
0015 - Target: A cylinder placed at the origin along the z-axis. The default
0016 size of the cylinder is 0.5 cm radius and 1 cm length, and its default
0017 material is "CsI".
0018
0019 - Detector: A tube centered at the origin along the z-axis, with inner
0020 radius matching the radius of the target. The default thickness of the
0021 tube is 2 cm and it is 5 cm long. The default material is "Germanium".
0022
0023 The user can change the target/detector size and material, using the
0024 commands in the directory
0025 ```
0026 /rdecay02/det
0027 ```
0028
0029 ## PHYSICS
0030
0031 The following physics processes are included by default:
0032
0033 - Standard electromagnetic
0034 - Decay
0035 - Radioactive Decay
0036 - Hadronic processes
0037
0038 By default radioactive decay is applied through out the geometry.
0039 The user can limit it to just the target by commands:
0040 ```
0041 /process/had/rdm/noVolumes
0042 /process/had/rdm/selectVolume Target
0043 ```
0044
0045
0046 ## EVENT: THE PRIMARY GENERATOR
0047
0048 The primary kinematic is a single particle or ion shooted at the
0049 centre of the target. The type of the particle and its energy are set in
0050 PrimaryGeneratorAction, and can be changed via the G4 build-in commands of
0051 ParticleGun class (see the macros provided with this example).
0052 Default is Ne24, at rest.
0053
0054 ## DETECTOR RESPONSE
0055
0056 The relevant informations are collected in TrackingAction or
0057 SteppingAction. These include:
0058
0059 - Emission particles in the RadioactiveDecay process:
0060 - particle PDGcode,
0061 - particle kinetic energy,
0062 - particle creation time,
0063 - particle weight.
0064 Note: the residual nuclei is not considered as an emitted particle.
0065
0066 - Radio-Isotopes. All the radioactive isotopes produced in the simulation:
0067 - isotope PDGcode,
0068 - isotope creation time,
0069 - isotope weight.
0070
0071 - Energy depositions in the target and detector by products of the
0072 RadioactiveDecay process:
0073 - energy depostion (positive value for target and negative for detector),
0074 - time,
0075 - weight.
0076
0077
0078 ## HISTOGRAMS
0079
0080 The test contains 7 built-in 1D histograms, which are managed by
0081 G4AnalysisManager and its Mmssenger. The histos can be individually
0082 activated with the command :
0083 ```
0084 /analysis/h1/set id nbBins valMin valMax unit
0085 ```
0086 where unit is the desired unit for the histo (MeV or keV, etc..)
0087 (see the macros xxxx.mac).
0088
0089 - histogram 0: The Pulse Height Spectrum (PHS) of the target.
0090 - histogram 1: The PHS of the detector.
0091 - histogram 2: The combined PHS of the target and detector.
0092 - histogram 3: The anti-coincidece PHS of the target.
0093 - histogram 4: The anti-coincidece PHS of the detector.
0094 - histogram 5: The coincidece PHS between the target and detector.
0095 - histogram 6: The emitted particle energy spectrum.
0096
0097 It is assumed the detector and target pulses both have an integration time
0098 of 1 microsecond, and the gate is 2 microsecond for the coincidence spectrum.
0099 The target and detctor have a threshold of 10 keV in the anti-/coincidence
0100 modes.
0101
0102 HistoManager includes also 4 ntuples whose contents are described in the above paragraphe
0103 (detector response)
0104 The ntuples can be activated with the command
0105 ```
0106 /analysis/ntuple/setActivation
0107 ```
0108
0109 One can control the name of the analysis file with the command:
0110 ```
0111 /analysis/setFileName name (default rdecay02)
0112 ```
0113
0114 It is possible to choose the format of the histogram file : root (default),
0115 xml, csv, by using namespace in HistoManager.hh
0116
0117 It is also possible to print selected histograms on an ascii file:
0118 ```
0119 /analysis/h1/setAscii id
0120 ```
0121 All selected histos will be written on a file name.ascii (default rdecay02)
0122
0123 Initially, all histograms but histogram 6 are inactive. They can all be turned on
0124 with the command
0125 ```
0126 /analysis/h1/setActivationToAll true
0127 ```
0128 or specific histograms can be turned on with the command
0129 ```
0130 /analysis/h1/setActivation i true
0131 ```
0132 where i is the histogram index (0,... n).
0133 To turn off, set the final argument to false
0134
0135 ## VISUALIZATION
0136
0137 The Visualization Manager is set in the main().
0138 The initialisation of the drawing is done via the commands
0139 /vis/... in the macro vis.mac. To get visualisation:
0140 ```
0141 > /control/execute vis.mac
0142 ```
0143 The tracks are drawn at the end of event, and erased at the end of run.
0144 - gamma green
0145 - neutron yellow
0146 - negative particles (e-, ...) red
0147 - positive particles (e+, ions, ...) blue
0148
0149 ## HOW TO START ?
0150
0151 Execute rdecay02 in 'batch' mode from macro files :
0152 ```
0153 % ./rdecay02 run.mac
0154 ```
0155
0156 Execute rdecay02 in 'interactive mode' with visualization :
0157 ```
0158 % ./rdecay02
0159 Idle> control/execute debug.mac
0160 ....
0161 Idle> type your commands
0162 ....
0163 Idle> exit
0164 ```
0165
0166 - run.mac : decay of Ne24. A run of 1000 events
0167 - debug.mac: interactively. One Ne24 decay,
0168 with visualization and tracking/verbose
0169
0170 ## FURTHER EXAMPLES
0171
0172 There are a number of macros files in the ./macros subdirectory, to show
0173 the features of the G4RadioactiveDecay process. Most of them will lead to
0174 the creation of a Root file in the same name of the macro file.
0175
0176 u238c.mac: shows the decays of the U238 chain in analogue MC mode.
0177
0178 th234c-b.mac: shows the decays of Th234 in variance reduction MC mode.
0179 All its secondaies in along the decay chains are generated. The default
0180 source profile and decay biasing schemes are used to determine the decay
0181 times and weights of the secondaries.
0182
0183 proton.mac: simulation of 1 GeV protons incident on a lead target.
0184 The decays of the radio-siotopes created in the proton-lead interactions
0185 are simulated with RadioactiveDecay in analogue MC mode.
0186
0187 proton-beam.mac: same as proton.mac, but the decays of the radio-siotopes
0188 created in the proton-lead interactions are simulated with
0189 RadioactiveDecay in variance reduction MC mode. The isotopes and those
0190 along the decay chains are forced to decay in the time windows specified
0191 by the user in file measures.data, and the weights of the decay products
0192 are determined by the beam profile as defined in the beam.data file and
0193 their decay times.
0194
0195 neutron.mac: macrofile to show the incident of low energy neutrons on an
0196 user specified NaI target and the decays of the induced radio-isotopes.
0197
0198 ne24.mac: this shows the decays of Ne-24 to Na-24 in variance reduction MC
0199 mode. Further decays of Na-24 are not simulated by applying the
0200 nucleuslimits in RadioactiveDecay. Two runs are carried out.
0201 One with the bracjing ratio biasing applied and one without.
0202
0203 isotopes.mac: to show the decays of a number of different isotopes in a
0204 single macro file.
0205
0206 f24.mac: to show the different treatments one can apply to the decays of F24.
0207 i) the complete decay chain from F24 to Mg24, in analogue mode;
0208 ii) the complete chain, but in variance reduction mode;
0209 iii) restrict to the decay of F24 only in analogue mode; iv) restrict to
0210 the decay of F24 only but in variance reduction mode.
0211
0212 as74.mac: The decays of As74 which has a rather complicated decay scheme.
0213 i) in analogue MC mode;
0214 ii) in variance reduction MC mode.
0215
0216 UserRadDataPb210Test.mac: show how the user can define its own radioactive
0217 decay datafile
0218
0219 UserEvapDataBiTest.mac: show how the user can define its own
0220 photo-evaporation datafile
0221
0222 No252Test.mac: show how to simulate Radoactive decay for nuclei with Z>100
0223 based on user datafile