Changeset 9524 for trunk/ippScripts/scripts/phase0_exp.pl

Timestamp:

Oct 12, 2006, 2:00:09 PM (20 years ago)

Author:

Paul Price

Message:

Using correct statistics: mean is mean of the background, stdev is the global stdev (\sum s_i^{2) / N, mean stdev is the stdev of the means.}

File:

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

trunk/ippScripts/scripts/phase0_exp.pl

@@ -100 +100 @@
 }
 my @backgrounds;                # Array of backgrounds
-my @stdevs;                        # Array of background standard deviations\
+my @variances;                        # Array of background standard deviations\
 my $obsType;                        # Observation type
   
@@ -124 +124 @@
     
     my $stdev = get_value($imfile, PHASE0_BG_MEAN_STDEV());
-    push @stdevs, $stdev;
+    push @variances, $stdev**2;
 }
 
@@ -147 +147 @@
         my $stats = Statistics::Descriptive::Sparse->new; # Statistics calculator
         $stats->add_data(@backgrounds);
-        push @command, ( '-' . PHASE0_BG() ), $stats->mean(), ( '-' . PHASE0_BG_STDEV() );
-        if (scalar @backgrounds == 1) {
-            push @command, 1;
-        } else {
-            push @command, $stats->standard_deviation();
-        }
+        push @command, ( '-' . PHASE0_BG() ), $stats->mean();
+        push @command, ( '-' . PHASE0_BG_MEAN_STDEV() ), ( $stats->standard_deviation() || 0 )
     }
     {
         my $stats = Statistics::Descriptive::Sparse->new; # Statistics calculator
-        $stats->add_data(@stdevs);
-        push @command, ( '-' . PHASE0_BG_MEAN_STDEV() ), $stats->mean();
+        $stats->add_data(@variances);
+        push @command, ( '-' . PHASE0_BG_STDEV() ), sqrt( $stats->mean() );
     }