footprint.sh

Go to the documentation of this file.

#!/bin/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:-./}
set_variable  FIT                 ${ACOUSTICS_FIT:-2}
set_variable  FORMULA             "[0] * exp(-0.5 * (x-[1])*(x-[1]) / ([2]*[2])) * exp(-0.5 * (y-[3])*(y-[3]) / ([4]*[4]))"

if ( do_usage $* ); then
    usage "$script <detector file> <Tx> <Ty> (input file)+"
fi

if (( $# < 4 )); then
    fatal "Wrong number of arguments."
fi

set_variable  DETECTOR            $argv[1]
set_variable  TX                  $argv[2]
set_variable  TY                  $argv[3]
set_array     INPUT_FILES         $argv[4,-1]

eval `JPrintDetector -a $DETECTOR -O IDENTIFIER`

source JAcoustics.sh   $DETECTOR_ID

CHECK_EXIT_CODE

set_variable  TMPDIR              $WORKDIR/.$$
mkdir -p     $TMPDIR

cp \
    $WORKDIR/tripod.txt                       \
    $TMPDIR/tripod.txt

JConvertDetectorFormat \
    -a $DETECTOR                              \
    -o $TMPDIR/detector.detx --!

JFootprint.sh \
    $TMPDIR/detector.detx                     \
    $TMPDIR/tripod.txt                        \
    $TX                                       \
    $TY

JKatoomba \
    -a $TMPDIR/detector.detx                  \
    -f "$INPUT_FILES[*]"                      \
    -o $TMPDIR/katoomba.root                  \
    -T $WORKDIR/tripod.txt                    \
    -V $WORKDIR/sound_velocity.txt            \
    -@ $WORKDIR/acoustics_fit_parameters.txt  \
    -@ "stdev      = 50"                      \
    -@ "mestimator = 1"                       \
    -F $FIT                                   \
    -u                                        \
    -d $DEBUG --!

JParramatta \
    -f $TMPDIR/katoomba.root                  \
    -o $TMPDIR/parramatta.root

set_array     X  `JPrintResult -f $TMPDIR/parramatta.root:% -F "GetMean(1)" -F "0.5*GetRMS(1)"`
set_array     Y  `JPrintResult -f $TMPDIR/parramatta.root:% -F "GetMean(2)" -F "0.5*GetRMS(2)"`

let XMIN="$X[1] - 0.5"
let XMAX="$X[1] + 0.5"
let YMIN="$Y[1] - 0.5"
let YMAX="$Y[1] + 0.5"

JFit2D \
    -f $TMPDIR/parramatta.root:%              \
    -o $TMPDIR/fit.root                       \
    -F "$FORMULA"                             \
    -x "$XMIN $XMAX"                          \
    -y "$YMIN $YMAX"                          \
    -@ "p0 = GetMaximum"                      \
    -@ "p1 = $X[1]"                           \
    -@ "p2 = $X[2]"                           \
    -@ "p3 = $Y[1]"                           \
    -@ "p4 = $Y[2]"                           \
    -d $DEBUG --!

JPlot2D \
    -f $TMPDIR/parramatta.root:%              \
    -O COLZ                                   \
    -\> "Tx [mrad]"                           \
    -\< "Ty [mrad]"                           \
    -N "X 505"                                \
    -N "Y 505"                                \
    -T ""                                     \
    -o footprint.gif

set_array    CHI2 `JPrintResult -f $TMPDIR/katoomba.root:chi2 -F "GetMean" -F "GetRMS / sqrt(GetEntries)"`
set_array    TX   `JPrintFit -f $TMPDIR/fit.root:% -@ p1`
set_array    TY   `JPrintFit -f $TMPDIR/fit.root:% -@ p3`

TX[1]=$(($TX[1]*1.0e-3))
TY[1]=$(($TY[1]*1.0e-3))

printf "Tx = %9.6f Ty = %9.6f chi2 = %7.3f %7.3f\n" $TX[1] $TY[1] ${CHI2[*]}

rm -rf $TMPDIR

---

Generated by   1.8.5