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")
