ó
4™xTc           @   sC  d  d l  m Z d  d l m Z d  d l Z d  d l Z d  d l m	 Z	 e	 j
 d d d d d	 d
 d d d d d d d d d d d d d d d d d d d d d g ƒ Z e	 j
 d  d! d" d# d$ d% d& d' d( d) d* d+ d, d- d. d/ d0 d1 d2 d3 d4 d5 d6 d7 d8 d9 d: g ƒ Z e	 j
 e j ƒ  d; ƒ Z d< „  Z e d= k r?e ƒ  n  d S(>   iÿÿÿÿ(   t   division(   t   xray_scatteringN(   t   flexg        gš™™™™™©?gš™™™™™¹?g333333Ã?gš™™™™™É?g      Ð?g333333Ó?gffffffÖ?gš™™™™™Ù?gÍÌÌÌÌÌÜ?g      à?g333333ã?gffffffæ?gš™™™™™é?gÍÌÌÌÌÌì?g      ð?g333333ó?gffffffö?gš™™™™™ù?gÍÌÌÌÌÌü?g       @g      @g      @g      @g      @g      @g      @g      $@gÌîÉÃB-#@g‹lçû©!@gŽuqà@g46<½Ò@gPüs×’@g_)ËÇ:@g'1¬Z@g’ËH¿}@gú~j¼t“@g;pÎˆÒ^@g×£p=
×þ?g¼t“Vû?gh³êsµù?gÊ2Ä±.n÷?g™*•Ô	ö?gŒÛh o�ó?gýöuàœñ?g°rh‘í|í?g€H¿}8é?gèÙ¬ú\må?g:#J{ƒ/Ü?gíž<,ÔšÒ?g~Œ¹k	ùÈ?g±áé•²Á?g o�Å�±?gà-� ø¡?g-Cëâ6
?c          C   sw   t  j d ƒ }  t j j j t t t |  j	 ƒ  ƒ } d Gt
 j | j ƒ  ƒ GHt j j j d d d |  j ƒ  d | ƒ d  S(   Ns   O2-s
   max error:t	   plots_dirt	   rez_plotst   labelt   gaussian_fit(   R   t   wk1995t   scitbxt   matht   gaussiant   fitt   table_2_stolt   table_2_o2minust   table_2_sigmast   fetchR   t   maxt   significant_relative_errorst   cctbxt   eltbxR   t   write_plotsR   (   t   wkR   (    (    sS   /public2/ftp/pub/ccp4/6.5/ccp4-6.5/lib/py2/cctbx/eltbx/development/rez_rez_grant.pyt   runE   s    t   __main__(   t
   __future__R    t   cctbx.eltbxR   t   scitbx.math.gaussian_fitR   t   cctbx.eltbx.gaussian_fitR   t   cctbx.array_familyR   t   doubleR   R   t   sizeR   R   t   __name__(    (    (    sS   /public2/ftp/pub/ccp4/6.5/ccp4-6.5/lib/py2/cctbx/eltbx/development/rez_rez_grant.pyt   <module>   s€   	