Jpp 19.3.0-rc.1
the software that should make you happy
Loading...
Searching...
No Matches
software/JAcoustics/plot-canberra.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 ./
13set_variable: XMIN CANBERRA_XMIN -1.0e-3
14set_variable: XMAX CANBERRA_XMAX +1.0e-3
15set_variable: ZMIN CANBERRA_ZMIN -0.5
16set_variable: ZMAX CANBERRA_ZMAX +0.5
17set_variable+ DISABLE CANBERRA_DISABLE -! $WORKDIR/disable.txt
18set_variable+ REVERT CANBERRA_REVERT -r
19set_variable: OPTION CANBERRA_OPTION -O HIST
20set_variable: FORMAT GRAPHICS_FORMAT gif
21set_variable+ BATCH GRAPHICS_BATCH -B
22
23if do_usage $*; then
24 usage "$script <detector file> (input file)+"\
25 "\nInput file corresponds to the output of JKatoomba[.sh]."
26fi
27
28if (( $# < 2 )); then
29 fatal "Wrong number of arguments."
30fi
31
32set_variable DETECTOR $argv[1]
33set_array INPUT_FILES $argv[2,-1]
34
35eval `JPrintDetector -a $DETECTOR -O IDENTIFIER`
36eval `JPrintDetector -a $DETECTOR -O SUMMARY`
37
38JAcoustics.sh $DETECTOR_ID
39
40JCanberra.sh $DETECTOR $INPUT_FILES[*] $WORKDIR/canberra\[%\].root
41
42source JAcousticsToolkit.sh
43
44CHECK_EXIT_CODE
45
46typeset -A EMITTERS
47
48get_tripods+ $WORKDIR/tripod.txt EMITTERS
49get_transmitters+ $WORKDIR/transmitter.txt EMITTERS
50
51
52typeset -Z 4 STRING
53typeset -Z 2 FLOOR
54
55for STRING in $STRINGS[*]; do
56
57 echo -n "Creating graphics for string $STRING.."
58
59 for (( FLOOR = $FIRST_FLOOR; $FLOOR <= $LAST_FLOOR; FLOOR += 1 )); do
60
61 eval MODULE=\${STRING_${(l:4::0::0:)STRING}\[$(($FLOOR))\]:--1}
62
63 INPUT_FILES=(`ls $WORKDIR/canberra\[*\].root`)
64
65 JPlot1D \
66 -f${^INPUT_FILES}:${MODULE} \
67 -> "#Deltat [s]" \
68 -\^ "number of events [a.u.]" \
69 -x "$XMIN $XMAX" \
70 $OPTION \
71 -N "X 505" \
72 -L TR \
73 -T "(${STRING},${FLOOR})" \
74 -o ${TMPDIR:-/tmp}/toa_${STRING}_${FLOOR}.$FORMAT $BATCH
75 done
76
77 montage \
78 -tile 6x3 \
79 -geometry +0+0 \
80 ${TMPDIR:-/tmp}/toa_${STRING}_{01..18}.$FORMAT \
81 $WORKDIR/toa_${STRING}${REVERT}.$FORMAT >& /dev/null
82
83 mv ${TMPDIR:-/tmp}/toa_${STRING}_00.$FORMAT \
84 $WORKDIR/toa_${STRING}_00${REVERT}.$FORMAT
85
86 echo " -> $WORKDIR/toa_${STRING}${REVERT}.$FORMAT"
87
88 rm -f ${TMPDIR:-/tmp}/toa_${STRING}_*.$FORMAT
89done
90
91
92for EMITTER in ${(@k)EMITTERS}; do
93
94 JPlot2D \
95 -f $WORKDIR/canberra\[${EMITTER}\].root:mean \
96 -> "string" \
97 -< "floor" \
98 -\^ "#Deltat [ms]" \
99 -z "$ZMIN $ZMAX" \
100 -N "Y 10" \
101 -T "[${EMITTER}]" \
102 -O COLZ \
103 -p 57 \
104 -o $WORKDIR/mean_${EMITTER}${REVERT}.$FORMAT $BATCH
105
106 JPlot2D \
107 -f $WORKDIR/canberra\[${EMITTER}\].root:sigma \
108 -> "string" \
109 -< "floor" \
110 -\^ "#sigma [ms]" \
111 -z "0.0 0.1" \
112 -N "Y 10" \
113 -T "[${EMITTER}]" \
114 -O COLZ \
115 -p 57 \
116 -o $WORKDIR/sigma_${EMITTER}${REVERT}.$FORMAT $BATCH
117done
118
119typeset -i NX
120typeset -i NY
121
122let "NY = sqrt(${#EMITTERS})"
123let "NX = (${#EMITTERS} + $NY - 1)/ $NY"
124
125for KEY in mean sigma; do
126
127 montage \
128 -tile ${NX}x${NY} \
129 -geometry +0+0 \
130 $WORKDIR/${KEY}_*${REVERT}.$FORMAT \
131 $WORKDIR/${KEY}${REVERT}.$FORMAT >& /dev/null
132done
133
134
135rm -f $WORKDIR/canberra\[*\].root