software/JCompass/JAHRSMonitor.sh

Go to the documentation of this file.

#!/bin/zsh
#
#
# \author mdejong
#
version=1.0
script=${0##*/}
 
if [ -z $JPP_DIR ]; then
    echo "Variable JPP_DIR undefined."
    exit
fi
 
source $JPP_DIR/setenv.sh $JPP_DIR
 
set_variable:  DEBUG        COMPASS_DEBUG          2
set_variable:  WORKDIR      COMPASS_WORKDIR        ./
set_variable:  FORMAT       GRAPHICS_FORMAT        gif
set_variable+  BATCH        GRAPHICS_BATCH         -B
 
if do_usage $*; then
    usage "$script <detector file> <run>"
fi
 
if (( $# != 2 )); then
    fatal "Wrong number of parameters."
fi
 
set_variable     DETECTOR                $argv[1]
set_variable     RUN                     $argv[2]
 
set_variable     AHRS_CALIBRATION_FILE   $WORKDIR/ahrs_calibration.txt
 
eval `JPrintDetector -a $DETECTOR -O IDENTIFIER`
 
set_variable     DETOID                  `getDetector -D $DETECTOR_ID -O string`
 
typeset -a  INPUT_FILES
 
if [[ ! "$RUN" == *"-"* ]]; then
    set_variable  INPUT_FILE   $WORKDIR/KM3NeT_${(l:8::0::0:)DETECTOR_ID}_${(l:8::0::0:)RUN}_ahrs.root
else
    set_variable  INPUT_FILE   $WORKDIR/KM3NeT_${(l:8::0::0:)DETECTOR_ID}_${(l:8::0::0:)RUN%%-*}-${(l:8::0::0:)RUN##*-}_ahrs.root
fi
 
if [[ ! -f $INPUT_FILE ]]; then
 
    echo -n "Writing $INPUT_FILE... "
    
    JConvertDB \
        -q "ahrs"                      \
        -@ "detid  = $DETECTOR_ID"     \
        -@ "minrun = ${RUN%%-*}"       \
        -@ "maxrun = ${RUN##*-}"       \
        -o $INPUT_FILE                 \
        -d $DEBUG --
 
    if (( $? == 0 )); then
        echo "OK"
    else
        fatal "No AHRS data"
    fi
fi
 
if [[ ! -f $AHRS_CALIBRATION_FILE ]]; then
 
    JAHRSCalibration \
        -D $DETECTOR_ID                \
        -o $AHRS_CALIBRATION_FILE      \
        -d $DEBUG
fi
 
JAHRSMonitor \
    -a  $DETECTOR                                \
    -c  $AHRS_CALIBRATION_FILE                   \
    -f  $INPUT_FILE                              \
    -o  $WORKDIR/ahrs.root                       \
    -q "0.1 0.2 0.3 0.4 0.5"                     \
    -Q  1                                        \
    -d  $DEBUG --!
 
 
JPlot1D \
    -f $WORKDIR/ahrs.root:h0                    \
    -Y                                          \
    -> "|A|"                                   \
    -\^ "number of entries"                     \
    -s 110010                                   \
    -T  "$DETOID"                               \
    -o A.$FORMAT                                $BATCH
 
JPlot1D \
    -f $WORKDIR/ahrs.root:h1                    \
    -Y                                          \
    -> "|H|"                                   \
    -\^ "number of entries"                     \
    -s 110010                                   \
    -T  "$DETOID"                               \
    -o H.$FORMAT                                $BATCH
 
JPlot2D \
    -f $WORKDIR/ahrs.root:h2                    \
    -> "|A|"                                   \
    -< "|H|"                                   \
    -O colz                                     \
    -T  "$DETOID"                               \
    -o AH.$FORMAT                               $BATCH
 
JPlot1D \
    -f $WORKDIR/ahrs.root:hn                    \
    -Y                                          \
    -> "number of entries"                     \
    -\^ "number of AHRSs"                       \
    -s 110010                                   \
    -T  "$DETOID"                               \
    -o ahrs.$FORMAT                             $BATCH
 
typeset -A BUFFER
 
BUFFER[ha]=all
BUFFER[hb]=valid
BUFFER[hc]=calibration
BUFFER[hd]=unique
 
for H2 LABEL in ${(kv)BUFFER}; do
 
    JPlot2D \
        -f $WORKDIR/ahrs.root:$H2                   \
        -> "string"                                \
        -< "floor"                                 \
        -O COLZ                                     \
        -T "${DETOID}/${LABEL}"                     \
        -o $H2.$FORMAT                              $BATCH
done
 
montage \
    -tile 4x1                                   \
    -geometry +0+0                              \
    ha.$FORMAT                                  \
    hb.$FORMAT                                  \
    hc.$FORMAT                                  \
    hd.$FORMAT                                  \
    H2.$FORMAT                                  >& /dev/null