| View previous topic :: View next topic |
| Author |
Message |
kumar pandey
Joined: 26 Nov 2008 Posts: 3 Affiliation: Inter University Centre for Astronomy and Astrophysics
|
|
| Back to top |
|
 |
Paniez Paykari
Joined: 08 Jun 2007 Posts: 14 Affiliation: CEA
|
Posted: December 15 2010 |
|
|
Hi All,
Does the s in the planck.dat file have to be necessarily smooth? Or can we have a simulated from a CMB map? I appreciate the help.
Thanks
Paniez |
|
| Back to top |
|
 |
Antony Lewis
Joined: 23 Sep 2004 Posts: 675 Affiliation: University of Sussex
|
Posted: December 15 2010 |
|
|
| You can use a simulation, just don't expect the posteriors to peak at exactly the simulation input values. |
|
| Back to top |
|
 |
Paniez Paykari
Joined: 08 Jun 2007 Posts: 14 Affiliation: CEA
|
Posted: December 15 2010 |
|
|
Thank you. So, that means I should expect some bias? I can trust the errors though?
Paniez |
|
| Back to top |
|
 |
Alex Conley
Joined: 08 Feb 2005 Posts: 11
|
Posted: April 23 2011 |
|
|
I'm a bit confused – why, in the example code, does
one write CTT+Nl and CEE+noiseFac*Nl, but
does not include any noise term when writing out
CTE? |
|
| Back to top |
|
 |
Antony Lewis
Joined: 23 Sep 2004 Posts: 675 Affiliation: University of Sussex
|
Posted: April 23 2011 |
|
|
| Just assuming the simple case where there is no noise correlation between Q/U and T. |
|
| Back to top |
|
 |
Alessio Notari
Joined: 26 Nov 2005 Posts: 7 Affiliation: McGill U.
|
Posted: July 11 2011 |
|
|
Hello,
I am trying to use the all_l_exact format on CosmoMC with the above procedure with some C_l table produced by the Healpix package anafast, and it works well :) .
One simple question: does it work only for full-sky data or should I do something more complicated (such as using planck_like.f90) when I put a mask with anafast on my sky maps?
And by the way, what is the planck_like.f90 routine doing exactly and how do I use it? (I cannot find any instructions on the web) ...
thanks a lot for your help!
Alessio |
|
| Back to top |
|
 |
Antony Lewis
Joined: 23 Sep 2004 Posts: 675 Affiliation: University of Sussex
|
Posted: July 12 2011 |
|
|
Planck_like implements the approximation of 0801.0554 as well as simple schemes like full-sky and Gaussian. You probably should use something that accounts for the covariance when using cut-sky data.
For sample data files see the bottom of
http://cosmologist.info/weightmixer/
or the python exact simulation script supplied with most recent cosmomc. |
|
| Back to top |
|
 |
Alessio Notari
Joined: 26 Nov 2005 Posts: 7 Affiliation: McGill U.
|
Posted: July 12 2011 |
|
|
Thanks for the answer.
Though, I am not sure if I understand the way to preceed. How do I tell to cosmomc that there is nontrivial covariance due to cut-sky?
Moreover, in the all_lexact format why is there an fsky column if this is valid only for full sky?
cheers |
|
| Back to top |
|
 |
Alessio Notari
Joined: 26 Nov 2005 Posts: 7 Affiliation: McGill U.
|
|
| Back to top |
|
 |
Jason Dossett
Joined: 19 Mar 2010 Posts: 21 Affiliation: The University of Texas at Dallas
|
Posted: February 28 2012 |
|
|
| Antony Lewis wrote: | Planck_like implements the approximation of 0801.0554 as well as simple schemes like full-sky and Gaussian. You probably should use something that accounts for the covariance when using cut-sky data.
For sample data files see the bottom of
http://cosmologist.info/weightmixer/
or the python exact simulation script supplied with most recent cosmomc. |
Hi,
I am trying to use the makePerfectForcast.py python script to create a new dataset for some future simulated data. I understand how to produce the input files:
lensedTotClFileRoot and lensPotential.
I am however unsure how exactly to construct the reconNoise input file. Any help would be greatly appreciated.
Thanks,
Jason |
|
| Back to top |
|
 |
Antony Lewis
Joined: 23 Sep 2004 Posts: 675 Affiliation: University of Sussex
|
Posted: February 28 2012 |
|
|
| You only need reconNoise if you want to include CMB lensing reconstruction (in which case it is the reconstruction noise variance, which you can calculate for isotropic full sky case using standard result involving 3j sum). |
|
| Back to top |
|
 |
Jason Dossett
Joined: 19 Mar 2010 Posts: 21 Affiliation: The University of Texas at Dallas
|
Posted: March 03 2012 |
|
|
| Thanks Antony. |
|
| Back to top |
|
 |
|