ó
¼™xTc           @   st   d  d l  m Z d  d l m Z d  d l m Z d  d l m Z d  d l Z d „  Z	 e
 d k rp e	 ƒ  d GHn  d S(	   iÿÿÿÿ(   t   division(   t   null_out(   t   easy_pickle(   t   StringIONc          C   sÙ  d d l  m }  d d l m } |  d d d t ƒ \ } } t ƒ  } | j d | | g d | ƒ } | j d	 d
 d t ƒ  ƒ t	 j
 | ƒ } | j ƒ  } t	 j
 | ƒ } t ƒ  }	 | j d |	 ƒ |	 j ƒ  j ƒ  }
 d |
 k sâ t ‚ d |
 k sô t ‚ d |
 k st ‚ | j d | | d g d t ƒ  ƒ } t ƒ  } | j ƒ  j d | ƒ | j ƒ  j ƒ  }
 d |
 k smt ‚ t ƒ  } | j d | | d g d | ƒ } | j ƒ  j ƒ  }
 d |
 k s»t ‚ t j | ƒ t j | ƒ d  S(   Niÿÿÿÿ(   t   generate_calcium_inputs(   t   find_peaks_holest	   file_baset   tst_find_peaks_holest	   anonymizet   argst   outt	   file_names   tst_fph_peaks.pdbt   logs'     anomalous H2O (anomalous > 3):      1s'     anomalous non-water atoms:          0s'     mFo-DFc >  9:                       0s   filter_peaks_by_2fofc=1.0s'     anomalous > 3:                      0s   include_peaks_near_model=Trues'     mFo-DFc >  9:                       1(   t)   mmtbx.regression.make_fake_anomalous_dataR   t   mmtbx.command_lineR   t   TrueR   t   runt   save_pdb_fileR   R   t   dumpst   get_summaryt   showt   getvaluet
   splitlinest   AssertionErrort   ost   remove(   R   R   t   mtz_filet   pdb_fileR
   t   peaks_holest   pt   st   spt   out2t   linest   out3(    (    sS   /public2/ftp/pub/ccp4/6.5/ccp4-6.5/lib/py2/mmtbx/regression/tst_find_peaks_holes.pyt   exercise   sB    									t   __main__t   OK(   t
   __future__R    t   libtbx.utilsR   t   libtbxR   t	   cStringIOR   R   R#   t   __name__(    (    (    sS   /public2/ftp/pub/ccp4/6.5/ccp4-6.5/lib/py2/mmtbx/regression/tst_find_peaks_holes.pyt   <module>   s   	#