Exception: parameter name not found: omeganuh2 for Planck Chains Using GetDist

Use of Healpix, camb, CLASS, cosmomc, compilers, etc.
Post Reply
Elimboto Yohana
Posts: 3
Joined: May 31 2018
Affiliation: TSSFL
Contact:

Exception: parameter name not found: omeganuh2 for Planck Ch

Post by Elimboto Yohana » May 31 2018

Hi CosmoCoffee members,

I am using GetDist to plot contours for cosmological parameters from Planck chains for base and w, wa parameters. The specific paths to the Planck data that I load are:

./COM_CosmoParams_fullGrid_R2.00/base/plikHM_TTTEEE_lowTEB_lensing/base_plikHM_TTTEEE_lowTEB_lensing_post_BAO_H070p6_JLA'

./COM_CosmoParams_fullGrid_R2.00/base_w_wa/plikHM_TTTEEE_lowTEB_BAO_H070p6_JLA/base_w_wa_plikHM_TTTEEE_lowTEB_BAO_H070p6_JLA_post_lensing'

Every time I run my code, I get this error:

Exception: parameter name not found: omeganuh2

But, the parameter omeganuh2 seems to be included in the chains.

Can you help on how to proceed to solve this problem?

Thank you,

Elimboto

Antony Lewis
Posts: 1364
Joined: September 23 2004
Affiliation: University of Sussex
Contact:

Re: Exception: parameter name not found: omeganuh2 for Planc

Post by Antony Lewis » May 31 2018

What code? The neutrino mass is fixed in those chains.

Elimboto Yohana
Posts: 3
Joined: May 31 2018
Affiliation: TSSFL
Contact:

Re: Exception: parameter name not found: omeganuh2 for Planc

Post by Elimboto Yohana » June 01 2018

Antony Lewis wrote:What code? The neutrino mass is fixed in those chains.

Hi Antony, thank you for reply.

Here is the code:

Code: Select all

#!/usr/bin/env python3

from getdist import loadMCSamples, plots

#Specify location of Planck chains 
path_planck_base_w_wa_TTTEEE_lowTEB_lensing_post_BAO_H0_JLA = '/home/tssfl/Desktop/COM_CosmoParams_fullGrid_R2.00/base_w_wa/plikHM_TTTEEE_lowTEB_BAO_H070p6_JLA/base_w_wa_plikHM_TTTEEE_lowTEB_BAO_H070p6_JLA_post_lensing'


#Load Planck chains
samples_base_w_wa_TTTEEE_lowTEB_lensing_post_BAO_H0_JLA = loadMCSamples("{}".format(path_planck_base_w_wa_TTTEEE_lowTEB_lensing_post_BAO_H0_JLA))

getplot = plots.getSubplotPlotter()
getplot.triangle_plot([samples_base_w_wa_TTTEEE_lowTEB_lensing_post_BAO_H0_JLA],
                  ['omegabh2', 'omegach2', 'w', 'wa', 'logA', 'H0', 'ns', 'omeganuh2'],
                  filled=True, 
                  legend_labels = ['Samples_base_w_wa_TTTEEE_lowTEB_lensing_post_BAO_H0_JLA'])

getplot.export('Samples_base_w_wa_TTTEEE_lowTEB_lensing_post_BAO_H0_JLA.png')
When I run this code, I get the error message: Exception: parameter name not found: omeganuh2


But, If I do not parse parameter omeganuh2, here is the output:


Image

omeganuh2 is listed in the parameters file and denoted by omeganuh2* (just like H0) to indicate that it is derived, so what is wrong with what I am doing?

Antony Lewis
Posts: 1364
Joined: September 23 2004
Affiliation: University of Sussex
Contact:

Re: Exception: parameter name not found: omeganuh2 for Planc

Post by Antony Lewis » June 01 2018

omeganuh2 is not varied in the chain because the neutrino mass is fixed.

Elimboto Yohana
Posts: 3
Joined: May 31 2018
Affiliation: TSSFL
Contact:

Re: Exception: parameter name not found: omeganuh2 for Planc

Post by Elimboto Yohana » June 01 2018

Antony Lewis wrote:omeganuh2 is not varied in the chain because the neutrino mass is fixed.
Does this mean it is not possible to plot omeganuh2 altogether with the rest of the parameters or there is a workaround to do so?

Post Reply