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   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> <Tx> <Ty> (input file)+"\
          "\nThe tilt angles should be given in units [rad]."
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`

JAcoustics.sh   $DETECTOR_ID

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

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

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

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

JKatoomba \
    -a $TMPDIR/detector.datx                  \
    -f "$INPUT_FILES[*]"                      \
    -o $TMPDIR/katoomba.root                  \
    -T $TMPDIR/tripod.txt                     \
    -V $WORKDIR/sound_velocity.txt            \
    -@ $WORKDIR/acoustics_fit_parameters.txt  \
    -C -JEvent                                \
    -u                                        \
    -d $DEBUG --!

JFootprint \
    -f $TMPDIR/katoomba.root                  \
    -o $TMPDIR/footprint.root

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

JProject2D -f $TMPDIR/footprint.root:% -P X -y "-1e3 +1e3" -o $TMPDIR/px.root
JProject2D -f $TMPDIR/footprint.root:% -P Y -x "-1e3 +1e3" -o $TMPDIR/py.root

set_array     X  `JPrintQuantiles -f $TMPDIR/px.root:% -Q "0.45 0.55"`
set_array     Y  `JPrintQuantiles -f $TMPDIR/py.root:% -Q "0.45 0.55"`

let "XMIN = $X[1] - $((2.0 * ($X[2] - $X[1])))"
let "XMAX = $X[2] + $((2.0 * ($X[2] - $X[1])))"
let "YMIN = $Y[1] - $((2.0 * ($X[2] - $X[1])))"
let "YMAX = $Y[2] + $((2.0 * ($X[2] - $X[1])))"

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

JPlot2D \
    -f $TMPDIR/fit.root:%                     \
    -x "$XMIN $XMAX"                          \
    -y "$YMIN $YMAX"                          \
    -\> "Tx [mrad]"                           \
    -\< "Ty [mrad]"                           \
    -O CONT3                                  \
    -N "X 505"                                \
    -N "Y 505"                                \
    -T ""                                     \
    -o fit.$FORMAT                            $BATCH

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