JOperationsRunQuality.sh

Go to the documentation of this file.

#!/bin/zsh
#
# \author mdejong
#
# Auxiliary script to produce QA/QC data for operations (near-real-time) monitoring.
#
# The data is stored int he <qaqc file> and can be viewed with the script JOperationsRunMonitor.sh.
version=1.0
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       QAQC_DEBUG           0
set_variable:  WORKDIR     QAQC_WORKDIR         ${TMPDIR:-/tmp}/.$$
 
if do_usage $*; then
    usage "$script <detector file> (run file(s))+ <qaqc file> [working directory]"
fi
 
case $# in
    4)  set_variable WORKDIR     $4;&
    3)  set_variable DETECTOR    $1;
        set_array    INPUT_FILES $2;
        set_variable QAQC_FILE   $3;;
    *)  fatal "Wrong number of arguments."
esac
 
TRIGGER_PARAMETERS=(
    trigger3DMuon.enabled
    trigger3DShower.enabled
    triggerMXShower.enabled
    triggerNB.enabled
    writeL0.prescale
    writeL1.prescale
    writeL2.prescale
    writeSN.prescale)
 
TRIGGERS=(`JPrintTriggerBits`)
 
DAQ=(
    JDAQTimeslice
    JDAQTimesliceL0
    JDAQTimesliceL1
    JDAQTimesliceL2
    JDAQTimesliceSN
    JDAQSummaryslice
    JDAQEvent)
 
typeset -a PARAMETERS
 
PARAMETERS=(
    GIT detector run livetime_s UTCMin_s UTCMax_s
    ${TRIGGER_PARAMETERS[*]}
    ${DAQ[*]}
    JTriggerReprocessor
    ${TRIGGERS[*]}
    in-sync out-sync DAQ WR HRV FIFO
    PMTs MEAN_Rate_Hz RMS_Rate_Hz
)
 
mkdir -p $WORKDIR
 
if [[ ! -s $QAQC_FILE ]]; then
    echo $PARAMETERS[*]                                                  >  $QAQC_FILE
fi
 
for INPUT_FILE in ${INPUT_FILES[*]}; do
 
    # check presence of trigger parameters
    
    JPrintTriggerParameters -f ${INPUT_FILE} -d 0
    CHECK_EXIT_CODE
 
    # determine values
    
    JTriggerMonitor -v |& awk '/version/ { print $2 }'            | read     GIT_VERSION
    CHECK_EXIT_CODE
 
    JPrintDAQHeader -f ${INPUT_FILE} -d 0                         | read -A  DAQ_HEADER
    CHECK_EXIT_CODE
 
    
 
    set_variable  DETECTOR_ID    ${DAQ_HEADER[1]}
    set_variable  RUN            ${DAQ_HEADER[2]}
    
    JDetectorDB -D ${DETECTOR_ID} -r ${RUN} -o ${DETECTOR}               # install detector file as-is for JTriggerReprocessor
    CHECK_EXIT_CODE
 
 
 
 
    getLivetime     -f ${INPUT_FILE} -d 0                         | read -A  LIVETIME_S
    CHECK_EXIT_CODE
 
    getUTC          -f ${INPUT_FILE} -d 0                         | read -A  UTC_S
    CHECK_EXIT_CODE
 
    { JTriggerMonitor -f ${INPUT_FILE} -d ${DEBUG} -Q 3 3>&1 1>&4            | read -A  TRIGGER_MONITOR } 4>&1
    CHECK_EXIT_CODE
 
    set_variable  OUTPUT_FILE    ${WORKDIR}/KM3NeT_${(l:8::0::0:)DETECTOR_ID}_${(l:8::0::0:)RUN}_trigger.root
 
    JTriggerReprocessor \
        -a ${DETECTOR}                          \
        -f ${INPUT_FILE}                        \
        -o ${OUTPUT_FILE}                       \
        -U                                      \
        -C -\.\*                                \
        -d ${DEBUG}                                  
    CHECK_EXIT_CODE
 
    JPrintTree -f ${OUTPUT_FILE} -@ "type = JDAQEvent" -k number_of_entries  | read     TRIGGER_REPROCESSOR
    
    { JTurbot \
        -f ${INPUT_FILE}                \
        -a ${DETECTOR}                  \
        -o /dev/null                    \
        -d ${DEBUG}                     \
        -Q 3 3>&1 1>&4                                                   | read -A  TURBOT } 4>&1
    CHECK_EXIT_CODE
 
    { JSummaryMonitor -f ${INPUT_FILE} -d ${DEBUG} -Q 3 3>&1 1>&4        | read -A  SUMMARY_MONITOR } 4>&1
    CHECK_EXIT_CODE
 
    set_variable  RUN            $DAQ_HEADER[2]
    
    # print values
 
    printf " %32s"   ${GIT_VERSION}                                      >> $QAQC_FILE
    printf " %8d"    ${DAQ_HEADER[1,2]}                                  >> $QAQC_FILE
    printf " %9.1f"  ${LIVETIME_S[1]}                                    >> $QAQC_FILE
    printf " %12.1f" ${UTC_S[*]}                                         >> $QAQC_FILE
    
    for KEY in ${TRIGGER_PARAMETERS[*]}; do
 
        JPrintTriggerParameters -f ${INPUT_FILE} -k $KEY -d 0                   | read VALUE
 
        printf " %4d" ${VALUE-0}                                         >> $QAQC_FILE
    done
 
    for KEY in ${DAQ[*]}; do
        
        JPrintTree -f ${INPUT_FILE} -@ "type = $KEY" -k number_of_entries       | read VALUE
 
        printf " %8d" ${VALUE-0}                                         >> $QAQC_FILE
    done
    
    printf " %8d"   ${TRIGGER_REPROCESSOR:-0}                            >> $QAQC_FILE
    printf " %8d"   ${TRIGGER_MONITOR[*]}                                >> $QAQC_FILE
    printf " %4d"   ${TURBOT[1,2]}                                       >> $QAQC_FILE
    printf " %7.5f" ${SUMMARY_MONITOR[1,4]}                              >> $QAQC_FILE
    printf " %8.1f" ${SUMMARY_MONITOR[5]}                                >> $QAQC_FILE
    printf " %7.0f" ${SUMMARY_MONITOR[6,7]}                              >> $QAQC_FILE
    printf "\n"                                                          >> $QAQC_FILE
done
 
rm -rf $WORKDIR