JSirene.sh

Go to the documentation of this file.

#!/usr/bin/env zsh
#
# \author mdejong
#
script=${0##*/}
 
# ------------------------------------------------------------------------------------------
#
#                            Auxiliary 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+  WRITE_EMSHOWERS  SIRENE_WRITE_EMSHOWERS                   -s
set_variable:  NUMBER_OF_HITS   SIRENE_NUMBER_OF_HITS                    1
set_variable:  SEED             SIRENE_SEED                              0
set_variable:  FACTOR           SIRENE_FACTOR                            1.0
set_variable:  CDF              SIRENE_CDF                               $JPP_DATA/I%p.dat
 
if do_usage $*; then
    usage "$script <detector file> (input file)+ <output file>"\
          "\nAuxiliary script for simulation of detector response."
fi
 
if (( $# < 3 )); then
    fatal "Wrong number of arguments."
fi
 
set_variable DETECTOR     $argv[1]
set_variable INPUT_FILE   $argv[2,-2]
set_variable OUTPUT_FILE  $argv[-1]
            
# Sanity checks
 
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
 
timer_start
 
JSirene \
    -F $CDF              \
    -a $DETECTOR         \
    -f $INPUT_FILE       \
    -o $OUTPUT_FILE      \
    -@ "Tmax_ns  = 2.0;"     \
    -@ "Nmax_PMT = 1000000;" \
    -U $FACTOR           \
    $WRITE_EMSHOWERS     \
    -N $NUMBER_OF_HITS   \
    -S $SEED             \
    -d $DEBUG --!
 
timer_stop
timer_print
 
JPrintMeta \
    -f $OUTPUT_FILE
 
JPrintSirene \
    -f $OUTPUT_FILE
 
JPrintTree \
    -f $OUTPUT_FILE