ó ÷™xTc@ sYddlmZddlZddlmZmZddlmZddd„ƒYZdS(iÿÿÿÿ(tdivisionN(t GoniometertKappaGoniometer(t!cbf_gonio_to_effective_axis_fixedtgoniometer_factorycB seZdZd„Zed„ƒZed„ƒZed„ƒZed„ƒZed„ƒZ ed„ƒZ ed„ƒZ ed „ƒZ RS( sA factory class for goniometer objects, which will encapsulate some standard goniometer designs to make it a little easier to get started with all of this - for cases when we are not using a CBF. When we have a CBF just use that factory method and everything will be peachy.cC sdS(N((tself((sD/public2/ftp/pub/ccp4/6.5/ccp4-6.5/lib/py2/dxtbx/model/goniometer.pyt__init__scC s+tttt|ƒƒttt|ƒƒƒS(N(Rttupletmaptfloat(t rotation_axistfixed_rotation((sD/public2/ftp/pub/ccp4/6.5/ccp4-6.5/lib/py2/dxtbx/model/goniometer.pytmake_goniometerscC s=tt|ƒt|ƒt|ƒt|ƒt|ƒt|ƒƒS(N(RR tstr(talphatomegatkappatphit directiont scan_axis((sD/public2/ftp/pub/ccp4/6.5/ccp4-6.5/lib/py2/dxtbx/model/goniometer.pytmake_kappa_goniometer"s     c C sd}d}tj||ƒS(sUConstruct a single axis goniometer which is canonical in the CBF reference frame.ii(iii( iiiiiiiii(RR (taxistfixed((sD/public2/ftp/pub/ccp4/6.5/ccp4-6.5/lib/py2/dxtbx/model/goniometer.pyt single_axis,sc C sd}d}tj||ƒS(soConstruct a single axis goniometer which is canonical in the CBF reference frame, but reversed in rotation.iõÿÿÿii(iõÿÿÿii( iiiiiiiii(RR (RR((sD/public2/ftp/pub/ccp4/6.5/ccp4-6.5/lib/py2/dxtbx/model/goniometer.pytsingle_axis_reverse6sc C s+t|ƒdkst‚d}t||ƒS(s‰Return an goniometer instance for a known rotation axis, assuming that nothing is known about the fixed element of the rotation axis.iii( iiiiiiiii(tlentAssertionErrorR(RR((sD/public2/ftp/pub/ccp4/6.5/ccp4-6.5/lib/py2/dxtbx/model/goniometer.pyt known_axis@scC stj||||||ƒS(sUReturn a kappa goniometer where omega is the primary axis (i,e. aligned with X in the CBF coordinate frame) and has the kappa arm with angle alpha attached to it, aligned with -z, +y, +z or -y at omega = 0, that being the direction, which in turn has phi fixed to it which should initially be coincident with omega. We also need to know which axis is being used for the scan i.e. phi or omega. All angles should be given in degrees. This will work by first constructing the rotation axes and then composing them to the scan axis and fixed component of the rotation.(RR(RRRRRR((sD/public2/ftp/pub/ccp4/6.5/ccp4-6.5/lib/py2/dxtbx/model/goniometer.pyRKs cC s]tjƒ}|j|tjƒ|jƒ}t|ƒ\}}|j|ƒ~tj||ƒS(s2Initialize a goniometer model from an imgCIF file.( tpycbftcbf_handle_structt read_filet MSG_DIGESTtconstruct_goniometerRt__swig_destroy__RR (tcif_filet cbf_handlet cbf_gonioRR((sD/public2/ftp/pub/ccp4/6.5/ccp4-6.5/lib/py2/dxtbx/model/goniometer.pytimgCIFZs   cC s>|jƒ}t|ƒ\}}|j|ƒ~tj||ƒS(svInitialize a goniometer model from an imgCIF file handle, where it is assumed that the file has already been read.(R RR!RR (R#R$RR((sD/public2/ftp/pub/ccp4/6.5/ccp4-6.5/lib/py2/dxtbx/model/goniometer.pytimgCIF_Hjs   ( t__name__t __module__t__doc__Rt staticmethodR RRRRRR%R&(((sD/public2/ftp/pub/ccp4/6.5/ccp4-6.5/lib/py2/dxtbx/model/goniometer.pyRs     (( t __future__RRtdxtbx_model_extRRtgoniometer_helpersRR(((sD/public2/ftp/pub/ccp4/6.5/ccp4-6.5/lib/py2/dxtbx/model/goniometer.pyts