Jpp  19.1.0-rc.1
the software that should make you happy
JParramatta.sh
Go to the documentation of this file.
1 #!/bin/zsh
2 script=${0##*/}
3 
4 if [ -z $JPP_DIR ]; then
5  echo "Variable JPP_DIR undefined."
6  exit
7 fi
8 
9 source $JPP_DIR/setenv.sh $JPP_DIR >& /dev/null
10 
11 set_variable: DEBUG ACOUSTICS_DEBUG 2
12 set_variable: WORKDIR ACOUSTICS_WORKDIR ./
13 set_variable: ZMAX_M ACOUSTICS_ZMAX_M 0.0
14 set_variable: TIMESTAMP GRAPHICS_TIMESTAMP utc
15 set_variable: FORMAT GRAPHICS_FORMAT gif
16 set_variable: SIZE GRAPHICS_SIZE 0.5
17 set_variable+ BATCH GRAPHICS_BATCH -B
18 
19 if do_usage $*; then
20  usage "$script (input file)+"\
21  "\nInput files correspond to the output of JKatoomba[.sh]."
22 fi
23 
24 if (( $# < 1 )); then
25  fatal "Wrong number of arguments."
26 fi
27 
28 set_variable INPUT_FILE $argv[1,-1]
29 
30 JParramatta \
31  -f "$INPUT_FILE[*]" \
32  -o $WORKDIR/parramatta.root \
33  -Z $ZMAX_M \
34  -d $DEBUG --!
35 
36 echo "Creating graphics output in directory $WORKDIR"
37 
38 JPlot1D \
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 
48 JPlot1D \
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 
59 JPlot1D \
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 
70 JPlot1D \
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 
78 typeset -Z4 STRING
79 
80 for 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
96 done
97 
98 if (( 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
127 fi