import math

LoadISISNexus(Filename='OFFSPEC00044684.nxs', OutputWorkspace='TOF_44684')
LoadISISNexus(Filename='OFFSPEC00044683.nxs', OutputWorkspace='TRANS_44683')
LoadISISNexus(Filename='OFFSPEC00044685.nxs', OutputWorkspace='TOF_44685')


def LambdaToQ(Wavelength, ThetaInRadians):
    return ((4 * math.pi)  * math.sin(ThetaInRadians)) / Wavelength

def QToLambda(Q, ThetaInRadians):
    return ((4 * math.pi)  * math.sin(ThetaInRadians)) / Q


#FirstWavelengthMin = 1.2
#FirstWavelengthMax = 12
#FirstWavelengthMin = 1.2
#FirstWavelengthMax = 12
#FirstMomentumTransferMin = LambdaToQ(Wavelength=WavelengthMax, ThetaInRadians= math.radians(0.3))  #0.0028800000000000002
#FirstMomentumTransferMax = LambdaToQ(Wavelength=WavelengthMin,ThetaInRadians=  math.radians(0.3))  #0.039320000000000001
#SecondMomentumTransferMin = LambdaToQ(Wavelength=WavelengthMax, ThetaInRadians= math.radians(0.5)) #0.0060200000000000002
#SecondMomentumTransferMax = LambdaToQ(Wavelength=WavelengthMin, ThetaInRadians= math.radians(0.5)) # 0.061289999999999997

FirstMomentumTransferMin = 0.0028800000000000002
FirstMomentumTransferMax = 0.039320000000000001
SecondMomentumTransferMin = 0.0060200000000000002
SecondMomentumTransferMax = 0.061289999999999997

FirstWavelengthMin = QToLambda(Q=FirstMomentumTransferMax, ThetaInRadians=math.radians(0.3))
FirstWavelengthMax = QToLambda(Q=FirstMomentumTransferMin, ThetaInRadians=math.radians(0.3))
SecondWavelengthMin = QToLambda(Q=SecondMomentumTransferMax, ThetaInRadians=math.radians(0.5))
SecondWavelengthMax = QToLambda(Q=SecondMomentumTransferMin, ThetaInRadians=math.radians(0.5))

ReflectometryReductionOneAuto(InputWorkspace='TOF_44684',
                                                    ProcessingInstructions='390-410',
                                                    ThetaIn=0.29999999999999999,
                                                    DetectorCorrectionType='RotateAroundSample',
                                                    WavelengthMin=FirstWavelengthMin,
                                                    WavelengthMax=FirstWavelengthMax,
                                                    I0MonitorIndex=1,
                                                    MonitorBackgroundWavelengthMin=15,
                                                    MonitorBackgroundWavelengthMax=20,
                                                    MonitorIntegrationWavelengthMin=2,
                                                    MonitorIntegrationWavelengthMax=14,
                                                    FirstTransmissionRun='TRANS_44683',
                                                    MomentumTransferMin=FirstMomentumTransferMin,
                                                    MomentumTransferStep=0.02,
                                                    MomentumTransferMax=FirstMomentumTransferMax,
                                                    ScaleFactor=1,
                                                    OutputWorkspaceBinned='IvsQ_binned_44684_q_range',
                                                    OutputWorkspace='IvsQ_44684',
                                                    OutputWorkspaceWavelength='IvsLam_44684_q_range', Diagnostics=1)

ReflectometryReductionOneAuto(InputWorkspace='TOF_44685',
                                                    ProcessingInstructions='390-410',
                                                    ThetaIn=0.5,
                                                    DetectorCorrectionType='RotateAroundSample',
                                                    WavelengthMin=SecondWavelengthMin,
                                                    WavelengthMax=SecondWavelengthMax,
                                                    I0MonitorIndex=1,
                                                    MonitorBackgroundWavelengthMin=15,
                                                    MonitorBackgroundWavelengthMax=20,
                                                    MonitorIntegrationWavelengthMin=2,
                                                    MonitorIntegrationWavelengthMax=14,
                                                    FirstTransmissionRun='TRANS_44683',
                                                    MomentumTransferMin=SecondMomentumTransferMin,
                                                    MomentumTransferStep=0.02,
                                                    MomentumTransferMax=SecondMomentumTransferMax,
                                                    ScaleFactor=1,
                                                    OutputWorkspaceBinned='IvsQ_binned_44685_q_range',
                                                    OutputWorkspace='IvsQ_44685_q_range',
                                                    OutputWorkspaceWavelength='IvsLam_44685_q_range', Diagnostics=1)

Stitch1DMany(InputWorkspaces='IvsQ_binned_44684_q_range,IvsQ_binned_44685_q_range',
                       OutputWorkspace='IvsQ_44684_44685_q_range',
                       Params='-0.02',
                       OutScaleFactors='0.962947')


def extra():
    LoadISISNexus(Filename='/home/ejb/Downloads/SampleData-ISIS/OFFSPEC00044683.nxs', OutputWorkspace='TRANS_44683')
    CreateTransmissionWorkspaceAuto(FirstTransmissionRun='TRANS_44683', ProcessingInstructions='390-410', WavelengthMin=1.2, WavelengthMax=12, I0MonitorIndex=1, MonitorBackgroundWavelengthMin=15, MonitorBackgroundWavelengthMax=20, MonitorIntegrationWavelengthMin=2, MonitorIntegrationWavelengthMax=14, StartOverlap=10, EndOverlap=12, OutputWorkspace='test3_normalised')
    RenameWorkspace(InputWorkspace='test3_normalised', OutputWorkspace='test3_trans')
    RebinToWorkspace(WorkspaceToRebin='test3_trans', WorkspaceToMatch='_LS_2_norm_monitor', OutputWorkspace='test3_trans_rebin')
    Divide(LHSWorkspace='_LS_2_norm_monitor', RHSWorkspace='test3_trans_rebin', OutputWorkspace='test3_norm_trans')
    CropWorkspace(InputWorkspace='test3_norm_trans', OutputWorkspace='test3_cropped', XMax=11.789999999999999)
    ConvertUnits(InputWorkspace='test3_cropped', OutputWorkspace='test3_converted', Target='MomentumTransfer')