Jpp 19.3.0-rc.1
the software that should make you happy
Loading...
Searching...
No Matches
tripod-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
11set_variable: DEBUG ACOUSTICS_DEBUG 2
12set_variable: WORKDIR ACOUSTICS_WORKDIR ./
13
14if do_usage $*; then
15 usage "$script <tripod identifier> <detector file> (input file)+"
16fi
17
18if (( $# < 3 )); then
19 fatal "Wrong number of arguments."
20fi
21
22set_variable TRIPOD $argv[1]
23set_variable DETECTOR $argv[2]
24set_array INPUT_FILES $argv[3,-1]
25
26eval `JPrintDetector -a $DETECTOR -O IDENTIFIER`
27
28JAcoustics.sh $DETECTOR_ID
29
30source JAcousticsToolkit.sh
31
32CHECK_EXIT_CODE
33
34typeset -A TRIPODS
35
36get_tripods $WORKDIR/tripod.txt TRIPODS
37
38if [[ -z "${TRIPODS[$TRIPOD]}" ]]; then
39 fatal "Invalid tripod $TRIPOD; "\
40 "possible values: ${(@k)TRIPODS}"
41fi
42
43set_variable TRIPOD_TXT $WORKDIR/tripodz_${TRIPOD}.txt
44set_variable TMPDIR $WORKDIR/.$$
45mkdir -p $TMPDIR
46
47for (( Z = -2.0; $Z <= 2.0; Z += 0.1 )); do
48
49 JEditTripod -f $WORKDIR/tripod.txt -T "$TRIPOD add 0.0 0.0 $Z" -o $TMPDIR/tripod.txt
50
51 rm -f $TMPDIR/katoomba.root
52
53 JKatoomba \
54 -a $DETECTOR \
55 -f "$INPUT_FILES[*]" \
56 -o $TMPDIR/katoomba.root \
57 -T $TMPDIR/tripod.txt \
58 -V $WORKDIR/sound_velocity.txt \
59 -M $WORKDIR/mechanics.txt \
60 -@ $WORKDIR/acoustics_fit_parameters.txt \
61 -! $WORKDIR/disable.txt \
62 -d $DEBUG --!
63
64 if (( `JPrintResult -f $TMPDIR/katoomba.root:chi2 -F "GetSumOfWeights"` > 1.0 )); then
65
66 set_array RESULT `JPrintResult -f $TMPDIR/katoomba.root:chi2 -F "GetMean" -F "GetRMS / sqrt(GetSumOfWeights)"`
67
68 printf "%7.3f %7.3f %7.3f\n" $Z $RESULT[*] >> $TRIPOD_TXT
69 else
70 error "No result."
71 fi
72done
73
74rm -rf $TMPDIR