JARCAEnergy.sh

Go to the documentation of this file.

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

# ------------------------------------------------------------------------------------------
#
#                         Utility script to run JEnergy.
#
# ------------------------------------------------------------------------------------------

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/gandalf.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  ENERGY_CORRECTION   energy_correction_arca.txt

#   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
set_variable       NUMBER_OF_PREFITS  1

timer_start

JEnergy \
    -a $DETECTOR                 \
    -f $INPUT_FILE               \
    -P $PDF                      \
    -E $ENERGY_CORRECTION        \
    -o $OUTPUT_FILE              \
    -M $LORENTZIAN               \
    -N $NUMBER_OF_PREFITS        \
    -d $DEBUG --!

timer_stop
timer_print