Changeset 581 in Sophya for trunk/Poubelle/archTOI.old/starmatcher.cc

Timestamp:

Nov 16, 1999, 5:10:08 PM (26 years ago)

Author:

ansari

Message:

sst

File:

• trunk/Poubelle/archTOI.old/starmatcher.cc (modified) (7 diffs)

trunk/Poubelle/archTOI.old/starmatcher.cc

@@ -12 +12 @@
 #define STARDUMP
 
+#include <math.h>
+
 extern "C" {
 #include "aa_hadec.h"
 #define NRANSI
 #include "nrutil.h"
+
+#ifndef M_PI
+#define M_PI 3.1415926535
+#endif
 
 void lfit(double x[], double y[], double sig[], int ndat, double a[], int ia[],
@@ -235 +241 @@
     double dmin = delta - dist; if (dmin<-90) dmin=-90;
     double dmax = delta + dist; if (dmax> 90) dmax= 90;
-    double amin = alpha - dist / cos(delta * 3.1415926/180) / 15.;
+    double amin = alpha - dist / cos(delta * M_PI/180) / 15.;
     if (amin<0) amin += 24;
-    double amax = alpha + dist / cos(delta * 3.1415926/180) / 15.;
+    double amax = alpha + dist / cos(delta * M_PI/180) / 15.;
     if (amax>24) amax -= 24;
   
@@ -256 +262 @@
     for (int i=imin; i<=imax; i++) {
       if (stars[i].ra >= amin && stars[i].ra <= amax) {
-        double ha = (ts/3600. - stars[i].ra) * 15. * 3.1415926/180.;
+        double ha = (ts/3600. - stars[i].ra) * 15. * M_PI/180.;
         double elv, azim;
-        hadec_aa(lat * 3.1415926/180., ha, stars[i].dec * 3.1415926/180.,
+        hadec_aa(lat * M_PI/180., ha, stars[i].dec * M_PI/180.,
                  &elv, &azim);
-        elv  *= 180/3.1415926;
-        azim *= 180/3.1415926;
+        elv  *= 180/M_PI;
+        azim *= 180/M_PI;
         if (azim<0) azim += 360;
       
         double da = azim-az; if (da>360) da -= 360;
         if (da < -0.6 || da > 0.4) continue;
-        double elv0 = elv - 1.41/45. * lastStar.NoDiode;
+        double elv0 = elv - GondolaGeom::sstPixelHeight * lastStar.NoDiode;
         if (fabs(elv0-GondolaGeom::elevSST0) > 0.25) continue; // Might be too strong
         
@@ -375 +381 @@
     if (ndata > 0 && az - oldAz >  180) az -= 360;
     if (ndata > 0 && az - oldAz < -180) az += 360;
-    fitElvAz.addData(s1.SN-sn0, s1.elvGSC - s1.nDiode*1.41/45., az);
+    fitElvAz.addData(s1.SN-sn0, s1.elvGSC - s1.nDiode*GondolaGeom::sstPixelHeight, az);
     oldAz = az;
     ndata++;
@@ -408 +414 @@
       matchStar sss = (*it1);
       if (!sss.ok) continue;
-      if (fabs(fitElvAz.valueY(sss.SN-sn0)-(sss.elvGSC - sss.nDiode*1.41/45.)) > dcutElv) {
+      if (fabs(fitElvAz.valueY(sss.SN-sn0)-
+               (sss.elvGSC - sss.nDiode*GondolaGeom::sstPixelHeight)) > dcutElv) {
         (*it1).ok = false;
         //logstream << "  kill " << sss.seq << " " << setprecision(11) << sss.SN << " "
@@ -483 +490 @@
     ndata++;
     snmean += st.SN;
-    azi[ndata]  = st.azGSC * 3.1415926/180;
-    elv0[ndata] = st.elvGSC - st.nDiode*1.41/45.;
+    azi[ndata]  = st.azGSC * M_PI/180;
+    elv0[ndata] = st.elvGSC - st.nDiode*GondolaGeom::sstPixelHeight;
     sig[ndata]  = 0.01;
   }
@@ -516 +523 @@
 
   double ampl = sqrt(cc*cc+ss*ss);
-  double phase = atan2(cc,ss)/(3.1415926/180);
+  double phase = atan2(cc,ss)/(M_PI/180);