examples/JTrigger/JTriggerEfficiencyRunByRun.sh

Go to the documentation of this file.

#!/bin/zsh
#
# \author mdejong
#
version=1.0
script=${0##*/}

# ------------------------------------------------------------------------------------------
#
#                Utility script to test JTriggerEfficiency in run-by-run mode.
#
# ------------------------------------------------------------------------------------------


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            ${TMPDIR:-/tmp}/                                  
set_variable   INPUT_FILE         $JPP_DATA/sirene.root                  
set_variable   DETECTOR           $JPP_DATA/km3net_reference.detx        
set_variable   TRIGGER            $JPP_DATA/trigger_parameters_arca.txt  
set_variable   BACKGROUND_HZ      5.0e3                                  
set_variable   FRAME_DURATION_NS  100.0e6                                
set_variable   NUMBER_OF_PMTS     31                                     


if do_usage $*; then
    usage "$script [input file [working directory]]"
fi

case $# in
    2) set_variable  WORKDIR     $2;&
    1) set_variable  INPUT_FILE  $1;;
esac

set_variable  DETECTOR_OUT_OF_SYNC                     $WORKDIR/detector_out_of_sync.detx
set_variable  DAQ_FILE_OUT_OF_SYNC                     $WORKDIR/trigger_efficiency_out_of_sync.root
set_variable  OUTPUT_FILE                              $WORKDIR/trigger_efficiency.root
set_variable  TRIGGEREFFICIENCY_TURBOT                 $WORKDIR/turbot.detx
set_variable  TRIGGEREFFICIENCY_TRIGGERED_EVENTS_ONLY  Y


function clean()
{
    rm -f $DETECTOR_OUT_OF_SYNC
    rm -f $DAQ_FILE_OUT_OF_SYNC
    rm -f $OUTPUT_FILE
    rm -f $TRIGGEREFFICIENCY_TURBOT
}


clean

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

set_variable  STRING         ${STRINGS[1]}
set_variable  MODULE         `getModule -a $DETECTOR -L "$STRING 10"`


# distort detector file

JEditDetector \
    -a $DETECTOR                            \
    -o $DETECTOR_OUT_OF_SYNC                \
    -M "$MODULE add $FRAME_DURATION_NS"

# generate fake raw data file, with distorted detector file.
# make sure to write summary data for JTriggerEfficiency and
# L2 data for JTurbot (inside JTriggerEfficiencyRunByRun.sh).

JTriggerEfficiency \
    -a ${DETECTOR_OUT_OF_SYNC}        \
    -f ${INPUT_FILE}                  \
    -o ${DAQ_FILE_OUT_OF_SYNC}        \
    -@ ${TRIGGER}                     \
    -@ "writeL0         = 0;"         \
    -@ "writeL1         = 0;"         \
    -@ "writeL2         = 1;"         \
    -@ "writeSN         = 0;"         \
    -@ "writeSummary    = 1;"         \
    -d ${DEBUG}                       \
    -B "$BACKGROUND_HZ"               \
    --!

rm -f $OUTPUT_FILE

$JPP_BIN/JTriggerEfficiencyRunByRun.sh \
    $DETECTOR                         \
    $DETECTOR                         \
    $INPUT_FILE                       \
    $OUTPUT_FILE                      \
    $DAQ_FILE_OUT_OF_SYNC             

for (( PMT = 0; $PMT != $NUMBER_OF_PMTS; ++PMT )); do

    set_variable  RESULT  `getPMTStatus -a $TRIGGEREFFICIENCY_TURBOT -P "$MODULE $PMT"  | sed -n 's/.*OUT_OF_SYNC *\([0-1]\)/\1/p'`

    echo "PMT ($MODULE,$PMT) has OUT_OF_SYNC $RESULT"

    if (( $RESULT != 1 )); then
        fatal "PMT ($MODULE,$PMT) has OUT_OF_SYNC $RESULT != 1"
    fi
done

set_array  RESULT  `$JPP_DIR/examples/JTrigger/JDAQHit -f $OUTPUT_FILE | grep "^ *$MODULE "`

printf "Module %10d: %8d %8d\n" $MODULE  $RESULT[2]   $RESULT[3] 

if (( ${#RESULT} >= 3 )); then

    if (( $RESULT[2] >= 20 )); then
        fatal "Module $MODULE has too many trigger hits $RESULT[2]"
    fi

    if (( $RESULT[3] <= 20 )); then
        fatal "Module $MODULE has too few snapshot hits $RESULT[3]"
    fi
else

    fatal "Module $MODULE has no data."
fi

clean

---

Generated by   1.8.5