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

[borisbolliet]

Hi,

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

CosmoMC/source/szcounts.f90

Line 2014 in d34a8d3

palpha*(cosmopar%alpha-1.789)**2./(2.*0.084**2.) + psigma*(cosmopar%sigmaM-0.075)**2./(2.*0.01**2.) !cosmomc_sz

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:

CosmoMC/source/szcounts.f90

Line 1282 in d34a8d3

arg=((lny-mu)/(sqrt(2.)*sigmaM))

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

If this needs to be corrected, then these lines:

CosmoMC/batch2/SZ_plus_CMB.ini

Line 14 in d34a8d3

param[scatter_SZ]=0.0750000 0.0220000 0.122000 0.00500000 0.00250000

CosmoMC/batch2/SZ_plus_priors.ini

Line 15 in d34a8d3

param[scatter_SZ]=0.0750000 0.0220000 0.122000 0.00500000 0.00250000

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

[cmbant]

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

[abonaldi]

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 @abonaldi<https://github.com/abonaldi> who perhaps can comment — You are receiving this because you were mentioned. Reply to this email directly, view it on GitHub<#34 (comment)>, or unsubscribe<https://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.

[cmbant]

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