4if [ -z $JPP_DIR ]; then
5 echo "Variable JPP_DIR undefined."
9source $JPP_DIR/setenv.sh $JPP_DIR >& /dev/null
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])/([5]*[5]) + 0.5*(y-[3])*(y-[3])/([6]*[6]) + [4]*(x-[2])*(y-[3])/([5]*[6]))"
14set_variable FORMULA "[0] + [1] * ((x-[2])*(x-[2])/([5]*[5]) + (y-[3])*(y-[3])/([6]*[6]) + [4]*(x-[2])*(y-[3])/([5]*[6]))"
15set_variable: FORMAT GRAPHICS_FORMAT gif
16set_variable+ BATCH GRAPHICS_BATCH -B
19 usage "$script <string identifier> [<detector file>]"
22if (( $# == 0 || $# > 2 )); then
23 fatal "Wrong number of arguments."
26set_variable STRING $argv[1]
28set_variable DETECTOR_TXT $WORKDIR/detectorzmul_${STRING}.txt
29set_variable DETECTOR_ROOT $WORKDIR/detectorzmul_${STRING}.root
30set_variable G1_TXT $WORKDIR/g1.txt
31set_variable G1_ROOT $WORKDIR/g1.root
33if [[ ! -f $DETECTOR_TXT ]] then
34 fatal "No detector data file $DETECTOR_TXT."
41sort -gr -k 3,3 $DETECTOR_TXT | tail -1 | read X Y Z RMS
56 END { print XMIN, XMAX, YMIN, YMAX }
58 if ( $1 < XMIN ) { XMIN = $1 }
59 if ( $1 > XMAX ) { XMAX = $1 }
60 if ( $2 < YMIN ) { YMIN = $2 }
61 if ( $2 > YMAX ) { YMAX = $2 }
63 ' $DETECTOR_TXT | read XMIN XMAX YMIN YMAX
65if (( $X > $XMIN && $X < $XMAX && $Y > $YMIN && $Y < $YMAX )); then
69 let "YMIN = $Y - 0.005"
70 let "YMAX = $Y + 0.005"
73 -f ${DETECTOR_ROOT}:\.\* \
74 -o $WORKDIR/fit.root \
85 -R "p2 = $XMIN $XMAX" \
86 -R "p3 = $YMIN $YMAX" \
90 set_array P2 `JPrintFit -f $WORKDIR/fit.root:user -@p2`
91 set_array P3 `JPrintFit -f $WORKDIR/fit.root:user -@p3`
96 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 "`
97 printf "optimal %4d %7.3f %7.3f %7.3f\n" $STRING $X $Y `JPrintFit -f $WORKDIR/fit.root:user -x "$X $Y"`
99 mv $WORKDIR/fit.root $DETECTOR_ROOT
103 echo "Optimum value at limit; -> no fit."
104 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`
105 printf "optimal %4d %7.3f %7.3f %7.3f\n" $STRING $X $Y $Z
109if (( $# == 2 )); then
111 set_variable DETECTOR $argv[2]
112 set_variable MODULE `getModule -a $DETECTOR -L "$STRING 0"`
116 -s "$STRING add 0.0 0.0 $X" \
117 -s "$STRING mul $Y" \
124 -f ${DETECTOR_ROOT}:\.\* \
133 -o $WORKDIR/detectorzmul_${STRING}.$FORMAT $BATCH
135rm -f $DETECTOR_ROOT $G1_TXT $G1_ROOT