1# awk program to scan loopstats files and report errors/statistics
2#
3# usage: awk -f loop.awk loopstats
4#
5# format of loopstats record
6#  MJD    sec   time (s)  freq (ppm)  poll
7# 49235  3.943  0.000016   22.4716    6
8#
9# format of output dataset (time values in milliseconds, freq in ppm)
10# loopstats.19960706
11# loop 1180, 0+/-11.0, rms 2.3, freq -24.45+/-0.045, var 0.019
12#
13BEGIN {
14	loop_tmax = loop_fmax = -1e9
15	loop_tmin = loop_fmin = 1e9
16}
17#
18# scan all records in file
19#
20{
21	if (NF >= 5) {
22		loop_count++
23		if ($3 > loop_tmax)
24			loop_tmax = $3
25		if ($3 < loop_tmin)
26			loop_tmin = $3
27		if ($4 > loop_fmax)
28			loop_fmax = $4
29		if ($4 < loop_fmin)
30			loop_fmin = $4
31		loop_time += $3
32		loop_time_rms += $3 * $3
33		loop_freq += $4
34		loop_freq_rms += $4 * $4
35	}
36} END {
37	if (loop_count > 0) {
38		loop_time /= loop_count
39                loop_time_rms = sqrt(loop_time_rms / loop_count - loop_time * loop_time)
40		loop_freq /= loop_count
41		loop_freq_rms = sqrt(loop_freq_rms / loop_count - loop_freq * loop_freq)
42		printf "loop %d, %.0f+/-%.1f, rms %.1f, freq %.2f+/-%0.3f, var %.3f\n", loop_count, (loop_tmax + loop_tmin) / 2 * 1e6, (loop_tmax - loop_tmin) / 2 * 1e6, loop_time_rms * 1e6, (loop_fmax + loop_fmin) / 2, (loop_fmax - loop_fmin) / 2, loop_freq_rms
43	}
44}
45
46