7 # ------------------------------------------------------------------------------------------
9 # Script to run the JPP Track Reconstruction Chain for ORCA.
11 # ------------------------------------------------------------------------------------------
13 if [ -z $JPP_DIR ]; then
14 echo "Variable JPP_DIR undefined."
18 source $JPP_DIR/setenv.sh $JPP_DIR
19 source $JPP_DIR/externals/km3net-dataformat/definitions/reconstruction.sh
21 set_variable: DEBUG RECONSTRUCTION_DEBUG 1
22 set_variable: PDF RECONSTRUCTION_PDF $JPP_DATA/J%p.dat
23 set_variable: SUFFIX RECONSTRUCTION_SUFFIX jmuon
24 set_variable: NUMBER_OF_EVENTS RECONSTRUCTION_NUMBER_OF_EVENTS $((2**31))
25 set_variable+ OVERWRITE RECONSTRUCTION_OVERWRITE 1
26 set_variable: TMAX_S RECONSTRUCTION_TMAX_S 100
29 usage "$script <detector file> <input file> <output file> <PMT parameters file> [(calibration file)+]"
33 fatal "Wrong number of arguments."
36 set_variable DETECTOR $argv[1]
37 set_variable INPUT_FILE $argv[2]
38 set_variable OUTPUT_FILE $argv[3]
39 set_variable PMT $argv[4]
40 set_array CALIBRATION $argv[5,-1]
42 if [[ ! "$OUTPUT_FILE" == *"%"* ]]; then
43 fatal "The output file must have the wildcard % in the name, e.g. %.root"
47 typeset -a APPLICATIONS
52 if [[ "$OVERWRITE" == 1 ]]; then
54 ${OUTPUT_FILE/\%/jprefit} \
55 ${OUTPUT_FILE/\%/jsimplex} \
56 ${OUTPUT_FILE/\%/jstart} \
57 ${OUTPUT_FILE/\%/jgandalf} \
58 ${OUTPUT_FILE/\%/jstart2} \
59 ${OUTPUT_FILE/\%/jenergy} \
60 ${OUTPUT_FILE/\%/jhistory} \
61 ${OUTPUT_FILE/\%/${SUFFIX}} \
62 ${OUTPUT_FILE/\%/${SUFFIX}.offline}
66 if (( ${#CALIBRATION} != 0 )); then
67 getMechanics.sh $DETECTOR
71 # ------------------------------------------------------------------------------------------
73 # ------------------------------------------------------------------------------------------
95 -o ${OUTPUT_FILE/\%/jprefit} \
96 -n ${NUMBER_OF_EVENTS} \
97 -@ `make_equation PARAMETERS` \
102 APPLICATIONS+=$JMUONPREFIT
108 # ------------------------------------------------------------------------------------------
110 # ------------------------------------------------------------------------------------------
125 -f ${OUTPUT_FILE/\%/jprefit} \
126 -o ${OUTPUT_FILE/\%/jsimplex} \
127 -@ `make_equation PARAMETERS` \
132 rm -vf ${OUTPUT_FILE/\%/jprefit}
134 APPLICATIONS+=$JMUONSIMPLEX
140 # ------------------------------------------------------------------------------------------
142 # ------------------------------------------------------------------------------------------
159 -f ${OUTPUT_FILE/\%/jsimplex} \
160 -o ${OUTPUT_FILE/\%/jstart} \
162 -@ `make_equation PARAMETERS` \
167 rm -vf ${OUTPUT_FILE/\%/jsimplex}
169 APPLICATIONS+=$JMUONSTART
175 # ------------------------------------------------------------------------------------------
177 # ------------------------------------------------------------------------------------------
196 -f ${OUTPUT_FILE/\%/jstart} \
197 -o ${OUTPUT_FILE/\%/jgandalf} \
199 -@ `make_equation PARAMETERS` \
204 rm -vf ${OUTPUT_FILE/\%/jstart}
206 APPLICATIONS+=$JMUONGANDALF
212 # ------------------------------------------------------------------------------------------
214 # ------------------------------------------------------------------------------------------
231 -f ${OUTPUT_FILE/\%/jgandalf} \
232 -o ${OUTPUT_FILE/\%/jstart2} \
234 -@ `make_equation PARAMETERS` \
239 rm -vf ${OUTPUT_FILE/\%/jgandalf}
241 APPLICATIONS+=$JMUONSTART
247 # ------------------------------------------------------------------------------------------
249 # ------------------------------------------------------------------------------------------
269 -f ${OUTPUT_FILE/\%/jstart2} \
270 -o ${OUTPUT_FILE/\%/jenergy} \
272 -@ `make_equation PARAMETERS` \
277 rm -vf ${OUTPUT_FILE/\%/jstart2}
279 APPLICATIONS+=$JMUONENERGY
285 # ------------------------------------------------------------------------------------------
287 # ------------------------------------------------------------------------------------------
292 -f ${OUTPUT_FILE/\%/jenergy} \
293 -o ${OUTPUT_FILE/\%/jhistory} \
294 -A "$APPLICATIONS[*]" \
299 rm -vf ${OUTPUT_FILE/\%/jenergy}
305 # ------------------------------------------------------------------------------------------
306 # rename outputfile and remove intermediate files
307 # ------------------------------------------------------------------------------------------
312 ${OUTPUT_FILE/\%/jhistory} \
313 ${OUTPUT_FILE/\%/${SUFFIX}}
319 # ------------------------------------------------------------------------------------------
321 # ------------------------------------------------------------------------------------------
326 -f ${OUTPUT_FILE/\%/${SUFFIX}} \
327 -o ${OUTPUT_FILE/\%/${SUFFIX}.offline} \