ó
�™xTc           @   sÄ   d  d l  m Z d  d l m Z d  d l m Z d  d l m Z m Z d  d l	 m
 Z
 d  d l Z d  d l m Z d  d l Z d  d l Z e d „ Z e d	 k rÀ e d
 d e j k ƒ d GHn  d S(   iÿÿÿÿ(   t   division(   t   merging_statistics(   t   flex(   t   approx_equalt   Exception_expected(   t   SorryN(   t   StringIOc         C   s`  t  j j d ƒ s d GHd  St  j j d d d t j j ƒ } | d d d g } |  rr | j d	 ƒ d
 j | ƒ GHn  t	 ƒ  } t
 j | d | ƒ} |  r¤ | j ƒ  GHn  d | j ƒ  k s¼ t ‚ d | j ƒ  k sÔ t ‚ d | j ƒ  k sì t ‚ | j ƒ  } d | k s
t ‚ t t | d ƒ | j j ƒ s/t ‚ t t | d ƒ | j j ƒ sTt ‚ t t | d ƒ | j j ƒ syt ‚ t t j | d ƒ d d d d d d d d d d g
 ƒ s¹t ‚ d | k sËt ‚ d | k sÝt ‚ | j d  d! ƒ j ƒ  } d" | k st ‚ d# | k st ‚ t | d$  ƒ d% d& g } t	 ƒ  } t
 j | d | ƒ} |  re| j ƒ  GHn  d' | j ƒ  k s}t ‚ t | d$  ƒ d( d) d* g } t	 ƒ  } t
 j | d | ƒ} |  rÌ| j ƒ  GHn  d+ | j ƒ  k sät ‚ d, | j ƒ  k süt ‚ t | d$  ƒ d- d. g } y t
 j | d | ƒ} Wn t k
 rD} n Xt ‚ t | d$  ƒ d( d/ g } y t
 j | d | ƒ} Wn t k
 r“} n Xt ‚ t  j j d d0 d t j j ƒ } t  j j d d1 d t j j ƒ }	 | d2 |	 d3 g } t	 ƒ  } t
 j | d | ƒ} |  r!| j ƒ  GHn  x8 d4 j ƒ  D]* }
 |
 | j ƒ  k s.t | j ƒ  ƒ ‚ q.Wd  S(5   Nt   phenix_regressions+   phenix_regression not configured, skipping.t   relative_paths&   phenix_regression/wizards/p9_se_w2.scat   tests   space_group=I4s)   unit_cell=113.949,113.949,32.474,90,90,90s   loggraph=Trues
   debug=Truet    t   outs   R-merge: 0.073s   R-meas:  0.079sX     1.81   1.74  12528   2073    6.04  97.10    1449.2     5.5  0.252  0.275  0.110  0.968t   _reflns_shells   _reflns.pdbx_Rpim_I_obss   _reflns.phenix_cc_stars   _reflns.phenix_cc_1/2s!   _reflns_shell.number_measured_obsiy=  ip=  i4=  i<  i”:  i³9  iK6  ií4  i–3  ið0  s   _reflns_shell.phenix_cc_stars   _reflns_shell.phenix_cc_1/2t
   wavelengthgKÈ=›Uï?s3   REMARK 200  <I/SIGMA(I)> FOR SHELL         : 5.4942s3   REMARK 200  WAVELENGTH OR RANGE        (A) : 0.9792iÿÿÿÿs   high_resolution=2.5s   low_resolution=15s   Resolution: 14.96 - 2.50s   high_resolution=1.5s   low_resolution=100s   --extend_d_max_mins   Resolution: 100.00 - 1.50s_     1.55   1.50      0      0    0.00   0.00       0.0     0.0   None   None   None  0.000  0.000s   high_resolution=15s   low_resolution=2.5s   low_resolution=1.6s)   phenix_regression/harvesting/unmerged.scas'   phenix_regression/harvesting/refine.mtzs   symmetry_file="%s"s   --estimate_cutoffss3    resolution of all data          :   2.000
  based on CC(1/2) >= 0.33        :   2.000
  based on mean(I/sigma) >= 2.0   :   2.155
  based on R-merge < 0.5          :   2.372
  based on R-meas < 0.5           :   2.520
  based on completeness >= 90%    :   2.000
  based on completeness >= 50%    :   2.000(   t   libtbxt   envt
   has_modulet   find_in_repositoriest   ost   patht   isfilet   appendt   joinR   R   t   runt   getvaluet   AssertionErrort   as_cif_blockR   t   floatt   overallt   r_pimt   cc_start   cc_one_halfR   t   intt   as_remark_200t
   splitlinest   listR   R   (   t   debugt   hkl_filet   argsR   t   resultt	   cif_blockt
   remark_200t   args2t   st	   symm_filet   line(    (    sU   /public2/ftp/pub/ccp4/6.5/ccp4-6.5/lib/py2/iotbx/regression/tst_merging_statistics.pyt   exercise   s–    		%%%-	
			t   __main__R$   s   --debugt   OK(   t
   __future__R    t   iotbx.command_lineR   t   cctbx.array_familyR   t   libtbx.test_utilsR   R   t   libtbx.utilsR   t   libtbx.load_envR   t	   cStringIOR   R   t   syst   FalseR.   t   __name__t   argv(    (    (    sU   /public2/ftp/pub/ccp4/6.5/ccp4-6.5/lib/py2/iotbx/regression/tst_merging_statistics.pyt   <module>   s   ]