Fix Mike's buggy code

Bug in sigma clipping calibration stars fixed.
young438 · Aug 13, 2018 · 4482401 · 4482401
1 parent b8dd2e4
commit 4482401
Showing 1 changed file with 5 additions and 4 deletions.
diff --git a/maphot/maphot.py b/maphot/maphot.py
@@ -306,10 +306,11 @@
 # Calculate magnitude calibration factor
 magCalibArray = (finalCat[FILTER + 'MeanPSFMag_CFHT']
                  - finalCat[magKeyName])
-dmagCalibration = np.std(magCalibArray)
-sigmaclip = [np.abs(magCalibArray) < 3 * dmagCalibration]
-magCalibration = np.median(magCalibArray[sigmaclip])
-dmagCalibration = np.std(magCalibArray[sigmaclip])
+dmagCalibration = np.nanstd(magCalibArray)
+magCalibration = np.nanmedian(magCalibArray)
+sigmaclip = [np.abs(magCalibArray - magCalibration) < 3 * dmagCalibration]
+magCalibration = np.nanmedian(magCalibArray[sigmaclip])
+dmagCalibration = np.nanstd(magCalibArray[sigmaclip])
 
 # Correct the TNO magnitude and zero point
 finalTNOphotCFHT = (TNOPhot.magnitude - lineAperCorr + magCalibration,