ó ű™xTc@ s{ddlmZddlmZdefd„ƒYZedkrwddlZx&ejdD]Zej eƒGHq\WndS(i˙˙˙˙(tdivision(tFormatPYtFormatPYmultitilecB sYeZed„ƒZd„Zd„Zd„Zd„Zdd„Z d„Z d„Z RS( cC s¤y1ddl}tj|dƒ}|j|ƒ}Wntk rEtSXddg}x|D]}||krYtSqYWddg}x|D]}||kr†tSq†WtS(Ni˙˙˙˙trbtTILESt METROLOGYtSIZE1tSIZE2(tcPickleRt open_filetloadtIOErrortFalsetTrue(t image_filetpickletstreamtdatatwanted_header_itemst header_itemtunwanted_header_items((sL/public2/ftp/pub/ccp4/6.5/ccp4-6.5/lib/py2/dxtbx/format/FormatPYmultitile.pyt understand s        cC s)|j|ƒst‚tj||ƒdS(s3Initialise the image structure from the given file.N(RtAssertionErrorRt__init__(tselfR((sL/public2/ftp/pub/ccp4/6.5/ccp4-6.5/lib/py2/dxtbx/format/FormatPYmultitile.pyR!scC sQddlm}||jƒ|_|jjƒ|jj|_|jj|_dS(Ni˙˙˙˙(t CSPadDetector(t"xfel.cftbx.detector.cspad_detectorRt _image_filet detectorbaset readHeadert_metrology_paramst_tiles(RR((sL/public2/ftp/pub/ccp4/6.5/ccp4-6.5/lib/py2/dxtbx/format/FormatPYmultitile.pyt_start)s   cC s |jjƒS(N(t_goniometer_factoryt single_axis(R((sL/public2/ftp/pub/ccp4/6.5/ccp4-6.5/lib/py2/dxtbx/format/FormatPYmultitile.pyt _goniometer=sc C sddlm}ddlm}ddlm}ddlm}m}|j j }|||j ƒj ƒ||j ƒ|dd|j j dfƒƒd}g|_|ƒ}xe|jD]Z} |||| j ƒj ƒ|| j ƒƒd} x"| jD]} | ||| j ƒj ƒ|| j ƒƒd} xß| jD]Ô} | ||| j ƒj ƒ|| j ƒƒd}|| j| jƒd}|||d ƒ}|||d| jddfƒ|}|||| jdddfƒ|}|j| j| j| jf}|jƒ}|jd ƒ|jd |ƒ|jg|jdd !D]}|d ^qFg|jdd !D]}|d ^qjg|jdd !D]}|d ^qŽƒ|jg| jD]}|d ^qľƒ|j| jƒ|jd| jfƒ|jj|j |ƒq5WqöWqˇW|S(säThe _detector() function returns a model for a CSPAD detector as used at LCLS's CXI and XPP endstations. It converts the metrology information in the pure Python object extracted from the image pickle to DXTBX-style transformation vectors. Only ASIC:s are considered, since DXTBX metrology is not concerned with hierarchies. Merged from xfel.cftbx.detector.cspad_detector.readHeader() and xfel.cftbx.detector.metrology.metrology_as_dxtbx_vectors(). i˙˙˙˙(tSimplePxMmStrategy(tHierarchicalDetector(tcol(t _transformtget_projection_matrixigüŠńŇMbP?itPADs %d:%d:%d:%dig@@(iii(!t dxtbx.modelR$tdxtbx.model.detectorR%t scitbx.matrixR&txfel.cftbx.detector.metrologyR'R(Rtdetectort orientationt normalizet translationtdistancet _raw_datatpaneltsensortasict pixel_sizet dimensiontserialt add_paneltset_typetset_nametset_local_frametelemstset_pixel_sizetset_image_sizetset_trusted_ranget saturationtappendR(RR$R%R&R'R(tdtTb_dR.tptTb_ptstTb_statTb_atPbtorigintfasttslowtkeyR4tt((sL/public2/ftp/pub/ccp4/6.5/ccp4-6.5/lib/py2/dxtbx/format/FormatPYmultitile.pyt _detectorAsP   (  ((  $$('#cC s"|dk r|j|S|jdS(s…If @p index is @c None, return the flex array for the first ASIC. Otherwise return the array for the ASIC at index @p index. iN(tNoneR3(Rtindex((sL/public2/ftp/pub/ccp4/6.5/ccp4-6.5/lib/py2/dxtbx/format/FormatPYmultitile.pyt get_raw_data”s  cC s|jj|jjƒS(s#Return a simple model for the beam.(t _beam_factorytsimpleRt wavelength(R((sL/public2/ftp/pub/ccp4/6.5/ccp4-6.5/lib/py2/dxtbx/format/FormatPYmultitile.pyt_beamžsc C sœddlm}ddlm}|jjd d}|jjd}|||dƒƒt|ƒ}|jjddd d |jj gd ddi|d 6ƒS(s+Return the scan information for this image.i˙˙˙˙(ttimegm(tstrptimeitUTCis%Y-%m-%dT%H:%M%Zt image_rangeitexposure_timesgV瞯Ňs Ž