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

Timestamp:

May 5, 2011, 7:25:18 PM (14 years ago)

Author:

ansari

Message:

modification rapport maxi a appliquer lors des corrections de beams, Reza 05/05/2011

File:

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

trunk/Cosmo/RadioBeam/calcpk2.cc

@@ -7 +7 @@
     R. Ansari , C. Magneville - Juin 2010 
 
-Usage: calcpk2 [-t -g] InMapLSS convFacLSS InMapSync convFacSync InMapRadioSource convFacRsc OutPkFile 
+Usage: calcpk2 [-t -g -mxr val] InMapLSS convFacLSS InMapSync convFacSync InMapRadioSource convFacRsc OutPkFile 
                [PixNoiseLevel] [Diameter/Four2DRespTableFile] [TargetBeamArcmin] [NSigSrcThr]
 ---------------------------------------------------------------  */
@@ -44 +44 @@
 {
   if ( (narg<6)||((narg>1)&&(strcmp(arg[1],"-h")==0)) ) {
-    cout << " Usage: [-t -g] calcpk2 InMapLSS convFacLSS InMapFgnd convFacFgnd OutPkFile \n" 
+    cout << " Usage: [-t -g -mxr val] calcpk2 InMapLSS convFacLSS InMapFgnd convFacFgnd OutPkFile \n" 
          << "        [PixNoiseLevel] [D_Dish/Four2DRespTableFile CorBeamDiam] \n"
          << "        [NSigSrcThr] [P2/P1] [RecMapFile] " << endl;
     if ((narg>1)&&(strcmp(arg[1],"-h")==0)) {
       cout << "-t -g : Triangular / gaussian beam shape (def=gaussian) \n" 
+           << "-mxr val: Max beam correction factor (default=10.) \n "
            << "- InMapLSS: Input 3D LSS cube (PPF file name) \n "
            << "- convFacLSS: LSS cube conversion factor to mK (milliKelvin) \n"
@@ -78 +79 @@
     // decodage argument optionnel 
     bool fgoptarg=true;
+    double maxratio=10.;
     while (fgoptarg) {
       string fbo = arg[1];
       if (fbo=="-t")  { fggaussian=false; arg++; narg--; }
       else if (fbo=="-g")  { fggaussian=true; arg++; narg--; }
+      else if (fbo=="-mxr")  { arg++; maxratio=atof(arg[1]); arg++; narg-=2; }
       else fgoptarg=false;
     }
@@ -197 +200 @@
     Four2DResponse tbeam(typcb, DoL, DoL );
 
-    ForegroundCleaner  cleaner(*arep_p, tbeam, skycube); 
+    ForegroundCleaner  cleaner(*arep_p, tbeam, skycube, maxratio); 
     if (fgcorrbeam) {
       cout << "calcpk2[3.b] : calling cleaner.BeamCorrections() for target beam Diameter=" << tbeamDiam