JAcousticsMonitor_short.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       ACOUSTICS_DEBUG        2
set_variable:  WORKDIR     ACOUSTICS_WORKDIR      ./
set_variable:  FORMAT      GRAPHICS_FORMAT        gif
set_variable   LOGBOOK     \"Individual Logbooks\"
set_variable   AUTHOR      \"Gatius C\"
set_variable   FACET       "ELcode"

if do_usage $*; then
    usage "$script <detector file> minrun maxrun report"\
    "\nIn case of failures, if report = 1 an elog entry is created, if report =! 1 no elog entry is created."
fi

if (( $# == 4 )); then
    set_variable  DETECTOR        $argv[1]
    set_variable  MINRUN          $argv[2]
    set_variable  MAXRUN          $argv[3]
    set_variable  REPORT          $argv[4]
else
    fatal "Wrong number of arguments."
fi

set_variable  USER            km3net
set_variable  PASSWORD        pyrosoma

JCookie.sh 

eval `JPrintDetector -a $DETECTOR -O IDENTIFIER`
eval `JPrintDetector -a $DETECTOR -O SUMMARY`

JAcoustics.sh   $DETECTOR_ID

source JAcousticsToolkit.sh

CHECK_EXIT_CODE

for ((RUN=$MINRUN; RUN<=$MAXRUN; RUN++)); do

    set_variable  INPUT_FILE        $WORKDIR/KM3NeT_${(l:8::0::0:)DETECTOR_ID}_${(l:8::0::0:)RUN}_toashort.root
    set_variable  OUTPUT_FILE       $WORKDIR/KM3NeT_${(l:8::0::0:)DETECTOR_ID}_${(l:8::0::0:)RUN}_event.root
    set_variable  MONITOR_FILE      $WORKDIR/KM3NeT_${(l:8::0::0:)DETECTOR_ID}_${(l:8::0::0:)RUN}_monitor.root
    set_variable  SUMMARY_FILE      $WORKDIR/summary_${(l:8::0::0:)RUN}.txt
    set_variable  TEST_FILE         $WORKDIR/KM3NeT_${(l:8::0::0:)DETECTOR_ID}_${(l:8::0::0:)RUN}_test.root
    set_variable  PARAMETERS_FILE   $JPP_DIR/examples/JAcoustics/acoustics_monitor_${(l:8::0::0:)DETECTOR_ID}.txt

    if [[ ! -f $INPUT_FILE ]]; then

        JConvertDB \
        -q toashort                                  \
        -@ "detid  = `getDetector -D $DETECTOR_ID`"  \
        -@ "minrun = $RUN"                           \
        -@ "maxrun = $RUN"                           \
        -o $INPUT_FILE                               \
        -d $DEBUG --!                                
    fi

    if [[ ! -f $OUTPUT_FILE ]]; then

        JAcousticsEventBuilder \
        -a $DETECTOR                                 \
        -f $INPUT_FILE                               \
        -T $WORKDIR/tripod.txt                       \
        -V $WORKDIR/sound_velocity.txt               \
        -W $WORKDIR/waveform.txt                     \
        -@ $WORKDIR/acoustics_trigger_parameters.txt \
        -H $WORKDIR/hydrophone.txt                   \
        -o $OUTPUT_FILE                              \
        -d $DEBUG --!                                
    fi

    set_array time_range $(JRuns -D "$DETECTOR_ID" -@ "run = ${RUN}" -F UNIXJOBSTART -F UNIXJOBEND --)
   
    set_variable lifetime_s $(( ($time_range[2]-$time_range[1])/1000 ))
    
    #if lifetime is unknown, assume 6hr
    if [ -z "$lifetime_s" ]; then
        lifetime_s=21600
    
    #avoid 0 sec lifetime
    elif (($lifetime_s == 0)); then
        lifetime_s=1
    fi
    
    $JPP_DIR/examples/JAcoustics/JAcousticsMonitor_short \
    -a $DETECTOR                                         \
    -f $OUTPUT_FILE                                      \
    -l $lifetime_s                                       \
    -o $MONITOR_FILE                                     \
    -d $DEBUG                                            

    read -r NUMBER_OF_FAILURES <<< "$($JPP_DIR/examples/JAcoustics/JAcousticsMonitorTest  \
    -f $MONITOR_FILE                                    \
    -a $DETECTOR                                        \
    -P $PARAMETERS_FILE                                 \
    -F $FACET                                           \
    -w $SUMMARY_FILE                                    \
    -t $TEST_FILE                                       \
    -r $RUN                                             \
    -d "1")"

    typeset -A EMITTERS

    get_tripods+     $WORKDIR/tripod.txt       EMITTERS
    get_transmitters+ $WORKDIR/transmitter.txt EMITTERS

    ATTACHMENTS=""

    for EMITTER in ${(@k)EMITTERS}; do

        JPlot2D \
        -f "${TEST_FILE}:H\[${EMITTER}\].rate-test"        \
        -> "string"                                       \
        -< "floor"                                        \
        -\^ "Rate outside range"                           \
        -z "0 1"                                           \
        -O COLZ                                            \
        -T "[${EMITTER}]"                                  \
        -o $WORKDIR/rate-test\[${EMITTER}\]_${RUN}.$FORMAT \
        -B

        ATTACHMENTS+=" -f ${WORKDIR}/rate-test\[${EMITTER}\]_${RUN}.${FORMAT}"
    done

    #report to e-log 
    if [[ ${NUMBER_OF_FAILURES[-1]} -gt 0 && $REPORT == 1 ]]; then

        eval elog                \
        -h elog.km3net.de        \
        -p 3000                  \
        -u $USER $PASSWORD       \
        -l $LOGBOOK              \
        -a Subject=\"Warning, detector ID ${DETECTOR_ID}\"  \
        -a Type=\"Acoustics monitoring\"         \
        -a Author=$AUTHOR        \
        -n 0 -v                  \
        -m $SUMMARY_FILE         \
        $ATTACHMENTS
    fi

    CHECK_EXIT_CODE

    #remove temporary files
    rm -f $SUMMARY_FILE
    rm -f $WORKDIR/rate-test*.$FORMAT
    
done