Jpp  master_rocky-43-ge265d140c
the software that should make you happy
JPlotPD0.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 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.0
25 set_variable EAST +1.0
26 set_variable SOUTH 0.0
27 set_variable WEST -1.0
28 
29 set_variable WORKDIR ${TMPDIR:-/tmp}
30 set_variable E 10
31 set_variable R 25
32 set_variable YMAX 2e-3
33 set_variable YMIN 1e-6
34 
35 set_variable PDF $JPP_DATA/J%p.dat
36 set_variable HISTOGRAM "510 -20 +150"
37 
38 set_variable: FORMAT GRAPHICS_FORMAT gif
39 set_variable+ BATCH GRAPHICS_BATCH -B
40 
41 if do_usage $*; then
42  usage "$script [energy [distance]]"
43 fi
44 
45 case $# in
46  2) set_variable R $2;&
47  1) set_variable E $1;;
48 esac
49 
50 
51 if (( 1 )); then
52 
53  for dir in NORTH EAST SOUTH WEST; do
54 
55  echo "Generating $dir."
56 
57  eval D=\$$dir
58 
59  JPlotPD0 \
60  -f ${PDF/\%/24} \
61  -R $R -E $E -c ${D} \
62  -H "${HISTOGRAM}" \
63  -o $WORKDIR/F\[${dir}\].root
64 
65  done
66 
67 fi
68 
69 
70 if (( 1 )); then
71 
72  JPlot1D \
73  -f $WORKDIR/F\[NORTH\].root:h0 \
74  -f $WORKDIR/F\[EAST\].root:h0 \
75  -f $WORKDIR/F\[SOUTH\].root:h0 \
76  -f $WORKDIR/F\[WEST\].root:h0 \
77  -y "$YMIN $YMAX" -Y \
78  -> "#Deltat [ns]" \
79  -\^ "dP/dt [npe/ns]" \
80  -L TR \
81  -T "E = ${E} GeV#; R = ${R} m" \
82  -o $WORKDIR/PDG.$FORMAT $BATCH
83 
84 fi