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

Timestamp:

Sep 24, 1999, 9:15:56 AM (26 years ago)

Author:

ansari

Message:

alpha delta a partir croisement galax

File:

• trunk/Poubelle/archTOI.old/templocator.cc (modified) (1 diff)

trunk/Poubelle/archTOI.old/templocator.cc

@@ -85 +85 @@
   double azCr1, azCr2;
   int rc = PlGalCross(ts/3600., lat, (90. - 41.), alpG, delG, azCr1, azCr2);
+  if (rc != 0) return; // pas deux points d'intersection
   
+  // Il faut determiner le croisement principal, ie le plus proche
+  // du centre galactique. Pendant le vol de Trapani, c'etait celui
+  // le plus proche de 220° d'azimut.
+  
+  double azCross = azCr1;
+  if (fabs(azCr2-220) < fabs(azCr1-220)) azCross = azCr2;
+  
+  double rotSpeed = 360./(crossings[icross] - crossings[icross-1]); // °/sample
+  
+  azimBolGC = azCross + (sampleNum - crossings[icross-1])*rotSpeed;
+  if (azimBolGC > 360) azimBolGC -= 360;
 }
+
+void TempLocator::getAltAzBolo(int sampleNum, int ibolo, double& elv, double& az) {
+  findGeomFromGC(sampleNum);
+  double delElv = 0;
+  double delAz = 0;  // relative to ch3-bolo6
+  elv = -99999;
+  az  = -99999;
+  switch (ibolo) {
+    case 11:
+      delElv =   0;
+      delAz  =   0.78;
+      break;
+    case 8:
+      delElv =   0.78 * sqrt(3.)/2.;
+      delAz  =   0.78 * 3./2.;
+      break;
+    case 13:
+      delElv = - 0.78 * sqrt(3.)/2.;
+      delAz  =   0.78 * 3./2.;
+      break;
+    case 9:
+      delElv =   0.78 * sqrt(3.)/2.;
+      delAz  =   0.78 * 1./2.;
+      break;
+    case 4:
+      delElv =   0.;
+      delAz  =   0.78 * 2.;
+      break;
+    case 15:
+      delElv =   0.;
+      delAz  =   0.;
+      break;
+    default:
+      return;
+  }
+  delAz /= cos(41 * M_PI/180);
+  elv = 41. + delElv;
+  az  = azimBolGC + delAz;
+  return;
+}
+
+double TempLocator::getAzimutBolo(int sampleNum, int ibolo) {
+  double elv, az;
+  getAltAzBolo(sampleNum, ibolo, elv, az);
+  return az;
+}
+
+double TempLocator::getElvBolo(int sampleNum, int ibolo) {
+  double elv, az;
+  getAltAzBolo(sampleNum, ibolo, elv, az);
+  return elv;
+}
+
+double TempLocator::getAlphaBolo(int sampleNum, int ibolo) {
+  double elv, az;
+  getAltAzBolo(sampleNum, ibolo, elv, az);
+  double ha;
+  aa_hadec (lat * M_PI/180, elv * M_PI/180, az * M_PI/180, &ha, &dec);
+  ra = - (ha * 180. / M_PI / 15) + (ts/3600.);
+  dec = dec * 180. / M_PI;
+  return ra;
+}
+
+double TempLocator::getDeltaBolo(int sampleNum, int ibolo) {
+  double elv, az;
+  getAltAzBolo(sampleNum, ibolo, elv, az);
+  double ha;
+  aa_hadec (lat * M_PI/180, elv * M_PI/180, az * M_PI/180, &ha, &dec);
+  ra = - (ha * 180. / M_PI / 15) + (ts/3600.);
+  dec = dec * 180. / M_PI;
+  return dec;
+}
+