pynbody.plot.stars.behroozi#
- pynbody.plot.stars.behroozi(xmasses, z, alpha=-1.412, Kravtsov=False)[source]#
Based on Behroozi+ (2013) return what stellar mass corresponds to the halo mass passed in.
Usage
>>> from pynbody.plot.stars import moster >>> xmasses = np.logspace(np.log10(min(totmasshalos)),1+np.log10(max(totmasshalos)),20) >>> ystarmasses, errors = moster(xmasses,halo_catalog._halos[1].properties['z']) >>> plt.fill_between(xmasses,np.array(ystarmasses)/np.array(errors), y2=np.array(ystarmasses)*np.array(errors), facecolor='#BBBBBB',color='#BBBBBB')