In [1]:

    

%matplotlib inline
import numpy as np
import sys
import corner
import pylab as plt


    

In [2]:

    

path = "/home/desvignes/projets/1906/new_data/results/inc-180/results-inc-om-fixed/"
data = np.loadtxt(path+"/chains-post_equal_weights.dat", usecols=(0,1,2,3,4))


    

In [ ]:

    

#figure = corner.corner(data, labels=[r"$\alpha$", r"$\delta$", r"$\phi_0$", r"$\psi_{00}$", r"$\Omega$"],
#                        quantiles=[0.16, 0.5, 0.84], bins=75, plot_datapoints=False,
#                        show_titles=True, label_kwargs={"fontsize": 14},
#                       levels=1.0 - np.exp(-0.5 * np.arange(1, 3.1, 1.) ** 2),
#                       hist2d_kwargs={"fontsize": 14}#
#                        #levels=np.array([0.682689492137086, 0.954499736103642, 0.997300203936740])
#                                          )

figure = corner.corner(data, labels=[r"$\alpha$", r"$\delta$", r"$\Phi_0$", r"$\psi_{00}$", r"$\phi_0$", r"$\phi_1$", r"$\phi_{13}$"], 
                         #quantiles=[0.16, 0.5, 0.84],
                         show_titles=False, title_args={"fontsize": 12},
                          bins=64, plot_datapoints=False, levels=1.0 - np.exp(-0.5 * np.arange(1, 3.1, 1.) ** 2))

args = {}
args['fontsize'] = 14

#figure = corner.corner(data, labels=[r"$\alpha$", r"$\delta$", r"$\Phi_0$", r"$\psi_{00}$"], 
#                         #quantiles=[0.16, 0.5, 0.84],
#                         show_titles=False, title_args={"fontsize": 12},
#                           range=[(98.7,99.8), (50,130), (215,315), (-50,70)],
#                          bins=96, plot_datapoints=False, levels=1.0 - np.exp(-0.5 * np.arange(1, 3.1, 1.) ** 2),
#                         label_kwargs=args)

#figure = triangle.corner(data, labels=[r"$\alpha$", r"$\delta$", r"$\phi_0$", r"$\psi_{00}$", r"$\Omega$"],
#                         quantiles=[0.16, 0.5, 0.84],
#                        show_titles=True, title_args={"fontsize": 12})


    

In [11]:

    

figure.savefig('globalRVMnest.eps', format='eps')


    

In [ ]:

    




    

In [ ]:

    

import read_TN
read_TN.makesubplot2d(data[:,4], data[:,5], bins=184, xlab=r"Precession rate, $\Omega$ (deg yr^${-1}$)", ylab="Inclination angle, $i$ (deg)")
plt.xlim(0,5)
plt.ylim(45,100)
plt.show()


    

In [ ]: