Photometry







 The package has a number of programs that can be used to make the photometry of a list of objects. Basically the photometry is done by fittng a PSF model to a fixed position in the subtracted image. This PSF model is generated by convolving with the kernel solution a model of the PSF in the reference image. The reference image can be divided in a number of areas (as for subtraction), and within each of these area a PSF model is constructed by median stacking of bright stars. The program which build this PSF map for the reference is called Bphot. The program which uses this PSF map and convolves it with the local kernel is called Cphot. This program calculates also the magnitude by fitting the profile, and write the data to the light curve. All these programs, as well as the image subtraction itself are managed by the shell script phot.csh. This script is in the sub-directory register of the distribution. Before you can run this script, you need to edit the configuration file of the photometric programs, which name is phot_config.
 

The most important parameter in phot_config the photometric radius, radphot. All the pixels within a distance radphot of the center of the object will be used to estimate the flux of the object by fitting the model of PSF. The other parameters do not need to be changed most of the time, or only changed slightly. Anyway be careful to setup the saturation parameter properly.
 

Before you run phot.csh, you must get accurate coordinates for the objects you want to make photometry. The program "stack" you already used to build the reference frame can do that for you. Just provide the program with a crude estimate of the coordinate in phot.data, and the program will return accurate coordinates. First you need to create a "temp.data" file which will contain the coordinates of the objects for which you want to have photometry.An example of such file format is:

 Xi    Yi
 250 100

 Xi and Yi are the coordinates of the maximum pixel of the object (integer coordinates). If you don't kow exactly what the coordinates of the maximum are, you can give approximate coordinates, and the program will find the nearest local maximum. You should then generate a few subtracted image in which your object has large flux variations from the reference, stack these images and get the position of your object. Note that
 

../bin/stack <image list> -o <ref name>  -i  temp.data

 The default output file is "phot.data", you need this file to do the photometry. Check phot.data, it should look like that:

   X            Y            Xi    Yi     File
 150.134 150.245 150 150  lc0.data

 X and Y are the accurate coordinates of the object and Xi and Yi are the coordinates of its maximum. File will be the name of the file containing the light curve of this object.

Do not forget to indicate the name of the reference file you want to use in the process_config file.
 
 

 Once you are done, you can run the photometry package, just type:

 ./phot.csh
 

 If it does not work, or does not do what you expect, first check that your file "process_config" is correct, and does contains the name of all the necessary files.