[807] | 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 | |
---|