def test_calcAnglularDistance(self): src1 = Source(ra = 0.0 , dec = 0.0) src2 = Source(ra = 0.0 , dec = 0.0) sc = SourceConflicts() d = sc.calcAngularDistance(src1, src2) self.assertAlmostEqual(0.0, d, 3) # from Carl's report For GBT12B-364: # source NGC 1266 src1.ra = 0.85526951111895499 src1.dec = -0.042364958710075701 # source NGC1266, GBT10A-014 src2.ra = 0.85526223891373798 #sexHrs2rad('03:16:00.0') src2.dec = -0.0423630194553513 #sexDeg2rad('-02:25:37.0') d = sc.calcAngularDistance(src1, src2) self.assertAlmostEqual(0.0258519896867, rad2arcMin(d), 5) #print "deltaRa: ", rad2arcMin(abs(src2.ra - src1.ra)) #print "deltaDec: ", rad2arcMin(abs(src2.dec - src1.dec)) # source NGC 3665 src1.ra = 2.9876837023847602 src1.dec = 0.67653858425479396 # source SDSS J112346 ...; GBT08A-033 src2.ra = 2.9835247282213602 #sexHrs2rad('11:23:46') src2.dec = 0.67338971939598802 #sexDeg2rad('38:34:56') d = sc.calcAngularDistance(src1, src2) self.assertAlmostEqual(15.549314042, rad2arcMin(d), 5)