Changeset 3705 in Sophya for trunk/AddOn/TAcq/brproc.cc

Timestamp:

Dec 9, 2009, 6:17:05 PM (16 years ago)

Author:

ansari

Message:

Ajout possibilite verification timetag, Reza 9/12/2009

File:

• trunk/AddOn/TAcq/brproc.cc (modified) (5 diffs)

trunk/AddOn/TAcq/brproc.cc

@@ -75 +75 @@
   cout << " BRVisibilityCalculator/Info NMean= " << nmean_ << " JF1=" << jf1_ 
        << " JF2=" << jf2_ << " DJF=" << djf_ << endl;
+
+  fgcktt_=false;
 }
 
@@ -87 +89 @@
   poc << chanum_;
 
+  if (fgcktt_) {
+    cout << " BRVisibilityCalculator -  Check TimeTag done: " << endl;
+    for(size_t fib=0; fib<(size_t)memgr_.NbFibres(); fib++)  {
+      cout << " BRTTCheck-Fiber[" << fib << "] NBadTT=" << vbadtt_[fib] 
+           << "  NotSameTT=" << vnsamett_[fib] << endl; 
+    }
+    POutPersist pott(outpath_+"ttfcmtx.ppf");
+    pott << PPFNameTag("FC") << fcmtx_; 
+    pott << PPFNameTag("TT") << ttmtx_;     
+  }
   delete[] xnt_;
+}
+
+/* --Methode-- */
+int BRVisibilityCalculator::ActivateTimeTagCheck(uint_8 maxnpaq)
+{
+  cout << " BRVisibilityCalculator::ActivateTimeTagCheck() - TT/Fc matrix NCols=" << maxnpaq << endl;
+  mindeltatt_=memgr_.PaqSize()/2;
+  if (mindeltatt_<1) mindeltatt_=1;
+  fcmtx_.SetSize(memgr_.NbFibres(), maxnpaq);
+  ttmtx_.SetSize(memgr_.NbFibres(), maxnpaq);
+  vlasttt_.resize(memgr_.NbFibres(), 0);
+  vbadtt_.resize(memgr_.NbFibres(), 0);
+  vnsamett_.resize(memgr_.NbFibres(), 0);
+  fgcktt_=true;
+  return 0;
 }
 
@@ -134 +161 @@
   
   moyfc_ += curfc_[0];
-  moytt_ += vpaq_[0].TimeTag();
+  moytt_ += (vpaq_[0].TimeTag()-ttfirst_[0]);
+  if (fgcktt_)  CheckTimeTag();
   totnbpaq_++;
   return 0;
@@ -160 +188 @@
         r_4 moyreal=0.;
         r_4 moyimag=0.;
-        for(sa_size_t jjf=jf; jjf<jf+djf_; jjf++) {
+        sa_size_t jjfmx=jf+djf_;
+        if (jjfmx > vismtx_.NCols()) jjfmx=vismtx_.NCols();
+        for(sa_size_t jjf=jf; jjf<jjfmx; jjf++) {
           moyreal+=vismtx_(rv,jjf).real();
           moyimag+=vismtx_(rv,jjf).imag();
@@ -170 +200 @@
         visdt_.AddRow(xnt_);
       }
+    }
+  }
+  return 0;
+}
+
+/* --Methode-- */
+int BRVisibilityCalculator::CheckTimeTag()
+{
+  if (totnbpaq_==0) {
+    for(size_t fib=0; fib<(size_t)memgr_.NbFibres(); fib++)  {
+      vlasttt_[fib]=ttfirst_[fib];
+      if (ttmtx_.NCols()>0) {
+        fcmtx_(fib,totnbpaq_) = curfc_[fib]; 
+        ttmtx_(fib,totnbpaq_) = vlasttt_[fib]; 
+      }
+    }
+    return 0;
+  }
+  for(size_t fib=0; fib<(size_t)memgr_.NbFibres(); fib++) {
+    uint_8 ld = vpaq_[fib].TimeTag()-vlasttt_[fib];
+    uint_8 fd = vpaq_[fib].TimeTag()-ttfirst_[fib]-vpaq_[0].TimeTag()+ttfirst_[0];
+    if ( (ld < mindeltatt_) || (fd<-5) || (fd>5)) { vbadtt_[fib]++; vnsamett_[fib]++; }
+    else {
+      if (fd!=0)  vnsamett_[fib]++;
+    }
+    vlasttt_[fib]=vpaq_[fib].TimeTag();
+    if (totnbpaq_<ttmtx_.NCols()) {
+      fcmtx_(fib,totnbpaq_) = curfc_[fib]; 
+      ttmtx_(fib,totnbpaq_) = vlasttt_[fib]; 
     }
   }