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