2,652
edits
INTEGRATE: Difference between revisions
m
no edit summary
mNo edit summary |
|||
| Line 26: | Line 26: | ||
XDS does not adjust the integration boxes such as to center them individually on the observed reflections: it only tries to minimize the deviations between observed and calculated spot coordinates by adjusting about a dozen diffraction parameters (those given by REFINE(INTEGRATE)) for the reflections in a certain range of frames (DELPHI). | XDS does not adjust the integration boxes such as to center them individually on the observed reflections: it only tries to minimize the deviations between observed and calculated spot coordinates by adjusting about a dozen diffraction parameters (those given by REFINE(INTEGRATE)) for the reflections in a certain range of frames (DELPHI). | ||
== How does INTEGRATE treat overlaps? == | |||
First of all, there is no conceptual difference in XDS between overlap on a frame (due to too close detector, or smeared spots), and overlap by rotation (due to too large delta-phi, or high mosaicity). | |||
The integration algorithm is: | |||
# each pixel of the detector is assigned to its nearest reflection in reciprocal space ("pixel-labelling", see [http://dx.doi.org/10.1107/S0021889888007903]) | |||
# some of these pixels will mostly allow the background estimation, others will mostly contribute to the integration area (but as they are transformed into a local coordinate system [http://dx.doi.org/10.1107/S0021889888007903] there is not a 1:1 relationship). At this step, pixels which should be background but are higher than expected (due to overlap) are rejected. | |||
# for each reflection, the background is estimated, and the 3D profile is assembled from the pixels contributing to it | |||
# a comparison is made: for a reflection, is the percentage of its observed profile assembled in 3. larger than some constant (called "MINPK" in XDS.INP)? If the answer is no, this reflection will be discarded (too much "overlap"). Otherwise, the observed intensity (from the incomplete profile) is scaled up, using the inverse of the observed fraction (this relies on the accuracy of the average profile). | |||
Among other things, this means that: | |||
# the program does _not_ look around each reflection to detect an overlap situation, it just tries to gather the pixels for each reflection | |||
# as a user, when your crystal-detector distance was chosen too low, or the reflections are very broad, or if the crystal has a high mosaicity (all of which result in many overlaps), you may try reducing MINPK down to 70, 65, 60, 55 or even 50. This will result in more completeness, but you should monitor the quality of the resulting data. Conversely, if you raise MINPK over its default of 75 you will discard more reflections, but the resulting dataset will be a bit cleaner. | |||