Post CTE Indirect Analysis Interface

[Michael Whitty]

Create interface for post-convert-to-energy steps for indirect instruments. From Spencer's scripts: fury, elwin, slice, msdfit, absorption correction

[Michael Whitty]

(In [5829]) re #1508 - adding interface for Fury, Elwin, Slice, MSDFit and Absorption stages of Indirect analysis steps extrapolated from Spencer Howell's scripts.

[Nick Draper]

fill in missing resolution

[Michael Whitty]

(In [5894]) re #1508 - fixing typo in python file

[Michael Whitty]

Instructions for testing:

NB: running some of these under debug mode seems to crash Python for some reason. test in release mode or on a recent installer/rpm build (>r5918).

Ok, leave the initial settings to Instrument: ISIS/IRIS, Analyser: Graphite, Reflection: 002 - ignore the other boxes on the first tab.

"Fury" - takes an "Icon" file (the output of Indirect ConvertToEnergy) and a RES file, which is pretty much the same thing but on a calibration run with a couple of post steps. Select these files (icon: irs21360_graphite002_ipg.nxs, res: irs26173_res.nxs) ELow/EWidth/EHigh are the rebin parameters. I tend to use "-0.5, 0.005, 0.5"

"Elwin" - takes irs21360_graphite002_ipg.nxs also. Window Start: -0.1, Window End: 0.1 - the result should be a workspace with a single histogram and the number of bins equal to the number of histograms in the input workspace, with the X-axis being the Q^{2 which is based off the detectors twoTheta value.}

"Slice" - takes a raw file - I'd suggest irs21360.raw which is in the test data directory. It takes also a calibration file generated in the ConvertToEnergy interface - irs26173_graphite002_calib.nxs attached. Suggested input values (taken from Spencer's scripts) for TOF ranges are: [TOF1A: 62500, TOF1B: 65000, TOF2A: 59000, TOF2B: 61500] This should IIRC give you a single-binned workspace which when you plot the bin, resembles in shape that retrieved from Elwin's plotted spectra.

"MSDFit" - takes an Elwin file (which you'll have generated by now, should be a nxs file ending with '_elw'). StartX: 0.3, EndX 1.3, check 'Plot Output', you'll get the Elwin plot and a line fitted to it, with some information in the title of the plot.

"Absorption" - takes the above-mentioned _ipg.nxs file. Select either flat plate from the combo box. Enter values of: attenuation=1.0, scattering=5.0, density=0.1 height=1.0, width=1.0, thickness=0.1, element size=0.5

Just waiting for mantidplot to build so I can create some test data, which will be attached shortly.

[Michael Whitty]

Testing file

[Michael Whitty]

Testing file

[Michael Whitty]

Testing file

[Michael Whitty]

That's the testing files added. Good luck. :)

[Stuart Campbell]

(In [5963]) Added test files for analysis testing. refs #1508

[Stuart Campbell]

All seem to run.

The only one I am not sure about is the "Slice"

It produces a workspace of 51 histograms of length 1 bin. Is this correct ? I can't figure out how to plot it ?

[Stuart Campbell]

I ended up using a different program to plot it and it looks ok.

[Stuart Campbell]

This ticket has been transferred to github issue ​2355