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;
Index: /trunk/ippScripts/scripts/detrend_reject_imfile.pl
===================================================================
--- /trunk/ippScripts/scripts/detrend_reject_imfile.pl	(revision 9457)
+++ /trunk/ippScripts/scripts/detrend_reject_imfile.pl	(revision 9458)
@@ -64,5 +64,5 @@
 # This measures how much variation there is in each imfile
 use constant REJECT_IMFILE_STDEV => {
-    'bias' => 0,		# Should be fairly flat; some CRs
+    'bias' => 15,		# Should be fairly flat; some CRs
     'dark' => 0,		# Lots of CRs
     'shutter' => undef,		# Can be significant structure
@@ -82,5 +82,5 @@
 # This measures how much variation there is across the imfiles
 use constant REJECT_EXPOSURE_STDEV => {
-    'bias' => 0,		# Should be little variation between chips
+    'bias' => 15,		# Should be little variation between chips
     'dark' => 0,		# Could be some glow on some chips
     'shutter' => undef,		# Can be significant structure
