JGizmo/JHistogram1D.sh

Go to the documentation of this file.

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

# ------------------------------------------------------------------------------------------
#
#                         Utility script to test JHistogram1D, JFit and JPlot1D.
#
# ------------------------------------------------------------------------------------------


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      /tmp/

if ( do_usage $* ); then
    usage "$script [working directory]"
fi

case $# in
    1) WORKDIR=$1;;
esac


if (( 1 )); then

    set_variable FORMULA      "[0] * exp(-0.5 * (x-[1])*(x-[1]) / ([2]*[2]))"
    set_variable OUTPUT_FILE  histogram.root
    
    JHistogram1D \
        -o $WORKDIR/$OUTPUT_FILE               \
        -F "$FORMULA"                          \
        -="p0 = 1.0; p1 = 0.0; p2 = 2.0"       \
        -x "100 -10 +10"                       \
        -n 10000
    
    
    JFit \
        -f $WORKDIR/${OUTPUT_FILE}:h0          \
        -o $WORKDIR/fit.root                   \
        -F "$FORMULA"                          \
        -@"p0 = GetMaximum"                    \
        -@"p1 = GetMean"                       \
        -@"p2 = GetRMS"                        \
        -d $DEBUG
    
    
    JPlot1D \
        -f $WORKDIR/fit.root:h0                

fi


if (( 1 )); then

    FORMULA="[0] * pow(10,[1]*x)"
    
    OUTPUT_FILE=histogram.root
    
    JHistogram1D \
        -o $WORKDIR/$OUTPUT_FILE               \
        -F "$FORMULA"                          \
        -="p0 = 1.0; p1 = -2"                  \
        -x "16 3 7"                            \
        -n 10000000
    
    
    JFit \
        -f $WORKDIR/${OUTPUT_FILE}:h0          \
        -o $WORKDIR/fit.root                   \
        -F "$FORMULA"                          \
        -@"p0 = GetMaximum * 1.0e6"            \
        -@"p1 = -2"                            \
        -d3
    
    
    JPlot1D \
        -f $WORKDIR/fit.root:h0                \
        -y "1e-1 1e8"                          \
        -XXY                                   \
        -G XY                                  
    
fi