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]]]]"\
          "\nAuxiliary script for simulation of detector response."\
          "\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