Jpp  master_rocky
the software that should make you happy
detector-XY:fit.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 FORMULA "[0] + [1] * log(1.0 + 0.5*(x-[2])*(x-[2])/([4]*[4]) + 0.5*(y-[3])*(y-[3])/([5]*[5]))"
14 set_variable FORMULA "[0] + [1] * ((x-[2])*(x-[2])/([4]*[4]) + (y-[3])*(y-[3])/([5]*[5]))"
15 set_variable: FORMAT GRAPHICS_FORMAT gif
16 set_variable+ BATCH GRAPHICS_BATCH -B
17 
18 if do_usage $*; then
19  usage "$script <string identifier> [<detector file>]"
20 fi
21 
22 if (( $# == 0 || $# > 2 )); then
23  fatal "Wrong number of arguments."
24 fi
25 
26 set_variable STRING $argv[1]
27 
28 set_variable DETECTOR_TXT $WORKDIR/detectorxy_${STRING}.txt
29 set_variable DETECTOR_ROOT $WORKDIR/detectorxy_${STRING}.root
30 set_variable G1_TXT $WORKDIR/g1.txt
31 set_variable G1_ROOT $WORKDIR/g1.root
32 
33 if [[ ! -f $DETECTOR_TXT ]] then
34  fatal "No detector data file $DETECTOR_TXT."
35 fi
36 
37 JGraph2D \
38  -f $DETECTOR_TXT \
39  -o $DETECTOR_ROOT
40 
41 sort -gr -k 3,3 $DETECTOR_TXT | tail -1 | read X Y Z RMS
42 
43 echo $X $Y > $G1_TXT
44 
45 JGraph \
46  -f $G1_TXT \
47  -o $G1_ROOT
48 
49 awk '
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 
65 if (( $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 
99 else
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 
105 fi
106 
107 if (( $# == 2 )); then
108 
109  set_variable DETECTOR $argv[2]
110 
111  JEditDetector -a $DETECTOR -S "$STRING add $X $Y 0.0" -o $DETECTOR
112 fi
113 
114 typeset -Z 4 STRING
115 
116 JPlot2D \
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 
128 rm -f $DETECTOR_ROOT $G1_TXT $G1_ROOT
129 
130