| 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 |
|
|---|
