Changeset 34980

Timestamp:

Jan 18, 2013, 11:43:34 AM (13 years ago)

Author:

bills

Message:

Correct the calculation of the magnitude limit for getstar
for skycalibration runs by using fixed zero point of 25
specified in new recipe value DVO.GETSTAR.FIXED.ZEROPT

Location:

tags/ipp-20121218

Files:

• ippconfig/gpc1/psastro.config (modified) (3 diffs, 1 prop)
• ippconfig/recipes/psastro.config (modified) (3 diffs, 1 prop)
• psastro/src/psastroLoadRefstars.c (modified) (2 diffs, 1 prop)

tags/ipp-20121218/ippconfig/gpc1/psastro.config

@@ -382 +382 @@
    PSASTRO.USE.MODEL           BOOL FALSE
    PSASTRO.MOSAIC.MODE         BOOL FALSE
+   DVO.GETSTAR.FIXED.ZEROPT    F32   25.0
 
    PSASTRO.CATDIR                       STR     PS1.REF.20120524
@@ -394 +395 @@
 END
 
-# initial setup for MD@100 inputs, should work for 50-500 input
+STATICSKY_CAL   METADATA
+    REFSTAR_MASK                BOOL FALSE
+    PSASTRO.SAVE.REFMATCH       BOOL FALSE
+    PSASTRO.FIX.CHIPS           BOOL FALSE
+    PSASTRO.USE.MODEL           BOOL FALSE
+    PSASTRO.MOSAIC.MODE         BOOL FALSE
+    DVO.GETSTAR.FIXED.ZEROPT    F32   25.0
+
+    PSASTRO.CATDIR                       STR     PS1.REF.20120524
+    ZERO.POINT.USE.MEAN                  BOOL    TRUE
+
+    ## for skycell, radius shouldn't need to be as large as 90 pixel default. deeper stacks, crowded fields can be problems
+    PSASTRO.MATCH.RADIUS.N0 F32    15
+    PSASTRO.MATCH.RADIUS.N1 F32    10
+    PSASTRO.MATCH.RADIUS.N2 F32    5
+    PSASTRO.MATCH.RADIUS.N3 F32    1
+    PSASTRO.MATCH.RADIUS.N4 F32    0.5
+    PSASTRO.MATCH.FIT.NITER S32    5
+
+    PSASTRO.MAX.ERROR      F32      0.2  # in arcsec units, if larger than 0.2" on skycell then a problem..
+
+END
+
+
+# initial setup for MD@10 inputs, should work for 2-30
+## MEH still testing with MD09 GR0 nightly 
+STATICSKY_NIGHTCAL   METADATA
+   REFSTAR_MASK                BOOL FALSE
+   PSASTRO.SAVE.REFMATCH       BOOL FALSE
+   PSASTRO.FIX.CHIPS           BOOL FALSE
+   PSASTRO.USE.MODEL           BOOL FALSE
+   PSASTRO.MOSAIC.MODE         BOOL FALSE
+   DVO.GETSTAR.FIXED.ZEROPT    F32   25.0
+
+   PSASTRO.CATDIR                       STR     PS1.REF.20120524
+   ZERO.POINT.USE.MEAN                  BOOL    TRUE
+
+   ## for skycell, radius shouldn't need to be as large as 90 pixel default. deeper stacks can be problems
+   PSASTRO.MATCH.RADIUS.N0 F32    15
+   PSASTRO.MATCH.RADIUS.N1 F32    10
+   PSASTRO.MATCH.RADIUS.N2 F32    5
+   PSASTRO.MATCH.RADIUS.N3 F32    1
+   PSASTRO.MATCH.RADIUS.N4 F32    0.5
+   PSASTRO.MATCH.FIT.NITER S32    5
+
+   PSASTRO.MAX.ERROR      F32      0.2  # in arcsec units
+   PSASTRO.MIN.NSTAR      S32      10   # edges may not have many and require special treatment 
+
+END
+
+# initial setup for MD@100 inputs, should work for 30-500 input
 STATICSKY_DEEPCAL   METADATA
    REFSTAR_MASK                BOOL FALSE
@@ -401 +452 @@
    PSASTRO.USE.MODEL           BOOL FALSE
    PSASTRO.MOSAIC.MODE         BOOL FALSE
+   DVO.GETSTAR.FIXED.ZEROPT    F32   25.0
 
    PSASTRO.CATDIR                       STR     PS1.REF.20120524

tags/ipp-20121218/ippconfig/recipes/psastro.config

@@ -190 +190 @@
 DVO.GETSTAR.MIN.MAG.INST      F32     -15.0
 
+# if non-zero use as the zero point for computing magnitude limits for getstar
+# by default use the ZEROPT in the PHOTCODE.DATA
+DVO.GETSTAR.FIXED.ZEROPT      F32       0.0     
+
 ZERO.POINT.BY.EXPOSURE          BOOL     TRUE
 ZERO.POINT.USE.MEAN             BOOL     TRUE
@@ -268 +272 @@
 
 STATICSKY_CAL   METADATA
-    REFSTAR_MASK                BOOL FALSE
-    PSASTRO.SAVE.REFMATCH       BOOL FALSE
-    PSASTRO.FIX.CHIPS           BOOL FALSE
-    PSASTRO.USE.MODEL           BOOL FALSE
-    PSASTRO.MOSAIC.MODE         BOOL FALSE
 END
 
@@ -278 +277 @@
 END
 
+STATICSKY_NIGHTCAL   METADATA
+END
+
 STATICSKY_DEEPCAL   METADATA
 END

tags/ipp-20121218/psastro/src/psastroLoadRefstars.c

@@ -363 +363 @@
     if (!status) ESCAPE ("missing DVO.GETSTAR.MIN.MAG.INST");
 
+    // if non zero override the zero point for the photcode (for stacks)
+    float fixedzeropt = psMetadataLookupF32(&status, recipe, "DVO.GETSTAR.FIXED.ZEROPT");
+
     // PHOTCODE.DATA is a multi of metadata items
     psListIterator *iter = psListIteratorAlloc(item->data.list, PS_LIST_HEAD, false);
@@ -391 +394 @@
             continue;
         }
+        if (fixedzeropt != 0.0) {
+            // override the recipe's zero point with the fixed value (for stacks)
+            zeropt = fixedzeropt;
+        }
 
         // convert the minInst to a calibrated minimum magnitude