Changeset 3933 in Sophya for trunk/Cosmo/RadioBeam/repicon.cc

Timestamp:

Dec 24, 2010, 1:22:04 AM (15 years ago)

Author:

ansari

Message:

correction/amelioration calcul reponse multi-dish ds plan (u,v), Reza 24/12/2010

File:

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

trunk/Cosmo/RadioBeam/repicon.cc

@@ -52 +52 @@
        << "    -D DishDiameter (default=5 m) \n"
        << "    -lmax array extension (default=100 m ) for response calculation kmax \n" 
+       << "    -repnxy Nx,Ny (default=200,200) ResponseTable binning \n"
+       << "    -beamnxy Nx,Ny (default=200,200) Beam sampling \n"
        << "    -rot ThetaMaxDeg,NTheta,PhiMaxDeg,NPhi (default NO-rotate/pointing -> 23,10,30,15 ) \n" 
        << "    -prt PrtLev,PrtModulo (default=0,10) \n" 
@@ -87 +89 @@
   int nphi=15;
   bool fgpoint=false;
-
+  int NRX=200;
+  int NRY=200;
+  int NBX=200;
+  int NBY=200;
+  
   int ka=1;
   while (ka<(narg-1)) {
@@ -101 +107 @@
     else if (strcmp(arg[ka],"-lmax")==0) {
       LMAX=atof(arg[ka+1]);  fgLMAXfixed=true;  ka+=2;
+    }
+    else if (strcmp(arg[ka],"-repnxy")==0) {
+      sscanf(arg[ka+1],"%d,%d",&NRX,&NRY);   ka+=2;
+    }
+    else if (strcmp(arg[ka],"-beamnxy")==0) {
+      sscanf(arg[ka+1],"%d,%d",&NBX,&NBY);   ka+=2;
     }
     else if (strcmp(arg[ka],"-rot")==0) {
@@ -124 +136 @@
   try {  // exception handling try bloc at top level
     cout << " ==== repicon.cc program , (u,v) plane response ==== " << endl;
+
+    H21Conversions conv;
+    double LAMBDA=0.357 ;  // 21 cm at z=0.7 
+    conv.setRedshift(z_Redshift);
+    LAMBDA = conv.getLambda();
     
     double Eta=0.95;
     double cylW=12.;   // Largeur des cylindres 
-    double cylRL=0.5;  // Longeur des elements de reception le long du cylindre
+    double cylRL=2*LAMBDA;  // Longeur des elements de reception le long du cylindre
     double etaW=0.95;   // Efficacite de couverture en largeur
     double etaRL=0.9;   // Efficacite de couverture le long du cylindre
 
-    int cnt=0;
-
     vector<Dish> vdishes;
     cout << " repicon[1] : creating dish vector/ MultiDish for configuration: " << config << endl;
@@ -146 +161 @@
     }
     else if (config=="f3cyl") {
-      cylW=10.;   cylRL=0.5; 
-      vdishes=CreateFilledCylConfig(3, 128, cylW, cylRL, etaW, etaRL, true);
+      cylW=10.;   cylRL=2*LAMBDA; 
+      vdishes=CreateFilledCylConfig(3, 92, cylW, cylRL, etaW, etaRL, true);
     }
     else if (config=="f3cylp") {
-      cylW=10.;   cylRL=0.5; 
-      vdishes=CreateFilledCylConfig(3, 128, cylW, cylRL, etaW, etaRL, false);
+      cylW=10.;   cylRL=2*LAMBDA; 
+      vdishes=CreateFilledCylConfig(3, 92, cylW, cylRL, etaW, etaRL, false);
     }
     else if (config=="f8cyl") {
-      cylW=12.;   cylRL=0.5; 
-      vdishes=CreateFilledCylConfig(8, 192, cylW, cylRL, etaW, etaRL, true);
+      cylW=12.;   cylRL=2*LAMBDA; 
+      vdishes=CreateFilledCylConfig(8, 160, cylW, cylRL, etaW, etaRL, true);
     }
     else if (config=="f8cylp") {
-      cylW=12.;   cylRL=0.5; 
-      vdishes=CreateFilledCylConfig(8, 192, cylW, cylRL, etaW, etaRL, false);
+      cylW=12.;   cylRL=2*LAMBDA; 
+      vdishes=CreateFilledCylConfig(8, 160, cylW, cylRL, etaW, etaRL, false);
     }
     else if (config=="confA") {
@@ -190 +205 @@
       return 99;
     }
-    
-    H21Conversions conv;
-    double LAMBDA=0.357 ;  // 21 cm at z=0.7 
-    conv.setRedshift(z_Redshift);
-    LAMBDA = conv.getLambda();
-    
+        
     double Dol = LMAX/LAMBDA;
     bool fgnoauto = true;
-    int NRX=200;
-    int NRY=200;
     
     cout << " repicon[1] : Lambda=" << LAMBDA << " LMAX= " << LMAX << " NDishes=" << vdishes.size() 
@@ -207 +215 @@
     mdish.SetPrtLevel(prtlev,prtmod);
     mdish.SetRespHisNBins(NRX,NRY); 
-
+    mdish.SetBeamNSamples(NBX,NBY);
     if (fgpoint)  {
       cout << " repicon[1.b] : activating pointing , ThetaMaxDeg=" << thetamxdeg << " NTheta=" << nteta
-           << "  PhiMaxDeg= " << phimxdeg << " NPhi" << nphi << endl;
+           << "  PhiMaxDeg= " << phimxdeg << " NPhi=" << nphi << endl;
       mdish.SetThetaPhiRange(Angle(thetamxdeg,Angle::Degree).ToRadian(),nteta,Angle(phimxdeg,Angle::Degree).ToRadian(), nphi);
     }