JSirene.sh

Go to the documentation of this file.

#!/bin/zsh
#
# \author mdejong
#
# utility script for Sirene Monte Carlo.
script=${0##*/}

# ------------------------------------------------------------------------------------------
#
#                            Example script to run JSirene.
#
# ------------------------------------------------------------------------------------------

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

source $JPP_DIR/setenv.sh $JPP_DIR

set_variable  DEBUG          ${SIRENE_DEBUG:-1}
set_variable  WORKDIR        ./
set_variable  CDF            $JPP_DATA/I%p.dat
set_variable  INPUT_FILE     $JPP_DATA/genhen.km3net_wpd_V2_0.evt.gz
set_variable  OUTPUT_FILE    $WORKDIR/sirene.root
set_variable  DETECTOR       $JPP_DATA/km3net_reference.detx
set_variable  KEEP           ${SIRENE_KEEP+-k}
set_variable  SEED           ${SIRENE_SEED:-0}
set_variable  FACTOR         ${SIRENE_FACTOR:-1.0}

if ( do_usage $* ); then
    usage "$script [detector file [input file [output file [CDF file descriptor]]]]"\
          "\nNote that if more than one input file is specified, all other arguments must be provided."
fi

case $# in
    4)  set_variable CDF          $4;&
    3)  set_variable OUTPUT_FILE  $3;&
    2)  set_variable INPUT_FILE   $2;&
    1)  set_variable DETECTOR     $1;;
    0)  ;;
    *)  set_variable DETECTOR     $1
        set_variable INPUT_FILE   $argv[2,-3]
        set_variable OUTPUT_FILE  $argv[-2]
        set_variable CDF          $argv[-1];;
esac
            
# Sanity checks

print_variable     DETECTOR INPUT_FILE OUTPUT_FILE CDF

for TYPE in 1 2 5 6 13 14; do
    if [[ ! -f ${CDF/\%/${TYPE}} ]]; then
        error "Missing CDF files."
        error "Run: JMakePDF.sh"
        JMakePDF.sh -h  
        exit
    fi
done

if [[ ! -f $DETECTOR ]]; then
    error "Missing detector file."
    error "Run: JDetector.sh"
    JDetector.sh -h
    exit
fi

check_input_file   $DETECTOR $INPUT_FILE
check_output_file  $OUTPUT_FILE

timer_start

JSirene \
    -F $CDF              \
    -a $DETECTOR         \
    -f $INPUT_FILE       \
    -o $OUTPUT_FILE      \
    -@ "Tmax_ns=2.0;"    \
    -U $FACTOR           \
    $KEEP                \
    -S $SEED             \
    -d $DEBUG --!

timer_stop
timer_print

JPrintMeta \
    -f $OUTPUT_FILE

JPrintSirene \
    -f $OUTPUT_FILE

JPrintTree \
    -f $OUTPUT_FILE