Я студент-физик, поэтому к кодированию и статистике я все еще привыкаю. Мой учитель создал эту библиотеку, поэтому она может показаться незнакомой, но я строю кривую Гаусса для радиоактивного распада. Мне просто интересно, какое правильное значение хи-квадрат, потому что оно говорит мне о двух значениях. Также совершенство подгонки.
import numpy as np
import LT.box as B
import LT_Fit.parameters as P
import LT_Fit.gen_fit as G
import scipy.stats as stats
file_name='HC.data'
f = B.get_file(file_name)
C = B.get_data(f,'T')
mean = np.average(C)
stdDev = np.std(C)
low = min(C)
high = max(C)
Bw = 15
#76 122
h = B.histo(C, range =(70, 130), bins = Bw)
h.plot()
B.pl.show()
h.fit()
B.pl.xlabel("Count [Hz]", fontsize=12)
B.pl.ylabel("Frequency", fontsize=12)
B.pl.title("High Count", fontsize=16)
hx = h.bin_center
hy = h.bin_content
mean = P.Parameter(mean,'mu')
sigma = P.Parameter(stdDev,'sigma')
height = P.Parameter(10,'height')
binwidth = 6 / Bw
scale_factor = len(C) * binwidth
def gaussFunct(x):
value = stats.norm.pdf((x-mean())/sigma())*height()
return value
fit = G.genfit(gaussFunct,[mean,sigma, height], x = hx, y = hy)
B.plot_line(fit.xpl, fit.ypl, color='red', linewidth=2)
