# dz.load_pymc_database(database_address)
# dz.plot_triangle_histContours(['He_abud', 'T_e', 'n_e', 'c_Hbeta', 'Tau', 'Xi'], [r'$y^{+}$', r'$T_{e}\,(K)$', r'$n_{e}\,(cm^{-3})$', r'$c(H\beta)$', r'$\tau$', r'$\xi$'], pymc_database=dz.pymc_database)
#
# dz.display_fig()
#Loop through files
for i in range(len(FilesList)):
version = '_real3_2000_100'
Folder_database = '/home/vital/Desktop/db_Testing/'
Db_name = 'db_Continuum' + version
Db_global_file_name = 'db_Continuum_Global' + version
#Analyze file address
CodeName, FileName, FileFolder = dz.Analyze_Address(FilesList[i])
lineslog_frame = dz.load_object_frame(FileFolder, CodeName, AbundancesFileExtension)
#Get reduce spectrum data
Wave_T, Int_T, ExtraData_T = dz.File_to_data(FileFolder,'obj' + CodeName + '_' + Catalogue_Dic['Datatype'] + '.fits')
Wave_N, Int_N, ExtraData_T = dz.File_to_data(FileFolder, CodeName + nebular_exten)
Wave_S, Int_S, ExtraData_N = dz.File_to_data(FileFolder, CodeName + Stellar_ext)
Wave_E, Int_E, ExtraData_N = dz.File_to_data(FileFolder, CodeName + emitting_ext)
#Load initial conditions
Te = dz.GetParameter_ObjLog(CodeName, FileFolder, 'TSIII_pn')
nHeII_HII = dz.GetParameter_ObjLog(CodeName, FileFolder, 'HeII_HII_from_O_pn')
nHeIII_HII = dz.GetParameter_ObjLog(CodeName, FileFolder, 'HeIII_HII_from_O_pn')
Halpha_Flux = lineslog_frame.loc['H1_6563A', 'line_Flux']
Hbeta_Flux = lineslog_frame.loc['H1_4861A', 'line_Flux']
obs_idx = (3550 < Wave_E) & (Int_E < 7000)
#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,
TitleColumn=['', '', 'HII Galaxy', CodeName, ''])
#Loop through all the emission recordings
objlines = obj_lines_frame.index.values
# emission_log_st = '_' + Catalogue_Dic['Datatype'] + '_emission_LinesLog_v3.txt'
#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, TitleColumn = ['', '', 'HII Galaxy', CodeName, ''])
#Loop through all the emission recordings
objlines = obj_lines_frame.index.values
for j in range(len(objlines)):
#Load pyneb data
Atom_dict = load_pyneb()
#Locate files on hard drive
FilesList = dz.Folder_Explorer(log_extension, Catalogue_Dic['Obj_Folder'], CheckComputer=False)
#Load catalogue properties frame
catalogue_frame = dz.load_catalogue_frame(FilesList)
#Loop through files to load frames into a dictionary:
Frames_dict = OrderedDict()
cHbeta_dict = OrderedDict()
for i in range(len(FilesList)):
CodeName, FileName, FileFolder = dz.Analyze_Address(FilesList[i])
Frames_dict[CodeName] = dz.load_object_frame(FileFolder, CodeName, emission_log, chbeta_coef = cHbeta_type)
cHbeta_dict[CodeName] = dz.GetParameter_ObjLog(CodeName, FileFolder, cHbeta_type, Assumption='float')
# #Dictionary with the variables we want to compare
variables_dict = OrderedDict()
# variables_dict['y_plus_inf'] = 'HeII_HII'
# variables_dict['Te_inf'] = 'TSIII'
# variables_dict['ne_inf'] = 'nSII'
# variables_dict['cHbeta_inf'] = 'cHBeta_stellar'
#Locate files on hard drive
FilesList = OrganizedObject_List()
#Get the dictionary with the headers format and the data
dz.Catalogue_headers()
