Jpp test-rotations-new
the software that should make you happy
Loading...
Searching...
No Matches
JPlotPDF_L1.sh
Go to the documentation of this file.
1#!/usr/bin/env zsh
2#
3# \author mdejong
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 T_NS 10
33set_variable: YMAX PHYSICS_YMAX 2e-0
34set_variable: YMIN PHYSICS_YMIN 1e-5
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]]"
44fi
45
46case $# in
47 2) set_variable R $2;&
48 1) set_variable E $1;;
49esac
50
51
52if (( 1 )); then
53
54 for dir in NORTH EAST SOUTH WEST; do
55
56 echo "Generating $dir."
57
58 eval D=\$$dir
59
60 $JPP_DIR/examples/JPhysics/JPlotPDF_L1 \
61 -f $PDF \
62 -D "${D}" -R $R -E $E \
63 -H "$HISTOGRAM" \
64 -T $T_NS \
65 -o $WORKDIR/F\[${dir}\].root
66
67 done
68
69fi
70
71
72if (( 1 )); then
73
74 JPlot1D \
75 -f $WORKDIR/F\[NORTH\].root:h0 \
76 -f $WORKDIR/F\[EAST\].root:h0 \
77 -f $WORKDIR/F\[SOUTH\].root:h0 \
78 -f $WORKDIR/F\[WEST\].root:h0 \
79 -y "$YMIN $YMAX" -Y \
80 -> "#Deltat [ns]" \
81 -\^ "dP/dt [npe/ns]" \
82 -L TR \
83 -T "E = ${E} GeV#; R = ${R} m" \
84 -o PDF.$FORMAT $BATCH
85
86fi