JMonteCarloFileWriter.sh

Go to the documentation of this file.

#!/usr/bin/env zsh
#
# \author mdejong
#
# utility script for Sirene Monte Carlo.
version=1.0
script=${0##*/}
 
# ------------------------------------------------------------------------------------------
#
#                        Example script to test ASCII I/O of Monte Carlo files.
#
# ------------------------------------------------------------------------------------------
 
 
if [ -z $JPP_DIR ]; then
    echo "Variable JPP_DIR undefined."
    exit
fi
 
 
source $JPP_DIR/setenv.sh $JPP_DIR
 
set_variable   DEBUG             1                                     
set_variable   WORKDIR           ${TMPDIR:-/tmp}/                                 
set_variable   INPUT_FILE        $JPP_DATA/genhen.km3net_wpd_V2_0.evt  
set_variable   DETECTOR          $JPP_DATA/km3net_reference.detx       
set_variable   NUMBER_OF_EVENTS  100                                   
set_variable   SEED              12345678                              
    
 
if do_usage $*; then
    usage "$script [detector file [input file [working directory]]]"
fi
 
case $# in 
    3) set_variable  WORKDIR     $3;&
    2) set_variable  INPUT_FILE  $2;&
    1) set_variable  DETECTOR    $1;;
esac
 
 
if (( 1 )); then
 
    if [[ -f $WORKDIR/sirene.root ]]; then
        rm -i $WORKDIR/sirene.root
    fi
    
    if [[ ! -f $WORKDIR/sirene.root ]]; then
 
        JSirene \
            -F $JPP_DATA/I%p.dat     \
            -a $DETECTOR             \
            -f $INPUT_FILE           \
            -o $WORKDIR/sirene.root  \
            -n $NUMBER_OF_EVENTS     \
            -S $SEED                 \
            -d $DEBUG
    fi
fi
 
if (( 1 )); then
 
    if [[ -f $WORKDIR/sirene.evt ]]; then
        rm -i $WORKDIR/sirene.evt
    fi
    
    if [[ ! -f $WORKDIR/sirene.evt ]]; then
 
        JSirene \
            -F $JPP_DATA/I%p.dat     \
            -a $DETECTOR             \
            -f $INPUT_FILE           \
            -o $WORKDIR/sirene.evt   \
            -n $NUMBER_OF_EVENTS     \
            -S $SEED                 \
            -d $DEBUG            
    fi
fi
 
 
rm -f $WORKDIR/sirene_root.txt
rm -f $WORKDIR/sirene_evt.txt
 
echo "Print $WORKDIR/sirene.root"
JPrint \
    -f $WORKDIR/sirene.root          \
    -n $NUMBER_OF_EVENTS             \
    > $WORKDIR/sirene_root.txt
 
echo "Print $WORKDIR/sirene.evt"
JPrint \
    -f $WORKDIR/sirene.evt           \
    -n $NUMBER_OF_EVENTS             \
    > $WORKDIR/sirene_evt.txt
 
 
echo "Comparing $WORKDIR/sirene_root.txt and $WORKDIR/sirene_evt.txt"
diff \
    $WORKDIR/sirene_root.txt         \
    $WORKDIR/sirene_evt.txt