xTc3@ sFddlmZddlmZddlZddlmZddlmZddl m Z ddl m Z ddl m Z dd lmZdd lmZddlZddlZd Zd Zejd ddddddddddddddddddd d!d"d#d$d%d&d'd(d)d*d+d,d-d.d/d0d1d2d3d4d5d6d7d8d9d:d;d<d=d>g2Zejd?d@dAdBdCdDdEdFdGdHdIdJdKdLdMdNdOdPdQdRdSdTdUdVdWdXdYdZd[d\d]d^d_d`dadbdcdddedfdgdhdidjdkdldmdndodpg2ZdqejjfdrYZdsejjfdtYZdS(ui(tdivision(tabsolute_scalingN(t data_plots(tflex(t scale_curves(tsimplex(tchebyshev_polynome(tSorry(t table_utilsg4gW?gl|F?g5&lKۿgd?g P׿gUFd}?g t=տg N9t?g?gΔ${agnqc ?gwtFg 1ۧ?gW8gHC u?gf~5EgOLl\gS,y_?gӸ)uIg1gUb?gqg!s?g`j3ttg/V0^y?g }vg̽Xyp?gSkg}vq?g~)\ng/e?goVn0*agØXlkm?g0p;ggy4a?g~"Wgd&c?gA?gTh2'g<]ޭ?gzG|Ɲggio?g /Wlgr?gg>5!7?g :Sg>]J?g.gՙ?g gM^ ?g7G1j/gʼ?go^}g_$?gZوjYg.ɨ J?g•gmZ?gg-)?g^gD{|?gtl䗁g|j2jd|_?|j?sd t+j3|j2|jt+j3|j2|j}d}t+j,||gt+j,|d|dgt+j,|d|dgg|_6t7j8d|j6||_9ndS(Nit miller_arrayt auto_kerneltn_binsi-tn_termiigMbP?g?g@s6No reflections within required resolution range after s filtering.g@gg?g?g?ig0@(@tindicestall_eqtAssertionErrortis_xray_amplitude_arrayt f_as_f_sqt deep_copytobstcalctmindtmaxdtmtntlevelRtkernel_normalisationtTruetd_star_sq_arrayt bin_selectiontselecttlow_limthigh_limt mean_I_arrayt obs_d_star_sqtcalc_d_star_sqtmean_obst mean_calct var_I_arraytvar_obstvar_calcRtcurve_interpolatortinterpolt interpolatetobs_atobs_btcalc_atcalc_bRt mean_coefstsizetmean_ratio_enginet std_coefststd_ratio_engineRtdoubletxtlow_lim_for_scalingtcountRt as_doublet weight_arraytsumtmeantsigmast sart_simplexRt simplex_opttoptit get_solutiontabstscaletb_valuetmodify_weightstall_bad_z_scores(tselft miller_obst miller_calct min_d_star_sqt max_d_star_sqtn_pointsRt norma_obst norma_calcR#R$t sel_calc_obst sel_obs_calctselttmp_obs_d_star_sqttmp_calc_d_star_sqt selectiontstbtsol((sK/public2/ftp/pub/ccp4/6.5/ccp4-6.5/lib/py2/mmtbx/scaling/relative_wilson.pyt__init__s!                 0303     ,#   ,#cC stj|j|j|j|j}|jdkj}tj|j ||jj |}tj |jj d}|j |||j| |j |j|j|j|jg}td|d|jd|jS(Nit all_curvesRRG(RtexpR$RER&RDR%t iselectiontlogRR6R2t set_selectedtcurvet get_z_scorestsummaryRRG(RHti_scaledRRtratiotratio_tcurves((sK/public2/ftp/pub/ccp4/6.5/ccp4-6.5/lib/py2/mmtbx/scaling/relative_wilson.pyRa{s'( icC sP|j|j|j}tjttj||k}|j||_dS(N(R`RDRERR6tlisttboolR;(RHRtz_scoresRR((sK/public2/ftp/pub/ccp4/6.5/ccp4-6.5/lib/py2/mmtbx/scaling/relative_wilson.pyRFs$c C s&tj|j||j|}|jdk|dk@j}|jj||j|}|j|jj|}|jdst tj |}|j |jj|}|jj|}|j dst tj |||} tj|jjd} | j|| | S(Nii(RR[R$R&R%R\RR_tall_gtRR]tstdtall_neRCR6R2R^( RHRDRERbRRRcR=tvart d_star_sqtztz_((sK/public2/ftp/pub/ccp4/6.5/ccp4-6.5/lib/py2/mmtbx/scaling/relative_wilson.pyR`s!c C sid}t|d}|d}|dkr5d}n|dkrJd}ntj|j||j|}||j}|j|j}|tj|}tj|tj t j jkrt dn||}tj|j } | dkst| dkrItj|t j j| krIt ddtj|n||j }tj|}|S( Ng?iigi@gis*Result array exceeds floating-point limit.s1Weighted result array will exceed floating-point s limit: %e(RCRR[R$R&R%R_tmaxtmathtsqrttsyst float_infot OverflowErrorR;RR<( RHtvectortvRDRERbRcR_tresulttwmax((sK/public2/ftp/pub/ccp4/6.5/ccp4-6.5/lib/py2/mmtbx/scaling/relative_wilson.pyttargets,     ! $ + cC s|jj|}|S(N(R3tf(RHRmRx((sK/public2/ftp/pub/ccp4/6.5/ccp4-6.5/lib/py2/mmtbx/scaling/relative_wilson.pyR_scC s|jj|}|S(N(R5R{(RHRmRx((sK/public2/ftp/pub/ccp4/6.5/ccp4-6.5/lib/py2/mmtbx/scaling/relative_wilson.pyRjscC s|jjd|S(Ntout(Ratshow(RHR|((sK/public2/ftp/pub/ccp4/6.5/ccp4-6.5/lib/py2/mmtbx/scaling/relative_wilson.pyt show_summarys( t__name__t __module__RYRaRFR`RzR_RjR~(((sK/public2/ftp/pub/ccp4/6.5/ccp4-6.5/lib/py2/mmtbx/scaling/relative_wilson.pyR s `      RacB s2eZdedZdZdZdZRS(g@cC stjdddddddgddgd dgd ttd gdtdg|D]}t|^qU|_||_||_dS(NttitlesRelative Wilson plott column_labelssMax. resolutionslog(I_exp/I_obs)sReference curvesZ-scoret graph_namest graph_labelssHigh resolutiontt graph_columnsitx_is_inverse_d_mintdata(sHigh resolutionR(Rt table_dataRftrangeRttabletcutoffRG(RHRZRRGtarray((sK/public2/ftp/pub/ccp4/6.5/ccp4-6.5/lib/py2/mmtbx/scaling/relative_wilson.pyRYs     % cC s4|j\}}}}||jk}|jtS(N(tdata_as_flex_arraysRR9R(RHtsstrrtiitzztflagged((sK/public2/ftp/pub/ccp4/6.5/ccp4-6.5/lib/py2/mmtbx/scaling/relative_wilson.pyt n_outliersscC s&g|jjD]}tj|^q S(N(RRRR6(RHtcolumn((sK/public2/ftp/pub/ccp4/6.5/ccp4-6.5/lib/py2/mmtbx/scaling/relative_wilson.pyRscC s|j\}}}}||jk}|j|}|j|}|j|} |j|} |jd|jd|j|j|jr|jd|jn|jd|j|j g} t |dkrxet ||| | D]N\} } }}t j d| }| jd|d|d|d | gqWtjd d d d dgd| }|j|n |jddS(NsRelative Wilson plotsThe relative Wilson plot compares the mean intensity of the observed data with the mean intensity computed from the model, as a function of resolution. This curve is expected to fall off at low resolution if no contribution for bulk solvent is provided for the calculated intensities, because the presence of bulk solvent reduces the observed intensities at low resolution by reducing the contrast. At high resolution, the curve should be a straight line with a slope that reflects the difference in overall B-factor between the model and the data. Compared to the normal Wilson plot, the relative Wilson plot is more linear because the influence of favored distances between atoms, caused by bonding and secondary structure, is cancelled out. sAll resolution shells have Z-scores above %4.2f sigma. This is indicative of severe problems with the input data, including processing errors or ice rings. We recommend checking the logs for data processing and inspecting the raw images. sBAll relative wilson plot outliers above %4.2f sigma are reported. ig?s%8.2fs%9.3es%5.2ftcolumn_headerss d-spacingsObs. Log[ratio]sExpected Log[ratio]sZ-scoret table_rowss=The Relative wilson plot doesn't indicate any serious errors.(RRRtshow_sub_headert show_textt show_plotRRGtwarntnewlinetlentzipRqRrtappendRt simple_tablet show_tableR}(RHR|RRRRRtsel_sstsel_ztsel_rtsel_itrowsRVRntrtitsssR((sK/public2/ftp/pub/ccp4/6.5/ccp4-6.5/lib/py2/mmtbx/scaling/relative_wilson.pyt _show_impls6    (-    (RRtFalseRYRRR(((sK/public2/ftp/pub/ccp4/6.5/ccp4-6.5/lib/py2/mmtbx/scaling/relative_wilson.pyRas   (t __future__Rt mmtbx.scalingRtmmtbxtiotbxRtscitbx.array_familyRt scitbx.mathRtscitbxRRt libtbx.utilsRtlibtbxRRsRqR R!R6R1R4tscalingtxtriage_analysisR Ra(((sK/public2/ftp/pub/ccp4/6.5/ccp4-6.5/lib/py2/mmtbx/scaling/relative_wilson.pyts