define radius 1.49

macro read sdsscoords.sm

data "tileOutChunk20.par"
lines 12 0
read <tile_ra_in 3 tile_dec_in 4>
eqsdss tile_ra_in tile_dec_in tile_ra tile_dec 
set tile_dec = (tile_dec<180.) ? tile_dec : tile_dec-360.
define ntiles (dimen(tile_ra_in))

define table_type bintable
table "tiOutputTargetsChunk20.fit"
read table 'byname' {ra dec mask tileId} 
eqsdss ra dec ra_in dec_in 
set ra = ra_in
#
# set ra = (ra < 180) ? ra : ra - 360
#
set dec = (dec_in<180.) ? dec_in : dec_in-360.

macro   domissing {
        set raout = ra if (tileId == -1 && (mask&2) == 2)
        set decout = dec if (tileId == -1 && (mask&2) == 2)
        ctype green     
        ptype 4 1 
        expand 2
  points raout decout
        expand 1.0001
        ctype default
}

macro   docollided {
        set raout = ra if (tileId == -1 && (mask&2) == 0)
        set decout = dec if (tileId == -1 && (mask&2) == 0)
        ctype magenta   
        ptype 4 1 
        expand 0.999
  points raout decout
        expand 1.0001
        ctype default
}

macro circles {
  define x0 (tile_ra[j])
  define y0 (tile_dec[j])
  set k = 0, 40
  set x1 = $x0 + $radius*cos(k*2.*3.14159265/40.)/cos($y0*3.14159/180.)
  set y1 = $y0 + $radius*sin(k*2.*3.14159265/40.)
        ctype red
  connect x1 y1
        ctype default
}

#limits ra 148 182
limits ra -182 -163
erase
box 1 2 0 0
ptype 1 1
points ra dec
#docollided
domissing
ctype red
lweight 1
for j=0 j<$ntiles j=j+1 circles
lweight 1
ctype default