Commit 842a671a authored by Maria Chiara Storer's avatar Maria Chiara Storer
Browse files

Signature modification to Modeller.addExtraParticles()

Within Modeller.addExtraParticles(), add the ability to pass down
an option to modify compiled.matchResidueToTemplate() 's
ignoreExternalBonds parameter.

https://simtk.org/plugins/phpBB/viewtopicPhpbb.php?f=161&t=11264&p=0&start=0&view=&sid=e384beeb08f3bfd767dcfb546006ad63
parent 75c1fcb6
...@@ -1004,7 +1004,7 @@ class Modeller(object): ...@@ -1004,7 +1004,7 @@ class Modeller(object):
del context del context
return actualVariants return actualVariants
def addExtraParticles(self, forcefield): def addExtraParticles(self, forcefield, ignoreExternalBonds=False):
"""Add missing extra particles to the model that are required by a force """Add missing extra particles to the model that are required by a force
field. field.
...@@ -1090,7 +1090,7 @@ class Modeller(object): ...@@ -1090,7 +1090,7 @@ class Modeller(object):
signature = _createResidueSignature([atom.element for atom in residue.atoms()]) signature = _createResidueSignature([atom.element for atom in residue.atoms()])
if signature in forcefield._templateSignatures: if signature in forcefield._templateSignatures:
for t in forcefield._templateSignatures[signature]: for t in forcefield._templateSignatures[signature]:
if compiled.matchResidueToTemplate(residue, t, bondedToAtom, False) is not None: if compiled.matchResidueToTemplate(residue, t, bondedToAtom, ignoreExternalBonds) is not None:
matchFound = True matchFound = True
if matchFound: if matchFound:
# Just copy the residue over. # Just copy the residue over.
...@@ -1109,7 +1109,7 @@ class Modeller(object): ...@@ -1109,7 +1109,7 @@ class Modeller(object):
if signature in forcefield._templateSignatures: if signature in forcefield._templateSignatures:
for t in forcefield._templateSignatures[signature]: for t in forcefield._templateSignatures[signature]:
if t in templatesNoEP: if t in templatesNoEP:
matches = compiled.matchResidueToTemplate(residueNoEP, templatesNoEP[t], bondedToAtomNoEP, False) matches = compiled.matchResidueToTemplate(residueNoEP, templatesNoEP[t], bondedToAtomNoEP, ignoreExternalBonds)
if matches is not None: if matches is not None:
template = t; template = t;
# Record the corresponding atoms. # Record the corresponding atoms.
......
Markdown is supported
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment