1 | #ifndef BRPROC_H_SEEN
|
---|
2 | #define BRPROC_H_SEEN
|
---|
3 |
|
---|
4 | //----------------------------------------------------------------
|
---|
5 | // Projet BAORadio - (C) LAL/IRFU 2008-2010
|
---|
6 | // Classes de threads de traitememt donnees BAORadio
|
---|
7 | //----------------------------------------------------------------
|
---|
8 |
|
---|
9 |
|
---|
10 | #include <string>
|
---|
11 | #include <vector>
|
---|
12 | #include <iostream>
|
---|
13 |
|
---|
14 | #include "racqumem.h"
|
---|
15 | #include "brbaseproc.h"
|
---|
16 | #include "parlex.h"
|
---|
17 | #include "sopnamsp.h"
|
---|
18 | #include "tvector.h"
|
---|
19 | #include "ntuple.h"
|
---|
20 | #include "datatable.h"
|
---|
21 |
|
---|
22 | using namespace std;
|
---|
23 |
|
---|
24 |
|
---|
25 | //---------------------------------------------------------------------
|
---|
26 | // Classe de traitement - calcul de visibilite pour n fibres
|
---|
27 | //---------------------------------------------------------------------
|
---|
28 | class BRVisibilityCalculator : public BRBaseProcessor, ParallelTaskInterface {
|
---|
29 | public:
|
---|
30 | BRVisibilityCalculator(RAcqMemZoneMgr& memgr, string outpath, uint_4 nmean=1000, size_t nthr=1);
|
---|
31 | virtual ~BRVisibilityCalculator(); // Fait le nettoyage final , sauvegarde des objets ...
|
---|
32 |
|
---|
33 | // Pour definir le rang de l'objet lors d'un traitement avec plusieurs BRVisibilityCalculator en //
|
---|
34 | void DefineRank(uint_4 nbc, uint_4 cid);
|
---|
35 |
|
---|
36 | // Pour limiter le calcul des visibilites a certaines fibres
|
---|
37 | uint_4 SelectPairs(uint_4 pair1=0, uint_4 nbpairs=99999999);
|
---|
38 | // pour definir le binning en frequence pour la VisibDataTable
|
---|
39 | int SelectFreqBinning(uint_4 freq1=0, uint_4 freq2=0, uint_4 nbfreq=1);
|
---|
40 |
|
---|
41 | // Pour activer la verification des timetag
|
---|
42 | // Creation / sauvegarde d'une matrice de TimeTag / FrameCounter avec NCols=maxnpaq (nb de paquets/trigger)
|
---|
43 | int ActivateTimeTagCheck(uint_8 maxnpaq);
|
---|
44 |
|
---|
45 | // Methode d'execution standard du thread (on redefinit la methode de BRBaseProcessor)
|
---|
46 | virtual void run();
|
---|
47 | // Methode de la classe ParallelTaskInterface
|
---|
48 | virtual int execute(int tid);
|
---|
49 |
|
---|
50 | // Le nombre total d'operations effectuees
|
---|
51 | inline r_8 TotNbFLOP() { return nb_flop_; }
|
---|
52 |
|
---|
53 | protected:
|
---|
54 | virtual int Process();
|
---|
55 | virtual int FillVisibTable(double fcm, double ttm);
|
---|
56 | virtual int CheckTimeTag();
|
---|
57 |
|
---|
58 | ParallelExecutor paralex_; // Pour l'execution en parallele
|
---|
59 | size_t nparthr_;
|
---|
60 |
|
---|
61 | uint_4 nmean_; // Nombre de spectres pour le calcul des moyennes
|
---|
62 | string outpath_; // directory pour fichiers de sortie
|
---|
63 | uint_4 nbcalc_, calcid_; // Nombre total de threads/objets BRVisibilityCalculator parallèles, l'identificateur de l'objet
|
---|
64 |
|
---|
65 | uint_4 pairst_, nbpairs_; // Numero de la premiere paire, nombre de paires
|
---|
66 | uint_4 jf1_,jf2_,djf_; // binning en frequence pour la datatable
|
---|
67 |
|
---|
68 | vector<TwoByteComplex*> vpdata_;
|
---|
69 | bool* fgallfibok;
|
---|
70 |
|
---|
71 | double moyfc_, moytt_;
|
---|
72 | TMatrix< complex<r_4> > vismtx_;
|
---|
73 | TVector< uint_4 > chanum_;
|
---|
74 | // FitsInOutFile dtfos_;
|
---|
75 | DataTable visdt_;
|
---|
76 |
|
---|
77 | uint_4 numfile_;
|
---|
78 | uint_8 totnbpaq_;
|
---|
79 |
|
---|
80 | r_8 nb_flop_; // Nombre total d'operations effectuees
|
---|
81 |
|
---|
82 | // Variables/objets pour la verifications des timetag
|
---|
83 | bool fgcktt_; // true -> checktimetag
|
---|
84 | int_8 mindeltatt_; // difference minimum de deux timetag successifs
|
---|
85 | TMatrix< uint_8 > fcmtx_; // matrices des framecounters
|
---|
86 | TMatrix< uint_8 > ttmtx_; // matrices des timetag
|
---|
87 | vector< uint_8 > vlasttt_; // dernieres valeurs de TimeTag / fibre
|
---|
88 | vector< uint_8 > vbadtt_; // Number of bad TimeTag paquets per fiber
|
---|
89 | vector< uint_8 > vnsamett_; // Number of paquets having TimeTag <> TT_Fiber1 per fiber
|
---|
90 | vector< uint_8 > vndiff5tt_; // Number of paquets having |TimeTag-TT_Fiber1|>5 per fiber
|
---|
91 |
|
---|
92 | };
|
---|
93 |
|
---|
94 |
|
---|
95 | //-------------------------------------------------------------------------------
|
---|
96 | // Classe Groupe (ensemble) de Calculateur de Visibilites, tourant en parallele
|
---|
97 | //-------------------------------------------------------------------------------
|
---|
98 |
|
---|
99 | class BRVisCalcGroup {
|
---|
100 | public:
|
---|
101 | BRVisCalcGroup(size_t nbcalc, RAcqMemZoneMgr& memgr, string outpath, uint_4 nmean=1000, size_t nthr=1);
|
---|
102 | ~BRVisCalcGroup();
|
---|
103 | int SelectFreqBinning(uint_4 freq1=0, uint_4 freq2=0, uint_4 nbfreq=1);
|
---|
104 | inline int ActivateTimeTagCheck(uint_8 maxnpaq)
|
---|
105 | { return viscalcp_[0]->ActivateTimeTagCheck(maxnpaq); }
|
---|
106 |
|
---|
107 | void start();
|
---|
108 | void join();
|
---|
109 | protected:
|
---|
110 | vector<BRVisibilityCalculator *> viscalcp_;
|
---|
111 | Timer tm_;
|
---|
112 | };
|
---|
113 |
|
---|
114 | //---------------------------------------------------------------------
|
---|
115 | // Classe de traitement simple - calcul de spectres moyennes / voie
|
---|
116 | //---------------------------------------------------------------------
|
---|
117 | class BRMeanSpecCalculator : public BRBaseProcessor
|
---|
118 | {
|
---|
119 | public:
|
---|
120 | BRMeanSpecCalculator(RAcqMemZoneMgr& memgr, string outpath, uint_4 nmean=1000);
|
---|
121 | virtual ~BRMeanSpecCalculator(); // Fait le nettoyage final , sauvegarde des objets ...
|
---|
122 |
|
---|
123 | protected:
|
---|
124 | virtual int Process();
|
---|
125 |
|
---|
126 | uint_4 nmean_; // Nombre de spectres pour le calcul des moyennes
|
---|
127 | string outpath_; // directory pour fichiers de sortie
|
---|
128 | uint_4 numfile_;
|
---|
129 |
|
---|
130 | uint_8 totnbpaq_;
|
---|
131 | TMatrix< r_4 > mspecmtx_; // Matrice des spectres moyennees
|
---|
132 | };
|
---|
133 |
|
---|
134 | //-------------------------------------------------------------------------
|
---|
135 | // **** OBSOLETE : BRProcA2C et BRProcB4C
|
---|
136 | // Anciennes classes de traitement (Avant Nov2009) utilises par mcrd.cc
|
---|
137 | //-------------------------------------------------------------------------
|
---|
138 | //---------------------------------------------------------------------
|
---|
139 | // Classe thread de traitement avec 2 voies par frame (donnees raw/FFT
|
---|
140 | //---------------------------------------------------------------------
|
---|
141 |
|
---|
142 | /*-- Arguments du constructeur :
|
---|
143 | o mem : gestionnaire zones memoires avec contenant les "frames" avec 2 voies/frame
|
---|
144 | o path : Chemin et nom des fichiers PPF produits (on ajoute _numero.ppf )
|
---|
145 | o nmean : Nombre de frame (paquets) moyennes
|
---|
146 | o nmax : Nombre maxi de frames traites
|
---|
147 | o fghist = true -> fait aussi l'histo des valeurs des echantillons temps
|
---|
148 | o nfsmap>0 : fabrication de cartes (image 2D) des spectres,
|
---|
149 | avec nfsmap=nombre de frequences sommes
|
---|
150 | o fgnotrl = true -> fichiers fits / paquets sans trailer
|
---|
151 | o card : numero de carte , juste pour les impressions
|
---|
152 | */
|
---|
153 | class BRProcA2C : public ZThread {
|
---|
154 | public:
|
---|
155 | BRProcA2C(RAcqMemZoneMgr& mem, string& path, bool fgraw,
|
---|
156 | uint_4 nmean=10, uint_4 nmax=100, bool fgtshist=false,
|
---|
157 | uint_4 nfsmap=0, bool fgnotrl=false, int card=1);
|
---|
158 | virtual void run();
|
---|
159 | void Stop();
|
---|
160 | inline void STOP() { stop_ = true; }
|
---|
161 | protected:
|
---|
162 | RAcqMemZoneMgr& memgr;
|
---|
163 | bool stop_;
|
---|
164 | bool fgraw_; // true -> raw data
|
---|
165 | uint_4 nmax_; // Nombre maxi de blocs MemZone traites
|
---|
166 | uint_4 nmean_; // Nombre de spectres pour le calcul des moyennes
|
---|
167 | uint_4 nfsmap_; // Binning en frequence pour les cartes 2D temps-freq
|
---|
168 | string path_;
|
---|
169 | bool fghist_; // if true, fill and save histogram with time sample values
|
---|
170 | bool fgnotrl_; // if true, don't check packet trailer - when using pre-june fits files
|
---|
171 | int card_; // Numeros de voies 2*card_-1, 2*card_
|
---|
172 | };
|
---|
173 |
|
---|
174 | //---------------------------------------------------------------------
|
---|
175 | // Classe thread de traitement 2 x 2 voies/paquet (Apres BRProcA2C)
|
---|
176 | //---------------------------------------------------------------------
|
---|
177 | /*-- Arguments du constructeur :
|
---|
178 | o mem1,mem2 : gestionnairea zones memoires avec contenant les donnees des 4 voies (2 voies/frame)
|
---|
179 | o path : Chemin et nom des fichiers PPF produits (on ajoute _numero.ppf )
|
---|
180 | o nmean : Nombre de frame (paquets) moyennes
|
---|
181 | o nmax : Nombre maxi de frames traites
|
---|
182 | o fgnotrl = true -> fichiers fits / paquets sans trailer
|
---|
183 | o card : numero de carte , juste pour les impressions
|
---|
184 | */
|
---|
185 |
|
---|
186 | class BRProcB4C : public ZThread {
|
---|
187 | public:
|
---|
188 | BRProcB4C(RAcqMemZoneMgr& mem1, RAcqMemZoneMgr& mem2, string& path,
|
---|
189 | bool fgraw, uint_4 nmean=10, uint_4 nmax=100, bool fgnotrl=false);
|
---|
190 | virtual void run();
|
---|
191 | void Stop();
|
---|
192 | inline void STOP() { stop_ = true; }
|
---|
193 | protected:
|
---|
194 | RAcqMemZoneMgr& memgr1;
|
---|
195 | RAcqMemZoneMgr& memgr2;
|
---|
196 | bool stop_;
|
---|
197 | bool fgraw_; // true -> raw data
|
---|
198 | uint_4 nmax_; // Nombre maxi de blocs MemZone traites
|
---|
199 | uint_4 nmean_; // Nombre de blocs pour le calcul des moyennes
|
---|
200 | string path_;
|
---|
201 | bool fgnotrl_; // if true, don't check packet trailer - when using pre-june fits files
|
---|
202 | };
|
---|
203 |
|
---|
204 |
|
---|
205 |
|
---|
206 | #endif
|
---|