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]

Navigation