def calc_ebvs(ra, decl, w1, absw1): absw1 = absw1.mean() sig = w1 - absw1 DMmax = sig deltaDM = 0.1 DMarry = np.arange(deltaDM, DMmax + deltaDM, deltaDM) nsteps = len(DMarry) # number of steps in distance will vary dist_arry = 10**((DMarry + 5.) / 5.) / 1E3 #convert DM to kpc ra_arry = np.ones(nsteps) * ra dec_arry = np.ones(nsteps) * decl glats, glons = ext.eqToGal(ra_arry, dec_arry) ebvs = np.array(ext.getDustValue(glats, glons, dist_arry)) coeff = [0.34, 0.21, 0.15, 0.11, 0.09, 0.13] A_w1 = coeff[3] * 2.751 * ebvs A_w1_tests = sig - DMarry diff = abs(A_w1 - A_w1_tests) if len(diff) > 0: magic = np.where(diff == min(diff)) return ebvs[magic] else: return 0
def calc_ebvs(ra,decl,w1,absw1): absw1 = absw1.mean() sig = w1 - absw1 DMmax = sig deltaDM = 0.1 DMarry = np.arange(deltaDM,DMmax + deltaDM, deltaDM) nsteps = len(DMarry) # number of steps in distance will vary dist_arry = 10**((DMarry+5.)/5.)/1E3 #convert DM to kpc ra_arry = np.ones(nsteps)*ra dec_arry = np.ones(nsteps)*decl glats, glons = ext.eqToGal(ra_arry, dec_arry) ebvs = np.array(ext.getDustValue(glats, glons, dist_arry)) coeff = [0.34,0.21,0.15,0.11,0.09,0.13] A_w1 = coeff[3]*2.751*ebvs A_w1_tests = sig - DMarry diff = abs(A_w1 - A_w1_tests) if len(diff) > 0: magic = np.where(diff == min(diff)) return ebvs[magic] else: return 0
def return_ebv(ra,decl,DMarry): ra_arry = ra*np.pi/180. dec_arry = decl*np.pi/180. dist_arry = 10**((DMarry+5.)/5.)/1E3 #convert DM to kpc glats, glons = ext.eqToGal(ra_arry, dec_arry) ebvs = np.array(ext.getDustValue(glats, glons, dist_arry)) return ebvs
def return_ebv(ra, decl, DMarry): ra_arry = ra * np.pi / 180. dec_arry = decl * np.pi / 180. dist_arry = 10**((DMarry + 5.) / 5.) / 1E3 #convert DM to kpc glats, glons = ext.eqToGal(ra_arry, dec_arry) ebvs = np.array(ext.getDustValue(glats, glons, dist_arry)) return ebvs