Changeset 18611

Timestamp:

Jul 17, 2008, 5:19:07 PM (18 years ago)

Author:

Paul Price

Message:

Need to use Poisson errors in order to get realistic weight values.

File:

• trunk/psModules/src/imcombine/pmPSFEnvelope.c (modified) (3 diffs)

trunk/psModules/src/imcombine/pmPSFEnvelope.c

@@ -34 +34 @@
 #define PEAK_FLUX 1.0e4                 // Peak flux for each source
 #define SKY_VALUE 0.0e0                 // Sky value for fake image
-#define WEIGHT_VAL 10.0                 // Weighting for image
+#define WEIGHT_VAL 3.0                  // Weighting for image
 #define WEIGHT_FACTOR 10.0              // Factor to multiply image by to get weighting
 #define PSF_STATS PS_STAT_SAMPLE_MEDIAN | PS_STAT_SAMPLE_STDEV // Statistics options for measuring PSF
@@ -251 +251 @@
     // Don't assume Poisson errors
     pmPSFOptions *options = pmPSFOptionsAlloc(); // Options for fitting a PSF
-    options->poissonErrorsPhotLMM = false;
+    options->poissonErrorsPhotLMM = true;
     options->poissonErrorsPhotLin = false;
-    options->poissonErrorsParams = false;
+    options->poissonErrorsParams = true;
     options->stats = psStatsAlloc(PSF_STATS);
     options->radius = radius;
@@ -264 +264 @@
     options->psfFieldYo = 0;
 
-    pmSourceFitModelInit(50, 0.01, WEIGHT_VAL, false);
+    pmSourceFitModelInit(50, 0.01, WEIGHT_VAL, true);
 
     pmPSFtry *try = pmPSFtryModel(fakes, modelName, options, 0, 0xff);