Back to home page

EIC code displayed by LXR

 
 

    


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