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

Timestamp:

May 9, 2008, 11:03:01 AM (18 years ago)

Author:

eugene

Message:

use clipmean and clipstdev to ensemble stats; clearer output messages

File:

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

trunk/ippScripts/scripts/detrend_reject_exp.pl

@@ -78 +78 @@
    [   
        #          KEYWORD                 STATISTIC          CHIPTOOL FLAG
-       { name => "bg",             type => "mean",  flag => "ensMeanMean",       dtype => "float" },
-       { name => "bg",             type => "stdev", flag => "ensMeanStdev",      dtype => "float" },
-       { name => "bg_mean_stdev",  type => "mean",  flag => "ensMeanStdevMean",  dtype => "float" },
-       { name => "bg_mean_stdev",  type => "stdev", flag => "ensMeanStdevStdev", dtype => "float" },
-       { name => "bg_stdev",       type => "mean",  flag => "ensStdevMean",      dtype => "float" },
-       { name => "bg_stdev",       type => "stdev", flag => "ensStdevStdev",     dtype => "float" },
-       { name => "bg_skewness",    type => "mean",  flag => "ensSkewness",       dtype => "float" },
-       { name => "bg_kurtosis",    type => "mean",  flag => "ensKurtosis",       dtype => "float" },
+       { name => "bg",             type => "clipmean",  flag => "ensMeanMean",       dtype => "float" },
+       { name => "bg",             type => "clipstdev", flag => "ensMeanStdev",      dtype => "float" },
+       { name => "bg_mean_stdev",  type => "clipmean",  flag => "ensMeanStdevMean",  dtype => "float" },
+       { name => "bg_mean_stdev",  type => "clipstdev", flag => "ensMeanStdevStdev", dtype => "float" },
+       { name => "bg_stdev",       type => "clipmean",  flag => "ensStdevMean",      dtype => "float" },
+       { name => "bg_stdev",       type => "clipstdev", flag => "ensStdevStdev",     dtype => "float" },
+       { name => "bg_skewness",    type => "clipmean",  flag => "ensSkewness",       dtype => "float" },
+       { name => "bg_kurtosis",    type => "clipmean",  flag => "ensKurtosis",       dtype => "float" },
 
    ];
@@ -129 +129 @@
 }
 
+# each image has a mean mean (average of means over all chips)
+# a standard deviation mean (average of standard deviations over all chips)
+# and a 
+
+
+
 # we use the statistics of the ensemble to accept/reject exposurs
-my $ensMeanMean       = $rejstats->value_for_flag ("ensMeanMean");      
-my $ensMeanStdev      = $rejstats->value_for_flag ("ensMeanStdev");     
-my $ensMeanStdevMean  = $rejstats->value_for_flag ("ensMeanStdevMean"); 
-my $ensMeanStdevStdev = $rejstats->value_for_flag ("ensMeanStdevStdev");
-my $ensStdevMean      = $rejstats->value_for_flag ("ensStdevMean");     
-my $ensStdevStdev     = $rejstats->value_for_flag ("ensStdevStdev");    
+my $ensMeanMean       = $rejstats->value_for_flag ("ensMeanMean");       # average of all exposure means (in turn averaged over all chips)
+my $ensMeanStdev      = $rejstats->value_for_flag ("ensMeanStdev");      # standard deviation of all exposure means 
+my $ensMeanStdevMean  = $rejstats->value_for_flag ("ensMeanStdevMean");  # average over all exposures of the stdev of the means for each chip
+my $ensMeanStdevStdev = $rejstats->value_for_flag ("ensMeanStdevStdev"); # standard deviation over all exposures of the stdev of the means for each chip
+my $ensStdevMean      = $rejstats->value_for_flag ("ensStdevMean");      # average over all exposures of the sum of the squares of the stdevs for each chip 
+my $ensStdevStdev     = $rejstats->value_for_flag ("ensStdevStdev");     # standard deviation over all exposures of the sum of the squares of the stdevs for each chip 
 
 $ipprc->define_camera($camera);
@@ -151 +157 @@
 $ipprc->outroot_prepare($outroot);
 
-my $logName = $outroot . "log"; # Name for log
+my $logName = "$outroot.log"; # Name for log
 
 my $logFile;
 unless ($no_op) {
     $logFile = $ipprc->file_create_append( $logName );
-    print $logFile "Ensemble mean $ensMeanMean +/- $ensMeanStdevMean, stdev $ensStdevMean\n\n";
+    print $logFile "Ensemble mean $ensMeanMean +/- $ensMeanStdev\n";
+    print $logFile "Ensemble stdev $ensStdevMean +/- $ensStdevStdev\n";
+    print $logFile "Ensemble mean rms (over imfiles) $ensMeanStdevMean +/- $ensMeanStdevStdev\n\n"; 
 }