Changeset 1632 in Sophya for trunk

Timestamp:

Aug 10, 2001, 4:35:27 PM (24 years ago)

Author:

cmv

Message:

add cut lat_gal cmv+rz 10/8/01

File:

• trunk/SophyaProg/PrgMap/cremskfrsph.cc (modified) (8 diffs)

trunk/SophyaProg/PrgMap/cremskfrsph.cc

@@ -26 +26 @@
  -n : negate the condition
       (switch to: sphere_value<min || max<sphere_value)
+ -b bmin,bmax : on coupe en latitude entre [bmin,bmax] degres (def=no+cut)
  -v valmsk,inimsk:
     valmsk: value to be put into mask if ok (def=1.)
@@ -42 +43 @@
     <<" -n : negate the condition"<<endl
     <<"      (switch to: sphere_value<min || max<sphere_value)"<<endl
+    <<" -b bmin,bmax : on coupe en latitude entre [bmin,bmax] degres (def=no+cut)"<<endl
     <<" -v valmsk,inimsk:"<<endl
     <<"    valmsk: value to be put into mask if ok (def=1.)"<<endl
@@ -50 +52 @@
 int main(int narg, char* arg[])
 {
-double vmin=-1.e30,vmax=1.e30,vmask=1.,vmaskini=0.;
+double vmin=-1.e30,vmax=1.e30,vmask=1.,vmaskini=0., bmin=1.,bmax=-1.;
 bool tstmin=false,tstmax=false,negate=false;
 string dum;
 int c;
-while((c = getopt(narg,arg,"hnm:M:v:")) != -1) {
+while((c = getopt(narg,arg,"hnm:M:v:b:")) != -1) {
   switch (c) {
   case 'n' :
@@ -69 +71 @@
   case 'v' :
     sscanf(optarg,"%lf,%lf",&vmask,&vmaskini);
+    break;
+  case 'b' :
+    sscanf(optarg,"%lf,%lf",&bmin,&bmax);
     break;
   case 'h' :
@@ -88 +93 @@
     <<"  ...mask set value "<<vmask<<endl
     <<"     mask init value "<<vmaskini<<endl
+    <<"Cut latitude between "<<bmin<<" and "<<bmax<<" deg"<<endl
     <<"Condition on data to set the mask: "
     <<dum<<"( "<<vmin<<" <= V <= "<<vmax<<" )"<<endl;
@@ -104 +110 @@
 SphereHEALPix<r_8> sphm(sph,false);
 
-// Filling Mask Sphere
+// Filling Mask Sphere with cut on values
 cout<<"Filling Mask"<<endl;
+double thetamin = (90.-bmax) * M_PI/180.;
+double thetamax = (90.-bmin) * M_PI/180.;
 uint_4 nmask = 0;
 for(int_4 i=0;i<sph.NbPixels();i++) {
@@ -113 +121 @@
   if((tstmin && sph(i)<vmin) || (tstmax && sph(i)>vmax)) skp = !negate;
   if(skp) continue;   // Do nothing
+
+  if(bmin<bmax) {
+    double theta,phi;
+    sphm.PixThetaPhi(i,theta,phi);
+    if(theta>=thetamin && theta<=thetamax) continue;
+  }
 
   sphm(i) = vmask;
@@ -124 +138 @@
 {
 dum = "rm -f "; dum += sphmsk; system(dum.c_str());
-FitsOutFile swfits(sphmsk);
+FitsOutFile swfits(sphmsk,FitsFile::clear);
 cout<<"Writing Mask Sphere Fits file"<<endl;
 swfits<<sphm;
 }
 
-exit(0);
+return 0;
 }