7doxygen=JDATABASE::JRunQuality
9if [ -z $JPP_DIR ]; then
10 echo "Variable JPP_DIR undefined."
14source $JPP_DIR/setenv.sh $JPP_DIR >& /dev/null
16# Auxiliary script to produce QA/QC data.
18# The output data are uploaded to the database with a tool that is part of the km3pipe package.
19# The locally produced data as well as the uploaded data can be viewed with applications/scripts JDataMonitor[.sh] and JDataQuality[.sh].
21# Note that if the list of values is extended, the include file JRunQuality.hh should accordingly be updated.
25set_variable: DEBUG QAQC_DEBUG 0
26set_variable: VERSION QAQC_VERSION D_1.0.0
27set_variable: WORKDIR QAQC_WORKDIR ${TMPDIR:-/tmp}/.$$
31 trigger3DShower.enabled
32 triggerMXShower.enabled
39TRIGGERS=(`JPrintTriggerBits`)
57 GIT JPP nb_of_meta UUID
58 detector run livetime_s UTCMin_s UTCMax_s
59 ${TRIGGER_PARAMETERS[*]}
63 in_sync out_sync DAQ WR HRV FIFO
64 PMTs MEAN_Rate_Hz RMS_Rate_Hz
65 hrv_fifo_failures duplic_timeslices
72 "$script (input file)+ <QA/QC file>" \
73 "\nAuxiliary script to produce the following QA/QC data:\n" \
75 "\nSee also: $JPP_PAGES/ Search \"JRunQuality\""
79 fatal "Wrong number of arguments."
82set_array INPUT_FILES ${argv[1,-2]}
83set_variable QAQC_FILE ${argv[-1]:a}
85for (( i = 1; $i <= ${#INPUT_FILES}; i += 1 )); do
86 if [[ -f $INPUT_FILES[$i] ]]; then
87 INPUT_FILES[$i]=${INPUT_FILES[$i]:a}
94set_variable ACOUSTICS_DEBUG ${DEBUG}
95set_variable COMPASS_DEBUG ${DEBUG}
96set_variable DETECTOR detector.detx # detector used during data taking
101# Locally install input files.
103if [[ ! -f sound_velocity.txt ]]; then
104 cat>sound_velocity.txt<<EOF
1051539.50 -17e-3 -2000.00
109if [[ ! -s $QAQC_FILE ]]; then
110 echo ${PARAMETERS//\.*/} > $QAQC_FILE
113for INPUT_FILE in ${INPUT_FILES[*]}; do
115 debug "Processing $INPUT_FILE"
119 JPrintMeta -v |& awk '/version/ { print $2 }' | read GIT_VERSION
121 debug "GIT version $GIT_VERSION"
123 set_array DAQ_VERSION `JPrintMeta -f ${INPUT_FILE} -k GIT`
124 debug "DAQ version $DAQ_VERSION[*]"
126 echo ${(s/./)DAQ_VERSION[1]} | read DAQ_MAJOR DAQ_MINOR DAQ_PATCH
128 if [[ -z $DAQ_VERSION[1] ]]; then
129 JPrintMeta -f ${INPUT_FILE} -k SVN | read DAQ_VERSION
132 getUUID -f ${INPUT_FILE} -d 0 | read UUID
135 set_array DAQ_HEADER `JPrintDAQHeader -f ${INPUT_FILE} -d 0`
136 set_variable DETECTOR_ID ${DAQ_HEADER[1]}
137 set_variable RUN ${DAQ_HEADER[2]}
139 JDetectorDB -D ${DETECTOR_ID} -r ${RUN} -o ${DETECTOR} -WW -d ${DEBUG} # install detector file as-is for a.o. JTriggerReprocessor
142 if [[ "$DAQ_MAJOR" != [0-9][0-9]* ]] || (( $DAQ_MAJOR <= 12 )); then # Jpp-version dependence
145 -P "-1 -1 set UDP_COUNTER_DISABLE" \
146 -P "-1 -1 set UDP_TRAILER_DISABLE" \
151 JAHRSCalibration.sh ${DETECTOR_ID} # install AHRS calibration file
154 set_array LIVETIME_S `getLivetime -f ${INPUT_FILE} -d 0`
155 set_array UTC_S `getUTC -f ${INPUT_FILE} -d 0`
157 {{ JTriggerMonitor -f ${INPUT_FILE} -d ${DEBUG} -Q 3 3>&1 1>&4 } | read -A TRIGGER_MONITOR } 4>&1
160 set_variable OUTPUT_FILE KM3NeT_${(l:8::0::0:)DETECTOR_ID}_${(l:8::0::0:)RUN}_reprocessor.root
162 JTriggerReprocessor \
171 JPrintTree -f ${OUTPUT_FILE} -@ "type = JDAQEvent" -k number_of_entries | read TRIGGER_REPROCESSOR
178 -Q 3 3>&1 1>&4 } | read -A TURBOT } 4>&1
181 {{ JSummaryMonitor -f ${INPUT_FILE} -t 50000 -d ${DEBUG} -Q 3 3>&1 1>&4 } | read -A SUMMARY_MONITOR } 4>&1
184 JAcousticsTriggerProcessor.sh ${DETECTOR} ${RUN}
186 if (( $? == 0 )); then
188 set_variable OUTPUT_FILE KM3NeT_${(l:8::0::0:)DETECTOR_ID}_${(l:8::0::0:)RUN}_event.root
190 NUMBER_OF_ACOUSTICS_EVENTS=(`JPrintTree -f ${OUTPUT_FILE} -@ "type = JACOUSTICS::JEvent" -k number_of_entries`)
192 NUMBER_OF_ACOUSTICS_EVENTS=(0)
195 JAHRS.sh ${DETECTOR_ID} ${RUN}
197 if (( $? == 0 )); then
199 set_variable OUTPUT_FILE KM3NeT_${(l:8::0::0:)DETECTOR_ID}_${(l:8::0::0:)RUN}_ahrs.root
204 -c ahrs_calibration.txt \
207 -Q 3 3>&1 1>&4 } | read -A AHRS_MONITOR } 4>&1
212 set_variable OUTPUT_FILE KM3NeT_${(l:8::0::0:)DETECTOR_ID}_${(l:8::0::0:)RUN}_trigger.root
218 -@ "trigger3DMuon.enabled = 1" \
219 -@ "trigger3DMuon.numberOfHits = 3" \
220 -@ "trigger3DMuon.numberOfModules = 3" \
221 -@ "trigger3DMuon.gridAngle_deg = 10" \
222 -@ "trigger3DMuon.TMaxExtra_ns = 500.0e3" \
223 -@ "TMaxEvent_ns = 500.0e3" \
224 -@ "TMaxLocal_ns = 20.0" \
225 -@ "L2.numberOfHits = 4" \
226 -@ "L2.TMaxLocal_ns = 20.0" \
227 -@ "L2.ctMin = 0.0" \
237 -Q 3 3>&1 1>&4 } | read -A SQUID } 4>&1
244 -Q 3 3>&1 1>&4 } | read -A MERMAID } 4>&1
248 printf " %12s" ${GIT_VERSION} >> $QAQC_FILE
249 printf " %12s" ${DAQ_VERSION[1]:-\?} >> $QAQC_FILE
250 printf " %1d" ${#DAQ_VERSION} >> $QAQC_FILE
251 printf " %s" ${UUID:-\?} >> $QAQC_FILE
252 printf " %8d" ${DAQ_HEADER[1,2]} >> $QAQC_FILE
253 printf " %9.1f" ${LIVETIME_S[1]} >> $QAQC_FILE
254 printf " %12.1f" ${UTC_S[1,2]} >> $QAQC_FILE
256 for KEY in ${TRIGGER_PARAMETERS[*]}; do
258 JPrintTriggerParameters -f ${INPUT_FILE} -k $KEY -d 0 | read VALUE
260 printf " %4d" ${VALUE:-0} >> $QAQC_FILE
263 for KEY in ${DAQ[*]}; do
265 JPrintTree -f ${INPUT_FILE} -@ "type = $KEY" -k number_of_entries | read VALUE
267 printf " %8d" ${VALUE:-0} >> $QAQC_FILE
270 printf " %8d" ${TRIGGER_REPROCESSOR:-0} >> $QAQC_FILE
271 printf " %8d" ${TRIGGER_MONITOR[*]} >> $QAQC_FILE
272 printf " %4d" ${TURBOT[1,2]} >> $QAQC_FILE
273 printf " %7.5f" ${SUMMARY_MONITOR[1,4]} >> $QAQC_FILE
274 printf " %8.1f" ${SUMMARY_MONITOR[5]} >> $QAQC_FILE
275 printf " %7.0f" ${SUMMARY_MONITOR[6,7]} >> $QAQC_FILE
276 printf " %8d" ${SUMMARY_MONITOR[8]} >> $QAQC_FILE
277 printf " %8d" ${SUMMARY_MONITOR[9]} >> $QAQC_FILE
278 printf " %4d" ${NUMBER_OF_ACOUSTICS_EVENTS[1]} >> $QAQC_FILE
279 printf " %5d" ${AHRS_MONITOR[1]} >> $QAQC_FILE
280 printf " %8.3f" ${AHRS_MONITOR[2]} >> $QAQC_FILE
281 printf " %4d" ${SQUID[1,2]} >> $QAQC_FILE
282 printf " %8d" ${MERMAID[1]} >> $QAQC_FILE
283 printf "\n" >> $QAQC_FILE