Skip to content

Commit

Permalink
Remove old fury code from IDA.py
Browse files Browse the repository at this point in the history
Refs #10747
  • Loading branch information
DanNixon committed Dec 10, 2014
1 parent 72e8ac2 commit 0d0c9ec
Showing 1 changed file with 0 additions and 77 deletions.
77 changes: 0 additions & 77 deletions Code/Mantid/scripts/Inelastic/IndirectDataAnalysis.py
Expand Up @@ -156,83 +156,6 @@ def confitSeq(inputWS, func, startX, endX, ftype, bgd, temperature=None, specMin

EndTime('ConvFit')

##############################################################################
# Fury
##############################################################################

def furyPlot(inWS, spec):
graph = mp.plotSpectrum(inWS, spec)
layer = graph.activeLayer()
layer.setScale(mp.Layer.Left, 0, 1.0)

def fury(samWorkspaces, res_workspace, rebinParam, RES=True, Save=False, Verbose=False,
Plot=False):

StartTime('Fury')
workdir = config['defaultsave.directory']
samTemp = samWorkspaces[0]
nsam,npt = CheckHistZero(samTemp)
Xin = mtd[samTemp].readX(0)
d1 = Xin[1]-Xin[0]
if d1 < 1e-8:
raise RuntimeError('Data energy bin is zero')

d2 = Xin[npt-1]-Xin[npt-2]
dmin = min(d1,d2)
pars = rebinParam.split(',')
if float(pars[1]) <= dmin:
raise RuntimeError('EWidth = ' + pars[1] + ' < smallest Eincr = ' + str(dmin))

outWSlist = []
# Process RES Data Only Once
CheckAnalysers(samTemp, res_workspace, Verbose)
nres,nptr = CheckHistZero(res_workspace)
if nres > 1:
CheckHistSame(samTemp,'Sample', res_workspace, 'Resolution')

tmp_res_workspace = '__tmp_' + res_workspace
Rebin(InputWorkspace=res_workspace, OutputWorkspace=tmp_res_workspace, Params=rebinParam)
Integration(InputWorkspace=tmp_res_workspace, OutputWorkspace='res_int')
ConvertToPointData(InputWorkspace=tmp_res_workspace, OutputWorkspace=tmp_res_workspace)
ExtractFFTSpectrum(InputWorkspace=tmp_res_workspace, OutputWorkspace='res_fft', FFTPart=2)
Divide(LHSWorkspace='res_fft', RHSWorkspace='res_int', OutputWorkspace='res')
for samWs in samWorkspaces:
(direct, filename) = os.path.split(samWs)
(root, ext) = os.path.splitext(filename)
Rebin(InputWorkspace=samWs, OutputWorkspace='sam_data', Params=rebinParam)
Integration(InputWorkspace='sam_data', OutputWorkspace='sam_int')
ConvertToPointData(InputWorkspace='sam_data', OutputWorkspace='sam_data')
ExtractFFTSpectrum(InputWorkspace='sam_data', OutputWorkspace='sam_fft', FFTPart=2)
Divide(LHSWorkspace='sam_fft', RHSWorkspace='sam_int', OutputWorkspace='sam')
# Create save file name
savefile = getWSprefix(samWs) + 'iqt'
outWSlist.append(savefile)
Divide(LHSWorkspace='sam', RHSWorkspace='res', OutputWorkspace=savefile)
#Cleanup Sample Files
DeleteWorkspace('sam_data')
DeleteWorkspace('sam_int')
DeleteWorkspace('sam_fft')
DeleteWorkspace('sam')
# Crop nonsense values off workspace
bin = int(math.ceil(mtd[savefile].blocksize()/2.0))
binV = mtd[savefile].dataX(0)[bin]
CropWorkspace(InputWorkspace=savefile, OutputWorkspace=savefile, XMax=binV)
if Save:
opath = os.path.join(workdir, savefile+'.nxs') # path name for nxs file
SaveNexusProcessed(InputWorkspace=savefile, Filename=opath)
if Verbose:
logger.notice('Output file : '+opath)
# Clean Up RES files
DeleteWorkspace(tmp_res_workspace)
DeleteWorkspace('res_int')
DeleteWorkspace('res_fft')
DeleteWorkspace('res')
if Plot:
specrange = range(0,mtd[outWSlist[0]].getNumberHistograms())
furyPlot(outWSlist, specrange)
EndTime('Fury')
return outWSlist

##############################################################################
# FuryFit
##############################################################################
Expand Down

0 comments on commit 0d0c9ec

Please sign in to comment.