Added a new class in halo_occupation for Leauthaud HOD

Added a new Class in halo_occupation.py for stellar mass based Leauthaud HOD. Not throughly tested, but checked that Ncen and Nsat give reasonable values.
Added default stellar mass threshold to defaults.py