Jpp  17.1.1
the software that should make you happy
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Macros Pages
JORCAShowerReconstruction.sh
Go to the documentation of this file.
1 #!/bin/zsh
2 #
3 # \author adomi
4 #
5 script=${0##*/}
6 
7 
8 # ------------------------------------------------------------------------------------------
9 #
10 # Script to run the JPP Shower Reconstruction Chain for ORCA.
11 #
12 # ------------------------------------------------------------------------------------------
13 
14 
15 if [ -z $JPP_DIR ]; then
16  echo "Variable JPP_DIR undefined."
17  exit
18 fi
19 
20 source $JPP_DIR/setenv.sh $JPP_DIR
21 source $JPP_DIR/externals/km3net-dataformat/definitions/reconstruction.sh
22 
23 set_variable: DEBUG RECONSTRUCTION_DEBUG 1
24 set_variable: PDF RECONSTRUCTION_PDF $JPP_DATA/J%p.dat
25 set_variable: SUFFIX RECONSTRUCTION_SUFFIX jshower
26 set_variable: ENERGY_CORRECTION RECONSTRUCTION_ENERGY_CORRECTION energy_correction_showerfit_fullorca.txt
27 set_variable: NUMBER_OF_EVENTS RECONSTRUCTION_NUMBER_OF_EVENTS $((2**31))
28 set_variable+ OVERWRITE RECONSTRUCTION_OVERWRITE 1
29 set_variable: TMAX_S RECONSTRUCTION_TMAX_S 100
30 
31 if do_usage $*; then
32  usage "$script <detector file> <input file> <output file> <PMT parameters file> [(calibration file)+]"
33 fi
34 
35 set_variable DETECTOR $argv[1]
36 set_variable INPUT_FILE $argv[2]
37 set_variable OUTPUT_FILE $argv[3]
38 set_variable PMT $argv[4]
39 set_array CALIBRATION $argv[5,-1]
40 
41 if [[ ! "$OUTPUT_FILE" == *"%"* ]]; then
42  fatal "The output file must have the wildcard % in the name, e.g. %.root"
43 fi
44 
45 typeset -A PARAMETERS
46 typeset -a APPLICATIONS
47 
48 APPLICATIONS=()
49 
50 
51 if [[ "$OVERWRITE" == 1 ]]; then
52  rm -vf \
53  ${OUTPUT_FILE/\%/jshprefit} \
54  ${OUTPUT_FILE/\%/jshpointsimplex} \
55  ${OUTPUT_FILE/\%/jshpositionfit} \
56  ${OUTPUT_FILE/\%/jshenergyprefit} \
57  ${OUTPUT_FILE/\%/jshdirectionprefit} \
58  ${OUTPUT_FILE/\%/jshfit} \
59  ${OUTPUT_FILE/\%/jhistory} \
60  ${OUTPUT_FILE/\%/${SUFFIX}} \
61  ${OUTPUT_FILE/\%/${SUFFIX}.aanet}
62 fi
63 
64 
65 # ------------------------------------------------------------------------------------------
66 # JShowerPrefit
67 # ------------------------------------------------------------------------------------------
68 
69 PARAMETERS=(
70  factoryLimit 40
71  NMaxHits 20
72  sigma_ns 4
73  numberOfOutliers 4
74  TMaxLocal_ns 120
75  ctMin -0.7
76  roadWidth_m 50
77  numberOfPrefits 50)
78 
79 timer_start
80 
81 eval JShowerPrefit \
82  -a $DETECTOR \
83  -+${^CALIBRATION} \
84  -T $TMAX_S \
85  -f $INPUT_FILE \
86  -o ${OUTPUT_FILE/\%/jshprefit} \
87  -n ${NUMBER_OF_EVENTS} \
88  -@ `make_equation PARAMETERS` \
89  -d ${DEBUG} --!
90 
91 CHECK_EXIT_CODE
92 
93 APPLICATIONS+=$JSHOWERPREFIT
94 
95 timer_stop
96 timer_print
97 
98 
99 # ------------------------------------------------------------------------------------------
100 # JShowerPointSimplex
101 # ------------------------------------------------------------------------------------------
102 
103 PARAMETERS=(
104  numberOfPrefits 0
105  numberOfOutfits 100
106  minPrefitsSize 9
107  sigma_ns 1
108  TMaxLocal_ns 15
109  TMax_ns 60
110  TMin_ns -60
111  roadWidth_m 60
112  ctMin 0.2
113  TWindow_ns "-130 120"
114  pos_grid_m " -20 20"
115  pos_step_m 20
116  time_grid_ns " -90 90"
117  time_step_ns 30
118  simplex_step_m 2
119  simplex_step_ns 10
120  tukey_std_dev 20.0)
121 
122 timer_start
123 
124 eval JShowerPointSimplex \
125  -a $DETECTOR \
126  -+${^CALIBRATION} \
127  -T $TMAX_S \
128  -f ${OUTPUT_FILE/\%/jshprefit} \
129  -o ${OUTPUT_FILE/\%/jshpointsimplex} \
130  -@ `make_equation PARAMETERS` \
131  -d ${DEBUG} --!
132 
133 CHECK_EXIT_CODE
134 
135 rm -vf ${OUTPUT_FILE/\%/jshprefit}
136 
137 APPLICATIONS+=$JSHOWERPOINTSIMPLEX
138 
139 timer_stop
140 timer_print
141 
142 
143 # ------------------------------------------------------------------------------------------
144 # JShowerPositionFit
145 # ------------------------------------------------------------------------------------------
146 
147 PARAMETERS=(
148  numberOfPrefits 0
149  numberOfOutfits 5
150  TMax_ns 120
151  TMin_ns -120
152  roadWidth_m 80
153  E_GeV 10
154  R_Hz 10000
155  TTS_ns 2)
156 
157 timer_start
158 
159 eval JShowerPositionFit \
160  -a $DETECTOR \
161  -+${^CALIBRATION} \
162  -T $TMAX_S \
163  -f ${OUTPUT_FILE/\%/jshpointsimplex} \
164  -F $PDF \
165  -o ${OUTPUT_FILE/\%/jshpositionfit} \
166  -@ `make_equation PARAMETERS` \
167  -d ${DEBUG} --!
168 
169 CHECK_EXIT_CODE
170 
171 rm -vf ${OUTPUT_FILE/\%/jshpointsimplex}
172 
173 APPLICATIONS+=$JSHOWERPOSITIONFIT
174 
175 timer_stop
176 timer_print
177 
178 
179 # ------------------------------------------------------------------------------------------
180 # JShowerEnergyPrefit
181 # ------------------------------------------------------------------------------------------
182 
183 PARAMETERS=(
184  numberOfPrefits 0
185  numberOfOutfits 1
186  Emin_GeV 1
187  Emax_GeV 100
188  resolution 0.01
189  TMax_ns 500
190  TMin_ns -500
191  R_Hz 10000
192  mestimator 3
193  roadWidth_m 80)
194 
195 timer_start
196 
197 eval JShowerEnergyPrefit \
198  -a $DETECTOR \
199  -+${^CALIBRATION} \
200  -T $TMAX_S \
201  -f ${OUTPUT_FILE/\%/jshpositionfit} \
202  -F $PDF \
203  -o ${OUTPUT_FILE/\%/jshenergyprefit} \
204  -@ `make_equation PARAMETERS` \
205  -d $DEBUG --!
206 
207 CHECK_EXIT_CODE
208 
209 rm -vf ${OUTPUT_FILE/\%/jshpositionfit}
210 
211 APPLICATIONS+=$JSHOWERENERGYPREFIT
212 
213 timer_stop
214 timer_print
215 
216 
217 # ------------------------------------------------------------------------------------------
218 # JShowerDirectionPrefit
219 # ------------------------------------------------------------------------------------------
220 
221 PARAMETERS=(
222  numberOfPrefits 0
223  numberOfOutfits 4
224  TMax_ns 30
225  TMin_ns -30
226  roadWidth_m 80
227  DMin_m 80
228  E_GeV 13
229  R_Hz 10000
230  mestimator 3
231  MaxAngle_deg_lowE 40
232  scanAngle_deg_lowE 10
233  MaxAngle_deg_highE 20
234  scanAngle_deg_highE 5)
235 
236 timer_start
237 
238 eval JShowerDirectionPrefit \
239  -a $DETECTOR \
240  -+${^CALIBRATION} \
241  -T $TMAX_S \
242  -f ${OUTPUT_FILE/\%/jshenergyprefit} \
243  -F $PDF \
244  -o ${OUTPUT_FILE/\%/jshdirectionprefit} \
245  -@ `make_equation PARAMETERS` \
246  -d $DEBUG --!
247 
248 CHECK_EXIT_CODE
249 
250 rm -vf ${OUTPUT_FILE/\%/jshenergyprefit}
251 
252 APPLICATIONS+=$JSHOWERDIRECTIONPREFIT
253 
254 timer_stop
255 timer_print
256 
257 
258 # ------------------------------------------------------------------------------------------
259 # JShowerFit
260 # ------------------------------------------------------------------------------------------
261 
262 PARAMETERS=(
263  numberOfPrefits 4
264  TMax_ns 30
265  TMin_ns -30
266  roadWidth_m 80
267  R_Hz 10000
268  mestimator 3
269  fit_step 0.05)
270 
271 timer_start
272 
273 eval JShowerFit \
274  -a $DETECTOR \
275  -+${^CALIBRATION} \
276  -T $TMAX_S \
277  -f ${OUTPUT_FILE/\%/jshdirectionprefit} \
278  -F $PDF \
279  -o ${OUTPUT_FILE/\%/jshfit} \
280  -@ `make_equation PARAMETERS` \
281  -E $ENERGY_CORRECTION \
282  -d $DEBUG --!
283 
284 CHECK_EXIT_CODE
285 
286 rm -vf ${OUTPUT_FILE/\%/jshdirectionprefit}
287 
288 APPLICATIONS+=$JSHOWERCOMPLETEFIT
289 
290 timer_stop
291 timer_print
292 
293 
294 # ------------------------------------------------------------------------------------------
295 # JHistory
296 # ------------------------------------------------------------------------------------------
297 
298 timer_start
299 
300 JHistory \
301  -f ${OUTPUT_FILE/\%/jshfit} \
302  -o ${OUTPUT_FILE/\%/jhistory} \
303  -A "$APPLICATIONS[*]" \
304  -d $DEBUG --!
305 
306 CHECK_EXIT_CODE
307 
308 rm -vf ${OUTPUT_FILE/\%/jshfit}
309 
310 timer_stop
311 timer_print
312 
313 
314 # ------------------------------------------------------------------------------------------
315 # rename outputfile
316 # ------------------------------------------------------------------------------------------
317 
318 timer_start
319 
320 mv -v \
321  ${OUTPUT_FILE/\%/jhistory} \
322  ${OUTPUT_FILE/\%/${SUFFIX}}
323 
324 timer_stop
325 timer_print
326 
327 
328 # ------------------------------------------------------------------------------------------
329 # JConvertEvt
330 # ------------------------------------------------------------------------------------------
331 
332 timer_start
333 
334 JConvertEvt \
335  -f ${OUTPUT_FILE/\%/${SUFFIX}} \
336  -o ${OUTPUT_FILE/\%/${SUFFIX}.aanet} \
337  -a $DETECTOR \
338  -+${^CALIBRATION} \
339  -T $TMAX_S \
340  -P $PMT \
341  -d $DEBUG --!
342 
343 timer_stop
344 timer_print
$argv
then usage $script< detector identifier >< startdate\"YYYY-MM-DDHH:MM:SS\"><finaldate\"YYYY-MM-DDHH:MM:SS\"><QA/QCfile> fi case set_variable QAQC_TXT $argv[4]
Definition: JDataMonitor.sh:24
PDF
then echo Variable JPP_DIR undefined exit fi source $JPP_DIR setenv sh $JPP_DIR set_variable DEBUG set_variable PDF $JPP_DATA J p dat set_variable HISTOGRAM if do_usage *then usage $script< detector file >< inputfile >< outputfile > fi case set_variable PDF
Definition: JARCACalibrateMuon.sh:23
NUMBER_OF_EVENTS
then usage $script< input_file >< detector_file >< output_file(.root)> fi set_variable NUMBER_OF_TIMESLICES set_variable NUMBER_OF_SUMMARYSLICES set_variable NUMBER_OF_EVENTS set_variable ALL_PLOTS case set_variable NUMBER_OF_EVENTS
Definition: JRunAnalyzer.sh:18
E
then usage E
Definition: JMuonPostfit.sh:35
JDATABASE::PBS::PMT
static const JPBS_t PMT(3, 4, 2, 3)
PBS of photo-multiplier tube (PMT)
km3net
rm f XXX txt echo nlist all archive files that do have specific variants getArchive sh $ARCHIVE test km3net(mjg_2|mjg_1)" getFile.sh $ARCHIVE/test 100 101 km3net "(mjg_2|mjg_1)" XXX
name
then echo Enter input within $TIMEOUT_S seconds echo n User name
Definition: JCookie.sh:42
f
o $QUALITY_ROOT d $DEBUG!JPlot1D f
Definition: JDataQuality.sh:66
parameters
*fatal Wrong number of arguments esac JCookie sh typeset Z DETECTOR typeset Z SOURCE_RUN typeset Z TARGET_RUN set_variable PARAMETERS_FILE $WORKDIR parameters
Definition: diff-Tuna.sh:38
exit
exit
Definition: JPizza.sh:36
o
then fatal Wrong number of arguments fi JConvertDetectorFormat a o
Definition: compare-detector.sh:31
echo
then echo
Definition: JDOMDAQDriver.sh:56
JTOOLS::n
const int n
Definition: JPolint.hh:697
timer_start
then timer_start(JRandomTimesliceWriter-a ${DETECTOR}-o ${PIPE}-n ${NUMBER_OF_SLICES}-P ${PMT_FILE}-d ${DEBUG}-S ${SEED}-B"$BACKGROUND_HZ"-N"$RECYCLING"--!&&for((i=0;$i!=10;++i));do;echo 'EOF' > $PIPE;usleep 100000;done &&rm-f $PIPE)&JTimesliceReprocessor-a $
Definition: JTimesliceReprocessor.sh:57
set_variable
then echo Variable JPP_DIR undefined exit fi source $JPP_DIR setenv sh $JPP_DIR &dev null set_variable
Definition: JAcoustics.sh:21
T
do set_variable OUTPUT_DIRECTORY $WORKDIR T
Definition: JCalibrateHeight.sh:61
F
then awk F
Definition: getMIIseacurrent.sh:86
usage
* usage
Definition: JLigier-local.sh:34
APPLICATIONS
then fatal The output file must have the wildcard in the e g root fi typeset A PARAMETERS typeset a APPLICATIONS APPLICATIONS
Definition: JARCAMuonReconstruction.sh:50
$INPUT_FILE
then JFileTuna f $INPUT_FILE
Definition: file-Tuna.sh:66
OUTPUT_FILE
&set_variable OUTPUT_FILE
Definition: JARCACalibrateMuon.sh:29
INPUT_FILE
set_variable INPUT_FILE
Definition: JARCACalibrateMuon.sh:30
if
then fatal Wrong number of arguments fi set_variable DETECTOR $argv[1] set_variable TRIPOD $argv[2] set_variable TX $argv[3] set_variable TY $argv[4] if[[!-f $DETECTOR]]
Definition: JFootprint.sh:28
root
do JPlot2D f $WORKDIR detector root
Definition: JCompareDetector.sh:95
JRECONSTRUCTION::JHistory
JFIT::JHistory JHistory
Definition: JHistory.hh:354
d
then JMuonMCEvt f $INPUT_FILE o $INTERMEDIATE_FILE d
Definition: JMuonPath.sh:47
$script
then echo Variable JPP_DIR undefined exit fi source $JPP_DIR setenv sh $JPP_DIR set_variable DEBUG set_variable WORKDIR if do_usage *then usage $script(input file)+" fi if (( $ fatal "Invalid number of arguments." fi JCookie.sh set_array INPUT_FILES $argv[1
a
then JCalibrateToT a
Definition: JTuneHV.sh:116
$DETECTOR
then set_variable MODULE getModule a $DETECTOR L $STRING $FLOOR JEditDetector a $DETECTOR M $MODULE add $X o $DETECTOR else echo No update of detector $DETECTOR
Definition: module-Z:fit.sh:118
DETECTOR
set_variable DETECTOR
Definition: JARCACalibrateMuon.sh:31
detector
do set_variable DETECTOR_TXT $WORKDIR detector
Definition: detector-XY:fit1d.sh:38
JTOOLS::v
data_type v[N+1][M+1]
Definition: JPolint.hh:777
in
then fatal Wrong number of arguments fi set_variable DETECTOR $argv[1] set_variable INPUT_FILE $argv[2] eval JPrintDetector a $DETECTOR O IDENTIFIER eval JPrintDetector a $DETECTOR O SUMMARY JAcoustics sh $DETECTOR_ID source JAcousticsToolkit sh CHECK_EXIT_CODE typeset A EMITTERS get_tripods $WORKDIR tripod txt EMITTERS get_transmitters $WORKDIR transmitter txt EMITTERS for EMITTER in
Definition: JCanberra.sh:46
A
source $JPP_DIR setenv csh $JPP_DIR &dev null eval JShellParser o a A
Definition: JShellParser.csh:15
script
script
Definition: JAcoustics.sh:2
P
then $DIR JPlotNPE PDG P
Definition: JPlotNPE-PDG.sh:62
DEBUG
#define DEBUG(A)
Message macros.
Definition: JMessage.hh:62
Generated by   doxygen 1.8.5