Changeset 1443 in Sophya for trunk/ArchTOIPipe/ProcWSophya

Timestamp:

Mar 15, 2001, 4:46:21 PM (25 years ago)

Author:

ansari

Message:

debug de Deglitcher, amelioration FITSTOIWriter - Reza 15/3/2001

Location:

trunk/ArchTOIPipe/ProcWSophya

Files:

• rztoi.cc (modified) (1 diff)
• simtoipr.cc (modified) (14 diffs)
• simtoipr.h (modified) (2 diffs)

trunk/ArchTOIPipe/ProcWSophya/rztoi.cc

@@ -212 +212 @@
     cout << " ... Processing remaining Samples - SN=" 
          << snlast << " ->" << sne << endl;
-    for(ii=snlast; ii<sne; ii++) {
+    for(ii=snlast; ii<=sne; ii++) {
       val = getData(0, ii);
       putData(0, ii, val*_fact);

trunk/ArchTOIPipe/ProcWSophya/simtoipr.cc

@@ -17 +17 @@
   maxpoints = maxnpt;
   totnscount = glnscount = glcount = 0;
+  deglitchdone = false;
 }
 
 SimpleDeglitcher::~SimpleDeglitcher()
 {
+}
+
+void SimpleDeglitcher::PrintStatus(ostream & os)
+{
+  os << "\n ------------------------------------------------------ \n" 
+     << " SimpleDeglitcher::PrintStatus() - WindowSize=" << WSize() 
+     << "  NbSigmas=" << NbSigmas() << " MaxPoints=" << MaxPoints() << endl;
+  TOIProcessor::PrintStatus(os);
+  if (deglitchdone) os << " Deglitching performed " << endl;
+  else os << " NO deglitching done " << endl;
+  double nst = (TotalSampleCount() > 0) ? TotalSampleCount() : 1.;
+  os << " TotalSampleCount=" << TotalSampleCount() << " GlitchCount= " 
+     << GlitchCount() << " GlitchSampleCount=" << GlitchSampleCount() 
+     << "( " << (double)GlitchSampleCount()*100./nst << " % )" << endl;
+  os << " ------------------------------------------------------ " << endl;
 }
 
@@ -29 +45 @@
   declareOutput("mean");
   declareOutput("sigma");
+  declareOutput("incopie");
   name = "SimpleDeglitcher";
   //  upExtra = 1;   A quoi ca sert ?
@@ -42 +59 @@
   bool fgmean = checkOutputTOIIndex(1);
   bool fgsigma = checkOutputTOIIndex(2);
+  bool fgincopie = checkOutputTOIIndex(3);
   
   if (!checkInputTOIIndex(0)) {
@@ -48 +66 @@
     throw ParmError("SimpleDeglitcher::run() Input TOI (in) not connected!");
   }
-  if (!fgout && !fgmean && !fgsigma) {
+  if (!fgout && !fgmean && !fgsigma &&!fgincopie) {
     cerr << " SimpleDeglitcher::run() - No Output TOI connected! "
          << endl;
@@ -59 +77 @@
   }
 
-  cout << " SimpleDeglitcher::run() SN=" << snb << " - " << sne << endl; 
+  cout << " SimpleDeglitcher::run() SNRange=" << snb << " - " << sne << endl; 
   try {
+    Timer tm("SimpleDeglitcher::run()");
     // Le debut 
     int wsz2 = wsize/2;
@@ -67 +86 @@
     for(k=0; k<wsize; k++) 
       vin(k) = getData(0, k+snb);
-    double s = vin.Sum()/wsize;
+    double s = vin.Sum();
     double mean = s/wsize;
     double s2 = vin.SumX2();
@@ -75 +94 @@
     bool fgglitch = false;
     double valcur,valsub,valadd;
-    for(k=0;k<sne-snb;k++) {
+    for(k=0;k<=sne-snb;k++) {
       totnscount++;
-      if (k%10000 == 0) cout << " DBG: K=" << k << endl;
+//       if (k%10000 == 0) cout << " DBG: K=" << k << endl;
       // Calcul mean-sigma
       if ((k>wsz2) && (k<sne-snb-wsz2)) {
         valsub = vin((k+wsz2)%wsize);
-        vin((k+wsz2)%wsize) = valadd = getData(0, k+snb);
+        vin((k+wsz2)%wsize) = valadd = getData(0, k+wsz2+snb);
         s += (valadd-valsub);
-        s2 += (valadd*valadd-valsub-valsub);
+        s2 += (valadd*valadd-valsub*valsub);
         mean = s/wsize;
         sigma = sqrt(s2/wsize-mean*mean);
@@ -94 +113 @@
         putData(2, k+snb, sigma, 0);
 
+      valcur = vin(k%wsize);
+      if (fgincopie) 
+        putData(3, k+snb, valcur, 0);
+        
       if (!fgout)   continue;  // Pas de sortie out (deglitche)
 
-      valcur = vin(k%wsize);
+      
+//             if (k<100) { 
+//              cout << "DBG-A-Deglitch[" << k << "] mean=" 
+//                   << mean << " sigma=" << sigma << " valcur=" 
+//                   << valcur << " Kgl=" << kgl ;
+//              if (fgglitch) cout << " In Glitch" ;
+//              cout << endl;
+//            }
 
       if (valcur < mean+nsig*sigma) {  // inferieur au seuil
@@ -139 +169 @@
       }
 
-      if (k%5000 == 0) cout << " ---DBG2: K=" << k << " glcount=" 
-                            << glcount << " LastPut= " << lastput << endl;
-    }  // Boucle sur les num-sample
-    cout << " La fin lastput=" << lastput << " SNE=" << sne; 
-    for(k=lastput-snb; k<sne-snb; k++) 
-      putData(0, k+snb, vin(k%wsize), 0);
-    cout << " DBG3- OUT of try bloc ! " << endl;
-      
+//           if (k%5000 == 0) cout << " ---DBG2: K=" << k << " glcount=" 
+//                              << glcount << " LastPut= " << lastput << endl;
+    }  // Fin de Boucle sur les num-sample
+
+    //DBG    cout << " La fin lastput=" << lastput << " SNE=" << sne; 
+    //DBG    for(k=lastput-snb+1; k<sne-snb; k++) 
+    //DBG      putData(0, k+snb, vin(k%wsize), 0);
+    //    cout << " DBG3- OUT of try bloc ! " << endl;
+    cout << " SimpleDeglitcher::run() - End of processing " 
+    << " TotalSampleCount=" << TotalSampleCount() 
+    <<  " GlitchCount= " << GlitchCount() << endl;
+    if (fgout) deglitchdone = true;
   }  // Bloc try 
   catch (PException & exc) {
@@ -174 +208 @@
 }
 
+void SimpleFilter::PrintStatus(ostream & os)
+{
+  os << "\n ------------------------------------------------------ \n" 
+     << " SimpleFilter::PrintStatus() - WindowSize=" << WSize() 
+     << " FilterKind= " << Type() << endl;
+  TOIProcessor::PrintStatus(os);
+  os << " Coeff= " ; 
+  for(int k=0; k<wsize; k++) os << coef[k] << " " ;
+  os << endl;
+  os << " TotalSampleCount=" << TotalSampleCount() << endl;
+  os << " ------------------------------------------------------ " << endl;
+}
+
 void SimpleFilter::init() {
   cout << "SimpleFilter::init" << endl;
   declareInput("in");
   declareOutput("out");
+  declareOutput("incopie");
   name = "SimpleFilter";
   //  upExtra = 1;
@@ -187 +235 @@
   int sne = getMaxIn();
 
-  cout << " SimpleFilter::run(): snb=" << snb << " sne=" << sne << endl;
-
   bool fgout = checkOutputTOIIndex(0);
+  bool fgincopie = checkOutputTOIIndex(1);
 
   if (!checkInputTOIIndex(0)) {
@@ -201 +248 @@
     throw ParmError("SimpleFilter::run() No output TOI connected!");
   }
+
+  cout << " SimpleFilter::run() SNRange=" << snb << " - " << sne << endl; 
   
 
   try {
+    Timer tm("SimpleFilter::run()");
     // Le debut
     int wsz2 = wsize/2;
@@ -212 +262 @@
     double mean = vin.Sum()/wsize;
     for(k=0; k<wsize/2; k++) vin(k) = mean; 
-    for(k=0;k<sne-snb;k++) {
+    for(k=0;k<=sne-snb;k++) {
       totnscount++;
       // Calcul mean-sigma
@@ -224 +274 @@
       }
       putData(0,k+snb,sortie,0);
+      if (fgincopie) putData(1,k+snb,vin(k%wsize),0);
     }  // Boucle sur les num-sample
+    cout << " SimpleFilter::run() - End of processing " << endl;
   }  // Bloc try 
 

trunk/ArchTOIPipe/ProcWSophya/simtoipr.h

@@ -27 +27 @@
   inline int    GlitchSampleCount() const { return glnscount; }
   
+  virtual void  PrintStatus(ostream & os) ; // const plus tard
+  
 protected:
   int totnscount;   // Nombre total d'echantillon processe
   int glnscount;    // Nombre total de glitch
   int glcount;      // Nombre de glitch detecte
+  bool deglitchdone;  // Deglitch effectue
 
   int wsize;        // Taille de fenetre de travail
@@ -62 +65 @@
   inline int    TotalSampleCount() const { return totnscount; }
 
+  virtual void  PrintStatus(ostream & os) ; // const plus tard
+
   virtual void  init();  
   virtual void  run();