JAcoustics/JToAshortMonitor.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

zmodload zsh/mathfunc

set_variable:  DEBUG       ACOUSTICS_DEBUG        2
set_variable:  WORKDIR     ACOUSTICS_WORKDIR      ./
set_variable:  XMIN        TOASHORT_XMIN         -2.0e-4
set_variable:  XMAX        TOASHORT_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/JToAshortMonitor \
        -a $DETECTOR                                              \
        -f "$INPUT_FILES[*]"                                      \
        -o $WORKDIR/toashort.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/toashort.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/toashort.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/toashort.root
    rm -f $WORKDIR/T_${EMITTER}_*.$FORMAT

done

---

Generated by   1.8.5