[807] | 1 | ------------------------------Advanced Example--------------------------------- README FILE |
---|
| 2 | |
---|
| 3 | |
---|
| 4 | NB: The documentation for this example is in the process of being updated. |
---|
| 5 | Please go to http://www.ge.infn.it/geant4/examples/ for the most up-to-date |
---|
| 6 | description, manual, design documentation and users requirement/traceability. |
---|
| 7 | |
---|
| 8 | JUNE 2002: |
---|
| 9 | The underground physics example has been updated since the December Geant4.0 |
---|
| 10 | release to include: 1) Analysis, using AIDA 2.2 and tested against Anaphe 4.0 |
---|
| 11 | 2) Full lab geometry - important for neutron scattering |
---|
| 12 | 3) Updated macros/messengers |
---|
| 13 | 4) Correction of a few minor features |
---|
| 14 | |
---|
| 15 | The messengers can be seen from typing help within the UI. |
---|
| 16 | Additional file handling and time/energy cut messengers have been implemented. |
---|
| 17 | |
---|
| 18 | (1) and (2) are discussed more fully below |
---|
| 19 | |
---|
| 20 | Note: Due to the importation of data files during the initialisation stage of |
---|
| 21 | Geant4, load-time may be in excess of 5 minutes. |
---|
| 22 | |
---|
| 23 | UNDERGROUND PHYSICS |
---|
| 24 | |
---|
| 25 | An example of a underground dark matter experiment. |
---|
| 26 | |
---|
| 27 | Please see the UserRequirements.txt and related web-pages referred to |
---|
| 28 | at the end of that document. |
---|
| 29 | |
---|
| 30 | Over-view: |
---|
| 31 | A single liquid xenon cell is simulated within Geant4 and the scintillation |
---|
| 32 | light produced from interactions from various calibration species is recorded |
---|
| 33 | as PhotoMultiplier hits. The output is then written to an ASCII file for |
---|
| 34 | future off-line analysis. |
---|
| 35 | |
---|
| 36 | Geometry: |
---|
| 37 | Experimental set-up: |
---|
| 38 | |
---|
| 39 | G4double worldWidth = 470.0*cm + 2.*wallThick; // "x" |
---|
| 40 | G4double worldLength = 690.0*cm + 2.*wallThick; // "y" |
---|
| 41 | G4double worldHeight = 280.0*cm + 2.*wallThick; // "z" |
---|
| 42 | |
---|
| 43 | |
---|
| 44 | |
---|
| 45 | A "cavern" of dimensions 5.18m x 7.38m x 3.28m with concrete walls is defined |
---|
| 46 | as the World Volume. A laboratory geometry is incorporated included desks, |
---|
| 47 | cupboards, door and windows. For ease this is included in a separate ".icc" |
---|
| 48 | file which can be removed should the code seem cumbersome. In the centre of |
---|
| 49 | the cavern a steel vacuum vessel containing liquid and gaseous xenon is placed. |
---|
| 50 | The internal construction of the vessel accurately reproduces an existing |
---|
| 51 | prototype Dark Matter detector which allows experimental comparison. The active |
---|
| 52 | detector volume is defined by a series of metal rings, complemented by |
---|
| 53 | a cover mirror and a PMT immersed in the liquid. Two grids and a |
---|
| 54 | thermalising copper shield are also incorporated. The liquid/gas |
---|
| 55 | interface is located 6mm away from the mirror surface. A Am241 |
---|
| 56 | calibration source is suspended from one of the grids in the liquid |
---|
| 57 | phase, above the PMT. |
---|
| 58 | |
---|
| 59 | XXX================XXX mirror |
---|
| 60 | XXX________________XXX gas phase |
---|
| 61 | XXX XXX |
---|
| 62 | XXX XXX liquid phase |
---|
| 63 | XXX XXX |
---|
| 64 | XXX.......U........XXX grid + calibrator |
---|
| 65 | XXX................XXX grid |
---|
| 66 | XXX| |XXX |
---|
| 67 | | ___------___ | |
---|
| 68 | || PMT || |
---|
| 69 | || || |
---|
| 70 | |
---|
| 71 | |
---|
| 72 | Hits Output (file "hits.out"): |
---|
| 73 | An ASCII file containing the following information: |
---|
| 74 | Evt # : event number |
---|
| 75 | Etot, MeV : energy deposited in liquid xenon |
---|
| 76 | LXe hits : number of hits in liquid xenon |
---|
| 77 | LXeTime, ns : time of first hit in liquid xenon |
---|
| 78 | PMT hits : number of hits in PMT (photocathode) |
---|
| 79 | PmtTime, ns : average PMT hit time relative to LXeTime |
---|
| 80 | First hit : first particle to hit liquid xenon |
---|
| 81 | Flags : particles contributing to energy deposition |
---|
| 82 | Seeds : the initial seed values for given hit events |
---|
| 83 | |
---|
| 84 | Note: |
---|
| 85 | |
---|
| 86 | The time information PmtTime is erroneous when forcing nuclear |
---|
| 87 | decay with the RadioactiveDecay module due to the precision |
---|
| 88 | required to detect nanosecond scintillation times on a global |
---|
| 89 | time of 432 years (in the case of 241Am decay). |
---|
| 90 | |
---|
| 91 | |
---|
| 92 | Pmt Output (file "pmt.out"): |
---|
| 93 | Photon hit positions within the PMT face (overwritten every event): |
---|
| 94 | "Hit# X, mm Y, mm Z, mm" |
---|
| 95 | |
---|
| 96 | |
---|
| 97 | To Run: |
---|
| 98 | Either run the macro files interactively or in batch with the command |
---|
| 99 | DMX macro_name.mac. |
---|
| 100 | |
---|
| 101 | |
---|
| 102 | Macros: |
---|
| 103 | |
---|
| 104 | initInter.mac |
---|
| 105 | Initialisation macro for interactive mode. |
---|
| 106 | |
---|
| 107 | gamma.mac |
---|
| 108 | Shoots one 60 keV gamma upwards from the calibrator and traces the |
---|
| 109 | scintillation light produced in LXe to the PMT. All tracks are drawn |
---|
| 110 | with custom colours. PMT hits in photocathode are also shown. Event |
---|
| 111 | summary is writen to file "hits.out". |
---|
| 112 | |
---|
| 113 | gamma_1000.mac |
---|
| 114 | Similar to above, but 1000 gammas are emmited isotropically from the |
---|
| 115 | source. No hits or tracks are draw, and the verbosity is reduced. |
---|
| 116 | |
---|
| 117 | alpha.mac |
---|
| 118 | Shoots one 5.486 MeV alpha particle upwards from the calibrator and |
---|
| 119 | traces the scintillation light to the PMT. All tracks except |
---|
| 120 | optical photons are drawn. Event summary is written to file |
---|
| 121 | "hits.out" and PMT hits to file "pmt.out". |
---|
| 122 | |
---|
| 123 | alpha_1000.mac |
---|
| 124 | Similar to above, but 1000 alphas are emmited isotropically from the |
---|
| 125 | source. No tracks are stored, and the verbosity is reduced. Event |
---|
| 126 | summary is written to file "hits.out", PMT hits are not written out. |
---|
| 127 | |
---|
| 128 | neutron.mac |
---|
| 129 | Shoots one 2.48 MeV neutron inside the room aimed at the detector. |
---|
| 130 | All tracks except scintillation photons are drawn (custom colours). |
---|
| 131 | Gammas are not killed at the concrete wall. Event summery is written |
---|
| 132 | to file "hits.out". |
---|
| 133 | |
---|
| 134 | ambe_spectrum.mac |
---|
| 135 | Produces a spectrum of neutrons according to an approximation of a Am/Be |
---|
| 136 | neutron source. However, this uses the GPS and therefore will only work |
---|
| 137 | after compilation with the DMXENV_GPS_USE environment variable (see below) |
---|
| 138 | |
---|
| 139 | |
---|
| 140 | sourceAm241.mac |
---|
| 141 | Forces the decay of 241Am nuclei in the calibrator and tracks the |
---|
| 142 | resulting particles (237Np + alpha + gamma from 237Np |
---|
| 143 | de-excitation). All tracks are drawn except scintillation photons. |
---|
| 144 | Event summery is written to file "hits.out". |
---|
| 145 | -> To be added for a future release |
---|
| 146 | |
---|
| 147 | |
---|
| 148 | sourceAm241_1000.mac |
---|
| 149 | Similar to above but for 1000 events. No tracks are stored. |
---|
| 150 | -> To be added for a future release |
---|
| 151 | |
---|
| 152 | |
---|
| 153 | Note: |
---|
| 154 | The following environment variables need to be set: |
---|
| 155 | G4RADIOACTIVEDATA : points to Radioactive Decay Data files |
---|
| 156 | G4LEDATA : points to low energy data base |
---|
| 157 | G4LEVELGAMMADATA : points to PhotoEvaporation data |
---|
| 158 | NeutronHPCrossSections : points to neutron data files |
---|
| 159 | |
---|
| 160 | In addition if you require to use the full General Particle Source then the |
---|
| 161 | variable DMXENV_GPS_USE can be set. The DMX gun is still included in order to |
---|
| 162 | allow forward compatibility should the GPS change. |
---|
| 163 | |
---|
| 164 | |
---|
| 165 | ANALYSIS: |
---|
| 166 | |
---|
| 167 | In order to use the AIDA 3.2.1 compliant analysis set the environment variable |
---|
| 168 | G4ANALYSIS_USE before building. If a previous build exists a gmake clean has |
---|
| 169 | to be executed. |
---|
| 170 | |
---|
| 171 | setenv G4ANALYSIS_USE 1 |
---|
| 172 | |
---|
| 173 | then gmake |
---|
| 174 | |
---|
| 175 | The program should then produce hbook histogram files |
---|
| 176 | |
---|
| 177 | Using AIDA 3.2.1 interfaces three different analysis packages can be utilised: |
---|
| 178 | Java Analysis Studio (JAS) |
---|
| 179 | OpenScientist |
---|
| 180 | PI |
---|
| 181 | |
---|
| 182 | Please for more information look at http://aida.freehep.org |
---|
| 183 | |
---|
| 184 | SEEDS: |
---|
| 185 | |
---|
| 186 | The seeds of event hits are stored in the hit record file. These can be used |
---|
| 187 | to repeat events for visualisation, test crashes/idiosyncracies: |
---|
| 188 | |
---|
| 189 | /random/setDirectoryName ./seeds |
---|
| 190 | /random/resetEngineFrom currentEvent.rndm |
---|
| 191 | /random/saveThisEvent |
---|
| 192 | /random/setSavingFlag |
---|
| 193 | |
---|
| 194 | The file currentEvent.rndm should contain the two seeds which were outputed in |
---|
| 195 | the scintHit file. |
---|
| 196 | |
---|
| 197 | ERRORS: |
---|
| 198 | |
---|
| 199 | When running interactively the following error will be shown: |
---|
| 200 | |
---|
| 201 | BooleanProcessor::caseIE : unimplemented case |
---|
| 202 | BooleanProcessor::caseIE : unimplemented case |
---|
| 203 | BooleanProcessor::caseIE : unimplemented case |
---|
| 204 | BooleanProcessor::caseIE : unimplemented case |
---|
| 205 | BooleanProcessor::caseIE : unimplemented case |
---|
| 206 | BooleanProcessor::caseIE : unimplemented case |
---|
| 207 | BooleanProcessor: boolean operation failed |
---|
| 208 | BooleanProcessor::caseIE : unimplemented case |
---|
| 209 | BooleanProcessor::caseIE : unimplemented case |
---|
| 210 | BooleanProcessor::caseIE : unimplemented case |
---|
| 211 | BooleanProcessor::caseIE : unimplemented case |
---|
| 212 | BooleanProcessor: boolean operation failed |
---|
| 213 | |
---|
| 214 | This is a "feature" of the visualisation of boolean volumes, but does not |
---|
| 215 | affect functionality/performance so can be ignored. |
---|
| 216 | |
---|
| 217 | NB: |
---|
| 218 | If using explicit libraries (?) i.e. non-shared then compilation time with |
---|
| 219 | neutrons in physics list is very long (>5 minutes) - check this............. |
---|
| 220 | It is more efficient to use shared libraries that are loaded at run-time with |
---|
| 221 | increased initialisation time (at run-time). |
---|
| 222 | |
---|
| 223 | Also if using shared libraries the load time at run-time may be several minutes |
---|
| 224 | - this is partially due to the neutron implementation requires full data sets |
---|
| 225 | for each isotope being specified. |
---|
| 226 | |
---|
| 227 | |
---|
| 228 | -------------------------------------- |
---|
| 229 | |
---|
| 230 | If running on Redhat 7.0 or above set G4SYSTEM to Linux-g++, alternatively you |
---|
| 231 | can install backward compatibility to egcs, however, requires |
---|
| 232 | config/sys/Linux-egcs.gmk to be altered so that CXX is set to kgcc |
---|
| 233 | (compared to g++ in original file) |
---|
| 234 | |
---|
| 235 | -------------------------------------- |
---|
| 236 | |
---|
| 237 | |
---|
| 238 | Alex Howard, 29/11/01 |
---|
| 239 | updated 18/06/02 |
---|
| 240 | |
---|