Changeset 1369 in Sophya for trunk/ArchTOIPipe

Timestamp:

Jan 3, 2001, 5:59:25 PM (25 years ago)

Author:

aubourg

Message:

flag handling

Location:

trunk/ArchTOIPipe/Kernel

Files:

• fitstoiwtr.cc (modified) (6 diffs)
• fitstoiwtr.h (modified) (2 diffs)

trunk/ArchTOIPipe/Kernel/fitstoiwtr.cc

@@ -10 +10 @@
   fname = fn;
   fstatus = 0;
+  outFlags = false;
 
   // Open file
@@ -20 +21 @@
 
 FITSTOIWriter::~FITSTOIWriter() {
+}
+
+void FITSTOIWriter::setOutFlags(bool yn) {
+  outFlags = yn;
 }
 
@@ -38 +43 @@
   }
 
-  int ncols = inIx.size();
+  int ndata = inIx.size();
+  int ncols = inIx.size() * (outFlags ? 2 : 1) + 1;
 
   char** colnames = new (char*[ncols]);
@@ -44 +50 @@
   char** colunits = new (char*[ncols]);
 
+  colnames[0] = "samplenum";
+  coltypes[0] = "1D";
+  colunits[0] = "integer";
+
   for (map<string, int>::iterator i = inIx.begin(); i != inIx.end(); i++) {
     int j = (*i).second;
     string n = (*i).first;
-    colnames[j] = const_cast<char*>(n.c_str());
-    coltypes[j] = "1D";
-    colunits[j] = "?";
+    if (outFlags) j += j;
+    colnames[j+1] = const_cast<char*>(n.c_str());
+    coltypes[j+1] = "1D";
+    colunits[j+1] = "?";
+    if (outFlags) {
+      colnames[j+2] = const_cast<char*>((n+"_flg").c_str());
+      coltypes[j+2] = "1D";
+      colunits[j+2] = "?";
+    }
   }
 
@@ -66 +82 @@
 
   int fitsLine = 1;
-  TOIManager* mgr = TOIManager::getManager();
-  // int snb = mgr->getRequestedBegin();
-  //int sne = mgr->getRequestedEnd();
   int snb = getMinIn();
   int sne = getMaxIn();
@@ -74 +87 @@
   for (int sn = snb; sn <= sne; sn++) {
     try {
-      for (int i=0; i<ncols; i++) {
-        //      cout << "wtr req " << sn << endl;
+      fits_lock();
+      double xx = sn;
+      fits_write_col_dbl(fptr, 1, fitsLine, 1, 1, &xx, &fstatus);
+      fits_unlock();
+      for (int i=0; i<ndata; i++) {
         double x = getData(i, sn);
-        //      cout << "wtr got " << sn << endl;
         fits_lock();
-        fits_write_col_dbl(fptr, i+1, fitsLine, 1, 1, &x, &fstatus);
+        if (outFlags) {
+          fits_write_col_dbl(fptr, 2*i+2, fitsLine, 1, 1, &x, &fstatus);
+          x = getFlag(i, sn);
+          fits_write_col_dbl(fptr, 2*i+3, fitsLine, 1, 1, &x, &fstatus);
+        } else {
+          fits_write_col_dbl(fptr, i+2, fitsLine, 1, 1, &x, &fstatus);
+        }
         if (fstatus != 0) {
           cerr << "fitstoiwtr error sn = " << sn << " i = "<< i << endl;

trunk/ArchTOIPipe/Kernel/fitstoiwtr.h

@@ -18 +18 @@
   virtual void  addInput(string name, TOI* toi);
 
+  virtual void  setOutFlags(bool yn = true);
+
   //virtual void  init();  
   virtual void  run();  
@@ -26 +28 @@
   int fstatus;
   vector<TOI*> fwinputs;
-
+  bool outFlags;
 };