def create(filename):
if re.search('\.fits$',filename) or re.search('\.fit$',filename):
return Catalog.create_from_fits(filename)
elif re.search('(\.P|\.P\.gz)',filename):
return load_object(filename)
else:
raise Exception("Unrecognized suffix in catalog file: %s"%(filename))
def create(filename):
if re.search('\.fits$', filename) or re.search('\.fit$', filename):
return Catalog.create_from_fits(filename)
elif re.search('(\.P|\.P\.gz)', filename):
return load_object(filename)
else:
raise Exception("Unrecognized suffix in catalog file: %s" %
(filename))
if args.iso0 is not None:
iso0 = np.loadtxt(args.iso0,unpack=True)
else:
iso0 = iso
if args.iso1 is not None:
iso1 = np.loadtxt(args.iso1,unpack=True)
else:
iso1 = iso
iso_spline = UnivariateSpline(np.log10(iso[0]),iso[1],k=1,s=0)
iso0_spline = UnivariateSpline(np.log10(iso0[0]),iso0[1],k=1,s=0)
iso1_spline = UnivariateSpline(np.log10(iso1[0]),iso1[1],k=1,s=0)
obj0 = load_object(args.files[0])
obj1 = load_object(args.files[1])
for k in sorted(obj0):
m0 = obj0[k]['flux']
ecenter = obj0[k]['ecenter']
m1 = obj1[k]['flux']
exp0 = obj0[k]['bexp']
exp1 = obj1[k]['bexp']
ms0 = m0.integrate(latrange=[-5.0,5.0])
ms1 = m1.integrate(latrange=[-5.0,5.0])
if args.iso0 is not None:
iso0 = np.loadtxt(args.iso0, unpack=True)
else:
iso0 = iso
if args.iso1 is not None:
iso1 = np.loadtxt(args.iso1, unpack=True)
else:
iso1 = iso
iso_spline = UnivariateSpline(np.log10(iso[0]), iso[1], k=1, s=0)
iso0_spline = UnivariateSpline(np.log10(iso0[0]), iso0[1], k=1, s=0)
iso1_spline = UnivariateSpline(np.log10(iso1[0]), iso1[1], k=1, s=0)
obj0 = load_object(args.files[0])
obj1 = load_object(args.files[1])
for k in sorted(obj0):
m0 = obj0[k]['flux']
ecenter = obj0[k]['ecenter']
m1 = obj1[k]['flux']
exp0 = obj0[k]['bexp']
exp1 = obj1[k]['bexp']
ms0 = m0.integrate(latrange=[-5.0, 5.0])
ms1 = m1.integrate(latrange=[-5.0, 5.0])
if ms0 == 0: continue