Changeset 1498 in Sophya for trunk/ArchTOIPipe/TestPipes

Timestamp:

May 17, 2001, 5:43:42 PM (24 years ago)

Author:

ansari

Message:

Commit de la modif toi2map.cc .h par Cecile - Reza 17/5/2001

File:

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

trunk/ArchTOIPipe/TestPipes/tsttoi2map.cc

@@ -1 +1 @@
 #include <unistd.h>
 #include <stdexcept>
-
+#include <stdlib.h>
+#include <stdio.h>
 #include "toi.h"
 #include "toiprocessor.h"
@@ -17 +18 @@
 void usage(void) {
  cout<<"tsttoi2map [-h] [-p lp] [-s samplemin,samplemax] [-w data_window_size]"<<endl
-     <<"           [-a label_alpha] [-d label_delta] [-b label_bolomuv]"<<endl
-     <<"           [-n nlat] [-G] fitsin fitsphout fitsphwout"<<endl;
+     <<"           [-a label_coord1] [-d label_coord2] [-b label_bolomuv]"<<endl
+     <<"           [-n nlat] coord_ini coord_fin fitsin_bolo fitsin_point fitsphout fitsphwout"<<endl
+     <<" coord_ini[_fin] = G for Galactic coordinates"<<endl
+     <<"                 = Q for equatorial coordinates"<<endl
+     <<" coord1 = alpha or longitude ; coord2 = delta or latitude"<<endl;
  return;
 }
@@ -30 +34 @@
 int lp = 0, width = 8192;
 int nlat = 128;   // npixel = 12 * nlat^2
-char *label_alpha = "alpha", *label_delta = "delta", *label_bolomuv = "boloMuV";
-bool coorgal=false;
+char *label_coord1 = "coord1", *label_coord2 = "coord2", *label_bolomuv = "boloMuV";
 long sdeb,sfin;
 int c;
@@ -52 +55 @@
     break;
   case 'a' :
-    label_alpha = optarg;
+    label_coord1 = optarg;
     break;
   case 'd' :
-    label_delta = optarg;
+    label_coord2 = optarg;
     break;
   case 'b' :
@@ -64 +67 @@
     if(nlat<0) nlat=128;
     break;
-  case 'G' :
-    coorgal = true;
-    break;
   case 'h' :
     usage(); exit(-1);
@@ -74 +74 @@
   }
 }
-if(optind+1>=narg) {usage(); exit(-2);}
-char * fitsin          = arg[optind];
-string const fitsphout = arg[optind+1];
+if(optind+3>=narg) {usage(); exit(-2);}
+bool mcoorgal=true;
+bool fcoorgal=true;
+ 
+if ( strcmp(arg[optind],"Q") == 0) mcoorgal=false; 
+if ( strcmp(arg[optind+1],"Q") == 0) fcoorgal=false; 
+ 
+char * fitsin_bolo = arg[optind+2];
+char * fitsin_point = arg[optind+3];
+string const fitsphout = arg[optind+4];
 string fitsphwout = "";
-if(optind+2<narg) fitsphwout = arg[optind+2];
+if(optind+5<narg) fitsphwout = arg[optind+5];
 
 cout<<">>>> tsttoi2map:"<<endl
-    <<"Fits Infile(snum,alpha,delta,bolomuv)= "<<fitsin<<endl
-    <<"  ...label_alpha "<<label_alpha<<"  ,  label_delta "<<label_delta<<endl
+    <<"Fits Infile(snum,bolomuv)= "<<fitsin_bolo<<endl
+    <<"Fits Infile(snum,coord1,coord2)= "<<fitsin_point<<endl
+    <<"  ...label_coord1 "<<label_coord1<<"  ,  label_coord2 "<<label_coord2<<endl
+    <<"  with coordinates Gal "<<mcoorgal<<endl
     <<"  ...label_bolomuv "<<label_bolomuv<<endl
     <<"Fits Sphere Healpix"<<fitsphout<<endl
     <<"  ...nlat "<<nlat<<endl
+    <<"  with coordinates Gal "<<fcoorgal<<endl
     <<"Fits Sphere Healpix Error "<<fitsphwout<<endl;
 
@@ -96 +106 @@
 
  // FITS reader
- FITSTOIReader rfits(fitsin);
- int ncol = rfits.getNOut();
- cout<<"Number of columns in fits Infile : "<<ncol<<endl;
- if(ncol<3) exit(-3);
+ FITSTOIReader rfitsb(fitsin_bolo);
+ int ncolb = rfitsb.getNOut();
+ cout<<"Number of columns in fits Infile_bolo : "<<ncolb<<endl;
+ if(ncolb<1) exit(-3);
+
+ FITSTOIReader rfitsp(fitsin_point);
+ int ncolp = rfitsp.getNOut();
+ cout<<"Number of columns in fits Infile_point : "<<ncolp<<endl;
+ if(ncolp<2) exit(-3);
 
  // Creation de la sphere Healpix
@@ -117 +132 @@
  TOI2Map toi2m(sph,wsph);
  cout<<"TOI2Map created"<<endl;
- toi2m.SetCoorGal(coorgal,2000.);  // equinoxe de ref. 2000. (pour archtoi)
+ toi2m.SetCoorGal(mcoorgal,fcoorgal,2000.);  // equinoxe de ref. 2000. (pour archtoi)
 
  // Definition des tuyaux
- TOISeqBuffered * toialphain = new TOISeqBuffered("toi_alpha_in",width);
- if(lp) toialphain->setDebugLevel(1);
- rfits.addOutput(label_alpha,toialphain);
- toi2m.addInput("AlphaIn",toialphain);
-
- TOISeqBuffered * toideltain = new TOISeqBuffered("toi_delta_in",width);
- if(lp) toideltain->setDebugLevel(1);
- rfits.addOutput(label_delta,toideltain);
- toi2m.addInput("DeltaIn",toideltain);
-
+ TOISeqBuffered * toicoord1in = new TOISeqBuffered("toi_coord1_in",width);
+ if(lp) toicoord1in->setDebugLevel(1);
+ rfitsp.addOutput(label_coord1,toicoord1in);
+ toi2m.addInput("Coord1In",toicoord1in);
+
+ TOISeqBuffered * toicoord2in = new TOISeqBuffered("toi_coord2_in",width);
+ if(lp) toicoord2in->setDebugLevel(1);
+ rfitsp.addOutput(label_coord2,toicoord2in);
+ toi2m.addInput("Coord2In",toicoord2in);
+ 
  TOISeqBuffered * toibolin = new TOISeqBuffered("toi_bolo_in",width);
  if(lp) toibolin->setDebugLevel(1);
- rfits.addOutput(label_bolomuv,toibolin);
+ rfitsb.addOutput(label_bolomuv,toibolin);
  toi2m.addInput("BoloIn",toibolin);
 
  // Run
  cout<<"----- FITSReaderTOI::PrintStatus() : -----"<<endl;
- rfits.PrintStatus(cout);
+ rfitsp.PrintStatus(cout);
+ rfitsb.PrintStatus(cout);
 
  PrtTim("starting threads");
- rfits.start();
+ rfitsp.start();
+ rfitsb.start();
  toi2m.start();
 
  if(lp>1)
    for(int jjjj=0;jjjj<5;jjjj++) {
-     cout<<*toialphain; 
+     cout<<*toicoord1in; 
      cout<<*toibolin; 
      sleep(2);