Commit a80c519d authored by Plaszczynski Stephane's avatar Plaszczynski Stephane
Browse files

better dumps for clik files

parent c815ece6
......@@ -74,8 +74,8 @@ application testKlass -group=test -s=../src/camel test/testKlass.cc
application testCLHEP -group=test -s=../src/camel/test testCLHEP.cc
application testBAO -group=test -s=../src/camel/test testBAO.cc
application testHiLLiPOP -group=test -s=../src/camel/test testHiLLiPOP.cc
application test_clik -group=test -s=../src/camel test/test_clik.cc
application test_lensing -group=test -s=../src/camel test/test_lensing.cc
application dump_clik -group=test -s=../src/camel test/test_clik.cc
application dump_lensing -group=test -s=../src/camel test/test_lensing.cc
application test_jla -group=test -s=../src/camel test/test_jla.cc
application testPico -group=test -s=../src/camel/pico/exec testPico.cc
......
......@@ -102,11 +102,11 @@ ClikLensingLikelihood::computePPLikelihood(const std::vector<unsigned int>& l,
std::vector<double>& nuisance){
//bretelle si fid
if (_fid) return computePPLikelihood(l,clphiphi,nuisance);
int lmin=l[0];
const int lmax=l.back();
//const int lmax=l.back();
//planck_assert(_lmax >= lmax(),"too low l given in CLIKlensing likelihood");
vector<double> fullCl;
......@@ -201,7 +201,7 @@ ClikLensingLikelihood::computePPLikelihood(const std::vector<unsigned int>& l,
err= &_err;
int lmin=l[0];
const int lmax=l.back();
//const int lmax=l.back();
vector<double> fullCl;
cltype t=PP;
......
......@@ -32,7 +32,7 @@ public:
enum cltype {PP=0,TT,EE,BB,TE,TB,EB}; //ordre clik
//constructors
ClikLensingLikelihood(char* hdf5file,bool use_fiducial=true);
ClikLensingLikelihood(char* hdf5file,bool use_fiducial=false);
//accessors
inline int getLmax(cltype t) const {return _lmax[t];}
......@@ -60,6 +60,8 @@ public:
// destructor
~ClikLensingLikelihood();
void setFiducial(bool b){ _fid=b;}
private:
clik_lensing_object* clikid;
......
......@@ -18,9 +18,7 @@ int main(int argc,char** argv){
//CLIK NOW
ClikLikelihood clik(argv[1]);
vector<double> nui(clik.nuisanceNames().size(),1);
clik.dump();
return 0;
}
......@@ -5,6 +5,8 @@
#include"Parser.hh"
#include<string>
#include<fstream>
#include<iomanip>
#include"Class/ClassEngine.hh"
#include"Lensing/ClikLensingLikelihood.hh"
#include<vector>
#include<string>
......@@ -35,17 +37,15 @@ int main(int argc,char** argv){
//CLIK NOW
vector<double> nuisance(1);
nuisance[0]=1;
ClikLensingLikelihood* lensing= new ClikLensingLikelihood(argv[1],true);
double lkl=lensing->computePPLikelihood(ell,clpp,nuisance);
cout << ">>>>>>>>>CLIK LKL (fid). =" << lkl << " in " << timer.partial()<< endl;
//use TT
ClikLensingLikelihood* lensing= new ClikLensingLikelihood(argv[1]);
lensing->dump();
//call to PP only (ie use fidicuial, default Cls in clik)
double lkl=lensing->computePPLikelihood(ell,clpp,nuisance);
cout << fixed << setprecision(5) << ">>>>>>>>>CLIK LKL (fid). =" << lkl << " in " << timer.partial()<< "s" << endl;
delete lensing;
// lensing=new PlensLensingLikelihood(argv[1]);
// lkl=lensing->computePPLikelihood(ell,clpp,cltt,nuisance);
// cout << ">>>>>>>>>LKL PLENS=" << lkl << " in " << timer.partial()<< endl;
// delete lensing;
return 0;
}
Markdown is supported
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment