ó (™xTc@ sÍddlmZddlZejjdƒZddlTddlZ ddl Z ddl Z e j Z xºddgD];Ze jjdeƒZe jeƒrree jd#sRi( tmmdbtManagertselect_atoms_with_same_resnametSetFlagtreducetoperatortor_tresidue_groupst PutPDBStringtformat_atom_recordt RuntimeErrortGetErrorDescription(troottflagstremove_duplicatestmanagertselectorRtatomtrc((sE/public2/ftp/pub/ccp4/6.5/ccp4-6.5/lib/py2/ccp4io_adaptbx/__init__.pytto_mmdbs    cC sri}x6|jƒD](}|j|jƒjgƒj|ƒqWtg|jƒD]}t|ƒ|f^qOƒdS(Ni(Rt setdefaulttparenttresnametappendtmaxtvaluestlen(Rtatoms_intat_((sE/public2/ftp/pub/ccp4/6.5/ccp4-6.5/lib/py2/ccp4io_adaptbx/__init__.pyR 2s&s"secondary structure does not matchsstructures are too remotes$can't make graph for first structuresempty graph for first structures%can't make graph for second structures empty graph for second structurecC stjj|dƒS(Nsundocumented return code(tssmtERROR_DESCRIPTION_FORtget(R((sE/public2/ftp/pub/ccp4/6.5/ccp4-6.5/lib/py2/ccp4io_adaptbx/__init__.pyREstSecondaryStructureMatchingcB sDeZdZejjejjd„Zd„Z d„Z d„Z RS(s-SSM matching with iotbx.pdb.hierarchy objectsc C sØg|_g|_g|_g|_xÂ||gD]´}td|ƒ}|jƒ}|jd|dtjj dddtj j ƒ|j d|ƒdkrµt dt|jƒd ‚n|jj|ƒ|jj|ƒ|jj|ƒq1Wt|jƒd kst‚t|jƒd kst‚t|jƒd ks:t‚tjƒ|_|jjd |jdd |jd d |d|d|jdd|jd ƒ}|tjjkrÂt tjd|ƒ‚n|jjƒ|_dS(NRtselHndtselTypetcidRtselKeyis%Empty atom selection for structure %siitmanager1tmanager2t precisiont connectivitytselHnd1tselHnd2R(tchainstmanagersthandlestqvaluesRt NewSelectiontSelectR tSELECTION_TYPEtATOMt SELECTION_KEYtNEWt GetSelLengthRR#R tAssertionErrorR'tSSMAligntAlignt RETURN_CODEtOkRt GetQvalues( tselftmovingt referenceR1R2tchainRthandleR((sE/public2/ftp/pub/ccp4/6.5/ccp4-6.5/lib/py2/ccp4io_adaptbx/__init__.pyt__init__Ls@            cC s|jS(N(R8(RF((sE/public2/ftp/pub/ccp4/6.5/ccp4-6.5/lib/py2/ccp4io_adaptbx/__init__.pyRE€scC s­|t|jƒkrdGHdS|jjd|jdd|jddtjjdtjjd|j dd |j dd |ƒ}|tj j kr©t tj d |ƒ‚ndS( NsNot that many matches availableR/iR0iR1R2R3R4t nselectedR(R#R8R'tAlignSelectedMatchR6t PRECISIONtNormalt CONNECTIVITYtFlexibleR7RCRDRR(RFRLR((sE/public2/ftp/pub/ccp4/6.5/ccp4-6.5/lib/py2/ccp4io_adaptbx/__init__.pyRM„s        cC s |jjS(N(R'tt_matrix(RF((sE/public2/ftp/pub/ccp4/6.5/ccp4-6.5/lib/py2/ccp4io_adaptbx/__init__.pyt get_matrix—s( t__name__t __module__t__doc__R'RNRORPRQRKRERMRS(((sE/public2/ftp/pub/ccp4/6.5/ccp4-6.5/lib/py2/ccp4io_adaptbx/__init__.pyR*Is  /  t SSMAlignmentcB s5eZdZd„Zd„Zd„ZeeƒZRS(sSSM alignment from SSM matchc C sßtjƒ}|jd|jdd|jdd|jƒ|d|jdƒ}|d|jdƒ}g|_g|_xa|jƒD]S\\}}}d„} |jj| ||ƒ| ||ƒfƒ|jj|ƒq„WdS(NR/iR0it ssm_alignRIcS sR|dkrdS|j|j|jf}||ksJtdt|ƒƒ‚||S(Ns Id %s missing(tNonetchain_idtresseqtinscodeR@tstr(trgitindexert identifier((sE/public2/ftp/pub/ccp4/6.5/ccp4-6.5/lib/py2/ccp4io_adaptbx/__init__.pyR)­s  "( R't XAlignTexttXAlignR6R5tpairststatst get_blocksR ( RFtmatchR_taligntindexer1tindexer2tftsR%R)((sE/public2/ftp/pub/ccp4/6.5/ccp4-6.5/lib/py2/ccp4io_adaptbx/__init__.pyRKŸs        (c C sFd}d}d}d}d}d}d}xÞ|jD]Ó}t|ƒdkr4|d7}|ddkrs|d7}n|ddkr|d7}n|ddkr­|d7}n|ddkrÊ|d7}n|ddkrç|d7}n|ddkr|d7}qq4q4W|d } || || || || || || |fS( Nigiigð?iiigY@(RdR#( RFtgaplessalignlents5ts4ts3ts2ts1ts0tetglen((sE/public2/ftp/pub/ccp4/6.5/ccp4-6.5/lib/py2/ccp4io_adaptbx/__init__.pytGetSSMSequenceIdentityºs0       cC sd„}|d|d|ƒS(NcS s>tg|jƒD]'}|j|jƒ|jf|f^qƒS(N(tdictRtidt resseq_as_intticode(RIR((sE/public2/ftp/pub/ccp4/6.5/ccp4-6.5/lib/py2/ccp4io_adaptbx/__init__.pyR×sRfR_((tclsRfR_((sE/public2/ftp/pub/ccp4/6.5/ccp4-6.5/lib/py2/ccp4io_adaptbx/__init__.pyRÕs (RTRURVRKRuRt classmethod(((sE/public2/ftp/pub/ccp4/6.5/ccp4-6.5/lib/py2/ccp4io_adaptbx/__init__.pyRWœs    ('t __future__Rt boost.pythontboosttpythont import_exttextRtlibtbx.load_envtlibtbxRtostpathtopR&tenvt under_disttisfiletenvirontgetenvtjointwarningstwarntTrueRR R'RCtNoHitstNoSuperpositiontNoGrapht NoVerticestNoGraph2t NoVertices2R(RtobjectR*RW(((sE/public2/ftp/pub/ccp4/6.5/ccp4-6.5/lib/py2/ccp4io_adaptbx/__init__.pyts:        3%         S