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