#!/usr/bin/perl
#======================================================================
#                    S C A N B I N S 
#                    doc: Mon Jan 27 17:55:34 2003
#                    dlm: Thu Mar 17 07:45:37 2016
#                    (c) 2003 A.M. Thurnherr
#                    uE-Info: 20 22 NIL 0 0 72 2 2 4 NIL ofnI
#======================================================================

# Collect Per-Bin Stats
#	NB: currently broken

# HISTORY:
#	Jan 27, 2003: - created
#	Sep 19, 2007: - adapted to new [RDI_BB_Read.pl] (not tested)
#   Mar 17, 2016: - adapted to new Getopt library
#				  - updated ancient library name

$0 =~ m{(.*)/[^/]+}; 
require "$1/RDI_PD0_IO.pl";
use Getopt::Std;

die("Usage: $0 " .
	"" .
	"<RDI file>\n")
		unless (&getopts("") && @ARGV == 1);

print(STDERR "Reading $ARGV[0]...");
readData($ARGV[0],\%dta);									# read data
print(STDERR "done\n");

for ($e=0; $e<=$#{$dta{ENSEMBLE}}; $e++) {
	checkEnsemble(\%dta,$e);									# sanity checks
	$nens++;

	next unless (defined($dta{ENSEMBLE}[$e]->{VELOCITY}[0][0]) &&
				 defined($dta{ENSEMBLE}[$e]->{VELOCITY}[0][1]) &&
				 defined($dta{ENSEMBLE}[$e]->{VELOCITY}[0][2]) &&
				 defined($dta{ENSEMBLE}[$e]->{VELOCITY}[0][3]));
	$ngoodens++;
	
	for ($b=0; $b<$dta{N_BINS}; $b++) {					# collect stats
		my($ngood) = defined($dta{ENSEMBLE}[$e]->{VELOCITY}[$b][0])
				   + defined($dta{ENSEMBLE}[$e]->{VELOCITY}[$b][1])
				   + defined($dta{ENSEMBLE}[$e]->{VELOCITY}[$b][2])
				   + defined($dta{ENSEMBLE}[$e]->{VELOCITY}[$b][3]);
		if 	  ($ngood == 4) { $ngood4[$b]++; }
		elsif ($ngood == 3) { $ngood3[$b]++; }
		else				{ $nbad[$b]++; }

		for ($i=0; $i<4; $i++) {
			if (defined($dta{ENSEMBLE}[$e]->{VELOCITY}[$b][$i])) {
				$ngood[$b][$i]++;
				$sumcor[$b][$i] += $dta{ENSEMBLE}[$e]->{CORRELATION}[$b][$i];
				$sumamp[$b][$i] += $dta{ENSEMBLE}[$e]->{ECHO_AMPLITUDE}[$b][$i];
			}
		}
	}
}

printf("$ngoodens good ensembles out of $nens\n");
for ($b=0; $b<$dta{N_BINS}; $b++) {						# gen output
	printf("%2d: vels: %3d%% 4-bin, %3d%% 3-bin, %3d%% bad; ",
		$b+1,
		100*$ngood4[$b]/$ngoodens,100*$ngood3[$b]/$ngoodens,
		100*$nbad[$b]/$ngoodens);
	printf("mean corr: %3d/%3d/%3d/%3d; mean amp: %3d/%3d/%3d/%3d",
		$sumcor[$b][0]/$ngood[$b][0], $sumcor[$b][1]/$ngood[$b][1],
		$sumcor[$b][2]/$ngood[$b][2], $sumcor[$b][3]/$ngood[$b][3],
		$sumamp[$b][0]/$ngood[$b][0], $sumamp[$b][1]/$ngood[$b][1],
		$sumamp[$b][2]/$ngood[$b][2], $sumamp[$b][3]/$ngood[$b][3]);
	print("\n");
}

exit(0);