Changeset 1191 in Sophya for trunk/SigPredictor/lightptsrclevsinband.cc

Timestamp:

Sep 19, 2000, 2:28:59 PM (25 years ago)

Author:

ansari

Message:

cleaned up for maching LevelS upgrade ongoing

D.Y.

File:

• trunk/SigPredictor/lightptsrclevsinband.cc (modified) (4 diffs)

trunk/SigPredictor/lightptsrclevsinband.cc

@@ -28 +28 @@
         while (HFI)
         {       // Data in file in radian and Jansky
+        //      cout<<compteur<<'\t'<<HFI.fail()<<'\n'<<endl;;
                 HFI>>theta;
                 HFI>>phi;
@@ -44 +45 @@
                 {       HFI.clear();
                         HFI.getline(Line, 199); 
-                        if(HFI.eof())
+                        if(HFI.eof())           // Fin de lecture du fichier
                         {       cerr<< "Fin de lecture du fichier: "<<PtSrcfile<<endl;
                                 cerr<< "Point source number : "<<MapOfPtSrc.size()<<endl;
+                                HFI.close();
                                 break;
                         }
@@ -67 +69 @@
                         UneSource.flux857G=f857*1.e-26;
                         
-                        const pair<const UnitVector, LevSPtSrcData> UneValeur=make_pair(UnCoord, UneSource);
+                        const pair<const UnitVector,LevSPtSrcData>UneValeur=make_pair(UnCoord, UneSource);
                         
                         MapOfPtSrc.insert(MapOfPtSrc.end(),UneValeur);
@@ -76 +78 @@
         return;
 }
+
+#include "fitsspherehealpix.h"
+
+void LightPtSrcLevSInBand::DumpMap(LevSPanckBand band)
+{       // Write a Mollweide projection of the Point Sources Sky
+        // power in the LevSPanckBand band
+        SphereHEALPix <r_4> Sphere(128);
+
+        multimap <UnitVector,LevSPtSrcData,LevSPtSrcCmpUnitVec>::iterator iter;
+        multimap <UnitVector,LevSPtSrcData,LevSPtSrcCmpUnitVec>::iterator iterBegin=MapOfPtSrc.begin();
+        multimap <UnitVector,LevSPtSrcData,LevSPtSrcCmpUnitVec>::iterator iterEnd=MapOfPtSrc.end();
+        
+        long indexSphere;
+        double theta;
+        double phi;
+        double power;
+        
+        for(iter=iterBegin; !(iter==iterEnd); iter++)
+        {       theta=(*iter).first.Theta();
+                phi=(*iter).first.Phi();
+                power=(*iter).second.getPower(band);
+                
+                indexSphere=Sphere.PixIndexSph(theta,phi);
+                Sphere(indexSphere) += power;
+        }
+        
+        FITS_SphereHEALPix <r_4> FSphere(Sphere);
+        FSphere.Mollweide_picture_projection("SkyOfPtSrcInBandMollweide.fits");
+
+        return;
+}