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

Timestamp:

Jun 29, 2010, 7:26:01 PM (15 years ago)

Author:

ansari

Message:

Suite modif utilisation de reponse calcule (tabulee) de l'interferometre pour soustraction avant-plans, Reza 29/06/2010

File:

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

trunk/Cosmo/RadioBeam/applobe.cc

@@ -3 +3 @@
     R. Ansari , C. Magneville - Juin 2010 
 
-  Usage: applobe In3DPPFName Out3DPPFName [ResmapleFactor=0.5,0.333...] [TargetBeamDoL]
+  Usage: applobe Diameter/Four2DRespTableFile In3DPPFName Out3DPPFName 
+                 [TargetBeamDoL] [ResmapleFactor=0.5,0.333...] 
 ---------------------------------------------------------------  */
 
@@ -50 +51 @@
 
   if (narg < 3) {
-    cout << "Usage: applobe In3DPPFName Out3DPPFName [ResmapleFactor=0.5,0.333...] [TargetBeamDoL]\n" << endl;
+    cout << "Usage: applobe Diameter/Four2DRespTableFile In3DPPFName Out3DPPFName " 
+         << "             [TargetBeamDoL] [ResmapleFactor=0.5,0.333...] \n" << endl;
     return 1;
   }
 
   // decodage arguments 
-  string outname = arg[2];
-  string inname = arg[1];
-  bool fgresample=false;
-  double fresamp=1.;
-  if (narg>3) {
-    fresamp=atof(arg[3]);
-    if ((fabs(fresamp)>1.e-2)&&(fabs(fresamp-1.)>1.e-2))  fgresample=true; 
+  bool fgresptbl=true;
+  double DIAMETRE=100.;
+  string resptblname;
+  if (isdigit(*arg[1])) {
+    fgresptbl=false;
+    DIAMETRE=atof(arg[1]);
   }
+  else resptblname=arg[1];
+
+  string inname = arg[2];
+  string outname = arg[3];
   bool fgcorrbeam=false;
   double tbeamDoL=135;
@@ -68 +73 @@
     tbeamDoL=atof(arg[4]);
     if (tbeamDoL>1)  fgcorrbeam=true;
+  }
+  bool fgresample=false;
+  double fresamp=1.;
+  if (narg>5) {
+    fresamp=atof(arg[5]);
+    if ((fabs(fresamp)>1.e-2)&&(fabs(fresamp-1.)>1.e-2))  fgresample=true; 
   }
 
@@ -101 +112 @@
     conv.setFrequency(Freq0MHz);
     double lambda = conv.getLambda();
-    Four2DResponse fresp(2, InterfArrayDiametre/lambda, InterfArrayDiametre/lambda, lambda);    
-    BeamEffect beam(fresp);
+    
+    Four2DResponse fresp(2, DIAMETRE/lambda, DIAMETRE/lambda, lambda);
+    Four2DResponse* fresp_p=&fresp;
+    Four2DRespTable resptbl;
+    if (fgresptbl) {
+      cout << "applobe[2.b]: initializing Four2DRespTable from file" << resptblname << endl;
+      resptbl.readFromPPF(resptblname);
+      fresp_p=&resptbl;
+    }
+    else cout << " applobe[2.b]: Four2DResponse ( Diameter=" << DIAMETRE << " Lambda= " << lambda
+              << " DoL=" << DIAMETRE/lambda << " ) " << endl;
+    BeamEffect beam(*fresp_p);
 
     if (fgcorrbeam) {
       double DoL = tbeamDoL;
       double tbeamarcmin = RadianToDegree(1.22/DoL)*60.;
-      Four2DResponse tbeam(2, DoL, DoL );
+      int typcb = 2;
+      if (fgresptbl) typcb=22;
+      Four2DResponse tbeam(typcb, DoL, DoL );
       cout << "applobe[3]: calling Correct2RefLobe() with target beam D/Lambda=" << DoL  
-           << " -> arcmin " << tbeamarcmin << endl;
+           << " -> arcmin " << tbeamarcmin << " TypDishResp=" << typcb << endl;
       beam.Correct2RefLobe(tbeam, incube, dx, dy, Freq0MHz, dfreq);
     }