Changeset 4022 in Sophya for trunk/Cosmo/RadioBeam/applobe.cc

Timestamp:

Sep 28, 2011, 5:13:51 PM (14 years ago)

Author:

ansari

Message:

modifs pour pouvoir imposer la moyenne en temp des plans X,Y des cubes lors de l'extraction du signal HI, Reza 28/9/2011

File:

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

trunk/Cosmo/RadioBeam/applobe.cc

@@ -58 +58 @@
 
   if ((narg < 3)||(strcmp(arg[1],"-h")==0)) {
-    cout << "Usage: applobe [-t -g -fib -mxr val] Diameter/Four2DRespTableFile In3DPPFName Out3DPPFName \n" 
+    cout << "Usage: applobe [-t -g -fib -kzf/-nzf -mxr val] Diameter/Four2DRespTableFile In3DPPFName Out3DPPFName \n" 
          << "               [TargetBeamDiam] [ResmapleFactor=0.5,0.333...] \n" << endl;
     if ((narg>1)&&(strcmp(arg[1],"-h")==0)) {
       cout << "   -t -g : Triangular / gaussian beam shape (def=gaussian) \n" 
            << "   -fib :  Application of a fixed (freq.independent) lobe dish-triangle or gaussian \n" 
+           << "   -kzf -nzf :  Keep (default) or Not zero space frequency when applying lobes (BeamEffect class) \n" 
            << "   -mxr val: Max beam correction factor (default=10.) \n"
            << "   Diameter/Four2DRespTableFile : dish diameter or 2D response PPF file name\n" 
@@ -79 +80 @@
   bool fixedbeam=false;  // true -> apply freq. independent beam 
   double maxratio=10.;   // valeur max du rapport des lobes lors de la correction de lobe 
+  bool preservefreq0=true; // true -> keep zero frequency when appyling lobe 
 
   // decodage argument optionnel 
@@ -87 +89 @@
     else if (fbo=="-g")  { fggaussian=true; arg++; narg--; }
     else if (fbo=="-fib")  { fixedbeam=true; arg++; narg--; }
+    else if (fbo=="-kzf")  { preservefreq0=true; arg++; narg--; }
+    else if (fbo=="-nzf")  { preservefreq0=false; arg++; narg--; }
     else if (fbo=="-mxr")  { arg++; maxratio=atof(arg[1]); arg++; narg-=2; }
     else fgoptarg=false;
@@ -164 +168 @@
     else cout << " applobe[2.b]: Four2DResponse ( Diameter=" << DIAMETRE << " Lambda= " << lambda
               << " DoL=" << DIAMETRE/lambda << " ) " << endl;
-    BeamEffect beam(*fresp_p);
+    BeamEffect beam(*fresp_p,preservefreq0);
 
     if (fgcorrbeam) {