From bc80a0b94c1c4640680305dfb1dd989734dc3627 Mon Sep 17 00:00:00 2001 From: Chun Ly Date: Fri, 29 May 2020 20:02:33 -0700 Subject: [PATCH] analysis.attenuation.compute_EBV: Use filename_dict for input files (iss #10) --- Metallicity_Stack_Commons/analysis/attenuation.py | 7 ++++++- 1 file changed, 6 insertions(+), 1 deletion(-) 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