Jpp test-rotations-old
the software that should make you happy
Loading...
Searching...
No Matches
detector-XY:fit.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 ./
13#set_variable FORMULA "[0] + [1] * log(1.0 + 0.5*(x-[2])*(x-[2])/([4]*[4]) + 0.5*(y-[3])*(y-[3])/([5]*[5]))"
14set_variable FORMULA "[0] + [1] * ((x-[2])*(x-[2])/([4]*[4]) + (y-[3])*(y-[3])/([5]*[5]))"
15set_variable: FORMAT GRAPHICS_FORMAT gif
16set_variable+ BATCH GRAPHICS_BATCH -B
17
18if do_usage $*; then
19 usage "$script <string identifier> [<detector file>]"
20fi
21
22if (( $# == 0 || $# > 2 )); then
23 fatal "Wrong number of arguments."
24fi
25
26set_variable STRING $argv[1]
27
28set_variable DETECTOR_TXT $WORKDIR/detectorxy_${STRING}.txt
29set_variable DETECTOR_ROOT $WORKDIR/detectorxy_${STRING}.root
30set_variable G1_TXT $WORKDIR/g1.txt
31set_variable G1_ROOT $WORKDIR/g1.root
32
33if [[ ! -f $DETECTOR_TXT ]] then
34 fatal "No detector data file $DETECTOR_TXT."
35fi
36
37JGraph2D \
38 -f $DETECTOR_TXT \
39 -o $DETECTOR_ROOT
40
41sort -gr -k 3,3 $DETECTOR_TXT | tail -1 | read X Y Z RMS
42
43echo $X $Y > $G1_TXT
44
45JGraph \
46 -f $G1_TXT \
47 -o $G1_ROOT
48
49awk '
50 BEGIN {
51 XMIN = +1e3
52 XMAX = -1e3
53 YMIN = +1e3
54 YMAX = -1e3
55 }
56 END { print XMIN, XMAX, YMIN, YMAX }
57 $1 !~ /^#/ {
58 if ( $1 < XMIN ) { XMIN = $1 }
59 if ( $1 > XMAX ) { XMAX = $1 }
60 if ( $2 < YMIN ) { YMIN = $2 }
61 if ( $2 > YMAX ) { YMAX = $2 }
62 }
63 ' $DETECTOR_TXT | read XMIN XMAX YMIN YMAX
64
65if (( $X > $XMIN && $X < $XMAX && $Y > $YMIN && $Y < $YMAX )); then
66
67 let "XMIN = $X - 0.5"
68 let "XMAX = $X + 0.5"
69 let "YMIN = $Y - 0.5"
70 let "YMAX = $Y + 0.5"
71
72 JFit2D \
73 -f ${DETECTOR_ROOT}:\.\* \
74 -o $WORKDIR/fit.root \
75 -F "$FORMULA" \
76 -x "$XMIN $XMAX" \
77 -y "$YMIN $YMAX" \
78 -@ "p0 = $Z" \
79 -@ "p1 = 2.0" \
80 -@ "p2 = $X" \
81 -@ "p3 = $Y" \
82 -@ "p4 = 0.5" \
83 -@ "p5 = 0.5" \
84 -R "p2 = $XMIN $XMAX" \
85 -R "p3 = $YMIN $YMAX" \
86 -d $DEBUG
87
88 set_array P2 `JPrintFit -f $WORKDIR/fit.root:user -@p2`
89 set_array P3 `JPrintFit -f $WORKDIR/fit.root:user -@p3`
90
91 set_variable X $P2[1]
92 set_variable Y $P3[1]
93
94 printf "nominal %4d %7.3f %7.3f %7.3f\n" $STRING 0.0 0.0 `JPrintFit -f $WORKDIR/fit.root:user -x "0.0 0.0"`
95 printf "optimal %4d %7.3f %7.3f %7.3f\n" $STRING $X $Y `JPrintFit -f $WORKDIR/fit.root:user -x "$X $Y"`
96
97 mv $WORKDIR/fit.root $DETECTOR_ROOT
98
99else
100
101 echo "Optimum value at limit; -> no fit."
102 printf "nominal %4d %7.3f %7.3f %7.3f\n" $STRING 0.0 0.0 `awk -v X=0.0 -v Y=0.0 '{ if ($1 == X && $2 == Y) { print $3 }}' $DETECTOR_TXT`
103 printf "optimal %4d %7.3f %7.3f %7.3f\n" $STRING $X $Y $Z
104
105fi
106
107if (( $# == 2 )); then
108
109 set_variable DETECTOR $argv[2]
110
111 JEditDetector -a $DETECTOR -S "$STRING add $X $Y 0.0" -o $DETECTOR
112fi
113
114typeset -Z 4 STRING
115
116JPlot2D \
117 -f ${DETECTOR_ROOT}:\.\* \
118 -f ${G1_ROOT}:\.\* \
119 -> "#Deltax [m]" \
120 -< "#Deltay [m]" \
121 -\^ "#chi^{2}/NDF" \
122 -N "X 505" \
123 -N "Y 505" \
124 -O CONT3 \
125 -T "[$STRING]" \
126 -o $WORKDIR/detectorxy_${STRING}.$FORMAT $BATCH
127
128rm -f $DETECTOR_ROOT $G1_TXT $G1_ROOT
129
130