Changeset 612 in Sophya for trunk/Poubelle/archTOI.old/gondolageom.cc

Timestamp:

Nov 22, 1999, 10:43:44 AM (26 years ago)

Author:

ansari

Message:

fin de fichier, start

File:

• trunk/Poubelle/archTOI.old/gondolageom.cc (modified) (5 diffs)

trunk/Poubelle/archTOI.old/gondolageom.cc

@@ -4 +4 @@
 #include <math.h>
 #include "gondolageom.h"
+
+
 extern "C" {
 #include "aa_hadec.h"
@@ -27 +29 @@
             
 
-GondolaGeom::GondolaGeom() {
+GondolaGeom::GondolaGeom()
+: fit(200, 2)
+{
   azPend = 0;
   angPend = 0;
@@ -51 +55 @@
   if (nstars<0) {
     nstars = 0;
-    saz=staz=st=st2=0;
+    fit.clear();
   }
   
@@ -61 +65 @@
   if (azCor - az0 >  180) azCor -= 360;
   if (azCor - az0 < -180) azCor += 360;
-  saz  += azCor;
-  staz += azCor*deltasn;
-  st   += deltasn;
-  st2  += deltasn*deltasn;
+  
+  fit.addData(deltasn, azCor);
 }
 
@@ -70 +72 @@
 int GondolaGeom::solveStars() {
   if (nstars<2) return -1;
-  staz /= nstars;
-  st   /= nstars;
-  saz  /= nstars;
-  st2  /= nstars;
-  double a = (staz - st*saz) / (st2 - st*st);
-  double b = saz - a*st;
-  
-  nstars = -1;
-  
-  azimut = b;
+  fit.doFit();
+  azimut = fit.value(0);
   
   if (azimut > 360) azimut -= 360;
   if (azimut < 0)   azimut += 360;
 
+  nstars = -1;
+  
   return 0;
 }