CosmoMC icon indicating copy to clipboard operation
CosmoMC copied to clipboard

SZ cluster counts likelihood: natural vs base 10 log in szcounts.f90

Open borisbolliet opened this issue 3 years ago • 3 comments

Hi,

I think the prior on the scatter parameter sigmaM, in line:

https://github.com/cmbant/CosmoMC/blob/d34a8d34a3f61f58d8dd601519ed459b006323a0/source/szcounts.f90#L2014

should be

(cosmopar%sigmaM-0.173)**2./(2.*0.023**2.)

rather than

(cosmopar%sigmaM-0.075)**2./(2.*0.01**2.)

because I think that the routine uses the parameter sigmaM defined in natural log, as described in the 2015 Planck SZ cluster paper (table 1 of https://arxiv.org/pdf/1502.01597.pdf and footnote c of the table), as can be seen here:

https://github.com/cmbant/CosmoMC/blob/d34a8d34a3f61f58d8dd601519ed459b006323a0/source/szcounts.f90#L1282

(Or maybe I am missing a conversion from natural log to base 10 log...?)

If this needs to be corrected, then these lines:

https://github.com/cmbant/CosmoMC/blob/d34a8d34a3f61f58d8dd601519ed459b006323a0/batch2/SZ_plus_CMB.ini#L14

https://github.com/cmbant/CosmoMC/blob/d34a8d34a3f61f58d8dd601519ed459b006323a0/batch2/SZ_plus_priors.ini#L15

in the input files should also be updated to be with respect to the natural log value of 0.173\pm0.023.

borisbolliet avatar Jun 09 '21 22:06 borisbolliet

I think this code was contributed by @abonaldi who perhaps can comment

cmbant avatar Jun 09 '21 22:06 cmbant

Is there a way I can see the original question posed?

thanks

Anna

Sent from my iPhone

On 9 Jun 2021, at 23:36, Antony Lewis @.***> wrote:



I think this code was contributed by @abonaldihttps://github.com/abonaldi who perhaps can comment

— You are receiving this because you were mentioned. Reply to this email directly, view it on GitHubhttps://github.com/cmbant/CosmoMC/issues/34#issuecomment-858145509, or unsubscribehttps://github.com/notifications/unsubscribe-auth/AG67ED2N7LB3RI6I3KDXN6TTR7UATANCNFSM46M5IOAA.

The SKA Observatory is an inter-governmental organisation and the successor of SKA Organisation, a private limited company by guarantee registered in England and Wales with registered number 07881918, with a registered office of Jodrell Bank Observatory, Lower Withington, Macclesfield, Cheshire, England, SK11 9FT.

This message is intended solely for the addressee and may contain confidential information. If you have received this message in error, please inform the sender, and immediately and permanently delete the email. Do not use, copy or disclose the information contained in this message or in any attachment.

This email has been scanned for viruses and malware, and may have been automatically archived, by Mimecast Ltd. Although SKA Observatory and SKA Organisation have taken reasonable precautions to ensure no viruses are present in this email, neither SKA Observatory nor SKA Organisation accept responsibility for any loss or damage sustained as a result of computer viruses and the recipient must ensure that the email (and attachments) are virus free.

abonaldi avatar Jun 10 '21 15:06 abonaldi

Here on github. if you are looking at the github email, click "view it on github".

cmbant avatar Jun 10 '21 16:06 cmbant