Hi Martin,
Here's a fourth option (and the technique worth noting :p ):
print len( set( [(i.chain,i.resi,i.resn) for i in
cmd.get_model(selection).atom] ) )
Cheers,
Tsjerk
On Tue, Nov 29, 2011 at 4:26 PM, Jason Vertrees
wrote:
> Hi Martin,
>
> You get three options:
>
> (1) You can count alph
Hi Martin,
You get three options:
(1) You can count alpha carbons:
fetch 1rx1, async=0
count_atoms n. CA
(2) But, a protein could be missing alpha carbons then this is more
complete (and the technique worth noting):
n=0
select qq, polymer
select pp, None
python
while cmd.count_atoms("qq"):