linformat_df.lambda_theo = round(linformat_df.lambda_theo.values, 2) for objName in catalogue_df.loc[dz.idx_include].index: if objName in ['8', 'SHOC579']: local_reference = objName.replace('_', '-') quick_reference = catalogue_df.loc[objName].quick_index pdf_address = '/home/vital/Dropbox/Astrophysics/Thesis/images/{}_absEffect.png'.format( quick_reference) dz.create_pdfDoc(pdf_address, pdf_type='table') dz.pdf_insert_table(table_format='l' + 'ccc') group_dict = OrderedDict() # Make dict with all the objects lines dataframes ouput_folder = '{}{}/'.format(catalogue_dict['Obj_Folder'], objName) linelog_reducAddress = '{objfolder}{codeName}_WHT_linesLog_reduc.txt'.format( objfolder=ouput_folder, codeName=objName) linelog_emisAddress = '{objfolder}{codeName}_WHT_linesLog_emission_2nd.txt'.format( objfolder=ouput_folder, codeName=objName) reduc_linedf = dz.load_lineslog_frame(linelog_reducAddress) emission_linedf = dz.load_lineslog_frame(linelog_emisAddress) Hbeta_F = reduc_linedf.loc['H1_4861A'].line_Flux Hbeta_emisF = emission_linedf.loc['H1_4861A'].line_Flux
] #Import library object dz = Dazer() #Load observational data catalogue_df = dz.load_excel_DF( 'E:\\Dropbox\\Astrophysics\\Data\\WHT_observations\\WHT_Galaxies_properties.xlsx' ) #Define data to load pdf_address = 'E:\\Dropbox\\Astrophysics\\Papers\\Yp_BayesianMethodology\\source files\\tables\\reference_table_noPreamble' #Generate pdf #dz.create_pdfDoc(pdf_address, pdf_type='table') dz.pdf_insert_table(headers) #Quick indexing dz.quick_indexing(catalogue_df) # Sample objects excludeObjects = ['SHOC579', 'SHOC575_n2', '11', 'SHOC588', 'SDSS1', 'SHOC36' ] # SHOC579, SHOC575, SHOC220, SHOC588, SHOC592, SHOC036 sampleObjects = catalogue_df.loc[ dz.idx_include & ~catalogue_df.index.isin(excludeObjects)].index.values for i in range(sampleObjects.size): objName = sampleObjects[i] local_refenrence = objName.replace('_', '-')
del metals_list[metals_list.index('HeI_HI' + ext_data)] del metals_list[metals_list.index('Ymass_O' + ext_data)] del metals_list[metals_list.index('Ymass_S' + ext_data)] #Generate pdf # dz.create_pdfDoc(pdf_address, pdf_type='table') # dz.pdfDoc.packages.append(Package('nicefrac')) # dz.pdfDoc.packages.append(Package('pifont')) # dz.pdfDoc.append(NoEscape(r'\newcommand{\cmark}{\ding{51}}')) # dz.pdfDoc.append(NoEscape(r'\newcommand{\xmark}{\ding{55}}')) # ['OI_HI_emis2nd', 'NI_HI_emis2nd', 'SI_HI_emis2nd'] # catalogue_df[metal_x].notnull() #Set the pdf format dz.pdf_insert_table(headers_format) for objName in catalogue_df.loc[dz.idx_include].index: regressions_employed = [] for element in ['O', 'N', 'S']: validity_entry = catalogue_df.loc[objName, element + '_valid'] element_abundance_key = '{}I_HI_emis2nd'.format(element) element_abundance_check = isnull(catalogue_df.loc[objName, element_abundance_key]) print objName, element, element_abundance_check if element_abundance_check is False: if (validity_entry not in ['ignored', 'NO_excess', 'Wide Component']): regressions_employed.append(element) else: print 'Fallo', objName, element name_superscript = r'\textsuperscript{{{regrens}}}'.format(regrens = ', '.join(regressions_employed))
# Quick indexing dz.quick_indexing(catalogue_df) # Sample objects excludeObjects = [ 'SHOC579', 'SHOC575_n2', '11', 'SHOC588', 'SDSS3', 'SDSS1', 'SHOC36' ] # SHOC579, SHOC575, SHOC220, SHOC588, SHOC592, SHOC036 sampleObjects = catalogue_df.loc[ dz.idx_include & ~catalogue_df.index.isin(excludeObjects)].index.values # Generate pdf tableAddress = article_folder + 'modelParameters' # print('Creating table in {}'.format(tableAddress)) # dz.create_pdfDoc(tableAddress, pdf_type='table') # dz.pdfDoc.packages.append(Package('nicefrac')) dz.pdf_insert_table(headers_format, addfinalLine=False) dz.addTableRow(['HII Galaxy', '$(cm^{-3})$', '$(K)$', '$(K)$', r'$c(H\beta)$'], last_row=True) # Loop through the objects for i in range(sampleObjects.size): # Object references objName = sampleObjects[i] local_reference = objName.replace('_', '-') quick_reference = catalogue_df.loc[objName].quick_index print '- Treating object {}: {} {}'.format(i, objName, quick_reference) # Declare configuration file objectFolder = '{}{}/'.format( root_folder, objName) # '{}{}\\'.format(root_folder, objName)
dz = Dazer() #Read table data df = pd.read_excel( '/home/vital/Dropbox/Astrophysics/Thesis/notes/HII_galaxies_properties.xlsx', sheetname='Sheet1') #print df.iloc[5].Author.replace("'", "\textquotesingle") #Define data to load pdf_address = '/home/vital/Dropbox/Astrophysics/Thesis/tables/HII_galaxies_properties' headers = ['Parameters range'] #Generate pdf #dz.create_pdfDoc(pdf_address, pdf_type='table') dz.pdf_insert_table(headers, table_format='c') for reference in df.index: low_limit, parameter, upper_limit, unit = df.loc[reference, df.columns].values if not pd.isnull(unit): unit = [unit.split(',')[0], unit.split(',')[1] ] if ',' in unit else [unit, unit] else: unit = ['', ''] if not pd.isnull(upper_limit): entry = '${low}{unit0}$ < ${variable}$ < ${up}{unit1}$'.format( variable=parameter,
headers_dic['c'] = r'c' headers_dic['d'] = r'd' headers_dic['e'] = r'd' # Import functions dz = Dazer() specS = SpectraSynthesizer() # Declare data location article_folder = 'E:\\Dropbox\\Astrophysics\\Papers\\Yp_BayesianMethodology\\source files\\tables\\' # # Generate pdf tableAddress = article_folder + 'emissivityCoefficients' # dz.create_pdfDoc(tableAddress, pdf_type='table') # # dz.pdfDoc.packages.append(Package('nicefrac')) dz.pdf_insert_table(headers_dic.values()) default_lines = [ 'H1_4341A', 'O3_4363A', 'He1_4471A', 'He2_4686A', 'Ar4_4740A', 'O3_4959A', 'O3_5007A', 'He1_5876A', 'S3_6312A', 'N2_6548A', 'H1_6563A', 'N2_6584A', 'He1_6678A', 'S2_6716A', 'S2_6731A', 'Ar3_7136A', 'O2_7319A_b', 'S3_9069A', 'S3_9531A' ] # Loop through the objects for i in range(len(default_lines)): # Object references lineLabel = default_lines[i] lineCoeffs = specS.config[lineLabel + '_coeffs']
] for objName in catalogue_df.loc[dz.idx_include].index: local_reference = objName.replace('_', '-') quick_reference = catalogue_df.loc[objName].quick_index pdf_address = '/home/vital/Dropbox/Astrophysics/Thesis/tables/object_lines/{}_lineFluxes'.format( quick_reference) table_address = '/home/vital/Dropbox/Astrophysics/Papers/Yp_AlternativeMethods/supplementary material online/{}_lineFluxes.txt'.format( quick_reference) #dz.create_pdfDoc(pdf_address, pdf_type='table') dz.pdf_insert_table(table_format='l' + 'c' * (3 * obj_per_page)) group_dict = OrderedDict() # Make dict with all the objects lines dataframes ouput_folder = '{}{}/'.format(catalogue_dict['Obj_Folder'], objName) linelog_reducAddress = '{objfolder}{codeName}_WHT_linesLog_reduc.txt'.format( objfolder=ouput_folder, codeName=objName) linelog_emisAddress = '{objfolder}{codeName}_WHT_linesLog_emission_2nd.txt'.format( objfolder=ouput_folder, codeName=objName) reduc_linedf = dz.load_lineslog_frame(linelog_reducAddress) emission_linedf = dz.load_lineslog_frame(linelog_emisAddress) cHbeta = catalogue_df.loc[objName, cHbeta_type] dz.deredden_lines(reduc_linedf, reddening_curve=red_curve,