initial FITS file copy code

knodlseder
1 parent 9204bcf5
Showing 3 changed files with 49 additions and 2 deletions
pyext/ctselect.i
src/ctselect/ctselect.cpp
src/ctselect/ctselect.hpp
@@ -30,10 +30,14 @@
 class ctselect : public GApplication  {
 public:
     // Constructors and destructors
-    ctselect(void) : GApplication(CTSELECT_NAME, CTSELECT_VERSION) { return; }
+    ctselect(void);
+    ctselect(int argc, char *argv[]);
+    ~ctselect(void);
     // Methods
-    int run(void);
+    void run(void);
+    void get_parameters(void);
+    void select(void);
 };
@@ -105,6 +105,9 @@ void ctselect::run(void)
         log << std::endl;
     }
+    // Select events
+    select();
+
     // Return
     return;
 }
@@ -134,6 +137,45 @@ void ctselect::get_parameters(void)
 }
+/***********************************************************************//**
+ * @brief Select events
+ *
+ * @todo An interesting alternative is the fits_select_rows routine which copies
+ * the selected rows from one FITS table to another. That's in fact the thing
+ * we want to achieve here. The only way to use this routine would be to
+ * open 2 FITS files, one for the input data and a fresh one for the output
+ * data. The events table would then be copied over using the select routine,
+ * any other tables would be copied over directly. The code would look like:
+ *
+ * GFits infile(m_infile);
+ * GFits outfile(m_outfile);
+ * for (int extno = 0; extno < infile.size(); ++extno) {
+ *     GFitsHDU* inhdu  = infile.hdu(i);
+ *     if (inhdu->extname() == "EVENTS") {
+ *         GFitsHDU* outhdu = inhdu->select("expression");
+ *         outfile.append_hdu(*outhdu);
+ *     }
+ *     else
+ *         outfile.append_hdu(*inhdu);
+ * }
+ * outfile.save();
+ * outfile.close();
+ *
+ ***************************************************************************/
+void ctselect::select(void)
+{
+    // Open FITS file
+    GFits file(m_infile);
+    //std::cout << file << std::endl;
+
+    // Save file
+    file.saveto(m_outfile, true);
+
+    // Return
+    return;
+}
+
+
 /*==========================================================================
  =                                                                         =
  =                             Private methods                             =
@@ -42,6 +42,7 @@ public:
     // Methods
     void run(void);
     void get_parameters(void);
+    void select(void);
 protected:
     // Protected methods