Jpp
JPlotNPE1D.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 set_variable D 100
31 
32 set_variable YMAX 1e-3
33 set_variable YMIN 1e-8
34 
35 set_variable PDF $JPP_DATA/J%p.dat
36 
37 
38 if ( do_usage $* ); then
39  usage "$script [energy [distance]]"
40 fi
41 
42 case $# in
43  2) set_variable D $2;;
44  1) set_variable E $1;;
45 esac
46 
47 
48 for dir in NORTH EAST SOUTH WEST; do
49 
50  echo "Generating $dir."
51 
52  eval DIR=\$$dir
53 
54  $JPP_DIR/examples/JPhysics/JPlotNPE1D \
55  -f ${PDF/\%/13} \
56  -f ${PDF/\%/14} \
57  -D "${DIR}" \
58  -R "${D}" \
59  -x "1000000 -1.0 +1.0" \
60  -E "${E}" \
61  -o F\[${dir}\].root --!
62 done
63 
64 if (( 1 )); then
65 
66  JPlot1D \
67  -f F\[NORTH\].root:h0 \
68  -f F\[EAST\].root:h0 \
69  -f F\[SOUTH\].root:h0 \
70  -f F\[WEST\].root:h0 \
71  -y "$YMIN $YMAX" -Y \
72  -> "cos #theta_{0}" \
73  -\^ "P [npe]" \
74  -L TL \
75  -T "E = ${E} GeV#; D = ${D} m" \
76  -o NPE.gif
77 fi
78 
79 if (( 1 )); then
80 
81  JPlot1D \
82  -f F\[NORTH\].root:h0 \
83  -f F\[EAST\].root:h0 \
84  -f F\[SOUTH\].root:h0 \
85  -f F\[WEST\].root:h0 \
86  -x "0.7 0.8" \
87  -> "cos #theta_{0}" \
88  -\^ "P [npe]" \
89  -L TL \
90  -T "E = ${E} GeV#; D = ${D} m" \
91  -o NPE-zoom.gif
92 fi
Generated by   doxygen 1.8.16