Re: [Open Babel] GAFF atom assignments

2011-12-23 Thread Eric Jang
Thanks Geoff, I took some time to really examine the SMARTS language today and it's not so bad - I panicked when I saw things like aliphatic/aromatic. Instead of writing a parser, I am going to try and use a bunch of if/else clauses. Would it be logically sound to proceed from bottom of the list

Re: [Open Babel] GAFF atom assignments

2011-12-22 Thread Geoff Hutchison
You don't have to use SMARTS -- you could program the atom types with a set of if/else clauses. Essentially all of MMFF94 atom typing is done this way in Open Babel. It's probably more tedious, though. > Sorry, but my question remains as to how can element is assigned its atom > type based sole

Re: [Open Babel] GAFF atom assignments

2011-12-19 Thread Eric Jang
Hi Geoff, Don't get me wrong, I think that Babel is the best. But I am unable to use Babel for various reasons (I am intending to implement a web-based docking experiment, that is one of them). I've been using Babel to verify my calculations. I'm so close to finishing, and I have everything nearly

Re: [Open Babel] GAFF atom assignments

2011-12-19 Thread Geoff Hutchison
On Dec 19, 2011, at 7:48 PM, Eric Jang wrote: > 1. Is the assignment of atom types somehow based on the SMARTS string of a > molecule? It would make sense that each atom has some SMARTS string property > that is compared to the gaff.prm file. SMARTS is essentially a programming language. So a

Re: [Open Babel] GAFF atom assignments

2011-12-19 Thread Eric Jang
Thanks for pointing out parsmart.cpp! Sorry, let me rephrase my question (or separate my many related questions): I have a very basic atomic structure (not part of OpenBabel) that I am implementing and I would like to assign GAFF atom types to it. All I know are the element names and the types of

Re: [Open Babel] GAFF atom assignments

2011-12-19 Thread Geoff Hutchison
> So I guess the question becomes: how do I obtain the SMARTS string for each > atom in GAFF? Again, all I know is which atoms are connected to each other > and the bond type (single, double). I'm not sure I understand your question. SMARTS is a pattern matching system. So the Open Babel code

[Open Babel] GAFF atom assignments

2011-12-19 Thread Eric Jang
Hello, I am trying to figure out how the forcefieldgaff.cpp file reads the gaff.prm file to assign atom types. Given a molecular model where I know the bonds connected between atoms (and I know their basic element name C,N,O,H, etc.), how can I deduce their atom type? (e.g. n2,n1,cp,cd, etc.) wit