JORCAShowerFit.sh

Go to the documentation of this file.

#!/bin/zsh
#
# \author adomi
#
version=1.0
script=${0##*/}

# ----------------------------------------------------------------------------------
#
#                    Utility script to test JORCAShowerFit.
#
# ----------------------------------------------------------------------------------


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

source $JPP_DIR/setenv.sh $JPP_DIR

set_variable  DEBUG          0
set_variable  WORKDIR        $JPP_DATA
set_variable  INPUT_FILE     $WORKDIR/shower_positionfit.root
set_variable  PDF_FILE       $WORKDIR/PDFShowerFit%p.dat
set_variable  OUTPUT_FILE    $WORKDIR/shower_fit.root
set_variable  DETECTOR       $WORKDIR/km3net_reference.detx

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

case $# in
    5) DEBUG=$5;&
    4) OUTPUT_FILE=$4;&
    3) PDF_FILE=$3;&
    2) INPUT_FILE=$2;&
    1) DETECTOR=$1;;
esac

print_variable     DETECTOR INPUT_FILE PDF_FILE OUTPUT_FILE DEBUG
check_input_file   $DETECTOR $INPUT_FILE
check_output_file  $OUTPUT_FILE

if (( 1 )); then

    # default input values                                                                                         
    set_variable ORCA_SHOWER_PARAMETERS  $JPP_DATA/showerReconstruction_parameters_orca.txt

    timer_start

    JORCAShowerFit                        \
        -a $DETECTOR                      \
        -f $INPUT_FILE                    \
        -F $PDF_FILE                      \
        -o $OUTPUT_FILE                   \
        -@ ${ORCA_SHOWER_PARAMETERS}      \
        -d ${DEBUG} --!

    timer_stop
    timer_print

fi