hiro/fpraktikum
1
0
Fork 0
mirror of https://github.com/vale981/fpraktikum synced 2025-03-05 09:31:44 -05:00
Code Issues Projects Releases Packages Wiki Activity Actions

remove scaling

Browse source
This commit is contained in:
hiro98 2019-11-30 19:02:00 +01:00
parent cc98779c45
commit 987d9d0e06
2 changed files with 148 additions and 159 deletions
Show stats Download patch file Download diff file Expand all files Collapse all files

296
CS/auswertung/Untitled.ipynb
View file

File diff suppressed because one or more lines are too long

11
CS/auswertung/utility.py
View file

@ -60,7 +60,7 @@ def plot_spec(spec, is_relative=False, offset=0, ticks=10, save=None, **pyplot_a
ax.grid(which='minor', alpha=.3)
ax.grid(which='major', alpha=.8)
ax.set_xlabel('Kanal')
ax.set_ylabel(('Normierte' if is_relative else '') + 'Zaehlrate')
ax.set_ylabel(('Normierte' if is_relative else '') + 'Ereignisszahl')
ax.set_xlim(chans[0], chans[-1])
ax.set_ylim(spec.min(), spec.max())
@ -79,17 +79,16 @@ def gauss_offset(x, mu, sigma, A, O):
def calibrate_peak(spec, start, end, save=None):
slc = spec[start:end]
slc = slc-slc.min()
slc = slc/slc.max()
chan = channels(spec)[start:end]
opt, d_opt = curve_fit(gauss_offset, chan,
slc, p0=((start+end)/2, (end-start)/2, 1, 0),
slc, p0=((start+end)/2, (end-start)/2, 1*slc.max(), 0),
sigma=1/2*np.ones_like(chan), absolute_sigma=True,
bounds=([start, -np.inf, 1/2, 0], [end, np.inf, 1, 1/2]))
bounds=([start, -np.inf, 1/2*slc.max(), 0], [end, np.inf, 1*slc.max(), 1/2*slc.max()]))
d_opt = np.sqrt(np.diag(d_opt))
fig, ax = plot_spec(slc, is_relative=True, offset=start, label='Gemessen')
fig, ax = plot_spec(slc, offset=start, label='Gemessen')
ax.plot(chan, gauss_offset(chan, *opt), label='Peak Fit', linestyle='--')
ax.legend()