REFERENCE DATA SET

Sometimes one wants to use a REFERENCE_DATA_SET which was not itself written out by XDS or XSCALE.

For a file to be accepted as a REFERENCE_DATA_SET by XDS, it must minimally look like:

!FORMAT=XDS_ASCII
!END_OF_HEADER
   1   0   4  48445. 983.361   
   1   0   5 2124.59 283.825  
   1   0   6  17955. 493.356   
   1   0   7  19328. 510.718   
   1   0   8  87172. 3437.40   
...
...
!END_OF_DATA

The columns do not have a specified width, and are separated by a blank (Fortran's free format). Additional columns besides h k l I sigI will be ignored. Anomalous data in that file will be treated as required by the setting of FRIEDEL'S_LAW in XDS.INP.

Please note that the REFERENCE_DATA_SET will only be used if SPACE_GROUP_NUMBER > 0, and proper UNIT_CELL_CONSTANTS are specified in XDS.INP !

Such a file, which I would call "fake XDS_ASCII" can be obtained with a little effort. As an example, the following assumes existence of a MTZ file refine.mtz containing column F-model:

# square model amplitudes
sftools <<EOF
read refine.mtz
select col F-model = present
calc col I-model = col F-model col F-model *
write I-model.mtz col I-model
EOF
# dump to ASCII format:
mtzdump hklin I-model.mtz > temp.hkl <<EOF
nref 10000000
end
EOF
# prepare I-model.hkl:
echo \!FORMAT=XDS_ASCII > I-model.hkl
echo \!END_OF_HEADER   >> I-model.hkl
# pick reflection info from temp.hkl:
awk '/LIST OF REFLECTIONS/,/<B><FONT COLOR=/' temp.hkl | tail -n +4 | grep -v FONT | awk '{print $0,1}' >> I-model.hkl
echo \!END_OF_DATA   >> I-model.hkl
# I-model.hkl is now ready to be used as REFERENCE_DATA_SET
rm temp.hkl

For use with XSCALE, the file must look like e.g.:

!FORMAT=XDS_ASCII
!NUMBER_OF_ITEMS_IN_EACH_DATA_RECORD=   5
!ITEM_H=1
!ITEM_K=2
!ITEM_L=3
!ITEM_IOBS=4
!ITEM_SIGMA(IOBS)=5
!END_OF_HEADER
   1   0   4  48445. 983.361   
   1   0   5 2124.59 283.825  
   1   0   6  17955. 493.356   
   1   0   7  19328. 510.718   
   1   0   8  87172. 3437.40   
...
...
!END_OF_DATA