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

Timestamp:

Oct 24, 2002, 12:30:50 PM (23 years ago)

Author:

aubourg

Message:

constructeur pour pipe xml

File:

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

trunk/ArchTOIPipe/ProcWSophya/toi2map.cc

@@ -3 +3 @@
 //                               Christophe Magneville
 //                               Reza Ansari
-// $Id: toi2map.cc,v 1.30 2002-08-05 16:25:13 cmv Exp $
+// $Id: toi2map.cc,v 1.31 2002-10-24 10:30:50 aubourg Exp $
 
 #include "machdefs.h"
@@ -13 +13 @@
 #include "smathconst.h"
 
+#include "fitsspherehealpix.h"
+
+
+////////////////////////////////////////////////////////////////////////
+TOI2Map::TOI2Map(string mapFName, string wmapFName, int nside) 
+  : mMap(new SphereHEALPix<r_8>(nside)), mWMap(new SphereHEALPix<r_8>(nside)),
+    mWMapInternal(false), totnscount(0), usednscount(0),
+    mMapFName(mapFName), mWMapFName(wmapFName)
+{
+  commonConst();
+}
+
 ////////////////////////////////////////////////////////////////////////
 TOI2Map::TOI2Map(PixelMap<r_8>* map,PixelMap<r_8>* wmap)
-  : mMap(map), mWMap(wmap), mWMapInternal(false), totnscount(0), usednscount(0)
-{
+  : mMap(map), mWMap(wmap), mWMapInternal(false), totnscount(0), usednscount(0),
+    mMapFName(""), mWMapFName("")
+{
+  commonConst();
+}
+
+void TOI2Map::commonConst() {
  SetEquinox();
  SetCoorIn();
@@ -230 +247 @@
 #endif
 
+ if (mMapFName != "") {
+   cout << "saving maps" << endl;
+   {
+     FitsOutFile sfits(mMapFName,FitsFile::clear);
+     sfits << * (SphereHEALPix<r_8>*)mMap;
+   }
+   {
+     FitsOutFile sfits(mWMapFName,FitsFile::clear);
+     sfits << * (SphereHEALPix<r_8>*)mWMap;
+   }
+ } 
+
 //---------------------------------------------------------
 } catch (PException & exc) {