source: Sophya/trunk/SigPredictor/AutresASiGPred.cc@ 2929

#include <stdio.h>
#include <stdlib.h>
#include <iostream>
#include "makebolotimeline.h"
#include "bruit.h"
//#include "fitsio.h"

#ifdef __MWERKS__
   #include "macenvvariables.h"
   #define SetNlatLSource (256)         // version debug
   #define NbSAMPLE (2000)                              //      Pour le debug
#else
   # define ERROR_LABEL -32768
   #define SetNlatLSource (1024)                // version de production CXX
   #define NbSAMPLE (2684300)
#endif


#include "alllobe.h"
#include "alllightsources.h"
#include "sigcalctools.h"
#include "allfilter.h"

int MakeGalaxyMap(double Resolution);   // Resolution en degres

int MakeGalaxyMap(double Pas)
{       
// On passe aux cartes de fonds physiques
// Initialisations
   char filename[150]="";
     
//      cerr.setf(ios::scientific);
        cerr<< "C'est parti"<<endl;
        
        // Des Lobes
        LobeGaussien GLobeTrapChan1(8./60.,75.e9,600.e9);  //degrŽs, freqmin, Freqmax
        
#ifndef __MWERKS__
  char* PathInstrumentData=getenv("PathInstrumentData");
#endif

        // Des Filtres
        sprintf(filename, "%sfiltretrap99chan1.dat",PathInstrumentData);
        TabulatedFilter FiltreTrapChan1(filename);      

        // Les Processus Physiques
        LightDiffDust DiffDust(SetNlatLSource);

        // Bolo 1_1 soit bolo11.
        sprintf(filename, "ASigPred_MapGalaxie");
        MakeBoloTimeline MapGalaxie(filename,DiffDust,GLobeTrapChan1,
                FiltreTrapChan1,RayleighJeans);
                
        // On definit les timelines desirŽees
        double alphadeg, deltadeg;
        double alpharad, deltarad;
        
        double NbreDePoint=360*40/(Pas*Pas);
        cout<<"Nbre de point a calculer: "<<NbreDePoint<<endl;
        
        long compteur=0;
        for(alphadeg=0; alphadeg<360.; alphadeg+=Pas)
        {       for(deltadeg=-20.; deltadeg<20.; deltadeg+=Pas)
                {       alpharad= alphadeg/180.*M_PI;
                        deltarad= deltadeg/180.*M_PI;           
                        MapGalaxie.addToStream(alpharad,deltarad);
                        
                        compteur++;
                        if(compteur%5000==0) 
                        {       double pourcent=100.*compteur/NbreDePoint;
                                cout<<"Nbre points calcules: "<<compteur<<endl;
                                cout<<"Travail effectue a :"<<pourcent<<" %"<<endl;
                        }
                }       
        }
        return 0;
}