I am an undergraduate student currently learning to use CosmoMC for an internship. To be fair, this is my first contact with both Fortran and Python and I have a few questions that I hope are not too stupid.
1) My first goal was to plot the experimental constraints for Omegal vs Omegam using the LCDM model on SNLS data. In order to do so, I sampled the chains using Cosmomc then ran them through getdist to do the analysis and generate the plots. But when I tried to run the python script using "python my_file_2D.py", I had the following issues concerning the tick labels and tight layout parameters:
Code: Select all
Traceback (most recent call last):
File "python_plots/bin_run/SNLS_100k_bin20_sc2_2D.py", line 10, in <module>
g.plots_2d(roots,param_pairs=pairs,filled=True)
File "/net/dapserv1/home/manip/mnt/snls/snls/CosmoMC/cosmomc/python/GetDistPlots.py", line 540, in plots_2d
self.plot_2d(roots, param_pair=pair, filled=filled, add_legend_proxy=i == 0)
File "/net/dapserv1/home/manip/mnt/snls/snls/CosmoMC/cosmomc/python/GetDistPlots.py", line 333, in plot_2d
self.setAxes(param_pair, **ax_args)
File "/net/dapserv1/home/manip/mnt/snls/snls/CosmoMC/cosmomc/python/GetDistPlots.py", line 369, in setAxes
self.setAxisProperties(ax.xaxis, True, prune)
File "/net/dapserv1/home/manip/mnt/snls/snls/CosmoMC/cosmomc/python/GetDistPlots.py", line 361, in setAxisProperties
tick_params(axis='both', which='major', labelsize=self.settings.axes_fontsize)
NameError: global name 'tick_params' is not defined
Code: Select all
Traceback (most recent call last):
File "python_plots/bin_run/SNLS_100k_bin20_sc2_2D.py", line 10, in <module>
g.plots_2d(roots,param_pairs=pairs,filled=True)
File "/net/dapserv1/home/manip/mnt/snls/snls/CosmoMC/cosmomc/python/GetDistPlots.py", line 542, in plots_2d
self.finish_plot([legend_labels, roots][legend_labels is None], legend_ncol=legend_ncol)
File "/net/dapserv1/home/manip/mnt/snls/snls/CosmoMC/cosmomc/python/GetDistPlots.py", line 484, in finish_plot
else: tight_layout()
NameError: global name 'tight_layout' is not defined
Do you know how I can fix this to get the proper tick_labels to be processed?
2) Using only SNLS data, I managed to generate many different plots. I am trying to get smoothed contours in order to have a plot that can be compared to an ellipse, yet I find it hard to get the proper shape.
According to various tests, it appears to me that two things seem to influence the appearance of my plots : the effective number of samples (that is given when launching getdist on a chain) and the num_bins parameter in getdist (and the smoothing scale parameter related to it).
From what I get, the effective number of samples seems to vary randomly from one chain to another and to be only related to the dataset used. Am I correct here or is there a way to increase this number on a given dataset by adjusting one of the input parameters of cosmomc? Apparently, it seems to be unrelated to the number of samples in the chain.
As for the num_bins parameter, reducing it quite logically increases the smoothing while diminishing the accuracy of the result. Though I can get some decent-looking almost elliptical plots by adjusting this parameter, I feel this approach to be too artificial and I fear that it lacks proper scientific rigor.
So, am I limited here by the SNLS dataset or is there a way to improve the quality of my plots?
3) And now for something completly different... I would have liked to try to incorporate the Gallileon model in CosmoMC. This model uses the WMAP Distance Priors. I have looked around in the code related to WMAP but have not found any reference to them. Are they currently implemented in CosmoMC or should I add them myself?
Thank you very much for your time,
Jacques Laurens