footprint.sh

Go to the documentation of this file.

#!/usr/bin/env zsh
script=${0##*/}
 
if [ -z $JPP_DIR ]; then
    echo "Variable JPP_DIR undefined."
    exit
fi
 
source $JPP_DIR/setenv.sh $JPP_DIR >& /dev/null
 
set_variable:  DEBUG         ACOUSTICS_DEBUG        2
set_variable:  WORKDIR       ACOUSTICS_WORKDIR      `pwd`
set_variable+  OVERWRITE     ACOUSTICS_OVERWRITE    Y
set_variable   FORMULA       "[0] * exp(-0.5 * (x-[1])*(x-[1]) / ([2]*[2])) * exp(-0.5 * (y-[3])*(y-[3]) / ([4]*[4]))"  
set_variable:  FORMAT        GRAPHICS_FORMAT        gif
set_variable+  BATCH         GRAPHICS_BATCH         -B
 
if do_usage $*; then
    usage "$script <detector file> <tripod file> (input file)+"\
          "\nInput files correspond to the output of JKatoomba[.sh]."\
          "\nThe environment variable ACOUSTICS_OVERWRITE should be set to overwrite the detector and tripod file."
fi
 
if (( $# < 3 )); then
    fatal "Wrong number of arguments."
fi
 
set_variable  DETECTOR            $argv[1]
set_variable  TRIPOD              $argv[2]
set_array     INPUT_FILES         $argv[3,-1]
 
eval `JPrintDetector -a $DETECTOR -O IDENTIFIER`
eval `JPrintDetector -a $DETECTOR -O CAN`
 
JAcoustics.sh   $DETECTOR_ID
 
JFootprint \
    -f${^INPUT_FILES}                         \
    -x "500 -20.0 +20.0"                      \
    -y "500 -20.0 +20.0"                      \
    -Z $CAN_ZMAX_M                            \
    -o $WORKDIR/footprint.root                \
    -d $DEBUG
 
JPlot2D \
    -f $WORKDIR/footprint.root:%              \
    -> "Tx [mrad]"                           \
    -< "Ty [mrad]"                           \
    -O COLZ                                   \
    -N "X 505"                                \
    -N "Y 505"                                \
    -T ""                                     \
    -o footprint.$FORMAT                      $BATCH
 
JPrintMaximum2D -f $WORKDIR/footprint.root:%  |  read TITLE X Y
 
let "XMIN = $X - 2.0"
let "XMAX = $X + 2.0"
let "YMIN = $Y - 2.0"
let "YMAX = $Y + 2.0"
 
JFit2D \
    -f $WORKDIR/footprint.root:%              \
    -o $WORKDIR/fit.root                      \
    -F "$FORMULA"                             \
    -x "$XMIN $XMAX"                          \
    -y "$YMIN $YMAX"                          \
    -@ "p0 = GetMaximum"                      \
    -@ "p1 = $X"                              \
    -@ "p2 = 1.0"                             \
    -@ "p3 = $Y"                              \
    -@ "p4 = 1.0"                             \
    -E "p0 = 0.1"                             \
    -E "p1 = 0.02"                            \
    -E "p2 = 0.02"                            \
    -E "p3 = 0.02"                            \
    -E "p4 = 0.02"                            \
    -O "L"                                    \
    -d $DEBUG --!
 
JPlot2D \
    -f $WORKDIR/fit.root:%                    \
    -x "$XMIN $XMAX"                          \
    -y "$YMIN $YMAX"                          \
    -> "Tx [mrad]"                           \
    -< "Ty [mrad]"                           \
    -O COLZ                                   \
    -N "X 505"                                \
    -N "Y 505"                                \
    -T ""                                     \
    -o fit.$FORMAT                            $BATCH
 
set_array    TX     `JPrintFit    -f $WORKDIR/fit.root:% -@ p1`
set_array    TY     `JPrintFit    -f $WORKDIR/fit.root:% -@ p3`
set_variable RESULT `JPrintFit    -f $WORKDIR/fit.root:% -F "GetChisquare/GetNDF"`
 
TX[1]=$(($TX[1]*1.0e-3))
TY[1]=$(($TY[1]*1.0e-3))
 
printf "Tx = %9.6f Ty = %9.6f chi2/NDF = %7.3f\n" $TX[1] $TY[1] $RESULT
 
if [[ "$OVERWRITE" == "Y" ]]; then
    JFootprint.sh  $DETECTOR  $TRIPOD  $TX[1]  $TY[1]
fi