49 cpl_frameset * frameset)
51 cpl_propertylist * idp_plist = cpl_propertylist_new ();
58 cpl_table * oi_vis2_SC_allpol = NULL;
59 cpl_table * oi_T3_SC_allpol = NULL;
60 cpl_table * oi_wave_SC_allpol = NULL;
71 for (
int pol = 1; pol < npol_sc; pol++) {
77 for (
int pol = 0; pol < npol_sc; pol++) {
79 vis2err += this_pol_vis2err * this_pol_vis2err;
81 visphierr += this_pol_visphierr * this_pol_visphierr;
83 t3phierr += this_pol_t3phierr * this_pol_t3phierr;
87 sprintf (qc_name,
"VIS2ERR");
89 cpl_propertylist_set_comment (idp_plist, qc_name,
"Representative squared visibility error [%]");
91 sprintf (qc_name,
"VISPHERR");
93 cpl_propertylist_set_comment (idp_plist, qc_name,
"Representative visibility phase error [%]");
95 sprintf (qc_name,
"T3PHIERR");
97 cpl_propertylist_set_comment (idp_plist, qc_name,
"Representative closure phase error [deg]");
99 double min_uvcoord, max_uvcoord;
102 double min_eff_wave = cpl_table_get_column_min(oi_wave_SC_allpol,
"EFF_WAVE");
103 double max_eff_wave = cpl_table_get_column_max(oi_wave_SC_allpol,
"EFF_WAVE");
105 double base_max = max_uvcoord / min_eff_wave;
106 double base_min = min_uvcoord / max_eff_wave;
107 if(isnan(base_max) || isinf(base_max))
109 if(isnan(base_min) || isinf(base_min))
112 sprintf (qc_name,
"BASE_MAX");
114 cpl_propertylist_set_comment (idp_plist, qc_name,
"Maximum baseline / Minimum effective wavelenth");
116 sprintf (qc_name,
"BASE_MIN");
118 cpl_propertylist_set_comment (idp_plist, qc_name,
"Minimum baseline / Maximum effective wavelenth");
122 sprintf (qc_name,
"NUM_CHAN");
123 cpl_propertylist_update_int (idp_plist, qc_name, cpl_table_get_nrow(oi_wave_SC_allpol) / npol_sc );
124 cpl_propertylist_set_comment (idp_plist, qc_name,
"Number of wavelength channels");
126 sprintf (qc_name,
"WAVELMAX");
127 double avg_band = cpl_table_get_column_mean(oi_wave_SC_allpol,
"EFF_BAND");
129 cpl_propertylist_set_comment (idp_plist, qc_name,
"[nm] Maximum wavelength");
131 sprintf (qc_name,
"WAVELMIN");
133 cpl_propertylist_set_comment (idp_plist, qc_name,
"[nm] Minimum wavelength");
135 cpl_table_duplicate_column(oi_wave_SC_allpol,
"SPEC_RES", oi_wave_SC_allpol,
"EFF_WAVE");
136 cpl_table_divide_columns(oi_wave_SC_allpol,
"SPEC_RES",
"EFF_BAND");
138 sprintf (qc_name,
"SPEC_RES");
140 cpl_propertylist_set_comment (idp_plist, qc_name,
"Spectral resolution");
144 sprintf (qc_name,
"EXPTIME");
145 double mean_int_time = 0.0;
146 cpl_size valid_elem = 0;
147 cpl_size nrow = cpl_table_get_nrow (oi_vis2_SC_allpol);
148 double * inttime = cpl_table_get_data_double (oi_vis2_SC_allpol,
"INT_TIME");
149 for (cpl_size r=0; r<nrow;r++)
153 mean_int_time += inttime[r];
159 mean_int_time = mean_int_time / valid_elem;
161 cpl_propertylist_set_comment (idp_plist, qc_name,
"Exposure time");
165 if(!cpl_propertylist_has(
header,
"TEXPTIME"))
168 for (
int pol = 0; pol < npol_sc; pol++) {
170 double this_pol_texptime = 0;
171 cpl_size nsets = cpl_table_get_nrow(this_pol_table) /
GRAVI_NBASE ;
172 for (cpl_size set = 0 ; set < nsets ; set++)
176 this_pol_texptime += max_int_time;
177 cpl_table_delete(int_time_this_set);
179 texptime = fmax(texptime, this_pol_texptime);
182 sprintf (qc_name,
"TEXPTIME");
184 cpl_propertylist_set_comment (idp_plist, qc_name,
"Total exposure time");
189 cpl_propertylist_get_double(
header,
"TEXPTIME") );
190 cpl_propertylist_set_comment (idp_plist,
"TEXPTIME",
"Total exposure time");
195 cpl_propertylist_update_string (idp_plist,
"PRODCATG",
"SCIENCE.VISIBILITY.UNCALIBRATED");
196 cpl_propertylist_set_comment (idp_plist,
"PRODCATG",
"Data product category");
199 double mjd_obs_first = DBL_MAX;
202 const cpl_frame *frame;
204 cpl_frameset_iterator *it = cpl_frameset_iterator_new(science_frames);
205 while ((frame = cpl_frameset_iterator_get(it)) != NULL) {
206 cpl_propertylist * this_frame_header = cpl_propertylist_load(cpl_frame_get_filename(frame), 0);
208 if (mjd_obs < mjd_obs_first)
209 mjd_obs_first = mjd_obs;
210 cpl_frameset_iterator_advance(it, 1);
211 cpl_propertylist_delete(this_frame_header);
213 cpl_frameset_delete(science_frames);
214 cpl_frameset_iterator_delete(it);
216 if (mjd_obs_first == DBL_MAX)
219 cpl_propertylist_set_comment (idp_plist,
"MJD-OBS",
"Start of observation");
222 double mjd_obs_last = 0;
223 double exptime_last = 0;
226 const cpl_frame *frame;
228 cpl_frameset_iterator *it = cpl_frameset_iterator_new(science_frames);
229 while ((frame = cpl_frameset_iterator_get(it)) != NULL) {
230 cpl_propertylist * this_frame_header = cpl_propertylist_load(cpl_frame_get_filename(frame), 0);
232 if (mjd_obs > mjd_obs_last)
234 mjd_obs_last = mjd_obs;
235 exptime_last = cpl_propertylist_get_double(this_frame_header,
"EXPTIME");
237 cpl_frameset_iterator_advance(it, 1);
238 cpl_propertylist_delete(this_frame_header);
240 cpl_frameset_delete(science_frames);
241 cpl_frameset_iterator_delete(it);
243 if (mjd_obs_last == 0)
245 if (exptime_last == 0)
247 if ( cpl_propertylist_has(
header,
"EXPTIME") )
248 exptime_last = cpl_propertylist_get_double(
header,
"EXPTIME");
250 exptime_last = cpl_propertylist_get_double(idp_plist,
"EXPTIME");
254 mjd_obs_last + exptime_last / 86400.);
255 cpl_propertylist_set_comment (idp_plist,
"MJD-END",
"End of observation");
258 cpl_propertylist_erase (
header,
"MJD-OBS");
261 cpl_propertylist_update_int (idp_plist,
"OBID1",
262 cpl_propertylist_get_int(
header,
"ESO OBS ID"));
263 cpl_propertylist_set_comment (idp_plist,
"OBID1",
"Obseration Block ID");
269 cpl_size nscience = cpl_frameset_get_size(science_frames);
270 cpl_frameset_delete(science_frames);
271 if(!cpl_propertylist_has(
header,
"NCOMBINE"))
275 cpl_propertylist_update_int (idp_plist,
"NCOMBINE", nscience);
276 cpl_propertylist_set_comment (idp_plist,
"NCOMBINE",
"Number of raw science combined");
281 cpl_propertylist_update_int (idp_plist,
"NCOMBINE",
282 cpl_propertylist_get_int(
header,
"NCOMBINE") );
283 cpl_propertylist_set_comment (idp_plist,
"NCOMBINE",
"Number of raw science combined");
290 if(cpl_propertylist_has(
header,
"ESO DPR TECH"))
292 cpl_propertylist_update_string (idp_plist,
"OBSTECH",
293 cpl_propertylist_get_string(
header,
"ESO DPR TECH") );
294 cpl_propertylist_set_comment (idp_plist,
"OBSTECH",
"Observation technique");
296 else if(cpl_propertylist_has(
header,
"OBSTECH"))
298 cpl_propertylist_update_string (idp_plist,
"OBSTECH",
299 cpl_propertylist_get_string(
header,
"OBSTECH") );
300 cpl_propertylist_set_comment (idp_plist,
"OBSTECH",
"Observation technique");
304 cpl_propertylist_update_string (idp_plist,
"SPECSYS",
"TOPOCENT");
305 cpl_propertylist_set_comment (idp_plist,
"SPECSYS",
"Frame of reference for spectral coordinates");
308 cpl_propertylist_update_string (idp_plist,
"TIMESYS",
"UTC");
309 cpl_propertylist_set_comment (idp_plist,
"TIMESYS",
"Time system");
316 if ( !strcmp (resolution,
"HIGH") )
318 if ( !strcmp (resolution,
"MED") )
320 if ( !strcmp (resolution,
"LOW") )
323 cpl_propertylist_set_comment (idp_plist,
"SPEC_ERR",
"Statistical error in spectral coordinate");
328 cpl_propertylist_set_comment (idp_plist,
"SPEC_SYE",
"Systematic error in spectral coordinate");
333 const cpl_frame *frame;
335 char prov_keyword[8];
336 cpl_frameset_iterator *it = cpl_frameset_iterator_new(frameset);
337 while ((frame = cpl_frameset_iterator_get(it)) != NULL) {
343 snprintf(prov_keyword, 7,
"PROV%zu",i_prov);
344 const char * filename = cpl_frame_get_filename(frame);
345 const char * filename_no_path = strrchr(filename,
'/');
346 if (filename_no_path == NULL)
347 filename_no_path = filename;
349 filename_no_path += 1;
350 cpl_propertylist_update_string(idp_plist, prov_keyword, filename_no_path);
353 cpl_frameset_iterator_advance(it, 1);
355 cpl_frameset_iterator_delete(it);
361 if(!cpl_table_has_column(oi_array,
"FOV"))
363 cpl_table_new_column(oi_array,
"FOV", CPL_TYPE_DOUBLE);
364 cpl_table_new_column(oi_array,
"FOVTYPE", CPL_TYPE_STRING);
367 if (telname == NULL) {
368 cpl_msg_warning(cpl_func,
"Cannot get TELNAME, FOV is not determined");
371 if (telname[0] ==
'U')
377 cpl_table_fill_column_window_double(oi_array,
"FOV", 0, cpl_table_get_nrow(oi_array), fov);
378 cpl_table_fill_column_window_string(oi_array,
"FOVTYPE", 0, cpl_table_get_nrow(oi_array),
"RADIUS");
383 cpl_table_delete(oi_vis2_SC_allpol);
384 cpl_table_delete(oi_T3_SC_allpol);
385 cpl_table_delete(oi_wave_SC_allpol);
typedefCPL_BEGIN_DECLS struct _gravi_data_ gravi_data
#define gravi_data_get_oi_t3(data, type, pol, npol)
#define gravi_data_get_oi_vis2(data, type, pol, npol)
#define gravi_data_get_oi_wave(data, type, pol, npol)
#define gravi_data_get_oi_vis(data, type, pol, npol)
#define GRAVI_SINGLE_SCIENCE_RAW
#define GRAVI_DUAL_CALIB_RAW
#define GRAVI_SINGLE_CALIB_RAW
#define GRAVI_DUAL_SCIENCE_RAW
cpl_propertylist * gravi_idp_compute(gravi_data *vis_data, cpl_propertylist *header, cpl_frameset *frameset)
Create IDP keywords to satisfy standard.
cpl_propertylist * header
cpl_propertylist_update_double(header, "ESO QC MINWAVE SC", cpl_propertylist_get_double(plist, "ESO QC MINWAVE SC"))
#define GRAVI_OI_VIS2_EXT
#define GRAVI_OI_ARRAY_EXT
double gravi_table_get_column_flagged_mean(cpl_table *table, const char *name)
Function to compute the mean of a column table with arrays.
double gravi_table_get_column_flagged_max(cpl_table *table, const char *name)
Function to compute the maximum of a column table with arrays.
int gravi_data_has_extension(gravi_data *raw_calib, const char *ext_name)
Check if data has extension with given EXTNAME.
cpl_table * gravi_data_get_table(gravi_data *self, const char *extname)
Return a pointer on a table extension by its EXTNAME.
cpl_frameset * gravi_frameset_extract_fringe_data(cpl_frameset *frameset)
int gravi_pfits_get_pola_num(const cpl_propertylist *plist, int type_data)
double gravi_pfits_get_mjd(const cpl_propertylist *plist)
const char * gravi_pfits_get_resolution(const cpl_propertylist *plist)
const char * gravi_conf_get_telname(int gravi_beam, cpl_propertylist *header)
cpl_error_code gravi_data_get_minmax_uvcoord(const cpl_table *oi_vis2, double *min_uvcoord, double *max_uvcoord)
Compute the minimum and maximum values of sqrt(ucoord**2 + vcoord**2)