source: Sophya/trunk/SigPredictor/lightsrcmappowerinband.cc@ 1670

                        // Dominique YVON, CEA/DAPNIA/SPP 02/2000

#include "lightsrcmappowerinband.h"
#include "fitsspherehealpix.h"

double LightSrcMapPowerInband::powerDensAmpli(double theta,double phi)
{       return (*pMap)(theta,phi);
                // Return power density Amplidude at coordinates
        // (W/m2/Hz/st)
}

double LightSrcMapPowerInband::spectre(double freq)
{       if((freq>freqMax)||(freq<freqMin)) return 0.;
        else return SpectreAmpli;
}

LightSrcMapPowerInband::LightSrcMapPowerInband
        (char filename[],long nside, double MinFreq, double MaxFreq)
  :freqMax(MaxFreq), freqMin(MinFreq), nlat(nside)
{  
   resolution=2*3.14159/4./nlat; // Resolution angulaire de la carte lue pour ce fond.
                                         // Radians
   SpectreAmpli=1./(freqMax-freqMin);
   if(MinFreq>MaxFreq) cerr<<"Erreur freq Min>freq Max in LightSrcMapPowerInband"<<endl;
   
   try {        pMap=new SphereHEALPix<r_4> (nlat);     }
   catch(bad_alloc) {
      cerr<<"Memory booking error LightSrcMapPowerInband, nlat= "<<nlat<<endl;
      exit(-1); 
   }
   cout<<"On lit la carte "<<filename<<endl;  
   
  FITS_SphereHEALPix<r_4> FSphere(*pMap);
  FSphere.Read(filename);
}

LightSrcMapPowerInband::~LightSrcMapPowerInband()
{       delete pMap;
}