4 if [ -z $JPP_DIR ]; then
5 echo "Variable JPP_DIR undefined."
9 source $JPP_DIR/setenv.sh $JPP_DIR >& /dev/null
13 set_variable DIR $JPP_DIR/examples/JAcoustics/
14 set_variable WORKDIR ./
15 set_variable: DEBUG ACOUSTICS_DEBUG 2
16 set_variable+ HYDROPHONE ACOUSTICS_HYDROPHONE -H hydrophone.txt
17 set_variable+ TRANSMITTER ACOUSTICS_TRANSMITTER -Y transmitter.txt
18 set_variable: SIGMA_S ACOUSTICS_SIGMA_S 50.0E-6
19 set_variable: NUMBER_OF_EVENTS ACOUSTICS_NUMBER_OF_EVENTS 1000
20 set_variable: FORMAT GRAPHICS_FORMAT gif
21 set_variable+ BATCH GRAPHICS_BATCH -B
24 usage "$script <detector file> <tripod file>"
27 if (( $# != 2 )); then
28 fatal "Wrong number of arguments."
31 set_variable DETECTOR $argv[1]
32 set_variable TRIPOD $argv[2]
34 if [[ ! -f $DETECTOR ]]; then
35 fatal "No file $DETECTOR."
38 if [[ ! -f $TRIPOD ]]; then
39 fatal "No file $TRIPOD."
42 eval `JPrintDetector -a $DETECTOR -O IDENTIFIER`
43 eval `JPrintDetector -a $DETECTOR -O SUMMARY`
51 -o footprint.$FORMAT $BATCH
54 if ( ! reuse_file $WORKDIR/katoomba.root ); then
58 -o $WORKDIR/katoomba.root \
61 -n $NUMBER_OF_EVENTS \
62 -@ acoustics_fit_parameters.txt \
63 -@ "sigma_s = $SIGMA_S;" \
72 rm -f $WORKDIR/strings.txt
74 for STRING in ${STRINGS[*]}; do
76 RESULT=(`JPrintResult -f "$WORKDIR/katoomba.root:string\[$STRING\]" -F "GetRMS(1)" -F "GetRMS(2)"`)
78 printf "%4d %9.5f %9.5f [mrad]\n" $STRING $RESULT[*]
80 printf "%9.5f\n" $((sqrt($RESULT[1]*$RESULT[2] + $RESULT[1]*$RESULT[2]))) >> $WORKDIR/strings.txt
83 $JPP_DIR/examples/JGizmo/JHistogram1D \
85 -f $WORKDIR/strings.txt \
86 -o $WORKDIR/strings.root
89 -f $WORKDIR/strings.root:\.\* \
90 -> "#sigma_{T} [mrad]" \
91 -\^ "number of strings" \
93 -o strings.$FORMAT $BATCH
95 rm -f $WORKDIR/strings.{txt,root}
101 -f "$WORKDIR/katoomba.root:cpu" \
103 -\^ "number of events [a.u.]" \
105 -o cpu.$FORMAT $BATCH
108 -f "$WORKDIR/katoomba.root:chi2" \
111 -\^ "number of events [a.u.]" \
114 -o chi2.$FORMAT $BATCH
119 source JAcousticsToolkit.sh
123 get_tripods+ tripod.txt EMITTERS
124 get_transmitters+ transmitter.txt EMITTERS
126 for EMITTER in % ${(@k)EMITTERS}; do
129 -f "$WORKDIR/katoomba.root:emitter\[${EMITTER}\]" \
131 -\^ "number of events [a.u.]" \
135 -o TX${EMITTER}.$FORMAT $BATCH
142 -o $WORKDIR/ellipse.root \
144 -@ "line.color = 2" \
147 for STRING in % $STRINGS[*]; do
150 -f "$WORKDIR/katoomba.root:string\[${STRING}\]" \
151 -f $WORKDIR/ellipse.root:\.\* \
153 -> "#DeltaTx [mrad]" \
154 -< "#DeltaTy [mrad]" \
160 -o RX${STRING}.$FORMAT $BATCH