#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;
}