Changeset 3948 in Sophya for trunk/Cosmo/RadioBeam/pknoise.cc

Timestamp:

Feb 16, 2011, 3:20:54 PM (15 years ago)

Author:

ansari

Message:

modifications de config cylindres pour papier, Reza 16/02/2011

File:

• trunk/Cosmo/RadioBeam/pknoise.cc (modified) (6 diffs)

trunk/Cosmo/RadioBeam/pknoise.cc

@@ -43 +43 @@
        << "    -noise NoiseLevel (default=1.) \n"
        << "    -renmax MaxValue (default : Do NOT renormalize 2D response value \n"   
-       << "    -scut SCutValue (default=100.) \n"
        << "    -ngen NGen (default=0) number of noise fourier amp generations \n"
        << "       NGen=0 -> Call ComputeNoisePk(), else generate Fourier Amplitudes (random) \n"
        << "    -z redshift (default=0.7) \n"
+       << "    -scut SCutValue (default= -100.) \n"
+       << "       if SCutValue<0. ==> SCut=MinNoisePower*(-SCutValue) \n"
        << "    -prt PrtLev,PrtModulo (default=0,10) " << endl;
   return;
@@ -76 +77 @@
   double rmax=1.;
   int NMAX = 0;
-  double SCut=0.;
+  double SCut=-100.;
+  bool fgautoscut=true;
+  double FacSCut=-SCut;
   double z_Redshift=0.7 ;  // 21 cm at z=0.7 -> 0.357 m  
   int prtlev=0;
@@ -92 +95 @@
     else if (strcmp(arg[ka],"-scut")==0) {
       SCut=atof(arg[ka+1]);    ka+=2; 
+      if (SCut<0.) { FacSCut=-SCut;  fgautoscut=true; }
     }
     else if (strcmp(arg[ka],"-ngen")==0) {
@@ -154 +158 @@
     else cout << " pknoise[1]: Four2DResponse ( Diameter=" << DIAMETRE << " Lambda= " << lambda
               << " DoL=" << DIAMETRE/lambda << " ) " << endl;
-
+    Histo2D h2drep = arep_p->GetResponse();
+    double repmax= h2drep.VMax();
+    if (fgautoscut) {
+      SCut = FacSCut/repmax;
+      cout << " pknoise[1.b]: Four2DResponse.RepMax=" << repmax << " --> SCut=" << FacSCut << "/repmax=" 
+           << SCut << endl;
+    }
+    else cout << " pknoise[1.b]: Four2DResponse.RepMax=" << repmax << " , SCut=" << SCut << endl;
     
     cout << " pknoise[2]: Instanciating object type Four3DPk  " << endl;
@@ -167 +178 @@
       pkn.SetPrtLevel(prtlev,prtmod);
       HProf hpn = pkn.Compute();
+      cout << " pknoise[3.b]: writing hpn noise profile to " << outfile << endl;
       po << hpn;
     }
@@ -173 +185 @@
       DataTable dtnoise;
       HProf hpn = m3d.ComputeNoisePk(*(arep_p),fracmodok,dtnoise,SCut);
-      po << hpn;
-      string outfile2 = "x"+outfile;
-      POutPersist po2(outfile2);
       HProf h1dnoise=arep_p->GetProjNoiseLevel();
       HProf h1drep=arep_p->GetProjResponse();
-      po2 << PPFNameTag("dtnoise") << dtnoise;
-      po2 << PPFNameTag("hpnoise") << hpn;
-      po2 << PPFNameTag("fracmodok") << fracmodok;
-      po2 << PPFNameTag("h1dnoise") << h1dnoise;
-      po2 << PPFNameTag("h1drep") << h1drep;
+      cout << " pknoise[3.b]: writing dtnoise,hpn,h2rep... with tags to " << outfile << endl;
+      po << PPFNameTag("dtnoise") << dtnoise;
+      po << PPFNameTag("hpnoise") << hpn;
+      po << PPFNameTag("fracmodok") << fracmodok;
+      po << PPFNameTag("h1dnoise") << h1dnoise;
+      po << PPFNameTag("h1drep") << h1drep;
+      po << PPFNameTag("h2drep") << h2drep;
     }
     rc = 0;