Changeset 3930 in Sophya for trunk/Cosmo/RadioBeam/specpk.h

Timestamp:

Dec 23, 2010, 12:49:22 AM (15 years ago)

Author:

ansari

Message:

ajout config double hexagon + programme calcul pknoise a partir du resultat de repicon.cc, Reza 22/12/2010

File:

• trunk/Cosmo/RadioBeam/specpk.h (modified) (2 diffs)

trunk/Cosmo/RadioBeam/specpk.h

@@ -48 +48 @@
   inline void SetCellSize(double dkx=DeuxPI, double dky=DeuxPI, double dkz=DeuxPI) 
   { dkx_=dkx;  dky_=dky;  dkz_=dkz; } 
-  inline int SetPrtLevel(int lev=0) 
-  { int olev=prtlev_; prtlev_=lev; return olev; }
+  inline int SetPrtLevel(int lev=0, int prtmod=10) 
+  { int olev=prtlev_; prtlev_=lev; prtmodulo_=prtmod; return olev; }
   void ComputeFourierAmp(SpectralShape& pk);
   void ComputeNoiseFourierAmp(Four2DResponse& resp, bool crmask=false);
@@ -71 +71 @@
   double dkx_, dky_, dkz_;
   int prtlev_;
+  int prtmodulo_;
 };
+
+// --- PkNoiseCalculator : 
+// - Classe de calcul du spectre de bruit PNoise(k) determine par une reponse
+//   2D de l'instrument 
+class PkNoiseCalculator
+{
+public:
+  PkNoiseCalculator(Four3DPk& pk3, Four2DResponse& rep, double s2cut=100., int ngen=1, const char* tit="PkNoise");
+  
+  inline void SetS2Cut(double s2cut=100.)
+  {  S2CUT=s2cut; }
+  inline double GetS2Cut() { return S2CUT; }
+  HProf Compute();
+  inline int SetPrtLevel(int lev=0, int prtmod=10) 
+  { int olev=prtlev_; prtlev_=lev; prtmodulo_=prtmod; return olev; }
+
+protected:
+  Four3DPk& pkn3d;
+  Four2DResponse& frep;
+  double S2CUT;
+  int NGEN;
+  string title;
+  int prtlev_;
+  int prtmodulo_;
+};
+