like_lc.pl

Typing 
like_lc.pl -help 
will display this help information.

Example perl script that generates Fermi LAT light curves using 
likelihood analysis. This is not an official part of the Fermi analysis
software.

WARNING! The models and other parameters used here
are for DEMONSTRATION purposes only! 
You should carefully evaluate the choice of parameters
for your own particular needs!

To use this script you must be set up to 
run both the Fermi science tools and HEASOFT.
This script checks the environment variables LHEASOFT,
PFILES, and FERMI_DIR as a check on this.

In addition, if you are using the auto-generation option to automatically
make a model file from the 1FGL catalog, the environment variable MY_FERMI_DIR
needs to be set. This is where the 1FGL catalog and the
make1FGLxml.py script must be located if used.

The model file generated by make1FGLxml.py uses the parameters
in the 1FGL catalog. In this example only the parameters of the
source of interest are left free, others are held fixed at their
catalog values. Spectra are power laws. This may very well
not be optimum!

Among the inputs prompted for are two files: 
slist.dat
	This should contain a list of source names and their coordinates
	in the form:
	1FGLJ0240.5+6113 40.1317, 61.2281
	i.e. the 1FGL catalog name must be given because if you will use
	make1FGLxml.py to generate the model file.
	Note that there is _no_ space in the source name.
	i.e. \"1FGLJ0240.5+6113\" rather than \"1FGL J0240.5+6113\"

Source names should not be duplicated because the output files are
constructed from the source names.

Optionally, an entry in the slist.dat file my have the form:
A_Source 40.1317, 61.2281 super_model.xml
where the 4th entry is the name of the model file to use for that source.
The source name (e.g. \"A_Source\") must be contained in the model file specified.
	

plist.dat
	This contains a list of photon files, one file per line.
It is also necessary to provide a single spacecraft (FT2) file.

Several temporary files are created while runing this script.
At least most are deleted at the end.

The output file for each source is an ASCII file called
XXX_like_lc.dat, where \"XXX\" is the source name as given in slist.dat
This file contains:
time (MJD) flux (cts/cm^2/s), flux error, time bin half width (days), TS value, spectral index, index error