Changeset 3192 in Sophya for trunk/Cosmo/RadioBeam/mbeamcyl.h
- Timestamp:
- Mar 20, 2007, 11:48:15 AM (19 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/Cosmo/RadioBeam/mbeamcyl.h
r3191 r3192 31 31 // ns = nb d'echantillons en temps de chaque paquet 32 32 // posy = position spatiale en Y pour chaque cylindre (unite ou lambda=1 pour f=1) 33 MultiBeamCyl(int nr=1024, int ns=4096, double posx=0., double posy=0.); 33 // posx = position spatiale en X pour chaque cylindre (unite ou lambda=1 pour f=1) 34 MultiBeamCyl(int nr=1024, int ns=4096, double posx=0., double posy=0. ); 34 35 ~MultiBeamCyl(); 35 36 36 37 // Niveau de print de debug 37 inline void SetPrintLevel(int prl=0) { PrtLev =prl; }38 inline void SetPrintLevel(int prl=0) { PrtLev_=prl; } 38 39 39 40 // Reponse angulaire pour un recepteur (cylindre+antenne) … … 41 42 inline double AngResponse(double angx, double angy) { return 1.; } 42 43 43 inline double getCylinderYPos() { return posY; } 44 inline double getCylinderXPos() { return posX; } 44 inline double getCylinderYPos() { return posY_; } 45 inline double getCylinderXPos() { return posX_; } 46 45 47 // Specification de la frequence de base f0 et espacement des recepteurs 46 48 // freq_vrai = freq_BRSourceGen + f0 47 49 // pour f_vrai = 2, lambda = 0.5 et lambda/2 = 0.25 48 50 inline void SetBaseFreqDa(double f0=2., double a=0.25) 49 { freq0 = f0; Da= a; }51 { freq0_ = f0; Da_ = a; } 50 52 51 53 // frequences reduites (entre 0 ... 0.5) , angle en radian, amp ~ 1 … … 55 57 56 58 // Definition du sigma du bruit gaussien sur les echantillons 57 inline void SetNoiseSigma(double sig=0.) { signoise = sig; }59 inline void SetNoiseSigma(double sig=0.) { signoise_ = sig; } 58 60 // Definition du sigma du jitter d'horloge (typ 0.01) 59 inline void SetTimeJitter(double tjit=0.) { timejitter = tjit; }61 inline void SetTimeJitter(double tjit=0.) { timejitter_ = tjit; } 60 62 // Definition du sigma des offsets d'horloge entre recepteurs (typ 0.02) 61 inline void SetTimeOffsetSigma(double tsig=0.) { toffsig = tsig; }63 inline void SetTimeOffsetSigma(double tsig=0.) { toffsig_ = tsig; } 62 64 63 65 // Definition du gain et sigma de fluctuations de gain … … 79 81 void ReconstructSourcePlane(bool fgzerocentre=true); 80 82 81 inline TMatrix< complex<r_4> > & getRecSrcPlane() { return cmplx_srcplane ; }83 inline TMatrix< complex<r_4> > & getRecSrcPlane() { return cmplx_srcplane_; } 82 84 83 85 //-------------- Variables membres ------------ 84 int NR , NS; // nb recepteurs, nb d'echantillons ds chaque paquet85 double pos Y; // Position selon Y (E-O) de la ligne de recepteurs86 double pos X; // Position selon X (N-S) de la ligne de recepteurs86 int NR_, NS_; // nb recepteurs, nb d'echantillons ds chaque paquet 87 double posX_; // Position selon X (N-S) du premier recepteur du cylindre 88 double posY_; // Position selon Y (E-O) de la ligne de recepteurs 87 89 88 int PrtLev ; // Niveau de print de debug90 int PrtLev_; // Niveau de print de debug 89 91 90 double Da ; // distance entre recepteurs91 double freq0 ; // frequence de base (freqvrai=f+freq0)92 double Da_; // distance entre recepteurs 93 double freq0_; // frequence de base (freqvrai=f+freq0) 92 94 93 Vector texact ; // les temps exact94 Vector tjitt ; // temps avec jitter95 Vector toffset ; // Offset en temps entre recepteurs96 double timejitter ; // le sigma du jitter en temps des coups d'horloge97 double toffsig ; // sigma des offsets en temps95 Vector texact_; // les temps exact 96 Vector tjitt_; // temps avec jitter 97 Vector toffset_; // Offset en temps entre recepteurs 98 double timejitter_; // le sigma du jitter en temps des coups d'horloge 99 double toffsig_; // sigma des offsets en temps 98 100 99 BRSourceGen* src ; // Les sources100 bool adfg ; // if true, delete src101 BRSourceGen* src_; // Les sources 102 bool adfg_; // if true, delete src 101 103 102 double signoise ; // sigma du bruit additif (bruit ampli ... )104 double signoise_; // sigma du bruit additif (bruit ampli ... ) 103 105 104 TVector<r_4> gain ; // Gain de chaque recepteur106 TVector<r_4> gain_; // Gain de chaque recepteur 105 107 106 TVector<r_4> signal ; // signal sans jitter en temps107 TVector<r_4> sigjitt ; // signal avec jitter en temps108 TVector<r_4> signal_; // signal sans jitter en temps 109 TVector<r_4> sigjitt_; // signal avec jitter en temps 108 110 109 TVector< complex<r_4> > f_sig , f_sigjit; // TF des vecteurs signal , sigjitt111 TVector< complex<r_4> > f_sig_, f_sigjit_; // TF des vecteurs signal , sigjitt 110 112 111 TMatrix< complex<r_4> > cmplx_srcplane ; // composantes complexe du plan-source reconstruit112 // TMatrix<r_4> rec_srcplane ; // Matrice plan source colonnes->angX, lignes->freq113 TMatrix< complex<r_4> > cmplx_srcplane_; // composantes complexe du plan-source reconstruit 114 // TMatrix<r_4> rec_srcplane_; // Matrice plan source colonnes->angX, lignes->freq 113 115 }; 114 116
Note:
See TracChangeset
for help on using the changeset viewer.