I have changed one of the derived quantities ([tex]r_{10}[/tex]) in the GetDist.f90 file, writing, in the MapParameters subroutine the following:
Code: Select all
invars(2+20)=(931400*invars(2+6)*invars(2+2))/(invars(2+21)*invars(2+21))
It is to define the neutrino mass as
[tex]
m_{\nu} = 93.14*f_{\nu}*\omega_{dm} (eV).
[/tex]
The results I got are sowing a double-peaked distribution for the neutrino mass [tex]m_{\nu}[/tex] and the fraction of DM in form of neutrinos [tex]f_{\nu}[/tex] (figure), and I think it should not be like that.
Someone suggested that there could be a problem with the priors, but I am no sure about that.
Do you have eny suggestion about this?
Thanks a lot,
MarioA.