golicnik schnell stroberg example.ipynb

If you know the t_q, you can use this simple Km calculator, which seems to work OK even for quite noisy readouts:

ests = []
signal = golicnik+np.random.normal(0, 0.0001, size=(len(golicnik),))
dp_ = np.gradient(signal)

for d in np.arange(1, 100):
    i = d2p.argmax()-d
    j = d2p.argmax()+d
    dyi = dp_[i]
    dyj = dp_[j]
    si = signal[i]
    sj = signal[j]
    

    km_est = ((dyj-dyi)* si*sj) / ((dyi*sj) - (dyj*si))
    if km_est<0:
        continue
    print(km_est)
    ests.append(km_est)


plt.hist(ests,bins=55)
plt.axvline(np.sqrt(np.power(ests, 2).mean()))