7# ------------------------------------------------------------------------------------------
9# Script to run the JPP Track Reconstruction Chain for ARCA.
11# ------------------------------------------------------------------------------------------
13JPP_DIR=/localstore/km3net/mjg/GIT/storage/
15if [ -z $JPP_DIR ]; then
16 echo "Variable JPP_DIR undefined."
20source $JPP_DIR/setenv.sh $JPP_DIR
21source $JPP_DIR/externals/km3net-dataformat/definitions/reconstruction.sh
23set_variable: DEBUG RECONSTRUCTION_DEBUG 1
24set_variable: PDF RECONSTRUCTION_PDF $JPP_DATA/J%p.dat
25set_variable: ENERGY_CORRECTION RECONSTRUCTION_ENERGY_CORRECTION energy_correction_arca.txt
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 6.0e3
33 usage "$script <detector file> <input file> <output file> <PMT parameters file> [(calibration file)+]"
37 fatal "Wrong number of arguments."
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]
46if (( ${#CALIBRATION} != 0 )); then
47 getMechanics.sh $DETECTOR
54 prefit.gridAngle_deg 1
56 prefit.numberOfOutliers 3
57 prefit.numberOfPrefits 50
59 prefit.numberOfDZMax 1
60 prefit.TMaxLocal_ns 18
62 prefit.roadWidth_m 200.0)
66 simplex.TMaxLocal_ns 15.0
68 simplex.roadWidth_m 200.0
69 simplex.numberOfPrefits 0)
75 gandalf.TMax_ns +450.0
77 gandalf.roadWidth_m 175.0
78 gandalf.numberOfPrefits 0)
81 start.roadWidth_m 160.0
87 start.numberOfPrefits 1)
90 energy.roadWidth_m 300.0
97 energy.numberOfPrefits 1)
99eval $JPP_DIR/examples/JReconstruction/JARCAMuonReconstruction \
106 -E $ENERGY_CORRECTION \
107 -n ${NUMBER_OF_EVENTS} \
108 -N ${NUMBER_OF_THREADS} \
109 -@ `make_equation PARAMETERS` \