Scripts: flat.py

# Compute plane interpolating selected atoms and report distances from atoms to plane.
import chimera, numpy
from StructMeasure import Planes, plane

# Currently selected atoms
atoms = chimera.selection.currentAtoms()

# Draw the plane.
Planes.planeManager.createPlane('plane', atoms)

# Compute distances
coords = chimera.numpyArrayFromAtoms(atoms)
p = plane(coords)
distances = numpy.array([p.distance(a.coord()) for a in atoms])

# Print results to reply log
print 'Signed distances for %d points' % len(atoms)
print distances
print ('Unsigned average %.4g, max %.4g, min %.4g, standard deviation %.4g'
       % (abs(distances).mean(), distances.max(), distances.min(), distances.std()))