def head_is_close
import bend-prot@0.1.0.0/topology.bend as Topologyimport Base
import bend-prot@0.1.0.0/geom.bend as G
import bend-prot@0.1.0.0/protein.bend as P
32def head_is_close(h: P.Atom, +qp: G.Vec3, +min2: F32) -> Bool:
33 match h:
34 case P.Atom{serial, elem, +pos}:
35 F32.is_lt(G.Vec3.dist2(pos, qp), min2)