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

Timestamp:

Oct 10, 2011, 11:31:03 PM (14 years ago)

Author:

ansari

Message:

Modif programme pknoise.cc et classes Four3DPk, PkNoiseCalculator pour calcul spectre de bruit en faisant varier la reponse instrumentale avec la frequence, Reza 10/10/2011

File:

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

trunk/Cosmo/RadioBeam/specpk.h

@@ -51 +51 @@
   { int olev=prtlev_; prtlev_=lev; prtmodulo_=prtmod; return olev; }
   void ComputeFourierAmp(SpectralShape& pk);
-  void ComputeNoiseFourierAmp(Four2DResponse& resp, bool crmask=false);
+// angscale is a multiplicative factor converting transverse k (wave number) values to angular wave numbers 
+// typically = ComovRadialDistance 
+  void ComputeNoiseFourierAmp(Four2DResponse& resp, double angscale=1., bool crmask=false);
+  void ComputeNoiseFourierAmp(Four2DResponse& resp, double f0, double df, double angscale=1.);
+
 // Return the array size 
   inline sa_size_t NCells() { return fourAmp.Size(); }
@@ -65 +69 @@
   void  ComputePkCumul(HProf& hp, double s2cut=0.);
 
-  HProf ComputeNoisePk(Four2DResponse& resp, Histo& fracmodok, DataTable& dt, 
+// angscale is a multiplicative factor converting transverse k (wave number) values to angular wave numbers 
+// typically = ComovRadialDistance 
+  HProf ComputeNoisePk(Four2DResponse& resp, Histo& fracmodok, DataTable& dt, double angscale=1.,
                        double s2cut=0., int nbin=256, double kmin=0., double kmax=-1.);
 
@@ -85 +91 @@
   PkNoiseCalculator(Four3DPk& pk3, Four2DResponse& rep, double s2cut=100., int ngen=1, const char* tit="PkNoise");
   
+  inline void SetFreqRange(double freq0=835.,double dfreq=0.5)
+  { freq0_=freq0;  dfreq_=dfreq; }
+  inline void SetAngScaleConversion(double angscale=1.)
+  { angscale_=angscale; }
   inline void SetS2Cut(double s2cut=100.)
   {  S2CUT=s2cut; }
@@ -95 +105 @@
   Four3DPk& pkn3d;
   Four2DResponse& frep;
+  double freq0_,dfreq_;
+  double angscale_; 
   double S2CUT;
   int NGEN;