JMultiPMT.sh

Go to the documentation of this file.

#!/usr/bin/env zsh
#
# \author mdejong
#
version=1.0
script=${0##*/}
 
# ------------------------------------------------------------------------------------------
#
#                         Example script to study L0 and L0 hit probabilities.
#
# ------------------------------------------------------------------------------------------
 
 
if [ -z $JPP_DIR ]; then
    echo "Variable JPP_DIR undefined."
    exit
fi
 
 
source $JPP_DIR/setenv.sh $JPP_DIR
 
set_variable   WORKDIR     ${TMPDIR:-/tmp}   
set_variable   DEBUG       3
set_variable   E_GEV       1e1        
set_variable   CD          0.5
set_variable   PDF         $JPP_DATA/J%p.dat.gz
 
set_variable:  PARTICLE    PHYSICS_PARTICLE       muon   # shower
 
set_variable:  FORMAT      GRAPHICS_FORMAT        gif
set_variable+  BATCH       GRAPHICS_BATCH         -B
 
if do_usage $*; then
    usage "$script <E [GeV]> [PDF]"\
          "The particle can be \"muon\" or \"shower\"."
fi
 
case $# in
    2) set_variable  PDF    $argv[2];&
    1) set_variable  E_GEV  $argv[1];;
    *) fatal "Wrong number of arguments."
esac
 
DIRECTIONS=(NORTH EAST SOUTH WEST)
 
for DIR in $DIRECTIONS[*]; do
 
        OUTPUT_FILE=F\[${DIR}\].root
 
        if [[ ! -f $OUTPUT_FILE ]]; then
 
            echo "Generating histograms for ${PARTICLE} direction ${DIR}."
 
                $JPP_DIR/examples/JPhysics/JMultiPMT \
                    -F $PDF                              \
                    -o $OUTPUT_FILE                      \
                    -O KM3NeT                            \
                    -E $E_GEV                            \
                    -c $CD                               \
                    -D "${(P)DIR}"                       \
            -P $PARTICLE                         \
                    -d $DEBUG
        fi
done
 
JPlot1D \
        -fF\[${^DIRECTIONS}\].root:L00                       \
    -G XY                                                \
    -Y                                                   \
        -> "R [m]"                                          \
        -\^ "P"                                              \
    -L TR                                                \
        -T "$DETECTOR $TYPE $PARTICLE"                       \
        -o L00.$FORMAT                                       $BATCH
 
for TYPE in L0 L1; do
 
        JPlot1D \
                -fF\[${^DIRECTIONS}\].root:${TYPE}${PARTICLE[1]}     \
                -> "R [m]"                                          \
                -\^ "P"                                              \
        -L TR                                                \
                -T "$DETECTOR $TYPE $PARTICLE"                       \
                -o F.${TYPE}.$FORMAT                                 $BATCH
done