initial commit

tangostationcontrol/tangostationcontrol/beam/old_test.py

+import datetime
+from delays import *
+
+if __name__ == '__main__':
+    # # create a frame tied to the reference position
+    reference_itrf = [3826577.066, 461022.948, 5064892.786] # CS002LBA, in ITRF2005 epoch 2012.5
+    d = delay_calculator(reference_itrf)
+
+    # # set the timestamp to solve for
+    timestamp = datetime.datetime(2021,1,1,0,0,5)
+    d.set_measure_time(timestamp)
+
+    # compute the delays for an antennas w.r.t. the reference position
+    antenna_itrf = [[3826923.546, 460915.441, 5064643.489]] # CS001LBA, in ITRF2005 epoch 2012.5
+
+    # # obtain the direction vector for a specific pointing
+    direction = "J2000","0deg","0deg"
+
+    # calculate the delays based on the set reference position, the set time and now the set direction and antenna positions.
+    delays = d.convert(direction, antenna_itrf)
+
+    # print the delays
+    # pprint.pprint(delays)
+
+
+    #test changing the time
+
+    print(f"Changing timestamp test\nBase parametres: Direction: {direction}, position: {antenna_itrf}")
+    old_val = 0
+
+    for i in range(10):
+        # # set the timestamp to solve for
+        timestamp = datetime.datetime(2021,1,1,0,i,5)
+        d.set_measure_time(timestamp)
+
+        old_val = delays
+        delays = d.convert(direction, antenna_itrf)
+        diff = delays[0] - old_val
+
+        # print the delays
+        print(f"Timestamp: {timestamp}:   {delays[0]}, diff: {diff}")
+
+
+    # reset time
+    timestamp = datetime.datetime(2021, 1, 1, 0, 0, 5)
+    d.set_measure_time(timestamp)
+
+
+    #test changing the antenna position
+    print(f"Changing Antenna position test.\nBase parametres: Time: {timestamp} Direction: {direction}")
+    old_val = 0
+
+    for i in range(10):
+        antenna_itrf = [[3826577.066 + i, 461022.948, 5064892.786]]  # CS002LBA, in ITRF2005 epoch 2012.5
+
+        old_val = delays
+        delays = d.convert(direction, antenna_itrf)
+        diff = delays[0] - old_val
+
+        # print the delays
+        print(f"Antenna position: {antenna_itrf}:   {delays}, diff: {diff}")
+
+    # test changing the direction
+
+    antenna_itrf = [[3826923.546, 460915.441, 5064643.489]]  # CS001LBA, in ITRF2005 epoch 2012.5
+    # antenna_itrf = [[0, 0, 0]]  # CS001LBA, in ITRF2005 epoch 2012.5
+
+    print(f"Changing direction test.\nBase parametres: Time: {timestamp} , position: {antenna_itrf}")
+
+    old_val = 0
+
+    for i in range(12):
+        deg = i * 30
+        direction = "J2000", f"0deg", f"{deg}deg"
+
+        old_val = delays
+        delays = d.convert(direction, antenna_itrf)
+        diff = delays[0] - old_val
+
+        # print the delays
+        print(f"Direction: {direction}:  {delays}, diff: {diff}")

tangostationcontrol/tangostationcontrol/beam/test_delays.py

 from delays import *
 
+from tangostationcontrol.test import base
+
+import asyncio
+import mock
+
+
+
+
+
+class TestAttributeTypes(base.TestCase):
+    def setUp(self):
+        # # create a frame tied to the reference position
+        reference_itrf = [3826577.066, 461022.948, 5064892.786]  # CS002LBA, in ITRF2005 epoch 2012.5
+        d = delay_calculator(reference_itrf)
+
+        # # set the timestamp to solve for
+        timestamp = datetime.datetime(2021, 1, 1, 0, 0, 5)
+        d.set_measure_time(timestamp)
+
+        # compute the delays for an antennas w.r.t. the reference position
+        antenna_itrf = [[3826923.546, 460915.441, 5064643.489]]  # CS001LBA, in ITRF2005 epoch 2012.5
+
+        # # obtain the direction vector for a specific pointing
+        direction = "J2000", "0deg", "0deg"
+
+
+    def test_time_change(self):
+        # # set the timestamp to solve for
+        timestamp = datetime.datetime(2021, 1, 1, 0, i, 5)
+        d.set_measure_time(timestamp)
+
+        delays = d.convert(direction, antenna_itrf)
+
+        # print the delays
+        print(f"Timestamp: {timestamp}:   {delays}")
+
+    def test_direction_change(self):
+        # test changing the direction
+
+        antenna_itrf = [[3826923.546, 460915.441, 5064643.489]]  # CS001LBA, in ITRF2005 epoch 2012.5
+        print(f"Changing direction test.\nBase parametres: Time: {timestamp} , position: {antenna_itrf}")
+
+        for i in range(10):
+            direction = "J2000", f"{i}deg", "0deg"
+
+            delays = d.convert(direction, antenna_itrf)
+
+            # print the delays
+            print(f"Direction: {direction}:  {delays}")
+
+    def test_position_change(self):
+        # test changing the antenna position
+        print(f"Changing Antenna position test.\nBase parametres: Time: {timestamp} Direction: {direction}")
+        for i in range(10):
+            antenna_itrf = [[3826577.066 + i, 461022.948, 5064892.786]]  # CS002LBA, in ITRF2005 epoch 2012.5
+
+            delays = d.convert(direction, antenna_itrf)
+
+            # print the delays
+            print(f"Antenna position: {antenna_itrf}:   {delays}")
+
+
 if __name__ == '__main__':
     # # create a frame tied to the reference position
     reference_itrf = [3826577.066, 461022.948, 5064892.786] # CS002LBA, in ITRF2005 epoch 2012.5