| Geant4 Cross Reference |
1 =========================================
2 Geant4 - ICRP110_HumanPhan
3 =========================================
4
5 The ICRP110_HumanPhantoms example is developed
6 Centre For Medical Radiation Physics (CMRP), U
7 and University of Manchester, UK.
8
9 Contacts:
10 - susanna@uow.edu.au
11 - mjl970@uowmail.edu.au
12 - malaroda@uow.edu.au
13 - John.Allison@g4ai.org
14
15 The example is based on the extended/medical/D
16
17 The authors acknowledge that this application
18 the International Commission on Radiological P
19
20 ----------------------------------------------
21 --------------------------------------> Introd
22 ----------------------------------------------
23
24 This application models the ICRP110 reference
25 the dose in individual voxels and in entire or
26
27 The human male phantom, provided kindly by the
28 individual with height 176 cm and mass approxi
29 whole body CT images of a 43yr old individual
30 individuals laying supine and with arms restin
31 closely approximate the ICRP adult Reference M
32
33 [1] HG Menzel, C Clement, and P DeLuca. ICRP
34 an icrp/icru joint effort: A report of adult
35 ICRP, 39(2):1, 2009. URL: http://www.icrp.or
36
37 [2] Valetin J 2002 Basic anatomical and phys
38 reference values: ICRP Publication 89 Ann. I
39
40 [3] Valetin J 2007 The 2007 recommendations
41 protection Ann. ICRP vol 37 (Oxford: Elsevie
42
43 The table below summarises the key features of
44
45 PROPERTY AM AF
46 _____________________________________
47 Height (m) 1.76 1.63
48
49 Mass(Kg) 73.0 60.0
50
51 Slice Thickness(mm) 8.0 4.84
52
53 Voxel in-plane- 2.137 1.775
54 -resolution (mm)
55
56 Voxels along x 254 299
57 (i.e. columns)
58
59 Voxels along y 127 137
60 (i.e. rows)
61
62 Number of Slices 222 348
63 (i.e. along z)
64 ______________________________________
65
66 ----------------------------------------------
67 ------------------------------> Application Su
68 ----------------------------------------------
69
70 - '/src': where the source .cc files are stor
71
72 - '/include': where header .hh files are stor
73
74 - '/ICRPdata': where the phantom data files (
75 It is downloaded automatically from URL https
76 during the configuration via cmake.
77
78 Phantom data files containing the voxelisatio
79 containing the definitions of the phantom org
80 code can be found in the folder /ICRPdata.
81
82 All data files used for this phantom were obt
83 - https://www.icrp.org/publication.asp
84
85 ----------------------------------------------
86 ----------------------------------> ICRP110Pha
87 ----------------------------------------------
88
89 Within the '/ICRPdata' directory, the followin
90
91 -> /ICRPdata/ : conta
92 simul
93
94 -> /ICRPdata/ICRP110_g4dat/AM/ : conta
95
96 -> /ICRPdata/ICRP110_g4dat/AF/ : conta
97
98 -> /ICRPdata/ICRP110_g4dat/P110_data_V
99
100 The final directory contains the raw ICRP110 p
101 5 files within folders for the AM and AF phant
102 supplementary data's included README file.
103
104 The array of organ identification numbers (in
105 AM.dat
106 AF.dat
107
108 A list of individually segmented structures,
109 AM_organs.dat
110 AF_organs.dat
111
112 A list of the media, their elemental composi
113 the file names are:
114 AM_media.dat
115 AF_media.dat
116
117 The mass ratios of bone constituents (trabec
118 the file names are:
119 AM_spongiosa.dat
120 AF_spongiosa.dat
121
122 The mass ratios of blood in various body tis
123 AM_blood.dat
124 AF_blood.dat
125
126 The primary data files AM.dat and AF.dat conta
127 Each number respresents the organ associated w
128 are listed slice by slice, within each slice r
129 index changes fastest, then the row index, the
130 along x, then along y and finally along z. Sli
131 row numbers increase from front to back; and c
132
133 For use in this application, the original AM.d
134 all voxels of the phantom were sub-divided int
135 As such, each file represents a 2D phantom sli
136 of each voxel. This allows for subsections of
137 need to simulate the entire phantom every time
138 allow for reductions in the simulation time de
139 This feature was achieved via a code developed
140
141 The AM human phantom is voxelised in x,y,z wit
142 The AF human phantom is voxelised in x,y,z wit
143
144 ----------------------------------------------
145 ---------------------------------------> How t
146 ----------------------------------------------
147
148 - Create a build folder for the phantom run
149 % mkdir build/
150
151 - Navigate to inside the build folder and init
152 % cmake ../
153
154 The ICRP110 phantom data will be automatical
155
156 - Compile and link to generate the executable
157 % make
158 This should make two executables - ICRP110ph
159
160 - Execute the application in 'interactive' mod
161 % ./ICRP110phantoms
162
163 - Execute the "standalone" application in 'int
164 % ./ICRP110standalone
165 This allows you to visualise the phantom wit
166 Of course, you cannot run or visualise traje
167
168 - Execute the application in 'batch' mode from
169 % ./ICRP110phantoms female_head.in
170
171 -----------------------------
172 AVAILABLE MACRO FILES
173 -----------------------------
174 For the users convenience, macro files have be
175 and trunk phantoms for both the male and femal
176 mode when executing the application as specifi
177 male or female phantom, the macros male.in and
178
179 - male_head.in/female_head.in : Creates a p
180 - male_trunk.in/female_trunk.in : Creates a p
181 - male.in : Creates ful
182 if the user
183 - female.in : Creates ful
184 'ICRPdata/Fe
185 - openGLVis.mac : macro for v
186 - vis.mac (default) : Executed by
187 - primary.mac : Contains th
188
189 At the very top of the various '.in' macro fil
190 which define the sex and section of the phanto
191
192 o /phantom/setPhantomSex <option> : Passes s
193 o /phantom/setScoreWriterSex <option> : Pass
194
195 o /phantom/setPhantomSection <option> : Pass
196 o /phantom/setScoreWriterSection <option> P
197
198 Available options for the first 2 commands are
199 Avalable options for the last 2 commands are:
200
201 In the event that the macro called upon by the
202 does not contain these commands (default case)
203
204 WARNING: the phantom model can be chosen only
205 It cannot be changed during the run session. T
206
207 ----------------------------------------------
208 ----------------------------------> Creating a
209 ----------------------------------------------
210
211 If the user wishes to construct a customised s
212 he/she has to create a specific macro or edit
213 phantom is outlined as follows.
214
215 The user should edit the macro 'male.in' and t
216 'MaleData.dat'. Firstly, in 'FemaleData.dat',
217 select a custom range of phantom slices to sim
218
219 1. The very first entry of each Data.dat indic
220 Changing this number will determine the num
221
222 2. Further down in the Data.dat files (beginni
223 by successive slices. Changing the slice fi
224 phantom to be simulated. As an indication t
225 male phantom below.
226
227 --> AM_Slice1.g4dat to AM_Slice20.g4dat: F
228
229 --> AM_Slice21.g4dat to AM_Slice121.g4dat:
230
231 --> AM_Slice169.g4dat: Single chest slice
232 of lungs, ribs, hea
233
234 --> AM_Slice182.g4dat to AM_Slice222.g4dat
235
236 NOTE: o Always order phantom slices beg
237 in slice number going down the
238 o Always use consecutive/adjacent
239 o The default number of slices fo
240 and starts at the feet of each
241
242 Once the user customises the MaleData.dat/Fema
243 he/she has also to fix appropriately the scori
244
245 ----------------------------------------------
246 ------------------------------> Scoring Mesh a
247 ----------------------------------------------
248
249 The macro primary.mac defines the radiation be
250 General Particle Source should be used to chan
251 the /run/beamOn command and can call upon the
252 '/control/execute primary.mac'.
253
254 Within male.in and female.in, a scoring mesh i
255 phantom voxel. The size of the scoring mesh is
256 to match the constructed phantom dimensions (w
257
258 The mesh dimensions are defined as half-dimens
259 a scoring mesh spanning from -100mm to +100mm
260 male phantom which has dimensions along x,y,z
261 as 271.399 x 135.6995 x 888. mm. The number of
262 should match the number of phantom voxels in x
263 directory.
264
265 If the user edits the MaleData.dat or FemaleDa
266 the scoring mesh dimensions and number of bins
267 only have to edit lines 54 and 55 of the male.
268
269 After completion of a simulation run, the phan
270 "PhantomMesh_Dose.txt". This text file lists t
271
272 The output PhantomMesh_Dose.txt file is create
273 in the voxels is analysed and associated to or
274
275 A final output file "ICRP.out" is then created
276
277 ----------------------------------------------
278 ----------------------------------------> Furt
279 ----------------------------------------------
280
281 -------> ColourMap.dat <--------
282
283 This file located in the build directory assig
284 The user may edit these as they wish for visua
285
286 ----------> Physics <-----------
287
288 The QGSP_BIC_HP Physics List is adopted. The u
289 cut of production of secondary particles.
290
291 -----> Primary particles <------
292
293 The G4 General Particle Source (gps) is used t
294 Macro primary.mac contains the definition of t