next up previous contents index
Next: Operational Status and Monitoring Up: Manual Data Analysis Previous: Match Structures   Contents   Index


Relative Photometry

A match structure as described above is not yet finished. To minimize systematic errors in source intensities, it is helpful to perform a ``realtive photometry'' procedure to remove photometric mismeasurements due to such difficulties as the presence of thin haze over a portion of a frame, poor fringe subtraction (Section 7.3.3), etc. The idl program relphot3.pro will perform this procedure. For each object in each observation, relphot3.pro determines a small relative photometric correction, or in more extreme situations calculates a systematic error and sets flags to specify the problem and allow one to exclude the affected regions.

The program begins by selecting a set of good template sources in each matched list. To qualify as a template source, an object must be detected in more than 75% of the epochs at which the source location was imaged. Next, the median magnitude of all good observations of the template objects is calculated. Each image is then divided into subtiles, 200 pixels on a side, and a photometric offset are calculated for each template source, grouped by subtile, that has a statistical uncertainty in its magnitude measurement less than 0.1 mag. The number of template sources within a typical subtile is $\sim50$. The relative photometry map is then calculated to be the array of median offsets for each subtile. The RMS deviation of the offsets in each subtile, which for good subtiles is typically around 0.03 mag, is used to estimate the systematic error. In addition, a subtile is classified as ``bad'' if there are less than five template sources within it, or if the RMS deviation of the template offsets is greater than 0.1 mag. Each magnitude measurement of each object in each observation is then corrected by the offset map using a bilinear interpolation, and the new systematic error estimate is added in quadrature to the previously estimated systematic error for each object in the observation. If the subtile is flagged as bad, then the object's observations are likewise flagged (flags are defined in Table C.2).

The function call for this procedure is:

IDL$>$ relphot3, match, newmatch, rpmap, stat=stat, save=save, archive=archive, over=over,

init=init

The variable match contains the input match structure (or it may also be the file name of a FITS file containing a match structure), while newmatch is where idl will put the corrected match structure. If the match variable is a match structure, then the stats variable must be set to the array of cobj header structures (Section 8.3.2) that was created with that match structure. If the match variable is a file name, the array of structures in the second extension will be read into the stat variable. The variable rpmap is a structure that contains the map of offset values as well as several diagnostic parameters. This structure is defined in Section C.6, and it is saved as the third extension to the relmat FITS file (Section A.4).

The application of the save, archive, and over keywords are the same as with regmatch3_list.pro, above. If init is set, then the program will blindly apply the procedure. Otherwise, it will check to make sure the procedure has not already been applied to match, and exit with a warning if it seems that relative photometry is not necessary in this case.


next up previous contents index
Next: Operational Status and Monitoring Up: Manual Data Analysis Previous: Match Structures   Contents   Index
Rotse Pager 2003-05-20