Changeset 31157
- Timestamp:
- Apr 4, 2011, 1:14:37 PM (15 years ago)
- Location:
- trunk/ppSim
- Files:
-
- 5 edited
-
. (modified) (1 prop)
-
src/Makefile.am (modified) (1 diff)
-
src/ppSimDetectionLimits.c (modified) (1 diff)
-
src/ppSimDetections.c (modified) (3 diffs)
-
src/ppSimMakeGalaxies.c (modified) (2 diffs)
Legend:
- Unmodified
- Added
- Removed
-
trunk/ppSim
- Property svn:ignore
-
old new 14 14 missing 15 15 test 16 a.out.dSYM
-
- Property svn:ignore
-
trunk/ppSim/src/Makefile.am
r30618 r31157 48 48 ppSimLoop.c \ 49 49 ppSimMergeReadouts.c \ 50 ppSimDetections.c \51 50 ppSimMergeSources.c \ 52 51 ppSimMosaicChip.c \ -
trunk/ppSim/src/ppSimDetectionLimits.c
r18011 r31157 3 3 // compare injected sources (PPSIM.SOURCES) and measured fake sources (PPSIM.FAKE.SOURCES) 4 4 // to determine detection limits and recovered magnitude errors 5 // XXX this function is not used 5 6 bool ppSimDetectionLimits (pmConfig *config, pmFPAview *view) { 6 7 -
trunk/ppSim/src/ppSimDetections.c
r18011 r31157 1 1 # include "ppSim.h" 2 2 3 // XXX this function is not used 3 4 bool ppSimDetections (psImage *significance, psMetadata *recipe, psArray *sources) { 4 5 psAssert (sources, "programming error: ppSimDetections passed NULL sources"); … … 12 13 float SIGMA_SMTH = psMetadataLookupF32 (&status, recipe, "SIGMA_SMOOTH"); 13 14 psAssert (status, "SIGMA_SMOOTH missing: call psphotSignificanceImage first"); 14 float effArea = 4.0*M_PI*PS_SQR(SIGMA_SMTH);15 // float effArea = 4.0*M_PI*PS_SQR(SIGMA_SMTH); 15 16 16 17 int row0 = significance->row0; … … 23 24 psAssert (peak, "peak is not defined for the source"); 24 25 25 peak->value = significance->data.F32[peak->y-row0][peak->x-col0]; 26 peak->SN = sqrt(peak->value*effArea); 27 26 peak->detValue = significance->data.F32[peak->y-row0][peak->x-col0]; 28 27 } 29 28 return true; -
trunk/ppSim/src/ppSimMakeGalaxies.c
r30618 r31157 16 16 17 17 bool galaxyGrid = psMetadataLookupBool(&mdok, recipe, "GALAXY.GRID"); // Density of fakes 18 int galaxyGridDX = psMetadataLookupS32 (&mdok, recipe, "GALAXY.GRID.DX"); // Density of fakes19 int galaxyGridDY = psMetadataLookupS32 (&mdok, recipe, "GALAXY.GRID.DY"); // Density of fakes18 int galaxyGridDX = psMetadataLookupS32 (&mdok, recipe, "GALAXY.GRID.DX"); // Density of fakes 19 int galaxyGridDY = psMetadataLookupS32 (&mdok, recipe, "GALAXY.GRID.DY"); // Density of fakes 20 20 bool galaxyGridRandom = psMetadataLookupBool(&mdok, recipe, "GALAXY.GRID.RANDOM"); // Density of fakes 21 float galaxyGridPeak = psMetadataLookupF32 (&mdok, recipe, "GALAXY.GRID.PEAK"); // peak flux of fakes 21 22 22 23 float galaxyRmajorMax = psMetadataLookupF32(&mdok, recipe, "GALAXY.RMAJOR.MAX"); // Density of fakes … … 107 108 galaxy->x = ix; 108 109 galaxy->y = iy; 109 110 // galaxy->peak = 1000; 111 galaxy->peak = bright / 100.0; 112 // fprintf (stderr, "Io: %f\n", bright); 110 galaxy->peak = galaxyGridPeak; 113 111 114 112 // galaxyIndex from user should be for function of this form: exp(-r^(1/n))
Note:
See TracChangeset
for help on using the changeset viewer.
