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,
