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

Timestamp:

Jun 25, 2007, 2:34:23 PM (19 years ago)

Author:

eugene

Message:

passing fringe amplitude and residual from imfiles to exp

File:

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

trunk/ippScripts/scripts/detrend_reject_imfile.pl

@@ -212 +212 @@
 my @meanStdevs;                 # Array of mean stdevs
 my @names;                      # Array of names (class_id)
+my @fringe_means;
+my @fringe_vars;
+my @dfringe_means;
+my @dfringe_vars;
 
 # load the arrays from the imfile output table
@@ -221 +225 @@
         push @fluxes,           $file->{bg};
     }
-    push @meanStdevs,   $file->{bg_mean_stdev};                # stdev of the mean counts (for imfile components)
-    push @variances,    $file->{bg_stdev}*$file->{bg_stdev};   # total variance from all imfile components
-    push @binVariances, $file->{bin_stdev}*$file->{bin_stdev}; # total variance of binned images over all imfile components
-    push @names,        $file->{class_id};                     # name of the component
+    push @meanStdevs,        $file->{bg_mean_stdev};                # stdev of the mean counts (for imfile components)
+    push @variances,         $file->{bg_stdev}*$file->{bg_stdev};   # total variance from all imfile components
+    push @binVariances,      $file->{bin_stdev}*$file->{bin_stdev}; # total variance of binned images over all imfile components
+    push @names,             $file->{class_id};                     # name of the component
+    push @fringe_means,      $file->{fringe_0};                     # fringe amplitude mean for imfile
+    push @fringe_vars,       $file->{fringe_1}*$file->{fringe_1};   # fringe variance
+    push @dfringe_means,     $file->{user_1};                       # fringe residual mean for imfile
+    push @dfringe_vars,      $file->{user_2}*$file->{user_2};       # fringe residual variance
+    # push @fringe_mean_stdev, $file->{fringe_2};                     # fringe amplitude stdev for imfile
 }
 
@@ -357 +366 @@
 $binVarianceStats->add_data(@binVariances);
 
+# background stats
 my $mean = $meanStats->mean();  # Mean of the imfile means
 my $flux = $fluxStats->mean();  # Mean of the imfile means
@@ -368 +378 @@
 my $exp_sn = 0.0;
 if ($stdev > 0) { $exp_sn = $mean / $stdev; }
+
+# prepare fringe amplitude stats
+my $fringeMeanStats = Statistics::Descriptive::Sparse->new();   # Statistics calculator for means
+$fringeMeanStats->add_data(@fringe_means);
+my $fringeVarStats = Statistics::Descriptive::Sparse->new();    # Statistics calculator for means
+$fringeVarStats->add_data(@fringe_vars);
+
+# fringe amplitude stats
+my $fringe_mean = $fringeMeanStats->mean();
+my $fringe_err = sqrt($fringeVarStats->mean());
+my $fringe_mean_stdev = $fringeMeanStats->standard_deviation();
+if (not defined $fringe_mean_stdev) {
+    $fringe_mean_stdev = 0;
+}
+
+# prepare fringe residual stast 
+my $dfringeMeanStats = Statistics::Descriptive::Sparse->new();  # Statistics calculator for means
+$dfringeMeanStats->add_data(@dfringe_means);
+my $dfringeVarStats = Statistics::Descriptive::Sparse->new();   # Statistics calculator for means
+$dfringeVarStats->add_data(@dfringe_vars);
+
+# fringe amplitude stats
+my $dfringe_mean = $dfringeMeanStats->mean();
+my $dfringe_err = sqrt($dfringeVarStats->mean());
+my $dfringe_mean_stdev = $dfringeMeanStats->standard_deviation();
+if (not defined $dfringe_mean_stdev) {
+    $dfringe_mean_stdev = 0;
+}
 
 ## Reject based on the exposure ensemble stats
@@ -486 +524 @@
     $command .= " -recip $recipe1,$recipe2 -path_base $outputRoot ";
     $command .= " -bg $bg -bg_stdev $bg_stdev -bg_mean_stdev $bg_mean_stdev -bin_stdev $bin_stdev";
+    $command .= " -fringe_0 $fringe_mean -fringe_1 $fringe_err -fringe_2 $fringe_mean_stdev";
+    $command .= " -user_1 $dfringe_mean -user_2 $dfringe_err -user_3 $dfringe_mean_stdev";
     $command .= ' -reject' if $reject;
     $command .= " -dbname $dbname" if defined $dbname;