Jpp 20.0.0-195-g190c9e876
the software that should make you happy
Loading...
Searching...
No Matches
JORCAReconstruction.sh
Go to the documentation of this file.
1#!/usr/bin/env zsh
2#
3# \author mdejong
4#
5script=${0##*/}
6
7# ------------------------------------------------------------------------------------------
8#
9# Script to run the JPP Reconstruction Chain for ORCA.
10#
11# ------------------------------------------------------------------------------------------
12
13if [ -z $JPP_DIR ]; then
14 echo "Variable JPP_DIR undefined."
15 exit
16fi
17
18source $JPP_DIR/setenv.sh $JPP_DIR
19
20
21source reconstruction.sh
22
23set_variable: DEBUG RECONSTRUCTION_DEBUG 1
24set_variable: OPTION RECONSTRUCTION_OPTION --!
25set_variable: PDF RECONSTRUCTION_PDF $JPP_DATA/J%p.dat.gz
26set_variable: NUMBER_OF_EVENTS RECONSTRUCTION_NUMBER_OF_EVENTS $((2**31))
27set_variable: NUMBER_OF_THREADS RECONSTRUCTION_NUMBER_OF_THREADS 0
28set_variable+ OVERWRITE RECONSTRUCTION_OVERWRITE 1
29set_variable: TMAX_S RECONSTRUCTION_TMAX_S 100
30set_variable: RATE_HZ RECONSTRUCTION_RATE_HZ 7.0e3
31
32if do_usage $*; then
33 usage "$script <detector file> <input file> <output file> <PMT parameters file> [(calibration file)+]"
34fi
35
36if (( $# < 4 )); then
37 fatal "Wrong number of arguments."
38fi
39
40set_variable DETECTOR $argv[1]
41set_variable INPUT_FILE $argv[2]
42set_variable OUTPUT_FILE $argv[3]
43set_variable PMT $argv[4]
44set_array CALIBRATION $argv[5,-1]
45
46typeset -A PARAMETERS
47
48PARAMETERS+=(
49 muon.prefit.sigma_ns 5.0
50 muon.prefit.gridAngle_deg 4.0
51 muon.prefit.useL0 1
52 muon.prefit.numberOfOutliers 2
53 muon.prefit.numberOfPrefits 60
54 muon.prefit.numberOfPostfits 20
55 muon.prefit.TMaxLocal_ns 12.0
56 muon.prefit.ctMin 0.0
57 muon.prefit.roadWidth_m 50.0
58 muon.prefit.Qwatershed 1.5)
59
60PARAMETERS+=(
61 muon.simplex.sigma_ns 3.0
62 muon.simplex.TMaxLocal_ns 15.0
63 muon.simplex.ctMin 0.0
64 muon.simplex.roadWidth_m 50.0
65 muon.simplex.numberOfPrefits 0)
66
67PARAMETERS+=(
68 muon.start.roadWidth_m 50.0
69 muon.start.R_Hz $RATE_HZ
70 muon.start.TMin_ns -25.0
71 muon.start.TMax_ns +25.0
72 muon.start.Pmin 1.0e-8
73 muon.start.Pmin1 1.0e-3
74 muon.start.Pmin2 1.0e-2
75 muon.start.numberOfPrefits 0)
76
77PARAMETERS+=(
78 muon.gandalf.TTS_ns 1.0
79 muon.gandalf.E_GeV 10.0
80 muon.gandalf.TMin_ns -50.0
81 muon.gandalf.TMax_ns +450.0
82 muon.gandalf.ZMin_m -3.0
83 muon.gandalf.ZMax_m +9.0
84 muon.gandalf.R_Hz $RATE_HZ
85 muon.gandalf.roadWidth_m 50.0
86 muon.gandalf.cosLR 0.7
87 muon.gandalf.numberOfPrefits 0)
88
89PARAMETERS+=(
90 muon.start2.roadWidth_m 50.0
91 muon.start2.R_Hz $RATE_HZ
92 muon.start2.TMin_ns -25.0
93 muon.start2.TMax_ns +25.0
94 muon.start2.Pmin1 1.0e-3
95 muon.start2.Pmin2 1.0e-2
96 muon.start2.numberOfPrefits 1
97 muon.start2.numberOfPostfits 0)
98
99PARAMETERS+=(
100 muon.energy.roadWidth_m 80.0
101 muon.energy.R_Hz $RATE_HZ
102 muon.energy.EMin_log 0.0
103 muon.energy.EMax_log 7.0
104 muon.energy.TMin_ns -25.0
105 muon.energy.TMax_ns +25.0
106 muon.energy.ZMin_m -7.5
107 muon.energy.resolution 0.01
108 muon.energy.mestimator 3
109 muon.energy.numberOfPrefits 0)
110
111PARAMETERS+=(
112 muon.features.TTS_ns 1.0
113 muon.features.E_GeV 10.0
114 muon.features.TMin_ns -50.0
115 muon.features.TMax_ns +450.0
116 muon.features.ZMin_m -3.0
117 muon.features.ZMax_m +9.0
118 muon.features.R_Hz $RATE_HZ
119 muon.features.roadWidth_m 50.0
120 muon.features.numberOfPrefits 0)
121
122PARAMETERS+=(
123 shower.prefit.numberOfPrefits 200
124 shower.prefit.numberOfGrids 2
125 shower.prefit.factoryLimit 20
126 shower.prefit.sigma_ns 2.5
127 shower.prefit.numberOfOutliers 4
128 shower.prefit.TMaxLocal_ns 10.0
129 shower.prefit.TMaxExtra_ns 80.0
130 shower.prefit.ctMin -0.3
131 shower.prefit.DMax_m 50.0
132 shower.prefit.pos_grid_m 10.0
133 shower.prefit.pos_step_m 10.0
134 shower.prefit.time_grid_ns 50.0
135 shower.prefit.time_step_ns 50.0)
136
137PARAMETERS+=(
138 shower.simplex.numberOfPrefits 0
139 shower.simplex.sigma_ns 2.5
140 shower.simplex.TMaxLocal_ns 15.0
141 shower.simplex.TMax_ns +60.0
142 shower.simplex.TMin_ns -60.0
143 shower.simplex.DMax_m 50.0
144 shower.simplex.ctMin 0.0
145 shower.simplex.mestimator 1)
146
147PARAMETERS+=(
148 shower.position.numberOfPrefits 100
149 shower.position.TMax_ns +120.0
150 shower.position.TMin_ns -120.0
151 shower.position.DMax_m 80.0
152 shower.position.Emin_GeV 1.0
153 shower.position.Emax_GeV 200.0
154 shower.position.En 3.0
155 shower.position.R_Hz $RATE_HZ
156 shower.position.TTS_ns 2.0)
157
158PARAMETERS+=(
159 shower.direction.numberOfPrefits 1
160 shower.direction.TMax_ns +30.0
161 shower.direction.TMin_ns -30.0
162 shower.direction.DMax_m 60.0
163 shower.direction.R_Hz $RATE_HZ
164 shower.direction.scanAngle_deg 10.0
165 shower.direction.Emin_GeV 1
166 shower.direction.Emax_GeV 1000
167 shower.direction.En 100)
168
169PARAMETERS+=(
170 shower.fit.numberOfPrefits 4
171 shower.fit.TMax_ns +30.0
172 shower.fit.TMin_ns -30.0
173 shower.fit.DMax_m 60.0
174 shower.fit.DStep_m 5.0
175 shower.fit.R_Hz $RATE_HZ
176 shower.fit.mestimator 3)
177
178eval JORCAReconstruction \
179 -a $DETECTOR \
180 -+${^CALIBRATION} \
181 -T $TMAX_S \
182 -f ${INPUT_FILE} \
183 -o ${OUTPUT_FILE} \
184 -F $PDF \
185 -n ${NUMBER_OF_EVENTS} \
186 -N ${NUMBER_OF_THREADS} \
187 -@ `make_equation PARAMETERS` \
188 -P $PMT \
189 -d $DEBUG \
190 $OPTION
191
192CHECK_EXIT_CODE
Generated by doxygen 1.12.0