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]