def main_alignment_test() :
print 'Start test in main_alignment_test()'
#path_calib = '/reg/d/psdm/CXI/cxi80410/calib/CsPad::CalibV1/CxiDs1.0:Cspad.0' # 2011-05-25
#path_calib = '/reg/d/psdm/CXI/cxi37411/calib/CsPad::CalibV1/CxiDs1.0:Cspad.0' # 2011-08-10
#path_calib = '/reg/d/psdm/CXI/cxi35711/calib/CsPad::CalibV1/CxiDs1.0:Cspad.0' # 2011-10-18
#path_calib = '/reg/neh/home1/dubrovin/LCLS/CSPadAlignment-v01/calib-cxi43312-Dsd' # 2012-01-12
#path_calib = '/reg/neh/home1/dubrovin/LCLS/CSPadAlignment-v01/calib-cxi80410-r1150-Ds1' # 2012-01-18
#path_calib = '/reg/neh/home1/dubrovin/LCLS/CSPadAlignment-v01/calib-cxi39112-r0009-Ds1' # 2012-02-17
#path_calib = '/reg/d/psdm/CXI/cxi39112/calib/CsPad::CalibV1/CxiDs1.0:Cspad.0' # 2012-02-17
#path_calib = '/reg/neh/home1/dubrovin/LCLS/CSPadAlignment-v01/calib-xpp-2012-02-26' # 2012-02-26
#path_calib = '/reg/neh/home1/dubrovin/LCLS/CSPadAlignment-v01/calib-cxi49812-r0073-Ds1' # 2012-03-08
#path_calib = '/reg/d/psdm/CXI/cxi49812/calib/CsPad::CalibV1/CxiDs1.0:Cspad.0' # 2012-03-08
#path_calib = '/reg/neh/home1/dubrovin/LCLS/CSPadAlignment-v01/calib-cxi49012-r0020-Ds1' # 2012-03-14
#path_calib = '/reg/d/psdm/CXI/cxi49012/calib/CsPad::CalibV1/CxiDs1.0:Cspad.0' # 2012-03-14
#path_calib = '/reg/d/psdm/XPP/xppcom10/calib/CsPad::CalibV1/XppGon.0:Cspad.0' # 2012-03-23 check
#path_calib = '/reg/neh/home1/dubrovin/LCLS/CSPadAlignment-v01/calib-xpp-2013-01-24' # 2013-01-24
#path_calib = '/reg/neh/home1/dubrovin/LCLS/CSPadAlignment-v01/calib-cxi64813-r0058-Ds1' # 2013-01-31
#path_calib = '/reg/d/psdm/cxi/cxi64813/calib/CsPad::CalibV1/CxiDs1.0:Cspad.0'
path_calib = '/reg/neh/home1/dubrovin/LCLS/CSPadAlignment-v01/calib-xpp-2013-01-29' # 2013-01-29
#path_calib = '/reg/neh/home1/dubrovin/LCLS/CSPadAlignment-v01/calib-cxi-ds1-2013-01-31' # 2013-02-21
#path_calib = '/reg/neh/home1/dubrovin/LCLS/CSPadAlignment-v01/calib-cxi80410-r1458-Ds1' # 2013-02-21 identic to calib-cxi-ds1-2013-01-31
#fname, runnum = '/reg/d/psdm/CXI/cxi35711/hdf5/cxi35711-r0009.h5', 9
#fname, runnum = '/reg/d/psdm/CXI/cxi37411/hdf5/cxi37411-r0080.h5', 80
#fname, runnum = '/reg/d/psdm/CXI/cxi37411/hdf5/cxi37411-r0039.h5', 39
#fname, runnum = '/reg/d/psdm/CXI/cxi80410/hdf5/cxi80410-r1150.h5', 1150
#fname, runnum = '/reg/d/psdm/CXI/cxi39112/hdf5/cxi39112-r0009.h5', 9
#fname, runnum = '/reg/d/psdm/XPP/xppcom10/hdf5/xppcom10-r1437.h5', 1437
#fname, runnum = '/reg/d/psdm/CXI/cxi49812/hdf5/cxi49812-r0073.h5', 73
#fname, runnum = '/reg/d/psdm/CXI/cxi49012/hdf5/cxi49012-r0020-raw.h5', 20
#fname, runnum = '/reg/d/psdm/CXI/cxi80410/hdf5/cxi80410-r0628.h5', 628
#fname, runnum = '/reg/d/psdm/xpp/xppcom13/hdf5/xppcom13-r0066.h5', 66
#fname, runnum = '/reg/d/psdm/CXI/cxi64813/hdf5/cxi64813-r0058.h5', 58
#fname, runnum = '/reg/d/psdm/CXI/cxi80410/hdf5/cxi80410-r1458.h5', 1458
fname, runnum = '/reg/d/psdm/xpp/xpp66213/hdf5/xpp66213-r0150.h5', 150
dsname = '/Configure:0000/Run:0000/CalibCycle:0000/CsPad::ElementV2/XppGon.0:Cspad.0/data'
#dsname = '/Configure:0000/Run:0000/CalibCycle:0000/CsPad::ElementV2/CxiDsd.0:Cspad.0/data'
#dsname = '/Configure:0000/Run:0000/CalibCycle:0000/CsPad::ElementV2/CxiDs1.0:Cspad.0/data'
event = 0
print 'Load calibration parameters from', path_calib
#calibpars = calp.calibpars.setCalibParsForPath ( run=runnum, path=path_calib )
calibpars = calp.CalibPars( path=path_calib, run=runnum )
print 'center_global:\n', calibpars.getCalibPars('center_global')
#print 'offset:\n', calibpars.getCalibPars('offset')
#calp.calibpars.printCalibPars()
#calp.calibpars.printCalibFiles ()
#calp.calibpars.printListOfCalibTypes()
#cald.calibparsdefault.printListOfCalibTypes()
#cald.calibparsdefault.printCalibParsDefault()
#cald.calibparsdefault.printCalibParsDefault('center_global')
print 'Get raw CSPad event %d from file %s \ndataset %s' % (event, fname, dsname)
ds1ev = hm.getOneCSPadEventForTest( fname, dsname, event )
#ds1ev = hm.getAverageCSPadEvent( fname, dsname, event=200, nevents=500 )
#ds1ev = hm.getAverageCSPadEvent( fname, dsname, event, nevents=10 )
#print 'ds1ev.shape = ',ds1ev.shape # should be (32, 185, 388)
#print 'ds1ev = ',ds1ev[1,:]
#print 'Subtract pedestals'
#ped_fname = '/reg/neh/home1/dubrovin/LCLS/CSPadPedestals/cspad-pedestals-cxi49812-r0072.dat' # shape = (5920, 388)
#ped_fname = '/reg/neh/home1/dubrovin/LCLS/CSPadPedestals/cspad-pedestals-cxi49012-r0008.dat' # shape = (5920, 388)
#ped_fname = '/reg/neh/home1/dubrovin/LCLS/CSPadPedestals/cspad-pedestals-cxi49012-r0038.dat' # shape = (5920, 388)
#ped_fname = '/reg/neh/home1/dubrovin/LCLS/CSPadPedestals/cspad-pedestals-cxi49012-r0027.dat' # shape = (5920, 388)
#ped_fname = '/reg/neh/home1/dubrovin/LCLS/CSPadPedestals/cspad-pedestals-xppcom10-r1435.dat' # shape = (5920, 388) low gain
#ped_fname = '/reg/d/psdm/CXI/cxi49012/calib/CsPad::CalibV1/CxiDs1.0:Cspad.0/pedestals/9-37.data' # shape = (5920, 388)
#ped_fname = '/reg/neh/home1/dubrovin/LCLS/CSPadPedestals/cspad-pedestals-cxi80410-r1453.dat' # shape = (5920, 388)
ped_fname = '/reg/neh/home1/dubrovin/LCLS/CSPadPedestals/cspad-pedestals-xpp66213-r0149.dat' # shape = (5920, 388)
#ds1ev = gm.getCSPadArrayFromFile(ped_fname)
#ds1ev -= gm.getCSPadArrayFromFile('/reg/neh/home1/dubrovin/LCLS/CSPadPedestals/cspad-pedestals-cxi49012-r0027.dat')
#ds1ev -= gm.getCSPadArrayFromFile('/reg/neh/home1/dubrovin/LCLS/CSPadPedestals/cspad-pedestals-xppcom10-r1442.dat')
ds1ev -= gm.getCSPadArrayFromFile(ped_fname)
print 'ds1ev.shape = ',ds1ev.shape # should be (32, 185, 388)
print 'Make the CSPad image from raw array'
#cspadimg = cip.CSPadImageProducer(rotation=0, tiltIsOn=True)#, mirror=True)
cspadimg = cip.CSPadImageProducer(calibpars, rotation=0, tiltIsOn=True)#, mirror=True)
#cspadimg.printInputPars()
#cspadimg.printGeometryPars()
#arr = cspadimg.getImageArrayForCSPadElement( ds1ev )
arr = cspadimg.getCSPadImage( ds1ev )
print 'Plot CSPad image'
#AmpRange = (800, 1300) # for cxi
AmpRange = (0, 10) # for cxi
#AmpRange = (-10, 50) # for xpp
#gg.plotImage(arr,range=AmpRange,figsize=(1.16*12,12))
gg.plotImageLarge(arr,range=AmpRange,figsize=(1.15*12,12))
gg.move(200,100)
gg.plotSpectrum(arr,range=AmpRange)
gg.move(50,50)
#gg.plotImageAndSpectrum(arr,range=(1,2001))
print 'To EXIT the test click on "x" in the top-right corner of each plot window.'
gg.show()
