4 if [ -z $JPP_DIR ]; then
5 echo "Variable JPP_DIR undefined."
9 source $JPP_DIR/setenv.sh $JPP_DIR >& /dev/null
11 set_variable: DEBUG ACOUSTICS_DEBUG 2
12 set_variable: WORKDIR ACOUSTICS_WORKDIR `pwd`
13 set_variable+ OVERWRITE ACOUSTICS_OVERWRITE Y
14 set_variable FORMULA "[0] * exp(-0.5 * (x-[1])*(x-[1]) / ([2]*[2])) * exp(-0.5 * (y-[3])*(y-[3]) / ([4]*[4]))"
15 set_variable: FORMAT GRAPHICS_FORMAT gif
16 set_variable+ BATCH GRAPHICS_BATCH -B
19 usage "$script <detector file> <tripod file> (input file)+"\
20 "\nInput files correspond to the output of JKatoomba[.sh]."\
21 "\nThe environment variable ACOUSTICS_OVERWRITE should be set to overwrite the detector and tripod file."
25 fatal "Wrong number of arguments."
28 set_variable DETECTOR $argv[1]
29 set_variable TRIPOD $argv[2]
30 set_array INPUT_FILES $argv[3,-1]
32 eval `JPrintDetector -a $DETECTOR -O IDENTIFIER`
33 eval `JPrintDetector -a $DETECTOR -O CAN`
35 JAcoustics.sh $DETECTOR_ID
39 -x "500 -20.0 +20.0" \
40 -y "500 -20.0 +20.0" \
42 -o $WORKDIR/footprint.root \
46 -f $WORKDIR/footprint.root:% \
53 -o footprint.$FORMAT $BATCH
55 JPrintMaximum2D -f $WORKDIR/footprint.root:% | read TITLE X Y
63 -f $WORKDIR/footprint.root:% \
64 -o $WORKDIR/fit.root \
68 -@ "p0 = GetMaximum" \
82 -f $WORKDIR/fit.root:% \
93 set_array TX `JPrintFit -f $WORKDIR/fit.root:% -@ p1`
94 set_array TY `JPrintFit -f $WORKDIR/fit.root:% -@ p3`
95 set_variable RESULT `JPrintFit -f $WORKDIR/fit.root:% -F "GetChisquare/GetNDF"`
97 TX[1]=$(($TX[1]*1.0e-3))
98 TY[1]=$(($TY[1]*1.0e-3))
100 printf "Tx = %9.6f Ty = %9.6f chi2/NDF = %7.3f\n" $TX[1] $TY[1] $RESULT
102 if [[ "$OVERWRITE" == "Y" ]]; then
103 JFootprint.sh $DETECTOR $TRIPOD $TX[1] $TY[1]