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])/([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
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/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
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
73 -f ${DETECTOR_ROOT}:\.\* \
74 -o $WORKDIR/fit.root \
84 -R "p2 = $XMIN $XMAX" \
85 -R "p3 = $YMIN $YMAX" \
88 set_array P2 `JPrintFit -f $WORKDIR/fit.root:user -@p2`
89 set_array P3 `JPrintFit -f $WORKDIR/fit.root:user -@p3`
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"`
97 mv $WORKDIR/fit.root $DETECTOR_ROOT
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
107if (( $# == 2 )); then
109 set_variable DETECTOR $argv[2]
111 JEditDetector -a $DETECTOR -S "$STRING add $X $Y 0.0" -o $DETECTOR
117 -f ${DETECTOR_ROOT}:\.\* \
126 -o $WORKDIR/detectorxy_${STRING}.$FORMAT $BATCH
128rm -f $DETECTOR_ROOT $G1_TXT $G1_ROOT