The data to be fitted is I_DA(t). I am less concerned about the convolution with g(t) and the f terms, although any comments on these would be appreciated. The fit needs to determine the parameters a and b. c can be determined by normalisation.
The question is really how to deal with the intgral over R. I tried modifying this example: