def number_negatives(self, source): #sources = filter(lambda src: src.getTag(tag), psources) tolerance = numpy.deg2rad(self.negregion * self.bmaj) if self.phaserad: radius = numpy.deg2rad(self.phaserad * self.bmaj) ra, dec = source.pos.ra, source.pos.dec # in radians within = self.nmodel.getSourcesNear(ra, dec, tolerance) if len(within) >= self.negthresh: if self.phaserad: if dist(self.ra0, self.dec0, ra, dec)[0] > radius: source.setTag(self.dd_tag, True) else: source.setTag(self.dd_tag, True)
def number_negatives(self): pmodel = Tigger.load(self.poscatalog, verbose=self.loglevel) nmodel = Tigger.load(self.negcatalog, verbose=self.loglevel) psources = pmodel.sources sources = filter(lambda src: src.getTag(self.high_corr_tag), psources) tolerance = numpy.deg2rad(self.negdetec_region * self.bmaj_deg) if self.phasecenter_excl_radius: radius = numpy.deg2rad(self.phasecenter_excl_radius * self.bmaj_deg) for srs in sources: ra, dec = srs.pos.ra, srs.pos.dec within = nmodel.getSourcesNear(ra, dec, tolerance) if len(within) >= self.negatives_thresh: if self.phasecenter_excl_radius: if dist( self.ra0, dec, ra, self.dec0)[0] > radius: srs.setTag(self.dd_tag, True) else: srs.setTag(self.dd_tag, True) pmodel.save(self.poscatalog)