Geant4 Cross Reference

Cross-Referencing   Geant4
Geant4/examples/extended/medical/dna/jetcounter/

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 ]

Name Size       Last modified (GMT) Description
Back Parent directory       2024-12-05 15:16:16
Folder include/       2024-12-05 15:16:16
Folder src/       2024-12-05 15:16:16
File CMakeLists.txt 869 bytes       2024-12-05 15:16:16
File History 739 bytes       2024-12-05 15:16:16
File README 4182 bytes       2024-12-05 15:16:16
File energy.spectrum 12864 bytes       2024-12-05 15:16:16
File init_phys.in 155 bytes       2024-12-05 15:16:16
C++ file jetcounter.cc 3257 bytes       2024-12-05 15:16:16
File jetcounter.out 41430 bytes       2024-12-05 15:16:16
File plot.py 2640 bytes       2024-12-05 15:16:16
File primaries.in 88 bytes       2024-12-05 15:16:16
File run.mac 288 bytes       2024-12-05 15:16:16
File vis.mac 433 bytes       2024-12-05 15:16:16

  1 -------------------------------------------------------------------
  2 
  3      =========================================================
  4         Geant4 - Jet Counter nanodosemeter example
  5      =========================================================
  6 
  7                             README file
  8                       ----------------------
  9 
 10 Authors: M. Pietrzak, M. Mietelska, A. Bancer, A. Rucinski and B. Brzozowska
 11 For any question, please contact the developer:
 12 marcin.pietrzak@fuw.edu.pl or marcin.pietrzak@ncbj.gov.pl
 13 
 14 This example is provided by the Geant4-DNA collaboration
 15 Any report or published results obtained using the Geant4-DNA software  
 16 shall cite the following Geant4-DNA collaboration publications:
 17 [1] Med. Phys. 37 (2010) 4692-4708
 18 [2] Phys. Med. 31 (2015) 861-874
 19 In addition, for this example, please cite:
 20 [3] M. Pietrzak et al 2021 Phys. Med. Biol. 66 225008 (doi: 10.1088/1361-6560/ac33eb).
 21 
 22 ----> 0. INTRODUCTION
 23 
 24 The jetcounter example provides a setup for simulation of a typical experiment with the Jet Counter nanodosemeter.
 25 
 26 This example allows the calculation of the ionisation cluster size distribution (ICSD) (the number of ionisations per projectile)
 27 in a real sized cylindrical interaction volume (IV) of the Jet Counter device.
 28 
 29 The results of each of the simulation are provided in the plain text format.
 30 Each thread produces its own output file.
 31 
 32 Each row in the output file consist information from a single event, where single primary particle (projectile) was generated.
 33 Each row contains 4 numbers (columns):
 34     1. Ionisation cluster size.
 35     2. Initial energy of the projectile.
 36     3. Interaction energy of the projectile.
 37     4. Final energy of the projectile.
 38 
 39 Geometry visualization is available through vis.mac file.
 40 
 41 Script written in Python is provided for visualization of the results.
 42 
 43 
 44 ----> 1. GEOMETRY SET-UP
 45 
 46 - The height of the interaction volume (IV) is equal to 20 mm and its diameter is equal to 10 mm.
 47 - The beam axis enters the IV perpendicularly to its axis of rotation, 5 mm from the top.
 48 - The IV is filled with liquid water material (G4_WATER), but its density is set to 0.45 µg/cm³,
 49 which corresponds to the nitrogen gas density in real experimental conditions. For details see [3].
 50 - The beam divergence is restricted by a cylindrical collimator (3 mm in diameter, 23 mm in length) adjacent to the source.
 51 - 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).
 52 - The wall is composed of two layers:
 53     * outer layer is made of Mylar, 0.96 mg/cm² in thickness;
 54     * inner layer is made of liquid water (G4_WATER), 0.02 mg/cm² in thickness.
 55 - There is a silicon (trigger) detector in the form of a thin disk placed behind the IV.
 56 
 57 
 58 ----> 2. PRIMARY PARTICLES
 59 
 60 If the beam energy is set to any positive value, the beam is mono-energetic with the given energy.
 61 If the beam energy is set to -1 (default value), the energy spectrum file is used.
 62 The energy spectrum of the beam is based on experimental data. Its mean energy is 4.6 MeV,
 63 which is the energy of alpha particles from gold-plated americium (241-Am) source.
 64 This source is chosen in the primaries.in macro file.
 65 
 66 
 67 ----> 3. PHYSICS LIST
 68 
 69 The PhysicsList used in this example show how to combine Geant4-DNA physics and standard physics in different regions
 70 using the G4EmDNAPhysicsActivator.
 71 1. G4EmDNAPhysics is active in the IV and inner layer of the wall (both made of G4_WATER).
 72 2. G4EmLivermorePhysics is active in other parts of the setup.
 73 
 74 The G4DNAPhysics_option is chosen in the init_phys.in macro file. Hence, it can be changed without the need of recompilation.
 75 
 76 
 77 ---->4. SET UP
 78 
 79 Make sure G4LEDATA points to the low energy electromagnetic data files.
 80 
 81 ---->5. HOW TO RUN THE EXAMPLE
 82 
 83  For the visualization:
 84     ./jetcounter vis
 85 
 86  For the complete simulation:
 87     ./jetcounter run
 88 
 89  In run.mac macro file will read init_phys.in and primaries.in.
 90 
 91  However, it is recommended to first increase the number of simulated events in the run.mac file to at least 4000.
 92 
 93  To produce figures from simulation results
 94  (requires numpy and matplotlib Python libraries):
 95     python3 plot.py