JARCAGandalf.sh

Go to the documentation of this file.

#!/bin/zsh
#
# \author mdejong
#
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               1
set_variable  WORKDIR             $JPP_DATA
set_variable  DETECTOR            $WORKDIR/km3net_reference.detx
set_variable  INPUT_FILE          $WORKDIR/simplex.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              1e3
set_variable R_HZ               6e3
set_variable ROADWIDTH_M        175

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}            \
    -d ${DEBUG} --!

timer_stop
timer_print