Jpp
JPlotNPE.sh
Go to the documentation of this file.
1 #!/bin/zsh
2 #
3 # \author lquinn
4 #
5 version=1.0
6 script=${0##*/}
7 
8 # ------------------------------------------------------------------------------------------
9 #
10 # Utility script to plot tabulated PDFs.
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 NORTH "0.00 0.00"
25 set_variable EAST "1.57 0.00"
26 set_variable SOUTH "3.14 0.00"
27 set_variable WEST "1.57 3.14"
28 
29 set_variable E 0
30 
31 set_variable YMAX 2e+2
32 set_variable YMIN 1e-9
33 
34 set_variable PDF $JPP_DATA/J%p.dat
35 
36 
37 if ( do_usage $* ); then
38  usage "$script [energy]"
39 fi
40 
41 case $# in
42  1) set_variable E $1;;
43 esac
44 
45 
46 for dir in NORTH EAST SOUTH WEST; do
47 
48  echo "Generating $dir."
49 
50  eval D=\$$dir
51 
52  if (( 1 )); then
53 
54  JPlotNPE \
55  -f ${PDF/\%/13} \
56  -f ${PDF/\%/14} \
57  -D "${D}" \
58  -E "${E}" \
59  -o F\[${dir}\].root --!
60 
61  fi
62 
63  if (( 1 )); then
64 
65  JPlot2D \
66  -f F\[${dir}\].root:h0 \
67  -z "$YMIN $YMAX" -Z \
68  -T "${dir} E = ${E} GeV" \
69  -O "colz" \
70  -> "D [m]" \
71  -< "cos #theta_{0}" \
72  -o NPE-${dir}-${E}GeV.gif
73  fi
74 done
Generated by   doxygen 1.8.16