from paida import *

analysisFactory = IAnalysisFactory.create()
treeFactory = analysisFactory.createTreeFactory()
tree = treeFactory.create()

functionFactory = analysisFactory.createFunctionFactory(tree)
func = functionFactory.createFunctionFromScript('G2', 2, 'amplitude * exp(-0.5 * (((x[0] - meanX) / sigmaX)**2 + ((x[1] - meanY) / sigmaY)**2))', 'amplitude,meanX,sigmaX,meanY,sigmaY', 'DoubleGaussian')
functionFactory.catalog().add('2DGaussian', func)

func.setParameter('amplitude', 8.0)
func.setParameter('meanX', 0.0)
func.setParameter('meanY', 0.0)
func.setParameter('sigmaX', 3.0)
func.setParameter('sigmaY', 3.0)

plotterFactory = analysisFactory.createPlotterFactory()
plotter = plotterFactory.create('PAIDA Plotter')
region = plotter.createRegion()

region.setTitle('User created function (cataloged as 2DGaussian)')
region.style().setParameter('showTitle', 'true')
region.plot(func, 'minX=-8.0;maxX=8.0;minY=-8.0;maxY=8.0;minZ=-3.0;maxZ=10.0')

### Wait.
dummy = raw_input('Hit any key.')