Changeset 423 in Sophya for trunk/Poubelle/archTOI.old/templocator.cc

Timestamp:

Sep 23, 1999, 4:26:26 PM (26 years ago)

Author:

ansari

Message:

galcross

File:

• trunk/Poubelle/archTOI.old/templocator.cc (modified) (2 diffs)

trunk/Poubelle/archTOI.old/templocator.cc

@@ -4 +4 @@
 #include "aa_hadec.h"
 }
+
+#include "fitsio.h"
+#include "plgalcross.h"
 
 #ifndef M_PI
@@ -14 +17 @@
 {
   lon = lat = ts = 0;
+  ra  = dec = -99999;
+  xSampleNum = -99999;
+  fitsfile* fptr;
+  int status=0;
+  fits_open_file(&fptr, "samplenum_gal_cross.fits", READONLY, &status);
+  int simple, bitpix, naxis;
+  long naxes;
+  long pcount, gcount;
+  int extend;
+  fits_read_imghdr(fptr, 1, &simple, &bitpix, &naxis, &naxes, &pcount, &gcount, &extend, &status);
+  nGalCross = naxes;
+  crossings = new long[nGalCross];
+  int anynul;
+  fits_read_img_lng(fptr, 0, 1, nGalCross, 0, crossings, &anynul, &status);
+  fits_close_file(fptr, &status);
+  fits_report_error(stderr, status);  /* print out any error messages */
+}
+           
+           
+void TempLocator::setEarthPos(double lon, double lat) {
+  if (this->lon == lon && this->lat == lat) return;
+  this->lon = lon;
+  this->lat = lat;
+  ra  = dec = -99999; xSampleNum = -99999;
 }
 
-void TempLocator::setEarthPos(double lon, double lat) {
-  this->lon = lon;
-  this->lat = lat;
+void TempLocator::setTSid(double ts) {
+  if (this->ts == ts) return;
+  this->ts  = ts;
+  ra  = dec = -99999; xSampleNum = -99999;
 }
-void TempLocator::setTSid(double ts) {
-  this->ts  = ts;
+  
+void   TempLocator::ComputeZenith() {
+  double ha;
+  aa_hadec (lat * M_PI/180, .5 * M_PI, 0, &ha, &dec);
+  ra = - (ha * 180. / M_PI / 15) + (ts/3600.);
+  dec = dec * 180. / M_PI;
 }
   
 double TempLocator::getAlphaZenith() {
-  double ha, dec;
-  aa_hadec (lat * M_PI/180, .5 * M_PI, 0, &ha, &dec);
-  double ra = - (ha * 180. / M_PI / 15) + (ts/3600.);
-  //dec = dec * 180. / M_PI;
+  if (ra < -100) ComputeZenith();
   return ra;
 }
 
 double TempLocator::getDeltaZenith() {
-  double ha, dec;
-  aa_hadec (lat * M_PI/180, .5 * M_PI, 0, &ha, &dec);
-  //double ra = - (ha * 180. / M_PI / 15) + (ts/3600.);
-  dec = dec * 180. / M_PI;
+  if (dec < -100) ComputeZenith();
   return dec;
 }
 
+void TempLocator::findGeomFromGC(int sampleNum) // pour le bolo qui voit les xing
+{
+  if (sampleNum == xSampleNum) return;
+  if (dec < -100) ComputeZenith();
+  
+  azimBolGC = -9999;
+  
+  // On trouve les croisements juste avant et juste apres notre sampleNum
+  int icross;
+  for (icross=0; icross<nGalCross; icross++) {
+    if (crossings[icross] > sampleNum) break;
+  }
+  if (icross == 0 || icross >= nGalCross) return;
+  
+  
+}
 
+
+                        //       ;
+                        //       ;                            /\         4     2
+                        //       ;                  elevation ||                       positive scanning, clockwise
+                        //       ;                            ||       6    1    5      ---------->
+                        //       ;                            ||                      == positive azimut
+                        //       ;                            ||         x     3
+            //                                         |----|      = 0.78 deg / cos(elev)
+            // bol 1 = 41° elevation
+