Changeset 3947 in Sophya for trunk/Cosmo/RadioBeam/mdish.h

Timestamp:

Feb 14, 2011, 12:58:29 AM (15 years ago)

Author:

ansari

Message:

Amelioration et modifs diverses lors de la redaction du papier, Reza 13/02/2011

File:

• trunk/Cosmo/RadioBeam/mdish.h (modified) (8 diffs)

trunk/Cosmo/RadioBeam/mdish.h

@@ -42 +42 @@
   inline void setLambda(double lambda=1.)
   { lambda_ = lambda;   lambda_ratio_ = lambda_/lambdaref_; }   
+
+  inline double getLambdaRef()   { return lambdaref_; }
+  inline double getLambda()  { return lambda_; } 
   
   // Return the 2D response for wave vector (kx,ky)
@@ -48 +51 @@
   { return Value(kx, ky); } 
   virtual Histo2D GetResponse(int nx=255, int ny=255);  
+  // Retourne le niveau moyen du bruit projete 1D en fonction (sqrt(u^2+v^2)
+  HProf GetProjNoiseLevel(int nbin=128, bool fgnorm1=true);
+  // Retourne la reponse moyenne projetee 1D en fonction (sqrt(u^2+v^2)
+  HProf GetProjResponse(int nbin=128, bool fgnorm1=true);
+
   inline double D() { return dx_; } ; 
   inline double Dx() { return dx_; } ; 
@@ -95 +103 @@
   // Circular dish 
   Dish(int id, double x, double y, double diam)
-    :  id_(id), X(x), Y(y), D(diam), Dx(D), Dy(D), fgcirc_(true)   {   } 
+    :  id_(id), X(x), Y(y), D(diam), Dx(D), Dy(D), fgcirc_(true), gain_(1.)   {   } 
   // Receiver with rectangular type answer in kx,ky plane 
   Dish(int id, double x, double y, double dx, double dy)
-    :  id_(id), X(x), Y(y), D(sqrt(dx*dy)), Dx(dx), Dy(dy), fgcirc_(false)   {   } 
+    :  id_(id), X(x), Y(y), D(sqrt(dx*dy)), Dx(dx), Dy(dy), fgcirc_(false), gain_(1.)    {   } 
   Dish(Dish const& a) 
-    :  id_(a.id_), X(a.X), Y(a.Y), D(a.D), Dx(a.Dx), Dy(a.Dy), fgcirc_(a.fgcirc_)     {   } 
+    :  id_(a.id_), X(a.X), Y(a.Y), D(a.D), Dx(a.Dx), Dy(a.Dy), fgcirc_(a.fgcirc_), gain_(a.gain_)  {   } 
+  inline void setGain(double gain) { gain_=gain; return; }
   inline bool isCircular() { return fgcirc_; }
   inline int ReflectorId() { return id_; }
@@ -106 +115 @@
   inline double DiameterX() { return Dx; }
   inline double DiameterY() { return Dy; }
+  inline double Gain() { return gain_; }
 
   int id_;   // numero de reflecteur
@@ -111 +121 @@
   double Dx, Dy;
   bool fgcirc_;  // false -> rectangular dish 
+  double gain_;
 };
 
@@ -132 +143 @@
   { beamnx_=nx; beamny_=ny; }
 
+  // Calcul la reponse ds le plan 2D  (u,v) = (kx,ky)
+  void ComputeResponse();
+  // Retourne la reponse 2D ds le plan (u,v) = (kx,ky) sous forme d'histo 2D 
   Histo2D GetResponse();
+
+  // Retourne le niveau moyen du bruit projete 1D en fonction (sqrt(u^2+v^2)
+  HProf GetProjNoiseLevel(int nbin=128, bool fgnorm1=true);
+  // Retourne la reponse moyenne projetee 1D en fonction (sqrt(u^2+v^2)
+  HProf GetProjResponse(int nbin=128, bool fgnorm1=true);
 
   double CumulResp(Four2DResponse& rd, double theta=0., double phi=0.);
@@ -143 +162 @@
   double AddToHisto(double kx0, double ky0, double x, double y, double w, bool fgfh);
 
-  double lambda_, dmax_;
+  double lambda_, dmax_, kmax_;
   vector<Dish> dishes_;
   bool fgnoauto_;
@@ -153 +172 @@
   //   Histo2D h2w_, h2cnt_;
   QHis2D h2w_;
+  bool fgcomputedone_;
   int mcnt_;
   int prtlev_,prtmodulo_;