1 | ------------------------------------------------------------------- |
---|
2 | $Id: README,v 1.24 2007/06/06 19:15:06 pia Exp $ |
---|
3 | ------------------------------------------------------------------- |
---|
4 | |
---|
5 | ========================================================= |
---|
6 | Geant4 - Radioprotection example |
---|
7 | ========================================================= |
---|
8 | |
---|
9 | README |
---|
10 | --------------------- |
---|
11 | |
---|
12 | 0. INTRODUCTION |
---|
13 | |
---|
14 | The Radioprotection example derives from a Geant4 application |
---|
15 | ( look www.ge.infn.it/geant4/space/remsim for more details ) whose scope |
---|
16 | is to evaluate the dose in astronauts, in vehicle concepts and |
---|
17 | Moon surface habitat configurations, in a defined interplanetary space |
---|
18 | radiation environment. |
---|
19 | |
---|
20 | 1. GEOMETRY |
---|
21 | |
---|
22 | The user can calculate the dose in the astronaut (phantom) in the |
---|
23 | following set-ups: |
---|
24 | |
---|
25 | - Vehicle configuration |
---|
26 | |
---|
27 | | ||sh| | | | | |
---|
28 | |S||ie| |SPE | | | |
---|
29 | ----> |I||ld| |she | |phantom| |
---|
30 | beam |H||in| |lter| | | |
---|
31 | | ||g | | | | | |
---|
32 | |
---|
33 | ---------------------------------------> |
---|
34 | Z axis |
---|
35 | |
---|
36 | - The SIH is the Simplified Inflatable Habitat. |
---|
37 | |
---|
38 | - The shielding is a layer of water, its scope it to protect the astronaut |
---|
39 | from Galactic Cosmic Rays (GCR). The user can add, delete this element |
---|
40 | in the geometrical configuration, change its thickness through UI comands. |
---|
41 | |
---|
42 | - The SPE shelter is a water layer (thickness = 75.cm along Z axis), its scope |
---|
43 | it to protect the astronaut from Solar Particle Events (SPE). |
---|
44 | The user can add, delete this element in the geometrical configuration |
---|
45 | through UI comands. |
---|
46 | |
---|
47 | - The phantom is the astronaut model; the energy deposit is collected in this |
---|
48 | geometrical component. The phantom is a box of water, |
---|
49 | it is 30. cm wide along Z axis, it is voxelised in 30 slices along Z axis. |
---|
50 | The energy deposit of primary and secondary particles is collected in |
---|
51 | each voxel. |
---|
52 | |
---|
53 | - Moon Habitat configuration |
---|
54 | _______________________________ |
---|
55 | | | |
---|
56 | / |Moon Surface | |
---|
57 | / | | |
---|
58 | | x | _________ | |
---|
59 | |<->|| _____ | | |
---|
60 | ---->| || |Phan | | <---shelter | |
---|
61 | beam | || |thom | | | |
---|
62 | | || |_____| | | |
---|
63 | \ ||_________| | |
---|
64 | \ | | |
---|
65 | pyramid | |
---|
66 | log| | |
---|
67 | |______________________________| |
---|
68 | |
---|
69 | -------------------------------------------> |
---|
70 | Z axis |
---|
71 | |
---|
72 | - The astronaut/phantom is set in the astronaut habitat (shelter). |
---|
73 | |
---|
74 | - The pyramid log is made of moon soil and protects the astronaut from |
---|
75 | GCR and SPE. The user can add, delete this element in the geometrical |
---|
76 | configuration, change its thickness (x) through UI comands. |
---|
77 | |
---|
78 | - The phantom is the astronaut model; the energy deposit is collected in this |
---|
79 | geometrical component. The phantom is a box of water, |
---|
80 | it is 30. cm wide along Z axis, it is voxelised in 30 slices along Z axis. |
---|
81 | The energy deposit of primary and secondary particles is collected in |
---|
82 | each voxel. |
---|
83 | |
---|
84 | 1.1 UI |
---|
85 | |
---|
86 | - The user can change the geometry set-up with the following UI commands: |
---|
87 | /configuration/choose vehicle -> choose the Vehicle configuration |
---|
88 | /configuration/choose moon -> choose the Moon Habitat configuration |
---|
89 | The user can not switch between these two configurations interactively. |
---|
90 | |
---|
91 | - The user can select in the vehicle configuration: |
---|
92 | /configuration/AddShielding On -> set the shielding water layer |
---|
93 | /configuration/AddShielding Off -> destroy the shielding water layer |
---|
94 | /shielding/thickness 30.cm -> set the thickness of the shielding layer |
---|
95 | /configuration/AddSPE On -> set the SPE shelter |
---|
96 | /configuration/AddSPE Off -> destroy the SPE shelter |
---|
97 | |
---|
98 | - The user can select in the Moon surface habitat configuration: |
---|
99 | /configuration/AddRoof On -> set the pyramid log |
---|
100 | /configuration/AddRoof Off -> destroy the pyramid log |
---|
101 | /roof/thickness 1. m -> set the height (x) of the pyramid log |
---|
102 | |
---|
103 | 2. PHYSICS LIST |
---|
104 | |
---|
105 | The user can select the physics processes to activate interactively as |
---|
106 | shown in the macro vis.mac. |
---|
107 | The example is provided of: |
---|
108 | - Low Energy electromagnetic processes for photons, e- |
---|
109 | - Standard electromagnetic processes for e+ |
---|
110 | - Low Energy electromagnetic processes with ICRU parameterisation |
---|
111 | for p, alpha and ions |
---|
112 | - Muon electromagnetic processes |
---|
113 | - Decay |
---|
114 | - Hadronic processes for p and alpha particles as primary particles. |
---|
115 | |
---|
116 | The user can choose the Bertini cascade approach or the |
---|
117 | binary cascade approach for the modeling of proton hadronic physics. |
---|
118 | |
---|
119 | 3. PRIMARY PARTICLES |
---|
120 | |
---|
121 | Primary particles are generated according spectra derived from the |
---|
122 | differential flux (CREME 96). |
---|
123 | gcr_min_z=1.txt and gcr_min_z=2.txt contain the differential flux |
---|
124 | of galactic cosmic protons and alpha particles with respect to the |
---|
125 | energy (MeV/nucl). |
---|
126 | These files are read by the primary particle component of the application |
---|
127 | and the spectra are derived. |
---|
128 | Primary particles are generated from a point set in |
---|
129 | the position (0., 0., -25. m), with a direction (0., 0., 1.) by default. |
---|
130 | The user can change these parameters interactively. |
---|
131 | |
---|
132 | 4. STEPPING |
---|
133 | |
---|
134 | Available UI command: |
---|
135 | |
---|
136 | /step/hadronicVerbose On -> print the hadronic processes undertaken by |
---|
137 | particles during the run |
---|
138 | |
---|
139 | /step/hadronicVerbose Off -> switch off the verbose level |
---|
140 | |
---|
141 | 5. ANALYSIS |
---|
142 | |
---|
143 | if ANALYSIS_USE = 1 in the variable environment, the output of |
---|
144 | the simulation is remsim.hbk or remsim.xml. |
---|
145 | The default output file is remsim.hbk. |
---|
146 | The user can change the format of the output file interactively as shown |
---|
147 | in section 5.1 |
---|
148 | |
---|
149 | The file contains histograms: |
---|
150 | - 10 (1) Energy Deposit (MeV)in the phantom (astronaut) versus |
---|
151 | the depth along Z axis |
---|
152 | |
---|
153 | - 20 (1) Initial energy per nucleon (MeV) of primary particles |
---|
154 | |
---|
155 | - 30 (1) Energy Deposit (MeV) in the phantom given by secondaries |
---|
156 | versus the depth along Z axis |
---|
157 | |
---|
158 | - 40 (1) Initial energy (MeV) of primaries reaching the phantom |
---|
159 | |
---|
160 | - 50 (1) Initial energy (MeV) of primaries ougoing the phantom |
---|
161 | |
---|
162 | - 60 (1) Energy (MeV) of primaries reaching the phantom |
---|
163 | |
---|
164 | - 70 (1) Energy (MeV) of primaries outgoing the phantom |
---|
165 | |
---|
166 | - 80 (2) Project the coordinate of the hits on the plane xy |
---|
167 | |
---|
168 | - 90 (2) Project the energy deposit of the hits on the plane xy |
---|
169 | |
---|
170 | - 100 (1) Secondary particles produced in the phantom |
---|
171 | |
---|
172 | - 110 (1) Energy of secondary p produced in the phantom |
---|
173 | |
---|
174 | - 120 (1) Energy of secondary n produced in the phantom |
---|
175 | |
---|
176 | - 130 (1) Energy of secondary pions produced in the phantom |
---|
177 | |
---|
178 | - 140 (1) Energy of secondary alpha produced in the phantom |
---|
179 | |
---|
180 | - 150 (1) Energy of secondary e+ produced in the phantom |
---|
181 | |
---|
182 | - 160 (1) Energy of secondary e- produced in the phantom |
---|
183 | |
---|
184 | - 170 (1) Energy of secondary gamma produced in the phantom |
---|
185 | |
---|
186 | - 180 (1) Energy of secondary mu produced in the phantom |
---|
187 | |
---|
188 | - 190 (1) Energy of secondary other particles produced in the phantom |
---|
189 | |
---|
190 | - 200 (1) Phantom Slice where secondary protons are produced |
---|
191 | |
---|
192 | - 210 (1) Phantom Slice where secondary neutrons are produced |
---|
193 | |
---|
194 | - 220 (1) Phantom Slice where secondary pions are produced |
---|
195 | |
---|
196 | - 230 (1) Phantom Slice where secondary alpha are produced |
---|
197 | |
---|
198 | - 240 (1) Phantom Slice where secondary positrons are produced |
---|
199 | |
---|
200 | - 250 (1) Phantom Slice where secondary electrons are produced |
---|
201 | |
---|
202 | - 260 (1) Phantom Slice where secondary gamma are produced |
---|
203 | |
---|
204 | - 270 (1) Phantom Slice where secondary muons are produced |
---|
205 | |
---|
206 | - 280 (1) Phantom Slice where secondary other particles are produced |
---|
207 | |
---|
208 | - 300 (1) secondary particles reaching the phantom |
---|
209 | |
---|
210 | - 310 (1) Energy of secondary p reaching the phantom |
---|
211 | |
---|
212 | - 320 (1) Energy of secondary n reaching the phantom |
---|
213 | |
---|
214 | - 330 (1) Energy of secondary pions reaching in the phantom |
---|
215 | |
---|
216 | - 340 (1) Energy of secondary alpha reaching the phantom |
---|
217 | |
---|
218 | - 350 (1) Energy of secondary e+ reaching the phantom |
---|
219 | |
---|
220 | - 360 (1) Energy of secondary e- reaching the phantom |
---|
221 | |
---|
222 | - 370 (1) Energy of secondary gamma reaching the phantom |
---|
223 | |
---|
224 | - 380 (1) Energy of secondary mu reaching the phantom |
---|
225 | |
---|
226 | - 390 (1) Energy of secondary other particles reaching the phantom |
---|
227 | |
---|
228 | - 400 (1) Energy of secondary neutrinos produced in the phantom |
---|
229 | |
---|
230 | - 500 (1) secondary particles in the vehicle |
---|
231 | |
---|
232 | 6.SET-UP |
---|
233 | |
---|
234 | - a standard Geant4 example GNUmakefile is provided |
---|
235 | |
---|
236 | setup with: |
---|
237 | compiler = gcc-3.2.3 |
---|
238 | G4SYSTEM = linux-g++ |
---|
239 | |
---|
240 | The following environment variables need to be set for the physics packages: |
---|
241 | G4LEDATA points to low energy data base - |
---|
242 | G4LEVELGAMMADATA points to PhotoEvaporation data |
---|
243 | G4RADIOACTIVEDATA points to Radioactive Decay data |
---|
244 | NeutronHPCrossSections points to neutron data - |
---|
245 | |
---|
246 | Setup for analysis: AIDA 3.2.1, PI 1.3.8 |
---|
247 | |
---|
248 | Users can download the analysis tools from: |
---|
249 | |
---|
250 | http://aida.freehep.org/ |
---|
251 | http://www.cern.ch/PI |
---|
252 | |
---|
253 | 7. HOW TO RUN THE EXAMPLE |
---|
254 | |
---|
255 | example macros are provided: |
---|
256 | |
---|
257 | - vehicle1.mac, vehicle2.mac are examples of simulation in the |
---|
258 | vehicle configuration |
---|
259 | - moon.mac is an example of simulation in the Moon habitat configuration |
---|
260 | |
---|
261 | - batch mode: |
---|
262 | $G4WORDIR/bin/Linux-g++/remsim vehicle1.mac |
---|
263 | $G4WORDIR/bin/Linux-g++/remsim vehicle2.mac |
---|
264 | $G4WORDIR/bin/Linux-g++/Brachy moon.mac |
---|
265 | |
---|
266 | - Interative mode: |
---|
267 | $G4WORDIR/bin/Linux-g++/remsim |
---|
268 | -> the vis.mac is loaded automatically |
---|
269 | |
---|
270 | ------------------------------------------------------------------------ |
---|
271 | |
---|
272 | |
---|