4 if [ -z $JPP_DIR ]; then
5 echo "Variable JPP_DIR undefined."
9 source $JPP_DIR/setenv.sh $JPP_DIR >& /dev/null
11 set_variable: DEBUG ACOUSTICS_DEBUG 2
12 set_variable: WORKDIR ACOUSTICS_WORKDIR ./
13 set_variable: FORMAT GRAPHICS_FORMAT gif
14 set_variable LOGBOOK \"Individual Logbooks\"
15 set_variable AUTHOR \"GatiusOliver C\"
16 set_variable FACET "ELcode"
19 usage "$script <detector file> minrun maxrun report"\
20 "\nIn case of failures, if report = 1 an elog entry is created, if report =! 1 no elog entry is created."
23 if (( $# == 4 )); then
24 set_variable DETECTOR $argv[1]
25 set_variable MINRUN $argv[2]
26 set_variable MAXRUN $argv[3]
27 set_variable REPORT $argv[4]
29 fatal "Wrong number of arguments."
32 set_variable USER km3net
33 set_variable PASSWORD pyrosoma
37 eval `JPrintDetector -a $DETECTOR -O IDENTIFIER`
38 eval `JPrintDetector -a $DETECTOR -O SUMMARY`
40 JAcoustics.sh $DETECTOR_ID
42 source JAcousticsToolkit.sh
46 for ((RUN=$MINRUN; RUN<=$MAXRUN; RUN++)); do
48 set_variable INPUT_FILE $WORKDIR/KM3NeT_${(l:8::0::0:)DETECTOR_ID}_${(l:8::0::0:)RUN}_toashort.root
49 set_variable OUTPUT_FILE $WORKDIR/KM3NeT_${(l:8::0::0:)DETECTOR_ID}_${(l:8::0::0:)RUN}_event.root
50 set_variable MONITOR_FILE $WORKDIR/KM3NeT_${(l:8::0::0:)DETECTOR_ID}_${(l:8::0::0:)RUN}_monitor.root
51 set_variable SUMMARY_FILE $WORKDIR/summary_${(l:8::0::0:)RUN}.txt
52 set_variable TEST_FILE $WORKDIR/KM3NeT_${(l:8::0::0:)DETECTOR_ID}_${(l:8::0::0:)RUN}_test.root
53 set_variable PARAMETERS_FILE $JPP_DIR/examples/JAcoustics/acoustics_monitor_${(l:8::0::0:)DETECTOR_ID}.txt
55 if [[ ! -f $INPUT_FILE ]]; then
59 -@ "detid = `getDetector -D $DETECTOR_ID`" \
66 if [[ ! -f $OUTPUT_FILE ]]; then
68 JAcousticsEventBuilder \
71 -T $WORKDIR/tripod.txt \
72 -Y $WORKDIR/transmitter.txt \
73 -V $WORKDIR/sound_velocity.txt \
74 -W $WORKDIR/waveform.txt \
75 -@ $WORKDIR/acoustics_trigger_parameters.txt \
76 -H $WORKDIR/hydrophone.txt \
81 set_array time_range $(JRuns -D "$DETECTOR_ID" -@ "run = ${RUN}" -F UNIXJOBSTART -F UNIXJOBEND --)
83 set_variable lifetime_s $(( ($time_range[2]-$time_range[1])/1000 ))
85 #if lifetime is unknown, assume 6hr
86 if [ -z "$lifetime_s" ]; then
90 elif (($lifetime_s == 0)); then
94 $JPP_DIR/examples/JAcoustics/JAcousticsMonitor_short \
101 read -r NUMBER_OF_FAILURES <<< "$($JPP_DIR/examples/JAcoustics/JAcousticsMonitorTest \
104 -P $PARAMETERS_FILE \
113 if [[ ${NUMBER_OF_FAILURES[-1]} -gt 0 && $REPORT == 1 ]]; then
117 get_tripods+ $WORKDIR/tripod.txt EMITTERS
118 get_transmitters+ $WORKDIR/transmitter.txt EMITTERS
122 for EMITTER in ${(@k)EMITTERS}; do
125 -f "${MONITOR_FILE}:H\[${EMITTER}\].event" \
131 -T "Emitter ${EMITTER}" \
132 -o $WORKDIR/rate\[${EMITTER}\]_${RUN}.$FORMAT \
136 -f "${TEST_FILE}:H\[${EMITTER}\].rate-test" \
141 -T "Emitter ${EMITTER}" \
142 -o $WORKDIR/rate-test\[${EMITTER}\]_${RUN}.$FORMAT \
146 set_variable PLOT_ACOUS_RATE $WORKDIR/KM3NeT_${(l:8::0::0:)DETECTOR_ID}_${RUN}_acous_rate
147 set_variable PLOT_ACOUS_RATE_TEST $WORKDIR/KM3NeT_${(l:8::0::0:)DETECTOR_ID}_${RUN}_acous_rate-test
152 -title "Average acoustic rate \n Run: ${RUN}" \
153 $WORKDIR/rate\[*\]_${RUN}.$FORMAT \
154 $WORKDIR/${PLOT_ACOUS_RATE}.$FORMAT >& /dev/null
159 -title "In red modules measuring a rate out of the accepted rate range: 10% < rate < 102% of the expected rate \n Run: ${RUN}" \
160 $WORKDIR/rate-test\[*\]_${RUN}.$FORMAT \
161 $WORKDIR/${PLOT_ACOUS_RATE_TEST}.$FORMAT >& /dev/null
163 ATTACHMENTS+=" -f $WORKDIR/${PLOT_ACOUS_RATE}.$FORMAT"
164 ATTACHMENTS+=" -f $WORKDIR/${PLOT_ACOUS_RATE_TEST}.$FORMAT"
171 -a Subject=\"Warning detector ID ${DETECTOR_ID}\" \
172 -a Type=\"Acoustics monitoring\" \
178 #remove temporary files
179 rm -f $WORKDIR/rate*_${RUN}.$FORMAT
181 rm -f $WORKDIR/${PLOT_ACOUS_RATE}.$FORMAT
182 rm -f $WORKDIR/${PLOT_ACOUS_RATE_TEST}.$FORMAT