Changeset 2050 in Sophya for trunk/ArchTOIPipe/TestPipes/tsttoi2map.cc

Timestamp:

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

Author:

ansari

Message:

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

File:

• trunk/ArchTOIPipe/TestPipes/tsttoi2map.cc (modified) (5 diffs)

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();