examples/JDB/JAHRSMonitor.sh

Go to the documentation of this file.

#!/usr/bin/env 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       2                       
set_variable   WORKDIR     `pwd`                   
set_variable   DIR         $JPP_DIR/examples/JDB/  
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]
 
eval `JPrintDetector -a $DETECTOR -O IDENTIFIER`
eval `JPrintDetector -a $DETECTOR -O SUMMARY`
 
set_variable     AHRS_CALIBRATION_FILE   $WORKDIR/ahrs_calibration.txt
set_variable     INPUT_FILE              $WORKDIR/KM3NeT_${(l:8::0::0:)DETECTOR_ID}_${(l:8::0::0:)RUN}_ahrs.root
 
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 --
 
    echo "OK"
fi
 
if [[ ! -f $AHRS_CALIBRATION_FILE ]]; then
 
    JAHRSCalibration \
        -D $DETECTOR_ID                \
        -o $AHRS_CALIBRATION_FILE      \
        -d $DEBUG
fi
 
$DIR/JAHRSMonitor \
    -f $INPUT_FILE                   \
    -o $WORKDIR/ahrs_old.root
 
$DIR/JAHRSMonitor \
    -f $INPUT_FILE                   \
    -A $AHRS_CALIBRATION_FILE        \
    -o $WORKDIR/ahrs_new.root
 
 
 
typeset -Z 4 STRING
typeset -Z 2 FLOOR
 
typeset -A X_LABEL
 
X_LABEL[A]="A"
X_LABEL[B]="H"
 
for key value in "${(@kv)X_LABEL}"; do
 
    for STRING in $STRINGS[*]; do
 
        echo -n "Creating $key graphics for string $STRING.."
 
        set_variable  FIRST_FLOOR  1
 
        for (( FLOOR = $FIRST_FLOOR; $FLOOR <= $LAST_FLOOR; FLOOR += 1 )); do
 
            eval MODULE=\${STRING_${(l:4::0::0:)STRING}\[$(($FLOOR))\]:--1}
 
            JPlot1D \
                -f "$WORKDIR/ahrs_old.root:${MODULE}.${key}"          \
                -f "$WORKDIR/ahrs_new.root:${MODULE}.${key}"          \
                -> "${value}"                                        \
                -\^ "number of events [a.u.]"                         \
                -s 110110                                             \
                -T "(${STRING},${FLOOR})"                             \
                -o ${TMPDIR:-/tmp}/${key}_${STRING}_${FLOOR}.$FORMAT             $BATCH
        done
 
        montage \
            -tile 6x3                            \
            -geometry +0+0                       \
            ${TMPDIR:-/tmp}/${key}_${STRING}_*.$FORMAT      \
            $WORKDIR/${key}_${STRING}.$FORMAT    >&  /dev/null
 
        echo " -> $WORKDIR/${key}_${STRING}.$FORMAT"
 
    done
done