Changeset 1725 in Sophya for trunk/ArchTOIPipe/Kernel/fitstoirdr.cc

Timestamp:

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

Author:

aubourg

Message:

fits reader pour fichers sans samplenum, level s planck

File:

• trunk/ArchTOIPipe/Kernel/fitstoirdr.cc (modified) (5 diffs)

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;