JORCAGandalf.sh

Go to the documentation of this file.

#!/bin/zsh
#
# \author lquinn, bofearraigh
#
script=${0##*/}

# ------------------------------------------------------------------------------------------
#
#                         Utility script to run JGandalf.
#
# ------------------------------------------------------------------------------------------

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

source $JPP_DIR/setenv.sh $JPP_DIR

# default values for testing

set_variable  DEBUG               2
set_variable  WORKDIR             $JPP_DATA
set_variable  DETECTOR            $WORKDIR/km3net_reference.detx
set_variable  INPUT_FILE          $WORKDIR/prefit.root
set_variable  OUTPUT_FILE         $WORKDIR/gandalf.root
set_variable  PDF                 $WORKDIR/J%p.dat

if ( do_usage $* ); then
    usage "$script [detector file [input file [output file [PDF file descriptor]]]]"
fi

case $# in
    4) set_variable PDF         $4;&
    3) set_variable OUTPUT_FILE $3;&
    2) set_variable INPUT_FILE  $2;&
    1) set_variable DETECTOR    $1;;
esac

print_variable     DETECTOR INPUT_FILE OUTPUT_FILE PDF
check_input_file   $DETECTOR $INPUT_FILE
check_output_file  $OUTPUT_FILE

#   default input values

set_variable TTS_NS             2.0
set_variable E_GEV              10.0
set_variable R_HZ               10e3
set_variable ROADWIDTH_M        50
set_variable ZRANGE_M           0 9
set_variable NUMBER_OF_PREFITS  0

timer_start

JGandalf \
    -a $DETECTOR               \
    -f $INPUT_FILE             \
    -P $PDF                    \
    -o $OUTPUT_FILE            \
    -T $TTS_NS                 \
    -E $E_GEV                  \
    -B $R_HZ                   \
    -R $ROADWIDTH_M            \
    -z $ZRANGE_M               \
    -N $NUMBER_OF_PREFITS      \
    -d $DEBUG --!

timer_stop
timer_print