//#define TOISEQBUFFERED #include #include #include #include "toi.h" #include "toiprocessor.h" #include "fitstoirdr.h" #include "fitstoiwtr.h" #include "toimanager.h" #ifdef TOISEQBUFFERED #include "toiseqbuff.h" #else #include "toisegment.h" #endif #include "sambainit.h" #include "flag2map.h" #include "fitsspherehealpix.h" #include "timing.h" #include void usage(void); void usage(void) { cout<<"tstflag2map [-h] [-p lp] [-s samplemin,samplemax] [-w data_window_size]"<= vmin"<=sdeb) mgr->setRequestedSample(sdeb,sfin); //else {cout<<"Bad sample interval "<=narg) {usage(); exit(3);} char * fitsin_point = arg[optind]; char * fitsin_bolo = arg[optind+1]; string const fitsphout = arg[optind+2]; if(optind+3>>> tstflag2map:"<* sph = new SphereHEALPix(nlat); cout<<"SphereHEALPix: Type de map : "<TypeOfMap()<NbPixels()<SizeIndex()<* wsph = NULL; if(fitsphwout.size()>0) { wsph = new SphereHEALPix; cout<<"SphereHEALPix Weight Created"<>>> Using TOISeqBuffered"<>>> Using TOISegmented"<setDebugLevel(1); rfitsp.addOutput(label_coord1,toicoord1in); flag2m.addInput("Coord1In",toicoord1in); #ifdef TOISEQBUFFERED TOISeqBuffered * toicoord2in = new TOISeqBuffered("toi_coord2_in",width); #else TOISegmented * toicoord2in = new TOISegmented("toi_coord2_in",width); #endif // toicoord2in->setDebugLevel(1); rfitsp.addOutput(label_coord2,toicoord2in); flag2m.addInput("Coord2In",toicoord2in); #ifdef TOISEQBUFFERED TOISeqBuffered * toibolin = new TOISeqBuffered("toi_bolo_in",width); #else TOISegmented * toibolin = new TOISegmented("toi_bolo_in",width); #endif // toibolin->setDebugLevel(1); rfitsb.addOutput(label_bolomuv,toibolin); flag2m.addInput("BoloIn",toibolin); // Run cout<<"----- FITSReaderTOI::PrintStatus() : -----"<1) // for(int jjjj=0;jjjj<5;jjjj++) { // cout<<*toicoord1in; // cout<<*toibolin; // sleep(2); // } // Affichage de l'avancement des TOIProcessors ProcSampleCounter stats(flag2m); stats.InfoMessage() = "tstflag2map/Info"; stats.PrintStats(); // Fin des traitements des TOIProcessors mgr->joinAll(); PrtTim("End threads"); // Ecriture de la sphere Healpix sur fits { FitsOutFile sfits(fitsphout,FitsFile::clear); cout<<"tstflag2map: Creating sphere fits file "<