############# Get 2MASS magnitude ############ print filt try: if filt == "J": ourMag = float(eval(settingsDict["2MASS_J"])[0]) elif filt == "H": ourMag = float(eval(settingsDict["2MASS_H"])[0]) elif filt == "Ks": ourMag = float(eval(settingsDict["2MASS_Ks"])[0]) else: ### We should get other catalogs in here! ourMag = "NA" except: ourMag = "NA" try: kepMag = kepler.kepMagCiardi(filt,ourMag) except: kepMag = "NA" ########################################################################### print "\n\n",annuli,skyBkgDict,"\n\n" ### Create output file limMagFile=ao.limMagFile(obj,filt,instrUsed) g = open(limMagFile,'w') print >>g,"File:",fitsFile #print >>g,"Camera-mode:",cameraMode print >>g,"Peak-counts:",peakPixelCount print >>g,"Sigma:",nSigma
def tableForAllStars(distCutoffArcsec=closeDist, filters =["J","Ks"], extra=""):
print "Making table of all stars within:",distCutoffArcsec," arcsec"
g = open(ao.compFileSorted(distCutoffArcsec, extra+".tsv"),"w")
gTex = open(ao.compFileSorted(distCutoffArcsec, extra+".tex"),"w")
## options for ARIES/Ks only
delimiter = "\t&"
print >>gTex, delimiter.join(["KNNNNN", "KeplerID", "KepMag", "Jmag", "Kmag","Filt","Star","Dist","PA","$\Delta_Mag$", "$Kep_{Dwarf}$", "\\\\"])
delimiter = "\t"
print >>g, delimiter.join(["KNNNNN", "KeplerID", "KepMag", "Jmag", "Kmag", "Filt", "Star","Dist","PA","Delta_Mag", "Kep_Dwarf"])
for obj in useObjects:
firstLineForObj = True # switch this when the first line is output
## read in the summary file
ss = open(ao.starsFile(obj,instr=instrUsed),"r")
lines = ss.readlines()
ss.close()
## find each filter
filtersFound=[]
for ii,line in enumerate(lines):
elems = line.split("\t")
if elems[0] in filters:
filtersFound.append(elems[0])
## read in the detected stars for each filter
fwhm = lines[1].rstrip().split()[2]
for line in lines:
elems = line.rstrip().split()
if elems[0] == 'J':
filt = 'J'
elif elems[0] == 'Ks':
filt = 'Ks'
## get rid of headers, also star 0
if elems[0] not in ["0","J","Ks","Star","object"]:
# print line, "xx", elems,"xx"
num = elems[0]
dist = elems[1]
deltaMag = elems[2]
xCoord = elems[3]
yCoord = elems[4]
pa = elems[5]
## estimate the kepler mag from K
try:
kMag = catalogs.getObjMag2MASS(obj,"Ks")[0]
jMag = catalogs.getObjMag2MASS(obj,"J")[0]
except:
jMag = "--"
kMag = "--"
# print "DIAG:",obj,kMag,jMag
try:
kepID = ksasDict[obj]["kepid"]
kepMag = ksasDict[obj]["kpmag"]
### There used to be a bug here (pre April 2013)! This should return the Kepler magnitude of the COMPANION, not the TARGET
if filt == "Ks":
kepMagDwarf = str(round(kepler.kepMagCiardi("Ks", eval(kMag)+eval(deltaMag)),1))
elif filt == "J":
kepMagDwarf = str(round(kepler.kepMagCiardi("J", eval(jMag)+eval(deltaMag)),1))
else:
kepMagDwarf = "--"
### Should allow for J AND K, since the transform is better, but matching stars is tricksy
except:
kepID = "--"
kepMag = "--"
kepMagDwarf = "--"
### Print info
if eval(dist) <= distCutoffArcsec:
if (firstLineForObj):
output = [obj, kepID, kepMag, jMag, kMag, filt, num,dist, pa, deltaMag, kepMagDwarf]
firstLineForObj = False ## for next star
else:
output = ["--", "--", "--","--", "--", filt, num, dist, pa, deltaMag, kepMagDwarf]
delimiter = "\t&"
print >>gTex, delimiter.join(output), "\\\\"
delimiter = "\t"
print >>g, delimiter.join(output)
gTex.close()
g.close()
