# D_T for nine blocks when dae playing up in cycle 09/1 or 09/2
wksp="SANS2D00001674"
wksp="5976_sans_nxs"
dimdet=192
spec1=8+1+192*192
tag="frt"
#spec1=8+1
#tag="rear"
#print speclist
i=0
namelist=["A TL","D TC","G TR","B ML","E MC","H MR","C BL","F BC","I BR"]
masklist=["h0>h127,v64>v191","h0>h127,v0>v63,v128>v191","h0>h127,v0>v127","h0>h63,h128>h191,v64>v191","h0>h63,h128>h191,v0>v63,v128>v191","h0>h63,h128>h191,v0>v127","h63>h191,v64>v191","h63>h191,v0>v63,v128>v191","h63>h191,v0>v127"]
for name in namelist:
print "i= "+str(i)
# guess "orientation"
list = SANSUtility.ConvertToSpecList(masklist[i], spec1, dimdet,'0')
CropWorkspace(wksp, OutputWorkspace = "D_T "+name+tag, StartWorkspaceIndex=(spec1-1), EndWorkspaceIndex=str(spec1+192*192-2))
SANSUtility.MaskBySpecNumber("D_T "+name+tag, list)
SumSpectra("D_T "+name+tag,"D_T "+name+tag)
i=i+1
print "done"
for wext in namelist:
wsp="D_T "+wext
if wext==namelist[0]:
plot2=plotSpectrum(wsp+tag,0)
layer=plot2.activeLayer()
layer.setTitle("D_T modules")
else:
mergePlots(plot2,plotSpectrum(wsp+tag,0))
layer.setCurveTitle(i,name)
#
