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)+"
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/JToAshortMonitor \
        -a $DETECTOR                                               \
        -f "$INPUT_FILES[*]"                                       \
        -o $WORKDIR/toashort.root                                  \
        -E $EMITTER                                                \
        -T "$RANGE_S"                                              \
        -x "$HISTOGRAM"                                            \
        -d 1

    JFit \
        -f "$WORKDIR/toashort.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