Commit 07e2f643 authored by Plaszczynski Stephane's avatar Plaszczynski Stephane
Browse files

add access to z/rs rec in classengine

parent 703f62a5
......@@ -127,7 +127,11 @@ public:
bool hasPk() const {return _hasPk;}
//for BAO
//recombination
inline double z_rec() const {return th.z_rec;}
inline double rs_rec() const {return th.rs_rec;}
//baryon drag
inline double z_drag() const {return th.z_d;}
inline double rs_drag() const {return th.rs_d;}
......
......@@ -59,6 +59,10 @@ public:
//forward
bool hasPk() const {return klass->hasPk();}
double get_Pk(double k, double z) { return klass->get_Pk(k,z);}
double get_PkNL(double k, double z) { return klass->get_PkNL(k,z);}
inline double z_drag() const {return klass->z_drag();}
inline double rs_drag() const {return klass->rs_drag();}
double get_Dv(double z) {return klass->get_Dv(z);}
......@@ -74,9 +78,9 @@ public:
double theta() const {return klass->theta();}
double YHe() const {return klass->YHe();}
double Omega_m() const {return klass->Omega_m();}
double get_Pk(double k,double z) {return klass->get_Pk(k,z);}
double get_PkNL(double k,double z) {return klass->get_PkNL(k,z);}
double z_rec() const {return klass->z_rec();}
double rs_rec() const {return klass->rs_rec();}
void writeCls(std::ostream &o) {klass->writeCls(o);}
......
......@@ -69,6 +69,9 @@ public:
virtual double get_Pk(double k, double z) { return undef("getpk");}
virtual double get_PkNL(double k, double z) { return undef("getpknl");}
virtual double z_rec() const {return undef("zrec");}
virtual double rs_rec() const {return undef("rs(rec)");}
//access thtought keyword: convention is that kw is name of previous methods without "get_"
// z values can be specified by appending them in parenthesis, eg. "H(0.57)" for get_H(0.57).
virtual double get(const std::string& kw);
......
......@@ -77,6 +77,7 @@ int main(int argc,char** argv){
//YHE
cout << " YHe=" << e->YHe() << endl << endl;;
cout << "recombination: zrec=" << e->z_rec() << " rs(zrec)=" << e->rs_rec() << endl;
//BAO
double rfid=149.28 ; // from BOSS DR11
double rd=e->rs_drag();
......
......@@ -106,6 +106,7 @@ int main(int argc,char** argv){
//YHE
cout << " YHe=" << e->YHe() << endl << endl;;
cout << "recombination: zrec=" << e->z_rec() << " rs(zrec)=" << e->rs_rec() << endl;
//BAO
cout << "BAO: zd=" << e->z_drag() << " rs(zd)=" << e->rs_drag() << endl;
......
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