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. 45, (2018) e722-e739
0029 Phys. Med. 31 (2015) 861-874
0030 Med. Phys. 37 (2010) 4692-4708
0031 Int. J. Model. Simul. Sci. Comput. 1 (2010) 157–178
0032
0033 ---->1. GEOMETRY SET-UP.
0034
0035 The geometry is a spherical nucleus surrounded by a spherical shell
0036 representing the cytoplasm. Nucleus radius and shell thickness can be selected
0037 from the provided macro file svalue.in, as well as materials
0038 (G4_WATER or G4_Galactic). The world is a sphere with radius 1e3 larger than
0039 the radius of the nucleus
0040
0041 Particles are shot randomly inside the cytoplasm.
0042
0043 Particle type and energy can be controlled by the
0044 svalue.in macro file.
0045
0046 The PrimaryGeneratorAction class is adapted (G4 state dependent)
0047 in order to enable generic physics list usage
0048 (empty modular physics list).
0049
0050 An alternative MyPrimaryGeneratorActionFromFile class is provided as an example
0051 for reading input data from a file in MT mode for the generation of primaries.
0052 The MyFileReader class is provided as well. They are document in:
0053 https://twiki.cern.ch/twiki/bin/view/Geant4/QuickMigrationGuideForGeant4V10
0054
0055 ---->2. SET-UP
0056
0057 Make sure G4LEDATA points to the low energy electromagnetic data files.
0058
0059 Set the variable MYFILE as 1 in MyFile.hh if you wish to use a spectrum file
0060 of incident energies. By default, this variable is not defined. The file name and
0061 number of lines to read can be specified in MyFile.cc. A spectrum.txt file is provided
0062 as example.
0063
0064 The code can be compiled with cmake.
0065
0066 It works in MT mode.
0067
0068 ---->3. HOW TO RUN THE EXAMPLE
0069
0070 In interactive mode, run:
0071
0072 ./svalue svalue.in
0073
0074 The svalue.in macro allows a full control of the simulation.
0075
0076 Two alternative macros, svalue-iodine125.in and svalue-iodine131.in
0077 are provided for the simulation of radioactive sources.
0078
0079 The svalue-spectrum.in macro shows how to shoot particles
0080 from the spectrum.txt file of energies (unit: eV), which has been created
0081 using the spectrum.C sample ROOT macro.
0082
0083 ---->4. PHYSICS
0084
0085 You can select Geant4-DNA physics constructor in svalue.in.
0086
0087 A tracking cut can be applied if requested.
0088
0089 ---->5. SIMULATION OUTPUT AND RESULT ANALYSIS
0090
0091 The output results consist in a text file (s.txt), containing :
0092 - the radius of the nucleus (in nm)
0093 - the thickness of the cytoplasm (in nm)
0094 - the energy of incident particles (in eV)
0095 - the S value for the cytoplasm (in Gy/Bq.s)
0096 - the rms on S value for the cytoplasm (in Gy/Bq.s)
0097 - the S value for the nucleus (in Gy/Bq.s)
0098 - the rms on S value for the nucleus (in Gy/Bq.s)
0099
0100 One can use the plot.C ROOT macro file to display results.
0101
0102 Note: rms values correspond to standard deviation.
