JAcousticsEventBuilder.sh

Go to the documentation of this file.

#!/bin/zsh
script=${0##*/}
 
if [ -z $JPP_DIR ]; then
    echo "Variable JPP_DIR undefined."
    exit
fi

source $JPP_DIR/setenv.sh $JPP_DIR >& /dev/null

set_variable  DEBUG               2
set_variable  DIR                 $JPP_DIR/examples/JAcoustics/
set_variable  WORKDIR             ./
set_variable  QUANTILE            0.0
set_variable  TMAX_S              0.1
set_variable  NUMBER_OF_HITS      15

if ( do_usage $* ); then
    usage "$script <detector> <run>"
fi

case $# in
    2)  set_variable RUN          $argv[2];
        set_variable DETECTOR     $argv[1];;
    *)  fatal "Wrong number of arguments.";;
esac

JCookie.sh

# create input files

if [[ ! -f $WORKDIR/detector.detx ]]; then

    JDetectorDB \
        -D $DETECTOR                          \
        -r $RUN                               \
        -o $WORKDIR/detector.detx --!
fi

if [[ ! -f $WORKDIR/tripod.txt ]]; then

    cat>$WORKDIR/tripod.txt<<EOF
# identifier UTM position (x,y,z)
#
  12    587888 4017079 -3450.0
  14    587436 4016992 -3450.0
  16    587843 4016727 -3450.0
EOF
fi 

if [[ ! -f $WORKDIR/base.txt ]]; then

    $DIR/JBaseModuleWriter \
        -a $WORKDIR/detector.detx                 \
        -o $WORKDIR/base.txt
fi

if (( 1 )); then

    JDrawDetector2D \
        -a $WORKDIR/detector.detx                 \
        -T $WORKDIR/tripod.txt                    \
        -S 0.8                                    \
        -L
fi

# create intemediate ROOT file for fast access

if [[ ! -f $WORKDIR/toashort.root ]]; then

    JConvertDB \
        -q toashort                               \
        -@ "detid  = `getDetector -D $DETECTOR`"  \
        -@ "minrun = $RUN"                        \
        -@ "maxrun = $RUN"                        \
        -o $WORKDIR/toashort.root                 \
        -d  $DEBUG --!
fi

JAcousticsEventBuilder \
    -a $WORKDIR/detector.detx                 \
    -f $WORKDIR/toashort.root                 \
    -@ "Q            = $QUANTILE;"            \
    -@ "TMax_s       = $TMAX_S;"              \
    -@ "numberOfHits = $NUMBER_OF_HITS;"      \
    -T $WORKDIR/tripod.txt                    \
    -B $WORKDIR/base.txt                      \
    -o $WORKDIR/event.root                    \
    -d  $DEBUG --!

JPrintAcousticsEvent \
    -f $WORKDIR/event.root                    \
    -n 3