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)
