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

GlobalCurve

This is script-driven global non-linear least squares parameter estimation written in Fortran 90. The software was originally written to analyze time resolved fluorescence data (decay). The observed fluorescence decay is the intrinsic decay convolved with the instrument resonse function (IRF).

Quirks

  • One must have global.in file (a .conf or .ini file) in the working directory of the .ans file. This controls the nature of the analysis, like the target change in chi^2 for stopping the analysis. For this reason, the .in file must be local.
  • When data needs to be pre- or post-processed in some way during the anslysis, this is specified in the data file itself. Specification for processing the data is in the header of the data file.

Limitations

  • This is compiled softare. Additional models must be compiled.
  • This was written using the Digital/Compac/Fujitsu/Intel F90 compiler, which uses digital extensions of the F90 standared. This restricts one to development on DOS.
  • Srcripts (usually called .ans files) are written in a custom format. These days, one would write in YAML or similar.
  • Software is very unforgiving of errors in data format or script files.

Installation

  • install globalcurve in c:\bin\
  • install pgplot plotting library (dependency) in c:\

Set system Path (environement variable) to recognize the globalcurve executable. Value on my computer is

%SystemRoot%\system32;%SystemRoot%;%SystemRoot%\System32\Wbem;%SYSTEMROOT%\System32\WindowsPowerShell\v1.0\;C:\bin\GlobalCurve

pgplot

install at c:

Usage

In windows, use PowerShell. To access documents, which are shared using VMWARE, cd z:. cd to directory of interest.

> globalcurve

Analysis examples

2 exponential fitting of TCSPC data

Data file

  3.4100000E-02          50         700          80
            0.0            1.0            1.0
            0.0            2.0            1.4
            1.0            3.0            1.7
            0.0            3.0            1.7
            0.0            1.0            1.0
            0.0            3.0            1.7
            0.0            3.0            1.7
  • data format is that of a .prn file (can export from excel)
  • windows formated EOL: CR+LF
  • usually given a .glo extension.
  • Format: Columns 1 and 2 are IRF and #counts in the donor channel, respectively. Rows are channels 0, 1, 2, .... Column 3, which specifies the precision of the data (sd) is optional. When omitted, the analysis assumed Poisson noise, where sd = sqrt(n). The header contains, successively,
    • time/channel (ns/chan)
    • start channel of the convolution
    • stop channel of chi^2 calculation
    • start channel of the chi^2 calculation

Ans file (analysis script)

Simulations

  • in default.in set simulation to 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
0.016,55,350,103
  • enter the fwhm of the impulse resonse curve. Here it is modeled as a gaussian.
0.05
  • add nose?
y

Clone this wiki locally