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
