Warning, /geant4/examples/extended/medical/dna/svalue/README is written in an unsupported language. File is not indexed.
0001 =========================================================
0002 Geant4 - svalue example
0003 =========================================================
0004
0005 README file
0006 ----------------------
0007
0008 CORRESPONDING AUTHOR
0009
0010 S. Incerti (a, *)
0011 a. LP2i, IN2P3 / CNRS / Bordeaux University, 33175 Gradignan, France
0012 * e-mail: incerti@lp2ib.in2p3.fr
0013
0014 ---->0. INTRODUCTION.
0015
0016 The svalue example shows how to calculate S values in liquid water
0017 using the Geant4-DNA physics processes and models.
0018
0019 It is adapted from TestEm12.
0020
0021 This example is provided by the Geant4-DNA collaboration.
0022
0023 These processes and models are further described at:
0024 http://geant4-dna.org
0025
0026 Any report or published results obtained using the Geant4-DNA software shall
0027 cite the following Geant4-DNA collaboration publications:
0028 Med. Phys. 51 (2024) 5873–5889
0029 Med. Phys. 45 (2018) e722-e739
0030 Phys. Med. 31 (2015) 861-874
0031 Med. Phys. 37 (2010) 4692-4708
0032 Int. J. Model. Simul. Sci. Comput. 1 (2010) 157–178
0033
0034 ---->1. GEOMETRY SET-UP.
0035
0036 The geometry is a spherical nucleus surrounded by a spherical shell
0037 representing the cytoplasm. Nucleus radius and shell thickness can be selected
0038 from the provided macro file svalue.in, as well as materials
0039 (G4_WATER or G4_Galactic). The world is a sphere with radius 1e3 larger than
0040 the radius of the nucleus
0041
0042 Particles are shot randomly inside the cytoplasm.
0043
0044 Particle type and energy can be controlled by the
0045 svalue.in macro file.
0046
0047 The PrimaryGeneratorAction class is adapted (G4 state dependent)
0048 in order to enable generic physics list usage
0049 (empty modular physics list).
0050
0051 An alternative MyPrimaryGeneratorActionFromFile class is provided as an example
0052 for reading input data from a file in MT mode for the generation of primaries.
0053 The MyFileReader class is provided as well. They are document in:
0054 https://twiki.cern.ch/twiki/bin/view/Geant4/QuickMigrationGuideForGeant4V10
0055
0056 ---->2. SET-UP
0057
0058 Make sure G4LEDATA points to the low energy electromagnetic data files.
0059
0060 Set the variable MYFILE as 1 in MyFile.hh if you wish to use a spectrum file
0061 of incident energies. By default, this variable is not defined. The file name and
0062 number of lines to read can be specified in MyFile.cc. A spectrum.txt file is provided
0063 as example.
0064
0065 The code can be compiled with cmake.
0066
0067 It works in MT mode.
0068
0069 ---->3. HOW TO RUN THE EXAMPLE
0070
0071 In interactive mode, run:
0072
0073 ./svalue svalue.in
0074
0075 The svalue.in macro allows a full control of the simulation.
0076
0077 Two alternative macros, svalue-iodine125.in and svalue-iodine131.in
0078 are provided for the simulation of radioactive sources.
0079
0080 The svalue-spectrum.in macro shows how to shoot particles
0081 from the spectrum.txt file of energies (unit: eV), which has been created
0082 using the spectrum.C sample ROOT macro.
0083
0084 ---->4. PHYSICS
0085
0086 You can select Geant4-DNA physics constructor in svalue.in.
0087
0088 A tracking cut can be applied if requested.
0089
0090 ---->5. SIMULATION OUTPUT AND RESULT ANALYSIS
0091
0092 The output results consist in a text file (s.txt), containing :
0093 - the radius of the nucleus (in nm)
0094 - the thickness of the cytoplasm (in nm)
0095 - the energy of incident particles (in eV)
0096 - the S value for the cytoplasm (in Gy/Bq.s)
0097 - the rms on S value for the cytoplasm (in Gy/Bq.s)
0098 - the S value for the nucleus (in Gy/Bq.s)
0099 - the rms on S value for the nucleus (in Gy/Bq.s)
0100
0101 One can use the plot.C ROOT macro file to display results.
0102
0103 Note: rms values correspond to standard deviation.
