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