Changeset 10538 for trunk/ippScripts/scripts/detrend_reject_imfile.pl

Timestamp:

Dec 7, 2006, 2:23:25 PM (19 years ago)

Author:

Paul Price

Message:

Adding rejection on the imfile mean stdevs.

File:

• trunk/ippScripts/scripts/detrend_reject_imfile.pl (modified) (4 diffs)

trunk/ippScripts/scripts/detrend_reject_imfile.pl

@@ -77 +77 @@
 my @means;                      # Array of means
 my @variances;                  # Array of variances
+my @meanStdevs;                 # Array of mean stdevs
 open my $list1File, '>' . $list1Name;
 open my $list2File, '>' . $list2Name;
@@ -83 +84 @@
     print $list2File ($ipprc->convert_filename_absolute( $file->{b2_uri} ) . "\n");
     push @means, $file->{bg};
+    push @meanStdevs, $file->{bg_mean_stdev};
     ## calculate the root-mean-square of the bd_stdevs
     push @variances, $file->{bg_stdev}*$file->{bg_stdev};
+
 }
 close $list1File;
@@ -111 +114 @@
 my $expected = rejection_limit( 'EXPECTED', $det_type, $filter ); # Expected mean
 # Rejection thresholds
-my $reject_imfile_mean   = rejection_limit( 'IMFILE.MEAN',   $det_type, $filter );
-my $reject_imfile_stdev  = rejection_limit( 'IMFILE.STDEV',  $det_type, $filter );
-my $reject_imfile_sn     = rejection_limit( 'IMFILE.SN',     $det_type, $filter );
-my $reject_exp_mean      = rejection_limit( 'EXP.MEAN',      $det_type, $filter );
-my $reject_exp_stdev     = rejection_limit( 'EXP.STDEV',     $det_type, $filter );
-my $reject_exp_meanstdev = rejection_limit( 'EXP.MEANSTDEV', $det_type, $filter );
-my $reject_exp_sn        = rejection_limit( 'EXP.SN',        $det_type, $filter );
+my $reject_imfile_mean      = rejection_limit( 'IMFILE.MEAN',      $det_type, $filter );
+my $reject_imfile_stdev     = rejection_limit( 'IMFILE.STDEV',     $det_type, $filter );
+my $reject_imfile_meanstdev = rejection_limit( 'IMFILE.MEANSTDEV', $det_type, $filter );
+my $reject_imfile_sn        = rejection_limit( 'IMFILE.SN',        $det_type, $filter );
+my $reject_exp_mean         = rejection_limit( 'EXP.MEAN',         $det_type, $filter );
+my $reject_exp_stdev        = rejection_limit( 'EXP.STDEV',        $det_type, $filter );
+my $reject_exp_meanstdev    = rejection_limit( 'EXP.MEANSTDEV',    $det_type, $filter );
+my $reject_exp_sn           = rejection_limit( 'EXP.SN',           $det_type, $filter );
 
 # Reject based on the stats of the imfiles
@@ -146 +150 @@
     } else {
         print "no rejection for imfile stdev\n";
+    }
+    if ($reject_imfile_meanstdev) {
+        if ($meanStdevs[$i] > $reject_imfile_meanstdev) {
+            print "Rejecting exposure based on bad imfile mean stdev for imfile $i: " .
+                $meanStdevs[$i] . " vs " . $reject_imfile_meanstdev . "\n";
+            $reject = 1;
+            last;
+        }
+    } else {
+        print "no rejection for imfile mean stdev\n";
     }
     if ($reject_imfile_sn) {