billabong-random:plot.sh

Go to the documentation of this file.

#!/usr/bin/env zsh
script=${0:t}
 
source $JPP_DIR/setenv.sh $JPP_DIR >& /dev/null
 
zmodload zsh/mathfunc
 
set_variable   DEBUG           2
set_variable:  FORMAT          GRAPHICS_FORMAT        gif
set_variable+  BATCH           GRAPHICS_BATCH         -B
 
if do_usage $*; then
    usage "$script <input file>" 
fi
 
if (( $# != 1 )); then
    fatal "Wrong number of arguments."
fi
 
set_variable  GRAPH_TXT    $1
set_variable  GRAPH_ROOT   graph.root
 
JGraph \
    -f $GRAPH_TXT   \
    -o $GRAPH_ROOT
 
let "XMIN = $(awk '{ print $1 }' $GRAPH_TXT | sort -g | head -1) - 0.05"
let "XMAX = $(awk '{ print $1 }' $GRAPH_TXT | sort -g | tail -1) + 0.50"
 
if (( 1 )) && [[ $GRAPH_TXT == *noise* ]]; then
    
    JFit \
        -f ${GRAPH_ROOT}:\.\*              \
        -F "[0] * pow(x,[1])"              \
        -@ "p0 =  1.0e-6"                  \
        -@ "p1 = 10.0"                     \
        -o fit.root
 
    JPlot1D \
        -@ "leftMargin  = 0.2"             \
        -@ "titleOffset = 1.3"             \
        -f fit.root:\.\*                   \
        -> "noise rate [Hz]"              \
        -\^ "event rate [Hz]"              \
        -x "$XMIN $XMAX"  -X               \
        -y "1.0e-6 1.0e0" -Y               \
        -T ""                              \
        -o noise.$FORMAT                   $BATCH
fi
 
if (( 0 )) && [[ $GRAPH_TXT == *noise* ]]; then
 
    JPlot1D \
        -@ "leftMargin  = 0.2"             \
        -@ "titleOffset = 1.3"             \
        -f ${GRAPH_ROOT}:\.\*              \
        -> "noise rate [Hz]"              \
        -\^ "event rate [Hz]"              \
        -x "$XMIN $XMAX"  -X               \
        -y "1.0e-6 1.0e0" -Y               \
        -T ""                              \
        -o noise.$FORMAT                   $BATCH
fi
 
if [[ $GRAPH_TXT == *cpu* ]]; then
 
    JFit \
        -f ${GRAPH_ROOT}:\.\*              \
        -F "[0] * pow(x,[1])"              \
        -@ "p0 = 2.0"                      \
        -@ "p1 = 3.5"                      \
        -o fit.root
 
    JPlot1D \
        -@ "leftMargin  = 0.2"             \
        -@ "titleOffset = 1.3"             \
        -f fit.root:\.\*                   \
        -> "noise rate [Hz]"              \
        -\^ "cpu"                          \
        -x "$XMIN $XMAX"  -X               \
        -y "1.00 5.0e3" -Y                 \
        -G XY                              \
        -T ""                              \
        -o cpu.$FORMAT                     $BATCH
fi
 
rm -f $GRAPH_ROOT >& /dev/null
rm -f fit.root    >& /dev/null