Changeset 1466 in Sophya

Timestamp:

Apr 11, 2001, 2:08:19 PM (24 years ago)

Author:

cmv

Message:

adapt getData cmv 11/4/01

Location:

trunk/ArchTOIPipe/ProcWSophya

Files:

• map2toi.cc (modified) (1 diff)
• toi2map.cc (modified) (3 diffs)

trunk/ArchTOIPipe/ProcWSophya/map2toi.cc

@@ -54 +54 @@
 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.) continue;
+  if(alpha<0. || alpha>=24.) flg=1;
   double delta = getData(1,k); // delta entre [-90,90] en degres
-  if(delta<-90. || delta>90.) continue;
+  if(delta<-90. || delta>90.) flg=1;
+  double bolo = -99999.;
 
-  double phi   = alpha *M_PI/12.; // sphere phi entre [0,2*Pi] en radian
-  double teta  = (90.-delta) *M_PI/180.; // sphere teta entre [0,Pi] en radian
-
-  double bolo  = mSph(teta,phi);
+  if(!flg) {
+    double phi   = alpha *M_PI/12.; // 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);
+  }
 
   if(connected_AlphaOut) putData(0,k,alpha,0);
   if(connected_DeltaOut) putData(1,k,delta,0);
-  putData(2,k,bolo,0);
+  putData(2,k,bolo,flg);
 }
 

trunk/ArchTOIPipe/ProcWSophya/toi2map.cc

@@ -72 +72 @@
 // Remplissage des spheres
 for(int s=snb;s<=sne;s++) {
-  int_4 fgbolo = 0; 
+  int_8 fgbolo = 0; 
+  double bolo;
 
   double alpha = getData(0,s); // alpha entre [0,24[ en heures
   double delta = getData(1,s); // delta entre [-90,90] en degres
-  double bolo  = getData(2,s);
-  //fgbolo = getFlag(2,s);
+  getData(2,s,bolo,fgbolo);
 
   if(delta<-90. || delta>90.)  fgbolo=1;
@@ -91 +91 @@
       EqtoGal(mjd,alpha,delta,&alpha,&delta);
       phi   = alpha * M_PI/180.;
-    } else phi   = alpha * M_PI/12.;
+    } else phi = alpha * M_PI/12.;
     theta = (90.-delta)*M_PI/180.;
     int_4 ipix = mSph->PixIndexSph(theta,phi);
@@ -98 +98 @@
     mNSnFill++;
   }
-
 }