Version:
[ ReleaseNotes ] [ 1.0 ] [ 1.1 ] [ 2.0 ] [ 3.0 ] [ 3.1 ] [ 3.2 ] [ 4.0 ] [ 4.0.p1 ] [ 4.0.p2 ] [ 4.1 ] [ 4.1.p1 ] [ 5.0 ] [ 5.0.p1 ] [ 5.1 ] [ 5.1.p1 ] [ 5.2 ] [ 5.2.p1 ] [ 5.2.p2 ] [ 6.0 ] [ 6.0.p1 ] [ 6.1 ] [ 6.2 ] [ 6.2.p1 ] [ 6.2.p2 ] [ 7.0 ] [ 7.0.p1 ] [ 7.1 ] [ 7.1.p1 ] [ 8.0 ] [ 8.0.p1 ] [ 8.1 ] [ 8.1.p1 ] [ 8.1.p2 ] [ 8.2 ] [ 8.2.p1 ] [ 8.3 ] [ 8.3.p1 ] [ 8.3.p2 ] [ 9.0 ] [ 9.0.p1 ] [ 9.0.p2 ] [ 9.1 ] [ 9.1.p1 ] [ 9.1.p2 ] [ 9.1.p3 ] [ 9.2 ] [ 9.2.p1 ] [ 9.2.p2 ] [ 9.2.p3 ] [ 9.2.p4 ] [ 9.3 ] [ 9.3.p1 ] [ 9.3.p2 ] [ 9.4 ] [ 9.4.p1 ] [ 9.4.p2 ] [ 9.4.p3 ] [ 9.4.p4 ] [ 9.5 ] [ 9.5.p1 ] [ 9.5.p2 ] [ 9.6 ] [ 9.6.p1 ] [ 9.6.p2 ] [ 9.6.p3 ] [ 9.6.p4 ] [ 10.0 ] [ 10.0.p1 ] [ 10.0.p2 ] [ 10.0.p3 ] [ 10.0.p4 ] [ 10.1 ] [ 10.1.p1 ] [ 10.1.p2 ] [ 10.1.p3 ] [ 10.2 ] [ 10.2.p1 ] [ 10.2.p2 ] [ 10.2.p3 ] [ 10.3 ] [ 10.3.p1 ] [ 10.3.p2 ] [ 10.3.p3 ] [ 10.4 ] [ 10.4.p1 ] [ 10.4.p2 ] [ 10.4.p3 ] [ 10.5 ] [ 10.5.p1 ] [ 10.6 ] [ 10.6.p1 ] [ 10.6.p2 ] [ 10.6.p3 ] [ 10.7 ] [ 10.7.p1 ] [ 10.7.p2 ] [ 10.7.p3 ] [ 10.7.p4 ] [ 11.0 ] [ 11.0.p1 ] [ 11.0.p2 ] [ 11.0.p3, ] [ 11.0.p4 ] [ 11.1 ] [ 11.1.1 ] [ 11.1.2 ] [ 11.1.3 ] [ 11.2 ] [ 11.2.1 ] [ 11.2.2 ] [ 11.3.0 ]
1 =========================================================
2 Geant4 - svalue example
3 =========================================================
4
5 README file
6 ----------------------
7
8 CORRESPONDING AUTHOR
9
10 S. Incerti (a, *)
11 a. LP2i, IN2P3 / CNRS / Bordeaux University, 33175 Gradignan, France
12 * e-mail: incerti@lp2ib.in2p3.fr
13
14 ---->0. INTRODUCTION.
15
16 The svalue example shows how to calculate S values in liquid water
17 using the Geant4-DNA physics processes and models.
18
19 It is adapted from TestEm12.
20
21 This example is provided by the Geant4-DNA collaboration.
22
23 These processes and models are further described at:
24 http://geant4-dna.org
25
26 Any report or published results obtained using the Geant4-DNA software shall
27 cite the following Geant4-DNA collaboration publications:
28 Med. Phys. 51 (2024) 5873–5889
29 Med. Phys. 45 (2018) e722-e739
30 Phys. Med. 31 (2015) 861-874
31 Med. Phys. 37 (2010) 4692-4708
32 Int. J. Model. Simul. Sci. Comput. 1 (2010) 157–178
33
34 ---->1. GEOMETRY SET-UP.
35
36 The geometry is a spherical nucleus surrounded by a spherical shell
37 representing the cytoplasm. Nucleus radius and shell thickness can be selected
38 from the provided macro file svalue.in, as well as materials
39 (G4_WATER or G4_Galactic). The world is a sphere with radius 1e3 larger than
40 the radius of the nucleus
41
42 Particles are shot randomly inside the cytoplasm.
43
44 Particle type and energy can be controlled by the
45 svalue.in macro file.
46
47 The PrimaryGeneratorAction class is adapted (G4 state dependent)
48 in order to enable generic physics list usage
49 (empty modular physics list).
50
51 An alternative MyPrimaryGeneratorActionFromFile class is provided as an example
52 for reading input data from a file in MT mode for the generation of primaries.
53 The MyFileReader class is provided as well. They are document in:
54 https://twiki.cern.ch/twiki/bin/view/Geant4/QuickMigrationGuideForGeant4V10
55
56 ---->2. SET-UP
57
58 Make sure G4LEDATA points to the low energy electromagnetic data files.
59
60 Set the variable MYFILE as 1 in MyFile.hh if you wish to use a spectrum file
61 of incident energies. By default, this variable is not defined. The file name and
62 number of lines to read can be specified in MyFile.cc. A spectrum.txt file is provided
63 as example.
64
65 The code can be compiled with cmake.
66
67 It works in MT mode.
68
69 ---->3. HOW TO RUN THE EXAMPLE
70
71 In interactive mode, run:
72
73 ./svalue svalue.in
74
75 The svalue.in macro allows a full control of the simulation.
76
77 Two alternative macros, svalue-iodine125.in and svalue-iodine131.in
78 are provided for the simulation of radioactive sources.
79
80 The svalue-spectrum.in macro shows how to shoot particles
81 from the spectrum.txt file of energies (unit: eV), which has been created
82 using the spectrum.C sample ROOT macro.
83
84 ---->4. PHYSICS
85
86 You can select Geant4-DNA physics constructor in svalue.in.
87
88 A tracking cut can be applied if requested.
89
90 ---->5. SIMULATION OUTPUT AND RESULT ANALYSIS
91
92 The output results consist in a text file (s.txt), containing :
93 - the radius of the nucleus (in nm)
94 - the thickness of the cytoplasm (in nm)
95 - the energy of incident particles (in eV)
96 - the S value for the cytoplasm (in Gy/Bq.s)
97 - the rms on S value for the cytoplasm (in Gy/Bq.s)
98 - the S value for the nucleus (in Gy/Bq.s)
99 - the rms on S value for the nucleus (in Gy/Bq.s)
100
101 One can use the plot.C ROOT macro file to display results.
102
103 Note: rms values correspond to standard deviation.
