Changeset 1148 in Sophya for trunk/SigPredictor/sigcalchackingtools.cc

Timestamp:

Aug 29, 2000, 12:03:24 PM (25 years ago)

Author:

ansari

Message:

mise a jour

File:

• trunk/SigPredictor/sigcalchackingtools.cc (modified) (2 diffs)

trunk/SigPredictor/sigcalchackingtools.cc

@@ -5 +5 @@
 #include <fstream>
 #ifdef __MWERKS__
-   #include "mwerksmath.h"
+    
    #include "unixmac.h"
 #endif
@@ -12 +12 @@
 //#include "integ.h"
  
+ int kmg_eulerRad(double ai, double bi, int select, double *ao, double *bo) {
+  // All coordinates are in Radian.
+  /* kmg_euler.c
+   *
+   * Converts between different coordinate systems. 
+   *
+   * SELECT   From           To         |   SELECT    From          To
+   *    1     RA-Dec (2000)  Galactic   |     4       Ecliptic      RA-Dec    
+   *    2     Galactic       RA-DEC     |     5       Ecliptic
+Galactic  
+   *    3     RA-Dec         Ecliptic   |     6       Galactic
+Ecliptic  
+   *
+   *    Date     Programmer Remarks
+   * ----------- ---------- -------
+   * 08-Aug-1999 K. Ganga   First version. Copied and modified EULER from 
+   *                        the IDL Astrolib. 
+   * May 2000, D. Yvon          Change coordinates units to Radians
+   */
+
+  /* Local Declarations */
+  double       a, b, sb, cb, cbsa;
+  long         i;
+  const double TWOPI   = 2.0*M_PI;
+  const double FOURPI  = 4.0*M_PI;
+  const double DEG2RAD = 180.0/M_PI;
+
+  /* J2000 coordinate conversions are based on the following constants
+   * eps    =  23.4392911111  Obliquity of the ecliptic
+   * alphaG = 192.85948       Right Ascension of Galactic North Pole
+   * deltaG =  27.12825       Declination of Galactic North Pole
+   * lomega =  32.93192       Galactic longitude of celestial equator
+   * alphaE = 180.02322       Ecliptic longitude of Galactic North Pole
+   * deltaE =  29.811438523   Ecliptic latitude of Galactic North Pole
+   * Eomega =   6.3839743     Galactic longitude of ecliptic equator
+   */
+
+  const double psi[6]    = {0.57477043300, 4.9368292465 ,
+                            0.00000000000, 0.0000000000 , 
+                            0.11142137093, 4.71279419371};
+  const double stheta[6] = {0.88998808748,-0.88998808748,
+                            0.39777715593,-0.39777715593,
+                            0.86766622025,-0.86766622025};    
+  const double ctheta[6] = {0.45598377618, 0.45598377618, 
+                            0.91748206207, 0.91748206207,
+                            0.49714719172, 0.49714719172};
+  const double phi[6]    = {4.9368292465 , 0.57477043300,
+                            0.0000000000 , 0.00000000000, 
+                            4.71279419371, 0.11142137093};
+
+  i    = select - 1;
+  a    = ai - phi[i];
+  b    = bi;
+  sb   = sin(b);
+  cb   = cos(b);
+  cbsa = cb*sin(a);
+  b    = -stheta[i]*cbsa + ctheta[i]*sb;
+  b    = ( (b >  1.0) ?  1.0 : b );
+  b    = ( (b < -1.0) ? -1.0 : b );
+  *bo  = asin(b);
+
+  a    = atan2( ctheta[i] * cbsa + stheta[i] * sb, cb * cos(a) );
+  *ao  = fmod(a + psi[i] + FOURPI, TWOPI);
+
+  /* Later */
+  return 0;
+}
+
+
 template <class T> void addToSkyMap(PixelMap<T>& Map, SigCalcTool& Tool) 
 {