FilesList = dz.Folder_Explorer(Pattern, Catalogue_Dic['Obj_Folder'], CheckComputer=False) Hbeta_values, Flux_values, names, sn_values, z_values = [], [], [], [], [] g_mags, r_mags = [], [] Declination_values = [] for i in range(len(FilesList)): #Get the frame row CodeName, FileName, FileFolder = dz.Analyze_Address(FilesList[i], verbose=False) #Load the observational data Hbeta_Flux = dz.GetParameter_ObjLog(CodeName, FileFolder, 'SDSS_Flux_Hbeta', Assumption='float') Hbeta_Ew = dz.GetParameter_ObjLog(CodeName, FileFolder, 'SDSS_Eqw_Hbeta', Assumption='float') SN_median = dz.GetParameter_ObjLog(CodeName, FileFolder, 'SDSS_SNmedian', Assumption='float') z_SDSS = dz.GetParameter_ObjLog(CodeName, FileFolder, 'SDSS_z', Assumption='float') g_mag = dz.GetParameter_ObjLog(CodeName, FileFolder,
#Get file list FilesList = dz.Folder_Explorer(emission_log, Catalogue_Dic['Obj_Folder'], CheckComputer=False) #Get the dictionary with the headers format and the data dz.EmissionLinesLog_header() #Generate list of objects (Dazer should have a method for this) for i in range(len(FilesList)): CodeName, FileName, FileFolder = dz.Analyze_Address(FilesList[i]) #load object data cHbeta = dz.GetParameter_ObjLog(CodeName, FileFolder, cHbeta_type, Assumption='float') obj_lines_frame = dz.load_object_frame(FileFolder, CodeName, emission_log, chbeta_coef=cHbeta_type) # obj_lines_frame_star = dz.load_object_frame(FileFolder, CodeName, emission_log_st, chbeta_coef = cHbeta_type) Hbeta_Flux = obj_lines_frame['line_Flux']['H1_4861A'] Hbeta_Int = obj_lines_frame['line_Int']['H1_4861A'] #Generate object row of data data_keys = dz.header_dict.keys() Object_row = ['None'] * len(dz.header_dict) #Establish table format dz.latex_header(table_address=FileFolder + CodeName + Table_Name,
Hbeta_label = 'H_beta' objects_which_failed = [] for i in range(len(FilesList)): 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")