JRadiation.sh

Go to the documentation of this file.

#!/usr/bin/env zsh
#
# \author mdejong
#
version=1.0
script=${0##*/}
 
# ------------------------------------------------------------------------------------------
#
#        Utility script to plot radiation cross sections, shower energy, range and b(E)
#
# ------------------------------------------------------------------------------------------
 
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}/  
set_variable:  FORMAT      GRAPHICS_FORMAT        gif
set_variable+  BATCH       GRAPHICS_BATCH         -B
 
if do_usage $*; then
    usage "$script [working directory]"
fi
 
case $# in
    1) set_variable  WORKDIR  $1;;
esac
 
if (( 1 )); then
 
    $JPP_DIR/examples/JPhysics/JRadiation \
        -o $WORKDIR/lambda.root                \
        -O length                              \
        -d $DEBUG
 
    for TYPE in eerad Brems gnrad; do
        
        JPlot1D \
            -f "$WORKDIR/lambda.root:${TYPE}"  \
            -y "1e-1 1e6"                      \
            -XXY                               \
            -\^ "\lambda [m]"                  \
            -> "E [GeV]"                      \
            -L TR -T "interaction length"      \
            -O "HIST]["                        \
            -G Y                               \
            -o lamdba-$TYPE.$FORMAT            $BATCH
    done
 
fi
 
 
if (( 1 )); then
 
    $JPP_DIR/examples/JPhysics/JRadiation \
        -o $WORKDIR/radiation.root             \
        -O energy                              \
        -n 40000                               \
        -d $DEBUG
 
    for TYPE in eerad Brems gnrad; do
        
        JPlot1D \
            -f "$WORKDIR/radiation.root:$TYPE" \
            -y "0.1 1e9"                       \
            -XXY                               \
            -\^ "<E_{s}> [GeV]"                \
            -> "E [GeV]"                      \
            -L TL -T "shower energy"           \
            -G Y                               \
            -o Es-$TYPE.$FORMAT                $BATCH
    done
 
fi
 
 
if (( 1 )); then
 
    $JPP_DIR/examples/JPhysics/JRadiation \
        -o $WORKDIR/radiation.root             \
        -O range                               \
        -n 1000                                \
        -d $DEBUG
 
    JPlot1D \
        -f $WORKDIR/radiation.root:R\.\*       \
        -y "0 40"                              \
        -XX                                    \
        -> "E [GeV]"                          \
        -\^ "R [km]"                           \
        -L TL -T "muon range"                  \
        -o range.$FORMAT                       $BATCH
 
fi
 
 
if (( 1 )); then
 
    $JPP_DIR/examples/JPhysics/JRadiation \
        -o $WORKDIR/radiation.root             \
        -O eloss                               \
        -n 4000000                             \
        -d $DEBUG
 
    JPlot1D   \
        -f $WORKDIR/radiation.root:hb          \
        -y "0 0.5"                             \
        -XX                                    \
        -> "E [GeV]"                          \
        -\^ "b [km^{-1}]"                      \
        -T "Energy loss"                       \
        -o b.$FORMAT                           $BATCH
 
fi
 
rm -f $WORKDIR/lambda.root                     >& /dev/null
rm -f $WORKDIR/radiation.root                  >& /dev/null