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

Timestamp:

May 17, 2010, 5:18:10 PM (15 years ago)

Author:

ansari

Message:

Ajout/modif classes (brproc.h,.cc) pour analyse de donnees BAORadio, Reza 17/05/2010

File:

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

trunk/AddOn/TAcq/brproc.cc

@@ -131 +131 @@
       if (stop_) break;
       if (memgr_.GetRunState() == MemZR_Stopped) break;
-      int mid = memgr_.FindMemZoneId(MemZA_ProcA);
+      int mid = memgr_.FindMemZoneId(mmact_);  // (MemZA_ProcA);
       Byte* buffg = memgr_.GetMemZone(mid);
       if (buffg == NULL) {
@@ -141 +141 @@
         fbuff_[fib] = memgr_.GetMemZone(mid,fib);
         if (fbuff_[fib] == NULL) { // cela ne devrait pas arriver 
-          cout << "BRBaseProcessor::run()/ERROR memgr.GetMemZone(" << mid << "," << fib << ") -> NULL" << endl;
+          cout << "BRVisibilityCalculator::run()/ERROR memgr.GetMemZone(" << mid << "," << fib << ") -> NULL" << endl;
           setRC(9);       fgrun=false;          
           break;
@@ -196 +196 @@
       if (rcpex!=0)  cout << " BRVisibilityCalculator::run() / Error Rc[paralex_.execute()]=" << rcpex << endl;
 
-      memgr_.FreeMemZone(mid, MemZS_ProcA);
+      memgr_.FreeMemZone(mid, mmsta_);  // (MemZS_ProcA);
     } // Fin de boucle sur les zones a traiter 
     //------------------------------------
@@ -207 +207 @@
   }
   catch (std::exception& exc) {
-    cout << " BRBaseProcessor::run()/catched std::exception " << exc.what() << endl;
+    cout << " BRVisibilityCalculator::run()/catched std::exception " << exc.what() << endl;
     setRC(98);  
     return; 
   }
   catch(...) {
-    cout << " BRBaseProcessor::run()/catched unknown ... exception " << endl;
+    cout << " BRVisibilityCalculator::run()/catched unknown ... exception " << endl;
     setRC(99);  
     return; 
@@ -397 +397 @@
   return 0;
 }
+
+
+//---------------------------------------------------------------------
+// Classe de traitement simple - calcul de spectres moyennes / voie
+//---------------------------------------------------------------------
+/* --Methode-- */
+BRMeanSpecCalculator::BRMeanSpecCalculator(RAcqMemZoneMgr& memgr, string outpath, uint_4 nmean)
+  : BRBaseProcessor(memgr), outpath_(outpath), nmean_(nmean)
+{
+  BRPaquet paq(memgr_.PaqSize());
+  mspecmtx_.SetSize(2*memgr_.NbFibres(), paq.DataSize()/4); 
+  numfile_=0;
+  totnbpaq_=0;
+}
+
+/* --Methode-- */
+BRMeanSpecCalculator::~BRMeanSpecCalculator()
+{
+}
+
+
+/* --Methode-- */
+int BRMeanSpecCalculator::Process()
+{
+      
+  if (totnbpaq_%nmean_ == 0) {
+    if (totnbpaq_ > 0) {
+      mspecmtx_.Info()["NPAQSUM"] = nmean_;
+      mspecmtx_ /= (double)nmean_;
+      char nfile[32];
+      sprintf(nfile,"mspecmtx%d.ppf",numfile_);
+      string flnm=outpath_+nfile;
+      POutPersist po(flnm);
+      po << mspecmtx_;
+      cout << numfile_ << "-BRMeanSpecCalculator::Process() NPaqProc="
+           << totnbpaq_ << "  -> Mean spectra Matrix in " << flnm << endl;
+      numfile_++;
+    }
+    mspecmtx_ = (r_8)(0.);
+  }
+
+  sa_size_t k=0;
+  for(size_t i=0; i<(size_t)2*memgr_.NbFibres(); i++) {
+    TwoByteComplex* zp=vpaq_[i/2].Data1C();
+    if (i%2==1)  zp=vpaq_[i/2].Data2C();
+    TVector< r_4 > spec = mspecmtx_.Row(k);   k++;
+    for(sa_size_t f=1; f<spec.Size(); f++) {
+      spec(f) += zp[f].module2F();
+      }
+    }
+
+  totnbpaq_++;
+  return 0;
+}
+
 
 //---------------------------------------------------------------