00001
00002
00003
00004
00005
00006
00007
00008
00009
00010
00011
00012
00013
00014
00015
00016
00017
00018
00019
00020
00021
00022
00023
00024
00025
00026
00027
00028 #ifdef HAVE_CONFIG_H
00029 #include <config.h>
00030 #endif
00031
00032
00033
00034
00035
00036 #include "crires_recipe.h"
00037
00038
00039
00040
00041
00042 #define RECIPE_STRING "crires_util_genstd"
00043
00044
00045
00046
00047
00048 static int crires_util_genstd_save(cpl_table *, const cpl_parameterlist *,
00049 cpl_frameset *);
00050
00051 static char crires_util_genstd_description[] =
00052 "This recipe is used to generate standard star photospheric flux tables.\n"
00053 "The sof consists of file names tagged with "CRIRES_UTIL_GENSTD_RAW".\n"
00054 "The specified files are named after the standard star they represent\n"
00055 "(e.g. HIP61007.txt).\n"
00056 "The first line of the file must contain the RA and DEC (hh mm ss).\n"
00057 "(e.g. # 13 20 35.818 -36 42 44.26).\n"
00058 "The rest of the file must contain two columns:\n"
00059 "1st: Wavelengths in increasing order (the unit is corrected by\n"
00060 " the factor option to obtain nanometers).\n"
00061 "2nd: The atmospheric emission.\n"
00062 "The file is generated using the ASCII files in the catalogs/stdstar\n"
00063 "directory of the CRIRES source-code distribution."
00064 "\n"
00065 "This recipe produces 1 file for each input file:\n"
00066 "First product: the table with the photospheric flux of the std.\n"
00067 " (PRO TYPE = "CRIRES_PROTYPE_PHO_FLUX")\n" ;
00068
00069 CRIRES_RECIPE_DEFINE(crires_util_genstd,
00070 CRIRES_PARAM_PLOT,
00071 "Generate standard star FITS tables",
00072 crires_util_genstd_description) ;
00073
00074
00075
00076
00077
00078 static struct {
00079
00080 int display ;
00081
00082 } crires_util_genstd_config ;
00083
00084
00085
00086
00087
00088
00097
00098 static int crires_util_genstd(
00099 cpl_frameset * framelist,
00100 const cpl_parameterlist * parlist)
00101 {
00102 FILE * in ;
00103 char line[1024] ;
00104 cpl_bivector * bivec_ref ;
00105 cpl_bivector * bivec ;
00106 int ra1, ra2, dec1, dec2;
00107 double ra3, dec3;
00108 double ra, dec ;
00109 char isign;
00110 cpl_frame * cur_frame ;
00111 const char * cur_fname ;
00112 cpl_table * tab ;
00113 int nvals, nframes, nvals_ref ;
00114 double * pwave ;
00115 double * pemiss ;
00116 double * pwave_ref ;
00117 cpl_array * array ;
00118 int i, j ;
00119
00120
00121 crires_util_genstd_config.display = crires_parameterlist_get_bool(
00122 parlist, RECIPE_STRING, CRIRES_PARAM_PLOT) ;
00123
00124
00125 if (crires_dfs_set_groups(framelist, NULL)) {
00126 cpl_msg_error(__func__, "Cannot identify RAW and CALIB frames") ;
00127 return -1 ;
00128 }
00129 nframes = cpl_frameset_get_size(framelist) ;
00130
00131
00132 cur_frame = cpl_frameset_get_frame(framelist, 0) ;
00133 cur_fname = cpl_frame_get_filename(cur_frame) ;
00134 if ((bivec_ref = cpl_bivector_read(cur_fname))==NULL){
00135 cpl_msg_error(__func__, "Cannot load the file in the bivector") ;
00136 return -1 ;
00137 }
00138 pwave_ref = cpl_bivector_get_x_data(bivec_ref) ;
00139 nvals_ref = cpl_bivector_get_size(bivec_ref) ;
00140
00141
00142 tab = cpl_table_new(nframes+1) ;
00143 cpl_table_new_column(tab, CRIRES_COL_STDNAME, CPL_TYPE_STRING) ;
00144 cpl_table_new_column(tab, CRIRES_COL_RA, CPL_TYPE_DOUBLE) ;
00145 cpl_table_new_column(tab, CRIRES_COL_DEC, CPL_TYPE_DOUBLE) ;
00146 cpl_table_new_column_array(tab, CRIRES_COL_PHOTOFLUX, CPL_TYPE_DOUBLE,
00147 nvals_ref) ;
00148
00149
00150 cpl_table_set_string(tab, CRIRES_COL_STDNAME, 0, "WAVE") ;
00151 cpl_table_set_double(tab, CRIRES_COL_RA, 0, -1.0) ;
00152 cpl_table_set_double(tab, CRIRES_COL_DEC, 0, -1.0) ;
00153 array = cpl_array_wrap_double(pwave_ref, nvals_ref) ;
00154 cpl_table_set_array(tab, CRIRES_COL_PHOTOFLUX, 0, array) ;
00155 cpl_array_unwrap(array) ;
00156
00157
00158 for (i=0 ; i<nframes ; i++) {
00159
00160
00161 cur_frame = cpl_frameset_get_frame(framelist, i) ;
00162 cur_fname = cpl_frame_get_filename(cur_frame) ;
00163
00164
00165 if ((in = fopen(cur_fname, "r")) == NULL) {
00166 cpl_msg_error(__func__, "Could not open %s", cur_fname) ;
00167 cpl_table_delete(tab) ;
00168 cpl_bivector_delete(bivec_ref) ;
00169 return -1 ;
00170 }
00171 if (fgets(line, 1024, in) == NULL) {
00172 fclose(in) ;
00173 cpl_table_delete(tab) ;
00174 cpl_bivector_delete(bivec_ref) ;
00175 return -1 ;
00176 }
00177 if (sscanf(line, "#%d %d %lg %c%d %d %lg ", &ra1, &ra2, &ra3, &isign,
00178 &dec1, &dec2, &dec3) != 7) {
00179 cpl_msg_error(__func__, "Invalid first line in file %s", cur_fname);
00180 fclose(in) ;
00181 cpl_table_delete(tab) ;
00182 cpl_bivector_delete(bivec_ref) ;
00183 return -1 ;
00184 }
00185 fclose(in) ;
00186 ra = crires_ra_hms2deg(ra1, ra2, ra3) ;
00187 dec = crires_dec_hms2deg(dec1, dec2, dec3) ;
00188 if (isign == '-') dec *= -1.0 ;
00189
00190
00191 if ((bivec = cpl_bivector_read(cur_fname))==NULL){
00192 cpl_msg_error(__func__, "Cannot load the file in the bivector") ;
00193 cpl_bivector_delete(bivec_ref) ;
00194 cpl_table_delete(tab) ;
00195 return -1 ;
00196 }
00197 pwave = cpl_bivector_get_x_data(bivec) ;
00198 pemiss = cpl_bivector_get_y_data(bivec) ;
00199 nvals = cpl_bivector_get_size(bivec) ;
00200
00201
00202 if (nvals != nvals_ref) {
00203 cpl_msg_error(__func__, "Invalid file size: %s", cur_fname) ;
00204 cpl_bivector_delete(bivec_ref) ;
00205 cpl_bivector_delete(bivec) ;
00206 cpl_table_delete(tab) ;
00207 return -1 ;
00208 }
00209
00210
00211 for (j=0 ; j<nvals ; j++) {
00212 if (pwave[j] != pwave_ref[j]) {
00213 cpl_msg_error(__func__, "Invalid bins in %s", cur_fname) ;
00214 cpl_bivector_delete(bivec_ref) ;
00215 cpl_bivector_delete(bivec) ;
00216 cpl_table_delete(tab) ;
00217 return -1 ;
00218 }
00219 }
00220
00221
00222 if (crires_util_genstd_config.display) {
00223 irplib_bivector_plot(
00224 "set grid;set xlabel 'Wavelength (nm)';set ylabel 'Flux (jy)';",
00225 "t 'Photospheric flux' w lines", "", bivec);
00226 }
00227
00228
00229 cpl_table_set_string(tab, CRIRES_COL_STDNAME, i+1,
00230 crires_get_root_name(crires_get_base_name(cur_fname))) ;
00231
00232
00233 cpl_table_set_double(tab, CRIRES_COL_RA, i+1, ra) ;
00234 cpl_table_set_double(tab, CRIRES_COL_DEC, i+1, dec) ;
00235
00236
00237 array = cpl_array_wrap_double(pemiss, nvals) ;
00238 cpl_table_set_array(tab, CRIRES_COL_PHOTOFLUX, i+1, array) ;
00239 cpl_array_unwrap(array) ;
00240
00241 cpl_bivector_delete(bivec) ;
00242 }
00243 cpl_bivector_delete(bivec_ref) ;
00244
00245
00246 cpl_msg_info(__func__, "Save the table") ;
00247 if (crires_util_genstd_save(tab, parlist, framelist) == -1) {
00248 cpl_msg_error(__func__, "Cannot write the table") ;
00249 cpl_table_delete(tab) ;
00250 return -1 ;
00251 }
00252 cpl_table_delete(tab) ;
00253
00254 return 0 ;
00255 }
00256
00257
00265
00266 static int crires_util_genstd_save(
00267 cpl_table * out_table,
00268 const cpl_parameterlist * parlist,
00269 cpl_frameset * set)
00270 {
00271 cpl_propertylist * plist ;
00272
00273 plist = cpl_propertylist_new();
00274 cpl_propertylist_append_string(plist, "INSTRUME", "CRIRES") ;
00275 cpl_propertylist_append_string(plist, CPL_DFS_PRO_CATG,
00276 CRIRES_CALPRO_STD_PHOTOFLUX) ;
00277 cpl_propertylist_append_string(plist, CPL_DFS_PRO_TYPE,
00278 CRIRES_PROTYPE_PHO_FLUX) ;
00279
00280 if (cpl_dfs_save_table(set, NULL, parlist, set, NULL, out_table,
00281 NULL, "crires_util_genstd", plist, NULL,
00282 PACKAGE "/" PACKAGE_VERSION,
00283 "crires_util_genstd.fits") != CPL_ERROR_NONE) {
00284 cpl_msg_error(__func__, "Cannot save the table") ;
00285 return -1 ;
00286 }
00287 cpl_propertylist_delete(plist) ;
00288
00289
00290 return 0 ;
00291 }