Jpp test-rotations-old
the software that should make you happy
Loading...
Searching...
No Matches
hydrophone-phi:run.sh
Go to the documentation of this file.
1#!/usr/bin/env zsh
2script=${0##*/}
3
4if [ -z $JPP_DIR ]; then
5 echo "Variable JPP_DIR undefined."
6 exit
7fi
8
9source $JPP_DIR/setenv.sh $JPP_DIR >& /dev/null
10
11zmodload zsh/mathfunc
12
13set_variable: DEBUG ACOUSTICS_DEBUG 2
14set_variable: WORKDIR ACOUSTICS_WORKDIR ./
15set_variable: QUALITY_FACTOR ACOUSTICS_QUALITYFACTOR 2.0e3
16
17if do_usage $*; then
18 usage "$script <string identifier> <detector file> (input file)+"
19fi
20
21if (( ${#} < 3 )); then
22 fatal "Wrong number of arguments."
23fi
24
25set_variable STRING $argv[1]
26set_variable DETECTOR $argv[2]
27set_array INPUT_FILES $argv[3,-1]
28
29eval `JPrintDetector -a $DETECTOR -O IDENTIFIER`
30eval `JPrintDetector -a $DETECTOR -O SUMMARY`
31
32JAcoustics.sh $DETECTOR_ID
33
34source JAcousticsToolkit.sh
35
36CHECK_EXIT_CODE
37
38typeset -A TRIPODS
39
40get_tripods $WORKDIR/tripod.txt TRIPODS
41
42if (( ${#TRIPODS} < 2 )); then
43 fatal "Not enough tripods " ${#TRIPODS}
44fi
45
46if [[ -z "${STRINGS[(r)$STRING]}" ]]; then
47 fatal "Invalid string $STRING; "\
48 "possible values: ${STRINGS[*]}"
49fi
50
51set_variable MODULE `getModule -a $DETECTOR -L "$STRING 0"`
52set_variable HYDROPHONE_TXT $WORKDIR/hydrophonephi_${STRING}.txt
53set_array QUANTILES 0.333 0.666
54set_variable PI 3.1415927
55set_variable TMPDIR $WORKDIR/.$$
56mkdir -p $TMPDIR
57
58JKatoomba.sh $DETECTOR $INPUT_FILES[*] $TMPDIR/katoomba.root
59
60for (( PHI = -$PI; $PHI <= +$PI; PHI += 0.02 )); do
61
62 JEditHydrophone -f $WORKDIR/hydrophone.txt -S "$STRING rot $PHI" -o $TMPDIR/hydrophone.txt
63
64 $JPP_DIR/examples/JAcoustics/JHydrophone \
65 -a $DETECTOR \
66 -f $TMPDIR/katoomba.root \
67 -o $TMPDIR/hydrophone.root \
68 -H $TMPDIR/hydrophone.txt \
69 -T $WORKDIR/tripod.txt \
70 -V $WORKDIR/sound_velocity.txt \
71 -W $WORKDIR/waveform.txt \
72 -Q $QUALITY_FACTOR \
73 -d $DEBUG --!
74
75 XMEAN=(0.0 0.0)
76 STDEV=(0.0 0.0)
77 TOTAL=0.0
78
79 for TRIPOD in ${(@k)TRIPODS}; do
80
81 JPrintQuantiles \
82 -f "$TMPDIR/hydrophone.root:${MODULE}\.${TRIPOD}" \
83 -Q "$QUANTILES[*]" \
84 -d 0 | read -A RESULT
85
86 if (( ${#RESULT} == ${#QUANTILES} )); then
87
88 VALUE=$((0.5e3 * ($RESULT[-1] + $RESULT[1]))) # average [ms]
89 TOTAL=$(($TOTAL + 1.0))
90
91 if (( $TOTAL == 1.0 )); then
92 XMEAN[1]=$VALUE
93 STDEV[1]=0.0
94 else
95 XMEAN[2]=$(($XMEAN[1] + ($VALUE - $XMEAN[1]) / $TOTAL))
96 STDEV[2]=$(($STDEV[1] + ($VALUE - $XMEAN[1]) * ($VALUE - $XMEAN[2])))
97 XMEAN[1]=$XMEAN[2]
98 STDEV[1]=$STDEV[2]
99 fi
100 fi
101 done
102
103 printf "%7.3f %9.5f\n" $PHI $((sqrt($TOTAL * $STDEV[1] / ($TOTAL * ($TOTAL - 1.0))))) >> $HYDROPHONE_TXT
104done
105
106rm -rf $TMPDIR