ó ÷™xTc@ sgddlmZddlZddlmZd„Zd„Zd„Zdd d„ƒYZd „Z dS( iÿÿÿÿ(tdivisionN(tmatrixcC sÚttt|dƒjƒjƒƒ}t|ƒdks<t‚t|dƒ}|dd!\}}|dd!}|d}|dd!}tt|dd !ƒ\}} |d d !\} } |d } |d d !\} }|d d !|d d!}}|dd!}t|dƒ}|dd!}|dd!|dd!|dd!}}}i|d6|d6|d6|d6|d6|d6|d6| d6| d6| d6| d6| d6|d 6|d!6|d"6|d#6|d$6|d%6|d&6|d'6|d(6}|S()stParse the XDS XPARM file, which contains a description of the detector and experimental geometry, to a dictionary.tri*iiiiii i iiiiiiii!i$i'tstarting_framet phi_startt phi_widthtaxist wavelengthtbeamtnxtnytpxtpytdistancetoxtoytxtytnormalt spacegrouptcelltatbtc(tmaptfloattopentreadtsplittlentAssertionErrortint(txds_xparm_filetdataRRRRRRR R R R R RRRRRRRRRRtresults((sJ/public2/ftp/pub/ccp4/6.5/ccp4-6.5/lib/py2/dxtbx/model/detector_helpers.pytread_xds_xparms0$      )##&cC s‘|\}}}|\}}}tj|j|j|jƒjƒ}|jƒddks`t‚tj|j|j|jƒjƒ} | jƒddks¢t‚|j|ƒtjrÙ|j |ƒ} |j|ƒ} nLtj |j|ƒtjƒdkr|j |ƒ} tj} n |} d} | j | ƒ} |} |j| |ƒ}| j | ƒ}|| }tj |jƒdƒdkst‚|S(s»Compute reference frame rotation to rotate from the original frame given by original = (x, y, z) to the to reference frame given by final = (_x, _y, _z). Returns M where M.x = _x etc.gð?gH¯¼šò×z>g( Rtsqrtelemst transposet determinantRtangletmathtpitcrosstfabst$axis_and_angle_as_r3_rotation_matrix(toriginaltfinalRRtzt_xt_yt_ztOtFt_ra_xt_a_xt_m_xt_ra_zt_a_zt_m_zt_m((sJ/public2/ftp/pub/ccp4/6.5/ccp4-6.5/lib/py2/dxtbx/model/detector_helpers.pytcompute_frame_rotation6s*&&%  %cC s$|jdƒ|jdƒ|jƒS(s:Given a cbf handle, get the value for the undefined pixel.tarray_intensitiestundefined_value(t find_categoryt find_columntget_doublevalue(t cbf_handle((sJ/public2/ftp/pub/ccp4/6.5/ccp4-6.5/lib/py2/dxtbx/model/detector_helpers.pytfind_undefined_valueas  tdetector_helper_sensorscB sDeZdZdZdZdZdZed„ƒZed„ƒZ RS(sÀA helper class which allows enumeration of detector sensor technologies which should help in identifying specific detectors when needed. These are currently limited to IMAGE_PLATE CCD PAD.t SENSOR_CCDt SENSOR_PADtSENSOR_IMAGE_PLATEtSENSOR_UNKNOWNcC s,|tjtjtjtjgkr(tStS(N(RERFRGRHRItTruetFalse(t sensor_type((sJ/public2/ftp/pub/ccp4/6.5/ccp4-6.5/lib/py2/dxtbx/model/detector_helpers.pyt check_sensorrs  cC stjtjtjgS(N(RERFRGRH(((sJ/public2/ftp/pub/ccp4/6.5/ccp4-6.5/lib/py2/dxtbx/model/detector_helpers.pytall{s( t__name__t __module__t__doc__RFRGRHRIt staticmethodRMRN(((sJ/public2/ftp/pub/ccp4/6.5/ccp4-6.5/lib/py2/dxtbx/model/detector_helpers.pyREhs c C s<ddlm}|j|jƒƒ\}}|jtt|ƒƒƒ}|j|ƒ|}x‡|D]}|j|jƒƒ} | |j|jƒƒ|d|j|j ƒƒ|d} |j d|jƒd|j ƒd| ƒq`W|j|jƒƒ\}}|j|ƒ|jtt|ƒƒƒj ƒdks8t ‚d S( s_ detector and beam are dxtbx objects, mosflm_beam_centre is a tuple of mm coordinates. iÿÿÿÿ(Riit fast_axist slow_axistorigingíµ ÷Æ°>N( tscitbxRtget_ray_intersectiontget_s0tcolttupletreversedt get_origint get_fast_axist get_slow_axistset_local_frametlengthR( tdetectorRtmosflm_beam_centreRtpanel_idtold_beam_centretnew_beam_centret origin_shifttpanelt old_origint new_origin((sJ/public2/ftp/pub/ccp4/6.5/ccp4-6.5/lib/py2/dxtbx/model/detector_helpers.pytset_mosflm_beam_centre‚s    (( t __future__RR)RVRR#R=RDRERj(((sJ/public2/ftp/pub/ccp4/6.5/ccp4-6.5/lib/py2/dxtbx/model/detector_helpers.pyts  ( +