- cctbxbb - phenix-online.org

Atom selections for distance (embarrassingly simple)
by Morten Grøftehauge 18 Feb '15

18 Feb '15

Hi everyone, I'm trying to do some distance extraction of hetero atoms in pdbs. I've got for model in pdb_obj.hierarchy.models(): for chain in model.chains(): for rg in chain.residue_groups(): for ag in rg.atom_groups(): for atom in ag.atoms(): if atom.hetero == True: print atom.distance(atom) Which prints out 0.0 as expected. But I can't figure out how to, say, measure the distance to an atom in the previous residue (if there is a previous residue). I also noticed something weird in this tutorial: http://cctbx.sourceforge.net/sbgrid2008/tutorial.html They truncate a non-standard amino acid to the beta-carbon but this only works because there's an alternate conformation with a standard amino acid (Cys). So it wouldn't work with a seleno-met for instance (I know there's a convert-to-ala function in mmtbx - I said it was weird, not terrible). Anyway, I hope someone can explain the distances thing for me. I'm sure there are multiple ways of doing it (should I just do a selection of "hetatm" and iterate over them and try to get some distanes?). Cheers, Morten

2 1

Partitioning of Grid/Unit Cell
by Nicholas Pearce 15 Feb '15

15 Feb '15

Hi Everyone, I’m trying to partition a grid/space based on smallest distance to atomic sites - as in a voronoi diagram (http://en.wikipedia.org/wiki/Voronoi_diagram <http://en.wikipedia.org/wiki/Voronoi_diagram>). Is there any code in cctbx that could be used to accomplish this? Thanks, Nick Nicholas Pearce, DPhil Student, Structural Genomics Consortium, University of Oxford

3 2

iotbx xray.reconstructed_amplitude and changed SIGF
by Huw Jenkins 12 Feb '15

12 Feb '15

Hi, This is my first post to this list and I hope I can explain my confusion properly. I have been using iotbx to manipulate data stored in mtz files and found that values in columns that should not have been altered by my code are different in the output mtz file. This seems only to occur when the input mtz file contains DANO, SIGDANO columns. When I run the code below F and DANO in input.mtz and output.mtz are the same but SIGF and SIGDANO are different (shown for 5 reflections listed). Is this expected behaviour? Given that in my understanding that code should just read a mtz file, convert the data to a miller array then write the array back out to a new mtz I would expect the values in the input and output to be identical (within rounding error). Thanks, Huw from __future__ import division import sys from iotbx import mtz from iotbx import reflection_file_reader from iotbx import reflection_file_utils def run(args): mtz_file = reflection_file_reader.any_reflection_file("input.mtz") mtz_as_miller_array = mtz_file.as_miller_arrays() f_obs = mtz_as_miller_array[0] mtz_out = f_obs.as_mtz_dataset(column_root_label="F_iotbx") mtz_out.mtz_object().write("output.mtz") if (__name__ == "__main__"): run(sys.argv[1:]) some reflections from input.mtz are: H K L F_ccp4 SIGF_ccp4 DANO_ccp4 SIGDANO_ccp4 ISYM_ccp4 1 2 2 4.90461 2.31536 0 3.27441 0 1 2 3 186.172 1.86381 119.745 3.84862 0 1 2 4 289.31 3.06709 -5.42197 6.46531 0 1 2 5 378.928 2.7118 87.3381 5.56933 0 1 2 6 540.866 4.38513 46.4155 8.8053 0 The same reflections in output.mtz are: H K L F_iotbx SIGF_iotbx DANOF_iotbx SIGDANOF_iotbx ISYMF_iotbx 1 2 2 4.90461 2.00516 0 4.01032 0 1 2 3 186.172 1.8943 119.745 3.7886 0 1 2 4 289.31 3.15096 -5.42197 6.30192 0 1 2 5 378.928 2.74848 87.3381 5.49695 0 1 2 6 540.866 4.3939 46.4155 8.7878 0

2 2

Jenkins spam imminent (but not now)
by markus.gerstel＠diamond.ac.uk 03 Feb '15

03 Feb '15

Hello everyone, We are currently building a continuous integration setup for cctbx using Jenkins at Diamond. This means Jenkins will run all tests for every source code change. If any test (or the build itself) fails, it will send a notification email to this list. This could result in a number of mails. However they follow a pretty obvious format, so should be easy to filter (and we can fine-tune the spamminess if needed). Mails will start soonishly. -Markus -- This e-mail and any attachments may contain confidential, copyright and or privileged material, and are for the use of the intended addressee only. If you are not the intended addressee or an authorised recipient of the addressee please notify us of receipt by returning the e-mail and do not use, copy, retain, distribute or disclose the information in or attached to the e-mail. Any opinions expressed within this e-mail are those of the individual and not necessarily of Diamond Light Source Ltd. Diamond Light Source Ltd. cannot guarantee that this e-mail or any attachments are free from viruses and we cannot accept liability for any damage which you may sustain as a result of software viruses which may be transmitted in or with the message. Diamond Light Source Limited (company no. 4375679). Registered in England and Wales with its registered office at Diamond House, Harwell Science and Innovation Campus, Didcot, Oxfordshire, OX11 0DE, United Kingdom

1 0

Best place to get ccp4io, boost, etc?
by Luc Bourhis 29 Jan '15

29 Jan '15

Hi, I noticed that the nightly build get their dependencies from http://cci.lbl.gov/repositories/ccp4io.gz and the like. Is this the preferred way to make sure one is in synch with trunk then? Best wishes, Luc

2 1

directory of cctbx tutorials
by David Waterman 14 Jan '15

14 Jan '15

Hi all, I've just spent a little time with Google collecting together tutorials for CCTBX and related open-source projects. I'm sure this is not a complete list, so I thought I'd share in case anyone knows of some other sources to add to it. I'm collecting this because the DIALS project is involved in workshops and teaching. My main interest is in diffraction data processing, but all tutorials that use core modules and types from cctbx are of interest, for potential updating/cannibalisation. General/developer - Tour of the cctbx http://cctbx.sourceforge.net/current/tour.html - Making compiled extensions. cctbx_introduction_3 at https://sourceforge.net/p/dials/code/HEAD/tree/trunk/doc/developer/ Data processing - Various tutorials on free-electron laser data processing (cctbx.xfel) at http://viper.lbl.gov/cctbx.xfel/index.php/Tutorials - CCTBX for diffraction calculations. cctbx_introduction_{1,2} at https://sourceforge.net/p/dials/code/HEAD/tree/trunk/doc/developer/ - Powder diffraction unit cell refinement ( <http://cctbx.sourceforge.net/siena2005/>Siena 2005) http://cctbx.sourceforge.net/siena2005/unit_cell_refinement.html - Data processing with DIALS. 4 tutorials at http://dials.sourceforge.net/doc/index.html Coordinate handling - PDB truncate to Ala (SBGrid 2008) (iotbx.pdb) at http://cctbx.sourceforge.net/sbgrid2008/tutorial.html Phasing - Direct methods (Siena 2005) http://cctbx.sourceforge.net/siena2005/direct_methods_light.html Model refinement - Rigid body refinement (IUCr 2008) at http://cctbx.sourceforge.net/iucr2008/rigid_body_refinement_core.html -- David

2 1

Boost Import Error
by Garrett Thompson 12 Jan '15

12 Jan '15

Hello, I recently looking at the cctbxbb mailing archive and I came across your email response a few months ago : http://phenix-online.org/pipermail/cctbxbb/2014-June/000933.html I was wondering if you could elaborate a bit on how to use the wrappers cctbx.python? I consider myself to be quite a beginner in this area and can't seem to import the same module, import scitbx.array_family.flex *ImportError: No module named boost.optional* When I do I receive the same error about having boost.optional being an import dependency. If you are able to, please be as explicit as possible, I am pretty bad with common conventions and terminology, as I have never studied computer science or any related area formally. Any help would be much appreciated! Sincerely, Garrett Thompson

2 1

e density from atomic model?
by TJ Lane 06 Jan '15

06 Jan '15

Crystallography pros, Field-outsider asking for advice here. Is there a standard way to estimate the electron density on a grid from an atomic model (e.g. PDB coordinates of atoms)? One thing that occurred to me was to simply FT the Cromer-Mann atom-centered Gaussians commonly used to model atomic form factors. Perhaps there is a more standard way, though? And possible code to mooch...? Thanks much! TJ

2 1

building on Mac OSX
by Roger Martin 05 Jan '15

05 Jan '15

Hi, I'm attempting to build on Mac. Via jenkins. The scons finds /usr/bin/c++ but my builds are based on gcc located at /opt/gcc/current/bin. A which c++ returns /opt/gcc/current/bin/c++. I tried to tell it by making a site_init.py according to http://www.scons.org/doc/production/HTML/scons-user.html#idm12919064 --------------- cd build mkdir -p site_scons cat <<EOF >./site_scons/site_init.py env = DefaultEnvironment(tools = ['gcc', 'gnulink'], CC = '/opt/gcc/current/bin/gcc', CXX = '/opt/gcc/currentbin/g++') EOF -------------- Do you have a way to tell it which compiler version to find that is compatible with the source code level? Regards, Roger

2 3

building from cctbx SVN repository
by Roger Martin 11 Dec '14

11 Dec '14

Hi, http://cctbx.sourceforge.net/current/installation.html shows svn export svn://svn.code.sf.net/p/cctbx/code/trunk/libtbx/development/cctbx_svn_getti… ./cctbx_svn_getting_started.csh where cctbx_svn_getting_started.csh starts with a curl: curl http://cci.lbl.gov/cctbx_build/results/current/cctbx_bundle.selfx -o cctbx_bundle.selfx yet currently http://cci.lbl.gov/cctbx_build/results/current/cctbx_bundle.selfx doesn't exist. http://cci.lbl.gov/cctbx_build/results/current/ and has the Nov 13 bundle. Has the build changed or this url http://cci.lbl.gov/cctbx_build/results/current/cctbx_bundle.selfx really missing? Thank you, I'm wanting to be at svn HEAD for development in case I find something to add or patch. Roger

1 0