Changeset 3796 in Sophya for trunk/Cosmo/RadioBeam/repicon.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/repicon.cc (modified) (5 diffs)

trunk/Cosmo/RadioBeam/repicon.cc

@@ -36 +36 @@
 {
   if (((narg>1)&&(strcmp(arg[1],"-h")==0))||(narg<3))  {
-    cout << " Usage: repicon configId OutPPFName [Lambda] " 
+    cout << " Usage: repicon configId OutPPFName [Lambda=0.357] [RenormalizeMax]" 
          << " configs: f4x4 , f8x8 , f20x20 Filled array of nxn dishes \n" 
          << "          confA , confB, confC : semi-filled array of dishes \n"
@@ -54 +54 @@
   if (outfile==".")  outfile = "repicon.ppf";
   double LAMBDA=0.357 ;  // 21 cm at z=0.7 
-  if (narg>3) LAMBDA = atof(arg[2]);
-
+  if (narg>3) LAMBDA = atof(arg[3]);
+  bool fgrenorm=false;
+  double rmax=1.;
+  if (narg>4) {
+    rmax=atof(arg[4]);
+    fgrenorm=true;
+  }
   //-- end command line arguments
   
@@ -61 +66 @@
   try {  // exception handling try bloc at top level
 
-    cout << " 2.c/ Creating MultiDish Filled Array " << endl;
     double Ddish = 5.;
     double Ddish2 = 7.5;
@@ -122 +126 @@
     double LMAX = D;
     bool fgnoauto = true;
-    int NRX=100;
-    int NRY=100;
+    int NRX=200;
+    int NRY=200;
 
     MultiDish mdish(LAMBDA, LMAX, vdishes, fgnoauto);
@@ -134 +138 @@
     PrtTim("Apres mdish.GetResponse()");
 
+    Four2DRespTable mdresp(hrep, Dol, LAMBDA);
+    if (fgrenorm) {
+      cout << " repicon[2.b] call to mdresp.renormalize(" << rmax << ")"; 
+      double omax=mdresp.renormalize(rmax);
+      cout << " Old Max=" << omax << endl;
+    }
     cout << " repicon[3] : saving Four2DRespTable for config " << config << " to PPF file " << outfile << endl;
-    Four2DRespTable mdresp(hrep, Dol, LAMBDA);
     mdresp.writeToPPF(outfile);