This is a REFMAC bug; now fixed I believe.. at least in version 82 and
higher..
You can download it from Garibs web site..
Eleanor
Jan Abendroth wrote:
> Hi all,
> I am trying to follow good practices and keep my set of free reflections
> between data sets, eg. in this case between an in-house low resolution and a
> synchrotron high resolution data set. High resolution data from hkl2000 were
> imported through the ccp4i task, keeping the low resolution FreeRs. This mtz
> file contains both unit cells, see below. The refined data set contains only
> one unit cell description, unfortunately the one from the FreeR
> (refmac5.5.72 and refmac5.5.70). As the two unit cells are sufficiently
> different, coot displays the model towards the edge of the density, real
> space refinement pulls the model back in the middle, refmac then starts with
> really high Rs and pulls the model back "out".
>
> When using rather ancient refmac5.2.0019, the mtz file has the correct unit
> cell description.
>
> Btw. both refinements look about the same. The only difference is a rather
> annoying shift of the electron density that is displayed in coot based on
> different unit cell in the mtz file.
>
> Is there a way to tell refmac which of the two unit cells to put in the
> output mtz file? Intuitively, it should be the one from which the amplitudes
> originate?
>
> Cheers
> Jan
>
>
> *mtz file after import*
> 1 myprotein
> high_reso
> synchrotron
> 79.0610 79.0610 311.7950 90.0000 90.0000 90.0000
> 1.00000
> 2 myprotein
> low_resol
> rotating_anode
> 78.5860 78.5860 311.1900 90.0000 90.0000 90.0000
> 1.54180
>
>
> *refmac5.5.0072*
> 2 myprotein
> low_reso
> rotating_anode
> 78.5860 78.5860 311.1900 90.0000 90.0000 90.0000
> 1.54180
>
> *refmac5.2.0019:*
> 1 myprotein
> high_reso
> synchrotron
> 79.0610 79.0610 311.7950 90.0000 90.0000 90.0000
> 1.00000
>
>
>
|