get_hkl: manipulate reflection data
Synopsis
get_hkl -i reflections.hkl -y pointgroup [options] ... -o output.hkl
get_hkl --help
Description
get_hkl performs various manipulations on reflection lists. Possible manipulations include: expanding to a point group of lower symmetry, 'twinning' to a point group of higher symmetry, adding noise, restricting the list to contain only reflections included in another 'template' list, and some less common miscellaneous manipulations.
The input filename should be specified with -i filename or --input=filename. The output filename should be specified with -o filename or --output=filename. Beyond these parameters, you can choose one of the manipulations described below.
The behaviour if multiple manipulations are requested is deliberately left undefined: it's much clearer to perform multiple manipulations in explicit separate steps.
- -p unitcell.cell
-p unitcell.pdb
--pdb=unitcell.pdb
Specify the name of the file containing the unit cell information, in PDB or CrystFEL format.
Choosing the output format
- --output-format=format
- The output file will be written in format, which can be mtz, mtz-bij or xds. Use mtz-bij to put Bijvoet pairs together in the same row, suitable for anomalous phasing. Otherwise, use mtz. If you omit this option, the output will be in the usual CrystFEL reflection list format.
Expanding reflections into a point group of lower symmetry
- -e pointgroup
--expand=pointgroup
- The reflections will be expanded out, according to the symmetry of the input reflections, to fill the asymmetric unit of pointgroup.
Reducing reflections into a point group of higher symmetry
- -w pointgroup
--twin=pointgroup
--no-need-all-parts
- Reflections in the list which are equivalent according to pointgroup will have their intensities summed. The output reflection list will contain the summed intensities in the asymmetric unit for pointgroup. Reflections for which any of the 'twin mates' are missing will not be written out, unless you use --no-need-all-parts
Adding noise
- --poisson
--noise
- Add either 10% flat random noise (with --noise) or generate Poisson noise (with --poisson). If Poisson noise is to be generated, the reflection intensities will be assumed to be measured in arbitrary units, and the conversion factor to photons must be given using --adu-per-photon=n, where n is the conversion factor.
- --adu-per-photon=n
- See --poisson.
Removing duplicated reflections
- --trim-centrics
- Reflections which are duplicated, according to the symmetry of the input reflection list, will be removed. This is sometimes useful when importing data from other programs.
Restricting reflections according to a template
- -t filename
--template=filename
- Only reflections in the input list which also appear in filename will be written to the output.
Multiplicity
- --multiplicity
- The intensities of the reflections will be multiplied by their symmetric multiplicities according to the point group of the input list.
Applying a resolution cutoff
- --cutoff-angstroms=n
--cutoff-angstroms=n1,n2,n3
--highres=n
--lowres=n
- In the first form, reflections with d (=lamba/2*sin(theta)) < n will be removed.
In the second form, anisotropic truncation will be performed with separate resolution limits n1, n2 and n3 along a*, b* and c* respectively. You must also specify -p or --pdb.
The option --lowres will remove reflections with d > n. --highres is a synonym for --cutoff-angstroms.
Reindexing reflections
- --reindex=op
- Reindex the reflections according to op. Example: --reindex=k,h,-l.
Overriding the input symmetry
- -y pg
--symmetry=pg
- Use pg for the symmetry of the input data. In almost all cases, you should omit this option which will cause the program to use the point group recorded in the headers of the input data.
Author
This page was written by Thomas White.
Reporting bugs
Report bugs to taw@physics.org
Copyright and Disclaimer
Copyright © 2012-2020 Deutsches Elektronen-Synchrotron DESY, a research centre of the Helmholtz Association.
Please read the AUTHORS file in the CrystFEL source code distribution for a full list of contributions and contributors.
CrystFEL is free software: you can redistribute it and/or modify it under the terms of the GNU General Public License as published by the Free Software Foundation, either version 3 of the License, or (at your option) any later version.
CrystFEL is distributed in the hope that it will be useful, but WITHOUT ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public License for more details.
You should have received a copy of the GNU General Public License along with CrystFEL. If not, see http://www.gnu.org/licenses/.
See also
If CrystFEL is installed on your computer, you can read this manual page offline using the command man get_hkl.