Changeset 2058 in Sophya for trunk/ArchTOIPipe/TestPipes

Timestamp:

Jun 18, 2002, 2:21:09 PM (23 years ago)

Author:

ansari

Message:

Ajout methode FITSTOIReader::setFlagFile() pour fichiers flag separes
de Level2 - Reza 18/6/2002

Location:

trunk/ArchTOIPipe/TestPipes

Files:

• toistat.cc (modified) (5 diffs)
• tsttoi2map.cc (modified) (7 diffs)

trunk/ArchTOIPipe/TestPipes/toistat.cc

@@ -45 +45 @@
     cout << "\n Usage : toistat [-intoi toiname] [-start snb] [-end sne] \n" 
          << "         [-wtoi sz] [-wclean wsz,nbw] [-range min,max] [-cleannsig nsig] \n"
-         << "         [-outppf] [-noprstat] [-useseqbuff] \n"
+         << "         [-sepflg sepFlagFile] [-outppf] [-noprstat] [-useseqbuff] \n"
          << "         inFitsName outFileName \n" 
          << "   -start snb : sets the start sample num \n"
@@ -52 +52 @@
          << "              default= -16000,16000\n"
          << "   -intoi toiName : select input TOI name (def bolo)\n"
+         << "   -sepflg sepFlagFileName: sets separate flag file (Level2)\n" 
          << "   -wtoi sz : sets TOISeqBuff buffer size (def= 8192)\n"
          << "   -wclean wsz,nbw : sets cleaner window parameters (256,5) \n"
@@ -90 +91 @@
   int clean_nbw = 5;
   double clean_nsig = 999999.;
+
+// Fichier de flag separe / Level2 / Reza 18/6/2002
+  string sepflagfile;   
+  bool sepflg = false;  
 
   // File names 
@@ -130 +135 @@
       intoi = arg[ia+1]; ia++;
     }    
+    else if (strcmp(arg[ia],"-sepflg") == 0) {
+      if (ia == narg-1) Usage(true);  
+      sepflagfile = arg[ia+1]; 
+      sepflg = true; ia++;
+    }    
     else if (strcmp(arg[ia],"-outppf") == 0)  fgoutppf = true;
 
@@ -166 +176 @@
     cout << "> Creating FITSTOIReader object - InFile=" << infile << endl;
     FITSTOIReader r(infile); 
+    if (sepflg) {
+      cout << " Setting separate flag file for InTOI_bolo File=" << sepflagfile 
+           << " (Flags=FlgToiSpike, FlgToiSource)" << endl;
+      vector<FlagToiDef> flgcol;
+      flgcol.push_back(FlgToiSpike);
+      flgcol.push_back(FlgToiSource);
+      r.setFlagFile(sepflagfile, flgcol);
+    }
 
     cout << "> Creating SimpleCleaner() " << endl;

trunk/ArchTOIPipe/TestPipes/tsttoi2map.cc

@@ -27 +27 @@
      <<"           [-a label_coord1] [-d label_coord2] [-b label_bolomuv]"<<endl
      <<"           [-n nlat] [-i c,h] [-o c,h]"<<endl
-     <<"           [-m vmin] [-M vmax] [-f flag]"<<endl
+     <<"           [-m vmin] [-M vmax] [-f flag] [-F sepFlagFileName]"<<endl
      <<"           fitsin_point fitsin_bolo fitsphout [fitsphwout]"<<endl
      <<" -p lp : print level (def=0)"<<endl
@@ -44 +44 @@
      <<" -f flag : samples are bad if match flag"<<endl
      <<" -N nbfiles : Number of fitsin_point files fitsin_point%d.fits 0..nb-1"<<endl
+     <<" -c calibcoeff : Coefficient de calibration a appliquer (def=1.)"<<endl
      <<" -I : sampleNum are implicit in fits files (def=no)"<<endl
+     <<" -F sepFlagFileName : separate flag file name for levelS"<<endl
+     <<"                      (FlagColName: Glitch PtSrc)"<<endl
      <<" fitsin_point : fits file for pointing"<<endl
      <<" fitsin_bolo : fits file for bolo values"<<endl
@@ -70 +73 @@
 bool fgprstat = true;
 int nbpointfiles = 0;
+double coeffcalib = 1.;
+string sepflagfile;   // Fichier de flag separe / Level2 / Reza 18/6/2002
+bool sepflg = false;  //   "     "     "      "       "       "       "
+
 int c;
-while((c = getopt(narg,arg,"hIp:s:w:a:d:b:n:i:o:m:M:f:e:N:")) != -1) {
+while((c = getopt(narg,arg,"hIp:s:w:a:d:b:n:i:o:m:M:f:e:N:c:F:")) != -1) {
   switch (c) {
   case 's' :
@@ -125 +132 @@
     nbpointfiles = atoi(optarg);
     break;
+  case 'c' :
+    coeffcalib = atof(optarg);
+    break;
   case 'I' :
     snimplicit = true;
+    break;
+  case 'F' :
+    sepflagfile = optarg;
+    sepflg = true;
     break;
   case 'h' :
@@ -154 +168 @@
     <<"  ...... ctype="<<tcoormap<<endl;
 cout<<"Equinoxe "<<equi<<" years"<<endl;
-
+cout<<"CoeffCalib "<<coeffcalib<<endl;
+ 
 SophyaInit();
 InitTim();
@@ -170 +185 @@
  if(ncolb<1) exit(-4);
 
+
+ sepflagfile = optarg;
+ sepflg = true;
+ if (sepflg) {
+   cout << " Setting separate flag file for InTOI_bolo File=" << sepflagfile 
+        << " (Flags=FlgToiSpike, FlgToiSource)" << endl;
+   vector<FlagToiDef> flgcol;
+   flgcol.push_back(FlgToiSpike);
+   flgcol.push_back(FlgToiSource);
+   rfitsb.setFlagFile(sepflagfile, flgcol);
+ }
  string pointfileI = fitsin_point;
  MuTyV numf=0;
@@ -212 +238 @@
  toi2m.SetTestMin(tmin,vmin);
  toi2m.SetTestMax(tmax,vmax);
+ toi2m.SetCalibrationFactor(coeffcalib);
  toi2m.Print(cout);