# 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) #