Changeset 10459

Timestamp:

Dec 4, 2006, 4:12:42 PM (20 years ago)

Author:

Paul Price

Message:

Rejection on mean stdev is relative to the mean (testing relative variation between components).

File:

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

trunk/ippScripts/scripts/detrend_reject_exp.pl

@@ -59 +59 @@
 my @means;                      # Array of means
 my @stdevs;                     # Array of standard deviations
-my @meanStdevs;                 # Array of standard deviations of the mean
+my @meanStdevs;                 # Array of standard deviations of the mean (normalised by the mean)
 my @variances;                  # Array of variances
 my @accept;                     # Array of accept flags
@@ -73 +73 @@
     push @means, $exposure->{bg};
     push @stdevs, $exposure->{bg_stdev};
-    push @meanStdevs, $exposure->{bg_mean_stdev};
+    push @meanStdevs, $exposure->{bg_mean_stdev} / $exposure->{bg};
     push @variances, $exposure->{bg_stdev}**2;
     push @accept, $exposure->{accept};
@@ -138 +138 @@
         if ($nSigma > $reject_meanstdev) {
             print "Rejecting $expTag based on outlier mean stdev: " .
-                $meanStdevs[$i] . " --> $nSigma sigma vs " . $reject_stdev . "\n";
+                $meanStdevs[$i] . " --> $nSigma sigma vs " . $reject_meanstdev . "\n";
             $reject = 1;
             goto UPDATE;