Constraints in number density per cubic parsec

[micheledoro]

I have to still have a look at the code, but I just ask, is there a way already coded to plot the constraints in number density per cubic parsec? In some cases this unit is useful and interesting too. Apologies if it was already there and I did not spot it.

[bradkav]

There isn't currently a way to plot the constraints in terms of number density per cubic parsec.

In principle, you could convert the bounds from f_PBH = Omega_PBH/Omega_DM, using n_PBH = f_PBH*<rho_DM>/M_PBH, with the mean cosmic DM density <rho_DM> ~ 3e-8 M_sun/pc^3 (if I've done the conversion right). However, I think this is probably misleading, because the constraints are generally local, so the constraint on the the number density of BHs would be proportional to the density of DM (in whatever region the constraint is derived from). Unfortunately, you would have to look up the relevant DM density for each constraints...

[micheledoro]

While it is probably true that constraints are local, following the evolution of the clustering of massive matter, I believe that even an average or effective number is interesting. All in all, as I understand it, even f_PBH is an average number over the sampled volume, isn't it?

PS: I asked in case you were intersted in updating the code to write different y-axis, all in all, that's ok with me, you can close this issue. Thanks