JPlotNPE-PDG.sh

Go to the documentation of this file.

#!/usr/bin/env zsh
#
# \author lquinn
#
version=1.0
script=${0##*/}
 
# ------------------------------------------------------------------------------------------
#
#                         Utility script to plot integrals of tabulated PDFs.
#
# ------------------------------------------------------------------------------------------
 
 
if [ -z $JPP_DIR ]; then
    echo "Variable JPP_DIR undefined."
    exit
fi
 
source $JPP_DIR/setenv.sh $JPP_DIR
 
zmodload zsh/mathfunc
 
set_variable   WORKDIR     ${TMPDIR:-/tmp}                        
set_variable   DEBUG       2                           
set_variable   TYPE        14                          
set_variable   E           1.0                         
 
set_variable:  YMAX        PHYSICS_YMAX                1e+2
set_variable:  YMIN        PHYSICS_YMIN                1e-5
set_variable:  PDF         PHYSICS_PDF                 $JPP_DATA/J%p.dat.gz
set_variable:  CDF         PHYSICS_CDF                 $JPP_DATA/I%p.dat.gz
 
set_variable:  FORMAT      GRAPHICS_FORMAT        gif
set_variable+  BATCH       GRAPHICS_BATCH         -B
 
if do_usage $*; then
    usage "$script [energy]"
fi
 
case $# in
    1) set_variable E  $1;;
    0) ;;
    *) fatal "Wrong number of arguments."
esac
 
DIRECTIONS=(NORTH EAST SOUTH WEST)
 
for DIR in $DIRECTIONS[*]; do
 
    echo "Generating $DIR."
 
    $JPP_DIR/examples/JPhysics/JPlotNPE-PDG \
        -P ${PDF/\%/${TYPE}}                   \
        -C ${CDF/\%/${TYPE}}                   \
        -o $WORKDIR/npe\[${DIR}\].root         \
        -D "${(P)DIR}"                         \
            -E "${E}"                              \
            -d $DEBUG
 
        for FUNCTION in pdf npe cdf; do
 
        JPlot2D \
                    -f $WORKDIR/npe\[${DIR}\].root:${FUNCTION}      \
                    -z "$YMIN $YMAX" -Z                             \
                    -T "${DIR} E = ${E} GeV"                        \
                    -O "colz"                                       \
                    -> "D [m]"                                     \
                    -< "cos #theta_{0}"                            \
                    -o ${FUNCTION}-${DIR}-${E}GeV.$FORMAT           $BATCH
        done
 
        set_variable  YMIN   0.90
        set_variable  YMAX   1.10
 
        for FUNCTION in npe cdf; do
 
        JOpera2D \
                    -f $WORKDIR/npe\[${DIR}\].root:${FUNCTION}        \
                    -f $WORKDIR/npe\[${DIR}\].root:pdf                \
            -u Divide                                         \
                    -o $WORKDIR/ratio-${FUNCTION}-${DIR}.root
        
        JPlot2D \
                    -f $WORKDIR/ratio-${FUNCTION}-${DIR}.root:Divide  \
                    -z "$YMIN $YMAX"                                  \
                    -T "${FUNCTION}/pdf ${DIR}"                       \
                    -> "D [m]"                                       \
                    -< "cos #theta_{0}"                              \
                    -O COLZ                                           \
                    -o ratio-${FUNCTION}-${DIR}.$FORMAT               $BATCH
        done
done
 
rm -f $WORKDIR/npe\[${^DIRECTIONS}\].root
rm -f $WORKDIR/ratio-npe-${^DIRECTIONS}.root
rm -f $WORKDIR/ratio-cdf-${^DIRECTIONS}.root