Jpp test-rotations-old
the software that should make you happy
Loading...
Searching...
No Matches
JParramatta.sh
Go to the documentation of this file.
1#!/usr/bin/env zsh
2script=${0##*/}
3
4if [ -z $JPP_DIR ]; then
5 echo "Variable JPP_DIR undefined."
6 exit
7fi
8
9source $JPP_DIR/setenv.sh $JPP_DIR >& /dev/null
10
11set_variable: DEBUG ACOUSTICS_DEBUG 2
12set_variable: WORKDIR ACOUSTICS_WORKDIR ./
13set_variable: ZMAX_M ACOUSTICS_ZMAX_M 0.0
14set_variable: TIMESTAMP GRAPHICS_TIMESTAMP utc
15set_variable: FORMAT GRAPHICS_FORMAT gif
16set_variable: SIZE GRAPHICS_SIZE 0.5
17set_variable+ BATCH GRAPHICS_BATCH -B
18
19if do_usage $*; then
20 usage "$script (input file)+"\
21 "\nInput files correspond to the output of JKatoomba[.sh]."
22fi
23
24if (( $# < 1 )); then
25 fatal "Wrong number of arguments."
26fi
27
28set_variable INPUT_FILE $argv[1,-1]
29
30JParramatta \
31 -f "$INPUT_FILE[*]" \
32 -o $WORKDIR/parramatta.root \
33 -Z $ZMAX_M \
34 -d $DEBUG --!
35
36echo "Creating graphics output in directory $WORKDIR"
37
38JPlot1D \
39 -w 1200x600 \
40 -f $WORKDIR/parramatta.root:g1 \
41 -\^ "#chi^{2}/NDF" \
42 -N "X 505" \
43 -t "$TIMESTAMP" \
44 -S $SIZE \
45 -T "" \
46 -o $WORKDIR/chi2-time.$FORMAT $BATCH
47
48JPlot1D \
49 -w 1200x600 \
50 -f "$WORKDIR/parramatta.root:G\[[0-9][0-9]*\].amplitude" \
51 -\^ "amplitude [mrad]" \
52 -N "X 505" \
53 -t "$TIMESTAMP" \
54 -S $SIZE \
55 -T "" \
56 -L TR \
57 -o $WORKDIR/amplitude.$FORMAT $BATCH
58
59JPlot1D \
60 -w 1200x600 \
61 -f "$WORKDIR/parramatta.root:G\[[0-9][0-9]*\].stretching" \
62 -\^ "stretching [%]" \
63 -N "X 505" \
64 -t "$TIMESTAMP" \
65 -S $SIZE \
66 -T "" \
67 -L BR \
68 -o $WORKDIR/stretching.$FORMAT $BATCH
69
70JPlot1D \
71 -f "$WORKDIR/parramatta.root:hx" \
72 -f "$WORKDIR/parramatta.root:hy" \
73 -\^ "#DeltaTx/y [mrad]" \
74 -> "string" \
75 -T "" \
76 -o $WORKDIR/txy.$FORMAT $BATCH
77
78typeset -Z4 STRING
79
80for H1 in `JPrintName -f "$WORKDIR/parramatta.root:^\[[0-9][0-9]*\]\.tiltdeviation"`; do
81
82 set_variable STRING `echo $H1 | sed 's/.*\[\‍([0-9]*\‍)\].*/\1/'`
83
84 JPlot2D \
85 -f "$WORKDIR/parramatta.root:^\[${STRING}\]\.tiltdeviation" \
86 -< "Ty - <Ty> [mrad]" \
87 -> "Tx - <Tx> [mrad]" \
88 -O COLZ \
89 -G X \
90 -G Y \
91 -N "X 505" \
92 -N "Y 505" \
93 -s 101110 \
94 -T "[${STRING}]" \
95 -o $WORKDIR/tiltdeviation_${STRING}.$FORMAT $BATCH
96done
97
98if (( 0 )); then
99
100 set_variable FORMULA "[0] + [1]*x + [2]*x*x"
101
102 for H1 in `JPrintName -f "$WORKDIR/parramatta.root:stretching\[[0-9][0-9]*\]"`; do
103
104 set_variable STRING `echo $H1 | sed 's/.*\[\‍([0-9]*\‍)\].*/\1/'`
105
106 H1=${H1/\[/\\\[}
107 H1=${H1/\]/\\\]}
108
109 JFit \
110 -f "$WORKDIR/parramatta.root:${H1}" \
111 -o ${TMPDIR:-/tmp}/fit.root \
112 -F "$FORMULA" \
113 -@ "p0 = 0.0" \
114 -@ "p1 = 1.0e-5" \
115 -@ "p2 = 1.0e-6" \
116 -d $DEBUG
117
118 JPlot1D \
119 -f ${TMPDIR:-/tmp}/fit.root:\.\* \
120 -> "|T| [mrad]" \
121 -\^ "stretching [%]" \
122 -y "-0.25 +0.25" \
123 -N "X 505" \
124 -T "$STRING" \
125 -o $WORKDIR/profile_${STRING}.$FORMAT $BATCH
126 done
127fi