1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
|
#!/usr/bin/env python
import os.path, sys
from fityk import Fityk
class GaussianFitter(Fityk):
def __init__(self, filename):
Fityk.__init__(self)
if not os.path.isfile(filename):
raise ValueError("File `%s' not found." % filename)
self.filename = filename
self.execute("@0 < '%s'" % filename)
print "Data info:", self.get_info("data", 0)
def run(self):
self.execute("guess Gaussian")
gauss = self.all_functions()[-1] # the last function (just created)
print "initial Gaussian center: %g" % gauss.get_param_value("center")
print "Fitting %s ..." % self.filename
self.execute("fit")
print "WSSR=", self.get_wssr()
print "Gaussian center: %g" % gauss.get_param_value("center")
def save_session(self, filename):
self.execute("info state >'%s'" % filename)
f = Fityk()
print f.get_info("version")
print "ln(2) =", f.calculate_expr("ln(2)")
del f
g = GaussianFitter(os.path.join(os.path.dirname(sys.argv[0]), "nacl01.dat"))
g.run()
g.save_session("tmp_save.fit")
|
