Changeset 3789 in Sophya for trunk/Cosmo/RadioBeam/calcpk2.cc

Timestamp:

Jun 27, 2010, 4:57:36 PM (15 years ago)

Author:

ansari

Message:

Corrections diverses, Reza 27/06/2010

File:

• trunk/Cosmo/RadioBeam/calcpk2.cc (modified) (4 diffs)

trunk/Cosmo/RadioBeam/calcpk2.cc

@@ -61 +61 @@
     if (narg>6) {
       pixsignoise=atof(arg[6]);
-      fgaddnoise=true;
+      if (pixsignoise>1.e-6)  fgaddnoise=true;
     }
     
     bool fgcorrbeam=true;
-    double tbeamarcmin=30.;
+    double tbeamarcmin=15.;
     if (narg>7) {
       tbeamarcmin=atof(arg[7]);
@@ -119 +119 @@
     conv.setFrequency(Freq0MHz);
     double lambda = conv.getLambda();
-    Four2DResponse arep(2, InterfArrayDiametre/lambda, InterfArrayDiametre/lambda);
+    Four2DResponse arep(2, InterfArrayDiametre/lambda, InterfArrayDiametre/lambda, lambda);
     double DoL = 1.22/ArcminToRadian(tbeamarcmin);
     Four2DResponse tbeam(2, DoL, DoL );
@@ -128 +128 @@
       cleaner.BeamCorrections();
     }
+    cout << " calcpk2[3.b] : calling cleaner.CleanNegatives() ... " << endl;
+    cleaner.CleanNegatives();
     if (fgclnsrc) {
-      cout << "calcpk2[3.b] : calling cleaner.CleanPointSources() with threshold NSigma=" << nsigsrc  << endl;
+      cout << "calcpk2[3.c] : calling cleaner.CleanPointSources() with threshold NSigma=" << nsigsrc  << endl;
       cleaner.CleanPointSources(nsigsrc);
     }
@@ -150 +152 @@
     RandomGenerator rg;
     Four3DPk pkc(four3d, rg);
+    double dkxmpc = DeuxPI/(double)exlss.SizeX()/XCellSizeMpc;
+    double dkympc = DeuxPI/(double)exlss.SizeY()/YCellSizeMpc;
+    double dkzmpc = DeuxPI/(double)exlss.SizeZ()/ZCellSizeMpc;
+    pkc.SetCellSize(dkxmpc, dkympc, dkzmpc);
     
     HProf hp = pkc.ComputePk(0.,256);