Simulations

Jump to bottom Edit New page

John Robinson edited this page Nov 1, 2015 · 8 revisions

Globalcurve can be used in simulation mode as well as curve fitting mode.

in default.in set simulation to t (true).
run globalcurve as usual
enter the .ans file name that contains the model you wish to simulate
give the tcal, lims1, lims2, lims3, which appear in the header of .glo files
enter the full width half maximum of the impulse response function (IRF), also called the lamp function (in time scale used for tcal, nanoseconds). The IRF is modeled as a Gaussian.
Add Poisson noise?

Enter answer filename: exp1.ans
 Enter simulation pearameters for expt #       1
 Enter tcal, lims1, lims2, lims3:
0.016,55,350,70
 Enter fwhm of the lamp function (in units of tcal):
0.1
 Simulation mode:
Add noise (t,f)? t

Figure 1. Plot of a simulated 1 exponential decay

Validate the simulation

change in default.in, change simulation to false.
change values in .ans file to give a different set of starting parametrs
you may need to change q-shift in the .glo file.
run globalcurve
recovered parameters are spot on
chi^2 should be 1. Autocorrelation should be flat with y scale +/- 0.1