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']