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

Timestamp:

Jun 14, 2001, 7:53:29 PM (24 years ago)

Author:

cmv

Message:

ameliorations et + de possibilites cmv 14/6/01

File:

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

trunk/ArchTOIPipe/TestPipes/tsttoi2map.cc

@@ -20 +20 @@
      <<"           [-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
      <<"           fitsin_point fitsin_bolo fitsphout [fitsphwout]"<<endl
      <<" -p lp : print level (def=0)"<<endl
@@ -29 +30 @@
      <<" -b label_bolomuv : label fits for bolo value (def=boloMuV)"<<endl
      <<" -n nlat : nlat for Healpix sphere (def=128)"<<endl
-     <<" -i c,h : coordIn caracteristics (c=G/E h=H/D) (def=G,D)"<<endl
+     <<" -i c,h : coordIn caracteristics (c=G/E h=H/D/R) (def=G,D)"<<endl
      <<" -o c,h : idem -i for coordOut"<<endl
+     <<" -m vmin : samples are good if sample value >= vmin"<<endl
+     <<" -M vmax : samples are good if sample value <= vmax"<<endl
+     <<" -f flag : samples are bad if match flag"<<endl
      <<" fitsin_point : fits file for pointing"<<endl
      <<" fitsin_bolo : fits file for bolo values"<<endl
@@ -40 +44 @@
 unsigned long typecoord(char typc,char hd)
 // typc : G=galactiques, E=equatoriales, autres=galactiques
-// hd : H=heure, D=degre, autres=(heure si typc==E, degre si typc==G)
+// hd : H=heure, D=degre, R=radian, autres=(heure si typc==E, degre si typc==G)
 {
   if(typc!='G' && typc!='E') typc='G';
-  if(hd!='H' && hd!='D') {if(typc=='E') hd='H'; else hd='D';}
+  if(hd!='H' && hd!='D' && hd!='R') {if(typc=='E') hd='H'; else hd='D';}
  unsigned long rc=TypCoordUndef;
   if(typc=='G') rc |= TypCoordGal;
     else        rc |= TypCoordEq;
-  if(hd=='D')   rc |= TypCoordDD;
-    else        rc |= TypCoordHD;
+  if(hd=='D')        rc |= TypCoordDD;
+    else if(hd=='R') rc |= TypCoordRR;
+      else           rc |= TypCoordHD;
   return rc;
 }
@@ -60 +65 @@
 int lp = 0, width = 8192;
 int nlat = 128;   // npixel = 12 * nlat^2
-char *label_coord1 = "coord1", *label_coord2 = "coord2", *label_bolomuv = "boloMuV";
+bool tflg=false, tmin=false, tmax=false;
+r_8 vmin=-1.e30, vmax=1.e30; int_8 badflg=0;
+char *label_coord1 = "coord1", *label_coord2 = "coord2"
+   , *label_bolomuv = "boloMuV";
 long sdeb,sfin;
 string fitsphwout = "";
 unsigned long tcoorin=typecoord(), tcoorout=typecoord();
 int c; char t=' ',h=' ';
-while((c = getopt(narg,arg,"hp:s:w:a:d:b:n:i:o:")) != -1) {
+while((c = getopt(narg,arg,"hp:s:w:a:d:b:n:i:o:m:M:f:")) != -1) {
   switch (c) {
   case 's' :
@@ -103 +111 @@
     tcoorout=typecoord(t,h);
     break;
+  case 'm' :
+    sscanf(optarg,"%lf",&vmin);
+    tmin = true;
+    break;
+  case 'M' :
+    sscanf(optarg,"%lf",&vmax);
+    tmax = true;
+    break;
+  case 'f' :
+    sscanf(optarg,"%ld",&badflg);
+    tflg = true;
+    break;
   case 'h' :
   default:
@@ -116 +136 @@
 if(optind+3<narg) fitsphwout = arg[optind+3];
 
-{
-unsigned long tg,te,hd,dd;
 cout<<">>>> tsttoi2map:"<<endl
     <<"Pipe Window Size "<<width<<endl
     <<"Fits Infile Bolo "<<fitsin_bolo<<endl
     <<"  ...label_bolomuv "<<label_bolomuv<<endl;
-tg = tcoorin&TypCoordGal; te = tcoorin&TypCoordEq;
-hd = tcoorin&TypCoordHD;  dd = tcoorin&TypCoordDD;
 cout<<"Fits Infile Pointing "<<fitsin_point<<endl
     <<"  ...label_coord1 "<<label_coord1<<endl
-    <<"  ...label_coord2 "<<label_coord2<<endl
-    <<"  ...... Gal="<<tg<<" Eq="<<te<<"   hour="<<hd<<" deg="<<dd<<endl;
-tg = tcoorout&TypCoordGal; te = tcoorout&TypCoordEq;
-hd = tcoorout&TypCoordHD;  dd = tcoorout&TypCoordDD; 
+    <<"  ...label_coord2 "<<label_coord2<<endl;
 cout<<"Fits Healpix Sphere "<<fitsphout<<endl
-    <<"  ...nlat "<<nlat<<endl
-    <<"  ...... Gal="<<tg<<" Eq="<<te<<"   hour="<<hd<<" deg="<<dd<<endl;
+    <<"  ...nlat "<<nlat<<endl;
 cout<<"Fits Healpix Weight Sphere "<<fitsphwout<<endl;
-}
 
 SophyaInit();
@@ -157 +168 @@
  SphereHEALPix<r_8>* sph = new SphereHEALPix<r_8>(nlat);
  cout<<"SphereHEALPix: Type de map : "<<sph->TypeOfMap()<<endl
-     <<"               Nombre de pixels : "<<sph->NbPixels()<<endl;
+     <<"               Nombre de pixels : "<<sph->NbPixels()<<endl
+     <<"               Nlat : "<<sph->SizeIndex()<<endl;
 
  // Creation de la sphere de poids Healpix
@@ -163 +175 @@
  if(fitsphwout.size()>0) {
    wsph = new SphereHEALPix<r_8>;
-   cout<<"SphereHEALPix Weight: Type de map : "<<wsph->TypeOfMap()<<endl
-       <<"              Nombre de pixels : "<<wsph->NbPixels()<<endl;
+   cout<<"SphereHEALPix Weight Created"<<endl;
  }
 
@@ -173 +184 @@
  toi2m.SetCoorIn((TypAstroCoord) tcoorin);
  toi2m.SetCoorOut((TypAstroCoord) tcoorout);
+ toi2m.SetTestFlag(tflg,badflg);
+ toi2m.SetTestMin(tmin,vmin);
+ toi2m.SetTestMax(tmax,vmax);
+ toi2m.Print(cout);
 
  // Definition des tuyaux