32 #include <fors_photometry_impl.h>
34 #include <fors_utils.h>
35 #include <test_simulate.h>
47 cpl_frameset_delete(frames); \
48 cpl_parameterlist_delete(parameters); \
49 cpl_propertylist_delete(header); \
57 cpl_parameterlist *parameters = cpl_parameterlist_new();
58 cpl_frameset *frames = cpl_frameset_new();
60 const char *filename[] = {
"photometry_aligned0.fits",
61 "photometry_aligned1.fits",
62 "photometry_aligned2.fits",
63 "photometry_aligned3.fits"};
64 cpl_propertylist *header = cpl_propertylist_new();
65 double exptime = 10.0;
66 cpl_frameset_insert(frames,
create_sky_flat(
"photometry_master_flat.fits",
73 CPL_FRAME_GROUP_CALIB));
77 cpl_propertylist_update_string(header,
"ORIGIN",
"ESO");
78 cpl_propertylist_update_double(header,
"MJD-OBS", 0.5);
82 for (i = 0; i < (int)(
sizeof(filename)/
sizeof(*filename)); i++)
85 airmass = 1.1 + 0.1*i;
88 cpl_propertylist_update_double(header,
"AIRMASS", airmass);
93 cpl_table *aligned = cpl_table_new(1);
95 cpl_frame *f = cpl_frame_new();
96 cpl_frame_set_tag(f, ALIGNED_PHOT);
97 cpl_frame_set_filename(f, filename[i]);
98 cpl_frameset_insert(frames, f);
100 cpl_table_new_column(aligned,
"INSTR_MAG", CPL_TYPE_DOUBLE);
101 cpl_table_new_column(aligned,
"DINSTR_MAG", CPL_TYPE_DOUBLE);
102 cpl_table_new_column(aligned,
"MAG", CPL_TYPE_DOUBLE);
103 cpl_table_new_column(aligned,
"DMAG", CPL_TYPE_DOUBLE);
104 cpl_table_new_column(aligned,
"CAT_MAG", CPL_TYPE_DOUBLE);
105 cpl_table_new_column(aligned,
"DCAT_MAG", CPL_TYPE_DOUBLE);
106 cpl_table_new_column(aligned,
"COLOR", CPL_TYPE_DOUBLE);
107 cpl_table_new_column(aligned,
"RA", CPL_TYPE_DOUBLE);
108 cpl_table_new_column(aligned,
"DEC", CPL_TYPE_DOUBLE);
109 cpl_table_new_column(aligned,
"X", CPL_TYPE_DOUBLE);
110 cpl_table_new_column(aligned,
"Y", CPL_TYPE_DOUBLE);
111 cpl_table_new_column(aligned,
"A", CPL_TYPE_DOUBLE);
112 cpl_table_new_column(aligned,
"B", CPL_TYPE_DOUBLE);
113 cpl_table_new_column(aligned,
"FWHM", CPL_TYPE_DOUBLE);
114 cpl_table_new_column(aligned,
"THETA", CPL_TYPE_DOUBLE);
115 cpl_table_new_column(aligned,
"CLASS_STAR", CPL_TYPE_DOUBLE);
116 cpl_table_new_column(aligned,
"OBJECT", CPL_TYPE_STRING);
117 cpl_table_new_column(aligned,
"USE_CAT", CPL_TYPE_INT);
119 cpl_table_set_double(aligned,
"INSTR_MAG", 0, -13.8 + (airmass-1)*0.1);
120 cpl_table_set_double(aligned,
"DINSTR_MAG", 0, 0.1);
121 cpl_table_set_double(aligned,
"MAG", 0, 14.4);
122 cpl_table_set_double(aligned,
"DMAG", 0, 0.004);
123 cpl_table_set_double(aligned,
"CAT_MAG", 0, 14.2);
124 cpl_table_set_double(aligned,
"DCAT_MAG", 0, 0.002);
125 cpl_table_set_double(aligned,
"COLOR", 0, 0.4);
127 cpl_table_set_double(aligned,
"RA", 0, 100);
128 cpl_table_set_double(aligned,
"DEC", 0, 1);
130 cpl_table_set_double(aligned,
"X", 0, 400);
131 cpl_table_set_double(aligned,
"Y", 0, 500);
132 cpl_table_set_double(aligned,
"A", 0, 2);
133 cpl_table_set_double(aligned,
"B", 0, 1);
134 cpl_table_set_double(aligned,
"FWHM", 0, 1.4);
135 cpl_table_set_double(aligned,
"THETA", 0, 0.1);
136 cpl_table_set_double(aligned,
"CLASS_STAR", 0, 0.1);
138 cpl_table_set_string(aligned,
"OBJECT", 0,
"MOBJEKT");
139 cpl_table_set_int(aligned,
"USE_CAT", 0, 1);
141 cpl_table_save(aligned, header, NULL, filename[i], CPL_IO_DEFAULT);
143 cpl_table_delete(aligned); aligned = NULL;
147 assure( !cpl_error_get_code(),
return, NULL );
149 fors_photometry_define_parameters(parameters);
150 cpl_parameter_set_int(cpl_parameterlist_find(parameters,
151 "fors.fors_photometry.degreef1"),
154 fors_photometry(frames, parameters);
static void test_photometry(void)
test
void create_standard_keys(cpl_propertylist *header, double exptime)
Write FORS standard keywords to simulated header.
int main(void)
Test of QC module.
cpl_frame * create_sky_flat(const char *filename, const char *tag, cpl_frame_group group, double exptime)
Simulate sky flat image.
cpl_frame * create_phot_table(const char *filename, const char *tag, cpl_frame_group group)
Create photometry table.