Open dkirkby opened 9 years ago
I was following the description in the joint BAO analysis of Aubourg et al., see equation 8 and the paragraph below.
It looks like the issue is the extra factor of g_c in Auborg et al to account for neutrino heating from e+e- annihilations, that is not present in Komatsu 2011. I seem to remember something about this in Turner & Kolb?
Can you check the (Nnu/3)^(1/4) factor used to calculate the neutrino temperature in massiveNuFunction[ ] ? I don't think it should be there, but probably removing it requires some compensating change somewhere else.