Xscale: Difference between revisions

621 bytes added ,  20 December 2007
m
no edit summary
No edit summary
mNo edit summary
Line 1: Line 1:
[http://www.mpimf-heidelberg.mpg.de/~kabsch/xds/html_doc/xscale_parameters.html XSCALE] ist the scaling program of the XDS suite.
[http://www.mpimf-heidelberg.mpg.de/~kabsch/xds/html_doc/xscale_parameters.html XSCALE] ist the scaling program of the XDS suite. At the XDS website, there is a short a long commented example of [http://www.mpimf-heidelberg.mpg.de/~kabsch/xds/html_doc/INPUT_templates/XSCALE.INP XSCALE.INP]


A minimal input file to combine two datasets into one file is:
A minimal input file to combine two datasets into one file is:
Line 6: Line 6:
  INPUT_FILE= ../fae-ip/xds_2/XDS_ASCII.HKL
  INPUT_FILE= ../fae-ip/xds_2/XDS_ASCII.HKL


The latest versions do not require
SPACE_GROUP_NUMBER=
UNIT_CELL_PARAMETERS=
in XSCALE.INP because they pick these parameters up from the header of the first input reflection file.


To enable radiation damage correction ([http://dx.doi.org/10.1107/S0907444903006516 K. Diederichs, S. McSweeney and R. B. G. Ravelli (2003) Zero-dose extrapolation as part of macromolecular synchrotron data reduction. Acta Cryst. (2003). D59, 903-909]) it suffices to use the CRYSTAL_NAME keyword. '''The CRYSTAL_NAME parameters of different datasets do not have to be different'''. If they are different, this results in more degrees of freedom (namely, the slopes of the reflection intensity as a function of dose) for the program to fit the observed changes of intensities which are induced by radiation damage. However, if the datasets are based on the same crystal, or the datasets are based on crystals from the same drop, it is reasonable to assume that the slopes are the same.  
== Radiation damage correction ==
 
To "switch on" radiation damage correction ([http://dx.doi.org/10.1107/S0907444903006516 K. Diederichs, S. McSweeney and R. B. G. Ravelli (2003) Zero-dose extrapolation as part of macromolecular synchrotron data reduction. ''Acta Cryst.'' '''D59''', 903-909]) it suffices to use the CRYSTAL_NAME keyword. '''The CRYSTAL_NAME parameters of different datasets do not have to be different'''. If they are different, this results in more degrees of freedom (namely, the slopes of the reflection intensity as a function of dose) for the program to fit the observed changes of intensities which are induced by radiation damage. However, if the datasets are based on the same crystal, or the datasets are based on crystals from the same drop, it is reasonable to assume that the slopes are the same.  
Example:
Example:
  OUTPUT_FILE=fae-merge.ahkl  
  OUTPUT_FILE=fae-merge.ahkl  
INPUT_FILE= ../fae-ip/xds_1/XDS_ASCII.HKL  !
  INPUT_FILE= ../fae-ip/xds_1/XDS_ASCII.HKL  !
CRYSTAL_NAME=ip
  CRYSTAL_NAME=ip
INPUT_FILE= ../fae-ip/xds_2/XDS_ASCII.HKL  ! same crystal, but translated along z
  INPUT_FILE= ../fae-ip/xds_2/XDS_ASCII.HKL  ! same crystal, but translated along z
  CRYSTAL_NAME=ip
  CRYSTAL_NAME=ip
This is the recommended way as it avoids overfitting.
 
If, however, the crystals represent different heavy atom soaks, it is advisable to give a different CRYSTAL_NAME to each dataset. Example:
  OUTPUT_FILE=hg.ahkl
  INPUT_FILE= ../xds-hg/XDS_ASCII.HKL  ! a mercury soak
  CRYSTAL_NAME=Hg
OUTPUT_FILE=pt.ahkl
  INPUT_FILE= ../xds-pt/XDS_ASCII.HKL  ! a platinum soak
  CRYSTAL_NAME=Pt


== A hint for long-time XSCALE users ==


''Much more to come''
The latest versions do not require
SPACE_GROUP_NUMBER=
UNIT_CELL_PARAMETERS=
in XSCALE.INP because they pick these parameters up from the header of the first input reflection file.
2,684

edits