class TestPeak(Test): def setUp(self): Test.setUp(self) self.p = GaussianPeakAndBackground() def test__init__(self): self.check__init__('peak', ('pos','offset','top', 'fwhm', 'residual'), 'pos') def test_process(self): pos = 15.014649472869898 offset = 0 top = 4.939879027681013 fwhm = 5.224729872484415 result = self.p._process(self.x, self.peak) expected = (pos, offset, top, fwhm) self.assert_(close(result, expected),"%s\n is not close to expected:\n%s"%(`result`,`expected`)) result = self.p._process(self.x, (self.peak + 123.)) expected = (pos, offset+123, top, fwhm) self.assert_(close(result, expected),"%s\n is not close to expected:\n%s"%(`result`,`expected`)) result = self.p._process(self.x, ((self.peak * -1) + 123)) expected = (pos, offset + 123, -top, fwhm) self.assert_(close(result, expected),"%s\n is not close to expected:\n%s"%(`result`,`expected`))
from gdascripts.analysis.datasetprocessor.oned.MinPositionAndValue import MinPositionAndValue from gdascripts.analysis.datasetprocessor.oned.CentreOfMass import CentreOfMass from gdascripts.analysis.datasetprocessor.oned.GaussianPeakAndBackground import GaussianPeakAndBackground from gdascripts.scan.process.ScanDataProcessor import ScanDataProcessor from gdascripts.scan.process.ScanDataProcessorResult import go from gdascripts.scan import specscans from gdascripts.scan import gdascans from gda.jython.commands.GeneralCommands import alias print "Setting up scan data processor, scan_processor" scan_processor = ScanDataProcessor([ MaxPositionAndValue(), MinPositionAndValue(), CentreOfMass(), GaussianPeakAndBackground() ], globals()) alias("go") print "Creating scan commands:" ascan = specscans.Ascan([scan_processor]) a2scan = specscans.A2scan([scan_processor]) a3scan = specscans.A3scan([scan_processor]) mesh = specscans.Mesh([scan_processor]) dscan = specscans.Dscan([scan_processor]) d2scan = specscans.D2scan([scan_processor]) d3scan = specscans.D3scan([scan_processor]) scan = gdascans.Scan([scan_processor]) rscan = gdascans.Rscan([scan_processor]) alias('ascan') print ascan.__doc__.split('\nUSAGE:\n\t\n ')[1]
from gdascripts.analysis.datasetprocessor.oned.MinPositionAndValue import MinPositionAndValue from gdascripts.analysis.datasetprocessor.oned.CentreOfMass import CentreOfMass from gdascripts.analysis.datasetprocessor.oned.GaussianPeakAndBackground import GaussianPeakAndBackground from gdascripts.analysis.datasetprocessor.oned.GaussianEdge import GaussianEdge from gdascripts.scan.process.ScanDataProcessor import ScanDataProcessor from gdascripts.scan import specscans from gdascripts.scan import gdascans from gda.jython.commands.GeneralCommands import alias print "Setting up scan data processor, scan_processor" ge = GaussianEdge() ge.smoothwidth = 5 ge.plotPanel = "Edge Fit Plot" gpab = GaussianPeakAndBackground() gpab.plotPanel = "Peak Fit Plot" scan_processor = ScanDataProcessor( [MaxPositionAndValue(), MinPositionAndValue(), gpab, ge], globals()) go = scan_processor.go alias("go") print "Creating spec-like commands:" ascan = specscans.Ascan([scan_processor]) #a2scan = specscans.A2scan([scan_processor]) #a3scan = specscans.A3scan([scan_processor]) #mesh = specscans.Mesh([scan_processor]) #dscan = specscans.Dscan([scan_processor]) #d2scan = specscans.D2scan([scan_processor]) #d3scan = specscans.D3scan([scan_processor])
from gdascripts.analysis.datasetprocessor.oned.MaxPositionAndValue import MaxPositionAndValue from gdascripts.analysis.datasetprocessor.oned.MinPositionAndValue import MinPositionAndValue from gdascripts.analysis.datasetprocessor.oned.CentreOfMass import CentreOfMass from gdascripts.analysis.datasetprocessor.oned.GaussianPeakAndBackground import GaussianPeakAndBackground from gdascripts.scan.process.ScanDataProcessor import ScanDataProcessor from gdascripts.scan.process.ScanDataProcessorResult import go from gdascripts.scan import specscans from gdascripts.scan import gdascans from gda.jython.commands.GeneralCommands import alias print "Setting up scan data processor, scan_processor" scan_processor = ScanDataProcessor( [MaxPositionAndValue(),MinPositionAndValue(),CentreOfMass(), GaussianPeakAndBackground() ], globals() ) alias("go") print "Creating scan commands:" ascan = specscans.Ascan([scan_processor]) a2scan = specscans.A2scan([scan_processor]) a3scan = specscans.A3scan([scan_processor]) mesh = specscans.Mesh([scan_processor]) dscan = specscans.Dscan([scan_processor]) d2scan = specscans.D2scan([scan_processor]) d3scan = specscans.D3scan([scan_processor]) scan=gdascans.Scan([scan_processor]) rscan=gdascans.Rscan([scan_processor]) alias('ascan'); print ascan.__doc__.split('\nUSAGE:\n\t\n ')[1] alias('a2scan');print a2scan.__doc__.split('\nUSAGE:\n\t\n ')[1] alias('a3scan');print a3scan.__doc__.split('\nUSAGE:\n\t\n ')[1] alias('mesh');print mesh.__doc__.split('\nUSAGE:\n\t\n ')[1] alias('dscan');print dscan.__doc__.split('\nUSAGE:\n\t\n ')[1] alias('d2scan');print d2scan.__doc__.split('\nUSAGE:\n\t\n ')[1]
def setUp(self): Test.setUp(self) self.p = GaussianPeakAndBackground()
def doGaussianTest(self, path): pwt = ProcessWithTimeout(path, GaussianPeakAndBackground()) pwt.start() pwt.join(10) self.assertFalse(pwt.isAlive(), "Processing still running after 10 seconds. WARNING: may effect future tests or test machine?")
def setUp(self): self.peak = GaussianPeakAndBackground()
from gdascripts.analysis.datasetprocessor.oned.MinPositionAndValue import MinPositionAndValue from gdascripts.analysis.datasetprocessor.oned.CentreOfMass import CentreOfMass from gdascripts.analysis.datasetprocessor.oned.GaussianPeakAndBackground import GaussianPeakAndBackground from gdascripts.analysis.datasetprocessor.oned.GaussianEdge import GaussianEdge from gdascripts.scan.process.ScanDataProcessor import ScanDataProcessor from gdascripts.scan import specscans from gdascripts.scan import gdascans from gda.jython.commands.GeneralCommands import alias print "Setting up scan data processor, scan_processor" ge=GaussianEdge() ge.smoothwidth=5 ge.plotPanel = "Edge Fit Plot" gpab = GaussianPeakAndBackground() gpab.plotPanel = "Peak Fit Plot" scan_processor = ScanDataProcessor( [MaxPositionAndValue(),MinPositionAndValue(),CentreOfMass(), gpab, ge ], globals() ) go = scan_processor.go alias("go") print "Creating spec-like commands:" ascan = specscans.Ascan([scan_processor]) #a2scan = specscans.A2scan([scan_processor]) #a3scan = specscans.A3scan([scan_processor]) #mesh = specscans.Mesh([scan_processor]) #dscan = specscans.Dscan([scan_processor]) #d2scan = specscans.D2scan([scan_processor]) #d3scan = specscans.D3scan([scan_processor]) alias('ascan');print ascan.__doc__.split('\n')[3] #alias('a2scan');print a2scan.__doc__.split('\n')[3]