DTREK2MTZ (CCP4: Supported Program)
NAMEdtrek2mtz - converts d*trek scalemerge output into MTZ format.
SYNOPSISdtrek2mtz HKLIN foo.dat HKLOUT foo.mtz
DESCRIPTIONThis program converts merged data from d*trek into an MTZ format. After the conversion use TRUNCATE to convert the Is into Fs. Finally, CAD (or sortmtz) should be run to sort the data and put it into the correct asymmetric unit for CCP4.
Assumption: if I+ and I- are on separate lines, then I- appears first. This assumption is checked and a count of the failure of the assumption is output AND I+/I- are swapped. BUT if this happens, you should be very wary of the result mtz file since d*TREK designates I+ and I- properly.
These are converted to the weighted means Imean and SigImean
Imean = (I+/SIGI+ + I-/SIGI-) / (1/SIGI+ + 1/SIGI-) SigImean = 1 / (1/SIGI+ + 1/SIGI-)
Data HarvestingData harvesting (and some other functionality) relies on the existence of appropriate project and dataset names, and these should be assigned in DTREK2MTZ using the PNAME and DNAME keywords. These will be written to the the output MTZ file, and be inherited by subsequent programs.
Available keywords are:
ANOMALOUS, CELL, DNAME (DATA), END, MONITOR, PNAME (PROJECT), REJECT, RESOLUTION, SCALE, SYMMETRY, TITLEThe CELL parameters and SYMMETRY are compulsory keyword.
PNAME <projectname>[default: 'unknown']
Project Name used for Data Harvesting. This keyword is not compulsory, but is advisable to obtain a sensible name.
A dataset, as listed in the MTZ header, is specified by a project-name/dataset-name pair. The project-name specifies a particular structure solution project, while the dataset-name specifies a particular dataset contributing to the structure solution. An entry in the PNAME keyword should therefore be accompanied by a corresponding entry in the DNAME keyword.
DNAME <datasetname>[default: 'unknown<ddmmyy>']
Dataset Name used for Data Harvesting. This keyword is not compulsory, but is advisable to obtain a sensible name rather than 'unknown<ddmmyy>' (where <ddmmyy> is the date, with no spaces).
A dataset, as listed in the MTZ header, is specified by a project-name/dataset-name pair. The project-name specifies a particular structure solution project, while the dataset-name specifies a particular dataset contributing to the structure solution. An entry in the DNAME keyword should therefore be accompanied by a corresponding entry in the PNAME keyword.
CELL <a> <b> <c> [ <alpha> <beta> <gamma> ]The cell dimensions in Angstroms and degrees. As the d*trek reflection list file does not include the cell dimensions this keyword is compulsory. The angles default to 90 degrees.
ANOMALOUS YES | NOSpecify whether or not input file contains anomalous data. Defaults to YES.
MONITOR <nmon>Monitors the <nmon>th reflection and writes it to the output stream. The default is every 1000th reflection.
REJECT <h k l>If specified then the specific hkl will not be written to the MTZ file.
RESOLUTION <maxres> <minres>If specified then the all reflections outside the range will be excluded from the MTZ file. The limits can be given in either order and in Angstrom or SineSquared(theta)/LamdaSquared units. The default is that all reflections are output.
SCALE <scale>All intensities and their sigmas will be scaled by <scale>. The default is 1.0.
SYMMETRY <SG name|SG number>This keyword is compulsory and can be given as the space group name or number.
TITLE <string>This title will become the title/history in the MTZ file.
ENDTerminates the keyword input. Must be last keyword.
SEE ALSOscalepack2mtz, cad, truncate, Data Harvesting.
AUTHORJ. W. Pflugrath (Dec. 1999) Modified from SCALEPACK2MTZ