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");
147 cpl_propertylist_set_comment (idp_plist, qc_name,
"Exposure time");
151 if(!cpl_propertylist_has(
header,
"TEXPTIME"))
153 cpl_size nsets = cpl_table_get_nrow(oi_vis2_SC_allpol) /
GRAVI_NBASE ;
155 for (cpl_size set = 0 ; set < nsets ; set++)
159 texptime += max_int_time;
160 cpl_table_delete(int_time_this_set);
163 sprintf (qc_name,
"TEXPTIME");
165 cpl_propertylist_set_comment (idp_plist, qc_name,
"Total exposure time");
170 cpl_propertylist_get_double(
header,
"TEXPTIME") );
171 cpl_propertylist_set_comment (idp_plist,
"TEXPTIME",
"Total exposure time");
176 cpl_propertylist_update_string (idp_plist,
"PRODCATG",
"SCIENCE.VISIBILITY.UNCALIBRATED");
177 cpl_propertylist_set_comment (idp_plist,
"PRODCATG",
"Data product category");
180 double mjd_obs_first = DBL_MAX;
183 const cpl_frame *frame;
185 cpl_frameset_iterator *it = cpl_frameset_iterator_new(science_frames);
186 while ((frame = cpl_frameset_iterator_get(it)) != NULL) {
187 cpl_propertylist * this_frame_header = cpl_propertylist_load(cpl_frame_get_filename(frame), 0);
189 if (mjd_obs < mjd_obs_first)
190 mjd_obs_first = mjd_obs;
191 cpl_frameset_iterator_advance(it, 1);
192 cpl_propertylist_delete(this_frame_header);
194 cpl_frameset_delete(science_frames);
195 cpl_frameset_iterator_delete(it);
197 if (mjd_obs_first == DBL_MAX)
200 cpl_propertylist_set_comment (idp_plist,
"MJD-OBS",
"Start of observation");
203 double mjd_obs_last = 0;
204 double exptime_last = 0;
207 const cpl_frame *frame;
209 cpl_frameset_iterator *it = cpl_frameset_iterator_new(science_frames);
210 while ((frame = cpl_frameset_iterator_get(it)) != NULL) {
211 cpl_propertylist * this_frame_header = cpl_propertylist_load(cpl_frame_get_filename(frame), 0);
213 if (mjd_obs > mjd_obs_last)
215 mjd_obs_last = mjd_obs;
216 exptime_last = cpl_propertylist_get_double(this_frame_header,
"EXPTIME");
218 cpl_frameset_iterator_advance(it, 1);
219 cpl_propertylist_delete(this_frame_header);
221 cpl_frameset_delete(science_frames);
222 cpl_frameset_iterator_delete(it);
224 if (mjd_obs_last == 0)
226 if (exptime_last == 0)
228 if ( cpl_propertylist_has(
header,
"EXPTIME") )
229 exptime_last = cpl_propertylist_get_double(
header,
"EXPTIME");
231 exptime_last = cpl_propertylist_get_double(idp_plist,
"EXPTIME");
235 mjd_obs_last + exptime_last / 86400.);
236 cpl_propertylist_set_comment (idp_plist,
"MJD-END",
"End of observation");
239 cpl_propertylist_erase (
header,
"MJD-OBS");
242 cpl_propertylist_update_int (idp_plist,
"OBID1",
243 cpl_propertylist_get_int(
header,
"ESO OBS ID"));
244 cpl_propertylist_set_comment (idp_plist,
"OBID1",
"Obseration Block ID");
250 cpl_size nscience = cpl_frameset_get_size(science_frames);
251 cpl_frameset_delete(science_frames);
252 if(!cpl_propertylist_has(
header,
"NCOMBINE"))
256 cpl_propertylist_update_int (idp_plist,
"NCOMBINE", nscience);
257 cpl_propertylist_set_comment (idp_plist,
"NCOMBINE",
"Number of raw science combined");
262 cpl_propertylist_update_int (idp_plist,
"NCOMBINE",
263 cpl_propertylist_get_int(
header,
"NCOMBINE") );
264 cpl_propertylist_set_comment (idp_plist,
"NCOMBINE",
"Number of raw science combined");
271 if(cpl_propertylist_has(
header,
"ESO DPR TECH"))
273 cpl_propertylist_update_string (idp_plist,
"OBSTECH",
274 cpl_propertylist_get_string(
header,
"ESO DPR TECH") );
275 cpl_propertylist_set_comment (idp_plist,
"OBSTECH",
"Observation technique");
277 else if(cpl_propertylist_has(
header,
"OBSTECH"))
279 cpl_propertylist_update_string (idp_plist,
"OBSTECH",
280 cpl_propertylist_get_string(
header,
"OBSTECH") );
281 cpl_propertylist_set_comment (idp_plist,
"OBSTECH",
"Observation technique");
285 cpl_propertylist_update_string (idp_plist,
"SPECSYS",
"TOPOCENT");
286 cpl_propertylist_set_comment (idp_plist,
"SPECSYS",
"Frame of reference for spectral coordinates");
289 cpl_propertylist_update_string (idp_plist,
"TIMESYS",
"UTC");
290 cpl_propertylist_set_comment (idp_plist,
"TIMESYS",
"Time system");
297 if ( !strcmp (resolution,
"HIGH") )
299 if ( !strcmp (resolution,
"MED") )
301 if ( !strcmp (resolution,
"LOW") )
304 cpl_propertylist_set_comment (idp_plist,
"SPEC_ERR",
"Statistical error in spectral coordinate");
309 cpl_propertylist_set_comment (idp_plist,
"SPEC_SYE",
"Systematic error in spectral coordinate");
314 const cpl_frame *frame;
316 char prov_keyword[8];
317 cpl_frameset_iterator *it = cpl_frameset_iterator_new(frameset);
318 while ((frame = cpl_frameset_iterator_get(it)) != NULL) {
322 snprintf(prov_keyword, 7,
"PROV%zu",i_prov);
323 const char * filename = cpl_frame_get_filename(frame);
324 const char * filename_no_path = strrchr(filename,
'/');
325 if (filename_no_path == NULL)
326 filename_no_path = filename;
328 filename_no_path += 1;
329 cpl_propertylist_update_string(idp_plist, prov_keyword, filename_no_path);
332 cpl_frameset_iterator_advance(it, 1);
334 cpl_frameset_iterator_delete(it);
340 if(!cpl_table_has_column(oi_array,
"FOV"))
342 cpl_table_new_column(oi_array,
"FOV", CPL_TYPE_DOUBLE);
343 cpl_table_new_column(oi_array,
"FOVTYPE", CPL_TYPE_STRING);
346 if (telname == NULL) {
347 cpl_msg_warning(cpl_func,
"Cannot get TELNAME, FOV is not determined");
350 if (telname[0] ==
'U')
356 cpl_table_fill_column_window_double(oi_array,
"FOV", 0, cpl_table_get_nrow(oi_array), fov);
357 cpl_table_fill_column_window_string(oi_array,
"FOVTYPE", 0, cpl_table_get_nrow(oi_array),
"RADIUS");
362 cpl_table_delete(oi_vis2_SC_allpol);
363 cpl_table_delete(oi_T3_SC_allpol);
364 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_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)