JF1.sh

Go to the documentation of this file.

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

# ------------------------------------------------------------------------------------------
#
#                         Utility script to test JF1.
#
# ------------------------------------------------------------------------------------------


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]))"
    
    JF1 \
        -o $WORKDIR/f1.root                    \
        -F "$FORMULA"                          \
        -x "-5.0 +5.0"                         \
        -@ "p0 = 1.0"                          \
        -@ "p1 = 0.0"                          \
        -@ "p2 = 1.0"                          \
        -d $DEBUG
    
    JPlot1D \
        -f $WORKDIR/f1.root:user              

fi


if (( 1 )); then

    set_variable FORMULA      "[0] * pow(10,[1]*x)"
    
    JF1 \
        -o $WORKDIR/f1.root                    \
        -F "$FORMULA"                          \
        -x "0.0 7.0"                           \
        -@ "p0 =  1.0"                         \
        -@ "p1 = -2.0"                         \
        -d $DEBUG
    
    JPlot1D \
        -f $WORKDIR/f1.root:user               \
        -y "1e-14 1.0"                         \
        -XXY                                   \
        -G XY                                  
    
fi