Changeset 10096 for trunk/psphot/src/psphotWeightBias.c

Timestamp:

Nov 19, 2006, 7:26:08 PM (20 years ago)

Author:

eugene

Message:

working on global linear fit with sky level

File:

• trunk/psphot/src/psphotWeightBias.c (modified) (1 diff)

trunk/psphot/src/psphotWeightBias.c

@@ -18 +18 @@
 
     // set fitting method to use non-poisson errors (local sky error)
-    pmSourceFitModelInit (15, 0.1, false);
+    float SKY_STDEV = psMetadataLookupF32 (&status, recipe, "SKY_STDEV");
+    if (!status) {
+      SKY_STDEV = 1.0;
+        psWarning("SKY_STDEV is not set --- defaulting to %f\n", SKY_STDEV);
+    }
+    // use poissonian errors or local-sky errors
+    bool POISSON_ERRORS = psMetadataLookupBool (&status, recipe, "POISSON_ERRORS");
+    if (!status) {
+        POISSON_ERRORS = true;
+        psWarning("POISSON_ERRORS is not set in the recipe --- defaulting to true.\n");
+    }
+    pmSourceFitModelInit (15, 0.1, PS_SQR(SKY_STDEV), POISSON_ERRORS);
 
     // option to limit analysis to a specific region