CosmoMC
CosmoMC copied to clipboard
SZ cluster counts likelihood: natural vs base 10 log in szcounts.f90
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.
I think this code was contributed by @abonaldi who perhaps can comment
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.
Here on github. if you are looking at the github email, click "view it on github".