Scripts: findIncomplete.py

File findIncomplete.py, 565 bytes (added by meng, 10 years ago)
Line 
1import chimera
2incomplete = []
3for m in chimera.openModels.list(modelTypes=[chimera.Molecule]):
4 for r in m.residues:
5 tmplRes = chimera.restmplFindResidue(r.type, False, False)
6 if not tmplRes:
7 continue
8 for aname in tmplRes.atomsMap.keys():
9 if aname[0] == "H":
10 continue
11 if aname not in r.atomsMap:
12 incomplete.append(r)
13 break
14if incomplete:
15 chimera.replyobj.info("Incomplete residues: "
16 + ", ".join([str(r) for r in incomplete]) + "\n")
17 chimera.selection.setCurrent(incomplete)
18else:
19 chimera.replyobj.info("No incomplete residues\n")