Changeset 24391 for trunk/extsrc/gpcsw/gpcsrc/analysis/libpscoords/pscoords.h

Timestamp:

Jun 12, 2009, 3:53:07 PM (17 years ago)

Author:

eugene

Message:

updates to burntool-and-related from camera svn

File:

• trunk/extsrc/gpcsw/gpcsrc/analysis/libpscoords/pscoords.h (modified) (6 diffs)

trunk/extsrc/gpcsw/gpcsrc/analysis/libpscoords/pscoords.h

@@ -2 +2 @@
 #define NINT(x) (x<0?(int)((x)-0.5):(int)((x)+0.5))
 
+#define OTA_CTE (5.0e-6)        /* Coefficient of thermal expansion of Si+Mo */
+#define CFFP_CTE (0.5e-6)       /* Coefficient of thermal expansion of CFFP */
+#define DELTA_T (-100.0)        /* Operating temp minus metrology temp */
+
+#define OTA_SCALE (1.0+OTA_CTE*DELTA_T) /* Room temp um to operating um of Si+Mo */
+#define CFFP_SCALE ((1.0+CFFP_CTE*DELTA_T)*25400) /* Room temp inches to um on CFFP */
+
 /* Silicon layout of OTAs (um) */
-#define PSC_PIXEL     10.0      /* Pixel size of an OTA */
-#define PSC_HCELL   5900.0      /* Horizontal cell size */
-#define PSC_VCELL   5980.0      /* Vertical cell size */
-#define PSC_HSTREET  120.0      /* Horizontal street between cells (HS) */
-#define PSC_VSTREET  180.0      /* Vertical street between cells (VS) */
-#define PSC_LBORDER  442.0      /* Left border (to left-most cell) (LB) */
-#define PSC_RBORDER  418.0      /* Right border (to right-most vertical street) (RB) */
-#define PSC_TBORDER  354.5      /* Top border (to top-most cell) (TB) */
-#define PSC_BBORDER  975.5      /* Bottom border (to lowest horizontal street) (BB) */
-#define PSC_HDIE   49500.0      /* Horizontal size of silicon die */
-#define PSC_VDIE   50130.0      /* Vertical size of silicon die */
+#define PSC_PIXEL     (10.0*OTA_SCALE)  /* Pixel size of an OTA */
+#define PSC_HCELL   (5900.0*OTA_SCALE)  /* Horizontal cell size */
+#define PSC_VCELL   (5980.0*OTA_SCALE)  /* Vertical cell size */
+#define PSC_HSTREET  (120.0*OTA_SCALE)  /* Horizontal street between cells (HS) */
+#define PSC_VSTREET  (180.0*OTA_SCALE)  /* Vertical street between cells (VS) */
+#define PSC_LBORDER  (442.0*OTA_SCALE)  /* Left border (to left-most cell) (LB) */
+#define PSC_RBORDER  (418.0*OTA_SCALE)  /* Right border (to right-most vertical street) (RB) */
+#define PSC_TBORDER  (354.5*OTA_SCALE)  /* Top border (to top-most cell) (TB) */
+#define PSC_BBORDER  (975.5*OTA_SCALE)  /* Bottom border (to lowest horizontal street) (BB) */
+#define PSC_HDIE   (49500.0*OTA_SCALE)  /* Horizontal size of silicon die */
+#define PSC_VDIE   (50130.0*OTA_SCALE)  /* Vertical size of silicon die */
 
 /* Mechanical layout of OTAs (um) */
-#define PSC_HMECH  (1.957*25400)        /* Horizontal mech spacing of OTA placement */
-#define PSC_VMECH  (2.025*25400)        /* Vertical mech spacing of OTA placement */
-#define PSC_VMOFF (PSC_VDIE-(1.568+0.375)*25400) /* Vert die offset between sides */
+#define PSC_HMECH  (1.957*CFFP_SCALE)   /* Horizontal mech spacing of OTA placement */
+#define PSC_VMECH  (2.025*CFFP_SCALE)   /* Vertical mech spacing of OTA placement */
+#define PSC_VMOFF (PSC_VDIE-(1.568+0.375)*CFFP_SCALE) /* Vert die offset between sides */
 
 #define PSC_REFRACT_CONST 60.0  /* Standard refraction ("/tanz) at STP */
@@ -123 +130 @@
 
 /* Overall focal plane offset and rotation (um, rad) */
-const PSC_OFFROT_T psc_fpoff;
+PSC_OFFROT_T psc_fpoff;
 
 /* Offsets and rotations of each OTA (um, millirad) */
-const PSC_OFFROT_T psc_otaoff[PSC_NX*PSC_NY];
-
+PSC_OFFROT_T psc_otaoff[PSC_NX*PSC_NY];
+
+
+/* This is from Run 3, 090516 */
+#define PS_scale 38.7932        /* Default PS plate scale [um/arcsec] */
+#define PS_d2  9.78e-7  /* Default quadratic PS distortion [arcsec^-1] */
+#define PS_d3  3.16e-11 /* Default cubic PS distortion [arcsec^-2] */
+
+
+#ifdef RUN_THREE_V1     /* This is from Run 3, 090228 */
 #define PS_scale 38.860 /* Default PS plate scale [um/arcsec] */
-#define PS_d3  1.49e-10 /* Default PS distortion [arcsec^-2] */
+#define PS_d3  1.49e-10 /* Default cubic PS distortion [arcsec^-2] */
+#endif
+
+
 #define PS_airdens 0.71 /* Default PS1 qir density (Haleakala) */
 
@@ -139 +157 @@
  * ---------------------------------------------------------------------------
  *   psc_tproject(a, d, a0, d0, pa, density, pi/2-alt, vertical, &x, &y);
- *   psc_psoptics(x/sec2rad, y/sec2rad, dx, dy, dpa, pscale, d3, &xfp, &yfp);
+ *   psc_psoptics2(x/sec2rad, y/sec2rad, dx, dy, dpa, pscale, d2, d3, &xfp, &yfp);
  *   psc_fp_to_pixel(xfp, yfp, &ota_xid, &ota_yid, &xota, &yota);
  *   psc_pixel_to_cell(xota, yota, &cell_xid, &cell_yid, &xcell, &ycell);
+ *
+ *  can also use
+ *   psc_psoptics(x/sec2rad, y/sec2rad, dx, dy, dpa, pscale, d3, &xfp, &yfp);
+ *
  *
  * To go from cell coords to a position (a,d) relative to (a0,d0,pa) pointing:
@@ -147 +169 @@
  *   psc_cell_to_pixel(cell_xid,cell_yid, xcell,ycell, &xota,&yota);
  *   psc_pixel_to_fp(ota_xid, ota_yid, xota, yota, &xfp, &yfp);
+ *   psc_invoptics2(xfp, yfp, dx, dy, dpa, pscale, d2, d3, &x, &y);
+ *   psc_tplonglat(x*sec2rad, y*sec2rad, a0, d0, pa, density,  pi/2-alt, vertical, &a, &d);
+ *
+ *  can also use
  *   psc_invoptics(xfp, yfp, dx, dy, dpa, pscale, d3, &x, &y);
- *   psc_tplonglat(x*sec2rad, y*sec2rad, a0, d0, pa, density,  pi/2-alt, vertical, &a, &d);
  */
 
@@ -155 +180 @@
 /**************/
 /* Enable application of chip offsets? */
-int psc_do_chipoff(int doit);   /* Apply chip offsets? (0/1, default=0) */
+int psc_do_chipoff(int doit);   /* Apply chip offsets? (0/1, default=1) */
 
 /* Convert a focal plane position to OTA otax,otay, Cell cellx,celly, Pixel */
@@ -202 +227 @@
                  double pscale/*38.86um/sec*/, double d3/*1.49e-10sec^-2*/, 
                  double *xfp_um, double *yfp_um);
+
+/* Evaluate an inverse optical model: microns -> arcsec notionally */
+/* Offset, scale with *small* distortion, rotate */
+int psc_invoptics2(double xfp_um, double yfp_um, double dx_sec, double dy_sec, double dpa_rad, 
+                   double pscale/*38.793um/sec*/,
+                   double d2/*9.8e-7sec^-2*/, double d3/*3.1e-11sec^-2*/, 
+                   double *x_sec, double *y_sec);
+
+/* Evaluate an optical model: arcsec -> microns notionally */
+/* Offset, scale with distortion, rotate */
+int psc_psoptics2(double x_sec, double y_sec, double dx_sec, double dy_sec, double dpa_rad, 
+                  double pscale/*38.793um/sec*/,
+                  double d2/*9.8e-7sec^-2*/, double d3/*3.1e-11sec^-2*/, 
+                  double *xfp_um, double *yfp_um);
+
+/* Enable application of chip offsets (um, mrad)? */
+int psc_load_otaoff(const char *fname);