1 | $Id: README,v 1.16 2010/06/03 23:42:16 cirrone Exp $ |
---|
2 | ------------------------------------------------------------------- |
---|
3 | |
---|
4 | ========================================================= |
---|
5 | Geant4 - human_phantom example |
---|
6 | ========================================================= |
---|
7 | |
---|
8 | README |
---|
9 | ----------------------- |
---|
10 | |
---|
11 | Authors: S. Guatelli (susanna.guatelli@uow.edu.au), G. Guerrieri |
---|
12 | and M. G. Pia (pia@ge.infn.it) |
---|
13 | |
---|
14 | Maintained by: S. Guatelli (susanna.guatelli@uow.edu.au) |
---|
15 | |
---|
16 | ------> Introduction |
---|
17 | |
---|
18 | The human_phantom example models anthropomorphic phantoms for |
---|
19 | Geant4 simulations. |
---|
20 | Two models are available: MIRD [1] and ORNL [2] (Male and Female for each approach). |
---|
21 | |
---|
22 | [1] W.S. Snyder, et al, "MIRD Pamphlet No. 5 Revised, Estimates of |
---|
23 | absorbed fractions for monoenergetic photon sources uniformly distributed |
---|
24 | in various organs of a heterogeneous phantom", |
---|
25 | J. Nucl. Med. Suppl., no. 3, pp. 5-52, 1969. |
---|
26 | |
---|
27 | [2] M. Cristy and K. F. Eckerman, "Specific absorbed fractions of energy |
---|
28 | at various ages from internal photon sources", ORNL/TM-8381/VI, Apr. 1987. |
---|
29 | |
---|
30 | The example is based on code developed by the undergraduate student |
---|
31 | G. Guerrieri. |
---|
32 | Note: this is a preliminary beta-version of the code; an improved |
---|
33 | version will be distributed in the next Geant4 public release, compliant |
---|
34 | with the design in a forthcoming publication, and subject to a design |
---|
35 | and code review. |
---|
36 | |
---|
37 | -----> Geometry |
---|
38 | |
---|
39 | The process of building a phantom is handled through the Builder |
---|
40 | design pattern. |
---|
41 | The creation of coherent models of the human phantom is handled through |
---|
42 | an Abstract Factory design pattern. |
---|
43 | |
---|
44 | The organs of the MIRD phantom are implemented in hard-code; the organs of |
---|
45 | the ORNL phantom are handled through GDML (Geometry Description Markup Language, |
---|
46 | www.cern.ch/gdml). The materials of the MIRD model are defined in the |
---|
47 | class G4HumanPhantomMaterial. The materials of the ORNL model are defined in |
---|
48 | the GDML files. |
---|
49 | |
---|
50 | *** Parameterized breast *** |
---|
51 | MIRD Female model: a breast is analytical and derives from the MIRD |
---|
52 | anthropomorphic phantom; the other breast derives from the model [3] and |
---|
53 | it is voxelised. |
---|
54 | |
---|
55 | [3] D.R. Dance and R. A. Hunt, "Voxel breast phantom to represent breasts |
---|
56 | of different sizes and glandularities for use with a Monte Carlo simulation |
---|
57 | program", Report RMTPC 02/1005. |
---|
58 | |
---|
59 | -----> Physics |
---|
60 | |
---|
61 | Particles: e-, e+, photon, geantino |
---|
62 | Physics list: electromagnetic processes are active for photons, e- and e+ |
---|
63 | The threshold of production of secondary particles is set to 0.1 * mm. |
---|
64 | |
---|
65 | -----> Primary particles |
---|
66 | |
---|
67 | The primary particles are geantino. |
---|
68 | UI commands: |
---|
69 | |
---|
70 | /gun/energy energy : to change the energy of primary particles |
---|
71 | /gun/particle particle_type : to change the type of primary particle |
---|
72 | /gun/setBeam beamAlongX : generate a primary particle beam along the X axis |
---|
73 | /gun/setBeam beamAlongY : generate a primary particle beam along the Y axis |
---|
74 | /gun/setBeam beamAlongZ : generate a primary particle beam along the Y axis |
---|
75 | /gun/setBeam isotropicFlux : generate a primary particle with random direction, |
---|
76 | from a random point. |
---|
77 | |
---|
78 | -----> Energy deposit |
---|
79 | |
---|
80 | The energy deposit is calculated in the organs of the phantom. |
---|
81 | At the end of the execution of the simulation the summary of the total energy deposit in |
---|
82 | each organ is print out. |
---|
83 | The user can set the sensitivity of each organ interactively through the UI command: |
---|
84 | /bodypart/addBodyPart organName yes |
---|
85 | (see adultFemale.mac as example) |
---|
86 | |
---|
87 | If the sensitivity is not set ( /bodypart/addBodyPart organName no ), |
---|
88 | the energy deposit is not calculated in the specific organ. |
---|
89 | |
---|
90 | The energy deposit is calculated in each voxel of the parameterised breast. |
---|
91 | |
---|
92 | -----> Analysis |
---|
93 | |
---|
94 | output file: g4humanphantom.hbk containing |
---|
95 | 1D histogram id=10: Particle Path Distribution (cm) |
---|
96 | 2D histogram id=20: Particle Projection XY (cm, cm) |
---|
97 | 2D histogram id=30: Particle Projection YZ (cm, cm) |
---|
98 | 2D histogram id=40: Particle Projection ZX (cm, cm) |
---|
99 | 2D histogram id=50: Energy Deposit (MeV, y axis) in Body Part (x axis) |
---|
100 | 2D histogram id=100: Edep(MeV) in the parameterized breast, x= slice, y= sector |
---|
101 | |
---|
102 | ----> Macro files: example of different human phantoms |
---|
103 | |
---|
104 | default.mac is executed by default in the simulation |
---|
105 | adultMIRDFemale.mac: example to define a MIRD female human phantom |
---|
106 | adultMIRDMale.mac: example to define a MIRD male human phantom |
---|
107 | adultHead.mac: example how to define one piece of the anatomy |
---|
108 | adultORNLFemale.mac: example to define a ORNL female human phantom |
---|
109 | adultORNLMale.mac: example to define a ORNL male human phantom |
---|
110 | adultMIXFemale.mac : example of MIRD human female phantom with |
---|
111 | parameterised breast |
---|
112 | |
---|
113 | -----> Setup |
---|
114 | |
---|
115 | A standard Geant4 example GNUmakefile is provided. |
---|
116 | |
---|
117 | The following environment variables need to be set: |
---|
118 | G4LEDATA : points to low energy data base |
---|
119 | set G4ANALYSIS_USE equal to 1 to activate the analysis |
---|
120 | |
---|
121 | - Requires the GDML module to be installed with Geant4: |
---|
122 | . |
---|
123 | Set G4LIB_BUILD_GDML to 1 |
---|
124 | Set XERCESCROOT variable defining the path to the Xercec-C library |
---|
125 | See also: http://cern.ch/gdml. |
---|
126 | |
---|
127 | - IMPORTANT!!!! Customize the '*.gdml' schema file in the directories |
---|
128 | gdmlData/Female and gdmlData/Male, specifying in the top line the |
---|
129 | right path where to retrieve the main GDML schema file. |
---|
130 | |
---|
131 | - Analysis set-up: AIDA 3.2.1, and iAIDA |
---|
132 | Users can download the analysis tools from: |
---|
133 | http://aida.freehep.org/ |
---|
134 | http://iaida.dynalias.net/ |
---|
135 | |
---|
136 | - Compile and link to generate an executable: |
---|
137 | % gmake |
---|
138 | |
---|
139 | - Execute the application: |
---|
140 | % phantom |
---|
141 | |
---|
142 | - Default macro: adultFemale.mac (MIRD, Female model) |
---|
143 | - Other macros: adultMale.mac (MIRD, Male model) |
---|
144 | adultORNLFemale.mac (ORNL, Female model) |
---|
145 | adultORNLMale.mac (ORNL, Male model) |
---|
146 | - Visualization macros: vrmlVis.mac, dawnVis.mac, openGLVis.mac, noVis.mac |
---|
147 | - Run particles: run.mac |
---|
148 | |
---|
149 | --------------------------------------------------------------------------- |
---|