Jpp
JGeanz.sh
Go to the documentation of this file.
1 #!/bin/zsh
2 #
3 # \author mdejong
4 #
5 version=1.0
6 script=${0##*/}
7 
8 # ------------------------------------------------------------------------------------------
9 #
10 # Utility script to plot longitudinal profile from EM showers.
11 #
12 # ------------------------------------------------------------------------------------------
13 
14 
15 if [ -z $JPP_DIR ]; then
16  echo "Variable JPP_DIR undefined."
17  exit
18 fi
19 
20 
21 source $JPP_DIR/setenv.sh $JPP_DIR
22 
23 
24 set_variable DEBUG 2
25 set_variable WORKDIR /tmp/
26 
27 if ( do_usage $* ); then
28  usage "$script [working directory]"
29 fi
30 
31 case $# in
32  1) WORKDIR=$1;;
33 esac
34 
35 
36 typeset -a INPUT_FILE
37 
38 for E_GEV in 1e0 1e1 1e2 1e3 1e4 1e5 1e6; do
39 
40  set_variable OUTPUT_FILE $WORKDIR/geanz\[${E_GEV}\ GeV\].root
41 
42  JGeanz \
43  -E "$E_GEV" \
44  -o $OUTPUT_FILE
45 
46  INPUT_FILE+=("$OUTPUT_FILE:Probability")
47  INPUT_FILE+=("$OUTPUT_FILE:Maximum")
48 
49 done
50 
51 
52 JPlot1D \
53  -f "$INPUT_FILE[*]" \
54  -g 2 \
55  -x "0 15" \
56  -> "z [m]" \
57  -\^ "dP/dz" \
58  -T "" -L TR -o geanz.gif