Changeset 37569

Timestamp:

Nov 8, 2014, 7:34:12 PM (12 years ago)

Author:

eugene

Message:

save fake qsos in ICRF format file

Location:

branches/eam_branches/ipp-20140904/Ohana/src/fakeastro/src

Files:

• ConfigInit.c (modified) (1 diff)
• make_2mass_measures.c (modified) (1 diff)
• make_fake_stars_catalog.c (modified) (1 diff)
• make_fakeqsos.c (modified) (2 diffs)

branches/eam_branches/ipp-20140904/Ohana/src/fakeastro/src/ConfigInit.c

@@ -55 +55 @@
   }
   if (!ScanConfig (config, "FAKEASTRO_REF_EPOCH", "%s", 0, FAKEASTRO_REF_EPOCH)) {
-    strcpy (FAKEASTRO_REF_EPOCH, "2011/05/11,00:00:00"); // epoch of reference catalog
+    strcpy (FAKEASTRO_REF_EPOCH, "2000/01/01,00:00:00"); // epoch of truth positions
   }
   if (!ScanConfig (config, "FAKEASTRO_2MASS_EPOCH", "%s", 0, FAKEASTRO_2MASS_EPOCH)) {
-    strcpy (FAKEASTRO_2MASS_EPOCH, "2000/01/01,00:00:00"); // epoch of reference catalog
+    strcpy (FAKEASTRO_2MASS_EPOCH, "2000/01/01,00:00:00"); // epoch of 2MASS astrometry
   }
 

branches/eam_branches/ipp-20140904/Ohana/src/fakeastro/src/make_2mass_measures.c

@@ -11 +11 @@
   Measure *measure = catalog->measure;
   
+  // we have starpar.R,D, which represent the true positions at the reference epoch, FAKEASTRO_REF_EPOCH
+  // the 2MASS stars are generated at locations for the 2MASS EPOCH
   time_t timeRef = ohana_date_to_sec(FAKEASTRO_REF_EPOCH);
   time_t tzero_2mass = ohana_date_to_sec(FAKEASTRO_2MASS_EPOCH);

branches/eam_branches/ipp-20140904/Ohana/src/fakeastro/src/make_fake_stars_catalog.c

@@ -20 +20 @@
   int Nstars = *nstars;
   
+  // we have starpar.R,D, which represent the true positions at the reference epoch, FAKEASTRO_REF_EPOCH
   time_t timeRef = ohana_date_to_sec(FAKEASTRO_REF_EPOCH);
   

branches/eam_branches/ipp-20140904/Ohana/src/fakeastro/src/make_fakeqsos.c

@@ -14 +14 @@
 
   CoordTransform *transback = InitTransform (COORD_CELESTIAL, COORD_GALACTIC);
+
+  FILE *f = NULL;
+  if (ICRF) {
+    f = fopen ("fake.icrf.dat", "w");
+    if (!f) {
+      fprintf (stderr, "unable to open fake.icrf.dat\n");
+    }
+  } else {
+    f = fopen ("fake.zero.dat", "w");
+    if (!f) {
+      fprintf (stderr, "unable to open fake.zero.dat\n");
+    }
+  } 
 
   int i;
@@ -68 +81 @@
     stars[i].starpar.uRA  = 0.0;
     stars[i].starpar.uDEC = 0.0;
+
+/*
+#                       hr mn seconds   deg mn seconds      mas    mas                   Jy     Jy      Jy     Jy      Jy     Jy     Jy     Jy      Jy     Jy
+C  2357+080 J0000+0816  00 00 07.031141 +08 16 45.05175    0.46   0.85   0.758     41  -1.00  -1.00   -1.00  -1.00    0.020 <0.014  -1.00  -1.00   -1.00  -1.00   X    rfc_2014c
+01234567890123456789012345678901234567890123456789012345678901234567890123456789012345678901234567890123456789012345678901234567890123456789012345678901234567890123456789
+0         1         2         3         4         5         6         7         8         9         0         1         2         3         4         5         6         
+*/
+
+    if (f) {
+      char Rline[32], Dline[32];
+      {
+        int flag = SIGN(R);
+        double hours = fabs(R / 15.0);  /* convert from degrees to hours */
+        int h = hours;
+        int m = 60.000001*(hours - h);
+        float s = 3600*(hours - h - m / 60.0);
+        if (flag > 0) {
+          snprintf (Rline, 32, " %02d %02d %09.6f  ", h, m, s);
+        } else {
+          snprintf (Rline, 32, "-%02d %02d %09.6f  ", h, m, s);
+        }       
+      }
+      {
+        int flag = SIGN(D);
+        double dec = fabs(D);
+        int d = dec;
+        int m = 60.000001*(dec - d);
+        float s = 3600*(dec - d - m / 60.0);
+        if (flag > 0) {
+          snprintf (Dline, 32, " %02d %02d %09.6f  ", d, m, s);
+        } else {
+          snprintf (Dline, 32, "-%02d %02d %09.6f  ", d, m, s);
+        }       
+      }
+
+      fprintf (f, "C  2357+080 J0000+0816 %16s %15s    0.46   0.85   0.758     41  -1.00  -1.00   -1.00  -1.00    0.020 <0.014  -1.00  -1.00   -1.00  -1.00   X    rfc_2014c\n", Rline, Dline);
+    }
   }
   fprintf (stderr, "\n");