#include #include #include "toi.h" #include "toiprocessor.h" #include "fitstoirdr.h" #include "fitstoiwtr.h" #include "toimanager.h" #include "toiseqbuff.h" #include "sambainit.h" #include "toi2map.h" #include "fitsspherehealpix.h" #include "timing.h" void usage(void); void usage(void) { cout<<"tsttoi2map [-h] [-p lp] [-s samplemin,samplemax] [-w data_window_size]"<=sdeb) mgr->setRequestedSample(sdeb,sfin); else {cout<<"Bad sample interval "<=narg) {usage(); exit(-2);} char * fitsin = arg[optind]; string const fitsphout = arg[optind+1]; string fitsphwout = ""; if(optind+2>>> tsttoi2map:"<* sph = new SphereHEALPix(nlat); cout<<"SphereHEALPix: Type de map : "<TypeOfMap()<NbPixels()<* wsph = NULL; if(fitsphwout.size()>0) { wsph = new SphereHEALPix; cout<<"SphereHEALPix Weight: Type de map : "<TypeOfMap()<NbPixels()<setDebugLevel(1); rfits.addOutput(label_alpha,toialphain); toi2m.addInput("AlphaIn",toialphain); TOISeqBuffered * toideltain = new TOISeqBuffered("toi_delta_in",width); if(lp) toideltain->setDebugLevel(1); rfits.addOutput(label_delta,toideltain); toi2m.addInput("DeltaIn",toideltain); TOISeqBuffered * toibolin = new TOISeqBuffered("toi_bolo_in",width); if(lp) toibolin->setDebugLevel(1); rfits.addOutput(label_bolomuv,toibolin); toi2m.addInput("BoloIn",toibolin); // Run cout<<"----- FITSReaderTOI::PrintStatus() : -----"<1) for(int jjjj=0;jjjj<5;jjjj++) { cout<<*toialphain; cout<<*toibolin; sleep(2); } mgr->joinAll(); PrtTim("End threads"); // Ecriture de la sphere Healpix sur fits { FitsOutFile sfits(fitsphout); cout<<"tsttoi2map: Creating sphere fits file "<