Jpp test-rotations-old
the software that should make you happy
Loading...
Searching...
No Matches
JPlotPDF.sh
Go to the documentation of this file.
1#!/usr/bin/env zsh
2#
3# \author mdejong, bofearraigh
4#
5version=1.0
6script=${0##*/}
7
8# ------------------------------------------------------------------------------------------
9#
10# Utility script to plot 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 NORTH "0.00 0.00"
25set_variable EAST "1.57 0.00"
26set_variable SOUTH "3.14 0.00"
27set_variable WEST "1.57 3.14"
28
29set_variable WORKDIR ${TMPDIR:-/tmp}
30set_variable E 1e3
31set_variable R 50
32set_variable z 60
33set_variable: YMAX PHYSICS_YMAX 2e-2
34set_variable: YMIN PHYSICS_YMIN 1e-7
35set_variable: HISTOGRAM PHYSICS_HISTOGRAM "1040 -20 +500"
36
37set_variable PDF $JPP_DATA/J%p.dat
38
39set_variable: FORMAT GRAPHICS_FORMAT gif
40set_variable+ BATCH GRAPHICS_BATCH -B
41
42if do_usage $*; then
43 usage "$script [energy [distance [z of PMT]]]"
44fi
45
46case $# in
47 3) set_variable z $3;&
48 2) set_variable R $2;&
49 1) set_variable E $1;;
50esac
51
52
53if (( 1 )); then
54
55 for dir in NORTH EAST SOUTH WEST; do
56
57 echo "Generating $dir."
58
59 eval D=\$$dir
60
61 JPlotPDF \
62 -f ${PDF/\%/1} \
63 -f ${PDF/\%/2} \
64 -f ${PDF/\%/3} \
65 -f ${PDF/\%/4} \
66 -f ${PDF/\%/5} \
67 -f ${PDF/\%/6} \
68 -D "${D}" -R $R -E $E -z $z \
69 -H "$HISTOGRAM" \
70 -o $WORKDIR/F\[${dir}\].root
71
72 done
73
74fi
75
76
77if (( 1 )); then
78
79 JPlot1D \
80 -f $WORKDIR/F\[NORTH\].root:h0 \
81 -f $WORKDIR/F\[EAST\].root:h0 \
82 -f $WORKDIR/F\[SOUTH\].root:h0 \
83 -f $WORKDIR/F\[WEST\].root:h0 \
84 -y "$YMIN $YMAX" -Y \
85 -> "#Deltat [ns]" \
86 -\^ "dP/dt [npe/ns]" \
87 -L TR \
88 -T "E = ${E} GeV#; R = ${R} m#; z = ${z} m" \
89 -o PDF.$FORMAT $BATCH
90
91fi