Changeset 1725 in Sophya

Timestamp:

Oct 29, 2001, 11:42:20 PM (24 years ago)

Author:

aubourg

Message:

fits reader pour fichers sans samplenum, level s planck

Location:

trunk/ArchTOIPipe/Kernel

Files:

• fitstoirdr.cc (modified) (5 diffs)
• fitstoirdr.h (modified) (2 diffs)
• toiprocessor.cc (modified) (1 diff)

trunk/ArchTOIPipe/Kernel/fitstoirdr.cc

@@ -13 +13 @@
   fptr = NULL;
   totnscount = 0;
+
+  implicitSN = false;
+  implicitSNStart = 0;
 }
 
 FITSTOIReader::~FITSTOIReader() {
+}
+
+void FITSTOIReader::setImplicitSN(int snStart) {
+  implicitSN = true;
+  implicitSNStart = snStart;
 }
 
@@ -61 +69 @@
 
   cout << "FITSTOIReader cols = " << ncols << " rows=" << nrows << endl;
-  int anyNul;
-  double y;
-  fits_read_col_dbl(fptr,1,1,1,1,0,&y,&anyNul,&fstatus);
-  firstSn = (int) (y+.1);
+  if (implicitSN) {
+    firstSn = implicitSNStart;
+  } else {
+    int anyNul;
+    double y;
+    fits_read_col_dbl(fptr,1,1,1,1,0,&y,&anyNul,&fstatus);
+    firstSn = (int) (y+.1);
+  }
   fits_unlock();  
 }
@@ -72 +84 @@
 
   fits_lock();
-  // Dans cette version, on s'attend a ce que la premiere colonne soit le samplenum
+  // si pas implicitSN, la premiere colonne est le sampleNum. 
+  // Sinon, le samplenum est la fitsline + offset.
   int itoi=-1;
-  for (int i=1; i<ncols; i++) {
+  int col1 = implicitSN ? 0 : 1;
+  for (int i=col1; i<ncols; i++) {
     char templt[10];
     sprintf(templt, "%d", i+1);
@@ -122 +136 @@
 }
 
+// run 1 : deprecated. NON MAINTENU. Incompatible avec implicit SN.
 void FITSTOIReader::run1() {
   // Il faudrait optimiser en fonction de ce qui a ete demande comme samplenum,
@@ -206 +221 @@
      int anyNul;
      fits_lock();
-     fits_read_col_dbl(fptr,1,ideb+1,1,n,0,samplenum,&anyNul,&fstatus);
+     if (implicitSN) {
+       for (long j=0; j<Buff_Sz; j++) {
+         samplenum[j] = j+ideb+implicitSNStart;
+       }
+     } else {
+       fits_read_col_dbl(fptr,1,ideb+1,1,n,0,samplenum,&anyNul,&fstatus);
+     }
      for(int k=0; k<getNOut(); k++) {
        if(colval[k]==NULL)  continue;

trunk/ArchTOIPipe/Kernel/fitstoirdr.h

@@ -15 +15 @@
   FITSTOIReader(string fn,int buff_sz=1000);
   ~FITSTOIReader();
+
+  void         setImplicitSN(int snStart=0);
 
   virtual void addFile(string fn);
@@ -42 +44 @@
   int       snBegin; // first file
   int       snEnd;   // last file
+
+  bool      implicitSN;
+  int       implicitSNStart;
   
   vector<string> allfn;

trunk/ArchTOIPipe/Kernel/toiprocessor.cc

@@ -246 +246 @@
   for (int i=0; i<n; i++) {
     TOI* toi = outTOIs[i];
-    toi->putDone();
+    if (toi) {
+      toi->putDone();
+    }
   }
 }