Changeset 12673 for trunk/ippScripts/scripts/detrend_reject_exp.pl

Timestamp:

Mar 29, 2007, 10:25:51 AM (19 years ago)

Author:

Paul Price

Message:

Write determinations to log file.

File:

• trunk/ippScripts/scripts/detrend_reject_exp.pl (modified) (7 diffs)

trunk/ippScripts/scripts/detrend_reject_exp.pl

@@ -121 +121 @@
 my $reject_meanstdev = rejection_limit( 'ENSEMBLE.MEANSTDEV', $det_type, $filter );
 
-print "Ensemble mean " . $meanStats->mean() . " +/- " . $meanStats->standard_deviation .
-    ", stdev " . $stdevStats->mean() . " +/- " . $stdevStats->standard_deviation() . "\n\n";
+$workdir = $ipprc->convert_filename_absolute( $workdir );
+
+my $logName = "$camera.$det_type.$det_id.$iter.detreject.log"; # Name for log
+$logName = File::Spec->catfile( $workdir, $logName );
+my $logFile;
+unless ($no_op) {
+    open $logFile, "> $logName" or &my_die("Unable to open log file $logName.\n", $det_id, $iter, $PS_EXIT_SYS_ERROR);
+    print $logFile "Ensemble mean " . $meanStats->mean() . " +/- " . $meanStats->standard_deviation .
+        ", stdev " . $stdevStats->mean() . " +/- " . $stdevStats->standard_deviation() . "\n\n";
+}
 
 # Go through again to do rejection, and update the database for each exposure
@@ -139 +147 @@
     if (not $accept[$i]) {
         # Rejected this at an earlier stage
-        print "Rejecting $expTag based on earlier determination.\n";
+        unless ($no_op) {
+            print $logFile "Rejecting $expTag based on earlier determination.\n";
+        }
         $reject = 1;
         goto UPDATE;
@@ -157 +167 @@
         my $dMean = abs($means[$i] - $meanStats->mean()) ;
         if ($dMean > ($reject_mean * $meanStats->standard_deviation())) {
-            print "Rejecting $expTag based on outlier mean value: " .
+            print $logFile "Rejecting $expTag based on ensemble mean value: " .
                 "$means[$i] --> $dMean vs " . $reject_mean * $meanStats->standard_deviation() . "\n";
             $reject = 1;
-            goto UPDATE;
+            #goto UPDATE;
         }
     } else {
-        print "no rejection for exposure mean\n";
+        print $logFile "No rejection for ensemble mean\n";
     }
 
@@ -169 +179 @@
         my $dMean = abs($stdevs[$i] - $stdevStats->mean());
         if ($dMean > ($reject_stdev * $stdevStats->standard_deviation())) {
-            print "Rejecting $expTag based on outlier stdev: " .
+            print $logFile "Rejecting $expTag based on ensemble stdev: " .
                 $stdevs[$i] . " --> $dMean sigma vs " . $reject_stdev . "\n";
             $reject = 1;
-            goto UPDATE;
+            #goto UPDATE;
         }
     } else {
-        print "no rejection for exposure stdev\n";
+        print $logFile "No rejection for ensemble stdev\n";
     }
     
@@ -181 +191 @@
         my $dMean = abs($meanStdevs[$i] - $meanStdevStats->mean());
         if ($dMean > ($reject_meanstdev * $meanStdevStats->standard_deviation())) {
-            print "Rejecting $expTag based on outlier mean stdev: " .
+            print $logFile "Rejecting $expTag based on ensemble mean stdev: " .
                 $meanStdevs[$i] . " --> $dMean sigma vs " . $reject_meanstdev . "\n";
             $reject = 1;
-            goto UPDATE;
+            #goto UPDATE;
         }
     } else {
-        print "no rejection for exposure mean stdev\n";
+        print $logFile "No rejection for ensemble mean stdev\n";
     }
 
@@ -198 +208 @@
     # Check for status changes
     if ((not $include[$i] and not $reject) or ($include[$i] and $reject)) {
-        print "Status of $expTag has changed.\n";
+        unless ($no_op) {
+            print $logFile "Status of $expTag has changed.\n";
+        }
         $numChanges++;
     }
@@ -228 +240 @@
 }
 
-print "Master: $master\n";
-print "Stop: $stop\n";
+unless ($no_op) {
+    print $logFile "Master: $master\n";
+    print $logFile "Stop: $stop\n";
+    close $logFile;
+}
 
 # Put results into the database