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