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
3 =========================================================
4 Geant4 - an Object-Oriented Toolkit for Simulation in HEP
5 =========================================================
6
7 icsd
8 -------
9
10 Authors: Sylvain Meylan, Yann Perrot and Carmen Villagrasa (IRSN, France)
11 Marcin Pietrzak (NCBJ, Poland)
12 For any question, please contact:
13 carmen.villagrasa@irsn.fr
14
15 This example is provided by the Geant4-DNA collaboration
16 Any report or published results obtained using the Geant4-DNA software
17 shall cite the following Geant4-DNA collaboration publications:
18 [1] Med. Phys. 37 (2010) 4692-4708
19 [2] Phys. Med. 31 (2015) 861-874
20 In addition, for this example using DNA-materials cross-sections, please cite:
21 [3] Rad. Phys and Chem. 130 (2017) 459-479
22 For this example using N2 cross-sections, please cite:
23 [4] Phys. Med. 102 (2022) 103-109
24
25 ---->0. INTRODUCTION.
26
27 The icsd name of this example stands for: Ionisation Cluster Size Distribution.
28 This example allows the calculation of the distribution concerning the number of ionisations per event in an small cylinder of 2.3 nanometers diameter and 3.4 nanometers heigh, typical dimensions of a 10 base pairs piece of chromatin.
29 It constitutes the first use of new DNA-like material's cross-sections that were experimentally obtained at PTB (Rad. Phys and Chem. 130 (2017) 459-479) for electrons and calculated for protons on THF, TMP, PY and PU materials.
30 The new physical models corresponding to the DNA materials a,d nitrogen have recently been implemented in Geant4-DNA thaks to new classes, that allow the use of other materials different from liquid water: G4VDNAModel, G4DNAModelInterface,G4DNADummyModel and other associated classes.
31
32 These new models are further described at:
33 http://geant4-dna.org
34
35 ---->1. GEOMETRY SET-UP AND PRIMARY PARTICLES
36
37 As indicated in the introduction, the geometry used in this exercise is a cylinder.
38 Two options are proposed: one nanometric cylinder made of DNA or a macroscopic volume with low nitrogen density.
39 The selection of the geometry is performed in the macro file:
40
41 --> /icsd/setGeom dna: a cylinder of 2.3 nanometers diameter and 3.4 nanometers height, typical dimensions of a 10 base pairs piece of chromatin that is centered in a cubic world of 10 nanometers side length filled with liquid water.
42 This cylinder is filled with THF material. The definition of the THF material must be done using the G4String="THF". Currently, other materials are available than also use the G4Strings to be recognized by the model classes "TMP", "PY" or "PU" for the DNA-like materials and "deoxyribose","adenine", "guanine", "thymine" or "cytosine" for the derived DNA materials.
43
44 --> /icsd/setGeom nanodosimeter: a cylinder of 10 mm nanometers diameter and 10 nanometers height filled with low pressure nitrogen (0.34 µg/cm3) mimicking the size of a short of segment od DNA of 2.3 nm in both diameter and height
45 The obtained results can be compared to experimental data from Jet Counter nanodosimeter as shown in Pietrazk et al.
46 Phys. Med. 102 (2022) 103-109
47
48 Primary particles are electrons starting at the cylinder's edge) and initial momentum (1,0,0)
49
50 ---->2. PHYSICS LIST
51
52 The PhysicsList used in this example show how to use the new classes that allow Geant4-DNA to use other materials than liquid water :
53 1/ Models are built: classical Geant4-DNA models for liquid water material and PTB models for DNA-like materials
54 2/ G4DNAModelInterface object is built for each type of interaction (elastic, ionisation, excitation)
55 3/ G4DNAModels are registered within the corresponding G4DNAModelInterface. For liquid water models, the particle name must be added to the RegisterModel function. For DNA-like material models and Vacuum no other parameter is needed.
56 4/ G4DNAProcesses objects are built and the corresponding G4DNAModelInterface is registered within. G4DNAModelInterface acts like a classical G4VEmModel
57
58 Remember that new DNA-like material classes are available for:
59
60 electrons from 12 eV-1 keV, including elastic, ionisation and excitation processes
61 protons from 70 keV to 10 MeV, only for ionisation processes following a modified HKS formalism
62
63 ---->3. SET UP
64
65 Make sure G4LEDATA points to the low energy electromagnetic data files.
66
67 ---->4. HOW TO RUN THE EXAMPLE
68
69 ./icsd will run in multi-threaded mode with 2 threads by default.
70 Two ntuples will be filled in the output file:
71 - Ntuple 1 contains the information for calculating the ionisation cluster size distribution per event
72 - Ntuple 2 contains interaction information at the step level
73
