location 4000 30000 4000 30000
lweight 2

data ../NGC5272R.XYVIQ
lines 1 99999
read {x 1 y 2 b 3 v 4 berr 7 verr 8 nb 10 nv 12 ffb 13 ffv 14 qb 15 qv 16 howb 17 howv 18 bmeta 22 vmeta 23}

set b_v = b - v 
set b_v_meta = bmeta - vmeta

window 1 1 1 1 
limits -3 4 -3 -15
limits -1 3 -3 -15
expand 0.8
ctype white
box 
xlabel B - V
ylabel B

ctype white
ptype 1 1
#poi b_v b 
poi b_v_meta bmeta
#poi b_v b if (nv==2 && nb==2)
#poi b_v_meta bmeta if (nv==2 && nb==2)
#poi b_v_meta bmeta if (nv==2 && nb==2 && qv<0.5 && qb<0.5)
#poi b_v_meta bmeta if (nv==2 && nb==2 && ffb<0.5 && ffr<0.5)
#poi b_v_meta bmeta if (nv==2 && nb==2 && qv<0.5 && qb<0.5 && ffb<0.5 && ffr<0.5)

ctype red
ptype 4 0
expand 2.2

# source 1; nothing inside 20 pix
set dr = ((x-1891.7)**2 + (y-3227.2)**2)**0.5
poi b_v_meta bmeta if (dr<20)

# source 2; nothing inside 20 pix
set dr = ((x-2119.9)**2 + (y-4214.6)**2)**0.5
poi b_v_meta bmeta if (dr<20)

# source 3
# SOMEWHAT BLUE STAR IN THE ERROR CIRCLE (0.3 mag blueward of MS) <10 pixels out (actually <2 pixels out)
set dr = ((x-2276.6)**2 + (y-3694.6)**2)**0.5
ctype blue
poi b_v_meta bmeta if (dr<10)
# IDing exact star...  this IS one H&P listed, but note that it's got several neighbors that didn't get fit.
poi b_v_meta bmeta if (dr<2)
ctype red

# source 4; nothing inside 20 pix
set dr = ((x-3277.7)**2 + (y-4871.2)**2)**0.5
poi b_v_meta bmeta if (dr<20)

# source 5; nothing inside 20 pix
set dr = ((x-3366.9)**2 + (y-3350.8)**2)**0.5
poi b_v_meta bmeta if (dr<20)

# source 6; nothing inside 20 pix
set dr = ((x-1904.3)**2 + (y-2394.8)**2)**0.5
poi b_v_meta bmeta if (dr<20)

# source 7
# this one has a SLIGHTLY blue star in/near error circle... not clear if significant; <10 pixels out
# and another slightly bluer one 15-20 pixels out...
set dr = ((x-2849.6)**2 + (y-3091.0)**2)**0.5
ctype green
poi b_v_meta bmeta if (dr<10)
# tracking down star... this IS one H&P identified... plausibly the counterpart, but only VERY SLIGHTLY blue...
set dr = ((x-2852)**2 + (y-3086)**2)**0.5
poi b_v_meta bmeta if (dr<2)
# tracking down other star that's further out...  DONE.  Not one H&P identified, and VERY crowded by brighter star in any case...
set dr = ((x-2850)**2 + (y-3108)**2)**0.5
poi b_v_meta bmeta if (dr<1)
ctype red

# source 8; nothing inside 20 pixels
set dr = ((x-3779.4)**2 + (y-4166.8)**2)**0.5
poi b_v_meta bmeta if (dr<20)

# source 9; nothing inside 20 pixels
set dr = ((x-1279.1)**2 + (y- 830.4)**2)**0.5
poi b_v_meta bmeta if (dr<20)

# source 10; nothing inside 20 pixels
set dr = ((x-2071.4)**2 + (y-4494.8)**2)**0.5
poi b_v_meta bmeta if (dr<20)

# source 11; nothing inside 20 pixels
set dr = ((x-2246.8)**2 + (y-1671.6)**2)**0.5
poi b_v_meta bmeta if (dr<20)

# source 12
# THIS ONE HAS A VERY RED STAR IN/NEAR THE ERROR CIRCLE (FOREGROUND M OR K STAR?), 10-15 pixels out
set dr = ((x-2273.9)**2 + (y-2805.6)**2)**0.5
ctype yellow
poi b_v_meta bmeta if (dr<15)
# tracking down exact star... this IS one of the stars H&P listed... WE HAVE AN ID! (not a very interesting one, though...)
set dr = ((x-2282)**2 + (y-2796)**2)**0.5
poi b_v_meta bmeta if (dr<2)
ctype red

# source 13
# this one is outside the FOV
set dr = ((x-2409.5)**2 + (y-5627.2)**2)**0.5
poi b_v_meta bmeta if (dr<20)

# source 14
# THIS ONE HAS A VERY BLUE STAR IN/NEAR THE ERROR CIRCLE, 10-15 pixels out
set dr = ((x-2594.7)**2 + (y-4316.8)**2)**0.5
ctype cyan
poi b_v_meta bmeta if (dr<15)
# locating exact star...  NOT one of the ones H&P thought could be the counterpart... (wrong direction after boresite correction)
set dr = ((x-2590)**2 + (y-4329)**2)**0.5
poi b_v_meta bmeta if (dr<2)
ctype red

# nothing after this produces any red squares... why? all out of FOV I guess?

# source 15
set dr = ((x-5541.6)**2 + (y- 189.6)**2)**0.5
poi b_v_meta bmeta if (dr<20)

# source 16
set dr = ((x-5647.6)**2 + (y-1188.0)**2)**0.5
poi b_v_meta bmeta if (dr<20)

# source 17
set dr = ((x-6230.8)**2 + (y-2509.6)**2)**0.5
poi b_v_meta bmeta if (dr<20)

# source 18
set dr = ((x-8427.1)**2 + (y- 587.7)**2)**0.5
poi b_v_meta bmeta if (dr<20)

# source 19
set dr = ((x+ 557.8)**2 + (y-5858.6)**2)**0.5
poi b_v_meta bmeta if (dr<20)

# source 20
set dr = ((x- 706.4)**2 + (y-5921.5)**2)**0.5
poi b_v_meta bmeta if (dr<20)

ctype white
expand 0.6
relocate (31000 800)
putlabel 4 $date

# Complete Chandra source list in Ata coords... (rounded to nearest 1/10th of a pixel)
# 
#  1891.7  3227.2   source1
#  2119.9  4214.6   source2
#  2276.6  3694.6   source3
#  3277.7  4871.2   source4
#  3366.9  3350.8   source5
#  1904.3  2394.8   source6
#  2849.6  3091.0   source7
#  3779.4  4166.8   source8
#  1279.1   830.4   source9
#  2071.4  4494.8   source10
#  2246.8  1671.6   source11
#  2273.9  2805.6   source12
#  2409.5  5627.2   source13
#  2594.7  4316.8   source14
#  5541.6   189.6   source15
#  5647.6  1188.0   source16
#  6230.8  2509.6   source17
#  8427.1   587.7   source18
#  -557.8  5858.6   source19
#   706.4  5921.5   source20