Changeset 2983 in Sophya

Timestamp:

Jun 21, 2006, 5:51:29 PM (19 years ago)

Author:

ansari

Message:

Ajout commande szidx2resol ds le module skymapmdule.cc , Reza 21/6/2006

File:

• trunk/SophyaPI/ProgPI/skymapmodule.cc (modified) (4 diffs)

trunk/SophyaPI/ProgPI/skymapmodule.cc

@@ -47 +47 @@
 
   void Resol2SizeIndex(vector<string>& tokens);
+  void SizeIndex2Resol(vector<string>& tokens);
   void TypeMap(vector<string>& tokens);
   void Map2Double(vector<string>& tokens);
@@ -105 +106 @@
 usage += "\n Usage: resol2szidx resol";
 mpiac->RegisterCommand(kw, usage, this, hgrp);
+kw = "szidx2resol";
+usage = "Compute resolution for a given SizeIndex (=nside for HEALPix)";
+usage += "\n Usage: szidx2resol szidx_m";
+mpiac->RegisterCommand(kw, usage, this, hgrp);
 
 kw = "typemap";
@@ -287 +292 @@
   }
   Resol2SizeIndex(tokens);
+} else if(kw == "szidx2resol") {
+  if(tokens.size()<1) {
+    cout<<"Usage: szidx2resol szidx_m"<<endl;
+    return(0);
+  }
+  SizeIndex2Resol(tokens);
 } else if(kw == "typemap") {
   if(tokens.size()<1) {
@@ -440 +451 @@
       << resrad << " radian  NPix ~= " << 4.*M_PI/(resrad*resrad) 
       << " \n  ===> nside=m_HEALPix= " << mh << "  m_ThetaPhi= " << mt << endl;
+ return;
+}
+
+void skymapmoduleExecutor::SizeIndex2Resol(vector<string>& tokens)
+{
+ int_4 m = atoi(tokens[0].c_str());
+ double resh = SphereHEALPix<r_4>::SizeIndexToResol(m);
+ double rest = SphereThetaPhi<r_4>::SizeIndexToResol(m);
+ cout << "SizeIndex2Resol: SizeIndex (m)= " << m
+      << "\n HEALPix:  " << Angle(resh).ToArcMin() << " arcmin "
+      << " NPix= " << 12*m*m 
+      << "\n ThetaPhi: " << Angle(rest).ToArcMin() << " arcmin "
+      << " NPix ~= " <<  4.*M_PI/(rest*rest) << endl;
  return;
 }