Exemplo n.º 1
0
    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)