JTriggerProfiler.sh

Go to the documentation of this file.

#!/usr/bin/env zsh
#
# \author mdejong
#
version=1.0
script=${0##*/}
 
# ------------------------------------------------------------------------------------------
#
#                  Utility script to plot JTriggerProcessor CPU requirements.
#
# ------------------------------------------------------------------------------------------
 
 
if [ -z $JPP_DIR ]; then
    echo "Variable JPP_DIR undefined."
    exit
fi
 
 
source $JPP_DIR/setenv.sh $JPP_DIR
 
 
set_variable   DEBUG             1                                      
set_variable   WORKDIR           ./                                     
set_variable   FRAMETIME_MS      100                                    
set_variable:  RATE_L1_HZ        PROFILER_RATE_L1_HZ    500                                    
set_variable:  LIVETIME          PROFILER_LIVETIME      1
set_variable   RECYCLING         5 100e3
set_variable:  FORMAT            GRAPHICS_FORMAT        gif
set_variable+  BATCH             GRAPHICS_BATCH         -B
 
if do_usage $*; then
    usage "$script <detector file> <trigger file> [working directory]"
fi
 
case $# in
    3) set_variable  WORKDIR   $3;&
    2) set_variable  TRIGGER   $2;
       set_variable  DETECTOR  $1;;
    *) fatal "Wrong number of arguments."
esac
 
 
set_variable   PIPE         $WORKDIR/pipe.dat
set_variable   CPU_TXT      $WORKDIR/cpu.txt
set_variable   CPU_ROOT     $WORKDIR/cpu.root
set_variable   RATE_TXT     $WORKDIR/rate.txt
set_variable   RATE_ROOT    $WORKDIR/rate.root
 
 
typeset -A SETUPS  # PMT rate -> number of slices
 
SETUPS[5e3]=$(($LIVETIME * 100))
SETUPS[7e3]=$(($LIVETIME * 100))
SETUPS[9e3]=$(($LIVETIME * 100))
SETUPS[11e3]=$(($LIVETIME * 100))
SETUPS[13e3]=$(($LIVETIME * 50))
SETUPS[15e3]=$(($LIVETIME * 50))
SETUPS[17e3]=$(($LIVETIME * 50))
SETUPS[19e3]=$(($LIVETIME * 20))
 
if (( 1 )); then
 
    echo "Detector: $DETECTOR"
    echo "Trigger parameters:"
    cat  $TRIGGER
 
    for RATE_L0_HZ NUMBER_OF_SLICES in ${(kv)SETUPS}; do
 
        LOGGER_FILE=$WORKDIR/trigger_${RATE_L0_HZ}.log
 
        if [[ ! -f $LOGGER_FILE ]]; then
            
            rm -f  $PIPE
            mkfifo $PIPE
            
            echo "Generating random background at ${RATE_L0_HZ} Hz."
            
            nohup JRandomTimesliceWriter \
                  -a $DETECTOR                   \
                  -o $PIPE                       \
                  -B "$RATE_L0_HZ $RATE_L1_HZ"   \
                  -n $NUMBER_OF_SLICES           \
                  -N "$RECYCLING"                \
                  -d $DEBUG                      > /dev/null 2>&1 &
            
            echo Processing random data.
            
            JTriggerProcessor \
                -a $DETECTOR                     \
                -f $PIPE                         \
                -o "/dev/null"                   \
                -@ $TRIGGER                      \
                -C JDAQTimesliceL0               \
                -P                               \
                -d $DEBUG                        > $LOGGER_FILE 2>&1
 
            rm -f  $PIPE
        fi
    done
fi
 
 
rm -f $CPU_TXT  $RATE_TXT   >&  /dev/null
rm -f $CPU_ROOT $RATE_ROOT  >&  /dev/null
 
for RATE_L0_HZ in `echo ${(k)SETUPS} | tr ' ' '\n' | sort -g`; do
    
    LOGGER_FILE=$WORKDIR/trigger_${RATE_L0_HZ}.log
 
    if [[ -f $LOGGER_FILE ]]; then
 
        let "TOTAL_CPU = 0.0"
    
        for KEY in "Calibration" "L0" "L1" "L2" "Timeslice router" "Trigger" "Trigger router"; do
        
            ELAPSED_TIME=`awk -v key="$KEY" '$0 == key,/CPU/ { if ($3 == "elapsed") print $1 }' $LOGGER_FILE`
        
            let "TOTAL_CPU += $ELAPSED_TIME / $FRAMETIME_MS"
        done
 
        TRIGGER_RATE_HZ=`awk '/Trigger rate \[Hz\]/ { print $4 }' $LOGGER_FILE`
 
        echo "$(($RATE_L0_HZ * 1e-3)) $TOTAL_CPU"       >> $CPU_TXT
        echo "$(($RATE_L0_HZ * 1e-3)) $TRIGGER_RATE_HZ" >> $RATE_TXT
    fi
done
 
 
JGraph \
    -f $CPU_TXT               \
    -o $CPU_ROOT
 
JPlot1D \
    -f $CPU_ROOT:\.\*         \
    -x "0 20"                 \
    -y "0 60"                 \
    -> "Rate [kHz]"          \
    -\^ "CPU [cores]"         \
    -O CP                     \
    -T ""                     \
    -o cpu.$FORMAT            $BATCH
 
 
JGraph \
    -f $RATE_TXT              \
    -o $RATE_ROOT
 
JPlot1D \
    -f $RATE_ROOT:\.\*        \
    -x "0 20"                 \
    -y "0 400"                \
    -> "Rate [kHz]"          \
    -\^ "Trigger rate [Hz]"   \
    -O CP                     \
    -T ""                     \
    -o rate.$FORMAT           $BATCH