ó ,§xTc@sÏddlZddlmZddlZddlmZddd„ƒYZd„Zd„Ze dkrËej d d krœej j d ƒej d ƒnx,ej d D]Zd eGHeeƒGHqªWndS(iÿÿÿÿN(t OrderedDict(t put_errortMtzMetacBs#eZd„Zd„Zd„ZRS(cCs»t|dƒtkst‚||_|ri|dd!\|_|_|_|d\|_|_|_ n||_ ||_ ||_ ||_ ||ksŸt‚t|ƒ|_||_dS(Nii(ttypettupletAssertionErrortcelltatbtctalphatbetatgammatsymmetryt sg_numbertdmintdmaxRtcolumnstfilename(tselfRR RRRRR((s4/public2/ftp/pub/ccp4/6.5/ccp4-6.5/lib/py2/c4/mtz.pyt__init__s     cCs d|jS(Ns…cell: %(cell)s symmetry: "%(symmetry)s" (symmetry group no. %(sg_number)d) resolution range: %(dmin)s - %(dmax)s columns: %(columns)s(t__dict__(R((s4/public2/ftp/pub/ccp4/6.5/ccp4-6.5/lib/py2/c4/mtz.pyt__str__scCst||jkr6td||jfƒtjdƒn|j|}||krptd||j||fƒtStS(NsColumn '%s' not found in %sis/Column '%s' in %s has type '%s' (expected '%s')(RRRtsystexittFalsetTrue(Rtlabelt expected_typetcol_type((s4/public2/ftp/pub/ccp4/6.5/ccp4-6.5/lib/py2/c4/mtz.pytcheck_col_types  (t__name__t __module__RRR(((s4/public2/ftp/pub/ccp4/6.5/ccp4-6.5/lib/py2/c4/mtz.pyRs  cCstjdd|gdtjdtjƒ}|jddƒ\}}|jƒ}|rgtd|ƒ‚n|jƒ}x_t|ƒD]Q\}}|jdƒs¡q€n|jd ƒrõy(t d „||d j ƒDƒƒ}WqÑt k rñd}qÑXq€|jd ƒr?|j d ƒdj ƒ} t|j ƒdjdƒƒ} q€|jdƒr…||d} t| dd!ƒ} t| dd!ƒ} q€|jdƒr«||dj ƒ}q€|jdƒr€||dj ƒ}q€q€Wt||ƒ}t|d| d| d| d| d|d|ƒS( sFfor now using mtzdump, directly calling libccp4/mtzlib would be bettertmtzdumptHKLINtstdintstdouttinputsHEAD ENDsmtzdump of %s failed.s * s0 * Dataset ID, project/crystal/dataset names, cecss|]}t|ƒVqdS(N(tfloat(t.0tx((s4/public2/ftp/pub/ccp4/6.5/ccp4-6.5/lib/py2/c4/mtz.pys 9sis * Space group = t'iiÿÿÿÿt)s * Resolution Rangeiii(i)i5s * Column Labelss * Column TypesR RRRRRN(t subprocesstPopentPIPEt communicatetpollt RuntimeErrort splitlinest enumeratet startswithRtsplitt ValueErrortNonetstriptinttrstripR&tzipR(thklintpt stdoutdatat stderrdatatretcodetlinestntlineRR Rtres_linet lower_resolt upper_resolt column_namest column_typesR((s4/public2/ftp/pub/ccp4/6.5/ccp4-6.5/lib/py2/c4/mtz.pyt read_metadata+s<  (  "cCsÂddg}t|ƒ}|j|jkrItd|j|jfƒ‚n|j|jkrztd|j|jfƒ‚nx1|D])}||jkr|j|dƒ|SqWtd|ƒ‚dS(Nt FreeR_flagtFREEs&free-R-flags dmax: %g (should be < %g)s'free-R-flags dmin: %g (should be >= %g)tIsfree-R column not found in %s(RHRR5RRR(tfree_mtzt data_mtz_metatnamest rfree_metatname((s4/public2/ftp/pub/ccp4/6.5/ccp4-6.5/lib/py2/c4/mtz.pytcheck_freerflags_columnMs   t__main__iisNo filenames. isFile: %s((R+t collectionsRRtc4.utilsRRRHRQRtargvtstderrtwriteRtarg(((s4/public2/ftp/pub/ccp4/6.5/ccp4-6.5/lib/py2/c4/mtz.pyts  $ "