JGizmo/JHistogram2D.sh

Go to the documentation of this file.

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

# ------------------------------------------------------------------------------------------
#
#                         Utility script to test JHistogram2D, JFit2D and JPlot2D.
#
# ------------------------------------------------------------------------------------------


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

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

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


if (( 1 )); then

    set_variable FORMULA      "[0] * exp(-0.5 * (x-[1])*(x-[1]) / ([2]*[2])) * exp(-0.5 * (y-[1])*(y-[1]) / ([2]*[2]))"
    
    $JPP_DIR/examples/JGizmo/JHistogram2D \
        -o $WORKDIR/histogram.root             \
        -F "$FORMULA"                          \
        -@ "p0 = 1.0; p1 = 0.0; p2 = 2.0"      \
        -x "100 -10 +10"                       \
        -y "100 -10 +10"                       \
        -n 100000

    JPrintMaximum2D -f $WORKDIR/histogram.root:h0
    
    JFit2D \
        -f $WORKDIR/histogram.root:h0          \
        -o $WORKDIR/fit.root                   \
        -F "$FORMULA"                          \
        -@ "p0 = GetMaximum"                   \
        -@ "p1 = GetMean"                      \
        -@ "p2 = GetRMS"                       \
        -d $DEBUG

    echo -n "chi2/NDF: "
    JPrintFit \
        -f $WORKDIR/fit.root:h0                \
        -F "GetChisquare/GetNDF"
    
    JPlot2D \
        -f $WORKDIR/fit.root:h0                \
        -O COLZ

fi


if (( 1 )); then

    set_variable FORMULA      "[0] * pow(10,[1]*x)"

    $JPP_DIR/examples/JGizmo/JHistogram2D \
        -o $WORKDIR/histogram.root             \
        -F "$FORMULA"                          \
        -@ "p0 = 1.0; p1 = -2"                 \
        -x "100  0  7"                         \
        -y "100 -1 +1"                         \
        -n 10000000


    JFit2D \
        -f $WORKDIR/histogram.root:h0          \
        -o $WORKDIR/fit.root                   \
        -F "$FORMULA"                          \
        -@ "p0 = GetMaximum"                   \
        -@ "p1 = -2"                           \
        -d3

    echo -n "chi2/NDF: "
    JPrintFit \
        -f $WORKDIR/fit.root:h0                \
        -F "GetChisquare/GetNDF"

    JPlot2D \
        -f $WORKDIR/fit.root:h0                \
        -z "1e-1 1e5"                          \
        -XXZ                                   \
        -O COLZ

fi


if (( 1 )); then

    set_variable FORMULA      "[0] * pow(10,[1]*y)"

    $JPP_DIR/examples/JGizmo/JHistogram2D \
        -o $WORKDIR/histogram.root             \
        -F "$FORMULA"                          \
        -@ "p0 = 1.0; p1 = -2"                 \
        -x "100 -1 +1"                         \
        -y "100  0  7"                         \
        -n 10000000


    JFit2D \
        -f $WORKDIR/histogram.root:h0          \
        -o $WORKDIR/fit.root                   \
        -F "$FORMULA"                          \
        -@ "p0 = GetMaximum"                   \
        -@ "p1 = -2"                           \
        -d3

    echo -n "chi2/NDF: "
    JPrintFit \
        -f $WORKDIR/fit.root:h0                \
        -F "GetChisquare/GetNDF"

    JPlot2D \
        -f $WORKDIR/fit.root:h0                \
        -z "1e-1 1e5"                          \
        -YYZ                                   \
        -O COLZ

fi