Ejemplos de Fisher.cont_from_fisher2d en Python

Ejemplo n.º 1

Archivo: matrices_inflation.py Proyecto: jchamilton75/MySoft

####### Inflation : from Dodelson, Kinney, Kolb, astro-ph/9702166v1 
clf()
xlim(0.8,1.1)
ylim(0,0.6)
xlabel('$n_s$',fontsize=20)
ylabel('$r$',fontsize=20)
al = 0.3
aa = fill_between([0.8, 1], [0,0], [0.45, 0],color='b', alpha=al)
text(0.82,0.18,'Small Field',fontsize=20)
bb = fill_between([0.8, 1], [0.45, 0], [1.4,0],color='r', alpha=al)
text(0.85,0.4,'Large Field',fontsize=20)
cc = fill_between([0.8, 1.1], [1.4,-0.7], [1e8,0],color='y', alpha=al)
text(1.02,0.3,'Hybrid',fontsize=20)
nsval = params['scalar_index']
rval = params['tensor_ratio']
a0=Fisher.cont_from_fisher2d(fm_nsr_bicep, [nsval, rval] ,color='k', onesigma=True)
a1=Fisher.cont_from_fisher2d(fm_nsr_qubic, [nsval, rval] ,color='b', onesigma=True)
a2=Fisher.cont_from_fisher2d(fm_nsr_qubicnew, [nsval, rval] ,color='g', onesigma=True)
a3=Fisher.cont_from_fisher2d(fm_nsr_qubicnew6, [nsval, rval] ,color='r', onesigma=True)
legend([a1,a0,a2,a3],['Planck $n_s$ + BICEP2', 'Planck $n_s$ + QUBIC 1 year, $f_{sky}=0.1$', 
	'Planck $n_s$ + QUBIC 1 year, $f_{sky}=0.4$', 
	'Planck $n_s$ + QUBIC 6 modules, $f_{sky}=0.4$'], fontsize=15, frameon=0)
savefig('inflationtype.png')




####### Cross-Check with NewFisher.py
from __future__ import division

import healpy as hp

Ejemplo n.º 2

Archivo: satellites.py Proyecto: jchamilton75/MySoft

####### Inflation : from Dodelson, Kinney, Kolb, astro-ph/9702166v1 
clf()
xlim(0.8,1.1)
ylim(0,0.6)
xlabel('$n_s$',fontsize=20)
ylabel('$r$',fontsize=20)
al = 0.3
aa = fill_between([0.8, 1], [0,0], [0.45, 0],color='b', alpha=al)
text(0.82,0.18,'Small Field',fontsize=20)
bb = fill_between([0.8, 1], [0.45, 0], [1.4,0],color='r', alpha=al)
text(0.85,0.4,'Large Field',fontsize=20)
cc = fill_between([0.8, 1.1], [1.4,-0.7], [1e8,0],color='y', alpha=al)
text(1.02,0.3,'Hybrid',fontsize=20)
nsval = params['scalar_index']
rval = params['tensor_ratio']
a0=Fisher.cont_from_fisher2d(fm_nsr_lo, [nsval, rval] ,color='g', onesigma=True)
a1=Fisher.cont_from_fisher2d(fm_nsr_s4, [nsval, rval] ,color='b', onesigma=True)
a2=Fisher.cont_from_fisher2d(fm_nsr_los4, [nsval, rval] ,color='r', onesigma=True)
a3=Fisher.cont_from_fisher2d(fm_nsr_sat, [nsval, rval] ,color='k', onesigma=True)
a4=Fisher.cont_from_fisher2d(fm_nsr_satplus, [nsval, rval] ,color='m', onesigma=True)
legend([a1,a0,a2,a4],['CMBS4', 'LRS', 
	'CMBS4+LRS', 
	'HRS','HRS++'], fontsize=15, frameon=0)
savefig('inflationtype_satellites._r{0:.2g}'.format(rvalue)+'.png')


clf()
xlim(0.8,1.1)
ylim(rvalue-0.01,rvalue+0.01)
xlabel('$n_s$',fontsize=20)
ylabel('$r$',fontsize=20)