Jpp test-rotations-new
the software that should make you happy
Loading...
Searching...
No Matches
detector-Zmul: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 <string identifier> <detector file> (input file)+"
16fi
17
18if (( $# < 3 )); then
19 fatal "Wrong number of arguments."
20fi
21
22set_variable STRING $argv[1]
23set_variable DETECTOR $argv[2]
24set_array INPUT_FILES $argv[3,-1]
25
26eval `JPrintDetector -a $DETECTOR -O IDENTIFIER`
27eval `JPrintDetector -a $DETECTOR -O SUMMARY`
28
29set_variable MODULE `getModule -a $DETECTOR -L "$STRING 0"`
30
31JAcoustics.sh $DETECTOR_ID
32
33CHECK_EXIT_CODE
34
35if [[ -z "${STRINGS[(r)$STRING]}" ]]; then
36 fatal "Invalid string $STRING; "\
37 "possible values: ${STRINGS[*]}"
38fi
39
40set_variable DETECTOR_TXT $WORKDIR/detectorzmul_${STRING}.txt
41set_variable TMPDIR $WORKDIR/.$$
42mkdir -p $TMPDIR
43
44for (( Z = -9.0; $Z <= +1.0; Z += 0.1 )); do
45 for (( MUL = -0.01; $MUL <= 0.03; MUL += 0.001 )); do
46
47 JEditDetector \
48 -a $DETECTOR \
49 -k "$STRING" \
50 -s "$STRING add 0.0 0.0 $Z" \
51 -s "$STRING mul $MUL" \
52 -o $TMPDIR/detector.detx
53
54 rm -f $TMPDIR/katoomba.root
55
56 JKatoomba \
57 -a $TMPDIR/detector.detx \
58 -f "$INPUT_FILES[*]" \
59 -o $TMPDIR/katoomba.root \
60 -T $WORKDIR/tripod.txt \
61 -V $WORKDIR/sound_velocity.txt \
62 -M $WORKDIR/mechanics.txt \
63 -@ $WORKDIR/acoustics_fit_parameters.txt \
64 -! $WORKDIR/disable.txt \
65 -d $DEBUG --!
66
67 if (( `JPrintResult -f $TMPDIR/katoomba.root:chi2 -F "GetSumOfWeights"` > 1.0 )); then
68
69 set_array RESULT `JPrintResult -f $TMPDIR/katoomba.root:chi2 -F "GetMean" -F "GetRMS / sqrt(GetSumOfWeights)"`
70
71 printf "%7.3f %7.3f %7.3f %7.4f\n" $Z $MUL $RESULT[*] >> $DETECTOR_TXT
72 else
73 error "No result."
74 fi
75 done
76done
77
78rm -rf $TMPDIR