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

trunk/Cosmo/RadioBeam/calcpk2.cc

@@ -7 +7 @@
     R. Ansari , C. Magneville - Juin 2010 
 
-Usage: calcpk2 InMapLSS convFacLSS InMapSync convFacSync InMapRadioSource convFacRsc OutPk 
-               [PixNoiseLevel] [TargetBeamArcmin] [NSigSrcThr]
+Usage: calcpk2 InMapLSS convFacLSS InMapSync convFacSync InMapRadioSource convFacRsc OutPkFile 
+               [PixNoiseLevel] [Diameter/Four2DRespTableFile] [TargetBeamArcmin] [NSigSrcThr]
 ---------------------------------------------------------------  */
 
@@ -44 +44 @@
 {
   if (narg<6) {
-    cout << " Usage: calcpk2 InMapLSS convFacLSS InMapSync convFacSync  OutPk \n" 
-         << "        [PixNoiseLevel] [TargetBeamArcmin] [NSigSrcThr]" << endl;
+    cout << " Usage: calcpk2 InMapLSS convFacLSS InMapSync convFacSync  OutPkFile \n" 
+         << "        [PixNoiseLevel] [Diameter/Four2DRespTableFile] [TargetBeamArcmin] [NSigSrcThr]" << endl;
     return 1;
   }
@@ -65 +65 @@
     
     bool fgcorrbeam=true;
+
+    bool fgresptbl=false;
+    double DIAMETRE=100.;
+    string resptblname;
+    if (narg>7) {
+      if (isdigit(*arg[7])) {
+        fgresptbl=false;
+        DIAMETRE=atof(arg[7]);
+      }
+      else { 
+        resptblname=arg[7];
+        fgresptbl=true;
+      }
+    }
     double tbeamDoL=135;
-    if (narg>7) {
-      tbeamDoL=atof(arg[7]);
+    if (narg>8) {
+      tbeamDoL=atof(arg[8]);
       if (tbeamDoL<1.)  fgcorrbeam=false;
     }
     bool fgclnsrc=true;
     double nsigsrc=5.;
-    if (narg>8) {
-      nsigsrc=atof(arg[8]);
+    if (narg>9) {
+      nsigsrc=atof(arg[9]);
       if (nsigsrc<1.e-6)  fgclnsrc=false;
     }
@@ -119 +133 @@
     conv.setFrequency(Freq0MHz);
     double lambda = conv.getLambda();
-    Four2DResponse arep(2, InterfArrayDiametre/lambda, InterfArrayDiametre/lambda, lambda);
+    Four2DResponse arep(2, DIAMETRE/lambda, DIAMETRE/lambda, lambda);
+    Four2DResponse* arep_p=&arep;
+    Four2DRespTable resptbl;
+    if (fgresptbl) {
+      cout << "calcpk2[3.a]: initializing Four2DRespTable from file" << resptblname << endl;
+      resptbl.readFromPPF(resptblname);
+      arep_p=&resptbl;
+    }
+    else cout << " calcpk2[3.a]: Four2DResponse ( Diameter=" << DIAMETRE << " Lambda= " << lambda
+              << " DoL=" << DIAMETRE/lambda << " ) " << endl;
+
     double DoL = tbeamDoL;
     double tbeamarcmin = RadianToDegree(1.22/DoL)*60.;
-    Four2DResponse tbeam(2, DoL, DoL );
-    ForegroundCleaner  cleaner(arep, tbeam, skycube); 
+    int typcb = 2;
+    //    if (fgresptbl) typcb=22;
+    Four2DResponse tbeam(typcb, DoL, DoL );
+
+    ForegroundCleaner  cleaner(*arep_p, tbeam, skycube); 
     if (fgcorrbeam) {
-      cout << "calcpk2[3.a] : calling cleaner.BeamCorrections() for target beam D/Lambda=" << DoL 
-           << "  -> arcmin " << tbeamarcmin << endl;
+      cout << "calcpk2[3.b] : calling cleaner.BeamCorrections() for target beam D/Lambda=" << DoL 
+           << "  -> arcmin " << tbeamarcmin << " TypDishResp=" << typcb << endl;
       cleaner.BeamCorrections();
     }
-    cout << " calcpk2[3.b] : calling cleaner.CleanNegatives() ... " << endl;
+    cout << " calcpk2[3.c] : calling cleaner.CleanNegatives() ... " << endl;
     cleaner.CleanNegatives();
     if (fgclnsrc) {
-      cout << "calcpk2[3.c] : calling cleaner.CleanPointSources() with threshold NSigma=" << nsigsrc  << endl;
+      cout << "calcpk2[3.d] : calling cleaner.CleanPointSources() with threshold NSigma=" << nsigsrc  << endl;
       cleaner.CleanPointSources(nsigsrc);
     }