diff --git a/Metallicity_Stack_Commons/analysis/attenuation.py b/Metallicity_Stack_Commons/analysis/attenuation.py index a5a4f2d..81d45f8 100644 --- a/Metallicity_Stack_Commons/analysis/attenuation.py +++ b/Metallicity_Stack_Commons/analysis/attenuation.py @@ -4,6 +4,7 @@ from os.path import join from .. import k_dict, line_name_short +from ..column_names import filename_dict HgHb_CaseB = 0.468 # Hg/Hb ratio for zero reddening HaHb_CaseB = 2.86 # Ha/Hb ratio for zero reddening @@ -15,7 +16,7 @@ k_HGAMMA = k_dict[HG] -def compute_EBV(fitspath, combine_asc): +def compute_EBV(fitspath): """ Purpose: Determines E(B-V) from Hg/Hb flux ratio @@ -24,6 +25,10 @@ def compute_EBV(fitspath, combine_asc): :param combine_asc: Astropy table containing emission-line flux """ + combine_file = join(fitspath, filename_dict['bin_fit_rev']) + print("Reading : " + combine_file) + combine_asc = asc.read(combine_file) + ID = combine_asc['ID'] HBETA = combine_asc[HB+'_Flux_Observed'].data HGAMMA = combine_asc[HG+'_Flux_Observed'].data