### find stars if os.path.exists(frame + ".fullcoo"): # worked = aries.readFullCoo(frame,useFWHM,useBkg,useThresh) worked = ao.readFullCooAllStars(frame, useFWHM, useBkg, useThresh) else: ao.findStars(frame, useFWHM, useBkg, useThresh) # worked = aries.readFullCoo(frame,useFWHM,useBkg,useThresh) worked = ao.readFullCooAllStars(frame, useFWHM, useBkg, useThresh) while worked[0] != True: useFWHM, useBkg, useThresh = worked os.system("rm -f " + frame + ".fullcoo") ao.findStars(frame, useFWHM, useBkg, useThresh) # worked = aries.readFullCoo(frame,useFWHM,useBkg,useThresh) worked = ao.readFullCooAllStars(frame, useFWHM, useBkg, useThresh) ### Now that we know where the stars are, we can check the FWHM estimate measFWHM, measX, measY = ao.getFWHM(frame + "//.red.sub", frame + ".coo", fwhmType=useFwhm) print frame + " Measured: X, Y (daophot); X, Y, FWHM (imexam)", worked[1], worked[2], measX, measY, measFWHM if 2 < measFWHM < 50: useFWHM = measFWHM #### NOTE: a seeing criteria is used in step6 # shiftListCoords[frame]=[worked[1],worked[2],measFWHM] shiftListCoords[frame] = [measX, measY, measFWHM] ### Now make shiftlist shiftListFile = "Shiftlist_" + properName[nn] + "_" + filt + ".txt" allObjFile = "allObj_" + properName[nn] + "_" + filt + "_use.txt" g = open(shiftListFile, "w") gAll = open(allObjFile, "w") try: refFrame = aries.refFrameForObj[properName[nn], useNight]
## The shiftlist, object list, and final image are saved in the main folder print "Copying shiftfile" os.system("cp Shiftlist_"+obj+"_"+filt+".txt "+objFilterDir) print "Copying list of frames" os.system("cp "+ff+" "+objFilterDir+ff) print "Copying final image" os.system("cp "+obj+"_"+filt+"_mp.fits "+objFilterDir+obj+"_"+filt+".fits") ## The intermediate steps get saved to a subdirectory print "Copying saved info..." os.system("cp "+obj+"_"+filt+"* "+savedDir) ## Also we should print the FWHM for our target in the combined image finalFrame = obj+"_"+filt+"_mp.fits" os.system("rm -f "+finalFrame+".fullcoo") print "Finding the FWHM for the final frame and exporting to a summary file" ao.findStars(finalFrame, aries.useFWHM, aries.useBkg ,aries.useThresh, extra="") worked = ao.readFullCooAllStars(finalFrame, aries.useFWHM, aries.useBkg, aries.useThresh, interactive=True) while worked[0] != True: useFWHM,useBkg,useThresh, useX, useY, useEdge = worked os.system("rm -f "+finalFrame+".fullcoo") ao.findStars(finalFrame, useFWHM, useBkg, useThresh,extra="") worked = ao.readFullCooAllStars(finalFrame, useFWHM, useBkg, useThresh, numXpixels=useX, numYpixels=useY, ignoreEdgePixels=useEdge, interactive=True) measFWHM, measX, measY = ao.getFWHM(finalFrame, finalFrame+".coo", fwhmType=aries.useFWHM) print obj,measX, measY, measFWHM, eval(measFWHM)*0.02 print >>g, obj,measX, measY, measFWHM, eval(measFWHM)*0.02 g.close()