Exception: parameter name not found: omeganuh2 for Planck Chains Using GetDist
-
- Posts: 3
- Joined: May 31 2018
- Affiliation: TSSFL
- Contact:
Exception: parameter name not found: omeganuh2 for Planck Ch
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
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
-
- Posts: 1943
- Joined: September 23 2004
- Affiliation: University of Sussex
- Contact:
Re: Exception: parameter name not found: omeganuh2 for Planc
What code? The neutrino mass is fixed in those chains.
-
- Posts: 3
- Joined: May 31 2018
- Affiliation: TSSFL
- Contact:
Re: Exception: parameter name not found: omeganuh2 for Planc
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')
But, If I do not parse parameter omeganuh2, here is the output:
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?
-
- Posts: 1943
- Joined: September 23 2004
- Affiliation: University of Sussex
- Contact:
Re: Exception: parameter name not found: omeganuh2 for Planc
omeganuh2 is not varied in the chain because the neutrino mass is fixed.
-
- Posts: 3
- Joined: May 31 2018
- Affiliation: TSSFL
- Contact:
Re: Exception: parameter name not found: omeganuh2 for Planc
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?Antony Lewis wrote:omeganuh2 is not varied in the chain because the neutrino mass is fixed.