diff options
Diffstat (limited to 'src/libcf/cf_active_region.c')
| -rw-r--r-- | src/libcf/cf_active_region.c | 90 |
1 files changed, 90 insertions, 0 deletions
diff --git a/src/libcf/cf_active_region.c b/src/libcf/cf_active_region.c new file mode 100644 index 0000000..8a3fd6e --- /dev/null +++ b/src/libcf/cf_active_region.c @@ -0,0 +1,90 @@ +/***************************************************************************** + * Johns Hopkins University + * Center For Astrophysical Sciences + * FUSE + ***************************************************************************** + * + * Synopsis: cf_active_region(fitsfile *header, long nevents, float *xfarf, + * float *yfarf, unsigned char *locflags) + * + * Description: Flags events outside of detector active area. + * + * Note: Both cf_ttag_init and cf_hist_init set LOCATION_SHLD + * when initializing IDF, but they look only at XRAW and use + * limits that are less restrictive than those used here. + * + * Arguments: fitsfile *header Pointer to FITS file containing the + * header of the intermediate data file + * long nevents The number of events + * float *xfarf An array of event X positions + * float *yfarf An array of event Y positions + * unsigned char *locflags Location flags of each event + * + * Calls: + * + * Return: 0 on success + * + * History: 11/12/02 1.1 peb Begin and finish work + * 03/11/03 1.2 wvd Changed locflags to unsigned char + * 05/20/03 1.3 rdr Added call to cf_proc_check + * 07/29/03 1.4 wvd If cf_proc_check fails, return errflg. + * 04/07/07 1.5 wvd Clean up compiler warnings. + * + ****************************************************************************/ + +#include <string.h> +#include <stdlib.h> +#include "calfuse.h" + +int +cf_active_region(fitsfile *header, long nevents, float *xfarf, float *yfarf, + unsigned char *locflags) +{ + char CF_PRGM_ID[] = "cf_active_region"; + char CF_VER_NUM[] = "1.5"; + + char elecfile[FLEN_VALUE]={'\0'}; + char keyword[FLEN_KEYWORD]={'\0'}, detector[FLEN_VALUE]={'\0'}; + int errflg=0, status=0, active_l, active_r, active_b, active_t; + long j; + fitsfile *elecfits=NULL; + + cf_error_init(CF_PRGM_ID, CF_VER_NUM, stderr); + cf_timestamp(CF_PRGM_ID, CF_VER_NUM, "Begin Processing"); + + if ((errflg = cf_proc_check(header, CF_PRGM_ID))) return errflg; + + FITS_read_key(header, TSTRING, "DETECTOR", detector, NULL, &status); + /* + * Read keywords from ELEC_CAL file + */ + FITS_read_key(header, TSTRING, "ELEC_CAL", elecfile, NULL, &status); + FITS_open_file(&elecfits, cf_cal_file(elecfile), READONLY, &status); + sprintf(keyword, "ACTVL_%s", detector); + FITS_read_key(elecfits, TINT, keyword, &active_l, NULL, &status); + sprintf(keyword, "ACTVR_%s", detector); + FITS_read_key(elecfits, TINT, keyword, &active_r, NULL, &status); + sprintf(keyword, "ACTVB_%s", detector); + FITS_read_key(elecfits, TINT, keyword, &active_b, NULL, &status); + sprintf(keyword, "ACTVT_%s", detector); + FITS_read_key(elecfits, TINT, keyword, &active_t, NULL, &status); + FITS_close_file(elecfits, &status); + cf_verbose(3, "Active area limits: X from %d to %d, Y from %d to %d", + active_l, active_r, active_b, active_t); + /* + * Apply active region flags to each event. + */ + for(j=0; j<nevents; j++) { + /* + * Flag events that fall outside of detector active region. + */ + if (xfarf[j] < active_l || xfarf[j] > active_r || + yfarf[j] < active_b || yfarf[j] > active_t) { + locflags[j] |= LOCATION_SHLD; + } + } + + cf_proc_update(header, CF_PRGM_ID, "COMPLETE"); + cf_timestamp(CF_PRGM_ID, CF_VER_NUM, "Done processing"); + return 0; +} |