Changeset 3192 in Sophya for trunk/Cosmo/RadioBeam/multicyl.h

Timestamp:

Mar 20, 2007, 11:48:15 AM (19 years ago)

Author:

legoff

Message:

input file read with DVcards by treccyl.cc

and by multiCylinders(char* filename)

using real units (meter and GHz)
using _ for member variables (e.g. x_)

File:

• trunk/Cosmo/RadioBeam/multicyl.h (modified) (6 diffs)

trunk/Cosmo/RadioBeam/multicyl.h

@@ -28 +28 @@
   //   longueur @ f=2 ~ 64 (256*lambda/2 = 256*0.25)
  MultiCylinders(int nr=256, int ns=1024);
+ MultiCylinders(char* filename);
  ~MultiCylinders();
 
  // Niveau de print de debug
- inline void SetPrintLevel(int prl=0) { PrtLev=prl; }
+ inline void SetPrintLevel(int prl=0) { PrtLev_=prl; }
 
- // Ajout d'un cylindre, en position posY
+ // Ajout d'un cylindre, en position posX, posY
  inline int AddCylinder(double posx, double posy)
-   {  mCyl.push_back( new MultiBeamCyl(NR, NS, posx, posy) ); return mCyl.size(); }
+   {  mCyl_.push_back( new MultiBeamCyl(NR_, NS_, posx, posy) ); return mCyl_.size(); }
 
  MultiBeamCyl & GetCylinder(int i);
@@ -41 +42 @@
  // Specification de la frequence de base f0 et espacement des recepteurs
  inline void SetBaseFreqDa(double f0=2., double a=0.25) 
-   {  freq0 = f0;   Da = a;  }
+   {  freq0_ = f0;   Da_ = a;  }
 
  // frequences reduites (entre 0 ... 0.5) , angle en radian, amp ~ 1
@@ -49 +50 @@
 
  // Definition du sigma du bruit gaussien sur les echantillons
- inline void SetNoiseSigma(double sig=0.) {  signoise = sig; }
+ inline void SetNoiseSigma(double sig=0.) {  signoise_ = sig; }
  // Definition du sigma du jitter d'horloge (typ 0.01)
- inline void SetTimeJitter(double tjit=0.) { timejitter = tjit; } 
+ inline void SetTimeJitter(double tjit=0.) { timejitter_ = tjit; } 
  // Definition du sigma des offsets d'horloge entre recepteurs (typ 0.02)
- inline void SetTimeOffsetSigma(double tsig=0.) { toffsig = tsig; }
+ inline void SetTimeOffsetSigma(double tsig=0.) { toffsig_ = tsig; }
 
  // Definition du gain et sigma de fluctuations de gain 
  // nzerogain: nb de recepteurs (choisis au hasard) avec gain mis a zero
  inline void SetGains(double g=1., double sigg=0., int nzerogain=0)
-   {  gain=g;  siggain=sigg;  ngainzero = nzerogain; }
+   {  gain_=g;  siggain_=sigg;  ngainzero_ = nzerogain; }
 
 
@@ -72 +73 @@
  // @f = 2 , lambda = 0.5 ===> posY <~ lambda/(2 sin resol) 
  //   ===> posY < ~ 20   
- void ReconstructSourceBox(int halfny=10, double stepangy=M_PI/300);
+// void ReconstructSourceBox(int halfny=10, double stepangy=M_PI/300);
+ void ReconstructSourceBox(int halfny, double stepangy, int ny, double stepangx);
 
 
@@ -78 +80 @@
  // avec la moyenne des modules entre kfmin <= k (selon z) <= kfmax 
  TMatrix< r_4 > getRecXYSlice(int kfmin, int kfmax);
+ TMatrix< r_4 > getRecYXSlice(int kfmin, int kfmax);
 
  // Acces a la boite 3D de sources reconstruite
- inline TArray< complex<r_4> > & getRecSrcBox() { return cmplx_srcbox; }
+ inline TArray< complex<r_4> > & getRecSrcBox() { return cmplx_srcbox_; }
 
  // Configure chaque cylindre a partir des parametres de la classe 
@@ -86 +89 @@
  void ConfigureCylinders();
 
+private:
  //-------------- Variables membres ------------ 
- int NR, NS;  // nb recepteurs, nb d'echantillons ds chaque paquet
+ int NR_, NS_;  // nb recepteurs, nb d'echantillons ds chaque paquet
 
- int PrtLev;   // Niveau de print de debug
+ int PrtLev_;   // Niveau de print de debug
 
- double Da; // distance entre recepteurs
- double freq0;  // frequence de base (freqvrai=f+freq0)
+ double Da_; // distance entre recepteurs
+ double freq0_;  // frequence de base (freqvrai=f+freq0)
  
- double timejitter; // le sigma du jitter en temps des coups d'horloge
- double toffsig;    // sigma des offsets en temps
+ double timejitter_; // le sigma du jitter en temps des coups d'horloge
+ double toffsig_;    // sigma des offsets en temps
  
- BRSourceGen* src;  // Les sources
- bool adfg; // if true, delete src 
+ BRSourceGen* src_;  // Les sources
+ bool adfg_; // if true, delete src 
 
- double signoise; // sigma du bruit additif (bruit ampli ... )
- double gain, siggain;
- int ngainzero;
+ double signoise_; // sigma du bruit additif (bruit ampli ... )
+ double gain_, siggain_;
+ int ngainzero_;
 
- vector<MultiBeamCyl *> mCyl;  // Les differents cylindres 
+ vector<MultiBeamCyl *> mCyl_;  // Les differents cylindres 
 
- TArray< complex<r_4> > cmplx_srcbox; // boite3D des sources (angX,Y,freq) reconstruit
+ TArray< complex<r_4> > cmplx_srcbox_; // boite3D des sources (angX,Y,freq) reconstruit
  };