Jpp test-rotations-old
the software that should make you happy
Loading...
Searching...
No Matches
JPlotNPE-PDG.sh
Go to the documentation of this file.
1#!/usr/bin/env zsh
2#
3# \author lquinn
4#
5version=1.0
6script=${0##*/}
7
8# ------------------------------------------------------------------------------------------
9#
10# Utility script to plot integrals of tabulated PDFs.
11#
12# ------------------------------------------------------------------------------------------
13
14
15if [ -z $JPP_DIR ]; then
16 echo "Variable JPP_DIR undefined."
17 exit
18fi
19
20
21source $JPP_DIR/setenv.sh $JPP_DIR
22
23
24set_variable DIR $JPP_DIR/examples/JPhysics
25
26set_variable NORTH "0.00 0.00"
27set_variable EAST "1.57 0.00"
28set_variable SOUTH "3.14 0.00"
29set_variable WEST "1.57 3.14"
30
31set_variable WORKDIR ${TMPDIR:-/tmp}
32set_variable DEBUG 2
33set_variable TYPE 14
34set_variable E 1.0
35set_variable YMAX 1e+2
36set_variable YMIN 1e-5
37
38set_variable PDF $JPP_DATA/J%p.dat
39set_variable CDF $JPP_DATA/I%p.dat
40
41set_variable: FORMAT GRAPHICS_FORMAT gif
42set_variable+ BATCH GRAPHICS_BATCH -B
43
44if do_usage $*; then
45 usage "$script [energy]"
46fi
47
48case $# in
49 1) set_variable E $1;;
50esac
51
52
53for dir in NORTH EAST SOUTH WEST; do
54
55 echo "Generating $dir."
56
57 eval D=\$$dir
58
59 if (( 1 )); then
60
61 $DIR/JPlotNPE-PDG \
62 -P ${PDF/\%/${TYPE}} \
63 -C ${CDF/\%/${TYPE}} \
64 -D "${D}" \
65 -E "${E}" \
66 -o $WORKDIR/npe\[${dir}\].root \
67 -d $DEBUG --!
68
69 fi
70
71 if (( 0 )); then
72
73 for FUNCTION in pdf npe cdf; do
74
75 JPlot2D \
76 -f $WORKDIR/npe\[${dir}\].root:${FUNCTION} \
77 -z "$YMIN $YMAX" -Z \
78 -T "${dir} E = ${E} GeV" \
79 -O "colz" \
80 -> "D [m]" \
81 -< "cos #theta_{0}" \
82 -o ${FUNCTION}-${dir}-${E}GeV.$FORMAT $BATCH
83 done
84 fi
85
86 if (( 1 )); then
87
88 set_variable YMIN 0.90
89 set_variable YMAX 1.10
90
91 for FUNCTION in npe cdf; do
92
93 JOpera2D \
94 -f $WORKDIR/npe\[${dir}\].root:${FUNCTION} \
95 -f $WORKDIR/npe\[${dir}\].root:pdf \
96 -u Divide \
97 -o $WORKDIR/ratio-${FUNCTION}-${dir}.root
98
99 JPlot2D \
100 -f $WORKDIR/ratio-${FUNCTION}-${dir}.root:Divide \
101 -z "$YMIN $YMAX" \
102 -T "${FUNCTION}/pdf ${dir}" \
103 -> "D [m]" \
104 -< "cos #theta_{0}" \
105 -O COLZ \
106 -o ratio-${FUNCTION}-${dir}.$FORMAT $BATCH
107 done
108 fi
109done