# HG changeset patch
# User Henry S. Thompson <ht@inf.ed.ac.uk>
# Date 1621977540 14400
# Node ID 672621ab4db422a944ea92fc474fc096d730be32
# Parent  56508a6033a9a6cbb7ef71b05bf83cb9cec4f18a
trim points we check to a circle per radius

diff -r 56508a6033a9 -r 672621ab4db4 python/cluster.py
--- a/python/cluster.py	Tue May 25 14:01:26 2021 -0400
+++ b/python/cluster.py	Tue May 25 17:19:00 2021 -0400
@@ -29,15 +29,20 @@
     print(l,file=outfile)
     ff=PPAT.split(l)
     (nr,ox,oy,oz)=intsMaybe(ff)
+    home=(float(ox),float(oy),float(oz))
     et=ff[9]
     l=infile.readline().rstrip()
     print(l,file=outfile)
     (orad,ymin,ymax)=intsMaybe(PPAT.split(l))
-    print(nr,ox,oy,oz,et,orad,ymin,ymax)
     _=infile.readline()
     for l in infile:
       found=False
-      q=[float(i) for i in l.split('\t')[2].split(',')]
+      (_,td,qd)=l.rstrip().split('\t')
+      q=[float(i) for i in qd.split(',')]
+      td=float(td)
+      if td>orad:
+        # Enforce a circular region, not square
+        break
       for c in cc:
         for p in c:
           if d(p,q)<=n: