uves_physmod_stability_check.c

00001 /*                                                                              *
00002  *   This file is part of the ESO UVES Pipeline                                 *
00003  *   Copyright (C) 2004,2005 European Southern Observatory                      *
00004  *                                                                              *
00005  *   This library is free software; you can redistribute it and/or modify       *
00006  *   it under the terms of the GNU General Public License as published by       *
00007  *   the Free Software Foundation; either version 2 of the License, or          *
00008  *   (at your option) any later version.                                        *
00009  *                                                                              *
00010  *   This program is distributed in the hope that it will be useful,            *
00011  *   but WITHOUT ANY WARRANTY; without even the implied warranty of             *
00012  *   MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the              *
00013  *   GNU General Public License for more details.                               *
00014  *                                                                              *
00015  *   You should have received a copy of the GNU General Public License          *
00016  *   along with this program; if not, write to the Free Software                *
00017  *   Foundation, 51 Franklin St, Fifth Floor, Boston, MA  02111-1307  USA       *
00018  *                                                                              */
00019 
00020 /*
00021  * $Author: amodigli $
00022  * $Date: 2007/06/06 08:17:33 $
00023  * $Revision: 1.10 $
00024  * $Name: uves-3_3_1 $
00025  * $Log: uves_physmod_stability_check.c,v $
00026  * Revision 1.10  2007/06/06 08:17:33  amodigli
00027  * replace tab with 4 spaces
00028  *
00029  * Revision 1.9  2006/11/06 15:19:41  jmlarsen
00030  * Removed unused include directives
00031  *
00032  * Revision 1.8  2006/08/24 06:36:10  amodigli
00033  * fixed doxygen warning
00034  *
00035  * Revision 1.7  2006/08/22 14:17:00  amodigli
00036  * fixed bug on computing MedDY
00037  *
00038  * Revision 1.6  2006/08/08 15:41:19  amodigli
00039  * fixed bug in aligning actual and ref tables
00040  *
00041  * Revision 1.5  2006/08/02 07:02:05  amodigli
00042  * stability check more robust: size of actual and ref table may differ
00043  *
00044  * Revision 1.4  2006/07/31 06:29:26  amodigli
00045  * added QC on stability test
00046  *
00047  * Revision 1.3  2006/07/28 14:51:26  amodigli
00048  * fixed some bugs on improper table selection
00049  *
00050  * Revision 1.2  2006/06/28 13:28:29  amodigli
00051  * improved output
00052  *
00053  * Revision 1.1  2006/02/03 07:46:30  jmlarsen
00054  * Moved recipe implementations to ./uves directory
00055  *
00056  * Revision 1.4  2006/01/20 10:36:25  amodigli
00057  *
00058  * Fixed warings from doxigen
00059  *
00060  * Revision 1.3  2006/01/19 08:47:24  jmlarsen
00061  * Inserted missing doxygen end tag
00062  *
00063  * Revision 1.2  2006/01/16 13:52:58  jmlarsen
00064  * Removed memory leak
00065  *
00066  * Revision 1.1  2006/01/16 08:02:50  amodigli
00067  *
00068  * Added
00069  *
00070  * Revision 1.4  2006/01/09 14:05:42  amodigli
00071  * Fixed doxigen warnings
00072  *
00073  * Revision 1.3  2005/12/20 08:11:44  jmlarsen
00074  * Added CVS  entry
00075  *
00076  */
00077 
00078 /*----------------------------------------------------------------------------*/
00082 /*----------------------------------------------------------------------------*/
00084 #ifdef HAVE_CONFIG_H
00085 #  include <config.h>
00086 #endif
00087 
00088 
00089 /*-----------------------------------------------------------------------------
00090                                 Includes
00091  -----------------------------------------------------------------------------*/
00092 #include <uves_physmod_stability_check.h>
00093 #include <uves_utils_wrappers.h>
00094 #include <uves_msg.h>
00095 #include <uves_error.h>
00096 
00097 /*-----------------------------------------------------------------------------
00098                                 Defines
00099  -----------------------------------------------------------------------------*/
00100 /*-----------------------------------------------------------------------------
00101                             Functions prototypes
00102  ----------------------------------------------------------------------------*/
00103 /*-----------------------------------------------------------------------------
00104                             Static variables
00105  -----------------------------------------------------------------------------*/
00106 
00107 /*-----------------------------------------------------------------------------
00108                             Functions code
00109  -----------------------------------------------------------------------------*/
00110 
00111 /*----------------------------------------------------------------------------*/
00127 /*----------------------------------------------------------------------------*/
00128 
00129 int uves_physmod_stability_check(cpl_table* m_tbl, 
00130                                  cpl_table* r_tbl, 
00131                                  double* med_dx,
00132                                  double* med_dy,
00133                                  double* avg_dx,
00134                                  double* avg_dy)
00135 
00136 {
00137 
00138 
00139   cpl_table* tmp_tbl1=NULL;
00140   cpl_table* tmp_tbl2=NULL;
00141   cpl_table* tmp_tbl3=NULL;
00142   cpl_table* tmp_tbl4=NULL;
00143   cpl_table* tmp_tbl5=NULL;
00144   int msz=0;
00145   int rsz=0;
00146   int ref=0;
00147   int status=0;
00148   /* note those are swapped */
00149   check_nomsg(cpl_table_duplicate_column(m_tbl,"X",m_tbl,"XMES"));
00150   check_nomsg(cpl_table_duplicate_column(m_tbl,"Y",m_tbl,"YMES"));
00151   check_nomsg(cpl_table_duplicate_column(m_tbl,"S",m_tbl,"STATUS"));
00152   check_nomsg(cpl_table_duplicate_column(m_tbl,"O",m_tbl,"ORDER"));
00153   check_nomsg(cpl_table_duplicate_column(m_tbl,"ID",m_tbl,"IDENT"));
00154 
00155   check_nomsg(cpl_table_duplicate_column(r_tbl,"X",r_tbl,"XMES"));
00156   check_nomsg(cpl_table_duplicate_column(r_tbl,"Y",r_tbl,"YMES"));
00157   check_nomsg(cpl_table_duplicate_column(r_tbl,"S",r_tbl,"STATUS"));
00158   check_nomsg(cpl_table_duplicate_column(r_tbl,"O",r_tbl,"ORDER"));
00159   check_nomsg(cpl_table_duplicate_column(r_tbl,"ID",r_tbl,"IDENT"));
00160 
00161   check_nomsg(msz=cpl_table_get_nrow(m_tbl));
00162   check_nomsg(rsz=cpl_table_get_nrow(r_tbl));
00163 
00164   if(msz > rsz) {
00165     check_nomsg(cpl_table_set_size(r_tbl,msz));
00166     uves_msg_warning("Actual and reference tables has different sizes");
00167   } else if (msz < rsz) {
00168     check_nomsg(cpl_table_set_size(m_tbl,rsz));
00169     uves_msg_warning("Actual and reference tables has different sizes");
00170   } 
00171   /* copy values measured on actual frame into table from reference frame */
00172   check_nomsg(cpl_table_duplicate_column(r_tbl,"Xm",m_tbl,"X"));
00173   check_nomsg(cpl_table_duplicate_column(r_tbl,"Ym",m_tbl,"Y"));
00174   check_nomsg(cpl_table_duplicate_column(r_tbl,"Sm",m_tbl,"S"));
00175   check_nomsg(cpl_table_duplicate_column(r_tbl,"Om",m_tbl,"O"));
00176   check_nomsg(cpl_table_duplicate_column(r_tbl,"IDm",m_tbl,"ID"));
00177   check_nomsg(cpl_table_duplicate_column(r_tbl,"WAVEm",m_tbl,"WAVE"));
00178 
00179   /* initialize difference columns with values from actual frame measures */
00180   check_nomsg(cpl_table_duplicate_column(r_tbl,"DX",r_tbl,"Xm"));
00181   check_nomsg(cpl_table_duplicate_column(r_tbl,"DY",r_tbl,"Ym"));
00182   check_nomsg(cpl_table_duplicate_column(r_tbl,"DO",r_tbl,"Om"));
00183   check_nomsg(cpl_table_duplicate_column(r_tbl,"DID",r_tbl,"IDm"));
00184   check_nomsg(cpl_table_duplicate_column(r_tbl,"DW",r_tbl,"WAVEm"));
00185 
00186   /* subtract values measured on reference frame */
00187   check_nomsg(cpl_table_subtract_columns(r_tbl,"DX","X"));
00188   check_nomsg(cpl_table_subtract_columns(r_tbl,"DY","Y"));
00189   check_nomsg(cpl_table_subtract_columns(r_tbl,"DO","O"));
00190   check_nomsg(cpl_table_subtract_columns(r_tbl,"DID","ID"));
00191   check_nomsg(cpl_table_subtract_columns(r_tbl,"DW","WAVE"));
00192 
00193   /* WE GET THE ABSOLUTE VALUE: sqrt(X^2)*/
00194   check(cpl_table_power_column(r_tbl,"DW",2.),"Error computing power column");
00195   check(cpl_table_power_column(r_tbl,"DW",0.5),"Error computing power column");
00196 
00197   /* select proper values */
00198   check(tmp_tbl1=uves_extract_table_rows(r_tbl,"DO",CPL_EQUAL_TO,
00199           0),"Error selecting DO");
00200   check(tmp_tbl2=uves_extract_table_rows(tmp_tbl1,"DW",CPL_LESS_THAN,
00201           0.001),"Error selecting DW");
00202   check(tmp_tbl3=uves_extract_table_rows(tmp_tbl2,"DID",CPL_LESS_THAN,
00203           0.001),"Error selecting DID");
00204   check(tmp_tbl4=uves_extract_table_rows(tmp_tbl3,"S",CPL_EQUAL_TO,
00205           0),"Error selecting S");
00206   check(tmp_tbl5=uves_extract_table_rows(tmp_tbl4,"Sm",CPL_EQUAL_TO,
00207           0),"Error selecting Sm");
00208   check_nomsg(ref=cpl_table_get_nrow(tmp_tbl5)/2);
00209 
00210   check_nomsg(*med_dx=cpl_table_get_column_median(tmp_tbl5,"DX"));
00211   check_nomsg(*med_dy=cpl_table_get_column_median(tmp_tbl5,"DY"));
00212   check_nomsg(*avg_dx=cpl_table_get_column_mean(tmp_tbl5,"DX"));
00213   check_nomsg(*avg_dy=cpl_table_get_column_mean(tmp_tbl5,"DY"));
00214   /*
00215   uves_msg("Stability check results: Median DX = %5.3f Median DY = %5.3f",
00216            *med_dx,*med_dy);
00217        */       
00218   uves_msg("Stability check results: Mean DX = %5.3f Mean DY = %5.3f",
00219            *avg_dx,*avg_dy);
00220   check_nomsg(uves_sort_table_1(tmp_tbl5,"DX",0));
00221   check_nomsg(*med_dx=cpl_table_get_double(tmp_tbl5,"DX",ref,&status));
00222   check_nomsg(uves_sort_table_1(tmp_tbl5,"DY",0));
00223   check_nomsg(*med_dy=cpl_table_get_double(tmp_tbl5,"DY",ref,&status));
00224   /* median computed as in UVES-MIDAS */
00225   uves_msg("Stability check results: Median DX = %5.3f Median DY = %5.3f",
00226            *med_dx,*med_dy);
00227 
00228  cleanup:
00229   uves_free_table(&tmp_tbl1);
00230   uves_free_table(&tmp_tbl2);
00231   uves_free_table(&tmp_tbl3);
00232   uves_free_table(&tmp_tbl4);
00233   uves_free_table(&tmp_tbl5);
00234 
00235  return 0;
00236 
00237 
00238 }

Generated on Tue Jun 19 14:39:16 2007 for UVES Pipeline Reference Manual by  doxygen 1.4.6