Difference between revisions of "Get systematic absences from .sca file"

From CCP4 wiki
Jump to navigationJump to search
(New page: Scalepack outputs only systematic absences, and not the rest of reflections. The best way to identify screw axis is probably to look at the whole set of reflections, so that you can compa...)
 
 
(2 intermediate revisions by the same user not shown)
Line 4: Line 4:
  
 
1.  Run scalepack in the spacegroup with no screw axis (e.g. P2, P222, I422, etc).
 
1.  Run scalepack in the spacegroup with no screw axis (e.g. P2, P222, I422, etc).
 +
 
2.  Use the following python script to extract specific reflection types.
 
2.  Use the following python script to extract specific reflection types.
  
 
The syntax is  
 
The syntax is  
  
./gethkl <file-name> 'absence code'
+
  ./gethkl <file-name> 'absence code'
  
 
Absence code should be in quotes but has a human(e) format, such as (0,0,l), (0,2k,0), etc.  Other examples are (h,0,0), (h,1,2), (h,2k,4) and so on.  
 
Absence code should be in quotes but has a human(e) format, such as (0,0,l), (0,2k,0), etc.  Other examples are (h,0,0), (h,1,2), (h,2k,4) and so on.  
Line 15: Line 16:
  
 
Click here for the code [http://docs.google.com/Doc?id=d64gxmg_0gvpgt9zq]
 
Click here for the code [http://docs.google.com/Doc?id=d64gxmg_0gvpgt9zq]
 +
 +
Back to [[Useful scripts (aka smart piece of code)]]

Latest revision as of 15:05, 16 December 2010

Scalepack outputs only systematic absences, and not the rest of reflections. The best way to identify screw axis is probably to look at the whole set of reflections, so that you can compare systematic absences with present reflections. Another relevant situation is when you need to identify which axis are screw in, say, P2221. You have to run scalepack at least twice to do it.

So, the general rule is:

1. Run scalepack in the spacegroup with no screw axis (e.g. P2, P222, I422, etc).

2. Use the following python script to extract specific reflection types.

The syntax is

 ./gethkl <file-name> 'absence code'

Absence code should be in quotes but has a human(e) format, such as (0,0,l), (0,2k,0), etc. Other examples are (h,0,0), (h,1,2), (h,2k,4) and so on.

Enjoy. (Tested on Ubuntu Gutsy, hence python 2.5).

Click here for the code [1]

Back to Useful scripts (aka smart piece of code)