Changeset 34260 for trunk/Ohana/src/libdvo/src/dvo_catalog.c

Timestamp:

Jul 31, 2012, 4:02:00 PM (14 years ago)

Author:

eugene

Message:

add PS1_V4 schema; add set mean fluxes based on stacks; define the primary image ID (based on distance from center or boundary tree) and set mean flux from stacks based on that info; fix (uniquify) psps ID for stack detections; mean fluxes are in Jy; detection fluxes are instrumental counts per second

File:

• trunk/Ohana/src/libdvo/src/dvo_catalog.c (modified) (3 diffs)

trunk/Ohana/src/libdvo/src/dvo_catalog.c

@@ -55 +55 @@
   if (!strcasecmp (catformat, "PS1_V2"))          return (DVO_FORMAT_PS1_V2);
   if (!strcasecmp (catformat, "PS1_V3"))          return (DVO_FORMAT_PS1_V3);
+  if (!strcasecmp (catformat, "PS1_V4"))          return (DVO_FORMAT_PS1_V4);
   if (!strcasecmp (catformat, "PS1_REF"))         return (DVO_FORMAT_PS1_REF);
   return (DVO_FORMAT_UNDEF);
@@ -66 +67 @@
   if (!strcasecmp (catmode, "SPLIT")) return (DVO_MODE_SPLIT);
   return (DVO_MODE_UNDEF);
+}
+
+// init all data, or just catalog data
+void dvo_average_init (Average *average) {
+  average->R               = 0;
+  average->D               = 0;
+  average->dR              = 0;
+  average->dD              = 0;
+
+  average->uR              = 0;
+  average->uD              = 0;
+  average->duR             = 0;
+  average->duD             = 0;
+  average->P               = 0;
+  average->dP              = 0;
+
+  average->ChiSqAve        = 0.0;
+  average->ChiSqPM         = 0.0;
+  average->ChiSqPar        = 0.0;
+  average->Tmean           = 0;
+  average->Trange          = 0;
+
+  average->Xp              = 0;
+  average->Npos            = 0;
+
+  average->Nmeasure        = 0;
+  average->Nmissing        = 0;
+  average->Nextend         = 0;
+  average->measureOffset   = -1;
+  average->missingOffset   = -1;
+  average->extendOffset    = -1;
+
+  average->flags           = 0;
+  average->photFlagsUpper  = 0;
+  average->photFlagsLower  = 0;
+
+  average->objID           = 0;
+  average->catID           = 0;
+  average->extID           = 0;
+}
+
+// init all data, or just catalog data
+void dvo_secfilt_init (SecFilt *secfilt) {
+  secfilt->M           = NAN;
+  secfilt->Map         = NAN;
+  secfilt->Mkron       = NAN;
+  secfilt->dMkron      = NAN;
+  secfilt->dM          = NAN;
+  secfilt->Xm          = NAN_S_SHORT;
+
+  secfilt->FluxPSF     = NAN;
+  secfilt->dFluxPSF    = NAN;
+  secfilt->FluxKron    = NAN;
+  secfilt->dFluxKron   = NAN;
+
+  secfilt->flags       = 0;
+  secfilt->Ncode       = 0;
+  secfilt->Nused       = 0;
+
+  secfilt->M_20        = NAN_S_SHORT;
+  secfilt->M_80        = NAN_S_SHORT;
+
+  secfilt->ubercalDist = 1000;
+  secfilt->Mstdev      = NAN_S_SHORT;
+  secfilt->stackID     = 0;
+
+  secfilt->dummy       = 0;
+}
+
+// init all data, or just catalog data
+void dvo_measure_init (Measure *measure) {
+ measure->M         = NAN;
+ measure->dR        = NAN;
+ measure->dD        = NAN;
+ measure->M         = NAN;
+ measure->Mcal      = NAN;
+ measure->Map       = NAN;
+ measure->Mkron     = NAN;
+ measure->dMkron    = NAN;
+ measure->dM        = NAN;
+ measure->dMcal     = NAN;
+ measure->dt        = NAN;
+
+ measure->FluxPSF   = NAN;
+ measure->dFluxPSF  = NAN;
+ measure->FluxKron  = NAN;
+ measure->dFluxKron = NAN;
+
+ measure->airmass   = NAN;
+ measure->az        = NAN;
+
+ measure->Xccd      = NAN;
+ measure->Yccd      = NAN;
+
+ measure->Sky       = NAN;
+ measure->dSky      = NAN;
+
+ measure->t         = 0;
+ measure->t_msec    = 0;
+ measure->averef    = 0;
+
+ measure->detID     = 0;
+ measure->imageID   = 0;
+ measure->objID     = 0;
+ measure->catID     = 0;
+ measure->extID     = 0;
+
+ measure->psfQual   = NAN;
+ measure->psfChisq  = NAN;
+ measure->psfNdof   = 0;
+ measure->psfNpix   = 0;
+ measure->crNsigma  = NAN;
+ measure->extNsigma = NAN;
+
+ measure->FWx       = 0;
+ measure->FWy       = 0;
+ measure->theta     = 0;
+
+ measure->Mxx       = 0;
+ measure->Mxy       = 0;
+ measure->Myy       = 0;
+
+ measure->dXccd     = 0;
+ measure->dYccd     = 0;
+ measure->dRsys     = 0;
+
+ measure->posangle  = 0;
+ measure->pltscale  = NAN;
+
+ measure->photcode  = 0;
+ measure->dbFlags   = 0;
+ measure->photFlags = 0;
 }
 
@@ -600 +733 @@
 }
 
+