JThetaMCS.sh

Go to the documentation of this file.

#!/bin/zsh
script=${0##*/}
 
if [ -z $JPP_DIR ]; then
    echo "Variable JPP_DIR undefined."
    exit
fi
 
source $JPP_DIR/setenv.sh $JPP_DIR >& /dev/null
 
set_variable:  DEBUG       PHYSICS_DEBUG          2
set_variable:  WORKDIR     PHYSICS_WORKDIR        ${TMPDIR:-/tmp}
set_variable:  FORMAT      GRAPHICS_FORMAT        gif
set_variable+  BATCH       GRAPHICS_BATCH         -B
 
if do_usage $*; then
    usage "$script"
fi
 
$JPP_DIR/examples/JPhysics/JThetaMCS \
    -o $WORKDIR/mcs.root          \
    -d $DEBUG
 
JPlot2D \
    -f $WORKDIR/mcs.root:h2       \
    -z "1e-5 1e-1" -Z             \
    -XX                           \
    -> "E [GeV]"                 \
    -< "x [m]"                   \
    -T "#theta [rad]"             \
    -O COLZ                       \
    -o mcs.$FORMAT                $BATCH
 
for X in 1 2 3 4 5 6; do
 
    let "XMIN = $X - 0.005"
    let "XMAX = $X + 0.005"
 
    JPlot1D \
        -f mcs.root:h2                          \
        -PY                                     \
        -x "$XMIN $XMAX"                        \
        -z "1e-6 2e-2" -Z                       \
        -> "x [m]"                             \
        -\^ "#theta_{0}"                        \
        -T "E = 10^{$X} [GeV]"                  \
        -o ${TMPDIR:-/tmp}/mcs_${X}.$FORMAT     $BATCH
done
 
montage \
    -tile 3x2                                   \
    -geometry +0+0                              \
    ${TMPDIR:-/tmp}/mcs_{1..6}.$FORMAT          \
    $WORKDIR/MCS.$FORMAT                        >& /dev/null