Changeset 1809 in Sophya for trunk/ArchTOIPipe/ProcWSophya

Timestamp:

Dec 3, 2001, 9:38:03 PM (24 years ago)

Author:

cmv

Message:

• map2toi et toi2map avec nouvelle interface sys coor et unites.
• les modifs map2toi (LocalMap) de EA sont introduites + certains bugs corriges.

cmv 3/12/01

Location:

trunk/ArchTOIPipe/ProcWSophya

Files:

• map2toi.cc (modified) (7 diffs)
• map2toi.h (modified) (3 diffs)
• toi2map.cc (modified) (10 diffs)
• toi2map.h (modified) (5 diffs)

trunk/ArchTOIPipe/ProcWSophya/map2toi.cc

@@ -3 +3 @@
 //                               Christophe Magneville
 //                               Reza Ansari
-// $Id: map2toi.cc,v 1.9 2001-11-13 16:22:47 aubourg Exp $
+// $Id: map2toi.cc,v 1.10 2001-12-03 20:38:02 cmv Exp $
 
 #include "toimanager.h"
@@ -9 +9 @@
 #include "ctimer.h"
 #include "map2toi.h"
+// La valeur "Pi" doit etre celle de smathconst.h a cause du test sur theta
+#include "smathconst.h"
 
 ////////////////////////////////////////////////////////////////////////
@@ -14 +16 @@
 : mSph(sph)
 {
+ SetBad();
  SetEquinox();
  SetCoorIn();
@@ -26 +29 @@
 void Map2TOI::Print(::ostream & os)
 {
-  os<<"Map2TOI::Print -- Sphere NLat = "<<mSph.SizeIndex()<<endl
-
-    <<"   - Equinoxe="<<mActualYear<<endl
-
-    <<"   - TypCoorIn: Gal("<<((mTypCoorIn&TypCoordGal)?1:0)
-    <<") Eq("<<((mTypCoorIn&TypCoordEq)?1:0)
-    <<") Deg("<<((mTypCoorIn&TypCoordDD)?1:0)
-    <<") Hour("<<((mTypCoorIn&TypCoordHD)?1:0)
-    <<") Rad("<<((mTypCoorIn&TypCoordRR)?1:0)<<")"<<endl
-
-    <<"   - TypCoorMap: Gal("<<((mTypCoorMap&TypCoordGal)?1:0)
-    <<") Eq("<<((mTypCoorMap&TypCoordEq)?1:0)
-    <<") Deg("<<((mTypCoorMap&TypCoordDD)?1:0)
-    <<") Hour("<<((mTypCoorMap&TypCoordHD)?1:0)
-    <<") Rad("<<((mTypCoorMap&TypCoordRR)?1:0)<<")"<<endl;
+ os<<"Map2TOI::Print -- Sphere NLat = "<<mSph.SizeIndex()<<endl
+   <<"   - Equinoxe:   "<<mActualYear<<" y"<<endl
+   <<"   - TypCoorIn:  "<<mTypCoorIn<<" = "<<DecodeTypAstro(mTypCoorIn)<<endl
+   <<"   - TypCoorMap: "<<mTypCoorMap<<" = "<<DecodeTypAstro(mTypCoorMap)<<endl;
 }
 
@@ -89 +81 @@
 try {
 
-uint_4 mSnRead=0, mSnFilled=0;
 double mjd = MJDfrYear(mActualYear);
+cout<<"Map2TOI::run() - modified Julian day "<<mjd<<endl;
+
+uint_4 mSnRead=0, mSnFilled=0, BadCoorRange=0;
 
 for(int k=snb;k<=sne;k++) {
@@ -100 +94 @@
   mSnRead++;
 
-  // Convert CoordIn to Standard
-  double phi=-1.,theta;
-  CoordConvertToStd(mTypCoorIn,c1,c2);
+  // Conversion de CoordIn to Standard
+  double phi=-1.;
+  CoordConvertToStd(mTypCoorIn,&c1,&c2);
+
+  // Conversion dans le systeme astronomique approprie
   if(mTypCoorIn&TypCoordEq && mTypCoorMap&TypCoordGal) { // CIn=Eq  CMap=Gal
     EqtoGal(mjd,c1,c2,&c1,&c2);
-    phi   = c1 * M_PI/180.;
+    phi   = c1 * Pi/180.;
   } else if(mTypCoorIn&TypCoordGal && mTypCoorMap&TypCoordEq) { // CIn=Gal CMap=Eq
     GaltoEq(mjd,c1,c2,&c1,&c2);
-    phi   = c1 * M_PI/12.;
+    phi   = c1 * Pi/12.;
   } else if(mTypCoorMap&TypCoordGal) { // CIn=Gal CMap=Gal
-    phi   = c1 * M_PI/180.;
+    phi   = c1 * Pi/180.;
   } else if(mTypCoorMap&TypCoordEq) { // CIn=Eq CMap=Eq
-    phi   = c1 * M_PI/12.;
+    phi   = c1 * Pi/12.;
   }
-  theta = (90.-c2) * M_PI/180.;
-  if(phi<0.   || phi>=2*M_PI) flg=FlgToiOut;
-  if(theta<0. || theta>=M_PI) flg=FlgToiOut;
+  ToCoLat(&c2,TypUniteD);
+  double theta = c2 * Pi/180.;
+  if(phi<0. || phi>=2*Pi || theta<0. || theta>Pi)
+               {BadCoorRange++; flg = mBadFlag; }
 
   if(!flg) {
@@ -128 +125 @@
 }
 
-cout<<"TOI2Map::run: Samples Read "<<mSnRead<<" Filled "<<mSnFilled<<endl;
+cout<<"TOI2Map::run: Samples Read "<<mSnRead
+    <<" Filled "<<mSnFilled
+    <<" BadCoorRange="<<BadCoorRange<<endl;
 
 //---------------------------------------------------------

trunk/ArchTOIPipe/ProcWSophya/map2toi.h

@@ -5 +5 @@
 //                               Christophe Magneville
 //                               Reza Ansari
-// $Id: map2toi.h,v 1.6 2001-11-13 16:22:47 aubourg Exp $
+// $Id: map2toi.h,v 1.7 2001-12-03 20:38:02 cmv Exp $
 
 #ifndef MAP2TOI_H
@@ -42 +42 @@
   virtual void  run(void);
 
-  // Pour savoir si on a de HourDeg,DegDeg,RadRad
+  // Flag pour les mauvais samples
+  inline void SetBad(unsigned long long flg=FlgToiOut)
+                    {mBadFlag = flg;}
+
+  // Pour definir la partie astrometrie
   inline void SetEquinox(double actualyear=2000.)
               {mActualYear = actualyear;}
-  inline void SetCoorIn(TypAstroCoord mfg=TypCoordGalStd)
+
+  inline void SetCoorIn(unsigned long mfg=TypCoordGalStd)
               {mTypCoorIn = mfg;}
-  inline void SetCoorMap(TypAstroCoord mfg=TypCoordGalStd)
+  inline void SetCoorIn(char const *ctype)
+              {mTypCoorIn = DecodeTypAstro(ctype);}
+
+  inline void SetCoorMap(unsigned long mfg=TypCoordGal)
               {mTypCoorMap = mfg;}
+  inline void SetCoorMap(char const *ctype)
+              {mTypCoorMap = DecodeTypAstro(ctype);}
+
   // Print
   void Print(::ostream & os);
@@ -54 +65 @@
 protected:
   SphereHEALPix<r_8>& mSph;
-  TypAstroCoord mTypCoorIn,mTypCoorMap;
+  unsigned long mTypCoorIn,mTypCoorMap;
+  unsigned long long mBadFlag;
   double mActualYear;
 };

trunk/ArchTOIPipe/ProcWSophya/toi2map.cc

@@ -3 +3 @@
 //                               Christophe Magneville
 //                               Reza Ansari
-// $Id: toi2map.cc,v 1.23 2001-12-01 13:35:15 aubourg Exp $
+// $Id: toi2map.cc,v 1.24 2001-12-03 20:38:02 cmv Exp $
 
 #include "machdefs.h"
@@ -10 +10 @@
 #include "ctimer.h"
 #include "toi2map.h"
+// La valeur "Pi" doit etre celle de smathconst.h a cause du test sur theta
+#include "smathconst.h"
 
 ////////////////////////////////////////////////////////////////////////
 TOI2Map::TOI2Map(PixelMap<r_8>* map,PixelMap<r_8>* wmap)
-  : mMap(map), mWMap(map), mWMapInternal(false)
+  : mMap(map), mWMap(wmap), mWMapInternal(false)
 {
  SetEquinox();
  SetCoorIn();
- SetCoorOut();
+ SetCoorMap();
  SetTestFlag();
  SetTestMin();
  SetTestMax();
-
- mIsColat = false;
 
  if(mMap->NbPixels()<1) {
@@ -31 +31 @@
  mMap->SetPixels(0.);
  int nlat = mMap->SizeIndex();
+ string typmap = mMap->TypeOfMap();
 
  if(mWMap==NULL) {
    // We would need a cloning function in maps. $$TBD$$ $CHECK$
-   string typmap = mMap->TypeOfMap();
    if (typmap == "LOCAL") {
      mWMap = new LocalMap<r_8>(*(LocalMap<r_8>*)mMap);
@@ -47 +47 @@
    mWMapInternal = false;
    if(nlat != mWMap->SizeIndex()) {
-     cout << "TOI2Map::TOI2Map() Bad size for sphere mWMap, does not "
-          << "correspond to mMap : " << mMap->SizeIndex() << ", " << mWMap->SizeIndex() << endl;
-     throw ParmError("TOI2Map::TOI2Map() - Different sizes for map and wmap");
+     cout<<"TOI2Map::TOI2Map() Bad size for sphere mWMap, does not "
+         <<"correspond to mMap : "<<mMap->SizeIndex()<<", "<<mWMap->SizeIndex()<<endl;
+     if (typmap == "LOCAL") {
+       throw ParmError("TOI2Map::TOI2Map() - Different sizes for map and wmap");
+     } else if(typmap == "RING") {
+       ((SphereHEALPix<r_8> *)mWMap)->Resize(nlat);
+       cout<<"Resize have been done..."<<endl;
+     } else {
+       cout << "TOI2Map::TOI2Map() cannot handle map of type " << typmap << endl;
+       throw ParmError("TOI2Map::TOI2Map() - bad type of map");
+     }
    }
  }
@@ -68 +76 @@
 void TOI2Map::Print(::ostream & os)
 {
-  os<<"TOI2Map::Print -- Map type " << mMap->TypeOfMap() << " SizeIndex = "<<mMap->SizeIndex()<<endl
-
-    <<"   - Equinoxe="<<mActualYear<<endl
-
-    <<"   - TypCoorIn: Gal("<<((mTypCoorIn&TypCoordGal)?1:0)
-    <<") Eq("<<((mTypCoorIn&TypCoordEq)?1:0)
-    <<") Deg("<<((mTypCoorIn&TypCoordDD)?1:0)
-    <<") Hour("<<((mTypCoorIn&TypCoordHD)?1:0)
-    <<") Rad("<<((mTypCoorIn&TypCoordRR)?1:0)<<")"<<endl
-
-    <<"   - TypCoorOut: Gal("<<((mTypCoorOut&TypCoordGal)?1:0)
-    <<") Eq("<<((mTypCoorOut&TypCoordEq)?1:0)
-    <<") Deg("<<((mTypCoorOut&TypCoordDD)?1:0)
-    <<") Hour("<<((mTypCoorOut&TypCoordHD)?1:0)
-    <<") Rad("<<((mTypCoorOut&TypCoordRR)?1:0)<<")"<<endl
-
-    <<"   - Tests: Flag("<<mTFlag<<") bad="<<mBadFlag
-    <<" / Value Min("<<mTMin<<")="<<mValMin
-    <<" , Max("<<mTMax<<")="<<mValMax<<endl;
+ os<<"TOI2Map::Print -- Map type " << mMap->TypeOfMap() << " SizeIndex = "<<mMap->SizeIndex()<<endl
+   <<"   - Equinoxe="<<mActualYear<<endl
+   <<"   - TypCoorIn:  "<<mTypCoorIn<<" = "<<DecodeTypAstro(mTypCoorIn)<<endl
+   <<"   - TypCoorMap: "<<mTypCoorMap<<" = "<<DecodeTypAstro(mTypCoorMap)<<endl
+   <<"  - Tests: Flag("<<mTFlag<<") bad="<<mBadFlag
+   <<" / Value Min("<<mTMin<<")="<<mValMin
+   <<" , Max("<<mTMax<<")="<<mValMax<<endl;
 }
 
@@ -116 +112 @@
   throw ParmError("TOI2Map::run() - Input Coordinates not Eq or Gal!");
 }
-if( !(mTypCoorOut&TypCoordEq || mTypCoorOut&TypCoordGal) ) {
+if( !(mTypCoorMap&TypCoordEq || mTypCoorMap&TypCoordGal) ) {
   cout<<"TOI2Map::run() - Output Coordinates not Eq or Gal! "<<endl;
   throw ParmError("TOI2Map::run() - Output Coordinates not Eq or Gal!");
@@ -126 +122 @@
 
 int ii;
-uint_4 mNSnFill=0, mNpixFill=0, NFill[NFILL];
+uint_4 mNSnFill=0, mNpixFill=0, NFill[NFILL], BadCoorRange=0;
 for(ii=0;ii<NFILL;ii++) NFill[ii]=0;
 double mjd = MJDfrYear(mActualYear);
@@ -175 +171 @@
   // sphere phi   entre [0,2*Pi] en radians
   // sphere theta entre [0,Pi]   en radians
-  double phi=-1.,theta;
-  CoordConvertToStd(mTypCoorIn,coord1,coord2);
-
-  if(mTypCoorIn&TypCoordEq && mTypCoorOut&TypCoordGal) { // Eq -> Gal
+  double phi=-1.;
+  CoordConvertToStd(mTypCoorIn,&coord1,&coord2);
+
+  if(mTypCoorIn&TypCoordEq && mTypCoorMap&TypCoordGal) { // Eq -> Gal
     EqtoGal(mjd,coord1,coord2,&coord1,&coord2);
-    phi   = coord1 * M_PI/180.;
-  } else if(mTypCoorIn&TypCoordGal && mTypCoorOut&TypCoordEq) { // Gal -> Eq
+    phi   = coord1 * Pi/180.;
+  } else if(mTypCoorIn&TypCoordGal && mTypCoorMap&TypCoordEq) { // Gal -> Eq
     GaltoEq(mjd,coord1,coord2,&coord1,&coord2);
-    phi   = coord1 * M_PI/12.;
-  } else if(mTypCoorOut&TypCoordGal) { // Gal -> Gal
-    phi   = coord1 * M_PI/180.;
-  } else if(mTypCoorOut&TypCoordEq) { // Eq -> Eq
-    phi   = coord1 * M_PI/12.;
+    phi   = coord1 * Pi/12.;
+  } else if(mTypCoorMap&TypCoordGal) { // Gal -> Gal
+    phi   = coord1 * Pi/180.;
+  } else if(mTypCoorMap&TypCoordEq) { // Eq -> Eq
+    phi   = coord1 * Pi/12.;
   }
-
-  if (phi<0) phi += 2*M_PI;
-  if (phi>=2*M_PI) phi -= 2*M_PI;
-  if (mIsColat) {
-    theta = coord2 * M_PI/180;
-  } else {
-    theta = (90.-coord2) * M_PI/180.;     
-  }
-
-  if(phi<0.   || phi>=2*M_PI) continue;
-  if(theta<0. || theta>=M_PI) continue;
+  ToCoLat(&coord2,TypUniteD);
+  double theta = coord2 * Pi/180.;
+
+  if(phi<0. | phi>=2*Pi || theta<0. || theta>Pi)
+                     {BadCoorRange++; continue;}
 
   int_4 ipix = mMap->PixIndexSph(theta,phi);
@@ -212 +202 @@
  for(int_4 i=0;i<mMap->NbPixels();i++) {
    r_8 wf = (*mWMap)(i);
-   if( wf > 0. ) {
-     mNpixFill++;
-     (*mMap)(i) /= wf;
-   }
+   if(wf>0.) {mNpixFill++; (*mMap)(i) /= wf;}
    int_4 nf = int_4(wf);
    if(nf>=NFILL) nf=NFILL-1; NFill[nf]++;
@@ -227 +214 @@
  for(ii=0;ii<NFILL;ii++) {cout<<NFill[ii]<<" "; if(ii%10==9) cout<<endl;}
  cout<<endl;
+ cout<<"  BadCoorRange="<<BadCoorRange<<endl;
 
 #ifndef SANS_BUFFER

trunk/ArchTOIPipe/ProcWSophya/toi2map.h

@@ -5 +5 @@
 //                               Christophe Magneville
 //                               Reza Ansari
-// $Id: toi2map.h,v 1.13 2001-12-01 13:34:09 aubourg Exp $
+// $Id: toi2map.h,v 1.14 2001-12-03 20:38:03 cmv Exp $
 
-#ifndef TOI2GMAP_H
-#define TOI2GMAP_H
+#ifndef TOI2MAP_H
+#define TOI2MAP_H
 
 #include "toiprocessor.h"
@@ -33 +33 @@
 //                     soit Galactiques  (GLong,GLat)
 // Sortie sur une sphere en coordonnees Equatoriales ou Galactiques
-// Si les CoordIn et les CoordOut sont equatoriales -> meme equinoxe!
+// Si les CoordIn et les CoordMap sont equatoriales -> meme equinoxe!
 
 class TOI2Map : public TOIProcessor {
 public:
                  TOI2Map(PixelMap<r_8>* map,PixelMap<r_8>* wmap=NULL);
+                 //TOI2Map(SphereHEALPix<r_8>* map,SphereHEALPix<r_8>* wmap=NULL);
   virtual       ~TOI2Map();
 
@@ -46 +47 @@
   inline void SetEquinox(double actualyear=2000.)
               {mActualYear = actualyear;}
-  inline void SetCoorIn(TypAstroCoord mfg=TypCoordGalStd)
+
+  inline void SetCoorIn(unsigned long mfg=TypCoordGalStd)
               {mTypCoorIn = mfg;}
-  inline void SetCoorOut(TypAstroCoord mfg=TypCoordGalStd)
-              {mTypCoorOut = mfg;}
+  inline void SetCoorIn(const char *ctype)
+              {mTypCoorIn = DecodeTypAstro(ctype);}
+
+  inline void SetCoorMap(unsigned long mfg=TypCoordGal)
+              {mTypCoorMap = mfg;}
+  inline void SetCoorMap(const char *ctype)
+              {mTypCoorMap = DecodeTypAstro(ctype);}
+
   // Test on flag value ? if yes, BAD sample have flag matching mBadFlag
   inline void SetTestFlag(bool tflg=false, uint_8 badflg=FlgToiAll)
@@ -64 +72 @@
   void Print(::ostream & os);
 
-  void SetIsColat(bool b=true) {mIsColat=b;} // latitude is a colatitude, 0-Pi
 protected:
   PixelMap<r_8>* mMap;
@@ -70 +77 @@
   bool mWMapInternal;
 
-  TypAstroCoord mTypCoorIn, mTypCoorOut;
+  unsigned long mTypCoorIn, mTypCoorMap;
   double mActualYear;
 
   bool mTFlag,mTMin,mTMax;
-  bool mIsColat; // temporary fix
   uint_8 mBadFlag;
   r_8 mValMin,mValMax;