FindDetectorsOutsideLimits v1¶

import numpy as np
ws = CreateSampleWorkspace(BankPixelWidth=10,NumBanks=1)

#create dome dead and noisy detectors
deadDetArray=[0.0] * ws.blocksize()
noisyDetArray= [100.0] * ws.blocksize()

for i in range(0,ws.getNumberHistograms(),5):
    ws.setY(i,np.array(deadDetArray))
    ws.setY(i+1,np.array(noisyDetArray))

print "With just the default LowThreshold of 0"
(wsOut,NumberOfFailures)=FindDetectorsOutsideLimits(ws)
print "%i spectra were outside the limits." % NumberOfFailures
print

print "With a High and LowThreshold, as well as restricting the XRange to consider"
(wsOut2,NumberOfFailures)=FindDetectorsOutsideLimits(ws, HighThreshold=1000,
    LowThreshold=0, RangeLower=200, RangeUpper=10000)
print "%i spectra were outside the limits." % NumberOfFailures

mtd.clear()

ws = CreateSimulationWorkspace('MARI','0,1,10')
nh = ws.getNumberHistograms()
for ind in xrange(nh):
    y = ws.dataY(ind)
    if ind>=100 and ind < 300:
        y.fill(100)
    else:
        y.fill(1)

mws1,nMasked1 = FindDetectorsOutsideLimits(ws,100)
mws2,nMasked2 = FindDetectorsOutsideLimits(ws,100,startWorkspaceIndex = 200)

print "****************************************"
print "full mask ws has {0} masked detectors".format(nMasked1)
print "part mask ws  has {0} masked detectors".format(nMasked2)
print "****************************************"
selected_spec = [99,100,199,200,299,300]
for spec in selected_spec:
    print "full mask ws Spec N{0} is masked: {1}".format(spec,mws1.readY(spec)[0]>0.5)
    print "part mask ws Spec N{0} is masked: {1}".format(spec,mws2.readY(spec)[0]>0.5)
print "****************************************"