Jpp test-rotations-old
the software that should make you happy
Loading...
Searching...
No Matches
JF2.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 test JF2.
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 DEBUG 2
25set_variable WORKDIR ${TMPDIR:-/tmp}/
26
27if do_usage $*; then
28 usage "$script [working directory]"
29fi
30
31case $# in
32 1) set_variable WORKDIR $1;;
33esac
34
35
36if (( 1 )); then
37
38 set_variable FORMULA "[0] * exp(-0.5 * (x-[1])*(x-[1]) / ([2]*[2])) * exp(-0.5 * (y-[1])*(y-[1]) / ([2]*[2]))"
39
40 JF2 \
41 -o $WORKDIR/f2.root \
42 -F "$FORMULA" \
43 -@ "p0 = 1.0; p1 = 0.0; p2 = 1.0" \
44 -x "-5 +5" \
45 -y "-5 +5" \
46 -N 1000
47
48 JPlot2D \
49 -f $WORKDIR/f2.root:user \
50 -z "1e-10 2" -Z \
51 -O COLZ
52
53fi
54
55
56if (( 1 )); then
57
58 set_variable FORMULA "[0] * pow(x,[1])"
59
60 JF2 \
61 -o $WORKDIR/f2.root \
62 -F "$FORMULA" \
63 -@ "p0 = 1.0; p1 = -2.0" \
64 -x "1e0 1e7" \
65 -y "-1 +1"
66
67 JPlot2D \
68 -f $WORKDIR/f2.root:user \
69 -z "1e-14 1.0" \
70 -XZ \
71 -O COLZ
72
73fi
74
75
76if (( 1 )); then
77
78 set_variable FORMULA "[0] * pow(10,[1]*x)"
79
80 JF2 \
81 -o $WORKDIR/f2.root \
82 -F "$FORMULA" \
83 -@ "p0 = 1.0; p1 = -2.0" \
84 -x " 0 7" \
85 -y "-1 +1"
86
87 JPlot2D \
88 -f $WORKDIR/f2.root:user \
89 -z "1e-14 1.0" \
90 -Z \
91 -O COLZ
92
93 JPlot2D \
94 -f $WORKDIR/f2.root:user \
95 -z "1e-14 1.0" \
96 -XXZ \
97 -O COLZ
98
99fi
100
101
102if (( 1 )); then
103
104 set_variable FORMULA "[0] * pow(y,[1])"
105
106 JF2 \
107 -o $WORKDIR/f2.root \
108 -F "$FORMULA" \
109 -@ "p0 = 1.0; p1 = -2.0" \
110 -x "-1 +1" \
111 -y "1e0 1e7"
112
113 JPlot2D \
114 -f $WORKDIR/f2.root:user \
115 -z "1e-14 1.0" \
116 -YZ \
117 -O COLZ
118
119fi
120
121
122if (( 1 )); then
123
124 set_variable FORMULA "[0] * pow(10,[1]*y)"
125
126 JF2 \
127 -o $WORKDIR/f2.root \
128 -F "$FORMULA" \
129 -@ "p0 = 1.0; p1 = -2.0" \
130 -x "-1 +1" \
131 -y " 0 7"
132
133 JPlot2D \
134 -f $WORKDIR/f2.root:user \
135 -z "1e-14 1.0" \
136 -Z \
137 -O COLZ
138
139 JPlot2D \
140 -f $WORKDIR/f2.root:user \
141 -z "1e-14 1.0" \
142 -YYZ \
143 -O COLZ
144
145fi