CosmoCoffee Forum Index CosmoCoffee

 FAQFAQ   SearchSearch  MemberlistSmartFeed   MemberlistMemberlist    RegisterRegister 
   ProfileProfile   Log inLog in 
Arxiv New Filter | Bookmarks & clubs | Arxiv ref/author:

Planck 2015 likelihood: Extracting covariance
Post new topic   Reply to topic    CosmoCoffee Forum Index -> Computers and software
View previous topic :: View next topic  
Author Message
Luke Hart

Joined: 13 Jul 2015
Posts: 30
Affiliation: University of Manchester

PostPosted: January 16 2018  Reply with quote

Hello everyone,

So I've spent the last few days trying to work out how to use the code supplied in plc_2.0 (COM_Likelihood_Code...) in the Planck 2015 Likelihood code with the plik data to construct a data covariance matrix. This is including the off-diagonal components due to multipole mixing as well.

Does anyone have any idea where to start? I've been looking at cliklike.F90 in the CosmoMC code (most recent version) as well, but there doesn't seem to be any clear links. I've only seen covariance codes set up for the foreground routines, in the cmb_only component of the code and as I say I don't really understand how to use it.

To be clear: trying to create a covariance matrix between the Cls including binned-ls and the mixing between different multipoles.

I've been looking at the clik header files and the computational aspects included in clik_compute but I cannot find the routines where the covariance is calculated??

Please any help would be appreciated :)

Back to top
View user's profile  
Graeme Addison

Joined: 17 Jul 2014
Posts: 17
Affiliation: Johns Hopkins University

PostPosted: February 14 2018  Reply with quote

Not sure this is still useful but the Cl-Cl' covariance can be accessed by reading in the criterion_gauss_mat file that lives in the plik/plik_dx11dr2_HM_v18_TT.clik/clik/lkl_0/ folder. In python you can do this for example by using

This will give you an array of floats that can be reshaped into a 765x765 covariance matrix - on each side the first 136 entries are the bandpower bins for 100x100 GHz, the next 199 are 143x143, then 215 for 143x217, and finally another 215 for 217x217 GHz.

This matrix has the signal (cosmic variance) and noise contributions included as described in the Planck 2015 likelihood paper. Note that the dependence of the CMB cosmic variance contribution on the actual cosmological parameters is ignored in the Planck likelihood, meaning the covariance matrix is fixed rather than having to be updated as the values of cosmological params change.

EDIT: Oops, sorry, this file contains the *inverse* covariance matrix.
Back to top
View user's profile  
Display posts from previous:   
Post new topic   Reply to topic    CosmoCoffee Forum Index -> Computers and software All times are GMT + 5 Hours
Page 1 of 1

Jump to:  
You cannot post new topics in this forum
You cannot reply to topics in this forum
You cannot edit your posts in this forum
You cannot delete your posts in this forum
You cannot vote in polls in this forum

Powered by phpBB © 2001, 2005 phpBB Group. Sponsored by WordWeb online dictionary and dictionary software.