Changeset 1454 in Sophya for trunk/ArchTOIPipe/TestPipes

Timestamp:

Apr 9, 2001, 2:14:21 PM (24 years ago)

Author:

ansari

Message:

Modifs effectuees a St Pierre de Charetreuse - Reza 9/4/2001

File:

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

trunk/ArchTOIPipe/TestPipes/simtst.cc

@@ -20 +20 @@
   }
   else {
-    cout << "\n Usage : tstrztoi [-sproc] [-rwtoi] [-dbg] [-start snb] [-end sne] \n" 
+    cout << "\n Usage : simtst [-sproc] [-rwtoi] [-dbg] [-start snb] [-end sne] \n" 
          << "         [-w1 sz] [-w2 sz] inFitsName outFitsName ppfFileName \n" 
          << "   -sproc : Run RzSimpleTOIProc (default RzTOIProc) \n"
@@ -27 +27 @@
          << "   -start snb : sets the start sample num \n"
          << "   -end sne : sets the end sample num \n"
+         << "   -range min,max : sets the acceptable range for intoi \n"
+         << "              default= -16000,16000\n"
+         << "   -intoi toiName : select input TOI name (def boloMuV_11)\n"
          << "   -w1 sz : sets main data window size (def= 8192)\n"
          << "   -w2 sz : sets secondary window size (def= 512) \n"
@@ -44 +47 @@
   bool fgdbg = false;
   int w1 = 8192;
-  int w2 = 512;
+  int w2 = 256;
   int nmax = 10;
   int istart = 104121000+w1*5;
   int iend = 0;
+  double range_min = -16000;
+  double range_max = 16000.;
   string infile;
   string outfile;
   string ppffile;
+  string intoi = "boloMuV_11";
 
   if (narg < 4) Usage(true);
@@ -72 +78 @@
       w2 = atoi(arg[ia+1]); ia++;
     }    
+    else if (strcmp(arg[ia],"-range") == 0) {
+      if (ia == narg-1) Usage(true);
+      sscanf(arg[ia+1],"%lf,%lf",&range_min, &range_max);
+      ia++;
+    }    
+    else if (strcmp(arg[ia],"-intoi") == 0) {
+      if (ia == narg-1) Usage(true);  
+      intoi = arg[ia+1]; ia++;
+    }    
     else if (strcmp(arg[ia],"-sproc") == 0)  fgsimple = true;
     else if (strcmp(arg[ia],"-rwtoi") == 0)  fgrwtoi = true;
@@ -93 +108 @@
   cout << ">>>> Window Size W1= " << w1 << " W2= " << w2
        << "  iStart= " << istart << " iEnd= " << iend << endl;
+  cout << ">>>> InTOIName= " << intoi << endl;
   try {
     TOIManager* mgr = TOIManager::getManager();
@@ -109 +125 @@
 
 
-    char * colname[5] = {"MJD", "UTC","boloMuV_11","magnFlux","pivot"};
-    char * toiname[5] = {"MJD", "UTC","bolo11","magneto","pivot"};
+    //    char * colname[5] = {"MJD", "UTC","boloMuV_11","magnFlux","pivot"};
+    //    char * toiname[5] = {"MJD", "UTC","bolo11","magneto","pivot"};
     TOISeqBuffered * toiin = new TOISeqBuffered("f2in", w1);
     if (fgdbg) toiin->setDebugLevel(1);
@@ -129 +145 @@
     //    }
 
-    int kk = 2;
-    r.addOutput(colname[kk], toiin);
+    //    int kk = 2;
+    r.addOutput(intoi, toiin);
     //  toi->dbg=true;
     //  r.addOutput("boloMuV_11", toi);
@@ -138 +154 @@
 
     cout << " Creating SimpleDeglitcher() " << endl;
-    SimpleDeglitcher degl;
+    SimpleDeglitcher degl(w2);
+    cout << " Setting Range for deglitcher: " << range_min 
+         << " - " << range_max << endl;
+    degl.SetRange(range_min, range_max);
     degl.addInput("in", toiin);
     degl.addOutput("out", toiout);
@@ -144 +163 @@
     degl.addOutput("sigma", toisig);
     degl.addOutput("incopie", toiincopie);
-    
-
+ 
+    /*   
+    cout << " Creating a GaussianFilter SimplFilterObject " << endl;
+    SimpleFilter filt(..)
+    TOISeqBuffered * toiincopie = new TOISeqBuffered("incopie", w1);
+    if (fgdbg) toiincopie->setDebugLevel(1);
+    */
 
     w.setOutFlags(true);
@@ -162 +186 @@
     degl.start();
     w.start();
-    
-    //    for(int jj=0; jj<3; jj++) {
-    //      sleep(2);
-    //      cout << *toi; 
-    //    }
+
+    /*    
+    for(int jj=0; jj<3; jj++) {
+      cout << *toiin; 
+      cout << *toiout; 
+      sleep(1);
+    }
+    */
 
     mgr->joinAll();