Ticket #4222: findpeakssx.py

File findpeakssx.py, 1.2 KB (added by Owen Arnold, 9 years ago)
Line 
1
2vanadium_raw_file = "WISH00018589.RAW"
3sample_raw_file = "WISH00017986.RAW"
4
5#load workspaces
6LoadRaw(Filename=vanadium_raw_file, OutputWorkspace="vanadium");
7SmoothNeighbours(InputWorkspace="vanadium", OutputWorkspace="smoothed_vanadium",Radius=3,WeightedSum='Flat',NumberOfNeighbours=35)
8DeleteWorkspace(Workspace='vanadium')
9
10LoadRaw(Filename=sample_raw_file, OutputWorkspace="sample");
11SmoothNeighbours(InputWorkspace="sample", OutputWorkspace="smoothed_sample",Radius=3,WeightedSum='Flat',NumberOfNeighbours=35)
12DeleteWorkspace(Workspace='sample')
13
14SmoothData(InputWorkspace="smoothed_vanadium", OutputWorkspace="smoothed_vanadium",NPoints=50);
15
16Divide(LHSWorkspace="smoothed_sample", RHSWorkspace="smoothed_vanadium",OutputWorkspace="normalised_sample");
17
18DeleteWorkspace(Workspace='smoothed_vanadium')
19DeleteWorkspace(Workspace='smoothed_sample')
20
21ReplaceSpecialValues(InputWorkspace='normalised_sample',OutputWorkspace='normalised_sample',NaNValue='0',InfinityValue='10000',InfinityError='10000',BigNumberThreshold='10000',BigNumberValue='10000',BigNumberError='10000')
22FindSXPeaks(InputWorkspace='product',SignalBackground='100',Resolution='0.20000000000000001',OutputWorkspace='csp79590_1_peaks')
23