Changeset 4022 in Sophya for trunk/Cosmo/RadioBeam/calcpk2.cc

Timestamp:

Sep 28, 2011, 5:13:51 PM (14 years ago)

Author:

ansari

Message:

modifs pour pouvoir imposer la moyenne en temp des plans X,Y des cubes lors de l'extraction du signal HI, Reza 28/9/2011

File:

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

trunk/Cosmo/RadioBeam/calcpk2.cc

@@ -44 +44 @@
 {
   if ( (narg<6)||((narg>1)&&(strcmp(arg[1],"-h")==0)) ) {
-    cout << " Usage: [-t -g -mxr val] calcpk2 InMapLSS convFacLSS InMapFgnd convFacFgnd OutPkFile \n" 
+    cout << " Usage: [-t -g -mxr val -fmt T0,alpha] calcpk2 InMapLSS convFacLSS InMapFgnd convFacFgnd OutPkFile \n" 
          << "        [PixNoiseLevel] [D_Dish/Four2DRespTableFile CorBeamDiam] \n"
          << "        [NSigSrcThr] [P2/P1] [RecMapFile] " << endl;
     if ((narg>1)&&(strcmp(arg[1],"-h")==0)) {
       cout << "-t -g : Triangular / gaussian beam shape (def=gaussian) \n" 
+           << "-fmt T0,alpha: fix mean slice temp. according to power law T0,alpha \n "
            << "-mxr val: Max beam correction factor (default=10.) \n "
            << "- InMapLSS: Input 3D LSS cube (PPF file name) \n "
@@ -79 +80 @@
     double maxratio=10.;   // valeur max du rapport des lobes lors de la correction de lobe
 
+    bool fgfix_mXYtemp=false;
+    double T0_pl=1.; 
+    double alpha_pl=0.;
     // decodage argument optionnel 
     bool fgoptarg=true;
@@ -86 +90 @@
       else if (fbo=="-g")  { fggaussian=true; arg++; narg--; }
       else if (fbo=="-mxr")  { arg++; maxratio=atof(arg[1]); arg++; narg-=2; }
+      else if (fbo=="-fmt")  
+        { fgfix_mXYtemp=true; arg++; sscanf(arg[1],"%lg,%lg",&T0_pl,&alpha_pl); arg++; narg-=2; }
       else fgoptarg=false;
     }
@@ -202 +208 @@
 
     ForegroundCleaner  cleaner(*arep_p, tbeam, skycube, maxratio); 
+    if (fgfix_mXYtemp) {
+      cout << "calcpk2[3.b] : calling cleaner.FixMeanXYTemp(T0=" << T0_pl << ",alpha=" << alpha_pl << ")" << endl;
+      cleaner.FixMeanXYTemp(T0_pl,alpha_pl);
+    }
     if (fgcorrbeam) {
-      cout << "calcpk2[3.b] : calling cleaner.BeamCorrections() for target beam Diameter=" << tbeamDiam  
+      cout << "calcpk2[3.c] : calling cleaner.BeamCorrections() for target beam Diameter=" << tbeamDiam  
            << " D/Lambda=" << DoL  << "  -> arcmin " << tbeamarcmin << " TypDishResp=" << typcb << endl;
       cleaner.BeamCorrections();
     }
-    cout << " calcpk2[3.c] : calling cleaner.CleanNegatives() ... " << endl;
+    cout << " calcpk2[3.d] : calling cleaner.CleanNegatives() ... " << endl;
     cleaner.CleanNegatives();
     if (fgclnsrc) {
-      cout << "calcpk2[3.d] : calling cleaner.CleanPointSources() with threshold NSigma=" << nsigsrc  << endl;
+      cout << "calcpk2[3.e] : calling cleaner.CleanPointSources() with threshold NSigma=" << nsigsrc  << endl;
       cleaner.CleanPointSources(nsigsrc);
     }