Changeset 3993 in Sophya for trunk/AddOn/TAcq/brprocGain.cc

Timestamp:

May 13, 2011, 12:37:14 AM (14 years ago)

Author:

ansari

Message:

Ajout possibilite de faire un DataTable (NTuple) avec la puissance spectrale dans plusieurs bandes en fonction du temps, Reza 13/05/2011

File:

• trunk/AddOn/TAcq/brprocGain.cc (modified) (2 diffs)

trunk/AddOn/TAcq/brprocGain.cc

@@ -185 +185 @@
     medfiltspecmtx_.Info()["NPAQSUM"] = nbtot_specwin_;
     medfiltspecmtx_.Info()[buff] = nbtot_specwin_;
-    //JEC timeStamp in the day @ millisec: faudra gerer le passage a minuit!
-    medfiltspecmtx_.Info()["TIMEWIN"] = (windowTime_.SecondsPart())*1000.;
-
-
     if ( nbtot_specwin_ > 0) {
       medfiltspecmtx_.Row(ir)  /= (r_4)nbtot_specwin_;
     }
   }//eo loop on channels
+  //JEC timeStamp in the day @ millisec: faudra gerer le passage a minuit!
+  //Reza : getObsTimeSeconds() fournit le temps ecoule en seconde a partir de 0h00 de la date du premier fichier
+  medfiltspecmtx_.Info()["TIMEWIN"] = (windowTime_.SecondsPart())*1000.;
+  if (nbpmoyttfc_>1)  {
+    moyfc_/=(double)nbpmoyttfc_;
+    moytt_/=(double)nbpmoyttfc_;
+  }
+  string ikey,ikdesc;    
+  ikey="MeanFC";  ikdesc="Mean FrameCounter";
+  medfiltspecmtx_.Info().SetI(ikey,moyfc_); 
+  medfiltspecmtx_.Info().SetComment(ikey,ikdesc);
+  ikey="MeanTT";  ikdesc="Mean TimeTag";
+  medfiltspecmtx_.Info().SetD(ikey,moytt_);
+  medfiltspecmtx_.Info().SetComment(ikey,ikdesc);
+
   char nfile[64];
   string flnm;
@@ -206 +217 @@
   //increment the file index
   nummedianfile_++;
+
+  if (dtpms_!=NULL)  FillPwrTmDTable(mspecmtx_);
     
   //reset the matirx
   medfiltspecmtx_ = (r_4)(0.);
-
   //reset counter
   nbtot_specwin_ = 0;
+  moyfc_=moytt_=0.;  
+  nbpmoyttfc_=0;
 
   return;