examples/JReconstruction/compare-detector.sh

Go to the documentation of this file.

#!/usr/bin/env 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