Source code :: scintilator1
[Return]
[Download]from numpy import *
from matplotlib.pyplot import *
# clear
hh = 5.5
h = 25 #pocet hladin, timto cislem vynasobim nejvyssi hodnotu a pak budu zaokrouhlovat, tim se vytvori hladiny blizke hodnoty se zaokrouhli na jednu hladinu
y = loadtxt('HXR_') #load as 1D array
z = zeros(h * amax(y) + 1) #1D array
u = h * y
v = round(u)
for k in xrange(1, 50000): #iterate from 1 to 49999
if (u[k] > u[k-1]) && (u[k+1] <= u[k]):
if v[k]>= 1:
z[v[k]] += 1
figure()
w = linspace(hh * len(z), len(z))
for j in xrange(len(w)):
a(j) = (662/345) * w(j);
end
r = bar(a, z);
set(r,'EdgeColor',[0 1 0]);
hold off;
xlabel('E [keV]');
ylabel('count of pulses');
print('-dpng', '-r300', 'spectrum.png');
Z = sum(z);
exit[Return]