| 1 | # Basic parameters |
|---|
| 2 | filename = "TOPAZ_3131_event.nxs" |
|---|
| 3 | #Name of the workspaces to create |
|---|
| 4 | ws = "TOPAZ_3131" |
|---|
| 5 | |
|---|
| 6 | # Load the original data |
|---|
| 7 | LoadEventNexus(Filename=filename,OutputWorkspace=ws) |
|---|
| 8 | |
|---|
| 9 | # Convert to reciprocal space, in the sample frame |
|---|
| 10 | ConvertToDiffractionMDWorkspace(InputWorkspace=ws,OutputWorkspace=ws+'_MD', |
|---|
| 11 | OutputDimensions='Q (lab frame)',LorentzCorrection='1') |
|---|
| 12 | |
|---|
| 13 | # Find peaks |
|---|
| 14 | FindPeaksMD(InputWorkspace=ws+'_MD',MaxPeaks='50',OutputWorkspace=ws+'_peaks') |
|---|
| 15 | |
|---|
| 16 | # Find the UB matrix using the peaks and known lattice parameters |
|---|
| 17 | FindUBUsingLatticeParameters(PeaksWorkspace=ws+'_peaks',a='10.3522',b='6.0768',c='4.7276', |
|---|
| 18 | alpha='90',beta='90',gamma='90', NumInitial='20', Tolerance='0.12') |
|---|
| 19 | # And index to HKL |
|---|
| 20 | IndexPeaks(PeaksWorkspace=ws+'_peaks', Tolerance='0.12') |
|---|
| 21 | |
|---|
