Jpp test-rotations-new
the software that should make you happy
Loading...
Searching...
No Matches
hydrophone-Z: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/hydrophonez_${STRING}.txt
53set_array QUANTILES 0.333 0.666
54set_variable TMPDIR $WORKDIR/.$$
55mkdir -p $TMPDIR
56
57JKatoomba.sh $DETECTOR $INPUT_FILES[*] $TMPDIR/katoomba.root
58
59for (( Z = -7.5; $Z <= +7.5; Z += 0.1 )); do
60
61 JEditDetector -a $DETECTOR -o $TMPDIR/detector.datx -M "$MODULE add 0.0 0.0 $Z"
62
63 $JPP_DIR/examples/JAcoustics/JHydrophone \
64 -a $TMPDIR/detector.datx \
65 -f $TMPDIR/katoomba.root \
66 -o $TMPDIR/hydrophone.root \
67 -H $WORKDIR/hydrophone.txt \
68 -T $WORKDIR/tripod.txt \
69 -V $WORKDIR/sound_velocity.txt \
70 -W $WORKDIR/waveform.txt \
71 -Q $QUALITY_FACTOR \
72 -d $DEBUG --!
73
74 XMEAN=(0.0 0.0)
75 STDEV=(0.0 0.0)
76 TOTAL=0.0
77
78 for TRIPOD in ${(@k)TRIPODS}; do
79
80 JPrintQuantiles \
81 -f "$TMPDIR/hydrophone.root:${MODULE}\.${TRIPOD}" \
82 -Q "$QUANTILES[*]" \
83 -d 0 | read -A RESULT
84
85 if (( ${#RESULT} == ${#QUANTILES} )); then
86
87 VALUE=$((0.5e3 * ($RESULT[-1] + $RESULT[1]))) # average [ms]
88 TOTAL=$(($TOTAL + 1.0))
89
90 if (( $TOTAL == 1.0 )); then
91 XMEAN[1]=$VALUE
92 STDEV[1]=0.0
93 else
94 XMEAN[2]=$(($XMEAN[1] + ($VALUE - $XMEAN[1]) / $TOTAL))
95 STDEV[2]=$(($STDEV[1] + ($VALUE - $XMEAN[1]) * ($VALUE - $XMEAN[2])))
96 XMEAN[1]=$XMEAN[2]
97 STDEV[1]=$STDEV[2]
98 fi
99 fi
100 done
101
102 printf "%7.3f %9.5f\n" $Z $((sqrt($TOTAL * $STDEV[1] / ($TOTAL * ($TOTAL - 1.0))))) >> $HYDROPHONE_TXT
103done
104
105rm -rf $TMPDIR