Index: trunk/ippScripts/scripts/detrend_reject_exp.pl
===================================================================
--- trunk/ippScripts/scripts/detrend_reject_exp.pl	(revision 9457)
+++ trunk/ippScripts/scripts/detrend_reject_exp.pl	(revision 9458)
@@ -98,6 +98,6 @@
 my $meanStats = Statistics::Descriptive::Sparse->new(); # Statistics calculator
 $meanStats->add_data(@means);
-my $variancestats = Statistics::Descriptive::Sparse->new(); # Statistics calculator
-$variancestats->add_data(@variances);
+my $varianceStats = Statistics::Descriptive::Sparse->new(); # Statistics calculator
+$varianceStats->add_data(@variances);
 my $meanStdevStats = Statistics::Descriptive::Sparse->new(); # Statistics calculator
 $meanStdevStats->add_data(@meanStdevs);
@@ -112,10 +112,13 @@
 
 # rejections based on comparison with ensemble statistics
+## these should probably be 3-sigma clipped statistics themselves...
 my $mean = $meanStats->mean();
 my $meanStdev = $meanStats->standard_deviation();
 if (not defined $meanStdev) { $meanStdev = 0; }
 
-my $var = $variances->mean();
-my $varStdev = $variances->standard_deviation();
+my $var = $varianceStats->mean();
+my $varStdev = $varianceStats->standard_deviation();
+
+print "Ensemble mean $mean +/- $meanStdev, stdev " . sqrt($var) . " +/- " . sqrt($varStdev) . "\n\n";
 
 # Go through again to do rejection, and update the database for each exposure
@@ -143,9 +146,9 @@
 	print "no rejection for exposure mean\n";
     }
-    if (defined REJECT_STDEV->{$det_type} && ($varStdev >0)) {
-	my $nSigma = abs($var[$i] - $var) / $varStdev;
+    if (REJECT_STDEV->{$det_type} && ($varStdev > 0)) {
+	my $nSigma = abs($variances[$i] - $var) / $varStdev;
 	if ($nSigma > REJECT_STDEV->{$det_type}) {
 	    print "Rejecting $expTag based on outlier stdev: " .
-		sqrt($variances[$i]) " is $nSigma vs " . REJECT_STDEV->{$det_type} . "\n";
+		sqrt($variances[$i]) . " is $nSigma vs " . REJECT_STDEV->{$det_type} . "\n";
 	    $reject = 1;
 	    goto UPDATE;