Python readchains 예제들

예제 #1

파일: analyse_mcmc_lya.py 프로젝트: jchamilton75/MySoft

#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,')

예제 #2

파일: test_sn1a.py 프로젝트: jchamilton75/MySoft

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)

예제 #3

파일: analyse_mcmc_lya_Halone.py 프로젝트: jchamilton75/MySoft

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)

예제 #4

파일: analyse_mcmc_lyaFev2014.py 프로젝트: jchamilton75/MySoft

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)