def calc_pos_angle(RA_DEC_hms_pri=['00:42.5 +41:12:00'],
RA_DEC_hms_sec=['00:41.5 +41:12:00']):
decimaldeg_pri = SkyCoord(RA_DEC_hms_pri,
unit=(units.hourangle, units.deg))
decimaldeg_sec = SkyCoord(RA_DEC_hms_sec,
unit=(units.hourangle, units.deg))
pos_angle = pyasl.positionAngle(decimaldeg_pri.ra.degree,
+decimaldeg_pri.dec.degree,
decimaldeg_sec.ra.degree,
+decimaldeg_sec.dec.degree)
return pos_angle
def calculate_PA(dataset, im):
PA = pyasl.positionAngle(dataset['%s_core'%im][0], dataset['%s_core'%im][1], dataset['%s_jet'%im][0], dataset['%s_jet'%im][1])
return PA
def position_angle(dataset, im):
PA = pyasl.positionAngle(dataset['%s_core'%im][0], dataset['%s_core'%im][1], dataset['%s_jet'%im][0], dataset['%s_jet'%im][1])
return PA
