forked from WCSim/WCSim
-
Notifications
You must be signed in to change notification settings - Fork 0
/
vis.mac
167 lines (134 loc) · 6.51 KB
/
vis.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
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
# Visualization and setup macro
# This file is run by default if you just type WCSim and nothing else.
/run/verbose 0
/tracking/verbose 0
/hits/verbose 0
## select the geometry
# Default config if you do nothing is currently SuperK
#
# The tube size is fixed for SK to 20"
# These are fixed geometries for validation
#/WCSim/WCgeom SuperK
# Some other SuperK options:
#/WCSim/WCgeom SuperK_20inchPMT_20perCent # Note: the actual coverage is 20.27%
#/WCSim/WCgeom SuperK_20inchBandL_20perCent # Note: the actual coverage is 20.27%
#/WCSim/WCgeom SuperK_12inchBandL_15perCent # Note: the actual coverage is 14.59%
#/WCSim/WCgeom SuperK_20inchBandL_14perCent # Note: the actual coverage is 13.51%
# Cylindrical detector with a height of 60 m and a diameter of 74 m
#/WCSim/WCgeom Cylinder_60x74_20inchBandL_14perCent # Note: the actual coverage is 13.51%
#/WCSim/WCgeom Cylinder_60x74_20inchBandL_40perCent
# Generic cylindrical detector with a height of 100m and a
# diameter of 69m with 12" HPD and 14.59% photocoverage
#/WCSim/WCgeom Cylinder_12inchHPD_15perCent # Note: the actual coverage is 14.59%
#/WCSim/WCgeom HyperK #default length is 49500 mm unless changed below.
#Added for detector length option (ONLY FOR HYPER-K GEOMETRY!!) 2014/07/29. This simulates a segment of the Hyper-K detector of a given length.
#/WCSim/HyperK/waterTank_Length 49500. mm # Equivalent length for 5 partitions(ie. will simulate a detector 1/10 of the HK length)
#/WCSim/HyperK/waterTank_Length 61875. mm # Equivalent length for 4 partitions
#/WCSim/HyperK/waterTank_Length 82500. mm # Equivalent length for 3 partitions
#/WCSim/HyperK/waterTank_Length 123750. mm # Equivalent length for 2 partitions
#/WCSim/HyperK/waterTank_Length 247500. mm # Equivalent length for 1 partition
#Added for the PMT QE option 08/17/10 (XQ)
# 1. Stacking only mean when the photon is generated
# the QE is applied to reduce the total number of photons
# 2. Stacking and sensitivity detector
# In the stacking part, the maximum QE is applied to reduce
# the total number of photons
# On the detector side, the rest of QE are applied according to QE/QE_max
# distribution. This option is in particular important for the WLS
# 3. The last option means all the QE are applied at the detector
# Good for the low energy running.
/WCSim/PMTQEMethod Stacking_Only
#/WCSim/PMTQEMethod Stacking_And_SensitiveDetector
#/WCSim/PMTQEMethod SensitiveDetector_Only
# turn on or off the collection efficiency
/WCSim/PMTCollEff on
#/WCSim/Construct # This must be uncommented to construct a new detector configuration.
# command to choose save or not save the pi0 info 07/03/10 (XQ)
#/WCSim/SavePi0 true
/WCSim/SavePi0 false
#choose the Trigger & Digitizer type (and options)
/DAQ/Digitizer SKI
/DAQ/Trigger NDigits
#generic digitizer options (defaults are class-specific. Can be overridden here)
# how long is the digitizer dead before a new hit on the same PMT can be recorded (ns)?
#/DAQ/DigitizerOpt/DeadTime 0
# how long does the digitizer integrate for (ns)?
#/DAQ/DigitizerOpt/IntegrationWindow 200
#generic trigger options (defaults are class-specfic. Can be overrideen here)
# allow the number of digits per PMT per trigger to be > 1?
#/DAQ/MultiDigitsPerTrigger false
#trigger algorithm development options (defaults are set in WCSimWCDAQMessenger. Can be overidden here)
# do we want to save only triggered events (mode 0), both triggered events & failed events (mode 1), or only failed events (mode 2)?
#/DAQ/TriggerSaveFailures/Mode 0
# what trigger time to use for the failed events (ns)?
#/DAQ/TriggerSaveFailures/TriggerTime 100
# the digits in which range around the trigger time (ns) should be saved with the event
#/DAQ/TriggerSaveFailures/PreTriggerWindow -400
#/DAQ/TriggerSaveFailures/PostTriggerWindow +950
#options for NDigits-like triggers (defaults are class-specific. Can be overridden here)
# control the NDigits trigger threshold
#/DAQ/TriggerNDigits/Threshold 25
# the trigger counts digits (looking for threshold) in this window (ns)
#/DAQ/TriggerNDigits/Window 200
# automatically adjust the NDigits threshold depending on the average noise occupancy during the trigger window
#/DAQ/TriggerNDigits/AdjustForNoise true
# the digits in which range around the trigger time (ns) should be saved with the event
#/DAQ/TriggerNDigits/PreTriggerWindow -400
#/DAQ/TriggerNDigits/PostTriggerWindow +950
# default dark noise frequency is 0 kHz
/DarkRate/SetDarkRate 4.2 kHz #This is the value for SKI set in SKDETSIM.
#/DarkRate/SetDarkRate 8.4 kHz #For 20 inch HPDs and Box and Line PMTs, based on High QE 20" R3600 dark rate from EGADS nov 2014
#/DarkRate/SetDarkRate 3.0 kHz #For 12 inch HPDs and Box and Line PMTs, based on High QE 20" R3600 dark rate from EGADS nov 2014
# command to multiply the dark rate.
# Convert dark noise frequency before digitization to after digitization by setting suitable factor
/DarkRate/SetConvert 1.367 #For Normal PMT
#/DarkRate/SetConvert 1.119 #For HPDs
#/DarkRate/SetConvert 1.126 #For Box and Line PMTs
# Select which time window(s) to add dark noise to
#/DarkRate/SetDarkMode 0 to add dark noise to a time window starting at
#/DarkRate/SetDarkLow to /DarkRate/SetDarkHigh [time in ns]
#/DarkRate/SetDarkMode 1 adds dark noise hits to a window of
#width /DarkRate/SetDarkWindow [time in ns] around each hit
#i.e. hit time ± (/DarkRate/SetDarkWindow) / 2
/DarkRate/SetDarkMode 1
/DarkRate/SetDarkHigh 100000
/DarkRate/SetDarkLow 0
/DarkRate/SetDarkWindow 4000
/vis/scene/create
/vis/open OGLSX
# You can also set the size and position of the window if you like
#/vis/open OGLSX 1000x1000-0+0
/vis/ogl/set/printSize 1000 1000
# Vis Settings for SK
/vis/viewer/zoom 1.2
#/vis/viewer/set/viewpointThetaPhi 45 45 deg
/vis/viewer/set/viewpointThetaPhi 110 0 deg
/vis/viewer/set/upVector 0 0 1
# Vis settings for HK
#/vis/viewer/set/upVector 0 1 0
#/vis/viewer/set/viewpointThetaPhi 50 165 deg
#/vis/viewer/zoom 1.6
#/vis/ogl/set/printMode pixmap
#/vis/ogl/printEPS
#/vis/open RayTracer
#/vis/viewer/set/upVector 0 0 1
#/vis/rayTracer/eyePosition 70 0 70
#/vis/rayTracer/trace
/vis/scene/add/hits
/vis/scene/add/trajectories
/vis/scene/endOfEventAction accumulate
## select the input nuance-formatted vector file
## you can of course use your own
#/mygen/vecfile inputvectorfile
#/mygen/vecfile h2o.2km.001-009x3_G4.kin
#/mygen/vecfile mu+.out
# Or you can use the G4 Particle Gun
/mygen/generator normal
/gun/particle mu+
/gun/energy 1000 MeV
/gun/direction 0 1 0
/gun/position 0 0 0
## change the name of the output root file, default = wcsim.root
/WCSimIO/RootFile wcsim.root
/run/beamOn 100
#exit