Changeset 1520 in Sophya for trunk/ArchTOIPipe/ProcWSophya/map2toi.cc

Timestamp:

Jun 12, 2001, 7:40:20 PM (24 years ago)

Author:

cmv

Message:

Re-shape de map2toi cmv 12/6/01

File:

• trunk/ArchTOIPipe/ProcWSophya/map2toi.cc (modified) (4 diffs)

trunk/ArchTOIPipe/ProcWSophya/map2toi.cc

@@ -8 +8 @@
 : mSph(sph)
 {
+ SetEquinox();
+ SetCoorIn();
+ SetCoorMap();
 }
 
@@ -17 +20 @@
 void Map2TOI::init() {
   cout << "Map2TOI::init" << endl;
-  declareInput("AlphaIn");     // input  index 0
-  declareInput("DeltaIn");     // input  index 1
-  declareOutput("AlphaOut");   // output index 0
-  declareOutput("DeltaOut");   // output index 1
-  declareOutput("BoloOut");    // output index 2
+  declareInput("Coord1In");     // input  index 0
+  declareInput("Coord2In");     // input  index 1
+  declareOutput("Coord1Out");   // output index 0
+  declareOutput("Coord2Out");   // output index 1
+  declareOutput("BoloOut");     // output index 2
 }
 
@@ -39 +42 @@
 }
 if(!checkInputTOIIndex(0) || !checkInputTOIIndex(1)) {
-  cout<<"Map2TOI::run() - Input TOI (AlphaIn or DeltaIn) not connected! "<<endl;
-  throw ParmError("Map2TOI::run() Output TOI (AlphaIn or DeltaIn) not connected!");
+  cout<<"Map2TOI::run() - Input TOI (Coord1In or Coord2In) not connected! "<<endl;
+  throw ParmError("Map2TOI::run() Output TOI (Coord1In or Coord2In) not connected!");
 }
-bool connected_AlphaOut = checkOutputTOIIndex(0);
-bool connected_DeltaOut = checkOutputTOIIndex(1);
+bool connected_Coord1Out = checkOutputTOIIndex(0);
+bool connected_Coord2Out = checkOutputTOIIndex(1);
 if(!checkOutputTOIIndex(2)) {
   cout<<"Map2TOI::run() - Output TOI (BoloOut) not connected! "<<endl;
   throw ParmError("Map2TOI::run() Output TOI (BoloOut) not connected!");
+}
+if( !(mTypCoorIn&TypCoordEq || mTypCoorIn&TypCoordGal) ) {
+  cout<<"Map2TOI::run() - CoordIn Coordinates not Eq or Gal! "<<endl;
+  throw ParmError("Map2TOI::run() - CoordIn Coordinates not Eq or Gal!");
+}
+if( !(mTypCoorMap&TypCoordEq || mTypCoorMap&TypCoordGal) ) {
+  cout<<"Map2TOI::run() - Sphere Coordinates not Eq or Gal! "<<endl;
+  throw ParmError("Map2TOI::run() - Sphere Coordinates not Eq or Gal!");
 }
 
@@ -52 +63 @@
 try {
 
+uint_4 mSnRead=0, mSnFilled=0;
+double mjd = MJDfrYear(mActualYear);
+
 for(int k=snb;k<=sne;k++) {
 
   int_8 flg = 0;
-  //  double alpha = getData(0,k); // alpha entre [0,24[ en heures
-  //  if(alpha<0. || alpha>=24.) flg=1;
-  // On passe l'angle en longitude en degres - de [0,360[ deg
-  double alpha = getData(0,k); // alpha entre [0,360[ en heures
-  if(alpha<0. || alpha>=360.) flg=1;
-  double delta = getData(1,k); // delta entre [-90,90] en degres
-  if(delta<-90. || delta>90.) flg=1;
-  double bolo = -99999.;
+  double coord1 = getData(0,k);
+  double coord2 = getData(1,k);
+  double c1 = coord1, c2=coord2, bolo = -99999.;
+  mSnRead++;
+
+  // Convert CoordIn to Standard
+  double phi=-1.,theta;
+  CoordConvertToStd(mTypCoorIn,c1,c2);
+  if(mTypCoorIn&TypCoordEq && mTypCoorMap&TypCoordGal) { // CIn=Eq  CMap=Gal
+    EqtoGal(mjd,c1,c2,&c1,&c2);
+    phi   = c1 * M_PI/180.;
+  } else if(mTypCoorIn&TypCoordGal && mTypCoorMap&TypCoordEq) { // CIn=Gal CMap=Eq
+    GaltoEq(mjd,c1,c2,&c1,&c2);
+    phi   = c1 * M_PI/12.;
+  } else if(mTypCoorMap&TypCoordGal) { // CIn=Gal CMap=Gal
+    phi   = c1 * M_PI/180.;
+  } else if(mTypCoorMap&TypCoordEq) { // CIn=Eq CMap=Eq
+    phi   = c1 * M_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;
 
   if(!flg) {
-    //    double phi   = alpha *M_PI/12.; // sphere phi entre [0,2*Pi] en radian
-    double phi   = alpha *M_PI/180.; // sphere phi entre [0,2*Pi] en radian
-    double teta  = (90.-delta) *M_PI/180.; // sphere teta entre [0,Pi] en radian
-    bolo = mSph(teta,phi);
+    bolo = mSph(theta,phi);
+    mSnFilled++;
   }
 
-  if(connected_AlphaOut) putData(0,k,alpha,0);
-  if(connected_DeltaOut) putData(1,k,delta,0);
+  if(connected_Coord1Out) putData(0,k,coord1,flg);
+  if(connected_Coord2Out) putData(1,k,coord2,flg);
   putData(2,k,bolo,flg);
 }
+
+cout<<"TOI2Map::run: Samples Read "<<mSnRead<<" Filled "<<mSnFilled<<endl;
 
 //---------------------------------------------------------