JGizmo/JHistogram3D.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]) / ([4]*[4])) * exp(-0.5 * (y-[2])*(y-[2]) / ([4]*[4])) * exp(-0.5 * (z-[3])*(z-[3]) / ([4]*[4]))"
    
    $JPP_DIR/examples/JGizmo/JHistogram3D \
        -o $WORKDIR/histogram.root                                \
        -F "$FORMULA"                                             \
        -@ "p0 = 1.0; p1 = 0.0; p2 = 2.0; p3 = 4.0; p4 = 1.0"     \
        -x "100 -10 +10"                                          \
        -y "100 -10 +10"                                          \
        -z "100 -10 +10"                                          \
        -n 100000
    
    
    JFit3D \
        -f $WORKDIR/histogram.root:h0                   \
        -o $WORKDIR/fit.root                            \
        -F "$FORMULA"                                   \
        -@ "p0 = GetMaximum"                            \
        -@ "p1 = GetMean(1)"                            \
        -@ "p2 = GetMean(2)"                            \
        -@ "p3 = GetMean(3)"                            \
        -@ "p4 = GetRMS(1)"                             \
        -d $DEBUG --!
 
    echo -n "chi2/NDF: "
    JPrintFit \
        -f $WORKDIR/fit.root:h0                         \
        -F "GetChisquare/GetNDF"
 
 
    for PROJECTION in xy xz zy; do
 
        echo "Plotting ${PROJECTION}-projection"
    
        JPlot2D \
            -f $WORKDIR/fit.root:h0                         \
            -P $PROJECTION                                  \
            -O COLZ
    done
fi
 
 
if (( 1 )); then
 
    set_variable FORMULA      "[0] * pow(10,[1]*x)"
 
    $JPP_DIR/examples/JGizmo/JHistogram3D \
        -o $WORKDIR/histogram.root                      \
        -F "$FORMULA"                                   \
        -@ "p0 = 1.0; p1 = -2;"                         \
        -x "100  0  7"                                  \
        -y "100 -1 +1"                                  \
        -z "100 -1 +1"                                  \
        -n 10000000
 
 
    JFit3D \
        -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"
 
    for PROJECTION in xy xz; do
 
        echo "Plotting ${PROJECTION}-projection"        
        
        JPlot2D \
            -f $WORKDIR/fit.root:h0                     \
            -z "1e-1 1e5"                               \
            -YYZ                                        \
            -P $PROJECTION                              \
            -O COLZ
    done
fi
 
 
if (( 1 )); then
 
    set_variable FORMULA      "[0] * pow(10,[1]*y)"
 
    $JPP_DIR/examples/JGizmo/JHistogram3D \
        -o $WORKDIR/histogram.root                     \
        -F "$FORMULA"                                  \
        -@ "p0 = 1.0; p1 = -2"                         \
        -x "100 -1 +1"                                 \
        -y "100  0  7"                                 \
        -z "100 -1 +1"                                 \
        -n 10000000
 
 
    JFit3D \
        -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"
 
    for PROJECTION in yx yz; do
 
        echo "Plotting ${PROJECTION}-projection"
        
        JPlot2D \
            -f $WORKDIR/fit.root:h0                    \
            -z "1e-1 1e5"                              \
            -YYZ                                       \
            -P $PROJECTION                             \
            -O COLZ
    done
fi
 
 
if (( 1 )); then
 
    set_variable FORMULA      "[0] * pow(10,[1]*z)"
 
    $JPP_DIR/examples/JGizmo/JHistogram3D \
        -o $WORKDIR/histogram.root                     \
        -F "$FORMULA"                                  \
        -@ "p0 = 1.0; p1 = -2"                         \
        -x "100 -1 +1"                                 \
        -y "100 -1 +1"                                 \
        -z "100  0  7"                                 \
        -n 10000000
 
 
    JFit3D \
        -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"
 
    for PROJECTION in zx zy; do
 
        echo "Plotting ${PROJECTION}-projection"        
 
        JPlot2D \
            -f $WORKDIR/fit.root:h0                    \
            -z "1e-1 1e5"                              \
            -YYZ                                       \
            -P $PROJECTION                             \
            -O COLZ
    done
fi