JF2.sh

Go to the documentation of this file.

#!/usr/bin/env zsh
#
# \author mdejong
#
version=1.0
script=${0##*/}
 
# ------------------------------------------------------------------------------------------
#
#                         Utility script to test JF2.
#
# ------------------------------------------------------------------------------------------
 
 
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]))"
    
    JF2 \
        -o $WORKDIR/f2.root                    \
        -F "$FORMULA"                          \
        -@ "p0 = 1.0; p1 = 0.0; p2 = 1.0"      \
        -x "-5 +5"                             \
        -y "-5 +5"                             \
        -N 1000
    
    JPlot2D \
        -f $WORKDIR/f2.root:user               \
        -z "1e-10 2" -Z                        \
        -O COLZ
 
fi
 
 
if (( 1 )); then
 
    set_variable FORMULA      "[0] * pow(x,[1])"
 
    JF2 \
        -o $WORKDIR/f2.root                    \
        -F "$FORMULA"                          \
        -@ "p0 = 1.0; p1 = -2.0"               \
        -x "1e0 1e7"                           \
        -y "-1 +1"                             
 
    JPlot2D \
        -f $WORKDIR/f2.root:user               \
        -z "1e-14 1.0"                         \
        -XZ                                    \
        -O COLZ
 
fi
 
 
if (( 1 )); then
 
    set_variable FORMULA      "[0] * pow(10,[1]*x)"
 
    JF2 \
        -o $WORKDIR/f2.root                    \
        -F "$FORMULA"                          \
        -@ "p0 = 1.0; p1 = -2.0"               \
        -x " 0  7"                             \
        -y "-1 +1"                             
 
    JPlot2D \
        -f $WORKDIR/f2.root:user               \
        -z "1e-14 1.0"                         \
        -Z                                     \
        -O COLZ
 
    JPlot2D \
        -f $WORKDIR/f2.root:user               \
        -z "1e-14 1.0"                         \
        -XXZ                                   \
        -O COLZ
 
fi
 
 
if (( 1 )); then
 
    set_variable FORMULA      "[0] * pow(y,[1])"
 
    JF2 \
        -o $WORKDIR/f2.root                    \
        -F "$FORMULA"                          \
        -@ "p0 = 1.0; p1 = -2.0"               \
        -x "-1 +1"                             \
        -y "1e0 1e7"                           
 
    JPlot2D \
        -f $WORKDIR/f2.root:user               \
        -z "1e-14 1.0"                         \
        -YZ                                    \
        -O COLZ
 
fi
 
 
if (( 1 )); then
 
    set_variable FORMULA      "[0] * pow(10,[1]*y)"
 
    JF2 \
        -o $WORKDIR/f2.root                    \
        -F "$FORMULA"                          \
        -@ "p0 = 1.0; p1 = -2.0"               \
        -x "-1 +1"                             \
        -y " 0  7"                             
 
    JPlot2D \
        -f $WORKDIR/f2.root:user               \
        -z "1e-14 1.0"                         \
        -Z                                     \
        -O COLZ
 
    JPlot2D \
        -f $WORKDIR/f2.root:user               \
        -z "1e-14 1.0"                         \
        -YYZ                                   \
        -O COLZ
 
fi