from rsf.proj import *

Flow('sin',None,'math n1=90000 label1=Sample d1=1 output="sin(0.0001*x1)" ')
Plot('sin','window j1=100 | graph title=Signal')

Flow('noi','sin','noise seed=2011 var=0.01')
Plot('noi','window j1=100 | graph title=Signal+Noise')

Result('noi','sin noi','OverUnderAniso')

pars = {
    'default': '',
    'variance': 'var=0.11 mean=1',
    'uniform':  'type=n',
    'range':    'type=n mean=0.5 range=2'
}

for case in pars.keys():
    Flow(case,'sin','noise rep=y seed=2011  %s ' % pars[case])
    Plot(case,'window j1=100 | graph title="%s"' % pars[case])

    hist = case+'-hist'
    Plot(hist,case,
         '''
         histogram n1=61 o1=-3 d1=0.1 |
         dd type=float | scale axis=1 |
         bargraph title="Normalized Distribution" label1=Value 
         ''')

    Result(case,[case,hist],'OverUnderAniso')

End()