From 16cc90a5f1a00d9dd7a5d1bb98b2e30909ebfe3e Mon Sep 17 00:00:00 2001 From: Abhijeet Anand Date: Tue, 25 Feb 2025 10:09:52 -0800 Subject: [PATCH] lam_edge_corrected --- qsoabsfind/absorberutils.py | 6 +++--- 1 file changed, 3 insertions(+), 3 deletions(-) diff --git a/qsoabsfind/absorberutils.py b/qsoabsfind/absorberutils.py index 65cbf28..595940c 100644 --- a/qsoabsfind/absorberutils.py +++ b/qsoabsfind/absorberutils.py @@ -547,7 +547,7 @@ def absorber_search_window(wavelength, residual, err_residual, zqso, absorber, m absorber (str): Options 'CIV', 'MgII'. min_wave (float): minimum wavelength edge (in Ang) max_wave (float): maximum wavelength edge (in Ang) - lam_edge_sep (float): separation from minimum wavelength, i.e. lam_min + lam_edge_sep, this is just to make sure that we avoid the very edge of the spectrum + lam_edge_sep (float): separation from minimum/maximum wavelength, i.e. lam_min +/- lam_edge_sep, this is just to make sure that we avoid the very edge of the spectrum verbose (bool, optional): If True will print time info. Default is False. Returns: @@ -559,12 +559,12 @@ def absorber_search_window(wavelength, residual, err_residual, zqso, absorber, m lam_CIV = lines['CIV_1549'] * (1 + zqso + lines['dz_start']) #redshifted from CIV emission lines lam_MgII = lines['MgI_2799'] * (1 + zqso - lines['dz_end']) #blueshifted MgII emission lines lam_start = max(min_wave, lam_CIV) + lam_edge_sep - lam_end = min(max_wave, lam_MgII) + lam_end = min(max_wave, lam_MgII) - lam_edge_sep elif absorber == 'CIV': dz = (lines['dv'] / speed_of_light) * (1 + zqso) lam_CIV = lines['CIV_1549'] * (1 + zqso + dz) lam_start = max(min_wave, 1310 * (1 + zqso)) + lam_edge_sep # This is from Cooksey et al 2013 - lam_end = min(lam_CIV, max_wave) + lam_end = min(lam_CIV, max_wave) - lam_edge_sep else: raise ValueError("Absorber must be 'CIV' or 'MgII'")