from mantid.simpleapi import LoadDiffCal, mtd, LoadEmptyInstrument, CalculateDIFC, MaskBTP, CloneWorkspace, SaveDiffCal, ApplyCalibration
import matplotlib.pyplot as plt
import tube
import numpy as np
import numpy.ma as ma
tube.readCalibrationFile(
'CalibTable', '/SNS/users/rwp/corelli/tube_calibration/CalibTableNew.txt')
corelli = LoadEmptyInstrument(InstrumentName='CORELLI')
corelli = CalculateDIFC(corelli)
difc0 = corelli.extractY().flatten()
ApplyCalibration('corelli', 'CalibTable')
corelli = CalculateDIFC(corelli)
difc = corelli.extractY().flatten()
plt.plot(difc / difc0)
plt.show()
LoadDiffCal(
Filename='../cal_Si_C60/cal_Si2_47327-47334_TubeCal_sum16_mask_lt_2.cal',
InstrumentName='CORELLI',
WorkspaceName='si')
MaskBTP(Workspace='si_mask', Pixel="1-16,241-256")
LoadDiffCal(
Filename='../cal_Si_C60/cal_C60_2_47367-47382_TubeCal_sum16_mask_lt_2.h5',
InstrumentName='CORELLI',
WorkspaceName='c60')
MaskBTP(Workspace='c60_mask', Pixel="1-16,241-256")
from mantid.simpleapi import LoadDiffCal, mtd, LoadEmptyInstrument, CalculateDIFC, MaskBTP, CloneWorkspace, SaveDiffCal, ApplyCalibration
# import matplotlib.pyplot as plt
import tube
import numpy as np
import numpy.ma as ma
tube.readCalibrationFile(
'CalibTable', '/SNS/users/rwp/corelli/tube_calibration/CalibTableNew.txt')
corelli = LoadEmptyInstrument(InstrumentName='CORELLI')
corelli = CalculateDIFC(corelli)
difc0 = corelli.extractY().flatten()
ApplyCalibration('corelli', 'CalibTable')
corelli = CalculateDIFC(corelli)
difc = corelli.extractY().flatten()
# plt.plot(difc/difc0)
# plt.show()
LoadDiffCal(Filename='../cal_Si_C60/cal_Si2_47327-47334_TubeCal_sum16.h5',
InstrumentName='CORELLI',
WorkspaceName='si')
MaskBTP(Workspace='si_mask', Pixel="1-16,241-256")
LoadDiffCal(Filename='../cal_Si_C60/cal_C60_2_47367-47382_TubeCal_sum16.h5',
InstrumentName='CORELLI',
WorkspaceName='c60')
MaskBTP(Workspace='c60_mask', Pixel="1-16,241-256")
si = mtd['si_cal']
c60 = mtd['c60_cal']
