Jan wrote:

Bob Hanson wrote:

Is this a "feature" or a "bug"? Should "within" restrict itself to a given model, or is it interesting or important that one be able to select "within" a distance across multiple model sets?

it is an important *feature* that one is able to use within across multiple model sets.

right, that's my point. It seems unnatural to me as well. That's why I was surprised and initially thought it was a bug. Personally I can't imagine it, either, unless....

...I suppose some day we will be loading two molecules and maneuvering them independently (like a drug and an active site). In that case they will (probably) be in two different frames because they will have been loaded from two different files, and the interesting thing might be to highlight what atoms in one molecule are near atoms in another.

this case is called "molecule" in RasMol2.7. and I don't want to load them as two separate PDB-MODELs in Jmol "animation frames". They could be addressed for maneuvering independently but should be in the same MODEL allowing measurements between molecules.


Jan, I wonder if you are confusing "model" and "molecule" perhaps? Your last statement there -- that you do NOT want to load them as two separate models suggests that you do NOT think that "it is an important feature to be able to use within across multiple models in multiple-model sets. Am I right that you are interested in this specifically across MOLECULES within the same model?

model set -- a set of frames, one model per frame
molecule -- as set of covalently bonded atoms, perhaps one per frame, but usually (as in crystals) one or more per model (per frame)




So maybe someday there will be a flag that can be set -- only within a model, only between models -- still I find it hard to believe one would ever want the current "within or between models".

I have a feeling that the logic would get totally unrulely if we tried to restrict "within" to a particular model. Because this has to result in a linear set of bits.

This is easy
select within(4.2, atomno=3 AND */1) AND */1


Right. Can you give an example of where you would want those */1 out? In other words:

select within(4.2, atomno=3)

? Question is whether you see a use for cross-MODEL within. (I think I do, but my point is simply that it shouldn't the the default.)

By the way, in the prototype I have added "within molecule" (not model -- real, covalent set molecules) and will be expanding on this.

select within(molecule, atomno=35)
#all atoms of molecule containing atomno=35 (in all frames)

select within(4.2, within(molecule, atomno=35)) and not within(molecule, atomno=35)
#all atoms in nearby molecules

color atoms molecule
#color from red to blue in molecule 1,2,3,..... by molecule

As a bonus, this calculates the molecular formula of the molecule. But I don't think I have that displayable yet.
I need "moleculeinfo" I think.

You can see this at work at http://fusion.stolaf.edu/chemistry/jmol/xtalx
Clicking the radio buttons should do some interesting molecule-based selections. You can see what those buttons are doing by using http://fusion.stolaf.edu/chemistry/jmol/xtalx/?NOAPPLET

highlight "first molecule"

select *; color atoms cpk;select within(molecule,atomno=1);color blue

highlight "by molecule"

select *; color atoms molecule



Bob


-------------------------------------------------------
Using Tomcat but need to do more? Need to support web services, security?
Get stuff done quickly with pre-integrated technology to make your job easier
Download IBM WebSphere Application Server v.1.0.1 based on Apache Geronimo
http://sel.as-us.falkag.net/sel?cmd=lnk&kid=120709&bid=263057&dat=121642
_______________________________________________
Jmol-users mailing list
[email protected]
https://lists.sourceforge.net/lists/listinfo/jmol-users

Reply via email to