diff options
author | Aaron LI <aaronly.me@outlook.com> | 2017-02-20 12:26:17 +0800 |
---|---|---|
committer | Aaron LI <aaronly.me@outlook.com> | 2017-02-20 12:26:17 +0800 |
commit | 4ea7a05ea9a7352602f1f48a860fd81c36e8bc04 (patch) | |
tree | beab7ec18d48c3e2093cd35fd8c79bd66f604a03 /mass_profile/dump_fit_qdp.cpp | |
parent | 9cec16d87f6dc2e0b34b605d88d0837a4a48d18c (diff) | |
download | chandra-acis-analysis-4ea7a05ea9a7352602f1f48a860fd81c36e8bc04.tar.bz2 |
Rename mass_profile to src; Add install & uninstall to Makefile
Diffstat (limited to 'mass_profile/dump_fit_qdp.cpp')
-rw-r--r-- | mass_profile/dump_fit_qdp.cpp | 55 |
1 files changed, 0 insertions, 55 deletions
diff --git a/mass_profile/dump_fit_qdp.cpp b/mass_profile/dump_fit_qdp.cpp deleted file mode 100644 index 85307d1..0000000 --- a/mass_profile/dump_fit_qdp.cpp +++ /dev/null @@ -1,55 +0,0 @@ -#include "dump_fit_qdp.hpp" - -namespace opt_utilities -{ - const static double kpc=3.086E21; - void dump_sbp_beta(std::ostream& os,fitter<double,double,std::vector<double>,double,std::string>& f,double cm_per_pixel,const std::vector<double>& r,const std::vector<double>& y,const std::vector<double>& ye) - { - os<<"read serr 1 2"<<std::endl; - os<<"skip single"<<std::endl; - os<<"la x \"radius (kpc)\""<<std::endl; - os<<"la y \"surface brightness (cts s\\u-1\\d pixel\\u-2\\d)\""<<std::endl; - os<<"li on 2"<<std::endl; - for(size_t i=1;i<r.size();++i) - { - os<<(r[i]+r[i-1])/2*cm_per_pixel/kpc<<"\t"<<(r[i]-r[i-1])/2*cm_per_pixel/kpc<<"\t"<<y[i-1]<<"\t"<<ye[i-1]<<std::endl; - } - os<<"no no no"<<std::endl; - std::vector<double> p=f.get_all_params(); - for(size_t i=1;i<r.size();++i) - { - double x=(r[i]+r[i-1])/2; - os<<x*cm_per_pixel/kpc<<"\t"<<0<<"\t"<<f.eval_model_raw(x,p)<<"\t"<<0<<std::endl; - } - } - - void dump_rho_beta(std::ostream& os,fitter<std::vector<double>,std::vector<double>,std::vector<double>,double,std::string>& f,double cm_per_pixel,const std::vector<double>& r,const std::vector<double>& sbps,const std::vector<double>& sbpe) - { - os<<"read serr 1 2"<<std::endl; - os<<"skip single"<<std::endl; - os<<"la x \"radius (kpc)\""<<std::endl; - os<<"la y \"density (cm\\u-3\\d)\""<<std::endl; - os<<"li on 2"<<std::endl; - - for(size_t i=1;i<r.size();++i) - { - double x=(r[i]+r[i-1])/2; - double y=sbps[i-1]; - double ye=sbpe[i-1]; - os<<x*cm_per_pixel/kpc<<"\t0\t"<<y<<"\t"<<ye<<std::endl; - } - - os<<"no no no"<<std::endl; - std::vector<double> p=f.get_all_params(); - std::vector<double> mv=f.eval_model_raw(r,p); - for(size_t i=1;i<r.size();++i) - { - double x=(r[i]+r[i-1])/2; - double y=mv[i-1]; - double ye=0; - os<<x*cm_per_pixel/kpc<<"\t0\t"<<y<<"\t"<<ye<<std::endl; - } - - } - -}; |