Jpp 21.0.0-rc.1
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.TMin_ns -25.0
70 muon.start.TMax_ns +25.0
71 muon.start.Pmin 1.0e-8
72 muon.start.Pmin1 1.0e-3
73 muon.start.Pmin2 1.0e-2
74 muon.start.numberOfPrefits 0)
75
76PARAMETERS+=(
77 muon.gandalf.TTS_ns 1.0
78 muon.gandalf.E_GeV 10.0
79 muon.gandalf.TMin_ns -50.0
80 muon.gandalf.TMax_ns +450.0
81 muon.gandalf.ZMin_m -3.0
82 muon.gandalf.ZMax_m +9.0
83 muon.gandalf.R_Hz $RATE_HZ
84 muon.gandalf.roadWidth_m 50.0
85 muon.gandalf.cosLR 0.7
86 muon.gandalf.numberOfPrefits 0)
87
88PARAMETERS+=(
89 muon.start2.roadWidth_m 50.0
90 muon.start2.TMin_ns -25.0
91 muon.start2.TMax_ns +25.0
92 muon.start2.Pmin1 1.0e-3
93 muon.start2.Pmin2 1.0e-2
94 muon.start2.numberOfPrefits 1
95 muon.start2.numberOfPostfits 0)
96
97PARAMETERS+=(
98 muon.energy.roadWidth_m 80.0
99 muon.energy.R_Hz $RATE_HZ
100 muon.energy.EMin_log 0.0
101 muon.energy.EMax_log 7.0
102 muon.energy.TMin_ns -25.0
103 muon.energy.TMax_ns +25.0
104 muon.energy.ZMin_m -7.5
105 muon.energy.resolution 0.01
106 muon.energy.mestimator 3
107 muon.energy.numberOfPrefits 0)
108
109PARAMETERS+=(
110 muon.features.TTS_ns 1.0
111 muon.features.E_GeV 10.0
112 muon.features.TMin_ns -50.0
113 muon.features.TMax_ns +450.0
114 muon.features.ZMin_m -3.0
115 muon.features.ZMax_m +9.0
116 muon.features.R_Hz $RATE_HZ
117 muon.features.roadWidth_m 50.0
118 muon.features.numberOfPrefits 0)
119
120PARAMETERS+=(
121 shower.prefit.numberOfPrefits 200
122 shower.prefit.numberOfGrids 2
123 shower.prefit.factoryLimit 20
124 shower.prefit.sigma_ns 2.5
125 shower.prefit.numberOfOutliers 4
126 shower.prefit.TMaxLocal_ns 10.0
127 shower.prefit.TMaxExtra_ns 80.0
128 shower.prefit.ctMin -0.3
129 shower.prefit.DMax_m 50.0
130 shower.prefit.pos_grid_m 10.0
131 shower.prefit.pos_step_m 10.0
132 shower.prefit.time_grid_ns 50.0
133 shower.prefit.time_step_ns 50.0)
134
135PARAMETERS+=(
136 shower.simplex.numberOfPrefits 0
137 shower.simplex.sigma_ns 2.5
138 shower.simplex.TMaxLocal_ns 15.0
139 shower.simplex.TMax_ns +60.0
140 shower.simplex.TMin_ns -60.0
141 shower.simplex.DMax_m 50.0
142 shower.simplex.ctMin 0.0
143 shower.simplex.mestimator 1)
144
145PARAMETERS+=(
146 shower.position.numberOfPrefits 100
147 shower.position.TMax_ns +120.0
148 shower.position.TMin_ns -120.0
149 shower.position.DMax_m 80.0
150 shower.position.Emin_GeV 1.0
151 shower.position.Emax_GeV 200.0
152 shower.position.En 3.0
153 shower.position.R_Hz $RATE_HZ
154 shower.position.TTS_ns 2.0)
155
156PARAMETERS+=(
157 shower.direction.numberOfPrefits 1
158 shower.direction.TMax_ns +30.0
159 shower.direction.TMin_ns -30.0
160 shower.direction.DMax_m 60.0
161 shower.direction.R_Hz $RATE_HZ
162 shower.direction.scanAngle_deg 10.0
163 shower.direction.Emin_GeV 1
164 shower.direction.Emax_GeV 1000
165 shower.direction.En 100)
166
167PARAMETERS+=(
168 shower.fit.numberOfPrefits 4
169 shower.fit.TMax_ns +30.0
170 shower.fit.TMin_ns -30.0
171 shower.fit.DMax_m 60.0
172 shower.fit.DStep_m 5.0
173 shower.fit.R_Hz $RATE_HZ
174 shower.fit.mestimator 3)
175
176eval JORCAReconstruction \
177 -a $DETECTOR \
178 -+${^CALIBRATION} \
179 -T $TMAX_S \
180 -f ${INPUT_FILE} \
181 -o ${OUTPUT_FILE} \
182 -F $PDF \
183 -n ${NUMBER_OF_EVENTS} \
184 -N ${NUMBER_OF_THREADS} \
185 -@ `make_equation PARAMETERS` \
186 -P $PMT \
187 -d $DEBUG \
188 $OPTION
189
190CHECK_EXIT_CODE
Generated by doxygen 1.12.0