Changeset 1483 in Sophya for trunk/ArchTOIPipe/TestPipes

Timestamp:

Apr 30, 2001, 10:52:01 AM (24 years ago)

Author:

ansari

Message:

Petites amelioration SimpleFourierFilter - Reza 30/4/2001

File:

• trunk/ArchTOIPipe/TestPipes/simtst.cc (modified) (10 diffs)

trunk/ArchTOIPipe/TestPipes/simtst.cc

@@ -34 +34 @@
     cout << "\n Usage : simtst [-dbg] [-start snb] [-end sne] \n" 
          << "         [-wtoi sz] [-wdegli sz] [-wfft sz] [-keepfft n] \n"
-         << "         inFitsName outFitsName  \n" 
+         << "         inFitsName outFitsName  outPPFName \n" 
          << "   -dbg : sets TOISeqBuffered debug level to 1 \n"
          << "   -start snb : sets the start sample num \n"
@@ -57 +57 @@
 
   bool fgdbg = false;
+  bool fgsetstart = false;
   int wtoi = 8192;
   int wdegli = 512;
@@ -62 +63 @@
   int keepfft = 0;
   int nmax = 10;
-  int istart = 104121000+wtoi*5;
+  int istart = 0;
   int iend = 0;
   double range_min = -16000;
@@ -68 +69 @@
   string infile;
   string outfile;
+  string outppfname;
   string intoi = "boloMuV_27";
   bool fg_f_filt = false;
@@ -78 +80 @@
       if (ia == narg-1) Usage(true);  // -start est suivi d'un argument 
       istart = atoi(arg[ia+1]); ia++;
+      fgsetstart = true;
     }    
     else if (strcmp(arg[ia],"-end") == 0) {
@@ -115 +118 @@
 
   if (iend < istart) iend = istart+wtoi*(nmax+5);
-  if ((narg-ko) < 2)  Usage(true); 
+  if ((narg-ko) < 3)  Usage(true); 
   infile = arg[ko];
   outfile = arg[ko+1];
+  outppfname = arg[ko+2];
 
   cout << " Initializing SOPHYA ... " << endl;
@@ -134 +138 @@
     //  mgr->setRequestedSample(11680920,11710584);
     //  mgr->setRequestedSample(104121000, 104946120);
-
-    mgr->setRequestedSample(istart, iend);
+    if (fgsetstart) 
+      mgr->setRequestedSample(istart, iend);
   
     //    FITSTOIReader r("/data/Archeops/bolo11.fits); 
@@ -146 +150 @@
 
     int w1 = wtoi;
-    int w2 = (fg_f_filt && (wfft > w1)) ? wfft : w1; 
+    int w2 = (fg_f_filt) ? 2*wfft+w1 : w1; 
 
     //    char * colname[5] = {"MJD", "UTC","boloMuV_11","magnFlux","pivot"};
@@ -158 +162 @@
     TOISeqBuffered * toisig = new TOISeqBuffered("sigma", w2);
     if (fgdbg) toisig->setDebugLevel(1);
-    TOISeqBuffered * toiincopie = new TOISeqBuffered("incopie", w1);
+    TOISeqBuffered * toiincopie = new TOISeqBuffered("incopie", w2);
     if (fgdbg) toiincopie->setDebugLevel(1);
     
@@ -234 +238 @@
     cout << " Creating Fourier Filter ... " << endl;
     SimpleFourierFilter sfft(fcoef);
-    sfft.KeepSpectra(keepfft);
+    sfft.KeepSpectra(outppfname, keepfft);
     TOISeqBuffered * toifft = NULL;
     TOISeqBuffered * toifiltcopie = NULL;