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 ./
13set_variable FORMULA "[0] + [1] * (sqrt(1.0 + 0.5*(x-[2])*(x-[2])/([4]*[4]) + 0.5*(y-[3])*(y-[3])/([5]*[5])) - 1.0)"
14set_variable: FORMAT GRAPHICS_FORMAT gif
15set_variable+ BATCH GRAPHICS_BATCH -B
18 usage "$script <string identifier> [<hydrophone file>]"
21if (( $# == 0 || $# > 2 )); then
22 fatal "Wrong number of arguments."
25set_variable STRING $argv[1]
27set_variable HYDROPHONE_TXT $WORKDIR/hydrophonexy_${STRING}.txt
28set_variable HYDROPHONE_ROOT $WORKDIR/hydrophonexy_${STRING}.root
29set_variable G1_TXT $WORKDIR/g1.txt
30set_variable G1_ROOT $WORKDIR/g1.root
32if [[ ! -f $HYDROPHONE_TXT ]] then
33 fatal "No hydrophone data file $HYDROPHONE_TXT."
40sort -gr -k 3,3 $HYDROPHONE_TXT | tail -1 | read X Y Z RMS
55 END { print XMIN, XMAX, YMIN, YMAX }
57 if ( $1 < XMIN ) { XMIN = $1 }
58 if ( $1 > XMAX ) { XMAX = $1 }
59 if ( $2 < YMIN ) { YMIN = $2 }
60 if ( $2 > YMAX ) { YMAX = $2 }
62 ' $HYDROPHONE_TXT | read XMIN XMAX YMIN YMAX
64if (( $X > $XMIN && $X < $XMAX && $Y > $YMIN && $Y < $YMAX )); then
72 -f ${HYDROPHONE_ROOT}:\.\* \
73 -o $WORKDIR/fit.root \
83 -R "p2 = $XMIN $XMAX" \
84 -R "p3 = $YMIN $YMAX" \
87 set_array P2 `JPrintFit -f $WORKDIR/fit.root:user -@p2`
88 set_array P3 `JPrintFit -f $WORKDIR/fit.root:user -@p3`
93 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"`
94 printf "optimal %4d %7.3f %7.3f %7.3f\n" $STRING $X $Y `JPrintFit -f $WORKDIR/fit.root:user -x "$X $Y"`
96 mv $WORKDIR/fit.root $HYDROPHONE_ROOT
100 echo "Optimum value at limit; -> no fit."
101 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 }}' $HYDROPHONE_TXT`
102 printf "optimal %4d %7.3f %7.3f %7.3f\n" $STRING $X $Y $Z
106if (( $# == 2 )); then
108 set_variable HYDROPHONE $argv[2]
110 JEditHydrophone -f $HYDROPHONE -S "$STRING set $X $Y 0.0" -o $HYDROPHONE
116 -f ${HYDROPHONE_ROOT}:\.\* \
125 -o $WORKDIR/hydrophonexy_${STRING}.$FORMAT $BATCH
127rm -f $HYDROPHONE_ROOT $G1_TXT $G1_ROOT