diff options
Diffstat (limited to 'src/fuv/cf_flux_calibrate.c')
| -rw-r--r-- | src/fuv/cf_flux_calibrate.c | 173 |
1 files changed, 173 insertions, 0 deletions
diff --git a/src/fuv/cf_flux_calibrate.c b/src/fuv/cf_flux_calibrate.c new file mode 100644 index 0000000..2ad92e6 --- /dev/null +++ b/src/fuv/cf_flux_calibrate.c @@ -0,0 +1,173 @@ +/***************************************************************************** + * Johns Hopkins University + * Center For Astrophysical Sciences + * FUSE + ***************************************************************************** + * + * Synopsis: cf_flux_calibrate options intermediate_file + * + * Description: Flux calibrate each photon assigned to a channel. + * + * ** NEXT STEPS COMMENTED OUT FOR NOW ** + * For each photon in the target aperture, + * apply astigmatism correction to XFARF array and use the + * resulting coordinate system to apply flat-field and + * worm corrections to photon weights. + * + * By default all corrections are performed. Command line + * options allow one or more corrections to be omitted. + * + * Arguments: input_file Intermediate data file + * + * Calibration files: AEFF_CAL + * + * Returns: 0 on successful completion + * + * Calls: cf_effective_area, cf_astig_farf, cf_flat_field, + * cf_worm_correction + * + * History: 12/06/02 wvd 1.1 Begin work + * 02/12/03 wvd 1.2 Read ERGCM2S from IDF + * Move cf_convert_to_ergs to end of + * module + * 03/11/03 wvd 1.3 Change channel to type char + * 03/12/03 wvd 1.4 If ERGCM2S[i] is undefined, + * set channel[i] = 0 + * 03/19/03 peb 1.5 Added verbose_level option. + * 05/06/03 rdr 1.6 Read from column ERGCM2 rather + * than ERGCM2S + * 05/07/03 wvd 1.7 Change ergcm2s to ergcm2 + * 06/11/03 wvd 1.8 Pass datatype to cf_read_col and + * cf_write_col. + * Comment out call to cf_astig_farf(). + * 08/25/03 wvd 1.9 Change coltype from string to int in + * cf_read_col and cf_write_col. + * 09/17/03 wvd 1.10 Add some documentation. + * 09/17/03 wvd 1.11 Change channel to unsigned char. + * 05/20/05 wvd 1.12 Clean up i/o. + * + ****************************************************************************/ + +#include <unistd.h> +#include <stdlib.h> +#include "calfuse.h" + +static char CF_PRGM_ID[] = "cf_flux_calibrate"; +static char CF_VER_NUM[] = "1.12"; + +int +main(int argc, char *argv[]) +{ + unsigned char *channel=NULL; + int effective_area=1, astig_farf=1, flat_field=1, worm_correction=1; + int status=0, hdutype, optc; + long nevents, nseconds; + float *time=NULL, *ttime=NULL, *xfarf=NULL, *yfarf=NULL, *weight=NULL; + float *lambda=NULL, *ycentl=NULL, *ycents=NULL, *ergcm2=NULL; + fitsfile *header; + + char opts[] = "heafwv:"; + char usage[] = + "Usage:\n" + " cf_flux_calibrate [-heafw] [-v level] idf_file\n"; + char option[] = + "Options:\n" + " -h: this help message\n" + " -v: verbosity level (default is 1; 0 is silent)\n" + " -e: no effective-area calibration\n" + " -a: no astigmatism correction of XFARF array\n" + " -f: no flat-field correction\n" + " -w: no worm correction\n"; + + verbose_level = 1; + + while ((optc = getopt(argc, argv, opts)) != -1) { + switch(optc) { + + case 'h': + printf("%s\n%s", usage, option); + return 0; + case 'v': + verbose_level = atoi(optarg); + break; + case 'e': + effective_area = 0; + break; + case 'a': + astig_farf = 0; + break; + case 'f': + flat_field = 0; + break; + case 'w': + worm_correction = 0; + break; + } + } + cf_error_init(CF_PRGM_ID, CF_VER_NUM, stderr); + if (argc != optind+1) { + printf("%s", usage); + cf_if_error("Incorrect number of command-line arguments"); + } + cf_timestamp(CF_PRGM_ID, CF_VER_NUM, "Begin Processing"); + + FITS_open_file(&header, argv[optind], READWRITE, &status); + + FITS_movabs_hdu(header, 2, &hdutype, &status); + nevents = cf_read_col(header, TFLOAT, "TIME", (void **) &time); + nevents = cf_read_col(header, TFLOAT, "WEIGHT", (void **) &weight); + nevents = cf_read_col(header, TFLOAT, "XFARF", (void **) &xfarf); + nevents = cf_read_col(header, TFLOAT, "YFARF", (void **) &yfarf); + nevents = cf_read_col(header, TBYTE, "CHANNEL", (void **) &channel); + nevents = cf_read_col(header, TFLOAT, "LAMBDA", (void **) &lambda); + nevents = cf_read_col(header, TFLOAT, "ERGCM2", (void **) &ergcm2); + + FITS_movabs_hdu(header, 4, &hdutype, &status) ; + nseconds = cf_read_col(header, TFLOAT, "TIME", (void **) &ttime) ; + nseconds = cf_read_col(header, TFLOAT, "YCENT_LIF", (void **) &ycentl) ; + nseconds = cf_read_col(header, TFLOAT, "YCENT_SIC", (void **) &ycents) ; + + FITS_movabs_hdu(header, 1, &hdutype, &status); + + /*************************************************************************** + * Since we currently have no flat-field or worm correction, + * these subroutines are commented out. + * + if (astig_farf) + cf_astig_farf(header, nevents, xfarf, yfarf, channel, time, + nseconds, ttime, ycentl, ycents); + if (flat_field) + cf_flat_field(header, nevents, time, weight, xfarf, yfarf, channel, + nseconds, ttime, ycentl, ycents); + + if (worm_correction) + cf_worm_correction(header, nevents, time, weight, xfarf, yfarf, channel, + nseconds, ttime, ycentl, ycents); + ***************************************************************************/ + + if (effective_area) + cf_convert_to_ergs(header, nevents, weight, ergcm2, channel, lambda); + + cf_error_init(CF_PRGM_ID, CF_VER_NUM, stderr); + + FITS_movabs_hdu(header, 2, &hdutype, &status); + cf_write_col(header, TFLOAT, "WEIGHT", (void *) weight, nevents); + cf_write_col(header, TBYTE, "CHANNEL", (void *) channel, nevents); + cf_write_col(header, TFLOAT, "ERGCM2", (void *) ergcm2, nevents); + + FITS_close_file(header, &status); + + free(time); + free(weight); + free(xfarf); + free(yfarf); + free(channel); + free(lambda); + free(ergcm2); + free(ttime); + free(ycentl); + free(ycents); + + cf_timestamp(CF_PRGM_ID, CF_VER_NUM, "Done processing"); + return 0; +} |