JShadowPDF.sh

Go to the documentation of this file.

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

# ------------------------------------------------------------------------------------------
#
#                         Utility script to plot PDFs.
#
# ------------------------------------------------------------------------------------------


if [ -z $JPP_DIR ]; then
    echo "Variable JPP_DIR undefined."
    exit
fi


source $JPP_DIR/setenv.sh $JPP_DIR


set_variable  NORTH  "0.00 0.00"
set_variable  EAST   "1.57 0.00"
set_variable  SOUTH  "3.14 0.00"
set_variable  WEST   "1.57 3.14"

set_variable  E      1e2
set_variable  R      10
set_variable  YMAX   2e-0
set_variable  YMIN   1e-7

set_variable  HISTOGRAM "100 -5 +50"

if ( do_usage $* ); then
    usage "$script [energy [distance]]"
fi

case $# in
    2) R=$2;&
    1) E=$1;;
esac


for dir in NORTH EAST SOUTH WEST; do

    echo "Generating $dir."

    eval D=\$$dir

    JDrawPDF \
        -F1 -F3 -F5                    \
        -D "${D}" -R $R -E $E          \
        -H "${HISTOGRAM}"              \
        -o fa\[${dir}\].root --!
    
    JDrawPDF \
        -@ "MODULE_RADIUS_M=0.00;"     \
        -F2 -F4 -F6                    \
        -D "${D}" -R $R -E $E          \
        -H "${HISTOGRAM}"              \
        -o fb\[${dir}\].root --!

    JDrawPDF \
        -@ "MODULE_RADIUS_M=0.25;"     \
        -F2 -F4 -F6                    \
        -D "${D}" -R $R -E $E          \
        -H "${HISTOGRAM}"              \
        -o fc\[${dir}\].root --!
        
    JPlot1D \
        -f fa\[${dir}\].root:h0        \
        -f fb\[${dir}\].root:h0        \
        -f fc\[${dir}\].root:h0        \
        -y "$YMIN $YMAX" -Y            \
        -\> "#Deltat [ns]"             \
        -\^ "dP/dt [npe/ns]"           \
        -L TR                          \
        -T "E = ${E} GeV#; R = ${R} m" \
        -o shadow-${R}-${dir}.gif -B

done