original_docstring: Experiment in progress, hence the underscore; expect changes.
docstring: Experiment in progress, hence the underscore; expect changes.
docstring_tokens: Experiment in progress hence the underscore ; expect changes .
code: def _plot_stars(catalog, observer, project, ax, mag1, mag2, margin=1.25):
      """Experiment in progress, hence the underscore; expect changes."""
  
      art = []
  
      # from astropy import wcs
      # w = wcs.WCS(naxis=2)
      # w.wcs.crpix = [-234.75, 8.3393]
      # w.wcs.cdelt = np.array([-0.066667, 0.066667])
      # w.wcs.crval = [0, -90]
      # w.wcs.ctype = ["RA---AIR", "DEC--AIR"]
      # w.wcs.set_pv([(2, 1, 45.0)])
  
      # import matplotlib.pyplot as plt
  
      # plt.subplot(projection=wcs)
      # #plt.imshow(hdu.data, vmin=-2.e-5, vmax=2.e-4, origin='lower')
      # plt.grid(color='white', ls='solid')
      # plt.xlabel('Galactic Longitude')
      # plt.ylabel('Galactic Latitude')
  
      xmin, xmax = ax.get_xlim()
      ymin, ymax = ax.get_xlim()
      lim = max(abs(xmin), abs(xmax), abs(ymin), abs(ymax)) * margin
      lims = (-lim, lim)
      ax.set_xlim(lims)
      ax.set_ylim(lims)
      ax.set_aspect('equal')
  
      o = observer[0]
  
      # Dim stars: points of with varying gray levels.
  
      c = catalog
      c = c[c['magnitude'] > mag1]
      c = c[c['magnitude'] <= mag2]
      #print('Second star group:', len(c))
      c = c.sort_values('magnitude', ascending=False)
      s = Star(ra_hours=c.ra_hours, dec_degrees=c.dec_degrees)
      spos = o.observe(s)
      x, y = project(spos)
      m = (mag2 - c['magnitude']) / (mag2 - mag1)
      # Note that "gray_r" is white for 0.0 and black for 1.0
      art.append(ax.scatter(
          x, y, s=1.0,
          c=1 - 0.8 * m, cmap='gray_r', vmin=0.0, vmax=1.0,
      ))
  
      # Bright stars: black circles of varying radius, surrounded by a
      # white gap in case stars are touching.  Draw the brightest stars
      # first to stop them from completely occluding smaller companions.
  
      def mag_to_radius(m):
          return (mag1 - m) * scale + 1.0
  
      c = catalog
      c = c[c['magnitude'] <= mag1]
      c = c.sort_values('magnitude', ascending=True)
      #print('First star group:', len(c))
      s = Star(ra_hours=c.ra_hours, dec_degrees=c.dec_degrees)
      spos = o.observe(s)
      x, y = project(spos)
      scale = 1.5
      radius = mag_to_radius(c['magnitude'])
  
      x2 = np.repeat(x, 2)
      y2 = np.repeat(y, 2)
      radius2 = (radius[:,None] + (3.0,0.0)).flatten()
      c2 = ('w', 'k')
      c2 = ('k', 'w')
  
      art.append(ax.scatter(x2, y2, s=radius2 ** 2.0, c=c2))
  
      return art, mag_to_radius
code_tokens: def _plot_stars ( catalog , observer , project , ax , mag1 , mag2 , margin = 1.25 ) : art = [ ] xmin , xmax = ax . get_xlim ( ) ymin , ymax = ax . get_xlim ( ) lim = max ( abs ( xmin ) , abs ( xmax ) , abs ( ymin ) , abs ( ymax ) ) * margin lims = ( - lim , lim ) ax . set_xlim ( lims ) ax . set_ylim ( lims ) ax . set_aspect ( 'equal' ) o = observer [ 0 ] c = catalog c = c [ c [ 'magnitude' ] > mag1 ] c = c [ c [ 'magnitude' ] <= mag2 ] c = c . sort_values ( 'magnitude' , ascending = False ) s = Star ( ra_hours = c . ra_hours , dec_degrees = c . dec_degrees ) spos = o . observe ( s ) x , y = project ( spos ) m = ( mag2 - c [ 'magnitude' ] ) / ( mag2 - mag1 ) art . append ( ax . scatter ( x , y , s = 1.0 , c = 1 - 0.8 * m , cmap = 'gray_r' , vmin = 0.0 , vmax = 1.0 , ) ) def mag_to_radius ( m ) : return ( mag1 - m ) * scale + 1.0 c = catalog c = c [ c [ 'magnitude' ] <= mag1 ] c = c . sort_values ( 'magnitude' , ascending = True ) s = Star ( ra_hours = c . ra_hours , dec_degrees = c . dec_degrees ) spos = o . observe ( s ) x , y = project ( spos ) scale = 1.5 radius = mag_to_radius ( c [ 'magnitude' ] ) x2 = np . repeat ( x , 2 ) y2 = np . repeat ( y , 2 ) radius2 = ( radius [ : , None ] + ( 3.0 , 0.0 ) ) . flatten ( ) c2 = ( 'w' , 'k' ) c2 = ( 'k' , 'w' ) art . append ( ax . scatter ( x2 , y2 , s = radius2 ** 2.0 , c = c2 ) ) return art , mag_to_radius
short_docstring: Experiment in progress, hence the underscore; expect changes.
short_docstring_tokens: Experiment in progress hence the underscore ; expect changes .
comment: """Experiment in progress, hence the underscore; expect changes."""$SEP$# from astropy import wcs$SEP$# w = wcs.WCS(naxis=2)$SEP$# w.wcs.crpix = [-234.75, 8.3393]$SEP$# w.wcs.cdelt = np.array([-0.066667, 0.066667])$SEP$# w.wcs.crval = [0, -90]$SEP$# w.wcs.ctype = ["RA---AIR", "DEC--AIR"]$SEP$# w.wcs.set_pv([(2, 1, 45.0)])$SEP$# import matplotlib.pyplot as plt$SEP$# plt.subplot(projection=wcs)$SEP$# #plt.imshow(hdu.data, vmin=-2.e-5, vmax=2.e-4, origin='lower')$SEP$# plt.grid(color='white', ls='solid')$SEP$# plt.xlabel('Galactic Longitude')$SEP$# plt.ylabel('Galactic Latitude')$SEP$# Dim stars: points of with varying gray levels.$SEP$#print('Second star group:', len(c))$SEP$# Note that "gray_r" is white for 0.0 and black for 1.0$SEP$# Bright stars: black circles of varying radius, surrounded by a$SEP$# white gap in case stars are touching.  Draw the brightest stars$SEP$# first to stop them from completely occluding smaller companions.$SEP$#print('First star group:', len(c))