Warning, /geant4/examples/extended/medical/dna/jetcounter/README is written in an unsupported language. File is not indexed.
0001 -------------------------------------------------------------------
0002
0003 =========================================================
0004 Geant4 - Jet Counter nanodosemeter example
0005 =========================================================
0006
0007 README file
0008 ----------------------
0009
0010 Authors: M. Pietrzak, M. Mietelska, A. Bancer, A. Rucinski and B. Brzozowska
0011 For any question, please contact the developer:
0012 marcin.pietrzak@fuw.edu.pl or marcin.pietrzak@ncbj.gov.pl
0013
0014 This example is provided by the Geant4-DNA collaboration
0015 Any report or published results obtained using the Geant4-DNA software
0016 shall cite the following Geant4-DNA collaboration publications:
0017 [1] Med. Phys. 37 (2010) 4692-4708
0018 [2] Phys. Med. 31 (2015) 861-874
0019 In addition, for this example, please cite:
0020 [3] M. Pietrzak et al 2021 Phys. Med. Biol. 66 225008 (doi: 10.1088/1361-6560/ac33eb).
0021
0022 ----> 0. INTRODUCTION
0023
0024 The jetcounter example provides a setup for simulation of a typical experiment with the Jet Counter nanodosemeter.
0025
0026 This example allows the calculation of the ionisation cluster size distribution (ICSD) (the number of ionisations per projectile)
0027 in a real sized cylindrical interaction volume (IV) of the Jet Counter device.
0028
0029 The results of each of the simulation are provided in the plain text format.
0030 Each thread produces its own output file.
0031
0032 Each row in the output file consist information from a single event, where single primary particle (projectile) was generated.
0033 Each row contains 4 numbers (columns):
0034 1. Ionisation cluster size.
0035 2. Initial energy of the projectile.
0036 3. Interaction energy of the projectile.
0037 4. Final energy of the projectile.
0038
0039 Geometry visualization is available through vis.mac file.
0040
0041 Script written in Python is provided for visualization of the results.
0042
0043
0044 ----> 1. GEOMETRY SET-UP
0045
0046 - The height of the interaction volume (IV) is equal to 20 mm and its diameter is equal to 10 mm.
0047 - The beam axis enters the IV perpendicularly to its axis of rotation, 5 mm from the top.
0048 - The IV is filled with liquid water material (G4_WATER), but its density is set to 0.45 µg/cm³,
0049 which corresponds to the nitrogen gas density in real experimental conditions. For details see [3].
0050 - The beam divergence is restricted by a cylindrical collimator (3 mm in diameter, 23 mm in length) adjacent to the source.
0051 - There is a wall around the IV in the form of a cylinder opened on both ends. The wall thickness is 0.98 mg/cm² (7 µm).
0052 - The wall is composed of two layers:
0053 * outer layer is made of Mylar, 0.96 mg/cm² in thickness;
0054 * inner layer is made of liquid water (G4_WATER), 0.02 mg/cm² in thickness.
0055 - There is a silicon (trigger) detector in the form of a thin disk placed behind the IV.
0056
0057
0058 ----> 2. PRIMARY PARTICLES
0059
0060 If the beam energy is set to any positive value, the beam is mono-energetic with the given energy.
0061 If the beam energy is set to -1 (default value), the energy spectrum file is used.
0062 The energy spectrum of the beam is based on experimental data. Its mean energy is 4.6 MeV,
0063 which is the energy of alpha particles from gold-plated americium (241-Am) source.
0064 This source is chosen in the primaries.in macro file.
0065
0066
0067 ----> 3. PHYSICS LIST
0068
0069 The PhysicsList used in this example show how to combine Geant4-DNA physics and standard physics in different regions
0070 using the G4EmDNAPhysicsActivator.
0071 1. G4EmDNAPhysics is active in the IV and inner layer of the wall (both made of G4_WATER).
0072 2. G4EmLivermorePhysics is active in other parts of the setup.
0073
0074 The G4DNAPhysics_option is chosen in the init_phys.in macro file. Hence, it can be changed without the need of recompilation.
0075
0076
0077 ---->4. SET UP
0078
0079 Make sure G4LEDATA points to the low energy electromagnetic data files.
0080
0081 ---->5. HOW TO RUN THE EXAMPLE
0082
0083 For the visualization:
0084 ./jetcounter vis
0085
0086 For the complete simulation:
0087 ./jetcounter run
0088
0089 In run.mac macro file will read init_phys.in and primaries.in.
0090
0091 However, it is recommended to first increase the number of simulated events in the run.mac file to at least 4000.
0092
0093 To produce figures from simulation results
0094 (requires numpy and matplotlib Python libraries):
0095 python3 plot.py
