Changeset 37321 for trunk/psphot/src

Timestamp:

Aug 27, 2014, 12:20:33 PM (12 years ago)

Author:

bills

Message:

In psphotFullForce, add code to read the warp to chip transformations stored in the
warp headers and identify the chip that generated the pixels for each source and
the chip coordinates of the peak.
Save this information in PS1_V5 cmf if the lensing parameters are also being saved.
(full force only)

Location:

trunk/psphot/src

Files:

• Makefile.am (modified) (2 diffs)
• psphot.h (modified) (1 diff)
• psphotChipParams.c (added)
• psphotFullForceReadout.c (modified) (1 diff)
• psphotSetThreads.c (modified) (1 diff)

trunk/psphot/src/Makefile.am

@@ -97 +97 @@
 # generalized forced photometry (including Kron, Petro, and Models) of specified positions given a specified psf
 psphotFullForce_SOURCES = \
-        psphotFullForce.c             \
-        psphotFullForceArguments.c    \
+        psphotFullForce.c          \
+        psphotFullForceArguments.c \
         psphotParseCamera.c        \
         psphotImageLoop.c          \
@@ -259 +259 @@
         psphotEfficiency.c             \
         psphotSetNFrames.c             \
-        psphotSourceMemory.c
+        psphotSourceMemory.c           \
+        psphotChipParams.c
 
 # not currently used

trunk/psphot/src/psphot.h

@@ -571 +571 @@
 bool psphotSourceMemoryReadout(pmConfig *config, const pmFPAview *view, const char *filerule, int index);
 
+bool psphotChipParams(pmConfig *config, const pmFPAview *view, const char *filerule);
+bool psphotChipParams_Threaded(psThreadJob *job);
+
 const char * psphotGetFilerule(const char *baseRule);
 extern bool psphotINpsphotStack;
 
+
 #endif

trunk/psphot/src/psphotFullForceReadout.c

@@ -106 +106 @@
     psphotKronFlux (config, view, filerule);
 
+    psphotChipParams (config, view, filerule);
+
     // measure petro fluxes
     psphotPetroFlux (config, view, filerule);

trunk/psphot/src/psphotSetThreads.c

@@ -87 +87 @@
     psFree(task);
 
+    task = psThreadTaskAlloc("PSPHOT_CHIP_PARAMS", 2);
+    task->function = &psphotChipParams_Threaded;
+    psThreadTaskAdd(task);
+    psFree(task);
+
     return true;
 }