Changeset 5958 for branches/eam_rel9_p0/psModules/src/objects/pmPSF.c

Timestamp:

Jan 9, 2006, 6:46:03 PM (21 years ago)

Author:

magnier

Message:

updated objects code with ApTrend concept (was in eam_rel9_b0)

File:

• branches/eam_rel9_p0/psModules/src/objects/pmPSF.c (modified) (3 diffs)

branches/eam_rel9_p0/psModules/src/objects/pmPSF.c

@@ -6 +6 @@
  *  @author EAM, IfA
  *
- *  @version $Revision: 1.3 $ $Name: not supported by cvs2svn $
- *  @date $Date: 2005-12-12 21:14:38 $
+ *  @version $Revision: 1.3.4.1 $ $Name: not supported by cvs2svn $
+ *  @date $Date: 2006-01-10 04:46:03 $
  *
  *  Copyright 2004 Maui High Performance Computing Center, University of Hawaii
@@ -69 +69 @@
  X-center
  Y-center
- ???: Sky background value?
- ???: Zo?
+ Sky background value
+ Object Normalization
  *****************************************************************************/
 pmPSF *pmPSFAlloc (pmModelType type)
@@ -83 +83 @@
     psf->dApResid = 0.0;
     psf->skyBias  = 0.0;
+
+    // the ApTrend components are (x, y, rflux)
+    psf->ApTrend = psPolynomial3DAlloc (2, 2, 1, PS_POLYNOMIAL_ORD);
+
+    // we do not allow any rflux vs X,Y cross-terms
+    psf->ApTrend->mask[0][1][1] = 1;  // rflux x^0 y^1
+    psf->ApTrend->mask[0][2][1] = 1;  // rflux x^0 y^2
+    psf->ApTrend->mask[1][0][1] = 1;  // rflux x^1 y^0
+    psf->ApTrend->mask[1][1][1] = 1;  // rflux x^1 y^1
+    psf->ApTrend->mask[1][2][1] = 1;  // rflux x^1 y^2
+    psf->ApTrend->mask[2][0][1] = 1;  // rflux x^2 y^0
+    psf->ApTrend->mask[2][1][1] = 1;  // rflux x^2 y^1
+    psf->ApTrend->mask[2][2][1] = 1;  // rflux x^2 y^2
+
+    // only 2nd order terms, no such combinations
+    psf->ApTrend->mask[2][2][0] = 1;  // x^2 y^2
+    psf->ApTrend->mask[2][1][0] = 1;  // x^2 y^1
+    psf->ApTrend->mask[1][2][0] = 1;  // x^1 y^2
+
+    // total free parameters = 18 - 11 = 7
 
     Nparams = pmModelParameterCount (type);