software/JCalibrate/JHobbit.sh

Go to the documentation of this file.

#!/bin/zsh
#
# \author mdejong
#
script=${0##*/}
 
# ------------------------------------------------------------------------------------------
#
#                         Script to run JHobbit.
#
# ------------------------------------------------------------------------------------------
 
if [ -z $JPP_DIR ]; then
    echo "Variable JPP_DIR undefined."
    exit
fi
source $JPP_DIR/setenv.sh $JPP_DIR
 
set_variable   DEBUG        2
set_variable   WORKDIR      ./
set_variable:  FITRANGE     CALIBRATE_FITRANGE      -7.5 +5.0
set_variable:  FUNCTION     CALIBRATE_FUNCTION      Gauss
set_variable:  PALETTE      GRAPHICS_PALETTE        57
set_variable:  FORMAT       GRAPHICS_FORMAT         gif
set_variable+  BATCH        GRAPHICS_BATCH          -B
 
if do_usage $*; then
    usage "$script <detector file> (input file)+"
fi
 
if (( $# < 2 )); then
    fatal "Wrong number of arguments."
fi
 
set_variable DETECTOR        $argv[1]
set_variable INPUT_FILES     $argv[2,-1]
 
if ( ! reuse_file $WORKDIR/hobbit.root ); then
 
    JHobbit \
        -a $DETECTOR                  \
        -f${^INPUT_FILES}             \
        -o $WORKDIR/hobbit.root       \
        -F $FUNCTION                  \
        -O L                          \
        -d $DEBUG
fi
 
eval `JPrintDetector -a km3net_reference.detx -O SUMMARY`
 
 
if (( ${#STRINGS} <= 15 )); then
    let "X = 600"
    let "Y = 600"
else
    let "X = 600 + 6 * (${#STRINGS} - 15)"
    let "Y = 600"
fi
 
JPlot1D \
    -w ${X}x${Y}                    \
    -f $WORKDIR/hobbit.root:hc      \
    -> "index"                     \
    -\^ "\chi^{2}/NDF"              \
    -T ""                           \
    -o chi2.$FORMAT                 $BATCH    
 
JPlot2D \
    -w ${X}x${Y}                    \
    -f $WORKDIR/hobbit.root:hi      \
    -> "string"                    \
    -< "floor"                     \
    -\^ "#Deltat [ns]"              \
    -T ""                           \
    -O COLZ                         \
    -p $PALETTE                     \
    -o detector.$FORMAT             $BATCH