Changeset 26788 for branches/eam_branches/20091201/psphot/src/psphotOutput.c

Timestamp:

Feb 5, 2010, 1:38:43 PM (16 years ago)

Author:

eugene

Message:

updates to psphot APIs to enable stack photometry

Location:

branches/eam_branches/20091201/psphot

Files:

• . (modified) (1 prop)
• src/psphotOutput.c (modified) (4 diffs)

branches/eam_branches/20091201/psphot/src/psphotOutput.c

@@ -126 +126 @@
 }
 
-bool psphotAddPhotcodeReadout (pmConfig *config, const pmFPAview *view, const char *filerule, int index) {
-
-    pmFPAfile *file = pmFPAfileSelectSingle(config->files, filerule, index); // File of interest
-    PS_ASSERT (file, false);
-
-    pmReadout  *readout = pmFPAviewThisReadout (view, file->fpa);
-
-    // determine PHOTCODE from fpa & view, overwrite in readout->analysis
-    char *photcode = pmConceptsPhotcodeForView (file, view);
-    PS_ASSERT (photcode, false);
-
-    psMetadataAddStr (readout->analysis, PS_LIST_TAIL, "PHOTCODE", PS_META_REPLACE, "photcode from FPA concepts", photcode);
-    psLogMsg ("psphot", 3, "PHOTCODE is %s", photcode);
-
-    psFree (photcode);
-    return true;
-}
-
 bool psphotAddPhotcode (pmConfig *config, const pmFPAview *view) {
 
@@ -161 +143 @@
 }
 
-bool psphotSetHeaderNstars (psMetadata *recipe, psArray *sources) {
+bool psphotAddPhotcodeReadout (pmConfig *config, const pmFPAview *view, const char *filerule, int index) {
+
+    pmFPAfile *file = pmFPAfileSelectSingle(config->files, filerule, index); // File of interest
+    PS_ASSERT (file, false);
+
+    pmReadout  *readout = pmFPAviewThisReadout (view, file->fpa);
+
+    // determine PHOTCODE from fpa & view, overwrite in readout->analysis
+    char *photcode = pmConceptsPhotcodeForView (file, view);
+    PS_ASSERT (photcode, false);
+
+    psMetadataAddStr (readout->analysis, PS_LIST_TAIL, "PHOTCODE", PS_META_REPLACE, "photcode from FPA concepts", photcode);
+    psLogMsg ("psphot", 3, "PHOTCODE is %s", photcode);
+
+    psFree (photcode);
+    return true;
+}
+
+bool psphotSetHeaderNstars (psMetadata *header, psArray *sources) {
 
     int nSrc = 0;
@@ -190 +190 @@
     }
 
-    psMetadataAddS32 (recipe, PS_LIST_TAIL, "NSTARS",   PS_META_REPLACE, "Number of sources", nSrc);
-    psMetadataAddS32 (recipe, PS_LIST_TAIL, "NDET_EXT", PS_META_REPLACE, "Number of extended sources", nEXT);
-    psMetadataAddS32 (recipe, PS_LIST_TAIL, "NDET_CR",  PS_META_REPLACE, "Number of cosmic rays", nCR);
-    return true;
-}
-
-// these values are saved in an output header stub - they are added to either the
-// PHU header (CMP) or the MEF table header (CMF)
-// XXX these are currently carried by the recipe -- this will not work in a Stack context
-// XXX they should be place in the readout->analysis of the given image
-psMetadata *psphotDefineHeader (psMetadata *recipe) {
+    psMetadataAddS32 (header, PS_LIST_TAIL, "NSTARS",   PS_META_REPLACE, "Number of sources", nSrc);
+    psMetadataAddS32 (header, PS_LIST_TAIL, "NDET_EXT", PS_META_REPLACE, "Number of extended sources", nEXT);
+    psMetadataAddS32 (header, PS_LIST_TAIL, "NDET_CR",  PS_META_REPLACE, "Number of cosmic rays", nCR);
+    return true;
+}
+
+// these values are saved in an output header stub - they are added to either the PHU header
+// (CMP) or the MEF table header (CMF).  before the header is created, each readout has these
+// values stored on readout->analysis
+psMetadata *psphotDefineHeader (psMetadata *analysis) {
 
     bool status = true;
@@ -207 +206 @@
 
     // write necessary information to output header
-    psMetadataItemSupplement (&status, header, recipe, "ZERO_PT");
-    psMetadataItemSupplement (&status, header, recipe, "PHOTCODE");
-
-    psMetadataItemSupplement (&status, header, recipe, "APMIFIT");
-    psMetadataItemSupplement (&status, header, recipe, "DAPMIFIT");
-    psMetadataItemSupplement (&status, header, recipe, "NAPMIFIT");
+    psMetadataItemSupplement (&status, header, analysis, "ZERO_PT");
+    psMetadataItemSupplement (&status, header, analysis, "PHOTCODE");
+
+    psMetadataItemSupplement (&status, header, analysis, "APMIFIT");
+    psMetadataItemSupplement (&status, header, analysis, "DAPMIFIT");
+    psMetadataItemSupplement (&status, header, analysis, "NAPMIFIT");
 
     // PSF model parameters (shape values for image center)
-    psMetadataItemSupplement (&status, header, recipe, "NPSFSTAR");
-    psMetadataItemSupplement (&status, header, recipe, "APLOSS");
-
-    psMetadataItemSupplement (&status, header, recipe, "FWHM_MAJ");
-    psMetadataItemSupplement (&status, header, recipe, "FW_MJ_SG");
-    psMetadataItemSupplement (&status, header, recipe, "FW_MJ_LQ");
-    psMetadataItemSupplement (&status, header, recipe, "FW_MJ_UQ");
-
-    psMetadataItemSupplement (&status, header, recipe, "FWHM_MIN");
-    psMetadataItemSupplement (&status, header, recipe, "FW_MN_SG");
-    psMetadataItemSupplement (&status, header, recipe, "FW_MN_LQ");
-    psMetadataItemSupplement (&status, header, recipe, "FW_MN_UQ");
-
-    psMetadataItemSupplement (&status, header, recipe, "ANGLE");
+    psMetadataItemSupplement (&status, header, analysis, "NPSFSTAR");
+    psMetadataItemSupplement (&status, header, analysis, "APLOSS");
+
+    psMetadataItemSupplement (&status, header, analysis, "FWHM_MAJ");
+    psMetadataItemSupplement (&status, header, analysis, "FW_MJ_SG");
+    psMetadataItemSupplement (&status, header, analysis, "FW_MJ_LQ");
+    psMetadataItemSupplement (&status, header, analysis, "FW_MJ_UQ");
+
+    psMetadataItemSupplement (&status, header, analysis, "FWHM_MIN");
+    psMetadataItemSupplement (&status, header, analysis, "FW_MN_SG");
+    psMetadataItemSupplement (&status, header, analysis, "FW_MN_LQ");
+    psMetadataItemSupplement (&status, header, analysis, "FW_MN_UQ");
+
+    psMetadataItemSupplement (&status, header, analysis, "ANGLE");
 
     // Image Quality measurements
-    psMetadataItemSupplement (&status, header, recipe, "IQ_NSTAR");
-
-    psMetadataItemSupplement (&status, header, recipe, "IQ_FW1");
-    psMetadataItemSupplement (&status, header, recipe, "IQ_FW1_E");
-    psMetadataItemSupplement (&status, header, recipe, "IQ_FW2");
-    psMetadataItemSupplement (&status, header, recipe, "IQ_FW2_E");
-
-    psMetadataItemSupplement (&status, header, recipe, "IQ_M2");
-    psMetadataItemSupplement (&status, header, recipe, "IQ_M2_ER");
-    psMetadataItemSupplement (&status, header, recipe, "IQ_M2_LQ");
-    psMetadataItemSupplement (&status, header, recipe, "IQ_M2_UQ");
-
-    psMetadataItemSupplement (&status, header, recipe, "IQ_M2C");
-    psMetadataItemSupplement (&status, header, recipe, "IQ_M2C_E");
-    psMetadataItemSupplement (&status, header, recipe, "IQ_M2C_L");
-    psMetadataItemSupplement (&status, header, recipe, "IQ_M2C_U");
-
-    psMetadataItemSupplement (&status, header, recipe, "IQ_M2S");
-    psMetadataItemSupplement (&status, header, recipe, "IQ_M2S_E");
-    psMetadataItemSupplement (&status, header, recipe, "IQ_M2S_L");
-    psMetadataItemSupplement (&status, header, recipe, "IQ_M2S_U");
-
-    psMetadataItemSupplement (&status, header, recipe, "IQ_M3");
-    psMetadataItemSupplement (&status, header, recipe, "IQ_M3_ER");
-    psMetadataItemSupplement (&status, header, recipe, "IQ_M3_LQ");
-    psMetadataItemSupplement (&status, header, recipe, "IQ_M3_UQ");
-
-    psMetadataItemSupplement (&status, header, recipe, "IQ_M4");
-    psMetadataItemSupplement (&status, header, recipe, "IQ_M4_ER");
-    psMetadataItemSupplement (&status, header, recipe, "IQ_M4_LQ");
-    psMetadataItemSupplement (&status, header, recipe, "IQ_M4_UQ");
+    psMetadataItemSupplement (&status, header, analysis, "IQ_NSTAR");
+
+    psMetadataItemSupplement (&status, header, analysis, "IQ_FW1");
+    psMetadataItemSupplement (&status, header, analysis, "IQ_FW1_E");
+    psMetadataItemSupplement (&status, header, analysis, "IQ_FW2");
+    psMetadataItemSupplement (&status, header, analysis, "IQ_FW2_E");
+
+    psMetadataItemSupplement (&status, header, analysis, "IQ_M2");
+    psMetadataItemSupplement (&status, header, analysis, "IQ_M2_ER");
+    psMetadataItemSupplement (&status, header, analysis, "IQ_M2_LQ");
+    psMetadataItemSupplement (&status, header, analysis, "IQ_M2_UQ");
+
+    psMetadataItemSupplement (&status, header, analysis, "IQ_M2C");
+    psMetadataItemSupplement (&status, header, analysis, "IQ_M2C_E");
+    psMetadataItemSupplement (&status, header, analysis, "IQ_M2C_L");
+    psMetadataItemSupplement (&status, header, analysis, "IQ_M2C_U");
+
+    psMetadataItemSupplement (&status, header, analysis, "IQ_M2S");
+    psMetadataItemSupplement (&status, header, analysis, "IQ_M2S_E");
+    psMetadataItemSupplement (&status, header, analysis, "IQ_M2S_L");
+    psMetadataItemSupplement (&status, header, analysis, "IQ_M2S_U");
+
+    psMetadataItemSupplement (&status, header, analysis, "IQ_M3");
+    psMetadataItemSupplement (&status, header, analysis, "IQ_M3_ER");
+    psMetadataItemSupplement (&status, header, analysis, "IQ_M3_LQ");
+    psMetadataItemSupplement (&status, header, analysis, "IQ_M3_UQ");
+
+    psMetadataItemSupplement (&status, header, analysis, "IQ_M4");
+    psMetadataItemSupplement (&status, header, analysis, "IQ_M4_ER");
+    psMetadataItemSupplement (&status, header, analysis, "IQ_M4_LQ");
+    psMetadataItemSupplement (&status, header, analysis, "IQ_M4_UQ");
 
     // XXX these need to be defined from elsewhere
     psMetadataAdd (header, PS_LIST_TAIL, "FSATUR",   PS_DATA_F32 | PS_META_REPLACE, "SATURATION MAG",      0.0);
     psMetadataAdd (header, PS_LIST_TAIL, "FLIMIT",   PS_DATA_F32 | PS_META_REPLACE, "COMPLETENESS MAG",    0.0);
-    psMetadataItemSupplement (&status, header, recipe, "NSTARS");
-
-    psMetadataItemSupplement (&status, header, recipe, "NDET_EXT");
-    psMetadataItemSupplement (&status, header, recipe, "NDET_CR");
+    psMetadataItemSupplement (&status, header, analysis, "NSTARS");
+
+    psMetadataItemSupplement (&status, header, analysis, "NDET_EXT");
+    psMetadataItemSupplement (&status, header, analysis, "NDET_CR");
 
     // sky background model statistics
-    psMetadataItemSupplement (&status, header, recipe, "MSKY_MN");
-    psMetadataItemSupplement (&status, header, recipe, "MSKY_SIG");
-    psMetadataItemSupplement (&status, header, recipe, "MSKY_MIN");
-    psMetadataItemSupplement (&status, header, recipe, "MSKY_MAX");
-    psMetadataItemSupplement (&status, header, recipe, "MSKY_NX");
-    psMetadataItemSupplement (&status, header, recipe, "MSKY_NY");
+    psMetadataItemSupplement (&status, header, analysis, "MSKY_MN");
+    psMetadataItemSupplement (&status, header, analysis, "MSKY_SIG");
+    psMetadataItemSupplement (&status, header, analysis, "MSKY_MIN");
+    psMetadataItemSupplement (&status, header, analysis, "MSKY_MAX");
+    psMetadataItemSupplement (&status, header, analysis, "MSKY_NX");
+    psMetadataItemSupplement (&status, header, analysis, "MSKY_NY");
 
     psMetadataAddF32 (header, PS_LIST_TAIL, "DT_PHOT", PS_META_REPLACE, "elapsed psphot time", psTimerMark ("psphotReadout"));