JORCAEnergy.sh

Go to the documentation of this file.

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

# ------------------------------------------------------------------------------------------
#
#                          Utility script to run JEnergy for ORCA
#
# ------------------------------------------------------------------------------------------

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/start.root
set_variable  OUTPUT_FILE         $WORKDIR/energy.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       ROADWIDTH_M        50
set_variable       TIMEWINDOW_NS      -25 +25
set_variable       START              1e-3 1e-2
set_variable       K40                10e3 500 50 5 0.5
set_variable       EWINDOW            0 2
set_variable       NUMBER_OF_PREFITS  1

#   choice of M-Estimator

set_variable       NORMAL             0
set_variable       LORENTZIAN         1   # Preferred for ARCA
set_variable       LINEAR             2
set_variable       NONE               3   # Preferred for ORCA

timer_start

JEnergy                          \
    -a $DETECTOR                 \
    -f $INPUT_FILE               \
    -P $PDF                      \
    -o $OUTPUT_FILE              \
    -R $ROADWIDTH_M              \
    -T $TIMEWINDOW_NS            \
    -s $START                    \
    -x $EWINDOW                  \
    -B $K40                      \
    -N $NUMBER_OF_PREFITS        \
    -M $NONE                     \
    -d $DEBUG --!

timer_stop
timer_print