Changeset 1516 in Sophya for trunk/ArchTOIPipe/ProcWSophya/toi2map.cc

Timestamp:

Jun 12, 2001, 6:00:35 PM (24 years ago)

Author:

cmv

Message:

Re-shape de toi2map cmv 12/6/01

File:

• trunk/ArchTOIPipe/ProcWSophya/toi2map.cc (modified) (5 diffs)

trunk/ArchTOIPipe/ProcWSophya/toi2map.cc

@@ -4 +4 @@
 #include "ctimer.h"
 #include "toi2map.h"
-#include "xastropack.h"
 
 ////////////////////////////////////////////////////////////////////////
 TOI2Map::TOI2Map(SphereHEALPix<r_8>* sph,SphereHEALPix<r_8>* wsph)
-  : mSph(sph), mWSph(wsph), mWSphInternal(false), mTypCoor(false), fTypCoor(false), mActualYear(2001.)
+  : mSph(sph), mWSph(wsph), mWSphInternal(false)
 {
+ SetEquinox();
+ SetCoorIn();
+ SetCoorOut();
+
  if(mSph->NbPixels()<1) {
   cout<<"TOI2Map::TOI2Map() Bad number of pixels in sphere mSph "
@@ -44 +47 @@
   declareInput("Coord1In");     // input index 0
   declareInput("Coord2In");     // input index 1
-  declareInput("BoloIn");      // input index 2
+  declareInput("BoloIn");       // input index 2
 }
 
@@ -61 +64 @@
   throw ParmError("TOI2Map::run() Output TOI (Coord1In or Coord2In or BoloIn) not connected!");
 }
+if( !(mTypCoorIn&TypCoordEq || mTypCoorIn&TypCoordGal) ) {
+  cout<<"TOI2Map::run() - Input Coordinates not Eq or Gal! "<<endl;
+  throw ParmError("TOI2Map::run() - Input Coordinates not Eq or Gal!");
+}
+if( !(mTypCoorOut&TypCoordEq || mTypCoorOut&TypCoordGal) ) {
+  cout<<"TOI2Map::run() - Output Coordinates not Eq or Gal! "<<endl;
+  throw ParmError("TOI2Map::run() - Output Coordinates not Eq or Gal!");
+}
 
 //---------------------------------------------------------
@@ -69 +80 @@
 uint_4 mNSnFill=0, mNpixFill=0, NFill[NFILL];
 for(ii=0;ii<NFILL;ii++) NFill[ii]=0;
+double mjd = MJDfrYear(mActualYear);
 
 // Remplissage des spheres
 for(int s=snb;s<=sne;s++) {
-  int_8 fgbolo = 0; 
+  int_8 fgbolo = 0;
   double bolo;
+  //              Equatoriales   /   Galactiques
+  // coord1,2 =   alpha,delta    /   gLon,gLat
+  double coord1 = getData(0,s);
+  double coord2 = getData(1,s);
 
-  double coord1 = getData(0,s); // gLat ou delta entre [-90,90] en degres
-  double coord2 = getData(1,s); // gLon entre [0,360[ en degres ou alpha entre [0,24[ en heures
   getData(2,s,bolo,fgbolo);
+  if(bolo<-32767.) continue; // Bidouille (A CORRIGER QUAND FGBOLO DISPO)
 
-  if(coord2<-90. || coord2>90.) fgbolo=1; 
-  if((coord1<0.) || (!mTypCoor && coord1>=24.) || (mTypCoor && coord1>=360.) ) {fgbolo=1; 
-  cout << "!!!!!!!!" <<coord1 << endl;
+  // sphere phi   entre [0,2*Pi] en radians
+  // sphere theta entre [0,Pi]   en radians
+  double phi,theta;
+  CoordConvertToStd(mTypCoorIn,coord1,coord2);
+  if(mTypCoorIn&TypCoordEq && mTypCoorOut&TypCoordGal) { // Eq -> Gal
+    EqtoGal(mjd,coord1,coord2,&coord1,&coord2);
+    phi   = coord1 * M_PI/180.;
+  } else if(mTypCoorIn&TypCoordGal && mTypCoorOut&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.;
   }
-  
-  
-  
-  if(bolo<-32767.) fgbolo=1;   // Bidouille Archeops
-  
-  if(!fgbolo) {
-    // sphere phi entre [0,2*Pi] en radians
-    // sphere phi entre [0,2*Pi] en radians
-    // sphere theta entre [0,Pi] en radians
-    double phi,theta;
-    if(fTypCoor && !mTypCoor) { //on a alpha,delta et on veut l,b
-      double mjd = MJDfrYear(mActualYear);
-      EqtoGal(mjd,coord1,coord2,&coord1,&coord2);
-    } 
+  theta = (90.-coord2) * M_PI/180.;
+  if(phi<0.   || phi>=2*M_PI) continue;
+  if(theta<0. || theta>=M_PI) continue;
 
-    if(!fTypCoor && mTypCoor) { //on a l,b et on veut alpha,delta
-      double mjd = MJDfrYear(mActualYear);
-      GaltoEq(mjd,coord1,coord2,&coord1,&coord2);
-    }
-    
-    if(fTypCoor) phi = coord1 * M_PI/180.;
-    else  phi = coord1 * M_PI/12.;
-
-    theta = (90.-coord2)*M_PI/180.;
-
-    int_4 ipix = mSph->PixIndexSph(theta,phi);
-    (*mSph)(ipix) += bolo;    
-    ((*mWSph)(ipix))++;
-    mNSnFill++;
-  }
+  int_4 ipix = mSph->PixIndexSph(theta,phi);
+  (*mSph)(ipix) += bolo;    
+  ((*mWSph)(ipix))++;
+  mNSnFill++;
 }
 
@@ -129 +133 @@
      <<"  mNpixFill="<<mNpixFill
      <<"  mNSnFill="<<mNSnFill<<endl
-     <<" --> FracSky="<<mNpixFill*100./(double)mSph->NbPixels()<<"%"<<endl
-     <<"NFill["<<NFILL<<"] = "<<endl;
- for(ii=0;ii<NFILL;ii++) cout<<NFill[ii]<<" ";
+     <<"  --> FracSky="<<mNpixFill*100./(double)mSph->NbPixels()<<"%"
+     <<"  NFill["<<NFILL<<"] ="<<endl;
+ for(ii=0;ii<NFILL;ii++) {cout<<NFill[ii]<<" "; if(ii%10==9) cout<<endl;}
  cout<<endl;