#xterm -e python ~/Python/Boss/McMc/mcmc_launcher.py owcdm BAO_DR7 & #xterm -e python ~/Python/Boss/McMc/mcmc_launcher.py owcdm BAO_Anderson & #xterm -e python ~/Python/Boss/McMc/mcmc_launcher.py owcdm BAO & #xterm -e python ~/Python/Boss/McMc/mcmc_launcher.py owcdm LyaDR11+BAO & #xterm -e python ~/Python/Boss/McMc/mcmc_launcher.py owcdm LyaDR11 & #xterm -e python ~/Python/Boss/McMc/mcmc_launcher.py owcdm planck & #xterm -e python ~/Python/Boss/McMc/mcmc_launcher.py owcdm planck+BAO & #xterm -e python ~/Python/Boss/McMc/mcmc_launcher.py owcdm planck+LyaDR11 & #xterm -e python ~/Python/Boss/McMc/mcmc_launcher.py owcdm planck+BAO+LyaDR11 & rep='/Users/hamilton/SDSS/LymanAlpha/JCMC_Chains/' ext='.db' ######################## Ok ################################ model='olambdacdm' planck=mcmc.readchains(rep+model+'-'+'planck'+ext,add_extra=True) planck_bao=mcmc.readchains(rep+model+'-'+'planck+BAO'+ext,add_extra=True) planck_bao_lya=mcmc.readchains(rep+model+'-'+'planck+BAO+LyaDR11'+ext,add_extra=True) planck_lya=mcmc.readchains(rep+model+'-'+'planck+LyaDR11'+ext,add_extra=True) lya=mcmc.readchains(rep+model+'-'+'LyaDR11'+ext,add_extra=True) beutler=mcmc.readchains(rep+model+'-'+'BAO_Beutler'+ext,add_extra=True) dr7=mcmc.readchains(rep+model+'-'+'BAO_DR7'+ext,add_extra=True) anderson=mcmc.readchains(rep+model+'-'+'BAO_Anderson'+ext,add_extra=True) bao=mcmc.readchains(rep+model+'-'+'BAO'+ext,add_extra=True) lya_bao=mcmc.readchains(rep+model+'-'+'LyaDR11+BAO'+ext,add_extra=True) ################################################### clf() plot(lya['om_ol'],lya['rssqrtomh2'],'k,')
import scipy import pickle xterm -e python ~/Python/Boss/McMc/mcmc_launcher.py flatwcdm planck & xterm -e python ~/Python/Boss/McMc/mcmc_launcher.py flatwcdm planck+BAO & xterm -e python ~/Python/Boss/McMc/mcmc_launcher.py flatwcdm planck+BAO+hRiess & xterm -e python ~/Python/Boss/McMc/mcmc_launcher.py flatwcdm planck+BAO+LyaDR11+hRiess & xterm -e python ~/Python/Boss/McMc/mcmc_launcher.py flatwcdm planckHp & xterm -e python ~/Python/Boss/McMc/mcmc_launcher.py flatwcdm planckHp+BAO & xterm -e python ~/Python/Boss/McMc/mcmc_launcher.py flatwcdm planckHp+BAO+hRiess & ############################ Flat w ################################### model='flatwcdm' rep='/Users/hamilton/SDSS/LymanAlpha/JCMC_Chains/' ext='.db' planck=mcmc.readchains(rep+model+'-'+'planck'+ext) planck_bao=mcmc.readchains(rep+model+'-'+'planck+BAO'+ext) planck_lya=mcmc.readchains(rep+model+'-'+'planck+LyaDR11'+ext) planck_bao_lya=mcmc.readchains(rep+model+'-'+'planck+BAO+LyaDR11'+ext) planck_bao_hriess=mcmc.readchains(rep+model+'-'+'planck+BAO+hRiess'+ext) planck_bao_lya_hriess=mcmc.readchains(rep+model+'-'+'planck+BAO+LyaDR11+hRiess'+ext) planckHp=mcmc.readchains(rep+model+'-'+'planckHp'+ext) planckHp_bao=mcmc.readchains(rep+model+'-'+'planckHp+BAO'+ext) planckHp_bao_hriess=mcmc.readchains(rep+model+'-'+'planckHp+BAO+hRiess'+ext) ################################################### clf() sm=6 xlim(0.,0.6) ylim(-2.,-0.3)
from McMc import cosmo_utils import scipy import pickle #### run in a shell xterm -e python ~/Python/Boss/McMc/mcmc_launcher.py owcdm BAO & xterm -e python ~/Python/Boss/McMc/mcmc_launcher.py owcdm LyaDR11_Halone+BAO & xterm -e python ~/Python/Boss/McMc/mcmc_launcher.py owcdm LyaDR11_Halone & rep='/Users/hamilton/SDSS/LymanAlpha/JCMC_Chains/' ext='.db' ######################## OwCDM ################################ model='owcdm' lya=mcmc.readchains(rep+model+'-'+'LyaDR11_Halone'+ext,add_extra=True) bao=mcmc.readchains(rep+model+'-'+'BAO'+ext,add_extra=True) lya_bao=mcmc.readchains(rep+model+'-'+'LyaDR11_Halone+BAO'+ext,add_extra=True) ################################################### reload(mcmc) clf() vars=['omega_M_0', 'omega_lambda_0','c_H0rs'] limits=[[0.,0.5],[0,1.2],[25,35]] doit=[True,True,True] a0=mcmc.matrixplot(bao,vars,'yellow',8,limits=limits,doit=doit) a1=mcmc.matrixplot(lya,vars,'purple',4,limits=limits,doit=doit) a2=mcmc.matrixplot(lya_bao,vars,'blue',4,limits=limits,doit=doit) subplot(3,3,3)
xterm -e python ~/Python/Boss/McMc/mcmc_launcher.py olambdacdm LyaDR11_HRiessPlanck_obh2BBN & xterm -e python ~/Python/Boss/McMc/mcmc_launcher.py flatwcdm LyaDR11_HRiessPlanck_obh2BBN & xterm -e python ~/Python/Boss/McMc/mcmc_launcher.py owcdm LyaDR11_HRiessPlanck_obh2BBN & xterm -e python ~/Python/Boss/McMc/mcmc_launcher.py olambdacdm BAO_HRiessPlanck_obh2BBN & xterm -e python ~/Python/Boss/McMc/mcmc_launcher.py flatwcdm BAO_HRiessPlanck_obh2BBN & xterm -e python ~/Python/Boss/McMc/mcmc_launcher.py owcdm BAO_HRiessPlanck_obh2BBN & rep='/Users/hamilton/SDSS/LymanAlpha/JCMC_Chains/' ext='.db' repfig = '/Users/hamilton/SDSS/LymanAlpha/JCMC_Chains/Fev2014/' ######################## Ok ################################ model='olambdacdm' lya_hp1_obh2p1=mcmc.readchains(rep+model+'-'+'LyaDR11_HPlanck1s_obh2Planck1s'+ext,add_extra=True) lya_hr2_obh2p1=mcmc.readchains(rep+model+'-'+'LyaDR11_HRiess2s_obh2Planck1s'+ext,add_extra=True) lya_hrp_obh2p1=mcmc.readchains(rep+model+'-'+'LyaDR11_HRiessPlanck_obh2Planck1s'+ext,add_extra=True) bao_hrp_obh2p1=mcmc.readchains(rep+model+'-'+'BAO_HRiessPlanck_obh2Planck1s'+ext,add_extra=True) lya_bao_hrp_obh2p1=mcmc.readchains(rep+model+'-'+'LyaDR11_BAO_HRiessPlanck_obh2Planck1s'+ext,add_extra=True) lya_bao=mcmc.readchains(rep+model+'-'+'LyaDR11+BAO'+ext,add_extra=True) ################################################### reload(mcmc) clf() limits=[[0,1],[0,1.5],[0.6,0.8],[0.021,0.023]] vars=['omega_M_0','omega_lambda_0','h','obh2'] a2=mcmc.matrixplot(lya_hrp_obh2p1,vars,'red',8,limits=limits,alpha=0.5) a0=mcmc.matrixplot(lya_hp1_obh2p1,vars,'blue',8,limits=limits,alpha=0.5) a1=mcmc.matrixplot(lya_hr1_obh2p1,vars,'green',8,limits=limits,alpha=0.5)
