4if [ -z $JPP_DIR ]; then
5 echo "Variable JPP_DIR undefined."
9source $JPP_DIR/setenv.sh $JPP_DIR >& /dev/null
13set_variable: DEBUG ACOUSTICS_DEBUG 2
14set_variable: WORKDIR ACOUSTICS_WORKDIR ./
15set_variable: FORMAT GRAPHICS_FORMAT gif
16set_variable+ BATCH GRAPHICS_BATCH -B
17set_variable RANGE_S 0 60
20 usage "$script <detector file> (input file)+"\
21 "\nInput files correspond to the output of JAcousticsEventBuilder[.sh]."
25 fatal "Wrong number of arguments."
28set_variable DETECTOR $argv[1]
29set_array INPUT_FILES $argv[2,-1]
31eval `JPrintDetector -a $DETECTOR -O SUMMARY`
32eval `JPrintDetector -a $DETECTOR -O IDENTIFIER`
34JAcoustics.sh $DETECTOR_ID
36source JAcousticsToolkit.sh
42get_tripods+ $WORKDIR/tripod.txt EMITTERS
43get_transmitters+ $WORKDIR/transmitter.txt EMITTERS
45set_variable HISTOGRAM "4000 -1.0e-1 +1.0e-1"
46set_variable FORMULA "[0] * exp(-0.5 * (x-[1])*(x-[1]) / ([2]*[2])) + [3]"
50for EMITTER in -1 ${(@k)EMITTERS}; do
52 echo Processing $EMITTER...
54 $JPP_DIR/examples/JAcoustics/JToAMonitor \
56 -f "$INPUT_FILES[*]" \
57 -o $WORKDIR/toa.root \
64 -f "$WORKDIR/toa.root:\[%\]\.t1" \
65 -o ${TMPDIR:-/tmp}/fit.root \
67 -x "-4.0e-4 +4.0e-4" \
68 -@ "p0 = GetMaximum * 0.8" \
74 P0=($(JPrintFit -f ${TMPDIR:-/tmp}/fit.root:\.\* -@ p0))
75 P3=($(JPrintFit -f ${TMPDIR:-/tmp}/fit.root:\.\* -@ p3))
77 if (( $? == 0 )); then
78 BUFFER[$EMITTER]=$(($P0[1] / $P3[1]))
82 -f ${TMPDIR:-/tmp}/fit.root:\.\* \
83 -x "-1.0e-3 +1.0e-3" \
86 -\^ "number of events [a.u.]" \
90 -o $WORKDIR/T_${EMITTER}.$FORMAT $BATCH
93for KEY VALUE in ${(kv)BUFFER}; do
94 printf "Emitter %2d S/N %7.1f\n" $KEY $VALUE