code: def get_starsep_RaDecDeg(ra1_deg, dec1_deg, ra2_deg, dec2_deg):
      
      sep = deltaStarsRaDecDeg(ra1_deg, dec1_deg, ra2_deg, dec2_deg)
      sgn, deg, mn, sec = degToDms(sep)
      if deg != 0:
          txt = '%02d:%02d:%06.3f' % (deg, mn, sec)
      else:
          txt = '%02d:%06.3f' % (mn, sec)
      return txt
docs: Calculate separation.