Geant4 Cross Reference |
1 ========================================= 1 ========================================================= 2 Geant4 - wvalue example 2 Geant4 - wvalue example 3 ========================================= 3 ========================================================= 4 4 5 README file 5 README file 6 -------------------- 6 ---------------------- 7 7 8 CORRESPONDING AUTHO << 8 CORRESPONDING AUTHOR 9 9 10 S. Incerti (a, *) << 10 S. Incerti et al. (a, *) 11 a. LP2i, IN2P3 / CNRS / Bordeaux University, 3 << 11 a. Centre d'Etudes Nucleaires de Bordeaux-Gradignan 12 * e-mail: incerti@lp2ib.in2p3.fr << 12 (CENBG), IN2P3 / CNRS / Bordeaux University, 33175 Gradignan, France 13 << 13 * e-mail:incerti@cenbg.in2p3.fr 14 ---->0. INTRODUCTION. << 15 14 >> 15 ---->0. INTRODUCTION. >> 16 16 The wvalue example shows how to calculate w in 17 The wvalue example shows how to calculate w in liquid water 17 for e- using the Geant4-DNA physics processes << 18 for e- using the Geant4-DNA physics processes and models. 18 19 19 w is computed as the ratio of the incident par 20 w is computed as the ratio of the incident particle energy 20 and the total number of ionisations. 21 and the total number of ionisations. 21 22 22 It is adapted from the svalue example. << 23 It is adapted from svalue. >> 24 >> 25 *** It is a preliminary version which might be buggy *** 23 26 24 This example is provided by the Geant4-DNA col 27 This example is provided by the Geant4-DNA collaboration. 25 28 26 These processes and models are further describ 29 These processes and models are further described at: 27 http://geant4-dna.org 30 http://geant4-dna.org 28 31 29 Any report or published results obtained using << 32 Any report or published results obtained using the Geant4-DNA software shall 30 cite the following Geant4-DNA collaboration pu 33 cite the following Geant4-DNA collaboration publications: 31 Med. Phys. 51 (2024) 5873–5889 << 32 Med. Phys. 45 (2018) e722-e739 << 33 Phys. Med. 31 (2015) 861-874 34 Phys. Med. 31 (2015) 861-874 34 Med. Phys. 37 (2010) 4692-4708 35 Med. Phys. 37 (2010) 4692-4708 35 Int. J. Model. Simul. Sci. Comput. 1 (2010) 15 << 36 36 37 This example is presented in the following pap 37 This example is presented in the following paper, which shall also be cited: 38 Med. Phys. 42 (2015) 3870-3876 << 38 Med. Phys. 42 (2015) 3870-3876 39 39 40 ---->1. GEOMETRY SET-UP. 40 ---->1. GEOMETRY SET-UP. 41 << 41 42 The geometry is a 1 m radius sphere of liquid 42 The geometry is a 1 m radius sphere of liquid water (G4_WATER 43 material). Particles are shot randomly from th 43 material). Particles are shot randomly from the sphere centre. 44 44 45 Radius of the sphere, physics constructor and << 45 Radius of the sphere, physics constructor and energy can be 46 controlled by the wvalue.in macro file. 46 controlled by the wvalue.in macro file. 47 47 48 The PrimaryGeneratorAction class is adapted (G << 48 The PrimaryGeneratorAction class is adapted (G4 state dependent) 49 in order to enable generic physics list usage << 49 in order to enable generic physics list usage 50 (empty modular physics list). 50 (empty modular physics list). 51 51 52 ---->2. SET-UP << 52 ---->2. SET-UP 53 << 53 54 Make sure G4LEDATA points to the low energy el 54 Make sure G4LEDATA points to the low energy electromagnetic data files. 55 55 56 The code can be compiled with cmake. 56 The code can be compiled with cmake. 57 57 58 It works in MT mode. 58 It works in MT mode. 59 59 60 ---->3. HOW TO RUN THE EXAMPLE << 60 ---->3. HOW TO RUN THE EXAMPLE 61 61 62 In interactive mode, run: 62 In interactive mode, run: 63 63 64 ./wvalue wvalue.in 64 ./wvalue wvalue.in 65 65 66 The wvalue.in macro allows a full control of t 66 The wvalue.in macro allows a full control of the simulation. 67 67 68 ---->4. PHYSICS 68 ---->4. PHYSICS 69 69 70 You can select Geant4-DNA physics constructor 70 You can select Geant4-DNA physics constructor in wvalue.in. 71 71 72 A tracking cut can be applied if requested. 72 A tracking cut can be applied if requested. 73 73 74 ---->5. SIMULATION OUTPUT AND RESULT ANALYSIS << 74 ---->5. SIMULATION OUTPUT AND RESULT ANALYSIS 75 75 76 The output results consist in a text file (wva << 76 The output results consist in a text file (wvalue.txt), containing : 77 - the energy of incident particles (in eV) 77 - the energy of incident particles (in eV) 78 - the mean number of ionisations 78 - the mean number of ionisations 79 - its rms 79 - its rms 80 - the w value (in eV) 80 - the w value (in eV) 81 - its rms (in eV) 81 - its rms (in eV) 82 82 83 Note: rms values correspond to standard deviat 83 Note: rms values correspond to standard deviation. 84 84 85 In addition, another macro (histo.in) is also << 85 In addition, another macro (histo.in) is also provided including 86 a series of histograms: << 86 a series of histograms : 87 - histogram #1 : nb of ionisation interactions 87 - histogram #1 : nb of ionisation interactions per event 88 - histogram #2 : total energy deposited in abs 88 - histogram #2 : total energy deposited in absorber 89 - histogram #3 : true track length of the prim 89 - histogram #3 : true track length of the primary particle 90 - histogram #4 : true step size of the primary 90 - histogram #4 : true step size of the primary particle 91 - histogram #5 : projected range of the primar 91 - histogram #5 : projected range of the primary particle 92 - histogram #6 : true track length of charged 92 - histogram #6 : true track length of charged secondaries 93 - histogram #7 : true track length of charged 93 - histogram #7 : true track length of charged secondaries >> 94 >> 95 --------------------------------------------------------------------------- >> 96 >> 97 Should you have any enquiry, please do not hesitate to contact: >> 98 incerti@cenbg.in2p3.fr >> 99