Skip to content

Commit

Permalink
Added changes. This refs #4303
Browse files Browse the repository at this point in the history
  • Loading branch information
JeanBilheux committed Jan 13, 2012
1 parent 9ceb8aa commit 233e090
Showing 1 changed file with 12 additions and 15 deletions.
Original file line number Diff line number Diff line change
Expand Up @@ -250,21 +250,18 @@
OutputWorkspace='NormWks')

#collapse data to 1 spectrum (peak region)
mt_NormWks = mtd['NormWks']
_tof_axis = mt_NormWks.readX(0)[:]
counts_vs_tof = zeros(len(_tof_axis))
from_norm_peak = list_norm_peak[0][0]
to_norm_peak = list_norm_peak[0][1]
NormPeakRange = arange(to_norm_peak-from_norm_peak+1) + from_norm_peak
for x in NormPeakRange:
print x
counts_vs_tof += mt_NormWks.readY(int(x))[:]
CreateWorkspace('NormWks', DataX=_tof_axis, DataY=counts_vs_tof, DataE=counts_vs_tof, Nspec=1)

# #replace NaN values
# ReplaceSpecialValues(InputWorkspace='NormWks',
# Outputworkspace='NormWks',
# NaNValue=0)
Integration('NormWks','IntegratedNormWks',)

# mt_NormWks = mtd['NormWks']
# _tof_axis = mt_NormWks.readX(0)[:]
# counts_vs_tof = zeros(len(_tof_axis))
# from_norm_peak = list_norm_peak[0][0]
# to_norm_peak = list_norm_peak[0][1]
# NormPeakRange = arange(to_norm_peak-from_norm_peak+1) + from_norm_peak
# for x in NormPeakRange:
# print x
# counts_vs_tof += mt_NormWks.readY(int(x))[:]
# CreateWorkspace('NormWks', DataX=_tof_axis, DataY=counts_vs_tof, DataE=counts_vs_tof, Nspec=1)

#### divide data by normalize histo workspace
Divide(LHSWorkspace='DataWks',
Expand Down

0 comments on commit 233e090

Please sign in to comment.