Changeset 36619

Timestamp:

Mar 29, 2014, 2:19:18 PM (12 years ago)

Author:

eugene

Message:

adding lensing parameters (pmSourceLensing, psphotLensing)

Location:

branches/eam_branches/ipp-20140226

Files:

• psModules/src/camera/pmFPAfileIO.c (modified) (1 diff)
• psModules/src/camera/pmReadoutFake.c (modified) (1 diff)
• psModules/src/extras/pmThreadTools.c (modified) (1 diff)
• psModules/src/extras/pmVisual.c (modified) (1 diff)
• psModules/src/imcombine/pmPSFEnvelope.c (modified) (1 diff)
• psModules/src/imcombine/pmSubtractionSimple.c (modified) (1 diff)
• psModules/src/imcombine/pmSubtractionStamps.c (modified) (1 diff)
• psModules/src/objects (modified) (1 prop)
• psModules/src/objects/Makefile.am (modified) (5 diffs)
• psModules/src/objects/mksource.pl (modified) (1 diff)
• psModules/src/objects/models/pmModel_DEV.c (modified) (1 diff)
• psModules/src/objects/models/pmModel_EXP.c (modified) (1 diff)
• psModules/src/objects/models/pmModel_GAUSS.c (modified) (1 diff)
• psModules/src/objects/models/pmModel_PGAUSS.c (modified) (1 diff)
• psModules/src/objects/models/pmModel_PS1_V1.c (modified) (1 diff)
• psModules/src/objects/models/pmModel_QGAUSS.c (modified) (1 diff)
• psModules/src/objects/models/pmModel_RGAUSS.c (modified) (1 diff)
• psModules/src/objects/models/pmModel_SERSIC.c (modified) (1 diff)
• psModules/src/objects/models/pmModel_TRAIL.c (modified) (1 diff)
• psModules/src/objects/pmDetEff.c (modified) (1 diff)
• psModules/src/objects/pmGrowthCurve.c (modified) (1 diff)
• psModules/src/objects/pmGrowthCurveGenerate.c (modified) (1 diff)
• psModules/src/objects/pmModelUtils.c (modified) (1 diff)
• psModules/src/objects/pmMoments.h (modified) (1 diff)
• psModules/src/objects/pmPCM_MinimizeChisq.c (modified) (1 diff)
• psModules/src/objects/pmPCMdata.c (modified) (1 diff)
• psModules/src/objects/pmPSF.c (modified) (1 diff)
• psModules/src/objects/pmPSF_IO.c (modified) (1 diff)
• psModules/src/objects/pmPSFtry.c (modified) (1 diff)
• psModules/src/objects/pmPSFtryFitEXT.c (modified) (1 diff)
• psModules/src/objects/pmPSFtryFitPSF.c (modified) (1 diff)
• psModules/src/objects/pmPSFtryMakePSF.c (modified) (1 diff)
• psModules/src/objects/pmPSFtryMetric.c (modified) (1 diff)
• psModules/src/objects/pmPSFtryModel.c (modified) (1 diff)
• psModules/src/objects/pmPhotObj.c (modified) (1 diff)
• psModules/src/objects/pmSource.c (modified) (4 diffs)
• psModules/src/objects/pmSource.h (modified) (1 diff)
• psModules/src/objects/pmSourceContour.c (modified) (1 diff)
• psModules/src/objects/pmSourceExtendedPars.h (modified) (1 diff)
• psModules/src/objects/pmSourceFitModel.c (modified) (1 diff)
• psModules/src/objects/pmSourceFitPCM.c (modified) (1 diff)
• psModules/src/objects/pmSourceFitSet.c (modified) (1 diff)
• psModules/src/objects/pmSourceGroups.c (modified) (1 diff)
• psModules/src/objects/pmSourceIO.c (modified) (7 diffs)
• psModules/src/objects/pmSourceIO.h (modified) (1 diff)
• psModules/src/objects/pmSourceIO_CFF.c (modified) (1 diff)
• psModules/src/objects/pmSourceIO_CMF.c.in (modified) (2 diffs)
• psModules/src/objects/pmSourceIO_CMP.c (modified) (1 diff)
• psModules/src/objects/pmSourceIO_MatchedRefs.c (modified) (1 diff)
• psModules/src/objects/pmSourceIO_OBJ.c (modified) (1 diff)
• psModules/src/objects/pmSourceIO_PS1_CAL_0.c (modified) (1 diff)
• psModules/src/objects/pmSourceIO_PS1_DEV_0.c (modified) (1 diff)
• psModules/src/objects/pmSourceIO_PS1_DEV_1.c (modified) (1 diff)
• psModules/src/objects/pmSourceIO_RAW.c (modified) (1 diff)
• psModules/src/objects/pmSourceIO_SMPDATA.c (modified) (1 diff)
• psModules/src/objects/pmSourceIO_SX.c (modified) (1 diff)
• psModules/src/objects/pmSourceLensing.c (added)
• psModules/src/objects/pmSourceLensing.h (added)
• psModules/src/objects/pmSourceMatch.c (modified) (1 diff)
• psModules/src/objects/pmSourceMoments.c (modified) (2 diffs)
• psModules/src/objects/pmSourceOutputs.c (modified) (1 diff)
• psModules/src/objects/pmSourcePhotometry.c (modified) (1 diff)
• psModules/src/objects/pmSourcePlotApResid.c (modified) (1 diff)
• psModules/src/objects/pmSourcePlotMoments.c (modified) (1 diff)
• psModules/src/objects/pmSourcePlotPSFModel.c (modified) (1 diff)
• psModules/src/objects/pmSourceSky.c (modified) (1 diff)
• psModules/src/objects/pmSourceUtils.c (modified) (1 diff)
• psModules/src/objects/pmSourceVisual.c (modified) (1 diff)
• psModules/src/psmodules.h (modified) (1 diff)
• psphot/src/Makefile.am (modified) (1 diff)
• psphot/src/psphot.h (modified) (1 diff)
• psphot/src/psphotLensing.c (added)
• psphot/src/psphotReadout.c (modified) (2 diffs)

branches/eam_branches/ipp-20140226/psModules/src/camera/pmFPAfileIO.c

@@ -39 +39 @@
 #include "pmSourceDiffStats.h"
 #include "pmSourceSatstar.h"
+#include "pmSourceLensing.h"
 #include "pmSource.h"
 #include "pmSourceFitModel.h"

branches/eam_branches/ipp-20140226/psModules/src/camera/pmReadoutFake.c

@@ -27 +27 @@
 #include "pmSourceDiffStats.h"
 #include "pmSourceSatstar.h"
+#include "pmSourceLensing.h"
 #include "pmSource.h"
 #include "pmSourceFitModel.h"

branches/eam_branches/ipp-20140226/psModules/src/extras/pmThreadTools.c

@@ -37 +37 @@
 #include "pmSourcePhotometry.h"
 #include "pmSourceSatstar.h"
+#include "pmSourceLensing.h"
 #include "pmSource.h"
 

branches/eam_branches/ipp-20140226/psModules/src/extras/pmVisual.c

@@ -41 +41 @@
 // #include "pmSourceDiffStats.h"
 #include "pmSourceSatstar.h"
+// #include "pmSourceLensing.h"
 // #include "pmSource.h"
 // #include "pmSourceFitModel.h"

branches/eam_branches/ipp-20140226/psModules/src/imcombine/pmPSFEnvelope.c

@@ -27 +27 @@
 #include "pmSourceDiffStats.h"
 #include "pmSourceSatstar.h"
+#include "pmSourceLensing.h"
 #include "pmSource.h"
 #include "pmSourceFitModel.h"

branches/eam_branches/ipp-20140226/psModules/src/imcombine/pmSubtractionSimple.c

@@ -34 +34 @@
 #include "pmSourceSatstar.h"
 
+#include "pmSourceLensing.h"
 #include "pmSource.h"
 

branches/eam_branches/ipp-20140226/psModules/src/imcombine/pmSubtractionStamps.c

@@ -26 +26 @@
 #include "pmSourceDiffStats.h"
 #include "pmSourceSatstar.h"
+#include "pmSourceLensing.h"
 #include "pmSource.h"
 

branches/eam_branches/ipp-20140226/psModules/src/objects

@@ -12 +12 @@
 pmSourceIO_CMF_PS1_V1.v1.c
 pmSourceIO_CMF_PS1_V4.c
+pmSourceIO_CMF_PS1_V5.c
 pmSourceIO_CMF_PS1_SV1.c
 pmSourceIO_CMF_PS1_SV2.c

@@ -12 +12 @@
 pmSourceIO_CMF_PS1_V1.v1.c
 pmSourceIO_CMF_PS1_V4.c
+pmSourceIO_CMF_PS1_V5.c
 pmSourceIO_CMF_PS1_SV1.c
 pmSourceIO_CMF_PS1_SV2.c

branches/eam_branches/ipp-20140226/psModules/src/objects/Makefile.am

@@ -49 +49 @@
         pmSourceIO_CMF_PS1_V3.c \
         pmSourceIO_CMF_PS1_V4.c \
+        pmSourceIO_CMF_PS1_V5.c \
         pmSourceIO_CMF_PS1_SV1.c \
         pmSourceIO_CMF_PS1_SV2.c \
@@ -76 +77 @@
         pmGrowthCurve.c \
         pmSourceMatch.c \
+        pmSourceLensing.c \
         pmDetEff.c \
         pmSourceGroups.c \
@@ -125 +127 @@
         pmGrowthCurveGenerate.h \
         pmSourceMatch.h \
+        pmSourceLensing.h \
         pmDetEff.h \
         pmSourceGroups.h \
@@ -145 +148 @@
 pmSourceIO_CMF_PS1_V3.c \
 pmSourceIO_CMF_PS1_V4.c \
+pmSourceIO_CMF_PS1_V5.c \
 pmSourceIO_CMF_PS1_DV1.c \
 pmSourceIO_CMF_PS1_DV2.c \
@@ -163 +167 @@
         mksource.pl pmSourceIO_CMF.c.in PS1_V4 pmSourceIO_CMF_PS1_V4.c
 
+pmSourceIO_CMF_PS1_V5.c : pmSourceIO_CMF.c.in mksource.pl
+        mksource.pl pmSourceIO_CMF.c.in PS1_V5 pmSourceIO_CMF_PS1_V5.c
+
 pmSourceIO_CMF_PS1_DV1.c : pmSourceIO_CMF.c.in mksource.pl
         mksource.pl pmSourceIO_CMF.c.in PS1_DV1 pmSourceIO_CMF_PS1_DV1.c

branches/eam_branches/ipp-20140226/psModules/src/objects/mksource.pl

@@ -19 +19 @@
                "PS1_V3", 3,
                "PS1_V4", 4,
+               "PS1_V5", 5,
     );
 %cmfmodes_dv = ("PS1_DV1", 1,

branches/eam_branches/ipp-20140226/psModules/src/objects/models/pmModel_DEV.c

@@ -43 +43 @@
 #include "pmSourceDiffStats.h"
 #include "pmSourceSatstar.h"
+#include "pmSourceLensing.h"
 #include "pmSource.h"
 #include "pmSourceFitModel.h"

branches/eam_branches/ipp-20140226/psModules/src/objects/models/pmModel_EXP.c

@@ -40 +40 @@
 #include "pmSourceDiffStats.h"
 #include "pmSourceSatstar.h"
+#include "pmSourceLensing.h"
 #include "pmSource.h"
 #include "pmSourceFitModel.h"

branches/eam_branches/ipp-20140226/psModules/src/objects/models/pmModel_GAUSS.c

@@ -40 +40 @@
 #include "pmSourceDiffStats.h"
 #include "pmSourceSatstar.h"
+#include "pmSourceLensing.h"
 #include "pmSource.h"
 #include "pmSourceFitModel.h"

branches/eam_branches/ipp-20140226/psModules/src/objects/models/pmModel_PGAUSS.c

@@ -40 +40 @@
 #include "pmSourceDiffStats.h"
 #include "pmSourceSatstar.h"
+#include "pmSourceLensing.h"
 #include "pmSource.h"
 #include "pmSourceFitModel.h"

branches/eam_branches/ipp-20140226/psModules/src/objects/models/pmModel_PS1_V1.c

@@ -42 +42 @@
 #include "pmSourceDiffStats.h"
 #include "pmSourceSatstar.h"
+#include "pmSourceLensing.h"
 #include "pmSource.h"
 #include "pmSourceFitModel.h"

branches/eam_branches/ipp-20140226/psModules/src/objects/models/pmModel_QGAUSS.c

@@ -42 +42 @@
 #include "pmSourceDiffStats.h"
 #include "pmSourceSatstar.h"
+#include "pmSourceLensing.h"
 #include "pmSource.h"
 #include "pmSourceFitModel.h"

branches/eam_branches/ipp-20140226/psModules/src/objects/models/pmModel_RGAUSS.c

@@ -41 +41 @@
 #include "pmSourceDiffStats.h"
 #include "pmSourceSatstar.h"
+#include "pmSourceLensing.h"
 #include "pmSource.h"
 #include "pmSourceFitModel.h"

branches/eam_branches/ipp-20140226/psModules/src/objects/models/pmModel_SERSIC.c

@@ -50 +50 @@
 #include "pmSourceDiffStats.h"
 #include "pmSourceSatstar.h"
+#include "pmSourceLensing.h"
 #include "pmSource.h"
 #include "pmSourceFitModel.h"

branches/eam_branches/ipp-20140226/psModules/src/objects/models/pmModel_TRAIL.c

@@ -40 +40 @@
 #include "pmSourceDiffStats.h"
 #include "pmSourceSatstar.h"
+#include "pmSourceLensing.h"
 #include "pmSource.h"
 #include "pmSourceFitModel.h"

branches/eam_branches/ipp-20140226/psModules/src/objects/pmDetEff.c

@@ -24 +24 @@
 #include "pmSourceDiffStats.h"
 #include "pmSourceSatstar.h"
+#include "pmSourceLensing.h"
 #include "pmSource.h"
 #include "pmSourceFitModel.h"

branches/eam_branches/ipp-20140226/psModules/src/objects/pmGrowthCurve.c

@@ -37 +37 @@
 #include "pmSourceDiffStats.h"
 #include "pmSourceSatstar.h"
+#include "pmSourceLensing.h"
 #include "pmSource.h"
 #include "pmSourceFitModel.h"

branches/eam_branches/ipp-20140226/psModules/src/objects/pmGrowthCurveGenerate.c

@@ -44 +44 @@
 #include "pmSourceDiffStats.h"
 #include "pmSourceSatstar.h"
+#include "pmSourceLensing.h"
 #include "pmSource.h"
 #include "pmSourceFitModel.h"

branches/eam_branches/ipp-20140226/psModules/src/objects/pmModelUtils.c

@@ -39 +39 @@
 #include "pmSourceDiffStats.h"
 #include "pmSourceSatstar.h"
+#include "pmSourceLensing.h"
 #include "pmSource.h"
 #include "pmSourceFitModel.h"

branches/eam_branches/ipp-20140226/psModules/src/objects/pmMoments.h

@@ -44 +44 @@
     float Myyyy;   ///< fourth moment
 
+  // float wSum;    ///< window-weighted sum (NOT needed by lensing)
+
     float Sum;    ///< Pixel sum above sky (background).
     float Peak;   ///< Peak counts above sky.

branches/eam_branches/ipp-20140226/psModules/src/objects/pmPCM_MinimizeChisq.c

@@ -38 +38 @@
 #include "pmSourceDiffStats.h"
 #include "pmSourceSatstar.h"
+#include "pmSourceLensing.h"
 #include "pmSource.h"
 #include "pmSourceFitModel.h"

branches/eam_branches/ipp-20140226/psModules/src/objects/pmPCMdata.c

@@ -38 +38 @@
 #include "pmSourceDiffStats.h"
 #include "pmSourceSatstar.h"
+#include "pmSourceLensing.h"
 #include "pmSource.h"
 #include "pmSourceFitModel.h"

branches/eam_branches/ipp-20140226/psModules/src/objects/pmPSF.c

@@ -44 +44 @@
 #include "pmSourceDiffStats.h"
 #include "pmSourceSatstar.h"
+#include "pmSourceLensing.h"
 #include "pmSource.h"
 #include "pmSourceFitModel.h"

branches/eam_branches/ipp-20140226/psModules/src/objects/pmPSF_IO.c

@@ -54 +54 @@
 #include "pmSourceDiffStats.h"
 #include "pmSourceSatstar.h"
+#include "pmSourceLensing.h"
 #include "pmSource.h"
 #include "pmSourceFitModel.h"

branches/eam_branches/ipp-20140226/psModules/src/objects/pmPSFtry.c

@@ -36 +36 @@
 #include "pmSourceDiffStats.h"
 #include "pmSourceSatstar.h"
+#include "pmSourceLensing.h"
 #include "pmSource.h"
 #include "pmSourceFitModel.h"

branches/eam_branches/ipp-20140226/psModules/src/objects/pmPSFtryFitEXT.c

@@ -36 +36 @@
 #include "pmSourceDiffStats.h"
 #include "pmSourceSatstar.h"
+#include "pmSourceLensing.h"
 #include "pmSource.h"
 #include "pmSourceUtils.h"

branches/eam_branches/ipp-20140226/psModules/src/objects/pmPSFtryFitPSF.c

@@ -34 +34 @@
 #include "pmSourceDiffStats.h"
 #include "pmSourceSatstar.h"
+#include "pmSourceLensing.h"
 #include "pmSource.h"
 #include "pmSourceFitModel.h"

branches/eam_branches/ipp-20140226/psModules/src/objects/pmPSFtryMakePSF.c

@@ -35 +35 @@
 #include "pmSourceDiffStats.h"
 #include "pmSourceSatstar.h"
+#include "pmSourceLensing.h"
 #include "pmSource.h"
 #include "pmSourceFitModel.h"

branches/eam_branches/ipp-20140226/psModules/src/objects/pmPSFtryMetric.c

@@ -35 +35 @@
 #include "pmSourceDiffStats.h"
 #include "pmSourceSatstar.h"
+#include "pmSourceLensing.h"
 #include "pmSource.h"
 #include "pmSourceFitModel.h"

branches/eam_branches/ipp-20140226/psModules/src/objects/pmPSFtryModel.c

@@ -36 +36 @@
 #include "pmSourceDiffStats.h"
 #include "pmSourceSatstar.h"
+#include "pmSourceLensing.h"
 #include "pmSource.h"
 #include "pmSourceFitModel.h"

branches/eam_branches/ipp-20140226/psModules/src/objects/pmPhotObj.c

@@ -36 +36 @@
 #include "pmSourceDiffStats.h"
 #include "pmSourceSatstar.h"
+#include "pmSourceLensing.h"
 #include "pmSource.h"
 

branches/eam_branches/ipp-20140226/psModules/src/objects/pmSource.c

@@ -41 +41 @@
 #include "pmSourcePhotometry.h"
 #include "pmSourceSatstar.h"
+#include "pmSourceLensing.h"
 #include "pmSource.h"
 
@@ -68 +69 @@
     psFree(tmp->galaxyFits);
     psFree(tmp->radialAper);
+    psFree(tmp->lensingOBJ);
+    psFree(tmp->lensingPSF);
     psTrace("psModules.objects", 10, "---- end ----\n");
 }
@@ -166 +169 @@
     source->diffStats = NULL;
     source->galaxyFits = NULL;
+    source->lensingOBJ = NULL;
+    source->lensingPSF = NULL;
     source->radialAper = NULL;
     source->parent = NULL;
@@ -249 +254 @@
 
     source->region           = in->region;
+
+    // XXX I am not copying the pointers to things like the blends, satstar profile, galaxyFits, etc
 
     return(source);

branches/eam_branches/ipp-20140226/psModules/src/objects/pmSource.h

@@ -120 +120 @@
     pmSourceDiffStats *diffStats;       ///< extra parameters for difference detections
     pmSourceGalaxyFits *galaxyFits;     ///< fits to galaxy models (psphotFullForce only)
+    pmSourceLensing *lensingOBJ;        ///< lensing moments parameters (per object)
+    pmSourceLensing *lensingPSF;        ///< lensing moments parameters (psf, interpolated)
     psArray *radialAper;                ///< radial flux in circular apertures
     pmSource *parent;                   ///< reference to the master source from which this is derived

branches/eam_branches/ipp-20140226/psModules/src/objects/pmSourceContour.c

@@ -40 +40 @@
 #include "pmSourceDiffStats.h"
 #include "pmSourceSatstar.h"
+#include "pmSourceLensing.h"
 #include "pmSource.h"
 

branches/eam_branches/ipp-20140226/psModules/src/objects/pmSourceExtendedPars.h

@@ -118 +118 @@
 pmSourceGalaxyFits *pmSourceGalaxyFitsAlloc (void);
 
-
 /// @}
 # endif /* PM_SOURCE_H */

branches/eam_branches/ipp-20140226/psModules/src/objects/pmSourceFitModel.c

@@ -40 +40 @@
 #include "pmSourceDiffStats.h"
 #include "pmSourceSatstar.h"
+#include "pmSourceLensing.h"
 #include "pmSource.h"
 #include "pmSourcePhotometry.h"

branches/eam_branches/ipp-20140226/psModules/src/objects/pmSourceFitPCM.c

@@ -38 +38 @@
 #include "pmSourceDiffStats.h"
 #include "pmSourceSatstar.h"
+#include "pmSourceLensing.h"
 #include "pmSource.h"
 #include "pmSourcePhotometry.h"

branches/eam_branches/ipp-20140226/psModules/src/objects/pmSourceFitSet.c

@@ -39 +39 @@
 #include "pmSourceDiffStats.h"
 #include "pmSourceSatstar.h"
+#include "pmSourceLensing.h"
 #include "pmSource.h"
 #include "pmSourcePhotometry.h"

branches/eam_branches/ipp-20140226/psModules/src/objects/pmSourceGroups.c

@@ -25 +25 @@
 #include "pmSourceDiffStats.h"
 #include "pmSourceSatstar.h"
+#include "pmSourceLensing.h"
 #include "pmSource.h"
 #include "pmSourceFitModel.h"

branches/eam_branches/ipp-20140226/psModules/src/objects/pmSourceIO.c

@@ -47 +47 @@
 #include "pmSourceDiffStats.h"
 #include "pmSourceSatstar.h"
+#include "pmSourceLensing.h"
 #include "pmSource.h"
 #include "pmSourceFitModel.h"
@@ -597 +598 @@
             PM_SOURCES_WRITE("PS1_V3",    CMF_PS1_V3);
             PM_SOURCES_WRITE("PS1_V4",    CMF_PS1_V4);
+            PM_SOURCES_WRITE("PS1_V5",    CMF_PS1_V5);
             PM_SOURCES_WRITE("PS1_SV1",   CMF_PS1_SV1);
             PM_SOURCES_WRITE("PS1_SV2",   CMF_PS1_SV2);
@@ -1116 +1118 @@
             PM_SOURCES_READ_PSF("PS1_V3",    CMF_PS1_V3);
             PM_SOURCES_READ_PSF("PS1_V4",    CMF_PS1_V4);
+            PM_SOURCES_READ_PSF("PS1_V5",    CMF_PS1_V5);
             PM_SOURCES_READ_PSF("PS1_SV1",   CMF_PS1_SV1);
             PM_SOURCES_READ_PSF("PS1_SV2",   CMF_PS1_SV2);
@@ -1389 +1392 @@
         PM_SOURCES_READ_XSRC("PS1_V3",    CMF_PS1_V3);
         PM_SOURCES_READ_XSRC("PS1_V4",    CMF_PS1_V4);
+        PM_SOURCES_READ_XSRC("PS1_V5",    CMF_PS1_V5);
         PM_SOURCES_READ_XSRC("PS1_SV1",   CMF_PS1_SV1);
         PM_SOURCES_READ_XSRC("PS1_SV2",   CMF_PS1_SV2);
@@ -1428 +1432 @@
         PM_SOURCES_READ_XFIT("PS1_V3",    CMF_PS1_V3);
         PM_SOURCES_READ_XFIT("PS1_V4",    CMF_PS1_V4);
+        PM_SOURCES_READ_XFIT("PS1_V5",    CMF_PS1_V5);
         PM_SOURCES_READ_XFIT("PS1_SV1",   CMF_PS1_SV1);
         PM_SOURCES_READ_XFIT("PS1_SV2",   CMF_PS1_SV2);
@@ -1466 +1471 @@
         PM_SOURCES_READ_XRAD("PS1_V3",    CMF_PS1_V3);
         PM_SOURCES_READ_XRAD("PS1_V4",    CMF_PS1_V4);
+        PM_SOURCES_READ_XRAD("PS1_V5",    CMF_PS1_V5);
         PM_SOURCES_READ_XRAD("PS1_SV1",   CMF_PS1_SV1);
         PM_SOURCES_READ_XRAD("PS1_SV2",   CMF_PS1_SV2);
@@ -1504 +1510 @@
         PM_SOURCES_READ_XGAL("PS1_V3",    CMF_PS1_V3);
         PM_SOURCES_READ_XGAL("PS1_V4",    CMF_PS1_V4);
+        PM_SOURCES_READ_XGAL("PS1_V5",    CMF_PS1_V5);
         PM_SOURCES_READ_XGAL("PS1_SV1",   CMF_PS1_SV1);
         PM_SOURCES_READ_XGAL("PS1_SV2",   CMF_PS1_SV2);

branches/eam_branches/ipp-20140226/psModules/src/objects/pmSourceIO.h

@@ -37 +37 @@
 MK_PROTO(CMF_PS1_V3);
 MK_PROTO(CMF_PS1_V4);
+MK_PROTO(CMF_PS1_V5);
 MK_PROTO(CMF_PS1_SV1);
 MK_PROTO(CMF_PS1_SV2);

branches/eam_branches/ipp-20140226/psModules/src/objects/pmSourceIO_CFF.c

@@ -44 +44 @@
 #include "pmSourceDiffStats.h"
 #include "pmSourceSatstar.h"
+#include "pmSourceLensing.h"
 #include "pmSource.h"
 #include "pmSourceFitModel.h"

branches/eam_branches/ipp-20140226/psModules/src/objects/pmSourceIO_CMF.c.in

@@ -44 +44 @@
 #include "pmSourceDiffStats.h"
 #include "pmSourceSatstar.h"
+#include "pmSourceLensing.h"
 #include "pmSource.h"
 #include "pmSourceFitModel.h"
@@ -178 +179 @@
         @>PS1_V2,PS1_SV?@         psMetadataAdd (row, PS_LIST_TAIL, "MOMENTS_M4C",      PS_DATA_F32, "fourth momemt cos theta",                    moments.M_c4);
         @>PS1_V2,PS1_SV?@         psMetadataAdd (row, PS_LIST_TAIL, "MOMENTS_M4S",      PS_DATA_F32, "fourth momemt sin theta",                    moments.M_s4);
+
+        // Lensing parameters:
+        @>PS1_V4@                 psMetadataAdd (row, PS_LIST_TAIL, "LENS_X11_SM_OBJ",  PS_DATA_F32, "smear polarizability element (objects)",     source->lensingOBJ->smear->X11);
+        @>PS1_V4@                 psMetadataAdd (row, PS_LIST_TAIL, "LENS_X12_SM_OBJ",  PS_DATA_F32, "smear polarizability element (objects)",     source->lensingOBJ->smear->X12);
+        @>PS1_V4@                 psMetadataAdd (row, PS_LIST_TAIL, "LENS_X22_SM_OBJ",  PS_DATA_F32, "smear polarizability element (objects)",     source->lensingOBJ->smear->X22);
+        @>PS1_V4@                 psMetadataAdd (row, PS_LIST_TAIL, "LENS_E1_SM_OBJ",   PS_DATA_F32, "smear polarizability element (objects)",     source->lensingOBJ->smear->e1);
+        @>PS1_V4@                 psMetadataAdd (row, PS_LIST_TAIL, "LENS_E2_SM_OBJ",   PS_DATA_F32, "smear polarizability element (objects)",     source->lensingOBJ->smear->e2);
+
+        @>PS1_V4@                 psMetadataAdd (row, PS_LIST_TAIL, "LENS_X11_SH_OBJ",  PS_DATA_F32, "shear polarizability element (objects)",     source->lensingOBJ->shear->X11);
+        @>PS1_V4@                 psMetadataAdd (row, PS_LIST_TAIL, "LENS_X12_SH_OBJ",  PS_DATA_F32, "shear polarizability element (objects)",     source->lensingOBJ->shear->X12);
+        @>PS1_V4@                 psMetadataAdd (row, PS_LIST_TAIL, "LENS_X22_SH_OBJ",  PS_DATA_F32, "shear polarizability element (objects)",     source->lensingOBJ->shear->X22);
+        @>PS1_V4@                 psMetadataAdd (row, PS_LIST_TAIL, "LENS_E1_SH_OBJ",   PS_DATA_F32, "shear polarizability element (objects)",     source->lensingOBJ->shear->e1); 
+        @>PS1_V4@                 psMetadataAdd (row, PS_LIST_TAIL, "LENS_E2_SH_OBJ",   PS_DATA_F32, "shear polarizability element (objects)",     source->lensingOBJ->shear->e2); 
+
+        @>PS1_V4@                 psMetadataAdd (row, PS_LIST_TAIL, "LENS_X11_SM_PSF",  PS_DATA_F32, "smear polarizability element (objects)",     source->lensingPSF->smear->X11);
+        @>PS1_V4@                 psMetadataAdd (row, PS_LIST_TAIL, "LENS_X12_SM_PSF",  PS_DATA_F32, "smear polarizability element (objects)",     source->lensingPSF->smear->X12);
+        @>PS1_V4@                 psMetadataAdd (row, PS_LIST_TAIL, "LENS_X22_SM_PSF",  PS_DATA_F32, "smear polarizability element (objects)",     source->lensingPSF->smear->X22);
+        @>PS1_V4@                 psMetadataAdd (row, PS_LIST_TAIL, "LENS_E1_SM_PSF",   PS_DATA_F32, "smear polarizability element (objects)",     source->lensingPSF->smear->e1); 
+        @>PS1_V4@                 psMetadataAdd (row, PS_LIST_TAIL, "LENS_E2_SM_PSF",   PS_DATA_F32, "smear polarizability element (objects)",     source->lensingPSF->smear->e2); 
+
+        @>PS1_V4@                 psMetadataAdd (row, PS_LIST_TAIL, "LENS_X11_SH_PSF",  PS_DATA_F32, "shear polarizability element (objects)",     source->lensingPSF->shear->X11);
+        @>PS1_V4@                 psMetadataAdd (row, PS_LIST_TAIL, "LENS_X12_SH_PSF",  PS_DATA_F32, "shear polarizability element (objects)",     source->lensingPSF->shear->X12);
+        @>PS1_V4@                 psMetadataAdd (row, PS_LIST_TAIL, "LENS_X22_SH_PSF",  PS_DATA_F32, "shear polarizability element (objects)",     source->lensingPSF->shear->X22);
+        @>PS1_V4@                 psMetadataAdd (row, PS_LIST_TAIL, "LENS_E1_SH_PSF",   PS_DATA_F32, "shear polarizability element (objects)",     source->lensingPSF->shear->e1); 
+        @>PS1_V4@                 psMetadataAdd (row, PS_LIST_TAIL, "LENS_E2_SH_PSF",   PS_DATA_F32, "shear polarizability element (objects)",     source->lensingPSF->shear->e2); 
 
         @>PS1_V2,PS1_SV?,>PS1_DV1@ psMetadataAdd (row, PS_LIST_TAIL, "MOMENTS_R1",       PS_DATA_F32, "first radial moment",                        moments.Mrf);

branches/eam_branches/ipp-20140226/psModules/src/objects/pmSourceIO_CMP.c

@@ -44 +44 @@
 #include "pmSourceDiffStats.h"
 #include "pmSourceSatstar.h"
+#include "pmSourceLensing.h"
 #include "pmSource.h"
 #include "pmSourceFitModel.h"

branches/eam_branches/ipp-20140226/psModules/src/objects/pmSourceIO_MatchedRefs.c

@@ -44 +44 @@
 #include "pmSourceDiffStats.h"
 #include "pmSourceSatstar.h"
+#include "pmSourceLensing.h"
 #include "pmSource.h"
 #include "pmSourceFitModel.h"

branches/eam_branches/ipp-20140226/psModules/src/objects/pmSourceIO_OBJ.c

@@ -44 +44 @@
 #include "pmSourceDiffStats.h"
 #include "pmSourceSatstar.h"
+#include "pmSourceLensing.h"
 #include "pmSource.h"
 #include "pmSourceFitModel.h"

branches/eam_branches/ipp-20140226/psModules/src/objects/pmSourceIO_PS1_CAL_0.c

@@ -44 +44 @@
 #include "pmSourceDiffStats.h"
 #include "pmSourceSatstar.h"
+#include "pmSourceLensing.h"
 #include "pmSource.h"
 #include "pmSourceFitModel.h"

branches/eam_branches/ipp-20140226/psModules/src/objects/pmSourceIO_PS1_DEV_0.c

@@ -44 +44 @@
 #include "pmSourceDiffStats.h"
 #include "pmSourceSatstar.h"
+#include "pmSourceLensing.h"
 #include "pmSource.h"
 #include "pmSourceFitModel.h"

branches/eam_branches/ipp-20140226/psModules/src/objects/pmSourceIO_PS1_DEV_1.c

@@ -44 +44 @@
 #include "pmSourceDiffStats.h"
 #include "pmSourceSatstar.h"
+#include "pmSourceLensing.h"
 #include "pmSource.h"
 #include "pmSourceFitModel.h"

branches/eam_branches/ipp-20140226/psModules/src/objects/pmSourceIO_RAW.c

@@ -44 +44 @@
 #include "pmSourceDiffStats.h"
 #include "pmSourceSatstar.h"
+#include "pmSourceLensing.h"
 #include "pmSource.h"
 #include "pmSourceFitModel.h"

branches/eam_branches/ipp-20140226/psModules/src/objects/pmSourceIO_SMPDATA.c

@@ -44 +44 @@
 #include "pmSourceDiffStats.h"
 #include "pmSourceSatstar.h"
+#include "pmSourceLensing.h"
 #include "pmSource.h"
 #include "pmSourceFitModel.h"

branches/eam_branches/ipp-20140226/psModules/src/objects/pmSourceIO_SX.c

@@ -44 +44 @@
 #include "pmSourceDiffStats.h"
 #include "pmSourceSatstar.h"
+#include "pmSourceLensing.h"
 #include "pmSource.h"
 #include "pmSourceFitModel.h"

branches/eam_branches/ipp-20140226/psModules/src/objects/pmSourceMatch.c

@@ -25 +25 @@
 #include "pmSourceDiffStats.h"
 #include "pmSourceSatstar.h"
+#include "pmSourceLensing.h"
 #include "pmSource.h"
 #include "pmSourceFitModel.h"

branches/eam_branches/ipp-20140226/psModules/src/objects/pmSourceMoments.c

@@ -42 +42 @@
 #include "pmSourceDiffStats.h"
 #include "pmSourceSatstar.h"
+#include "pmSourceLensing.h"
 #include "pmSource.h"
 
@@ -435 +436 @@
         }
     }
+    // NOT needed : source->moments->wSum = Sum;
+
     source->moments->Mxx = XX/Sum;
     source->moments->Mxy = XY/Sum;

branches/eam_branches/ipp-20140226/psModules/src/objects/pmSourceOutputs.c

@@ -33 +33 @@
 #include "pmSourceDiffStats.h"
 #include "pmSourceSatstar.h"
+#include "pmSourceLensing.h"
 #include "pmSource.h"
 #include "pmSourceFitModel.h"

branches/eam_branches/ipp-20140226/psModules/src/objects/pmSourcePhotometry.c

@@ -40 +40 @@
 #include "pmSourceDiffStats.h"
 #include "pmSourceSatstar.h"
+#include "pmSourceLensing.h"
 #include "pmSource.h"
 #include "pmSourceFitModel.h"

branches/eam_branches/ipp-20140226/psModules/src/objects/pmSourcePlotApResid.c

@@ -42 +42 @@
 #include "pmSourceDiffStats.h"
 #include "pmSourceSatstar.h"
+#include "pmSourceLensing.h"
 #include "pmSource.h"
 #include "pmSourceFitModel.h"

branches/eam_branches/ipp-20140226/psModules/src/objects/pmSourcePlotMoments.c

@@ -45 +45 @@
 #include "pmSourceDiffStats.h"
 #include "pmSourceSatstar.h"
+#include "pmSourceLensing.h"
 #include "pmSource.h"
 #include "pmSourceFitModel.h"

branches/eam_branches/ipp-20140226/psModules/src/objects/pmSourcePlotPSFModel.c

@@ -46 +46 @@
 #include "pmSourceDiffStats.h"
 #include "pmSourceSatstar.h"
+#include "pmSourceLensing.h"
 #include "pmSource.h"
 #include "pmSourceFitModel.h"

branches/eam_branches/ipp-20140226/psModules/src/objects/pmSourceSky.c

@@ -41 +41 @@
 #include "pmSourceDiffStats.h"
 #include "pmSourceSatstar.h"
+#include "pmSourceLensing.h"
 #include "pmSource.h"
 

branches/eam_branches/ipp-20140226/psModules/src/objects/pmSourceUtils.c

@@ -41 +41 @@
 #include "pmSourceDiffStats.h"
 #include "pmSourceSatstar.h"
+#include "pmSourceLensing.h"
 #include "pmSource.h"
 

branches/eam_branches/ipp-20140226/psModules/src/objects/pmSourceVisual.c

@@ -23 +23 @@
 #include "pmSourceDiffStats.h"
 #include "pmSourceSatstar.h"
+#include "pmSourceLensing.h"
 #include "pmSource.h"
 #include "pmSourceFitModel.h"

branches/eam_branches/ipp-20140226/psModules/src/psmodules.h

@@ -135 +135 @@
 #include <pmSourceSatstar.h>
 #include <pmSourceDiffStats.h>
+#include <pmSourceLensing.h>
 #include <pmSource.h>
 #include <pmSourceFitModel.h>

branches/eam_branches/ipp-20140226/psphot/src/Makefile.am

@@ -208 +208 @@
         psphotMakeGrowthCurve.c        \
         psphotMagnitudes.c             \
+        psphotLensing.c        \
         psphotSetMaskBits.c            \
         psphotSkyReplace.c             \

branches/eam_branches/ipp-20140226/psphot/src/psphot.h

@@ -144 +144 @@
 bool            psphotMagnitudesReadout(pmConfig *config, psMetadata *recipe, const pmFPAview *view, pmReadout *readout, psArray *sources, pmPSF *psf);
 bool            psphotMagnitudes_Threaded (psThreadJob *job);
+
+bool            psphotLensing (pmConfig *config, const pmFPAview *view, const char *filerule);
+bool            psphotLensingReadout(pmConfig *config, const pmFPAview *view, const char *filerule, int index, psMetadata *recipe);
+bool            psphotLensingPSFtrendsReadout(pmConfig *config, const pmFPAview *view, const char *filerule, int index, psMetadata *recipe);
 
 bool            psphotEfficiency (pmConfig *config, const pmFPAview *view, const char *filerule);

branches/eam_branches/ipp-20140226/psphot/src/psphotReadout.c

@@ -340 +340 @@
     // calculate source magnitudes
     if (!psphotMagnitudes(config, view, filerule)) { // pass 1 (detections->allSources)
-      psErrorStackPrint(stderr, "Unable to do magnitudes.");
+        psErrorStackPrint(stderr, "Unable to do magnitudes.");
         psErrorClear();
     }
+
+    // calculate lensing parameters
+    if (!psphotLensing(config, view, filerule)) {
+        psErrorStackPrint(stderr, "Unable to do lensing parameters.");
+        psErrorClear();
+    }
+
     if (!psphotEfficiency(config, view, filerule)) { // pass 1
         psErrorStackPrint(stderr, "Unable to determine detection efficiencies from fake sources");
@@ -353 +360 @@
     // replace background in residual image
     if (!psphotSkyReplace (config, view, filerule)) { // pass 1
-      psErrorStackPrint(stderr, "Unable to replace sky");
-      psErrorClear();
+        psErrorStackPrint(stderr, "Unable to replace sky");
+        psErrorClear();
     }
 
     // drop the references to the image pixels held by each source
     if (!psphotSourceFreePixels (config, view, filerule)) { // pass 1
-      psErrorStackPrint(stderr, "Unable to free source pixels");
-      psErrorClear();
+        psErrorStackPrint(stderr, "Unable to free source pixels");
+        psErrorClear();
     }
 

@@ -12 +12 @@
 pmSourceIO_CMF_PS1_V1.v1.c
 pmSourceIO_CMF_PS1_V4.c
+pmSourceIO_CMF_PS1_V5.c
 pmSourceIO_CMF_PS1_SV1.c
 pmSourceIO_CMF_PS1_SV2.c