7# ------------------------------------------------------------------------------------------
9# Script to run the JPP Track Reconstruction Chain for ARCA.
11# ------------------------------------------------------------------------------------------
13if [ -z $JPP_DIR ]; then
14 echo "Variable JPP_DIR undefined."
18source $JPP_DIR/setenv.sh $JPP_DIR
19source $JPP_DIR/externals/km3net-dataformat/definitions/reconstruction.sh
21set_variable: DEBUG RECONSTRUCTION_DEBUG 1
22set_variable: PDF RECONSTRUCTION_PDF $JPP_DATA/J%p.dat
23set_variable: SUFFIX RECONSTRUCTION_SUFFIX jmuon
24set_variable: ENERGY_CORRECTION RECONSTRUCTION_ENERGY_CORRECTION energy_correction_arca.txt
25set_variable: NUMBER_OF_EVENTS RECONSTRUCTION_NUMBER_OF_EVENTS $((2**31))
26set_variable: NUMBER_OF_THREADS RECONSTRUCTION_NUMBER_OF_THREADS 0
27set_variable+ OVERWRITE RECONSTRUCTION_OVERWRITE 1
28set_variable: TMAX_S RECONSTRUCTION_TMAX_S 100
29set_variable: RATE_HZ RECONSTRUCTION_RATE_HZ 6.0e3
32 usage "$script <detector file> <input file> <output file> <PMT parameters file> [(calibration file)+]"
36 fatal "Wrong number of arguments."
39set_variable DETECTOR $argv[1]
40set_variable INPUT_FILE $argv[2]
41set_variable OUTPUT_FILE $argv[3]
42set_variable PMT $argv[4]
43set_array CALIBRATION $argv[5,-1]
45if [[ ! "$OUTPUT_FILE" == *"%"* ]]; then
46 fatal "The output file must have the wildcard % in the name, e.g. %.root"
50typeset -a APPLICATIONS
55if [[ "$OVERWRITE" == 1 ]]; then
57 ${OUTPUT_FILE/\%/jprefit} \
58 ${OUTPUT_FILE/\%/jsimplex} \
59 ${OUTPUT_FILE/\%/jgandalf} \
60 ${OUTPUT_FILE/\%/jstart} \
61 ${OUTPUT_FILE/\%/jenergy} \
62 ${OUTPUT_FILE/\%/jhistory} \
63 ${OUTPUT_FILE/\%/${SUFFIX}} \
64 ${OUTPUT_FILE/\%/${SUFFIX}.offline}
67if (( ${#CALIBRATION} != 0 )); then
68 getMechanics.sh $DETECTOR
72# ------------------------------------------------------------------------------------------
74# ------------------------------------------------------------------------------------------
95 -o ${OUTPUT_FILE/\%/jprefit} \
96 -n ${NUMBER_OF_EVENTS} \
97 -N ${NUMBER_OF_THREADS} \
98 -@ `make_equation PARAMETERS` \
103APPLICATIONS+=$JMUONPREFIT
109# ------------------------------------------------------------------------------------------
111# ------------------------------------------------------------------------------------------
126 -f ${OUTPUT_FILE/\%/jprefit} \
127 -o ${OUTPUT_FILE/\%/jsimplex} \
128 -N ${NUMBER_OF_THREADS} \
129 -@ `make_equation PARAMETERS` \
134rm -vf ${OUTPUT_FILE/\%/jprefit}
136APPLICATIONS+=$JMUONSIMPLEX
142# ------------------------------------------------------------------------------------------
144# ------------------------------------------------------------------------------------------
161 -f ${OUTPUT_FILE/\%/jsimplex} \
162 -o ${OUTPUT_FILE/\%/jgandalf} \
164 -N ${NUMBER_OF_THREADS} \
165 -@ `make_equation PARAMETERS` \
170rm -vf ${OUTPUT_FILE/\%/jsimplex}
172APPLICATIONS+=$JMUONGANDALF
178# ------------------------------------------------------------------------------------------
180# ------------------------------------------------------------------------------------------
197 -f ${OUTPUT_FILE/\%/jgandalf} \
198 -o ${OUTPUT_FILE/\%/jstart} \
200 -N ${NUMBER_OF_THREADS} \
201 -@ `make_equation PARAMETERS` \
206rm -vf ${OUTPUT_FILE/\%/jgandalf}
208APPLICATIONS+=$JMUONSTART
214# ------------------------------------------------------------------------------------------
216# ------------------------------------------------------------------------------------------
234 -f ${OUTPUT_FILE/\%/jstart} \
235 -o ${OUTPUT_FILE/\%/jenergy} \
237 -E $ENERGY_CORRECTION \
238 -@ `make_equation PARAMETERS` \
239 -N ${NUMBER_OF_THREADS} \
244rm -vf ${OUTPUT_FILE/\%/jstart}
246APPLICATIONS+=$JMUONENERGY
252# ------------------------------------------------------------------------------------------
254# ------------------------------------------------------------------------------------------
259 -f ${OUTPUT_FILE/\%/jenergy} \
260 -o ${OUTPUT_FILE/\%/jhistory} \
261 -A "$APPLICATIONS[*]" \
266rm -vf ${OUTPUT_FILE/\%/jenergy}
272# ------------------------------------------------------------------------------------------
273# rename outputfile and remove intermediate files
274# ------------------------------------------------------------------------------------------
279 ${OUTPUT_FILE/\%/jhistory} \
280 ${OUTPUT_FILE/\%/${SUFFIX}}
286# ------------------------------------------------------------------------------------------
288# ------------------------------------------------------------------------------------------
293 -f ${OUTPUT_FILE/\%/${SUFFIX}} \
294 -o ${OUTPUT_FILE/\%/${SUFFIX}.offline} \