def _build_src_dicts_from_peaks(self, peaks, maps, src_dict_template):
tsmap = maps['ts']
amp = maps['amplitude']
src_dicts = []
names = []
for p in peaks:
o = fit_error_ellipse(tsmap, (p['ix'], p['iy']), dpix=2)
skydir = o['skydir']
p['fit_loc'] = o
p['fit_skydir'] = o['skydir']
p.update(o)
name = utils.create_source_name(skydir)
src_dict = copy.deepcopy(src_dict_template)
norm_par = get_function_norm_par_name(
src_dict_template['SpectrumType'])
src_dict.update({
norm_par: amp.data[p['iy'], p['ix']],
'ra': skydir.icrs.ra.deg,
'dec': skydir.icrs.dec.deg
})
src_dict['glon_err'] = o['glon_err']
src_dict['glat_err'] = o['glat_err']
src_dict['ra_err'] = o['ra_err']
src_dict['dec_err'] = o['dec_err']
src_dict['pos_err'] = o['pos_err']
src_dict['pos_err_semimajor'] = o['pos_err_semimajor']
src_dict['pos_err_semiminor'] = o['pos_err_semiminor']
src_dict['pos_r68'] = o['pos_r68']
src_dict['pos_r95'] = o['pos_r95']
src_dict['pos_r99'] = o['pos_r99']
src_dict['pos_angle'] = o['pos_angle']
src_dict['pos_gal_cov'] = o['pos_gal_cov']
src_dict['pos_gal_corr'] = o['pos_gal_corr']
src_dict['pos_cel_cov'] = o['pos_cel_cov']
src_dict['pos_cel_corr'] = o['pos_cel_corr']
self.logger.info('Found source\n' + 'name: %s\n' % name +
'ts: %f' % p['amp']**2)
names.append(name)
src_dicts.append(src_dict)
return names, src_dicts
def _build_src_dicts_from_peaks(self, peaks, maps, src_dict_template):
tsmap = maps['ts']
amp = maps['amplitude']
src_dicts = []
names = []
for p in peaks:
o = fit_error_ellipse(tsmap, (p['ix'], p['iy']), dpix=2)
skydir = o['skydir']
p['fit_loc'] = o
p['fit_skydir'] = o['skydir']
p.update(o)
name = utils.create_source_name(skydir)
src_dict = copy.deepcopy(src_dict_template)
norm_par = get_function_norm_par_name(
src_dict_template['SpectrumType'])
src_dict.update({norm_par: amp.counts[p['iy'], p['ix']],
'ra': skydir.icrs.ra.deg,
'dec': skydir.icrs.dec.deg})
src_dict['glon_err'] = o['glon_err']
src_dict['glat_err'] = o['glat_err']
src_dict['ra_err'] = o['ra_err']
src_dict['dec_err'] = o['dec_err']
src_dict['pos_err'] = o['pos_err']
src_dict['pos_err_semimajor'] = o['pos_err_semimajor']
src_dict['pos_err_semiminor'] = o['pos_err_semiminor']
src_dict['pos_r68'] = o['pos_r68']
src_dict['pos_r95'] = o['pos_r95']
src_dict['pos_r99'] = o['pos_r99']
src_dict['pos_angle'] = o['pos_angle']
src_dict['pos_gal_cov'] = o['pos_gal_cov']
src_dict['pos_gal_corr'] = o['pos_gal_corr']
src_dict['pos_cel_cov'] = o['pos_cel_cov']
src_dict['pos_cel_corr'] = o['pos_cel_corr']
self.logger.info('Found source\n' +
'name: %s\n' % name +
'ts: %f' % p['amp'] ** 2)
names.append(name)
src_dicts.append(src_dict)
return names, src_dicts
def _build_src_dicts_from_peaks(self, peaks, maps, src_dict_template):
tsmap = maps['ts']
amp = maps['amplitude']
src_dicts = []
names = []
for p in peaks:
o, skydir = fit_error_ellipse(tsmap, (p['ix'], p['iy']), dpix=2)
p['fit_loc'] = o
p['fit_skydir'] = skydir
p.update(o)
if o['fit_success']:
skydir = p['fit_skydir']
else:
skydir = p['skydir']
name = utils.create_source_name(skydir)
src_dict = copy.deepcopy(src_dict_template)
src_dict.update({
'Prefactor': amp.counts[p['iy'], p['ix']],
'ra': skydir.icrs.ra.deg,
'dec': skydir.icrs.dec.deg
})
src_dict['pos_sigma'] = o['sigma']
src_dict['pos_sigma_semimajor'] = o['sigma_semimajor']
src_dict['pos_sigma_semiminor'] = o['sigma_semiminor']
src_dict['pos_r68'] = o['r68']
src_dict['pos_r95'] = o['r95']
src_dict['pos_r99'] = o['r99']
src_dict['pos_angle'] = np.degrees(o['theta'])
self.logger.info('Found source\n' + 'name: %s\n' % name +
'ts: %f' % p['amp']**2)
names.append(name)
src_dicts.append(src_dict)
return names, src_dicts