double theta= Math.PI;
double phi = Math.PI;
double radius = Math.toRadians(0.2/3600.); // One arcse
long nside = pt.GetNSide(radius);
System.out.println(" calculated nside="+nside);
long cpix = pt.ang2pix_ring(nside,theta,phi);
Vector3d vc = pt.pix2vect_ring(nside, cpix);
ArrayList pixlist;
pixlist = pt.query_disc(nside, vc, radius, nest, inclusive);
int nlist = pixlist.size();