Changeset 21365 for trunk/ppSim/src/ppSimPhotomReadoutForce.c
- Timestamp:
- Feb 5, 2009, 4:47:00 PM (17 years ago)
- File:
-
- 1 edited
-
trunk/ppSim/src/ppSimPhotomReadoutForce.c (modified) (7 diffs)
Legend:
- Unmodified
- Added
- Removed
-
trunk/ppSim/src/ppSimPhotomReadoutForce.c
r21183 r21365 12 12 } 13 13 14 # if 0 14 # if 0 15 15 // set the photcode for this image 16 16 if (!psphotAddPhotcode (recipe, config, view, "PPSIM.FAKE.CHIP")) { … … 22 22 // *** in this section, perform the photometry for real + force sources on PPSIM.FORCE.CHIP *** 23 23 24 // find the currently selected readout. 24 // find the currently selected readout. 25 25 // we always perform photometry on the mosaiced chip 26 26 pmReadout *readout = pmFPAfileThisReadout (config->files, view, "PPSIM.FORCE.CHIP"); … … 35 35 psAssert (forceSources, "failed to load force photometry sources"); 36 36 37 // Generate the mask and weightimages, including the user-defined analysis region of interest38 psphotSetMaskAnd Weight(config, readout, recipe);37 // Generate the mask and variance images, including the user-defined analysis region of interest 38 psphotSetMaskAndVariance (config, readout, recipe); 39 39 40 40 // XXX need to define the source pixels … … 46 46 psArray *realSources = psArrayAlloc (realMeasuredSources->n); 47 47 for (int i = 0; i < realMeasuredSources->n; i++) { 48 realSources->data[i] = pmSourceCopy (realMeasuredSources->data[i]);48 realSources->data[i] = pmSourceCopy (realMeasuredSources->data[i]); 49 49 } 50 50 51 51 // load the psf model, if suppled. FWHM_X,FWHM_Y,etc are saved in the recipe 52 // this function uses PSPHOT.PSF.LOAD as the pmFPAfile 52 // this function uses PSPHOT.PSF.LOAD as the pmFPAfile 53 53 pmPSF *psf = psphotLoadPSF (config, view, recipe); 54 54 assert (psf); 55 55 56 // remove all sources 56 // remove all sources 57 57 psphotRemoveAllSources (realSources, recipe); 58 58 … … 69 69 70 70 // XXX fake sources should measure peak->x,y, force sources should not 71 psImageMaskType maskVal = 0xff; 71 psImageMaskType maskVal = 0xff; 72 72 psImage *significance = psphotSignificanceImage (readout, recipe, 1, maskVal); 73 73 ppSimDetections (significance, recipe, forceSources); … … 83 83 psphotGuessModels (config, readout, realSources, psf); 84 84 psphotGuessModels (config, readout, forceSources, psf); 85 85 86 86 // linear fit to real + force sources 87 87 psArray *sources = ppSimMergeSources (realSources, forceSources); … … 106 106 pmReadout *forceReadout = pmFPAfileThisReadout (config->files, view, "PPSIM.FORCE.SOURCES"); 107 107 if (!forceReadout) { 108 forceReadout = pmReadoutAlloc (forceCell);109 psFree (forceReadout); // there is a copy on 'cell' as well108 forceReadout = pmReadoutAlloc (forceCell); 109 psFree (forceReadout); // there is a copy on 'cell' as well 110 110 } 111 111 psAssert (forceReadout, "no forceReadout?");
Note:
See TracChangeset
for help on using the changeset viewer.
