Changeset 13925

Timestamp:

Jun 20, 2007, 12:43:07 PM (19 years ago)

Author:

Paul Price

Message:

Adding code to parse fringe amplitudes and errors.

File:

• trunk/PS-IPP-Config/lib/PS/IPP/Metadata/Stats.pm (modified) (5 diffs)

trunk/PS-IPP-Config/lib/PS/IPP/Metadata/Stats.pm

@@ -1 +1 @@
 # Copyright (c) 2006  Paul Price, Joshua Hoblitt
 #
-# $Id: Stats.pm,v 1.13 2007-02-09 03:55:22 price Exp $
+# $Id: Stats.pm,v 1.14 2007-06-20 22:43:07 price Exp $
 
 package PS::IPP::Metadata::Stats;
@@ -24 +24 @@
 
 use base qw( Class::Accessor::Fast );
-__PACKAGE__->mk_accessors( qw( bg_mean bg_stdev bg_mean_stdev constants variables bg_data bg_stdev_data ) );
+__PACKAGE__->mk_accessors( qw( bg_mean bg_stdev bg_mean_stdev fringe fringe_err
+                               constants variables bg_data bg_stdev_data ) );
 
 #$::RD_TRACE = 1;
@@ -38 +39 @@
                  bg_stdev => undef, # Standard deviation of the mean backgrounds
                  bg_mean_stdev => undef, # Mean of the standard deviation of the backgrounds
+                 fringe => [],  # Fringe amplitudes
+                 fringe_err => [], # Fringe amplitude errors
                  constants => $constants, # Array of values that should be constant through the FPA
                  variables => $variables, # Array of values that may be variable through the FPA
                  bg_data => [], # Array of background values
-                 bg_stdev_data => [],# Array of background standard deviations
+                 bg_stdev_data => [], # Array of background standard deviations
+                 fringe_data => [], # Fringe amplitudes each component
+                 fringe_err_data => [], # Fringe errors for each component
                  data => {}        # The data
                  };
@@ -128 +133 @@
     } else {
         $self->bg_stdev(undef);
+    }
+
+    # Get fringe measurements
+    if (scalar @{$self->{fringe_data}} > 0) {
+        foreach my $array (@{$self->{fringe_data}}) {
+            my $stats = Statistics::Descriptive::Sparse->new(); # Statistics calculator
+            $stats->add_data(@$array);
+            push @{$self->{fringe}}, $stats->mean();
+        }
+    }
+    if (scalar @{$self->{fringe_err_data}} > 0) {
+        foreach my $array (@{$self->{fringe_err_data}}) {
+            my @variances;
+            foreach my $value (@$array) {
+                push @variances, $value**2;
+            }
+            my $stats = Statistics::Descriptive::Sparse->new(); # Statistics calculator
+            $stats->add_data(@variances);
+            push @{$self->{fringe_err}}, sqrt($stats->mean());
+        }
     }
 
@@ -207 +232 @@
                 push @{$self->{bg_data}}, $entry->{value};
             }
+        } elsif ($entry->{name} =~ /^FRINGE/) {
+            my ($num) = $entry->{name} =~ /_(\d+)$/; # Component number
+            if ($entry->{name} =~ /_ERR_/) {
+                my $arrayRef = ${$self->{fringe_err_data}}[$num];
+                unless (defined $arrayRef) {
+                    my @array;
+                    $arrayRef = \@array;
+                    ${$self->{fringe_err_data}}[$num] = $arrayRef;
+                }
+                push @$arrayRef, $entry->{value};
+            } else {
+                my $arrayRef = ${$self->{fringe_data}}[$num];
+                unless (defined $arrayRef) {
+                    my @array;
+                    $arrayRef = \@array;
+                    ${$self->{fringe_data}}[$num] = $arrayRef;
+                }
+                push @$arrayRef, $entry->{value};
+            }
         } else {
             $self->_check_values($entry);