Changeset 92

Timestamp:

Jul 23, 2003, 3:36:00 PM (23 years ago)

Author:

eugene

Message:

added abs to sqrt to clear rounding error problems

File:

• trunk/Ohana/src/perl/src/mktrans (modified) (3 diffs)

trunk/Ohana/src/perl/src/mktrans

@@ -10 +10 @@
 foreach $code (@photcodes) {
 
+    print STDERR "filter: $code\n";
     @list = `photsearch -trange $start $stop -photcode $code`;
 
     # the lines from photsearch look like this:
-    # V 2001/02/03,15:38:53  26.249  0.017 elixir
+    # filter date            zp      dzp   N T label (N - Nmeas, T - Ntime)
+    # V 2001/02/03,15:38:53  26.249  0.017 3 1 elixir
 
     @dates = ();
@@ -29 +31 @@
         # add to the stats 
         $sumN{$date} ++;
+        $sumM{$date} += $words[5];
         $sum1{$date} += $words[2];
         $sum2{$date} += $words[2]*$words[2];
-
-        # print STDERR "date: $date\n";
-
     }
 
@@ -41 +41 @@
     foreach $date (@sdates) {
         
+        $Nm  = $sumM{$date};
         $Np  = $sumN{$date};
         $zp  = $sum1{$date} / $Np;
-        $dzp = sqrt ($sum2{$date} / $Np - $zp*$zp);
+        $dzp = sqrt (abs($sum2{$date} / $Np - $zp*$zp));
 
-        printf STDOUT "transreg %s %6.3f %6.4f %s 1d\n", $code, $zp, $dzp, $date;
+        printf STDOUT "photreg -trans -photcode $code -zp %6.3f -dzp %6.4f -trange $date 1d -Nmeas $Nm -Ntime $Np\n", $zp, $dzp;
     }
 
+    %sumM = ();
     %sumN = ();
     %sum1 = ();