JEMShower.sh

Go to the documentation of this file.

#!/bin/zsh
#
# \author mdejong
#
version=1.0
script=${0##*/}
 
# ------------------------------------------------------------------------------------------
#
#             Utility script to plot number of photons - GeV for EM-showers
#
# ------------------------------------------------------------------------------------------
 
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     ./  
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
 
set_variable    GRAPH_TXT    $WORKDIR/graph.txt
set_variable    GRAPH_ROOT   $WORKDIR/graph\[Geant4\].root
set_variable    GEANC_ROOT   $WORKDIR/geanc\[JSirene\].root
 
# number of photons per GeV
# from Geant4 simulation by Daniel Lopez
 
cat>$GRAPH_TXT<<EOF
  1.0e-3    90.96
  2.0e-3   277.36
  3.0e-3   485.82
  4.0e-3   692.83
  5.0e-3   890.01
  6.0e-3  1098.53
  7.0e-3  1285.47
  8.0e-3  1502.86
  9.0e-3  1687.15
 10.0e-3  1856.97
100.0e-3  1.905e4
  1.0     1.889e5
 10.0     1.875e6
100.0     1.881e7
EOF
 
JGraph \
    -f $GRAPH_TXT   \
    -o $GRAPH_ROOT  
 
JEMShower \
    -o $GEANC_ROOT
 
JPlot1D \
    -f ${GRAPH_ROOT}:\.\*         \
    -f ${GEANC_ROOT}:h0           \
    -x "-4 +4"                    \
    -XX                           \
    -y "1e1 1e9"                  \
    -Y                            \
    -> "E [GeV]"                 \
    -\^ "number of photons"       \
    -T ""                         \
    -o $WORKDIR/geant4.$FORMAT    $BATCH
 
 
JPlot1D \
    -f ${GEANC_ROOT}:h1          \
    -XX                          \
    -y "0.0 1.1"                 \
    -> "E [GeV]"                \
    -\^ "Geant4/JSirene"         \
    -T ""                        \
    -o $WORKDIR/ratio.$FORMAT    $BATCH