Jpp test-rotations-old
the software that should make you happy
Loading...
Searching...
No Matches
JAcousticsMonitorRateSummary.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: QUALITY_FACTOR ACOUSTICS_QUALITYFACTOR 2.0e3
14set_variable: FORMAT GRAPHICS_FORMAT gif
15set_variable LOGBOOK \"Individual Logbooks\"
16set_variable FACET "ELcode"
17
18if do_usage $*; then
19 usage "$script <detector file> <minrun> <maxrun> <acoustics_monitor_file>"\
20 "\nApply test criteria to the acoustic rate per emitter for each DOM."
21fi
22
23if (( $# == 4 )); then
24 set_variable DETECTOR $argv[1]
25 set_variable MINRUN $argv[2]
26 set_variable MAXRUN $argv[3]
27 set_variable PARAMETERS_FILE $argv[4]
28else
29 fatal "Wrong number of arguments."
30fi
31
32set_variable USER km3net
33set_variable PASSWORD pyrosoma
34
35JCookie.sh
36
37eval `JPrintDetector -a $DETECTOR -O IDENTIFIER`
38eval `JPrintDetector -a $DETECTOR -O SUMMARY`
39
40JAcoustics.sh $DETECTOR_ID
41
42source JAcousticsToolkit.sh
43
44CHECK_EXIT_CODE
45
46for ((RUN=$MINRUN; RUN<=$MAXRUN; RUN++)); do
47
48 set_variable INPUT_FILE $WORKDIR/KM3NeT_${(l:8::0::0:)DETECTOR_ID}_${(l:8::0::0:)RUN}_toa.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 STATUS_FILE $WORKDIR/status_rates_${(l:8::0::0:)RUN}.txt
53
54 echo "Using Parameter file " $PARAMETERS_FILE
55 if [[ ! -f $INPUT_FILE ]]; then
56
57 JToA \
58 -@ "detid = `getDetector -D $DETECTOR_ID`" \
59 -@ "minrun = $RUN" \
60 -@ "maxrun = $RUN" \
61 -@ "QUALITYFACTOR >= $QUALITY_FACTOR" \
62 -o $INPUT_FILE \
63 -d $DEBUG --!
64 fi
65
66 if [[ ! -f $OUTPUT_FILE ]]; then
67
68 JAcousticsEventBuilder \
69 -a $DETECTOR \
70 -f $INPUT_FILE \
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 \
77 -o $OUTPUT_FILE \
78 -d $DEBUG --!
79 fi
80
81 set_array time_range $(JRuns -D "$DETECTOR_ID" -@ "run = ${RUN}" -F UNIXJOBSTART -F UNIXJOBEND --)
82
83 set_variable lifetime_s $(( ($time_range[2]-$time_range[1])/1000 ))
84
85 #if lifetime is unknown, assume 6hr
86 if [ -z "$lifetime_s" ]; then
87 lifetime_s=21600
88
89 #avoid 0 sec lifetime
90 elif (($lifetime_s == 0)); then
91 lifetime_s=1
92 fi
93
94 $JPP_DIR/examples/JAcoustics/JAcousticsMonitorRateSummary \
95 -a $DETECTOR \
96 -f $OUTPUT_FILE \
97 -l $lifetime_s \
98 -o $MONITOR_FILE \
99 -d $DEBUG
100
101 $JPP_DIR/examples/JAcoustics/JAcousticsMonitorRateTest \
102 -i $MONITOR_FILE \
103 -a $DETECTOR \
104 -@ $PARAMETERS_FILE \
105 -F $FACET \
106 -s $SUMMARY_FILE \
107 -o $STATUS_FILE \
108 -r $RUN \
109 -d $DEBUG
110
111done