Changeset 2228 in Sophya for trunk

Timestamp:

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

Author:

aubourg

Message:

constructeur pour pipe xml

Location:

trunk/ArchTOIPipe/ProcWSophya

Files:

• toi2map.cc (modified) (3 diffs)
• toi2map.h (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) {

trunk/ArchTOIPipe/ProcWSophya/toi2map.h

@@ -5 +5 @@
 //                               Christophe Magneville
 //                               Reza Ansari
-// $Id: toi2map.h,v 1.17 2002-06-18 20:55:21 ansari Exp $
+// $Id: toi2map.h,v 1.18 2002-10-24 10:30:50 aubourg Exp $
 
 #ifndef TOI2MAP_H
@@ -38 +38 @@
 public:
                  TOI2Map(PixelMap<r_8>* map,PixelMap<r_8>* wmap=NULL);
+                 TOI2Map(string mapFName, string wmapFName, int nside);
                  //TOI2Map(SphereHEALPix<r_8>* map,SphereHEALPix<r_8>* wmap=NULL);
   virtual       ~TOI2Map();
@@ -81 +82 @@
 
 protected:
+  void commonConst();
+
   PixelMap<r_8>* mMap;
   PixelMap<r_8>* mWMap;
   bool mWMapInternal;
+
+  string mMapFName, mWMapFName;
 
   unsigned long mTypCoorIn, mTypCoorMap;