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

Timestamp:

Apr 18, 2011, 5:30:44 PM (14 years ago)

Author:

ansari

Message:

Corrections diverses: choix lobe gaussien/triangle et specif DishDiameter au lieu de DoL ds applobe/calcpk2, possibilite application lobe freq.independante ds applobe, Reza 18/04/2011

File:

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

trunk/Cosmo/RadioBeam/applobe.cc

@@ -3 +3 @@
     R. Ansari , C. Magneville - Juin 2010 
 
-  Usage: applobe Diameter/Four2DRespTableFile In3DPPFName Out3DPPFName 
-                 [TargetBeamDoL] [ResmapleFactor=0.5,0.333...] 
+  Usage: applobe [-g -t -fib]  Diameter/Four2DRespTableFile In3DPPFName Out3DPPFName 
+                 [TargetBeamDiam] [ResmapleFactor=0.5,0.333...] 
+   o -t : triangle beam shape in k-space
+   o -g : gaussian beam shape in k-space 
+   o -fib: application d'un lobe fixe (independant de la frequence) 
+   o Diameter/Four2DRespTableFile : diametre dish ou nom de fichier PPF reponse 2D
+   o In3DPPFName Out3DPPFName : Noms des fichiers (TArray<r_4> 3D) entree-sortie 
+   o TargetBeamDiam : Correction à un beam fiduciel, independant de la frequence 
+     avec specification de la valeur du diametre pour la frequence la plus basse
+   o ResmapleFactor : reechantillonnage du cube (2 -> surechantillonnage, 0.5 : 1 sur 2)
 ---------------------------------------------------------------  */
 
@@ -48 +56 @@
 
   InitTim();   // Initializing the CPU timer
+
+  if ((narg < 3)||(strcmp(arg[1],"-h")==0)) {
+    cout << "Usage: applobe [-t -g -fib] 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" 
+           << "   Diameter/Four2DRespTableFile : dish diameter or 2D response PPF file name\n" 
+           << "   In3DPPFName Out3DPPFName: Input/Output PPF file names  (TArray<r_4> 3D) \n" 
+           << "   TargetBeamDiam: Corrected beam target diameter (D/Lambda at lowest frequency) \n" 
+           << "      DoL = 100 --> beam ~ 35 arcmin (D=30m @ z~0.5 Lambda~30cm)) \n"
+           << "   ResmapleFactor : Resampling (2-> oversampling, 0.5 : 1/2 undersampling) \n" << endl;    
+    }  
+    return 1;
+  }
+
   Timer tm("applobe");
 
+  // decodage arguments 
+  bool fggaussian=true;  // true -> gaussian beam 
+  bool fixedbeam=false;  // true -> apply freq. independent beam 
+
+  // decodage argument optionnel 
+  bool fgoptarg=true;
+  while (fgoptarg) {
+    string fbo = arg[1];
+    if (fbo=="-t")  { fggaussian=false; arg++; narg--; }
+    else if (fbo=="-g")  { fggaussian=true; arg++; narg--; }
+    else if (fbo=="-fib")  { fixedbeam=true; arg++; narg--; }
+    else fgoptarg=false;
+  }
   if (narg < 3) {
-    cout << "Usage: applobe Diameter/Four2DRespTableFile In3DPPFName Out3DPPFName " 
-         << "             [TargetBeamDoL] [ResmapleFactor=0.5,0.333...] \n" << endl;
-    return 1;
-  }
-
-  // decodage arguments 
+    cout << " applobe/error arguments , applobe -h for help " << endl;
+    return 2;
+  }
+
   bool fgresptbl=true;
   double DIAMETRE=100.;
@@ -64 +99 @@
     DIAMETRE=atof(arg[1]);
   }
-  else resptblname=arg[1];
-
+  else  resptblname=arg[1];
+  
   string inname = arg[2];
   string outname = arg[3];
   bool fgcorrbeam=false;
-  double tbeamDoL=135;
+  double tbeamDiam=50.;
   if (narg>4) {
-    tbeamDoL=atof(arg[4]);
-    if (tbeamDoL>1)  fgcorrbeam=true;
+    tbeamDiam=atof(arg[4]);
+    if (tbeamDiam>1.)  fgcorrbeam=true;
   }
   bool fgresample=false;
@@ -83 +118 @@
   int rc = 91;
 
-  cout << " ====== applobe :   Input skycube  name= " << inname << " OutName=" << outname;
+  cout << " ====== applobe :   Input skycube  name= " << inname << " OutName=" << outname << endl;
+  cout << ((fggaussian)?" Gaussian ":" Triangular") << ((fixedbeam)?" Fixed (freq.independent)":"") << " beams" << endl;
   bool fginmap=true;
   try {
@@ -113 +149 @@
     double lambda = conv.getLambda();
     
-    Four2DResponse fresp(2, DIAMETRE/lambda, DIAMETRE/lambda, lambda);
+    int typbeam=(fggaussian)?1:2;
+    Four2DResponse fresp(typbeam, DIAMETRE/lambda, DIAMETRE/lambda, lambda);
     Four2DResponse* fresp_p=&fresp;
     Four2DRespTable resptbl;
@@ -119 +156 @@
       cout << "applobe[2.b]: initializing Four2DRespTable from file" << resptblname << endl;
       resptbl.readFromPPF(resptblname);
+      resptbl.renormalize(1.);
       fresp_p=&resptbl;
     }
@@ -126 +164 @@
 
     if (fgcorrbeam) {
-      double DoL = tbeamDoL;
+      double DoL = tbeamDiam/lambda;
       double tbeamarcmin = RadianToDegree(1.22/DoL)*60.;
-      int typcb = 2;
+      int typcb = (fggaussian)?1:2;
       //      if (fgresptbl) typcb=22;
       Four2DResponse tbeam(typcb, DoL, DoL );
-      cout << "applobe[3]: calling Correct2RefLobe() with target beam D/Lambda=" << DoL  
-           << " -> arcmin " << tbeamarcmin << " TypDishResp=" << typcb << endl;
+      cout << "applobe[3]: calling Correct2RefLobe() with target beam Diameter=" << tbeamDiam 
+           << " D/Lambda=" << DoL << " -> arcmin " << tbeamarcmin << " TypDishResp=" << typcb << endl;
       beam.Correct2RefLobe(tbeam, incube, dx, dy, Freq0MHz, dfreq);
+    }
+    else if (fixedbeam) {
+      cout << "applobe[3]: calling ApplyLobe() (freq.independent beam) ... " << endl;
+      beam.ApplyLobe(incube, dx, dy, Freq0MHz);
     }
     else {