examples/JReconstruction/compare-detector.sh

Go to the documentation of this file.

#!/bin/zsh

if [ -z $JPP_DIR ]; then
    echo "Variable JPP_DIR undefined."
    exit
fi

source $JPP_DIR/setenv.sh $JPP_DIR

zmodload zsh/mathfunc

set_variable   DEBUG       2
set_variable   WORKDIR     ${TMPDIR:-/tmp}/
set_variable:  PRECISION   RECONSTRUCTION_PRECISION    1.0
set_variable:  FORMAT      GRAPHICS_FORMAT             gif
set_variable+  BATCH       GRAPHICS_BATCH              -B

if do_usage $*; then
    usage "$script <detector file> <detector file>"
fi

if (( $# != 2 )); then
    fatal "Wrong number of arguments."
fi

DETECTOR=($argv[1] $argv[2])

JConvertDetectorFormat \
    -a $DETECTOR[2]                   \
    -o $WORKDIR/detector.datx

# correct global time offset

let "T  = 0.0"

for (( i = 0; $i != 2; ++i )); do

    JCompareDetector \
        -a $DETECTOR[1]                   \
        -b $WORKDIR/detector.datx         \
        -o $WORKDIR/compare.root          \
        -d0                               >& /dev/null 

    JPrintResult -f $WORKDIR/compare.root:T2 -F "GetMinimum"  |  read ZMIN
    JPrintResult -f $WORKDIR/compare.root:T2 -F "GetMaximum"  |  read ZMAX

    let "Z  = 0.5 * ($ZMIN + $ZMAX)"
    let "T += $Z"

    JEditDetector \
        -a $WORKDIR/detector.datx         \
        -o $WORKDIR/detector.datx         \
        -s "-1 add $Z"                    \
        -d 0                              >& /dev/null 
done

printf "Correct global time offset %9.6f [ns].\n" $T
printf "Test with precision %9.6f [ns].\n" $PRECISION

JCompareDetector \
    -a $DETECTOR[1]                   \
    -b $WORKDIR/detector.datx         \
    -o $WORKDIR/compare.root          \
    -p "tcal = $PRECISION"            \
    -d1

JPlot2D \
    -w 1200x600                       \
    -f $WORKDIR/compare.root:T2       \
    -\> "string"                      \
    -\< "floor"                       \
    -O COLZ                           \
    -T "#Deltat [ns]"                 \
    -o detector.$FORMAT               $BATCH