Changeset 33736 for branches/eam_branches/ipp-20120405/Ohana/src/relphot/src/StarOps.c

Timestamp:

Apr 9, 2012, 2:06:33 PM (14 years ago)

Author:

eugene

Message:

add Mstdev to secfilt; add some test code to MosaicOps.c

File:

• branches/eam_branches/ipp-20120405/Ohana/src/relphot/src/StarOps.c (modified) (2 diffs)

branches/eam_branches/ipp-20120405/Ohana/src/relphot/src/StarOps.c

@@ -376 +376 @@
         catalog[i].secfilt[Nsecfilt*j+Nsec].M  = stats.mean;
         catalog[i].secfilt[Nsecfilt*j+Nsec].Map = apstats.mean;
-        catalog[i].secfilt[Nsecfilt*j+Nsec].dM = MAX (stats.error, stats.sigma);
+        catalog[i].secfilt[Nsecfilt*j+Nsec].dM = stats.error;
+        catalog[i].secfilt[Nsecfilt*j+Nsec].Mstdev = 1000.0*stats.sigma; // Mstdev is in millimags (not enough space for more precision)
         catalog[i].secfilt[Nsecfilt*j+Nsec].Xm = (stats.Nmeas > 1) ? 100.0*log10(stats.chisq) : NAN_S_SHORT;
         catalog[i].secfilt[Nsecfilt*j+Nsec].Ncode = Ncode;
@@ -488 +489 @@
         /* use sigma or error in dM for output? */
         catalog[i].secfilt[Nsecfilt*j+Ns].M  = stats.mean;
-        catalog[i].secfilt[Nsecfilt*j+Ns].dM = MAX (stats.error, stats.sigma);
+        catalog[i].secfilt[Nsecfilt*j+Ns].dM = stats.error;
+        catalog[i].secfilt[Nsecfilt*j+Ns].Mstdev = 1000.0*stats.sigma; // Mstdev is in millimags (not enough space for more precision)
         catalog[i].secfilt[Nsecfilt*j+Ns].Xm = (stats.Nmeas > 1) ? 100.0*log10(stats.chisq) : NAN_S_SHORT;
         catalog[i].secfilt[Nsecfilt*j+Ns].Ncode = N;