deletions | additions       diff --git a/res_fitting.md b/res_fitting.md  index a769002..f376b66 100644  --- a/res_fitting.md  +++ b/res_fitting.md 
...
\[  F_i(t) = b(\mathbf{q}_i, t) + \sum_{k=0}^{m} f(\mathbf{p}_{i,k},t) + W + R  \]  , where \(b\) is a n-th order polynomial with \(\mathbf{q}_i\) being the coefficients for the *i*-th pixel, summation is performed over all *m* events in the pixel, \(W\) consists of the noise and R \(R\)  is remaining residual not captured in the baseline nor events. Ideally, \(R=0\), but in achieving this we are limited by the accuracy of the event region detection (we cannot fit what we do not detect) and whether or not our fitting function is general enough to be able to approximate various types of events.  Because it is not know which part of the signal is the event and which is the baseline the first fit also has to estimate the baseline properties.  Signal in the candidate region is fitted with an extended fit function (Figure \ref{fig:fig2}) that also depends on relaxation baseline \(B\) and baseline offset \(C\). The \(C\) parameter allows for the possibility of an elevated background before the release event.