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