Changeset 2050 in Sophya for trunk/ArchTOIPipe/TestPipes

Timestamp:

Jun 10, 2002, 4:23:56 PM (23 years ago)

Author:

ansari

Message:

Modifs lecture multiple fichiers de pointage - Reza 10/6/2002

Location:

trunk/ArchTOIPipe/TestPipes

Files:

• aksj02.cc (modified) (4 diffs)
• tsttoi2map.cc (modified) (5 diffs)

trunk/ArchTOIPipe/TestPipes/aksj02.cc

@@ -44 +44 @@
          << "         [-soff2 wsz,nptfit,degpol] [-soff2nt PPFName] \n"
          << "         [-bgal toiname] [-bgalfile pointFitsName] [-bgalnf N]\n"
-         << "         [-bgal toiname] [-bgalfile pointFitsName] \n"
          << "         [-bgalcut bmin,bmax]  [-bgal2cut bmin,bmax] \n"
          << "         [-gfilt wsz,sigma] \n" 
@@ -74 +73 @@
          << "   -wrtms : Write mean/sigma TOI's from cleaner \n"
          << "   -wrtifo : Write incopie flagged  TOI's from cleaner/ offset from SimpleOffsetEstimator \n"
-         << "   -prstat : PrintStat with ProcSampleCounter \n"
+         << "   -noprstat : DO NOT PrintStat with ProcSampleCounter \n"
          << "   -useseqbuff : Use TOISeqBuffered TOI's (default: TOISegmented) \n"
          << endl;
@@ -90 +89 @@
   bool fgsetstart = false;
 
-  bool fgprstat = false;
+  bool fgprstat = true;
   bool fgsegmented = true;
 
@@ -250 +249 @@
 
     else if (strcmp(arg[ia],"-prstat") == 0)  fgprstat = true;
+    else if (strcmp(arg[ia],"-noprstat") == 0)  fgprstat = false;
     else if (strcmp(arg[ia],"-useseqbuff") == 0)  fgsegmented = false;
 

trunk/ArchTOIPipe/TestPipes/tsttoi2map.cc

@@ -43 +43 @@
      <<" -M vmax : samples are good if sample value <= vmax"<<endl
      <<" -f flag : samples are bad if match flag"<<endl
+     <<" -N nbfiles : Number of fitsin_point files fitsin_point%d.fits 0..nb-1"<<endl
      <<" -I : sampleNum are implicit in fits files (def=no)"<<endl
      <<" fitsin_point : fits file for pointing"<<endl
@@ -67 +68 @@
 long sdeb,sfin;
 bool snimplicit = false;
+bool fgprstat = true;
+int nbpointfiles = 0;
 int c;
-while((c = getopt(narg,arg,"hIp:s:w:a:d:b:n:i:o:m:M:f:e:")) != -1) {
+while((c = getopt(narg,arg,"hIp:s:w:a:d:b:n:i:o:m:M:f:e:N:")) != -1) {
   switch (c) {
   case 's' :
@@ -118 +121 @@
     sscanf(optarg,"%ul",&badflg);
     tflg = true;
+    break;
+  case 'N' :
+    nbpointfiles = atoi(optarg);
     break;
   case 'I' :
@@ -139 +145 @@
     <<"Fits Infile Bolo "<<fitsin_bolo<<endl
     <<"  ...label_bolomuv "<<label_bolomuv<<endl;
-cout<<"Fits Infile Pointing "<<fitsin_point<<endl
+cout<<"Fits Infile Pointing "<<fitsin_point<<" NbFiles=" <<nbpointfiles<<endl
     <<"  ...label_coord1 "<<label_coord1<<endl
     <<"  ...label_coord2 "<<label_coord2<<endl
@@ -164 +170 @@
  if(ncolb<1) exit(-4);
 
- FITSTOIReader rfitsp(fitsin_point);
+ string pointfileI = fitsin_point;
+ MuTyV numf=0;
+ if(nbpointfiles > 0) 
+   pointfileI += (string)numf + ".fits";
+ cout << " Creating FITSTOIReader for pointing file " << pointfileI << endl;
+ FITSTOIReader rfitsp(pointfileI);
+ if(nbpointfiles > 0) 
+   for(int kff=1; kff<nbpointfiles; kff++) {
+     numf=kff;
+     pointfileI = fitsin_point;
+     pointfileI += (string)numf + ".fits";
+     cout << " pointing FITSTOIReader.addFile( " << pointfileI << " )" << endl;
+     rfitsp.addFile( pointfileI );
+   }
+ 
  if(snimplicit) rfitsp.setImplicitSN();
  int ncolp = rfitsp.getNOut();