JAcoustics/JToAshortMonitor.sh

Go to the documentation of this file.

#!/usr/bin/env zsh
script=${0##*/}
 
if [ -z $JPP_DIR ]; then
    echo "Variable JPP_DIR undefined."
    exit
fi
 
source $JPP_DIR/setenv.sh $JPP_DIR >& /dev/null
 
zmodload zsh/mathfunc
 
set_variable:  DEBUG       ACOUSTICS_DEBUG        2
set_variable:  WORKDIR     ACOUSTICS_WORKDIR      ./
set_variable:  XMIN        TOA_XMIN               -2.0e-4
set_variable:  XMAX        TOA_XMAX               +2.0e-4
set_variable:  FORMAT      GRAPHICS_FORMAT        gif
set_variable+  BATCH       GRAPHICS_BATCH         -B
set_variable:  OPTION      GRAPHICS_OPTION
set_variable   RANGE_S     0 60
 
if do_usage $*; then
    usage "$script <detector file> (input file)+"\
          "\nInput files correspond to the output of JAcousticsEventBuilder[.sh]."
fi
 
if (( $# < 2 )); then
    fatal "Wrong number of arguments."
fi
 
set_variable  DETECTOR            $argv[1]
set_array     INPUT_FILES         $argv[2,-1]
 
eval `JPrintDetector -a $DETECTOR -O SUMMARY`
eval `JPrintDetector -a $DETECTOR -O IDENTIFIER`
 
JAcoustics.sh   $DETECTOR_ID
 
source JAcousticsToolkit.sh
 
CHECK_EXIT_CODE
 
typeset -A EMITTERS
 
get_tripods+      $WORKDIR/tripod.txt       EMITTERS
get_transmitters+ $WORKDIR/transmitter.txt  EMITTERS
 
typeset -i NX
typeset -i NY
 
let "NY = sqrt(${#STRINGS}) * 0.75 + 0.25"
let "NX = (${#STRINGS} + $NY - 1)/ $NY"
 
for EMITTER in -1 ${(@k)EMITTERS}; do
 
    $JPP_DIR/examples/JAcoustics/JToAMonitor \
        -a $DETECTOR                                              \
        -f "$INPUT_FILES[*]"                                      \
        -o $WORKDIR/toa.root                                      \
        -E $EMITTER                                               \
        -x "$((int(2.5e6*($XMAX - XMIN)))) $XMIN $XMAX"           \
        -T "$RANGE_S"                                             \
        -d $DEBUG
 
    echo -n "Creating graphics for emitter $EMITTER.."
 
    JPlot1D \
        -f "$WORKDIR/toa.root:\[%\]\.t1"                          \
        -> "#Deltat [s]"                                         \
        -\^ "number of events [a.u.]"                             \
        -s 111110                                                 \
        $OPTION                                                   \
        -N "X 505"                                                \
        -T "[$EMITTER]"                                           \
        -o $WORKDIR/T_${EMITTER}s.$FORMAT                         $BATCH
 
    for STRING in ${STRINGS[*]}; do
 
        JPlot1D \
            -f "$WORKDIR/toa.root:\[\ *${STRING}\ \ *.*\]\.t1"        \
            -> "#Deltat [s]"                                         \
            -\^ "number of events [a.u.]"                             \
            -L TL                                                     \
            $OPTION                                                   \
            -N "X 505"                                                \
            -T "[$EMITTER]"                                           \
            -o $WORKDIR/T_${EMITTER}_${(l:4::0::0:)STRING}.$FORMAT    $BATCH
    done
 
    montage \
        -tile ${NX}x${NY}                        \
        -geometry +0+0                           \
        $WORKDIR/T_${EMITTER}_*.$FORMAT          \
        $WORKDIR/T_${EMITTER}.$FORMAT            >& /dev/null
 
    echo " -> $WORKDIR/T_${EMITTER}.$FORMAT"
 
    rm -f $WORKDIR/toa.root
    rm -f $WORKDIR/T_${EMITTER}_*.$FORMAT
 
done