forked from legend-exp/g4simple
-
Notifications
You must be signed in to change notification settings - Fork 0
/
Copy pathrun.mac
119 lines (94 loc) · 3.37 KB
/
run.mac
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
/run/verbose 1
/event/verbose 0
/tracking/verbose 0
/g4simple/setRandomSeed true
#/random/setSeed 12345678
#Need to set before Phyics List
/optics/fiberDetProb 0.4
# Need to set the physics list before we can do some of the other commands.
/g4simple/setReferencePhysList Shielding
# Set GDML file name
# The bool after the file name turns validation on / off
/g4simple/setDetectorGDML L200 false
# Set up output. Choose a format:
#/g4simple/setOutputFormat csv
#/g4simple/setOutputFormat xml
#/g4simple/setOutputFormat root
#/g4simple/setOutputFormat hdf5
# Uncomment to override an output's standard option
#/g4simple/setOutputOption stepwise
#/g4simple/setOutputOption eventwise
# Change the name of the output file (by default it's g4simple.[ext])
#/analysis/setFileName geCounterOut
# Steps are only written out if they are in a volume that matches a pattern set
# via the setVolID command
# Example setting a full volume name (exact match):
/g4simple/setVolID innerShroud 1
/g4simple/setVolID outerShroud 1 #dont distinguish for now
#the following volumes should NOT show up as detectors any more
/g4simple/setVolID wslrCopper 3 #TODO keep typo?
/g4simple/setVolID larVolume 4
# Example using a regular expression to match multiple volumes, extract an
# integer from the name, and set it as the volID
#/g4simple/setVolID .*Detector([0-9]*).* $1
# To get all steps in all volumes, comment out all setVolID commands (removes
# all volume IDs from output), or uncomment the following line to keep the
# volume IDs:
#/g4simple/recordAllSteps
/g4simple/verbose 0
/generator/verbose 0
#set geometry
/geometry/innerShroud/innerRadius 122.5 mm
/geometry/innerShroud/outerRadius 134.5 mm
/geometry/innerShroud/height 1300 mm
/geometry/outerShroud/innerRadius 283.0 mm
/geometry/outerShroud/outerRadius 295.0 mm
/geometry/outerShroud/height 1500 mm
/geometry/wlsr/radius 700 mm
/geometry/wlsr/height 3500 mm
/geometry/wlsr/tpbThickness 0.001 mm
/geometry/wlsr/cuThickness 0.03 mm
/geometry/wlsr/tetraTexThickness 0.01 mm
/geometry/cryostat/wallThickness 10 mm
/geometry/Ge/discHeight 100 mm
/geometry/Ge/discRadius 40 mm
/geometry/Ge/discGap 50 mm
/geometry/Ge/arrayRadius 200 mm
/geometry/Ge/discPerString 12
/geometry/Ge/nrStrings 14
#set optics
/optics/fiberAbsProb .2
/optics/lArAbsLength 100 cm
/optics/visAbsLength 10 m
/optics/lArScintWL 128 nm
/optics/tpbScintWL 450 nm
/optics/lArRayToggle false
/optics/setBlackWLSR false
#update geometry & optics
/update
/run/initialize
# If you want to see the list of available NIST materials (e.g. to help you
# build your gdml file) uncomment this line
#/material/nist/listMaterials
# If you want to see the list of all physical volume names (e.g. to set volume
# IDs above) uncomment this line
#/g4simple/listPhysVols
# Set up gps after /run/initialize
# This is for 208Tl events distributed throughout the "source" volume
#/gps/List
#/gps/particle opticalphoton
#/gps/energy 9.68626 eV
#/gps/direction 1 0 0
#/gps/pos/centre 0 0 0
#rad min: only 0 taken by program, even if different in macro
/generator/SetRadiusMax 700 mm
/generator/SetRadiusMin 0 mm
/generator/SetHeight 0 mm
/generator/SetBinWidth 10 mm
/generator/SetNParticles 100
/generator/SetCenterVector 0.0 0.0 0.0 mm
/generator/Set1D false
#should a voxel be aborted (and counted as zero) upon a single non-lar primary?
/generator/abortOnNonlar true
/write/filename temp.root
#/run/beamOn 10 NO SUCH STUFF HERE XXX