4if [ -z $JPP_DIR ]; then
5 echo "Variable JPP_DIR undefined."
9source $JPP_DIR/setenv.sh $JPP_DIR >& /dev/null
11set_variable: WORKDIR ACOUSTICS_WORKDIR ./
12set_variable: FORMAT GRAPHICS_FORMAT gif
15 usage "$script <detector file> minrun maxrun"
19 set_variable DETECTOR $argv[1]
20 set_variable MINRUN $argv[2]
21 set_variable MAXRUN $argv[3]
22 set_variable PARAMETERS_FILE $argv[4]
24 fatal "Wrong number of arguments."
29eval `JPrintDetector -a $DETECTOR -O IDENTIFIER`
30eval `JPrintDetector -a $DETECTOR -O SUMMARY`
32JAcoustics.sh $DETECTOR_ID
34source JAcousticsToolkit.sh
38set_variable MONITOR_FILE_WEEK $WORKDIR/KM3NeT_${(l:8::0::0:)DETECTOR_ID}_${MINRUN}_${MAXRUN}_monitor.root
39set_variable TEST_FILE_WEEK $WORKDIR/KM3NeT_${(l:8::0::0:)DETECTOR_ID}_${MINRUN}_${MAXRUN}_test.root
40set_variable PLOT_ACOUS_RATE $WORKDIR/KM3NeT_${(l:8::0::0:)DETECTOR_ID}_${MINRUN}_${MAXRUN}_acous_rate
41set_variable PLOT_ACOUS_RATE_TEST $WORKDIR/KM3NeT_${(l:8::0::0:)DETECTOR_ID}_${MINRUN}_${MAXRUN}_acous_rate-test
42set_variable NUM_RUNS $(($MAXRUN - $MINRUN + 1))
43set_array MONITOR_FILES
46for ((RUN=$MINRUN; RUN<=$MAXRUN; RUN++)); do
48 set_variable MONITOR_FILE $WORKDIR/KM3NeT_${(l:8::0::0:)DETECTOR_ID}_${(l:8::0::0:)RUN}_monitor.root
49 set_variable TEST_FILE $WORKDIR/KM3NeT_${(l:8::0::0:)DETECTOR_ID}_${(l:8::0::0:)RUN}_test.root
51 MONITOR_FILES+=($MONITOR_FILE)
52 TEST_FILES+=($TEST_FILE)
54 if [[ ! -f $MONITOR_FILE || ! -f $TEST_FILE ]]; then
55 $JPP_DIR/examples/JAcoustics/JAcousticsMonitorRateSummary.sh $DETECTOR $RUN $RUN $PARAMETERS_FILE 0
60hadd -f $MONITOR_FILE_WEEK $MONITOR_FILES
61hadd -f $TEST_FILE_WEEK $TEST_FILES
65get_tripods+ $WORKDIR/tripod.txt EMITTERS
66get_transmitters+ $WORKDIR/transmitter.txt EMITTERS
68for EMITTER in ${(@k)EMITTERS}; do
71 -f "${MONITOR_FILE_WEEK}:H\[${EMITTER}\].event" \
73 -o "scale.root" >& /dev/null
76 -f "scale.root:H\[${EMITTER}\].event" \
82 -T "Emitter ${EMITTER}" \
83 -o ${PLOT_ACOUS_RATE}_\[${EMITTER}\].$FORMAT -B
86 -f "${TEST_FILE_WEEK}:H\[${EMITTER}\].rate-test" \
88 -o "scale_test.root" >& /dev/null
91 -f "scale_test.root:H\[${EMITTER}\].rate-test" \
94 -\^ "fraction of runs outside the rate range" \
97 -T "Emitter ${EMITTER}" \
98 -o ${PLOT_ACOUS_RATE_TEST}_\[${EMITTER}\].$FORMAT -B
100 #remove temporary files
102 rm -f "scale_test.root"
105if [[ -f ${PLOT_ACOUS_RATE}.$FORMAT ]]; then
106 rm -f ${PLOT_ACOUS_RATE}.$FORMAT
109if [[ -f ${PLOT_ACOUS_RATE_TEST}.$FORMAT ]]; then
110 rm -f ${PLOT_ACOUS_RATE_TEST}.$FORMAT
113set_array start_time $(JRuns -D "$DETECTOR_ID" -@ "run = ${MINRUN}" -F UNIXJOBSTART -F UNIXJOBEND --)
114set_array end_time $(JRuns -D "$DETECTOR_ID" -@ "run = ${MAXRUN}" -F UNIXJOBSTART -F UNIXJOBEND --)
116set_variable start_time_h $(date '+%d-%m-%Y %H:%M' -d @$(($start_time[1]/1000)))
117set_variable end_time_h $(date '+%d-%m-%Y %H:%M' -d @$(($end_time[2]/1000)))
122 -title "Average acoustic rate \n Runs: ${MINRUN}-${MAXRUN}, Time: ${start_time_h} - ${end_time_h}" \
123 ${PLOT_ACOUS_RATE}_*.$FORMAT \
124 ${PLOT_ACOUS_RATE}.$FORMAT >& /dev/null
129 -title "Fraction of runs outside accepted rate range: 10% < rate < 102% of the expected rate \n Runs: ${MINRUN}-${MAXRUN}, Time: ${start_time_h} - ${end_time_h}" \
130 ${PLOT_ACOUS_RATE_TEST}_*.$FORMAT \
131 ${PLOT_ACOUS_RATE_TEST}.$FORMAT >& /dev/null
134rm -f ${PLOT_ACOUS_RATE}_*.$FORMAT
135rm -f ${PLOT_ACOUS_RATE_TEST}_*.$FORMAT