#Get the object SDSS parameters
Name, Catalogue, mjd, plate, fiberID, alfalfa_code = Candiates_frame.index[i], row['catalogue'].values[0], row['mjd'].values[0], row['plate'].values[0], row['fiber'].values[0], row['alfalfa_name'].values[0]
#Generate object folder
CodeName = pv.generate_catalogue_tree(Catalogue_Dic, obj = Name)
FileFolder = Catalogue_Dic['Obj_Folder'] + CodeName + '/'
#Store parameter in object log file
pv.SetLogFile(CodeName + pv.ObjectLog_extension, FileFolder)
#Query the object table
#----------------Sloan objects----------------------
if Catalogue == 'sloan':
mjd, plate, fiberID = int(mjd), int(plate), int(fiberID)
obj_table = SDSS.query_specobj(mjd = mjd, plate = plate, fiberID=fiberID)
if obj_table != None:
#print Name, str(ephem.hours(math.radians(obj_table['ra'][0]))), str(ephem.degrees(math.radians(obj_table['dec'][0]))), '\n'
SDSS_RA = obj_table['ra'][0]
SDSS_DEC = obj_table['dec'][0]
SDSS_RA_hours = str(ephem.hours(math.radians(obj_table['ra'][0])))
SDSS_DEC_hours = str(ephem.degrees(math.radians(obj_table['dec'][0])))
co = coords.SkyCoord(float(SDSS_RA), float(SDSS_DEC), unit="deg")
Obj_query = SDSS.query_crossid(co, photoobj_fields=['modelMag_u', 'modelMag_g', 'modelMag_r'])
mag_u, mag_g, mag_r = Obj_query['modelMag_u'][0], Obj_query['modelMag_g'][0], Obj_query['modelMag_r'][0]
website = "http://dr12.sdss3.org/spectrumDetail?mjd={mjd}&fiber={fiber}&plateid={plateid}".format(mjd = mjd, fiber = fiberID, plateid = plate)
print Name, website
try:
#Get the frame row
CodeName, FileName_Blue, FileFolder = dz.Analyze_Address(FilesList[i])
#Store parameter in object log file
dz.SetLogFile(CodeName + log_file_extension, FileFolder)
#Query the object table
#----------------Sloan objects----------------------
mjd = dz.GetParameter_ObjLog(CodeName, FileFolder, 'SDSS_MJD', Assumption = 'float')
fiberID = dz.GetParameter_ObjLog(CodeName, FileFolder, 'SDSS_Fiber', Assumption = 'float')
plate = dz.GetParameter_ObjLog(CodeName, FileFolder, 'SDSS_Plate', Assumption = 'float')
obj_table = SDSS.query_specobj(mjd = mjd, plate = plate, fiberID=fiberID)
if obj_table != None:
#print Name, str(ephem.hours(math.radians(obj_table['ra'][0]))), str(ephem.degrees(math.radians(obj_table['dec'][0]))), '\n'
SDSS_RA = obj_table['ra'][0]
SDSS_DEC = obj_table['dec'][0]
SDSS_RA_hours = str(ephem.hours(math.radians(obj_table['ra'][0])))
SDSS_DEC_hours = str(ephem.degrees(math.radians(obj_table['dec'][0])))
co = coords.SkyCoord(float(SDSS_RA), float(SDSS_DEC), unit="deg")
Obj_query = SDSS.query_crossid(co, photoobj_fields=['modelMag_u', 'modelMag_g', 'modelMag_r'])
mag_u, mag_g, mag_r = Obj_query['modelMag_u'][0], Obj_query['modelMag_g'][0], Obj_query['modelMag_r'][0]
website = "http://dr12.sdss3.org/spectrumDetail?mjd={mjd}&fiber={fiber}&plateid={plateid}".format(mjd = mjd, fiber = fiberID, plateid = plate)
