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      ./

if ( do_usage $* ); then
    usage "$script [working directory]"
fi

case $# in
    1) 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.gif


JPlot1D \
    -f ${GEANC_ROOT}:h1     \
    -XX                     \
    -y "0.0 1.1"            \
    -> "E [GeV]"           \
    -\^ "Geant4/JSirene"    \
    -T ""                   \
    -o $WORKDIR/ratio.gif