Changeset 1191 in Sophya for trunk/SigPredictor/lobelevshealpix.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/lobelevshealpix.cc (modified) (6 diffs)

trunk/SigPredictor/lobelevshealpix.cc

@@ -1 @@
-#include "fitsspherehealpix.h"
+#include <stdlib.h>
+#include <iostream>
 #include "lobelevshealpix.h"
-
+#include "lightptsrclevsinband.h"
 
 LobeLevSHealpix::LobeLevSHealpix(char FileOflobeMap[],double angMax_degre, double MeanFreq_Hz)
@@ -9 +10 @@
                 cosanglemax=cos(angleMax);
                 
-                sprintf(Name,"Lobe LevelS Healpix");
+                sprintf(Name,"Lobe_LevelS_Healpix");
+                
+        // Version 4
+                // On lit la Sphere Healpix
+                FitsInFile File(FileOflobeMap); 
+                FITS_SphereHEALPix <r_4> FSphere;
+                FSphere.Read(File);
+                cerr<<"Lecture des donnŽes du lobe LevelSHealpix"<<endl;
+                
+                // On stocke les donnees
+                pSphere=new SphereHEALPix<r_4> ((SphereHEALPix<r_4>) FSphere );
 
 /*              // Version 1            
@@ -16 +27 @@
                 FITS_SphereHEALPix<float> FSphere(*pSphere);
                 FSphere.Read(FileOflobeMap);
-*/
-/*              // Version 2
+
+                // Version 2
                 FITS_SphereHEALPix<float> FSphere(FileOflobeMap);       
                 cerr<<"Lecture des donnŽes du lobe LevelSHealpix"<<endl;
                 pSphere=new SphereHEALPix<r_4>;
                 *pSphere= (SphereHEALPix<r_4>) FSphere;
-*/              
-/*              // Version 3
+        
+                // Version 3
                 FITS_SphereHEALPix<float> FSphere(FileOflobeMap);       
                 cerr<<"Lecture des donnŽes du lobe LevelSHealpix"<<endl;
@@ -29 +40 @@
                 pSphere=&Sphere;
 */
+
+                 
+/*         pSphere= new SphereHEALPix<float>;
+           *pSphere= ( SphereHEALPix<float> ) FSphere;
+*/   
+           /*
+                SphereHEALPix <float>* pUNSphere=pFSphere->getObj();
+                long Nlat=pUNSphere->SizeIndex();
+                pSphere=new SphereHEALPix <float>(Nlat);
+                *pSphere=(SphereHEALPix <float>) (*pFSphere);
+                delete pFSphere;
+                */
 }
 
 LobeLevSHealpix::~LobeLevSHealpix()
 {       delete pSphere;
+        
 }
 
@@ -39 +63 @@
 {       double cosinusZ=VP*VInteg;
         if(cosinusZ<cosanglemax) return 0.;
+        if(cosinusZ==1.) return (*pSphere)(pSphere->PixIndexSph(0.,0.));
         else
         {       UnitVector VX=VP^VY;
@@ -47 +72 @@
                 
                 double norm=Vproj.Norm();
-                if(norm==0.)
-                {       phi=0.;}
-                else
-                {       Vproj/=norm;
-                        phi = scangle(VY*Vproj,VX*Vproj);       // Fonction de geometry.h
-                }                                                                                                       // double scangle(double sinus, double cosinus)         
+                Vproj/=norm;
+                phi = scangle(VY*Vproj,VX*Vproj);       // Fonction de geometry.h
+                                                                                                                // double scangle(double sinus, double cosinus)         
           // double PixNumb = pSphere->PixIndexSph(theta,phi);
            double weigthVal= (*pSphere)(pSphere->PixIndexSph(theta,phi));