JCompareSummaryslice.sh

Go to the documentation of this file.

#!/usr/bin/env zsh
#
# \author mdejong
#
script=${0##*/}
 
source  $JPP_DIR/setenv.sh  $JPP_DIR  >&  /dev/null
 
set_variable   DEBUG       2
set_variable:  FORMAT      GRAPHICS_FORMAT        gif
set_variable+  BATCH       GRAPHICS_BATCH         -B
 
if do_usage $*; then
    usage "$script <input file> <input file>"
fi
 
if (( $# == 2 )); then
    INPUT_FILES=($argv[1,2])
else
    fatal "Wrong number of arguments."
fi
 
$JPP_DIR/examples/JSummaryslice/JCompareSummaryslice \
    -f "$INPUT_FILES[1,2]" \
    -R "-0.1 +0.1"         \
    -o summary.root        \
    -d $DEBUG --!
 
JPlot1D \
    -f summary.root:h0               \
    -> "#DeltaR [kHz]"              \
    -\^ "number of events [a.u.]"    \
    -T ""                            \
    -s 110                           \
    -o h0.$FORMAT                    $BATCH
    
JPlot2D \
    -f summary.root:ha               \
    -> "R [kHz]"                    \
    -< "#DeltaR [kHz]"              \
    -X                               \
    -Z                               \
    -O colz                          \
    -T ""                            \
    -o ha.$FORMAT                    $BATCH
    
JPlot2D \
    -f summary.root:hb               \
    -> "R [kHz]"                    \
    -< "#Deltabit"                  \
    -y "-10 +10"                     \
    -X                               \
    -Z                               \
    -O colz                          \
    -T ""                            \
    -o hb.$FORMAT                    $BATCH
 
typeset -A H2
 
H2[h1]="R [kHz]"
H2[h2]="high-rate veto"
H2[h3]="FIFO (almost) full"
H2[h4]="UDP maximum"
H2[h5]="UDP count"
H2[h6]="UDP trailer"
 
JPrintRange2D -f summary.root:h1 | read -A BUFFER
 
JLine -p "$BUFFER[1] $BUFFER[2] $BUFFER[4] $BUFFER[5]" -@ "width = 1" -o line.root
 
for KEY VALUE in ${(kv)H2}; do
 
    if [[ "$VALUE" == *"Hz"* ]]; then
        OPTIONS=(-X -Y)
    else
        OPTIONS=()
 
        if (( $(JPrintResult -f summary.root:$KEY -F GetNbinsX) <= 2 )); then
            OPTIONS+=("-N X 2")
        fi
        if (( $(JPrintResult -f summary.root:$KEY -F GetNbinsY) <= 2 )); then
            OPTIONS+=("-N Y 2")
        fi
    fi
 
    if [[ "$KEY" == "h1" ]]; then
        OPTIONS+=("-f line.root:TLine")
    fi
 
    JPlot2D \
        -f summary.root:$KEY         \
        -> "$VALUE"                 \
        -< "$VALUE"                 \
        $OPTIONS[*]                  \
        -Z                           \
        -O colz                      \
        -T ""                        \
        -o $KEY.$FORMAT              $BATCH
done