Changeset 3905 in Sophya for trunk/AddOn/TAcq/brproc.h

Timestamp:

Oct 15, 2010, 11:02:36 PM (15 years ago)

Author:

ansari

Message:

implementation de la fenetre en temps pour les spectres ds BRMeanSpecCalculator, Reza 15/10/2010

File:

• trunk/AddOn/TAcq/brproc.h (modified) (4 diffs)

trunk/AddOn/TAcq/brproc.h

@@ -6 +6 @@
 // Classes de threads de traitememt donnees BAORadio 
 //  Principales classes utilisables 
-//  - BRMeanSpecCalculator (donnees FFT)
+//  - BRMeanSpecCalculator (donnees FFT ou RAW apres BRFFTCalculator)
 //  - BRFFTCalculator (donnees RAW) 
 //----------------------------------------------------------------
@@ -27 +27 @@
 using namespace std;
 
-
  
 //---------------------------------------------------------------------
@@ -40 +39 @@
   virtual ~BRMeanSpecCalculator(); // Fait le nettoyage final , sauvegarde des objets ...
 
+// Definit la taille de la fenetre temporelle (winsz = nb de paquets) des spectres a traiter
+// Cette zone est etendu de wszext paquet en debut et en fin de paquet 
+  void SetSpectraWindowSize(uint_4 winsz=15, uint_4 wszext=3);
+// Taille de la fenetre en temps (nb de paquets a traiter en bloc
+  inline sa_size_t GetSpectraWindowSize() { return spec_window_.SizeZ()-2*spwin_ext_sz_; }
+// Taille de la zone d'extension de la fenetre de chaque cote 
+  inline sa_size_t GetSpecWinExtensionSize() { return spwin_ext_sz_; }
+// La taille de recouvrement est le double de la taille de la zone d'extension 
+  inline sa_size_t GetSpecWinOverlapSize() { return 2*spwin_ext_sz_; }
+// Retourne l'index (selon Z) dans le tableau des spectres en fonction du numero de paquet
+  inline sa_size_t PaqNumToArrayIndex(uint_8 numpaq)  { return numpaq%spec_window_.SizeZ() ; }
+
+// Pour definir le nombre maximum de fichiers de fenetres de spectres 
+  inline void SetMaxNbSepcWinFiles(uint_4 nmax=0)  {  nmaxfiles_specw_=nmax; }
+
   inline void SetVarianceLimits(double vmin=0., double vmax=9.e99)
-     { varmin_=vmin;  varmax_=vmax; }
-
-  void ReadGainFitsFile(string filename);
-
-protected:
-  virtual int Process(); 
-  virtual void FlagBadPackets();
-  virtual void SaveSpectra();
+    { varmin_=vmin;  varmax_=vmax; }
+
+// lecture des gains en fonction de la frequence depuis un fichier FITS
+// flag fgapp: Active ou desactive la renormalisation de coefficients de Fourier par les gains=f(frequence)  
+  void ReadGainFitsFile(string filename, bool fgapp=true);
+
+protected:
+  virtual int Process();   // appele paquet par paquet 
+// appele par Process() pour traiter un groupe de paquets correspondant a une fenetre en temps
+  virtual void ProcSpecWin(uint_8 numpaqstart, uint_8 numpaqend);  // appele fenetre par fenetre  
+  virtual void FlagBadPackets(uint_8 numpaqstart, uint_8 numpaqend);   
+  virtual void SaveMeanSpectra();  // Pour sauver les spectres moyennes ( + sigma ) dans un fichier 
+  virtual void SaveSpectraWindow();  // Pour sauver les spectres de la fenetre temporel ds un fichier 
 
 
@@ -56 +75 @@
   bool fgsinglechannel_;
   uint_4 numfile_;
-  uint_8 nbpaq4mean_;
-
+  vector<uint_8> nbpaq4mean_;
   uint_8 totnbpaq_;
+
   TMatrix< r_4 > mspecmtx_;    // Matrice des spectres moyennees 
   TMatrix< r_4 > sigspecmtx_;    // Matrice des sigmas des spectres  
 
+  TArray< r_4 > spec_window_;    // fenetre en temps sur les spectres  
+  uint_4 spwin_ext_sz_;          // Taille de la zone d'extension en debut et en fin 
+  TArray< uint_2 > clnflg_;      // flag de nettoyage resultant de la routine   FlagBadPackets() 
+  uint_8 paqnum_w_start;     // numero de paquet du debut de la fenetre en temps (inclusif) 
+  uint_8 paqnum_w_end;       // numero de paquet final de la fenetre en temps (non inclusif) 
+
+  uint_4 nbtot_specwin_;    // nombre total de fenetre de spectres traites 
+  uint_4 nmaxfiles_specw_;   // nombre maximum de fichier de fenetre_spectre 
+  
   TMatrix< r_4 > sgain_;   // gain en fonction de la frequence  
-
-  vector<int> clnflg_;      // resultat de la routine "cleaner" 
+  bool fg_apply_gains_; 
+
   double varmin_, varmax_;  // limites pour la coupure sur la variance
   vector<uint_8> nbadpaq_;