diff options
Diffstat (limited to 'src/analysis/extract_jitter.c')
| -rw-r--r-- | src/analysis/extract_jitter.c | 198 |
1 files changed, 198 insertions, 0 deletions
diff --git a/src/analysis/extract_jitter.c b/src/analysis/extract_jitter.c new file mode 100644 index 0000000..a4ceb4d --- /dev/null +++ b/src/analysis/extract_jitter.c @@ -0,0 +1,198 @@ + +/************************************************************************* + * Johns Hopkins University + * Center For Astrophysical Sciences + * FUSE + ************************************************************************* + * + * + * Usage: + * extract_jitter [-h] [-v level] input_file + * + * + * + * Arguments: + * input_file : jitter file + * + * + * + * + * + * + * Options: + * -h: this help message + * -v: verbosity level (=1; 0 is silent) + * + * + * + * + * History: 11/04/2004 v1.0 bjg + * 06/03/2005 v1.1 wvd Fix bug that set expstart=expend. + * Delete unused variables. + * + ***************************************************************************/ + + +#include <stdlib.h> +#include <stdio.h> +#include <string.h> +#include <time.h> +#include <unistd.h> +#include <math.h> + +#include "calfuse.h" + + +static char CF_PRGM_ID[]= "extract_jitter"; +static char CF_VER_NUM[]= "1.0"; + + +int main(int argc,char *argv[]){ + + + + char input_filename[FLEN_CARD]; + char output_filename_x[FLEN_CARD]; + char output_filename_y[FLEN_CARD]; + + FILE * output_file_x; + FILE * output_file_y; + + fitsfile * infits; + + + + int intnull=0,anynull; + int status=0; + int hdutype=0; + + long i, N; + int ncol; + + double * hdu2_time; + double * hdu2_dx; + double * hdu2_dy; + int * hdu2_trkflg; + + + + double mjd,dx,dy; + + int trkflg; + + double expstart, expend; + + int optc; + + char opts[] = "hfv:"; + char usage[] = + "Usage:\n" + " extract_jitter [-h] [-v level] input_file\n\n" + "Arguments:\n" + " input_file : jitter file (FITS)\n\n"; + char option[] = + "Options:\n" + " -h: this help message \n" + " -v: verbosity level (=1; 0 is silent) \n"; + + verbose_level = 1; + + + /* Check number of options and arguments */ + while ((optc = getopt(argc, argv, opts)) != -1) { + switch(optc) { + case 'h': + printf("%s\n%s", usage, option); + return EXIT_SUCCESS; + case 'v': + verbose_level = atoi(optarg); + break; + } + } + + /* Initialize error checking. */ + cf_error_init(CF_PRGM_ID, CF_VER_NUM, stderr); + if (argc != optind+1) + cf_if_error("%s\nIncorrect number of program arguments", usage); + + + cf_timestamp(CF_PRGM_ID, CF_VER_NUM, "Started execution."); + + strcpy(input_filename,argv[optind]); + + + + strcpy(output_filename_x,input_filename); + strcpy(output_filename_y,input_filename); + strcat(output_filename_x,".jtx"); + strcat(output_filename_y,".jty"); + + + /* Create the output file */ + output_file_x = fopen(output_filename_x,"w"); + output_file_y = fopen(output_filename_y,"w"); + + + + FITS_open_file(&infits,input_filename,READONLY,&status); + + FITS_read_key(infits,TDOUBLE,"EXPSTART",&expstart,NULL,&status); + FITS_read_key(infits,TDOUBLE,"EXPEND",&expend,NULL,&status); + + FITS_movabs_hdu(infits,2,&hdutype,&status); + + FITS_read_key(infits,TLONG,"NAXIS2",&N,NULL,&status); + + hdu2_time=(double *)malloc(N*sizeof(double)); + hdu2_dx=(double *)malloc(N*sizeof(double)); + hdu2_dy=(double *)malloc(N*sizeof(double)); + hdu2_trkflg=(int *)malloc(N*sizeof(int)); + + FITS_get_colnum(infits, TRUE, "TIME", &ncol, &status); + FITS_read_col(infits, TDOUBLE, ncol, 1, 1, N, &intnull, + hdu2_time, &anynull, &status); + + FITS_get_colnum(infits, TRUE, "DX", &ncol, &status); + FITS_read_col(infits, TDOUBLE, ncol, 1, 1, N, &intnull, + hdu2_dx, &anynull, &status); + + FITS_get_colnum(infits, TRUE, "DY", &ncol, &status); + FITS_read_col(infits, TDOUBLE, ncol, 1, 1, N, &intnull, + hdu2_dy, &anynull, &status); + + FITS_get_colnum(infits, TRUE, "TRKFLG", &ncol, &status); + FITS_read_col(infits, TINT, ncol, 1, 1, N, &intnull, + hdu2_trkflg, &anynull, &status); + + + for (i=0;i<N-100;i++){ + mjd=expstart+(hdu2_time[i])/(3600.0*24.0); + dx=hdu2_dx[i]; + dy=hdu2_dy[i]; + trkflg=hdu2_trkflg[i]; + + if ((mjd>=expstart) && (mjd<=expend) && (trkflg==5)){ + fprintf(output_file_x,"%15lf %15lf\n",mjd,dx); + fprintf(output_file_y,"%15lf %15lf\n",mjd,dy); + } + } + + + + free(hdu2_time); + free(hdu2_dx); + free(hdu2_dy); + free(hdu2_trkflg); + + + + FITS_close_file(infits,&status); + + fclose(output_file_x); + fclose(output_file_y); + + + cf_timestamp(CF_PRGM_ID, CF_VER_NUM, "Finished execution."); + return EXIT_SUCCESS; + +} |