toashort_monitor.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:  FORMAT      GRAPHICS_FORMAT        gif
set_variable+  BATCH       GRAPHICS_BATCH         -B
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
 
set_variable HISTOGRAM    "4000 -1.0e-1 +1.0e-1"
set_variable FORMULA      "[0] * exp(-0.5 * (x-[1])*(x-[1]) / ([2]*[2])) + [3]"
 
typeset -A BUFFER
 
for EMITTER in -1 ${(@k)EMITTERS}; do
 
    echo Processing $EMITTER...
 
    $JPP_DIR/examples/JAcoustics/JToAMonitor \
        -a $DETECTOR                                               \
        -f "$INPUT_FILES[*]"                                       \
        -o $WORKDIR/toa.root                                       \
        -E $EMITTER                                                \
        -T "$RANGE_S"                                              \
        -x "$HISTOGRAM"                                            \
        -d 1
 
    JFit \
        -f "$WORKDIR/toa.root:\[%\]\.t1"                           \
        -o  ${TMPDIR:-/tmp}/fit.root                               \
        -F "$FORMULA"                                              \
        -x "-4.0e-4 +4.0e-4"                                       \
        -@ "p0 = GetMaximum * 0.8"                                 \
        -= "p1 =  0.0"                                             \
        -@ "p2 = 30.0e-6"                                          \
        -@ "p3 = 10.0"                                             \
        -d $DEBUG
    
    P0=($(JPrintFit -f ${TMPDIR:-/tmp}/fit.root:\.\* -@ p0))
    P3=($(JPrintFit -f ${TMPDIR:-/tmp}/fit.root:\.\* -@ p3))
 
    if (( $? == 0 )); then
        BUFFER[$EMITTER]=$(($P0[1] / $P3[1]))
    fi
 
    JPlot1D \
        -f  ${TMPDIR:-/tmp}/fit.root:\.\*                          \
        -x "-1.0e-3 +1.0e-3"                                       \
        -Y                                                         \
        -> "#Deltat [s]"                                          \
        -\^ "number of events [a.u.]"                              \
        -N "X 505"                                                 \
        -s 1000000000                                              \
        -T "[$EMITTER]"                                            \
        -o $WORKDIR/T_${EMITTER}.$FORMAT                           $BATCH
done
 
for KEY VALUE in ${(kv)BUFFER}; do
    printf "Emitter %2d  S/N %7.1f\n" $KEY $VALUE
done