#!/usr/bin/env python
#------------------------------
"""
Class :py:class:`SegGeometryStore` is a factory class/method
============================================================
Switches between different device-dependent segments/sensors
to access their pixel geometry using :py:class:`SegGeometry` interface.
Usage::
from PSCalib.SegGeometryStore import sgs
sg = sgs.Create('SENS2X1:V1', pbits=0377)
sg2= sgs.Create('EPIX100:V1', pbits=0377)
sg3= sgs.Create('PNCCD:V1', pbits=0377)
sg4= sgs.Create('ANDOR3D:V1', pbits=0377)
sg5= sgs.Create('JUNGFRAU:V1',pbits=0377)
sg.print_seg_info(pbits=0377)
size_arr = sg.size()
rows = sg.rows()
cols = sg.cols()
shape = sg.shape()
pix_size = sg.pixel_scale_size()
area = sg.pixel_area_array()
mask = sg.pixel_mask(mbits=0377)
sizeX = sg.pixel_size_array('X')
sizeX, sizeY, sizeZ = sg.pixel_size_array()
X = sg.pixel_coord_array('X')
X,Y,Z = sg.pixel_coord_array()
xmin = sg.pixel_coord_min('X')
ymax = sg.pixel_coord_max('Y')
xmin, ymin, zmin = sg.pixel_coord_min()
xmax, ymax, zmax = sg.pixel_coord_mas()
...
See:
* :py:class:`GeometryObject`,
* :py:class:`SegGeometry`,
* :py:class:`SegGeometryCspad2x1V1`,
* :py:class:`SegGeometryEpix100V1`,
* :py:class:`SegGeometryJungfrauV1`,
* :py:class:`SegGeometryMatrixV1`,
* :py:class:`SegGeometryStore`
For more detail see `Detector Geometry <https://confluence.slac.stanford.edu/display/PSDM/Detector+Geometry>`_.
This software was developed for the SIT project.
If you use all or part of it, please give an appropriate acknowledgment.
Created: 2013-03-08 by Mikhail Dubrovin
"""
#------------------------------
from PSCalib.SegGeometryCspad2x1V1 import cspad2x1_one
from PSCalib.SegGeometryEpix100V1 import epix2x2_one
from PSCalib.SegGeometryMatrixV1 import SegGeometryMatrixV1, segment_one, matrix_pars
from PSCalib.SegGeometryJungfrauV1 import jungfrau_one
#------------------------------
[docs]class SegGeometryStore() :
"""Factory class for SegGeometry-base objects of different detectors"""
#------------------------------
def __init__(sp) :
pass
#------------------------------
[docs] def Create(sp, segname='SENS2X1:V1', pbits=0 ) :
""" Factory method returns device dependent SINGLETON object with interface implementation
"""
if segname=='SENS2X1:V1' : return cspad2x1_one # SegGeometryCspad2x1V1(use_wide_pix_center=False)
if segname=='EPIX100:V1' : return epix2x2_one # SegGeometryEpix100V1(use_wide_pix_center=False)
if segname=='PNCCD:V1' : return segment_one # SegGeometryMatrixV1()
if segname[:4]=='MTRX' :
rows, cols, psize_row, psize_col = matrix_pars(segname)
return SegGeometryMatrixV1(rows, cols, psize_row, psize_col,\
pix_size_depth=100,\
pix_scale_size=min(psize_row, psize_col))
if segname=='JUNGFRAU:V1': return jungfrau_one # SegGeometryJungfrauV1()
#if segname=='ANDOR3D:V1' : return seg_andor3d # SegGeometryMatrixV1()
return None
#------------------------------
sgs = SegGeometryStore()
#------------------------------
#----------- TEST -------------
#------------------------------
def test_seggeom() :
import sys
from time import time
t0_sec = time()
if len(sys.argv)==1 : print 'For test(s) use command: python', sys.argv[0], '<test-number=1-5>'
elif(sys.argv[1]=='1') :
sg = sgs.Create('SENS2X1:V1', pbits=0377)
sg.print_seg_info(pbits=0377)
elif(sys.argv[1]=='2') :
sg = sgs.Create('EPIX100:V1', pbits=0377)
sg.print_seg_info(pbits=0377)
elif(sys.argv[1]=='3') :
sg = sgs.Create('PNCCD:V1', pbits=0377)
sg.print_seg_info(pbits=0377)
elif(sys.argv[1]=='4') :
sg = sgs.Create('MTRX:512:512:54:54', pbits=0377)
print 'Consumed time for MTRX:512:512:54:54 (sec) =', time()-t0_sec
sg.print_seg_info(pbits=0377)
elif(sys.argv[1]=='5') :
sg = sgs.Create('JUNGFRAU:V1', pbits=0377)
sg.print_seg_info(pbits=0377)
else : print 'Non-expected arguments: sys.argv=', sys.argv, ' use 0,1,2,...'
#------------------------------
if __name__ == "__main__" :
test_seggeom()
print 'End of test.'
#------------------------------