Coordinate system used for healpy.query_disc()
Solved. Used the output of ang2vec() to give the vector.
The healpy.query_disc() function takes an argument vec which is a three-component unit vector defining the center of the disc. What coordinate system is being used here - why is there a third dimension for a 2-d projection? What point is the "tail" of this vector llocated at?