[chimera-dev] python for finding clashes
Eric Pettersen
pett at cgl.ucsf.edu
Tue Oct 28 11:17:45 PDT 2008
Hi JD,
You should use DetectClash.detectClash(). The only mandatory
argument is a list of atoms. So if your selection was in a variable
named 's', you might call the function like so:
from DetectClash import detectClash
clashes = detectClash(s.atoms())
The return value is a dictionary. The keys of the dictionary are
atoms that have clashes. The values are dictionaries. These
dictionaries have keys which are the atoms that clash with the
original key atom. The values are the amount of the clash. So to get
the number of clashes, you would use something like this:
sum = 0
for clashDict in clashes.values():
sum += len(clashDict)
numClashes = sum / 2
Now, there are a lot of optional keyword arguments to the
detectClash() function. They are pretty well documented in the doc
string of the function. So you should look through the doc string and
post any questions that you might still have...
--Eric
Eric Pettersen
UCSF Computer Graphics Lab
http://www.cgl.ucsf.edu
On Oct 28, 2008, at 7:15 AM, Jean Didier Pie Marechal wrote:
> Hi everyone,
>
> I am writing a script for which I need the number of clashes for a
> given selection for continuing.
>
> My first thought was to simply run in my script:
> runCommand("findclash selection ignoreIntraRes false")
> after having selected the region I was interested in.
> However, I realize that I come to the recurrent problem: having the
> information I am interested in printed in the reply log. I don't
> know how to have the info passed into an array or any kind of
> variable in the python environment. Could you tell me how to do so
> please?
>
> On the other side, I would be pleased if you could tell me what is
> the correct arguments I have to pass and syntax I have to use for
> the functions detectClash and/or cmdDetectClash (which I guess is
> the basis of the same findclash).
>
> All the best,
> JD
>
> Dr. Jean-Didier Maréchal
> Lecturer
> Computational Bioorganic and Bioionorganic Chemistry @ Transmet
> Unitat de Química Física
> Departament de Química
> Universitat Autònoma de Barcelona
> Edifici C.n.
> 08193 Cerdanyola (Barcelona)
> Tel: +34.935814936
> e-mail: JeanDidier.Marechal at uab.es
>
>
>
> _______________________________________________
> Chimera-dev mailing list
> Chimera-dev at cgl.ucsf.edu
> http://www.cgl.ucsf.edu/mailman/listinfo/chimera-dev
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://plato.cgl.ucsf.edu/pipermail/chimera-dev/attachments/20081028/bce05d00/attachment.html>
More information about the Chimera-dev
mailing list