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