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

Timestamp:

Jun 4, 2007, 2:25:37 PM (19 years ago)

Author:

eugene

Message:

calculate and supply the binned stdev

File:

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

trunk/ippScripts/scripts/detrend_reject_imfile.pl

@@ -120 +120 @@
 my @means;                      # Array of means
 my @variances;                  # Array of variances
+my @binVariances;               # Array of binned variances
 my @meanStdevs;                 # Array of mean stdevs
 my @names;                      # Array of names (class_id)
@@ -129 +130 @@
     ## calculate the root-mean-square of the bd_stdevs
     push @variances, $file->{bg_stdev}*$file->{bg_stdev};
+    push @binVariances, $file->{bin_stdev}*$file->{bin_stdev};
     push @names, $file->{class_id};
 }
@@ -180 +182 @@
 my $reject_imfile_meanstdev = rejection_limit( 'IMFILE.MEANSTDEV', $det_type, $filter );
 my $reject_imfile_sn        = rejection_limit( 'IMFILE.SN',        $det_type, $filter );
+my $reject_imfile_bin_sn    = rejection_limit( 'IMFILE.BIN.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_exp_bin_sn       = rejection_limit( 'EXP.BIN.SN',       $det_type, $filter );
 
 my $logName = caturi( $workdir, "$exp_tag.detreject.$det_id.$iter.log" ); # Name for log
@@ -199 +203 @@
     $mean -= $expected;
     my $stdev = sqrt($variances[$i]);   # Stdev for this imfile
+    my $binStdev = sqrt($binVariances[$i]);     # Stdev for this imfile
     my $name = $names[$i];
 
@@ -247 +252 @@
     } else {
         print $logFile "No rejection on imfile S/N for $name\n";
+    }
+    if ($reject_imfile_bin_sn) {
+        if ($mean < $binStdev * $reject_imfile_bin_sn) {
+            print $logFile "Rejecting exposure based on bad imfile binned S/N for $name: ";
+            $reject = 1;
+        } else {
+            print $logFile "Imfile binned S/N for $name meets requirements: ";
+        }
+        print $logFile "mean: $mean vs binStdev*SNlimit: " . $binStdev * $reject_imfile_bin_sn . "\n";
+    } else {
+        print $logFile "No rejection on imfile binned S/N for $name\n";
     }
 }
@@ -255 +271 @@
 my $varianceStats = Statistics::Descriptive::Sparse->new(); # Statistics calculator for variances
 $varianceStats->add_data(@variances);
+my $binVarianceStats = Statistics::Descriptive::Sparse->new(); # Statistics calculator for variances
+$binVarianceStats->add_data(@binVariances);
 
 my $mean = $meanStats->mean();  # Mean of the imfile means
@@ -263 +281 @@
 }
 my $stdev = sqrt($varianceStats->mean()); # Root-Mean-Square of the imfile stdevs (root mean of variances)
+my $binStdev = sqrt($binVarianceStats->mean()); # Root-Mean-Square of the imfile stdevs (root mean of variances)
 my $exp_sn = 0.0;
 if ($stdev > 0) { $exp_sn = $mean / $stdev; }
@@ -318 +337 @@
         print $logFile "No rejection for exp S/N\n";
     }
+    # reject if the signal-to-noise is insufficient
+    if ($reject_exp_bin_sn) {
+        if (abs($mean) < abs($binStdev * $reject_exp_bin_sn)) {
+            print $logFile "Rejecting exposure based on poor binned S/N: \n";
+            $reject = 1;
+        } else {
+            print $logFile "Exposure binned S/N meets requirements: \n";
+        }
+        print $logFile "signal: $mean vs noise: $binStdev (s/n limit is: $reject_exp_bin_sn)\n";
+    } else {
+        print $logFile "No rejection for exp binned S/N\n";
+    }
 
     close $logFile;
@@ -326 +357 @@
 my $bg_stdev = $stdev;
 my $bg_mean_stdev = $meanStdev;
+my $bin_stdev = $binStdev;
 
 unless ($no_update) {
     my $command = "$dettool -addresidexp -det_id $det_id -iteration $iter -exp_tag $exp_tag";
     $command .= " -recip $recipe1,$recipe2 -path_base $outputRoot ";
-    $command .= " -bg $bg -bg_stdev $bg_stdev -bg_mean_stdev $bg_mean_stdev";
+    $command .= " -bg $bg -bg_stdev $bg_stdev -bg_mean_stdev $bg_mean_stdev -bin_stdev $bin_stdev";
     $command .= ' -reject' if $reject;
     $command .= " -dbname $dbname" if defined $dbname;